Evolution of the sheep industry and genetic research in the United States: opportunities for convergence in the twenty-first century.

The continuous development and application of technology for genetic improvement is a key element for advancing sheep production in the United States. The US sheep industry has contracted over time but appears to be at a juncture where a greater utilization of technology can facilitate industry expansion to new markets and address inefficiencies in traditional production practices. Significant transformations include the increased value of lamb in relation to wool, and a downtrend in large-scale operations but a simultaneous rise in small flocks. Additionally, popularity of hair breeds not requiring shearing has surged, particularly in semi-arid and subtropical US environments. A variety of domestically developed composite breeds and newly established technological approaches are now widely available for the sheep industry to use as it navigates these ongoing transformations. These genetic resources can also address long-targeted areas of improvement such as growth, reproduction and parasite resistance. Moderate progress in production efficiency has been achieved by producers who have employed estimated breeding values, but widespread adoption of this technology has been limited. Genomic marker panels have recently shown promise for reducing disease susceptibility, identifying parentage and providing a foundation for marker-assisted selection. As the ovine genome is further explored and genomic assemblies are improved, the sheep research community in the USA can capitalize on new-found information to develop and apply genetic technologies to improve the production efficiency and profitability of the sheep industry.

Genetic structure and admixture in sheep from terminal breeds in the United States.

Selection for performance in diverse production settings has resulted in variation across sheep breeds worldwide. Although sheep are an important species to the United States, the current genetic relationship among many terminal sire breeds is not well characterized. Suffolk, Hampshire, Shropshire and Oxford (terminal) and Rambouillet (dual purpose) sheep (n = 248) sampled from different flocks were genotyped using the Applied Biosystems Axiom Ovine Genotyping Array (50K), and additional Shropshire sheep (n = 26) using the Illumina Ovine SNP50 BeadChip. Relationships were investigated by calculating observed heterozygosity, inbreeding coefficients, eigenvalues, pairwise Wright's FST estimates and an identity by state matrix. The mean observed heterozygosity for each breed ranged from 0.30 to 0.35 and was consistent with data reported in other US and Australian sheep. Suffolk from two different regions of the United States (Midwest and West) clustered separately in eigenvalue plots and the rectangular cladogram. Further, divergence was detected between Suffolk from different regions with Wright's FST estimate. Shropshire animals showed the greatest divergence from other terminal breeds in this study. Admixture between breeds was examined using admixture, and based on cross-validation estimates, the best fit number of populations (clusters) was K = 6. The greatest admixture was observed within Hampshire, Suffolk, and Shropshire breeds. When plotting eigenvalues, US terminal breeds clustered separately in comparison with sheep from other locations of the world. Understanding the genetic relationships between terminal sire breeds in sheep will inform us about the potential applicability of markers derived in one breed to other breeds based on relatedness.

Genome-wide association study to identify genetic loci associated with gastrointestinal nematode resistance in Katahdin sheep.

Resistance to gastrointestinal nematodes has previously been shown to be a moderately heritable trait in some breeds of sheep, but the mechanisms of resistance are not well understood. Selection for resistance currently relies upon faecal egg counts (FEC), blood packed cell volumes and FAMACHA visual indicator scores of anaemia. Identifying genomic markers associated with disease resistance would potentially improve the selection process and provide a more reliable means of classifying and understanding the biology behind resistant and susceptible sheep. A GWAS was conducted to identify possible genetic loci associated with resistance to Haemonchus contortus in Katahdin sheep. Forty animals were selected from the top and bottom 10% of estimated breeding values for FEC from a total pool of 641 sires and ram lambs. Samples were genotyped using Applied Biosystems™ Axiom™ Ovine Genotyping Array (50K) consisting of 51 572 SNPs. Following quality control, 46 268 SNPs were included in subsequent analyses. Analyses were conducted using a linear regression model in plink v1.90 and a single-locus mixed model in snp and variation suite. Genome-wide significance was determined by a Bonferroni correction for multiple testing. Using linear regression, loci on chromosomes 2, 3, 16, 23 and 24 were significantly associated at the genome level with FEC estimated breeding values, and we identified a region on chromosome 2 that was significant using both statistical analyses. We suggest a potential role for the gene DIS3L2 for gastrointestinal nematode resistance in Katahdin sheep, although further research is needed to validate these findings.

Identification of polymorphisms influencing feed intake and efficiency in beef cattle.

Feed efficiency is an economically important trait in beef cattle. Net feed efficiency, measured as residual feed intake (RFI), is the difference between actual feed intake and the predicted feed intake required for maintenance and gain of the animal. SNPs that show associations with RFI may be useful quantitative trait nucleotides for marker-assisted selection. This study identified associations between SNPs underlying five RFI QTL on five bovine chromosomes (BTA2, 5, 10, 20 and 29) with measures of dry matter intake (DMI), RFI and feed conversion ratio (FCR) in beef cattle. Six SNPs were found to have effects on RFI (P < 0.05). The largest single SNP allele substitution effect for RFI was -0.25 kg/day located on BTA2. The combined effects of the SNPs found significant in this experiment explained 6.9% of the phenotypic variation of RFI. Not all the RFI SNPs showed associations with DMI and FCR even though these traits are highly correlated with RFI (r = 0.77 and r = 0.62 respectively). This shows that these SNPs may be affecting the underlying biological mechanisms of feed efficiency beyond feed intake control and weight gain efficiency. These SNPs can be used in marker-assisted selection but first it will be important to verify these effects in independent populations of cattle.

Polymorphisms and haplotypes in the bovine neuropeptide Y, growth hormone receptor, ghrelin, insulin-like growth factor 2, and uncoupling proteins 2 and 3 genes and their associations with measures of growth, performance, feed efficiency, and carcass merit in beef cattle.

Genes that regulate metabolism and energy partitioning have the potential to influence economically important traits in farm animals, as do polymorphisms within these genes. In the current study, SNP in the bovine neuropeptide Y (NPY), growth hormone receptor (GHR), ghrelin (GHRL), uncoupling proteins 2 and 3 (UCP2 and UCP3), IGF2, corticotrophin-releasing hormone (CRH), cocaine and amphetamine regulated transcript (CART), melanocortin-4 receptor (MC4R), proopiomelanocortin (POMC), and GH genes were evaluated for associations with growth, feed efficiency, and carcass merit in beef steers. In total, 24 SNP were evaluated for associations with these traits and haplotypes were constructed within each gene when 2 or more SNP showed significant associations. An A/G SNP located in intron 4 of the GHR gene had the largest effects on BW of the animals (dominance effect P < 0.01) and feed efficiency (allele substitution effect P < 0.05). Another A/G SNP located in the promoter region of GHR had similar effects but the haplotypes of these 2 SNP reduced the effects of the SNP located in intron 4. Three SNP in the NPY gene showed associations to marbling (P < 0.001) as well as with ADG, BW, and feed conversion ratio (FCR; P < 0.05). The combination of these 3 SNP into haplotypes generally improved the association or had a similar scale of association as each single SNP. Only 1 SNP in UCP3, an A/G SNP in intron 3, was associated with ADG (P = 0.025), partial efficiency of growth, and FCR (P < 0.01). Three SNP in UCP2 gene were in almost complete linkage disequilibrium and showed associations with lean meat yield, yield grade, DMI, and BW (P < 0.05). Haplo-types between the SNP in UCP3 and UCP2 generally reduced the associations seen individually in each SNP. An A/G SNP in the GHRL gene tended to show effects on residual feed intake, FCR, and partial efficiency of growth (P < 0.10). The IGF2 SNP most strongly affected LM area (P < 0.01), back fat, ADG, and FCR (P < 0.05). The SNP in the CART, MC4R, POMC, GH, and CRH genes did not show associations at P < 0.05 with any of the traits. Although most of the SNP that showed associations do not cause amino acid changes, these SNP could be linked to other yet to be detected causative mutations or nearby QTL. It will be very important to verify these results in other cattle populations.

Construction of bovine whole-genome radiation hybrid and linkage maps using high-throughput genotyping.

High-density whole-genome maps are essential for ordering genes or markers and aid in the assembly of genome sequence. To increase the density of markers on the bovine radiation hybrid map, and hence contribute to the assembly of the bovine genome sequence, an Illumina BeadStation was used to simultaneously type large numbers of markers on the Roslin-Cambridge 3000 rad bovine-hamster whole-genome radiation hybrid panel (WGRH3000). In five multiplex reactions, 6738 sequence tagged site (STS) markers were successfully typed on the WGRH3000 panel DNA. These STSs harboured SNPs that were developed as a result of the bovine genome sequencing initiative. Typically, the most time consuming and expensive part of creating high-density radiation hybrid (RH) maps is genotyping the markers on the RH panel with conventional approaches. Using the method described in this article, we have developed a high-density whole-genome RH map with 4690 loci and a linkage map with 2701 loci, with direct comparison to the bovine whole-genome sequence assembly (Btau_2.0) in a fraction of the time it would have taken with conventional typing and genotyping methods.