Equine arteritis virus (EAV) is the causative agent of equine viral arteritis (EVA), a respiratory, systemic, and reproductive disease of horses and other equid species. Following natural infection, 10–70% of the infected stallions can become persistently infected and continue to shed EAV in their semen for periods ranging from several months to life. Recently, we reported that some stallions possess a subpopulation(s) of CD3+ T lymphocytes that are susceptible to in vitro EAV infection and that this phenotypic trait is associated with long-term carrier status following exposure to the virus. In contrast, stallions not possessing the CD3+ T lymphocyte susceptible phenotype are at less risk of becoming long-term virus carriers. A genome wide association study (GWAS) using the Illumina Equine SNP50 chip revealed that the ability of EAV to infect CD3+ T lymphocytes and establish long-term carrier status in stallions correlated with a region within equine chromosome 11. Here we identified the gene and mutations responsible for these phenotypes. Specifically, the work implicated three allelic variants of the equine orthologue of CXCL16 (EqCXCL16) that differ by four non-synonymous nucleotide substitutions (XM_00154756; c.715 A → T, c.801 G → C, c.804 T → A/G, c.810 G → A) within exon 1. This resulted in four amino acid changes with EqCXCL16S (XP_001504806.1) having Phe, His, Ile and Lys as compared to EqCXL16R having Tyr, Asp, Phe, and Glu at 40, 49, 50, and 52, respectively. Two alleles (EqCXCL16Sa, EqCXCL16Sb) encoded identical protein products that correlated strongly with long-term EAV persistence in stallions (P<0.000001) and are required for in vitro CD3+ T lymphocyte susceptibility to EAV infection. The third (EqCXCL16R) was associated with in vitro CD3+ T lymphocyte resistance to EAV infection and a significantly lower probability for establishment of the long-term carrier state (viral persistence) in the male reproductive tract. EqCXCL16Sa and EqCXCL16Sb exert a dominant mode of inheritance. Most importantly, the protein isoform EqCXCL16S but not EqCXCL16R can function as an EAV cellular receptor. Although both molecules have equal chemoattractant potential, EqCXCL16S has significantly higher scavenger receptor and adhesion properties compared to EqCXCL16R.

Skeletal atavism in Shetland ponies is a heritable disorder characterized by abnormal growth of the ulna and fibula that extend the carpal and tarsal joints, respectively. This causes abnormal skeletal structure and impaired movements, and affected foals are usually killed. In order to identify the causal mutation we subjected six confirmed Swedish cases and a DNA pool consisting of 21 control individuals to whole genome resequencing. We screened for polymorphisms where the cases and the control pool were fixed for opposite alleles and observed this signature for only 25 SNPs, most of which were scattered on genome assembly unassigned scaffolds. Read depth analysis at these loci revealed homozygosity or compound heterozygosity for two partially overlapping large deletions in the pseudoautosomal region (PAR) of chromosome X/Y in cases but not in the control pool. One of these deletions removes the entire coding region of the SHOX gene and both deletions remove parts of the CRLF2 gene located downstream of SHOX. The horse reference assembly of the PAR is highly fragmented, and in order to characterize this region we sequenced bacterial artificial chromosome (BAC) clones by single-molecule real-time (SMRT) sequencing technology. This considerably improved the assembly and enabled size estimations of the two deletions to 160−180 kb and 60−80 kb, respectively. Complete association between the presence of these deletions and disease status was verified in eight other affected horses. The result of the present study is consistent with previous studies in humans showing crucial importance of SHOX for normal skeletal development.

Previous studies in our laboratory have identified equine CXCL16 (EqCXCL16) to be a candidate molecule and possible cell entry receptor for equine arteritis virus (EAV). In horses, the CXCL16 gene is located on equine chromosome 11 (ECA11) and encodes a glycosylated, type I transmembrane protein with 247 amino acids. Stable transfection of HEK-293T cells with plasmid DNA carrying EqCXCL16 (HEK-EqCXCL16 cells) increased the proportion of the cell population permissive to EAV infection from <3% to almost 100%. The increase in permissiveness was blocked either by transfection of HEK-EqCXCL16 cells with small interfering RNAs (siRNAs) directed against EqCXCL16 or by pretreatment with guinea pig polyclonal antibody against EqCXCL16 protein (Gp anti-EqCXCL16 pAb). Furthermore, using a virus overlay protein-binding assay (VOPBA) in combination with far-Western blotting, gradient-purified EAV particles were shown to bind directly to the EqCXCL16 protein in vitro. The binding of biotinylated virulent EAV strain Bucyrus at 4°C was significantly higher in HEK-EqCXCL16 cells than nontransfected HEK-293T cells. Finally, the results demonstrated that EAV preferentially infects subpopulations of horse CD14؉ monocytes expressing EqCXCL16 and that infection of these cells is significantly reduced by pretreatment with Gp anti-EqCXCL16 pAb. The collective data from this study provide confirmatory evidence that the transmembrane form of EqCXCL16 likely plays a major role in EAV host cell entry processes, possibly acting as a primary receptor molecule for this virus. IMPORTANCEOutbreaks of EVA can be a source of significant economic loss for the equine industry from high rates of abortion in pregnant mares, death in young foals, establishment of the carrier state in stallions, and trade restrictions imposed by various countries. Similar to other arteriviruses, EAV primarily targets cells of the monocyte/macrophage lineage, which, when infected, are believed to play a critical role in EVA pathogenesis. To this point, however, the host-specified molecules involved in EAV binding and entry into monocytes/macrophages have not been identified. Identification of the cellular receptors for EAV may provide insights to design antivirals and better prophylactic reagents. In this study, we have demonstrated that EqCXCL16 acts as an EAV entry receptor in EAV-susceptible cells, equine monocytes. These findings represent a significant advance in our understanding of the fundamental mechanisms associated with the entry of EAV into susceptible cells.

No abstract

Chondrodysplastic dwarfism in Miniature horses appeared to be a recessive genetic trait based on the occurrence of affected offspring by normal parents. Dwarf phenotypes vary and range from abnormal abortuses to viable offspring with evidence of skeletal dysplasia. A genome-wide association study implicated a region of ECA1 with dwarfism in Miniature horses. Aggrecan (ACAN) was a candidate gene in that region, and exons were sequenced to compare DNA sequences for dwarf and non-dwarf horses. Sequencing led to the discovery of variants in exons 2, 6, 7 and 15 associated with dwarfism. The four variants are identified with reference to Ecab 3.0 (GCF_002863925.1) as g.95291270del (rs1095048841), g.95284530C>T (ERP107353), g.95282140C>G (rs1095048823) and g.95257480_95257500del (rs1095048839) and designated here as D1, D2, D3* and D4 respectively. A previous study at another laboratory reported dwarfism associated with homozygosity for D3*. Homozygotes for those variants and compound heterozygotes for any combination of those variants always expressed a dwarfism phenotype. However, eight additional horses with dwarfism were found, seven of which were heterozygotes for D2, D3* or D4, suggesting the existence of additional ACAN alleles causing dwarfism. Among Miniature horses, the combined frequency of D1, D2, D3* and D4 was 0.163, suggesting a carrier rate of 26.2% for alleles causing chondrodysplastic dwarfism.

Mushroom is a unique coat color phenotype in Shetland Ponies characterized by the dilution of the chestnut coat color to a sepia tone and is hypothesized to be a recessive trait. A genome wide association study (GWAS), utilizing the Affymetrix 670K array (MNEc670k) and a single locus mixed linear model analysis (EMMAX), identified a locus on ECA7 for further investigation (Pcorrected = 2.08 × 10−10). This locus contained a 3 Mb run of homozygosity in the 12 mushroom ponies tested. Analysis of high throughput Illumina sequencing data from one mushroom Shetland pony compared to 87 genomes from horses of various breeds, uncovered a frameshift variant, p.Asp201fs, in the MFSD12 gene encoding the major facilitator superfamily domain containing 12 protein. This variant was perfectly concordant with phenotype in 96 Shetland Ponies (P = 1.15 × 10−22), was identified in the closely related Miniature Horse for which the mushroom phenotype is suspected to occur (fmu = 0.02), and was absent in 252 individuals from seven additional breeds not reported to have the mushroom phenotype. MFSD12 is highly expressed in melanocytes and variants in this gene in humans, mice, and dogs impact pigmentation. Given the role of MFSD12 in melanogenesis, we propose that p.Asp201fs is causal for the dilution observed in mushroom ponies.

No abstract

Horse owners choose whether or not to train American Saddlebred horses (ASHs) to perform the 4-beat gaits called rack and slow gait. The rack and slow gait are similar to ambling gaits shown to be associated with variation in the DMRT3 gene in other breeds but are trained rather than naturally occurring gaits. A premature stop codon in the DMRT3 gene (DMRT3_Ser301STOP) caused by the Ch23:g.22999655C>A SNP has an effect on the pattern of locomotion in horses and allows for the pacing gait and strong association with performance of ambling gaits in diverse breeds. We used horse show records to identify ASHs trained to perform as 5-gaited horses and compared them to other Saddlebreds for the prevalence of the A allele of the Ch23:g.22999655C>A SNP of DMRT3. Genomic DNA was typed using a PCR-RFLP technique for 5-gaited ASHs (n = 55), 3-gaited or harness ASHs (n = 64) and ASHs not shown (n = 128). We did not observe differences in the frequencies of the Ch23:g.22999655C>A SNP of DMRT3 when comparing 5-gaited ASHs to other ASHs (P > 0.22). These results suggest that DMRT3 alleles do not affect the choice of breeders to train and show ASHs as 5-gaited horses.