- Research article
- Open Access
Transcription of a protein-coding gene on B chromosomes of the Siberian roe deer (Capreolus pygargus)
BMC Biologyvolume 11, Article number: 90 (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.
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)
For chromosome-specific cDNA selection we used a protocol based on the previously published selection of hybrids by affinity capture (SHAC)  technique with some modifications. A biotinylated library derived from degenerate oligonucleotide-primed polymerase chain reaction (DOP-PCR) amplification of flow-sorted B chromosomes was used as ‘target DNA’. We confirmed the B chromosomal origin of the biotinylated library by Fluorescence in situ hybridization (FISH) without Cot5 suppression. Signals were detected on all eight B chromosomes of the Siberian roe deer (CPY_d specimen, Table 1) (Figure 1). In addition, the subcentromeric heterochromatic blocks on all autosomes (but not on X chromosomes) were painted, suggesting the presence of shared heterochromatic sequences on all autosomes and B chromosomes. No signals on autosomal arms were detected. The ‘source’ cDNA was synthesized from total RNA isolated from the CPY_d fibroblast culture (passage 6); the cDNAs of the ‘source’ library were flanked by primers (provided in the library construction kit). We used an excess of CPY Cot5 DNA in hybridization experiments to reduce the presence of the repetitive DNA in the final B-specific cDNA library. The quality of the CBCL (in terms of chromosome specificity and depletion of the repetitive DNA) was confirmed by FISH (Figure 2). Noteworthy were the absence of signals on autosomal pericentric heterochromatic blocks and a high specificity of the resulting B-specific probe.
B-specific cDNA cloning and sequencing
To characterize the CBCL and estimate its content, we cloned the PCR products in plasmid vectors and sequenced 33 random clones. The sequences were aligned against cattle genome assemblies (UCSC genome browser on Cow, October 2011 (Baylor Btau_4.6.1/bosTau7)) using the BLAT search genome alignment tool and the following results were obtained (here we present only matches that had over 80% similarity of 150 bp): eight clones showed similarity to dispersed repetitive elements (mostly long interspersed nuclear elements (LINEs) from the L1 BT family); three clones were homologous to genes from the mitochondrial genome. Of the 22 clones that showed a high similarity to unique parts of the bovine genome (Bos taurus (BTA)), 18 clones were mapped to a 1.8-Mbp region of BTA3 (74.6 to 76.4 Mbp according to (UCSC genome browser on Cow, October 2011 (Baylor Btau_4.6.1/bosTau7))) (Table 2, Additional file 1) (GenBank: JN871269 to JN871285). The remaining four unique clones were homologous to regions on BTA 2, 7, 10 and 28 (data not shown). The cDNAs homologous to the BTA3 region were not randomly distributed across the 1.8 Mbp genomic segment; 2 areas of high match density were observed: 74.6 to 74.9 Mbp (7 matches) and 76.0 to 76.4 Mbp (11 matches), separated by a 1 Mbp region lacking any matches (Figure 3).
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.
All BACs produced specific signals on the Siberian roe deer chromosome 1 and bovine chromosome 3. These results are consistent with our data of cross-species chromosome painting that demonstrated complete homology of chromosomes BTA3 and CPY1 . Six bovine BAC clones (CH240-10H15, CH240-444I8, CH240-493Р4, CH240-515С3, CH240-454D22 and CH240-351I13) produced strong signals on B chromosomes of C. pygargus, confirming the presence of the CPY1 fragment (Figure 4). BAC clones CH240-131I21 and CH240-385G2 did not give any specific signal on B chromosomes, presumably due to the absence of the respective genomic loci on these elements, although this does not exclude the presence of a short flanking fragment from these BACs.
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 characterized the copy number of regions present on B chromosomes (Bs) in Siberian roe deer without B chromosomes (CPY_a), Siberian roe deer with four Bs (CPY_b) and Siberian roe deer with eight Bs (CPY_d) (Figure 5). The qPCR results showed that the regions of TNNI3K and LRRIQ3 genes are present in two copies per diploid genome in CPY_a (similar to PTGFR) whereas the region of FPGT gene is present in four copies per diploid genome in the same animal. This indicates a duplication of the region containing FPGT gene in CPY_a, although no extra chromosomes were found in the specimen.
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.
White MJD: Animal Cytology and Evolution. 1973, Cambridge, UK: Cambridge University Press
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Graphodatsky AS, Sharshov AA, Shutov VV: Karyotypic relationships between Cervidae. Zool Zh. 1990, 69: 101-114.
Danilkin AA, Baskevich MI: Karyotypes of the roe deer (Capreolus pygargus) in the Far East. Zool Zh. 1987, 66: 314-317.
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.
Danilkin AA: Deer (Cervidae). Mammals of Russia and Adjacent Territories [in Russian]. 1999, Moscow, Russia: GEOS Press
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.
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.
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.
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.
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.
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.
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.
Danilkin AA: Karyotypes of Eurasian roe deer (Capreolus Gray): a speciation hypothesis. Dokl Akad Nauk. 1985, 285: 1513-1516.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
NCBI: GenBank. http://www.ncbi.nlm.nih.gov/Genbank,
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.
The authors declare they have no competing interests.
VAT conceived the project and designed the experiments. VAT, PVD and PCMO’B performed the experiments. VAT, PVD, DML, PCMO’B, PLP, FY, MAF and ASG wrote the manuscript. All authors read and approved the final manuscript.
Vladimir A Trifonov, Polina V Dementyeva contributed equally to this work.