Genome-wide associations with longevity and reproductive traits in U.S. rangeland ewes.

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.

Model definition for genetic evaluation of purebred and crossbred lambs including heterosis.

Crossbreeding is a common practice among commercial sheep producers to improve animal performance. However, genetic evaluation of U.S. sheep is performed within breed type (terminal sire, semi-prolific, and western range). While incorporating crossbred records may improve assessment of purebreds, it requires accounting for heterotic and breed effects in the evaluation. The objectives of this study were to: 1) describe the development of a paternal composite (PC) line, 2) determine the effect of direct and maternal heterosis on growth traits of crossbred lambs, 3) estimate (co)variance components for direct and maternal additive, and uncorrelated maternal environmental, effects, and 4) provide an interpretation of the estimates of random effects of genetic groups, and to use those solutions to compare the genetic merit of founding breed subpopulations. Data included purebred and crossbred records on birth weight (BN; n = 14,536), pre-weaning weight measured at 39 or 84 d (WN; n = 9,362) depending on year, weaning weight measured at 123 d (WW; n = 9,297), and post-weaning weight measured at 252 d (PW; n = 1,614). Mean (SD) body weights were 5.3 (1.1), 16.8 (3.9) and 28.0 (7.6), 39.1 (7.2), and 54.2 (8.7) kg for BN, WN (at the two ages), WW, and PW, respectively. In designed experiments, the Siremax, Suffolk, Texel, Polypay, Columbia, Rambouillet, and Targhee breeds were compared within the same environment. Estimates of heterotic effects and covariance components were obtained using a multiple trait animal model. Genetic effects based on founders' breeds were significant and included in the model. Percent estimates of direct heterosis were 2.89 ± 0.61, 2.60 ± 0.65, 4.24 ± 0.56, and 6.09 ± 0.86, and estimates of maternal heterosis were 1.92 ± 0.87, 4.64 ± 0.80, 3.95 ± 0.66, and 4.04 ± 0.91, for BN, WN, WW, and PW, respectively. Correspondingly, direct heritability estimates were 0.17 ± 0.02, 0.13 ± 0.02, 0.17 ± 0.02, and 0.46 ± 0.04 for BN, WN, WW, and PW. Additive maternal effects accounted for trivial variation in PW. For BN, WN, and WW, respectively, maternal heritability estimates were 0.16 ± 0.02, 0.10 ± 0.02, and 0.07 ± 0.01. Uncorrelated maternal environmental effects accounted for little variation in any trait. Direct and maternal heterosis had considerable impact on growth traits, emphasizing the value of crossbreeding and the need to account for heterosis, in addition to breed effects, if crossbred lamb information is included in genetic evaluation.

Crossbreeding is common in commercial sheep enterprises. It allows breeds with different attributes to be combined to generate crossbred progeny tailored to production environments and customer preferences. Additionally, crossbreds often benefit from heterosis, performing at levels above the average of their parental breeds. Over two decades, body weights were collected at birth and at pre-weaning, weaning, and post-weaning ages on purebred and crossbred lambs from semi-prolific (Polypay), western range (Columbia, Rambouillet, Targhee), and terminal sire (Siremax, Suffolk, Texel) breeds at the U.S. Sheep Experiment Station. When combined, the value of direct heterosis­that due to a lamb being crossbred­and maternal heterosis­that due to the lamb's dam being crossbred­increased birth (5%) and post-natal (up to 10%) weights in crossbred lambs. This highlights the value of crossbreeding to the U.S. sheep industry, especially in western range production systems. Genetic variation between and within breeds also was detected for the purebred parental breeds. Such heterotic and breed effects must be accounted for if crossbred performance is to be incorporated in genetic evaluation of purebreds. Therefore, these results provide the foundation for utilizing crossbred information in the evaluation and selection of purebred sheep in the United States.

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.

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.

Association of TMEM8B and SPAG8 with Mature Weight in Sheep.

Signature of selection studies have identified many genomic regions with known functional importance and some without verified functional roles. Multiple studies have identified Transmembrane protein 8B (TMEM8B)rs426272889 as having been recently under extreme selection pressure in domesticated sheep, but no study has provided sheep phenotypic data clarifying a reason for extreme selection. We tested rs426272889 for production trait association in 770 U.S. Rambouillet, Targhee, Polypay, and Suffolk sheep. TMEM8Brs426272889 was associated with mature weight at 3 and 4 years (p < 0.05). This suggested selection for sheep growth and body size might explain the historical extreme selection pressure in this genomic region. We also tested Sperm-associated antigen 8 (SPAG8) rs160159557 encoding a G493C substitution. While this variant was associated with mature weights at ages 3 and 4, it was not as strongly associated as TMEM8Brs426272889. Transmembrane protein 8B has little functional information except as an inhibitor of cancer cell proliferation. To our knowledge, this is the first study linking TMEM8B to whole organism growth and body size under standard conditions. Additional work will be necessary to identify the underlying functional variant(s). Once identified, such variants could be used to improve sheep production through selective breeding.

Evaluation of Rambouillet, Polypay, and Romanov-White Dorper × Rambouillet ewes mated to terminal sires in an extensive rangeland production system: body weight and wool characteristics1,2.

Though lamb sales account for the majority of annual receipts on U.S. sheep operations, wool is an important income source for many Western flocks. Crossing fine-wool sheep with prolific or composite hair sheep breeds can increase lamb production, but fleece quality and marketability may be reduced by a greater content of nonwool fibers (e.g., med and kemp). The objectives of this study were to compare BW and wool characteristics of Rambouillet, Polypay, and Romanov-White Dorper × Rambouillet (RW-RA) ewes under extensive rangeland management conditions. Ewe BW was collected before mating (fall) and 30 d postlambing (spring) each year from 1 yr up to 4 yr of age. In spring and fall, Rambouillet and Polypay ewes were similar in BW (P ≥ 0.94). Spring BW did not differ (P = 0.13) between RW-RA and Polypay, but Rambouillet ewes were heavier than RW-RA in the spring (P = 0.02). Both Rambouillet and Polypay ewes were heavier (P < 0.07) than RW-RA in the fall. Greasy fleece weight and mid-side wool samples were collected from ewes at 1 and 4 yr of age. Clean fleece weights (CFW) were estimated from average laboratory scoured yield of mid-side wool samples composited within-breed. Average fiber diameter (AFD), SD of fiber diameter (SD-FD), and percentage med (%M), kemp (%K), and total medullated fibers (%T) were quantified on individual mid-side wool samples. There was no difference in 1-yr-old CFW among breed types (P ≥ 0.96). Four-yr-old Rambouillet ewes had heavier CFW (2.29 kg; P < 0.001) than 4-yr-old Polypay (1.83 kg) and RW-RA ewes (1.86 kg), which were not different (P > 0.99). Within 1- and 4-yr-olds, AFD differed among breed type (P < 0.001) and was the finest for Rambouillet (20.1 and 21.9 µm, respectively), intermediate for RW-RA (22.8 and 24.8 µm), and coarsest for Polypay (24.2 and 26.7 µm). Also within 1- and 4-yr olds, SD-FD was lowest in Rambouillet, intermediate in Polypay, and greatest in RW-RA (P < 0.01). Wool from RW-RA ewes had greater %M, %K, and %T (P < 0.001) than wool from Rambouillet and Polypay ewes, which were not different (P > 0.99). Results indicated superior wool production for Rambouillet compared with the coarser, more variable wool produced by Polypay and RW-RA. Still, past research reported greater lamb production in Polypay and RW-RA ewes which, under recent market conditions, would be associated with greater annual gross revenue for these breed types than for Rambouillet.

Effects of rearing triplet lambs on ewe productivity, lamb survival and performance, and future ewe performance.

Increasing prolificacy has been proposed to be the most effective way to increase the biological efficiency and profitability of sheep production. However, use of prolific breeds and genes with major effects on ovulation rate can increase prolificacy to levels that may not be desirable or sustainable in extensive rangeland production systems. This study thus evaluated effects of triplet births on ewe productivity and ewe and lamb performance. An initial study used 666 purebred Polypay litters to compare ewes with triplet litters that were required to raise all the lambs (Treatment A) with those whose triplet litters were reduced to 2 lambs (Treatment R). Adult Polypay ewes had an average litter size of 2.35 lambs per litter. The frequency of litters of 3 or more lambs was 43.2%; 56.0% of lambs were born in litters of 3 or more lambs. Ewes that had singles weaned fewer lambs and less body weight (BW) of lambs (P < 0.001; 0.94 lambs and 40.4 kg, respectively) than ewes that had twins or triplets. Ewes with triplet litters in Treatment A weaned more lambs (P < 0.01) and more BW of lambs (P < 0.05) than ewes that had triplets in Treatment R (2.13 lambs and 62.9 kg, respectively, vs. 1.79 lambs and 55.0 kg, respectively), and weaned more lambs than ewes that had twins (1.77 lambs; P < 0.01). However, neither group of triplet-bearing ewes weaned more BW of lambs than ewes that had twins (58.9 kg; P ≥ 0.34). In 2 sets of data involving 442 purebred Polypay litters and 987 litters from Polypay or Romanov‒White Dorper × Rambouillet ewes mated to terminal sires, ewes were required to raise all triplet-born lambs. Death losses for triplets in these studies (39.6 and 31.6%, respectively) were higher than those in Treatment A of the initial study (26.2%), resulting in greater numbers of lambs weaned for triplet, compared to twin, litters (1.79 vs. 1.68, respectively; P = 0.02) but no greater weight of lambs weaned (54.3 vs. 55.4 kg, respectively; P = 0.17). Based on these 3 sets of data, ewes that were required to rear triplet lambs weaned 0.20 more lambs per litter than ewes that had twins but also had 0.75 additional dead lambs per litter, and thus a lamb mortality overhead of 3.75 additional dead lambs for each additional weaned lamb. We conclude that there is an intermediate optimum prolificacy level for extensive rangeland production systems. If optimum prolificacy is exceeded, removal and artificial rearing of surplus lambs are necessary to avoid increased lamb death losses.

Ovar-DRB1 haplotypes *2001 and *0301 are associated with sheep growth and ewe lifetime prolificacy.

BACKGROUND: The major histocompatibility complex (MHC) is an organized cluster of tightly linked vertebrate genes with immunological and non-immunological functions. While the important MHC gene DRB1 has been examined in regard to many sheep infectious disease traits, only one study, based on microsatellite markers, has previously examined DRB1 and sheep production traits. Furthermore, to our knowledge no studies have examined DRB1 relationship with lifetime ewe prolificacy traits. Therefore, we analyzed association between the presence of DRB1 SNP haplotypes with internationally recognized standard names and production traits including growth and lifetime prolificacy in 370 Rambouillet, Columbia, and Polypay sheep. RESULTS: The DRB1 *2001 haplotype was associated with increased weaning and mature weights, as well as average daily gain (Sidák P<0.05; corrected for the number of haplotypes tested). Interestingly, the *2001 haplotype also showed a trend toward association with increased total number of lifetime lambs born (Sidák P=0.084) and number of lambs born alive (Sidák P=0.084). In contrast, the DRB1 *0301 haplotype was associated with decreased mature weight (Sidák P=0.01). CONCLUSIONS: Since the *2001 haplotype was present in all three breeds, these results suggest there is at least one functional mutation in the region that influences growth and prolificacy traits that may be broadly present across several breeds. Furthermore, combined use of the similar *2001 and *0301 multi-marker haplotypes that nonetheless have opposing directions of production trait associations will enhance mutation discovery in this region. If undesirable alleles for underlying mutations can be identified, selective pressure against one or a small number of undesirable alleles may improve production with limited impact on MHC genetic diversity and infectious disease susceptibility.

Classical scrapie prions are associated with peripheral blood monocytes and T-lymphocytes from naturally infected sheep.

BACKGROUND: Classical scrapie is a transmissible spongiform encephalopathy (TSE) that affects sheep and goats. Our previous bioassay studies in lambs revealed that scrapie prions could be detected in association with peripheral blood monocular cells (PBMC), B lymphocytes and platelet-rich plasma fractions. In the present study, bioassay in lambs was again used to determine if scrapie prions are associated with the other two subsets of PBMC, monocytes and T lymphocytes. RESULTS: PBMC, monocytes and T lymphocytes were isolated from two preclinically affected VRQ/VRQ sheep naturally infected with classical ovine scrapie and intravenously transfused into VRQ/VRQ lambs post-weaning. As determined using standard immunohistochemistry for scrapie, abnormal isoforms of prion protein were detected in lymphoid tissues of lambs inoculated with PBMC (4/4 recipient lambs), monocytes (2/5) and T lymphocytes (1/4). Prion protein misfolding activity was detected by serial protein misfolding cyclic amplification (sPMCA) in PBMC from monocyte and T lymphocyte recipient sheep in agreement with antemortem rectal biopsy results, but such prion protein misfolding activity was not detected from other recipients. CONCLUSIONS: These findings show that scrapie prions are associated with monocytes and T lymphocytes circulating in the peripheral blood of sheep naturally infected with classical scrapie. Combined with our previous findings, we can now conclude that all three major subsets of PBMC can harbor prions during preclinical disease and thus, present logical targets for development of a sensitive assay to detect scrapie prions. In this regard, we have also demonstrated that sPMCA can be used to detect scrapie prions associated with PBMC.

Genome-Wide Association Identifies SLC2A9 and NLN Gene Regions as Associated with Entropion in Domestic Sheep.

Entropion is an inward rolling of the eyelid allowing contact between the eyelashes and cornea that may lead to blindness if not corrected. Although many mammalian species, including humans and dogs, are afflicted by congenital entropion, no specific genes or gene regions related to development of entropion have been reported in any mammalian species to date. Entropion in domestic sheep is known to have a genetic component therefore, we used domestic sheep as a model system to identify genomic regions containing genes associated with entropion. A genome-wide association was conducted with congenital entropion in 998 Columbia, Polypay, and Rambouillet sheep genotyped with 50,000 SNP markers. Prevalence of entropion was 6.01%, with all breeds represented. Logistic regression was performed in PLINK with additive allelic, recessive, dominant, and genotypic inheritance models. Two genome-wide significant (empirical P<0.05) SNP were identified, specifically markers in SLC2A9 (empirical P = 0.007; genotypic model) and near NLN (empirical P = 0.026; dominance model). Six additional genome-wide suggestive SNP (nominal P<1x10(-5)) were identified including markers in or near PIK3CB (P = 2.22x10(-6); additive model), KCNB1 (P = 2.93x10(-6); dominance model), ZC3H12C (P = 3.25x10(-6); genotypic model), JPH1 (P = 4.68x20(-6); genotypic model), and MYO3B (P = 5.74x10(-6); recessive model). This is the first report of specific gene regions associated with congenital entropion in any mammalian species, to our knowledge. Further, none of these genes have previously been associated with any eyelid traits. These results represent the first genome-wide analysis of gene regions associated with entropion and provide target regions for the development of sheep genetic markers for marker-assisted selection.

A divergent Artiodactyl MYADM-like repeat is associated with erythrocyte traits and weight of lamb weaned in domestic sheep.

A genome-wide association study (GWAS) was performed to investigate seven red blood cell (RBC) phenotypes in over 500 domestic sheep (Ovis aries) from three breeds (Columbia, Polypay, and Rambouillet). A single nucleotide polymorphism (SNP) showed genome-wide significant association with increased mean corpuscular hemoglobin concentration (MCHC, P = 6.2×10(-14)) and genome-wide suggestive association with decreased mean corpuscular volume (MCV, P = 2.5×10(-6)). The ovine HapMap project found the same genomic region and the same peak SNP has been under extreme historical selective pressure, demonstrating the importance of this region for survival, reproduction, and/or artificially selected traits. We observed a large (>50 kb) variant haplotype sequence containing a full-length divergent artiodactyl MYADM-like repeat in strong linkage disequilibrium with the associated SNP. MYADM gene family members play roles in membrane organization and formation in myeloid cells. However, to our knowledge, no member of the MYADM gene family has been identified in development of morphologically variant RBCs. The specific RBC differences may be indicative of alterations in morphology. Additionally, erythrocytes with altered morphological structure often exhibit increased structural fragility, leading to increased RBC turnover and energy expenditure. The divergent artiodactyl MYADM-like repeat was also associated with increased ewe lifetime kilograms of lamb weaned (P = 2×10(-4)). This suggests selection for normal RBCs might increase lamb weights, although further validation is required before implementation in marker-assisted selection. These results provide clues to explain the strong selection on the artiodactyl MYADM-like repeat locus in sheep, and suggest MYADM family members may be important for RBC morphology in other mammals.

Genome-wide genetic diversity and differentially selected regions among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep.

Sheep are among the major economically important livestock species worldwide because the animals produce milk, wool, skin, and meat. In the present study, the Illumina OvineSNP50 BeadChip was used to investigate genetic diversity and genome selection among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds from the United States. After quality-control filtering of SNPs (single nucleotide polymorphisms), we used 48,026 SNPs, including 46,850 SNPs on autosomes that were in Hardy-Weinberg equilibrium and 1,176 SNPs on chromosome × for analysis. Phylogenetic analysis based on all 46,850 SNPs clearly separated Suffolk from Rambouillet, Columbia, Polypay, and Targhee, which was not surprising as Rambouillet contributed to the synthesis of the later three breeds. Based on pair-wise estimates of F(ST), significant genetic differentiation appeared between Suffolk and Rambouillet (F(ST) = 0.1621), while Rambouillet and Targhee had the closest relationship (F(ST) = 0.0681). A scan of the genome revealed 45 and 41 differentially selected regions (DSRs) between Suffolk and Rambouillet and among Rambouillet-related breed populations, respectively. Our data indicated that regions 13 and 24 between Suffolk and Rambouillet might be good candidates for evaluating breed differences. Furthermore, ovine genome v3.1 assembly was used as reference to link functionally known homologous genes to economically important traits covered by these differentially selected regions. In brief, our present study provides a comprehensive genome-wide view on within- and between-breed genetic differentiation, biodiversity, and evolution among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds. These results may provide new guidance for the synthesis of new breeds with different breeding objectives.

Genome-wide association identifies multiple genomic regions associated with susceptibility to and control of ovine lentivirus.

BACKGROUND: Like human immunodeficiency virus (HIV), ovine lentivirus (OvLV) is macrophage-tropic and causes lifelong infection. OvLV infects one quarter of U.S. sheep and induces pneumonia and body condition wasting. There is no vaccine to prevent OvLV infection and no cost-effective treatment for infected animals. However, breed differences in prevalence and proviral concentration have indicated a genetic basis for susceptibility to OvLV. A recent study identified TMEM154 variants in OvLV susceptibility. The objective here was to identify additional loci associated with odds and/or control of OvLV infection. METHODOLOGY/PRINCIPAL FINDINGS: This genome-wide association study (GWAS) included 964 sheep from Rambouillet, Polypay, and Columbia breeds with serological status and proviral concentration phenotypes. Analytic models accounted for breed and age, as well as genotype. This approach identified TMEM154 (nominal P=9.2×10(-7); empirical P=0.13), provided 12 additional genomic regions associated with odds of infection, and provided 13 regions associated with control of infection (all nominal P<1 × 10(-5)). Rapid decline of linkage disequilibrium with distance suggested many regions included few genes each. Genes in regions associated with odds of infection included DPPA2/DPPA4 (empirical P=0.006), and SYTL3 (P=0.051). Genes in regions associated with control of infection included a zinc finger cluster (ZNF192, ZSCAN16, ZNF389, and ZNF165; P=0.001), C19orf42/TMEM38A (P=0.047), and DLGAP1 (P=0.092). CONCLUSIONS/SIGNIFICANCE: These associations provide targets for mutation discovery in sheep susceptibility to OvLV. Aside from TMEM154, these genes have not been associated previously with lentiviral infection in any species, to our knowledge. Further, data from other species suggest functional hypotheses for future testing of these genes in OvLV and other lentiviral infections. Specifically, SYTL3 binds and may regulate RAB27A, which is required for enveloped virus assembly of human cytomegalovirus. Zinc finger transcription factors have been associated with positive selection for repression of retroviral replication. DLGAP1 binds and may regulate DLG1, a known regulator of HIV infectivity.

Ovine progressive pneumonia provirus levels associate with breed and Ovar-DRB1.

Previous studies initiated defining the role of host genetics in influencing the outcome of exposure to ovine progressive pneumonia virus. However, specific genes influencing host control of virus replication and disease progression have not been identified. This study, using 383 ewes of the Columbia, Polypay, and Rambouillet breeds, tested the hypothesis that host control of OPPV as measured by provirus levels in the peripheral blood associates with certain breeds and MHC class II Ovis aries (Ovar)-DRB1 expressed alleles. Rambouillet ewes were less likely to have measurable provirus levels as compared to Columbia ewes at ages 5 and 6 (P value < 0.02), and they exhibited lower provirus levels when compared to both Columbia and Polypay ewes of the same ages (P value < 0.05). The presence of DRB1*0403- or DRB1*07012-expressed alleles were significantly associated (P value = 0.019 and 0.0002, respectively) with lower OPP provirus levels but only were only found in 11% of the ewe flock. Analysis of each segregating amino acid in the beta1 domain of DR beta-chain revealed that amino acids Y31, T32, N37, T51, Q60, or N74 significantly associated (P value range = 0.0003-0.018) with lower OPP provirus levels, whereas amino acids H32, A38, or I67 associated (P value range = 0.013-0.043) with higher OPP provirus levels. These results suggest that Ovar-DRB1 contributes as one host genetic factor that controls OPP provirus levels, but does not fully account for the breed-specific OPP proviral differences.

Development and validation of an ovine progressive pneumonia virus quantitative PCR.

Ovine progressive pneumonia virus (OPPV) infects at least one sheep in 81% of U.S. sheep flocks, as determined by serology, and can cause viral mastitis, arthritis, dyspnea, and cachexia. Diagnostic tests that quantify OPPV proviral load in peripheral blood leukocytes (PBL) provide an additional method for identification of infected sheep and may help to further understanding of the pathogenesis of OPPV-induced disease. In this study, we compared a new OPPV real-time quantitative PCR (qPCR) assay specific for the transmembrane region of the envelope gene (tm) with a competitive inhibition enzyme-linked immunosorbent assay (cELISA) using 396 PBL samples and sera from Idaho sheep. The OPPV qPCR had a positive concordance of 96.2% +/- 2.3% and a negative concordance of 97.7% +/- 2.5% compared to the cELISA, with a kappa value of 0.93, indicating excellent agreement between the two tests. In addition, the presence of tm in the three OPPV qPCR-positive and cELISA-negative sheep and in 15 sheep with different OPPV proviral loads was confirmed by cloning and sequencing. These data indicate that the OPPV qPCR may be used as a supplemental diagnostic tool for OPPV infection and for measurement of viral load in PBLs of infected sheep.

Enhancement of equine infectious anemia virus virulence by identification and removal of suboptimal nucleotides.

Pathogenicity was reportedly restored to an avirulent molecular clone of equine infectious anemia virus (EIAV) by substitution of 3' sequences from the pathogenic variant strain (EIAV(PV)). However, the incidence of disease in horses/ponies was found to be significantly lower (P = 0.016) with the chimeric clone (EIAV(UK)) than with EIAV(PV). This was attributable to 3' rather than 5' regions of the proviral genome, where EIAV(UK) differs from the consensus EIAV(PV) sequence by having a 68-bp duplication in the 3' LTR and arginine (R(103)) rather than tryptophan (W(103)) at position 103 in the second exon of rev. In EIAV(UK) recipients the duplication was rapidly eliminated and R(103) replaced by W(103) in the viral population. Furthermore, removal of the 3' variant sequences from EIAV(UK) (EIAV(UK3)) resulted in an equivalent (P = 0.013) disease potential in Equus caballus to EIAV(PV). The 68-bp duplication and/or R(103) may limit peak viral RNA accumulation during acute infection.