An ordered BAC contig map of the equine major histocompatibility complex

Genetic diversity and signatures of selection in Icelandic horses and Exmoor ponies

BACKGROUND

The Icelandic horse and Exmoor pony are ancient, native breeds, adapted to harsh environmental conditions and they have both undergone severe historic bottlenecks. However, in modern days, the selection pressures on these breeds differ substantially. The aim of this study was to assess genetic diversity in both breeds through expected (HE) and observed heterozygosity (HO) and effective population size (Ne). Furthermore, we aimed to identify runs of homozygosity (ROH) to estimate and compare genomic inbreeding and signatures of selection in the breeds.

RESULTS

HO was estimated at 0.34 and 0.33 in the Icelandic horse and Exmoor pony, respectively, aligning closely with HE of 0.34 for both breeds. Based on genomic data, the Ne for the last generation was calculated to be 125 individuals for Icelandic horses and 42 for Exmoor ponies. Genomic inbreeding coefficient (FROH) ranged from 0.08 to 0.20 for the Icelandic horse and 0.12 to 0.27 for the Exmoor pony, with the majority of inbreeding attributed to short ROHs in both breeds. Several ROH islands associated with performance were identified in the Icelandic horse, featuring target genes such as DMRT3, DOCK8, EDNRB, SLAIN1, and NEURL1. Shared ROH islands between both breeds were linked to metabolic processes (FOXO1), body size, and the immune system (CYRIB), while private ROH islands in Exmoor ponies were associated with coat colours (ASIP, TBX3, OCA2), immune system (LYG1, LYG2), and fertility (TEX14, SPO11, ADAM20).

CONCLUSIONS

Evaluations of genetic diversity and inbreeding reveal insights into the evolutionary trajectories of both breeds, highlighting the consequences of population bottlenecks. While the genetic diversity in the Icelandic horse is acceptable, a critically low genetic diversity was estimated for the Exmoor pony, which requires further validation. Identified signatures of selection highlight the differences in the use of the two breeds as well as their adaptive trait similarities. The results provide insight into genomic regions under selection pressure in a gaited performance horse breed and various adaptive traits in small-sized native horse breeds. This understanding contributes to preserving genetic diversity and population health in these equine populations.

Newly identified variability of the antigen binding site coding sequences of the equine major histocompatibility complex class I and class II genes

The major histocompatibility complex (MHC) with its class I and II genes plays a crucial role in the immune response to pathogens by presenting oligopeptide antigens to various immune response effector cells. In order to counteract the vast variability of infectious agents, MHC class I and II genes usually retain high levels of SNPs mainly concentrated in the exons encoding the antigen binding sites. The aim of the study was to reveal new variability of selected MHC genes with a special focus on MHC class I physical haplotypes. Long-range NGS to was used to identify exon 2-exon 3 alleles in three genetically distinct horse breeds. A total of 116 allelic variants were found in the MHC class I genes Eqca-1, Eqca-2, Eqca-7 and Eqca-Ψ, 112 of which were novel. The MHC class II DRA locus was confirmed to comprise five exon 2 alleles, and no new sequences were observed. Additional variability in terms of 15 novel exon 2 alleles was identified in the DQA1 locus. Extensive overall variability across the entire MHC region was confirmed by an analysis of MHC-linked microsatellite loci. Both diversifying and purifying selection were detected within the MHC class I and II loci analyzed.

Identification of vaccine candidates against rhodococcus equi by combining pangenome analysis with a reverse vaccinology approach

Rhodococcus equi (R. equi) is a zoonotic opportunistic pathogen that can cause life-threatening infections. The rapid evolution of multidrug-resistant R. equi and the fact that there is no currently licensed effective vaccine against R. equi warrant the need for vaccine development. Reverse vaccinology (RV), which involves screening a pathogen's entire genome and proteome using various web-based prediction tools, is considered one of the most effective approaches for identifying vaccine candidates. Here, we performed a pangenome analysis to determine the core proteins of R. equi. We then used the RV approach to examine the subcellular localization, host and gut flora homology, antigenicity, transmembrane helices, physicochemical properties, and immunogenicity of the core proteins to select potential vaccine candidates. The vaccine candidates were then subjected to epitope mapping to predict the exposed antigenic epitopes that possess the ability to bind with major histocompatibility complex I/II (MHC I/II) molecules. These vaccine candidates and epitopes will form a library of elements for the development of a polyvalent or universal vaccine against R. equi. Sixteen R. equi complete proteomes were found to contain 6,238 protein families, and the core proteins consisted of 3,969 protein families (∼63.63% of the pangenome), reflecting a low degree of intraspecies genomic variability. From the pool of core proteins, 483 nonhost homologous membrane and extracellular proteins were screened, and 12 vaccine candidates were finally identified according to their antigenicity, physicochemical properties and other factors. These included four cell wall/membrane/envelope biogenesis proteins; four amino acid transport and metabolism proteins; one cell cycle control, cell division and chromosome partitioning protein; one carbohydrate transport and metabolism protein; one secondary metabolite biosynthesis, transport and catabolism protein; and one defense mechanism protein. All 12 vaccine candidates have an experimentally validated 3D structure available in the protein data bank (PDB). Epitope mapping of the candidates showed that 16 MHC I epitopes and 13 MHC II epitopes with the strongest immunogenicity were exposed on the protein surface, indicating that they could be used to develop a polypeptide vaccine. Thus, we utilized an analytical strategy that combines pangenome analysis and RV to generate a peptide antigen library that simplifies the development of multivalent or universal vaccines against R. equi and can be applied to the development of other vaccines.

High-Resolution Genotyping of Expressed Equine MHC Reveals a Highly Complex MHC Structure

The Major Histocompatibility Complex (MHC) genes play a key role in a number of biological processes, most notably in immunological responses. The MHCI and MHCII genes incorporate a complex set of highly polymorphic and polygenic series of genes, which, due to the technical limitations of previously available technologies, have only been partially characterized in non-model but economically important species such as the horse. The advent of high-throughput sequencing platforms has provided new opportunities to develop methods to generate high-resolution sequencing data on a large scale and apply them to the analysis of complex gene sets such as the MHC. In this study, we developed and applied a MiSeq-based approach for the combined analysis of the expressed MHCI and MHCII repertoires in cohorts of Thoroughbred, Icelandic, and Norwegian Fjord Horses. The approach enabled us to generate comprehensive MHCI/II data for all of the individuals (n = 168) included in the study, identifying 152 and 117 novel MHCI and MHCII sequences, respectively. There was limited overlap in MHCI and MHCII haplotypes between the Thoroughbred and the Icelandic/Norwegian Fjord horses, showcasing the variation in MHC repertoire between genetically divergent breeds, and it can be inferred that there is much more MHC diversity in the global horse population. This study provided novel insights into the structure of the expressed equine MHC repertoire and highlighted unique features of the MHC in horses.

Genetic diversity, evolution and selection in the major histocompatibility complex DRB and DQB loci in the family Equidae

Background

The mammalian Major Histocompatibility Complex (MHC) is a genetic region containing highly polymorphic genes with immunological functions. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. The MHC class II sub-region contains genes expressed in antigen presenting cells. The antigen binding site is encoded by the second exon of genes encoding antigen presenting molecules. The exon 2 sequences of these MHC genes have evolved under the selective pressure of pathogens. Interspecific differences can be observed in the class II sub-region. The family Equidae includes a variety of domesticated, and free-ranging species inhabiting a range of habitats exposed to different pathogens and represents a model for studying this important part of the immunogenome. While equine MHC class II DRA and DQA loci have received attention, the genetic diversity and effects of selection on DRB and DQB loci have been largely overlooked. This study aimed to provide the first in-depth analysis of the MHC class II DRB and DQB loci in the Equidae family.

Results

Three DRB and two DQB genes were identified in the genomes of all equids. The genes DRB2, DRB3 and DQB3 showed high sequence conservation, while polymorphisms were more frequent at DRB1 and DQB1 across all species analyzed. DQB2 was not found in the genome of the Asiatic asses Equus hemionus kulan and E. h. onager. The bioinformatic analysis of non-zero-coverage-bases of DRB and DQB genes in 14 equine individual genomes revealed differences among individual genes. Evidence for recombination was found for DRB1, DRB2, DQB1 and DQB2 genes. Trans-species allele sharing was identified in all genes except DRB1. Site-specific selection analysis predicted genes evolving under positive selection both at DRB and DQB loci. No selected amino acid sites were identified in DQB3.

Conclusions

The organization of the MHC class II sub-region of equids is similar across all species of the family. Genomic sequences, along with phylogenetic trees suggesting effects of selection as well as trans-species polymorphism support the contention that pathogen-driven positive selection has shaped the MHC class II DRB/DQB sub-regions in the Equidae.

Origin and evolution of the major histocompatibility complex class I region in eutherian mammals

Major histocompatibility complex (MHC) genes in vertebrates are vital in defending against pathogenic infections. To gain new insights into the evolution of MHC Class I (MHCI) genes and test competing hypotheses on the origin of the MHCI region in eutherian mammals, we studied available genome assemblies of nine species in Afrotheria, Xenarthra, and Laurasiatheria, and successfully characterized the MHCI region in six species. The following numbers of putatively functional genes were detected: in the elephant, four, one, and eight in the extended class I region, and κ and β duplication blocks, respectively; in the tenrec, one in the κ duplication block; and in the four bat species, one or two in the β duplication block. Our results indicate that MHCI genes in the κ and β duplication blocks may have originated in the common ancestor of eutherian mammals. In the elephant, tenrec, and all four bats, some MHCI genes occurred outside the MHCI region, suggesting that eutherians may have a more complex MHCI genomic organization than previously thought. Bat-specific three- or five-amino-acid insertions were detected in the MHCI α1 domain in all four bats studied, suggesting that pathogen defense in bats relies on MHCIs having a wider peptide-binding groove, as previously assayed by a bat MHCI gene with a three-amino-acid insertion showing a larger peptide repertoire than in other mammals. Our study adds to knowledge on the diversity of eutherian MHCI genes, which may have been shaped in a taxon-specific manner.

A high-density BAC physical map covering the entire MHC region of addax antelope genome

BACKGROUND

The mammalian major histocompatibility complex (MHC) harbours clusters of genes associated with the immunological defence of animals against infectious pathogens. At present, no complete MHC physical map is available for any of the wild ruminant species in the world.

RESULTS

The high-density physical map is composed of two contigs of 47 overlapping bacterial artificial chromosome (BAC) clones, with an average of 115 Kb for each BAC, covering the entire addax MHC genome. The first contig has 40 overlapping BAC clones covering an approximately 2.9 Mb region of MHC class I, class III, and class IIa, and the second contig has 7 BAC clones covering an approximately 500 Kb genomic region that harbours MHC class IIb. The relative position of each BAC corresponding to the MHC sequence was determined by comparative mapping using PCR screening of the BAC library of 192,000 clones, and the order of BACs was determined by DNA fingerprinting. The overlaps of neighboring BACs were cross-verified by both BAC-end sequencing and co-amplification of identical PCR fragments within the overlapped region, with their identities further confirmed by DNA sequencing.

CONCLUSIONS

We report here the successful construction of a high-quality physical map for the addax MHC region using BACs and comparative mapping. The addax MHC physical map we constructed showed one gap of approximately 18 Mb formed by an ancient autosomal inversion that divided the MHC class II into IIa and IIb. The autosomal inversion provides compelling evidence that the MHC organizations in all of the ruminant species are relatively conserved.

MHC haplotype diversity in Icelandic horses determined by polymorphic microsatellites

The Icelandic horse has been maintained as a closed population in its eponymous homeland for many generations, with no recorded introductions of new horses of any breed since the year 1000 CE. Here we determined the diversity of major histocompatibility complex (MHC) haplotypes in 156 Icelandic horses from two groups, based on a panel of 12 polymorphic intra-MHC microsatellites tested in families of various composition. We identified a total of 79 MHC haplotypes in these two groups, including one documented intra-MHC recombination event from a total of 147 observed meioses. None of these MHC haplotypes have been previously described in any other horse breed. Only one MHC homozygote was found in the entire population studied. These results indicate a very high level of MHC heterozygosity and haplotype diversity in the Icelandic horse. The environment in Iceland is remarkable for its lack of common agents of equine infectious disease, including equine herpesvirus type 1, influenza virus, and streptococcus equi. The driving forces for maintenance of MHC heterozygosity in Icelandic horses must thus be sought outside of these major horse pathogens. Based on our results, we propose that intra-MHC recombination may play a major role in the generation of novel haplotypes.

Copy number variations in Friesian horses and genetic risk factors for insect bite hypersensitivity

Background

Many common and relevant diseases affecting equine welfare have yet to be tested regarding structural variants such as copy number variations (CNVs). CNVs make up a substantial proportion of total genetic variability in populations of many species, resulting in more sequence differences between individuals than SNPs. Associations between CNVs and disease phenotypes have been established in several species, but equine CNV studies have been limited. Aim of this study was to identify CNVs and to perform a genome-wide association (GWA) study in Friesian horses to identify genomic loci associated with insect bite hypersensitivity (IBH), a common seasonal allergic dermatitis observed in many horse breeds worldwide.

Results

Genotypes were obtained using the Axiom® Equine Genotyping Array containing 670,796 SNPs. After quality control of genotypes, 15,041 CNVs and 5350 CNV regions (CNVRs) were identified in 222 Friesian horses. Coverage of the total genome by CNVRs was 11.2% with 49.2% of CNVRs containing genes. 58.0% of CNVRs were novel (i.e. so far only identified in Friesian horses). A SNP- and CNV-based GWA analysis was performed, where about half of the horses were affected by IBH. The SNP-based analysis showed a highly significant association between the MHC region on ECA20 and IBH in Friesian horses. Associations between the MHC region on ECA20 and IBH were also detected based on the CNV-based analysis. However, CNVs associated with IBH in Friesian horses were not often in close proximity to SNPs identified to be associated with IBH.

Conclusions

CNVs were identified in a large sample of the Friesian horse population, thereby contributing to our knowledge on CNVs in horses and facilitating our understanding of the equine genome and its phenotypic expression. A clear association was identified between the MHC region on ECA20 and IBH in Friesian horses based on both SNP- and CNV-based GWA studies. These results imply that MHC contributes to IBH sensitivity in Friesian horses. Although subsequent analyses are needed for verification, nucleotide differences, as well as more complex structural variations like CNVs, seem to contribute to IBH sensitivity. IBH should be considered as a common disease with a complex genomic architecture.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0657-0) contains supplementary material, which is available to authorized users.

Microsatellite markers for evaluating the diversity of the natural killer complex and major histocompatibility complex genomic regions in domestic horses

Genotyping microsatellite markers represents a standard, relatively easy, and inexpensive method of assessing genetic diversity of complex genomic regions in various animal species, such as the major histocompatibility complex (MHC) and/or natural killer cell receptor (NKR) genes. MHC-linked microsatellite markers have been identified and some of them were used for characterizing MHC polymorphism in various species, including horses. However, most of those were MHC class II markers, while MHC class I and III sub-regions were less well covered. No tools for studying genetic diversity of NKR complex genomic regions are available in horses. Therefore, the aims of this work were to establish a panel of markers suitable for analyzing genetic diversity of the natural killer complex (NKC), and to develop additional microsatellite markers of the MHC class I and class III genomic sub-regions in horses. Nine polymorphic microsatellite loci were newly identified in the equine NKC. Along with two previously reported microsatellites flanking this region, they constituted a panel of 11 loci allowing to characterize genetic variation in this functionally important part of the horse genome. Four newly described MHC class I/III-linked markers were added to 11 known microsatellites to establish a panel of 15 MHC markers with a better coverage of the class I and class III sub-regions. Major characteristics of the two panels produced on a group of 65 horses of 13 breeds and on five Przewalski's horses showed that they do reflect genetic variation within the horse species.

Insights into the ancestral organisation of the mammalian MHC class II region from the genome of the pteropid bat, Pteropus alecto

BACKGROUND

Bats are an extremely successful group of mammals and possess a variety of unique characteristics, including their ability to co-exist with a diverse range of pathogens. The major histocompatibility complex (MHC) is the most gene dense and polymorphic region of the genome and MHC class II (MHC-II) molecules play a vital role in the presentation of antigens derived from extracellular pathogens and activation of the adaptive immune response. Characterisation of the MHC-II region of bats is crucial for understanding the evolution of the MHC and of the role of pathogens in shaping the immune system.

RESULTS

Here we describe the relatively contracted MHC-II region of the Australian black flying-fox (Pteropus alecto), providing the first detailed insight into the MHC-II region of any species of bat. Twelve MHC-II genes, including one locus (DRB2) located outside the class II region, were identified on a single scaffold in the bat genome. The presence of a class II locus outside the MHC-II region is atypical and provides evidence for an ancient class II duplication block. Two non-classical loci, DO and DM and two classical, DQ and DR loci, were identified in P. alecto. A putative classical, DPB pseudogene was also identified. The bat's antigen processing cluster, though contracted, remains highly conserved, thus supporting its importance in antigen presentation and disease resistance.

CONCLUSIONS

This detailed characterisation of the bat MHC-II region helps to fill a phylogenetic gap in the evolution of the mammalian class II region and is a stepping stone towards better understanding of the immune responses in bats to viral, bacterial, fungal and parasitic infections.

Genomic structure of the horse major histocompatibility complex class II region resolved using PacBio long-read sequencing technology

The mammalian Major Histocompatibility Complex (MHC) region contains several gene families characterized by highly polymorphic loci with extensive nucleotide diversity, copy number variation of paralogous genes, and long repetitive sequences. This structural complexity has made it difficult to construct a reliable reference sequence of the horse MHC region. In this study, we used long-read single molecule, real-time (SMRT) sequencing technology from Pacific Biosciences (PacBio) to sequence eight Bacterial Artificial Chromosome (BAC) clones spanning the horse MHC class II region. The final assembly resulted in a 1,165,328 bp continuous gap free sequence with 35 manually curated genomic loci of which 23 were considered to be functional and 12 to be pseudogenes. In comparison to the MHC class II region in other mammals, the corresponding region in horse shows extraordinary copy number variation and different relative location and directionality of the Eqca-DRB, -DQA, -DQB and -DOB loci. This is the first long-read sequence assembly of the horse MHC class II region with rigorous manual gene annotation, and it will serve as an important resource for association studies of immune-mediated equine diseases and for evolutionary analysis of genetic diversity in this region.

Polymorphism at expressed DQ and DR loci in five common equine MHC haplotypes

The polymorphism of major histocompatibility complex (MHC) class II DQ and DR genes in five common equine leukocyte antigen (ELA) haplotypes was determined through sequencing of mRNA transcripts isolated from lymphocytes of eight ELA homozygous horses. Ten expressed MHC class II genes were detected in horses of the ELA-A3 haplotype carried by the donor horses of the equine bacterial artificial chromosome (BAC) library and the reference genome sequence: four DR genes and six DQ genes. The other four ELA haplotypes contained at least eight expressed polymorphic MHC class II loci. Next generation sequencing (NGS) of genomic DNA of these four MHC haplotypes revealed stop codons in the DQA3 gene in the ELA-A2, ELA-A5, and ELA-A9 haplotypes. Few NGS reads were obtained for the other MHC class II genes that were not amplified in these horses. The amino acid sequences across haplotypes contained locus-specific residues, and the locus clusters produced by phylogenetic analysis were well supported. The MHC class II alleles within the five tested haplotypes were largely non-overlapping between haplotypes. The complement of equine MHC class II DQ and DR genes appears to be well conserved between haplotypes, in contrast to the recently described variation in class I gene loci between equine MHC haplotypes. The identification of allelic series of equine MHC class II loci will aid comparative studies of mammalian MHC conservation and evolution and may also help to interpret associations between the equine MHC class II region and diseases of the horse.

Inability of FMDV replication in equine kidney epithelial cells is independent of integrin αvβ3 and αvβ6

Integrins can function as receptors for foot-and-mouth disease virus (FMDV) in epithelium. Horses are believed to be insusceptible to this disease, but the mechanism of resistance remains unclear. To detect whether FMDV can use integrin to attach to equine epithelial, we compared the utilities of αvβ3 and αvβ6 between bovine and equine kidney epithelial cells (KECs). Equine KECs showed almost equal efficiency to those of bovine. Further, the integrin αv, β3, and β6 subunits from bovine and equine were cloned and vectors were transfected into SW480 cells and COS-1 cells alone or together, and virus titers were used to determine the viral replication. In all cases, the virus reproduced successfully. Overall, FMDV can replicate in SW480 cells transfected with equine β3/β6 subunits and COS-1 cells transfected with equine αvβ3/αvβ6 integrins, but not in EKECs. These results indicated that failure of FMDV replication in EKECs was not attributed to integrin receptors.

Evolution and comparative analysis of the bat MHC-I region

Bats are natural hosts to numerous viruses and have ancient origins, having diverged from other eutherian mammals early in evolution. These characteristics place them in an important position to provide insights into the evolution of the mammalian immune system and antiviral immunity. We describe the first detailed partial map of a bat (Pteropus alecto) MHC-I region with comparative analysis of the MHC-I region and genes. The bat MHC-I region is highly condensed, yet relatively conserved in organisation, and is unusual in that MHC-I genes are present within only one of the three highly conserved class I duplication blocks. We hypothesise that MHC-I genes first originated in the β duplication block, and subsequently duplicated in a step-wise manner across the MHC-I region during mammalian evolution. Furthermore, bat MHC-I genes contain unique insertions within their peptide-binding grooves potentially affecting the peptide repertoire presented to T cells, which may have implications for the ability of bats to control infection without overt disease.

A comprehensive mapping of the structure and gene organisation in the sheep MHC class I region

Background

The major histocompatibility complex (MHC) is a chromosomal region that regulates immune responsiveness in vertebrates. This region is one of the most important for disease resistance because it has been associated with resistance or susceptibility to a wide variety of diseases and because the MHC often accounts for more of the variance than other loci. Selective breeding for disease resistance is becoming increasingly common in livestock industries, and it is important to determine how this will influence MHC polymorphism and resistance to diseases that are not targeted for selection. However, in sheep the order and sequence of the protein coding genes is controversial. Yet this information is needed to determine precisely how the MHC influences resistance and susceptibility to disease.

Methods

CHORI bacterial artificial chromosomes (BACs) known to contain sequences from the sheep MHC class I region were sub-cloned, and the clones partially sequenced. The resulting sequences were analysed and re-assembled to identify gene content and organisation within each BAC. The low resolution MHC class I physical map was then compared to the cattle reference genome, the Chinese Merino sheep MHC map published by Gao, et al. (2010) and the recently available sheep reference genome.

Results

Immune related class I genes are clustered into 3 blocks; beta, kappa and a novel block not previously identified in other organisms. The revised map is more similar to Bovidae maps than the previous sheep maps and also includes several genes previously not annotated in the Chinese Merino BAC assembly and others not currently annotated in the sheep reference chromosome 20. In particular, the organisation of nonclassical MHC class I genes is similar to that present in the cattle MHC. Sequence analysis and prediction of amino acid sequences of MHC class I classical and nonclassical genes was performed and it was observed that the map contained one classical and eight nonclassical genes together with three possible pseudogenes.

Conclusions

The comprehensive physical map of the sheep MHC class I region enhances our understanding of the genetic architecture of the class I MHC region in sheep and will facilitate future studies of MHC function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1992-4) contains supplementary material, which is available to authorized users.

The common equine class I molecule Eqca-1*00101 (ELA-A3.1) is characterized by narrow peptide binding and T cell epitope repertoires

Here we describe a detailed quantitative peptide-binding motif for the common equine leukocyte antigen (ELA) class I allele Eqca-1*00101, present in roughly 25 % of Thoroughbred horses. We determined a preliminary binding motif by sequencing endogenously bound ligands. Subsequently, a positional scanning combinatorial library (PSCL) was used to further characterize binding specificity and derive a quantitative motif involving aspartic acid in position 2 and hydrophobic residues at the C-terminus. Using this motif, we selected and tested 9- and 10-mer peptides derived from the equine herpesvirus type 1 (EHV-1) proteome for their capacity to bind Eqca-1*00101. PSCL predictions were very efficient, with an receiver operating characteristic (ROC) curve performance of 0.877, and 87 peptides derived from 40 different EHV-1 proteins were identified with affinities of 500 nM or higher. Quantitative analysis revealed that Eqca-1*00101 has a narrow peptide-binding repertoire, in comparison to those of most human, non-human primate, and mouse class I alleles. Peripheral blood mononuclear cells from six EHV-1-infected, or vaccinated but uninfected, Eqca-1*00101-positive horses were used in IFN-γ enzyme-linked immunospot (ELISPOT) assays. When we screened the 87 Eqca-1*00101-binding peptides for T cell reactivity, only one Eqca-1*00101 epitope, derived from the intermediate-early protein ICP4, was identified. Thus, despite its common occurrence in several horse breeds, Eqca-1*00101 is associated with a narrow binding repertoire and a similarly narrow T cell response to an important equine viral pathogen. Intriguingly, these features are shared with other human and macaque major histocompatibility complex (MHC) molecules with a similar specificity for D in position 2 or 3 in their main anchor motif.

Characterisation of major histocompatibility complex class I genes at the fetal-maternal interface of marsupials

Major histocompatibility complex class I molecules (MHC-I) are expressed at the cell surface and are responsible for the presentation of self and non-self antigen repertoires to the immune system. Eutherian mammals express both classical and non-classical MHC-I molecules in the placenta, the latter of which are thought to modulate the maternal immune response during pregnancy. Marsupials last shared a common ancestor with eutherian mammals such as humans and mice over 160 million years ago. Since, like eutherians, they have an intra-uterine development dependent on a placenta, albeit a short-lived and less invasive one, they provide an opportunity to investigate the evolution of MHC-I expression at the fetal-maternal interface. We have characterised MHC-I mRNA expression in reproductive tissues of the tammar wallaby (Macropus eugenii) from the time of placental attachment to day 25 of the 26.5 day pregnancy. Putative classical MHC-I genes were expressed in the choriovitelline placenta, fetus, and gravid endometrium throughout the whole of this period. The MHC-I classical sequences were phylogenetically most similar to the Maeu-UC (50/100 clones) and Maeu-UA genes (7/100 clones). Expression of three non-classical MHC-I genes (Maeu-UD, Maeu-UK and Maeu-UM) were also present in placental samples. The results suggest that expression of classical and non-classical MHC-I genes in extant marsupial and eutherian mammals may have been necessary for the evolution of the ancestral therian placenta and survival of the mammalian fetus at the maternal-fetal interface.

Identification of equine major histocompatibility complex haplotypes using polymorphic microsatellites

A system for identifying equine major histocompatibility complex (MHC) haplotypes was developed based on five polymorphic microsatellites located within the MHC region on ECA 20. Molecular signatures for 50 microsatellite haplotypes were recognized from typing 353 horses. Of these, 23 microsatellite haplotypes were associated with 12 established equine leucocyte antigen (ELA) haplotypes in Thoroughbreds and Standardbreds. Five ELA serotypes were associated with multiple microsatellite subhaplotypes, expanding the estimates of diversity in the equine MHC. The strong correlations between serological and microsatellite typing demonstrated a linkage to known MHC class I protein polymorphisms and validated this assay as a useful supplement to ELA serotyping, and in some applications, a feasible alternative method for MHC genotyping in horse families and in population studies.

A cytogenetic and comparative map of camelid chromosome 36 and the minute in alpacas

Recent advances in camelid genomics have provided draft sequence assemblies and the first comparative and gene maps for the dromedary (CDR) and the alpaca (LPA). However, no map information is currently available for the smallest camelid autosome-chr36. The chromosome is also of clinical interest because of its involvement in the minute chromosome syndrome (MCS) in infertile alpacas. Here, we developed molecular markers for camelid chr36 by direct sequencing CDR36 and LPA minute and by bioinformatics analysis of alpaca unplaced sequence scaffolds. We constructed a cytogenetic map for chr36 in the alpaca, llama, and dromedary and showed its homology to human chromosome 7 (HSA7) at 49.8-55.5 Mb. The chr36 map comprised seven markers, including two genes-ZPBP and WVC2. Comparative status of HSA7 was further refined by cytogenetic mapping of 16 HSA7 orthologs in camelid chromosomes 7 and 18 and by the analysis of HSA7-conserved synteny blocks across 11 vertebrate species. Finally, mapping chr36 markers in infertile alpacas confirmed that the minute chromosome was a derivative of chr36, but the small size was not a result of a large deletion or a translocation. Instead, cytogenetic mapping of 5.8S, 18S, and 28S rRNA genes (nucleolus organizer region (NOR)) revealed that the size difference between chr36 homologs in infertile alpacas was due to a heterozygous presence of NOR, whereas chr36 in fertile alpacas had no NOR. We theorized that the heterozygous NOR might affect chr36 pairing, recombination, and segregation in meiosis and, thus fertility.

The equine immune responses to infectious and allergic disease: a model for humans?

The modern horse, Equus caballus has historically made important contributions to the field of immunology, dating back to Emil von Behring's description of curative antibodies in equine serum over a century ago. While the horse continues to play an important role in human serotherapy, the mouse has replaced the horse as the predominant experimental animal in immunology research. Nevertheless, continuing efforts have led to an improved understanding of the equine immune response in a variety of infectious and non-infectious diseases. Based on this information, we can begin to identify specific situations where the horse may provide a unique immunological model for certain human diseases.

Immediate-early protein of equid herpesvirus type 1 as a target for cytotoxic T-lymphocytes in the Thoroughbred horse

Cytotoxic T-lymphocytes (CTLs) are associated with protective immunity against disease caused by equid herpesvirus type 1 (EHV-1). However, the EHV-1 target proteins for CTLs are poorly defined. This limits the development of vaccine candidates designed to stimulate strong CTL immunity. Here, classical CTL assays using lymphocytes from horses of three defined MHC class I types that experienced natural infection with EHV-1 and a modified vaccinia virus construct containing an EHV-1 gene encoding the immediate-early (IE) protein are reported. Horses homozygous for the equine leukocyte antigen (ELA)-A2 haplotype, but not the ELA-A5 haplotype, produced MHC-restricted CTL responses against the IE protein. Previously, horses homozygous for the ELA-A3 haplotype also mounted CTL responses against the IE protein. Both haplotypes are common in major horse breeds, including the Thoroughbred. Thus, the IE protein is an attractive candidate molecule for future studies of T-cell immunity to EHV-1 in the horse.

Genetic risk factors for insidious equine recurrent uveitis in Appaloosa horses

Appaloosa horses are predisposed to equine recurrent uveitis (ERU), an immune-mediated disease characterized by recurring inflammation of the uveal tract in the eye, which is the leading cause of blindness in horses. Nine genetic markers from the ECA1 region responsible for the spotted coat color of Appaloosa horses, and 13 microsatellites spanning the equine major histocompatibility complex (ELA) on ECA20, were evaluated for association with ERU in a group of 53 Appaloosa ERU cases and 43 healthy Appaloosa controls. Three markers were significantly associated (corrected P-value <0.05): a SNP within intron 11 of the TRPM1 gene on ECA1, an ELA class I microsatellite located near the boundary of the ELA class III and class II regions and an ELA class II microsatellite located in intron 1 of the DRA gene. Association between these three genetic markers and the ERU phenotype was confirmed in a second population of 24 insidious ERU Appaloosa cases and 16 Appaloosa controls. The relative odds of being an ERU case for each allele of these three markers were estimated by fitting a logistic mixed model with each of the associated markers independently and with all three markers simultaneously. The risk model using these markers classified ~80% of ERU cases and 75% of controls in the second population as moderate or high risk, and low risk respectively. Future studies to refine the associations at ECA1 and ELA loci and identify functional variants could uncover alleles conferring susceptibility to ERU in Appaloosa horses.

Comparative organization and gene expression profiles of the porcine pseudoautosomal region

The pseudoautosomal region (PAR) has important biological functions in spermatogenesis, male fertility and early development. Even though pig (Sus scrofa, SSC) is an agriculturally and biomedically important species, and its genome is sequenced, current knowledge about the porcine PAR is sparse. Here we defined the PAR in SSCXp/Yp by demarcating the sequence of the pseudoautosomal boundary at X:6,743,567 bp in intron 3-4 of SHROOM2 and showed that SHROOM2 is truncated in SSCY. Cytogenetic mapping of 20 BAC clones containing 15 PAR and X-specific genes revealed that the pig PAR is largely collinear with other mammalian PARs or Xp terminal regions. The results improved the current SSCX sequence assembly and facilitated distinction between the PAR and X-specific genes to study their expression in adult and embryonic tissues. A pilot analysis showed that the PAR genes are expressed at higher levels than X-specific genes during early development, whereas the expression of PAR genes was higher at day 60 compared to day 26, and higher in embryonic tissues compared to placenta. The findings advance the knowledge about the comparative organization of the PAR in mammals and suggest that the region might have important functions in early development in pigs.

A high-resolution radiation hybrid map of the river buffalo major histocompatibility complex and comparison with BoLA

The major histocompatibility complex (MHC) in mammals codes for antigen-presenting proteins. For this reason, the MHC is of great importance for immune function and animal health. Previous studies revealed this gene-dense and polymorphic region in river buffalo to be on the short arm of chromosome 2, which is homologous to cattle chromosome 23. Using cattle-derived STS markers and a river buffalo radiation hybrid (RH) panel (BBURH5000 ), we generated a high-resolution RH map of the river buffalo MHC region. The buffalo MHC RH map (cR5000 ) was aligned with the cattle MHC RH map (cR12000 ) to compare gene order. The buffalo MHC had similar organization to the cattle MHC, with class II genes distributed in two segments, class IIa and class IIb. Class IIa was closely associated with the class I and class III regions, and class IIb was a separate cluster. A total of 53 markers were distributed into two linkage groups based on a two-point LOD score threshold of ≥8. The first linkage group included 32 markers from class IIa, class I and class III. The second linkage group included 21 markers from class IIb. Bacterial artificial chromosome clones for seven loci were mapped by fluorescence in situ hybridization on metaphase chromosomes using single- and double-color hybridizations. The order of cytogenetically mapped markers in the region corroborated the physical order of markers obtained from the RH map and served as anchor points to align and orient the linkage groups.

Development and application of camelid molecular cytogenetic tools

Cytogenetic chromosome maps offer molecular tools for genome analysis and clinical cytogenetics and are of particular importance for species with difficult karyotypes, such as camelids (2n = 74). Building on the available human-camel zoo-fluorescence in situ hybridization (FISH) data, we developed the first cytogenetic map for the alpaca (Lama pacos, LPA) genome by isolating and identifying 151 alpaca bacterial artificial chromosome (BAC) clones corresponding to 44 specific genes. The genes were mapped by FISH to 31 alpaca autosomes and the sex chromosomes; 11 chromosomes had 2 markers, which were ordered by dual-color FISH. The STS gene mapped to Xpter/Ypter, demarcating the pseudoautosomal region, whereas no markers were assigned to chromosomes 14, 21, 22, 28, and 36. The chromosome-specific markers were applied in clinical cytogenetics to identify LPA20, the major histocompatibility complex (MHC)-carrying chromosome, as a part of an autosomal translocation in a sterile male llama (Lama glama, LGL; 2n = 73,XY). FISH with LPAX BACs and LPA36 paints, as well as comparative genomic hybridization, were also used to investigate the origin of the minute chromosome, an abnormally small LPA36 in infertile female alpacas. This collection of cytogenetically mapped markers represents a new tool for camelid clinical cytogenetics and has applications for the improvement of the alpaca genome map and sequence assembly.

Adaptive molecular evolution of the Major Histocompatibility Complex genes, DRA and DQA, in the genus Equus

BACKGROUND

Major Histocompatibility Complex (MHC) genes are central to vertebrate immune response and are believed to be under balancing selection by pathogens. This hypothesis has been supported by observations of extremely high polymorphism, elevated nonsynonymous to synonymous base pair substitution rates and trans-species polymorphisms at these loci. In equids, the organization and variability of this gene family has been described, however the full extent of diversity and selection is unknown. As selection is not expected to act uniformly on a functional gene, maximum likelihood codon-based models of selection that allow heterogeneity in selection across codon positions can be valuable for examining MHC gene evolution and the molecular basis for species adaptations.

RESULTS

We investigated the evolution of two class II MHC genes of the Equine Lymphocyte Antigen (ELA), DRA and DQA, in the genus Equus with the addition of novel alleles identified in plains zebra (E. quagga, formerly E. burchelli). We found that both genes exhibited a high degree of polymorphism and inter-specific sharing of allele lineages. To our knowledge, DRA allelic diversity was discovered to be higher than has ever been observed in vertebrates. Evidence was also found to support a duplication of the DQA locus. Selection analyses, evaluated in terms of relative rates of nonsynonymous to synonymous mutations (dN/dS) averaged over the gene region, indicated that the majority of codon sites were conserved and under purifying selection (dN <dS). However, the most likely evolutionary codon models allowed for variable rates of selection across codon sites at both loci and, at the DQA, supported the hypothesis of positive selection acting on specific sites.

CONCLUSIONS

Observations of elevated genetic diversity and trans-species polymorphisms supported the conclusion that balancing selection may be acting on these loci. Furthermore, at the DQA, positive selection was occurring at antigen binding sites, suggesting that a few selected residues may play a significant role in equid immune function. Future studies in natural equid populations will be valuable for understanding the functional significance of the uniquely diverse DRA locus and for elucidating the mechanism maintaining diversity at these MHC loci.

Major histocompatibility complex diversity in the endangered Ethiopian wolf (Canis simensis)

The major histocompatibility complex (MHC) influences immune response to infection and vaccination. In most species, MHC genes are highly polymorphic, but few wild canid populations have been investigated. In Ethiopian wolves, we identified four DLA (dog leucocyte antigen)-DRB1, two DLA-DQA1 and five DQB1 alleles. Ethiopian wolves, the world's rarest canids with fewer than 500 animals worldwide, are further endangered and threatened by rabies. Major rabies outbreaks in the Bale Mountains of southern Ethiopia (where over half of the Ethiopian wolf population is located) have killed over 75% of wolves in the affected sub-populations. In 2004, following a rabies outbreak, 77 wolves were vaccinated, and 19 were subsequently recaptured to monitor the effectiveness of the intervention. Pre- and post-vaccination rabies antibody titres were available for 18 animals, and all of the animals sero-converted after vaccination. We compared the haplotype frequencies of this group of 18 with the post-vaccination antibody titre, and showed that one haplotype was associated with a lower response (uncorrected P < 0.03). In general, Ethiopian wolves probably have an adequate amount of MHC variation to ensure the survival of the species. However, we sampled only the largest Ethiopian wolf population in Bale, and did not take the smaller populations further north into consideration.

Genetic diversity of the class II major histocompatibility DRA locus in European, Asiatic and African domestic donkeys

The major histocompatibility complex (MHC) genes coding for antigen presenting molecules are the most polymorphic genes in vertebrate genome. The MHC class II DRA gene shows only small variation in many mammalian species, but it exhibits relatively high level of polymorphism in Equidae, especially in donkeys. This extraordinary degree of polymorphism together with signatures of selection in specific amino acids sites makes the donkey DRA gene a suitable model for population diversity studies. The objective of this study was to investigate the DRA gene diversity in three different populations of donkeys under infectious pressure of protozoan parasites, Theileria equi and Babesia caballi. Three populations of domestic donkeys from Italy (N = 68), Jordan (N = 43), and Kenya (N = 78) were studied. A method of the donkey MHC DRA genotyping based on PCR-RFLP and sequencing was designed. In addition to the DRA gene, 12 polymorphic microsatellite loci were genotyped. The presence of Theileria equi and Babesia caballi parasites in peripheral blood was investigated by PCR. Allele and genotype frequencies, observed and expected heterozygosities and F(IS) values were computed as parameters of genetic diversity for all loci genotyped. Genetic distances between the three populations were estimated based on F(ST) values. Statistical associations between parasite infection and genetic polymorphisms were sought. Extensive DRA locus variation characteristic for Equids was found. The results showed differences between populations both in terms of numbers of alleles and their frequencies as well as variation in expected heterozygosity values. Based on comparisons with neutral microsatellite loci, population sub-structure characteristics and association analysis, convincing evidence of pathogen-driven selection at the population level was not provided. It seems that genetic diversity observed in the three populations reflects mostly effects of selective breeding and their different genetic origins.

Equine major histocompatibility complex class I molecules act as entry receptors that bind to equine herpesvirus-1 glycoprotein D

The endotheliotropism of equine herpesvirus-1 (EHV-1) leads to encephalomyelitis secondary to vasculitis and thrombosis in the infected horse central nervous system (CNS). To identify the host factors involved in EHV-1 infection of CNS endothelial cells, we performed functional cloning using an equine brain microvascular endothelial cell cDNA library. Exogenous expression of equine major histocompatibility complex (MHC) class I heavy chain genes conferred susceptibility to EHV-1 infection in mouse NIH3T3 cells, which are not naturally susceptible to EHV-1 infection. Equine MHC class I molecules bound to EHV-1 glycoprotein D (gD), and both anti-gD antibodies and a soluble form of gD blocked viral entry into NIH3T3 cells stably expressing the equine MHC class I heavy chain gene (3T3-A68 cells). Treatment with an anti-equine MHC class I monoclonal antibody blocked EHV-1 entry into 3T3-A68 cells, equine dermis (E. Derm) cells and equine brain microvascular endothelial cells. In addition, inhibition of cell surface expression of MHC class I molecules in E. Derm cells drastically reduced their susceptibility to EHV-1 infection. These results suggest that equine MHC class I is a functional gD receptor that plays a pivotal role in EHV-1 entry into equine cells.

Genetic diversity of the major histocompatibility complex class II in Alaskan caribou herds

We have sampled five different herds of caribou in Alaska to ascertain their major histocompatibility complex (MHC) class II diversity, and to assess whether the herds were significantly different in their MHC class II allele profiles. We complemented the MHC results with data from nine neutral microsatellite markers. The results indicate that while the microsatellites are diverse, there are no significant differences between the herds. However, for the MHC, we have shown that there is diversity at three of the four loci studied, the different herds have significantly different MHC class II allele profiles. It is also clear that although some of the herds have overlapping ranges, they are still different for their MHC class II alleles.

Maternal immune responses to trophoblast: the contribution of the horse to pregnancy immunology

The horse has proven to be a distinctively informative species in the study of pregnancy immunology for several reasons. First, unique aspects of the anatomy and physiology of the equine conceptus facilitate approaches that are not possible in other model organisms, such as non-surgical recovery of early stage embryos and conceptuses and isolation of pure trophoblast cell populations. Second, pregnant mares make strong cytotoxic antibody responses to paternal major histocompatibility complex class I antigens expressed by the chorionic girdle cells, permitting detailed evaluation of the antigenicity of these invasive trophoblasts and how they affect the maternal immune system. Third, there is abundant evidence for local maternal cellular immune responses to the invading trophoblasts in the pregnant mare. The survival of the equine fetus in the face of strong maternal immune responses highlights the complex immunoregulatory mechanisms that result in materno-fetal tolerance. Finally, the parallels between human and horse trophoblast cell types, their gene expression, and function make the study of equine pregnancy highly relevant to human health. Here, we review the most pertinent aspects of equine reproductive immunology and how studies of the pregnant mare have contributed to our understanding of maternal acceptance of the allogeneic fetus.

Polymorphism of sheep MHC Class IIb gene TAPASIN

The Major Histocompatibility Complex (MHC) is one of the most gene dense regions in the genome and studies in several species have shown significant associations between the MHC and disease. The endoplasmic reticular glycoprotein, tapasin, is involved in the MHC class I antigen presentation pathway. Sheep TAPASIN is located in the class IIb region of the MHC. Sheep TAPASIN was subcloned from BAC and cosmid genomic clones and DNA sequenced. TAPASIN is 9549bp in length and encodes a protein of 447 amino acids. The structure of sheep TAPASIN was similar to other mammals and consisted of eight exons with a distinctively larger intron between exon three and four. Sheep TAPASIN gene had high sequence identity with other mammalian TAPASINs. The TAPASIN gene sequence is conserved across many mammalian species and is possibly maintained through purifying selection with the average ratio of ds/dn of 3.9. Twenty-six SNPs in sheep TAPASIN were identified.

Analysis of MHC class I genes across horse MHC haplotypes

The genomic sequences of 15 horse major histocompatibility complex (MHC) class I genes and a collection of MHC class I homozygous horses of five different haplotypes were used to investigate the genomic structure and polymorphism of the equine MHC. A combination of conserved and locus-specific primers was used to amplify horse MHC class I genes with classical and nonclassical characteristics. Multiple clones from each haplotype identified three to five classical sequences per homozygous animal and two to three nonclassical sequences. Phylogenetic analysis was applied to these sequences, and groups were identified which appear to be allelic series, but some sequences were left ungrouped. Sequences determined from MHC class I heterozygous horses and previously described MHC class I sequences were then added, representing a total of ten horse MHC haplotypes. These results were consistent with those obtained from the MHC homozygous horses alone, and 30 classical sequences were assigned to four previously confirmed loci and three new provisional loci. The nonclassical genes had few alleles and the classical genes had higher levels of allelic polymorphism. Alleles for two classical loci with the expected pattern of polymorphism were found in the majority of haplotypes tested, but alleles at two other commonly detected loci had more variation outside of the hypervariable region than within. Our data indicate that the equine major histocompatibility complex is characterized by variation in the complement of class I genes expressed in different haplotypes in addition to the expected allelic polymorphism within loci.

Polymorphism and selection in the major histocompatibility complex DRA and DQA genes in the family Equidae

The major histocompatibility complex genes coding for antigen binding and presenting molecules are the most polymorphic genes in the vertebrate genome. We studied the DRA and DQA gene polymorphism of the family Equidae. In addition to 11 previously reported DRA and 24 DQA alleles, six new DRA sequences and 13 new DQA alleles were identified in the genus Equus. Phylogenetic analysis of both DRA and DQA sequences provided evidence for trans-species polymorphism in the family Equidae. The phylogenetic trees differed from species relationships defined by standard taxonomy of Equidae and from trees based on mitochondrial or neutral gene sequence data. Analysis of selection showed differences between the less variable DRA and more variable DQA genes. DRA alleles were more often shared by more species. The DQA sequences analysed showed strong amongst-species positive selection; the selected amino acid positions mostly corresponded to selected positions in rodent and human DQA genes.

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 map of the class III region of the sheep major histocompatibilty complex

Background

The central, or class III, region of the major histocompatibility complex (MHC) is an important gene rich sub-region of the MHC of mammals and contains many loci implicated in disease processes and potential productivity traits. As a prelude to identifying MHC loci associated with productivity traits in sheep, we have used BAC and cosmid libraries of genomic DNA to generate a physical map of the sheep MHC class III region. This map will facilitate association studies and provide insights into the distribution of recombination events in this chromosomal segment.

Results

Twenty eight sheep genes were identified in 10 BAC clones which spanned approximately 700 kbp of a chromosomal region adjacent to the class I region of the sheep MHC and which therefore covers most, if not all, of the class III of the sheep MHC. The relative positions of 17 of these genes was established as well as two additional groups of genes for which the intragroup order was not known. Cosmid mapping permitted a more detailed mapping of the complement genes present in the class III and showed a local inversion (relative to humans) of one pair of the duplicated complement C4 and CYP21 loci. A panel of 26 single nucleotide polymorphisms (SNPs) was identified in 10 loci, covering ≈600 kbp of the mapped region.

Conclusion

This report provides a physical map covering ≈700 kbp of the class III of the sheep MHC together with a SNP panel which will facilitate disease and productivity association studies. The presence of a local inversion (relative to humans) of one pair of the duplicated C4 and CYP21 loci and a previously described dinucleotide tandem repeat locus (BfMs) has been located within an intron of the SK12VL gene.

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.

Potential applications of equine genomics in dissecting diseases and fertility

Following the recent development of high-resolution gene maps and generation of several basic tools and resources to use them in analyzing traits that are economically important to horse owners, genome analysis in horses is witnessing a shift towards developing an ability to analyze complex traits. The likelihood of this happening in the very near future is great, mainly because of the recent availability of the whole genome sequence in the horse. The latter has triggered the development of novel tools like SNP-chip and expression arrays that will permit rapid genome-wide analysis. While these tools will be used for a range of multi-factorial disease traits, attempts are underway to develop focused tools that can target reproduction, fertility and sex determination. For this, a catalog of sex and reproduction related (SRR) genes is being developed in horses. A recently developed dense map of the horse Y chromosome will provide genes that are expressed exclusively in males and, therefore, have an impact on stallion fertility. Overall, these advances in equine genome analysis hold promise for improved diagnosis and treatment of various conditions in horses.

Comparative genomic structure of human, dog, and cat MHC: HLA, DLA, and FLA

Comparisons of the genomic structure of 3 mammalian major histocompatibility complexes (MHCs), human HLA, canine DLA, and feline FLA revealed remarkable structural differences between HLA and the other 2 MHCs. The 4.6-Mb HLA sequence was compared with the 3.9-Mb DLA sequence from 2 supercontigs generated by 7x whole-genome shotgun assembly and 3.3-Mb FLA draft sequence. For FLA, we confirm that 1) feline FLA was split into 2 pieces within the TRIM (member of the tripartite motif) gene family found in human HLA, 2) class II, III, and I regions were placed in the pericentromeric region of the long arm of chromosome B2, and 3) the remaining FLA was located in subtelomeric region of the short arm of chromosome B2. The exact same chromosome break was found in canine DLA structure, where class II, III, and I regions were placed in a pericentromeric region of chromosome 12 whereas the remaining region was located in a subtelomeric region of chromosome 35, suggesting that this chromosome break occurred once before the split of felid and canid more than 55 million years ago. However, significant differences were found in the content of genes in both pericentromeric and subtelomeric regions in DLA and FLA, the gene number, and amplicon structure of class I genes plus 2 other class I genes found on 2 additional chromosomes; canine chromosomes 7 and 18 suggest the dynamic nature in the evolution of MHC class I genes.

AnimalQTLdb: a livestock QTL database tool set for positional QTL information mining and beyond

The Animal Quantitative Trait Loci (QTL) database (AnimalQTLdb) is designed to house all publicly available QTL data on livestock animal species from which researchers can easily locate and compare QTL within species. The database tools are also added to link the QTL data to other types of genomic information, such as radiation hybrid (RH) maps, finger printed contig (FPC) physical maps, linkage maps and comparative maps to the human genome, etc. Currently, this database contains data on 1287 pig, 630 cattle and 657 chicken QTL, which are dynamically linked to respective RH, FPC and human comparative maps. We plan to apply the tool to other animal species, and add more structural genome information for alignment, in an attempt to aid comparative structural genome studies ().

A molecular approach to the identification of cytotoxic T-lymphocyte epitopes within equine herpesvirus 1

Equine herpesvirus 1 (EHV-1) causes respiratory and neurological disease and abortion in horses. Animals with high frequencies of cytotoxic T lymphocytes (CTL) show reduced severity of respiratory disease and frequency of abortion, probably by CTL-mediated control of cell-associated viraemia. This study aimed to identify CTL epitopes restricted by selected major histocompatibility complex (MHC) class I alleles expressed in the equine leukocyte antigen (ELA) A3 haplotype. Effector CTL were induced from EHV-1-primed ponies and thoroughbreds with characterized MHC class I haplotypes and screened against P815 target cells transfected with selected EHV-1 genes and MHC class I genes. Targets that expressed EHV-1 gene 64 and the MHC B2 gene were lysed by effector CTL in a genetically restricted manner. There was no T-cell recognition of targets expressing either the MHC B2 gene and EHV-1 genes 2, 12, 14, 16, 35, 63 or 69, or the MHC C1 gene and EHV-1 genes 12, 14, 16 or 64. A vaccinia virus vector encoding gene 64 (NYVAC-64) was also investigated. Using lymphocytes from ELA-A3 horses, the recombinant NYVAC-64 virus induced effector CTL that lysed EHV-1-infected target cells; the recombinant virus also supplied a functional peptide that was expressed by target cells and recognized in an MHC-restricted fashion by CTL induced with EHV-1. This construct may therefore be used to determine the antigenicity of EHV-1 gene 64 for other MHC haplotypes. These techniques are broadly applicable to the identification of additional CTL target proteins and their presenting MHC alleles, not only for EHV-1, but for other equine viruses.

YAC/BAC contig spanning the MHC class III region of cattle

A contig of the class III region of the bovine major histocompatibility complex (MHC) was established from bacterial and yeast artificial chromosomes using PCR and BAC-end sequencing. The marker content of individual clones was determined by gene and BAC-end specific PCR, and the location of genes and BAC-ends was confirmed analyzing somatic hybrid cells. A comparative analysis indicated that the content and order of MHC class III genes is strongly conserved between cattle and other mammalian species. Fluorescence in situ hybridization localized the bovine class III region to BTA23q21-->q22. The results show that the collection of sequenced BAC-ends is a powerful resource for generating high-resolution comparative chromosome maps.

Integration of hybridization-based markers (overgos) into physical maps for comparative and evolutionary explorations in the genus Oryza and in Sorghum

BACKGROUND

With the completion of the genome sequence for rice (Oryza sativa L.), the focus of rice genomics research has shifted to the comparison of the rice genome with genomes of other species for gene cloning, breeding, and evolutionary studies. The genus Oryza includes 23 species that shared a common ancestor 8-10 million years ago making this an ideal model for investigations into the processes underlying domestication, as many of the Oryza species are still undergoing domestication. This study integrates high-throughput, hybridization-based markers with BAC end sequence and fingerprint data to construct physical maps of rice chromosome 1 orthologues in two wild Oryza species. Similar studies were undertaken in Sorghum bicolor, a species which diverged from cultivated rice 40-50 million years ago.

RESULTS

Overgo markers, in conjunction with fingerprint and BAC end sequence data, were used to build sequence-ready BAC contigs for two wild Oryza species. The markers drove contig merges to construct physical maps syntenic to rice chromosome 1 in the wild species and provided evidence for at least one rearrangement on chromosome 1 of the O. sativa versus Oryza officinalis comparative map. When rice overgos were aligned to available S. bicolor sequence, 29% of the overgos aligned with three or fewer mismatches; of these, 41% gave positive hybridization signals. Overgo hybridization patterns supported colinearity of loci in regions of sorghum chromosome 3 and rice chromosome 1 and suggested that a possible genomic inversion occurred in this syntenic region in one of the two genomes after the divergence of S. bicolor and O. sativa.

CONCLUSION

The results of this study emphasize the importance of identifying conserved sequences in the reference sequence when designing overgo probes in order for those probes to hybridize successfully in distantly related species. As interspecific markers, overgos can be used successfully to construct physical maps in species which diverged less than 8 million years ago, and can be used in a more limited fashion to examine colinearity among species which diverged as much as 40 million years ago. Additionally, overgos are able to provide evidence of genomic rearrangements in comparative physical mapping studies.

Comparative analysis of the bovine MHC class IIb sequence identifies inversion breakpoints and three unexpected genes

The bovine major histocompatibility complex (MHC) or BoLA is organized differently from typical mammalian MHCs in that a large portion of the class II region, called class IIb, has been transposed to a position near the centromere on bovine chromosome 23. Gene mapping indicated that the rearrangement resulted from a single inversion, but the boundaries and gene content of the inverted segment have not been fully determined. Here, we report the genomic sequence of BoLA IIb. Comparative sequence analysis with the human MHC revealed that the proximal inversion breakpoint occurred approximately 2.5 kb from the 3' end of the glutamate-cysteine ligase, catalytic subunit (GCLC) locus and that the distal breakpoint occurred about 2 kb from the 5' end from a divergent class IIDRbeta-like sequence designated DSB. Gene content, order and orientation of BoLA IIb are consistent with the single inversion hypothesis when compared with the corresponding region of the human class II MHC (HLA class II). Differences with HLA include the presence of a single histone H2B gene located between the proteasome subunit, beta type, 9 (PSMB9) and DMB loci and a duplicated TAP2 with a variant splice site. BoLA IIb spans approximately 450 kb DNA, with 20 apparently intact genes and no obvious pseudogenes. The region contains 227 simple sequence repeats (SSRs) and approximately 167 kb of retroviral-related repetitive DNA. Nineteen of the 20 genes identified in silico are supported by bovine EST data indicating that the functional gene content of BoLA IIb has not been diminished because it has been transposed from the remainder of BoLA genes.

A BAC clone-based physical map of ovine major histocompatibility complex

An ovine bacterial artificial chromosome (BAC) library containing 190,000 BAC clones was constructed and subsequently screened to construct a BAC-based physical map for the ovine major histocompatibility complex (MHC). Two hundred thirty-three BAC clones were selected by 84 overgo probes designed on human, mouse, and swine MHC sequence homologies. Ninety-four clones were ordered by DNA fingerprinting to form contigs I, II, and III that correspond to ovine MHC class I-class III, class IIa, and class IIb. The minimum tiling paths of contigs I, II, and III are 15, 4, and 4 BAC clones, spanning approximately 1900, 400, and 300 kb, respectively. The order and orientation of most BAC clones in each contig were confirmed by BAC-end sequencing. An open gap exists between class IIa and class III. This work helps to provide a foundation for detailed study of ovine MHC genes and of evolution of MHCs in mammals.

The genomic sequence and analysis of the swine major histocompatibility complex

We describe the generation and analysis of an integrated sequence map of a 2.4-Mb region of pig chromosome 7, comprising the classical class I region, the extended and classical class II regions, and the class III region of the major histocompatibility complex (MHC), also known as swine leukocyte antigen (SLA) complex. We have identified and manually annotated 151 loci, of which 121 are known genes (predicted to be functional), 18 are pseudogenes, 8 are novel CDS loci, 3 are novel transcripts, and 1 is a putative gene. Nearly all of these loci have homologues in other mammalian genomes but orthologues could be identified with confidence for only 123 genes. The 28 genes (including all the SLA class I genes) for which unambiguous orthology to genes within the human reference MHC could not be established are of particular interest with respect to porcine-specific MHC function and evolution. We have compared the porcine MHC to other mammalian MHC regions and identified the differences between them. In comparison to the human MHC, the main differences include the absence of HLA-A and other class I-like loci, the absence of HLA-DP-like loci, and the separation of the extended and classical class II regions from the rest of the MHC by insertion of the centromere. We show that the centromere insertion has occurred within a cluster of BTNL genes located at the boundary of the class II and III regions, which might have resulted in the loss of an orthologue to human C6orf10 from this region.