Identification of genetic variants and individual genes associated with postpartum hypocalcemia in Holstein cows.

Periparturient hypocalcemia is a complex metabolic disorder that occurs at the onset of lactation because of a sudden irreversible loss of Ca incorporated into colostrum and milk. Some cows are unable to quickly adapt to this demand and succumb to clinical hypocalcemia, commonly known as milk fever, whereas a larger proportion of cows develop subclinical hypocalcemia. The main goal of this study was to identify causative mutations and candidate genes affecting postpartum blood calcium concentration in Holstein cows. Data consisted of blood calcium concentration measured in 2513 Holstein cows on the first three days after parturition. All cows had genotypic information for 79 k SNP markers. Two consecutive rounds of imputation were performed: first, the 2513 Holstein cows were imputed from 79 k to 312 k SNP markers. This imputation was performed using a reference set of 17,131 proven Holstein bulls with 312 k SNP markers. Then, the 2513 Holstein cows were imputed from 312 k markers to whole-genome sequence data. This second round of imputation used 179 Holstein animals from the 1000 Bulls Genome Project as a reference set. Three alternative phenotypes were evaluated: (1) total calcium concentration in the first 24 h postpartum, (2) total calcium concentration in the first 72 h postpartum calculated as the area under the curve; and (3) the recovery of total calcium concentration calculated as the difference in total calcium concentration between 72 and 24 h. The identification of genetic variants associated with these traits was performed using a two-step mixed model-based approach implemented in the R package MixABEL. The most significant variants were located within or near genes involved in calcium homeostasis and vitamin D transport (GC), calcium and potassium channels (JPH3 and KCNK13), energy and lipid metabolism (CA5A, PRORP, and SREBP1), and immune response (IL12RB2 and CXCL8), among other functions. This work provides the foundation for the development of novel breeding and management tools for reducing the incidence of periparturient hypocalcemia in dairy cattle.

Degree of leukochimerism with calf age in Angus-crossbred twin sets.

Twinning in cattle is infrequent and usually undesired. It can result in an increased occurrence of abortion and dystocia, reduced calf survival and a high likelihood of freemartinism in mixed-sex twins. Twin gestations are also commonly associated with the formation of placental vascular anastomoses (PVA) between twins. Through PVA they share blood, hormones (leading to freemartinism in mixed sex twins) and hematopoietic stem cells, which are the progenitors of white blood cells. The sharing of stem cells between twins can result in leukochimeric twin sets. These are twins that have white blood cells derived from both self and the co-twin owing to the fetal migration of hematopoietic stem cells from the extraembryonic mesoderm of the yolk sac to final sites like bone marrow and thymus. This study examined the degree to which this leukochimerism changes with age. DNA was extracted from hair bulbs containing mesenchymal dermal papilla to determine the individual's true genotype and blood samples were obtained at six time points from 1 week to 8 months of age to assess leukochimerism. Samples were genotyped using a medium density SNP chip, and quantitative estimates of allele frequency were determined using SNPs for which members of a twin set had alternative homozygous genotypes. The results indicate statistically significant changes in the proportion of self and co-twin with age and suggest that by 2-4 months of age the genotypic mix in white blood cells represents the hematopoetic stem cell population resident in the individual (i.e. permanently found in thymus and bone marrow).

Preantral follicle numbers and size in heifers carrying Trio, a bovine high fecundity allele.

In brief: The bovine high fecundity allele, Trio, results in the occurrence of multiple ovulations and is characterized by antral follicles that develop slower and acquire ovulatory capacity at smaller sizes. This study provides novel information on the effect of the Trio allele on early folliculogenesis. Abstract: The bovine high fecundity allele, Trio, causes overexpression in granulosa cells (GCs) of SMAD6, an inhibitor of BMP15-activated SMAD signalling. Furthermore, the Trio allele results in antral follicles that develop slower, acquire ovulatory capacity at smaller sizes, and have three-fold greater ovulation rate compared to half-sib non-carriers. The present study was designed to determine preantral follicle numbers and size in Trio carrier and non-carrier cattle testing the hypothesis that inhibition of SMAD signalling would alter preantral follicle activation and/or growth. Ovarian tissues from Trio carrier (n = 12) and non-carrier (n = 12) heifers were obtained by laparotomy after follicle wave synchronization. Follicle numbers and dimensions were determined for each stage of development (primordial, transitional, primary, and secondary) from paraffin-embedded sections. There were no differences in the number of primordial, transitional, or secondary follicles or in antral follicle count, circulating AMH, or ovarian volume between carriers and non-carriers. Trio carriers had ~2.5-fold greater (P < 0.01) number of primary follicles than non-carriers, and transitional and primary follicles were larger (~1.2-fold; P < 0.1) in Trio carriers. Oocyte volume of primordial and transitional follicles was not different between genotypes; however, oocytes were larger (P < 0.05) in primary (~1.3-fold) and secondary (~1.8-fold) follicles for Trio carriers. Granulosa cell numbers were not different (P > 0.3) between carriers and non-carriers, irrespective of the stage of development. These results suggest that, after primordial follicle activation, follicles in Trio carrier cattle have slower progression through the primary stage, hence the larger oocyte and greater number of primary follicles.

SMAD6 inhibits granulosa cell proliferation and follicle growth rate in carrier and noncarrier heifers of the Trio allele.

In brief: Follicle selection is a key event in monovular species. In this manuscript, we demonstrate the role of SMAD6 in promoting decreased granulosa cell proliferation and follicle growth rate in carriers vs noncarriers of the Trio allele and after vs before follicle deviation. Abstract: Cattle are generally considered a monovular species; however, recently, a bovine high fecundity allele, termed the Trio allele, was discovered. Carriers of Trio have an elevated ovulation rate (3-5), while half-sibling noncarriers are monovular. Carriers of the Trio allele have overexpression in granulosa cells of SMAD6, an inhibitor of oocyte-derived regulators of granulosa cell proliferation and differentiation. In experiment 1, follicle size was tracked for each follicle during a follicular wave. Follicle growth rate was greater before vs after follicle deviation in both carriers and noncarriers. Additionally, follicle growth rate was consistently less in carriers vs noncarriers. In experiment 2, we collected granulosa cells from follicles before and after deviation for evaluation of granulosa cell gene expression. Granulosa cell proliferation was less in carriers vs noncarriers and after vs before follicle deviation (decreased expression of cell cycle genes CCNB1 and CCNA2). The decreased granulosa cell proliferation in noncarriers after deviation was associated with increased SMAD6 expression. Similarly, in experiment 3, decreased expression of SMAD6 in granulosa cells of noncarriers cultured in vitro for 60 h was associated with increased expression of cell cycle genes. This suggests that SMAD6 may not just be inhibiting follicle growth rate in carriers of Trio but may also play a role in the decreased follicle growth after deviation in noncarriers. The hypotheses were supported that (1) follicle growth and granulosa cell proliferation decrease after deviation in both carriers and noncarriers and that (2) granulosa cell proliferation is reduced in carriers compared to noncarriers.

Identifying genetic variants and pathways influencing daughter averages for twinning in North American Holstein cattle and evaluating the potential for genomic selection.

Multiple birth in dairy cattle is a detrimental trait both economically for producers and for animal health. Genetics of twinning is complex and has led to several quantitative trait loci regions being associated with increased twinning. To identify variants associated with this trait, calving records from 2 time periods were used to estimate daughter averages for twinning for Holstein bulls. Multiple analyses were conducted and compared including GWAS, genomic prediction, and gene set enrichment analysis for pathway detection. Although pathway analysis did not yield many congruent pathways of interest between data sets, it did indicate two of interest. Both pathways have ties to the strong candidate region on BTA11 from the genome-wide association analysis across data sets. This region does not overlap with previously identified quantitative trait loci regions for twinning or ovulation rate in cattle. The strongest associated SNPs were upstream from 2 candidate genes LHCGR and FSHR, which are involved in folliculogenesis. Genomic prediction showed a moderate correlation accuracy (0.43) when predicting genomic breeding values for bulls with estimates from calving records from 2010 to 2016. Future analysis of the region on BTA11 and the relation of the candidate genes could improve this accuracy.

Selection of fewer dominant follicles in Trio carriers given GnRH antagonist and luteinizing hormone action replaced by nonpulsatile human chorionic gonadotropin†.

Studying selection of multiple dominant follicles (DFs) in monovulatory species can advance our understanding of mechanisms regulating selection of single or multiple DFs. Carriers of the bovine high fecundity Trio allele select multiple DFs, whereas half-sib noncarriers select a single DF. This study compared follicle selection during endogenous gonadotropin pulses versus during ablation of pulses with Acyline (GnRH antagonist) and luteinizing hormone (LH) action replaced with nonpulsatile human chorionic gonadotropin (hCG) treatment in Trio carriers (n = 28) versus noncarriers (n = 32). On Day 1.5 (Day 0 = ovulation), heifers were randomized: (1) Control, untreated; (2) Acyline, two i.m. doses (Days 1.5 and D3) of 3 µg/kg; (3) hCG, single i.m. dose of 50 IU hCG on Day 1.5 followed by daily doses of 100 IU; and (4) Acyline + hCG. Treatments with nonpulsatile hCG were designed to replace LH action in heifers treated with Acyline. Acyline treatment resulted in cessation of follicle growth on Day 3 with smaller (P < 0.0001) maximum follicle diameter in Trio carriers (6.6 ± 0.2 mm) than noncarriers (8.7 ± 0.4 mm). Replacement of LH action (hCG) reestablished follicle diameter deviation and maximum diameter of DFs in both genotypes (8.9 ± 0.3 mm and 13.1 ± 0.5 mm; P < 0.0001). Circulating follicle stimulating hormone (FSH) was greater in Acyline-treated than in controls. Finally, Acyline + hCG decreased (P < 0.0001) the number of DFs from 2.7 ± 0.2 to 1.3 ± 0.2 in Trio carriers, with most heifers having only one DF. This demonstrates the necessity for LH in acquisition of dominance in Trio carriers (~6.5 mm) and noncarriers (~8.5 mm) and provides evidence for a role of GnRH-induced FSH/LH pulses in selection of multiple DFs in Trio carriers and possibly other physiologic situations with increased ovulation rate.

Practical application of an impractical bovine genotype: creating bilateral twin pregnancies in Trio allele carriers.

Bovine twin birth is associated with detriments, including increased embryo/fetal losses, malpresentation, and dystocia. Incidence of these is lessened in bilateral compared with unilateral twin pregnancy. This study was undertaken to assess the use of follicular ablation by aspiration to create bilateral twin pregnancies in females with genetic potential for ~3.5 ovulations per cycle (Trio allele carriers). In experiment 1, carriers (n = 30) and noncarriers (n = 10) were synchronized for ovulation and timed artificial insemination (TAI). Follicles (>5 mm) in excess of one per ovary were aspirated ~16 h preceding TAI. Follicle count for females with follicles on only one ovary was reduced to two. Blood was sampled 2 wk post-TAI to assess progesterone (P4) concentrations; embryo count was determined by ultrasound 6 wk post-TAI. Circulating P4 concentration post-TAI was significantly (P < 0.001) associated with both genotype and subsequent pregnancy status (pregnant noncarriers: 7.06 ± 0.68 ng/mL; pregnant carriers: 5.54 ± 0.55 ng/mL; nonpregnant noncarriers: 5.22 ± 1.05 ng/mL; nonpregnant carriers: 3.13 ± 0.42 ng/mL). Experiment 2 was undertaken to offset the negative effects of follicular aspiration on subsequent P4 concentration observed in experiment 1. Carriers (n = 38) and noncarriers (n = 32) were submitted to TAI and follicle ablation as described for experiment 1. Additionally, accessory corpora lutea (CL) were induced in carriers by the administration of human chorionic gonadotropin (carriers) at day 6 post-TAI. Consequently, P4 concentration post-TAI was significantly (P < 0.05) associated with subsequent pregnancy status (pregnant: 8.48 ± 0.61 ng/mL; nonpregnant: 6.70 ± 0.63 ng/mL) but not with genotype (carrier: 8.01 ± 0.59 ng/mL; noncarrier: 7.17 ± 0.64 ng/mL). Embryo number was greater in carriers (exp. 1: 1.64 ± 0.81; exp 2: 1.45 ± 0.09) vs. noncarriers (1.00 ± 0.00, both experiments). Single, twin, and triplet pregnancies occurred in carriers in experiment 1, whereas multiples in experiment 2 were limited to twin pregnancies. Genotype effects on pregnancy rate were not significant (P > 0.10) in either experiment. Results suggest that follicular ablation to create bilateral twin pregnancies in Trio carriers is feasible but requires the induction of accessory CL to offset the negative effects of follicular aspiration on subsequent P4 concentration and associated fertility outcomes.

Genomic analysis of bovine respiratory disease and lung consolidation in preweaned Holstein calves using clinical scoring and lung ultrasound.

Bovine respiratory disease (BRD) is a leading cause of morbidity and mortality in dairy calves, with detrimental long-term effects that include stunted growth, increased age at first calving, and decreased milk production in first lactation. The objectives of this study were to establish a protocol for objective and efficient assessment of BRD phenotypes in preweaned dairy calves, develop a genomic reference population with well-defined clinical and subclinical phenotypes, identify chromosomal regions associated with BRD in a genome-wide association study, estimate genetic parameters of BRD, and predict genomic breeding values of dairy calves. A total of 1,107 Holstein calves from 6 dairy farms in southern Wisconsin were examined using clinical respiratory scoring and lung ultrasound at 3 and 6 wk of age. The clinical respiratory score was based on visual appraisal of eyes, nose, ears, cough, and temperature. Lung ultrasound scores were assigned based on the amount of consolidation present. Calves were genotyped with a commercially available SNP array and after quality control and imputation to higher density, 690,291 SNP markers and 1,014 individuals remained. Single-step genome-wide association study and single-step genomic best linear unbiased prediction were applied to binary 3- and 6-wk phenotypes considered as overall respiratory healthy versus affected (RESP3, RESP6) or as presence or absence of lung consolidation (CON3, CON6). Lung ultrasound combined with a clinical scoring system allowed for efficient and objective assessment for the prevalence of BRD. Proportions of variance attributed to 1-Mb non-overlapping windows suggested genomic regions that may contain putative candidate genes, most notably regions on Bos taurus autosomes 1, 6, 7, 10, 11, 12, 15, 17, 18, 27, and 28 that explained 0.70 to 1.45% of the genetic variance. Heritability estimates were higher at 3 wk (0.214 and 0.241 for CON3 and RESP3, respectively) than 6 wk (0.084 and 0.111 for CON6 and RESP6, respectively), and mean reliabilities of genomic estimated breeding vales for calves with genotypes and phenotypes ranged from 0.12 for CON6 to 0.30 for RESP3.

Short communication: Heritability of susceptibility to infection by Mycobacterium avium ssp. paratuberculosis in Holstein cattle.

Johne's disease in cattle is the result of infection of the small intestine by Mycobacterium avium ssp. paratuberculosis (MAP), leading to an incurable inflammatory bowel disease (Johne's disease or paratuberculosis). The disease is a concern both for its direct cost to dairy producers and for its zoonotic potential. The objective of this study was to estimate the heritability for susceptibility to infection of cattle by MAP using Johne's testing records (ELISA test for presence of antibodies to MAP in milk or blood) from US Holstein cattle from 2009 to 2016. Data sets were edited to include records from herds with 100 or more total records and sires with 50 or more daughters. Data sets were further edited to include (1) only herds with at least 1 positive test, (2) herds with at least 2.5% positive test results, and (3) herds with at least 5% positive test results to examine the effect of data from herds with higher proportions of positive tests, and presumably higher pathogen exposure, on heritability estimates. Two models were used in this study, a linear sire model and a binary threshold-probit sire model. Both were mixed models considering fixed effects of herd and age at test, the latter as a covariate accounting for linear and quadratic effects; random effects included sire and residual. Analyses were conducted using a restricted maximum likelihood method. Heritability estimates (±standard error) from the linear model were 0.041 ± 0.004, 0.050 ± 0.004, and 0.062 ± 0.007 for data from herds with at least 1 positive test, ≥2.5% positive tests, and ≥5% positive tests, respectively. Heritability estimates from the threshold model were 0.157 ± 0.014, 0.174 ± 0.016, and 0.186 ± 0.021 for data from herds with at least 1 positive test, ≥2.5% positive tests, and ≥5% positive tests, respectively. Heritability estimates from the linear model were affected by population incidence for positive tests, in contrast to estimates from the threshold model, likely accounting for the difference in magnitude of heritability estimates between models and suggesting that the threshold model analysis is the better choice. Heritability estimates increased as data were restricted to herds with presumed higher MAP exposure for both linear model and threshold model analyses. These estimates are similar to previous estimates in other dairy cattle populations and suggest the potential for selection to lessen susceptibility to MAP infection.

A positional candidate gene association analysis of susceptibility to paratuberculosis on bovine chromosome 7.

Mycobacterium avium subspecies paratuberculosis (MAP) is the causative pathogen for paratuberculosis, which is a chronic inflammation of the small intestine in ruminants and some wild animals. It affects negatively on the economics of dairy operations worldwide and has a zoonotic concern for its potential relationship to Crohn's disease in humans. In this study, we used different approaches to investigate genetic markers associated with MAP infection on bovine chromosome 7. The interleukin-2-inducible T-cell kinase (ITK) gene is a non-receptor tyrosine kinase expressed in T cells, which also is involved in regulation and development of immune cells and some inflammatory diseases. The gene was suggested in a previous GWAS analysis as a positional candidate gene located on BTA7 for MAP infection. In this study, bovine ITK was sequenced and seventeen identified SNPs were genotyped in 1419 Holstein and Jersey cows. Association analysis revealed no significantly associated SNP in the bovine ITK gene using either single marker or haplotype analyses. In a complementary analysis, Holstein genotypes were imputed from 50 K SNPs to the full genome sequence of BTA7. Single SNP association tests identified two SNPs at 15 Mb (p = 5 × 10-7) significantly associated with MAP infection. Our results suggest an additional region of BTA7 contributes to susceptibility to MAP infection in cattle, relative to our previous report, and further investigations are required.

Proteomic analysis of follicular fluid in carriers and non-carriers of the Trio allele for high ovulation rate in cattle.

This study was conducted to characterise differences in follicular fluid proteins between carriers and non-carriers of a bovine allele for high ovulation rate. A total of four non-carrier and five carrier females were used in an initial study with four and six additional non-carriers and carriers respectively used in a validation study. Emergence of the follicular wave was synchronised and the ovaries containing the dominant follicle(s) were extracted by ovariectomy for follicular fluid collection. A hexapeptide ligand library was used to overcome the masking effect of high-abundance proteins and to increase detection of low-abundance proteins in tandem mass spectrometry. After correcting for multiple comparisons, only two proteins, glia-derived nexin precursor (SERPINE2) and inhibin ß B chain precursor (INHBB), were significantly differentially expressed (false-discovery rate <0.05). In a replicate study of analogous design differential expression was confirmed (P<0.05). Joint analysis of results from the two studies indicated that three additional proteins were consistently differentially expressed between genotypes. For three of these five, previous studies have indicated that expression is increased by transforming growth factor-ß-bone morphogenetic protein signalling; their reduction in follicular fluid from carrier animals is consistent with the ~9-fold overexpression of SMAD family member 6 (SMAD6) in carriers that is inhibitory to this pathway.

Trio a novel bovine high-fecundity allele: II. Hormonal profile and follicular dynamics underlying the high ovulation rate.

The newly discovered Trio high-fecundity allele produces multiple ovulations in cattle. This study evaluated (1) size and growth rates of follicles in Trio carriers during a synchronized follicular wave, induced by follicle aspiration; (2) follicle-stimulating hormone (FSH) patterns associated with the follicular wave; (3) size of corpora lutea (CL) and circulating progesterone; and (4) intrafollicular estradiol concentrations prior to normal deviation. Trio carriers had mean dominant follicles that were significantly smaller in diameter and volume than noncarriers. Onset of diameter deviation occurred at ∼3 days after the last follicle aspiration in both genotypes despite Trio carriers having much smaller individual follicles. Follicles of Trio carriers grew at a slower rate than noncarrier follicles (∼65% in mm/day or ∼30% in mm3/day) resulting in much smaller individual dominant follicles (∼25% volume). However, total dominant follicle volume, calculated as the sum of all dominant follicles in each animal, was similar in carriers and noncarriers of Trio throughout the entire follicular wave. Circulating FSH was greater in Trio carriers during the 24 h encompassing deviation. Trio carriers had significantly more ovulations than noncarriers, and individual CL volume was smaller, although total luteal tissue volume and circulating P4 were not different. Thus, increased ovulation rate in Trio carriers relates to smaller individual follicles (one-third the volume) near the time of deviation due to slower follicle growth rate, although time of deviation is similar, with increased circulating FSH near deviation leading to selection of multiple dominant follicles in Trio carriers with similar total follicle volume.

Trio, a novel bovine high fecundity allele: III. Acquisition of dominance and ovulatory capacity at a smaller follicle size.

The acquisition of dominance and ovulatory capacity was evaluated in follicles from cows that were carriers or half-sibling noncarriers of the Trio allele. Follicle size at acquisition of follicular dominance was determined by evaluating whether follicles ovulate after GnRH challenge (ovulatory capacity-experiment 1) and by determination of intrafollicular concentrations of estradiol and free insulin like growth factor 1 (IGF1) and relative mRNA expression of cytochrome P450 family 19 subfamily A member 1 (CYP19A1), luteinizing hormone/choriogonadotropin receptor (LHCGR), and pappalysin 1 (PAPPA, previously known as pregnancy-associated plasma protein A, pappalysin 1) in granulosa cells from follicles of different sizes (experiment 2). Ovulatory capacity developed in follicles at 8.3 mm (50% ovulatory capacity) in noncarriers but at smaller sizes (5.5 mm) in Trio carriers. Similarly, in experiment 2, follicles of Trio carriers acquired a dominant phenotype, as determined by intrafollicular estradiol and CYP19A1, LHCGR, and PAPPA mRNA expression in granulosa cells, at significantly smaller sizes but at a similar time after wave emergence. Overall, dominance/ovulatory capacity was acquired when follicles of Trio carriers were ∼30% the size (volume basis) of follicles in noncarriers. In addition, follicles in Trio carriers appear to acquire dominance in a hierarchal manner, as demonstrated by the progressively greater number of follicles with a dominant phenotype between days 2 and 4 after wave emergence. Thus, results from this study provide further support for a physiological model in which selection of multiple follicles in Trio allele carriers is characterized by acquisition of dominance at a smaller follicle size but at a similar time in the follicular wave with multiple follicles acquiring dominance in a hierarchal sequence.

Short communication: Heritability of twinning rate in Holstein cattle.

Multiple births or twinning in cattle is a naturally occurring reproductive phenomenon. For dairy cattle, twinning is considered a detrimental trait as it can be harmful to cow and calf as well as costly to the producer. The objective of this study was to examine recent US calving records for the Holstein breed to determine a current estimate of heritability for twinning rate along with effects of season and parity. Two models were used in this study: a linear sire model and a binary threshold-logit sire model. Both were mixed models considering fixed effects and random effects. Analyses were conducted using a restricted maximum likelihood method. Heritability estimates were 0.0192 ± 0.0009 and 0.1420 ± 0.0069 for the linear and threshold models, respectively. Repeatabilities from the linear and threshold-logit models were 0.0443 ± 0.0012 and 0.2310 ± 0.0072, respectively. The nonzero estimates of heritability indicate the potential to select against this trait for genetic improvement of Holstein cattle.

Mechanisms regulating follicle selection in ruminants: lessons learned from multiple ovulation models.

Selection of a single dominant follicle from a cohort of growing follicles is a unique biological process, a key step in female reproductive function in monovular species, and lies at the core of reproductive technologies in cattle. Follicle growth and the number of follicles that ovulate are regulated by precise endocrine, paracrine, and autocrine mechanisms. Most of our current understanding about follicle selection focuses on the role of FSH, LH, and the IGF family in follicle growth and selection of the dominant follicle. However, more recently the role of members of the TGF-ß family has been highlighted, particularly in high fecundity genotypes in sheep. Intercellular signaling between the oocyte and granulosa cells (GC) regulates proliferation and differentiation due to actions of bone morphogenetic protein 15 (BMP15) and growth and differentiation factor 9 (GDF9) within the follicle. Mutations that either knockout or reduce the activity of BMP15 or GDF9 have been found to increase ovulation rate in heterozygotes and generally cause severe follicle abnormalities in homozygotes. A mutation in the intracellular kinase domain of the BMPR1B receptor (Booroola fecundity gene) increases ovulation rate in heterozygotes with further increases in ovulation in homozygotes. The physiological mechanisms linking these mutations to increased ovulation rates are still not well defined. A recently identified high fecundity bovine genotype, Trio, causes increased expression of SMAD6, an intracellular inhibitor of the BMP15/GDF9 signalling pathways. This bovine model has provided insights into the mechanisms associated with selection of multiple dominant follicles and multiple ovulations in carriers of fecundity alleles. The present review focuses on the mechanisms involved in follicle selection in ruminants with a special emphasis on the contribution made by multiple ovulation models in both cattle and sheep. The evaluation of multiple ovulation models in ruminants has allowed us to construct a new physiological model that relates changes in the BMP15/GDF9 signalling pathways to the physiological changes that result in selection of multiple dominant follicles. This model is characterized by acquisition of dominance at a smaller follicle size but at a similar time in the follicular wave with multiple follicles acquiring dominance in a hierarchal sequence, delaying FSH suppression and, thus allowing additional follicles to continue to grow and acquire dominance.

Trio, a novel high fecundity allele: I. Transcriptome analysis of granulosa cells from carriers and noncarriers of a major gene for bovine ovulation rate.

A major gene for bovine ovulation rate has been mapped to a 1.2 Mb region of chromosome 10. Screening of coding regions of positional candidate genes within this region failed to reveal a causative polymorphism, leading to the hypothesis that the phenotype results from differences in candidate gene expression rather than alteration of gene structure. This study tested differences in expression of positional candidate genes in granulosa cells between carriers and noncarriers of the high fecundity allele, as well as characterizing differences in the transcriptomic profile between genotypes. Five carriers and five noncarriers, female descendants of "Trio," a carrier of the high fecundity allele were initially used in an RNA-seq analysis of gene expression. Four of ten samples were contaminated with theca cells, so that six samples were used in the final analysis (three of each genotype). Of 14 973 genes expressed, 143 were differentially expressed (false discovery rate P < 0.05) in carriers versus noncarriers. Among the positional candidate genes, SMAD6 was 6.6-fold overexpressed in the carriers compared to noncarriers (P < 5 × 10-5). This result was replicated in an independent group of 12 females (7 carriers and 5 noncarriers) using quantitative real-time PCR; SMAD6 was 9.3-fold overexpressed in carriers versus noncarriers (P = 1.17 × 10-6). Association of overexpression of SMAD6, an inhibitor of the BMP/SMAD signaling pathway, with high ovulation rate corresponds well with disabling mutations in ligands (BMP15 and GDF9) and a receptor (BMPR1B) of this pathway that cause increased ovulation rate in sheep.

T and B cell activation profiles from cows with and without Johne's disease in response to in vitro stimulation with Mycobacterium avium subspecies paratuberculosis.

Johne's disease (JD) is a chronic wasting disease of ruminants caused by infection with Mycobacterium avium subspecies paratuberculosis (MAP). JD is particularly problematic on US dairy farms: estimates show that over 50% of farms are MAP-contaminated and as many as 91% of dairy herds could be infected. Although estimates vary widely, JD may cost the dairy industry between $200 million and $1.5 billion every year. One major obstacle to JD management is that JD is difficult to detect in many animals, in part due to the variable immunity against MAP demonstrated by JD+ cattle. To characterize the diversity of immune responses against MAP, peripheral blood mononuclear cells from 154 JD test negative and 96 JD test positive cows from the same dairy herds were stimulated with MAP in vitro. The activation of CD4+, CD8+ and γδ T cells and surface IgM+ B cells was measured using flow cytometry. CD4+CD45R0+ T cells, γδ+MHCII+ and γδ+MHCII- T cells and SIgM+ B cells from JD test positive cows all exhibited increased proportions expressing CD25 after MAP stimulation, while CD8+ T cells did not demonstrate increased CD25 expression in response to MAP.

Ovulation rate, antral follicle count, and circulating anti-Müllerian hormone in Trio allele carriers, a novel high fecundity bovine genotype.

High fecundity genotypes in sheep are a valuable model to study the physiological mechanisms underlying follicle selection and the control of ovulation rate. Similar genotypes in cattle had not been described until the recent identification of a major bovine allele, termed Trio, which had a large effect on ovulation rate. The present study was designed to evaluate ovulation rate, antral follicle count (AFC), circulating ant-müllerian hormone (AMH), and the association among these measures in unstimulated and superstimulated Trio carrier cattle. We hypothesized that AFC and AMH would be variable among individual cows but would be similar between Trio carriers and non-carrier control cows and that there would be no association between these measures of follicle numbers and ovulation rate. In experiment 1, ovulation rate was determined during 4 consecutive estrous cycles in Trio carriers (n = 34) and non-carrier controls (n = 27). Ovulation rate, on average, was greater (P < 0.01) in Trio carriers (3.5 ± 0.2) compared to non-carrier controls (1.1 ± 0.1) with ∼70% of carrier cycles (n = 136) having 3-4 ovulations while only ∼5% had single ovulations. In contrast, non-carrier cycles (n = 108) were mostly single ovulation (89%) with none having more than two ovulations. In experiment 2, AFC, determined at wave emergence, was not different (P = 0.54) between Trio carriers (24.5 ± 1.3; n = 45) and non-carrier controls (23.1 ± 0.9; n = 37), and no correlation was found between AFC and mean ovulation rate in either genotype (r = -0.009 and r = -0.07; P > 0.70, respectively). In Experiment 3, circulating AMH was also not different between genotypes (P = 0.65) while correlations were found between AFC and AMH in Trio carriers (r = 0.43; P = 0.05; n = 27) and non-carrier controls (r = 0.78; P < 0.01; n = 19). In experiment 4, AFC and AMH were determined in Trio-carriers (n = 9) in relation to a synchronized follicular wave which was unstimulated or stimulated with exogenous FSH. Stimulation with FSH increased ovulation rate, compared to unstimulated Trio carriers, however no association was found between AFC or AMH and ovulation rate regardless of whether superstimulation with exogenous FSH was used. In conclusion, the novel high fecundity bovine genotype Trio, results in consistent multiple ovulations despite having similar AFC and AMH. Therefore, our results suggest that differences in antral follicle numbers during the final stages of follicle development are not a key component of the mechanism underlying multiple ovulations in Trio carriers.

Follicular waves and hormonal profiles during the estrous cycle of carriers and non-carriers of the Trio allele, a major bovine gene for high ovulation and fecundity.

A high fecundity bovine genotype has recently been discovered and genetic mapping indicates evidence for segregation of a major gene with influence on ovulation rate located on bovine chromosome 10. Cattle carrying the high fecundity allele, referred to as the Trio allele, have multiple ovulations while half-sibling, non-carriers generally have single ovulations. The present study was designed to evaluate follicle wave patterns and associated circulating hormones during the estrous cycle of Trio allele carriers (n = 7) and non-carrier half-sib controls (n = 5). We hypothesized that Trio allele carriers would exhibit multiple smaller dominant follicles and greater circulating FSH than non-carrier controls. The proportion of Trio carrier and non-carrier cows with 2 or 3-wave patterns was not different between genotypes with the majority (>70%) exhibiting 3-wave patterns. Trio carriers had greater (P < 0.01) number of ovulations (∼4 vs ∼1 ovulations) and smaller preovulatory follicles (8.9 vs. 14.9 mm; P < 0.01) than non-carrier controls. However, total luteal tissue volume and circulating progesterone, normalized to the initial ovulation or to the onset of luteolysis, were not different between genotypes (P > 0.10). Follicular waves were found to be associated with an FSH surge in both genotypes. Peak FSH concentration at each follicular wave (3-wave patterns) was not different (P > 0.05) between genotypes, but circulating FSH during the decline and nadir, encompassing the day of deviation, was greater (P < 0.05) in Trio carriers. Despite a difference (P = 0.032) in the length of the estrous cycle (23.0 vs. 25.2 d; Trio carrier and non-carriers respectively), the pattern of follicle growth, such as day of wave emergence, day of follicle deviation, and day of maximum diameter of the dominant follicle, when normalized to the FSH surge of each follicular wave were similar in Trio carriers compared to non-carriers although Trio carriers consistently had much smaller-sized follicles (P < 0.05). Thus, decreased follicle size and greater circulating FSH are key components of the mechanism that produces multiple ovulations in cattle that are carriers of the Trio high fecundity allele.

An across-breed genome wide association analysis of susceptibility to paratuberculosis in dairy cattle.

Paratuberculosis is a chronic disease of ruminants caused by Mycobacterium avium subspecies paratuberculosis (MAP). It occurs worldwide and causes a significant loss in the animal production industry. There is no cure for MAP infection and vaccination is problematic. Identification of genetics of susceptibility could be a useful adjunct for programs that focus on management, testing and culling of diseased animals. A case-control, genome-wide association study (GWAS) was conducted using Holstein and Jersey cattle in a combined analysis in order to identify markers and chromosomal regions associated with susceptibility to MAP infection across-breed. A mixed-model method (GRAMMAR-GC) implemented in the GenABEL R package and a Bayes C analysis implemented in GenSel software were used as alternative approaches to conduct GWAS analysis focused on single SNPs and chromosomal segments, respectively. After conducting quality control, 22 406 SNPs from 2157 individuals were available for the GRAMMAR-GC (Bayes C) analysis and 45 640 SNPs from 2199 individuals were available for the Bayes C analysis. One SNP located on BTA27 (8·6 Mb) was identified as moderately associated (P < 5 × 10-5, FDR = 0·44) in the GRAMMAR-GC analysis of the combined breed data. Nine 1 Mb windows located on BTA 2, 3 (3 windows), 6, 8, 25, 27 and 29 each explained ≥1% of the total proportion of genetic variance in the Bayes C analysis. In an analysis ignoring differences in linkage phase, two moderately significantly associated SNPs were identified; ARS-BFGL-NGS-19381 on BTA23 (32 Mb) and Hapmap40994-BTA-46361 on BTA19 (61 Mb). New common genomic regions and candidate genes have been identified from the across-breed analysis that might be involved in the immune response and susceptibility to MAP infection.