Transcription of a protein-coding gene on B chromosomes of the Siberian roe deer (Capreolus pygargus)
- Vladimir A Trifonov†1Email author,
- Polina V Dementyeva†1,
- Denis M Larkin2,
- Patricia CM O’Brien3,
- Polina L Perelman1,
- Fengtang Yang4,
- Malcolm A Ferguson-Smith3 and
- Alexander S Graphodatsky1
© Trifonov et al.; licensee BioMed Central Ltd. 2013
Received: 9 April 2013
Accepted: 16 July 2013
Published: 6 August 2013
Most eukaryotic species represent stable karyotypes with a particular diploid number. B chromosomes are additional to standard karyotypes and may vary in size, number and morphology even between cells of the same individual. For many years it was generally believed that B chromosomes found in some plant, animal and fungi species lacked active genes. Recently, molecular cytogenetic studies showed the presence of additional copies of protein-coding genes on B chromosomes. However, the transcriptional activity of these genes remained elusive. We studied karyotypes of the Siberian roe deer (Capreolus pygargus) that possess up to 14 B chromosomes to investigate the presence and expression of genes on supernumerary chromosomes.
Here, we describe a 2 Mbp region homologous to cattle chromosome 3 and containing TNNI3K (partial), FPGT, LRRIQ3 and a large gene-sparse segment on B chromosomes of the Siberian roe deer. The presence of the copy of the autosomal region was demonstrated by B-specific cDNA analysis, PCR assisted mapping, cattle bacterial artificial chromosome (BAC) clone localization and quantitative polymerase chain reaction (qPCR). By comparative analysis of B-specific and non-B chromosomal sequences we discovered some B chromosome-specific mutations in protein-coding genes, which further enabled the detection of a FPGT-TNNI3K transcript expressed from duplicated genes located on B chromosomes in roe deer fibroblasts.
Discovery of a large autosomal segment in all B chromosomes of the Siberian roe deer further corroborates the view of an autosomal origin for these elements. Detection of a B-derived transcript in fibroblasts implies that the protein coding sequences located on Bs are not fully inactivated. The origin, evolution and effect on host of B chromosomal genes seem to be similar to autosomal segmental duplications, which reinforces the view that supernumerary chromosomal elements might play an important role in genome evolution.
KeywordsB chromosomes Segmental duplications Gene duplications Karyotype evolution Cervidae
B chromosomes are traditionally classified as a special class of eukaryotic chromosomes different from both autosomes and sex chromosomes . A characteristic trait of B chromosomes is their dispensability; that is, there may be some individuals lacking them in the population . The number of B chromosomes may vary between cells, tissues, and individuals within a population, which is due to their irregular meiotic and mitotic behavior . B chromosomal meiotic drive has been detected in many plant and animal species (reviewed in ) and it plays a key role in the maintaining of B chromosomes through generations. The molecular mechanisms responsible for this drive remained unknown for a long period of time, but recently a non-disjunction of B chromatids accompanied by centromere activity was demonstrated in the male gametophyte of rye . Although it is widely accepted that B chromosomes are highly heterogeneous and may have different properties in different organisms, until recently the concept of the totally heterochromatic and transcriptionally inert nature of B chromosomes has prevailed. Indeed, supernumerary chromosomes of various organisms often contain tandemly arranged  and dispersed  repetitive elements, interstitial telomeric sequences [7, 8], ribosomal DNA clusters [9, 10], histone genes  and so on. The first protein-coding genes were found in B chromosomes of the fungus Nectria haematococca. Recently, a growing body of evidence has accumulated on the presence of coding genes [13–17] and even organellar genome sequences  on B chromosomes of different animal and plant species. This has corroborated the idea of a complex mosaic structure of supernumerary chromosomes, representing autonomous blocks of segmental duplications and heterochromatin . Transcription of ribosomal genes has been detected in some plants  and animals [10, 20]. However, no transcription of B-specific protein-coding genes has been detected in vertebrates so far.
The Siberian roe deer (Capreolus pygargus (CPY), 2n = 70 + 1–14 Bs) is a widely distributed cervid species inhabiting wide areas of Asia from the Volga River to the Pacific coast. The Siberian roe deer is closely related to the European roe deer (Capreolus capreolus, 2n = 70), but these two species do have some morphological and karyological differences and are well resolved in molecular phylogenies [21, 22]. The only derived karyotypic character of the Siberian roe deer is the presence of B chromosomes , the number of supernumerary elements usually varies both between individuals and between cells within the same individual albeit some individuals may possess a stable number of Bs . Recently, we have studied the standard karyotype of the Siberian roe deer by comparative chromosome painting and found it to be highly conserved and similar to that of the cervid ancestor . Although some authors proposed the classification of several subspecies of C. pygargus[26, 27], our studies of the mitochondrial DNA control region in the populations of extant and ancient roe deer failed to reveal any subspecies structure .
Here, we report on three protein-coding genes that are located within the duplicated autosomal regions on the B chromosomes of the Siberian roe deer (Capreolus pygargus) and demonstrate the transcription of a B-specific protein-coding sequence in fibroblast culture.
Generation of Siberian roe deer B chromosome-specific cDNA library (CBCL)
Characteristics of animal samples used in the study
No. of B chromosomes
Yakutia region, Russia
Novosibirsk region, Russia
Cervus elaphus sibiricus
Altai region, Russia
Novosibirsk region, Russia
Altai region, Russia
Altai region, Russia
Altai region, Russia
Leningrad Oblast, Russia
Berlin Zoo, Germany
B-specific cDNA cloning and sequencing
The characteristics of sequenced B chromosome-specific cDNA clones
Name of cDNA clone
Bos taurus(BTA) chromosome
Homology to BTA segment, %
Location of homologous region on BTA3, bp
GenBank accession no.
74,572,819 to 74,573,252
74,666,209 to 74,666,514
74,727,542 to 74,727,927
TNNI3K (exon) FPGT (exons)
74,729,423 to 74,782,291
74,855,578 to 74,856,167
74,856,714 to 74,857,042
74,866,986 to 74,867,372
76,073,726 to 76,074,320
76,143,925 to 76,144,402
76,215,946 to 76,216,107
See Additional file 1
76,256,747 to 76,257,369
76,257,204 to 76,257,458
76,260,507 to 76,261,080
76,261,257 to 76,261,587
76,261,288 to 76,261,613
76,402,615 to 76,403,200
76,402,676 to 76,403,051
76,448,848 to 76,449,833
The mapping of B chromosomes with bovine bacterial artificial chromosome (BAC) clones
To confirm the presence of a BTA3-homologous region on B chromosomes of C. pygargus we performed FISH of bovine BAC clones from the respective genomic segments. Based on the distribution of B chromosome-specific cDNA clones we selected sequenced bovine BAC clones from the CHORI-240 library for FISH mapping (Additional file 2: Table S1). FISH was performed on the metaphase chromosomes of C. pygargus with four Bs (CPY_b) and with eight Bs (CPY_d). As a control we also mapped the same BAC clones via FISH on bovine metaphase chromosomes.
Some BAC clones produced signals of varying intensity on different B chromosomes suggesting molecular heterogeneity among supernumerary elements, for example, the BAC clone CH240-444I8 gave more intensive signals on six out of eight B chromosomes in the CPY_d, whereas all four B chromosomes were painted in equal intensity in the CPY_b. CH240-454D22 gave signals of different intensities on B chromosomes of CPY_b: two out of four B chromosomes consistently showed signals of two to three times stronger intensity than the remaining ones. This could reflect a different degree of regional amplification.
The localization of BAC clones further confirmed the results of cDNA mapping. Assuming that the cDNAs and BACs demonstrate the presence of a continuous segment homologous to BTA3, full copies and parts of some genes must be present in the proximal part of this segment (TNNI3K, LRRIQ3, FPGT, non-protein coding FGD5-AS1, and a pseudogene RNPS1), while the distal part of the segment seems to be gene sparse (it contains only a pseudogene of Arpc2) and represents a large gene desert (Figure 3).
The mapping of B chromosomes by PCR
To refine the boundaries of the amplified regions on B chromosomes we used PCR-assisted mapping on the flow-sorted B chromosome library. We designed 40 primer pairs from BTA3 based on conserved regions of the bovine genome (UCSC genome browser on Cow, October 2011 (Baylor Btau_4.6.1/bosTau7)), taking into account the distribution of B chromosome-specific cDNAs and BAC mapping results (Figure 3). The PCR was carried out using the DOP-PCR library, derived from flow-sorted B chromosomes of the Siberian roe deer, as template. Parallel amplification was conducted with bovine and total roe deer genomic DNAs. The PCR products of appropriate size resulting from the B chromosome-specific library, and genomic DNA of both cattle and Siberian roe deer indicated the presence of these fragments on supernumerary chromosomes of C. pygargus. The absence of several closely located products in the B chromosome-specific library indicated deletions or segment boundaries. We confirmed the specificity of several PCR products by direct sequencing. This method allowed us to define precisely the boundaries of the B chromosome-specific segmental duplication. The PCR-assisted mapping extended the size of the segment on supernumerary chromosomes to approximately 2 Mbp in total, with some possible deletions in the central part of the region (alternatively, these might be due to incomplete coverage of B-specific libraries). To determine the upstream boundary of the segment we used five primer pairs. This boundary was found to be located between positions 74,538 kbp and 74,600 kbp. Using three primer pairs we mapped the downstream boundary of the segment to the interval between 76,468 kbp and 76,502 kbp.
Although the homology of the roe deer B chromosomes to each other has been demonstrated by reverse painting of probes derived from single-copy microdissection experiments (data not shown) and was implicated by the data on BAC localization, we did not have any evidence on the general structural homogeneity of B chromosomes.
The estimation of copy number variations of the regions localized on the B chromosomes of the Siberian roe deer
To estimate the copy numbers of amplified B genes we applied qualitative PCR (qPCR) using conserved parts including the sequence of the fourth exon from FPGT and LRRIQ3 and a promoter region of TNNI3K (Additional file 3: Figure S1). As a single copy reference we used a conserved segment of PTGFR gene located on bovine chromosome 3 outside the region amplified on B chromosomes of the Siberian roe deer (approximately 4.6 Mbp downstream). To achieve equal efficiency of amplification in both cattle and Siberian roe deer we designed primers based on identical short sequences in both species. The short sequences in the Siberian roe deer were obtained by targeted sequencing.
We found that the regions of FPGT and LRRIQ3 genes are present in six copies per diploid genome in CPY_b (four Bs), suggesting a direct correlation between the number of extra copies of the gene fragment and the number of B chromosomes. We showed that the region of TNNI3K gene is present in four copies per diploid genome in CPY_b, suggesting that the copies of TNNI3K gene segment are localized only on two out of four B chromosomes. We found the region of the TNNI3K gene is present in ten copies per diploid genome in the specimen CPY_d (eight Bs) (suggesting each of B chromosomes carries a single copy of this segment). The regions of FPGT and LRRIQ3 genes were found to be represented by eight copies per diploid genome in CPY_d (suggesting that at least two Bs lack the segment). Generally these findings correlate well with the results of cattle BAC clone localization and demonstrate that not all B chromosomes contain full copies of the proposed amplified fragment. Although the number of segment copies depends on the number of B chromosomes, this correlation is not necessarily direct, as some B chromosomes may lack the full fragment. Our results demonstrate some variation in the composition of B chromosomes despite overall similarity and general homology.
DNA sequence comparison of non-B chromosomal and B chromosome-specific fragments
We used targeted sequencing to compare some homologous regions from B chromosomes and non-B chromosomes. We sequenced several regions (from 84 bp to 1,443 bp), comprising the promoter region of the TNNI3K gene (GenBank: JN871289 to JN871291), the first three exons of FPGT (Additional file 1) and the second exon of LRRIQ3 gene (GenBank: JN871294, JN871295) in the library of flow-sorted B chromosomes as well as in the total genomic DNA of the European roe deer and several specimens of the Siberian roe deer (CPY_a-CPY_d) (GenBank: JN871289 to JN871291, JN871294, JN871295). In addition, we sequenced the fourth exon of FPGT gene (GenBank: JN871287, JN871288) and the first exon of TNNI3K gene (Additional file 1) in the European roe deer and CPY_d. We found nucleotide substitutions specific to B chromosomes in two cases. The B chromosomal second exon of LRRIQ3 gene carries a transition (G > A) (valine141 to isoleucine141); the comparison was based on both the roe deer without B chromosomes and cattle, in the DOP-library of B chromosomes (GenBank: JN871295), in B chromosome-specific cDNA and in a fragment from total cDNA derived from CPY_d (obtained sequences were identical to GenBank: JN871295).
The B chromosomal second exon of FPGT carries a transition A > G (lysine83 to arginine83) (Additional file 1). The same B-specific mutation was found in the second exon of FPGT in a cDNA clone (described below).
It is intriguing that the first triplet of FPGT encodes for valine both in B chromosomes and in the genomic DNA of the Siberian roe deer with four and eight B chromosomes (CPY_b, CPY_d) (due to the use of total genomic DNA as template we could not reveal the underrepresented adenosines). A heterozygous condition (valine/methionine) was revealed in the genomes of two individuals of the European roe deer, while we found only methionine in the FPGT first position in the Siberian roe deer without B chromosomes (specimen CPY_a) (Additional file 1). Additional experiments are needed to test if this A1 > G1 transition might be characteristic for a pseudogene in the European roe deer, while the valine triplet (GTG) is associated with the segmental duplication of FPGT located on B chromosomes, since the GTG is known to be a less effective starting codon and may result in the loss of expression or N-terminally truncated proteins [30, 31]. However, it is possible that an alternative upstream exon with a conserved ATG codon may be used. Interestingly the same A/G single nucleotide variation at the first position has been described in human FPGT.
Evidence of FPGT transcription from B chromosomes
The finding of a B chromosome-specific mutation (A > G in the second exon of FPGT) in the processed FPGT-TNNI3K cDNA (GenBank: JN871272) prompted us to study this case in detail. First, we sequenced the second exon of FPGT in CPY_a (without Bs), in red deer (Cervus elaphus), gray brocket deer (Mazama gouazoubira), and Eurasian elk (Alces alces), and found the ‘wild type A’ condition to be conserved (sequencing of B-specific DOP-PCR library, genomic DNA of CPY_d-specific and B-specific cDNA revealed the derived ‘G’ condition). Comparison with many other sequenced vertebrates including opossum, platypus, rooster and anole lizard further corroborated the high conservation of the position in various phylogenetic lineages. To exclude the possibility that the mutation occurred on the autosomal copies of CPY with B chromosomes, we conducted comparative sequencing of autosomal and B chromosomal segments from the same animal. To achieve this, we microdissected separately chromosomes CPY1 and Bs (from the specimens CPY_d and CPY_c) in separate tubes and sequenced the region of interest. Four sequences derived from non-B chromosomes gave ‘a wild-type A’ condition, while three B chromosomal sequences gave ‘a derived G’ condition. To exclude a possible polymorphism among B chromosomes we sequenced ten FPGT exon 2 carrying clones derived from two different flow sorting fractions. All these clones carried the ‘G’ mutation. Based on our discovery of a processed transcript carrying the derived state (in CPY_d) we concluded that this region is transcribed not from autosomes but from B chromosomes.
Siberian roe deer belong to the cetartiodactyl family Cervidae, where B chromosomes have also been described in five Mazama species . It is interesting that both Capreolus and Mazama genera likely possess identical standard karyotypes (2n = 70 + Bs) and have GC-rich small B chromosomes, which may suggest a common ancestry of supernumerary elements. However, our studies on the brown brocket deer (Mazama gouazoubira) did not reveal any homologous unique segments between B chromosomes of these species (unpublished data), suggesting a fairly recent and independent origin of Bs in Capreolus and Mazama. However, there could be a particular predisposition in the genome organization that underlies the higher frequency of B chromosomes appearance in Capreolinae.
Concerning the general polymorphism of B chromosomes, we should emphasize that a Siberian roe deer specimen carrying no B chromosomes (CPY_a) was described here for the first time. The finding of this animal was quite unexpected, considering the many previous karyotyping studies of Siberian roe deer, describing from 1 to 14 copies of B chromosomes per diploid genome [23, 24, 34, 35]. The presence of an animal lacking Bs raised the question of whether this could be a misidentified C. capreolus individual. The analysis of the mitochondrial DNA control region allowed us to unambiguously describe the species as C. pygargus (the haplotype Ns99 is identical to the ancient haplotype DC6 (, GenBank: GU811828). This finding indicates that the amount of B chromosomes in the Siberian roe deer is highly variable between individuals (0 to 14) and that the presence of B chromosomes is not essential for animal physiology and cannot be used as a species discriminating character. It is of interest that some animals were reported as lacking the intraindividual mosaicism of Bs . In this study, we have chosen animals with stable karyotypes from Altai.
The same method of cDNA-selection-based mapping of B chromosomes described here has also been applied to the red fox (Vulpes vulpes). In this species we mapped cDNAs to seven different autosomal regions, and these results were further confirmed by subsequent canine BAC mapping . It seems that the roe deer B chromosomes contain a large segment from a single genomic locus, while canid B chromosomes have a complex structure comprising segments from multiple regions. B chromosomes of the red fox and the raccoon dog share some segments (homologous to canine chromosome 13) while most other segmental duplications (five in the red fox and three in the Chinese raccoon dog) are species specific.
It is worth mentioning that the fragment identified on Siberian roe deer B chromosomes is located near the evolutionary breakpoint on cattle chromosome 3 and contains a large gene desert conserved in both human and cattle. Moreover, we found several large segmental duplications in the cattle genomic region. This might have led to a general genomic instability of the locus and its involvement in B chromosome formation.
The finding of coding genes on B chromosomes of the Siberian roe deer further emphasizes the importance of studying regions previously believed to be non-coding. Representing a variation in number of coding and non-coding sequences these elements may influence many aspects of genome activity. Finding a significant amount of segmental duplications in many vertebrate species suggests that B chromosomes act as autonomously segregating karyotypic elements, carrying some genomic regions that could potentially influence overall genome plasticity and contribute to novel adaptations.
The genes that were identified on the Siberian roe deer B chromosomes, TNNI3K (partial), FPGT and LRRIQ3, may be expressing and influencing metabolism of the cells. Partial sequencing of these genes on B chromosomes and in the Siberian roe deer without any Bs has demonstrated a high sequence conservation and very low number of B-specific substitutions, which might reflect a relatively recent origin of B chromosomes in the Siberian roe deer from a large segment derived from CPY1. The B-specific transcript (GenBank: JN871272) described here represents an alternatively spliced variant, consisting of two exons of FPGT (the first and the second) and a large part of the second exon of TNNI3K. Multiple transcripts composed of different exons due to alternative splicing of FPGT-TNNI3K readthrough transcript have been described in human . However, the combination of exons in the roe deer transcript is different from previously published. TNNI3K seems to be very important in cardiomyogenesis, since its expression promotes the differentiation of cardiomyocytes and enhances cardiac performance [38, 39], while FPGT (the full copy of the gene is present on roe deer Bs) participates in reutilization of l-fucose arising from cellular metabolism of glycoproteins and glycolipids . The locus including both TNNI3K and FPGT genes was found to be responsible for susceptibility to viral myocarditis in mice . Transcription of truncated copies might generally affect the expression of autosomal copies, but more detailed research is needed here. However, these genes maybe partly or completely epigenetically inactivated. Thus, the finding of the B-derived transcript corroborates the hypothesis of transcriptional activity of segmental duplications. Further studies will determine if transcription of B chromosomal copies of the gene LRRIQ3 also takes place. At least we found a putative B chromosome-specific mutation that might be used as a marker of B-derived transcripts, but have not yet excluded a possible case of autosomal heterozygosity of this mutation. It should be emphasized that the transcription of the B chromosomal gene was demonstrated only in fibroblasts and additional experiments are required to test the expression of this gene in other cells and tissues of the organism. It is also possible that cultured fibroblasts might manifest different properties and have different activity of some genes.
Application of novel molecular technologies has allowed the localization of unique sequences, comprising protein-coding genes, on B chromosomes of several vertebrate species [13, 14, 16, 36]. Use of high-throughput sequencing methods to obtain a large amount of information about the structure of whole genomes looks very promising for studying species with supernumerary chromosomes. However, our study shows that it should be supplemented by individual chromosome copy isolation and analysis to differentiate the non-B chromosomal and B-specific components. New data on supernumerary elements of the genome clearly indicate that the role of B chromosomes in evolution and their influence on the host are largely underestimated.
Considering that mammalian B chromosomes carrying additional copies of expressed genes do not have any pronounced adverse effect (at least none described so far), their use for the future development of artificial mammalian chromosomes seems promising: our studies demonstrate that they can be isolated easily from the rest of the genome by chromosome sorting, and studies of other authors show that flow-sorted chromosomes can be inserted into oocytes . Employing newly developed methods of zinc finger nuclease technology  it would be of interest to use animals carrying B chromosomes as models for transgenic studies.
Ethics statement and sample collection
All Siberian roe deer, European roe deer (Leningrad oblast), Eurasian elk and red deer tissue samples were provided by local hunters; animals were not killed for the purpose of this study. Fibroblast tissue cultures of European roe deer (UK), cattle and brown brocket were taken for DNA isolation from the collection of the Comparative genomics group (Cambridge). Table 1 lists the number of specimens and their characteristics.
Cell culture and chromosome preparation
Fibroblast cell lines were established from the ear cartilage tissue of Siberian roe deer individuals with different numbers of B chromosomes as described previously . The cell lines were cultured at 37°C in Dulbecco’s modified Eagle medium (DMEM) and αMEM (Gibco, Invitrogen, Paisley, UK) enriched with 10% fetal bovine serum (Gibco), penicillin (100 units/ml) and streptomycin (100 μg/ml) or kanamycin (100 μg/ml). Before harvest, the cells were treated with 1 μg/ml of ethidium bromide (final concentration) for 2 h and treated with 0.05 mg/ml colchicine (final concentration) for 45 minutes. Preparation of metaphase chromosomes was made according to standard procedures that included a 45-minute hypotonic treatment in 0.56% KCl, fixation in 3:1 methanol/glacial acetic acid, followed by slide preparation and air drying. The number of B chromosomes was analyzed in all C. pygargus samples (at least 50 cells per sample) and we detected no variation in the number of B chromosomes between cells.
B chromosomes of the Siberian Roe deer (CPY_d) have been flow sorted on a FACStar Plus, Becton Dickinson, USA, as described previously [45, 46]. The peak containing B chromosomes was very well resolved and distinct from non-B chromosomal peaks (Additional file 4: Figure S2). We collected B chromosomes into 4 tubes containing 20 μl of water (1,000 chromosomes per tube).
Generation of B chromosome-specific libraries
Chromosome-specific libraries were made by DOP-PCR amplification of flow-sorted chromosomes . DOP-PCR-amplified chromosome-specific DNA was labeled during the secondary PCR by incorporating biotin-16-dUTP (Roche, Switzerland) or Cy3-dUTP (GE Healthcare, UK) . Chromosome specificity of the library was confirmed by FISH (Figure 1).
Selection of B chromosome-specific cDNAs
Total RNA was isolated from several millions of fibroblasts from the sixth passage of the Siberian roe deer primary culture (CPY_d) using the MasterPure RNA purification Kit (Epicentre, UK) according to the manufacturer’s protocol. A double stranded total cDNA library was constructed using the BD SMART™ PCR cDNA Synthesis Kit (BD Biosciences) according to the manufacturer’s protocol. Two rounds of B chromosome-specific cDNA selection were conducted using a method based on previously published selection of hybrids by affinity capture (SHAC) . Briefly, 50 μl of cDNA PCR library (driver) and 1 μl of biotin-labeled B chromosome-specific library mixed with 5 μg of C. pygargus Cot5 DNA (tracer) were denatured for 3 minutes at 95°C. After denaturing the products were adjusted with sodium dodecyl sulfate (SDS) (up to 0.1%) and saline-sodium citrate (SSC) (up to 1 ×) and prehybridized separately at 65°C for 30 minutes to suppress the repetitive elements. Both solutions were mixed together and hybridized overnight at 65°C. The resulting DNA duplexes were captured by streptavidin-coated magnetic beads, Dynabeads® M-280 Streptavidin (Invitrogen, UK) for 20 minutes at 65°C. Magnetic beads were washed twice in 2 × SSC, 0.1% SDS for 5 minutes at 65°C, five times in 0.2 × SSC, 0.1% SDS for 5 minutes at 65°C and once in 0.2 × SSC at 30°C. The cDNAs were released by denaturing the beads in Tris-ethylenediaminetetra-acetic acid (TE) buffer for 5 minutes at 95°C. Selected cDNAs were amplified by PCR with specific primer (provided by BD Bioscience). cDNAs, after the first round of selection were subjected to the second round, to increase the chromosome specificity. The quality of the final CBCLs was controlled by FISH (Figure 2).
Cloning of cDNA fragments
The B-specific cDNA library was cloned using TOPO® TA Cloning® Kit (Invitrogen, UK) according to the manufacturer’s instructions. In all, 50 single colonies were cultured for plasmid DNA isolation; 33 clones with BD primer flanked inserts were selected for sequencing.
Fluorescence in situhybridization
Probes for fluorescence in situ hybridization were labeled either by direct incorporation of labeled dUTP in PCR (with DOP or BD primers) or using a nick translation kit (Invitrogen, UK) following the manufacturer’s instructions (for bovine BACs). Fluorescence in situ hybridization was performed using a standard protocol [49, 50].
Images were captured and processed using the CytoVision Genus system (Applied Imaging, USA) and a Cohu CCD camera mounted on an Olympus BX-60 microscope, or using Videotest 2.0 Image Analysis System (Saint Petersburg, Russia) and a Baumer Optronics CCD Camera mounted on Axioscope 2 plus microscope (Carl Zeiss, Germany).
Isolation of DNA from bovine BAC clones
DNA was isolated from bovine BAC clones (Additional file 2: Table S1) of the genomic BAC library CHORI-240 using the QIAGEN Large-Construct Kit (QIAGEN, UK).
Isolation of DNA from animal tissues and cell cultures
DNA was isolated from all samples listed in Table 1 using DNeasy Blood & Tissue Kit (QIAGEN) according to the manufacturer’s instructions.
PCR-assisted mapping of the flow-sorted CPY B chromosome-specific library
PCR-assisted mapping was conducted using primers designed from cattle chromosome 3 region (UCSC genome browser on Cow, October 2011 (Baylor Btau_4.6.1/bosTau7)) listed in Additional file 5: Table S2. A total of 20 ng DNA of the Siberian roe deer with B chromosomes (CPY_d), flow sorting-derived CPY B chromosome-specific library and cattle genomic DNA were used as templates. The PCR program included: 94°С for 2 minutes; 30 cycles of 94°С (15 s), Х°С (30 s), 72°С (30 s); and a final extension at 72°С for 5 minutes, where Х corresponds to the respective annealing temperature for each pair of primers. PCR products were analyzed on a 1.5% agarose gel.
The fragment of the PTGFR gene amplified using the NBF and NBR primers (Additional file 6: Table S3) on cattle genomic DNA was used as a reference fragment for the real time PCR calibration. The primer sequences were designed from the genomic regions conserved between cattle and the Siberian roe deer. Real time PCR was performed using standard curve methods  using a C1000 Thermal Cycler (BioRad). C1000 manager software was utilized for the analysis of the results.
The copy number of regions was estimated in the genomes of cattle and Siberian roe deer with different numbers of B chromosomes. Each primer pair was tested by non-template PCR to check the primer-dimer formation in the reactions. A final dissociation step was always performed at the end of each PCR to identify the unique and specific amplification of the target sequence. Additional file 6: Table S3 lists the primers used for real time PCR in the present study. The amplification reactions were run three times. Reactions contained 10 μl of DNA (0.2 μg), 5 μl of primers mix (final primer concentration of 1 mM each), and 10 μl of Power SYBR Green PCR Master Mix (Sintol, Moscow, Russia), in a final volume of 25 μl. Amplifications were carried out in a Thermal Cycler C1000 (BioRad). The reactions started with an initial step at 95°C for 5 minutes. The reaction proceeded with 40 cycles of 95°C for 15 s, 60°C for 20 s and 72°C for 30 s. A final dissociation step was always performed to control the amplicons melting curves.
B-specific cDNA clones were sequenced with standard M13 primers using DTCS Quick Start Master Mix and CEQ 2000 sequencer (Beckman-Coulter) in the Department of Veterinary Medicine (Cambridge, UK). PCR products for sequence analysis of FPGT, TNNI3K and LRRIQ3 genes were generated using primers based on cattle chromosome 3 region (UCSC genome browser on Cow, October 2011 (Baylor Btau_4.6.1/bosTau7)) (listed in Additional file 7: Table S4). For comparative analysis of non-B chromosomal and B chromosome-specific copies we amplified the segments from B chromosome-specific libraries, DNA of individuals lacking Bs and microdissected autosomes. Sequencing was carried out in the Interinstitutional center of DNA sequencing at the Siberian Branch of Russian Academy of Sciences (Novosibirsk, Russia) using an ABI3130×1 Genetic Analyzer (Applied Biosystems Inc., CA, USA) with ABI Big Dye kit according to a standard protocol. Sequence analysis was accomplished using Sequence scanner V1.0 software (Applied Biosystems), Blast and Blat alignment tools [52, 53].
We microdissected several copies of the largest autosome (CPY1) and B chromosomes as described earlier  with some modifications: after incubation in the collection drop with proteinase K and SDS, the sample was transferred to water, denatured and adjusted with all reagents necessary for PCR with T2 primers (Additional file 5: Table S2). In all, 40 cycles were performed (controlled with samples containing other autosomes and total genomic DNA) revealing fragments of interest. The fragments were directly sequenced after ExoSap treatment.
The GenBank  accession numbers for DNA sequences are JN871269 to JN871295.
Bacterial artificial chromosome
B chromosome-specific cDNA library
Fluorescence in situ hybridization
Selection of hybrids by affinity capture.
We thank Drs NV Rubtsova (ICG SB RAS), SA Romanenko (IMCB SB RAS) and LS Biltueva (IMCB SB RAS) for help in establishing tissue cultures, AA Sharshov (ICG SB RAS) for tissue samples of the roe deer, Jesus Aguirre-Hernandez (Cambridge University) and Daria Grafodatskaya (Toronto University) for help in cDNA sequencing, and Willem Rens for helpful discussion. This work was funded by a Wellcome Trust Grant to VAT, MAF-S and PCMO’B. VAT, PVD, PLP, NVV and ASG were funded by MCB, RFBR and SB RAS programs. We thank the three anonymous reviewers whose suggestions helped to improve this manuscript.
- White MJD: Animal Cytology and Evolution. 1973, Cambridge, UK: Cambridge University PressGoogle Scholar
- Camacho JPM, Sharbel TF, Beukeboom LW: B-chromosome evolution. Philos Trans R Soc Lond B Biol Sci. 2000, 355: 163-178. 10.1098/rstb.2000.0556.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones N, Houben A: B chromosomes in plants: escapees from the A chromosome genome?. Trends Plant Sci. 2003, 8: 417-423. 10.1016/S1360-1385(03)00187-0.View ArticlePubMedGoogle Scholar
- Banaei-Moghaddam AM, Schubert V, Kumke K, Weiss O, Klemme S, Nagaki K, Macas J, Gonzalez-Sanchez M, Heredia V, Gomez-Revilla D, Gonzalez-Garcia M, Vega JM, Puertas MJ, Houben A: Nondisjunction in favor of a chromosome: the mechanism of rye B chromosome drive during pollen mitosis. Plant Cell. 2012, 24: 4124-4134. 10.1105/tpc.112.105270.PubMed CentralView ArticlePubMedGoogle Scholar
- Potapov VA, Solov’ev VV, Romashchenko AG, Sosnovtsev SV, Ivanov SV: Features of the structure and evolution of complex, tandemly organized Bsp-repeats in the fox genome. I. Structure and internal organization of the BamHI-dimer. Mol Biol (Mosk). 1990, 24: 1649-1665.Google Scholar
- Peppers JA, Wiggins LE, Baker RJ: Nature of B chromosomes in the harvest mouse Reithrodontomys megalotis by fluorescence in situ hybridization (FISH). Chromosome Res. 1997, 5: 475-479. 10.1023/A:1018421114607.View ArticlePubMedGoogle Scholar
- Wurster-Hill DH, Ward OG, Davis BH, Park JP, Moyzis RK, Meyne J: Fragile sites, telomeric DNA sequences, B chromosomes, and DNA content in raccoon dogs, Nyctereutes procyonoides, with comparative notes on foxes, coyote, wolf, and raccoon. Cytogenet Cell Genet. 1988, 49: 278-281. 10.1159/000132677.View ArticlePubMedGoogle Scholar
- Szczerbal I, Switonski M: B chromosomes of the Chinese raccoon dog (Nyctereutes procyonoides procyonoides Gray) contain inactive NOR-like sequences. Caryologia. 2003, 56: 213-216.View ArticleGoogle Scholar
- Stitou S, de la Guardia RD, Jimenez R, Burgos M: Inactive ribosomal cistrons are spread throughout the B chromosomes of Rattus rattus (Rodentia, Muridae). Implications for their origin and evolution. Chromosome Res. 2000, 8: 305-311. 10.1023/A:1009227427575.View ArticlePubMedGoogle Scholar
- Rubtsov NB, Karamysheva TV, Andreenkova OV, Bochkaerev MN, Kartavtseva IV, Roslik GV, Borissov YM: Comparative analysis of micro and macro B chromosomes in the Korean field mouse Apodemus peninsulae (Rodentia, Murinae) performed by chromosome microdissection and FISH. Cytogenet Genome Res. 2004, 106: 289-294. 10.1159/000079301.View ArticlePubMedGoogle Scholar
- Teruel M, Cabrero J, Perfectti F, Camacho JP: B chromosome ancestry revealed by histone genes in the migratory locust. Chromosoma. 2010, 119: 217-225. 10.1007/s00412-009-0251-3.View ArticlePubMedGoogle Scholar
- Miao VP, Covert SF, VanEtten HD: A fungal gene for antibiotic resistance on a dispensable (“B”) chromosome. Science. 1991, 254: 1773-1776. 10.1126/science.1763326.View ArticlePubMedGoogle Scholar
- Graphodatsky AS, Kukekova AV, Yudkin DV, Trifonov VA, Vorobieva NV, Beklemisheva VR, Perelman PL, Graphodatskaya DA, Trut LN, Yang F, Ferguson-Smith MA, Acland GM, Aquirre GD: The proto-oncogene C-KIT maps to canid B-chromosomes. Chromosome Res. 2005, 13: 113-122. 10.1007/s10577-005-7474-9.View ArticlePubMedGoogle Scholar
- Yudkin DV, Trifonov VA, Kukekova AV, Vorobieva NV, Rubtsova NV, Yang F, Acland GM, Ferguson-Smith MA, Graphodatsky AS: Mapping of KIT adjacent sequences on canid autosomes and B chromosomes. Cytogenet Genome Res. 2007, 116: 100-103. 10.1159/000097424.View ArticlePubMedGoogle Scholar
- Lamatsch DK, Trifonov V, Schories S, Epplen JT, Schmid M, Schartl M: Isolation of a cancer-associated microchromosome in the sperm-dependent parthenogen Poecilia formosa. Cytogenet Genome Res. 2011, 135: 135-142. 10.1159/000331271.View ArticlePubMedGoogle Scholar
- Yoshida K, Terai Y, Mizoiri S, Aibara M, Nishihara H, Watanabe M, Kuroiwa A, Hirai H, Hirai Y, Matsuda Y, Okada N: B chromosomes have a functional effect on female sex determination in Lake Victoria cichlid fishes. Plos Genet. 2011, 7: e1002203-10.1371/journal.pgen.1002203.PubMed CentralView ArticlePubMedGoogle Scholar
- Martis MM, Klemme S, Banaei-Moghaddam AM, Blattner FR, Macas J, Schmutzer T, Scholz U, Gundlach H, Wicker T, Šimková H, Novák P, Neumann P, Kubaláková M, Bauer E, Haseneyer G, Fuchs J, Doležel J, Stein N, Mayer KF, Houben A: Selfish supernumerary chromosome reveals its origin as a mosaic of host genome and organellar sequences. Proc Natl Acad Sci USA. 2012, 109: 13343-13346. 10.1073/pnas.1204237109.PubMed CentralView ArticlePubMedGoogle Scholar
- Trifonov VA, Dementyeva PV, Beklemisheva VR, Yudkin DV, Vorobieva NV, Graphodatsky AS: Supernumerary chromosomes, segmental duplications, and evolution. Russian Journal of Genetics. 2010, 46: 1234-1236.View ArticleGoogle Scholar
- Leach CR, Houben A, Field B, Pistrick K, Demidov D, Timmis JN: Molecular evidence for transcription of genes on a B chromosome in Crepis capillaris. Genetics. 2005, 171: 269-278. 10.1534/genetics.105.043273.PubMed CentralView ArticlePubMedGoogle Scholar
- Ruiz-Estevez M, Lopez-Leon MD, Cabrero J, Camacho JPM: B chromosome ribosomal DNA is functional in the grasshopper Eyprepocnemis plorans. PLoS ONE. 2012, 7: e36600-10.1371/journal.pone.0036600.PubMed CentralView ArticlePubMedGoogle Scholar
- Randi E, Pierpaolo M, Danilkin A: Mitochondrial DNA polymorphism in populations of Siberian and European roe deer (Capreolus pygargus and C. capreolus). Heredity. 1998, 80: 429-437. 10.1046/j.1365-2540.1998.00318.x.View ArticlePubMedGoogle Scholar
- Xiao C, Xhang M, Fu Y, Koh H: Mitochondrial DNA distinction of northeastern China roe deer, Siberian roe deer, and European roe deer, to clarify the taxonomic status of Northeastern China roe deer. Biochem Genet. 2007, 45: 93-102. 10.1007/s10528-006-9052-z.View ArticlePubMedGoogle Scholar
- Graphodatsky AS, Sharshov AA, Shutov VV: Karyotypic relationships between Cervidae. Zool Zh. 1990, 69: 101-114.Google Scholar
- Danilkin AA, Baskevich MI: Karyotypes of the roe deer (Capreolus pygargus) in the Far East. Zool Zh. 1987, 66: 314-317.Google Scholar
- Dementyeva PV, Trifonov VA, Kulemzina AI, Graphodatsky AS: Reconstruction of the putative cervidae ancestral karyotype by chromosome painting of Siberian roe deer (Capreolus pygargus) with dromedary probes. Cytogenet Genome Res. 2010, 128: 228-235. 10.1159/000298878.View ArticlePubMedGoogle Scholar
- Danilkin AA: Deer (Cervidae). Mammals of Russia and Adjacent Territories [in Russian]. 1999, Moscow, Russia: GEOS PressGoogle Scholar
- Sheremetyeva IN, Sheremetyev IS: Skull variation in the Siberian roe deer Capreolus pygargus from the Far East: a revision of the distribution of the subspecies. Eur J Wildlife Res. 2008, 54: 557-569. 10.1007/s10344-008-0180-0.View ArticleGoogle Scholar
- Vorobieva NV, Sherbakov DY, Druzhkova AS, Stanyon R, Tsybankov AA, Vasil’ev SK, Shunkov MV, Trifonov VA, Graphodatsky AS: Genotyping of Capreolus pygargus fossil DNA from Denisova cave reveals phylogenetic relationships between ancient and modern populations. PLoS ONE. 2011, 6: e24045-10.1371/journal.pone.0024045.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen-Liu LW, Huang BC, Scalzi JM, Hall BK, Sims KR, Davis LM, Siebert PD, Hozier JC: Selection of hybrids by affinity capture (SHAC): a method for the generation of cDNAs enriched in sequences from a specific chromosome region. Genomics. 1995, 30: 388-392. 10.1006/geno.1995.0038.View ArticlePubMedGoogle Scholar
- Reddy P, Peterkofsky A, McKenney K: Translational efficiency of the Escherichia coli adenylate cyclase gene: mutating the UUG initiation codon to GUG or AUG results in increased gene expression. Proc Natl Acad Sci USA. 1985, 82: 5656-5660. 10.1073/pnas.82.17.5656.PubMed CentralView ArticlePubMedGoogle Scholar
- Révay T, Villagómez DAF, Brewer D, Chenier T, King WA: GTG mutation in the start codon of the androgen receptor gene in a family of horses with 64, XY disorder of sex development. Sex Dev. 2012, 6: 108-116. 10.1159/000334049.View ArticlePubMedGoogle Scholar
- NCBI dbSNP: Reference SNP(refSNP) Cluster Report: rs149824098. http://www.ncbi.nlm.nih.gov/SNP/snp_ref.cgi?rs=149824098,
- Duarte JMB, Jorge W: Chromosomal polymorphism in several populations of deer (genus Mazama) from Brazil. Arch Zootec. 1996, 45: 281-287.Google Scholar
- Sokolov VE, Orlov VN, Chudinovskaya GA, Danilkin AA: Differences in chromosomes between two subspecies Capreolus capreolus capreolus L. and C. c. pygargus Pall. Zool Zhurn. 1978, 57: 1109-1112.Google Scholar
- Danilkin AA: Karyotypes of Eurasian roe deer (Capreolus Gray): a speciation hypothesis. Dokl Akad Nauk. 1985, 285: 1513-1516.Google Scholar
- Duke Becker SE, Thomas R, Trifonov VA, Wayne RK, Graphodatsky AS, Breen M: Anchoring the dog to its relatives reveals new evolutionary breakpoints across 11 species of the Canidae and provides new clues for the role of B chromosomes. Chromosome Res. 2011, 19: 685-708. 10.1007/s10577-011-9233-4.View ArticleGoogle Scholar
- NCBI: FPGT-TNNI3K FPGT-TNNI3K readthrough. http://www.ncbi.nlm.nih.gov/gene?cmd=Retrieve&dopt=Graphics&list_uids=100526835,
- Wheeler FC, Tang H, Marks OA, Hadnott TN, Chu PL, Mao L, Rockman HA, Marchuk DA: TNNI3K modifies disease progression in murine models of cardiomyopathy. PloS Genet. 2009, 5: e1000647-10.1371/journal.pgen.1000647.PubMed CentralView ArticlePubMedGoogle Scholar
- Lai ZF, Chen YZ, Feng LP, Meng XM, Ding JF, Wang LY, Ye J, Li P, Cheng XS, Kitamoto Y, Monzen K, Komuro I, Sakaguchi N, Kim-Mitsuyama S: Overexpression of TNNI3K, a cardiac-specific MAP kinase, promotes P19CL6-derived cardiac myogenesis and prevents myocardial infarction-induced injury. Am J Physiol Heart Circ Physiol. 2008, 295: H708-H716. 10.1152/ajpheart.00252.2008.View ArticlePubMedGoogle Scholar
- Pastuszak I, Ketchum C, Hermanson G, Sjoberg EJ, Drake R, Elbein AD: GDP-L-fucose pyrophosphorylase. Purification, cDNA cloning, and properties of the enzyme. J Biol Chem. 1998, 273: 30165-30174. 10.1074/jbc.273.46.30165.View ArticlePubMedGoogle Scholar
- Wiltshire SA, Leiva-Torres GA, Vidal SM: Quantitative trait locus analysis, pathway analysis, and consomic mapping show genetic variants of TNNI3K, FPGT, or H28 control susceptibility to viral myocarditis. J Immunol. 2011, 186: 6398-6405. 10.4049/jimmunol.1100159.View ArticlePubMedGoogle Scholar
- Co DO, Borowski AH, Leung JD, van der Kaa J, Hengst S, Platenburg GJ, Pieper FR, Perez CF, Jirik FR, Drayer JI: Generation of transgenic mice and germline transmission of a mammalian artificial chromosome introduced into embryos by pronuclear microinjection. Chromosome Res. 2000, 8: 183-191. 10.1023/A:1009206926548.View ArticlePubMedGoogle Scholar
- Le Provost F, Lillico S, Passet B, Young R, Whitelaw B, Vilotte JL: Zinc finger nuclease technology heralds a new era in mammalian transgenesis. Trends Biotechnol. 2010, 28: 134-141. 10.1016/j.tibtech.2009.11.007.View ArticlePubMedGoogle Scholar
- Trifonov VA, Stanyon R, Nesterenko AI, Fu B, Perelman PL, O’Brien PC, Stone G, Rubtsova NV, Houck ML, Robinson TJ, Ferguson-Smith MA, Dobigny G, Graphodatsky AS, Yang F: Multidirectional cross-species painting illuminates the history of karyotypic evolution in Perissodactyla. Chromosome Res. 2008, 16: 89-107. 10.1007/s10577-007-1201-7.View ArticlePubMedGoogle Scholar
- Yang F, Carter NP, Shi L, Ferguson-Smith MA: A comparative study of karyotypes of muntjacs by chromosome painting. Chromosoma. 1995, 103: 642-652. 10.1007/BF00357691.View ArticlePubMedGoogle Scholar
- Yang F, OBrien PCM, Wienberg J, Neitzel H, Lin CC, Ferguson-Smith MA: Chromosomal evolution of the Chinese muntjac (Muntiacus reevesi). Chromosoma. 1997, 106: 37-43. 10.1007/s004120050222.View ArticlePubMedGoogle Scholar
- Telenius H, Pelmear AH, Tunnacliffe A, Carter NP, Behmel A, Ferguson-Smith MA, Nordenskjöld M, Pfragner R, Ponder BA: Cytogenetic analysis by chromosome painting using DOP-PCR amplified flow-sorted chromosomes. Genes Chromosomes Cancer. 1992, 4: 257-263. 10.1002/gcc.2870040311.View ArticlePubMedGoogle Scholar
- Yang F, Fu B, O’Brien PCM, Robinson TJ, Ryder OA, Ferguson-Smith MA: Karyotypic relationships of horses and zebras: results of cross-species chromosome painting. Cytogenet Genome Res. 2003, 102: 235-243. 10.1159/000075755.View ArticlePubMedGoogle Scholar
- Yang F, O’Brien PC, Milne BS, Graphodatsky AS, Solanky N, Trifonov V, Rens W, Sargan D, Ferguson-Smith MA: A complete comparative chromosome map for the dog, red fox, and human and its integration with canine genetic maps. Genomics. 1999, 62: 189-202. 10.1006/geno.1999.5989.View ArticlePubMedGoogle Scholar
- Graphodatsky AS, Yang F, O’Brien PC, Serdukova N, Milne BS, Trifonov V, Ferguson-Smith MA: A comparative chromosome map of the Arctic fox, red fox and dog defined by chromosome painting and high resolution G-banding. Chromosome Res. 2000, 8: 253-263. 10.1023/A:1009217400140.View ArticlePubMedGoogle Scholar
- Rutledge RG, Cote C: Mathematics of quantitative kinetic PCR and the application of standard curves. Nucleic Acids Res. 2003, 31: e93-10.1093/nar/gng093.PubMed CentralView ArticlePubMedGoogle Scholar
- NCBI: BLAST. http://blast.ncbi.nlm.nih.gov/Blast.cgi,
- UCSC Genome Bioinformatics: FAQ. http://genome.ucsc.edu/FAQ/FAQblat.html,
- Yang F, Trifonov VA, Ng BL, Kosyakova N, Carter NP: Generation of paint probes by flow-sorted and microdissected chromosomes. Fluorescence in situ hybridization (FISH) - Application Guide. Edited by: Liehr T. 2009, Berlin Heidelberg: Springer, 35-52.View ArticleGoogle Scholar
- NCBI: GenBank. http://www.ncbi.nlm.nih.gov/Genbank,
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.