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

Unlocking Horse Y Chromosome Diversity

The present equine genetic variation mirrors the deep influence of intensive breeding programs during the last 200 years. Here, we provide a comprehensive current state of knowledge on the trends and prospects on the variation in the equine male-specific region of the Y chromosome (MSY), which was assembled for the first time in 2018. In comparison with the other 12 mammalian species, horses are now the most represented, with 56 documented MSY genes. However, in contrast to the high variability in mitochondrial DNA observed in many horse breeds from different geographic areas, modern horse populations demonstrate extremely low genetic Y-chromosome diversity. The selective pressures employed by breeders using pedigree data (which are not always error-free) as a predictive tool represent the main cause of this lack of variation in the Y-chromosome. Nevertheless, the detailed phylogenies obtained by recent fine-scaled Y-chromosomal genotyping in many horse breeds worldwide have contributed to addressing the genealogical, forensic, and population questions leading to the reappraisal of the Y-chromosome as a powerful genetic marker to avoid the loss of biodiversity as a result of selective breeding practices, and to better understand the historical development of horse breeds.

Horse Clinical Cytogenetics: Recurrent Themes and Novel Findings

Clinical cytogenetic studies in horses have been ongoing for over half a century and clearly demonstrate that chromosomal disorders are among the most common non-infectious causes of decreased fertility, infertility, and congenital defects. Large-scale cytogenetic surveys show that almost 30% of horses with reproductive or developmental problems have chromosome aberrations, whereas abnormal karyotypes are found in only 2-5% of the general population. Among the many chromosome abnormalities reported in the horse, most are unique or rare. However, all surveys agree that there are two recurrent conditions: X-monosomy and SRY-negative XY male-to-female sex reversal, making up approximately 35% and 11% of all chromosome abnormalities, respectively. The two are signature conditions for the horse and rare or absent in other domestic species. The progress in equine genomics and the development of molecular tools, have qualitatively improved clinical cytogenetics today, allowing for refined characterization of aberrations and understanding the underlying molecular mechanisms. While cutting-edge genomics tools promise further improvements in chromosome analysis, they will not entirely replace traditional cytogenetics, which still is the most straightforward, cost-effective, and fastest approach for the initial evaluation of potential breeding animals and horses with reproductive or developmental disorders.

Ten years of the horse reference genome: insights into equine biology, domestication and population dynamics in the post-genome era

The horse reference genome from the Thoroughbred mare Twilight has been available for a decade and, together with advances in genomics technologies, has led to unparalleled developments in equine genomics. At the core of this progress is the continuing improvement of the quality, contiguity and completeness of the reference genome, and its functional annotation. Recent achievements include the release of the next version of the reference genome (EquCab3.0) and generation of a reference sequence for the Y chromosome. Horse satellite‐free centromeres provide unique models for mammalian centromere research. Despite extremely low genetic diversity of the Y chromosome, it has been possible to trace patrilines of breeds and pedigrees and show that Y variation was lost in the past approximately 2300 years owing to selective breeding. The high‐quality reference genome has led to the development of three different SNP arrays and WGSs of almost 2000 modern individual horses. The collection of WGS of hundreds of ancient horses is unique and not available for any other domestic species. These tools and resources have led to global population studies dissecting the natural history of the species and genetic makeup and ancestry of modern breeds. Most importantly, the available tools and resources, together with the discovery of functional elements, are dissecting molecular causes of a growing number of Mendelian and complex traits. The improved understanding of molecular underpinnings of various traits continues to benefit the health and performance of the horse whereas also serving as a model for complex disease across species.

Satellite DNA at the Centromere is Dispensable for Segregation Fidelity

The typical vertebrate centromeres contain long stretches of highly repeated DNA sequences (satellite DNA). We previously demonstrated that the karyotypes of the species belonging to the genus Equus are characterized by the presence of satellite-free and satellite-based centromeres and represent a unique biological model for the study of centromere organization and behavior. Using horse primary fibroblasts cultured in vitro, we compared the segregation fidelity of chromosome 11, whose centromere is satellite-free, with that of chromosome 13, which has similar size and a centromere containing long stretches of satellite DNA. The mitotic stability of the two chromosomes was compared under normal conditions and under mitotic stress induced by the spindle inhibitor, nocodazole. Two independent molecular-cytogenetic approaches were used—the interphase aneuploidy analysis and the cytokinesis-block micronucleus assay. Both assays were coupled to fluorescence in situ hybridization with chromosome specific probes in order to identify chromosome 11 and chromosome 13, respectively. In addition, we tested if the lack of centromeric satellite DNA affected chromatid cohesion under normal and stress conditions. We demonstrated that, in our system, the segregation fidelity of a chromosome is not influenced by the presence of long stretches of tandem repeats at its centromere. To our knowledge, the present study is the first analysis of the mitotic behavior of a natural satellite-free centromere.

Birth, evolution, and transmission of satellite-free mammalian centromeric domains

Mammalian centromeres are associated with highly repetitive DNA (satellite DNA), which has so far hindered molecular analysis of this chromatin domain. Centromeres are epigenetically specified, and binding of the CENPA protein is their main determinant. In previous work, we described the first example of a natural satellite-free centromere on Equus caballus Chromosome 11. Here, we investigated the satellite-free centromeres of Equus asinus by using ChIP-seq with anti-CENPA antibodies. We identified an extraordinarily high number of centromeres lacking satellite DNA (16 of 31). All of them lay in LINE- and AT-rich regions. A subset of these centromeres is associated with DNA amplification. The location of CENPA binding domains can vary in different individuals, giving rise to epialleles. The analysis of epiallele transmission in hybrids (three mules and one hinny) showed that centromeric domains are inherited as Mendelian traits, but their position can slide in one generation. Conversely, centromere location is stable during mitotic propagation of cultured cells. Our results demonstrate that the presence of more than half of centromeres void of satellite DNA is compatible with genome stability and species survival. The presence of amplified DNA at some centromeres suggests that these arrays may represent an intermediate stage toward satellite DNA formation during evolution. The fact that CENPA binding domains can move within relatively restricted regions (a few hundred kilobases) suggests that the centromeric function is physically limited by epigenetic boundaries.

Mitochondrial DNA insertions in the nuclear horse genome

The insertion of mitochondrial DNA in the nuclear genome generates numts, nuclear sequences of mitochondrial origin. In the horse reference genome, we identified 82 numts and showed that the entire horse mitochondrial DNA is represented as numts without gross bias. Numts were inserted in the horse nuclear genome at random sites and were probably generated during the repair of DNA double-strand breaks. We then analysed 12 numt loci in 20 unrelated horses and found that null alleles, lacking the mitochondrial DNA insertion, were present at six of these loci. At some loci, the null allele is prevalent in the sample analysed, suggesting that, in the horse population, the number of numt loci may be higher than 82 present in the reference genome. Contrary to humans, the insertion polymorphism of numts is extremely frequent in the horse population, supporting the hypothesis that the genome of this species is in a rapidly evolving state.

Centromere sliding on a mammalian chromosome

The centromere directs the segregation of chromosomes during mitosis and meiosis. It is a distinct genetic locus whose identity is established through epigenetic mechanisms that depend on the deposition of centromere-specific centromere protein A (CENP-A) nucleosomes. This important chromatin domain has so far escaped comprehensive molecular analysis due to its typical association with highly repetitive satellite DNA. In previous work, we discovered that the centromere of horse chromosome 11 is completely devoid of satellite DNA; this peculiar feature makes it a unique model to dissect the molecular architecture of mammalian centromeres. Here, we exploited this native satellite-free centromere to determine the precise localization of its functional domains in five individuals: We hybridized DNA purified from chromatin immunoprecipitated with an anti CENP-A antibody to a high resolution array (ChIP-on-chip) of the region containing the primary constriction of horse chromosome 11. Strikingly, each individual exhibited a different arrangement of CENP-A binding domains. We then analysed the organization of each domain using a single nucleotide polymorphism (SNP)-based approach and single molecule analysis on chromatin fibres. Examination of the ten instances of chromosome 11 in the five individuals revealed seven distinct ‘positional alleles’, each one extending for about 80–160 kb, were found across a region of about 500 kb. Our results demonstrate that CENP-A binding domains are autonomous relative to the underlying DNA sequence and are characterized by positional instability causing the sliding of centromere position. We propose that this dynamic behaviour may be common in mammalian centromeres and may determine the establishment of epigenetic alleles.

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.

Whole genome comparative analysis of channel catfish (Ictalurus punctatus) with four model fish species

BACKGROUND

Comparative mapping is a powerful tool to study evolution of genomes. It allows transfer of genome information from the well-studied model species to non-model species. Catfish is an economically important aquaculture species in United States. A large amount of genome resources have been developed from catfish including genetic linkage maps, physical maps, BAC end sequences (BES), integrated linkage and physical maps using BES-derived markers, physical map contig-specific sequences, and draft genome sequences. Application of such genome resources should allow comparative analysis at the genome scale with several other model fish species.

RESULTS

In this study, we conducted whole genome comparative analysis between channel catfish and four model fish species with fully sequenced genomes, zebrafish, medaka, stickleback and Tetraodon. A total of 517 Mb draft genome sequences of catfish were anchored to its genetic linkage map, which accounted for 62% of the total draft genome sequences. Based on the location of homologous genes, homologous chromosomes were determined among catfish and the four model fish species. A large number of conserved syntenic blocks were identified. Analysis of the syntenic relationships between catfish and the four model fishes supported that the catfish genome is most similar to the genome of zebrafish.

CONCLUSION

The organization of the catfish genome is similar to that of the four teleost species, zebrafish, medaka, stickleback, and Tetraodon such that homologous chromosomes can be identified. Within each chromosome, extended syntenic blocks were evident, but the conserved syntenies at the chromosome level involve extensive inter-chromosomal and intra-chromosomal rearrangements. This whole genome comparative map should facilitate the whole genome assembly and annotation in catfish, and will be useful for genomic studies of various other fish species.

Genomic analysis and mRNA expression of equine type I interferon genes

This study aimed at identifying all of the type I interferon (IFN) genes of the horse and at monitoring their expression in equine cells on in vitro induction. We identified 32 putative type I IFN loci on horse chromosome 23 and an unplaced genomic scaffold. A phylogentic analysis characterized these into 8 different type I IFN classes, that is, putative functional genes for 6 IFN-α, 4 IFN-β, 8 IFN-ω (plus 4 pseudogenes), 3 IFN-δ (plus 1 pseudogene), 1 IFN-κ and 1 IFN-ε, plus 1 IFN-ν pseudogene, and 3 loci belonging to what has previously been called IFN-αω. Our analyses indicate that the IFN-αω genes are quite distinct from both IFN-α and IFN-ω, and we refer to this type I IFN as IFN-μ. Results from cell cultures showed that leukocytes readily expressed IFN-α, IFN-β, IFN-δ, IFN-μ, and IFN-ω mRNA on induction with, for example, live virus; while fibroblasts only expressed IFN-β mRNA on stimulation. IFN-κ or IFN-ε expression was not consistently induced in these cell cultures. Thus, the equine type I IFN family comprised 8 classes, 7 of which had putative functional genes, and mRNA expression of 5 was induced in vitro. Moreover, a relatively low number of IFN-α subtypes was found in the horse compared with other eutherian mammals.

Comparative genomic analysis of catfish linkage group 8 reveals two homologous chromosomes in zebrafish and other teleosts with extensive inter-chromosomal rearrangements

Background

Comparative genomics is a powerful tool to transfer genomic information from model species to related non-model species. Channel catfish (Ictalurus punctatus) is the primary aquaculture species in the United States. Its existing genome resources such as genomic sequences generated from next generation sequencing, BAC end sequences (BES), physical maps, linkage maps, and integrated linkage and physical maps using BES-associated markers provide a platform for comparative genomic analysis between catfish and other model teleost fish species. This study aimed to gain understanding of genome organizations and similarities among catfish and several sequenced teleost genomes using linkage group 8 (LG8) as a pilot study.

Results

With existing genome resources, 287 unique genes were identified in LG8. Comparative genome analysis indicated that most of these 287 genes on catfish LG8 are located on two homologous chromosomes of zebrafish, medaka, stickleback, and three chromosomes of green-spotted pufferfish. Large numbers of conserved syntenies were identified. Detailed analysis of the conserved syntenies in relation to chromosome level similarities revealed extensive inter-chromosomal and intra-chromosomal rearrangements during evolution. Of the 287 genes, 35 genes were found to be duplicated in the catfish genome, with the vast majority of the duplications being interchromosomal.

Conclusions

Comparative genome analysis is a powerful tool even in the absence of a well-assembled whole genome sequence. In spite of sequence stacking due to low resolution of the linkage and physical maps, conserved syntenies can be identified although the exact gene order and orientation are unknown at present. Through chromosome-level comparative analysis, homologous chromosomes among teleosts can be identified. Syntenic analysis should facilitate annotation of the catfish genome, which in turn, should facilitate functional inference of genes based on their orthology.

Construction and preliminary characterization analysis of Wuzhishan miniature pig bacterial artificial chromosome library with approximately 8-fold genome equivalent coverage

Bacterial artificial chromosome (BAC) libraries have been invaluable tools for the genome-wide genetic dissection of complex organisms. Here, we report the construction and characterization of a high-redundancy BAC library from a very valuable pig breed in China, Wuzhishan miniature pig (Sus scrofa), using its blood cells and fibroblasts, respectively. The library contains approximately 153,600 clones ordered in 40 superpools of 10 × 384-deep well microplates. The average insert size of BAC clones was estimated to be 152.3 kb, representing approximately 7.68 genome equivalents of the porcine haploid genome and a 99.93% statistical probability of obtaining at least one clone containing a unique DNA sequence in the library. 19 pairs of microsatellite marker primers covering porcine chromosomes were used for screening the BAC library, which showed that each of these markers was positive in the library; the positive clone number was 2 to 9, and the average number was 7.89, which was consistent with 7.68-fold coverage of the porcine genome. And there were no significant differences of genomic BAC library from blood cells and fibroblast cells. Therefore, we identified 19 microsatellite markers that could potentially be used as genetic markers. As a result, this BAC library will serve as a valuable resource for gene identification, physical mapping, and comparative genomics and large-scale genome sequencing in the porcine.

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.

Construction and characterization of the BAC library for common carp Cyprinus carpio L. and establishment of microsynteny with zebrafish Danio rerio

A bacterial artificial chromosome (BAC) library of common carp Cyprinus carpio L. was constructed as a part of ongoing common carp genome project, which is aiming assembly of common carp genome. The library, containing a total of 92,160 BAC clones with an average insert size of 141 kb, was constructed into the restriction site of Hind III on BAC vector CopyControl pCC1BAC, covering 7.7 X haploid genome equivalents. Three dimension pools and superpools of the BAC library were established and 23 positive clones of 14 targets were identified from one-fifth of the BAC library. Pilot project of BAC end sequencing was conducted on 2,688 BAC ends from 1,344 clones and harvested 2,522 high-quality Q20 sequences with average length of 677 bp. The sequencing success rate was 93.8% and pair-end success rate was 92.3%. A total of 212 microsyntenies had been established between common carp and zebrafish genomes as a trial for genome-wide comparative genomics in these two closely related species.

Construction of an American mink bacterial artificial chromosome (BAC) library and sequencing candidate genes important for the fur industry

Background

Bacterial artificial chromosome (BAC) libraries continue to be invaluable tools for the genomic analysis of complex organisms. Complemented by the newly and fast growing deep sequencing technologies, they provide an excellent source of information in genomics projects.

Results

Here, we report the construction and characterization of the CHORI-231 BAC library constructed from a Danish-farmed, male American mink (Neovison vison). The library contains approximately 165,888 clones with an average insert size of 170 kb, representing approximately 10-fold coverage. High-density filters, each consisting of 18,432 clones spotted in duplicate, have been produced for hybridization screening and are publicly available. Overgo probes derived from expressed sequence tags (ESTs), representing 21 candidate genes for traits important for the mink industry, were used to screen the BAC library. These included candidate genes for coat coloring, hair growth and length, coarseness, and some receptors potentially involved in viral diseases in mink. The extensive screening yielded positive results for 19 of these genes. Thirty-five clones corresponding to 19 genes were sequenced using 454 Roche, and large contigs (184 kb in average) were assembled. Knowing the complete sequences of these candidate genes will enable confirmation of the association with a phenotype and the finding of causative mutations for the targeted phenotypes.

Additionally, 1577 BAC clones were end sequenced; 2505 BAC end sequences (80% of BACs) were obtained. An excess of 2 Mb has been analyzed, thus giving a snapshot of the mink genome.

Conclusions

The availability of the CHORI-321 American mink BAC library will aid in identification of genes and genomic regions of interest. We have demonstrated how the library can be used to identify specific genes of interest, develop genetic markers, and for BAC end sequencing and deep sequencing of selected clones. To our knowledge, this is the first report of 454 sequencing of selected BAC clones in mammals and re-assures the suitability of this technique for obtaining the sequence information of genes of interest in small genomics projects. The BAC end sequences described in this paper have been deposited in the GenBank data library [HN339419-HN341884, HN604664-HN604702]. The 454 produced contigs derived from selected clones are deposited with reference numbers [GenBank: JF288166-JF288183 &JF310744].

Comparative genomics in teleost species: Knowledge transfer by linking the genomes of model and non-model fish species

Comparative genomics is a powerful tool to transfer knowledge coming from model fish species to non-model fish species of economic or/and evolutionary interest. Such transfer is of importance as functional studies either are difficult to perform with most non-model species. The first comparative map constructed using the human and the chimpanzee genome allowed the identification of putative orthologues. Although comparative mapping in teleosts is still in its infancy, five model teleost genomes from different orders have been fully sequenced to date and the sequencing of several commercially important species are also underway or near completion. The accessibility of these whole genome sequences and rapid developments in genomics of fish species are paving the way towards new and valuable research in comparative genetics and genomics. With the accumulation of information in model species, the genetic and genomic characterization of non-model, but economically, physiologically or evolutionary important species is now feasible. Furthermore, comparison of low coverage gene maps of non-model fish species against fully sequenced fish species will enhance the efficiency of candidate gene identification projected for quantitative trait loci (QTL) scans for traits of special interest.

[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.

Comparative analysis of catfish BAC end sequences with the zebrafish genome

Background

Comparative mapping is a powerful tool to transfer genomic information from sequenced genomes to closely related species for which whole genome sequence data are not yet available. However, such an approach is still very limited in catfish, the most important aquaculture species in the United States. This project was initiated to generate additional BAC end sequences and demonstrate their applications in comparative mapping in catfish.

Results

We reported the generation of 43,000 BAC end sequences and their applications for comparative genome analysis in catfish. Using these and the additional 20,000 existing BAC end sequences as a resource along with linkage mapping and existing physical map, conserved syntenic regions were identified between the catfish and zebrafish genomes. A total of 10,943 catfish BAC end sequences (17.3%) had significant BLAST hits to the zebrafish genome (cutoff value ≤ e^-5), of which 3,221 were unique gene hits, providing a platform for comparative mapping based on locations of these genes in catfish and zebrafish. Genetic linkage mapping of microsatellites associated with contigs allowed identification of large conserved genomic segments and construction of super scaffolds.

Conclusion

BAC end sequences and their associated polymorphic markers are great resources for comparative genome analysis in catfish. Highly conserved chromosomal regions were identified to exist between catfish and zebrafish. However, it appears that the level of conservation at local genomic regions are high while a high level of chromosomal shuffling and rearrangements exist between catfish and zebrafish genomes. Orthologous regions established through comparative analysis should facilitate both structural and functional genome analysis in catfish.

Interleukin 4 receptor alpha (IL4R) and calcium-activated chloride channel 1 (CLCA1) genes map to donkey chromosome

The results obtained in the present study enabled the physical map of the donkey genome to be extended with markers associated with recurrent airway obstruction (RAO), a major performance-limiting disease of Equidae. The equine BAC clone containing the IL4R and CLCA1 genes were localized to EAS 14q13 and EAS 6q15 respectlivy by fluorescent in situ hybridization. Identification of their locus confirmed the distribution of syntenic regions between the domestic horse and the domestic donkey within the chromosomes analysed.

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.

BAC library construction and BAC end sequencing of five Drosophila species: the comparative map with the D. melanogaster genome

We constructed and characterized arrayed bacterial artificial chromosome (BAC) libraries of five Drosophila species (D. melanogaster, D. simulans, D. sechellia, D. auraria, and D. ananassae), which are genetically well characterized in the studies of meiosis, evolution, population genetics, and developmental biology. The BAC libraries comprise 8,000 to 12,500 clones for each species, estimated to cover the most of the genomes. We sequenced both ends of most of these BAC clones with a success rate of 91%. Of these, 53,701 clones consisting of non-repetitive BAC end sequences (BESs) were mapped with reference of the public D. melanogaster genome sequences. The BES mapping estimated that the BAC libraries of D. auraria and D. ananassae covered 47% and 57% of the D. melanogaster genome, respectively, and those of D. melanogaster, D. sechellia, and D. simulans covered 94-97%. The low coverage by BESs of D. auraria and D. ananassae may be due to the high sequence divergence with D. melanogaster. From the comparative BES mapping, 111 possible breakpoints of chromosomal rearrangements were identified in these four species. The breakpoints of the major chromosome rearrangement between D. simulans and D. melanogaster on the third chromosome were determined within 20 kb in 84E and 30 kb in 93E/F. Corresponding breakpoints were also identified in D. sechellia. The BAC clones described here will be an important addition to the Drosophila genomic resources.

Allelic heterogeneity at the equine KIT locus in dominant white (W) horses

White coat color has been a highly valued trait in horses for at least 2,000 years. Dominant white (W) is one of several known depigmentation phenotypes in horses. It shows considerable phenotypic variation, ranging from approximately 50% depigmented areas up to a completely white coat. In the horse, the four depigmentation phenotypes roan, sabino, tobiano, and dominant white were independently mapped to a chromosomal region on ECA 3 harboring the KIT gene. KIT plays an important role in melanoblast survival during embryonic development. We determined the sequence and genomic organization of the approximately 82 kb equine KIT gene. A mutation analysis of all 21 KIT exons in white Franches-Montagnes Horses revealed a nonsense mutation in exon 15 (c.2151C>G, p.Y717X). We analyzed the KIT exons in horses characterized as dominant white from other populations and found three additional candidate causative mutations. Three almost completely white Arabians carried a different nonsense mutation in exon 4 (c.706A>T, p.K236X). Six Camarillo White Horses had a missense mutation in exon 12 (c.1805C>T, p.A602V), and five white Thoroughbreds had yet another missense mutation in exon 13 (c.1960G>A, p.G654R). Our results indicate that the dominant white color in Franches-Montagnes Horses is caused by a nonsense mutation in the KIT gene and that multiple independent mutations within this gene appear to be responsible for dominant white in several other modern horse populations.

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.

Using comparative genomics to reorder the human genome sequence into a virtual sheep genome

BACKGROUND

Is it possible to construct an accurate and detailed subgene-level map of a genome using bacterial artificial chromosome (BAC) end sequences, a sparse marker map, and the sequences of other genomes?

RESULTS

A sheep BAC library, CHORI-243, was constructed and the BAC end sequences were determined and mapped with high sensitivity and low specificity onto the frameworks of the human, dog, and cow genomes. To maximize genome coverage, the coordinates of all BAC end sequence hits to the cow and dog genomes were also converted to the equivalent human genome coordinates. The 84,624 sheep BACs (about 5.4-fold genome coverage) with paired ends in the correct orientation (tail-to-tail) and spacing, combined with information from sheep BAC comparative genome contigs (CGCs) built separately on the dog and cow genomes, were used to construct 1,172 sheep BAC-CGCs, covering 91.2% of the human genome. Clustered non-tail-to-tail and outsize BACs located close to the ends of many BAC-CGCs linked BAC-CGCs covering about 70% of the genome to at least one other BAC-CGC on the same chromosome. Using the BAC-CGCs, the intrachromosomal and interchromosomal BAC-CGC linkage information, human/cow and vertebrate synteny, and the sheep marker map, a virtual sheep genome was constructed. To identify BACs potentially located in gaps between BAC-CGCs, an additional set of 55,668 sheep BACs were positioned on the sheep genome with lower confidence. A coordinate conversion process allowed us to transfer human genes and other genome features to the virtual sheep genome to display on a sheep genome browser.

CONCLUSION

We demonstrate that limited sequencing of BACs combined with positioning on a well assembled genome and integrating locations from other less well assembled genomes can yield extensive, detailed subgene-level maps of mammalian genomes, for which genomic resources are currently limited.

A chromosome inversion near the KIT gene and the Tobiano spotting pattern in horses

Tobiano is a white spotting pattern in horses caused by a dominant gene, Tobiano(TO). Here, we report TO associated with a large paracentric chromosome inversion on horse chromosome 3. DNA sequences flanking the inversion were identified and a PCR test was developed to detect the inversion. The inversion was only found in horses with the tobiano pattern, including horses with diverse genetic backgrounds, which indicated a common genetic origin thousands of years ago. The inversion does not interrupt any annotated genes, but begins approximately 100 kb downstream of the KIT gene. This inversion may disrupt regulatory sequences for the KIT gene and cause the white spotting pattern.

Genetic marker anchoring by six-dimensional pools for development of a soybean physical map

BACKGROUND

Integrated genetic and physical maps are extremely valuable for genomic studies and as important references for assembling whole genome shotgun sequences. Screening of a BAC library using molecular markers is an indispensable procedure for integration of both physical and genetic maps of a genome. Molecular markers provide anchor points for integration of genetic and physical maps and also validate BAC contigs assembled based solely on BAC fingerprints. We employed a six-dimensional BAC pooling strategy and an in silico approach to anchor molecular markers onto the soybean physical map.

RESULTS

A total of 1,470 markers (580 SSRs and 890 STSs) were anchored by PCR on a subset of a Williams 82 BstY I BAC library pooled into 208 pools in six dimensions. This resulted in 7,463 clones (approximately 1x genome equivalent) associated with 1470 markers, of which the majority of clones (6,157, 82.5%) were anchored by one marker and 1106 (17.5%) individual clones contained two or more markers. This contributed to 1184 contigs having anchor points through this 6-D pool screening effort. In parallel, the 21,700 soybean Unigene set from NCBI was used to perform in silico mapping on 80,700 Williams 82 BAC end sequences (BES). This in silico analysis yielded 9,835 positive results anchored by 4152 unigenes that contributed to 1305 contigs and 1624 singletons. Among the 1305 contigs, 305 have not been previously anchored by PCR. Therefore, 1489 (78.8%) of 1893 contigs are anchored with molecular markers. These results are being integrated with BAC fingerprints to assemble the BAC contigs. Ultimately, these efforts will lead to an integrated physical and genetic map resource.

CONCLUSION

We demonstrated that the six-dimensional soybean BAC pools can be efficiently used to anchor markers to soybean BACs despite the complexity of the soybean genome. In addition to anchoring markers, the 6-D pooling method was also effective for targeting BAC clones for investigating gene families and duplicated regions in the genome, as well as for extending physical map contigs.

Evolutionary history of chromosome 11 featuring four distinct centromere repositioning events in Catarrhini

Panels of BAC clones used in FISH experiments allow a detailed definition of chromosomal marker arrangement and orientation during evolution. This approach has disclosed the centromere repositioning phenomenon, consisting in the activation of a novel, fully functional centromere in an ectopic location, concomitant with the inactivation of the old centromere. In this study, appropriate panels of BAC clones were used to track the chromosome 11 evolutionary history in primates and nonprimate boreoeutherian mammals. Chromosome 11 synteny was found to be highly conserved in both primate and boreoeutherian mammalian ancestors. Amazingly, we detected four centromere repositioning events in primates (in Old World monkeys, in gibbons, in orangutans, and in the Homo-Pan-Gorilla (H-P-G) clade ancestor), and one in Equidae. Both H-P-G and Lar gibbon novel centromeres were flanked by large duplicons with high sequence similarity. Outgroup species analysis revealed that this duplicon was absent in phylogenetically more distant primates. The chromosome 11 ancestral centromere was probably located near the HSA11q telomere. The domain of this inactivated centromere, in humans, is almost devoid of segmental duplications. An inversion occurred in chromosome 11 in the common ancestor of H-P-G. A large duplicon, again absent in outgroup species, was found located adjacent to the inversion breakpoints. In Hominoidea, almost all the five largest duplicons of this chromosome appeared involved in significant evolutionary architectural changes.

A BAC-based physical map of the channel catfish genome

Catfish is the major aquaculture species in the United States. To enhance its genome studies involving genetic linkage and comparative mapping, a bacterial artificial chromosome (BAC) contig-based physical map of the channel catfish (Ictalurus punctatus) genome was generated using four-color fluorescence-based fingerprints. Fingerprints of 34,580 BAC clones (5.6x genome coverage) were generated for the FPC assembly of the BAC contigs. A total of 3307 contigs were assembled using a cutoff value of 1x10(-20). Each contig contains an average of 9.25 clones with an average size of 292 kb. The combined contig size for all contigs was 0.965 Gb, approximately the genome size of the channel catfish. The reliability of the contig assembly was assessed by both hybridization of gene probes to BAC clones contained in the fingerprinted assembly and validation of randomly selected contigs using overgo probes designed from BAC end sequences. The presented physical map should greatly enhance genome research in the catfish, particularly aiding in the identification of genomic regions containing genes underlying important performance traits.

Fluorescent in situ hybridization mapping of the epidermal growth factor receptor gene in donkey

The physical localization of the epidermal growth factor receptor (EGFR) gene was performed on donkey chromosomes. Bacterial artificial chromosome DNA containing the equine EGFR gene was used to map this gene by fluorescent in situ hybridization on donkey metaphase chromosomes. The gene was mapped on donkey 1q21.1 region.

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.

A radiation hybrid map of sheep chromosome 23 based on ovine BAC-end sequences

More than 375,000 BAC-end sequences (BES) of the CHORI-243 ovine BAC library have been deposited in public databases. blastn searches with these BES against HSA18 revealed 1806 unique and significant hits. We used blastn-anchored BES for an in silico prediction of gene content and chromosome assignment of comparatively mapped ovine BAC clones. Ovine BES were selected at approximately 1.3-Mb intervals of HSA18 and incorporated into a human-sheep comparative map. An ovine 5000-rad whole-genome radiation hybrid panel (USUoRH5000) was typed with 70 markers, all of which mapped to OAR23. The resulting OAR23 RH map included 43 markers derived from BES with high and unique BLAST hits to the sequence of the orthologous HSA18, nine EST-derived markers, 16 microsatellite markers taken from the ovine linkage map and two bovine microsatellite markers. Six new microsatellite markers derived from the 43 mapped BES and the two bovine microsatellite markers were linkage-mapped using the International Mapping Flock (IMF). Thirteen additional microsatellite markers were derived from other ovine BES with high and unique BLAST hits to the sequence of the orthologous HSA18 and also positioned on the ovine linkage map but not incorporated into the OAR23 RH map. This resulted in 24 markers in common and in the same order between the RH and linkage maps. Eight of the BES-derived markers were mapped using fluorescent in situ hybridization (FISH), to thereby align the RH and cytogenetic maps. Comparison of the ovine chromosome 23 RH map with the HSA18 map identified and localized three major breakpoints between HSA18 and OAR23. The positions of these breakpoints were equivalent to those previously shown for syntenic BTA24 and HSA18. This study presents evidence for the usefulness of ovine BES when constructing a high-resolution comprehensive map for a single sheep chromosome. The comparative analysis confirms and refines knowledge about chromosomal conservation and rearrangements between sheep, cattle and human. The constructed RH map demonstrates the resolution and utility of the newly constructed ovine RH panel.

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 scaffold for the Chinese hamster genome

Chinese hamster ovary (CHO) cells are a prevalent tool in biological research and are among the most widely used host cell lines for production of recombinant therapeutic proteins. While research in other organisms has been revolutionized through the development of DNA sequence-based tools, the lack of comparable genomic resources for the Chinese hamster has impeded similar work in CHO cell lines. A comparative genomics approach, based upon the completely sequenced mouse genome, can facilitate genomic work in this important organism. Using chromosome synteny to define regions of conserved linkage between Chinese hamster and mouse chromosomes, a working scaffold for the Chinese hamster genome has been developed. Mapping CHO and Chinese hamster sequences to the mouse genome creates direct access to relevant information in public databases. Additionally, mapping gene expression data onto a chromosome scaffold affords the ability to interpret information in a genomic context, potentially revealing important structural and regulatory features in the Chinese hamster genome. Further development of this genomic scaffold will provide opportunities to use biomolecular tools for research in CHO cell lines today and will be an asset to future efforts to sequence the Chinese hamster genome.