- Research article
- Open Access
A genetic linkage map and comparative mapping of the prairie vole (Microtus ochrogaster) genome
BMC Genetics volume 12, Article number: 60 (2011)
The prairie vole (Microtus ochrogaster) is an emerging rodent model for investigating the genetics, evolution and molecular mechanisms of social behavior. Though a karyotype for the prairie vole has been reported and low-resolution comparative cytogenetic analyses have been done in this species, other basic genetic resources for this species, such as a genetic linkage map, are lacking.
Here we report the construction of a genome-wide linkage map of the prairie vole. The linkage map consists of 406 markers that are spaced on average every 7 Mb and span an estimated ~90% of the genome. The sex average length of the linkage map is 1707 cM, which, like other Muroid rodent linkage maps, is on the lower end of the length distribution of linkage maps reported to date for placental mammals. Linkage groups were assigned to 19 out of the 26 prairie vole autosomes as well as the X chromosome. Comparative analyses of the prairie vole linkage map based on the location of 387 Type I markers identified 61 large blocks of synteny with the mouse genome. In addition, the results of the comparative analyses revealed a potential elevated rate of inversions in the prairie vole lineage compared to the laboratory mouse and rat.
A genetic linkage map of the prairie vole has been constructed and represents the fourth genome-wide high-resolution linkage map reported for Muroid rodents and the first for a member of the Arvicolinae sub-family. This resource will advance studies designed to dissect the genetic basis of a variety of social behaviors and other traits in the prairie vole as well as our understanding of genome evolution in the genus Microtus.
Genomes evolve by a number of molecular mechanisms including chromosomal rearrangements . The genomes of rodents, in particular Muroid rodents (super-family: Muroidea) have been shown to have elevated rates of chromosomal rearrangements compared to most other placental mammals (e.g. ). Though the genome sequence of just two Muroid rodents, the laboratory mouse (Mus musculus) and rat (Rattus norvegicus), are currently available [3, 4], interspecies comparisons based on cross-species comparative chromosome painting has provided a low-resolution view of the similarities and differences in genome organization in a number of species within this super-family (e.g. ).
The genus Microtus is comprised of 62 species of voles and is one of the most, if not the most, speciose Muroid genus [6, 7]. Microtus is particularly interesting to study with respect to the process of genome evolution because the rate of speciation in this genus is estimated to be 20-fold higher than the average mammalian lineage  and because their genomes have been associated with rapid rates of evolution [8, 9]. Comparisons of G-banded karyotypes, cross-species chromosome painting, and multi-color banding has yielded a low-resolution view of how the karyotypes within the Microtus genus differ from one another, the likely number and type of large-scale chromosomal rearrangements that have led to those observed differences, as well as the reconstruction of a proposed ancestral Microtus genome [10–15]. In addition, cross-species chromosome painting studies have been used to establish synteny maps between Microtus genomes and that of the mouse [14, 15]. However, a higher resolution map, such as a genetic linkage map, that could reveal additional insights into the rates and patterns of genome evolution within the Microtus genus has not reported.
The North American prairie vole (Microtus ochrogaster) is an emerging model for studying the genetic and molecular bases of social behavior and how it evolves . Breeding colonies of prairie voles have been used in a laboratory research setting for more than 40 years (e.g. ) and a karyotype for the prairie vole was first reported in 1974 . However, despite the history of breeding prairie voles in captivity and interest in the genetic basis of inherited traits in this species [17, 18], no genetic linkage map has been constructed for the prairie vole. Previously we reported a low-resolution comparative cytogenetic map between the prairie vole and the laboratory mouse . With the goal of developing additional genetic resources for this species and for facilitating studies of genome evolution in this lineage, here we describe the construction of a linkage map of the prairie vole (2N = 54) and comparative analyses of this genome with respect to the laboratory mouse (2N = 40) and rat (2N = 42), as well as other Microtus genomes.
Genotyping and SNP features
A total of 624 SNPs were genotyped on a panel of 353 prairie voles. After applying quality control measures and other filters to the genotyping results, the final data set used for the linkage mapping included the genotypes of 431 SNPs from 285 individuals (see Materials and Methods). Most (n = 392) of the prairie vole loci tagged by the filtered set of SNPs could be assigned an orthologous position in the mouse genome, i.e., Type I markers (e.g. ). The average physical spacing of Type I markers along each mouse chromosome varied from 4.8 to 16.7 Mb, for an average spacing across the genome of one marker every 6.9 Mb. Note that 284 of the Type I markers are within genes (Additional File 1).
Linkage groups and chromosomal assignments
Thirty-five linkage groups that included a total of 406 SNPs were identified (Figure 1 and 2, and see Table 1 for a summary of the number of informative meioses). The total length of the sex-averaged genetic map is 1707 cM. The majority of the chromosome and linkage groups (24 of 36) are longer in females than males, and the total length of the female and male autosomal maps are 1885 and 1575 cM, respectively. In order to anchor the linkage groups to chromosomes, we integrated mapping information from the Type I markers in the linkage groups with a previously published low-resolution comparative cytogenetic map of the prairie vole genome . By using the integrated mapping data, we were able to confidently assign 22 linkage groups to 19 of the 26 prairie vole autosomes and one linkage group to the X chromosome (Figure 1 and 2 and Additional File 2). Prairie vole chromosomes 11, 12, 20, 23 and 25 could not be associated with a linkage group due to an absence of cytogenetic data. Linkage groups could not be assigned to chromosomes 3, 9 and 13 due to uninformative or conflicting genetic and cytogenetic data (see Additional File 2). Conversely, thirteen linkage groups could not be associated with a specific chromosome because of a lack of definitive cytogenetic data (see Additional File 2).
The physical linkage, and to a lesser extent the order of genes tend to be conserved between mammalian genomes (e.g. ). Of the Type I markers in the prairie vole linkage map, 89% (345/387) fall within 61 blocks of conserved synteny with mouse that combined span 1.7 Gb (Figure 1, 2 and 3). These blocks of conserved synteny range in length from 46 kb to 143 Mb and are on average 28 Mb. Within the regions of conserved synteny are 91 segments in which the marker order is identical in the prairie vole linkage map and sequenced mouse genome. These segments of conserved gene order range in size from 20 kb to 76 Mb and average 11 Mb (Figure 3).
The differences in linkage and the order of loci between the prairie vole and mouse genomes are a reflection of chromosomal rearrangements that have occurred since the divergence of the lineages leading to these species. Based on a pairwise comparison of the marker order in the prairie vole linkage map to that in the sequenced mouse genome, the GRIMM algorithm  estimated a total of 177 rearrangements between these genomes. Similarly, 198 rearrangements were estimated to have occurred between the prairie vole and rat genomes. In order to reconstruct the history of the chromosomal rearrangements that led to the differences between these rodent genomes, we applied the MGR algorithm  to a four-way comparison between the prairie vole, mouse, rat, and human genomes (Figure 4). Significantly more rearrangements are predicted to have occurred in the prairie vole versus the mouse/rat lineages (Figure 4, χ2 test, p < 0.012). This difference is primarily due to an estimated higher rate of inversions in the prairie vole versus the mouse/rat lineages (χ2 test, p < 0.0003, Figure 4A), whereas no significant difference is detectable with respect to the number of inter-chromosomal rearrangements, i.e. translocations, fissions and fusions (χ2 test, p > 0.64, Figure 4B).
Genome-wide genetic linkage maps including hundreds of markers have been constructed for a number of placental mammals [24–38]. Among those species, Muroid rodents tend to have shorter genetic linkage maps than most other placental mammals, which could be indicative of lower genome-wide recombination rates within this super-family. The length calculated for the sex-averaged prairie vole linkage map is 1707 cM. The prairie vole linkage map is therefore longer than those of the previously sampled Muroid rodents, i.e, 1111-1647 cM [24, 27, 28], but shorter than the map lengths reported in 13 out of 15 other placental mammals, i.e. 2048-4370 cM [25–36], the exceptions being the American mink and silver fox, which have map lengths of 1340 cM and 1480 cM, respectively [37, 38]. Based on the syntenic locations of the Type I markers in the genome of the laboratory mouse, we estimate the linkage map reported here spans ~89% of the prairie vole genome. Thus we believe our observed map length for the prairie vole is a fairly accurate reflection of the true genetic map length. In addition, the length of the prairie vole linkage map is consistent with the previous observation that Muroid rodent genomes tend to have shorter genetic linkage maps than most placental mammals that have been sampled to date. This observation is likely because rodents tend to possess a larger number of acrocentric chromosomes resulting in fewer cross-over events relative to higher order mammals .
Our comparative analyses revealed 61 blocks of conserved synteny between the prairie vole and mouse that together span the majority of those genomes. Not surprisingly, the largest block of conserved synteny was associated with the X chromosome. Other chromosomes, in particular MMU3, 4, 6, 7, 9, 12, 14 and 16, were also associated with blocks of conserved synteny > 50 Mb corresponding to prairie vole chromosomes 21, LG5/10, 14, 22, 5, 1, 15 and 17, respectively, including an 76 Mb segment of MMU14 in which the order of the markers was the same as on prairie vole chromosome 17. We further estimate that 177 to 201 chromosomal rearrangements have occurred in the 24 million years since the most recent common ancestor of the prairie vole and laboratory mouse . It should be noted, however, that our comparative analyses did not have unlimited resolution, as would be the case when comparing a pair of sequenced and assembled genomes. As such, there are most certainly other large chromosomal rearrangements between the genomes of the prairie vole and laboratory mouse that would not have been detected in our study. For example, the number of chromosomal rearrangements between the genomes of the laboratory mouse and rat estimated by GRIMM and MGR based solely on our Type I markers included in the prairie vole linkage map are 71 to 81, respectively, and are thus lower than the estimated 100 to 107 rearrangements between those rodents that were based on alignments of the genome sequence assemblies . Therefore, with the caveat that the prairie vole linkage map is likely to contain some errors with respect to the true maker order in the genome (see Discussion below), our estimates of the number of chromosomal rearrangements between the genomes of the prairie vole and laboratory mouse and rat are likely reasonable underestimates of the approximate number of large chromosomal rearrangements that distinguish the karyotypes of these species.
The rate of chromosomal rearrangements varies across the mammalian phylogeny . From our comparison of the prairie vole linkage map to that of sequenced mouse, rat and human genomes, we were able to detect a significant difference in the rate of chromosomal rearrangements in the prairie vole lineage compared to that of the laboratory mouse/rat (Murinae) lineage. In particular, the number of inversions in the prairie vole lineage was much higher than that observed in the Murinae lineage (Figure 4). This observation could therefore be a true indicator of an increase in the rate of inversions in the prairie vole lineage, a relative decrease in the rate of inversions in the Murinae lineage, or both. Alternatively, while Microtus genomes have previously been reported to be rapidly evolving [8, 9], it is possible that the number of inversions estimated in this study is artificially inflated due to the limited power to accurately order closely linked loci with the pedigree-based linkage map. Future comparative analyses of the pending genome assembly of the prairie vole http://www.genome.gov/10002154 along with the genome of another member of the Cricetidae family, the deer mouse (Peromyscus maniculatus), should be able to address this question further and better resolve when shifts in the rate of genome evolution may have occurred.
Comparative analyses of G-banded karyotypes from a number of Arvicolid rodents, including the prairie vole and other members of the Microtus genus, identified derived characteristics of the prairie vole genome . More recently, comparative analysis of several Microtus genomes other than the prairie vole by cross-species chromosome painting has led to the reconstruction of a putative ancestral Microtus karyotype . While we cannot directly compare the results of our comparative mapping data with those previously reported due to differences in methodology, an attempt to integrate and synthesize the genome mapping results of the prairie vole and other Microtus genomes with respect to our higher-resolution comparative mapping results is of potential value. For example, by using synteny maps with the laboratory mouse reported here for the prairie vole and previously for other voles as a common point of reference , we predict that the putative ancestral Microtus autosomes 2, 5, 6, 7, 12, 13, 14, 19, 21, 23 and 24 proposed in  are likely orthologous to prairie vole chromosomes 2, 6, 7, 8, 18, 16, 17, 15, LG2, 22 and 21, respectively. In addition, prairie vole chromosomes 1 and 4 likely represent derived fissions of ancestral chromosomes 12 and 16, and 11 and 18, respectively. Though this strategy could not assign all prairie vole chromosomes to an ancestral Microtus chromosome, or vice versa, it does suggest that at the macro-level the prairie vole genome differs from the ancestral Microtus karyotype by the aforementioned pair of fission events, and presumably at least two fusion events which would be needed to maintain the diploid karyotype of 2n = 54 present in both the prairie vole and ancestral state. However, it is important to note that future improvement of the vole linkage map may alter some of the comparisons we report here. Future comparative mapping studies of other Microtus genomes that can now be undertaken by leveraging the genomic resources being developed for the prairie vole will provide further insights into the rates and mechanisms by which genomes within this genus have evolved.
The genetic linkage map of the prairie vole will provide an important resource towards our understanding of genome evolution in the genus Microtus. Further, it will provide a beneficial resource for furthering our understanding of the genetic basis of social behaviors and other traits while giving insight into how these traits evolve.
SNP discovery, selection and genotyping
SNPs were discovered by either re-sequencing loci corresponding to previously reported bacterial artificial chromosome (BAC)-end sequences from the prairie vole CHORI-232 library , or using data from an ongoing unpublished transcriptome sequencing project (see GenBank SRA Accessions: SRX018685, SRX018510, SRX018516, SRX018515, SRX018511, SRX018514, SRX018513, SRX018512). In the case of the BAC-end re-sequencing, loci-specific M13F and M13R tailed primers were designed with PRIMER3  and used to PCR amplify and then directly sequence (Sanger) each locus in three individuals from our local colony of prairie voles. Variant calls from PolyPhred  were manually validated and a prioritized list of SNPs was selected for genotyping. In particular, SNPs linked to loci that shared homology with the mouse genome were preferentially selected for genotyping using the SNPStream® (Beckman Coulter) platform. Transcriptome sequencing data generated with the Roche 454 sequencing platform  was scanned for SNPs by mapping the individual reads back to the contigs in which they were assembled by Newbler (Roche) using the GS Reference Mapper algorithm (Roche). The candidate SNPs were then filtered to select only those with > 30% (minimum of 3) reads supporting the alternative allele and were > 60-bp from the nearest flanking variable site, repetitive sequence (RepeatMasker, http://repeatmasker.org/), or predicted intron-exon boundary defined based on alignment of the prairie vole transcripts to the mouse genome.
The chromosome locations for the mouse orthologs of the 2773 prairie vole transcripts in which at least one SNP was identified were used as a proxy to select a uniformly-spaced and genome-wide set of 384 SNPs for genotyping on the GoldenGate (Illumina) platform. The names, sequences, and other accessory information for each genotyped locus are described in Additional File 1.
A total of 353 prairie voles from our local colony were genotyped. After excluding individuals with questionable parentage, i.e. children for which genotypes at > 3% of the makers were inconsistent with the parental genotypes, two pedigrees were selected for constructing the linkage map: a 3-4 generation pedigree derived from interbreeding descendents of three founder breeding pairs (n = 263 individuals, Additional File 3), and a single nuclear family (n = 22 individuals, Additional File 4). Euthanasia and collection of liver samples for DNA extractions were performed as per guidelines that were reviewed and approved by the Emory Institutional Animal Care and Use Committee and were conducted in accordance with the Guide for Care and Use of Laboratory Animals published by the National Research Council.
Genotypes of X-linked loci were re-coded to account for hemizygosity in males and markers that had high failure rates or indicated the presence of a paralogous sequence variant, i.e., all heterozygous genotypes (n = 47), were monomorphic (n = 109), or associated with an inferred genotype error rate in excess of 7% (n = 37) were excluded from the linkage analyses. Inconsistent genotypes, which represented 0.5% (n = 638) of the genotype calls for the remaining 438 markers, were re-coded as no data. The final data set used for the linkage analyses included 122,071 genotype calls and 2,759 missing data points.
Linkage map construction
A modified version of CRI-MAP , v2.503 (kindly provided by J. F. Maddox), was used to estimate two-point LOD scores between all markers. The markers were then placed into linkage groups based on a LOD threshold of ≥ 2.5, but note that only 5 markers in the final linkage groups had a maximum two-point LOD score of < 3.0. When necessary, the initial linkage groups were further subdivided by increasing the LOD threshold in order to yield groups of markers that limited the number of potentially spurious linkages, i.e., those that were at odds with the published comparative cytogenetic map . The order of markers within an individual linkage group was determined by establishing a framework map with the BUILD option and then adding and iteratively re-ordering markers with the FLIPS3/4 option until no order with a better likelihood score could be found. Due to the size of the first pedigree this procedure proved to be computationally impractical for prairie vole chromosomes 1, 2, 4 and 16, and linkage groups 1, 3 and 4. To overcome this problem we split the large pedigree into smaller overlapping sub-families, optimized the marker order as above, and when possible then calculated the likelihood map scores and cM positions for the best map using the original pedigree structure.
An inherent methodological limitation of the likelihood-based method used by CRI-MAP is that it does not necessarily explore all possible marker orders, and thus the marker order found to have the best likelihood score starting from a given framework map may not represent the true optimal order of the markers. To address that limitation, when the linkage groups included more than 2 loci that mapped to the same orthologous mouse chromosome and the order of the markers was not the same as that in mouse we also calculated the likelihood score of each group of markers assuming an order equivalent to that of the mouse. For prairie vole chromsomes 6, 10, 15 and 21, the marker order predicted by the orthologous positions in the mouse yielded a map with a better likelihood score than the synteny naïve method, in which case the mouse order was used as starting point to improve the ordering of the prairie vole markers with the FLIPS option. Additionally, several markers were uniquely placed to the same position on the map Pair-wise LOD scores for these loci are indicated in Additional File 5. The linkage maps were plotted using MapChart2.2 .
Integration of linkage and cytogenetic maps
The prairie vole-mouse comparative cytogenetic map described in  was integrated with the linkage mapping data based on the orthologous position of the markers in the mouse genome. Specifically, the orthologous mouse positions of the FISH mapped BAC clones in  were converted from mm8 to mm9 coordinates. Both the cytogenetic and linkage markers were then sorted based on their position/order along the prairie vole chromosomes and orthologous position in the mouse genome (Additional File 2).
Loci in the prairie vole linkage map with a known orthologous position in the mouse genome (Additional File 1) were used to construct a comparative map with the mouse (genome assembly version mm9). The corresponding orthologous position and orientation in the human (genome assembly version hg19) and rat genomes (genome assembly version rn4) of those loci was also inferred using the UCSC Genome Browser LiftOver tool http://genome.ucsc.edu/hgLiftOver. GRIMM  was used to infer the orientation of the loci in prairie vole linkage map via a comparison to the mouse genome. MGR  was then applied to the four-way comparison between the marker order in prairie vole linkage map and the mouse, rat and human genomes to estimate the number and types of rearrangements that had occurred across the phylogeny. Note that the relative order of prairie vole markers that mapped to the identical location in the linkage map were assigned randomly, or when applicable by maximizing conserved synteny with the mouse (Supplementary Tables 2 and 3). A conserved block of synteny was defined as 2 or more consecutive markers (ignoring markers without a known orthologous position in the mouse genome) that were present as uninterrupted strings of loci, independent of order, in the prairie vole linkage map and mouse genome. The size of a conserved syntenic block was defined as the distance between the first and last loci within the block based on the mouse chromosomal position and order.
Coghlan A, Eichler EE, Oliver SG, Paterson AH, Stein L: Chromosome evolution in eukaryotes: a multi-kingdom perspective. Trends in Genetics. 2005, 21: 673-682. 10.1016/j.tig.2005.09.009.
Murphy WJ, Larkin DM, Everts-van der Wind A, Bourque G, Tesler G, Auvil L, Beever JE, Chowdhary BP, Galibert F, Gatzke L: Dynamics of mammalian chromosome evolution inferred from multispecies comparative maps. Science. 2005, 309 (5734): 613-617. 10.1126/science.1111387.
Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, Agarwal P, Agarwala R, Ainscough R, Alexandersson M, An P: Initial sequencing and comparative analysis of the mouse genome. Nature. 2002, 420 (6915): 520-562. 10.1038/nature01262.
Gibbs RA, Weinstock GM, Metzker ML, Muzny DM, Sodergren EJ, Scherer S, Scott G, Steffen D, Worley KC, Burch PE: Genome sequence of the Brown Norway rat yields insights into mammalian evolution. Nature. 2004, 428 (6982): 493-521.
Romanenko SA, Perelman PL, Serdukova NA, Trifonov VA, Biltueva LS, Wang J, Li T, Nie W, O'Brien PC, Volobouev VT: Reciprocal chromosome painting between three laboratory rodent species. Mamm Genome. 2006, 17 (12): 1183-1192. 10.1007/s00335-006-0081-z.
Nowak RM: Walker's mammals of the world. 1999, Baltimore: The Johns Hopkins University Press, 6
Wilson DE, Reeder DM, eds: Mammal Species of the World. A Taxonomic and Geographic Reference. 2005, Baltimore: Johns Hopkins University Press, 3
Triant DA, Dewoody JA: Accelerated molecular evolution in Microtus (Rodentia) as assessed via complete mitochondrial genome sequences. Genetica. 2006, 128 (1-3): 95-108. 10.1007/s10709-005-5538-6.
Maruyama T, Imai HT: Evolutionary rate of the mammalian karyotype. J Theor Biol. 1981, 90 (1): 111-121. 10.1016/0022-5193(81)90125-9.
Modi W: Phylogenetic analyses of the chromosomal banding patterns among the Neartic Arviclidae. Systematic Zoology. 1987, 36: 109-136. 10.2307/2413264.
Li T, Wang J, Su W, Yang F: Chromosomal mechanisms underlying the karyotype evolution of the oriental voles (Muridae, Eothenomys). Cytogenet Genome Res. 2006, 114 (1): 50-55. 10.1159/000091928.
Lemskaya NA, Romanenko SA, Golenishchev FN, Rubtsova NV, Sablina OV, Serdukova NA, O'Brien PC, Fu B, Yigit N, Ferguson-Smith MA: Chromosomal evolution of Arvicolinae (Cricetidae, Rodentia). III. Karyotype relationships of ten Microtus species. Chromosome Res. 2010, 18 (4): 459-471. 10.1007/s10577-010-9124-0.
Romanenko SA, Sitnikova NA, Serdukova NA, Perelman PL, Rubtsova NV, Bakloushinskaya IY, Lyapunova EA, Just W, Ferguson-Smith MA, Yang F: Chromosomal evolution of Arvicolinae (Cricetidae, Rodentia). II. The genome homology of two mole voles (genus Ellobius), the field vole and golden hamster revealed by comparative chromosome painting. Chromosome Res. 2007, 15 (7): 891-897. 10.1007/s10577-007-1171-9.
Sitnikova NA, Romanenko SA, O'Brien PC, Perelman PL, Fu B, Rubtsova NV, Serdukova NA, Golenishchev FN, Trifonov VA, Ferguson-Smith MA: Chromosomal evolution of Arvicolinae (Cricetidae, Rodentia). I. The genome homology of tundra vole, field vole, mouse and golden hamster revealed by comparative chromosome painting. Chromosome Res. 2007, 15 (4): 447-456. 10.1007/s10577-007-1137-y.
Trifonov VA, Kosyakova N, Romanenko SA, Stanyon R, Graphodatsky AS, Liehr T: New insights into the karyotypic evolution in muroid rodents revealed by multicolor banding applying murine probes. Chromosome Res. 2010, 18 (2): 265-275. 10.1007/s10577-010-9110-6.
McGraw LA, Young LJ: The prairie vole: an emerging model organism for understanding the social brain. Trends Neurosci. 2010, 33 (2): 103-109. 10.1016/j.tins.2009.11.006.
Richmond M, Conaway CH: Management, breeding, and reproductive performance of the vole, Microtus ochrogaster, in a laboratory colony. Lab Anim Care. 1969, 19 (1): 80-87.
Hartke GT, Leipold HW, Huston K, Cook JE, Saperstein G: Three mutations and the karyotype of the prairie vole. White spotting, polydipsia, and muscular dystrophy in Microtus ochrogaster. J Hered. 1974, 65 (5): 301-307.
McGraw LA, Davis JK, Lowman JJ, ten Hallers BF, Koriabine M, Young LJ, de Jong PJ, Rudd MK, Thomas JW: Development of genomic resources for the prairie vole (Microtus ochrogaster): construction of a BAC library and vole-mouse comparative cytogenetic map. BMC Genomics. 2010, 11: 70-
O'Brien SJ: Mammalian genome mapping: lessons and prospects. Curr Opin Genet Develop. 1991, 1: 105-111. 10.1016/0959-437X(91)80050-V.
O'Brien SJ, Menotti-Raymond M, Murphy WJ, Nash WG, Wienberg J, Stanyon R, Copeland NG, Jenkins NA, Womack J, Graves JA: The promise of comparative genomics in mammals. Science. 1999, 286 (5439): 458-481. 10.1126/science.286.5439.458.
Tesler G: GRIMM: genome rearrangements web server. Bioinformatics. 2002, 18 (3): 492-493. 10.1093/bioinformatics/18.3.492.
Lin CH, Zhao H, Lowcay SH, Shahab A, Bourque G: webMGR: an online tool for the multiple genome rearrangement problem. Bioinformatics. 2010, 26 (3): 408-410. 10.1093/bioinformatics/btp689.
Okuizumi H, Ohsumi T, Sakaki N, Imoto H, Mizuno Y, Hanami T, Yamashita H, Kamiya M, Takada S, Kitamura A: Linkage map of Syrian hamster with restriction landmark genomic scanning. Mamm Genome. 1997, 8 (2): 121-128. 10.1007/s003359900370.
Wong AK, Ruhe AL, Dumont BL, Robertson KR, Guerrero G, Shull SM, Ziegle JS, Millon LV, Broman KW, Payseur BA: A comprehensive linkage map of the dog genome. Genetics. 2010, 184 (2): 595-605. 10.1534/genetics.109.106831.
Menotti-Raymond M, David VA, Schaffer AA, Tomlin JF, Eizirik E, Phillip C, Wells D, Pontius JU, Hannah SS, O'Brien SJ: An autosomal genetic linkage map of the domestic cat, Felis silvestris catus. Genomics. 2009, 93 (4): 305-313. 10.1016/j.ygeno.2008.11.004.
Brown DM, Matise TC, Koike G, Simon JS, Winer ES, Zangen S, McLaughlin MG, Atkinson OS, Hudson JR, Chakravarti A: An integrated genetic linkage map of the laboratory rat. Mammalian Genome. 1998, 9: 521-530. 10.1007/s003359900812.
Cox A, Ackert-Bicknell CL, Dumont BL, Ding Y, Bell JT, Brockmann GA, Wergedal JE, Bult C, Paigen B, Flint J: A new standard genetic map for the laboratory mouse. Genetics. 2009, 182 (4): 1335-1344. 10.1534/genetics.109.105486.
Poissant J, Hogg JT, Davis CS, Miller JM, Maddox JF, Coltman DW: Genetic linkage map of a wild genome: genomic structure, recombination and sexual dimorphism in bighorn sheep. BMC Genomics. 2010, 11: 524-10.1186/1471-2164-11-524.
Rogers J, Mahaney MC, Witte SM, Nair S, Newman D, Wedel S, Rodriguez LA, Rice KS, Slifer SH, Perelygin A: A genetic linkage map of the baboon (Papio hamadryas) genome based on human microsatellite polymorphisms. Genomics. 2000, 67: 237-247. 10.1006/geno.2000.6245.
Rogers J, Garcia R, Shelledy W, Kaplan J, Arya A, Johnson Z, Bergstrom M, Novakowski L, Nair P, Vinson A: An initial genetic linkage map of the rhesus macaque (Macaca mulatta) genome using human microsatellite loci. Genomics. 2006, 87 (1): 30-38. 10.1016/j.ygeno.2005.10.004.
Jasinska AJ, Service S, Levinson M, Slaten E, Lee O, Sobel E, Fairbanks LA, Bailey JN, Jorgensen MJ, Breidenthal SE: A genetic linkage map of the vervet monkey (Chlorocebus aethiops sabaeus). Mamm Genome. 2007, 18 (5): 347-360. 10.1007/s00335-007-9026-4.
Matise TC, Chen F, Chen W, De La Vega FM, Hansen M, He C, Hyland FC, Kennedy GC, Kong X, Murray SS: A second-generation combined linkage physical map of the human genome. Genome Res. 2007, 17 (12): 1783-1786. 10.1101/gr.7156307.
Ihara N, Takasuga A, Mizoshita K, Takeda H, Sugimoto M, Mizoguchi Y, Hirano T, Itoh T, Watanabe T, Reed KM: A comprehensive genetic map of the cattle genome based on 3802 microsatellites. Genome Res. 2004, 14 (10A): 1987-1998. 10.1101/gr.2741704.
Swinburne JE, Boursnell M, Hill G, Pettitt L, Allen T, Chowdhary B, Hasegawa T, Kurosawa M, Leeb T, Mashima S: Single linkage group per chromosome genetic linkage map for the horse, based on two three-generation, full-sibling, crossbred horse reference families. Genomics. 2006, 87 (1): 1-29. 10.1016/j.ygeno.2005.09.001.
Rohrer GA, Alexander LJ, Hu Z, Smith TP, Keele JW, Beattie CW: A comprehensive map of the porcine genome. Genome Res. 1996, 6 (5): 371-391. 10.1101/gr.6.5.371.
Kukekova AV, Trut LN, Oskina IN, Johnson JL, Temnykh SV, Kharlamova AV, Shepeleva DV, Gulievich RG, Shikhevich SG, Graphodatsky AS: A meiotic linkage map of the silver fox, aligned and compared to the canine genome. Genome Res. 2007, 17 (3): 387-399. 10.1101/gr.5893307.
Anistoroaei R, Ansari S, Farid A, Benkel B, Karlskov-Mortensen P, Christensen K: An extended anchored linkage map and virtual mapping for the American mink genome based on homology to human and dog. Genomics. 2009, 94 (3): 204-210. 10.1016/j.ygeno.2009.05.014.
Coop G, Przeworski M: An evolutionary view of human recombination. Nature Reviews Genetics. 2007, 8: 23-34.
Steppan S, Adkins R, Anderson J: Phylogeny and divergence-date estimates of rapid radiations in muroid rodents based on multiple nuclear genes. Syst Biol. 2004, 53 (4): 533-553. 10.1080/10635150490468701.
Bourque G, Pevzner PA, Tesler G: Reconstructing the genomic architecture of ancestral mammals: lessons from human, mouse, and rat genomes. Genome Res. 2004, 14 (4): 507-516. 10.1101/gr.1975204.
Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz SA, Misener S. 2000, Totowa: Humana Press, 365-386.
Nickerson DA, Tobe VO, Taylor SL: PolyPhred: automating the detection and genotyping of single nucleotide substitutions using fluorescence-based resequencing. Nucleic Acids Research. 1997, 25: 2745-2751. 10.1093/nar/25.14.2745.
Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437 (7057): 376-380.
Green P, Falls K, Crooks SM: Documentation for CRI-MAP, version 2.4. 1990
Voorrips RE: MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002, 93 (1): 77-78. 10.1093/jhered/93.1.77.
Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics. Genome Res. 2009, 19 (9): 1639-1645. 10.1101/gr.092759.109.
The authors wish to acknowledge Lanikea King for his assistance with the re-sequencing and SNP discovery, Greg K. Tharp for computational support, the Emory University School of Medicine Biomarker core for performing the genotyping, and Jillian F. Maddox (University of Melbourne, Department of Veterinary Science) for the updated version of CRI-MAP. This work was supported by National Institutes of Health grants 1R21MH082225 to JWT, JKD, and LJY and 1F32MH079661 to LAM.
JKD performed the re-sequencing of the BES for SNP discovery. LAM, LJY, and JWT conceived and oversaw the project. LAM and JWT performed the analyses and wrote the manuscript. All authors approved the final version of the manuscript.
Electronic supplementary material
Additional file 1:Supplementary Table 1. Marker information. Locus name refers to either the BAC-end locus or the orthologos mouse gene name. SNP names were asssigned similarly where names beginning in "FI" correspond to vole BAC-end loci names and names beginning in "NM" correspond to the accession number of the orthologos mouse gene. (XLS 282 KB)
Additional file 2:Supplementary Table 2. Integration of the prairie vole cytogenetic and linkage maps. Markers that did not have an orthologus position in the mouse genome were omitted from this table. Concordant markers were those in which the cytogenetic and genetic linkage maps agreed with respect to blocks of prairie vole-mouse synteny. (XLS 107 KB)
Additional file 3:FamilyPedigree1. Pedigree of large multi-generational family used for constructing the linkage map. (PDF 193 KB)
Additional file 5:Supplementary Table 3. Pair-wise LOD scores for markers that were uniquely placed to the same position on the map. (XLS 24 KB)
About this article
Cite this article
McGraw, L.A., Davis, J.K., Young, L.J. et al. A genetic linkage map and comparative mapping of the prairie vole (Microtus ochrogaster) genome. BMC Genet 12, 60 (2011). https://doi.org/10.1186/1471-2156-12-60
- Linkage Group
- Bacterial Artificial Chromosome
- Chromosomal Rearrangement
- Mouse Genome
- Placental Mammal