Assessment of genetic diversity and population structure of U.S. Polypay sheep from breed origins to future genomic selection.

Knowledge of past and present genetic diversity within a breed is critical for the design and optimization of breeding programs as well as the development of strategies for the conservation of genetic resources. The Polypay sheep breed was developed at the U.S. Sheep Experiment Station (USSES) in 1968 with the goal of improving productivity in Western U.S. range flocks. It has since flourished in the more intensively managed production systems throughout the U.S. The genetic diversity of the breed has yet to be documented. Therefore, the primary objective of this study was to perform a comprehensive evaluation of the genetic diversity and population structure of U.S. Polypay sheep using both pedigree- and genomic-based methods. Pedigree data from 193 Polypay flocks participating in the National Sheep Improvement Program (NSIP) were combined with pedigree records from USSES (n = 162,997), tracing back to the breed's origin. A subset of these pedigreed sheep from 32 flocks born from 2011 to 2023 were genotyped with the GGP Ovine 50K BeadChip containing 51,867 single nucleotide polymorphisms (SNPs). Four subgroups were used for the pedigree-based analyses: 1) the current generation of animals born in 2020-2022 (n = 20,701), 2) the current generation with a minimum of four generations of known ancestors (n = 12,685), 3) only genotyped animals (n = 1,856), and 4) the sires of the current generation (n = 509). Pedigree-based inbreeding for the full population was 2.2%, with a rate of inbreeding of 0.22% per generation. Pedigree-based inbreeding, Wright's inbreeding, and genomic inbreeding based on runs of homozygosity were 2.9%, 1.3%, and 5.1%, respectively, for the genotyped population. The effective population size ranged from 41 to 249 for the pedigree-based methods and 118 for the genomic-based estimate. Expected and observed heterozygosity levels were 0.409 and 0.403, respectively. Population substructure was evident based on the fixation index (FST), principal component analysis, and model-based population structure. These analyses provided evidence of differentiation from the foundation flock (USSES). Overall, the Polypay breed exhibited substantial genetic diversity and the presence of a population substructure that provides a basis for the implementation of genomic selection in the breed.

Puberty classifications in beef heifers are moderately to highly heritable and associated with candidate genes related to cyclicity and timing of puberty.

Introduction: Pubertal attainment is critical to reproductive longevity in heifers. Previously, four heifer pubertal classifications were identified according to attainment of blood plasma progesterone concentrations > 1 ng/ml: 1) Early; 2) Typical; 3) Start-Stop; and 4) Non-Cycling. Early and Typical heifers initiated and maintained cyclicity, Start-Stop started and then stopped cyclicity and Non-Cycling never initiated cyclicity. Start-Stop heifers segregated into Start-Stop-Discontinuous (SSD) or Start-Stop-Start (SSS), with SSD having similar phenotypes to Non-Cycling and SSS to Typical heifers. We hypothesized that these pubertal classifications are heritable, and loci associated with pubertal classifications could be identified by genome wide association studies (GWAS). Methods: Heifers (n = 532; 2017 - 2022) genotyped on the Illumina Bovine SNP50 v2 or GGP Bovine 100K SNP panels were used for variant component estimation and GWAS. Heritability was estimated using a univariate Bayesian animal model. Results: When considering pubertal classifications: Early, Typical, SSS, SSD, and Non-Cycling, pubertal class was moderately heritable (0.38 ± 0.08). However, when heifers who initiated and maintained cyclicity were compared to those that did not cycle (Early+Typical vs. SSD+Non-Cycling) heritability was greater (0.59 ± 0.19). A GWAS did not identify single nucleotide polymorphisms (SNPs) significantly associated with pubertal classifications, indicating puberty is a polygenic trait. A candidate gene approach was used, which fitted SNPs within or nearby a set of 71 candidate genes previously associated with puberty, PCOS, cyclicity, regulation of hormone secretion, signal transduction, and methylation. Eight genes/regions were associated with pubertal classifications, and twenty-two genes/regions were associated with whether puberty was attained during the trial. Additionally, whole genome sequencing (WGS) data on 33 heifers were aligned to the reference genome (ARS-UCD1.2) to identify variants in FSHR, a gene critical to pubertal attainment. Fisher's exact test determined if FSHR SNPs segregated by pubertal classification. Two FSHR SNPs that were not on the bovine SNP panel were selected for additional genotyping and analysis, and one was associated with pubertal classifications and whether they cycled during the trial. Discussion: In summary, these pubertal classifications are moderately to highly heritable and polygenic. Consequently, genomic tools to inform selection/management of replacement heifers would be useful if informed by SNPs associated with cyclicity and early pubertal attainment.

Daily Eicosapentaenoic Acid Infusion in IUGR Fetal Lambs Reduced Systemic Inflammation, Increased Muscle ADRß2 Content, and Improved Myoblast Function and Muscle Growth.

Intrauterine growth-restricted (IUGR) fetuses exhibit systemic inflammation that contributes to programmed deficits in myoblast function and muscle growth. Thus, we sought to determine if targeting fetal inflammation improves muscle growth outcomes. Heat stress-induced IUGR fetal lambs were infused with eicosapentaenoic acid (IUGR+EPA; n = 9) or saline (IUGR; n = 8) for 5 days during late gestation and compared to saline-infused controls (n = 11). Circulating eicosapentaenoic acid was 42% less (p < 0.05) for IUGR fetuses but was recovered in IUGR+EPA fetuses. The infusion did not improve placental function or fetal O2 but resolved the 67% greater (p < 0.05) circulating TNFα observed in IUGR fetuses. This improved myoblast function and muscle growth, as the 23% reduction (p < 0.05) in the ex vivo differentiation of IUGR myoblasts was resolved in IUGR+EPA myoblasts. Semitendinosus, longissimus dorsi, and flexor digitorum superficialis muscles were 24-39% lighter (p < 0.05) for IUGR but not for IUGR+EPA fetuses. Elevated (p < 0.05) IL6R and reduced (p < 0.05) ß2 adrenoceptor content in IUGR muscle indicated enhanced inflammatory sensitivity and diminished ß2 adrenergic sensitivity. Although IL6R remained elevated, ß2 adrenoceptor deficits were resolved in IUGR+EPA muscle, demonstrating a unique underlying mechanism for muscle dysregulation. These findings show that fetal inflammation contributes to IUGR muscle growth deficits and thus may be an effective target for intervention.

An autosomal recessive variant in PYGM causes myophosphorylase deficiency in Red Angus composite cattle.

BACKGROUND: Between 2020 and 2022, eight calves in a Nebraska herd (composite Simmental, Red Angus, Gelbvieh) displayed exercise intolerance during forced activity. In some cases, the calves collapsed and did not recover. Available sire pedigrees contained a paternal ancestor within 2-4 generations in all affected calves. Pedigrees of the calves' dams were unavailable, however, the cows were ranch-raised and retained from prior breeding seasons, where bulls used for breeding occasionally had a common ancestor. Therefore, it was hypothesized that a de novo autosomal recessive variant was causative of exercise intolerance in these calves. RESULTS: A genome-wide association analysis utilizing SNP data from 6 affected calves and 715 herd mates, followed by whole-genome sequencing of 2 affected calves led to the identification of a variant in the gene PYGM (BTA29:g.42989581G > A). The variant, confirmed to be present in the skeletal muscle transcriptome, was predicted to produce a premature stop codon (p.Arg650*). The protein product of PYGM, myophosphorylase, breaks down glycogen in skeletal muscle. Glycogen concentrations were fluorometrically assayed as glucose residues demonstrating significantly elevated glycogen concentrations in affected calves compared to cattle carrying the variant and to wild-type controls. The absence of the PYGM protein product in skeletal muscle was confirmed by immunohistochemistry and label-free quantitative proteomics analysis; muscle degeneration was confirmed in biopsy and necropsy samples. Elevated skeletal muscle glycogen persisted after harvest, resulting in a high pH and dark-cutting beef, which is negatively perceived by consumers and results in an economic loss to the industry. Carriers of the variant did not exhibit differences in meat quality or any measures of animal well-being. CONCLUSIONS: Myophosphorylase deficiency poses welfare concerns for affected animals and negatively impacts the final product. The association of the recessive genotype with dark-cutting beef further demonstrates the importance of genetics to not only animal health but to the quality of their product. Although cattle heterozygous for the variant may not immediately affect the beef industry, identifying carriers will enable selection and breeding strategies to prevent the production of affected calves.

Daily Injection of the ß2 Adrenergic Agonist Clenbuterol Improved Muscle Glucose Metabolism, Glucose-Stimulated Insulin Secretion, and Hyperlipidemia in Juvenile Lambs Following Heat-Stress-Induced Intrauterine Growth Restriction.

Stress-induced fetal programming diminishes ß2 adrenergic tone, which coincides with intrauterine growth restriction (IUGR) and lifelong metabolic dysfunction. We determined if stimulating ß2 adrenergic activity in IUGR-born lambs would improve metabolic outcomes. IUGR lambs that received daily injections of saline or the ß2 agonist clenbuterol from birth to 60 days were compared with controls from pair-fed thermoneutral pregnancies. As juveniles, IUGR lambs exhibited systemic inflammation and robust metabolic dysfunction, including greater (p < 0.05) circulating TNFα, IL-6, and non-esterified fatty acids, increased (p < 0.05) intramuscular glycogen, reduced (p < 0.05) circulating IGF-1, hindlimb blood flow, glucose-stimulated insulin secretion, and muscle glucose oxidation. Daily clenbuterol fully recovered (p < 0.05) circulating TNFα, IL-6, and non-esterified fatty acids, hindlimb blood flow, muscle glucose oxidation, and intramuscular glycogen. Glucose-stimulated insulin secretion was partially recovered (p < 0.05) in clenbuterol-treated IUGR lambs, but circulating IGF-1 was not improved. Circulating triglycerides and HDL cholesterol were elevated (p < 0.05) in clenbuterol-treated IUGR lambs, despite being normal in untreated IUGR lambs. We conclude that deficient ß2 adrenergic regulation is a primary mechanism for several components of metabolic dysfunction in IUGR-born offspring and thus represents a potential therapeutic target for improving metabolic outcomes. Moreover, benefits from the ß2 agonist were likely complemented by its suppression of IUGR-associated inflammation.

A de novo mutation in CACNA1A is associated with autosomal dominant bovine familial convulsions and ataxia in Angus cattle.

Bovine familial convulsions and ataxia (BFCA) is considered an autosomal dominant syndrome with incomplete penetrance. Nine Angus calves from the same herd were diagnosed with BFCA within days of birth. Necropsy revealed cerebellar and spinal cord lesions associated with the condition. Parentage testing confirmed that all affected calves had a common sire. The sire was then bred to 36 cows across two herds using artificial insemination, producing an additional 14 affected calves. The objective of this investigation was to identify hypothesized dominant genetic variation underlying the condition. Whole-genome sequencing was performed on the sire, six affected and seven unaffected paternal half-sibling calves and combined with data from 135 unrelated controls. The sire and five of the six affected calves were heterozygous for a nonsense variant (Chr7 g.12367906C>T, c.5073C>T, p.Arg1681*) in CACNA1A. The other affected calves (N = 8) were heterozygous for the variant but it was absent in the other unaffected calves (N = 7) and parents of the sire. This variant was also absent in sequence data from over 6500 other cattle obtained via public repositories and collaborator projects. The variant in CACNA1A is expressed in the cerebellum of the ataxic calves as detected in the transcriptome and was not differentially expressed compared with controls. The CACNA1A protein is part of a highly expressed cerebellar calcium voltage gated channel. The nonsense variant is proposed to cause haploinsufficiency, preventing proper transmission of neuronal signals through the channel and resulting in BFCA.

A recessive CLN3 variant is responsible for delayed-onset retinal degeneration in Hereford cattle.

Thirteen American Hereford cattle were reported blind with presumed onset when ~12-mo-old. All blind cattle shared a common ancestor through both the maternal and paternal pedigrees, suggesting a recessive genetic origin. Given the pedigree relationships and novel phenotype, we characterized the ophthalmo-pathologic changes associated with blindness and identified the responsible gene variant. Ophthalmologic examinations of 5 blind cattle revealed retinal degeneration. Histologically, 2 blind cattle had loss of the retinal photoreceptor layer. Whole-genome sequencing (WGS) of 7 blind cattle and 9 unaffected relatives revealed a 1-bp frameshift deletion in ceroid lipofuscinosis neuronal 3 (CLN3; chr25 g.26043843del) for which the blind cattle were homozygous and their parents heterozygous. The identified variant in exon 16 of 17 is predicted to truncate the encoded protein (p. Pro369Argfs*8) battenin, which is involved in lysosomal function necessary for photoreceptor layer maintenance. Of 462 cattle genotyped, only blind cattle were homozygous for the deletion. A query of WGS data of > 5,800 animals further revealed that the variant was only observed in related Hereford cattle. Mutations in CLN3 are associated with human juvenile neuronal ceroid lipofuscinosis (JNCL), or Batten disease, which results in early-onset retinal degeneration and lesions similar to those observed in our cases. Our data support the frameshift variant of CLN3 as causative of blindness in these Hereford cattle, and provide additional evidence of the role of this gene in retinal lesions, possibly as a model for human non-syndromic JNCL.

A Comprehensive Allele Specific Expression Resource for the Equine Transcriptome.

Background: Allele-specific expression (ASE) analysis provides a nuanced view of cis-regulatory mechanisms affecting gene expression. Results: In this work, we introduce and highlight the significance of an equine ASE analysis, containing integrated long- and short-read RNA sequencing data, along with insight from histone modification data, from four healthy Thoroughbreds (2 mares and 2 stallions) across 9 tissues. Conclusions: This valuable publicly accessible resource is poised to facilitate investigations into regulatory variation in equine tissues and foster a deeper understanding of the impact of allelic imbalance in equine health and disease at the molecular level.