A high-resolution physical map of equine homologs of HSA19 shows divergent evolution compared with other mammals

Subchromosomal karyotype evolution in Equidae

Equidae is a small family which comprises horses, African and Asiatic asses, and zebras. Despite equids having diverged quite recently, their karyotypes underwent rapid evolution which resulted in extensive differences among chromosome complements in respective species. Comparative mapping using whole-chromosome painting probes delineated genome-wide chromosome homologies among extant equids, enabling us to trace chromosome rearrangements that occurred during evolution. In the present study, we performed subchromosomal comparative mapping among seven Equidae species, representing the whole family. Region-specific painting and bacterial artificial chromosome probes were used to determine the orientation of evolutionarily conserved segments with respect to centromere positions. This allowed assessment of the configuration of all fusions occurring during the evolution of Equidae, as well as revealing discrepancies in centromere location caused by centromere repositioning or inversions. Our results indicate that the prevailing type of fusion in Equidae is centric fusion. Tandem fusions of the type telomere-telomere occur almost exclusively in the karyotype of Hartmann's zebra and are characteristic of this species' evolution. We revealed inversions in segments homologous to horse chromosomes 3p/10p and 13 in zebras and confirmed inversions in segments 4/31 in African ass, 7 in horse and 8p/20 in zebras. Furthermore, our mapping results suggested that centromere repositioning events occurred in segments homologous to horse chromosomes 7, 8q, 10p and 19 in the African ass and an element homologous to horse chromosome 16 in Asiatic asses. Centromere repositioning in chromosome 1 resulted in three different chromosome types occurring in extant species. Heterozygosity of the centromere position of this chromosome was observed in the kiang. Other subtle changes in centromere position were described in several evolutionary conserved chromosomal segments, suggesting that tiny centromere repositioning or pericentric inversions are quite frequent in zebras and asses.

Karyotypic relationships in Asiatic asses (kulan and kiang) as defined using horse chromosome arm-specific and region-specific probes

Cross-species chromosome painting has been applied to most of the species making up the numerically small family Equidae. However, comparative mapping data were still lacking in Asiatic asses kulan (Equus hemionus kulan) and kiang (E. kiang). The set of horse arm-specific probes generated by laser microdissection was hybridized onto kulan (E. hemionus kulan) and kiang (E. kiang) chromosomes in order to establish a genome-wide chromosomal correspondence between these Asiatic asses and the horse. Moreover, region-specific probes were generated to determine fusion configuration and orientation of conserved syntenic blocks. The kulan karyotype (2n = 54) was ascertained to be almost identical to the previously investigated karyotype of onager E. h. onager (2n = 56). The only difference is in fusion/fission of chromosomes homologous to horse 2q/3q, which are involved in chromosome number polymorphism in many Equidae species. E. kiang karyotype differs from the karyotype of E. hemionus by two additional fusions 8q/15 and 7/25. Chromosomes equivalent to 2q and 3q are not fused in kiang individuals with 2n = 52. Several discrepancies in centromere positions among kulan, kiang and horse chromosomes have been described. Most of the chromosome fusions in Asiatic asses are of centromere-centromere type. Comparative chromosome painting in kiang completed the efforts to establish chromosomal homologies in all representatives of the family Equidae. Application of region-specific probes allows refinement comparative maps of Asiatic asses.

A high resolution RH map of the bovine major histocompatibility complex

Background

The cattle MHC is termed the bovine leukocyte antigen (BoLA) and, along with the MHCs of other ruminants, is unique in its genomic organization. Consequently, correct and reliable gene maps and sequence information are critical to the study of the BoLA region. The bovine genome sequencing project has produced two assemblies (Btau_3.1 and 4.0) that differ substantially from each other and from conventional gene maps in the BoLA region. To independently compare the accuracies of the different sequence assemblies, we have generated a high resolution map of BoLA using a 12,000_rad radiation hybrid panel. Seventy-seven unique sequence tagged site (STS) markers chosen at approximately 50 kb intervals from the Btau 2.0 assembly and spanning the IIa-III-I and IIb regions of the bovine MHC were mapped on a 12,000_rad bovine radiation hybrid (RH) panel to evaluate the different assemblies of the bovine genome sequence.

Results

Analysis of the data generated a high resolution RH map of BoLA that was significantly different from the Btau_3.1 assembly of the bovine genome but in good agreement with the Btau_4.0 assembly. Of the few discordancies between the RH map and Btau_4.0, most could be attributed to closely spaced markers that could not be precisely ordered in the RH panel. One probable incorrectly-assembled sequence and three missing sequences were noted in the Btau_4.0 assembly. The RH map of BoLA is also highly concordant with the sequence-based map of HLA (NCBI build 36) when reordered to account for the ancestral inversion in the ruminant MHC.

Conclusion

These results strongly suggest that studies using Btau_3.1 for analyses of the BoLA region should be reevaluated in light of the Btau_4.0 assembly and indicate that additional research is needed to produce a complete assembly of the BoLA genomic sequences.

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.

Glycogen synthase (GYS1) mutation causes a novel skeletal muscle glycogenosis

Polysaccharide storage myopathy (PSSM) is a novel glycogenosis in horses characterized by abnormal glycogen accumulation in skeletal muscle and muscle damage with exertion. It is unlike glycogen storage diseases resulting from known defects in glycogenolysis, glycolysis, and glycogen synthesis that have been described in humans and domestic animals. A genome-wide association identified GYS1, encoding skeletal muscle glycogen synthase (GS), as a candidate gene for PSSM. DNA sequence analysis revealed a mutation resulting in an arginine-to-histidine substitution in a highly conserved region of GS. Functional analysis demonstrated an elevated GS activity in PSSM horses, and haplotype analysis and allele age estimation demonstrated that this mutation is identical by descent among horse breeds. This is the first report of a gain-of-function mutation in GYS1 resulting in a glycogenosis.

Construction of radiation hybrid panels

Whole-genome radiation hybrid (RH) mapping has proven to be a powerful tool for mapping genes and comparing genome architecture. We describe a protocol for constructing RH panels by rescuing irradiated fibroblast donor cells of any mammalian species by polyethylene glycol fusion to a thymidine kinase-deficient hamster cell line. Characterization and expansion of a panel of 90-100 cell lines can be used to map virtually any PCR-based marker that can be distinguished from the recipient hamster genome. The described procedure has been used successfully to create RH panels from diverse mammalian species such as macaques, elephants, alpacas, and armadillos, and may be applicable to nonmammalian vertebrates as well.

Karyotypic relationships among Equus grevyi, Equus burchelli and domestic horse defined using horse chromosome arm-specific probes

Using laser microdissection we prepared a set of horse chromosome arm-specific probes. Most of the probes were generated from horse chromosomes, some of them were derived from Equus zebra hartmannae. The set of probes were hybridized onto E. grevyi chromosomes in order to establish a genome-wide chromosomal correspondence between this zebra and horse. The use of arm-specific probes provided us with more information on the mutual arrangement of the genomes than we could obtain by means of whole-chromosome paints generated by flow sorting, even if we used reciprocal painting with probe sets from both species. By comparison of our results and results of comparative mapping in E. burchelli, we also established the chromosomal correspondence between E. grevyi and E. burchelli, providing evidence for a very close karyotypic relationship between these two zebra species. Establishment of the comparative map for E. grevyi contributes to the knowledge of the karyotypic phylogeny in the Equidae family.

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.

Construction of a river buffalo (Bubalus bubalis) whole-genome radiation hybrid panel and preliminary RH mapping of chromosomes 3 and 10

The buffalo (Bubalus bubalis) is a source of milk and meat, and also serves as a draft animal. In this study, a 5000-rad whole-genome radiation hybrid (RH) panel for river buffalo was constructed and used to build preliminary RH maps for BBU3 and BBU10 chromosomes. The preliminary maps contain 66 markers, including coding genes, cattle expressed sequence tags (ESTs) and microsatellite loci. The RH maps presented here are the starting point for mapping additional loci that will allow detailed comparative maps between buffalo, cattle and other species whose genomes may be mapped in the future. A large quantity of DNA has been prepared from the cell lines forming the river buffalo RH panel and will be made publicly available to the international community both for the study of chromosome evolution and for the improvement of traits important to the role of buffalo in animal agriculture.

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.

A 1.5-Mb-resolution radiation hybrid map of the cat genome and comparative analysis with the canine and human genomes

We report the construction of a 1.5-Mb-resolution radiation hybrid map of the domestic cat genome. This new map includes novel microsatellite loci and markers derived from the 2X genome sequence that target previous gaps in the feline-human comparative map. Ninety-six percent of the 1793 cat markers we mapped have identifiable orthologues in the canine and human genome sequences. The updated autosomal and X-chromosome comparative maps identify 152 cat-human and 134 cat-dog homologous synteny blocks. Comparative analysis shows the marked change in chromosomal evolution in the canid lineage relative to the felid lineage since divergence from their carnivoran ancestor. The canid lineage has a 30-fold difference in the number of interchromosomal rearrangements relative to felids, while the felid lineage has primarily undergone intrachromosomal rearrangements. We have also refined the pseudoautosomal region and boundary in the cat and show that it is markedly longer than those of human or mouse. This improved RH comparative map provides a useful tool to facilitate positional cloning studies in the feline model.

Design factors that influence PCR amplification success of cross-species primers among 1147 mammalian primer pairs

Background

Cross-species primers have been used with moderate success to address a variety of questions concerning genome structure, evolution, and gene function. However, the factors affecting their success have never been adequately addressed, particularly with respect to producing a consistent method to achieve high throughput. Using 1,147 mammalian cross-species primer pairs (1089 not previously reported), we tested several factors to determine their influence on the probability that a given target will amplify in a given species under a single amplification condition. These factors included: number of mismatches between the two species (the index species) used to identify conserved regions to which the primers were designed, GC-content of the gene and amplified region, CpG dinucleotides in the primer region, degree of encoded protein conservation, length of the primers, and the degree of evolutionary distance between the target species and the two index species.

Results

The amplification success rate for the cross-species primers was significantly influenced by the number of mismatches between the two index species (6–8% decrease per mismatch in a primer pair), the GC-content within the amplified region (for the dog, GC ≥ 50%, 56.9% amplified; GC<50%, 74.2% amplified), the degree of protein conservation (R² = 0.14) and the relatedness of the target species to the index species. For the dog, 598 products of 930 primer pairs (64.3%) (excluding primers in which dog was an index species) were sequenced and shown to be the expected product, with an additional three percent producing the incorrect sequence. When hamster DNA was used with the single amplification condition in a microtiter plate-based format, 510 of 1087 primer pairs (46.9%) produced amplified products. The primer pairs are spaced at an average distance of 2.3 Mb in the human genome and may be used to produce up to several hundred thousand bp of species-specific sequence.

Conclusion

The most important factors influencing the proportion of successful amplifications are the number of index species mismatches, GC-richness of the target amplimer, and the relatedness of the target species to the index species, at least under the single PCR condition used. The 1147 cross-species primer pairs can be used in a high throughput manner to generate data for studies on the genetics and genomics of non-sequenced mammalian genomes.

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.

Comparative genome maps of the pangolin, hedgehog, sloth, anteater and human revealed by cross-species chromosome painting: further insight into the ancestral karyotype and genome evolution of eutherian mammals

To better understand the evolution of genome organization of eutherian mammals, comparative maps based on chromosome painting have been constructed between human and representative species of three eutherian orders: Xenarthra, Pholidota, and Eulipotyphla, as well as between representative species of the Carnivora and Pholidota. These maps demonstrate the conservation of such syntenic segment associations as HSA3/21, 4/8, 7/16, 12/22, 14/15 and 16/19 in Eulipotyphla, Pholidota and Xenarthra and thus further consolidate the notion that they form part of the ancestral karyotype of the eutherian mammals. Our study has revealed many potential ancestral syntenic associations of human chromosomal segments that serve to link the families as well as orders within the major superordinial eutherian clades defined by molecular markers. The HSA2/8 and 7/10 associations could be the cytogenetic signatures that unite the Xenarthrans, while the HSA1/19p could be a putative signature that links the Afrotheria and Xenarthra. But caution is required in the interpretation of apparently shared syntenic associations as detailed analyses also show examples of apparent convergent evolution that differ in breakpoints and extent of the involved segments.

A human-horse comparative map based on equine BAC end sequences

In an effort to increase the density of sequence-based markers for the horse genome we generated 9473 BAC end sequences (BESs) from the CHORI-241 BAC library with an average read length of 677 bp. BLASTN searches with the BESs revealed 4036 meaningful hits (E <or= 10(-5)) in the human genome that provide useful markers for the human-horse comparative map. The 4036 BLASTN hits allowed the anchoring of 3079 BAC clones to the human genome, on average one corresponding equine BAC clone per megabase of human DNA. We used the BLASTN anchored BESs for an in silico prediction of the gene content and chromosome assignment of comparatively mapped equine BAC clones. As a first verification of our in silico mapping strategy we placed 19 equine BESs with matches to HSA6 onto the RH map. All markers were assigned to the predicted localizations on ECA10, ECA20, and ECA31, respectively.

A 1.3-Mb interval map of equine homologs of HSA2

A comparative approach that utilizes information from more densely mapped or sequenced genomes is a proven and efficient means to increase our knowledge of the structure of the horse genome. Human chromosome 2 (HSA2), the second largest human chromosome, comprising 243 Mb, and containing 1246 known genes, corresponds to all or parts of three equine chromosomes. This report describes the assignment of 140 new markers (78 genes and 62 microsatellites) to the equine radiation hybrid (RH) map, and the anchoring of 24 of these markers to horse chromosomes by FISH. The updated equine RH maps for ECA6p, ECA15, and ECA18 resulting from this work have one, two, and three RH linkage groups, respectively, per chromosome/chromosome-arm. These maps have a three-fold increase in the number of mapped markers compared to previous maps of these chromosomes, and an increase in the average marker density to one marker per 1.3 Mb. Comparative maps of ECA6p, ECA15, and ECA18 with human, chimpanzee, dog, mouse, rat, and chicken genomes reveal blocks of conserved synteny across mammals and vertebrates.