Comprehensive genomic characterization of five canine lymphoid tumor cell lines
BMC Veterinary Research volume 12, Article number: 207 (2016)
Leukemia/lymphoma cell lines have been critical in the investigation of the pathogenesis and therapy of hematological malignancies. While human LL cell lines have generally been found to recapitulate the primary tumors from which they were derived, appropriate characterization including cytogenetic and transcriptional assessment is crucial for assessing their clinical predictive value.
In the following study, five canine LL cell lines, CLBL-1, Ema, TL-1 (Nody-1), UL-1, and 3132, were characterized using extensive immunophenotyping, karyotypic analysis, oligonucleotide array comparative genomic hybridization (oaCGH), and gene expression profiling. Genome-wide DNA copy number data from the cell lines were also directly compared with 299 primary canine round cell tumors to determine whether the cell lines represent primary tumors, and, if so, what subtype each most closely resembled.
Based on integrated analyses, CLBL-1 was classified as B-cell lymphoma, Ema and TL-1 as T-cell lymphoma, and UL-1 as T-cell acute lymphoblastic leukemia. 3132, originally classified as a B-cell lymphoma, was reclassified as a histiocytic sarcoma based on characteristic cytogenomic properties. In combination, these data begin to elucidate the clinical predictive value of these cell lines which will enhance the appropriate selection of in vitro models for future studies of canine hematological malignancies.
Hematological diseases in humans are widely heterogeneous including numerous molecular subtypes with wide ranging prognoses and therapeutic responses. Oncogenesis is well understood for some subtypes, while the molecular changes leading to other subtypes remain unknown . A similar level of molecular heterogeneity likely exists in the >1000 leukemia/lymphoma (LL) cell lines that have been described to date .
Although the use of cell lines provides numerous advantages, including ease of handling and manipulation, high homogeneity, provision of a continuous source of sample material, and accessibility to the scientific community , detailed characterization is crucial before they are used as an in vitro preclinical cancer models. Nearly 100 % of examined human LL cell lines carry stable genetic alterations and karyotypic changes that maintain the major features of the original cells  while a number of other human LL cell line studies have provided evidence that these cell lines have a high clinical predictive value that may translate into a favorable response rates in Phase II clinical trials [4, 5].
Cell line characterization minimally requires that a comprehensive set of immunophenotyping and cytogenetic data have been published . However, with the recent explosion of advances in the genomics field, it is now possible to provide a more in-depth assessment of cytogenetic and transcriptional characteristics of cell lines that can provide further insight into biological processes including chromosomal translocations, signaling pathways, mutational analysis, gene dysregulation, and RNAi gene silencing . Of the >1000 human LL cell lines described, ~40 % have been characterized in sufficient detail for accurate classification as discrete LL subtypes .
Spontaneously occurring lymphoid malignancies in dogs share the same histopathological and clinical features of their human couterparts, in addition to evolutionarily conserved chromosome aberrations and mutations, indicating shared pathogenesis across species [7–9]. There are only a small number of established canine LL cell lines, six of which have already been characterized at the genomic level [10, 11]. Similar to human LL cell lines, the importance of detailed genomic and phenotypic characterization in an effort to define a framework to assess their clinical predictive value was recently emphasized . Five additional canine LL cell lines, CLBL-1, Ema, UL-1, TL-1 (Nody-1), and 3132, have been used previously in in vitro studies [13–21]. These cell lines have varying levels of characterization, none of which includes an in-depth genomic and transcriptomic approach [22–24].
We present a comprehensive characterization of five canine LL cell lines starting with an extended panel of immunophenotyping. High resolution oligonucleotide array comparative genomic hybridization (oaCGH) was performed to assess genome-wide copy number status, and multicolor fluorescence in situ hybridization (FISH) analysis was used to further identify copy number imbalances and structural changes in karyotype architecture. Transcription status of each cell line was investigated using high-density array based gene expression profiling (GEP) and quantitative reverse transcriptase polymerase chain reaction (qRT-PCR). Additionally, genome wide copy number data of each cell line were compared with data from primary canine round cell tumors to further confirm their classification and relevance as in vitro preclinical models of lymphoid neoplasia for canine and comparative medicine.
Canine LL cell lines
Five previously established canine LL cell lines with varying levels of initial characterization were included in this study: CLBL-1  (kind gift from Dr. Barbara Rutgen, University of Veterinary Medicine Vienna, Austria, Ema  (kind gift from Dr. Takuya Mizuno, Yamaguchi University, Japan), TL-1  (Nody-1, kind gift from Dr. Yasuhiko Okamura, Iwate University, Japan), UL-1  (kind gift from Dr. Hajime Tsujimoto, University of Tokyo, Japan), and 3132 [24, 25] (kind gift from Dr. Mark Holmes, University of Cambridge, UK). All cell lines were maintained at 37 °C and 5 % CO2 in RPMI-1640 culture medium (Mediatech, Hendon, VA) supplemented with 10 % fetal bovine serum (FBS, Mediatech), 2 mM Glutamax (Life Technologies, Grand Island, NY), and 100 μg/ml Primocin (Invivogen, San Diego, CA) and tested negative using a PCR Mycoplasma test kit (Applichem, Cheshire, CT).
Immunophenotyping of each line was completed using flow cytometry as previously described [11, 26] at the UC Davis Leukocyte Antigen Biology Laboratory using a panel of monoclonal antibodies reactive with canine leukocyte antigens, including CD1a, CD1c, CD3, CD4, CD5, CD8α, CD8β, CD11a, CD11b, CD11c, CD11d, CD14, CD18, CD21, CD34, CD45, CD45RA, CD49a, CD54, CD79α, CD80, CD86, DM5, MHC-II, 5G2, AG5, TCRαβ, TCRγδ, TCRαβ, TCRCCγδ, and Thy-1 (CD90). Analysis at the NCSU Clinical Immunology Laboratory was also performed using a smaller number of antibodies, including CD3, CD4, CD5, CD8α, CD21, CD34, and CD79α/β.
Polymerase chain reaction for antigen receptor rearrangement (PARR) was completed as previously described [27, 28] to assess clonality and, possibly lineage. PCR products were separated using capillary gel electrophoresis (QIAxcel Electrophoresis System, Qiagen, Valencia, CA). A clonal sample was determined if one or more discrete bands were seen on the gel, and a polyclonal sample was determined if multiple bands or a smear of amplicons were seen. A negative sample was determined if no bands were seen.
Isolation of cell line DNA and RNA and generation of metaphase preparations
Aliquots of 1 × 107cells were removed from the same culture flask at the same time point for isolation of DNA and RNA and preparation of metaphase chromosome preparations to ensure consistency in downstream analyses. DNA was isolated using the DNeasy Blood and Tissue kit (Qiagen, Valencia, CA) and manufacturer’s protocol, and quantity and quality of DNA were evaluated using spectrophotometry (260/280 > 1.8) and gel electrophoresis. RNA was isolated using the RNeasy Plus Mini Kit (Qiagen) and manufacturer’s protocol, and assessed using the 2100 Bioanalyzer RNA 6000 Nano Kit (Agilent Technologies, Santa Clara, CA) to confirm an RNA integrity number (RIN) >9.0. Metaphase chromosome preparations were prepared from each cell line as previously described  using conventional techniques of colcimid arrest (final concentration of 50 ng/ml for 1 h), hypotonic treatment, and methanol-glacial acetic acid fixation prior to being dropped onto glass slides.
Identification of copy number aberrations (CNAs) using oaCGH
oaCGH was completed as previously described using a 180,000 feature canine oligonucleotide array (Agilent Technologies) with repeat-masked 60mer oligonucleotides spaced ~13 kb across the genome . An equimolar pool of DNA from 25 clinically healthy female dogs was used as a common reference for all cell lines. Cell line and reference DNA was labeled with Cyanine 3-dUTP and Cyanine 5-dUTP, respectively, using the Agilent Enzymatic Labeling Kit, and probe hybridization, array washing, and scanning was performed as described elsewhere .
Scan data were processed using Feature Extraction v.10.10 software (Agilent Technologies) and imported into Nexus Copy Number v7.5 (Biodiscovery, Hawthorne, CA). Raw data were evaluated to identify and exclude probes displaying non-uniform hybridization or signal saturation, and copy number calls were made using Biodiscovery’s FASST2 segmentation algorithm. Copy number calls were based on a minimum of three consecutive probes per segment, and mean log2 cell line:reference thresholds of +/− 0.2 were used to define gain and loss, respectively. Genes within the defined intervals were identified using the UCSC canine genome browser (CanFam2 assembly; http://genome.ucsc.edu/) and the NCBI gene database (http://www.ncbi.nlm.nih.gov/gene). Genes previously associated with cancer were based on those reported in the Cancer Gene Census (http://cancer.sanger.ac.uk/cosmic/census) .
Further statistical analyses using the Feature Extraction data were performed using R . The signals (rProcessedSignal and gProcessedSignal) were normalized using the following equations:
Where, Processed ratio is the centered and normalized ratio of the Agilent processed fluorescent signals. Segmentation was performed across all chromosomes using circular binary segmentation . Data were further dichotomized as gain (1), no change (0), or loss (−1), based on segments that were +/− 3 MAD (mean absolute deviation) from the median of each sample’s response across all chromosomes. Hierarchical clustering of the five cell lines was performed using dichotomous data using Euclidean distance and Ward’s method. Additionally, hierarchical clustering of the cell lines with 299 canine primary round cell tumors including 123 leukemias , 106 lymphomas (Thomas et al., in preparation), and 70 histiocytic malignancies (Kennedy et al., in preparation) was performed using segmented data using Euclidean distance and Ward’s method. Tumor type and subtype were both annotated on the heatmap.
FISH of all cells lines was performed as previously described  to evaluate structural changes and verify oaCGH copy number data using panels of clones from the CHORI-82 dog bacterial artificial chromosome (BAC) library (www.chori.org). Initially, two clones (326 K03 and 330E21), previously determined to hybridize to the centromeric regions of canine autosomes , were fluorescently labeled and hybridized to metaphase preparations of each cell line. Centromeric signals were used to properly orient the chromosomes, aid in confirming modal chromosome count, and identify bi-armed chromosomes.
Twenty additional BAC clones were selected to contain known oncogenes and tumor suppressor genes that met two of the three following criteria: (1) located in a region of CNA in at least one of the five cells lines, (2) displayed differential expression between cell lines based on the GEP, and (3) have been associated with human and/or canine lymphoid malignancies in prior studies [8, 9, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46]. The BAC clones selected represented the following genes: BCL11B, IGH, VEGFA, CCNC, FOXO3A, CDKN2A, MYC, KIT, CDK6, EZH2, MYCBP2, FLT3, PTEN, HEY1, E2F5, NFKB2, ERG, MLLT2, CD83, and DEK (Table 1). The 20 clones were divided into four panels of five for multicolor FISH, and each panel was hybridized to metaphase spreads of each cell line and healthy dog controls. Copy number status of each probe was scored in at least 50 cells of each cell line and normal controls.
Total RNA from each cell line was used to perform gene expression profiling (GEP) as described elsewhere [11, 47] using the GeneChip Canine Genome 2.0 array (Affymetrix, Santa Clara, CA) which is comprised of 18,000 Canis familiaris mRNA transcripts and over 20,000 non-redundant predicted genes. Additionally, total RNA was isolated as described from lymph nodes harvested from six healthy mixed breed dogs that showed no evidence of lymphoid neoplasia at necropsy . Microarrays were processed by the Lineberger Functional Genomic Core Facility at the University of North Carolina Chapel Hill. Total RNA (1 μg) was processed for microarray hybridization using the MessageAmp II-Biotin Enhanced Kit (Ambion, Life Technologies, Grand Island, NY) and hybridization was performed according to Affymetrix technical protocols. GEP analysis was performed using GeneSpring GX v12 (Agilent Technologies). Expression array data were normalized using the GC-RMA procedure  and signals were median-centered across all arrays. Data were filtered to remove probe sets with limited variation (standard deviation <2) across all arrays, and fold change analysis was performed for each cell line using the averaged expression data from the normal lymph node controls as baseline. Additionally, unsupervised hierarchical clustering analysis was performed across all cell lines and controls using the filtered probe set. Functional analysis was performed by evaluating for enrichment in genes that were up- or down-regulated in each cell line by >5-fold compared to normal lymph node controls. Enrichment analysis in Gene Ontology (GO) biological processes and Kyoto Encyclopedia of Genes and Genome (KEGG) Pathways was completed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) v 6.7 [49, 50].
One ug of total RNA from each cell line and two normal lymph node controls was used to perform qRT-PCR to validate GEP as previously described . qRT-PCR was performed using Applied Biosystems OneStepPlus Real-Time PCR system (Life Technologies) and a cycling protocol with an initial denaturation at 95 °C for 3 min; followed by 40 cycles of 95 °C for 3 s, 62 °C for 20 s, and 72 °C for 15 s; with a final extension at 72 °C for 5 min, followed by a melt curve analysis. All assays were performed in triplicate.
Primers were designed using NCBI Primer-BLAST as previously described  for MYC, KIT, FLT3, PTEN, and RPL32 (Table 2). Relative quantification using the comparative Ct method (ΔΔCT) was performed as described previously using normal lymph node as the baseline for comparisons. RPL32 was used as the reference gene based on its stable expression across all samples in GEP analysis and previous identification as a suitable qRT-PCR reference gene for canine lymphoid neoplasia .
Flow cytometry data are presented in Table 3. All cell lines expressed CD45 and CD45RA, and 4/5 expressed CD18, which verifies a leukocytic origin. CLBL-1 displayed positive staining for CD1a, CD1c, CD11a, CD11b, CD54, CD79α, CD80, CD86, and MHC-II and no staining for CD3, CD4, CD5, CD8, and Thy-1 indicating a mature B-cell origin, which is further supported by a monoclonal product resulting from PCR analysis of the IGH gene. Ema was positive for Thy-1, a T-cell marker. TL-1 cells were negative for most antigens, however weak expression of MHCII was observed, which may indicate a T-cell phenotype as MHCII is expected to have high intensity on B cells and antigen presenting cells, with much weaker expression observed on T cells . TL-1 and Ema both were found to have rearrangement of TCRG gene (oligoclonal for TL-1), further supporting a mature T-cell phenotype for both cell lines. UL-1 expressed most antigens including CD34, which indicates an immature precursor phenotype. Positive expression for CD3, CD4, CD5, CD8α, CD8β, CD11d, and Thy-1 all support a T-cell phenotype. UL-1 is also positive for CD14, which is normally a marker for myeloid cells; however, expression is also demonstrated in B and T cell precursors but not mature cells , further supporting an immature phenotype. Additionally, PCR of TCRG indicated a monoclonal product further supporting T-cell phenotype. 3132 is most likely of dendritc/histiocytic origin based on the positive expression of CD1a and CD11c, which is characteristic of histiocytic diseases  and a combination not found in other cell lines of lymphoid origin in this or previous studies [10, 11]. Additionally, strong intensity for MHC-II further supports dendritic cell origin. However, a T-cell origin cannot be completely ruled out based on expression of CD4, CD8α, and a monoclonal product resulting from PCR of TCRG.
Chromosome enumeration and centromere localization via FISH of metaphase chromosomes was used to assess the gross karyotypic architecture of each cell line. The normal dog karyotype includes 38 pairs of acrocentric autosomes, a large sub-metacentric X chromosome and a small metacentric Y chromosome . Enumeration of chromosomes from 30 metaphase spreads of each cell line indicated varying levels of aneuploidy (Fig. 1). CLBL-1, Ema, and 3132 were all hypodiploid, while TL-1 and UL-1 contained normal modal chromosome counts (Table 4).
Genome wide copy number imbalance
oaCGH analysis revealed varying levels of genomic imbalance in each of the cell lines with detected CNAs ranging in size from 19.5 Mb to entire chromosomes (Fig. 2). Large differences in the number of aberrations (Fig. 3a) and the percent of the genome impacted by regions of imbalance (Fig. 3b) were also noted. CLBL-1, TL-1, and UL-1 have a similar and even distribution of the number of aberrant gains and losses; whereas Ema was skewed towards an increase in the number of losses, and 3132 was skewed towards regions of copy number gain. When the cumulative size of CNAs and the percent of the genome impacted by imbalance were assessed, Ema and TL-1 were minimally changed which is reiterated from the oaCGH profiles in which both cell lines had less visually obvious CNAs. CLBL-1 and UL-1 had similar total percent of the genome changed, but CLBL-1 had a greater percent of the genome with imbalances of loss and UL-1 had a greater percent of the genome gained. 3132 had the greatest percent of the genome changed which was expected given the number of aberrations detected and the extensive genomic complexity visualized in the oaCGH profile. Cancer-associated genes in regions of imbalance were identified and are detailed in Additional file 1: Table S1. Cell lines were clustered to determine their relatedness to each other based on genome-wide copy number assessment (Fig. 4). 3132 branched away from the other four, and Ema and TL-1, the two mature T-cell neoplasms based on immunophenotyping, were the most closely related.
FISH analysis was completed for 20 genes relevant to lymphoid neoplasia to verify oaCGH data, enumerate the level of imbalance at each locus, and investigate any structural aberrations associated with the selected genes of interest. The 20 genes were divided into four panels of five for multicolor FISH analysis. Each of the four FISH panels is identified in the insets of Fig. 5 with control chromosomes from clinically healthy dogs indicating the appropriate localization of each BAC clone. A representative metaphase spread with CNAs or structural changes in one of the cell lines is shown in each panel.
Structural changes were identified in each cell line. CLBL-1 contains a derivative 13 chromosome that appears to include two copies of Canis familiaris 13 (CFA 13) joined by the centromere resulting in the formation of a metacentric chromosome (Fig. 5b), in addition to a grossly normal copy of CFA 13 based on the location of the 2 probes on CFA 13. A similar abnormality occurred in UL-1 with chromosome 25, as is evidenced by 2 copies of FLT3 forming a metacentric chromosome in addition to an acrocentric chromosome with FLT3 placement in the expected location on CFA 25 (Fig. 5c). Both Ema and CLBL-1 contained a derivative metacentric chromosome resulting from the apparent fusion of CFA 28 and CFA 35 (Fig. 5d). Ema contained additional copies of CFA 28 and CFA 35, while CLBL-1 contained another CFA 28 and a heterozygous loss of CFA 35. Numerous structural aberrations occurred in 3132, and appeared to be random from cell to cell. Further detail of these structural changes was difficult to ascertain given the heterogeneity observed in the metaphase spreads.
A summary of the frequency of CNA at each of the 20 investigated loci is presented in Fig. 6. Ema and TL-1 showed grossly normal copy number changes, as was seen in the oaCGH; however, both were found to have a homozygous deletion of CDKN2A. Ema also had a heterozygous loss of MLLT2 in 100 % of cells and a heterozygous loss of FLT3 in 15 % of cells which was detected via oaCGH. Additionally, TL-1 had three copies of E2F5 and HEY1, both located on CFA 30, in 100 % of cells. UL-1 and CLBL-1 showed copy number changes in nine and 10 of the loci evaluated, respectively. CLBL-1 exhibited heterozygous losses of CCNC, FOXO3A, CDK6, MYCBP2, PTEN, CD83, and DEK, and a copy number of three for MYC and KIT (both located on CFA 13) in the majority of cells (>87 %). UL-1 had more genes gained, with three copies of BCL11B, IGH, VEGFA, CCNC, FOXO3A, FLT3, and NFKB2 in the vast majority (>93 %) of cells. Additionally, 50 % of cells gained an extra copy of CDK6. FISH analysis of 3132 showed copy number changes at 19/20 loci, many of which were strikingly heterogeneous. MYC, KIT, and CD83 each have proportions of cells with normal copy number and gains of 3, 4, and >4 copies. Losses of ERG and EZH2 were identified in all or nearly all 3132 cells, with proportions (10–23 %) of cells displaying homozygous deletion.
Gene expression analysis
GEP data were filtered to remove probe sets with limited variation (standard deviation <2) across the five cell lines and six non-neoplastic lymph nodes resulting in 1153 probe sets used for subsequent analyses. Unsupervised clustering resulted in immediate branching of 3132 from the other 10 samples. The remaining samples were then separated into two discrete groups including the four cell lines, and the six control lymph nodes (Fig. 7). Normal lymph nodes unsurprisingly display greater transcription conservation across biological replicates compared with the cell lines as indicated by their shorter connecting branches. The cell lines displaying the greatest transcriptional similarity were Ema and TL-1, the two mature T-cell lines, as reiterated from the oaCGH clustering analysis. UL-1, most likely an immature T-cell phenotype, is also more similar to Ema and TL-1 than CLBL-1, a B-cell line.
Fold change analysis was also completed by comparing each cell line with the mean expression of the non-neoplastic lymph nodes. Fold changes of known cancer related genes and the transcripts with the 50 largest absolute fold changes (up or down) in each cell line are listed in Additional file 1: Tables S2 and S3, respectively, as genes that may have functional relevance. Additionally, genes with a fold change >5 were further assessed in functional gene-annotation enrichment analysis. Due to limited functional annotation of canine genes, official gene symbols were used for enrichment analysis based on human gene annotations. Up to the top ten enriched GO biological processes and KEGG pathways are listed in Additional file 1: Tables S4 and S5, respectively. All five cell lines exhibited upregulation of genes associated with the GO biological processes of cell proliferation and division, and three of five exhibited upregulation of genes involved in metabolic processes (e.g. sterol biosynthetic process, glucose metabolic process, cholesterol metabolic process). All five cell lines also exhibited downregulation of genes involved in immune response and leukocyte activation. When the KEGG pathways were assessed, four of the cell lines exhibited upregulation of metabolic and biosynthesis pathways (e.g. steroid biosynthesis, pyruvate metabolism) as similarly identified in the GO term analysis. TL-1 was also found to exhibit upregulation of genes in both the ERBB and Jak-STAT signaling pathways; while genes involved in the MAPK signaling pathways and the NOD-like and Toll-like receptor signaling pathways were upregulated in 3132. All cell lines exhibited downregulation of cell adhesion molecules.
qRT-PCR was performed to verify GEP findings and analyze the relationship between copy number change and expression change in four well known cancer associated genes including MYC, KIT, FLT3, and PTEN (Fig. 8). Homozygous loss of PTEN was observed in CLBL-1 and 3132, which both show the greatest decreases in PTEN expression. Three copies of FLT3 were present in UL-1, and this gene was upregulated 17 fold (and was the only cell line to show an upregulation). Fold changes in MYC and KIT are more variable across the cell lines and have minimal correlation with copy number changes observed. CLBL-1 and Ema had slight downregulation of MYC, while the other three showed modest upregulation. KIT was largely downregulated in CLBL-1, Ema, and 3132 (two of which had copy number gains), and upregulated 8-fold in TL-1 in which the loci was copy number neutral.
Genome-wide CNA comparison with primary tumors
Finally, we compared the five cell lines with genome wide copy number data generated in our laboratory from 299 primary round cell tumors via clustering analysis to determine whether the cell lines represented primary tumors, and, if so, what subtype each most closely resembled (Fig. 9). All cell lines were found to cluster within the primary tumors, indicating they are more closely related to primary tumors than other cell lines. 3132 segregated within a large cluster composed of primarily histiocytic malignancies. CLBL-1 branches from a cluster of mature B-cell neoplasms consisting primarily of B-cell lymphomas with a smaller subset of B-cell chronic lymphocytic leukemias. TL-1 and Ema both clustered in a somewhat heterogeneous grouping comprised chiefly of T-cell lymphomas, and both these cell lines were most closely related to a T-cell lymphoma. UL-1 grouped in a cluster composed of acute leukemias, most of which were lymphoid in origin.
We present a detailed characterization of five canine LL cell lines using a genome-wide molecular approach including oaCGH and GEP analysis. The data presented here, combined with previously published canine LL cell line characterizations , provides the opportunity to more appropriately select canine LL cell lines for preclinical in vitro studies. Overall, our data suggests that CLBL-1 represents a mature B-cell lymphoma as previously reported , and Ema and TL-1 both represent T-cell lymphomas as previously reported . In contrast, UL-1 should be reclassified as T-cell acute lymphoblastic leukemia (instead of a T-cell lymphoma) and 3132 should be reclassified as a histiocytic sarcoma (instead of a B-cell lymphoma).
CLBL-1 immunophenotyping in this study agreed with the immunophenotyping published when the cell line was initially established. Positive expression of CD11a, CD79α, CD45, CD45RA, and MHCII and no expression of CD3, CD4, CD5, CD8, CD11d, CD14, CD21, CD34, and CD56 classifies it as a B-cell lymphoma. Additionally, B-cell PARR revealed a monoclonal rearrangement, which we also found . CLBL-1 was the only cell line in this dataset that had been previously karyotyped . It was found to be hypodiploid with a modal chromosome number of 70–71, which matches our findings. There was only one biarmed chromosome mentioned, derivative chromosome 13, which we also identified via FISH analysis. Since centromere visualization was not previously completed, and a non-standard canine chromosome nomenclature was used, correlation of further previously identified chromosome aberrations with our findings via DAPI banding, FISH, and oaCGH was not possible. Numerous CNAs identified in CLBL-1 via oaCGH are shared with primary canine B-cell lymphomas (B-LSA). Gain of CFA 13 has been identified in 25 % of B-LSA. Loss of CFA14 and loss of the proximal region of CFA 3 have also been found to occur in ~10 % of B-LSA . Additionally, CLBL-1 clustered with primary mature B-cell lymphoid neoplasia when clustered with 299 primary canine round cell tumors, further supporting the classification of CLBL-1 as a B-cell lymphoma.
Ema immunophenotyping, at the time of establishment, was positive for CD3, CD45, CD45RA, and Thy-1. In our assessment, Ema was also positive (intermediate signal) for CD11a and CD18 and negative for CD3. Similarly, TL-1 was originally positive for CD3, CD18, CD45, CD45RA, and MHC11 while our immunophenotyping indicated no expression of CD3 and CD18. The slight differences in the immunophenotyping data may be due to factors that can influence flow cytometry data including, the amount of antibody used, the call threshold of positive versus negative, instrumentation, and other reagents. Alternatively, it’s possible that the loss or gain of antigens on the cell surface occurred while in culture. Both Ema and TL-1 were previously T-cell PARR positive, which is consistent with our PARR data. Overall, the changes observed in immunophenotype did not impact the overall interpretation of the classification of both these cell lines as mature T-cell lymphoid neoplasia.
Both Ema and TL-1 were found to have a copy number loss in the region containing CDKN2A in oaCGH analysis, which was further confirmed to be a homozygous loss of the gene loci in FISH analysis. CDKN2A loss occurs in >55 % of T-cell lymphomas (T-LSA) and <2 % of other subtypes of canine leukemia and lymphoma [8, 34]. Ema was also found to have a loss of the proximal half of CFA 22, which occurs in 20 % of T-LSA, and TL-1 had a gain of CFA 29 previously identified in 40 % of T-LSA . Ema and TL-1 were also the cell lines most closely related in clustering analyses using both the oaCGH and GEP data, and both segregated with primary T-LSA when clustered with primary canine round cell tumors.
UL-1 immunophenotyping at the time of establishment was positive for CD8α, CD18, CD45, CD45RA. We found additional positive expression of CD3, CD4, CD5, CD14, CD21, CD34, Thy-1, and MHC-II. All were detected at an intermediate level therefore, as discussed earlier, differences in methods and data analysis could have resulted in previous negative findings. UL-1 was also T-cell PARR positive which agrees with our PARR findings. UL-1 was previously classified as a T cell lymphoma; however, based on expression of CD34, a surface glycoprotein expressed on hematopoietic stem cells, we suggest it is more representative of a T-cell acute leukemia (T-cell ALL). In veterinary medicine, CD34 expression can be used to suggest an acute leukemia since it is less commonly expressed in malignancies of more mature cells, such as lymphomas and chronic lymphocytic leukemias [26, 56]. Several regions of CNA are also shared with primary canine ALLs including the loss of the distal end of CFA 1 which occurs in 20 % of ALLs, gain of CFA 12 and 25 which occurs in 15 and 10 % of ALLs, respectively, and loss of CFA 35 which occurs in 15 % of ALLs. Of the CNAs identified in UL-1, more were in common with ALLs than other subtypes. UL-1 also segregated with a group of ALLs when clustered with primary round cell tumors. Together, these findings support the classification of UL-1 as a T-cell ALL.
3132 immunophenotyping in our laboratory, in conjunction with the original description of the cells (cellular pleomorphism, multinucleated giant cells, numerous mitotic figures, extreme variations in nuclear:cytoplasm ratio) suggests this cell line is representative of a disseminated histiocytic sarcoma (HS). No immunophentypic or genomic analysis of this cell line has been previously reported, although it was later described as a B-cell lymphoma based on detectable surface immunoglobulins , and it has since been cited in the literature as a B-cell lymphoma . Our immunophenotyping revealed expression of CD1, CD11c, CD18, and MHCII which is the characteristic immunophenotype of canine HS [53, 57]. Since the cells also exhibit low and intermediate expression of CD4 and TCRδγ, respectively, we could not completely rule out a T-cell lymphoma with aberrant expression of CD11c [52, 53], although subsequent genomic results further supported HS. 3132 cells were also T-cell PARR positive, although rearrangements of both T-cell receptor and immunoglobulin genes have been previously identified in human histiocytic sarcomas [1, 58] and in canine histiocytomas . Finally, our karyotypic, CNA, and GEP data also suggests that 3132 cells represent a histiocytic sarcoma. Metaphase spreads from primary cases contain highly variable chromosome numbers with a range of 42–53 and an abundance of aberrant bi-armed chromosomes which is comparable to karyotypic findings of 3132. oaCGH of 3132 cells was indicative of high genomic instability based on the number of CNAs identified, the percentage of the genome involved in regions of CNA, and the observation that several chromosomes include numerous of gains and losses across the chromosome, all of which have been found in primary HS . Several specific CNAs are also shared with primary HS, including loss of CFA 16 and loss of CFA 31. Loss of CFA 12, 14, and 36 and gain of CFA 13 are also conserved with CNAs in primary canine HS . When examining cluster analysis of GEP data with primary canine round cell tumors, 3132 segregated with a large group of histiocytic malignancies.
Dysregulation of several genes previously associated with canine diffuse large B-cell lymphoma in the NF-kB signaling and B-cell receptor signaling pathways shared conserved expression in CLBL-1 including the following: KRAS, NRAS, PIK3R5, PLCG2, TGFBR2, TNFAIP3, TRADD, BCL2A1, CAMK2D, NFATC2 , BUB1B, PRKCD, CD83, CXCL13, CD36, IL8, IL2, CD40LG, LCK, LTBR, and TNFSF11 . Similarities were also noted between CLBL-1 and other targeted B-cell lymphoma gene expression studies, including KIT (almost 400-fold based on qRT-PCR)  and ZAP70 downregulation (decreased 10 fold) .
There are less reported data regarding gene expression changes in other canine hematopoietic malignancies. KIT expression is decreased in canine T-cell LSA , and was decreased three fold in the TL-1 cell line. SYK expression is downregulated in canine T cell malignancies  and was downregulated 12–14 fold in all three T-cell cell lines (Ema, TL-1, and UL-1). MMP9 and TIMP1 expression is significantly upregulated in canine T-cell lymphomas , and was upregulated in both Ema and TL-1 cells. UL-1 was the only cell line to upregulate FLT3, and increased FLT3 expression has been previously identified in canine a subset of acute lymphocytic leukemias and the GL-1 cell line . Upregulation of VEGFA was also found in UL-1 and has been previously reported in canine ALL . Finally, upregulation of GTSF1, LUM, and PYPH and downregulation of CLEC12A and CD9 in primary canine histiocytic sarcomas [63, 64] were found to be dysregulated similarly in 3132 cells.
Cross contamination of cell lines is a common issue in the scientific community with ~19 % of human LL cell lines being cross contaminated . This problem is further evidenced by the identification of human cells in five recently characterized canine B-cell lymphoma cell lines . Our cytogenetic analysis of chromosomal architecture, in conjunction with oaCGH and FISH analysis, proves none of the examined cell lines are cross contaminated with other canine or human cells.
In summary, we present the comprehensive immunophenotypic and genomic characterization of 5 canine LL cell lines that confirms their cell of origin in 3 cell lines (CLBL-1, Ema, and TL-1) and refutes their cell of origin in 2 cell lines (UL-1 and 3132). These data provide valuable information that can be used to select cell lines for preclinical in vitro studies based on the presence or absence of particular immunophenotypic and/or genomic characteristics of interest, which will enhance their clinical predictive value. The generation of integrated molecular profiling of cell lines and comparison with primary tumors will allow further exploration into their biology and clinical utility in veterinary medicine and contribute to comparative and translational studies of hematopoietic malignancies in dogs and humans.
Copy number analysis
Fluorescent in situ hybridization
Gene expression profiling
Oligonucleotide array comparative genomic hybridization
Polymerase chain reaction for antigen receptor rearrangement
Swerdlow S, Campo E, Harris NL. WHO classification of tumours of haematopoietic and lymphoid tissues. France: IARC Press; 2008.
MacLeod RA, Nagel S, Scherr M, Schneider B, Dirks WG, Uphoff CC, Quentmeier H, Drexler HG. Human leukemia and lymphoma cell lines as models and resources. Curr Med Chem. 2008;15(4):339–59.
Ferreira D, Adega F, Chaves R. The importance of cancer cell lines as in vitro models in cancer methylome analysis and anticancer drugs testing. In: López-Camarillo C, Aréchaga-Ocampo E, editors. Oncogenomics and cancer proteomics - novel approaches in biomarkers discovery and therapeutic targets in cancer, InTech. 2013.
Rücker FG, Sander S, Döhner K, Döhner H, Pollack JR, Bullinger L. Molecular profiling reveals myeloid leukemia cell lines to be faithful model systems characterized by distinct genomic aberrations. Leukemia. 2006;20(6):994–1001.
Andersson A, Edén P, Lindgren D, Nilsson J, Lassen C, Heldrup J, Fontes M, Borg Å, Mitelman F, Johansson B. Gene expression profiling of leukemic cell lines reveals conserved molecular signatures among subtypes with specific genetic aberrations. Leukemia. 2005;19(6):1042–50.
Drexler HG, Matsuo Y, MacLeod RA. Continuous hematopoietic cell lines as model systems for leukemia–lymphoma research. Leuk Res. 2000;24(11):881–911.
Breen M, Modiano JF. Evolutionarily conserved cytogenetic changes in hematological malignancies of dogs and humans – man and his best friend share more than companionship. Chromosom Res. 2008;16(1):145–54.
Thomas R, Seiser EL, Motsinger-Reif A, Borst L, Valli VE, Kelley K, Suter SE, Argyle D, Burgess K, Bell J, et al. Refining tumor-associated aneuploidy through ‘genomic recoding’ of recurrent DNA copy number aberrations in 150 canine non-Hodgkin’s lymphomas. Leuk Lymphoma. 2011;52(7):1321–35.
Suter S, Small G, Seiser E, Thomas R, Breen M, Richards K. FLT3 mutations in canine acute lymphocytic leukemia. BMC Cancer. 2011;11(1):38.
Kisseberth WC, Nadella MVP, Breen M, Thomas R, Duke SE, Murahari S, Kosarek CE, Vernau W, Avery AC, Burkhard MJ. A novel canine lymphoma cell line: a translational and comparative model for lymphoma research. Leuk Res. 2007;31(12):1709–20.
Seiser E, Thomas R, Richards K, Kathryn Kelley M, Moore P, Suter S, Breen M. Reading between the lines: molecular characterization of five widely used canine lymphoid tumour cell lines. Vet Comp Oncol. 2011;11(1):30–50.
Comazzi S, Guscetti F, Marconato L. First meeting of the European canine lymphoma group. Workshop: state of the art and comparative aspects in canine lymphoma. CH-Lugano, 22 June 2013. Hematol Oncol. 2013;32(2):68–71.
Aricò A, Giantin M, Gelain ME, Riondato F, Comazzi S, Rütgen BC, Essler SE, Dacasto M, Castagnaro M, Aresu L. The role of vascular endothelial growth factor and matrix metalloproteinases in canine lymphoma: in vivo and in vitro study. BMC Vet Res. 2013;9(1):94.
Fujiwara-Igarashi A, Goto-Koshino Y, Mochizuki H, Sato M, Fujino Y, Ohno K, Tsujimoto H. Inhibition of p16 tumor suppressor gene expression via promoter hypermethylation in canine lymphoid tumor cells. Res Vet Sci. 2014;97(1):60–3.
Fujiwara N, Kawasaki H, Yabe R, Christensen DJ, Vitek MP, Mizuno T, Sato K, Ohama T. A potential therapeutic application of SET/I2PP2A Inhibitor OP449 for Canine T-cell lymphoma. J Vet Med Sci. 2013;75(3):349–54.
Kojima K, Fujino Y, Goto-Koshino Y, Ohno K, Tsujimoto H. Analyses on activation of NF-kB and effect of bortezomib in canine neoplastic lymphoid cell lines. J Vet Med Sci. 2012;75(6):727–31.
Shiomitsu K, Xia X, Waite K, Sehgal I, Li S. Evaluation of the aurora kinase inhibitor, ZM447439, in canine malignant lymphoid cells in vitro. Open J Vet Med. 2013;3:29.
Tomiyasu H, Goto-Koshino Y, Fujino Y, Ohno K, Tsujimoto H. Antitumour effect and modulation of expression of the ABCB1 gene by perifosine in canine lymphoid tumour cell lines. Vet J. 2014;199(1):103–9.
Tomiyasu H, Goto-Koshino Y, Fujino Y, Ohno K, Tsujimoto H. The regulation of the expression of ABCG2 Gene through mitogen-activated protein kinase pathways in canine lymphoid tumor cell lines. J Vet Med Sci. 2013;76(2):237–42.
Tomiyasu H, Goto-Koshino Y, Fujino Y, Ohno K, Tsujimoto H. Epigenetic regulation of the ABCB1 gene in drug-sensitive and drug-resistant lymphoid tumour cell lines obtained from canine patients. Vet J. 2014;199(1):103–9.
Mudaliar MA, Haggart RD, Miele G, Sellar G, Tan KA, Goodlad JR, Milne E, Vail DM, Kurzman I, Crowther D. Comparative gene expression profiling identifies common molecular signatures of NF-kB activation in canine and human diffuse large B cell lymphoma (DLBCL). PLoS One. 2013;8(9):e72591.
Rütgen BC, Hammer SE, Gerner W, Christian M, De Arespacochaga AG, Willmann M, Kleiter M, Schwendenwein I, Saalmüller A. Establishment and characterization of a novel canine B-cell line derived from a spontaneously occurring diffuse large cell lymphoma. Leuk Res. 2010;34(7):932–8.
Umeki S, Ema Y, Suzuki R, Kubo M, Hayashi T, Okamura Y, Yamazaki J, Tsujimoto H, Tani K, Hiraoka H. Establishment of five canine lymphoma cell lines and tumor formation in a xenotransplantation model. J Vet Med Sci. 2012;75(4):467–74.
Holmes MA, Duffus W, Gorman N. Natural cytotoxicity in the dog: description of two new allogeneic tumour targets. Vet Immunol Immunopathol. 1989;23(1):161–70.
Strandstrom H, Rimaila-Parnanen E. Canine atypical malignant lymphoma. Am J Vet Res. 1979;40(7):1033.
Vernau W, Moore P. An immunophenotypic study of canine leukemias and preliminary assessment of clonality by polymerase chain reaction. Vet Immunol Immunopathol. 1999;69(2):145–64.
Burnett R, Vernau W, Modiano J, Olver C, Moore P, Avery A. Diagnosis of canine lymphoid neoplasia using clonal rearrangements of antigen receptor genes. Vet Pathol. 2003;40(1):32–41.
Thalheim L, Williams L, Borst L, Fogle J, Suter S. Lymphoma immunophenotype of dogs determined by immunohistochemistry, flow cytometry, and polymerase chain reaction for antigen receptor rearrangements. J Vet Intern Med. 2013;27(6):1509–16.
Breen M, Hitte C, Lorentzen TD, Thomas R, Cadieu E, Sabacan L, Scott A, Evanno G, Parker HG, Kirkness EF, et al. An integrated 4249 marker FISH/RH map of the canine genome. BMC Genomics. 2004;5(1):65.
Thomas R, Borst L, Rotroff D, Motsinger-Reif A, Lindblad-Toh K, Modiano JF, Breen M. Genomic profiling reveals extensive heterogeneity in somatic DNA copy number aberrations of canine hemangiosarcoma. Chromosom Res. 2014;22(3):305–19.
Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, Rahman N, Stratton MR. A census of human cancer genes. Nat Rev Cancer. 2004;4(3):177–83.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013.
Olshen AB, Venkatraman E, Lucito R, Wigler M. Circular binary segmentation for the analysis of array‐based DNA copy number data. Biostatistics. 2004;5(4):557–72.
Roode SC, Rotroff D, Avery AC, Suter SE, Bienzle D, Schiffman JD, Motsinger-Reif A, Breen M. Genome-wide assessment of recurrent genomic imbalances in canine leukemia identifies evolutionarily conserved regions for subtype differentiation. Chromosom Res. 2015;23(4):681–708.
Thomas R, Duke SE, Karlsson EK, Evans A, Ellis P, Lindblad-Toh K, Langford CF, Breen M. A genome assembly-integrated dog 1 Mb BAC microarray: a cytogenetic resource for canine cancer studies and comparative genomic analysis. Cytogenet Genome Res. 2008;122(2):110–21.
Kapushesky M, Adamusiak T, Burdett T, Culhane A, Farne A, Filippov A, Holloway E, Klebanov A, Kryvych N, Kurbatova N. Gene expression atlas update—a value-added database of microarray and sequencing-based functional genomics experiments. Nucleic Acids Res. 2012;40(D1):1077–81.
Wall M, Campbell L. Cytogenetics of Lymphomas. In: Wiernik PH, Goldman J, Dutcher J, Kyle R, editors. Neoplastic Diseases of the Blood. 5th ed. New York: Springer Science; 2013. p. 945–84.
Zhao W-L, Mourah S, Mounier N, Leboeuf C, Daneshpouy ME, Legrès L, Meignin V, Oksenhendler E, Le Maignin C, Calvo F. Vascular endothelial growth factor-A is expressed both on lymphoma cells and endothelial cells in angioimmunoblastic T-cell lymphoma and related to lymphoma progression. Lab Investig. 2004;84(11):1512–9.
van Delft FW, Horsley S, Colman S, Anderson K, Bateman C, Kempski H, Zuna J, Eckert C, Saha V, Kearney L. Clonal origins of relapse in ETV6-RUNX1 acute lymphoblastic leukemia. Blood. 2011;117(23):6247–54.
Li Z, Zhang H, Chen Y, Fan L, Fang J. Forkhead transcription factor FOXO3a protein activates nuclear factor kB through B-cell lymphoma/leukemia 10 (BCL10) protein and promotes tumor cell survival in serum deprivation. J Biol Chem. 2012;287(21):17737–45.
Li H, Lahti JM, Valentine M, Saito M, Reed SI, Look AT, Kidd VJ. Molecular Cloning and Chromosomal Localization of the Human Cyclin C (CCNC) and Cyclin E (CCNE) genes: deletion of the CCNC gene in human tumors. Genomics. 1996;32(2):253–9.
Simon C, Chagraoui J, Krosl J, Gendron P, Wilhelm B, Lemieux S, Boucher G, Chagnon P, Drouin S, Lambert R. A key role for EZH2 and associated genes in mouse and human adult T-cell acute leukemia. Genes Dev. 2012;26(7):651–6.
Fujiwara S, Yamashita Y, Nakamura N, Choi Y, Ueno T, Watanabe H, Kurashina K, Soda M, Enomoto M, Hatanaka H. High-resolution analysis of chromosome copy number alterations in angioimmunoblastic T-cell lymphoma and peripheral T-cell lymphoma, unspecified, with single nucleotide polymorphism-typing microarrays. Leukemia. 2008;22(10):1891–8.
Costello R, Sanchez C, Le Treut T, Rihet P, Imbert J, Sébahoun G. Peripheral T‐cell lymphoma gene expression profiling and potential therapeutic exploitations. Br J Haematol. 2010;150(1):21–7.
Nagel S, Venturini L, Przybylski GK, Grabarczyk P, Meyer C, Kaufmann M, Battmer K, Schmidt CA, Drexler HG, Scherr M. NK-like homeodomain proteins activate NOTCH3-signaling in leukemic T-cells. BMC Cancer. 2009;9(1):371.
Abdelali RB, Asnafi V, Leguay T, Boissel N, Buzyn A, Chevallier P, Thomas X, Lepretre S, Huguet F, Vey N. Pediatric-inspired intensified therapy of adult T-ALL reveals the favorable outcome of NOTCH1/FBXW7 mutations, but not of low ERG/BAALC expression: a GRAALL study. Blood. 2011;118(19):5099–107.
Richards KL, Motsinger-Reif AA, Chen H-W, Fedoriw Y, Fan C, Nielsen DM, Small GW, Thomas R, Smith C, Dave SS. Gene profiling of canine B-cell lymphoma reveals germinal center and postgerminal center subtypes with different survival times, modeling human DLBCL. Cancer Res. 2013;73(16):5029–39.
Wu Z, Irizarry RA, Gentleman R, Murillo FM, Spencer F. A model based background adjustment for oligonucleotide expression arrays. In: Johns Hopkins University, Dept of Biostatistics Working Papers. 2004. Working Paper 1.
Dennis Jr G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4(5):3.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008;4(1):44–57.
Tsai P-C, Breen M. Array-based comparative genomic hybridization–guided identification of reference genes for normalization of real-time quantitative polymerase chain reaction assay data for lymphomas, histiocytic sarcomas, and osteosarcomas of dogs. Am J Vet Res. 2012;73(9):1335–43.
Ortolani C. Flow cytometry of hematological malignancies. Wiley; 2011.
Moore P. A review of histiocytic diseases of dogs and cats. Vet Pathol. 2014;51(1):167–84.
Breen M. Canine cytogenetics--from band to basepair. Cytogenet Genome Res. 2008;120(1–2):50–60.
Rütgen BC, Willenbrock S, Reimann-Berg N, Walter I, Fuchs-Baumgartinger A, Wagner S, Kovacic B, Essler SE, Schwendenwein I, Nolte I. Authentication of primordial characteristics of the CLBL-1 cell line prove the integrity of a canine B-cell lymphoma in a murine in vivo model. PLoS One. 2012;7(6):e40078.
Reggeti F, Bienzle D. Flow cytometry in veterinary oncology. Veterinary Pathology. 2011;48(1):223–35.
Affolter V, Moore P. Localized and disseminated histiocytic sarcoma of dendritic cell origin in dogs. Vet Pathol. 2002;39(1):74–83.
Hanson CA, Jaszcz W, Kersey JH, Astorga MG, Peterson BA, Gajl‐Peczalska KJ, Frizzera G. True histiocytic lymphoma: histopathologic, immunophenotypic and genotypic analysis. Br J Haematol. 1989;73(2):187–98.
Hedan B, Thomas R, Motsinger-Reif A, Abadie J, Andre C, Cullen J, Breen M. Molecular cytogenetic characterization of canine histiocytic sarcoma: A spontaneous model for human histiocytic cancer identifies deletion of tumor suppressor genes and highlights influence of genetic background on tumor behavior. BMC Cancer. 2011;11(1):201.
Giantin M, Aresu L, Aricò A, Gelain M, Riondato F, Comazzi S, Dacasto M. Evaluation of tyrosine-kinase receptor c-kit mutations, mRNA and protein expression in canine lymphoma: might c-kit represent a therapeutic target? Vet Immunol Immunopathol. 2013;154(3):153–9.
Mortarino M, Gelain M, Gioia G, Ciusani E, Bazzocchi C, Comazzi S. ZAP-70 and Syk expression in canine lymphoid cells and preliminary results on leukaemia cases. Vet Immunol Immunopathol. 2009;128(4):395–401.
Aricò A, Giantin M, Gelain M, Riondato F, Mortarino M, Comazzi S, Dacasto M, Castagnaro M, Aresu L. Matrix metalloproteinases and vascular endothelial growth factor expression in canine leukaemias. Vet J. 2013;196(2):260–2.
Boerkamp KM, van der Kooij M, van Steenbeek FG, van Wolferen ME, Koerkamp MJG, van Leenen D, Grinwis GC, Penning LC, Wiemer EA, Rutteman GR. Gene expression profiling of histiocytic sarcomas in a canine model: the predisposed flatcoated retriever dog. PLoS One. 2013;8(8):e71094.
Boerkamp KM, van Steenbeek FG, Penning LC, Koerkamp MJG, van Leenen D, Vos-Loohuis M, Grinwis GC, Rutteman GR. The two main forms of histiocytic sarcoma in the predisposed flatcoated retriever dog display variation in gene expression. PLoS One. 2014;9(6):e98258.
Zwingenberger AL, Vernau W, Sh C, Yan W, Chen X, Gordon IK, Kent MS. Corrigendum to “Development and characterization of 5 canine B-cell lymphoma cell lines” [Leuk. Res. 36 (5)(2011) 601–606]. Leuk Res. 2013;37(1):120.
We thank Dr. Rachael Thomas and Katie Kennedy for providing oaCGH data for primary canine lymphomas and histiocytic sarcomas, and Linda English (NCSU Clinical Immunology Laboratory) for flow cytometry data.
This study was funded in part by awards made by the North Carolina University Cancer Research Fund administered by the Lineberger Comprehensive Cancer Center and the Comparative Medicine Institute (SES), AKC-Canine Health Foundation (Grant number 01843, MB) and the NCSU Cancer Genomics Fund (MB). SCR was supported in part by an NCSU Comparative Biomedical Science DVM-PhD training fellowship and a NIH T35 training grant.
Availability of data and materials
GEP data has been deposited in the NCBI GEO database (GSE85720).
SRC, MB, and SES conceived, designed, and performed the experiments. DR, AMR, and KLR provided statistical analysis. PM performed the immunophenotypic analysis. SCR, KLR, YO, TM, HT, SES, and MB wrote and edited the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
IACUC approval or owner consent was not required for this study since no cells from living animals were utilized.
Known cancer-associated genes located in regions of genomic imbalance in each of the five cell lines. Table S2. Summary of fold changes of cancer-associated genes for each cell line when normalized expression level was compared with mean normalized expression levels of non-neoplastic lymph nodes. Table S3. List of 50 transcripts from each cell line with the largest absolute fold changes when normalized expression level was compared with mean normalized expression levels of non-neoplastic lymph nodes. Table S4. Gene ontology (GO) biological processes enriched for genes upregulated or downregulatued by >5 fold in each cell line when compared with non-neoplastic lymph nodes. The top ten terms associated with up and down regulated genes are listed. Table S5. KEGG pathways associatated with genes upregulated or downregulated by >5 fold in each cell line when compared with non-neoplastic lymph node. Up to the top ten pathways associated with up or down regulated genes are listed. (XLSX 113 kb)
About this article
Cite this article
Roode, S.C., Rotroff, D., Richards, K.L. et al. Comprehensive genomic characterization of five canine lymphoid tumor cell lines. BMC Vet Res 12, 207 (2016). https://doi.org/10.1186/s12917-016-0836-z