A genome-wide association study for quantitative trait loci of show-jumping in Hanoverian warmblood horses

Show-jumping is an economically important breeding goal in Hanoverian warmblood horses. The aim of this study was a genome-wide association study (GWAS) for quantitative trait loci (QTL) for show-jumping in Hanoverian warmblood horses, employing the Illumina equine SNP50 Beadchip. For our analyses, we genotyped 115 stallions of the National State stud of Lower Saxony. The show-jumping talent of a horse includes style and ability in free-jumping. To control spurious associations based on population stratification, two different mixed linear animal model (MLM) approaches were employed, besides linear models with fixed effects only and adaptive permutations for correcting multiple testing. Population stratification was explained best in the MLM considering Hanoverian, Thoroughbred, Trakehner and Holsteiner genes and the marker identity-by-state relationship matrix. We identified six QTL for show-jumping on horse chromosomes (ECA) 1, 8, 9 and 26 (-log(10) P-value >5) and further putative QTL with -log(10) P-values of 3-5 on ECA1, 3, 11, 17 and 21. Within six QTL regions, we identified human performance-related genes including PAPSS2 on ECA1, MYL2 on ECA8, TRHR on ECA9 and GABPA on ECA26 and within the putative QTL regions NRAP on ECA1, and TBX4 on ECA11. The results of our GWAS suggest that genes involved in muscle structure, development and metabolism are crucial for elite show-jumping performance. Further studies are required to validate these QTL in larger data sets and further horse populations.     

A Non-Synonymous HMGA2 Variant Decreases Height in Shetland Ponies and Other Small Horses

The identification of quantitative trait loci (QTL) such as height and their underlying causative variants is still challenging and often requires large sample sizes. In humans hundreds of loci with small effects control the heritable portion of height variability. In domestic animals, typically only a few loci with comparatively large effects explain a major fraction of the heritability. We investigated height at withers in Shetland ponies and mapped a QTL to ECA 6 by genome-wide association (GWAS) using a small cohort of only 48 animals and the Illumina equine SNP70 BeadChip. Fine-mapping revealed a shared haplotype block of 793 kb in small Shetland ponies. The HMGA2 gene, known to be associated with height in horses and many other species, was located in the associated haplotype. After closing a gap in the equine reference genome we identified a non-synonymous variant in the first exon of HMGA2 in small Shetland ponies. The variant was predicted to affect the functionally important first AT-hook DNA binding domain of the HMGA2 protein (c.83G>A; p.G28E). We assessed the functional impact and found impaired DNA binding of a peptide with the mutant sequence in an electrophoretic mobility shift assay. This suggests that the HMGA2 variant also affects DNA binding in vivo and thus leads to reduced growth and a smaller stature in Shetland ponies. The identified HMGA2 variant also segregates in several other pony breeds but was not found in regular-sized horse breeds. We therefore conclude that we identified a quantitative trait nucleotide for height in horses.     

Evidence that In Vitro Susceptibility of CD3+ T Lymphocytes to Equine Arteritis Virus Infection Reflects Genetic Predisposition of Naturally Infected Stallions To Become Carriers of the Virus

We investigated the correlation between in vitro susceptibility of CD3⁺ T lymphocytes to equine arteritis virus (EAV) infection and establishment of persistent infection among 14 stallions following natural infections. The data showed that carrier stallions with a CD3⁺ T lymphocyte susceptibility phenotype to in vitro EAV infection may be at higher risk of becoming carriers than those that lack this phenotype (P = 0.0002).     

A Genome-Wide Association Study Identifies Risk Loci to Equine Recurrent Uveitis in German Warmblood Horses

Equine recurrent uveitis (ERU) is a common eye disease affecting up to 3–15% of the horse population. A genome-wide association study (GWAS) using the Illumina equine SNP50 bead chip was performed to identify loci conferring risk to ERU. The sample included a total of 144 German warmblood horses. A GWAS showed a significant single nucleotide polymorphism (SNP) on horse chromosome (ECA) 20 at 49.3 Mb, with IL-17A and IL-17F being the closest genes. This locus explained a fraction of 23% of the phenotypic variance for ERU. A GWAS taking into account the severity of ERU, revealed a SNP on ECA18 nearby to the crystalline gene cluster CRYGA-CRYGF. For both genomic regions on ECA18 and 20, significantly associated haplotypes containing the genome-wide significant SNPs could be demonstrated. In conclusion, our results are indicative for a genetic component regulating the possible critical role of IL-17A and IL-17F in the pathogenesis of ERU. The associated SNP on ECA18 may be indicative for cataract formation in the course of ERU.     

Macrophages and CD4+ T lymphocytes from two multiply exposed, uninfected individuals resist infection with primary non-syncytium-inducing isolates of human immunodeficiency virus type 1.

Despite multiple, high-risk sexual exposures, some individuals remain uninfected with human immunodeficiency virus type 1 (HIV-1). CD4+ lymphocytes from these individuals are less susceptible to infection in vitro with some strains of HIV-1, suggesting that the phenotype of the virus may influence its ability to interact with certain CD4+ cells. In the present study, we examined the susceptibility of CD4+ T lymphocytes and macrophages from two exposed uninfected individuals (EU2 and EU3) to infection with a panel of biologically cloned isolates of HIV-1 having either a non-syncytium-inducing (NSI) or a syncytium-inducing (SI) phenotype. Our results indicate that CD4+ T lymphocytes from EU2 and EU3 are resistant to infection with NSI isolates of HIV-1 but are susceptible to infection with primary SI isolates. In addition, we found that macrophages from EU2 and EU3 are resistant to infection with both NSI and SI isolates. The latter finding was confirmed by using several uncloned NSI and SI isolates obtained from patients during acute HIV-1 infection. In further experiments, env clones encoding glycoproteins characteristic of NSI or SI viruses were used in single-cycle infectivity assays to evaluate infection of CD4+ lymphocytes and macrophages from EU2 and EU3. Consistent with our previous results, we found that macrophages from these individuals are resistant to infection with NSI and SI env-pseudotyped viruses, while CD4+ T lymphocytes are resistant to NSI, but not SI, pseudotyped viruses. Overall, our results demonstrate that CD4+ cells from two exposed uninfected individuals resist infection in vitro with primary, macrophage-tropic, NSI isolates of HIV-1, which is the predominant viral phenotype found following HIV-1 transmission. Furthermore, infection with NSI isolates was blocked in both CD4+ T lymphocytes and macrophages from these individuals, suggesting that there may be a common mechanism for resistance in both cell types.     

Fine mapping of a quantitative trait locus for osteochondrosis on horse chromosome 2

In this study, we refine a quantitative trait locus for equine osteochondrosis (OC) on horse chromosome (ECA) 2 to a genome-wide significant interval at 20.08-30.94 Mb. The marker set contained 27 newly developed microsatellites equidistantly distributed over ECA2 and 44 nucleotide polymorphisms, located in 16 positional candidate genes for OC. Genotyping was performed in 211 Hanoverian horses from 14 paternal half-sib groups. A NCDN-associated SNP and haplotype were significantly associated with OC in fetlock and/or hock joints. This study is a further step towards the identification of genes responsible for OC in horses.     

Concerted evolution of vertebrate CCR2 and CCR5 genes and the origin of a recombinant equine CCR5/2 gene

Chemokine receptors (CCRs) play an essential role in the initiation of an innate immune host response. Several of these receptors have been shown to modulate the outcome of viral infections. The recent availability of complete genome sequences from a number of species provides a unique opportunity to analyze the evolution of the CCR genes. A phylogenetic analysis revealed that the CCR2 gene evolved in concert with the paralogous CCR5 gene, but not with another paralogous gene, CCR3, in the opossum, platypus, rabbit, guinea pig, cat, and rodent lineages. In addition, evidence of concerted evolution of the CCR2 and CCR5 genes was observed in chicken and lizard genomes. A unique CCR5/2 gene that originated by unequal crossing over between the CCR2 and CCR5 genes was detected in the domestic horse. The CCR2, CCR5, and CCR5/2 genes were mapped to ECA16q21 using fluorescent in situ hybridization (FISH). Single-nucleotide polymorphisms identified in the equine CCR5 gene and characterized within 5 horse breeds provide haplotype markers for future case/control studies investigating the genetic bases of horse susceptibility to infectious diseases.     

Whole-genome scan identifies quantitative trait loci for chronic pastern dermatitis in German draft horses

Chronic pastern dermatitis (CPD), also known as chronic progressive lymphedema (CPL), is a skin disease that affects draft horses. This disease causes painful lower-leg swelling, nodule formation, and skin ulceration, interfering with movement. The aim of this whole-genome scan was to identify quantitative trait loci (QTL) for CPD in German draft horses. We recorded clinical data for CPD in 917 German draft horses and collected blood samples from these horses. Of these 917 horses, 31 paternal half-sib families comprising 378 horses from the breeds Rhenish German, Schleswig, Saxon-Thuringian, and South German were chosen for genotyping. Each half-sib family was constituted by only one draft horse breed. Genotyping was done for 318 polymorphic microsatellites evenly distributed on all equine autosomes and the X chromosome with a mean distance of 7.5 Mb. An across-breed multipoint linkage analysis revealed chromosome-wide significant QTL on horse chromosomes (ECA) 1, 9, 16, and 17. Analyses by breed confirmed the QTL on ECA1 in South German and the QTL on ECA9, 16, and 17 in Saxon-Thuringian draft horses. For the Rhenish German and Schleswig draft horses, additional QTL on ECA4 and 10 and for the South German draft horses an additional QTL on ECA7 were found. This is the first whole-genome scan for CPD in draft horses and it is an important step toward the identification of candidate genes.     

An Overview of Ten Italian Horse Breeds through Mitochondrial DNA

BackgroundThe climatic and cultural diversity of the Italian Peninsula triggered, over time, the development of a great variety of horse breeds, whose origin and history are still unclear. To clarify this issue, analyses on phenotypic traits and genealogical data were recently coupled with molecular screening.MethodologyTo provide a comprehensive overview of the horse genetic variability in Italy, we produced and phylogenetically analyzed 407 mitochondrial DNA (mtDNA) control-region sequences from ten of the most important Italian riding horse and pony breeds: Bardigiano, Esperia, Giara, Lipizzan, Maremmano, Monterufolino, Murgese, Sarcidano, Sardinian Anglo-Arab, and Tolfetano. A collection of 36 Arabian horses was also evaluated to assess the genetic consequences of their common use for the improvement of some local breeds.ConclusionsIn Italian horses, all previously described domestic mtDNA haplogroups were detected as well as a high haplotype diversity. These findings indicate that the ancestral local mares harbored an extensive genetic diversity. Moreover, the limited haplotype sharing (11%) with the Arabian horse reveals that its impact on the autochthonous mitochondrial gene pools during the final establishment of pure breeds was marginal, if any. The only significant signs of genetic structure and differentiation were detected in the geographically most isolated contexts (i.e. Monterufolino and Sardinian breeds). Such a geographic effect was also confirmed in a wider breed setting, where the Italian pool stands in an intermediate position together with most of the other Mediterranean stocks. However, some notable exceptions and peculiar genetic proximities lend genetic support to historical theories about the origin of specific Italian breeds.     

Genetic divergence with emergence of novel phenotypic variants of equine arteritis virus during persistent infection of stallions

The persistently infected carrier stallion is the critical natural reservoir of equine arteritis virus (EAV), as venereal infection of mares frequently occurs after breeding to such stallions. Two Thoroughbred stallions that were infected during the 1984 outbreak of equine viral arteritis in central Kentucky subsequently became long-term EAV carriers. EAV genomes amplified from the semen of these two stallions were compared by sequence analysis of the six 3' open reading frames (ORFs 2 through 7), which encode the four known structural proteins and two uncharacterized glycoproteins. The major variants of the EAV population that sequentially arose within the reproductive tract of each carrier stallion varied by approximately 1% per year, and the heterogeneity of the viral quasispecies increased during the course of long-term persistent infection. The various ORFs of the dominant EAV variants evolved independently, and there was apparently strong selective pressure on the uncharacterized GP3 protein during persistent infection. Amino acid changes also occurred in the V1 variable region of the GL protein. This region has been previously identified as a crucial neutralization domain, and selective pressures exerted on the V1 region during persistent EAV infection led to the emergence of virus variants with distinct neutralization properties. Thus, evolution of the EAV quasispecies that occurs during persistent infection of the stallion clearly can influence viral phenotypic properties such as neutralization and perhaps virulence.     

Refinement of a quantitative trait locus on equine chromosome 5 responsible for fetlock osteochondrosis in Hanoverian warmblood horses

In this report, we provide 29 new informative microsatellites distributed over a region of 21 Mb on horse chromosome (ECA) 5 and refine a quantitative trait locus (QTL) for fetlock osteochondrosis dissecans (OCD) to a genome-wide significant interval between 78.03 and 90.23 Mb on ECA5. Genotyping was performed in 211 Hanoverian warmblood horses from 14 paternal half-sib groups. Within this OCD-QTL, collagen type XXIV alpha 1 was identified as a potential functional candidate gene for equine osteochondrosis. This report is a further step towards unravelling the genes that cause equine osteochondrosis.     

FISH mapping of the IGF2 gene in horse and donkey—detection of homoeology with HSA11

Three genomic subclones derived from a phage clone containing the equine IGF2 gene were used to FISH map the gene on horse (ECA) and donkey (EAS) metaphase chromosomes. The gene mapped on ECA 12q13 band and is the first locus mapped to this horse chromosome. In donkey the gene mapped very terminal on the long arm of one small submetacentric chromosome that shows almost identical DAPI-banding pattern with ECA12. This is the first locus mapped in donkey genome. Cross species chromosome painting of equine metaphase chromosomes with human Chromosome (Chr) 11-specific probe showed homoeology of this human chromosome with ECA12 and ECA7. The novel ECA12 comparative painting results are thus in accordance with the localization of the equine IGF2 gene. Comparison of the hitherto known physical locations of IGF2 in different species, viz. human, cattle, sheep, horse, and donkey, shows that this gene tends to maintain a terminal location on the chromosome arm. Received: 12 January 1997 / Accepted: 17 March 1997     

Effect of polymorphisms in the Slc11a1 coding region on resistance to brucellosis by macrophages in vitro and after challenge in two Bos breeds (Blanco Orejinegro and Zebu)

The resistance/susceptibility of selected cattle breeds to brucellosis was evaluated in an F1 population generated by crossing animals classified as resistant (R) and susceptible (S) (R x R, R x S, S x R, S x S) based on challenges in vitro and in vivo. The association between single nucleotide polymorphisms identified in the coding region of the Slc11a1 gene and resistance/susceptibility was estimated. The trait resistance or susceptibility to brucellosis, evaluated by a challenge in vitro, showed a high heritable component in terms of additive genetic variance (h(2) = 0.54 ± 0.11). In addition, there was a significant association (p < 0.05) between the control of bacterial survival and two polymorphisms (a 3'UTR and SNP4 located in exon 10). The antibody response of animals classified as resistant to infection by Brucella abortus differed significantly (p < 0.05) from that of susceptible animals. However, there was no significant association between single nucleotide polymorphisms located in the Slc11a1 gene and the antibody response stimulated by a challenge in vivo.     

Characterization of equine and other vertebrate TLR3, TLR7, and TLR8 genes

Toll-like receptors 3, 7, and 8 (TLR3, TLR7, and TLR8) were studied in the genomes of the domestic horse and several other mammals. The messenger RNA sequences and exon/intron structures of these TLR genes were determined. An equine bacterial artificial chromosome clone containing the TLR3 gene was assigned by fluorescent in situ hybridization to the horse chromosomal location ECA27q16–q17 and this map location was confirmed using an equine radiation hybrid panel. Direct sequencing revealed 13 single-nucleotide polymorphisms in the coding regions of the equine TLR 3, 7, and 8 genes. Of these polymorphisms, 12 were not previously reported. The allelic frequency was estimated for each single-nucleotide polymorphism from genotyping data obtained for 154 animals from five horse breeds. Some of these frequencies varied significantly among different horse breeds. Domain architecture predictions for the three equine TLR protein sequences revealed several conserved regions within the variable leucine-rich repeats between the corresponding horse and cattle TLR proteins. A phylogenetic analysis did not indicate that any significant exchanges had occurred between paralogous TLR7 and TLR8 genes in 20 vertebrate species analyzed. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Whole-genome scan for guttural pouch tympany in Arabian and German warmblood horses

Equine guttural pouch tympany (GPT) is a hereditary disease in foals of several breeds, including thoroughbreds, Arabian, Quarter and warmblood horses. We performed a whole-genome scan for GPT in 143 horses from five Arabian and five German warmblood families and genotyped 257 microsatellites. Chromosome-wide significant linkage was detected on ECA2 and ECA15 using multipoint non-parametric linkage analyses. Analyses stratified by sex revealed chromosome-wide significant linkage on ECA2 for fillies and chromosome-wide significant linkage on ECA15 for colts. For Arabian colts, the quantitative trait locus (QTL) on ECA15 was genome-wide significant. Haplotypes including two to four microsatellites within the QTL on ECA2 and 15 in fillies and colts, respectively, were significantly associated with GPT for both breeds. Thus, our analysis indicated sex-specific QTL, a fact which is in agreement with a two- to fourfold higher incidence of GPT in females. This is the first report of QTL for equine GPT and a first step towards identifying genes responsible for GPT.     

A genome-wide association study for racing performances in Thoroughbreds clarifies a candidate region near the MSTN gene

Using 1400 microsatellites, a genome-wide association study (GWAS) was performed to identify genomic regions associated with lifetime earnings and performance ranks, as determined by the Japan Racing Association (JRA). The minimum heritability (h(2) ) was estimated at 7-8% based on the quantitative trait model, suggesting that the racing performance is heritable. Following GWAS with microsatellites, fine mapping led to identification of three SNPs on ECA18, namely, g.65809482T>C (P=1.05E-18), g.65868604G>T (P=6.47E-17), and g.66539967A>G (P=3.35E-14) associated with these performance measures. The haplotype of these SNPs, together with a recently published nearby SNP, g.66493737C>T (P=9.06E-16) in strong linkage disequilibrium, also showed a very clear association with the performance (P<1E-05). The candidate genomic region contained eight genes annotated by ENSEMBL, including the myostatin gene (MSTN). These findings suggest the presence of a gene affecting the racing performance in Thoroughbred racehorses in this region on ECA18.     

Antibodies against equine herpesviruses and equine arteritis virus in Burchell's zebras (Equus burchelli ) from the Serengeti ecosystem

A total of 51 sera from a migratory population of Burchell's zebras (Equus burchelli) were collected in the Serengeti National Park (Tanzania) between 1999 and 2001 to assess levels of exposure to equine herpesvirus types 1, 2, 4, 9 (EHV-1, -2, -4, -9), EHV-1 zebra isolate T965, and equine arteritis virus (EAV). Using virus-specific neutralizing antibody tests, seroprevalence was high for EHV-9 (60% of 45), moderate for EAV (24% of 51), and lower for the EHV-1-related zebra isolate (17% of 41), EHV-1 (14% of 49), and EHV-4 (2% of 50). No evidence for exposure to EHV-2 was found (0% of 51). The high level of exposure to EHV-9 is interesting because evidence of infection with this virus has not been previously described in any wild equine population. Although the epidemiology of EHV-9 in Burchell's zebras is presently unknown, our results suggest that in East Africa, this species may be a natural host of EHV-9, a neuropathogenic virus that was only recently isolated from captive Thomson's gazelles (Gazella thomsoni) in Japan. There is currently no evidence that EHV-9 induced mortality in Burchell's zebras in the Serengeti, but because of the reported virulence of this virus for more susceptible species such as Thomson's gazelles, viral transmission from infected zebras to ungulates may result in mortality.     

The equine CD1 gene family is the largest and most diverse yet identified

The CD1 family is a group of non-polymorphic MHC class I-like molecules that present lipid-based antigens to T cells. Previous work in our laboratory demonstrated that cytotoxic T lymphocytes from immune adult horses recognize lipids from the cell wall of an important equine pathogen, Rhodococcus equi. These findings suggest an important role for the equine CD1 antigen presentation system in protective immune responses to microbial pathogens in the horse. In this study, we characterized and mapped the equine CD1 gene cluster. The equine genome was found to contain 13 complete CD1 genes; seven genes were classified as homologues of human CD1a, two CD1b, one CD1c, one CD1d, and two CD1e, making it the largest CD1 family to date. All but one of the eqCD1 molecules were expressed in all antigen-presenting cells investigated. The major amino acid differences between equine CD1 isoforms are located in the predicted antigen binding site, suggesting that a variety of lipid antigens can be presented. R. equi survives and replicates within professional phagocytes by arresting phagosome maturation between the early endosome and late phagosome. Based on the absence of a tyrosine sorting motif in all eqCD1a, CD1a molecules are predicted to co-localize with R. equi in the early endosome. Here, they could acquire lipid antigen and present it to T lymphocytes. The extraordinarily large number of CD1 molecules in the horse may reflect their crucial role in immunity to R. equi.     

The 5'' end of the equine arteritis virus replicase gene encodes a papainlike cysteine protease.

The presence of a papainlike cysteine protease (PCP) domain in the N-terminal region of the equine arteritis virus (EAV) replicase, which had been postulated on the basis of limited sequence similarities with cellular and viral thiol proteases, was confirmed by in vitro translation and mutagenesis studies. The EAV protease was found to direct an autoproteolytic cleavage at its C terminus which leads to the production of an approximately 30-kDa N-terminal replicase product (nsp1) containing the PCP domain. Amino acid residues Cys-164 and His-230 of the EAV replicase polyprotein were identified as the most likely candidates for the role of PCP catalytic residues. By means of N-terminal sequence analysis of a PCP cleavage product, derived from a bacterial expression system, it was shown that cleavage occurs between Gly-260 and Gly-261. No evidence for PCP-directed cleavages at other positions in the EAV replicase was obtained. In cotranslational and posttranslational trans-cleavage assays, neither EAV nsp1 nor its precursor was able to process the PCP cleavage site in trans.