Loci associated with spontaneous abortion in primiparous Holstein cattle

Introduction: Spontaneous abortion (SA) in cattle refers to pregnancy loss occurring between days 42 and 260 of gestation. SA is costly, inefficient, and often leads to premature culling of cows.

Methods: This study aimed to identify loci associated with SA in primiparous Holstein cows by conducting a genome-wide association analysis of cows bred via artificial insemination (AI; 679 controls and 69 cases that aborted) or that were embryo transfer (ET) recipients (236 controls and 41 cases) from a single dairy.

Results: In AI-bred cows, 86 loci (413 single nucleotide polymorphisms or SNPs) and 168 positional candidate genes were associated (FDR < 0.05) with SA, while in ET recipients, 4 loci (10 SNPs) and 16 positional candidate genes were associated (FDR < 0.05) with SA. No SA-associated loci were shared between AI-bred or ET recipient cows, but nine associated loci (FDR < 0.05) in AI-bred cows were shared with AI-bred heifers.

Discussion: The difference in loci associated with AI-bred and ET recipient cattle may be due to differences in mechanisms associated with the maintenance of pregnancy between in vivo and in vitro derived embryos, or a larger sample size may be needed to identify additional shared loci. Identifying loci associated with SA in AI-bred and ET recipient cows provides an opportunity to enhance selection for reproductive efficiency in Holstein cattle.

1 Introduction

Spontaneous abortion (SA) or fetal loss typically occurs in 4 to 24% of pregnancies in primiparous cows (1–3). Spontaneous abortion results in increased culling of animals and reduced milk production, resulting in significant economic losses for the dairy industry. Keshavarzi et al. (4) found that SA reduced milk yield by 19.4% in cows that initiated lactation due to SA and by 7.3% in cows that experienced SA during lactation. Animals experiencing SA have an increased risk of retaining fetal membranes and other health issues, which may further reduce reproductive performance and increase costs to the industry (5–7).

Genomic selection for conception rate, daughter pregnancy rate, calving ability/ease (the ability of calves to be born unassisted), and calving interval (the time between the birth of two calves for the same cow) has been utilized to improve reproductive efficiency (8–10). In Holstein cattle, the conception rate improved by approximately 5% in heifers and 18% in cows between 2010 and 2022 (11, 12). Although the use of fertility traits in genomic selection has accelerated genetic progress for fertility, especially for embryonic loss, there is limited information regarding SA. Specifically, the genetic impact on SA risk in cows bred by artificial insemination (AI) compared to those that are embryo transfer (ET) recipients is not well understood. Recent reports have indicated the genetic effects of heifer genomics on SA for those bred by AI and those that were ET recipients, but there is limited research on SA in primiparous cows (13).

Some loci associated with fertility traits may differ between heifers and cows due to the distinct biological demands experienced by the two parities. For example, there are differences in the nutritional demands of a growing heifer compared to those of a lactating primiparous cow. The energy demand for heifers is focused on growth and puberty, while primiparous cattle must recover from parturition, lactate, and return to estrus. Milk production typically peaks between 60 and 100 days post-calving, coinciding with the time when cows are being bred (14). This creates a significantly greater nutritional demand on primiparous cows as they support their energy needs, the energy required for milk production, and, after conception, the energy demands of a developing fetus. These differences necessitate evaluating primiparous cows separately from heifers. To better understand how the use of ET affects SA compared to AI-bred cattle, the loci and positional candidate genes associated with SA in primiparous cows must be identified. Additionally, it is important to determine whether selecting for reduced fetal loss in AI-bred cattle can also improve calving rates in ET recipients by further evaluating these loci. Therefore, the objective of this study was to identify loci and genes associated with SA in US Holstein primiparous cows bred by AI or serving as ET recipients using a genome-wide association analysis (GWAA). These results will pinpoint loci that could be utilized in genomic selection indices and provide identification of specific genes associated with SA in cattle bred by AI or those that are ET recipients, leading to a better understanding of the fetal loss process.

2 Materials and methods

2.1 Study population and phenotypes

Health and breeding records for 5,750 primiparous Holstein cows were obtained using Dairy Comp 305 (Valley Agricultural Software, Tulare, CA, USA) records from a single dairy in southern Idaho. Data were recorded daily for breeding, production, and disease. The dairy consists of a milking herd of 2,300 cows housed in a dry lot with shade and fed a total mixed ration. Only cows that were pregnant to the first AI (n = 748) or were confirmed pregnant to the first ET (n = 277) were included in the study to eliminate possible confounding of phenotypes for embryonic loss. Cows were bred following a double ovulation synchronization protocol or after standing estrus prior to AI. There was no difference (p < 0.05, ANOVA) in the frequency of SA in cows bred by AI during observed estrus or following synchronization. Cows that were ET recipients received a double ovulation synchronization protocol with embryo transfer on day 7 post-estrus. Pregnancy was determined using trans-rectal ultrasound at 30 days post-breeding. Phenotypes were based on whether a cow maintained a pregnancy to ≥ 260 days of gestation (control) or spontaneously aborted between days 42 and 260 of gestation (case). Cases were identified by the cows returning to estrus following a confirmed pregnancy. AI controls were compared to AI cases, while ET controls were compared to ET cases. Four AI cows and three ET recipient cows were removed from the study as they experienced mastitis, metritis, metabolic issues, lameness, or respiratory disease between approximately day 30 before breeding and/or before calving. After filtering for these events, 679 AI controls, 69 AI cases, 236 ET controls and 41 ET cases remained in the study. All cows were at least 25 months of age at the beginning of their first lactation.

2.2 Genotyping and imputation

Genotypes from Zoetis (Kalamazoo, MI, USA) were imputed to a higher density utilizing Beagle (V. 4.1) and an in-house reference population of 4,800 US Holsteins genotyped with the Illumina (San Diego, CA, USA) BovineHD BeadChip (15). The 636,042 imputed single nucleotide polymorphisms (SNPs) had a call accuracy of >95% across the genome. Accuracy was determined by reducing the genotyping density of animals in the reference population to approximately 50,000 SNPs. Imputation was then performed on the lower density genotypes to restore them to the BovineHD BeadChip level for assessing the accuracy of imputation. Genotypes from the imputed SNPs were compared with the SNPs obtained from the initial genotyping using the BovineHD BeadChip.

2.3 Quality control

Prior to the GWAA, imputed genotypes were filtered for quality. SNPs were removed if they had a call rate <0.90, had a minor allele frequency (MAF) <0.01 or if they failed (p < 1 x 10⁻¹⁰⁰) Hardy–Weinberg equilibrium testing. Primiparous cows bred by AI had 4,733 SNPs removed for call rate, 140,423 SNPs were removed for MAF, and 1,012 SNPs were removed for failing Hardy–Weinberg equilibrium testing. For primiparous cows receiving ET 4,658 SNPs were removed for call rate, and 151,385 SNPs were removed for low MAF. No SNPs failed Hardy–Weinberg equilibrium testing for cows that were ET recipients, nor were any cows removed for genotyping call rate (<0.90) or for differences between genotypic and phenotypic sex. After quality control was complete, a total of 1,025 cows remained for the analyses, with 489,927 SNPs remaining for cows bred by AI and 479,927 SNPs remaining for ET recipient cows.

2.4 Genome-wide association analysis

A GWAA was performed using the SNP and Variation Suite (SVS) software version 8.1 (70). An Efficient Mixed Model Associated eXpedited (EMMAX) statistical approach was used with an identity by state matrix (16). The general mixed model for EMMAX is expressed as γ = Xβ+Zμ+ε where y = the vector of observed phenotypic values, X = a matrix of fixed effects, β = regression coefficients, Z = a matrix containing the observed random effects, u = vector of random effects related to variants of allele substitutions in the population, and ε = residual effects (16). Associations with SA were identified when FDR <0.05 (17, 18). Dosage compensation for SNPs on the X chromosome were not accounted for in EMMAX.

Due to the unknown inheritance of SA in cattle, three inheritance models (additive, dominant, and recessive) were analyzed for AI and ET groups. The additive model identified associations where having two minor alleles (aa) was twice as likely to impact the phenotype (SA) compared to having no minor alleles (AA), and half as likely to have an impact when one minor allele (Aa) was present. The dominant model identified associations where having one or two minor alleles (Aa or aa) was likely to have more impact on the phenotype compared to not having a minor allele (AA). A recessive model identified associations where having two minor alleles (aa) was likely to have more impact on the phenotype compared to having one or no minor allele (Aa or AA).

A principal component analysis was performed to evaluate population stratification. Two distinct clusters within the population based on birth year were identified (Supplementary Figure S1). An ANOVA (p < 0.05) was used to test if birth year, service sire, and season of conception affected SA. The tested covariates with the associated p-values for each population are listed in Table 1. To ensure that the population stratification was accounted for, the genomic inflation factor lambda (λ_GC) was calculated in SVS for each model. For AI-bred primiparous cows, no covariates were included as none of the tested covariates showed a significant (p < 0.05) effect on SA (Number of Sires used). In contrast, ET recipient primiparous cows had a difference (p < 0.05) in SA by season, therefore, this covariate was included in the analysis. Additional covariates, not identified by ANOVA as significant, were not included in the analysis to minimize false negative errors due to over-correction of possible population stratification.

Table 1

Table 1. P-value of tested covariates for artificial insemination (AI) bred and embryo transfer (ET) recipient primiparous cows.

If multiple SNPs were associated with SA on a chromosome, loci were defined by a threshold of D'> 0.7 (19, 20). For each locus, positional candidate genes were identified within the average haplotype size of 30.5 kb, calculated using the method of Gabriel (21), 5′ or 3′ to the associated SNP based on the bovine ARS-UCD 1.2 genome assembly (22).

The proportion of variance explained by a SNP was calculated in SVS utilizing the notation of Further Optimization when Covariates are Present, as described by Segura (72) and Vilhjalmsson (73). The sum of the proportion of variance explained for all SNPs will exceed 100% due to SNPs within a locus not being independent. This calculation was used to evaluate the percentage of total variation observed in spontaneous abortion that can be attributed to these SNPs.

A measure of the strength of the association between having the risk allele and experiencing a SA, referred to as relative risk, was calculated utilizing the frequencies of the minor and major alleles (23). The following equations were utilized for calculating relative risk for SA: a(a+b)/c(c+d), where a is the minor allele frequency in cases, b is the minor allele frequency in controls, c is the major allele frequency in cases, and d is the major allele frequency in controls.

Heritability was estimated for SA using a genomic best linear unbiased predictor (GBLUP) analysis (24, 25) with the average information algorithm (AI-REML), which is a bivariate restricted maximum likelihood analysis (26, 71). The AI-REML GBLUP method calculates variance components that are then used to calculate heritability, as the pseudo-heritability estimated by EMMAX can over-inflate the estimate in limited sample sizes.

Literature and database searches were conducted to identify any shared genomic regions or genes associated with SA and other fertility and production traits. These traits included the fertility index in Nordic Red Cattle, which encompasses multiple traits such as the number of inseminations per conception, duration of the interval from calving to first insemination, days from first to last insemination, and 46-day non-return rate (27). Other fertility traits included cow conception rate (percentage of cows that are pregnant at each service), daughter pregnancy rate (percentage of cows that become pregnant during each 21-day period), days open, and early embryonic loss (20, 28–31). Additionally, a comparison was made using tissue expression data based on fertility status and expression in reproductive tissues (32, 33). Any shared loci or positional candidate genes will be listed and discussed in the Discussion section.

Additionally, the comparison of production traits was conducted in two ways: first, through a literature search for shared loci and positional candidate genes related to production traits such as milk yield and milk components like fat yield, fat percentage, protein yield, and protein percentage (34). Secondly, multiple GWAA were utilized on the population of animals in this study that had complete lactation records for the following three production traits: total milk, total fat, and total protein.

3 Results

The combined rate of SA for first-service Holstein primiparous cows in AI-bred and ET recipient groups was 10.7%. Cows bred by AI experienced an SA rate of 9.2%, while cows that were ET recipients had an SA rate of 14.8%. The SA rate of 9.2% for cows bred by AI is consistent with other reports where the rate of SA ranged from 2 to 12% (3, 6). The SA rate of 14.8% for ET recipient cows was slightly higher than previous reports of 10 to 13% SA associated with in vitro-produced embryos (35, 36).

3.1 Loci associated with spontaneous abortion in primiparous cows bred by artificial insemination

The highest incidence of SA for primiparous cows bred by AI occurred between days 42 and 91 of gestation. The single highest day for SA was day 45 of gestation, with 12 cows aborting (Figure 1). These gestational days when SA was noted represent the day the cow returned to estrus after being identified as pregnant by trans-rectal ultrasound at approximately day 30 of gestation. The 54 SA occurrences between days 42 and 91 of gestation accounted for 78% of the fetal loss in this population. The second greatest interval of gestation associated with SA was between days 168 and 260, during which 20% of fetal loss occurred.

Figure 1

Figure 1. Timing of spontaneous abortions during the fetal period (days 42 to 260 of gestation) for primiparous cows bred by artificial insemination. The x-axis lists the gestation day while the y-axis indicates the number of spontaneous abortions occurring on each day of gestation during the fetal period.

The λ_GC for SA in AI cows was 0.99 for the additive inheritance model, 1.01 for the dominant inheritance model, and 0.96 for the recessive model. The estimated heritability for SA in primiparous cows bred by AI was 0.03 (±0.053). As the standard error for the heritability estimate is larger than the estimate itself, this provides low confidence in the estimated heritability. Since the majority of the loci associated with SA were recessive, and heritability measures additive genetic variance, it is unsurprising that the heritability estimate for SA is low. This estimated heritability for SA is similar to that reported for primiparous cows at 0.09 (±0.08) in a previous study (3) and for other reproductive traits such as the number of services, days from calving to first service, days open, and calving interval, which have been reported to be 0.040 ± 0.017, 0.034 ± 0.011, 0.053 ± 0.019, and 0.056 ± 0.014, respectively (3, 37). Rahbar et al. (38) also reported similar heritability estimates for several fertility traits, such as days open (0.016) and gestation length (0.123). Two other traits with reported heritability estimates comparable to that calculated in this study include the interval from first calving to first insemination (0.013) and the interval from first to successful insemination for primiparous cows (0.038) (39).

In the recessive model, there were 85 loci associated (FDR <0.05) with SA (Figure 2, Supplementary Table S1). The most significant loci associated with SA in the recessive model were on BTA3 at 110 Mb (FDR = 2.9 × 10⁻⁵), BTA22 at 3 Mb (FDR = 4.7 × 10⁻⁵), BTA24 at 37 Mb (FDR = 3.5 × 10⁻⁵), and BTA28 at 9 Mb (FDR = 1 × 10⁻⁴) and 39 Mb (FDR = 3.07 × 10⁻⁵). Positional candidate genes at these loci include: SFPQ, ZMYM4, KIAA0319L, RBMS3, DLGAP1, MTR, TRNAE-UUC, LOC112444209, LOC112446145, LOC618220, and LOC107131937. An additional 156 positional candidate genes were identified as associated with SA, for a total of 168 positional candidate genes associated with SA in primiparous cows bred by AI. The positional candidate genes were grouped into 16 functional groups as shown in Supplementary Table S2.

Figure 2

Figure 2. Loci associated with spontaneous abortion in Holstein primiparous cows that were bred by artificial insemination. The associations for the additive, dominant and recessive models of inheritance are represented in (A–C), respectively. All plots have the Bos taurus chromosome on the x-axis, and the -log10FDR-value on the y-axis. The red line indicates the threshold for an association (FDR <0.05) with spontaneous abortion. The Y chromosome is absent as all animals were female.

The highest proportion of phenotypic variance explained for the primiparous cows bred by AI was 5.6% from the locus on BTA28 at 9 Mb, with two positional candidate genes, MTR and TRNAE-UUC. Additionally, 21 loci representing 45 positional candidate genes, listed in Table 2, had a proportion of variance explained >3%. The average relative risk for the most significant SNPs associated with SA in primiparous cows bred by AI was 2.20 (Supplementary Table S3). The four positional candidate genes of the most significant loci are further described below.

Table 2

Table 2. Loci and positional candidate genes associated with spontaneous abortion in cows bred by artificial insemination with a proportion of variance explained >3%.

Positional candidate genes associated with SA in primiparous cows bred by AI were differentially expressed in the endometrium and in different cell types of the bovine placenta (Supplementary Figure S2) during the embryonic and fetal periods. Twenty-eight positional candidate genes were differentially expressed in high fertility compared to sub-fertile heifers at day 17 of gestation (32). An additional 42 positional candidate genes were expressed in the reproductive tissues of cattle, based on single-cell data during the embryonic and fetal periods, which support their involvement in pregnancy (Table 3; Supplementary Figure S2) (33).

Table 3

Table 3. Positional candidate genes for spontaneous abortion shared with fertility traits.

3.2 Association of spontaneous abortion in primiparous cows that were embryo transfer recipients

Primiparous cows that were ET recipients also experienced the greatest number of SA on day 45 of gestation (Figure 3). Seventy-six percent of the SA in ET recipients occurred during the first trimester.

Figure 3

Figure 3. Timing of spontaneous abortions during the fetal period (days 42 to 260 of gestation) for primiparous cows that were embryo transfer recipients. The x-axis lists the gestation day while the y-axis indicates the number of spontaneous abortions occuring on each day of gestation during the fetal period.

The λ_GC for SA in cows that were ET recipients was 1.02, 1.02, and 1.03 for the additive, dominant, and recessive inheritance models, respectively. No loci were associated with SA in cows that were ET recipients in the additive or dominant inheritance models, but four loci were associated (FDR <0.05) with SA in the recessive inheritance model on BTA9 and BTA13 (Figure 4, Table 4). There were 14 positional candidate genes within the two loci on BTA13 (Table 4). All four loci associated with SA in ET recipient primiparous cows had a proportion of variance explained >3%. The smallest variance explained was 8.7%, belonging to the locus on BT 13 at 46 Mb, while the largest proportion of variance explained was 10.5% for the locus on BTA13 at 47 Mb. The positional candidate genes for these loci can be seen listed in Table 4.

Figure 4

Figure 4. Loci associated with spontaneous abortion in Holstein primiparous cows that were embryo transfer recipients. The associations for the additive, dominant and recessive models of inheritance are represented in (A–C), respectively. All plots have the Bos taurus chromosome on the x-axis, and the -log10FDR-value on the y-axis. The red line indicates the threshold for an association (FDR <0.05) with spontaneous abortion. The Y chromosome is absent as all animals were female.

Table 4

Table 4. Loci associated with spontaneous abortions in cows that were embryo transfer recipients.

Loci and positional candidate genes associated with SA were also compared with production traits to identify possible correlations to determine if selection for improved fertility would negatively impact selection for production traits. The loci and positional candidate genes associated with SA in AI bred and ET recipient primiparous cows were not associated (p <0.05) with dairy production traits. These traits included total milk, total fat, total protein, milk yield, fat yield, fat percentage, protein yield and protein percentage as mentioned in the materials and methods. This indicates that selection for production traits while selecting for reduced SA could be done simultaneously without reducing genetic progress in production traits.

4 Discussion

In cattle, ET has been found to be associated with abnormal morphology of the embryo, abnormal placentation, and elevated immune responses compared with pregnancies established by natural service or AI (40–43). Studies in mice have found that embryos resulting from ET have increased placental weight, a reduced fetal-to-placental weight ratio, and significantly larger tissues when compared to naturally conceived controls (42). ET also increases epigenetic perturbations that could potentially lead to adverse neonatal and long-term health outcomes in offspring (42). In vitro-produced embryos have a more fragile zona pellucida, reduced intracellular communication, a higher incidence of chromosome abnormalities, errors of imprinting, a slower growth rate, higher thermal sensitivity, lower inner cell mass/trophectoderm ratio, and differences in gene expression compared to in vivo embryos (43). The differences between in vivo-derived embryos and in vitro-produced embryos could be reflected in the different genes and loci associated with SA.

The identification of loci and positional candidate genes provides added information to pursue in understanding how these genomic regions are essential for the maintenance of pregnancy. Genomic mechanisms that influence pregnancy can result from modifications to positional candidate genes that alter gene function or expression, changes in regulatory element sites, or changes in methylation or chromatin accessibility that affect transcription.

4.1 Comparison of loci associated with spontaneous abortion in cows bred by artificial insemination and previous fertility studies

One locus on BTA17 at 51 Mb was associated (FDR = 0.049) with SA in the additive model with the positional candidate gene nuclear receptor corepressor 2 (NCOR2) (Figure 2, Supplementary Table S1). The NCOR2 gene codes for a nuclear corepressor that interacts with multiple transcription factors. This corepressor orchestrates proliferation and differentiation during B cell development and is critical for fetal development (44). There were six loci associated with SA in cows bred by AI that were previously identified as associated with fertility traits. For a locus to be determined to be shared, the associated SNPs were required to be within a haplotype (31 kb) of each other. Two of these six loci, located on BTA1 and BTA4, were shared with genomic regions identified by Cole et al. (28) for daughter pregnancy rate (DPR) in Holstein cows. The sharing of loci between DPR and SA may be due to similar physiological processes necessary for the initiation and maintenance of pregnancy in cattle. A locus on BTA16 associated with cow conception rate to the first service and number of times bred in Holstein cattle was also associated with SA (31). A study by Höglund et al. (45) identified two loci (BTA2 at 46 Mb and BTA29 at 48 Mb) that were associated with fertility traits in Nordic cattle and shared with SA in this study. The locus on BTA2 was also shared with a genomic region identified for the Nordic female fertility index and days from first to last insemination in Nordic Holstein, Nordic Red, and Jersey cows (45). The Nordic female fertility index is a multi-trait index based on data collected on dairy cattle in Denmark, Sweden, and Finland and includes the number of inseminations required per pregnancy in cows and heifers, the length in days of the interval from calving to first insemination in cows, days from first to last insemination in cows and heifers, 56-day non-return rate in cows and heifers, and heat strength in cows and heifers (45). A locus on BTA12 at 52 Mb was also associated with a fertility index in Nordic Red cattle and with SA in Holstein cattle (27). Finally, a locus on BTA29 that was associated with SA was also associated with the number of inseminations per conception and days from first to last insemination in cows and heifers in Nordic Holstein, Nordic Red, and Jersey cows (45). Additionally, there were nine loci (BTA1 at 128 Mb, BTA5 at 85 Mb, BTA7 at 6 Mb, BTA7 at 85 Mb, two on BTA13 at 75 Mb, BTA17 at 35 Mb, BTA18 at 9 Mb, and BTA23 at 43 Mb) associated with SA in heifers bred by AI (13) and primiparous cows bred by AI. This indicates that shared genomic regions across parities could be used for genomic selection to help lower the risk and occurrences of SA.

There were 168 positional candidate genes, representing 86 loci, associated with SA in cows bred by AI, of which 68 were previously identified as having a role in fertility (Table 3). Two positional candidate genes (FST5 on BTA17 and CDH13 on BTA18) have been identified as positional candidate genes associated with SA in a separate population of Holstein cattle and were associated with embryonic loss (3, 20). The cadherin 13 (CDH13) gene, identified as a positional candidate gene for SA in AI bred cows, encodes a calcium-dependent cell adhesion, cadherin protein that protects vascular endothelial cells from apoptosis due to oxidative stress (46, 47). CDH13 plays a role in cell migration and acts as a transcription repressor as a tumor suppressor gene (48–50). CDH13 is correlated with trophoblastic invasion anomalies that lead to fetal growth restriction in humans, which can cause fetal demise (51, 52). Other studies have found that CDH13 plays a role in early reproductive development and cholesterol biosynthesis (46). There were 23 positional candidate genes for SA in AI bred cattle that have previously been identified as enriched for embryonic loss and one (EYA2) that was enriched for SA in Holstein cattle (3, 20, 31). An additional 14 genes (SOX5, LOC101904842, LOC100847115, EYA2, LOC101905203, ZMYND8, LOC104973934, LOC112449409, CDH13, PHACTR1, AP1M1, LOC112447625, FAM32A, and TRNAE-UUC) were reported as positional candidate genes associated with SA in Holstein heifers bred by AI (13).

The four positional candidate genes of the most significant loci include the splicing factor proline and glutamine-rich (SFPQ) gene, which encodes an RNA-binding protein that is critical for maintaining transcriptional elongation of long genes and enabling histone deacetylase binding activity. Histone deacetylases are enzymes that catalyze the removal of acetyl functional groups from lysine residues of both nonhistone and histone proteins (53). The removal of epigenetic modifications regulates chromatin structure and transcription, while the deacetylation of non-histones controls a variety of cellular processes (53).

The second positional candidate gene from the most significant loci for SA is the RNA binding motif single-stranded interacting protein 3 (RBMS3) gene, which is involved in DNA replication, gene transcription, cell cycle progression, and apoptosis (54). This gene influences the regulation of microRNA expression or stabilization, inhibits the Wnt/β-catenin pathway, and affects other epithelial–mesenchymal transition-related transcription factors (55).

The third positional candidate gene is DLG associated protein 1 (DLGAP1), which is a part of the DLGAP family of proteins that act as scaffold proteins in the brain (56). These proteins play a vital role in synaptic scaling by regulating the turnover of both metabotropic and ionotropic glutamate receptors in response to synaptic activity (56). DLGAP1 is specifically involved in megakaryocyte biology and platelet function, potentially giving a proliferative advantage in hematopoietic cells (57).

Lastly, the fourth positional candidate gene is LOC107131937, a pseudogene of the small ribosomal subunit protein uS12 that is conserved across all domains of life (58). The ribosomal protein S12 is a component of the decoding center of the 30S ribosomal subunit and is involved in transfer RNA selection (59). These transfer RNAs are key in synthesizing proteins, serving as an adapter or link between the messenger RNA (mRNA) molecule and the growing chain of amino acids that make up the protein (60). Pseudogenes can be either functional or non-functional; the status of this pseudogene is currently unknown. If it is functional, it could impact the synthesis of important proteins that affect SA. Holstein infertility haplotypes have also been discovered, which are significant due to their association with infertility, embryonic loss, and stillbirth in Holstein cattle. This makes it important to determine if any SA-associated loci found in this study may also correspond to several identified Holstein infertility (HH) haplotypes. Several Holstein infertility haplotypes (HH) have been identified and are currently being tested to reduce infertility. These haplotypes include HH1 to HH6 (61) and are located on BTA5 at 63 Mb (HH1), on BTA1 at 94 Mb (HH2), on BTA8 at 95 Mb (HH3), on BTA1 at 1.3 Mb (HH4), on BTA9 at 92 Mb (HH5), and on BTA16 at 29 Mb (HH6) (62, 63). The nearest locus associated with SA in primiparous cows bred by AI was approximately 1.2 Mb away from HH5 on BTA 9. The next closest loci were four loci on BTA 8 that ranged from 8 Mb to 22 Mb away from HH3. Additionally, there was one locus that was 19 Mb away from HH2 on BTA 1, one locus that was 13 Mb away from HH6 on BTA16, and another that was 22 Mb away from HH1 on BTA5. Any other loci associated with SA were over 22 Mb away from these identified Holstein haplotypes, indicating independent regions.

4.2 Comparison of loci associated with spontaneous abortion in cows bred as embryo transfer recipients and associated or enriched with other fertility traits

Sixteen positional candidate genes were associated with SA in cows that were ET recipients, with five (SLC23A2, IDI1, PRNP, RASSF2, and GTPB4) previously identified as having roles in fertility. Solute carrier family 23 member 2 (SLC23A2) encodes a sodium-dependent vitamin C transporter linked to spontaneous preterm birth in women, with research suggesting that genetic variants may increase the risk of preterm delivery (64). IDI1 was previously reported as a leading-edge gene enriched with embryonic loss in Holstein heifers (20). Ras association domain family member 2 (RASSF2) is a tumor suppressor gene that regulates Ras signaling (65). The GTP binding protein 4 (GTPBP4) encodes a small GTP-binding protein located in the nucleus and is involved in signal transduction (66). There were no loci identified as associated with SA in primiparous cows receiving ET that were on the same chromosome as any of the six Holstein infertility haplotypes mentioned earlier.

4.3 Comparison of loci associated with spontaneous abortion in cows that were embryo transfer recipients and are differentially expressed

Three genes (RASSF2, GTPBP4, and IDI1) have been previously reported to be differentially expressed in the endometrium of pregnant beef heifers compared to non-pregnant, high-fertility beef heifers (32). IDI1, PRNP, and LARP4B were expressed in the reproductive tissues of cattle with pregnancies generated via embryo transfer, thus supporting their role in the function of a normal pregnancy (33).

Research is ongoing to determine the genomic differences between SA in AI-bred and ET-recipient cattle. There was no overlap of loci or genes associated with SA in cows bred by AI with those associated with SA in ET-recipient cows. This lack of overlap could be due to insufficient power (power = 0.2) to detect an association or due to distinct physiological, biological, and genomic regions essential for pregnancy establishment and maintenance in AI-bred and ET recipients. Lower statistical power also increases the risk of false negatives in the analysis; therefore, further studies need to be conducted with larger sample sizes to confirm these results. To the author's knowledge, there are no other association studies that have investigated loci associated with fetal loss in ET-recipient cattle, but there have been studies investigating embryonic loss in ET recipients. Similar results have been identified in loci associated with embryonic loss in cows and heifers bred by AI and those that were ET recipients (67) and are consistent with physiological studies that have identified differences between AI-bred and ET-recipient pregnancies.

For primiparous cows bred by AI and those receiving ET, the greatest number of SA (as indicated by cows returning to estrus) occurred on day 45 of gestation. This could indicate possible issues or deformities of the placenta that lead to SA. Functional insufficiency of the placenta often results in increased offspring mortality (68). The gestational interval between days 28 and 60 is a pivotal period for the development of the placenta (69). This period of gestation requires substantial growth and differentiation of embryonic and trophoblastic cells, with significant communication between the dam and the embryo occurring (69). After day 30 of gestation, a firm connection between fetal cotyledons and maternal caruncles occurs, with the fetus developing a distinguishable adult organ system by day 60. Fetal loss during this period can result from developmental deficiencies such as inadequate placentation, difficulty transitioning from amniotic to allantoic nutrition, changes in placental vascularization, and underdevelopment of the embryo or fetus (69). Twenty-five percent of pregnancies from in vitro-produced embryos are reported to be lost between days 28 and 60 of gestation (69). In these embryos, there is a clear underlying physiology of placental abnormalities that contribute to the high incidence of pregnancy loss during this period.

Other factors can impact SA, such as infectious diseases, nutritional status, and overall management. In this study, we accounted for infectious diseases and some nutritional issues such as ketosis by excluding animals from the study that were recorded as experiencing these events during their pregnancy. However, not everything can be carefully managed or recorded in the dairy's records (such as body condition score), so there is potential for confounding variables to go undetected and potentially impact the occurrence of SA.

A future direction for the analysis of SA with a greater number of cases would be to split the gestational window into two categories to include early losses (approximately 42 to 120 days) and late losses (121–260 days). As reported in the results, there was a greater number of SA that occurred between days 42 and 120 for both primiparous cows bred by AI and those that were ET recipients. Further analysis focusing on SA occurring between days 42 and 120 could also be beneficial, as early losses are more often linked to genetic factors than late-term losses. A future study examining these early and late losses could be very helpful in identifying genetic regions and biological components that impact SA losses at different times during gestation.

5 Conclusions

This study identified loci and positional candidate genes associated with SA in Holstein primiparous cows bred by AI and those that were ET recipients. While physiological differences between AI and ET services have previously been identified, this study further supports the existence of physiological differences in the maintenance of pregnancies from embryos produced by cows bred by AI compared to those that received ET, as evidenced by the different genomic regions identified in the GWAA. None of the loci identified for cows bred by AI and those receiving ET overlapped with reported loci associated with dairy production traits. This suggests that SA and dairy production traits can be selected for simultaneously without a negative correlated response. By identifying the loci and positional candidate genes associated with SA in AI bred and ET recipient cows, this study enhances the understanding of the underlying factors that contribute to SA in cattle.

There were six loci associated with SA in AI bred cows and five loci associated with SA in ET recipients that have previously been identified as associated with either SA or embryonic loss in cattle. These shared regions with embryonic loss could be due to the similarity of pregnancy loss in the SAs occurring at day 45 of gestation for both primiparous cows bred by AI and those that were ET recipients. Seven loci were shared with fertility traits in previous GWAA studies, and 18 positional candidate genes were related to pregnancy during the fetal period. Sixty-seven positional candidate genes for SA were differentially expressed in reproductive tissues during the fetal stage of pregnancy, supporting their role in pregnancy initiation or maintenance.

The occurrence of spontaneous abortion in Holstein dairy cattle results in higher culling rates, longer calving intervals, and decreased milk production. This makes SA a substantial economic cost to producers each year. By utilizing identified genomic regions associated with SA in AI bred and ET recipients, this study provides a basis for improving reproductive efficiency through genomic selection.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://osf.io/qym4e/, doi: 10.17605/OSF.IO/QYM4E.

Ethics statement

The animal studies were approved by Washington State University Institutional Animal Use and Care Committee. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author contributions

ES: Data curation, Investigation, Validation, Writing – review & editing, Visualization, Software, Writing – original draft, Formal analysis. VK: Formal analysis, Validation, Software, Writing – review & editing. JK: Resources, Validation, Writing – review & editing, Data curation, Methodology, Software. KD: Visualization, Resources, Writing – review & editing, Data curation. BM: Conceptualization, Methodology, Writing – review & editing, Funding acquisition. AH: Resources, Data curation, Writing – review & editing. HN: Supervision, Writing – review & editing, Investigation, Writing – original draft, Conceptualization, Methodology, Funding acquisition, Resources, Project administration, Visualization, Formal analysis.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This project was supported by Agriculture and Food Research Initiative Competitive Grant No. 2023-67015-39524 and Hatch Project WNP0007 (Accession #1025787) from the US Department of Agriculture National Institute of Food and Agriculture.

Acknowledgments

The authors would like to acknowledge the Idaho dairy for their assistance in collecting of records for phenotyping, as well as Zoetis for the delivery of genotypes.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2025.1599401/full#supplementary-material

References

1. Humblot P. Use of pregnancy specific proteins and progesterone assays to monitor pregnancy and determine the timing, frequencies and sources of embryonic mortality in ruminants. Theriogenology. (2001) 56:1417–33. doi: 10.1016/S0093-691X(01)00644-6

Crossref Full Text | Google Scholar

2. Atkins JA, Smith MF, MacNeil MD, Jinks EM, Abreu FM, Alexander LJ, et al. Pregnancy establishment and maintenance in cattle1. J Anim Sci. (2013) 91:722–33. doi: 10.2527/jas.2012-5368

PubMed Abstract | Crossref Full Text | Google Scholar

3. Oliver KF, Wahl A, Dick M, Toenges JA, Kiser JN, Galliou JM, et al. Genomic analysis of spontaneous abortion in holstein heifers and primiparous cows. Genes. (2019) 10:954. doi: 10.3390/genes10120954

PubMed Abstract | Crossref Full Text | Google Scholar

4. Keshavarzi H, Sadeghi-Sefidmazgi A, Ghorbani GR, Kowsar R, Razmkabir M, Amer P. Effect of abortion on milk production, health, and reproductive performance of holstein dairy cattle. Anim Reprod Sci. (2020) 217:106458. doi: 10.1016/j.anireprosci.2020.106458

PubMed Abstract | Crossref Full Text | Google Scholar

5. Gilbert RO. Symposium review: mechanisms of disruption of fertility by infectious diseases of the reproductive tract. J Dairy Sci. (2019) 102:3754–65. doi: 10.3168/jds.2018-15602

PubMed Abstract | Crossref Full Text | Google Scholar

6. Sigdel A, Bisinotto RS, Peñagaricano F. Genes and pathways associated with pregnancy loss in dairy cattle. Sci Reports. (2021) 11:1. doi: 10.1038/s41598-021-92525-0

PubMed Abstract | Crossref Full Text | Google Scholar

7. Wilson DJ, Grohn YT, Bennett GJ, González RN, Schukken YH, Spatz J. Milk production change following clinical mastitis and reproductive performance compared among J5 vaccinated and control dairy cattle. J Dairy Sci. (2008) 91:3869–79. doi: 10.3168/jds.2008-1405

PubMed Abstract | Crossref Full Text | Google Scholar

8. Jiang J, Ma L, Prakapenka D, VanRaden PM, Cole JW, Da Y. A large-scale genome-wide association study in U.S. Holstein cattle. Front Genet. (2019) 10:412. doi: 10.3389/fgene.2019.00412

PubMed Abstract | Crossref Full Text | Google Scholar

9. Temesgen MY, Assen AA, Gizaw TT, Minalu BA, Mersha AY. Factors affecting calving to conception interval (days open) in dairy cows located at Dessie and Kombolcha towns, Ethiopia. PLoS ONE. (2022) 17:e0264029. doi: 10.1371/journal.pone.0264029

PubMed Abstract | Crossref Full Text | Google Scholar

10. Probo M, Guadagnini M, Sala G, Amodeo P, Bolli A. Calving ease risk factors and subsequent survival, fertility and milk production in Italian holstein cows. Animals. (2022) 12:671. doi: 10.3390/ani12060671

PubMed Abstract | Crossref Full Text | Google Scholar

11. Council on Dairy Cattle Breeding. In Genetic Trend: Cow Conception Rate for Holstein. Bowie, MD: CDCB. (2024). Available online at: https://webconnect.uscdcb.com/#/summary-stats/genetic-trend (accessed August 7, 2024).

Google Scholar

12. Council on Dairy Cattle Breeding. In Genetic Trend: Heifer Conception Rate for Holstein. Bowie, MD: CDCB. (2024). Available online at: https://webconnect.uscdcb.com/#/summary-stats/genetic-trend (accessed May 7, 2024).

Google Scholar

13. Suarez EM, Kelson VC, Kiser JN, Davenport KM, Murdoch BM, Neibergs HL. Genomic regions associated with spontaneous abortion in Holstein heifers. Genes. (2024) 15:1498. doi: 10.3390/genes15121498

PubMed Abstract | Crossref Full Text | Google Scholar

14. Hutjens MF. Dairy farm management systems | dry lot dairy cow breeds. Encyclop Dairy Sci. (2011) 2011:52–8. doi: 10.1016/B978-0-12-374407-4.00117-5

Crossref Full Text | Google Scholar

15. Browning BL, Browning SR. Genotype imputation with millions of reference samples. Am J Human Genet. (2016) 98:116–26. doi: 10.1016/j.ajhg.2015.11.020

PubMed Abstract | Crossref Full Text | Google Scholar

16. Kang HM, Sul JH, Service SK, Zaitlen NA, Kong S, Freimer NB, et al. Variance component model to account for sample structure in genome-wide association studies. Nat Genet. (2010) 42:348–54. doi: 10.1038/ng.548

PubMed Abstract | Crossref Full Text | Google Scholar

17. Otowa T, Maher BS, Aggen SH, McClay JL, van den Oord EJ, Hettema JM. Genome-wide and gene-based association studies of anxiety disorders in European and African American samples. PLoS ONE. (2014) 9:e112559. doi: 10.1371/journal.pone.0112559

PubMed Abstract | Crossref Full Text | Google Scholar

18. Genovese CR. False discovery rate control. Elsevier eBooks. (2015) 1:501–7. doi: 10.1016/B978-0-12-397025-1.00323-7

Crossref Full Text | Google Scholar

19. Weiss KM, Clark AG. Linkage disequilibrium and the mapping of complex human traits. Trends in Genetics. (2002) 18:19–24. doi: 10.1016/S0168-9525(01)02550-1

Crossref Full Text | Google Scholar

20. Galliou JM, Kiser JN, Oliver KF, Seabury CM, Moraes JGN, Burns GW, et al. Identification of loci and pathways associated with Heifer conception rate in U.S. Holsteins. Genes. (2020) 11:767. doi: 10.3390/genes11070767

PubMed Abstract | Crossref Full Text | Google Scholar

21. Gabriel SB. The structure of haplotype blocks in the human genome. Science. (2002) 296:2225–9. doi: 10.1126/science.1069424

PubMed Abstract | Crossref Full Text | Google Scholar

22. Rosen BD, Bickhart DM, Schnabel RD, Koren S, Elsik CG, Tseng E, et al. De Novo assembly of the cattle reference genome with single-molecule sequencing. GigaScience. (2020) 9:21. doi: 10.1093/gigascience/giaa021

PubMed Abstract | Crossref Full Text | Google Scholar

23. Ahsan H, Hodge SE, Heiman GA, Begg MD, Susser ES. Relative risk for genetic associations: the case-parent triad as a variant of case-cohort design. Int J Epidemiol. (2002) 31:669–78. doi: 10.1093/ije/31.3.669

PubMed Abstract | Crossref Full Text | Google Scholar

24. VanRaden PM. Efficient methods to compute genomic predictions. J Dairy Sci. (2008) 91:4414–23. doi: 10.3168/jds.2007-0980

PubMed Abstract | Crossref Full Text | Google Scholar

25. Taylor JF. Implementation and accuracy of genomic selection. Aquaculture. (2014) 420–421:S8–S14. doi: 10.1016/j.aquaculture.2013.02.017

Crossref Full Text | Google Scholar

26. Lee SH, Yang J, Goddard ME, Visscher PM, Wray NR. Estimation of pleiotropy between complex diseases using single-nucleotide polymorphism-derived genomic relationships and restricted maximum likelihood. Bioinformatics. (2012) 28:2540–2. doi: 10.1093/bioinformatics/bts474

PubMed Abstract | Crossref Full Text | Google Scholar

27. Höglund JK, Buitenhuis B, Guldbrandtsen B, Lund MS, Sahana G. Genome-wide association study for female fertility in nordic red cattle. BMC Genet. (2015) 16:11. doi: 10.1093/med/9780199398348.003.0011

Crossref Full Text | Google Scholar

28. Cole JB, Wiggans GR, Ma L, Sonstegard TS, Lawlor TJ, Crooker BA, et al. Genome-wide association analysis of thirty one production, health, reproduction and body conformation traits in contemporary U.S. Holstein cows. BMC Genom. (2011) 12:408. doi: 10.1186/1471-2164-12-408

PubMed Abstract | Crossref Full Text | Google Scholar

29. Cochran SD, Cole JB, Null DJ, Hansen PJ. Discovery of single nucleotide polymorphisms in candidate genes associated with fertility and production traits in Holstein cattle. BMC Genet. (2013) 14:49. doi: 10.1186/1471-2156-14-49

PubMed Abstract | Crossref Full Text | Google Scholar

30. Ortega MS, Denicol AC, Cole JB, Null DJ, Hansen PJ. Use of single nucleotide polymorphisms in candidate genes associated with daughter pregnancy rate for prediction of genetic merit for reproduction in Holstein Cows. Anim Genet. (2016) 47:288–97. doi: 10.1111/age.12420

PubMed Abstract | Crossref Full Text | Google Scholar

31. Kiser JN, Clancey E, Moraes JGN, Dalton J, Burns GW, Spencer TE, Neibergs HL. Identification of loci associated with conception rate in primiparous Holstein cows. BMC Genom. (2019) 20:2. doi: 10.1186/s12864-019-6203-2

PubMed Abstract | Crossref Full Text | Google Scholar

32. Moraes JGN, Behura SK, Geary TW, Hansen PJ, Neibergs HL, Spencer TE. Uterine influences on conceptus development in fertility-classified animals. Proc Nat Acad Sci. (2018) 115:E1749–58. doi: 10.1073/pnas.1721191115

PubMed Abstract | Crossref Full Text | Google Scholar

33. Davenport KM, O'Neil EV, Ortega MS, Patterson A, Kelleher AM, Warren WC, et al. Single-cell insights into development of the bovine placenta. Biol Reprod. (2024) 110:169–84. doi: 10.1093/biolre/ioad123

PubMed Abstract | Crossref Full Text | Google Scholar

34. Bekele R, Taye M, Abebe G, Meseret S. Genomic regions and candidate genes associated with milk production traits in Holstein and its crossbred cattle: a review. Int J Genomics. (2023) 2023:1–18. doi: 10.1155/2023/8497453

PubMed Abstract | Crossref Full Text | Google Scholar

35. Kruip AM, den Daas JHG. In vitro produced and cloned embryos: effects on pregnancy, parturition and offspring. Theriogenology. (1997) 47:43–52. doi: 10.1016/S0093-691X(96)00338-X

Crossref Full Text | Google Scholar

36. Hasler JF. In-vitro production of cattle embryos: problems with pregnancies and parturition. Hum Reprod. (2000) 15:47–58. doi: 10.1093/humrep/15.suppl_5.47

PubMed Abstract | Crossref Full Text | Google Scholar

37. Guo G, Guo X, Wang Y, Zhang X, Zhang S, Li X, et al. Estimation of genetic parameters of fertility traits in chinese holstein cattle. Can J Anim Sci. (2014) 94:281–5. doi: 10.4141/cjas2013-113

Crossref Full Text | Google Scholar

38. Rahbar R, Aminafshar M, Abdullahpour R, Chamani M. Genetic analysis of fertility traits of Holstein dairy cattle in warm and temperate climate. Acta Scientiarum Animal Sci. (2016) 38:333. doi: 10.4025/actascianimsci.v38i3.31377

Crossref Full Text | Google Scholar

39. Otwinowska-Mindur A, Ptak E, Jagusiak W, Zarnecki A. Estimation of genetic parameters for female fertility traits in the polish holstein-friesian population. Animals. (2022) 12:1485. doi: 10.3390/ani12121485

PubMed Abstract | Crossref Full Text | Google Scholar

40. Demetrio DGB, Santos RM, Demetrio CGB, Vasconcelos JLM. Factors affecting conception rates following artificial insemination or embryo transfer in lactating Holstein cows. J Dairy Sci. (2007) 90:5073–82. doi: 10.3168/jds.2007-0223

PubMed Abstract | Crossref Full Text | Google Scholar

41. Sartori R, Gümen A, Guenther JN, Souza AH, Caraviello DZ, Wiltbank MC. Comparison of artificial insemination versus embryo transfer in lactating dairy cows. Theriogenology. (2006) 65:1311–21. doi: 10.1016/j.theriogenology.2005.05.055

PubMed Abstract | Crossref Full Text | Google Scholar

42. de Waal E, Vrooman LA, Fischer E, Ord T, Mainigi M, Coutifaris C, et al. The cumulative effect of assisted reproduction procedures on placental development and epigenetic perturbations in a mouse model. Human Mol Genet. (2015) 24(24), ddv400–ddv400. doi: 10.1093/hmg/ddv400

PubMed Abstract | Crossref Full Text | Google Scholar

43. Van Eetvelde M, Heras S, Leroy JLMR, Van Soom A, Opsomer G. The importance of the periconception period: immediate effects in cattle breeding and in assisted reproduction such as artificial insemination and embryo transfer. Adv Exp Med Biol. (2017) 1014:41–68. doi: 10.1007/978-3-319-62414-3_3

PubMed Abstract | Crossref Full Text | Google Scholar

44. Lee RD, Knutson TP, Munro SA, Miller JT, Heltemes-Harris LM, Mullighan CG, et al. Nuclear corepressors NCOR1/NCOR2 regulate B cell development, maintain genomic integrity and prevent transformation. Nat Immunol. (2022) 23:1763–76. doi: 10.1038/s41590-022-01343-7

PubMed Abstract | Crossref Full Text | Google Scholar

45. Höglund JK, Sahana G, Guldbrandtsen B, Lund MS. Validation of associations for female fertility traits in nordic holstein, nordic red and jersey dairy cattle. BMC Genet. (2014) 15:8. doi: 10.1186/1471-2156-15-8

PubMed Abstract | Crossref Full Text | Google Scholar

46. Coen S, Keogh K, Lonergan P, Fair S, Kenny DA. Early life nutrition affects the molecular ontogeny of testicular development in the young bull calf. Sci Reports. (2023) 13:1. doi: 10.1038/s41598-022-23743-3

PubMed Abstract | Crossref Full Text | Google Scholar

47. Gene (2024). Gene. Bethesda (MD): National Library of Medicine (US), National Centerfor Biotechnology Information (2004). Available online at: https://www.ncbi.nlm.nih.gov/gene/ (accessed March 25, 2024).

Google Scholar

48. Chen F, Huang T, Ren Y, Wei J, Lou Z, Wang X, et al. Clinical significance of CDH13 promoter methylation as a biomarker for bladder cancer: a meta-analysis. BMC Urol. (2016) 16:1. doi: 10.1186/s12894-016-0171-5

PubMed Abstract | Crossref Full Text | Google Scholar

49. Killen AC, Barber M, Paulin JJW, Ranscht B, Parnavelas JG, Andrews W. Protective role of cadherin 13 in interneuron development. Brain Struct Funct. (2017) 222:3567–85. doi: 10.1007/s00429-017-1418-y

PubMed Abstract | Crossref Full Text | Google Scholar

50. Ziegler GC, Ehlis A-C, Weber H, Vitale MR, Zöller JEM, Ku H-P, et al. A common CDH13 variant is associated with low agreeableness and neural responses to working memory tasks in ADHD. Genes. (2021) 12:1356. doi: 10.3390/genes12091356

PubMed Abstract | Crossref Full Text | Google Scholar

51. Teng MS, Hsu LA, Wu S, Sun Y, Juan SH, Ko YL. Association of CDH13 genotypes/haplotypes with circulating adiponectin levels, metabolic syndrome, and related metabolic phenotypes: the role of the suppression effect. PLoS ONE. (2015) 10(4):e0122664. doi: 10.1371/journal.pone.0122664

PubMed Abstract | Crossref Full Text | Google Scholar

52. Biyik I, Metineren H, Ozturk E, Simsek S, Oztas E, Guldur ME, et al. Placental T-cadherin correlates with trophoblastic invasion anomalies: placenta percreta and fetal growth restriction. Int J Gynecol Patho. (2022) 42:293–300. doi: 10.1097/PGP.0000000000000884

PubMed Abstract | Crossref Full Text | Google Scholar

53. Seto E, Yoshida M. Erasers of histone acetylation: the histone deacetylase enzymes. Cold Spring Harb Perspect Biol. (2014) 6:a018713–a018713. doi: 10.1101/cshperspect.a018713

PubMed Abstract | Crossref Full Text | Google Scholar

54. Sukhija N, Malik AA, Devadasan JM, Dash A, Bidyalaxmi K, Ravi Kumar D, et al. Genome-wide selection signatures address trait specific candidate genes in cattle indigenous to arid regions of india. Animal Biotechnol. (2023) 35:2290521. doi: 10.1080/10495398.2023.2290521

PubMed Abstract | Crossref Full Text | Google Scholar

55. Górnicki T, Lambrinow J, Mrozowska M, Podhorska-Okołów M, Dziegiel P, Grzegrzółka J. Role of RBMS3 novel potential regulator of the EMT phenomenon in physiological and pathological processes. Int J Mol Sci. (2022) 23:10875–10875. doi: 10.3390/ijms231810875

PubMed Abstract | Crossref Full Text | Google Scholar

56. Rasmussen AH, Rasmussen HB, Silahtaroglu A. The DLGAP family: neuronal expression, function and role in brain disorders. Mol Brain. (2017) 10:9. doi: 10.1186/s13041-017-0324-9

PubMed Abstract | Crossref Full Text | Google Scholar

57. Kwiatkowski BA, Burwick NR, Richard RE. DLGAP1 directs megakaryocytic growth and differentiation in an mpl dependent manner in hematopoietic cells. Biomarker Res. (2019) 7:13. doi: 10.1186/s40364-019-0165-z

PubMed Abstract | Crossref Full Text | Google Scholar

58. Datta M, Pillai M, Modak MJ, Liiv A, Khaja FT, Hussain T, et al. A mutation in the ribosomal protein uS12 reveals novel functions of its universally conserved PNSA loop. Mol Microbiol. (2021) 115:1292–308. doi: 10.1111/mmi.14675

PubMed Abstract | Crossref Full Text | Google Scholar

59. Vila-Sanjurjo A, Lu Y, Aragonez JL, Starkweather RE, Sasikumar M, O'Connor M. Modulation of 16S rRNA function by ribosomal protein S12. Gene Struct Expres. (2007) 1769:462–71. doi: 10.1016/j.bbaexp.2007.04.004

PubMed Abstract | Crossref Full Text | Google Scholar

60. Raina M, Ibba M. TRNAs as regulators of biological processes. Front Genet. (2014) 5. doi: 10.3389/fgene.2014.00171

PubMed Abstract | Crossref Full Text | Google Scholar

61. Council of Dairy Cattle Breeding. Haplotypes & Genetic conditions. Bowie, MD: CDCB. (2024). Available online at: https://uscdcb.com/haplotypes/ (accessed May 7, 2024).

Google Scholar

62. Cole JB, Null DJ, VanRaden PM. Phenotypic and genetic effects of recessive haplotypes on yield, longevity, and fertility. J Dairy Sci. (2016) 99:7274–88. doi: 10.3168/jds.2015-10777

PubMed Abstract | Crossref Full Text | Google Scholar

63. Khan YA, Omar AI, He Y, Chen S, Zhang S, Xiao W, et al. Prevalence of nine genetic defects in Chinese Holstein cattle. Vet Med Sci. (2021) 7:1728–35. doi: 10.1002/vms3.525

PubMed Abstract | Crossref Full Text | Google Scholar

64. Erichsen HC, Engel SA, Eck PK, Welch R, Yeager M, Levine M, et al. Genetic variation in the sodium-dependent vitamin C transporters, SLC23A1, and SLC23A2 and risk for preterm delivery. Am J Epidemiol. (2005) 163:245–54. doi: 10.1093/aje/kwj035

PubMed Abstract | Crossref Full Text | Google Scholar

65. Akino K, Toyota M, Suzuki H, Mita H, Sasaki Y, Ohe-Toyota M, et al. The Ras effector RASSF2 is a novel tumor-suppressor gene in human colorectal cancer. Gastroenterology. (2005) 129:156–69. doi: 10.1053/j.gastro.2005.03.051

PubMed Abstract | Crossref Full Text | Google Scholar

66. Wu J, Chen G, Wang W, Yang Y, Yuan Y, Shang A, et al. GTPBP4: a new therapeutic target gene promotes tumor progression in non-small cell lung cancer via EMT. J Oncol. (2022) 2022:1–11. doi: 10.1155/2022/2164897

PubMed Abstract | Crossref Full Text | Google Scholar

67. Kelson VC, Kiser JN, Davenport KM, Suarez EM, Murdoch BM, Neibergs HL. Identifying regions of the genome associated with conception rate to the first service in holstein heifers bred by artificial insemination and as embryo transfer recipients. Genes. (2024) 15:765–765. doi: 10.3390/genes15060765

PubMed Abstract | Crossref Full Text | Google Scholar

68. Tanner AR, Kennedy VC, Lynch CS, Hord TK, Winger QA, Rozance PJ, et al. (2022). In vivo investigation of ruminant placenta function and physiology—a review. J Animal Sci. 100:skac045. doi: 10.1093/jas/skac045

PubMed Abstract | Crossref Full Text | Google Scholar

69. Wiltbank MC, Baez GM, Garcia-Guerra A, Toledo MZ, Monteiro PLJ, Melo LF, et al. Pivotal periods for pregnancy loss during the first trimester of gestation in lactating dairy cows. Theriogenology. (2016) 86:239–53. doi: 10.1016/j.theriogenology.2016.04.037

PubMed Abstract | Crossref Full Text | Google Scholar

70. SNP and Variation Suite. SNP and Variation Suite (Version 9.1) [SVS]. Bozeman MT: Golden Helix, INC. (n.d.). Available online at: http://www.goldenhelix.com/products/SNP_Variation/index.html%E2%80%AF%E2%80%AF (accessed 15 December, 2023).

Google Scholar

71. Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: a tool for genome-wide complex trait analysis. Am J Human Genet. (2011) 88:76–82. doi: 10.1016/j.ajhg.2010.11.011

PubMed Abstract | Crossref Full Text | Google Scholar

72. Segura V, Vilhjalmsson BJ, Platt A, Korte A, Seren U, Long Q, et al. An efficient multi-locus mixed-model approach for genome-wide association studies in structured populations. Nat. Genet. (2012) 44:825–30. doi: 10.1038/ng.2314

PubMed Abstract | Crossref Full Text | Google Scholar

73. Vilhjalmsson BJ, Nordborg M. The nature of confounding in genome-wide association studies. Nat. Rev. Genet. (2013) 14:1–2. doi: 10.1038/nrg3382

PubMed Abstract | Crossref Full Text | Google Scholar