ABSTRACT
Telomere-to-telomere (T2T) assemblies reveal new insights into the structure and function of the previously 'invisible' parts of the genome and allow comparative analyses of complete genomes across entire clades. We present here an open collaborative effort, termed the 'Ruminant T2T Consortium' (RT2T), that aims to generate complete diploid assemblies for numerous species of the Artiodactyla suborder Ruminantia to examine chromosomal evolution in the context of natural selection and domestication of species used as livestock.
Subject(s)
Ruminants , Telomere , Telomere/genetics , Animals , Ruminants/genetics , Evolution, Molecular , Genome/genetics , Selection, Genetic , Phylogeny , DiploidyABSTRACT
BACKGROUND: Managing genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (Ne), runs of homozygosity (ROH) and Wright's fixation index (FST) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs). RESULTS: Katahdin sheep had the largest current Ne at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the FST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the FST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all FST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms. CONCLUSIONS: In this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity.
Subject(s)
Polymorphism, Single Nucleotide , Selection, Genetic , Animals , Sheep/genetics , Homozygote , Genetic Variation , United States , Inbreeding , Sheep, Domestic/genetics , Breeding/methodsABSTRACT
Skeletal muscle growth is an economically important trait in the cattle industry. Secreted muscle-derived proteins, referred to as myokines, have important roles in regulating the growth, metabolism, and health of skeletal muscle in human and biomedical research models. Accumulating evidence supports the importance of myokines in skeletal muscle and whole-body health, though little is known about the potential presence and functional significance of these proteins in cattle. This study evaluates and confirms that secreted proteins acidic and rich in cysteine (SPARC), fibroblast growth factor 21 (FGF-21), myostatin (MSTN), and decorin (DCN) are expressed and SPARC, FGF-21, and DCN are secreted by primary bovine satellite cells from 3- (BSC3; n = 3) and 11- (BSC11; n = 3) month -old commercial angus steers. Cells were cultured and collected at zero, 12, 24, and 48 hours to characterize temporal expression and secretion from undifferentiated and differentiated cells. The expression of SPARC was higher in the undifferentiated (p = 0.04) and differentiated (p = 0.07) BSC11 than BSC3. The same was observed with protein secretion from undifferentiated (p <0.0001) BSC11 compared to BSC3. Protein secretion of FGF-21 was higher in undifferentiated BSC11 (p < 0.0001) vs. BSC3. DCN expression was higher in differentiated BSC11 (p = 0.006) vs. BSC3. Comparing undifferentiated vs. differentiated BSC, MSTN expression was higher in differentiated BSC3 (p ≤ 0.001) for 0, 12, and 24 hours and in BSC11 (p ≤ 0.03) for 0, 12, 24, and 48 hours. There is also a change over time for SPARC expression (p ≤ 0.03) in undifferentiated and differentiated BSC and protein secretion (p < 0.0001) in undifferentiated BSC, as well as FGF-21 expression (p = 0.007) in differentiated BSC. This study confirms SPARC, FGF-21, and DCN are secreted, and SPARC, FGF-21, MSTN, and DCN are expressed in primary bovine muscle cells with age and temporal differences.
Subject(s)
Cell Differentiation , Decorin , Fibroblast Growth Factors , Osteonectin , Animals , Cattle , Osteonectin/metabolism , Osteonectin/genetics , Fibroblast Growth Factors/metabolism , Decorin/metabolism , Cells, Cultured , Male , Satellite Cells, Skeletal Muscle/metabolism , Satellite Cells, Skeletal Muscle/cytology , Aging/metabolism , Myostatin/metabolism , Muscle, Skeletal/metabolism , Muscle, Skeletal/cytologyABSTRACT
Introduction: Improving ewe longevity is an important breeding and management goal, as death loss and early culling of mature ewes are economic burdens in the sheep industry. Ewe longevity can be improved by selecting for positive reproductive outcomes. However, the breeding approaches for accomplishing this come with the challenge of recording a lifetime trait. Characterizing genetic factors underpinning ewe longevity and related traits could result in the development of genomic selection strategies to improve the stayability of sheep through early, informed selection of replacement ewes. Methods: Towards this aim, a genome-wide association study (GWAS) was performed to identify genetic markers associated with ewe longevity, reproductive, and production traits. Traits evaluated included longevity (i.e., length of time in the flock), parity and the lifetime number of lambs born, lambs born alive, lambs weaned, and weight of lambs weaned. Ewe records from previous studies were used. Specifically, Rambouillet (n = 480), Polypay (n = 404), Suffolk (n = 182), and Columbia (n = 64) breed ewes (N = 1,130) were analyzed against 503,617 SNPs in across-breed and within-breed GWAS conducted with the Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK) model in R. Results: The across-breed GWAS identified 25 significant SNPs and the within-breed GWAS for Rambouillet, Polypay, and Suffolk ewes identified an additional 19 significant SNPs. The most significant markers were rs411309094 (13:22,467,143) associated with longevity in across-breed GWAS (p-value = 8.3E-13) and rs429525276 (2:148,398,336) associated with both longevity (p-value = 6.4E-15) and parity (p-value = 4.8E-15) in Rambouillet GWAS. Significant SNPs were identified within or in proximity (±50 kb) of genes with known or proposed roles in reproduction, dentition, and the immune system. These genes include ALPL, ANOS1, ARHGEF26, ASIC2, ASTN2, ATP8A2, CAMK2D, CEP89, DISC1, ITGB6, KCNH8, MBNL3, MINDY4, MTSS1, PLEKHA7, PRIM2, RNF43, ROBO2, SLCO1A2, TMEM266, TNFRSF21, and ZNF804B. Discussion: This study proposes multiple SNPs as candidates for use in selection indices and suggests genes for further research towards improving understanding of the genetic factors contributing to longevity, reproductive, and production traits of ewes.
ABSTRACT
Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. We assembled complete and gapless telomere to telomere (T2T) Y chromosomes for these species. The pseudo-autosomal regions were similar in length, but the total chromosome size was substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity was accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 18MYA. The centromeres also differed dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosome have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies.
ABSTRACT
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.
Subject(s)
Gene Expression Profiling , Genomics , Cattle/genetics , Animals , Sequence Analysis, RNA , Transcriptome , Quantitative Trait Loci , RNA , Protein Isoforms , Molecular Sequence AnnotationABSTRACT
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.
Subject(s)
Feces , Haemonchiasis , Haemonchus , Parasite Egg Count , Sheep Diseases , Animals , Sheep , Sheep Diseases/parasitology , Sheep Diseases/immunology , Sheep Diseases/genetics , Haemonchiasis/veterinary , Haemonchiasis/parasitology , Haemonchiasis/immunology , Parasite Egg Count/veterinary , Feces/parasitology , Selective Breeding , Male , Female , Genetic Predisposition to Disease , BreedingABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
Subject(s)
Genome-Wide Association Study , Nerve Tissue Proteins , Animals , Male , Sheep/genetics , Introns/genetics , RNA-Binding Proteins/genetics , ExonsABSTRACT
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.
Subject(s)
Bone Development , Osteogenesis , Cattle/genetics , Animals , Female , Osteogenesis/genetics , Genotype , Bone and Bones , EstrogensABSTRACT
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.
Subject(s)
Haemonchiasis , Nematoda , Sheep Diseases , Animals , Sheep/genetics , Genome-Wide Association Study , Sheep Diseases/genetics , Sheep Diseases/parasitology , Haemonchiasis/genetics , Haemonchiasis/veterinary , Nematoda/genetics , Gastrointestinal TractABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
Subject(s)
Mycoplasma ovipneumoniae/physiology , Sheep, Domestic/genetics , Sheep, Domestic/microbiology , Animals , Genome-Wide Association Study , Genotype , Lung/microbiology , Sheep , Sheep, Domestic/immunologyABSTRACT
The ability of livestock to reproduce efficiently is critical to the sustainability of animal agriculture. Antral follicle count (AFC) and reproductive tract scores (RTS) can be used to estimate fertility in beef heifers, but the genetic mechanisms influencing variation in these measures are not well understood. Two genome-wide association studies (GWAS) were conducted to identify the significant loci associated with these traits. In total, 293 crossbred beef heifers were genotyped on the Bovine GGP 50K chip and genotypes were imputed to 836,121 markers. A GWAS was performed with the AFC phenotype for 217 heifers with a multi-locus mixed model, conducted using the year, age at time of sampling and principal component analysis groupings as the covariates. The RTS GWAS was performed with 289 heifers using an additive correlation/trend test comparing prepubertal to pubertal heifers. The loci on chromosomes 2, 3 and 23 were significant in the AFC GWAS and the loci on chromosomes 2, 8, 10 and 11 were significant in the RTS GWAS. The significant region on chromosome 2 was similar between both analyses. These regions contained genes associated with cell proliferation, transcription, apoptosis and development. This study proposes candidate genes for beef cattle fertility, although future research is needed to elucidate the precise mechanisms.