Genetic diversity in the modern horse illustrated from genome-wide SNP data.

Horses were domesticated from the Eurasian steppes 5,000-6,000 years ago. Since then, the use of horses for transportation, warfare, and agriculture, as well as selection for desired traits and fitness, has resulted in diverse populations distributed across the world, many of which have become or are in the process of becoming formally organized into closed, breeding populations (breeds). This report describes the use of a genome-wide set of autosomal SNPs and 814 horses from 36 breeds to provide the first detailed description of equine breed diversity. F(ST) calculations, parsimony, and distance analysis demonstrated relationships among the breeds that largely reflect geographic origins and known breed histories. Low levels of population divergence were observed between breeds that are relatively early on in the process of breed development, and between those with high levels of within-breed diversity, whether due to large population size, ongoing outcrossing, or large within-breed phenotypic diversity. Populations with low within-breed diversity included those which have experienced population bottlenecks, have been under intense selective pressure, or are closed populations with long breed histories. These results provide new insights into the relationships among and the diversity within breeds of horses. In addition these results will facilitate future genome-wide association studies and investigations into genomic targets of selection.

Drawing together the knowledge of forestry and pastoralism experts in the construction of a technical support tool for silvopastoralism.

A growing body of literature documents emphasizes the value of the hybridisation of different types of knowledge for environmental management. This is particularly true for silvopastoralism: although silvopastoralism appears to be a relevant way to manage large forests areas in the northern Mediterranean region, knowledge pertaining to silvopastoralism is distributed between foresters and pastoralists, with a lack of shared knowledge making joint action difficult. The consensus about the need of knowledge hybridisation, however, is not enough to make it happen and a whole field of research investigates the methods which can be used to create shared knowledge. This paper studies one of these methods - the participation of experts in modelling - used in a French research development program to design a technical support tool for silvopastoralism. After having described the modelling process conducted by forestry and pastoralists experts, it presents the conceptual silvopastoralism model they produced and analyses the shared knowledge - especially intermediary concepts and expert typologies - it contains. The paper finally outlines the conditions needed for the participation of experts in modelling to create shared knowledge, its limitations in generating quantitative expertise and discusses how the knowledge created can be shared with a wider audience.

Genome-wide analysis reveals selection for important traits in domestic horse breeds.

Intense selective pressures applied over short evolutionary time have resulted in homogeneity within, but substantial variation among, horse breeds. Utilizing this population structure, 744 individuals from 33 breeds, and a 54,000 SNP genotyping array, breed-specific targets of selection were identified using an F(ST)-based statistic calculated in 500-kb windows across the genome. A 5.5-Mb region of ECA18, in which the myostatin (MSTN) gene was centered, contained the highest signature of selection in both the Paint and Quarter Horse. Gene sequencing and histological analysis of gluteal muscle biopsies showed a promoter variant and intronic SNP of MSTN were each significantly associated with higher Type 2B and lower Type 1 muscle fiber proportions in the Quarter Horse, demonstrating a functional consequence of selection at this locus. Signatures of selection on ECA23 in all gaited breeds in the sample led to the identification of a shared, 186-kb haplotype including two doublesex related mab transcription factor genes (DMRT2 and 3). The recent identification of a DMRT3 mutation within this haplotype, which appears necessary for the ability to perform alternative gaits, provides further evidence for selection at this locus. Finally, putative loci for the determination of size were identified in the draft breeds and the Miniature horse on ECA11, as well as when signatures of selection surrounding candidate genes at other loci were examined. This work provides further evidence of the importance of MSTN in racing breeds, provides strong evidence for selection upon gait and size, and illustrates the potential for population-based techniques to find genomic regions driving important phenotypes in the modern horse.

A method for proteomic analysis of equine subchondral bone and epiphyseal cartilage.

Proteomic analyses of cartilage and, to a lesser extent, of bone have long been impaired because of technical challenges related to their structure and biochemical properties. We have developed a unified method based on phenol extraction, 2DE, silver staining, and subsequent LC-MS/MS. This method proved to be efficient to characterize the proteome of equine cartilage and bone samples collected in vivo. Since proteins from several cellular compartments could be recovered, our procedure is mainly suitable for in situ molecular physiology studies focused on the cellular content of chondrocytes, osteoblasts, and osteoclasts as well as that of the extracellular matrix, with the exception of proteoglycans. Our method alleviates some drawbacks of cell culture that can mask physiological differences, as well as reduced reproducibility due to fractionation. Proteomic comparative studies between cartilage and bone samples from healthy and affected animals were thus achieved successfully. This achievement will contribute to increasing knowledge on the molecular mechanisms underlying the physiopathology of numerous osteoarticular diseases in horses and in humans.

A high density SNP array for the domestic horse and extant Perissodactyla: utility for association mapping, genetic diversity, and phylogeny studies.

An equine SNP genotyping array was developed and evaluated on a panel of samples representing 14 domestic horse breeds and 18 evolutionarily related species. More than 54,000 polymorphic SNPs provided an average inter-SNP spacing of ∼43 kb. The mean minor allele frequency across domestic horse breeds was 0.23, and the number of polymorphic SNPs within breeds ranged from 43,287 to 52,085. Genome-wide linkage disequilibrium (LD) in most breeds declined rapidly over the first 50-100 kb and reached background levels within 1-2 Mb. The extent of LD and the level of inbreeding were highest in the Thoroughbred and lowest in the Mongolian and Quarter Horse. Multidimensional scaling (MDS) analyses demonstrated the tight grouping of individuals within most breeds, close proximity of related breeds, and less tight grouping in admixed breeds. The close relationship between the Przewalski's Horse and the domestic horse was demonstrated by pair-wise genetic distance and MDS. Genotyping of other Perissodactyla (zebras, asses, tapirs, and rhinoceros) was variably successful, with call rates and the number of polymorphic loci varying across taxa. Parsimony analysis placed the modern horse as sister taxa to Equus przewalski. The utility of the SNP array in genome-wide association was confirmed by mapping the known recessive chestnut coat color locus (MC1R) and defining a conserved haplotype of ∼750 kb across all breeds. These results demonstrate the high quality of this SNP genotyping resource, its usefulness in diverse genome analyses of the horse, and potential use in related species.

LIPH expression in skin and hair follicles of normal coat and Rex rabbits.

Natural mutations in the LIPH gene were shown to be responsible for hair growth defects in humans and for the rex short hair phenotype in rabbits. In this species, we identified a single nucleotide deletion in LIPH (1362delA) introducing a stop codon in the C-terminal region of the protein. We investigated the expression of LIPH between normal coat and rex rabbits during critical fetal stages of hair follicle genesis, in adults and during hair follicle cycles. Transcripts were three times less expressed in both fetal and adult stages of the rex rabbits than in normal rabbits. In addition, the hair growth cycle phases affected the regulation of the transcription level in the normal and mutant phenotypes differently. LIPH mRNA and protein levels were higher in the outer root sheath (ORS) than in the inner root sheath (IRS), with a very weak signal in the IRS of rex rabbits. In vitro transfection shows that the mutant protein has a reduced lipase activity compared to the wild type form. Our results contribute to the characterization of the LIPH mode of action and confirm the crucial role of LIPH in hair production.

Chromosomal assignment of six genes (EIF4G3, HSP90, RBBP6, IL8, TERT, and TERC) in four species of the genus Equus.

We mapped six genes (EIF4G3, HSP90, RBBP6, IL8, TERT, and TERC) on the chromosomes of Equus caballus, Equus asinus, Equus grevyi, and Equus burchelli by fluorescence in situ hybridization. Our results add six type I markers to the cytogenetic map of these species and provide new information on the comparative genomics of the genus Equus.

A deletion in exon 9 of the LIPH gene is responsible for the rex hair coat phenotype in rabbits (Oryctolagus cuniculus).

The fur of common rabbits is constituted of 3 types of hair differing in length and diameter while that of rex animals is essentially made up of amazingly soft down-hair. Rex short hair coat phenotypes in rabbits were shown to be controlled by three distinct loci. We focused on the "r1" mutation which segregates at a simple autosomal-recessive locus in our rabbit strains. A positional candidate gene approach was used to identify the rex gene and the corresponding mutation. The gene was primo-localized within a 40 cM region on rabbit chromosome 14 by genome scanning families of 187 rabbits in an experimental mating scheme. Then, fine mapping refined the region to 0.5 cM (Z = 78) by genotyping an additional 359 offspring for 94 microsatellites present or newly generated within the first defined interval. Comparative mapping pointed out a candidate gene in this 700 kb region, namely LIPH (Lipase Member H). In humans, several mutations in this major gene cause alopecia, hair loss phenotypes. The rabbit gene structure was established and a deletion of a single nucleotide was found in LIPH exon 9 of rex rabbits (1362delA). This mutation results in a frameshift and introduces a premature stop codon potentially shortening the protein by 19 amino acids. The association between this deletion and the rex phenotype was complete, as determined by its presence in our rabbit families and among a panel of 60 rex and its absence in all 60 non-rex rabbits. This strongly suggests that this deletion, in a homozygous state, is responsible for the rex phenotype in rabbits.

Gene expression profiling in equine polysaccharide storage myopathy revealed inflammation, glycogenesis inhibition, hypoxia and mitochondrial dysfunctions.

BACKGROUND: Several cases of myopathies have been observed in the horse Norman Cob breed. Muscle histology examinations revealed that some families suffer from a polysaccharide storage myopathy (PSSM). It is assumed that a gene expression signature related to PSSM should be observed at the transcriptional level because the glycogen storage disease could also be linked to other dysfunctions in gene regulation. Thus, the functional genomic approach could be conducted in order to provide new knowledge about the metabolic disorders related to PSSM. We propose exploring the PSSM muscle fiber metabolic disorders by measuring gene expression in relationship with the histological phenotype. RESULTS: Genotypying analysis of GYS1 mutation revealed 2 homozygous (AA) and 5 heterozygous (GA) PSSM horses. In the PSSM muscles, histological data revealed PAS positive amylase resistant abnormal polysaccharides, inflammation, necrosis, and lipomatosis and active regeneration of fibers. Ultrastructural evaluation revealed a decrease of mitochondrial number and structural disorders. Extensive accumulation of an abnormal polysaccharide displaced and partially replaced mitochondria and myofibrils. The severity of the disease was higher in the two homozygous PSSM horses.Gene expression analysis revealed 129 genes significantly modulated (p < 0.05). The following genes were up-regulated over 2 fold: IL18, CTSS, LUM, CD44, FN1, GST01. The most down-regulated genes were the following: mitochondrial tRNA, SLC2A2, PRKCalpha, VEGFalpha. Data mining analysis showed that protein synthesis, apoptosis, cellular movement, growth and proliferation were the main cellular functions significantly associated with the modulated genes (p < 0.05). Several up-regulated genes, especially IL18, revealed a severe muscular inflammation in PSSM muscles. The up-regulation of glycogen synthase kinase-3 (GSK3beta) under its active form could be responsible for glycogen synthase (GYS1) inhibition and hypoxia-inducible factor (HIF1alpha) destabilization. CONCLUSION: The main disorders observed in PSSM muscles could be related to mitochondrial dysfunctions, glycogenesis inhibition and the chronic hypoxia of the PSSM muscles.

Characterization of the equine glycogen debranching enzyme gene (AGL): Genomic and cDNA structure, localization, polymorphism and expression.

Glycogen debranching enzyme (AGL) is a multifunctional enzyme acting in the glycogen degradation pathway. In humans, the AGL activity deficiency causes a type III glycogen storage disease (Cori-Forbes disease). One particularity of AGL gene expression lies in the multiple alternative splicing in its 5' region. The AGL gene was localized on ECA5q14-q15. The sequence of the equine cDNA was determined to be 7.5 kb in length with an open reading frame of 4602 bp. The gene is 69 kb long and contains 35 exons. The equine AGL gene has an ubiquitous expression and presents five tissue-dependent cDNA variants arising from alternative splicing of the first exons. The equine skeletal muscle and heart contain four out of six variants previously described in humans and the equine liver express three of these four human variants. We identified a new alternative splicing variant expressed in equine skeletal and heart muscles. All these mRNA variants most probably encode only two different protein isoforms of 1533 and 1377 amino-acids. Four SNPs were detected in the mRNA. The equine in silico promoter sequence reveals a structure similar to those of other mammalian species. The disposition of the transcription factor biding sites does not correlate to the transcription start sites of tissue-specific variants.

Evolutionary movement of centromeres in horse, donkey, and zebra.

Centromere repositioning (CR) is a recently discovered biological phenomenon consisting of the emergence of a new centromere along a chromosome and the inactivation of the old one. After a CR, the primary constriction and the centromeric function are localized in a new position while the order of physical markers on the chromosome remains unchanged. These events profoundly affect chromosomal architecture. Since horses, asses, and zebras, whose evolutionary divergence is relatively recent, show remarkable morphological similarity and capacity to interbreed despite their chromosomes differing considerably, we investigated the role of CR in the karyotype evolution of the genus Equus. Using appropriate panels of BAC clones in FISH experiments, we compared the centromere position and marker order arrangement among orthologous chromosomes of Burchelli's zebra (Equus burchelli), donkey (Equus asinus), and horse (Equus caballus). Surprisingly, at least eight CRs took place during the evolution of this genus. Even more surprisingly, five cases of CR have occurred in the donkey after its divergence from zebra, that is, in a very short evolutionary time (approximately 1 million years). These findings suggest that in some species the CR phenomenon could have played an important role in karyotype shaping, with potential consequences on population dynamics and speciation.

cDNA sequence of the horse (Equus caballus) LAMA3 gene and characterization of two intronic SNP markers.

Laminins are large heterotrimeric basement membrane glycoproteins composed of alpha, beta and gamma chains. The Laminin 5 isoform has an alpha3beta3gamma2 composition and is essential for the adhesion of basal keratinocytes to the underlying epithelial basement membrane where it is mainly located. Mutations in the genes coding for the 3 chains have been associated with a severe skin blistering disease, Herlitz's junctional epidermolysis bullosa (JEB), observed in different species as man, dog, cat and horse. In this study, we report the sequence of the 5.2 kb horse laminin alpha 3 cDNA (LAMA3) as well as the detection of two intronic SNPs. These data will be useful to further identify causal mutations for the disease in this gene.

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.

A mutation in the LAMC2 gene causes the Herlitz junctional epidermolysis bullosa (H-JEB) in two French draft horse breeds.

Epidermolysis bullosa (EB) is a heterogeneous group of inherited diseases characterised by skin blistering and fragility. In humans, one of the most severe forms of EB known as Herlitz-junctional EB (H-JEB), is caused by mutations in the laminin 5 genes. EB has been described in several species, like cattle, sheep, dogs, cats and horses where the mutation, a cytosine insertion in exon 10 of the LAMC2 gene, was very recently identified in Belgian horses as the mutation responsible for JEB. In this study, the same mutation was found to be totally associated with the JEB phenotype in two French draft horse breeds, Trait Breton and Trait Comtois. This result provides breeders a molecular test to better manage their breeding strategies by genetic counselling.

A mutation in the MATP gene causes the cream coat colour in the horse.

In horses, basic colours such as bay or chestnut may be partially diluted to buckskin and palomino, or extremely diluted to cream, a nearly white colour with pink skin and blue eyes. This dilution is expected to be controlled by one gene and we used both candidate gene and positional cloning strategies to identify the "cream mutation". A horse panel including reference colours was established and typed for different markers within or in the neighbourhood of two candidate genes. Our data suggest that the causal mutation, a G to A transition, is localised in exon 2 of the MATP gene leading to an aspartic acid to asparagine substitution in the encoded protein. This conserved mutation was also described in mice and humans, but not in medaka.

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.

Molecular characterization of the equine testis-specific protein 1 (TPX1) and acidic epididymal glycoprotein 2 (AEG2) genes encoding members of the cysteine-rich secretory protein (CRISP) family.

The cysteine-rich secretory protein (CRISP) family consists of three members called acidic epididymal glycoprotein 1 (AEG1), AEG2, and testis-specific protein 1 (TPX1), which share 16 conserved cysteine residues at their C-termini. The CRISP proteins are primarily expressed in different sections of the male genital tract and are thought to mediate cell-cell interactions of male germ cells with other cells during sperm maturation or during fertilization. Therefore, their genes are of interest as candidate genes for inherited male fertility dysfunctions and as putative quantitative trait loci for male fertility traits. In this report, the cloning and DNA sequence of 137 kb of horse genomic DNA from equine chromosome 20q22 containing the closely linked equine TPX1 and AEG2 genes are described. The equine TPX1 gene consists of ten exons spanning 18 kb while the AEG2 gene consists of eight exons that are spread over 24 kb. The expression of these two genes was investigated in several tissues by reverse transcription polymerase chain reaction analysis and Western blotting. Comparative genome analysis between horse, human, and mouse indicates that all three CRISP genes are clustered on one chromosomal location, which shows conserved synteny between these species.

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.

A specific pattern of splicing for the horse alphaS1-Casein mRNA and partial genomic characterization of the relevant locus.

Mares' milk has a composition very different from that of cows' milk. It is much more similar to human milk, in particular in its casein fraction. This study reports on the sequence of a 994 bp amplified fragment corresponding to a horse alphaS1-Casein (alphaS1-Cn) cDNA and its comparison with its caprine, pig, rabbit and human counterparts. The alignment of these sequences revealed a specific pattern of splicing for this horse primary transcript. As in humans, exons 3', 6' and 13' are present whereas exons 5, 13 and 14 are absent in this equine mRNA sequence. BAC clones, screened from a horse BAC library, containing the alphaS1-Cn gene allowed the mapping of its locus by FISH on equine chromosome 3q22.2-q22.3 which is in agreement with the Zoo-FISH results. Genomic analysis of the alphaS1-Cn gene showed that the region from the second exon to the last exon is scattered within a nucleotide stretch nearly 15-kb in length which is quite similar in size to its ruminant and rabbit counterparts. The region between alphaS1- and beta-Cn genes, suspected to contain cis-acting elements involved in the expression of all clustered casein genes, is similar in size (ca 15-kb) to the caprine and mouse intergenic region.