Roles of M1 Macrophages and Their Extracellular Vesicles in Cancer Therapy.

Tumor-associated macrophages (TAMs) are inflammatory cells that are important components of the tumor microenvironment. TAMs are functionally heterogeneous and divided into two main subpopulations with distinct and opposite functions: M1 and M2 macrophages. The secretory function of TAMs is essential for combating infections, regulating immune responses, and promoting tissue repair. Extracellular vesicles (EVs) are nanovesicles that are secreted by cells. They play a crucial role in mediating intercellular information transfer between cells. EVs can be secreted by almost all types of cells, and they contain proteins, microRNAs, mRNAs, and even long non-coding RNAs (lncRNAs) that have been retained from the parental cell through the process of biogenesis. EVs can influence the function and behavior of target cells by delivering their contents, thus reflecting, to some extent, the characteristics of their parental cells. Here, we provide an overview of the role of M1 macrophages and their EVs in cancer therapy by exploring the impact of M1 macrophage-derived EVs (M1-EVs) on tumors by transferring small microRNAs. Additionally, we discuss the potential of M1-EVs as drug carriers and the possibility of reprogramming M2 macrophages into M1 macrophages for disease treatment. We propose that M1-EVs play a crucial role in cancer therapy by transferring microRNAs and loading them with drugs. Reprogramming M2 macrophages into M1 macrophages holds great promise in the treatment of cancers.

A Modified Tridecapeptide Probe for Imaging Cell Junction.

Cell junctions, which are typically associated with dynamic cytoskeletons, are essential for a wide range of cellular activities, including cell migration, cell communication, barrier function and signal transduction. Observing cell junctions in real-time can help us understand the mechanisms by which they regulate these cellular activities. This study examined the binding capacity of a modified tridecapeptide from Connexin 43 (Cx43) to the cell junction protein zonula occludens-1 (ZO-1). The goal was to create a fluorescent peptide that can label cell junctions. A cell-penetrating peptide was linked to the modified tridecapeptide. The heterotrimeric peptide molecule was then synthesized. The binding of the modified tridecapeptide was tested using pulldown and immunoprecipitation assays. The ability of the peptide to label cell junctions was assessed by adding it to fixed or live Caco-2 cells. The testing assays revealed that the Cx43-derived peptide can bind to ZO-1. Additionally, the peptide was able to label cell junctions of fixed cells, although no obvious cell junction labeling was observed clearly in live cells, probably due to the inadequate affinity. These findings suggest that labeling cell junctions using a peptide-based strategy is feasible. Further efforts to improve its affinity are warranted in the future.

Strategies for activity analysis of single nucleotide polymorphisms associated with human diseases.

Genome-wide association studies (GWAS) have identified a large number of single nucleotide polymorphism (SNP) sites associated with human diseases. In the annotation of human diseases, especially cancers, SNPs, as an important component of genetic factors, have gained increasing attention. Given that most of the SNPs are located in non-coding regions, the functional verification of these SNPs is a great challenge. The key to functional annotation for risk SNPs is to screen SNPs with regulatory activity from thousands of disease associated-SNPs. In this review, we systematically recapitulate the characteristics and functional roles of SNP sites, discuss three parallel reporter screening strategies in detail based on barcode tag classification, and recommend the common in silico strategies to help supplement the annotation of SNP sites with epigenetic activity analysis, prediction of target genes and trans-acting factors. We hope that this review will contribute to this exuberant research field by providing robust activity analysis strategies that can facilitate the translation of GWAS results into personalized diagnosis and prevention measures for human diseases.

Chromosomal instability (CIN) in HAP1 cell lines revealed by multiplex fluorescence in situ hybridisation (M-FISH).

BACKGROUND: HAP1, a near-haploid human leukemic cancer cell line is often used in combination with CRISPR-Cas9 gene editing technology for genetic screens. HAP1 carries the Philadelphia chromosome (Ph) and an additional ~ 30 Mb fragment of chromosome 15 inserted into chromosome 19. The potential use of an in vitro cell line as a model system in biomedical research studies depends on its ability to maintain genome stability. Being a cancer cell line with a near-haploid genome, HAP1 is prone to genetic instability, which is further compounded by its tendency to diploidise in culture spontaneously. Moreover, CRISPR-Cas9 gene editing coupled with prolonged in-vitro cell culturing has the potential to induce unintended 'off-target' cytogenetic mutations. To gain an insight into chromosomal instability (CIN) and karyotype heterogeneity, 19 HAP1 cell lines were cytogenetically characterised, 17 of which were near-haploids and two double-haploids, using multiplex fluorescence in situ hybridisation (M-FISH), at single cell resolution. We focused on novel numerical (N) and structural (S) CIN and discussed the potential causal factors for the observed instability. For each cell line we examined its ploidy, gene editing status and its length of in-vitro cell culturing. RESULTS: Sixteen of the 19 cell lines had been gene edited with passage numbers ranging from 10 to 35. Diploidisation in 17 near-haploid cell lines ranged from 4 to 35% and percentage of N- and S-CIN in [1n] and [2n] metaphases ranged from 7 to 50% with two cell lines showing no CIN. Percentage of cells with CIN in the two double-haploid cell lines were 96% and 100% respectively. The most common S-CIN observed was deletion followed by translocation of both types, non-reciprocal and Robertsonian. Interestingly, we observed a prevalence of S-CIN associated with chromosome 13 in both near-and double-haploid cell lines, with a high incidence of Robertsonian translocation involving chromosome 13. Furthermore, locus-specific BAC (bacterial artificial chromosome) FISH enabled us to show for the first time that the additional chromosome 15 fragment is inserted into the p-arm rather than the q-arm of chromosome 19 of the HAP1 genome. CONCLUSION: Our study revealed a high incidence of CIN leading to karyotype heterogeneity in majority of the HAP1 cell lines with the number of chromosomal aberrations varying between cell lines. A noteworthy observation was the high frequency of structural chromosomal aberrations associated with chromosome 13. We showed that CRISPR-Cas9 gene editing technology in combination with spontaneous diploidisation and prolonged in-vitro cell culturing is potentially instrumental in inducing further chromosomal rearrangements in the HAP1 cell lines with existing CIN. We highlight the importance of maintaining cell lines at low passage and the need for regular monitoring to prevent implications in downstream applications. Our study also established that the additional fragment of chromosome 15 in the HAP1 genome is inserted into chromosome 19p rather than 19q.

Werner Helicase Is a Synthetic-Lethal Vulnerability in Mismatch Repair-Deficient Colorectal Cancer Refractory to Targeted Therapies, Chemotherapy, and Immunotherapy.

Targeted therapies, chemotherapy, and immunotherapy are used to treat patients with mismatch repair-deficient (dMMR)/microsatellite instability-high (MSI-H) colorectal cancer. The clinical effectiveness of targeted therapy and chemotherapy is limited by resistance and drug toxicities, and about half of patients receiving immunotherapy have disease that is refractory to immune checkpoint inhibitors. Loss of Werner syndrome ATP-dependent helicase (WRN) is a synthetic lethality in dMMR/MSI-H cells. To inform the development of WRN as a therapeutic target, we performed WRN knockout or knockdown in 60 heterogeneous dMMR colorectal cancer preclinical models, demonstrating that WRN dependency is an almost universal feature and a robust marker for patient selection. Furthermore, models of resistance to clinically relevant targeted therapy, chemotherapy, and immunotherapy retain WRN dependency. These data show the potential of therapeutically targeting WRN in patients with dMMR/MSI-H colorectal cancer and support WRN as a therapeutic option for patients with dMMR/MSI-H cancers refractory to current treatment strategies. SIGNIFICANCE: We found that a large, diverse set of dMMR/MSI-H colorectal cancer preclinical models, including models of treatment-refractory disease, are WRN-dependent. Our results support WRN as a promising synthetic-lethal target in dMMR/MSI-H colorectal cancer tumors as a monotherapy or in combination with targeted agents, chemotherapy, or immunotherapy.This article is highlighted in the In This Issue feature, p. 1861.

Combinatorial CRISPR screen identifies fitness effects of gene paralogues.

Genetic redundancy has evolved as a way for human cells to survive the loss of genes that are single copy and essential in other organisms, but also allows tumours to survive despite having highly rearranged genomes. In this study we CRISPR screen 1191 gene pairs, including paralogues and known and predicted synthetic lethal interactions to identify 105 gene combinations whose co-disruption results in a loss of cellular fitness. 27 pairs influence fitness across multiple cell lines including the paralogues FAM50A/FAM50B, two genes of unknown function. Silencing of FAM50B occurs across a range of tumour types and in this context disruption of FAM50A reduces cellular fitness whilst promoting micronucleus formation and extensive perturbation of transcriptional programmes. Our studies reveal the fitness effects of FAM50A/FAM50B in cancer cells.

Correction to: PiggyBac mutagenesis and exome sequencing identify genetic driver landscapes and potential therapeutic targets of EGFR-mutant gliomas.

An amendment to this paper has been published and can be accessed via the original article.

PiggyBac mutagenesis and exome sequencing identify genetic driver landscapes and potential therapeutic targets of EGFR-mutant gliomas.

BACKGROUND: Glioma is the most common intrinsic brain tumor and also occurs in the spinal cord. Activating EGFR mutations are common in IDH1 wild-type gliomas. However, the cooperative partners of EGFR driving gliomagenesis remain poorly understood. RESULTS: We explore EGFR-mutant glioma evolution in conditional mutant mice by whole-exome sequencing, transposon mutagenesis forward genetic screening, and transcriptomics. We show mutant EGFR is sufficient to initiate gliomagenesis in vivo, both in the brain and spinal cord. We identify significantly recurrent somatic alterations in these gliomas including mutant EGFR amplifications and Sub1, Trp53, and Tead2 loss-of-function mutations. Comprehensive functional characterization of 96 gliomas by genome-wide piggyBac insertional mutagenesis in vivo identifies 281 known and novel EGFR-cooperating driver genes, including Cdkn2a, Nf1, Spred1, and Nav3. Transcriptomics confirms transposon-mediated effects on expression of these genes. We validate the clinical relevance of new putative tumor suppressors by showing these are frequently altered in patients' gliomas, with prognostic implications. We discover shared and distinct driver mutations in brain and spinal gliomas and confirm in vivo differential tumor suppressive effects of Pten between these tumors. Functional validation with CRISPR-Cas9-induced mutations in novel genes Tead2, Spred1, and Nav3 demonstrates heightened EGFRvIII-glioma cell proliferation. Chemogenomic analysis of mutated glioma genes reveals potential drug targets, with several investigational drugs showing efficacy in vitro. CONCLUSION: Our work elucidates functional driver landscapes of EGFR-mutant gliomas, uncovering potential therapeutic strategies, and provides new tools for functional interrogation of gliomagenesis.

Analysis pipelines for cancer genome sequencing in mice.

Mouse models of human cancer have transformed our ability to link genetics, molecular mechanisms and phenotypes. Both reverse and forward genetics in mice are currently gaining momentum through advances in next-generation sequencing (NGS). Methodologies to analyze sequencing data were, however, developed for humans and hence do not account for species-specific differences in genome structures and experimental setups. Here, we describe standardized computational pipelines specifically tailored to the analysis of mouse genomic data. We present novel tools and workflows for the detection of different alteration types, including single-nucleotide variants (SNVs), small insertions and deletions (indels), copy-number variations (CNVs), loss of heterozygosity (LOH) and complex rearrangements, such as in chromothripsis. Workflows have been extensively validated and cross-compared using multiple methodologies. We also give step-by-step guidance on the execution of individual analysis types, provide advice on data interpretation and make the complete code available online. The protocol takes 2-7 d, depending on the desired analyses.

Expanded potential stem cell media as a tool to study human developmental hematopoiesis in vitro.

Pluripotent stem cell (PSC) differentiation in vitro represents a powerful and tractable model to study mammalian development and an unlimited source of cells for regenerative medicine. Within hematology, in vitro PSC hematopoiesis affords novel insights into blood formation and represents an exciting potential approach to generate hematopoietic and immune cell types for transplantation and transfusion. Most studies to date have focused on in vitro hematopoiesis from mouse PSCs and human PSCs. However, differences in mouse and human PSC culture protocols have complicated the translation of discoveries between these systems. We recently developed a novel chemical media formulation, expanded potential stem cell medium (EPSCM), that maintains mouse PSCs in a unique cellular state and extraembryonic differentiation capacity. Herein, we describe how EPSCM can be directly used to stably maintain human PSCs. We further demonstrate that human PSCs maintained in EPSCM can spontaneously form embryoid bodies and undergo in vitro hematopoiesis using a simple differentiation protocol, similar to mouse PSC differentiation. EPSCM-maintained human PSCs generated at least two hematopoietic cell populations, which displayed distinct transcriptional profiles by RNA-sequencing (RNA-seq) analysis. EPSCM also supports gene targeting using homologous recombination, affording generation of an SPI1 (PU.1) reporter PSC line to study and track in vitro hematopoiesis. EPSCM therefore provides a useful tool not only to study pluripotency but also hematopoietic cell specification and developmental-lineage commitment.

Functional linkage of gene fusions to cancer cell fitness assessed by pharmacological and CRISPR-Cas9 screening.

Many gene fusions are reported in tumours and for most their role remains unknown. As fusions are used for diagnostic and prognostic purposes, and are targets for treatment, it is crucial to assess their function in cancer. To systematically investigate the role of fusions in tumour cell fitness, we utilized RNA-sequencing data from 1011 human cancer cell lines to functionally link 8354 fusion events with genomic data, sensitivity to >350 anti-cancer drugs and CRISPR-Cas9 loss-of-fitness effects. Established clinically-relevant fusions were identified. Overall, detection of functional fusions was rare, including those involving cancer driver genes, suggesting that many fusions are dispensable for tumour fitness. Therapeutically actionable fusions involving RAF1, BRD4 and ROS1 were verified in new histologies. In addition, recurrent YAP1-MAML2 fusions were identified as activators of Hippo-pathway signaling in multiple cancer types. Our approach discriminates functional fusions, identifying new drivers of carcinogenesis and fusions that could have clinical implications.

PiggyBac transposon tools for recessive screening identify B-cell lymphoma drivers in mice.

B-cell lymphoma (BCL) is the most common hematologic malignancy. While sequencing studies gave insights into BCL genetics, identification of non-mutated cancer genes remains challenging. Here, we describe PiggyBac transposon tools and mouse models for recessive screening and show their application to study clonal B-cell lymphomagenesis. In a genome-wide screen, we discover BCL genes related to diverse molecular processes, including signaling, transcriptional regulation, chromatin regulation, or RNA metabolism. Cross-species analyses show the efficiency of the screen to pinpoint human cancer drivers altered by non-genetic mechanisms, including clinically relevant genes dysregulated epigenetically, transcriptionally, or post-transcriptionally in human BCL. We also describe a CRISPR/Cas9-based in vivo platform for BCL functional genomics, and validate discovered genes, such as Rfx7, a transcription factor, and Phip, a chromatin regulator, which suppress lymphomagenesis in mice. Our study gives comprehensive insights into the molecular landscapes of BCL and underlines the power of genome-scale screening to inform biology.

Structural rearrangements generate cell-specific, gene-independent CRISPR-Cas9 loss of fitness effects.

BACKGROUND: CRISPR-Cas9 genome editing is widely used to study gene function, from basic biology to biomedical research. Structural rearrangements are a ubiquitous feature of cancer cells and their impact on the functional consequences of CRISPR-Cas9 gene-editing has not yet been assessed. RESULTS: Utilizing CRISPR-Cas9 knockout screens for 250 cancer cell lines, we demonstrate that targeting structurally rearranged regions, in particular tandem or interspersed amplifications, is highly detrimental to cellular fitness in a gene-independent manner. In contrast, amplifications caused by whole chromosomal duplication have little to no impact on fitness. This effect is cell line specific and dependent on the ploidy status. We devise a copy-number ratio metric that substantially improves the detection of gene-independent cell fitness effects in CRISPR-Cas9 screens. Furthermore, we develop a computational tool, called Crispy, to account for these effects on a single sample basis and provide corrected gene fitness effects. CONCLUSION: Our analysis demonstrates the importance of structural rearrangements in mediating the effect of CRISPR-Cas9-induced DNA damage, with implications for the use of CRISPR-Cas9 gene-editing in cancer cells.

ATM orchestrates the DNA-damage response to counter toxic non-homologous end-joining at broken replication forks.

Mutations in the ATM tumor suppressor gene confer hypersensitivity to DNA-damaging chemotherapeutic agents. To explore genetic resistance mechanisms, we performed genome-wide CRISPR-Cas9 screens in cells treated with the DNA topoisomerase I inhibitor topotecan. Thus, we here establish that inactivating terminal components of the non-homologous end-joining (NHEJ) machinery or of the BRCA1-A complex specifically confer topotecan resistance to ATM-deficient cells. We show that hypersensitivity of ATM-mutant cells to topotecan or the poly-(ADP-ribose) polymerase (PARP) inhibitor olaparib reflects delayed engagement of homologous recombination at DNA-replication-fork associated single-ended double-strand breaks (DSBs), allowing some to be subject to toxic NHEJ. Preventing DSB ligation by NHEJ, or enhancing homologous recombination by BRCA1-A complex disruption, suppresses this toxicity, highlighting a crucial role for ATM in preventing toxic LIG4-mediated chromosome fusions. Notably, suppressor mutations in ATM-mutant backgrounds are different to those in BRCA1-mutant scenarios, suggesting new opportunities for patient stratification and additional therapeutic vulnerabilities for clinical exploitation.

Shieldin complex promotes DNA end-joining and counters homologous recombination in BRCA1-null cells.

BRCA1 deficiencies cause breast, ovarian, prostate and other cancers, and render tumours hypersensitive to poly(ADP-ribose) polymerase (PARP) inhibitors. To understand the resistance mechanisms, we conducted whole-genome CRISPR-Cas9 synthetic-viability/resistance screens in BRCA1-deficient breast cancer cells treated with PARP inhibitors. We identified two previously uncharacterized proteins, C20orf196 and FAM35A, whose inactivation confers strong PARP-inhibitor resistance. Mechanistically, we show that C20orf196 and FAM35A form a complex, 'Shieldin' (SHLD1/2), with FAM35A interacting with single-stranded DNA through its C-terminal oligonucleotide/oligosaccharide-binding fold region. We establish that Shieldin acts as the downstream effector of 53BP1/RIF1/MAD2L2 to promote DNA double-strand break (DSB) end-joining by restricting DSB resection and to counteract homologous recombination by antagonizing BRCA2/RAD51 loading in BRCA1-deficient cells. Notably, Shieldin inactivation further sensitizes BRCA1-deficient cells to cisplatin, suggesting how defining the SHLD1/2 status of BRCA1-deficient tumours might aid patient stratification and yield new treatment opportunities. Highlighting this potential, we document reduced SHLD1/2 expression in human breast cancers displaying intrinsic or acquired PARP-inhibitor resistance.

The Origins and Vulnerabilities of Two Transmissible Cancers in Tasmanian Devils.

Transmissible cancers are clonal lineages that spread through populations via contagious cancer cells. Although rare in nature, two facial tumor clones affect Tasmanian devils. Here we perform comparative genetic and functional characterization of these lineages. The two cancers have similar patterns of mutation and show no evidence of exposure to exogenous mutagens or viruses. Genes encoding PDGF receptors have copy number gains and are present on extrachromosomal double minutes. Drug screening indicates causative roles for receptor tyrosine kinases and sensitivity to inhibitors of DNA repair. Y chromosome loss from a male clone infecting a female host suggests immunoediting. These results imply that Tasmanian devils may have inherent susceptibility to transmissible cancers and present a suite of therapeutic compounds for use in conservation.

Evolutionary routes and KRAS dosage define pancreatic cancer phenotypes.

The poor correlation of mutational landscapes with phenotypes limits our understanding of the pathogenesis and metastasis of pancreatic ductal adenocarcinoma (PDAC). Here we show that oncogenic dosage-variation has a critical role in PDAC biology and phenotypic diversification. We find an increase in gene dosage of mutant KRAS in human PDAC precursors, which drives both early tumorigenesis and metastasis and thus rationalizes early PDAC dissemination. To overcome the limitations posed to gene dosage studies by the stromal richness of PDAC, we have developed large cell culture resources of metastatic mouse PDAC. Integration of cell culture genomes, transcriptomes and tumour phenotypes with functional studies and human data reveals additional widespread effects of oncogenic dosage variation on cell morphology and plasticity, histopathology and clinical outcome, with the highest KrasMUT levels underlying aggressive undifferentiated phenotypes. We also identify alternative oncogenic gains (Myc, Yap1 or Nfkb2), which collaborate with heterozygous KrasMUT in driving tumorigenesis, but have lower metastatic potential. Mechanistically, different oncogenic gains and dosages evolve along distinct evolutionary routes, licensed by defined allelic states and/or combinations of hallmark tumour suppressor alterations (Cdkn2a, Trp53, Tgfß-pathway). Thus, evolutionary constraints and contingencies direct oncogenic dosage gain and variation along defined routes to drive the early progression of PDAC and shape its downstream biology. Our study uncovers universal principles of Ras-driven oncogenesis that have potential relevance beyond pancreatic cancer.

Genomic evidence reveals a radiation of placental mammals uninterrupted by the KPg boundary.

The timing of the diversification of placental mammals relative to the Cretaceous-Paleogene (KPg) boundary mass extinction remains highly controversial. In particular, there have been seemingly irreconcilable differences in the dating of the early placental radiation not only between fossil-based and molecular datasets but also among molecular datasets. To help resolve this discrepancy, we performed genome-scale analyses using 4,388 loci from 90 taxa, including representatives of all extant placental orders and transcriptome data from flying lemurs (Dermoptera) and pangolins (Pholidota). Depending on the gene partitioning scheme, molecular clock model, and genic deviation from molecular clock assumptions, extensive sensitivity analyses recovered widely varying diversification scenarios for placental mammals from a given gene set, ranging from a deep Cretaceous origin and diversification to a scenario spanning the KPg boundary, suggesting that the use of suboptimal molecular clock markers and methodologies is a major cause of controversies regarding placental diversification timing. We demonstrate that reconciliation between molecular and paleontological estimates of placental divergence times can be achieved using the appropriate clock model and gene partitioning scheme while accounting for the degree to which individual genes violate molecular clock assumptions. A birth-death-shift analysis suggests that placental mammals underwent a continuous radiation across the KPg boundary without apparent interruption by the mass extinction, paralleling a genus-level radiation of multituberculates and ecomorphological diversification of both multituberculates and therians. These findings suggest that the KPg catastrophe evidently played a limited role in placental diversification, which, instead, was likely a delayed response to the slightly earlier radiation of angiosperms.

Recurrent Rearrangements of Human Amylase Genes Create Multiple Independent CNV Series.

The human amylase gene cluster includes the human salivary (AMY1) and pancreatic amylase genes (AMY2A and AMY2B), and is a highly variable and dynamic region of the genome. Copy number variation (CNV) of AMY1 has been implicated in human dietary adaptation, and in population association with obesity, but neither of these findings has been independently replicated. Despite these functional implications, the structural genomic basis of CNV has only been defined in detail very recently. In this work, we use high-resolution analysis of copy number, and analysis of segregation in trios, to define new, independent allelic series of amylase CNVs in sub-Saharan Africans, including a series of higher-order expansions of a unit consisting of one copy each of AMY1, AMY2A, and AMY2B. We use fiber-FISH (fluorescence in situ hybridization) to define unexpected complexity in the accompanying rearrangements. These findings demonstrate recurrent involvement of the amylase gene region in genomic instability, involving at least five independent rearrangements of the pancreatic amylase genes (AMY2A and AMY2B). Structural features shared by fundamentally distinct lineages strongly suggest that the common ancestral state for the human amylase cluster contained more than one, and probably three, copies of AMY1.

Contrasting origin of B chromosomes in two cervids (Siberian roe deer and grey brocket deer) unravelled by chromosome-specific DNA sequencing.

BACKGROUND: B chromosomes are dispensable and variable karyotypic elements found in some species of animals, plants and fungi. They often originate from duplications and translocations of host genomic regions or result from hybridization. In most species, little is known about their DNA content. Here we perform high-throughput sequencing and analysis of B chromosomes of roe deer and brocket deer, the only representatives of Cetartiodactyla known to have B chromosomes. RESULTS: In this study we developed an approach to identify genomic regions present on chromosomes by high-throughput sequencing of DNA generated from flow-sorted chromosomes using degenerate-oligonucleotide-primed PCR. Application of this method on small cattle autosomes revealed a previously described KIT gene region translocation associated with colour sidedness. Implementing this approach to B chromosomes from two cervid species, Siberian roe deer (Capreolus pygargus) and grey brocket deer (Mazama gouazoubira), revealed dramatically different genetic content: roe deer B chromosomes consisted of two duplicated genomic regions (a total of 1.42-1.98 Mbp) involving three genes, while grey brocket deer B chromosomes contained 26 duplicated regions (a total of 8.28-9.31 Mbp) with 34 complete and 21 partial genes, including KIT and RET protooncogenes, previously found on supernumerary chromosomes in canids. Sequence variation analysis of roe deer B chromosomes revealed a high frequency of mutations and increased heterozygosity due to either amplification within B chromosomes or divergence between different Bs. In contrast, grey brocket deer B chromosomes were found to be more homogeneous and resembled autosomes in patterns of sequence variation. Similar tendencies were observed in repetitive DNA composition. CONCLUSIONS: Our data demonstrate independent origins of B chromosomes in the grey brocket and roe deer. We hypothesize that the B chromosomes of these two cervid species represent different stages of B chromosome sequences evolution: probably nascent and similar to autosomal copies in brocket deer, highly derived in roe deer. Based on the presence of the same orthologous protooncogenes in canids and brocket deer Bs we argue that genomic regions involved in B chromosome formation are not random. In addition, our approach is also applicable to the characterization of other evolutionary and clinical rearrangements.