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.
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.
Paratuberculosis, or Johne's disease, is a chronic, granulomatous, gastrointestinal tract disease of cattle and other ruminants caused by the bacterium Mycobacterium avium, subspecies paratuberculosis (MAP). Control of Johne's disease is based on programs of testing and culling animals positive for infection with MAP while concurrently modifying management to reduce the likelihood of infection. The current study is motivated by the hypothesis that genetic variation in host susceptibility to MAP infection can be dissected and quantifiable associations with genetic markers identified. For this purpose, a case-control, genome-wide association study was conducted using US Holstein cattle phenotyped for MAP infection using a serum ELISA and/or fecal culture test. Cases included cows positive for either serum ELISA, fecal culture or both. Controls consisted of animals negative for the serum ELISA test or both serum ELISA and fecal culture when both were available. Controls were matched by herd and proximal birth date with cases. A total of 856 cows (451 cases and 405 controls) were used in initial discovery analyses, and an additional 263 cows (159 cases and 104 controls) from the same herds were used as a validation data set. Data were analyzed in a single marker analysis controlling for relatedness of individuals (GRAMMAR-GC) and also in a Bayesian analysis in which multiple marker effects were estimated simultaneously (GenSel). For the latter, effects of non-overlapping 1 Mb marker windows across the genome were estimated. Results from the two discovery analyses were generally concordant; however, discovery results were generally not well supported in analysis of the validation data set. A combined analysis of discovery and validation data sets provided strongest support for SNPs and 1 Mb windows on chromosomes 1, 2, 6, 7, 17 and 29.