Enhanced bovine genome annotation through integration of transcriptomics and epi-transcriptomics datasets facilitates genomic biology.

BACKGROUND: The accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues. RESULTS: A total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network. CONCLUSIONS: These validated results show significant improvement over current bovine genome annotations.

Reducing fecal egg count through selective breeding alters dorper lamb response to Haemonchus contortus in an artificial challenge trial.

Infection by gastrointestinal nematodes (GIN), particularly Haemonchus contortus, can be detrimental to sheep health and performance. Genetic susceptibility to GIN varies between breeds, with those lacking high levels of natural resistance often requiring frequent anthelmintic treatment when facing parasitic challenge. Genetic technology can serve as a tool to decrease GIN susceptibility via selection for sheep with reduced fecal egg count (FEC) estimated breeding values (EBVs). However, the physiological changes that result from implementation of this strategy are not well described. Additionally, there is a need for comparison of animals from recent selective breeding against breeds with inherent GIN resistance. In this study we administered a challenge of H. contortus to Dorper x White Dorper (DWD; n = 92) lambs that have been genetically selected for either low (DWD-) or high (DWD+) FEC EBVs and Barbados Blackbelly x Mouflon (BBM; n = 19) lambs from a genetically resistant breed backgrounds. Lamb FEC, packed-cell volume (PCV) and serum IgG were measured at intermittent levels over 5 weeks. At day 21 and day 35, the selectively bred DWD- had a lower mean FEC compared to DWD+, but were higher than BBM. Reductions in both PCV and serum IgG from initial day 0 levels were observed in DWD lambs, but not in BBM. Furthermore, from a subset of lambs (n = 24) harvested at day 21, DWD- only tended (p = 0.056) to have lower mean worm counts than DWD+, with BBM having the lowest mean worm count. Differentially expressed genes (DEGs) identified via RNA-sequencing of abomasal tissue at day 21 indicate a more pronounced Th2 immune response and more rapid worm expulsion occurred in iBBM than iDWD- and iDWD+ lambs. However, gene expression in DWD- suggests an association between reduced FEC EBV and gastric acid secretion and the ability to limit worm fecundity. Ultimately, selection of Dorper sheep for low FEC EBV can reduce susceptibility to GIN, but it will likely require multiple generations with this trait as a breeding priority before presenting a similar resistance level to Caribbean breeds.

Myokines Produced by Cultured Bovine Satellite Cells Harvested from 3- and 11-Month-Old Angus Steers.

The myokines interleukin 6 (IL-6), interleukin 15 (IL-15), myonectin (CTRP15), fibronectin type III domain containing protein 5/irisin (FNDC5), and brain-derived neurotrophic factor (BDNF) are associated with skeletal muscle cell proliferation, differentiation, and muscle hypertrophy in biomedical model species. This study evaluated whether these myokines are produced by cultured bovine satellite cells (BSCs) harvested from 3- and 11-month-old commercial black Angus steers and if the expression and secretion of these targets change across 0, 12, 24, and 48 h in vitro. IL-6, IL-15, FNDC5, and BDNF expression were greater (p ≤ 0.05) in the differentiated vs. undifferentiated BSCs at 0, 12, 24, and 48 h. CTRP15 expression was greater (p ≤ 0.03) in the undifferentiated vs. differentiated BSCs at 24 and 48 h. IL-6 and CTRP15 protein from culture media were greater (p ≤ 0.04) in undifferentiated vs. differentiated BSCs at 0, 12, 24, and 48 h. BDNF protein was greater in the media of differentiated vs. undifferentiated BSCs at 0, 12, 24, and 48 h. IL-6, 1L-15, FNDC5, and BDNF are expressed in association with BSC differentiation, and CTRP15 appears to be expressed in association with BSC proliferation. This study also confirms IL-6, IL-15, CTRP15, and BDNF proteins present in media collected from primary cultures of BSCs.

Current challenges and future of agricultural genomes to phenomes in the USA.

Dramatic improvements in measuring genetic variation across agriculturally relevant populations (genomics) must be matched by improvements in identifying and measuring relevant trait variation in such populations across many environments (phenomics). Identifying the most critical opportunities and challenges in genome to phenome (G2P) research is the focus of this paper. Previously (Genome Biol, 23(1):1-11, 2022), we laid out how Agricultural Genome to Phenome Initiative (AG2PI) will coordinate activities with USA federal government agencies expand public-private partnerships, and engage with external stakeholders to achieve a shared vision of future the AG2PI. Acting on this latter step, AG2PI organized the "Thinking Big: Visualizing the Future of AG2PI" two-day workshop held September 9-10, 2022, in Ames, Iowa, co-hosted with the United State Department of Agriculture's National Institute of Food and Agriculture (USDA NIFA). During the meeting, attendees were asked to use their experience and curiosity to review the current status of agricultural genome to phenome (AG2P) work and envision the future of the AG2P field. The topic summaries composing this paper are distilled from two 1.5-h small group discussions. Challenges and solutions identified across multiple topics at the workshop were explored. We end our discussion with a vision for the future of agricultural progress, identifying two areas of innovation needed: (1) innovate in genetic improvement methods development and evaluation and (2) innovate in agricultural research processes to solve societal problems. To address these needs, we then provide six specific goals that we recommend be implemented immediately in support of advancing AG2P research.

Pedigree diversity and implications for genetic selection of Katahdin sheep.

The Katahdin hair breed gained popularity in the United States as low input and prolific, with a propensity to exhibit parasite resistance. With the introduction of genomically enhanced estimated breeding values (GEBV) to the Katahdin genetic evaluation, defining the diversity present in the breed is pertinent. Utilizing pedigree records (n = 92,030) from 1984 to 2019 from the National Sheep Improvement Program, our objectives were to (i) estimate the completeness and quality of the pedigree, (ii) calculate diversity statistics for the whole pedigree and relevant reference subpopulations and (iii) assess the impact of current diversity on genomic selection. Reference 1 was Katahdins born from 2017 to 2019 (n = 23,494), while reference 2 was a subset with at least three generations of Katahdin ancestry (n = 9327). The completeness of the whole pedigree, and the pedigrees of reference 1 and reference 2, were above 50% through the fourth, fifth and seventh generation of ancestors, respectively. Effective population size (Ne) averaged 111 animals with a range from 42.2 to 451.0. The average generation interval was 2.9 years for the whole pedigree and reference 1, and 2.8 years for reference 2. The mean individual inbreeding and average relatedness coefficients were 1.62% and 0.91%, 1.74% and 0.90% and 2.94% and 1.46% for the whole pedigree, reference 1, and reference 2, respectively. There were over 300 effective founders in the whole pedigree and reference 1, with 169 in reference 2. Effective number of ancestors were over 150 for the whole pedigree and reference 1, while there were 67 for reference 2. Prediction accuracies increased as the reference population grew from 1k to 7.5k and plateaued at 15k animals. Given the large number of founders and ancestors contributing to the base genetic variation in the breed, the Ne is sufficient to maintain diversity while achieving progress with selection. Stable low rates of inbreeding and relatedness suggest that incorporating genetic conservation in breeding decisions is currently not of high priority. Current Ne suggests that with limited genotyping, high levels of accuracy for genomic prediction can be achieved. However, intense selection on GEBV may cause loss of genetic diversity long term.

A computational framework for improving genetic variants identification from 5,061 sheep sequencing data.

BACKGROUND: Pan-genomics is a recently emerging strategy that can be utilized to provide a more comprehensive characterization of genetic variation. Joint calling is routinely used to combine identified variants across multiple related samples. However, the improvement of variants identification using the mutual support information from multiple samples remains quite limited for population-scale genotyping. RESULTS: In this study, we developed a computational framework for joint calling genetic variants from 5,061 sheep by incorporating the sequencing error and optimizing mutual support information from multiple samples' data. The variants were accurately identified from multiple samples by using four steps: (1) Probabilities of variants from two widely used algorithms, GATK and Freebayes, were calculated by Poisson model incorporating base sequencing error potential; (2) The variants with high mapping quality or consistently identified from at least two samples by GATK and Freebayes were used to construct the raw high-confidence identification (rHID) variants database; (3) The high confidence variants identified in single sample were ordered by probability value and controlled by false discovery rate (FDR) using rHID database; (4) To avoid the elimination of potentially true variants from rHID database, the variants that failed FDR were reexamined to rescued potential true variants and ensured high accurate identification variants. The results indicated that the percent of concordant SNPs and Indels from Freebayes and GATK after our new method were significantly improved 12%-32% compared with raw variants and advantageously found low frequency variants of individual sheep involved several traits including nipples number (GPC5), scrapie pathology (PAPSS2), seasonal reproduction and litter size (GRM1), coat color (RAB27A), and lentivirus susceptibility (TMEM154). CONCLUSION: The new method used the computational strategy to reduce the number of false positives, and simultaneously improve the identification of genetic variants. This strategy did not incur any extra cost by using any additional samples or sequencing data information and advantageously identified rare variants which can be important for practical applications of animal breeding.

Angular limb deformity associated with TSPAN18, NRG3 and NOVA2 in Rambouillet rams.

Angular limb deformity (ALD) affects many species of livestock and companion animals. The mechanisms of ALD development are not well understood, but previous research suggests the involvement of genetic risk factors. A case-control genome-wide association study (GWAS) was conducted with 40 ALD-affected and 302 unaffected Rambouillet rams and 40,945 single nucleotide polymorphisms (SNPs). Forelimbs of 6 ALD-affected rams were examined and diagnosed with osteochondrosis. Genome-wide or chromosome-wide significant SNPs were positioned exonic, intronic or within the 3'UTR of genes TSPAN18, NRG3 and NOVA2, respectively. These genes have previously described roles related to angiogenesis and osteoblast, osteoclast and chondrocyte proliferation and differentiation, which suggests the possibility for their involvement in the pathogenesis of osteochondrosis. Functional consequences of SNPs were evaluated through transcription factor binding site analysis, which predicted binding sites for transcription factors of known importance to bone growth, including SOX6, SOX9 and RUNX2. The identification of genetic risk factors for ALD may help to improve animal welfare and production in Rambouillet, a breed known to be at risk for ALD development. This study proposes genes TSPAN18, NRG3 and NOVA2 as targets for further research towards understanding the etiology of ALD in Rambouillet sheep.

Advanced Skeletal Ossification Is Associated with Genetic Variants in Chronologically Young Beef Heifers.

Osteogenesis is a developmental process critical for structural support and the establishment of a dynamic reservoir for calcium and phosphorus. Changes in livestock breeding over the past 100 years have resulted in earlier bone development and increased physical size of cattle. Advanced skeletal maturity is now commonly observed at harvest, with heifers displaying more mature bone than is expected at 30 months of age (MOA). We surmise that selection for growth traits and earlier reproductive maturity resulted in co-selection for accelerated skeletal ossification. This study examines the relationship of single nucleotide polymorphisms (SNPs) in 793 beef heifers under 30 MOA with USDA-graded skeletal maturity phenotypes (A-, B-, C- skeletal maturity). Further, the estrogen content of FDA-approved hormonal implants provided to heifers prior to harvest was evaluated in association with the identified SNPs and maturities. Association tests were performed, and the impact of the implants were evaluated as covariates against genotypes using a logistic regression model. SNPs from the ESR1, ALPL, PPARGC1B, SORCS1 genes, and SNPs near KLF14, ANKRD61, USP42, H1C1, OVCA2, microRNA mir-29a were determined to be associated with the advanced skeletal ossification phenotype in heifers. Higher dosage estrogen implants increased skeletal maturity in heifers with certain SNP genotypes.

Genome-Wide Analysis of Sheep Artificially or Naturally Infected with Gastrointestinal Nematodes.

The anthelmintic resistance of gastrointestinal nematodes (GINs) poses a significant threat to sheep worldwide, but genomic selection can serve as an alternative to the use of chemical treatment as a solution for parasitic infection. The objective of this study is to conduct genome-wide association studies (GWASs) to identify single nucleotide polymorphisms (SNPs) in Rambouillet (RA) and Dorper × White Dorper (DWD) lambs associated with the biological response to a GIN infection. All lambs were genotyped with a medium-density genomic panel with 40,598 markers used for analysis. Separate GWASs were conducted using fecal egg counts (FECs) from lambs (<1 year of age) that acquired their artificial infections via an oral inoculation of 10,000 Haemonchus contortus larvae (n = 145) or naturally while grazing on pasture (n = 184). A GWAS was also performed for packed cell volume (PCV) in artificially GIN-challenged lambs. A total of 26 SNPs exceeded significance and 21 SNPs were in or within 20 kb of genes such as SCUBE1, GALNT6, IGF1R, CAPZB and PTK2B. The ontology analysis of candidate genes signifies the importance of immune cell development, mucin production and cellular signaling for coagulation and wound healing following epithelial damage in the abomasal gastric pits via H. contortus during GIN infection in lambs. These results add to a growing body of the literature that promotes the use of genomic selection for increased sheep resistance to GINs.

AgAnimalGenomes: browsers for viewing and manually annotating farm animal genomes.

Current genome sequencing technologies have made it possible to generate highly contiguous genome assemblies for non-model animal species. Despite advances in genome assembly methods, there is still room for improvement in the delineation of specific gene features in the genomes. Here we present genome visualization and annotation tools to support seven livestock species (bovine, chicken, goat, horse, pig, sheep, and water buffalo), available in a new resource called AgAnimalGenomes. In addition to supporting the manual refinement of gene models, these browsers provide visualization tracks for hundreds of RNAseq experiments, as well as data generated by the Functional Annotation of Animal Genomes (FAANG) Consortium. For species with predicted gene sets from both Ensembl and RefSeq, the browsers provide special tracks showing the thousands of protein-coding genes that disagree across the two gene sources, serving as a valuable resource to alert researchers to gene model issues that may affect data interpretation. We describe the data and search methods available in the new genome browsers and how to use the provided tools to edit and create new gene models.

Response of Rambouillet Lambs to an Artificial Gastrointestinal Nematode Infection.

Gastrointestinal nematodes (GIN) threaten the productivity and health of sheep worldwide, prompting the need for genetic selection to reduce GIN susceptibility. Fecal egg count (FEC), packed-cell volume (PCV), and various production traits were examined in parasitized Rambouillet sheep and compared to sire FEC estimated breeding value (EBV). Rambouillet lambs (n = 77) were inoculated with 10,000 H. contortus L3 larvae. Subsequently, FEC, PCV, and body weight (BW) were captured at seven-day intervals for six weeks. Lambs were sired by one of two rams with post-weaning FEC EBV of −9% or +9%. Mean FEC differed (p = 0.0132) with lambs from the lower EBV sire ("Sire L") being reduced, versus those from the higher EBV sire ("Sire H"), being 2135 ± 211 vs. 2912 ± 207 eggs per gram, respectively. Males and females did not differ for FEC, but females exhibited a higher mean PCV than males, (33.74 vs. 29.65%, p < 0.0001). Lambs were shorn ~120 d post artificial infection and wool measurements were captured. A negative correlation between FEC and grease fleece weight was observed. Our results describe the response of Rambouillet lambs to artificial H. contortus infection and suggest FEC EBV can reduce susceptibility to GIN in this breed.

Single Nucleotide Polymorphism Effects on Lamb Fecal Egg Count Estimated Breeding Values in Progeny-Tested Katahdin Sires.

Estimated breeding values (EBV) for fecal egg counts (FEC) at 42-90 days of age (WFEC) and 91-150 days of age (PFEC) for 84 progeny-tested Katahdin sires were used to identify associations of deregressed EBV with single-nucleotide polymorphisms (SNP) using 388,000 SNP with minor-allele frequencies ≥0.10 on an Illumina high-density ovine array. Associations between markers and FEC EBV were initially quantified by single-SNP linear regression. Effects of linkage disequilibrium (LD) were minimized by assigning SNP to 2,535 consecutive 1-Mb bins and focusing on the effect of the most significant SNP in each bin. Bonferroni correction was used to define bin-based (BB) genome- and chromosome-wide significance. Six bins on chromosome 5 achieved BB genome-wide significance for PFEC EBV, and three of those SNP achieved chromosome-wide significance after Bonferroni correction based on the 14,530 total SNP on chromosome 5. These bins were nested within 12 consecutive bins between 59 and 71 Mb on chromosome 5 that reached BB chromosome-wide significance. The largest SNP effects were at 63, 67, and 70 Mb, with LD among these SNP of r 2 ≤ 0.2. Regional heritability mapping (RHM) was then used to evaluate the ability of different genomic regions to account for additive variance in FEC EBV. Chromosome-level RHM indicated that one 500-SNP window between 65.9 and 69.9 Mb accounted for significant variation in PFEC EBV. Five additional 500-SNP windows between 59.3 and 71.6 Mb reached suggestive (p < 0.10) significance for PFEC EBV. Although previous studies rarely identified markers for parasite resistance on chromosome 5, the IL12B gene at 68.5 Mb codes for the p40 subunit of both interleukins 12 and 23. Other immunoregulatory genes are also located in this region of chromosome 5, providing opportunity for additive or associative effects.

A high-density genome-wide association with absolute blood monocyte count in domestic sheep identifies novel loci.

Monocytes are a core component of the immune system that arise from bone marrow and differentiate into cells responsible for phagocytosis and antigen presentation. Their derivatives are often responsible for the initiation of the adaptive immune response. Monocytes and macrophages are central in both controlling and propagating infectious diseases such as infection by Coxiella burnetii and small ruminant lentivirus in sheep. Genotypes from 513 Rambouillet, Polypay, and Columbia sheep (Ovis aries) were generated using the Ovine SNP50 BeadChip. Of these sheep, 222 animals were subsequently genotyped with the Ovine Infinium® HD SNP BeadChip to increase SNP coverage. Data from the 222 HD genotyped sheep were combined with the data from an additional 258 unique sheep to form a 480-sheep reference panel; this panel was used to impute the low-density genotypes to the HD genotyping density. Then, a genome-wide association analysis was conducted to identify loci associated with absolute monocyte counts from blood. The analysis used a single-locus mixed linear model implementing EMMAX with age and ten principal components as fixed effects. Two genome-wide significant peaks (p < 5x10-7) were identified on chromosomes 9 and 1, and ten genome-wide suggestive peaks (p < 1x10-5) were identified on chromosomes 1, 2, 3, 4, 9, 10, 15, and 16. The identified loci were within or near genes including KCNK9, involved into cytokine production, LY6D, a member of a superfamily of genes, some of which subset monocyte lineages, and HMGN1, which encodes a chromatin regulator associated with myeloid cell differentiation. Further investigation of these loci is being conducted to understand their contributions to monocyte counts. Investigating the genetic basis of monocyte lineages and numbers may in turn provide information about pathogens of veterinary importance and elucidate fundamental immunology.

Variants Within Genes EDIL3 and ADGRB3 are Associated With Divergent Fecal Egg Counts in Katahdin Sheep at Weaning.

Gastrointestinal nematodes (GIN) pose a severe threat to sheep production worldwide. Anthelmintic drug resistance coupled with growing concern regarding potential environmental effects of drug use have demonstrated the necessity of implementing other methods of GIN control. The aim of this study was to test for genetic variants associated with resistance or susceptibility to GIN in Katahdin sheep to improve the current understanding of the genetic mechanisms responsible for host response to GIN. Linear regression and case-control genome-wide association studies were conducted with high-density genotype data and cube-root transformed weaning fecal egg counts (tFEC) of 583 Katahdin sheep. The case-control GWAS identified two significant SNPs (P-values 1.49e-08 to 1.01e-08) within introns of the gene adhesion G protein-coupled receptor B3 (ADGRB3) associated with lower fecal egg counts. With linear regression, four significant SNPs (P-values 7.82e-08 to 3.34e-08) were identified within the first intron of the gene EGF-like repeats and discoidin domains 3 (EDIL3). These identified SNPs were in very high linkage disequilibrium (r 2 of 0.996-1), and animals with alternate homozygous genotypes had significantly higher median weaning tFEC phenotypes compared to all other genotypes. Significant SNPs were queried through public databases to identify putative transcription factor binding site (TFBS) and potential lncRNA differences between reference and alternate alleles. Changes in TFBS were predicted at two SNPs, and one significant SNP was found to be within a predicted lncRNA sequence with greater than 90% similarity to a known lncRNA in the bovine genome. The gene EDIL3 has been described in other species for its roles in the inhibition and resolution of inflammation. Potential changes of EDIL3 expression mediated through lncRNA expression and/or transcription factor binding may impact the overall immune response and reduce the ability of Katahdin sheep to control GIN infection. This study lays the foundation for further research of EDIL3 and ADGRB3 towards understanding genetic mechanisms of susceptibility to GIN, and suggests these SNPs may contribute to genetic strategies for improving parasite resistance traits in sheep.

An improved ovine reference genome assembly to facilitate in-depth functional annotation of the sheep genome.

BACKGROUND: The domestic sheep (Ovis aries) is an important agricultural species raised for meat, wool, and milk across the world. A high-quality reference genome for this species enhances the ability to discover genetic mechanisms influencing biological traits. Furthermore, a high-quality reference genome allows for precise functional annotation of gene regulatory elements. The rapid advances in genome assembly algorithms and emergence of sequencing technologies with increasingly long reads provide the opportunity for an improved de novo assembly of the sheep reference genome. FINDINGS: Short-read Illumina (55× coverage), long-read Pacific Biosciences (75× coverage), and Hi-C data from this ewe retrieved from public databases were combined with an additional 50× coverage of Oxford Nanopore data and assembled with canu v1.9. The assembled contigs were scaffolded using Hi-C data with Salsa v2.2, gaps filled with PBsuitev15.8.24, and polished with Nanopolish v0.12.5. After duplicate contig removal with PurgeDups v1.0.1, chromosomes were oriented and polished with 2 rounds of a pipeline that consisted of freebayes v1.3.1 to call variants, Merfin to validate them, and BCFtools to generate the consensus fasta. The ARS-UI_Ramb_v2.0 assembly is 2.63 Gb in length and has improved continuity (contig NG50 of 43.18 Mb), with a 19- and 38-fold decrease in the number of scaffolds compared with Oar_rambouillet_v1.0 and Oar_v4.0. ARS-UI_Ramb_v2.0 has greater per-base accuracy and fewer insertions and deletions identified from mapped RNA sequence than previous assemblies. CONCLUSIONS: The ARS-UI_Ramb_v2.0 assembly is a substantial improvement in contiguity that will optimize the functional annotation of the sheep genome and facilitate improved mapping accuracy of genetic variant and expression data for traits in sheep.

Genetic association of wool quality characteristics in United States Rambouillet sheep.

Introduction: Fine wool production is an important source of revenue, accounting for up to 13% of total revenue in extensively managed wool sheep production systems of the United States. The Rambouillet are a predominant breed that excels in wool quality characteristics. Understanding the genetic basis of wool quality characteristics would aid in the development of genomic breeding strategies to facilitate genetic improvement. Methods: Wool characteristics and DNA were collected for rams enrolled in the North Dakota State University and University of Wyoming annual central performance ram tests over a three-year period (2019-2021, N = 313). The relationships of wool quality characteristics including grease fleece weight adjusted 365 days (wt. 365 adj.), clean fleece wt. 365 adj., staple length 365 adj., average fiber diameter, face wool cover, amount of skin wrinkles and belly wool were evaluated through genome-wide association studies (GWAS), Pearson correlation and ANOVA. Results: The GWAS identified four genome-wide significant genetic markers (p-value <1.19e-06) and five chromosome-wide significant markers (p-value <1.13e-05) on chromosomes 1, 2, 4, 15, and 19. Significant markers were associated with genes notable for relevant wool biological functions, including the gene ABCC8 which codes for SUR1, an ATP-sensitive potassium channel known to affect hair growth and 60S ribosomal protein L17-like, previously found to be expressed during follicle formation. The strongest Pearson correlation coefficients were identified between clean fleece wt. 365 adj. and grease fleece wt. 365 adj. (r = 0.83) and between clean fleece wt. 365 adj. and staple length 365 adj. (r = 0.53). Additionally, clean fleece wt. 365 adj. was correlated with final body weight (r = 0.35) and scrotal circumference (r = 0.16). Staple length 365 adj. (p-value = 5e-04), average fiber diameter (p-value = .0053) and clean fleece wt. 365 adj. (p-value = .014) were significantly associated with belly wool score. Discussion: The results of this study provide important insight into the relationships between wool quality characteristics and report specific markers that Rambouillet sheep producers may use to help inform selection and breeding decisions for improved wool quality.

Using whole genome sequence to compare variant callers and breed differences of US sheep.

As whole genome sequence (WGS) data sets have become abundant and widely available, so has the need for variant detection and scoring. The aim of this study was to compare the accuracy of commonly used variant calling programs, Freebayes and GATK HaplotypeCaller (GATK-HC), and to use U.S. sheep WGS data sets to identify novel breed-associated SNPs. Sequence data from 145 sheep consisting of 14 U.S. breeds were filtered and biallelic single nucleotide polymorphisms (SNPs) were retained for genotyping analyses. Genotypes from both programs were compared to each other and to genotypes from bead arrays. The SNPs from WGS were compared to the bead array data with breed heterozygosity, principal component analysis and identifying breed associated SNPs to analyze genetic diversity. The average sequence read depth was 2.78 reads greater with 6.11% more SNPs being identified in Freebayes compared to GATK-HC. The genotype concordance of the variant callers to bead array data was 96.0% and 95.5% for Freebayes and GATK-HC, respectively. Genotyping with WGS identified 10.5 million SNPs from all 145 sheep. This resulted in an 8% increase in measured heterozygosity and greater breed separation in the principal component analysis compared to the bead array analysis. There were 1,849 SNPs identified in only the Romanov sheep where all 10 rams were homozygous for one allele and the remaining 135 sheep from 13 breeds were homozygous for the opposite allele. Both variant calling programs had greater than 95% concordance of SNPs with bead array data, and either was suitably accurate for ovine WGS data sets. The use of WGS SNPs improved the resolution of PCA analysis and was critical for identifying Romanov breed-associated SNPs. Subsets of such SNPs could be used to estimate germplasm composition in animals without pedigree information.

Variation in type two taste receptor genes is associated with bitter tasting phenylthiocarbamide consumption in mature Targhee and Rambouillet rams.

Bitter taste perception in sheep can lead to avoidance of specific types of forage, such as sagebrush, which is present on many rangeland grazing systems in the Intermountain West. In humans, bitter taste perception is influenced by variation in several TAS2R genes, including more extensively studied TAS2R38 and TAS2R16. We hypothesize that variation in taste receptor genes in sheep is associated with bitter taste. Therefore, the objective of this study was to examine variation in TAS2R genes in relation to consumption of a bitter tasting compound phenylthiocarbamide (PTC) which determines bitter "taster" and "non-taster" status in humans. Rambouillet and Targhee rams (n = 26) were offered various concentrations of PTC solution (0.2-12.29 mM) and water in a side-by-side presentation during two experiments. Blood was collected for DNA isolation and sequencing. Nineteen TAS2R genes were amplified and sequenced with long read Oxford Nanopore MinION technology. A total of 1,049 single nucleotide polymorphisms (SNPs) and 26 haplotypes were identified in these genes. Of these, 24 SNPs and 11 haplotypes were significantly (P < 0.05) associated with PTC consumption in TAS2R3, TAS2R5, TAS2R8, TAS2R9, TAS2R16, TAS2R31-like, TAS2R38, TAS2R39, and TAS2R42-like. Over 50% of the SNPs resulted in a change in amino acid sequence and several resided in potential regulatory regions, which could have downstream functional consequences and influence bitter taste perception in sheep. Further research is needed to validate these associations and elucidate the mechanisms that link variation in TAS2R genes to bitter taste perception in sheep. This may enable producers to select sheep more likely to consume bitter forage such as sagebrush as a flock and rangeland management strategy.

Genes involved in immune, gene translation and chromatin organization pathways associated with Mycoplasma ovipneumoniae presence in nasal secretions of domestic sheep.

Mycoplasma ovipneumoniae contributes to polymicrobial pneumonia in domestic sheep. Elucidation of host genetic influences of M. ovipneumoniae nasal detection has the potential to reduce the incidence of polymicrobial pneumonia in sheep through implementation of selective breeding strategies. Nasal mucosal secretions were collected from 647 sheep from a large US sheep flock. Ewes of three breeds (Polypay n = 222, Rambouillet n = 321, and Suffolk n = 104) ranging in age from one to seven years, were sampled at three different times in the production cycle (February, April, and September/October) over four years (2015 to 2018). The presence and DNA copy number of M. ovipneumoniae was determined using a newly developed species-specific qPCR. Breed (P<0.001), age (P<0.024), sampling time (P<0.001), and year (P<0.001) of collection affected log10 transformed M. ovipneumoniae DNA copy number, where Rambouillet had the lowest (P<0.0001) compared with both Polypay and Suffolk demonstrating a possible genetic component to detection. Samples from yearlings, April, and 2018 had the highest (P<0.046) detected DNA copy number mean. Sheep genomic DNA was genotyped with the Illumina OvineHD BeadChip. Principal component analysis identified most of the variation in the dataset was associated with breed. Therefore, genome wide association analysis was conducted with a mixed model (EMMAX), with principal components 1 to 6 as fixed and a kinship matrix as random effects. Genome-wide significant (P<9x10-8) SNPs were identified on chromosomes 6 and 7 in the all-breed analysis. Individual breed analysis had genome-wide significant (P<9x10-8) SNPs on chromosomes 3, 4, 7, 9, 10, 15, 17, and 22. Annotated genes near these SNPs are part of immune (ANAPC7, CUL5, TMEM229B, PTPN13), gene translation (PIWIL4), and chromatin organization (KDM2B) pathways. Immune genes are expected to have increased expression when leukocytes encounter M. ovipneumoniae which would lead to chromatin reorganization. Work is underway to narrow the range of these associated regions to identify the underlying causal mutations.