Development and validation of a fast and automated DNA identification line.

The importance of DNA evidence for gaining investigative leads demands a fast workflow for forensic DNA profiling performed in large volumes. Therefore, we developed software solutions for automated DNA profile analysis, contamination check, major donor inference, DNA database (DDB) comparison and reporting of the conclusions. This represents the Fast DNA IDentification Line (FIDL) and this study describes its development, validation and implementation in criminal casework at the authors' institute. This first implementation regards single donor profiles and major contributors to mixtures. The validation included testing of the software components on their own and examination of the performance of different DDB search strategies. Furthermore, end-to-end testing was performed under three conditions: (1) testing of scenarios that can occur in DNA casework practice, (2) tests using three months of previous casework data, and (3) testing in a casework production environment in parallel to standard casework practices. The same DNA database candidates were retrieved by this automated line as by the manual workflow. The data flow was correct, results were reproducible and robust, results requiring manual analysis were correctly flagged, and reported results were as expected. Overall, we found FIDL valid for use in casework practice in our institute. The results from FIDL are automatically reported within three working days from receiving the trace sample. This includes the time needed for registration of the case, DNA extraction, quantification, polymerase chain reaction and capillary electrophoresis. FIDL itself takes less than two hours from intake of the raw CE data to reporting. Reported conclusions are one of five options: (1) candidate retrieved from DDB, (2) no candidate retrieved from DDB, (3) high evidential value with regards to reference within the case, (4) results require examination of expert, or (5) insufficient amount of DNA obtained to generate a DNA profile. In our current process, the automated report is sent within three working days and a complete report, with confirmation of the FIDL results, and signed by a reporting officer is sent at a later time. The signed report may include additional analyses regarding e.g. minor contributors. The automated report with first case results is quickly available to the police enabling them to act upon the DNA results prior to receiving the full DNA report. This line enables a uniform and efficient manner of handling large numbers of traces and cases and provides high value investigative leads in the early stages of the investigation.

Multi-laboratory validation of DNAxs including the statistical library DNAStatistX.

This study describes a multi-laboratory validation of DNAxs, a DNA eXpert System for the data management and probabilistic interpretation of DNA profiles [1], and its statistical library DNAStatistX to which, besides the organising laboratory, four laboratories participated. The software was modified to read multiple data formats and the study was performed prior to the release of the software to the forensic community. The first exercise explored all main functionalities of DNAxs with feedback on user-friendliness, installation and general performance. Next, every laboratory performed likelihood ratio (LR) calculations using their own dataset and a dataset provided by the organising laboratory. The organising laboratory performed LR calculations using all datasets. The datasets were generated with different STR typing kits or analysis systems and consisted of samples varying in DNA amounts, mixture ratios, number of contributors and drop-out level. Hypothesis sets had the correct, under- and over-assigned number of contributors and true and false donors as person of interest. When comparing the results between laboratories, the LRs were foremost within one unit on log10 scale. The few LR results that deviated more had differences for the parameters estimated by the optimizer within DNAStatistX. Some of these were indicated by failed iteration results, others by a failed model validation, since unrealistic hypotheses were included. When these results that do not meet the quality criteria were excluded, as is in accordance with interpretation guidelines, none of the analyses in the different laboratories yielded a different statement in the casework report. Nonetheless, changes in software parameters were sought that minimized differences in outcomes, which made the DNAStatistX module more robust. Overall, the software was found intuitive, user-friendly and valid for use in multiple laboratories.

Using online tools at the Bovine Genome Database to manually annotate genes in the new reference genome.

With the availability of a new highly contiguous Bos taurus reference genome assembly (ARS-UCD1.2), it is the opportune time to upgrade the bovine gene set by seeking input from researchers. Furthermore, advances in graphical genome annotation tools now make it possible for researchers to leverage sequence data generated with the latest technologies to collaboratively curate genes. For many years the Bovine Genome Database (BGD) has provided tools such as the Apollo genome annotation editor to support manual bovine gene curation. The goal of this paper is to explain the reasoning behind the decisions made in the manual gene curation process while providing examples using the existing BGD tools. We will describe the sources of gene annotation evidence provided at the BGD, including RNA-seq and Iso-Seq data. We will also explain how to interpret various data visualizations when curating gene models, and will demonstrate the value of manual gene annotation. The process described here can be applied to manual gene curation for other species with similar tools. With a better understanding of manual gene annotation, researchers will be encouraged to edit gene models and contribute to the enhancement of livestock gene sets.

DNAxs/DNAStatistX: Development and validation of a software suite for the data management and probabilistic interpretation of DNA profiles.

The data management, interpretation and comparison of sets of DNA profiles can be complex, time-consuming and error-prone when performed manually. This, combined with the growing numbers of genetic markers in forensic identification systems calls for expert systems that can automatically compare genotyping results within (large) sets of DNA profiles and assist in profile interpretation. To that aim, we developed a user-friendly software program or DNA eXpert System that is denoted DNAxs. This software includes features to view, infer and match autosomal short tandem repeat profiles with connectivity to up and downstream software programs. Furthermore, DNAxs has imbedded the 'DNAStatistX' module, a statistical library that contains a probabilistic algorithm to calculate likelihood ratios (LRs). This algorithm is largely based on the source code of the quantitative probabilistic genotyping system EuroForMix [1]. The statistical library, DNAStatistX, supports parallel computing which can be delegated to a computer cluster and enables automated queuing of requested LR calculations. DNAStatistX is written in Java and is accessible separately or via DNAxs. Using true and non-contributors to DNA profiles with up to four contributors, the DNAStatistX accuracy and precision were assessed by comparing the DNAStatistX results to those of EuroForMix. Results were the same up to rare differences that could be attributed to the different optimizers used in both software programs. Implementation of dye specific detection thresholds resulted in larger likelihood values and thus a better explanation of the data used in this study. Furthermore, processing time, robustness of DNAStatistX results and the circumstances under which model validations failed were examined. Finally, guidelines for application of the software are shared as an example. The DNAxs software is future-proof as it applies a modular approach by which novel functionalities can be incorporated.

Potential emigration of Siberian cattle germplasm on Chirikof Island, Alaska.

Feral cattle residing in Chirikof Island, Alaska, are relatively distinct from breeds used in commercial production in North America. However, preliminary evidence suggested that they exhibit substantial genetic relationship with cattle from Yakutian region of Siberia. Thus, our objective was to further elucidate quantify the origins, admixture and divergence of the Chirikof Island cattle relative to cattle from Siberia and USA. Subject animals were genotyped at 15 microsatellite loci. Compared with Turano-Mongolian and North American cattle, Chirikof Island cattle had similar variation, with slightly less observed heterozygosity, fewer alleles per locus and a positive fixation index. Analysis of the genetic distances revealed two primary clusters; one that contained the North American breeds and the Kazakh White head, and a second that contained the Yakutian and Kalmyk breeds, and the Chirikof population. Thus, it is suggested that Chirikof Island cattle may be a composite of British breeds emanating from North America and Turano-Mongolian cattle. A potential founder effect, consistent with historical records of the Russian-American period, may contribute to the adaptation of the Chirikof Island cattle to their harsh high-latitude environment. Further study of adaptive mechanisms manifest by these cattle is warranted.

Origins of cattle on Chirikof Island, Alaska, elucidated from genome-wide SNP genotypes.

Feral livestock may harbor genetic variation of commercial, scientific, historical or esthetic value. The origins and uniqueness of feral cattle on Chirikof Island, Alaska, are uncertain. The island is now part of the Alaska Maritime Wildlife Refuge and Federal wildlife managers want grazing to cease, presumably leading to demise of the cattle. Here we characterize the cattle of Chirikof Island relative to extant breeds and discern their origins. Our analyses support the inference that Yakut cattle from Russia arrived first on Chirikof Island, then ~120 years ago the first European taurine cattle were introduced to the island, and finally a large wave of Hereford cattle were introduced on average 40 years ago. In addition, this mixture of European and East-Asian cattle is unique compared with other North American breeds and we find evidence that natural selection in the relatively harsh environment of Chirikof Island has further impacted their genetic architecture. These results provide an objective basis for decisions regarding conservation of the Chirikof Island cattle.

Preovulatory estradiol and the establishment and maintenance of pregnancy in suckled beef cows.

In postpartum beef cows, GnRH-induced ovulation of small dominant follicles decreased pregnancy rates and increased late embryonic/fetal mortality. In Exp. 1, single ovulation reciprocal embryo transfer (ET) was used to examine the relationship between preovulatory serum concentrations of estradiol at GnRH-induced ovulation in donor and recipient cows and establishment and maintenance of pregnancy. Suckled beef cows (n = 1,164) were administered GnRH (GnRH1, 100 µg) on d -9 (GnRH1), PGF(2α) on d -2, and GnRH2 (GnRH2, 100 µg) on d 0 (CO-Synch protocol) either with (donors; n = 810) or without (recipients; n = 354) AI. Single embryos (n = 394) or oocytes (n = 45) were recovered from the donor cows (d 7; ET) and all live embryos were transferred into recipients. Serum concentration of estradiol at GnRH2 was positively correlated with follicle size at GnRH2 (r = 0.45, P < 0.01) and progesterone at ET (r =0.34, P < 0.01). Donor cows with greater estradiol at GnRH2 were more likely to yield an embryo than an unfertilized oocyte (P < 0.01). Donor and recipient cows were retrospectively divided into 4 groups [low estradiol (<8.4 pg/mL) or high estradiol (≥8.4 pg/mL)] based on serum concentration of estradiol at GnRH2. Pregnancy rate at d 27 for low-low (n = 78), low-high (n = 80), high-low (n = 91), and high-high (n = 101) groups (donor-recipient, respectively) was 45, 65, 43, and 61% respectively (P < 0.02). Because recipient cows with greater estradiol concentration at GnRH2 had greater pregnancy rates in Exp. 1, the objective of Exp. 2 was to evaluate the effect of estradiol supplementation on pregnancy rate. Ovulation was synchronized in suckled beef cows (n = 600) using the CO-Synch protocol with the insertion of a controlled internal drug release (CIDR; intravaginal progesterone supplement) from d -9 until d -2. Approximately one-half of the cows (n = 297) received an injection of estradiol cypionate (ECP; 0.5 mg intramuscularly) 24 h before AI. Compared with the no treatment (Control) cows, ECP treatment increased (P < 0.01) pregnancy rates of cows induced to ovulate smaller dominant follicles (<12.2 mm). In conclusion, GnRH-induced ovulation of small dominant follicles was associated with reduced serum estradiol, fertilization rate (donor cows), and pregnancy establishment (recipient cows). Furthermore, ECP supplementation during the preovulatory period increased pregnancy rates in cows induced to ovulate smaller dominant follicles.

Pregnancy establishment and maintenance in cattle.

A single ovulation, reciprocal embryo transfer study was used to investigate effects of oocyte competence and maternal environment on pregnancy establishment and maintenance in beef cows. Estrous cycles were synchronized in suckled beef cows and embryo donors were inseminated on d 0 (n = 810). Cows were classified on d 0 as having a small (<12.5 mm) or large (≥12.5 mm) ovulatory follicle and randomly chosen as donors or recipients to remove confounding effects of ovulatory follicle size on fertility. Embryos (n = 393) or oocytes (n = 44) were recovered on d 7, and all viable embryos were transferred into recipients (n = 354). All statistical analyses were conducted using the GLM procedure of SAS. Path analysis (with significance set at P < 0.10) was used to examine potential cause-effect relationships among the measured variables. Greater donor cow BW, circulating estradiol concentration at insemination, postpartum interval, and ovulatory follicle size directly increased (P < 0.10) fertilization success. Greater donor cow age was the only factor that directly decreased (P < 0.10) fertilization success. Viability of d-7 embryos was directly inhibited (P < 0.10) by rapid follicular growth rate from d -2 to 0 and heavier BW. Direct beneficial effects to embryo viability were increased serum progesterone concentration on d -2 and ovulatory follicle size. Pregnancy maintenance from d 7 to 27 was enhanced (P < 0.10) by increased serum estradiol concentration on d 0 and progesterone concentration on d 7 in the recipient cow. Increased follicular diameter in the recipient cow on d 0 was detrimental to pregnancy maintenance from d 7 to 27. This manuscript defines the complex interplay and relative contributions of endocrine and physical factors both prior and subsequent to fertilization that influence both oocyte competence and maternal environment and their roles in establishment and maintenance of pregnancy.

Effects of reduced panel, reference origin, and genetic relationship on imputation of genotypes in Hereford cattle.

The objective of this study was to investigate alternative methods of designing and using reduced SNP panels for imputing SNP genotypes. Two purebred Hereford populations, an experimental population known as Line 1 Hereford (L1, n = 240) and registered Hereford with American Hereford Association (AHA, n = 311), were used. Using different reference samples of 62 to 311 animals with 39,497 SNP on 29 autosomes and study samples of 57 or 62 animals for which genotypes were available for ~2,600 SNP (reduced panels), imputations were performed to predict the other ~36,900 loci that had been masked. An imputation package, including LinkPHASE and DAGPHASE, was used for imputation. Four reduced panels differing in minor allele frequency (MAF) and marker spacing were evaluated. Reduced panels included every 15th SNP across the genome (SNP_space), commercial Illumina Bovine3K Beadchip (SNP_3K), SNP with the highest MAF (SNP_MAF), and SNP with high MAF that were also evenly spaced across the genome (SNP_MS). Imputation accuracy was defined as the correlation of imputed genotypes and real genotypes. Reference samples were either from L1 or AHA. Among animals with genotypes, genetic relationships were estimated based on molecular marker genotypes or pedigree. Reduced panel design, number of animals in the reference sample, reference origin and genetic relationship between animals in the reference, and study samples all affected imputation accuracy (P < 0.001). Across genotyping schemes, imputed genotypes from SNP_MS had the greatest accuracy. A 0.1 increase in average pedigree relationship or average molecular relationship between reference and study samples increased imputation accuracy 10 to 20%. Using reference samples from the L1 population resulted in lower imputation accuracy than using reference samples from the admixed population AHA (P < 0.001). Increasing the number of animals in the reference panel by 100 individuals increased imputation accuracy by 8% when pedigree relationship was used as a covariate and 6% when molecular relationship was used as a covariate. We concluded that imputation accuracy would be increased through optimization of reduced panel design and genotyping strategy.

Using 50 k single nucleotide polymorphisms to elucidate genomic architecture of line 1 hereford cattle.

Hereford is a major beef breed in the USA, and a sub-population, known as Line 1 (L1), was established in 1934 using two paternal half-sib bulls and 50 unrelated females. L1 has since been maintained as a closed population and selected for growth to 1 year of age. Objectives were to characterize the molecular genetic architecture of L1 (n = 240) by comparing a cross-section of L1 with the general US. Hereford population (AHA, n = 311), estimating effects of imposed selection within L1 based on allele frequencies at 50 K SNP loci, and examining loci-specific effects of heterozygosity on the selection criterion. Animals were genotyped using the Illumina BovineSNP50 Beadchip, and SNP were mapped to UMD3.0 assembly of the bovine genome sequence. Average linkage disequilibrium (LD), measured by square of Pearson correlation, of adjacent SNP was 0.36 and 0.16 in L1 and AHA, respectively. Difference in LD between L1 and AHA decreased as SNP spacing increased. Persistence of phase between L1 and AHA decreased from 0.45 to 0.14 as SNP spacing increased from 50 to 5,000 kb. Extended haplotype homozygosity was greater in L1 than in AHA for 95.6% of the SNP. Knowledge of selection applied to L1 facilitated a novel approach to QTL discovery. Minor allele frequency was (FDR < 0.01) affected by cumulative selection differential at 191 out of 25,901 SNP. With the FDR relaxed to 0.05, 13 regions on BTA2, 5, 6, 9, 11, 14, 15, 18, 23, and 26 are co-located with previously identified QTL for growth. After adjustment of postweaning gain phenotypes for fixed effects and direct additive genetic effects, regression of residuals on genome-wide heterozygosity was -235.3 ± 91.6 kg. However, no SNP-specific loci where heterozygotes were significantly superior to the average of homozygotes were revealed (FDR ≥ 0.17). In conclusion, genome-wide SNP genotypes clarified effects of selection and inbreeding within L1 and differences in genomic architecture between the population segment L1 and the AHA population.

Quantitative trait loci with effects on feed efficiency traits in Hereford x composite double backcross populations.

Two half-sib families of backcross progeny were produced by mating F(1) Line 1 Hereford (L1) x composite gene combination (CGC) bulls with L1 and CGC cows. Feed intake and periodic weights were measured for 218 backcross progeny. These progenies were genotyped using 232 microsatellite markers that spanned the 29 BTA. Progeny from L1 and CGC females was analysed separately using composite interval mapping to find quantitative trait loci (QTL) affecting daily dry matter intake (DMI), average daily gain (ADG), feed conversion (FCR) and residual feed intake (RFI). Results from both backcrosses were pooled to find additional QTL. In the backcross to L1, QTL were detected for RFI and DMI on BTA11, FCR on BTA16, and ADG on BTA9. In the backcross to CGC, QTL were detected for RFI on BTA10, FCR on BTA12 and 16 and ADG on BTA15 and 17. After pooling, QTL were detected for RFI on BTA 2, 6, 7, 10, 11, 13 and 16; for FCR on BTA 9, 12, 16, 17 and 21; for ADG on BTA 9, 14, 15, 17; and for DMI on BTA 2, 5, 6, 9, 10, 11, 20 and 23.

A Limousin specific myostatin allele affects longissimus muscle area and fatty acid profiles in a Wagyu-Limousin F2 population.

A microsatellite-based genome scan of a Wagyu x Limousin F(2) cross population previously demonstrated QTL affecting LM area and fatty acid composition were present in regions near the centromere of BTA2. In this study, we used 70 SNP markers to examine the centromeric 24 megabases (Mb) of BTA2, including the Limousin-specific F94L myostatin allele (AB076403.1; 415C > A) located at approximately 6 Mb on the draft genome sequence of BTA2. A significant effect of the F94L marker was observed (F = 60.17) for LM area, which indicated that myostatin is most likely responsible for the effect. This is consistent with previous reports that the substitution of Leu for Phe at AA 94 of myostatin (caused by the 415C > A transversion) is associated with increased muscle growth. Surprisingly, several fatty acid trait QTL, which affected the amount of unsaturated fats, also mapped to or very near the myostatin marker, including the ratio of C16:1 MUFA to C16:0 saturated fat (F = 16.72), C18:1 to C18:0 (F = 18.88), and total content of MUFA (F = 17.12). In addition, QTL for extent of marbling (F = 14.73) approached significance (P = 0.05), and CLA concentration (F = 9.22) was marginally significant (P = 0.18). We also observed associations of SNP located at 16.3 Mb with KPH (F = 15.00) and for the amount of SFA (F = 12.01). These results provide insight into genetic differences between the Wagyu and Limousin breeds and may lead to a better tasting and healthier product for consumers through improved selection for lipid content of beef.

Detection of germline and somatic copy number variations in cattle.

As a complement to the Bovine HapMap Consortium project, we initiated a systematic study of the copy numbervariation (CNV) within the same cattle population using array comparative genomic hybridization (array CGH). Oligonucleotide CGH arrays were designed and fabricated to cover all chromosomes with an average interval of 6 kb using the latest bovine genome assembly. In the initial screening, three Holstein bulls were selected to represent major paternal lineages of the Holstein breed with some maternal linkages between these lines. Dual-label hybridizations were performed using either Hereford L1 Dominette 01449 or L1 Domino 99375 as reference. The CNVs were represented by gains and losses of normalized fluorescence intensities relative to the reference. The data presented here, for the first time, demonstrated that significant amounts of germline and fewer somatic CNVs exist in cattle, that many CNVs are common both across diverse cattle breeds and among individuals within a breed, and that array CGH is an effective tool to systematically detect bovine CNV. Selected CNVs have been confirmed by independent methods using real-time (RT) PCR. The strategy used in this study, based on genome higher-orderarchitecture variation, is a powerful approach to generating resources for the identification of novel genomic variation and candidate genes for economically important traits.

Comparative and physical mapping of 111 previously reported and 105 new porcine microsatellites.

Here we report radiation hybrid mapping of 105 new porcine microsatellite markers on the IMpRH(7000) radiation hybrid panel. In addition, we searched flanking sequences of these markers, as well as 673 previously reported RH-mapped microsatellite markers, for orthology to human sequences. Eighty-seven new and 111 previously mapped sequences exhibited orthology to human sequences. Using a stringent sequence alignment, 25 microsatellite-flanking sequences were found to be highly similar to genic sequences, whereas 173 were similar to non-genic sequences in the human genome. Five markers were located near known breakpoints of synteny between human and swine.

Quantitative trait loci with additive effects on palatability and fatty acid composition of meat in a Wagyu-Limousin F2 population.

A whole-genome scan was conducted on 328 F(2) progeny in a Wagyu x Limousin cross to identify quantitative trait loci (QTL) affecting palatability and fatty acid composition of beef at an age-constant endpoint. We have identified seven QTL on five chromosomes involved in lipid metabolism and tenderness. None of the genes encoding major enzymes involved in fatty acid metabolism, such as fatty acid synthase (FASN), acetyl-CoA carboxylase alpha (ACACA), solute carrier family 2 (facilitated glucose transporter) member 4 (SLC2A4), stearoyl-CoA desaturase (SCD) and genes encoding the subunits of fatty acid elongase, was located in these QTL regions. The present study may lead to a better-tasting and healthier product for consumers through improved selection for palatability and lipid content of beef.

Influence of seasonal forage quality on glucose kinetics of young beef cows.

Extensive range livestock production systems in the western United States rely heavily on rangeland forages to meet the nutritional needs of grazing livestock throughout the year. Interannual variation in the quantity and quality of rangeland forage in the Northern Great Plains, as well as throughout much of the western United States, may play a pivotal role in how well grazing ruminants sequester nutrients in their tissues. This variation in forage quality may influence the ability of a beef cow to utilize dietary nutrients via changes in tissue responsiveness to insulin. Identifying specific periods and production states in which this phenomenon is manifested will provide insight into the development and implementation of strategic and targeted supplementation practices that improve nutrient utilization during times of nutritional imbalance and may improve the lifetime productivity of grazing range beef cows. A 2-yr study was conducted to monitor serum metabolites, glucose kinetics during glucose tolerance tests, and forage chemical composition every 90 d in young cows (2 to 4 yr of age; n = 28). In yr 1 and 2, cows were managed on 4 pastures varying in size from 36 to 76 ha in yr 1 and 49 to 78 ha in yr 2. Regardless of year, cow age, or cow physiological status, the main factor influencing glucose half-life was season of the year (P = 0.02). Effects of season on glucose half-life closely followed assessments describing forage quality, with glucose half-lives of 46, 39, 43, and 51 +/- 3.9 min for May, August, December, and March, respectively. Elevated glucose half-life during seasons in which forage quality is of lower nutritive value indicated that tissue responsiveness to the actions of insulin followed seasonal changes in forage quality. Glucose half-life tended (P = 0.09) to decrease between May and August, increased (P = 0.04) between December and March, and showed a tendency (P = 0.10) to decrease in seasons of greater nutrient density (May and August) compared with seasons of lower nutrient density (December and March). Seasonal changes in serum metabolites were also observed and corresponded with changes in forage quality. The results support our hypothesis that as the season progresses and forage quality declines, maternal tissues become less responsive to insulin, indicating that targeted supplementation with glucogenic precursors during these seasons of nutritional stress may improve cow performance.

Quantitative trait loci with additive effects on growth and carcass traits in a Wagyu-Limousin F2 population.

A whole-genome scan for carcass traits [average daily gain during the pre-weaning, growth and finishing periods; birth weight; hot carcass weight and longissimus muscle area (LMA)] was performed on 328 F(2) progeny produced from Wagyu x Limousin-cross parents derived from eight founder Wagyu bulls. Nine significant (P

Genetic relationships between feral cattle from Chirikof Island, Alaska and other breeds.

The origin of cattle on Chirikof Island, off the coast of Alaska, is not well documented. We assessed genetic differentiation of cattle isolated on Chirikof Island from several breeds commonly used for commercial production in North America including breeds popularly believed to have contributed to the Chirikof Island population. A set of 34 microsatellite loci was used to genotype Angus, Charolais, Hereford, Highland, Limousin, Red Angus, Salers, Shorthorn, Simmental, Tarentaise and Texas Longhorn cattle sampled from North America and the Chirikof Island population. Resulting F(ST) statistics for these loci ranged from 0.06 to 0.22 and on average, 14% of total genetic variation was between breeds. Whether population structure was modelled as a bifurcating tree or genetic network, Chirikof Island cattle appeared to be unique and strongly differentiated relative to the other breeds that were sampled. Bayesian clustering for multiple-locus assignment to genetic groups indicated low levels of admixture in the Chirikof Island population. Thus, the Chirikof Island population may be a novel genetic resource of some importance for conservation and industry.

Genetic partitioning of variation in ovulatory follicle size and probability of pregnancy in beef cattle.

The objectives of this research were to partition variation in ovulatory follicle size into genetic and nongenetic components and to assess the utility of ovulatory follicle size as an indicator trait associated with reproductive success in beef cattle. Data were collected during the years 2002 to 2005 from 780 beef females that ranged in age from 1 to 12 yr (mean of 2.4 observations per female). Data were analyzed with a multiple trait Gibbs sampler for animal models to make Bayesian inferences from flat priors. A chain of 500,000 Gibbs samples was thinned to every 200th sample to produce a posterior distribution composed of 2,500 samples. Heritability estimates (posterior mean +/- SD) were 0.16 +/- 0.03 for follicle size and 0.07 +/- 0.02 and 0.02 +/- 0.01 for pregnancy rate as a trait of the female and service sire, respectively. Posterior means of genetic correlations were all <0.10, with 0.00 contained within the respective 90% probability density posterior intervals. Results indicate that whereas follicle size is of greater heritability than pregnancy rate, its usefulness to improve reproductive rate is greatest as an ancillary phenotype in multiple trait selection.