Construction and characterization of a genomic BAC library for the Mus m. musculus mouse subspecies (PWD/Ph inbred strain)
BMC Genomics volume 6, Article number: 161 (2005)
The genome of classical laboratory strains of mice is an artificial mosaic of genomes originated from several mouse subspecies with predominant representation (>90%) of the Mus m. domesticus component. Mice of another subspecies, East European/Asian Mus m. musculus, can interbreed with the classical laboratory strains to generate hybrids with unprecedented phenotypic and genotypic variations. To study these variations in depth we prepared the first genomic large insert BAC library from an inbred strain derived purely from the Mus m. musculus-subspecies. The library will be used to seek and characterize genomic sequences controlling specific monogenic and polygenic complex traits, including modifiers of dominant and recessive mutations.
A representative mouse genomic BAC library was derived from a female mouse of the PWD/Ph inbred strain of Mus m. musculus subspecies. The library consists of 144 768 primary clones from which 97% contain an insert of 120 kb average size. The library represents an equivalent of 6.7 × mouse haploid genome, as estimated from the total number of clones carrying genomic DNA inserts and from the average insert size. The clones were arrayed in duplicates onto eight high-density membranes that were screened with seven single-copy gene probes. The individual probes identified four to eleven positive clones, corresponding to 6.9-fold coverage of the mouse genome. Eighty-seven BAC-ends of PWD/Ph clones were sequenced, edited, and aligned with mouse C57BL/6J (B6) genome. Seventy-three BAC-ends displayed unique hits on B6 genome and their alignment revealed 0.92 single nucleotide polymorphisms (SNPs) per 100 bp. Insertions and deletions represented 0.3% of the BAC end sequences.
Analysis of the novel genomic library for the PWD/Ph inbred strain demonstrated coverage of almost seven mouse genome equivalents and a capability to recover clones for specific regions of PWD/Ph genome. The single nucleotide polymorphism between the strains PWD/Ph and C57BL/6J was 0.92/100 bp, a value significantly higher than between classical laboratory strains. The library will serve as a resource for dissecting the phenotypic and genotypic variations between mice of the Mus m. musculus subspecies and classical laboratory mouse strains.
PWD/Ph is a highly inbred strain currently at 81 generations of brother × sister matings. It originated from the Mus m. musculus mouse subspecies by systematic inbreeding of a pair of wild mice trapped in 1972 [1, 2]. The mouse subspecies M. m. musculus and M. m. domesticus diverged from their common ancestor about 300 thousand years  to 1 million years ago  and at present they display signs of incomplete reproductive isolation [5–7]. As a consequence of the interrupted gene flow between both subspecies, the mice of the PWD/Ph strain exhibit a high degree of DNA polymorphisms and a broad range of phenotypic differences when compared to classical laboratory strains [2, 8]. Because of this unique feature, the PWD/Ph inbred strain has been nominated among 15 mouse strains, the genomes of which are being resequenced using high-density oligonucleotide array technology by Perlegen Sciences, Inc. . Moreover, PWD/Ph serves as the chromosome donor strain in construction of a set of C57BL/6-ChrPWD chromosome substitution strains (Gregorova, Forejt and coworkers, in preparation).
Bacterial Artificial Chromosome (BAC) genomic libraries are source of large genomic DNA insert clones for sequencing projects, physical mapping and isolation of intact genes [10, 11]. Although BAC clones may carry large inserts of genomic DNA (up to 200 kb) they display low rate of de novo rearrangements and are easy to handle. These features are in strong favor of the BAC libraries over the Yeast Artificial Chromosome (YAC) libraries, which can contain up to 60% of chimeric clones . Transgenic mice can be generated using BAC clones to examine candidate genes in context of all regulatory DNA elements required for their function and the phenotype of a mutant mouse can be rescued by BAC transgenesis [13, 14]. Moreover a targeted modification at exact positions within a genomic BAC clone can be introduced by recombineering [15, 16].
Here we report construction and characterization of the PWD/Ph BAC library, the first genomic library of the Mus m. musculus mouse subspecies. This library together with the upcoming panel of chromosome substitution strains will serve as a tool for analysis of complex traits by taking advantage of the evolutionary divergence between the two closely related mouse subspecies.
Results and discussion
Construction of the PWD/Ph-BAC library
The BAC library was prepared by cloning the Eco RI-partially digested genomic DNA from the spleen of a PWD/Ph female mouse in the vector pBACe3.6. Female DNA was chosen to gain an unbiased representation of the X chromosome in the library. The primary clones were picked and arrayed in 377 individual 384-well plates. The library consists of two segments containing 54 144 and 90 624 clones, respectively. Together 144 768 primary clones were arrayed on eight high-density nylon membranes (18 342 clones in doublets per membrane). The high-density membranes were utilized in subsequent hybridization experiments.
Average insert size of the library
The average insert size of the library was determined on a set of 400 randomly selected BAC clones. DNA samples were prepared from 164 and 236 BAC clones from the library segments 1 and 2, respectively, and subjected to Not I restriction analysis. The products of the digestion reactions were resolved by pulsed-field gel electrophoresis (PFGE) along with the high molecular weight markers. The average insert size for the first and the second library segment was 101.1 kb (SD ± 21.4) and 129.5 kb (SD ± 14.7), respectively (Figure 1). In the first and second library segments 6% and 1.2% clones were observed without insert, corresponding to 97% of insert-containing clones for the entire library. Estimation of 6.7-fold redundancy of the library was based on the average insert size (120 kb) and 2.6 × 109 bp size of the mouse genome.
Library screening and BAC end sequencing
A probability to find any given unique sequence in the library is 99.85%, according to the published formula (P = 1 - eN.ln(1-I/GS), where P is probability, N is number of clones, I is insert size, and GS is size of genome) . To further assess the genome redundancy and possible cloning bias of the library experimentally we performed a screening of the library with 7 single-copy gene probes. The probes were designed to amplify PCR products on the PWD/Ph genomic DNA template (Table 1). Seven probes detected in total 48 positive clones by hybridization on 8 high-density library membranes, 4 to 11 clones for each individual probe. The average number of clones recognized by a single probe was 6.9, in good accordance with the assessment of the library redundancy based on the average insert size.
To characterize the inserts of the PWD/Ph BAC library at the DNA sequence level we sequenced and manually edited 87 BAC ends from 47 BAC clones (total 38,339 nucleotides). The BAC end sequences (BESs) were masked for repeats and aligned on the C57BL/6J mouse genome. BES pairs of 29 BAC clones mapped to unique positions in the B6 genome on the opposite DNA strands within the distance up to 200 kb (Additional file 1). The mapping allowed us to estimate the average insert size of the BAC clones based on their locations on the B6 genome as 127 kb, which was slightly higher estimate than the average insert size acquired by restriction analysis (120 kb). These values corresponded well with the average insert size calculated for another set of clones recovered by the library screening described above (Table 2). A BES pair belonging to the clone 307-9O mapped to two distinct chromosomes. Whether it represents a chimeric insert or a chromosomal rearrangement in the PWD/Ph genome remains to be determined by fluorescence in situ hybridization (FISH) analysis. For each of the additional 13 BAC clones we found unambiguous positions for only one BES of a pair. Mapping of remaining 14 BESs was prevented by a high content of repetitive elements.
Analysis of SNPs and DNA polymorphism
To find out the degree of nucleotide polymorphism between the PWD/Ph and C57BL/6J mouse strains, we aligned 73 uniquely mapped BESs (32,182 nucleotides) with their C57BL/6J genomic counterparts and found 297 single nucleotide substitutions. The calculated SNP rate 0.92 per 100 bp is significantly higher than SNP frequency between laboratory strains [18–20] and corresponds well to the rate between the closely related subspecies Mus m. molossinus and the C57BL/6J strain (0.96%) . The insertions and deletions (indels) were found with lower frequency than SNPs: single nucleotide indels occurred with frequency 0.19% while multinucleotide indels with only 0.08% frequency. All nucleotide changes observed in the alignments of 87 PWD/Ph BESs and their B6 counterparts are summarized in Additional files 2 and 3. The high number of SNPs of the PWD/Ph strain is reflected by a high frequency of genetic and phenotypic variations between PWD/Ph and B6 inbred mice. An initial study performed to compare behavior of the PWD/Ph inbred strain with the B6 revealed substantial behavioral differences between these two strains . Using dense SNP maps of various laboratory and wild-derived inbred strains [20, 22] it will be possible to map genes responsible for particular complex traits more efficiently. For ultimate validation of candidate genes genomic BAC libraries will be highly desirable.
The first genomic BAC library was constructed for the Mus m. musculus subspecies of the house mouse, represented by the PWD/Ph inbred strain. The quality of the PWD BAC library was verified by hybridization with seven unique probes that identified multiple positive clones. BAC end sequencing provided a new piece of evidence on the high incidence of SNPs (0.92/100 bp) between C57BL/6J and PWD/Ph inbred strains. The mouse PWD/Ph BAC library will serve as a tool for functional genomics of complex genetic traits with the ultimate goal to identify and clone responsible genes. The PWD BAC library will become accessible to the scientific community via RZPD, Berlin, Germany .
Mouse strain and DNA isolation
Mouse manipulation was in accordance with the Czech Animal Protection Act No. 246/92, 162/93, and decrees No. 311/97, fully compatible with the NIH Publication No. 85-23, revised 1985. A female mouse of the PWD/Ph inbred strain, derived from the Mus mus musculus subspecies  was used for high molecular weight DNA (HMW-DNA) preparation. The mouse was killed by cervical dislocation, spleen dissected and single cell suspension prepared in PBS using a glass homogenizer. The agarose-embeded HMW-DNA was prepared as described in detail elsewhere .
The agarose HMW-DNA plugs were subjected to pre-electrophoresis in a CHEF-DR-III apparatus (BioRad) in 1% agarose and 0.5 × TBE buffer for 12 hrs (4 V/cm, 5 s pulse, 14°C). Genomic DNA was partially digested with the mixture of Eco RI endonuclease and Eco RI methylase. The optimal ratio of the enzymes was determined by titration: usually 5–25 units of methylase per 1 unit of endonuclease were employed. DNA fragments were prepared by slight modification of an approach described before . Briefly, DNA fragments were separated from the digested agarose plugs in the CHEF-DR-III in 1% agarose and 0.5 × TBE buffer for 16 hrs (5 V/cm, 0.1 to 40 s pulse, 14°C). Subsequently, three stripes corresponding to fragment size between 150 kb and 200 kb were excised and subjected to another size selection by additional electrophoresis in 0.5 × TBE buffer for 12 hr (5 V/cm, 2.5 to 4.5 s pulse, 14°C). The second size selection effectively removed short fragments while keeping long fragments in the agarose strips. The appropriate fragments were isolated by electroelution and ligated to the Eco RI site of the pBACe3.6 vector . The ligation mixtures were dialyzed on ice in a well created by 0.5% agarose with 1 M glucose for 1 hr. The desalted ligation mixtures were electroporated into E. coli electrocompetent DH10B ElectroMax cells (Invitrogen) by a Gene Pulser apparatus (BioRad) in 0.1 cm cuvette with the following parameters: voltage 1.8 V, impedance 200 Ω, capacitance 25 μF, time constant between 3.5 to 4.5 s. The electroporated cells were diluted in 1 ml SOC medium and incubated in an orbital shaker at 37°C and 200 rpm for 1 hr. The titer of each electroporation reaction was determined by spreading an aliquot on selection agar plates (LB, 20 μg/ml chloramphenicol, 5% sucrose) as described . The remainder containing the primary clones was supplemented with glycerol to the final concentration of 15%, then quickly frozen in liquid nitrogen and stored at -80°C. The frozen stocks of the primary clones were recovered and spread on large selection plates. The colonies were picked with multi-functional robotical system Gene TacTM-G3 (Genomic Solutions) and arrayed in 377 individual 384-well dishes containing LB medium supplied by 7.5% glycerol and 20 μg/ml chloramphenicol. The clones were gridded using the robot on 8 nylon membranes (18 342 unique clones in duplicates per membrane). Afterwards, the bacterial colonies were lysed, their DNA denatured, and crosslinked to the membranes by standard methods .
Estimation of average insert size
One hundred sixty-two clones from the library segment 1 and 236 clones from the segment 2 were randomly picked and grown in 15 ml 2xYT medium (16 g/L tryptone, 10 g/L yeast extract, 5 g/L NaCl, 20 μg/ml chloramphenicol) at 37°C and 300 rpm for 18–20 h. BAC DNA was prepared using a modification of the standard protocol based on alkalic lysis. Briefly, 15 ml of the overnight culture was spun down, the bacterial pellet was resuspended in 300 μl of lysis buffer I (50 mM glucose, 25 mM Tris-HCl with pH = 8.0, 10 mM EDTA), then lysed with 600 μl of freshly made lysis buffer II (0.2 M NaOH, 1% SDS) and precipitated with 450 μl of lysis buffer III (3 M KOAc, pH = 4.8) followed by incubation on ice for 1 hour. The resulted precipitate was spun down in a microfuge for 15 min at maximum speed. The BAC DNA was further precipitated at room temperature with 0.6 volumes of isopropanol for 30 min and centrifuged at maximum speed for 10 min. The pellet was washed with 70% ethanol, air dried shortly and dissolved in 25 μl of TE. The BAC DNA was subjected to Not I (Fermentas) restriction overnight to achieve complete digestion. The reactions were resolved along with the mid range PFG marker I (New England Biolabs, cat # N3551S) in the CHEF-DR-III in 1% agarose and 0.5 × TBE buffer for 16 hrs (5 V/cm, 0.1 to 40 s pulse, 14°C). The insert size was estimated after ethidium bromide staining and the average insert size for both segments of the library was calculated.
Hybridization of high-density colony filters
Seven probes for single-copy mouse genes were used to screen the library on high-density membranes. Six of them were produced by PCR on 50 ng of HMW-DNA isolated from the brain of a female PWD/Ph mouse. The primers were designed using Oligo 6 (MBI) software and the GenBank mRNA sequences (Table 1). PCR was performed at standard conditions: 96°C for 1 min, 30 cycles at 96°C for 10 s, 58°C for 20 s, 72°C for 1 min, and final extension at 72°C for 3 min. The 1085-bp DNA fragment of Tbp (U63933) was derived from the tbp-1942 cDNA clone (kind gift from Dr. Trachtulec) . The resulted amplicons and Tbp fragment were purified and labeled by random priming with [α-32P]dCTP. Hybridization was done in Church buffer  with 15% formamide at 60°C over night. The membranes were washed in 0.2 × SSC, 0.1% SDS buffer at 42°C for 20 min and then in 0.2 × SSC, 0.1% SDS buffer at 60°C for 10 min. The membranes were then autoradiographed for 2 to 10 days, the positive clones identified and recovered from the frozen 384-well plates.
BAC end sequencing and analysis
BAC DNA was prepared from 60 clones as described above and purified using a QIAGEN kit following the manufacturer's instructions. The sequencing was performed using a Big Dye Terminator v3.1 cycle sequencing kit in an ABI 310 instrument (Applied Biosystems) with primers T7 (GGTCGAGCTTGACATTGTAG) and SP6 (GATCCTCCCGAATTGACTAGTG). Each DNA sample was sequenced twice. The sequences from the same BAC end were aligned and manually edited in order to obtain a consensus sequence. BESs were masked for mouse repeats using RepeatMasker  (sensitive settings) and aligned to the mouse genome sequence (mm5 assembly, May 2004, UCSC) [30, 31] using BLAT . The mouse genome sequence had already been soft-masked for repeats by UCSC and BLAT was set to produce all possible alignments (tile size = 10, minimum score = 0, minimum sequence identity = 0). The hits were filtered to keep only those with minimum alignment ratio = 0.8. After manual inspection, a list of BESs mapped to unambiguous positions in the genome was compiled. The corresponding genomic sequences were excised and aligned with the appropriate BESs (unmasked sequence) using SSEARCH  (standard settings). A Perl script was written to process the pair-wise alignments and enumerate the sequence polymorphisms (SNPs, insertions, deletions, etc). The visualization of DNA polymorphisms was made by TeXshade LaTeX package . All intermediate steps were performed using customized Perl scripts and utilities available from UCSC website .
Pavljukova H, Forejt J: Private communication. Mouse News Lett. 1981, 65: 44-
Gregorova S, Forejt J: PWD/Ph and PWK/Ph inbred mouse strains of Mus m. musculus subspecies--a valuable resource of phenotypic variations and genomic polymorphisms. Folia Biol (Praha). 2000, 46: 31-41.
Boursot P, Din W, Anand R, Darviche D, Dod B, VonDeimling F, Talwar GP, Bonhomme F: Origin and radiation of the house mouse: Mitochondrial DNA phylogeny. Journal of Evolutionary Biology. 1996, 9, 391-415.: 391-415. 10.1046/j.1420-9101.1996.9040391.x.
Silver LM: Mouse Genetics Concepts and Applications. Concepts and Applications. 1995, , Oxford University Press, 376-
Forejt J, Ivanyi P: Genetic studies on male sterility of hybrids between laboratory and wild mice (Mus musculus L.). Genet Res. 1974, 24: 189-206.
Forejt J: Hybrid sterility in the mouse. Trends Genet. 1996, 12: 412-417. 10.1016/0168-9525(96)10040-8.
Storchova R, Gregorova S, Buckiova D, Kyselova V, Divina P, Forejt J: Genetic analysis of X-linked hybrid sterility in the house mouse. Mamm Genome. 2004, 15: 515-524.
Fernandes C, Liu L, Paya-Cano JL, Gregorova S, Forejt J, Schalkwyk LC: Behavioral Characterization of Wild Derived Male Mice ( Mus musculus musculus ) of the PWD/Ph Inbred Strain: High Exploration Compared to C57BL/6J. Behav Genet. 2004, 34: 621-630. 10.1007/s10519-004-5589-y.
Perlegen Sciences, Inc. [http://www.perlegen.com/]
Osoegawa K, Zhu B, Shu CL, Ren T, Cao Q, Vessere GM, Lutz MM, Jensen-Seaman MI, Zhao S, de Jong PJ: BAC resources for the rat genome project. Genome Res. 2004, 14: 780-785. 10.1101/gr.2033904.
McPherson JD, Marra M, Hillier L, Waterston RH, Chinwalla A, Wallis J, Sekhon M, Wylie K, Mardis ER, Wilson RK, Fulton R, Kucaba TA, Wagner-McPherson C, Barbazuk WB, Gregory SG, Humphray SJ, French L, Evans RS, Bethel G, Whittaker A, Holden JL, McCann OT, Dunham A, Soderlund C, Scott CE, Bentley DR, Schuler G, Chen HC, Jang W, Green ED, Idol JR, Maduro VV, Montgomery KT, Lee E, Miller A, Emerling S, Kucherlapati, Gibbs R, Scherer S, Gorrell JH, Sodergren E, Clerc-Blankenburg K, Tabor P, Naylor S, Garcia D, de Jong PJ, Catanese JJ, Nowak N, Osoegawa K, Qin S, Rowen L, Madan A, Dors M, Hood L, Trask B, Friedman C, Massa H, Cheung VG, Kirsch IR, Reid T, Yonescu R, Weissenbach J, Bruls T, Heilig R, Branscomb E, Olsen A, Doggett N, Cheng JF, Hawkins T, Myers RM, Shang J, Ramirez L, Schmutz J, Velasquez O, Dixon K, Stone NE, Cox DR, Haussler D, Kent WJ, Furey T, Rogic S, Kennedy S, Jones S, Rosenthal A, Wen G, Schilhabel M, Gloeckner G, Nyakatura G, Siebert R, Schlegelberger B, Korenberg J, Chen XN, Fujiyama A, Hattori M, Toyoda A, Yada T, Park HS, Sakaki Y, Shimizu N, Asakawa S, Kawasaki K, Sasaki T, Shintani A, Shimizu A, Shibuya K, Kudoh J, Minoshima S, Ramser J, Seranski P, Hoff C, Poustka A, Reinhardt R, Lehrach H: A physical map of the human genome. Nature. 2001, 409: 934-941. 10.1038/35057157.
Banfi S, Ledbetter SA, Chinault AC, Zoghbi HY: An easy and rapid method for the detection of chimeric yeast artificial chromosome clones. Nucleic Acids Res. 1992, 20: 1814-
Antoch MP, Song EJ, Chang AM, Vitaterna MH, Zhao Y, Wilsbacher LD, Sangoram AM, King DP, Pinto LH, Takahashi JS: Functional identification of the mouse circadian Clock gene by transgenic BAC rescue. Cell. 1997, 89: 655-667. 10.1016/S0092-8674(00)80246-9.
Abe K, Hazama M, Katoh H, Yamamura K, Suzuki M: Establishment of an efficient BAC transgenesis protocol and its application to functional characterization of the mouse Brachyury locus. Exp Anim. 2004, 53: 311-320. 10.1538/expanim.53.311.
Lee EC, Yu D, Martinez de Velasco J, Tessarollo L, Swing DA, Court DL, Jenkins NA, Copeland NG: A highly efficient Escherichia coli-based chromosome engineering system adapted for recombinogenic targeting and subcloning of BAC DNA. Genomics. 2001, 73: 56-65. 10.1006/geno.2000.6451.
Gong S, Yang XW, Li C, Heintz N: Highly efficient modification of bacterial artificial chromosomes (BACs) using novel shuttle vectors containing the R6Kgamma origin of replication. Genome Res. 2002, 12: 1992-1998. 10.1101/gr.476202.
Clarke L, Carbon J: A colony bank containing synthetic Col El hybrid plasmids representative of the entire E. coli genome. Cell. 1976, 9: 91-99. 10.1016/0092-8674(76)90055-6.
Lindblad-Toh K, Winchester E, Daly MJ, Wang DG, Hirschhorn JN, Laviolette JP, Ardlie K, Reich DE, Robinson E, Sklar P, Shah N, Thomas D, Fan JB, Gingeras T, Warrington J, Patil N, Hudson TJ, Lander ES: Large-scale discovery and genotyping of single-nucleotide polymorphisms in the mouse. Nat Genet. 2000, 24: 381-386. 10.1038/74215.
Wiltshire T, Pletcher MT, Batalov S, Barnes SW, Tarantino LM, Cooke MP, Wu H, Smylie K, Santrosyan A, Copeland NG, Jenkins NA, Kalush F, Mural RJ, Glynne RJ, Kay SA, Adams MD, Fletcher CF: Genome-wide single-nucleotide polymorphism analysis defines haplotype patterns in mouse. Proc Natl Acad Sci U S A. 2003, 100: 3380-3385. 10.1073/pnas.0130101100.
Frazer KA, Wade CM, Hinds DA, Patil N, Cox DR, Daly MJ: Segmental phylogenetic relationships of inbred mouse strains revealed by fine-scale analysis of sequence variation across 4.6 mb of mouse genome. Genome Res. 2004, 14: 1493-1500. 10.1101/gr.2627804.
Abe K, Noguchi H, Tagawa K, Yuzuriha M, Toyoda A, Kojima T, Ezawa K, Saitou N, Hattori M, Sakaki Y, Moriwaki K, Shiroishi T: Contribution of Asian mouse subspecies Mus musculus molossinus to genomic constitution of strain C57BL/6J, as defined by BAC-end sequence-SNP analysis. Genome Res. 2004, 14: 2439-2447. 10.1101/gr.2899304.
Pletcher MT, McClurg P, Batalov S, Su AI, Barnes SW, Lagler E, Korstanje R, Wang X, Nusskern D, Bogue MA, Mural RJ, Paigen B, Wiltshire T: Use of a dense single nucleotide polymorphism map for in silico mapping in the mouse. PLoS Biol. 2004, 2: e393-10.1371/journal.pbio.0020393.
RZPD, Berlin, Germany. [http://www.rzpd.de/]
Osoegawa K, de Jong PJ, Frengen E, Ioannou PA: Construction of Bacterial Arteficial Chromosome (BAC/PAC) Librarries. Current Protocols in Molecular Biology. Edited by: Ausubel FM, Brent R, Kingston RE, Moore DD, Seidman JG, Smith JA and Struhl K. 2001, , John Wiley & Sons, Inc, 1: 5.9.1-5.9.33.
Frengen E, Weichenhan D, Zhao B, Osoegawa K, van Geel M, de Jong PJ: A modular, positive selection bacterial artificial chromosome vector with multiple cloning sites. Genomics. 1999, 58: 250-253. 10.1006/geno.1998.5693.
Ausubel FM, Brent R, Kingston RE, Moore DD, Seidman JG, Smith JA, Struhl K: Current Protocols in Molecular Biology. Edited by: Chanda VB. 2005, , John Wiley & Sons, Inc.
Trachtulec Z, Hamvas RM, Forejt J, Lehrach HR, Vincek V, Klein J: Linkage of TATA-binding protein and proteasome subunit C5 genes in mice and humans reveals synteny conserved between mammals and invertebrates. Genomics. 1997, 44: 1-7. 10.1006/geno.1997.4839.
Church GM, Gilbert W: Genomic sequencing. Proc Natl Acad Sci U S A. 1984, 81: 1991-1995.
Smit AF, Hubley R, Green P: RepeatMasker Open-3.0. [http://www.repeatmasker.org/]
Karolchik D, Baertsch R, Diekhans M, Furey TS, Hinrichs A, Lu YT, Roskin KM, Schwartz M, Sugnet CW, Thomas DJ, Weber RJ, Haussler D, Kent WJ: The UCSC Genome Browser Database. Nucleic Acids Res. 2003, 31: 51-54. 10.1093/nar/gkg129.
Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, Agarwal P, Agarwala R, Ainscough R, Alexandersson M, An P, Antonarakis SE, Attwood J, Baertsch R, Bailey J, Barlow K, Beck S, Berry E, Birren B, Bloom T, Bork P, Botcherby M, Bray N, Brent MR, Brown DG, Brown SD, Bult C, Burton J, Butler J, Campbell RD, Carninci P, Cawley S, Chiaromonte F, Chinwalla AT, Church DM, Clamp M, Clee C, Collins FS, Cook LL, Copley RR, Coulson A, Couronne O, Cuff J, Curwen V, Cutts T, Daly M, David R, Davies J, Delehaunty KD, Deri J, Dermitzakis ET, Dewey C, Dickens NJ, Diekhans M, Dodge S, Dubchak I, Dunn DM, Eddy SR, Elnitski L, Emes RD, Eswara P, Eyras E, Felsenfeld A, Fewell GA, Flicek P, Foley K, Frankel WN, Fulton LA, Fulton RS, Furey TS, Gage D, Gibbs RA, Glusman G, Gnerre S, Goldman N, Goodstadt L, Grafham D, Graves TA, Green ED, Gregory S, Guigo R, Guyer M, Hardison RC, Haussler D, Hayashizaki Y, Hillier LW, Hinrichs A, Hlavina W, Holzer T, Hsu F, Hua A, Hubbard T, Hunt A, Jackson I, Jaffe DB, Johnson LS, Jones M, Jones TA, Joy A, Kamal M, Karlsson EK, Karolchik D, Kasprzyk A, Kawai J, Keibler E, Kells C, Kent WJ, Kirby A, Kolbe DL, Korf I, Kucherlapati RS, Kulbokas EJ, Kulp D, Landers T, Leger JP, Leonard S, Letunic I, Levine R, Li J, Li M, Lloyd C, Lucas S, Ma B, Maglott DR, Mardis ER, Matthews L, Mauceli E, Mayer JH, McCarthy M, McCombie WR, McLaren S, McLay K, McPherson JD, Meldrim J, Meredith B, Mesirov JP, Miller W, Miner TL, Mongin E, Montgomery KT, Morgan M, Mott R, Mullikin JC, Muzny DM, Nash WE, Nelson JO, Nhan MN, Nicol R, Ning Z, Nusbaum C, O'Connor MJ, Okazaki Y, Oliver K, Overton-Larty E, Pachter L, Parra G, Pepin KH, Peterson J, Pevzner P, Plumb R, Pohl CS, Poliakov A, Ponce TC, Ponting CP, Potter S, Quail M, Reymond A, Roe BA, Roskin KM, Rubin EM, Rust AG, Santos R, Sapojnikov V, Schultz B, Schultz J, Schwartz MS, Schwartz S, Scott C, Seaman S, Searle S, Sharpe T, Sheridan A, Shownkeen R, Sims S, Singer JB, Slater G, Smit A, Smith DR, Spencer B, Stabenau A, Stange-Thomann N, Sugnet C, Suyama M, Tesler G, Thompson J, Torrents D, Trevaskis E, Tromp J, Ucla C, Ureta-Vidal A, Vinson JP, Von Niederhausern AC, Wade CM, Wall M, Weber RJ, Weiss RB, Wendl MC, West AP, Wetterstrand K, Wheeler R, Whelan S, Wierzbowski J, Willey D, Williams S, Wilson RK, Winter E, Worley KC, Wyman D, Yang S, Yang SP, Zdobnov EM, Zody MC, Lander ES: Initial sequencing and comparative analysis of the mouse genome. Nature. 2002, 420: 520-562. 10.1038/nature01262.
Kent WJ: BLAT--the BLAST-like alignment tool. Genome Res. 2002, 12: 656-664. 10.1101/gr.229202. Article published online before March 2002.
Pearson WR: Flexible sequence similarity searching with the FASTA3 program package. Methods Mol Biol. 2000, 132: 185-219.
Beitz E: TEXshade: shading and labeling of multiple sequence alignments using LATEX2 epsilon. Bioinformatics. 2000, 16: 135-139. 10.1093/bioinformatics/16.2.135.
UCSC Genome Browser, Jim Kent's Web Page. [http://www.soe.ucsc.edu/~kent/]
We thank Heinz Himmelbauer for invaluable help at the start of the project, Jaroslav Doležel for the access to the core facility of his laboratory, Jaroslav Janda, Michaela Libošvárová and Radka Tusková for help with picking and gridding the library, Vladana Fotopulosová for help with DNA sequencing, and Zdeněk Trachtulec and Šárka Takáčová for reading the manuscript. This work was supported by the Ministry of Education, Youth and Sports of the Czech Republic Grant No. 1M6837805002 – Center for Applied Genomics, by the Grant Agency of the Academy of Sciences Grant AVOZ50520514 and by NIH grant 1R01 HG00318. Jiří Forejt was supported as an International Scholar of the Howard Hughes Medical Institute by Grant 55000306.
PJ constructed the library, participated in the sequence alignment and drafted the manuscript with JF and PD. PD carried out the sequence alignments and the analysis of DNA polymorphism. JF conceived the study, participated in its design and coordination. PJ and JF wrote the manuscript and all authors read and approved its final version.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Jansa, P., Divina, P. & Forejt, J. Construction and characterization of a genomic BAC library for the Mus m. musculus mouse subspecies (PWD/Ph inbred strain). BMC Genomics 6, 161 (2005). https://doi.org/10.1186/1471-2164-6-161