Single linkage group per chromosome genetic linkage map for the horse, based on two three-generation, full-sibling, crossbred horse reference families

Exploration of fine-scale recombination rate variation in the domestic horse

Total genetic map length and local recombination landscapes typically vary within and across populations. As a first step to understanding the recombination landscape in the domestic horse, we calculated population recombination rates and identified likely recombination hotspots using approximately 1.8 million SNP genotypes for 485 horses from 32 distinct breeds. The resulting breed-averaged recombination map spans 2.36 Gb and accounts for 2939.07 cM. Recombination hotspots occur once per 23.8 Mb on average and account for ∼9% of the physical map length. Regions with elevated recombination rates in the entire cohort were enriched for genes in pathways involving interaction with the environment: immune system processes (specifically, MHC class I and class II genes), responses to stimuli, and serotonin receptor pathways. We found significant correlations between differences in local recombination rates and population differentiation quantified by F _ST Analysis of breed-specific maps revealed thousands of hotspot regions unique to particular breeds, as well as unique "coldspots," regions where a particular breed showed below-average recombination, whereas all other breeds had evidence of a hotspot. Finally, we identified relative enrichment (P = 5.88 × 10^-27) for the in silico-predicted recognition motif for equine PR/SET domain 9 (PRDM9) in recombination hotspots. These results indicate that selective pressures and PRDM9 function contribute to variation in recombination rates across the domestic horse genome.

Genetic variability of Akhal-Teke horses bred in Italy

Background

The Akhal-Teke horse (AKH) is native of the modern Turkmenistan area. It was introduced in Italy from 1991 to 2000 mainly as an endurance horse. This paper characterizes the genetic variability of the whole Italian AKH horse population and evaluates their inbreeding level by analyzing microsatellite markers and mitochondrial D-Loop sequences.

Methods

Seventeen microsatellite marker loci were genotyped on 95 DNA samples from almost all the AKH horses bred in Italy in the last 20 years. Standard genetic variability measures (H_o, H_e, F_IS) were compared against the same variables published on other eight AKH populations. In addition, 397 bp of mtDNA D-loop region were sequenced on a sub-group of 22 unrelated AKH out of the 95 sampled ones, and on 11 unrelated Arab horses. The haplotypes identified in the Italian population were aligned to sequences of AKH (56), Arab (five), Caspian Pony (13), Przewalskii (two) and Barb (15) horses available in GenBank. The Median Joining Network (MJN), Principal Component Analysis (PCA) and Neighbor-joining (NJ) tree were calculated on the total 126 sequences.

Results

Nucleic markers showed a high degree of polymorphism (H_o = 0.642; H_e = 0.649) and a low inbreeding level (F_IS = 0.016) in Italian horses, compared to other AKH populations (ranged from −0.103 AKH from Estonia to 0.114 AKH from Czech Republic). High variability was also recorded in the D-Loop region. 11 haplotypes were identified with haplotype diversity (hd), nucleotide diversity (π) and average number of nucleotide differences (k) of 0.938, 0.021 and 6.448, respectively. When all the 126 D-Loop sequences were compared, 51 haplotypes were found, and four were here found only in the Italian AKH horses. The 51 haplotypes were conformed to eight recognized mtDNA haplogroups (A, C, F, G, L, M, P and Q) and confirmed by MJN analysis, Italian horses being assigned to five haplogroups (A, C, G, L and M). Using a PCA approach to the same data, the total haplotypes were grouped into two clusters including A+C+M+P and G+F haplogroups, while L and Q haplogroups remained ungrouped. Finally, the NJ algorithm effectively discretizes only the L haplogroup. All the above data univocally indicate good genetic variability and accurate management of the Akhal-Teke population in Italy.

Variation and Evolution of the Meiotic Requirement for Crossing Over in Mammals

The segregation of homologous chromosomes at the first meiotic division is dependent on the presence of at least one well-positioned crossover per chromosome. In some mammalian species, however, the genomic distribution of crossovers is consistent with a more stringent baseline requirement of one crossover per chromosome arm. Given that the meiotic requirement for crossing over defines the minimum frequency of recombination necessary for the production of viable gametes, determining the chromosomal scale of this constraint is essential for defining crossover profiles predisposed to aneuploidy and understanding the parameters that shape patterns of recombination rate evolution across species. Here, I use cytogenetic methods for in situ imaging of crossovers in karyotypically diverse house mice (Mus musculus domesticus) and voles (genus Microtus) to test how chromosome number and configuration constrain the distribution of crossovers in a genome. I show that the global distribution of crossovers in house mice is thresholded by a minimum of one crossover per chromosome arm, whereas the crossover landscape in voles is defined by a more relaxed requirement of one crossover per chromosome. I extend these findings in an evolutionary metaanalysis of published recombination and karyotype data for 112 mammalian species and demonstrate that the physical scale of the genomic crossover distribution has undergone multiple independent shifts from one crossover per chromosome arm to one per chromosome during mammalian evolution. Together, these results indicate that the chromosomal scale constraint on crossover rates is itself a trait that evolves among species, a finding that casts light on an important source of crossover rate variation in mammals.

Developmental progress and current status of the Animal QTLdb

The Animal QTL Database (QTLdb; http://www.animalgenome.org/QTLdb) has undergone dramatic growth in recent years in terms of new data curated, data downloads and new functions and tools. We have focused our development efforts to cope with challenges arising from rapid growth of newly published data and end users' data demands, and to optimize data retrieval and analysis to facilitate users' research. Evidenced by the 27 releases in the past 11 years, the growth of the QTLdb has been phenomenal. Here we report our recent progress which is highlighted by addition of one new species, four new data types, four new user tools, a new API tool set, numerous new functions and capabilities added to the curator tool set, expansion of our data alliance partners and more than 20 other improvements. In this paper we present a summary of our progress to date and an outlook regarding future directions.

A genome-wide scan for tying-up syndrome in Japanese Thoroughbreds

Tying-up syndrome, also known as recurrent exertional rhabdomyolysis in Thoroughbreds, is a common muscle disorder for racehorses. In this study, we performed a multipoint linkage analysis using LOKI based on the Bayesian Markov chain Monte Carlo method using 5 half-sib families (51 affected and 277 nonaffected horses in total), and a genome-wide association study (GWAS) using microsatellites (144 affected and 144 nonaffected horses) to map candidate regions for tying-up syndrome in Japanese Thoroughbreds. The linkage analysis identified one strong L-score (82.45) between the loci UCDEQ411 and COR058 (24.9-27.9 Mb) on ECA12. The GWAS identified two suggestive genomic regions on ECA12 (24.9-27.8 Mb) and ECA20 (29.3-33.5 Mb). Based on both results, the genomic region between UCDEQ411 and TKY499 (24.9-27.8 Mb) on ECA12 was the most significant and was considered as a candidate region for tying-up syndrome in Japanese Thoroughbreds.

Linkage disequilibrium and historical effective population size in the Thoroughbred horse

Many genomic methodologies rely on the presence and extent of linkage disequilibrium (LD) between markers and genetic variants underlying traits of interest, but the extent of LD in the horse has yet to be comprehensively characterized. In this study, we evaluate the extent and decay of LD in a sample of 817 Thoroughbreds. Horses were genotyped for over 50,000 single nucleotide polymorphism (SNP) markers across the genome, with 34,848 autosomal SNPs used in the final analysis. Linkage disequilibrium, as measured by the squared correlation coefficient (r(2)), was found to be relatively high between closely linked markers (>0.6 at 5 kb) and to extend over long distances, with average r(2) maintained above non-syntenic levels for single nucleotide polymorphisms (SNPs) up to 20 Mb apart. Using formulae which relate expected LD to effective population size (N(e)), and assuming a constant actual population size, N(e) was estimated to be 100 in our population. Values of historical N(e), calculated assuming linear population growth, suggested a decrease in N(e) since the distant past, reaching a minimum twenty generations ago, followed by a subsequent increase until the present time. The qualitative trends observed in N(e) can be rationalized by current knowledge of the history of the Thoroughbred breed, and inbreeding statistics obtained from published pedigree analyses are in agreement with observed values of N(e). Given the high LD observed and the small estimated N(e), genomic methodologies such as genomic selection could feasibly be applied to this population using the existing SNP marker set.

Male horse meiosis: metaphase I chromosome configuration and chiasmata distribution

Chromosome configurations and chiasma frequency during the metaphase I stage of spermatogenesis in the male horse are characterized in this work. The genome-wide frequency and distribution of chiasmata was detected as 49.45 ± 2.07 for 14 fertile stallions. All X and Y chromosomes shared a single chiasma at their pseudoautosomal region, while 1-4 chiasmata were observed in autosomal chromosomes. The chiasma frequency and distribution were further studied for 8 different bivalents identified by FISH in 5 fertile stallions. Genetic length was calculated from chiasmata data for the whole genome as well as for these 8 chromosomes. The findings complement the genetic linkage data and provide insight into the genetic basis of spermatogenesis in normal stallions.

Applied equine genetics

Genome sequencing of the domestic horse and subsequent advancements in the field of equine genomics have led to an explosion in the development of tools for mapping traits and diseases and evaluating gene expression. The objective of this review is to discuss the current progress in the field of equine genomics, with specific emphasis on assembly and analysis of the reference sequence and subsequent sequencing of a Quarter Horse mare; the genomic tools currently available to researchers and their implications in genomic investigations in the horse; the genomics of Mendelian and non-Mendelian traits; the genomics of performance traits and considerations regarding genetic testing in the horse. The whole-genome sequencing of a Quarter Horse mare has provided additional variants within the equine genome that extend past single nucleotide polymorphisms to include insertions/deletions and copy number variants. Equine single nucleotide polymorphism arrays have allowed for the investigation of both simple and complex genetic traits while DNA microarrays have provided a tool for examining gene expression across various tissues and with certain disease conditions. Recently, next-generation sequencing has become more affordable and both whole-genome DNA sequencing and transcriptome-wide RNA sequencing are methodologies that are being applied to equine genomic research. Research in the field of equine genomics continues to expand rapidly as the cost of genotyping and sequencing decreases, resulting in a need for quality bioinformatics software and expertise to appropriately handle both the size and complexity of these data.

The utility of low-density genotyping for imputation in the Thoroughbred horse

Background

Despite the dramatic reduction in the cost of high-density genotyping that has occurred over the last decade, it remains one of the limiting factors for obtaining the large datasets required for genomic studies of disease in the horse. In this study, we investigated the potential for low-density genotyping and subsequent imputation to address this problem.

Results

Using the haplotype phasing and imputation program, BEAGLE, it is possible to impute genotypes from low- to high-density (50K) in the Thoroughbred horse with reasonable to high accuracy. Analysis of the sources of variation in imputation accuracy revealed dependence both on the minor allele frequency of the single nucleotide polymorphisms (SNPs) being imputed and on the underlying linkage disequilibrium structure. Whereas equidistant spacing of the SNPs on the low-density panel worked well, optimising SNP selection to increase their minor allele frequency was advantageous, even when the panel was subsequently used in a population of different geographical origin. Replacing base pair position with linkage disequilibrium map distance reduced the variation in imputation accuracy across SNPs. Whereas a 1K SNP panel was generally sufficient to ensure that more than 80% of genotypes were correctly imputed, other studies suggest that a 2K to 3K panel is more efficient to minimize the subsequent loss of accuracy in genomic prediction analyses. The relationship between accuracy and genotyping costs for the different low-density panels, suggests that a 2K SNP panel would represent good value for money.

Conclusions

Low-density genotyping with a 2K SNP panel followed by imputation provides a compromise between cost and accuracy that could promote more widespread genotyping, and hence the use of genomic information in horses. In addition to offering a low cost alternative to high-density genotyping, imputation provides a means to combine datasets from different genotyping platforms, which is becoming necessary since researchers are starting to use the recently developed equine 70K SNP chip. However, more work is needed to evaluate the impact of between-breed differences on imputation accuracy.

Distribution of MLH1 foci in horse male synaptonemal complex

Advances in molecular cytogenetics have provided the opportunity to study events during prophase I of meiosis. Immunofluorescent localization of different meiotic protein components were used to characterize the early stages of the first meiotic division in horse spermatocytes. The frequency and distribution of recombination events during prophase I were investigated using the mutL homolog 1 (MLH1) protein that is known to be associated with these events. The frequency and distribution of MLH1 foci were investigated in pachytene nuclei of 6 fertile stallions, and the average relative synaptonemal complex length was found to be highly correlated with the average number of MLH1 foci. The frequency and distribution of MLH1 foci were found to closely correspond to the frequency and distribution of chiasmata on metaphase I chromosomes, and genetic length, calculated from MLH1 foci data, for the whole genome was 2,505.5 cM.

Marsupial genome sequences: providing insight into evolution and disease

Marsupials (metatherians), with their position in vertebrate phylogeny and their unique biological features, have been studied for many years by a dedicated group of researchers, but it has only been since the sequencing of the first marsupial genome that their value has been more widely recognised. We now have genome sequences for three distantly related marsupial species (the grey short-tailed opossum, the tammar wallaby, and Tasmanian devil), with the promise of many more genomes to be sequenced in the near future, making this a particularly exciting time in marsupial genomics. The emergence of a transmissible cancer, which is obliterating the Tasmanian devil population, has increased the importance of obtaining and analysing marsupial genome sequence for understanding such diseases as well as for conservation efforts. In addition, these genome sequences have facilitated studies aimed at answering questions regarding gene and genome evolution and provided insight into the evolution of epigenetic mechanisms. Here I highlight the major advances in our understanding of evolution and disease, facilitated by marsupial genome projects, and speculate on the future contributions to be made by such sequences.

A genetic linkage map and comparative mapping of the prairie vole (Microtus ochrogaster) genome

BACKGROUND

The prairie vole (Microtus ochrogaster) is an emerging rodent model for investigating the genetics, evolution and molecular mechanisms of social behavior. Though a karyotype for the prairie vole has been reported and low-resolution comparative cytogenetic analyses have been done in this species, other basic genetic resources for this species, such as a genetic linkage map, are lacking.

RESULTS

Here we report the construction of a genome-wide linkage map of the prairie vole. The linkage map consists of 406 markers that are spaced on average every 7 Mb and span an estimated ~90% of the genome. The sex average length of the linkage map is 1707 cM, which, like other Muroid rodent linkage maps, is on the lower end of the length distribution of linkage maps reported to date for placental mammals. Linkage groups were assigned to 19 out of the 26 prairie vole autosomes as well as the X chromosome. Comparative analyses of the prairie vole linkage map based on the location of 387 Type I markers identified 61 large blocks of synteny with the mouse genome. In addition, the results of the comparative analyses revealed a potential elevated rate of inversions in the prairie vole lineage compared to the laboratory mouse and rat.

CONCLUSIONS

A genetic linkage map of the prairie vole has been constructed and represents the fourth genome-wide high-resolution linkage map reported for Muroid rodents and the first for a member of the Arvicolinae sub-family. This resource will advance studies designed to dissect the genetic basis of a variety of social behaviors and other traits in the prairie vole as well as our understanding of genome evolution in the genus Microtus.

Gametic phase disequilibrium between the syntenic multiallelic HTG4 and HMS3 markers widely used for parentage testing in Thoroughbred horses

Validation of parentage and horse breed registries through DNA typing relies on estimates of random match probabilities with DNA profiles generated from multiple polymorphic loci. Of the twenty-seven microsatellite loci recommended by the International Society for Animal Genetics for parentage testing in Thoroughbred horses, eleven are located on five chromosomes. An important aspect in determining combined exclusion probabilities is the ascertainment of the genetic linkage status of syntenic markers, which may affect reliable use of the product rule in estimating random match probabilities. In principle, linked markers can be in gametic phase disequilibrium (GD). We aimed at determining the extent, by frequency and strength, of GD between the HTG4 and HMS3 multiallelic loci, syntenic on chromosome 9. We typed the qualified offspring (n (1) = 27; n (2) = 14) of two Quarter Bred stallions (registered by the Brazilian Association of Quarter Horse Breeders) and 121 unrelated horses from the same breed. In the 41 informative meioses analyzed, the frequency of recombination between the HTG4 and HMS3 loci was 0.27. Consistent with genetic map distances, this recombination rate does not fit to the theoretical distribution for independently segregated markers. We estimated sign-based D' coefficients as a measure of GD, and showed that the HTG4 and HMS3 loci are in significant, yet partial and weak, disequilibrium, with two allele pairs involved (HTG4 M/HMS3 P, D'(+) = 0.6274; and HTG4 K/HMS3 P, D'(-) = -0.6096). These results warn against the inadequate inclusion of genetically linked markers in the calculation of combined power of discrimination for Thoroughbred parentage validation.

[Progress on horse genome project]

There is unique genetic information belonging to various kinds of living beings. Understanding of the formation process of organisms and a variety of vital movement is associated with the achievements of genome study. As horse has a notable health condition and great record of the genealogy in the world, thus it becomes a valuable model animal for studying life science. Despite of a late start, the map of the horse genome has undergone unprecedented expansion during the last few years. The current progresses of the horse genome, including genetic map, physical map, comparative genomic map, and functional genomics, were reviewed in this paper. The maps are currently used worldwide to discover genes associated with various traits of significance in horse including general health, disease resistance, reproduction, fertility, athletic performance, phenotypic characteristics like coat color, etc. The results are believed to provide new ideas and approaches for prevention, diagnostics, and therapeutic for horses, and also better foundation of breed selection and equine genetic breeding.

A whole-genome scan for recurrent airway obstruction in Warmblood sport horses indicates two positional candidate regions

Recurrent airway obstruction (RAO), or heaves, is a naturally occurring asthma-like disease that is related to sensitisation and exposure to mouldy hay and has a familial basis with a complex mode of inheritance. A genome-wide scanning approach using two half-sibling families was taken in order to locate the chromosome regions that contribute to the inherited component of this condition in these families. Initially, a panel of 250 microsatellite markers, which were chosen as a well-spaced, polymorphic selection covering the 31 equine autosomes, was used to genotype the two half-sibling families, which comprised in total 239 Warmblood horses. Subsequently, supplementary markers were added for a total of 315 genotyped markers. Each half-sibling family is focused around a severely RAO-affected stallion, and the phenotype of each individual was assessed for RAO and related signs, namely, breathing effort at rest, breathing effort at work, coughing, and nasal discharge, using an owner-based questionnaire. Analysis using a regression method for half-sibling family structures was performed using RAO and each of the composite clinical signs separately; two chromosome regions (on ECA13 and ECA15) showed a genome-wide significant association with RAO at P < 0.05. An additional 11 chromosome regions showed a more modest association. This is the first publication that describes the mapping of genetic loci involved in RAO. Several candidate genes are located in these regions, a number of which are interleukins. These are important signalling molecules that are intricately involved in the control of the immune response and are therefore good positional candidates.

A genome scan for positive selection in thoroughbred horses

Thoroughbred horses have been selected for exceptional racing performance resulting in system-wide structural and functional adaptations contributing to elite athletic phenotypes. Because selection has been recent and intense in a closed population that stems from a small number of founder animals Thoroughbreds represent a unique population within which to identify genomic contributions to exercise-related traits. Employing a population genetics-based hitchhiking mapping approach we performed a genome scan using 394 autosomal and X chromosome microsatellite loci and identified positively selected loci in the extreme tail-ends of the empirical distributions for (1) deviations from expected heterozygosity (Ewens-Watterson test) in Thoroughbred (n = 112) and (2) global differentiation among four geographically diverse horse populations (F(ST)). We found positively selected genomic regions in Thoroughbred enriched for phosphoinositide-mediated signalling (3.2-fold enrichment; P<0.01), insulin receptor signalling (5.0-fold enrichment; P<0.01) and lipid transport (2.2-fold enrichment; P<0.05) genes. We found a significant overrepresentation of sarcoglycan complex (11.1-fold enrichment; P<0.05) and focal adhesion pathway (1.9-fold enrichment; P<0.01) genes highlighting the role for muscle strength and integrity in the Thoroughbred athletic phenotype. We report for the first time candidate athletic-performance genes within regions targeted by selection in Thoroughbred horses that are principally responsible for fatty acid oxidation, increased insulin sensitivity and muscle strength: ACSS1 (acyl-CoA synthetase short-chain family member 1), ACTA1 (actin, alpha 1, skeletal muscle), ACTN2 (actinin, alpha 2), ADHFE1 (alcohol dehydrogenase, iron containing, 1), MTFR1 (mitochondrial fission regulator 1), PDK4 (pyruvate dehydrogenase kinase, isozyme 4) and TNC (tenascin C). Understanding the genetic basis for exercise adaptation will be crucial for the identification of genes within the complex molecular networks underlying obesity and its consequential pathologies, such as type 2 diabetes. Therefore, we propose Thoroughbred as a novel in vivo large animal model for understanding molecular protection against metabolic disease.

An autosomal genetic linkage map of the domestic cat, Felis silvestris catus

We report on the completion of an autosomal genetic linkage (GL) map of the domestic cat (Felis silvestris catus). Unlike two previous linkage maps of the cat constructed with a hybrid pedigree between the domestic cat and the Asian leopard cat, this map was generated entirely with domestic cats, using a large multi-generational pedigree (n=256) maintained by the Nestlé Purina PetCare Company. Four hundred eighty-three simple tandem repeat (STR) loci have been assigned to linkage groups on the cat's 18 autosomes. A single linkage group spans each autosome. The length of the cat map, estimated at 4370 cM, is long relative to most reported mammalian maps. A high degree of concordance in marker order was observed between the third-generation map and the 1.5 Mb-resolution radiation hybrid (RH) map of the cat. Using the cat 1.9x whole-genome sequence, we identified map coordinates for 85% of the loci in the cat assembly, with high concordance observed in marker order between the linkage map and the cat sequence assembly. The present version represents a marked improvement over previous cat linkage maps as it (i) nearly doubles the number of markers that were present in the second-generation linkage map in the cat, (ii) provides a linkage map generated in a domestic cat pedigree which will more accurately reflect recombination distances than previous maps generated in a hybrid pedigree, and (iii) provides single linkage groups spanning each autosome. Marker order was largely consistent between this and the previous maps, though the use of a hybrid pedigree in the earlier versions appears to have contributed to some suppression of recombination. The improved linkage map will provide an added resource for the mapping of phenotypic variation in the domestic cat and the use of this species as a model system for biological research.

The candidate gene XIRP2 at a quantitative gene locus on equine chromosome 18 associated with osteochondrosis in fetlock and hock joints of South German Coldblood horses

A whole-genome scan for radiological signs of osteochondrosis (OC) and osteochondrosis dissecans (OCD) in South German Coldblood (SGC) horses using 250 microsatellite markers identified a genome-wide significant quantitative trait locus (QTL) for fetlock OCD and a chromosome-wide QTL for hock OC on Equus caballus chromosome (ECA) 18 at a relative position of 45.9-78.2 cM. The aim of this study was to analyze associations of single-nucleotide polymorphisms (SNPs) in candidate genes for OC in this QTL region using 96 SGC horses. The OC-QTL on ECA18 could be confirmed and narrowed down to an interval of 13 Mb between GALNT13 and Xin actin-binding repeat containing 2 (XIRP2). SNPs in the XIRP2 gene were significantly associated with fetlock OC, fetlock OCD, and hock OC. The significant associations of SNPs in XIRP2 could be confirmed in linear animal models controlling for systematic environmental and residual quantitative genetic effects. The significant additive genetic effects of the intronic SNPs (AJ885515:g.159A>G, AJ885515:g.445T>C) in XIRP2 were 0.15 (P = 0.01) for fetlock OC, 0.27 (P = 0.01) for fetlock OCD, and 0.15-0.16 (P = 0.01-0.02) for hock OC. Homozygous (A/A or T/T) and heterozygous horses were at a 1.3- to 2.4-fold higher risk for fetlock and hock OC. These results suggest that dominant variants of XIRP2 may be involved in pathogenesis of equine OC.

A domestic cat X chromosome linkage map and the sex-linked orange locus: mapping of orange, multiple origins and epistasis over nonagouti

A comprehensive genetic linkage map of the domestic cat X chromosome was generated with the goal of localizing the genomic position of the classic X-linked orange (O) locus. Microsatellite markers with an average spacing of 3 Mb were selected from sequence traces of the cat 1.9x whole genome sequence (WGS), including the pseudoautosomal region 1 (PAR1). Extreme variation in recombination rates (centimorgans per megabase) was observed along the X chromosome, ranging from a virtual absence of recombination events in a region estimated to be >30 Mb to recombination frequencies of 15.7 cM/Mb in a segment estimated to be <0.3 Mb. This detailed linkage map was applied to position the X-linked orange gene, placing this locus on the q arm of the X chromosome, as opposed to a previously reported location on the p arm. Fine mapping placed the locus between markers at positions 106 and 116.8 Mb in the current 1.9x-coverage sequence assembly of the cat genome. Haplotype analysis revealed potential recombination events that could reduce the size of the candidate region to 3.5 Mb and suggested multiple origins for the orange phenotype in the domestic cat. Furthermore, epistasis of orange over nonagouti was demonstrated at the genetic level.

Polysaccharide storage myopathy phenotype in quarter horse-related breeds is modified by the presence of an RYR1 mutation

In this study we examined a family of Quarter Horses with Polysaccharide Storage Myopathy (PSSM) with a dominant mutation in the skeletal muscle glycogen synthase (GYS1) gene. A subset of horses within this family had a more severe and occasionally fatal PSSM phenotype. The purpose of this study was to identify a modifying gene(s) for the severe clinical phenotype. A genetic association analysis was used to identify RYR1 as a candidate modifying gene. A rare, known equine RYR1 mutation, associated with malignant hyperthermia (MH), was found to segregate in this GYS1 PSSM family. Retrospective analysis of patient records (n=179) demonstrated that horses with both the GYS1 and RYR1 mutations had a more severe clinical phenotype than horses with the GYS1 mutation alone. A treadmill trial (n=8) showed that serum creatine kinase activity was higher and exercise intolerance greater in horses with both mutations compared to the GYS1 mutation alone.

A 4,103 marker integrated physical and comparative map of the horse genome

A comprehensive second-generation whole genome radiation hybrid (RH II), cytogenetic and comparative map of the horse genome (2n = 64) has been developed using the 5000rad horse x hamster radiation hybrid panel and fluorescence in situ hybridization (FISH). The map contains 4,103 markers (3,816 RH; 1,144 FISH) assigned to all 31 pairs of autosomes and the X chromosome. The RH maps of individual chromosomes are anchored and oriented using 857 cytogenetic markers. The overall resolution of the map is one marker per 775 kilobase pairs (kb), which represents a more than five-fold improvement over the first-generation map. The RH II incorporates 920 markers shared jointly with the two recently reported meiotic maps. Consequently the two maps were aligned with the RH II maps of individual autosomes and the X chromosome. Additionally, a comparative map of the horse genome was generated by connecting 1,904 loci on the horse map with genome sequences available for eight diverse vertebrates to highlight regions of evolutionarily conserved syntenies, linkages, and chromosomal breakpoints. The integrated map thus obtained presents the most comprehensive information on the physical and comparative organization of the equine genome and will assist future assemblies of whole genome BAC fingerprint maps and the genome sequence. It will also serve as a tool to identify genes governing health, disease and performance traits in horses and assist us in understanding the evolution of the equine genome in relation to other species.

The horse genome derby: racing from map to whole genome sequence

The map of the horse genome has undergone unprecedented expansion during the past six years. Beginning from a modest collection of approximately 300 mapped markers scattered on the 31 pairs of autosomes and the X chromosome in 2001, today the horse genome is among the best-mapped in domestic animals. Presently, high-resolution linearly ordered gene maps are available for all autosomes as well as the X and the Y chromosome. The approximately 4350 mapped markers distributed over the approximately 2.68 Gbp long equine genome provide on average 1 marker every 620 kb. Among the most remarkable developments in equine genome analysis is the availability of the assembled sequence (EquCab2) of the female horse genome and the generation approximately 1.5 million single nucleotide polymorphisms (SNPs) from diverse breeds. This has triggered the creation of new tools and resources like the 60K SNP-chip and whole genome expression microarrays that hold promise to study the equine genome and transcriptome in ways not previously envisaged. As a result of these developments it is anticipated that, during coming years, the genetics underlying important monogenic traits will be analyzed with improved accuracy and speed. Of larger interest will be the prospects of dissecting the genetic component of various complex/multigenic traits that are of vital significance for equine health and welfare. The number of investigations recently initiated to study a multitude of such traits hold promise for improved diagnostics, prevention and therapeutic approaches for horses.

A BAC contig map over the proximal approximately 3.3 Mb region of horse chromosome 21

A total of 207 BAC clones containing 155 loci were isolated and arranged into a map of linearly ordered overlapping clones over the proximal part of horse chromosome 21 (ECA21), which corresponds to the proximal half of the short arm of human chromosome 19 (HSA19p) and part of HSA5. The clones form two contigs - each corresponding to the respective human chromosomes - that are estimated to be separated by a gap of approximately 200 kb. Of the 155 markers present in the two contigs, 141 (33 genes and 108 STS) were generated and mapped in this study. The BACs provide a 4-5x coverage of the region and span an estimated length of approximately 3.3 Mb. The region presently contains one mapped marker per 22 kb on average, which represents a major improvement over the previous resolution of one marker per 380 kb obtained through the generation of a dense RH map for this segment. Dual color fluorescence in situ hybridization on metaphase and interphase chromosomes verified the relative order of some of the BACs and helped to orient them accurately in the contigs. Despite having similar gene order and content, the equine region covered by the contigs appears to be distinctly smaller than the corresponding region in human (3.3 Mb vs. 5.5-6 Mb) because the latter harbors a host of repetitive elements and gene families unique to humans/primates. Considering limited representation of the region in the latest version of the horse whole genome sequence EquCab2, the dense map developed in this study will prove useful for the assembly and annotation of the sequence data on ECA21 and will be instrumental in rapid search and isolation of candidate genes for traits mapped to this region.

Livestock genomics: bridging the gap between mice and men

Dissecting the genetic control of variation in complex traits, such as disease resistance and agricultural-product quality, remains very challenging. Farm animals are now well placed to bridge the gap between human biology and traditional model species. Livestock species share with model species the benefits of controlled breeding, and their biology is often much closer to that of humans. Genetic research in model species focuses on differences between homogenous lines, whereas genetic research in humans focuses on genetic variation within populations. Livestock genetics has the strengths of both human and model-species genetics because researchers can exploit both the abundant genetic variation between divergent breeds and the variation that is segregating within breeds. Therefore, livestock genomics fills the void where the genetics of model species proves intractable or where model species are not a good proxy for human biology.

Genome-wide search for markers associated with osteochondrosis in Hanoverian warmblood horses

A genome-wide scan was performed to detect quantitative trait loci (QTLs) for osteochondrosis (OC) and osteochondrosis dissecans (OCD) in horses. The marker set comprised 260 microsatellites. We collected data from 211 Hanoverian warmblood horses consisting of 14 paternal half-sib families. Traits used were OC (fetlock and/or hock joints affected), OCD (fetlock and/or hock joints affected), fetlock OC, fetlock OCD, hock OC, and hock OCD. The first genome scan included 172 microsatellite markers. In a second step 88 additional markers were chosen to refine putative QTLs found in the first scan. Genome-wide significant QTLs were located on equine chromosomes 2, 4, 5, and 16. QTLs for fetlock OC and hock OC partly overlapped on the same chromosomes, indicating that these traits may be genetically related. QTLs reached the chromosome-wide significance level on eight different equine chromosomes: 2, 3, 4, 5, 15, 16, 19, and 21. This whole-genome scan was a first step toward the identification of candidate genome regions harboring genes responsible for equine OC. Further investigations are necessary to refine the map positions of the QTLs already identified for OC.

Mapping quantitative trait loci for osteochondrosis in fetlock and hock joints and palmar/plantar osseus fragments in fetlock joints of South German Coldblood horses

The aim of this study was to identify quantitative trait loci (QTL) for osteochondrosis (OC) and palmar/plantar osseous fragments (POF) in fetlock joints in a whole-genome scan of 219 South German Coldblood horses. Symptoms of OC and POF were checked by radiography in 117 South German Coldblood horses at a mean age of 17 months. The radiographic examination comprised the fetlock and hock joints of all limbs. The genome scan included 157 polymorphic microsatellite markers. All microsatellite markers were equally spaced over the 31 autosomes and the X chromosome, with an average distance of 17.7 cM and a mean polymorphism information content (PIC) of 63%. Sixteen chromosomes harbouring putative QTL regions were further investigated by genotyping the animals with 93 additional markers. QTL that had chromosome-wide significance by non-parametric Z-means and LOD scores were found on 10 chromosomes. This included seven QTL for fetlock OC and one QTL on ECA18 associated with hock OC and fetlock OC. Significant QTL for POF in fetlock joints were located on equine chromosomes 1, 4, 8, 12 and 18. This genome scan is an important step towards the identification of genes responsible for OC in horses.

Genome-wide search for microsatellite markers associated with radiologic alterations in the navicular bone of Hanoverian warmblood horses

The aim of this study was to identify quantitative trait loci (QTLs) for pathologic changes in the navicular bone in Hanoverian warmblood horses. Seventeen paternal half-sib groups comprising 192 individuals were analyzed in a whole-genome scan. These families included 144 progeny and grandchildren, which were randomly chosen from the Hanoverian warmblood. Three different traits were considered: deformed canales sesamoidales and radiographic changes in the contour and in the structure of the navicular bone. The genome scan included in total 214 highly polymorphic microsatellite markers. The putatively linked genomic regions on equine chromosomes (ECA) 2, 3, 10, and 15 were refined using 53 additional microsatellites. Chromosome-wide significant QTLs were located on five different equine chromosomes (ECA2, 3, 4, 10, and 26). Genome-wide significant QTLs were on ECA2 at 48 cM and on ECA10 from 45.5 to 49.8 cM. This study was a first step to get more insight into the molecular genetic determination of radiologic changes in the equine navicular bone.

Homozygosity mapping approach identifies a missense mutation in equine cyclophilin B (PPIB) associated with HERDA in the American Quarter Horse

Hereditary equine regional dermal asthenia (HERDA), a degenerative skin disease that affects the Quarter Horse breed, was localized to ECA1 by homozygosity mapping. Comparative genomics allowed the development of equine gene-specific markers which were used with a set of affected horses to detect a homozygous, identical-by-descent block spanning approximately 2.5 Mb, suggesting a recent origin for the HERDA mutation. We report a mutation in cyclophilin B (PPIB) as a novel, causal candidate gene for HERDA. A c.115G>A missense mutation in PPIB alters a glycine residue that has been conserved across vertebrates. The mutation was homozygous in 64 affected horses and segregates concordant with inbreeding loops apparent in the genealogy of 11 affected horses. Screening of control Quarter Horses indicates a 3.5% carrier frequency. The development of a test that can detect affected horses prior to development of clinical signs and carriers of HERDA will allow Quarter Horse breeders to eliminate this debilitating disease.

Improved resolution of the comparative horse–human map: Investigating markers with in silico and linkage mapping approaches

Genetic maps are extremely important tools for tracing the genes that govern economically significant traits, and microsatellites are a significant component of these. In this study, we isolated 2346 novel horse microsatellites as resources for the construction of high-density horse genetic maps. Of these 2346 markers, 339 (14.5%) horse sequences showed sequence homology to DNA sequences in the human genome, demonstrating that microsatellites as type II markers are valuable resources for developing linkage maps and that they have a potential equal to that of type I markers for developing comparative maps. Of the 339 markers, 206 (60.8%) were assigned to horse chromosomes using the Animal Health Trust (AHT) full-sib reference family, and 195 (94.6%) of these localized to the expected syntenic locations on the human genome. These results confirmed the high level of accuracy of in silico mapping. Thus, the 339 markers that exhibited homology to the human genome increased the density of markers on the horse-human comparative map. The resulting comparative map will facilitate the use of horse microsatellites as genetic markers for the identification of quantitative trait loci (QTL) that have been mapped on the human genome. In addition, although the in silico and linkage mapping data did not agree for the other 11 (5.4%) of the assigned 206 markers, these may represent new putative regions of horse-human synteny.

Whole-genome linkage disequilibrium screening for complex traits in horses

The identification of candidate genes for significant traits is crucial. In this study, we developed and tested effective and systematic methods based on linkage disequilibrium (LD) for the identification of candidate regions for genes with Mendelian inheritance and those associated with complex traits. Our approach entailed the combination of primary screening using pooled DNA samples based on DeltaTAC, secondary screening using an individual typing method and tertiary screening using a permutation test based on the differences in the haplotype frequency between two neighbouring microsatellites. This series of methods was evaluated using horse coat colour traits (chestnut/non-chestnut) as a simple Mendelian inheritance model. In addition, the methods were evaluated using a complex trait model constructed by mixing samples from chestnut and non-chestnut horses. Using both models, the methods could detect the expected regions for the horse coat colour trait. The results revealed that LD extends up to several centimorgans in horses, indicating that whole-genome LD screening in horses could be performed systematically and efficiently by combining the above-mentioned methods. Since genetic maps based on microsatellites have been constructed for many other species, the approaches present here could have wide applicability.

A high-resolution comparative radiation hybrid map of equine chromosome 4q12-q22

In this study, we present a comprehensive 5000-rad radiation hybrid map of a 40-cM region on equine chromosome 4 (ECA4) that contains quantitative trait loci for equine osteochondrosis. We mapped 29 gene-associated sequence tagged site markers using primers designed from equine expressed sequence tags or BAC clones in the ECA4q12-q22 region. Three blocks of conserved synteny, showing two chromosomal breakpoints, were identified in the segment of ECA4q12-q22. Markers from other segments of HSA7q mapped to ECA13p and ECA4p, and a region of HSA7p was homologous to ECA13p. Therefore, we have improved the resolution of the human-equine comparative map, which allows the identification of candidate genes underlying traits of interest.

Molecular characterization of the equine collagen, type IX, alpha 2 (COL9A2) gene on horse chromosome 2p16-->p15

The mammalian collagen, type IX, alpha 2 gene (COL9A2) encodes the alpha-2 chain of type IX collagen and is located on horse chromosome 2p16-->p14 harbouring a quantitative trait locus for osteochondrosis. We isolated a bacterial artificial chromosome (BAC) clone containing the equine COL9A2 gene and determined the complete genomic sequence of this gene. Cloning and characterization of equine COL9A2 revealed that the equine gene consists of 32 exons spanning approximately 15 kb. The COL9A2 transcript encodes a single protein of 688 amino acids. Thirty two single nucleotide polymorphisms (SNPs) equally distributed in the gene were detected in a mutation scan of eight unrelated Hanoverian warmblood stallions, including one SNP that affects the amino acid sequence of COL9A2. Comparative analyses between horse, human, mouse and rat indicate that the chromosomal location of equine COL9A2 is in agreement with known chromosomal synteny relationships. The comparison of the gene structure and transcript revealed a high degree of conservation towards the other mammalian COL9A2 genes. We chose three informative SNPs for association and linkage disequilibrium tests in three to five paternal half-sib families of Hanoverian warmblood horses consisting of 44 to 75 genotyped animals. The test statistics did not reach the significance threshold of 5% and so we could not show an association of COL9A2 with equine osteochondrosis.

High-resolution gene maps of horse chromosomes 14 and 21: additional insights into evolution and rearrangements of HSA5 homologs in mammals

High-resolution physically ordered gene maps for equine homologs of human chromosome 5 (HSA5), viz., horse chromosomes 14 and 21 (ECA14 and ECA21), were generated by adding 179 new loci (131 gene-specific and 48 microsatellites) to the existing maps of the two chromosomes. The loci were mapped primarily by genotyping on a 5000-rad horse x hamster radiation hybrid panel, of which 28 were mapped by fluorescence in situ hybridization. The approximately fivefold increase in the number of mapped markers on the two chromosomes improves the average resolution of the map to 1 marker/0.9 Mb. The improved resolution is vital for rapid chromosomal localization of traits of interest on these chromosomes and for facilitating candidate gene searches. The comparative gene mapping data on ECA14 and ECA21 finely align the chromosomes to sequence/gene maps of a range of evolutionarily distantly related species. It also demonstrates that compared to ECA14, the ECA21 segment corresponding to HSA5 is a more conserved region because of preserved gene order in a larger number of and more diverse species. Further, comparison of ECA14 and the distal three-quarters region of ECA21 with corresponding chromosomal segments in 50 species belonging to 11 mammalian orders provides a broad overview of the evolution of these segments in individual orders from the putative ancestral chromosomal configuration. Of particular interest is the identification and precise demarcation of equid/Perissodactyl-specific features that for the first time clearly distinguish the origins of ECA14 and ECA21 from similar-looking status in the Cetartiodactyls.

Exclusion of linkage of theRYR1, CACNA1S, andATP2A1genes to recurrent exertional rhabdomyolysis in Thoroughbreds

OBJECTIVE

To determine whether there was genetic linkage between the recurrent exertional rhabdomyolysis (RER) trait in Thoroughbred horse pedigrees and DNA markers in genes (the sarcoplasmic reticulum calcium release channel [RYR1] gene, the sarcoplasmic reticulum calcium ATPase [ATP2A1] gene, and the transverse tubule dihydropyridine receptor-voltage sensor [CACNA1S] gene) that are important in myoplasmic calcium regulation.

ANIMALS

34 horses in the University of Minnesota RER resource herd and 62 Thoroughbreds from 3 families of Thoroughbreds outside of the university in which RER-affected status was assigned after 2 or more episodes of ER had been observed.

PROCEDURES

Microsatellite DNA markers from the RYR1, ATP2A1, and CACNA1S gene loci on equine chromosomes 10, 13, and 30 were identified. Genotypes were obtained for all horses in the 4 families affected by RER, and data were used to test for linkage of these 3 loci to the RER phenotype.

RESULTS

Analysis of the RYR1, CACNA1S, and ATP2A1 microsatellites excluded a link between those markers and the RER trait.

CONCLUSIONS AND CLINICAL RELEVANCE

It is likely that the heritable alterations in muscle contractility that are characteristic of RER are caused by a gene that is not yet known to cause related muscle disease in other species.