The increasing number of immunotherapies developed in the last two decades presents the need for an appropriate animal model to evaluate the efficacy of these treatments. The spontaneous nature of cancer in dogs and the common features they share with human malignancies make the dog a favorable translational model. The major histocompatibility complex (MHC) molecules in dogs are referred to as dog leukocyte antigens (DLA). Here, we introduce two antibodies for the characterization of the DLA class I immunopeptidome from primary canine tumors. We show that up to 55% of the peptides presented by tumor DLA are identical to peptides reported from common HLA class I molecules, displaying striking similarity in length and anchoring positions. Intriguingly, hundreds of these tumor DLA peptides are derived from well-established cancer-associated antigens. In summary, we demonstrate that canine and human MHC class I molecules are highly homologous in their antigen presentation function and peptide repertoire. These findings exhibit promising implications for advancing cancer immunotherapies and their translation from dogs to humans.

Despite the widespread application of immunotherapy, treating immune-cold tumors remains a significant challenge in cancer therapy. Using multiomic spatial analyses and experimental validation, we identify MGAT1, a glycosyltransferase, as a pivotal factor governing tumor immune response. Overexpression of MGAT1 leads to immune evasion due to aberrant elevation of CD73 membrane translocation, which suppresses CD8+ T cell function, especially in immune-cold triple-negative breast cancer (TNBC). Mechanistically, addition of N-acetylglucosamine to CD73 by MGAT1 enables the CD73 dimerization necessary for CD73 loading onto VAMP3, ensuring membrane fusion. We further show that THBS1 is an upstream etiological factor orchestrating the MGAT1-CD73-VAMP3-adenosine axis in suppressing CD8+ T cell antitumor activity. Spatial transcriptomic profiling reveals spatially resolved features of interacting malignant and immune cells pertaining to expression levels of MGAT1 and CD73. In preclinical models of TNBC, W-GTF01, an inhibitor specifically blocked the MGAT1-catalyzed CD73 glycosylation, sensitizing refractory tumors to anti-PD-L1 therapy via restoring capacity to elicit a CD8+ IFNγ-producing T cell response. Collectively, our findings uncover a strategy for targeting the immunosuppressive molecule CD73 by inhibiting MGAT1.

Triple-negative breast cancer (TNBC) in obese patients remains challenging. Recent studies have linked obesity to an increased risk of TNBC and malignancies. Through multiomic analysis and experimental validation, a dysfunctional Eukaryotic Translation Initiation Factor 3 Subunit H (EIF3H)/Yes-associated protein (YAP) proteolytic axis is identified as a pivotal junction mediating the interplay between cancer-associated adipocytes and the response to anti-cancer drugs in TNBC. Mechanistically, cancer-associated adipocytes drive metabolic reprogramming resulting in an upregulated hexosamine biosynthetic pathway (HBP). This aberrant upregulation of HBP promotes YAP O-GlcNAcylation and the subsequent recruitment of EIF3H deubiquitinase, which stabilizes YAP, thus promoting tumor growth and chemotherapy resistance. It is found that Retatrutide, an anti-obesity agent, inhibits HBP and YAP O-GlcNAcylation leading to increased YAP degradation through the deprivation of EIF3H-mediated deubiquitylation of YAP. In preclinical models of obese TNBC, Retatrutide downregulates HBP, decreases YAP protein levels, and consequently decreases tumor size and enhances chemotherapy efficacy. This effect is particularly pronounced in obese mice bearing TNBC tumors. Overall, these findings reveal a critical interplay between adipocyte-mediated metabolic reprogramming and EIF3H-mediated YAP proteolytic control, offering promising therapeutic strategies to mitigate the adverse effects of obesity on TNBC progression.

The dog serves as a key translational model in cancer immunotherapy. Understanding the T cell receptor (TCR) repertoire is needed for various cancer immunotherapies. Compared to humans where >300 million TCRs have been identified, <100 canine TCRs are reported. To address this deficiency, we assembled >200,000 complete TCR complementarity-determining region 3 (CDR3) sequences from RNA-seq data published for ~2,000 canine samples of blood, lymph node, and other tissues, of which 613 are tumors. We collected 1,324 human RNA-seq samples to compare the similarities and differences in T-cell repertoires between humans and dogs. Notably, our analysis revealed distinct variable gene usage patterns between blood samples and solid tissues in both canine and human samples for TRA and TRB loci. Moreover, our investigation led to the discovery of novel V gene and allele candidates in the canine genome. Our findings also revealed that the canine CDR3 resembled human CDR3 in terms of length and motifs. Additionally, our study unveiled shared traits in cancer TCRs between dogs and humans, including longer lengths and higher hydrophobicity of private CDR3s. Our results indicated the diversity of canine to be more comparable to that of humans than mice. Our study provides an initial landscape of the canine TCR repertoire, highlighting both its similarities and differences with the human counterpart, thus laying the groundwork for future research in comparative immunology and vaccine development.

The genes of the Major Histocompatibility Complex class I (MHC-I) are among the most diverse in the mammalian genome, playing a crucial role in immunology. Understanding the diversity landscape of MHC-I is therefore of paramount importance. The dog is a key translational model in various biomedical fields. However, our understanding of the canine MHC-I diversity landscape lags significantly behind that of humans. To address this deficiency, we used our newly developed software, KPR de novo assembler and genotyper, to genotype 1,325 samples from 1,025 dogs with paired-end RNA-seq data from 43 BioProjects, after extensive quality control. Among 926 dogs that pass the QC, 591 dogs (64%) have at least one allele genotyped, and a total of 97 known alleles and 52 putative new alleles were identified. Further analysis reveals that DLA-I gene expression levels vary among the tissues, with lowest for testis and brain tissues and highest for blood, corpus luteum, and spleen. We identified dominant alleles in each of the 17 canine breeds, as well as among the entire canine population. Furthermore, our analysis also identifies breed-specific alleles and mutually co-occurred/exclusive alleles. Our study indicates that canine DLA-88 is as diversified as human HLA-A/B/C genes within the entire population, but less diversified within a breed than with HLA-A/B/C within an ethnic group. Lastly, we examined the hypervariable regions (HVR) within or across human/canine MHC-I alleles and found that 80% of the HVRs overlap between the two species. We further noted that 80% of the HVRs are within 4A contact with the peptides, and that the dog-human difference overlaps with only 20% HVRs. Our research offers valuable insights for immunological studies involving dogs.

The major histocompatibility complex class I (MHC-I) genes are highly polymorphic. MHC-I genotyping is required for determining the peptide epitopes available to an individual's T-cell repertoire. Current genotyping software tools do not work for the dog, due to very limited known canine alleles. To address this, we developed a Kmer-based paired-end read (KPR) de novo assembler and genotyper, which assemble paired-end RNA-seq reads from MHC-I regions into contigs, and then ge-notype each contig and estimate its expression level. KPR tools outperform other popular software examined in typing new alleles. We used KPR tools to success-fully genotype152 dogs from a published dataset. The study discovers 33 putative new alleles, finds dominant alleles in 4 dog breeds, and builds allele diversity and expression landscapes among the 152 dogs. Our software meets a significant need in biomedical research.

Fig. S1. Supplementary data for large scale genomic aberrations illustrated in Figure 2, including: A) the distribution of mapped read density of WGS; B) PCR confirmation of translocation and fusion junctions of a superamplicon; and C) PCR and RT-PCR amplification of MGAM, ZFAND3, and the ZFAND3-MGAM fusion gene. Fig. S2. The distribution of mapped read density of WES, supplementary data for Figure 3. Fig. S3. Differentially expressed genes and functions among complex, simple carcinomas and normal mammary glands, supplementary data for Figure 3. Fig. S4. IHC images of histone modifications in canine normal mammary glands at 400X magnification, supplementary data for Figure 4. Fig. S5. IHC images of histone modifications in simple carcinomas, complex carcinomas and normal mammary glands at 200X magnification, supplementary data for Figure 5. Fig. S6. Supplementary data for homology between canine MCs and human breast cancers described in Figure 6. These include differentially expressed genes between ER+ and ER- carcinomas at FDR ≤ 0.1 and an example of PAM50 clustering indicating canine and human tumor IDs.

Supplementary Table S3. Supplementary data for mutation discovery in canine MCs by WES and RNA-seq analyses, including read mapping information, mutated gene listed, base substitution type, etc.

Naturally occurring canine cancers have remarkable similarities to their human counterparts. To better understand these similarities, we investigated 671 client-owned dogs from 96 breeds with 23 common tumor types, including those whose mutation profile are unknown (anal sac carcinoma and neuroendocrine carcinoma) or understudied (thyroid carcinoma, soft tissue sarcoma and hepatocellular carcinoma). We discovered mutations in 50 well-established oncogenes and tumor suppressors, and compared them to those reported in human cancers. As in human cancer, TP53 is the most commonly mutated gene, detected in 22.5% of canine tumors overall. Canine tumors share mutational hotspots with human tumors in oncogenes including PIK3CA , KRAS , NRAS , BRAF , KIT and EGFR . Hotspot mutations with significant association to tumor type include NRAS G61R and PIK3CA H1047R in hemangiosarcoma, ERBB2 V659E in pulmonary carcinoma, and BRAF V588E (equivalent of V600E in humans) in urothelial carcinoma. Our findings better position canines as a translational model of human cancer to investigate a wide spectrum of targeted therapies.

Abstract Spontaneously occurring canine mammary cancer represents an excellent model of human breast cancer, but is greatly understudied. To better use this valuable resource, we performed whole-genome sequencing, whole-exome sequencing, RNA-seq, and/or high-density arrays on twelve canine mammary cancer cases, including seven simple carcinomas and four complex carcinomas. Canine simple carcinomas, which histologically match human breast carcinomas, harbor extensive genomic aberrations, many of which faithfully recapitulate key features of human breast cancer. Canine complex carcinomas, which are characterized by proliferation of both luminal and myoepithelial cells and are rare in human breast cancer, seem to lack genomic abnormalities. Instead, these tumors have about 35 chromatin-modification genes downregulated and are abnormally enriched with active histone modification H4-acetylation, whereas aberrantly depleted with repressive histone modification H3K9me3. Our findings indicate the likelihood that canine simple carcinomas arise from genomic aberrations, whereas complex carcinomas originate from epigenomic alterations, reinforcing their unique value. Canine complex carcinomas offer an ideal system to study myoepithelial cells, the second major cell lineage of the mammary gland. Canine simple carcinomas, which faithfully represent human breast carcinomas at the molecular level, provide indispensable models for basic and translational breast cancer research. Cancer Res; 74(18); 5045–56. ©2014 AACR.

Background About 20% of breast cancers in humans are basal-like, a subtype that is often triple-negative and difficult to treat. An effective translational model for basal-like breast cancer is currently lacking and urgently needed. To determine whether spontaneous mammary tumors in pet dogs could meet this need, we subtyped canine mammary tumors and evaluated the dog–human molecular homology at the subtype level. Methods We subtyped 236 canine mammary tumors from 3 studies by applying various subtyping strategies on their RNA-seq data. We then performed PAM50 classification with canine tumors alone, as well as with canine tumors combined with human breast tumors. We identified feature genes for human BLBC and luminal A subtypes via machine learning and used these genes to repeat canine-alone and cross-species tumor classifications. We investigated differential gene expression, signature gene set enrichment, expression association, mutational landscape, and other features for dog–human subtype comparison. Results Our independent genome-wide subtyping consistently identified two molecularly distinct subtypes among the canine tumors. One subtype is mostly basal-like and clusters with human BLBC in cross-species PAM50 and feature gene classifications, while the other subtype does not cluster with any human breast cancer subtype. Furthermore, the canine basal-like subtype recaptures key molecular features (e.g., cell cycle gene upregulation, TP53 mutation) and gene expression patterns that characterize human BLBC. It is enriched in histological subtypes that match human breast cancer, unlike the other canine subtype. However, about 33% of canine basal-like tumors are estrogen receptor negative (ER−) and progesterone receptor positive (PR+), which is rare in human breast cancer. Further analysis reveals that these ER−PR+ canine tumors harbor additional basal-like features, including upregulation of genes of interferon- γ response and of the Wnt-pluripotency pathway. Interestingly, we observed an association of PGR expression with gene silencing in all canine tumors and with the expression of T cell exhaustion markers (e.g., PDCD1 ) in ER−PR+ canine tumors. Conclusions We identify a canine mammary tumor subtype that molecularly resembles human BLBC overall and thus could serve as a vital translational model of this devastating breast cancer subtype. Our study also sheds light on the dog–human difference in the mammary tumor histology and the hormonal cycle.

Supplementary Table S7. Supplementary data indicating altered genes and their function information.

Pet dogs develop spontaneous cancers at a rate estimated to be five times higher than that of humans, providing a unique opportunity to study disease biology and evaluate novel therapeutic strategies in a model system that possesses an intact immune system and mirrors key aspects of human cancer biology. Despite decades of interest, effective utilization of pet dog cancers has been hindered by a limited repertoire of necessary cellular and molecular reagents for both in vitro and in vivo studies, as well as a dearth of information regarding the genomic landscape of these cancers. Recently, many of these critical gaps have been addressed through the generation of a highly annotated canine reference genome, the creation of several tools necessary for multi-omic analysis of canine tumours, and the development of a centralized repository for key genomic and associated clinical information from canine cancer patients, the Integrated Canine Data Commons. Together, these advances have catalysed multidisciplinary efforts designed to integrate the study of pet dog cancers more effectively into the translational continuum, with the ultimate goal of improving human outcomes. The current review summarizes this recent progress and provides a guide to resources and tools available for comparative study of pet dog cancers.

The major histocompatibility complex class I (MHC-I) genes are highly polymorphic among individuals. MHC-I genotyping is required for determining the antigen-binding specificity of each MHC-I molecule in an individual. Numerous tools have been developed for human MHC-I genotyping using deep sequencing data such as RNA-seq; however they do not work for the dog, due to very limited information for canine alleles. To address this issue, we developed a Kmer-based paired-end read (KPR) de novo assembler and genotyper, which first assemble paired-end RNA-seq reads mapped to the MHC-I regions into contigs de novo and then genotype each contig. Our KPR tools are validated by Sanger sequencing, simulation and published genotype data. Applying our KPR tools on the published RNA-seq data of 158 tumor and 64 normal samples from 158 dogs, we have achieved a genotyping success rate of 86%, which includes 133 tumor and 57 normal samples from 142 dogs. We have identified 39 known alleles and 83 new alleles of high confidence in these dogs, yielding a more comprehensive MHC-I allele diversity landscape for the dog.### Competing Interest StatementThe authors have declared no competing interest.

Spontaneous canine cancers are a valuable but relatively understudied and underutilized model in cancer research. To enhance their usage, we reanalyzed whole exome sequencing data published for 601 dogs with mammary cancer, osteosarcoma, oral melanoma, lymphoma, glioma or hemangiosarcoma from over 35 breeds, after rigorous quality control, including breed validation. Each cancer type harbors distinct molecular features, with major pathway alterations matching its human counterpart (e.g., PI3K for mammary cancer and p53 for osteosarcoma). On average, mammary cancer and glioma have lower mutation rates (median <0.5 mutation per Mb), whereas oral melanoma, osteosarcoma and hemangiosarcoma have higher mutation rates (median ≥1 mutation per Mb). Across cancer types and across breeds, the mutation rate is strongly associated with TP53 mutation, but not with PIK3CA mutation. The mutation rate is also associated with a mutation signature enriched in osteosarcoma of Golden Retrievers, independent of TP53 mutation. Finally, compared to other breeds examined, DNA repair genes appear to be less conserved in Golden Retriever which is predisposed to numerous cancers.

Naturally occurring canine cancers have remarkable similarities to their human counterparts. In order to determine whether these similarities occur at the molecular level, we investigated hotspot mutations in a variety of spontaneously arising canine cancers and found high concordance in oncogenic drivers between cancers in both species. These findings suggest that canines may present a powerful and complementary model for preclinical investigations for targeted cancer therapeutics. Through analysis of 708 client-owned dogs from 96 breeds (plus mixed breeds) with 23 common tumor types, we discovered mutations in 50 well-established oncogenes and tumor suppressors, and compared them to those reported in human cancers. TP53 is the most commonly mutated gene, detected in 30.81% of canine tumors overall and >40% in hemangiosarcoma and osteosarcoma. Canine tumors share mutational hotspots with human tumors in oncogenes including PIK3CA, KRAS, NRAS, BRAF, KIT and EGFR. Hotspot mutations with significant (P<0.0001) association to tumor type include NRAS G61R and PIK3CA H1047R in hemangiosarcoma, ERBB2 V659E in pulmonary carcinoma, and BRAF V588E in urothelial carcinoma. This work positions canines as excellent spontaneous models of human cancers that can help to investigate a wide spectrum of targeted therapies.

Madin-Darby canine kidney II (MDCKII) cells are widely used to study epithelial morphogenesis. To better understand this process, we performed time course RNA-seq analysis of MDCKII 3D cystogenesis, along with polarized 2D cells for comparison. Our study reveals a biphasic change in the transcriptome that occurs after the first cell cycle and coincides with lumen establishment. This change appears to be linked to translocation of β-catenin, supported by analyses with AVL9- and DENND5A-knockdown clones, and regulation by HNF1B, supported by ATAC-seq study. These findings indicate a qualitative change model for transcriptome remodeling during epithelial morphogenesis, leading to cell proliferation decrease and cell polarity establishment. Furthermore, our study reveals that active mitochondria are retained and chromatin accessibility decreases in 3D cysts but not in 2D polarized cells. This indicates that 3D culture is a better model than 2D culture for studying epithelial morphogenesis.

Spontaneous tumors in pet dogs represent a valuable but undercharacterized cancer model. To better use this resource, we performed an initial global comparison between proliferative and invasive colorectal tumors from 20 canine cases, and evaluated their molecular homology to human colorectal cancer (CRC). First, proliferative canine tumors harbor overactivated WNT/β-catenin pathways and recurrent CTNNB1 (β-catenin) mutations S45F/P, D32Y and G34E. Invasive canine tumors harbor prominent fibroblast proliferation and overactivated stroma. Both groups have recurrent TP53 mutations. We observed three invasion patterns in canine tumors: collective, crypt-like and epithelial–mesenchymal transition (EMT). We detected enriched Helicobacter bilis and Alistipes finegoldii in proliferative and crypt-like tumors, but depleted mucosa-microbes in the EMT tumor. Second, guided by our canine findings, we classified 79% of 478 human colon cancers from The Cancer Genome Atlas into four subtypes: primarily proliferative, or with collective, crypt-like or EMT invasion features. Their molecular characteristics match those of canine tumors. We showed that consensus molecular subtype 4 (mesenchymal) of human CRC should be further divided into EMT and crypt-like subtypes, which differ in TGF-β activation and mucosa-microbe content. Our canine tumors share the same pathogenic pathway as human CRCs. Dog-human integration identifies three CRC invasion patterns and improves CRC subtyping.

The majority of human cancers are derived from epithelial tissues. Understanding epithelial morphogenesis is hence important in cancer research. Madin-Darby canine kidney II (MDCKII) cells, cultured in 3D and 2D conditions, are used extensively as a model for studying cell polarity, epithelial morphogenesis and carcinogenesis. Like other cell differentiation systems, gene expression plays a key role in MDCKII epithelial morphogenesis.To better understand the change in the transcriptome, we performed a time course RNA-seq analysis of MDCKII 3D cytogenesis, along with fully polarized 2D cells for comparison. Surprisingly, our study reveals that the change is not linear, but rather clearly biphasic. Specifically, about 3,000 genes are significantly up- or down-regulated between 24 hour and day 3 after seeding, when the lumen is forming, compared to < 120 such genes during other time intervals. Because previous studies have shown that Rab11-coordinated intracellular trafficking plays an essential role in MDCKII lumen formation, we hypothesize this qualitative change in the transcriptome is linked to intracellular trafficking. To test this hypothesis, we used knockdown clones of AVL9 and DENND5A, both interacting with Rab11 and participating in trafficking. Knockdown cells have altered cell polarity and defective trafficking, but not as significant changes in transcriptome as wild type cells, supporting our hypothesis. We then focused on β-catenin to better understand the mechanism. Our study reveals that in wild type cells, following the first cell division, β-catenin is depleted from the nucleus and enriched in cell-cell junction. This results in MYC transcriptional silencing (supported by ATAC-seq analysis), which in turn initiated down-regulation of numerous MYC target genes. These observations are however not observed in knockdown cells.Our study supports the qualitative change (“switch on-off”) model, rather than the quantitative (gradual) change model, in transcriptome remodeling during epithelial cell differentiation. Moreover, our work supports that intracellular trafficking likely initiates this quantitative change in transcriptome.Citation Format: Tianfang Wang, Shaying Zhao, Shi-Yuan Cheng. A qualitative change in transcriptome during MDCKII 3D epithelial morphogenesis is linked to intracellular trafficking [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 5236.

Carcinomas have two broadly distinct modes of invasion: collective and individual. Collective invasion is characterized as groups of cells that migrate while retaining cell-cell contacts. Individual invasion occurs through mechanisms such as the acquiring of mesenchymal-like features by cancer cells. While these different modes of invasion have been well studied at the genomic and transcriptional levels, their metabolic alterations remain much less understood. To address this deficiency, we performed untargeted NMR metabolomics, in conjunction with RNA-seq, on two variants of MCF7, a widely used breast cancer cell line. We report herein several findings. First, the two MCF7 variants, one parental (referred to as MCF-7 hereafter) and another having undergone many passages (referred to as Augusta hereafter), differ morphologically, with MCF-7 being more epithelial and Augusta being more mesenchymal-like. Importantly, fibroblast co-culture experiments indicate that MCF-7 cells resemble collective invasion, while Augusta cells resemble individual invasion. Furthermore, both scratch-assay and gene set enrichment analysis of RNA-seq data indicate that Augusta cells are more invasive. Second, our NMR metabolomics analysis reveals that branched chain amino acids (BCAAs) represent a significant metabolic difference between MCF-7 and Augusta cells. Partial least squares discriminant analysis of the metabolomics data identified BCAAs as essential distinguishing features between the two cell lines. Specifically, compared to MCF-7 cells, BCAAs are significantly depleted inside Augusta cells, consistent with the upregulation of several key BCAA degradation genes. We hypothesize that Augusta cells can more efficiently utilize BCAAs, and are performing 13C-labeled BCAA flux analysis to test this hypothesis. Previous studies have shown that an aberrantly expressed BCAA metabolism gene sustains growth and enhances invasiveness of breast cancer cells. However, increased catabolism of BCAAs has not been reported in literature on this subject. Our simple cell line model allows for a deeper molecular understanding of BCAA metabolism in breast cancer cell invasion.Citation Format: Maxwell B. Colonna, Arthur S. Edison, Shaying Zhao. A potential role of branched chain amino acid metabolism in breast cancer cell invasiveness [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 5270.