Copy number variation in the horse genome.

We constructed a 400K WG tiling oligoarray for the horse and applied it for the discovery of copy number variations (CNVs) in 38 normal horses of 16 diverse breeds, and the Przewalski horse. Probes on the array represented 18,763 autosomal and X-linked genes, and intergenic, sub-telomeric and chrY sequences. We identified 258 CNV regions (CNVRs) across all autosomes, chrX and chrUn, but not in chrY. CNVs comprised 1.3% of the horse genome with chr12 being most enriched. American Miniature horses had the highest and American Quarter Horses the lowest number of CNVs in relation to Thoroughbred reference. The Przewalski horse was similar to native ponies and draft breeds. The majority of CNVRs involved genes, while 20% were located in intergenic regions. Similar to previous studies in horses and other mammals, molecular functions of CNV-associated genes were predominantly in sensory perception, immunity and reproduction. The findings were integrated with previous studies to generate a composite genome-wide dataset of 1476 CNVRs. Of these, 301 CNVRs were shared between studies, while 1174 were novel and require further validation. Integrated data revealed that to date, 41 out of over 400 breeds of the domestic horse have been analyzed for CNVs, of which 11 new breeds were added in this study. Finally, the composite CNV dataset was applied in a pilot study for the discovery of CNVs in 6 horses with XY disorders of sexual development. A homozygous deletion involving AKR1C gene cluster in chr29 in two affected horses was considered possibly causative because of the known role of AKR1C genes in testicular androgen synthesis and sexual development. While the findings improve and integrate the knowledge of CNVs in horses, they also show that for effective discovery of variants of biomedical importance, more breeds and individuals need to be analyzed using comparable methodological approaches.

Parascaris univalens--a victim of large-scale misidentification?

The equine ascarid parasite Parascaris equorum is well known as a ubiquitous parasite infecting foals. A sibling species, Parascaris univalens, was first described over 130 years ago, but very little attention has been given to its existence and possible implications for anthelmintic resistance, clinical disease, or host age spectrum. P. univalens only possesses one germ line chromosome pair as opposed to two for P. equorum, but the two species are otherwise considered morphologically identical. For the present study, live worms obtained from the University of Kentucky parasitology horse herd were dissected and identified using karyotyping techniques. With no exception, all specimens (n = 30) were identified to be P. univalens. Further, the karyotyping technique was adapted to ascarid eggs derived from fecal samples and carried out on samples collected from 25 Thoroughbred foals from three farms in Central Kentucky. P. equorum was not identified among these, whereas P. univalens was found in 17 samples, with the remaining being inconclusive. The mitochondrial genome was sequenced, assembled, and annotated from one male worm identified as P. univalens, and comparison with available sequence reads labeled as P. equorum revealed only 0.16% nucleotide differences. However, it is unlikely that the sequences available in public databases have been unequivocally identified to species level by karyotyping. Taken together, these data suggest that P. univalens is likely the main species now observed in equines and that perhaps the designation Parascaris spp. should be used unless cytological characterization has confirmed the species.

Detection of two equine trisomies using SNP-CGH.

Chromosomal aberrations in the horse are known to cause congenital abnormalities, embryonic loss, and infertility. While diagnosed mainly by karyotyping and FISH in the horse, the use of SNP array comparative genome hybridization (SNP-CGH) is becoming increasingly common in human diagnostics. Normalized probe intensities and allelic ratios are used to detect changes in copy number genome-wide. Two horses with suspected chromosomal abnormalities and six horses with FISH-confirmed aberrant karyotypes were chosen for genotyping on the Equine SNP50 array. Karyotyping of the first horse indicated mosaicism for an additional small, acrocentric chromosome, although the identity of the chromosome was unclear. The second case displayed a similar phenotype to human disease caused by a gene deletion and so was chosen for SNP-CGH due to the ability to detect changes at higher resolutions than those achieved with conventional karyotyping. The results of SNP-CGH analysis for the six horses with known chromosomal aberrations agreed completely with previous karyotype and FISH analysis. The first undiagnosed case showed a pattern of altered allelic ratios without a noticeable shift in overall intensity for chromosome 27, consistent with a mosaic trisomy. The second case displayed a more drastic change in both values for chromosome 30, consistent with a complete trisomy. These results indicate that SNP-CGH is a viable method for detection of chromosomal aneuploidies in the horse.

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.

The effect of age and telomere length on immune function in the horse.

Telomeres, specialized structures present at the ends of linear eukaryotic chromosomes, function to maintain chromosome stability and integrity. Telomeres shorten with each cell division eventually leading to replicative senescence, a process thought to be associated with age-related decline in immune function. We hypothesized that shortened PBMC telomere length is a factor contributing to immunosenescence of the aged horse. Telomere length was assessed in 19 horses ranging in age from 1 to 25 years. Mitogen-induced 3H-thymidine incorporation, total serum IgG, and pro-inflammatory cytokine expression was also determined for each horse. Relative telomere length (RTL) was highly correlated with overall age. RTL was positively correlated with 3H-thymidine incorporation and total IgG. Expression of pro-inflammatory cytokines was negatively correlated with RTL. These measures were also correlated with age, as expected. However, RTL was not correlated with immunosenescence and inflammaging in the oldest horse.

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.

Characterization of the equine 2'-5' oligoadenylate synthetase 1 (OAS1) and ribonuclease L (RNASEL) innate immunity genes.

BACKGROUND: The mammalian OAS/RNASEL pathway plays an important role in antiviral host defense. A premature stop-codon within the murine Oas1b gene results in the increased susceptibility of mice to a number of flaviviruses, including West Nile virus (WNV). Mutations in either the OAS1 or RNASEL genes may also modulate the outcome of WNV-induced disease or other viral infections in horses. Polymorphisms in the human OAS gene cluster have been previously utilized for case-control analysis of virus-induced disease in humans. No polymorphisms have yet been identified in either the equine OAS1 or RNASEL genes for use in similar case-control studies. RESULTS: Genomic sequence for equine OAS1 was obtained from a contig assembly generated from a shotgun subclone library of CHORI-241 BAC 100I10. Specific amplification of regions of the OAS1 gene from 13 horses of various breeds identified 33 single nucleotide polymorphisms (SNP) and two microsatellites. RNASEL cDNA sequences were determined for 8 mammals and utilized in a phylogenetic analysis. The chromosomal location of the RNASEL gene was assigned by FISH to ECA5p17-p16 using two selected CHORI-241 BAC clones. The horse genomic RNASEL sequence was assembled. Specific amplification of regions of the RNASEL gene from 13 horses identified 31 SNPs. CONCLUSION: In this report, two dinucleotide microsatellites and 64 single nucleotide polymorphisms within the equine OAS1 and RNASEL genes were identified. These polymorphisms are the first to be reported for these genes and will facilitate future case-control studies of horse susceptibility to infectious diseases.

Comparative analysis of vertebrate EIF2AK2 (PKR) genes and assignment of the equine gene to ECA15q24-q25 and the bovine gene to BTA11q12-q15.

The structures of the canine, rabbit, bovine and equine EIF2AK2 genes were determined. Each of these genes has a 5' non-coding exon as well as 15 coding exons. All of the canine, bovine and equine EIF2AK2 introns have consensus donor and acceptor splice sites. In the equine EIF2AK2 gene, a unique single nucleotide polymorphism that encoded a Tyr329Cys substitution was detected. Regulatory elements predicted in the promoter region were conserved in ungulates, primates, rodents, Afrotheria (elephant) and Insectifora (shrew). Western clawed frog and fugu EIF2AK2 gene sequences were detected in the USCS Genome Browser and compared to those of other vertebrate EIF2AK2 genes. A comparison of EIF2AK2 protein domains in vertebrates indicates that the kinase catalytic domains were evolutionarily more conserved than the nucleic acid-binding motifs. Nucleotide substitution rates were uniform among the vertebrate sequences with the exception of the zebrafish and goldfish EIF2AK2 genes, which showed substitution rates about 20% higher than those of other vertebrates. FISH was used to physically assign the horse and cattle genes to chromosome locations, ECA15q24-q25 and BTA11q12-15, respectively. Comparative mapping data confirmed conservation of synteny between ungulates, humans and rodents.

Genomic characterization of MHC class I genes of the horse.

The availability of a contig of bacterial artificial chromosome (BAC) clones spanning the equine major histocompatibility complex (MHC) made possible a detailed analysis of horse MHC class I genes. Prior to this study, only a single horse MHC class I gene had been sequenced at the genomic level. Although many ( approximately 60) MHC class I cDNA sequences had been determined and published, from this information, it was not possible to determine how many class I loci are expressed in horses or to assign individual sequences to allelic series. In this study, 15 MHC class I genes were identified in BAC subclones and fully sequenced. Because the BAC library donor horse had been bred for homozygosity at the MHC, these 15 genomic clones represent distinct MHC class I genes and pseudogenes and not alleles at a smaller number of loci. For five of the genes, cDNA sequences from these loci had previously been identified. Two additional expressed class I genes were discovered, bringing the known total of different equine MHC class I genes (loci) expressed as mRNA to seven. Expression of all seven loci was detected by reverse transcriptase-polymerase chain reaction in adult, fetal, and placental tissues. The remaining eight genes were designated as pseudogenes. This work resulted in moderate expansion of the horse MHC BAC contig length, and the remaining gap was shortened. The information contained in these equine MHC class I sequences will permit comparison of MHC class I genes expressed across different horse MHC haplotypes and between horses and other mammalian species.

The complete map of the Ig heavy chain constant gene region reveals evidence for seven IgG isotypes and for IgD in the horse.

This report contains the first map of the complete Ig H chain constant (IGHC) gene region of the horse (Equus caballus), represented by 34 overlapping clones from a new bacterial artificial chromosome library. The different bacterial artificial chromosome inserts containing IGHC genes were identified and arranged by hybridization using overgo probes specific for individual equine IGHC genes. The analysis of these IGHC clones identified two previously undetected IGHC genes of the horse. The newly found IGHG7 gene, which has a high homology to the equine IGHG4 gene, is located between the IGHG3 and IGHG4 genes. The high degree of conservation shared between the nucleotide sequences of the IGHG7 and IGHG4 genes is unusual for the IGHG genes of the horse and suggests that these two genes duplicated most recently during evolution of the equine IGHG genes. Second, we present the genomic nucleotide sequence of the equine IGHD gene, which is located downstream of the IGHM gene. Both the IGHG7 and IGHD genes were found to be expressed at the mRNA level. The order of the 11 IGHC genes in the IGH-locus of the horse was determined to be 5'-M-D-G1-G2-G3-G7-G4-G6-G5-E-A-3', confirming previous studies using lambda phage clones, with the exception that the IGHG5 gene was found to be the most downstream-located IGHG gene. Fluorescence in situ hybridization was used to localize the IGHC region to Equus caballus (ECA) 24qter, the horse chromosome corresponding to human chromosome 14, where the human IGH locus is found.

Natural killer cell receptors in the horse: evidence for the existence of multiple transcribed LY49 genes.

In rodents, the Ly49 family encodes natural killer (NK) receptors interacting with classical MHC class I molecules, whereas the corresponding receptors in primates are members of the killer cell immunoglobulin-like receptor (KIR) family. Recent evidence indicates that the cattle, domestic cat, dog, and pig have a single LY49 and multiple KIR genes, suggesting that predominant NK receptors in most non-rodent mammals might be KIR. Here, we show that the horse has at least six LY49 genes, five with an immunoreceptor tyrosine-based inhibition motif (ITIM) and one with arginine in the transmembrane region. Interestingly, none of the horse KIR-like cDNA clones isolated by library screening encoded molecules likely to function asNK receptors; four types of clones were KIR-Ig-like transcript (KIR-ILT) hybrids and contained premature stop codons and/or frameshift mutations, and two putative allelic sequences predicting KIR3DL molecules had mutated ITIM. To our knowledge, this is the first report suggesting that non-rodent mammals may use LY49 as NK receptors for classical MHC class I. We also show that horse spleen expresses ILT-like genes with unique domain organizations. Radiation hybrid mapping and fluorescence in situ hybridization localized horse LY49 and KIR/ILT genes to chromosomes 6q13 and 10p12, respectively.

Characterization of the beta2-microglobulin gene of the horse.

A clone containing beta(2)-microglobulin (beta(2)-m), the light chain of the major histocompatibility complex class I cell surface molecule, was isolated from an equine bacterial artificial chromosome library. This clone was used as a template for polymerase chain reaction (PCR) and unidirectional sequencing to elucidate the genomic sequence and intron/exon boundaries. We obtained 7,000 bases of sequence, extending from 1,100 nucleotides (nt) upstream of the coding region start through 1,698 nt downstream of the stop codon. The sequence contained regulatory elements in the region upstream of the coding sequence similar to those of the beta(2)-m gene of other species. The beta(2)-m gene was localized to horse chromosome ECA1q23-q25 by fluorescent in situ hybridization. This was confirmed by synteny mapping on a (horse x mouse) somatic cell hybrid panel. The sequence and intron/exon boundaries determined were used to design PCR primers to amplify and sequence the coding region of the beta(2)-m gene in other equids, including five breeds of domestic horse, one Przewalski's horse, five domestic donkeys and five zebras. A high degree of conservation was found among equids, illustrated by >98% (349/354) identity at the nucleotide level and 95% (113/118) at the amino acid level, because of non-synonymous nucleotide substitutions. The promoter detected in the region upstream of the coding sequence was subcloned and used in chloramphenicol acetyl transferase (CAT) assays to demonstrate the presence of a functional promoter. This study provides tools for the analysis of regulation of not only the horse beta(2)-m gene, but also for any genes dependent upon beta(2)-m for expression.

Cytogenetic localization of 136 genes in the horse: comparative mapping with the human genome.

The aim of this study was to increase the number of type I markers on the horse cytogenetic map and to improve comparison with maps of other species, thus facilitating positional candidate cloning studies. BAC clones from two different sources were FISH mapped: homologous horse BAC clones selected from our newly extended BAC library using consensus primer sequences and heterologous goat BAC clones. We report the localization of 136 genes on the horse cytogenetic map, almost doubling the number of cytogenetically mapped genes with 48 localizations from horse BAC clones and 88 from goat BAC clones. For the first time, genes were mapped to ECA13p, ECA29, and probably ECA30. A total of 284 genes are now FISH mapped on the horse chromosomes. Comparison with the human map defines 113 conserved segments that include new homologous segments not identified by Zoo-FISH on ECA7 and ECA13p.

The first-generation whole-genome radiation hybrid map in the horse identifies conserved segments in human and mouse genomes.

A first-generation radiation hybrid (RH) map of the equine (Equus caballus) genome was assembled using 92 horse x hamster hybrid cell lines and 730 equine markers. The map is the first comprehensive framework map of the horse that (1) incorporates type I as well as type II markers, (2) integrates synteny, cytogenetic, and meiotic maps into a consensus map, and (3) provides the most detailed genome-wide information to date on the organization and comparative status of the equine genome. The 730 loci (258 type I and 472 type II) included in the final map are clustered in 101 RH groups distributed over all equine autosomes and the X chromosome. The overall marker retention frequency in the panel is approximately 21%, and the possibility of adding any new marker to the map is approximately 90%. On average, the mapped markers are distributed every 19 cR (4 Mb) of the equine genome--a significant improvement in resolution over previous maps. With 69 new FISH assignments, a total of 253 cytogenetically mapped loci physically anchor the RH map to various chromosomal segments. Synteny assignments of 39 gene loci complemented the RH mapping of 27 genes. The results added 12 new loci to the horse gene map. Lastly, comparison of the assembly of 447 equine genes (256 linearly ordered RH-mapped and additional 191 FISH-mapped) with the location of draft sequences of their human and mouse orthologs provides the most extensive horse-human and horse-mouse comparative map to date. We expect that the foundation established through this map will significantly facilitate rapid targeted expansion of the horse gene map and consequently, mapping and positional cloning of genes governing traits significant to the equine industry.