- Research article
- Open Access
A physical map of a BAC clone contig covering the entire autosome insertion between ovine MHC Class IIa and IIb
BMC Genomics volume 13, Article number: 398 (2012)
The ovine Major Histocompatibility Complex (MHC) harbors genes involved in overall resistance/susceptibility of the host to infectious diseases. Compared to human and mouse, the ovine MHC is interrupted by a large piece of autosome insertion via a hypothetical chromosome inversion that constitutes ~25% of ovine chromosome 20. The evolutionary consequence of such an inversion and an insertion (inversion/insertion) in relation to MHC function remains unknown. We previously constructed a BAC clone physical map for the ovine MHC exclusive of the insertion region. Here we report the construction of a high-density physical map covering the autosome insertion in order to address the question of what the inversion/insertion had to do with ruminants during the MHC evolution.
A total of 119 pairs of comparative bovine oligo primers were utilized to screen an ovine BAC library for positive clones and the orders and overlapping relationships of the identified clones were determined by DNA fingerprinting, BAC-end sequencing, and sequence-specific PCR. A total of 368 positive BAC clones were identified and 108 of the effective clones were ordered into an overlapping BAC contig to cover the consensus region between ovine MHC class IIa and IIb. Therefore, a continuous physical map covering the entire ovine autosome inversion/insertion region was successfully constructed. The map confirmed the bovine sequence assembly for the same homologous region. The DNA sequences of 185 BAC-ends have been deposited into NCBI database with the access numbers HR309252 through HR309068, corresponding to dbGSS ID 30164010 through 30163826.
We have constructed a high-density BAC clone physical map for the ovine autosome inversion/insertion between the MHC class IIa and IIb. The entire ovine MHC region is now fully covered by a continuous BAC clone contig. The physical map we generated will facilitate MHC functional studies in the ovine, as well as the comparative MHC evolution in ruminants.
The mammalian Major Histocompatibility Complex (MHC) harbors genes involved in overall resistance/susceptibility of animals to infectious pathogens, including viral, bacterial, internal and external parasites. Pathogens serve as sources of selection pressure to their host animals, and the hosts are forced to develop effective strategies to fight against the pathogens in various environments. Such co-evolutionary struggles may have left distinct marks in the genome of each species involved, and mammalian MHC regions have been shaped into clusters of immunological gene families by such host-pathogen interactions, probably via functional gene duplications [1–3]. The implications of ovine MHC molecules in providing protection against pathogens [4–8] and the associated structures of the artiodactyl’s MHC region in general have led to a number of studies into the sheep MHC [9–15].
The ovine MHC, also called ovine leukocyte antigen (OLA), is located on the long arm of ovine chromosome 20 (OAR 20q15–20q23) with a similar structure and organization to that of human and other mammals . The literature shows that MHC genes play vital roles in resistance of animals to foot rot , parasites , and bovine leukemia virus . To date, the majority of studies on the structure and organization of the ovine MHC have focused on the gene content and polymorphism of the class II region [18–23]. Although most loci in the sheep MHC are found to be homologous to their counterparts in the human MHC [12, 21, 24, 25], there are significant differences. Examples of such differences include the DP loci in human being replaced by DY in sheep [19, 21, 26, 27], and the number of DQA loci varying significantly among sheep breeds [20, 22, 28].
Compared to human and mouse, the structure of the sheep MHC is interrupted by a piece of ~14 Mb autosome insertion, possibly via a hypothetical chromosome inversion (inversion/insertion) in the class II region, similar to that of cattle [24, 29–32]. The inversion/insertion constitutes ~25% of ovine chromosome 20, which spliced the MHC class II region into IIa and IIb. The significance of such an insertion in relation to the ovine MHC functions remains unknown. The evolutionary consequence of such an event is also worthy of attention, because some of the ovine-specific MHC loci like DY, and Dsb are located near the boundary region of the inversion/insertion. We previously constructed a physical map of BAC clone contigs covering the ovine MHC except the autosome insertion region [12, 13], and a high accuracy sequence map of sheep OLA was accordingly constructed .
With the initial release of sheep whole genome reference sequences by the International Sheep Genomic Consortium (ISGC), much more genome sequence information is now accessible for functional and comparative studies . Nevertheless, the sequence map would serve the research community even better if it is cross-referenced/checked for accuracy in DNA sequence and assembly, at least for some chromosome regions, by an alternative approach. In this regard, the detailed information is still not fully available for the gene structure, organization, and DNA sequence for the ovine chromosome region between OLA class IIa and IIb [12, 14, 27].
In this paper, we describe the construction of a BAC physical map covering the entire autosome insertion between ovine MHC class IIa and IIb. Because ovine and bovine species share the consensus structure and organization in the entire MHC region [24, 29–32], we used comparative approaches to screen a sheep BAC library with 119 bovine oligo nucleotide primers designed from the bovine genomic sequences for the consensus region. The order and overlapping relationship of the identified BAC clones were determined by DNA fingerprinting, BAC-end sequencing, and sequence-specific PCR. A total of 108 effective overlapping BAC clones were selected to fully cover the region between class IIa and IIb. The physical map we constructed will help to generate ovine MHC sequencing map with a high level of accuracy, which in turn will facilitate MHC functional and comparative MHC evolution studies in ruminants.
Comparative design of oligo primers
A BAC library was previously constructed using the genome DNA from a male Chinese merino sheep, with a total of 190,500 BAC clones and an average insert length of 133 kb [12, 13]. To screen the BAC library for positive clones in the target genome region between ovine MHC class IIa and IIb, we adapted a comparative strategy to design bovine oligo nucleotide primers using the bovine reference DNA sequences in the consensus genome region . At the time this study was conducted, no sheep genomic sequence was publicly available for the genome region of our concern. Bovine DNA sequences of homologous genes, exon, intron, or partial STS sequences were acquired from the NCBI website (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/genome/sts/). Primers were designed along the bovine MHC region between class IIa and IIb, approximately 80–160 kb apart between two neighbor loci using the software Prime Primer 5.0 (Biosoft International, CA). A total of 119 bovine primer pairs were designed for screening the sheep genomic BAC library (Table1).
BAC library organization and screening
To facilitate large scale PCR screening, all the 190,500 clones of the BAC library were organized into 3-dimensional BAC clone pools of plates, rows, and columns. Random BAC clones from each of 496 permanent 384-well storage plates were duplicated onto a Luria-Bertani (LB) agar plate for overnight growth at 37°C, using a 384-pin Multi-Blot Replicator as tool for BAC clone duplication (V & P Scientific, Inc., San Diego, CA). The overnight E. coli colonies were then harvested and pooled for plate (n = 496), row (n = 16), or column (n = 24). The standard alkaline lyses methodology was adapted for isolation of the pooled BAC plasmid DNA and the resulting DNA was assembled into super plates for routine PCR screening . The first dimension of the BAC clone pool consisted of 496 DNA samples, each representing one of 496 BAC plates (P001-P496). The second and third dimension consisted of 16 and 24 DNA samples, respectively, for the pooled 16 rows (R01-R16) and 24 columns (C01-C24) of the random BAC clones.
To screen the BAC library using each of 119 pairs of comparative oligo primer pairs, the diluted DNA from each well of the super pool plates was used as a DNA template. The individual PCR reaction was adapted in a total of 10 μl reaction volume with 50 μM of dNTPs, 1.5 mM Mg++, 0.2 μM of each primer pair, 1 × PCR buffer, and 0.1 unit of Tag DNA polymerase. The PCR products were resolved by 1.5% agarose gel electrophoresis and the specific PCR fragment band with the expected size indicated a potential positive BAC clone for the gene loci of oligo primers used. The exact location of the target clone in the BAC library was determined by sequential PCRs using the super row and super column DNA as templates, respectively.
DNA fingerprinting and contig assembling
DNA fingerprinting was performed to determine the overlapping relationship among the identified positive BAC clones . DNA from the positive BAC clone was purified from host E. coli by QIAGEN column and subjected for complete restriction enzyme digestion using Hin dIII. The enzyme digested products were analyzed on 1% TAE agarose gel electrophoresis for recoding of DNA fragment patterns. The fingerprinting images were captured with UVP Labworks System (UVP Inc., Upland, CA) for systematic analysis. Restriction fragment patterns were analyzed to identify overlapping BAC clones, which were then manually assembled into draft contigs based on the modified methods of Marra  and Soderlund .
BAC-end sequencing was performed for the selected clones to facilitate verification of the overlapping relationships of the BAC clones. The sequencing was performed on an ABI 3730X DNA analyzer at the core facilities of the Institute of Genetics and Developmental Biology, the Chinese Academy of Sciences. The oligo nucleotide primers used for the DNA sequencing were Copycontrol pCC1BAC vector-derived sequencing primer T7 (5’-TAATACGACTCACTATAGGG3’), pCC1/pEpiFOS RP-2 (abbr. RP2) (5’-TACGCCAAGCTATTTAGGTGAGA-3’), and pCC1/pEpiFOS RP-1(abbr. RP-1) (5'-CTCGTATGTTGTGTGGAATTGTGAGC-3'). The resulting sequences were analyzed for overlapping, and used as templates for oligo primer design. Based on the sequence data generated by BAC-end sequencing, PCR primers (Additional file 1: Table S1) were designed to amplify the common genetic loci in two overlapped BACs for confirmation. Sequence-Specific PCRs (SP-PCRs) were performed in 20 μl system including approximately 2 ng BAC DNA, 0.5 U Taq DNA polymerase, 0.1 mM dNTPs, 1.5 mM Mg++, 0.25 μM each primer, and 1× PCR buffer. When necessary, the PCR products were verified by cloning the fragments into a TA vector for verifying DNA sequencing.
Assemble of the BAC clone contig
A continuous BAC clone contig was eventually assembled based on the integrated results of DNA fingerprinting, BAC-end sequencing, and sequence specific PCR amplification of the common loci on the overlapping clones. Redundant BAC clones were removed from the assembly based on the necessity and the relative contribution of each overlapping BACs on the contig. Gaps in the contig were closed by the repeated cycles of PCR screening of BAC clones, DNA fingerprinting of additional BAC clones identified, BAC-end sequencing, and SP-PCR verification. Additional effort was made to link the existing BAC clone contig to the physical map constructed previously, for a complete physical map covering the entire ovine MHC including the autosome insertion between class IIa and IIb.
For comparison of the MHC structure and organization between sheep and other mammals, multiple comparisons were performed for the representative MHC and extended DNA sequences from human, chimpanzees, mouse, cattle, and sheep. Sequence data were downloaded from the NCBI database and other related public websites designated for the sheep genomic information.
Target BAC identification
We successfully identified a total of 368 positive BAC clones for ovine chromosome 20 between MHC class IIa and IIb, utilizing bovine primers designed from the consensus genome region (Table1). Out of 119 pairs of oligo primers designed, 92 pairs worked effectively to generate specific target gene fragments of the expected sizes. This approach resulted in the successful identification of positive ovine BAC clones in the target genome region, and the overall efficiency of comparative PCR reached 80%. The relatively high rate of success for the comparative SP-PCR not only facilitated our mapping efforts, but also helped to confirm the homologous nature of MHC regions between bovine and ovine species.
Organization of ~190,500 random ovine BAC clones into three dimensional super DNA pool of rows (n = 16), columns (n = 24), and plates (n = 496) significantly increased the efficiency of PCR screening of the sheep BAC library (Figure1). The whole BAC library of 8.4× genome equivalents was screened through with a maximum of 536 (=496 + 16 + 24) PCR reactions, and a positive BAC clone could be frequently identified by as few as 136 (=96 + 16 + 24) PCR reactions using the super pool DNA as templates. In addition, PCR-based BAC clone screening also helped to eliminate the need for hybridization-based screening using radioactive 32P labeling.
DNA fingerprinting and BAC-end sequencing
The initial order of the positive BAC clones identified was successfully determined by inferring the overlapping relationships among the clones via DNA fingerprinting, using Hin dIII for restriction enzyme digestion of the BAC clone DNAs (Figure2). Out of 368 positive BAC clones subjected for the DNA fingerprinting, 185 clones with their overlapping relationships were successfully determined. The resulting BAC contig covered the entire autosome insertion region between the MHC class IIa and IIb. After removing the redundant clones, a total of 108 effective BACs were ordered to form an overlapping BAC contig (Additional file 1: Table S1).
For cross-checking of the clone order, BAC-end sequencing was performed for all overlapping BAC clones, and the sequences generated were used to design BAC-end oligo primers (Additional file 1: Table S1) for further verification of overlapping relationships. The sequences of 185 BAC-ends have been deposited into the NCBI database with the access number HR309252 through HR309068, corresponding to dbGSS ID 30164010 through 30163826.
Cross verification and physical map assembling
For additional cross-verification of the BAC clone orders, a total of 108 pairs of BAC-end oligo primers were designed for amplification by PCR of the common loci in two overlapping BACs (Figure3). Verification PCR confirmed the results of DNA fingerprinting at a high level of accuracy. Out of the 108 primer pairs used, 103 produced the specific PCR products with the expected size, the overall success rate reached 95% (Additional file 1: Table S1). An overlapping relationship between two BACs was further verified if the common target loci were detected from both BACs in the overlapped region. A total of five pairs of oligo primers failed to generate the specific PCR band, or failed to produce the PCR fragment at the expected size.
A complete physical map of a BAC clone contig for the ovine MHC region between class IIa and IIb was successfully assembled (Figure4), based on the integrated results of DNA fingerprinting, BAC-end sequencing, and confirmation PCR of the BAC ends. The fully assembled physical map was composed of 108 effective ovine BAC clones organized into a continuous contig that covered the entire region between ovine MHC class IIa and IIb (Figure4). Based on the results of DNA fingerprinting, no gaps exist in the constructed BAC clone physical map which spans approximately 14 Mb genome region of ovine chromosome 20, indicating the even distribution of BAC clones in the library we previously constructed.
Using the comparative approaches, we successfully constructed a 14 Mb BAC clone contig map for a region in ovine chromosome 20 that harbors the MHC. Comparison between the identified ovine BAC contig and the orthologous bovine genomic region showed that the two species share essentially the same genomic structure and organization for the entire inversion/insertion between MHC class IIa and IIb (Figure5). For the available genetic loci generated via the SP-PCR and BAC-end sequencing, our results essentially confirmed the sheep genome sequence assembly presented by ISGC in the MHC region .
The physical map of ovine BAC contig we constructed helped to provide additional evidence to support the hypothesis that, there was an ancient chromosome rearrangement in the ancestor of ruminants which shaped the MHC structures currently observed in the ovine and bovine (Figure5). It is obvious that the MHC region in human, mouse and chimpanzees is continuous with no interruption, but in bovine and ovine it is interrupted by a large piece of autosome insertion which divided MHC class II into IIa and IIb subregions (Figure5). Given the fact of opposite loci order and orientation for the insertion region in ovine and bovine relative to those of human and mouse, it is highly possible that an event of genetic recombination occurred to the ancestor chromosome of ruminants, probably via chromosome looping and the subsequent crossover. This possibility was suggested by researchers previously [29, 38].
Examination of the bovine DNA sequence from the public database showed that the total length of bovine MHC is ~20 Mb, including the extended Class IIb region . However, the total length of the orthologous ovine MHC was ~14.3 Mb as determined in this study, which is approximately 5.7 Mb shorter than the MHC of bovine. On the other hand, the sequence of the same bovine region presented in the NCBI database is ~18 Mb in length (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/projects/mapview/maps.cgi?taxid=9913&chr=23). These discrepancies may not likely be resolved unless highly accurate sequence maps for the entire MHC regions become available.
The reliability of the ovine BAC contig map reported here is sufficiently high in theory, partially due to the fact that the DNA fingerprinting was utilized to infer the BAC clone orders, plus the results were cross-verified by both of the BAC-end sequencing and SP-PCR amplification of the target loci. However, it is not escaped from our attention that there are 5 out of the 108 overlapping locations in the BAC map where the SP-PCR failed to generate the expected PCR products between the overlapping BAC clones (data not shown). The significance of such failure in relation to the overall quality of the map remains to be determined. The possible explanations include the error in SP-PCR primer sequences, the high level of heterogeneity or polymorphism of the target locus involved, or the mistake in the interpretation of results of DNA fingerprinting.
Combined with our previous BAC physical map for the ovine MHC, we have now assembled a completed BAC clone physical map with the inversion/insertion region included (Additional file 2: Figure S1). The physical map will help to generate an ovine MHC sequencing map with a high level of accuracy, which in turn will facilitate MHC functional studies and comparative MHC evolution studies in ruminants. DNA sequencing of the BACs is currently underway.
We constructed a high-density physical map for the sheep genome region between MHC class IIa and IIb via comparative approaches. A total of 108 effective ovine BAC clones were selected to form a continuous BAC contig that covers the entire non-MHC insertion. The map spans approximately 14 Mb in length, constituting ~25% of ovine chromosome 20. The entire ovine MHC region, including the autosome insertion for which the physical map has been constructed, is now fully covered by a continuous BAC clone contig. The accuracy of DNA sequences play vital roles in detailed SNP and other functional studies of MHC genes, as well as for genome evolution studies. The physical map will help to generate ovine MHC sequencing map with a high level of accuracy, which in turn will facilitate MHC functional studies, as well as the comparative MHC evolution in ruminants.
Nei M, Rooney AP: Concerted and birth-and-death evolution of multigene families. Annu Rev Genet. 2005, 39: 121-152.
Spurgin LG, Richardson DS: How pathogens drive genetic diversity: MHC, mechanisms and misunderstandings. Proc Biol Sci. 2010, 277 (1684): 979-988.
Trowsdale J: The MHC, disease and selection. Immunol Lett. 2011, 137 (1–2): 1-8.
Bonneaud C, Richard M, Faivre B, Westerdahl H, Sorci G: An Mhc class I allele associated to the expression of T-dependent immune response in the house sparrow. Immunogenetics. 2005, 57 (10): 782-789.
Dukkipati VS, Blair HT, Garrick DJ, Murray A: Ovar-Mhc–ovine major histocompatibility complex: role in genetic resistance to diseases. N Z Vet J. 2006, 54 (4): 153-160.
Galindo RC, Ayoubi P, Garcia-Perez AL, Naranjo V, Kocan KM, Gortazar C, de la Fuente J: Differential expression of inflammatory and immune response genes in sheep infected with Anaplasma phagocytophilum. Vet Immunol Immunopathol. 2008, 126 (1–2): 27-34.
Konnai S, Takeshima SN, Tajima S, Yin SA, Okada K, Onuma M, Aida Y: The influence of ovine MHC class II DRB1 alleles on immune response in bovine leukemia virus infection. Microbiol Immunol. 2003, 47 (3): 223-232.
Mena A, Nichani AK, Popowych Y, Ioannou XP, Godson DL, Mutwiri GK, Hecker R, Babiuk LA, Griebel P: Bovine and ovine blood mononuclear leukocytes differ markedly in innate immune responses induced by Class A and Class B CpG-oligodeoxynucleotide. Oligonucleotides. 2003, 13 (4): 245-259.
Buitkamp J, Filmether P, Stear MJ, Epplen JT: Class I and class II major histocompatibility complex alleles are associated with faecal egg counts following natural, predominantly Ostertagia circumcincta infection. Parasitol Res. 1996, 82 (8): 693-696.
Dukkipati VS, Blair HT, Garrick DJ, Murray A: 'Ovar-Mhc' - ovine major histocompatibility complex: structure and gene polymorphisms. Genet Mol Res. 2006, 5 (4): 581-608.
Gruszczynska J, Brokowska K, Charon KM, Swiderek WP: Restriction fragment length polymorphism of exon 2 Ovar-DRB1 gene in Polish Heath Sheep and Polish Lowland Sheep. J Appl Genet. 2005, 46 (3): 311-314.
Liu H, Liu K, Wang J, Ma RZ: A BAC clone-based physical map of ovine major histocompatibility complex. Genomics. 2006, 88 (1): 88-95.
Liu K, Zhang P, Gao J, Liu H, Li G, Qiu Z, Zhang Y, Ren J, Tan P, Ma RZ: Closing a gap in the physical map of the ovine major histocompatibility complex. Anim Genet. 2011, 42 (2): 204-207.
Gao J, Liu K, Liu H, Blair HT, Li G, Chen C, Tan P, Ma RZ: A complete DNA sequence map of the ovine major histocompatibility complex. BMC Genomics. 2010, 11: 466-
Miltiadou D, Ballingall KT, Ellis SA, Russell GC, McKeever DJ: Haplotype characterization of transcribed ovine major histocompatibility complex (MHC) class I genes. Immunogenetics. 2005, 57 (7): 499-509.
Mahdy EA, Makinen A, Chowdhary BP, Andersson L, Gustavsson I: Chromosomal localization of the ovine major histocompatibility complex (OLA) by in situ hybridization. Hereditas. 1989, 111 (1): 87-90.
Escayg AP, Hickford JG, Bullock DW: Association between alleles of the ovine major histocompatibility complex and resistance to footrot. Res Vet Sci. 1997, 63 (3): 283-287.
Ballingall KT, Fardoe K, McKeever DJ: Genomic organisation and allelic diversity within coding and non-coding regions of the Ovar-DRB1 locus. Immunogenetics. 2008, 60 (2): 95-103.
Deverson EV, Wright H, Watson S, Ballingall K, Huskisson N, Diamond AG, Howard JC: Class II major histocompatibility complex genes of the sheep. Anim Genet. 1991, 22 (3): 211-225.
Escayg AP, Hickford JG, Montgomery GW, Dodds KG, Bullock DW: Polymorphism at the ovine major histocompatibility complex class II loci. Anim Genet. 1996, 27 (5): 305-312.
Scott PC, Choi CL, Brandon MR: Genetic organization of the ovine MHC class II region. Immunogenetics. 1987, 25 (2): 116-122.
Snibson KJ, Maddox JF, Fabb SA, Brandon MR: Allelic variation of ovine MHC class II DQA1 and DQA2 genes. Anim Genet. 1998, 29 (5): 356-362.
van der Poel JJ, Groenen MA, Dijkhof RJ, Ruyter D, Giphart MJ: The nucleotide sequence of the bovine MHC class II alpha genes: DRA, DOA, and DYA. Immunogenetics. 1990, 31 (1): 29-36.
Childers CP, Newkirk HL, Honeycutt DA, Ramlachan N, Muzney DM, Sodergren E, Gibbs RA, Weinstock GM, Womack JE, Skow LC: Comparative analysis of the bovine MHC class IIb sequence identifies inversion breakpoints and three unexpected genes. Anim Genet. 2006, 37 (2): 121-129.
Devilee P, Warnaar JN, Giphart MJ: MHC homology between various mammalian species at the DNA level: its relevance to MHC evolution. Exp Clin Immunogenet. 1984, 1 (2): 90-98.
Qin J, Mamotte C, Cockett NE, Wetherall JD, Groth DM: A map of the class III region of the sheep major histocompatibilty complex. BMC Genomics. 2008, 9: 409-
Wright H, Ballingall KT, Redmond J: The DY sub-region of the sheep MHC contains an A/B gene pair. Immunogenetics. 1994, 40 (3): 230-234.
Hickford JG, Ridgway HJ, Escayg AP: Evolution of the ovine MHC DQA region. Anim Genet. 2000, 31 (3): 200-205.
Amills M, Ramiya V, Norimine J, Lewin HA: The major histocompatibility complex of ruminants. Rev Sci Tech. 1998, 17 (1): 108-120.
Everts-van Der Wind A, Kata SR, Band MR, Rebeiz M, Larkin DM, Everts RE, Green CA, Liu L, Natarajan S, Goldammer T, et al: A 1463 gene cattle-human comparative map with anchor points defined by human genome sequence coordinates. Genome Res. 2004, 14 (7)): 1424-1437.
Maddox JF, Davies KP, Crawford AM, Hulme DJ, Vaiman D, Cribiu EP, Freking BA, Beh KJ, Cockett NE, Kang N, et al: An enhanced linkage map of the sheep genome comprising more than 1000 loci. Genome Res. 2001, 11 (7): 1275-1289.
McShane RD, Gallagher DS, Newkirk H, Taylor JF, Burzlaff JD, Davis SK, Skow LC: Physical localization and order of genes in the class I region of the bovine MHC. Anim Genet. 2001, 32 (5): 235-239.
Archibald AL, Cockett NE, Dalrymple BP, Faraut T, Kijas JW, Maddox JF, McEwan JC, Hutton Oddy V, Raadsma HW, Wade C, et al: The sheep genome reference sequence: a work in progress. Anim Genet. 2010, 41 (5): 449-453.
Elsik CG, Tellam RL, Worley KC, Gibbs RA, Muzny DM, Weinstock GM, Adelson DL, Eichler EE, Elnitski L, Guigo R, et al: The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science. 2009, 324 (5926): 522-528.
Osoegawa K, Woon PY, Zhao B, Frengen E, Tateno M, Catanese JJ, de Jong PJ: An improved approach for construction of bacterial artificial chromosome libraries. Genomics. 1998, 52 (1): 1-8.
Marra MA, Kucaba TA, Dietrich NL, Green ED, Brownstein B, Wilson RK, McDonald KM, Hillier LW, McPherson JD, Waterston RH: High throughput fingerprint analysis of large-insert clones. Genome Res. 1997, 7 (11): 1072-1084.
Soderlund C, Longden I, Mott R: FPC: a system for building contigs from restriction fingerprinted clones. Comput Appl Biosci. 1997, 13 (5): 523-535.
Lewin HA, Russell GC, Glass EJ: Comparative organization and function of the major histocompatibility complex of domesticated cattle. Immunol Rev. 1999, 167: 145-158.
The authors are very appreciative of the expert reviewers who helped to improve the quality of the manuscript significantly. This work was funded by research grants from National Natural Science Foundation of China (30125024; 30771148), Ministry of Science and Technology of China (2006DFB33750; 2010CB530204), and China Ministry of Agriculture (2009ZX08008-005B).
Authors declare no conflict of interests.
GL carried out BAC library organization and SP-PCR screening. KL carried out DNA fingerprinting and contig assembling. SJ and GL performed oligo primer design and BAC-end sequencing. HL constructed the sheep BAC library. HB carried out data analysis. XC carried out certain verification experiments. PT and PZ carried out data cross checking. RM and JG supervised the studies and wrote the manuscript. All authors read and approved the final version of the manuscript.
Gang Li, Ka Liu contributed equally to this work.
Electronic supplementary material
Additional file 2: Figure S1. A complete physical map of entire ovine MHC with the insertion region between class IIa and IIb included. Order and orientation of overlapping BAC clones were jointly determined by combinations of DNA fingerprinting, BAC-end sequencing, and sequence-specific PCR. Genes identified by BAC-end sequencing are marked with erect black lines, with their names listed above. A horizontal bar stands for individual BAC with its identification marked above. Red, purple and green color represent the MHC class I, class III, and class II, representatively. (DOC 86 KB)
About this article
Cite this article
Li, G., Liu, K., Jiao, S. et al. A physical map of a BAC clone contig covering the entire autosome insertion between ovine MHC Class IIa and IIb. BMC Genomics 13, 398 (2012) doi:10.1186/1471-2164-13-398
- Physical map
- Comparative mapping