Canine lymphoma (CL) is one of the most common malignant tumors in dogs. The cause of CL remains unclear. Genetic mutations that have been suggested as possible causes of CL are not fully understood. Whole-exome sequencing (WES) is a time- and cost-effective method for detecting genetic variants targeting only protein-coding regions (exons) that are part of the entire genome region. A total of eight patients with B-cell lymphomas were recruited, and WES analysis was performed on whole blood and lymph node aspirate samples from each patient. A total of 17 somatic variants (GOLIM4, ITM2B, STN1, UNC79, PLEKHG4, BRF1, ENSCAFG00845007156, SEMA6B, DSC1, TNFAIP1, MYLK3, WAPL, ADORA2B, LOXHD1, GP6, AZIN1, and NCSTN) with moderate to high impact were identified by WES analysis. Through a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of 17 genes with somatic mutations, a total of 16 pathways were identified. Overall, the somatic mutations identified in this study suggest novel candidate mutations for CL, and further studies are needed to confirm the role of these mutations.

Keywords:

Subject: Medicine and Pharmacology - Veterinary Medicine

1. Introduction

Advances in next-generation sequencing technologies have made it easy and inexpensive to generate large amounts of genomic data. Determining the sequence of the entire genome is called whole-genome sequencing (WGS), and determining the sequence of the entire exon is called whole-exome sequencing (WES) [1]. Almost all protein-coding genes have discontinuous structures. Protein-coding regions are fragmented into several species called exons [2]. Exomes represent only approximately 1% of the genome; therefore, WES is less expensive than WGS [3]. WES technology is a proven method for identifying functionally relevant genetic variants in diseases such as cancers [4].

Studies using WES in human lymphomas have shown that significantly mutated genes such as CD79B, TP53, CARD11, MYD88, and EZH2 are associated with large B-cell lymphoma (BCL) [5,6,7]). A study using WES in three breeds of dogs (Cocker Spaniel, Golden Retriever, and Boxer) predisposed to canine lymphoma (CL) confirmed mutations in TRAF3-MAP3K14, FBXW7, and POT1 [8]. The mutations identified using WGS studies in CL were in ST6GALNAC5, ENSCAFG00000007370, PPP2CB, TP53, SH2B3, ZNF503, SETD2, and COX18 [9].

The purpose of this study is to identify somatic mutations in dogs with BCL by pairwise WES of DNA from eight dogs to determine the impact of variants from lymph node (LN) aspirate samples compared with whole blood (WB) samples. In addition, a bioinformatic Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using somatic variant genes identified in this study was performed to elucidate their roles in CL.

2. Materials and Methods

2.1. Samples

From 2016 to 2020, among the samples requested for cytological and PCR for antigen receptor rearrangements (PARR) tests for the diagnosis of lymphoma at the Laboratory of Veterinary Laboratory Medicine, Chungbuk National University, samples containing WB were used for WES analysis. The inclusion criteria for the study were multicentric lymphoma with enlarged LNs and no previous or current cancer diagnosis other than lymphoma. Of the eight subjects enrolled in this study, two were Maltese, two were Welsh Corgis, one was a Cocker spaniel, one was a Shih Tzu, one was a White Terrier, and one was a mixed breed. The detailed information of the eight dog patients is presented in Table 1. All had BCL as determined by cytology and PARR.

2.2. Immunophenotyping

Immunophenotyping of all samples was determined by PARR. DNA was extracted from LN aspirates using a MagPurix® Tissue DNA Extraction Kit and MagPurix® 12s automated nucleic acid purification system (Zinexts Life Science Corp., Taiwan) according to the manufacturer’s instructions. DNA preparations were stored at -80 °C until use. For this extraction, fine-needle aspiration (FNA) needles were washed with 600 μL of phosphate-buffered saline. The collected samples were vortexed shortly, and 200 μL was used for extraction. The final volume of the DNA elution was 50 μl. PCR was performed using primers used for amplification of Cµ (positive control), TCRγ CDR3, and Ig CDR3 (Table 2) as previously described [10]. Using distilled water, a negative control was run on each sample to ensure that no contamination was present. The amplification was performed in TaKaRa PCR Thermal Cycler Dice® Touch (TaKaRa Bio, Japan). The PCR reaction conditions for each product (Cµ, IgH major, IgH minor, and TCRγ) are shown in Table 3. The PCR products were separated by electrophoresis and detected using the Qsep 100 automatic nucleic acid protein analysis system (BiOptic, Taiwan). The DNA samples were processed with a premade cartridge supplied with a DNA size marker and dye (Standard Cartridge Cat. No: C105201). According to the manual, all peaks and alleles were analyzed using Qsep 100 software.

2.3. DNA extraction and sequencing

DNA was extracted from WB and LN aspirates using the MagPurix Blood DNA Extraction kit and MagPurix Tissue DNA extraction kit (Zinexts Life Science Corp., Taiwan), respectively, according to the manufacturer’s instructions. The sample quantity and purity were assessed with a NanoDrop spectrophotometer (260/280: 1.6 - 2.3, 260/230 > 1.6), electrophoresis (No DNA degradation), and Qubit fluorometric quantitation (≥ 100 ng/μl). Extracted DNA samples were sent to Theragen Bio (Gyeonggi, Korea) for NovaSeq 6000 (Illumina Inc., San Diego, CA, USA) sequencing with qPCR and SureSelect XT Canine All Exon V2 kit (Agilent, Santa Clara, CA, USA).

2.4. Variant calling

The raw sequencing data in the FASTQ format were subjected to a quality check stage using FastQC v0.10.1 [11] program, and the adapter sequences were removed using Cutadapt v1.81 [12]. The processed data were aligned to a canine reference genome from the Boxer breed (CanFam3.1) using BWA v0.7.17 [13] and were sorted and marked for duplicates using Picard. Somatic variants were called using VarScan v2.3.9 using the alignment data from blood or oral samples as the control. Variant effects were predicted using VEP Ensembl web interface v107 [14] with ROS Cfam 1.0 as the reference. Variant effects such as LOF, disruptive insertions and deletions, stop gain and loss, and splice region variants were classified as high-impact variants, while missense variants, untranslated region (UTR) variants, in-frame insertions, and deletions were classified as moderate-impact variants. Synonymous, intron, and intergenic variants were classified as low-impact variants [9]. Variants with less than or equal to 5% variant frequency in blood samples and greater than or equal to 15% variant frequency in tumor samples were considered tumor-specific according to the VarScan guidelines [15], while variants with greater than or equal to 15% variant frequency in blood were classified as germline variants. A schematic diagram of the pipeline analysis is shown in Figure 1.

2.5. Sharing analysis

Identification of individual variants and genes targeted by said variants that were commonly seen among subjects was performed using an in-house Python 3.9 code that uses VEP files as input.

2.6. Protein-protein interaction (PPI) network construction ana analysis of modules

The PPI network of differentially expressed genes (DEGs) was obtained from the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING, Version 10.0; https://www.strind-db.org/) online database, and a confidence score ≥ 0.4 with FDR stringency of 5% was set as the cutoff criterion, which is the default setting. The PPI network was imported into Cytoscape 3.10.0 (https://cytoscape.org/) software for network visualization.

2.7. KEGG pathway analysis

The gene names that were obtained in this study were converted to KEGG gene IDs according to UniProt (https://uniprot.org/). Then, the KEGG pathway analysis was performed on a KEGG mapper (https://www.genome.jp/kegg/mapper/color.html) using the obtained KEGG gene IDs.

3. Results

3.1. Variant level analysis reveals highly shared GOLIM4 variant

Out of the myriad of somatic variants identified within each subject, variants with a higher likelihood of being directly associated with BCL were selected because the chances of a random mutation appearing at the exact location in the exact same way in multiple subjects are slim. Among 480 variants that were shared by at least two subjects, 19 were classified as moderate- or high-impact variants by VEP. By mapping these shared somatic variants onto the chromosome structure of canines, we identified multiple chromosomes with multiple somatic mutations shared by more than three subjects. In particular, chromosome 25 had consecutive highly shared variants near its edge, and chromosome 34 had a somatic variant shared by half of the subjects (Figure 2).

We additionally attempted to identify shared somatic variants with moderate to high predicted impact and the genes that they target. We found that the highly shared variant in chromosome 34 was a variant in the 3' UTR region of the gene Golgi integral membrane protein 4 (GOLIM4), which plays a role in the endosome to Golgi protein trafficking pathway. Other notable variants include frameshift variants in gene desmocollin1 (DSC1), lipoxygenase homology domains 1 (LOXHD1), and glycoprotein VI platelet (GP6), a start lost variant in gene Adenosine A2B receptor (ADORA2B), and a stop lost variant in the gene antizyme inhibitor 1 (AZIN1) (Table 4).

3.2. Gene level sharing analysis shows variant accumulation in a specific PPI network

After investigating the individual shared variants, we investigated genes that were commonly targeted by somatic variants in multiple subjects although they were not identical variants. We identified 2,131 genes that were commonly targeted, and more than 150 genes were shared by more than half of the subjects. Some examples of chromosomes containing many of the highly shared genes include chromosomes 1, 2, 13, and 15 (Figure 3A). By filtering out 48 genes that were commonly targeted by somatic variants with high or moderate impact, we observed that GOLIM4 was the target of many more 3' UTR variants that differed from the previously identified locus. Other genes that were highly shared include wings apart-like protein homolog (WAPL), which also had 3' UTR variants, and C-C motif chemokine ligand 23 (CCL23), which had missense and 3' UTR variants. (Figure 3B).

To identify PPI networks regarding the shared variants, we analyzed 333 highly shared annotated genes targeted by tumor variants that were shared by greater than or equal to half of the subjects and found that there was a large network of interconnected proteins with over 50 nodes that were affected by them as well as additional smaller networks (Figure 3C).

3.3. KEGG pathway analysis results

As a result of KEGG pathway analysis of 17 genes with somatic mutations, a total of 16 pathways (vascular smooth muscle contraction, calcium signaling pathway, platelet activation, ECM-receptor interaction, oxytocin signaling pathway, regulation of actin cytoskeleton, axon guidance, apelin signaling pathway, alcoholism, cGMP-PKG signaling pathway, gastric acid secretion, Alzheimer disease, focal adhesion, notch signaling pathway, neuroactive ligand-receptor interaction, and Rap1 signaling pathway) were identified (Table 5). After excluding pathways associated with only one gene, the following pathways were identified: vascular smooth muscle contraction, calcium signaling pathway, and platelet activation (Figure 4, Figure 5 and Figure 6).

3.4. Highly shared germline mutations in diffuse large BCL (DLBCL)-associated genes and pathways

Analysis of germline variants present in 12 known DLBCL-associated genes (MYC, BCL2, BCL6, EZH2, GNA13, PTEN, TP53, TNFAIP3, PRDM1, CD79B, CARD11, MYD88 [5,6,7,16,17,18,19], and their PPI networks annotated via string DB revealed several stereotypic moderate- to high-impact variants (Figure 7). Among these cases, TRRAP, EP300, USP7, and ATM showed variants that could greatly impact the function of proteins such as stop gains, stop losses, and frameshifts. In particular, for USP7, the same germline variant causing a stop gain was present in five of the eight subjects. Several missense variants including those present in TAX1BP1, PRMT5, and SYK were shared by all or seven out of the eight subjects involved, indicating a strong possibility that such mutations in the germline may be associated with DLBCL.

4. Discussion

WES analysis of eight patients with BCL in this study identified 17 shared somatic variants with moderate or high impact (GOLIM4, ITM2B, STN1, UNC79, PLEKHG4, BRF1, ENSCAFG00845007156, SEMA6B, DSC1, TNFAIP1, MYLK3, WAPL, ADORA2B, LOXHD1, GP6, AZIN1, and NCSTN). According to the literature, considering the functions of these genes and related diseases, GOLIM4, ITM2B, STN1, DSC1, TNFAIP1, WAPL, and NCSTN can be assumed to play roles as tumor suppressor genes. GOLIM4 (also known as GPP130) is a component of the Golgi transport complex that plays an important role in the transport of Golgi proteins [15]. Although the Golgi apparatus may be involved in tumor biological processes, the function of GOLIM4 during tumorigenesis remains unclear. Nevertheless, the Golgi apparatus and endosome dysfunction are involved in the progression of various tumors, and increased expression of GOLIM4 has been shown to inhibit cancer cell proliferation, promote apoptosis, and induce G1 phase arrest in human head and neck cancer cell lines [16]. Integral membrane protein 2B (ITM2B; also known as BRI2) is a type II transmembrane protein that is a substrate for regulated intramembrane proteolysis [17]. ITM2B induces apoptosis and inhibits proliferation [18]. Downregulation of ITM2B has been identified in human lung cancer tissues; therefore, ITM2B appears to play a role as a tumor suppressor gene [19]. STN1, along with CTC1 and TEN1, is a component of the CST (CTC1-STN1-TEN1) complex and is responsible for maintaining telomere and genome integrity [20]. CST was first identified as a telomere-binding protein complex and functions in telomere replication and protection. CST can mediate end protection at double-strand breaks, likely using a similar strategy to filling in of the telomeric C-strand. Supporting this observation, CST has been shown to promote ploy polymerase inhibitor sensitivity in BRCA1-deficient cancer cells. Given its essential roles in replication and DNA repair, CST is known to be important for genome stability [21]. DSC1 was predicted to encode a sodium channel based on a high sequence similarity with vertebrate and invertebrate sodium channel genes. In human medicine, decreased expression of DSC1 was related to the poor differentiation and prognosis of head and neck squamous cell carcinoma, lung cancer, melanoma, and colorectal carcinoma [22,23,24,25]. TNF-α-induced protein 1 (TNFAIP1; B12) gene is highly conserved gene in several species and is known as a tumor suppressor gene. TNFAIP1 is induced by TNF-α and interleukin-6, and it is mainly involved in DNA synthesis and repair and apoptosis [26,27,28]. In human medicine, breast cancer, gastric carcinomas, and lung cancer are known to be associated with TNFAIP1 mutations [29,30,31]. WAPL is important in regulating the level of aggregation of chromosomes by separating cohesive loops from chromatin. By separating cohesion from chromatin, WAPL is a regulator of the loading and unloading cycle. The loss of WAPL is known to result in p53-dependent cell cycle arrest [32,33]. The role of nicastrin (NCSTN) remains unknown. However, NCSTN is known to be related to AKT and p-AKT, which affect cell proliferation, growth, and differentiation [34]. In addition, incomplete expression of NCSTN is known to reduce the expression of miR-100-5p, which acts as a tumor suppressor involved in cell self-renewal and wound healing [35].

On the other hand, BRF1, SEMA6B, ADORA2B, GP6, and AZIN1 can be assumed to function as oncogenes. The BRF1 gene encodes the BRF1 protein in its zinc ribbon domain and directly participates in the process of protein synthesis. Deregulation of BRF1 is associated with cell proliferation, cell transformation, and tumorigenesis. BRF1 is overexpressed in hepatocellular carcinoma, breast cancer, gastric cancer, prostate cancer, and lung cancer in humans [36,37,38,39]. Semaphorin 6b (SEMA6B) is a member of the semaphoring axon-guidance family and was initially characterized as an axon guidance factor with axon navigation functions but has also been demonstrated to induce or inhibit tumor progression. The overexpression of this gene is related to colorectal cancer in humans [40,41,42]. ADORA2B encodes a protein belonging to the G protein-coupled receptor superfamily that plays a role in tissue distribution along with A1, A2A and A3. Abnormal expression of ADORA2B may play a pathophysiological role in some human cancers [43]. ADORA2B is highly expressed in oral cancer, lung adenocarcinoma, and prostate cancer and promotes proliferation and metastasis of carcinoma cells [43,44,45]. GP6 is a transmembrane protein that is the major signaling receptor for collagen on platelets and regulates several platelet functions, such as adhesion, aggregation, and procoagulant activity [46,47]. In addition, GP6 plays a role in supporting platelet adhesion to tumor cells, which is known to be involved in the metastasis of colorectal cancer and breast cancer [48]. Adenosine to inosine (A-to-I) RNA editing catalyzed by adenosine deaminases acting on RNA enzymes is a post-transcriptional modification that has emerged as a key player in tumorigenesis and cancer progression. AZIN1 has been identified as one of the most frequently occurring A-to-I RNA alterations in colorectal cancer and hepatocellular carcinoma and acts as an oncogene [49,50].

The function of three genes (UNC79, PLEKHG4B, and ENSCAFG00845007156) and their associations with cancers have not been clearly identified. UNC79 protein forms an NALCN complex with NALCN, FAM155, and UNC80 proteins, which are involved in voltage-gated sodium and calcium channels [51,52]. There is no research about the association between UNC79 mutations and cancers. Mutations in the pleckstrin homology domain-containing family g member 4b (PLEKHG4B; puratrophin-1) gene are associated with the hereditary neurological disorder autosomal dominant spinocerebellar ataxia. However, the biochemical function of this gene product has not been described [53]. Moreover, there is no research about the association between PLEKHG4B mutations and cancers. ENSCAFG00845007156 is similar to human Aldo-keto reductase family 1 member D1, but its function has not been identified, and protein analysis has not been performed to date.

For both genes (MYLK3 and LOXHD1), the function of each gene was not associated with cancers but was associated with diseases other than cancers. Myosin light chain kinase 3 (MYLK3) is a protein-coding gene that acts as a regulator of the actin cytoskeleton and immune response signaling. WES revealed that MYLK3 mutations are associated with dilated cardiomyopathy in humans [54]. LOXHD1 encodes a protein consisting of 15 polycystin lipoxygenase α-toxin repeats, which can bind lipids and proteins in other proteins [55]. Mutations in LOXHD1 cause progressive hearing loss (Grillet et al., 2009). As a result, mutations in tumor suppressor genes can affect cancers such as CL. Therefore, GOLIM4, ITM2B, STN1, DSC1, TNFAIP1, WAPL, and NCSTN mutations may be prognostic markers for patients with CL. However, the functions of UNC79, PLEKHG4, and ENSCAFG00845007156 and their associations with cancers have not yet been identified. In addition, MYLK3 and LOXHD1 mutations appear to have a very low association with CL.

As a result of constructing a PPI network for approximately330 genes targeting 480 shared somatic variants, one large network was formed. A total of seven genes (STN1, AZIN1, ITM2B, ADORA2B, SEMA6B, NCSTN, and DSC1) in the network were genes identified with moderate to high impact mutations. However, none of the many genes observed in the network appear to be directly related to CL.

KEGG pathway analysis revealed that three pathways are associated with at least two mutated genes. Among them, two pathways suspected to be closely related to BCL are the calcium signaling pathway and the platelet activation pathway. Processes such as cell proliferation and death and gene transcription are essential for regulating cellular functions, and tight regulation of calcium signaling is fundamental in this process [57]. Therefore, changes in calcium signaling can cause various diseases including tumors, and these changes have been confirmed in cancer cell lines [58]. Platelet function plays an important role not only in hemostasis but also in tumor metastasis [59]. Some research has demonstrated reduced tumor metastasis after experimentally inducing thrombocytopenia in mouse models [60]. Therefore, mutations in genes (ADORA2B, MYLK3, and GP6) involved in both pathways may be mutations of interest in BCL.

Overall, considering the function of each gene, GOLIM4, ITM2B, STN1, DSC1, TNFAIP1, WAPL, and NCSTN mutations were found to be highly associated with BCL, and ADORA2B, MYLK3, and GP6 mutations were suspected to be associated with BCL through KEGG pathway analysis. Therefore, GOLIM4, ITM2B, STN1, DSC1, TNFAIP1, WAPL, NCSTN, ADORA2B, MYLK3, and GP6 mutations are proposed as candidate mutations associated with BCL.

There are some limitations in this study. First, the patients with BCL in this study were various canine breeds. Different breeds of dogs exhibit different physical characteristics, and these characteristics may be related at a genetic level. Second, the number of patients used in this study was small. The inclusion of more patients will more clearly identify commonly observed somatic variants. Third, a number of germline mutations were observed, but these could not be analyzed due to the vast amount of data. A review of germline mutations in human DLBCL suggests a possible association between germline mutations and lymphoma [61]; therefore, further studies on germline mutations and CL will be needed.

5. Conclusions

This study demonstrates the utility of WES in identifying somatic mutations in dogs with BCL. The identification of shared variants and their associated genes contributes to the understanding of the molecular mechanisms underlying BCL development and progression in dogs. Further analysis using bioinformatic tools, such as PPI network construction and KEGG pathway analysis, may provide additional insights into the functional roles of these genes and their involvement in CL. These findings pave the way for future research focusing on targeted therapies and personalized medicine for dogs with BCL.

Author Contributions

Conceptualization, S.K. and KJ.N.; methodology, S.K., N.K. and KJ.N.; software, S.K., N.K. and KJ.N.; validation, S.K., N.K. and KJ.N.; formal analysis, S.K., N.K. and KJ.N.; investigation, S.K., HM.K., and KJ.N.; resources, S.K., HM.K., HJ. J. and KJ.N.; data curation, S.K., N.K. and KJ.N.; writing—original draft preparation, S.K., N.K. and KJ.N.; writing—review and editing, S.K., N.K. and KJ.N.; visualization, S.K. and N.K.; supervision, KJ.N.; project administration, KJ.N.; funding acquisition, KJ.N. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (NRF-2016R1D1A1B03932312) and Hyundai Car Chung Mong-Koo Foundation.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank aniDAP company and Haemaru Referral Animal Hospital for the provision of samples.

Conflicts of Interest

The authors declare no conflict of interest.

References

Bick, D.; Dimmock, D. Whole Exome and Whole Genome Sequencing. Curr Opin Pediatr 2011, 23, 594–600. [Google Scholar] [CrossRef] [PubMed]
Sastre, L. Exome Sequencing: What Clinicians Need to Know. Adv Genom Genetics 2014, 4, 15–27. [Google Scholar] [CrossRef]
Teer, J.K; Mullikin, J.C. Exome Seqeuncing: The Sweet Spot before Whole Genomes. Hum Mol Genet 2010, 19, 145–151. [Google Scholar] [CrossRef] [PubMed]
Kunstman, J.W.; Juhlin, C.C.; Goh, G.; Brown, T.C.; Stenman, A.; Healy, J.M.; Rubinstein, J.C.; Choi, M.; Kiss, N.; Nelson-Williams, C.; et al. Characterization of the Mutational Landscape of Anaplastic Thyroid Cancer via Whole-Exome Sequencing. Hum Mol Genet 2015, 24, 2318–2329. [Google Scholar] [CrossRef] [PubMed]
Davis, R.E.; Ngo, V.N.; Lenz, G.; Tolar, P.; Young, R.M.; Romesser, P.B.; Kohlhammer, H.; Lamy, L.; Zhao, H.; Yang, Y.; et al. Chronic Active B-Cell-Receptor Signalling in Diffuse Large B-Cell Lymphoma. Nature 2010, 463, 88–92. [Google Scholar] [CrossRef]
Ngo, V.N.; Young, R.M.; Schmitz, R.; Jhavar, S.; Xiao, W.; Lim, K.-H.; Kohlhammer, H.; Xu, W.; Yang, Y.; Zhao, H.; et al. Oncogenically Active MYD88 Mutations in Human Lymphoma. Nature 2011, 470, 115–119. [Google Scholar] [CrossRef]
Lenz, G.; Davis, R.E.; Ngo, V.N.; Lam, L.; George, T.C.; Wright, G.W.; Dave, S.S.; Zhao, H.; Xu, W.; Rosenwald, A.; et al. Oncogenic CARD11 Mutations in Human Diffuse Large B Cell Lymphoma. Science 2008, 319, 1676–1679. [Google Scholar] [CrossRef]
Elvers, I.; Turner-Maier, J.; Swofford, R.; Koltookian, M.; Johnson, J.; Stewart, C.; Zhang, C.-Z.; Schumacher, S.E.; Beroukhim, R.; Rosenberg, M.; et al. Exome Sequencing of Lymphomas from Three Dog Breeds Reveals Somatic Mutation Patterns Reflecting Genetic Background. Genome Res 2015, 25, 1634–1645. [Google Scholar] [CrossRef]
Sparks, A.; Woods, J.P.; Bienzle, D.; Wood, G.A.; Coomber, B.L. Whole Genome Sequencing Analysis of High Confidence Variants of B-Cell Lymphoma in Canis Familiaris. Plos One 2020, 15, e0238183. [Google Scholar] [CrossRef]
Burnett, R.C.; Vernau, W.; Modiano, J.F.; Olver, C.S.; Moore, P.F.; Avery, A.C. Diagnosis of Canine Lymphoid Neoplasia Using Clonal Rearrangements of Antigen Receptor Genes. Vet Pathol 2003, 40, 32–41. [Google Scholar] [CrossRef]
Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.
Martin, M. Cutadapt Removes Adapter Sequences from High-Throughput Sequencing Reads. EMBnetJ. 2011, 17, 10–12. [Google Scholar] [CrossRef]
Li, H.; Durbin, R. Fast and Accurate Short Read Alignment with Burrows–Wheeler Transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed]
McLaren, W.; Pritchard, B.; Rios, D.; Chen, Y.; Flicek, P.; Cunningham, F. Deriving the Consequences of Genomic Variants with the Ensembl API and SNP Effect Predictor. Bioinformatics 2010, 26, 2069–2070. [Google Scholar] [CrossRef] [PubMed]
Koboldt, D.C. Best Practices for Variant Calling in Clinical Sequencing. Genome Med. 2020, 12, 1–13. [Google Scholar] [CrossRef]
Künstner, A.; Witte, H.M.; Riedl, J.; Bernard, V.; Stölting, S.; Merz, H.; Olschewski, V.; Peter, W.; Ketzer, J.; Busch, Y.; et al. Muta-tional Landscape of High-Grade B-Cell Lymphoma with MYC-, BCL2 and/or BCL6 Rearrangements Characterized by Whole-Exome Sequencing. Haematologica 2021, 107, 1850–1863. [Google Scholar] [CrossRef]
Lohr, J.G.; Stojanov, P.; Lawrence, M.S.; Auclair, D.; Chapuy, B.; Sougnez, C.; Cruz-Gordillo, P.; Knoechel, B.; Asmann, Y.W.; Slager, S.L.; et al. Discovery and Prioritization of Somatic Mutations in Diffuse Large B-Cell Lymphoma (DLBCL) by Whole-Exome Sequencing. Proc. Natl. Acad. Sci. 2012, 109, 3879–3884. [Google Scholar] [CrossRef]
Park, H.Y.; Lee, S.; Yoo, H.; Kim, S.; Kim, W.; Kim, J.; Ko, Y. Whole-Exome and Transcriptome Sequencing of Re-fractory Diffuse Large B-Cell Lymphoma. Oncotarget 2016, 7, 86433–86445. [Google Scholar] [CrossRef]
Braggio, E.; Wier, S.V.; Ojha, J.; McPhail, E.; Asmann, Y.W.; Egan, J.; Silva, J.A. da; Schiff, D.; Lopes, M.B.; Decker, P.A.; et al. Ge-nome-Wide Analysis Uncovers Novel Recurrent Alterations in Primary Central Nervous System Lymphomas. Clin. Cancer Res. 2015, 21, 3986–3994. [Google Scholar] [CrossRef]
Witkos, T.M.; Lowe, M. Recognition and Tethering of Transport Vesicles at the Golgi Apparatus. Curr Opin Cell Biol 2017, 47, 16–23. [Google Scholar] [CrossRef]
Bai, Y.; Cui, X.; Gao, D.; Wang, Y.; Wang, B.; Wang, W. Golgi Integral Membrane Protein 4 Manipulates Cellular Proliferation, Apoptosis, and Cell Cycle in Human Head and Neck Cancer. Bioscience Rep 2018, 38, BSR20180454. [Google Scholar] [CrossRef]
Martin, L.; Fluhrer, R.; Reiss, K.; Kremmer, E.; Saftig, P.; Haass, C. Regulated Intramembrane Proteolysis of Bri2 (Itm2b) by ADAM10 and SPPL2a/SPPL2b*. J Biol Chem 2008, 283, 1644–1652. [Google Scholar] [CrossRef] [PubMed]
Baron, B.W.; Pytel, P. Expression Pattern of the BCL6 and ITM2B Proteins in Normal Human Brains and in Alzheimer Disease. Appl Immunohisto M M 2017, 25, 489–496. [Google Scholar] [CrossRef] [PubMed]
Zhou, J.; Yao, Z.; Zheng, Z.; Yang, J.; Wang, R.; Fu, S.; Pan, X.; Liu, Z.; Wu, K. G-MDSCs-Derived Exosomal MiRNA-143-3p Promotes Proliferation via Targeting of ITM2B in Lung Cancer. Oncotargets Ther 2020, 13, 9701–9719. [Google Scholar] [CrossRef]
Stewart, J.A.; Wang, Y.; Ackerson, S.M.; Schuck, P.L. Emerging Roles of CST in Maintaining Genome Stability and Human Disease. Frontiers Biosci Landmark Ed 2018, 23, 1564–1586. [Google Scholar] [CrossRef] [PubMed]
Lyu, X.; Sang, P.B.; Chai, W. CST in Maintaining Genome Stability: Beyond Telomeres. Dna Repair 2021, 102, 103104. [Google Scholar] [CrossRef]
Wang, Y.; Chen, C.; Wang, X.; Jin, F.; Liu, Y.; Liu, H.; Li, T.; Fu, J. Lower DSC1 Expression Is Related to the Poor Differentiation and Prognosis of Head and Neck Squamous Cell Carcinoma (HNSCC). J Cancer Res Clin 2016, 142, 2461–2468. [Google Scholar] [CrossRef]
Cui, T.; Chen, Y.; Yang, L.; Mireskandari, M.; Knösel, T.; Zhang, Q.; Kohler, L.H.; Kunze, A.; Presselt, N.; Petersen, I. Diagnostic and Prognostic Impact of Desmocollins in Human Lung Cancer. J Clin Pathol 2012, 65, 1100. [Google Scholar] [CrossRef]
Knösel, T.; Chen, Y.; Hotovy, S.; Settmacher, U.; Altendorf-Hofmann, A.; Petersen, I. Loss of Desmocollin 1-3 and Homeobox Genes PITX1 and CDX2 Are Associated with Tumor Progression and Survival in Colorectal Carcinoma. Int J Colorectal Dis 2012, 27, 1391–1399. [Google Scholar] [CrossRef]
Nikolaev, S.I.; Rimoldi, D.; Iseli, C.; Valsesia, A.; Robyr, D.; Gehrig, C.; Harshman, K.; Guipponi, M.; Bukach, O.; Zoete, V.; et al. Exome Sequencing Identifies Recurrent Somatic MAP2K1 and MAP2K2 Mutations in Melanoma. Nat Genet 2012, 44, 133–139. [Google Scholar] [CrossRef]
Wolf, F.W.; Marks, R.M.; Sarma, V.; Byers, M.G.; Katz, R.W.; Shows, T.B.; Dixit, V.M. Characterization of a Novel Tumor Necrosis Factor-Alpha-Induced Endothelial Primary Response Gene. J Biol Chem 1992, 267, 1317–1326. [Google Scholar] [CrossRef]
Link, C.D.; Taft, A.; Kapulkin, V.; Duke, K.; Kim, S.; Fei, Q.; Wood, D.E.; Sahagan, B.G. Gene Expression Analysis in a Transgenic Caenorhabditis Elegans Alzheimer’s Disease Model. Neurobiol Aging 2003, 24, 397–413. [Google Scholar] [CrossRef] [PubMed]
Yang, L.; Liu, N.; Hu, X.; Zhang, W.; Wang, T.; Li, H.; Zhang, B.; Xiang, S.; Zhou, J.; Zhang, J. CK2 Phosphorylates TNFAIP1 to Affect Its Subcellular Localization and Interaction with PCNA. Mol Biol Rep 2010, 37, 2967–2973. [Google Scholar] [CrossRef] [PubMed]
Grinchuk, O.V.; Motakis, E.; Kuznetsov, V.A. Grinchuk, O.V.; Motakis, E.; Kuznetsov, V.A. Complex Sense-Antisense Architecture of TNFAIP1/POLDIP2 on 17q11.2 Represents a Novel Transcriptional Structural-Functional Gene Module Involved in Breast Cancer Progression. Bmc Genomics 2010, 11, S9–S9.
Zhou, C.; Li, X.; Zhang, X.; Liu, X.; Tan, Z.; Yang, C.; Zhang, J. MicroRNA-372 Maintains Oncogene Characteristics by Targeting TNFAIP1 and Affects NFκB Signaling in Human Gastric Carcinoma Cells. Int J Oncol 2013, 42, 635–642. [Google Scholar] [CrossRef] [PubMed]
Cui, R.; Meng, W.; Sun, H.-L.; Kim, T.; Ye, Z.; Fassan, M.; Jeon, Y.-J.; Li, B.; Vicentini, C.; Peng, Y.; et al. MicroRNA-224 Promotes Tumor Progression in Nonsmall Cell Lung Cancer. Proc National Acad Sci 2015, 112, E4288–E4297. [Google Scholar] [CrossRef] [PubMed]
Tedeschi, A.; Wutz, G.; Huet, S.; Jaritz, M.; Wuensche, A.; Schirghuber, E.; Davidson, I.F.; Tang, W.; Cisneros, D.A.; Bhaskara, V.; et al. Wapl Is an Essential Regulator of Chromatin Structure and Chromosome Segregation. Nature 2013, 501, 564–568. [Google Scholar] [CrossRef]
Liu, N.Q.; Maresca, M.; Brand, T. van den; Braccioli, L.; Schijns, M.M.G.A.; Teunissen, H.; Bruneau, B.G.; Nora, E.P.; Wit, E. de WAPL Maintains a Cohesin Loading Cycle to Preserve Cell-Type Specific Distal Gene Regulation. Nat Genet 2021, 53, 100–109. [Google Scholar] [CrossRef]
He, Y.; Li, C.; Xu, H.; Duan, Z.; Liu, Y.; Zeng, R.; Li, M.; Wang, B. AKT-dependent Hyperproliferation of Keratinocytes in Familial Hidradenitis Suppurativa with a NCSTN Mutation: A Potential Role of Defective MiR-100-5p. Brit J Dermatol 2020, 182, 500–502. [Google Scholar] [CrossRef]
Jin, Y.; Tymen, S.D.; Chen, D.; Fang, Z.J.; Zhao, Y.; Dragas, D.; Dai, Y.; Marucha, P.T.; Zhou, X. MicroRNA-99 Family Targets AKT/MTOR Signaling Pathway in Dermal Wound Healing. Plos One 2013, 8, e64434. [Google Scholar] [CrossRef]
Fang, Z.; Yi, Y.; Shi, G.; Li, S.; Chen, S.; Lin, Y.; Li, Z.; He, Z.; Li, W.; Zhong, S. Role of Brf1 Interaction with ER α, and Significance of Its Overexpression, in Human Breast Cancer. Mol Oncol 2017, 11, 1752–1767. [Google Scholar] [CrossRef]
Zhong, Q.; Xi, S.; Liang, J.; Shi, G.; Huang, Y.; Zhang, Y.; Levy, D.; Zhong, S. The Significance of Brf1 Overexpression in Human Hepatocellular Carcinoma. Oncotarget 2015, 7, 6243–6254. [Google Scholar] [CrossRef] [PubMed]
Lin, M.; Huang, C.; Ren, W.; Chen, J.; Xia, N.; Zhong, S. Mitogen- and Stress-Activated Protein Kinase 1 Mediates Alcohol-Upregulated Transcription of Brf1 and TRNA Genes to Cause Phenotypic Alteration. Oxid Med Cell Longev 2020, 2067959. [Google Scholar] [CrossRef] [PubMed]
Loveridge, C.J.; Slater, S.; Campbell, K.J.; Nam, N.A.; Knight, J.; Ahmad, I.; Hedley, A.; Lilla, S.; Repiscak, P.; Patel, R.; et al. Correction: BRF1 Accelerates Prostate Tumourigenesis and Perturbs Immune Infiltration. Oncogene 2020, 39, 2450–2450. [Google Scholar] [CrossRef]
Li, T.; Yan, Z.; Wang, W.; Zhang, R.; Gan, W.; Lv, S.; Zeng, Z.; Hou, Y.; Yang, M. SEMA6B Overexpression Predicts Poor Prognosis and Correlates With the Tumor Immunosuppressive Microenvironment in Colorectal Cancer. Frontiers Mol Biosci 2021, 8, 687319. [Google Scholar] [CrossRef] [PubMed]
Műzes, G.; Sipos, F. Relation of Immune Semaphorin/Plexin Signaling to Carcinogenesis. Eur J Cancer Prev 2014, 23, 469–476. [Google Scholar] [CrossRef] [PubMed]
Neufeld, G.; Mumblat, Y.; Smolkin, T.; Toledano, S.; Nir-Zvi, I.; Ziv, K.; Kessler, O. The Role of the Semaphorins in Cancer. Cell Adhes Migr 2016, 10, 652–674. [Google Scholar] [CrossRef]
Sui, Y.; Liu, J.; Zhang, J.; Zheng, Z.; Wang, Z.; Jia, Z.; Meng, Z. Expression and Gene Regulation Network of Adenosine Receptor A2B in Lung Adenocarcinoma: A Potential Diagnostic and Prognostic Biomarker. Frontiers Mol Biosci 2021, 8, 663011. [Google Scholar] [CrossRef]
Vecchio, E.A.; Tan, C.Y.R.; Gregory, K.J.; Christopoulos, A.; White, P.J.; May, L.T. Ligand-Independent Adenosine A2B Receptor Constitutive Activity as a Promoter of Prostate Cancer Cell Proliferation. J Pharmacol Exp Ther 2016, 357, 36–44. [Google Scholar] [CrossRef]
Kasama, H.; Sakamoto, Y.; Kasamatsu, A.; Okamoto, A.; Koyama, T.; Minakawa, Y.; Ogawara, K.; Yokoe, H.; Shiiba, M.; Tanzawa, H.; et al. Adenosine A2b Receptor Promotes Progression of Human Oral Cancer. Bmc Cancer 2015, 15, 563. [Google Scholar] [CrossRef]
Hermans, C.; Wittevrongel, C.; Thys, C.; Smethurst, P.A.; Geet, C.V.; Freson, K. A Compound Heterozygous Mutation in Glycoprotein VI in a Patient with a Bleeding Disorder. J Thromb Haemost 2009, 7, 1356–1363. [Google Scholar] [CrossRef]
Induruwa, I.; Moroi, M.; Bonna, A.; Malcor, J. -D.; Howes, J. -M.; Warburton, E.A.; Farndale, R.W.; Jung, S.M. Platelet Collagen Receptor Glycoprotein VI-dimer Recognizes Fibrinogen and Fibrin through Their D-domains, Contributing to Platelet Adhesion and Activation during Thrombus Formation. J Thromb Haemost 2018, 16, 389–404.
Mammadova-Bach, E.; Gil-Pulido, J.; Sarukhanyan, E.; Burkard, P.; Shityakov, S.; Schonhart, C.; Stegner, D.; Remer, K.; Nurden, P.; Nurden, A.T.; et al. Platelet Glycoprotein VI Promotes Metastasis through Interaction with Cancer Cell-Derived Galectin-3. Blood 2020, 135, 1146–1160. [Google Scholar] [CrossRef] [PubMed]
Wei, Y.; Zhang, H.; Feng, Q.; Wang, S.; Shao, Y.; Wu, J.; Jin, G.; Lin, W.; Peng, X.; Xu, X. A Novel Mechanism for A-to-I RNA-Edited AZIN1 in Promoting Tumor Angiogenesis in Colorectal Cancer. Cell Death Dis 2022, 13, 294. [Google Scholar] [CrossRef] [PubMed]
Chen, L.; Li, Y.; Lin, C.H.; Chan, T.H.M.; Chow, R.K.K.; Song, Y.; Liu, M.; Yuan, Y.-F.; Fu, L.; Kong, K.L.; et al. Recoding RNA Editing of AZIN1 Predisposes to Hepatocellular Carcinoma. Nat Med 2013, 19, 209–216. [Google Scholar] [CrossRef]
Chua, H.C.; Wulf, M.; Weidling, C.; Rasmussen, L.P.; Pless, S.A. The NALCN Channel Complex Is Voltage Sensitive and Directly Modulated by Extracellular Calcium. Sci Adv 2020, 6, eaaz3154. [Google Scholar] [CrossRef] [PubMed]
Stephens, R.F.; Guan, W.; Zhorov, B.S.; Spafford, J.D. Selectivity Filters and Cysteine-Rich Extracellular Loops in Voltage-Gated Sodium, Calcium, and NALCN Channels. Front Physiol 2015, 6, 153. [Google Scholar] [CrossRef]
Gupta, M.; Kamynina, E.; Morley, S.; Chung, S.; Muakkassa, N.; Wang, H.; Brathwaite, S.; Sharma, G.; Manor, D. Plekhg4 Is a Novel Dbl Family Guanine Nucleotide Exchange Factor Protein for Rho Family GTPases*. J Biol Chem 2013, 288, 14522–14530. [Google Scholar] [CrossRef]
Tobita, T.; Nomura, S.; Morita, H.; Ko, T.; Fujita, T.; Toko, H.; Uto, K.; Hagiwara, N.; Aburatani, H.; Komuro, I. Identification of MYLK3 Mutations in Familial Dilated Cardiomyopathy. Sci Rep-uk 2017, 7, 17495.
Bateman, A.; Sandford, R. The PLAT Domain: A New Piece in the PKD1 Puzzle. Curr Biol 1999, 9, R588–S2. [Google Scholar] [CrossRef]
Grillet, N.; Schwander, M.; Hildebrand, M.S.; Sczaniecka, A.; Kolatkar, A.; Velasco, J.; Webster, J.A.; Kahrizi, K.; Najmabadi, H.; Kimberling, W.J.; et al. Mutations in LOXHD1, an Evolutionarily Conserved Stereociliary Protein, Disrupt Hair Cell Function in Mice and Cause Progressive Hearing Loss in Humans. Am J Hum Genetics 2009, 85, 328–337. [Google Scholar] [CrossRef]
Berridge, M.J.; Bootman, M.D.; Roderick, H.L. Calcium Signalling: Dynamics, Homeostasis and Remodelling. Nat Rev Mol Cell Bio 2003, 4, 517–529. [Google Scholar] [CrossRef] [PubMed]
Stewart, T.A.; Yapa, K.T.D.S.; Monteith, G.R. Altered Calcium Signaling in Cancer Cells. Biochimica Et Biophysica Acta Bba - Biomembr 2015, 1848, 2502–2511.
Borsig, L. The Role of Platelet Activation in Tumor Metastasis. Expert Rev Anticanc 2008, 8, 1247–1255. [Google Scholar] [CrossRef]
Karpatkin, S.; Pearlstein, E. Role of Platelets in Tumor Cell Metastases. Ann Intern Med 1981, 95, 636. [Google Scholar] [CrossRef] [PubMed]
Leeksma, O.C.; Miranda, N.F. de; Veelken, H. Germline Mutations Predisposing to Diffuse Large B-Cell Lymphoma. Blood Cancer J 2017, 7, e541. [Google Scholar] [CrossRef] [PubMed]

Figure 1. Standard bioinformatic analysis pipeline used in this study.

Figure 2. Analysis of somatic tumor variants show specific highly shared loci. Individual somatic variants shared between at least two subjects were plotted based on the degree of sharing (number of subjects sharing said variant) and the location of each variant on the genome.

Figure 3. Gene-level analysis showing a major PPI network being subjected to somatic mutations.: (A) Genes targeted by somatic mutations in multiple subjects were plotted based on the locus of the chromosome and the degree of sharing. (B) Genes targeted by somatic mutations with high or moderate impact in multiple subjects. (C) PPI network of the named genes targeted by shared somatic variants. Nodes with no connections were removed for enhanced visibility, and the built-in hierarchical layout of Cytoscape was used.

Figure 4. Vascular smooth muscle contraction pathway confirmed in the KEGG pathway analysis. The asterisk (*) indicates the genes associated with this pathway.

Figure 5. Calcium signaling pathway confirmed in the KEGG pathway analysis. The asterisk (*) indicates the genes associated with this pathway.

Figure 6. Platelet activation pathway confirmed in the KEGG pathway analysis. The asterisk (*) indicates the genes associated with this pathway.

Figure 7. Stereotypic moderate to high impact germline mutations found in lymphoma subjects. The mutation, position and chromosome based on CanFam3.1 is written within the bar plots with the VEP predicted effect shown as icons. Those with two or more icons have multiple effects assigned.

Table 1. Clinical characteristics of dogs used in this study.

Sample No.	Breed	Age (years)	Sex	Cytological results	PARR results (monoclonal)	Lymphoma WHO stage	Survival time (days)
1	Maltese	5	SF	Diffuse, large	IgH major	stage IV, substage a	Euthanasia
2	Welsh corgi	8	IF	Diffuse, large	IgH major	stage IV, substage b	30, loss
3	Cocker spaniel	11	IF	Diffuse, intermediate	IgH major	stage IV, substage b	Loss
4	Shih tzu	10	CM	Diffuse, intermediate	IgH major	stage V, substage b	406
5	Maltese	6	SF	Diffuse, intermediate	IgH minor	stage IV, substage b	407
6	White terrier	6	IF	Diffuse, large	IgH minor	stage IV, substage a	435, loss
7	Mixed breed	7	IM	Diffuse, large	IgH major	stage IV, substage a	-
8	Welsh corgi	9	CM	Diffuse, intermediate	IgH major	stage IV, substage a	232

PARR, PCR for antigen receptor rearrangement; WHO, World Health Organization; SF, spayed female; IF, intact female; CM, castrated male; IM, intact male.

Table 2. Primers used for amplification of Cµ (positive control), TCRγ CDR3, and Ig CDR3.

Reaction No.	Product	Primer names	Primer specificity	Primer sequence
1	Cµ	Sigmf1	Cµ	TTC CCC CTC ATC ACC TGT GA
1	Cµ	Srµ3	Cµ	GGT TGT TGA TTG CAC TGA GG
2	IgH major	CB1	VH	CAG CCT GAG AGC CGA GGA CAC
2	IgH major	CB2	JH	TGA GGA GAC GGT GAC CAG GGT
3	IgH minor	CB1	VJ	CAG CCT GAG AGC CGA GGA CAC
3	IgH minor	CB3	JH	TGA GGA CAC AAA GAG TGA GG
4	TCRγ	TCRγ1	JH	ACC CTG AGA ATT GTG CCA GG
		TCRγ2	JH	GTT ACT ATA AAC CTG GTA AC
		TCRγ3	VH	TCT GGG RTG TAY TAC TGT GCT GTC TGG

Table 3. PCR reaction conditions used for this study.

Reaction No.	Product	Initial denaturation	40 cycles			Final extension
Reaction No.	Product	Initial denaturation	Denaturation	Annealing	Extension	Final extension
1	Cµ	94°C, 15 s	94°C, 15 s	57°C, 15 s	72°C, 15 s	72°C, 15 s
2	IgH major	94°C, 15 min	94°C, 15 s	63°C, 15 s	72°C, 15 s	-
3	IgH minor	94°C, 15 min	94°C, 15 s	57°C, 15 s	72°C, 15 s	72°C, 1 min
4	TCRγ	94°C, 15 min	94°C, 15 s	52°C, 15 s	72°C, 15 s	-

Table 4. Moderate- or high-impact gene mutations found in patients with CL with existing, defined variant names.

Gene name	Chromosome	Variant type(Position reference base variant impact)	Degree of sharing
GOLIM4	34	32880086_C_T_3_prime_UTR_variant	4
ITM2B	22	3396316_CTGGGGGCGGGTGGG_C_5_prime_UTR_variant	2
STN1	28	16517367_CA_C_splice_polypyrimidine_tract_variant&intron_variant	2
PLEKHG4B	34	11903730_G_C_missense_variant&splice_region_variant	2
PLEKHG4B	34	11903747_G_T_splice_polypyrimidine_tract_variant&intron_variant	2
UNC79	8	63061314_A_G_splice_polypyrimidine_tract_variant&intron_variant	2
UNC79	8	63061318_G_A_splice_polypyrimidine_tract_variant&intron_variant	2
BRF1	8	72881347_G_A_missense_variant	2
ENSCAFG 00845007156	16	15352766_C_G_missense_variant	2
SEMA6B	20	55831203_G_A_3_prime_UTR_variant	2
DSC1	7	58326567_CT_C_frameshift_variant	2
TNFAIP1	9	43411403_T_G_missense_variant	2
MYLK3	15	8510405_G_A_splice_donor_5th_base_variant&intron_variant	2
WAPL	4	34552139_TTC_T_3_prime_UTR_variant	2
ADORA2B	5	40018212_CA_C_frameshift_variant&start_lost	2
LOXHD1	7	44757742_GAA_G_frameshift_variant	2
GP6	1	103278325_A_AG_frameshift_variant	2
AZIN1	13	4527111_A_G_stop_lost	2
NCSTN	38	21997405_T_G_missense_variant	2

GOLIM4, Golgi Integral Membrane Protein 4; ITM2B, Integral membrane protein 2B; PLEKHG4B, Pleckstrin homology domain-containing family g member 4b; SEMA6B, Semaphorin 6b; DSC1, Desmocollin1; TNFAIP1, TNF alpha induced protein 1; MYLK3, Myosin light chain kinase 3; WAPL, Wings apart-like protein homolog; ADORA2B, Adenosine A2B receptor; LOXHD1, Lipoxygenase homology domins 1; GP6, Glycoprotein VI platelet; AZIN1, Antizyme inhibitor 1; NCSTNj, Nicastrin.

Table 5. KEGG pathway analysis results using genes whose somatic mutations were confirmed through WES.

pathway ID	Name of pathways	No. of genes involved	Name of genes involved
hsa04270	vascular smooth muscle contraction	2	ADORA2B, MYLK3
hsa04020	calcium signaling pathway	2	ADORA2B, MYLK3
hsa04611	platelet activation	2	GP6, MYLK3
hsa04512	ECM-receptor interaction	1	GP6
hsa04921	oxytocin signaling pathway	1	MYLK3
hsa04810	regulation of actin cytoskeleton	1	MYLK3
hsa04360	axon guidance	1	SEMA6B
hsa04371	apelin signaling pathway	1	MYLK3
hsa05034	alcoholism	1	ADORA2B
hsa04022	cGMP-PKG signaling pathway	1	MYLK3
hsa04971	gastric acid secretion	1	MYLK3
hsa05010	Alzheimer disease	1	NCSTN
hsa04510	focal adhesion	1	MYLK3
hsa04330	notch signaling pathway	1	NCSTN
hsa04080	neuroactive ligand-receptor interaction	1	ADORA2B
hsa04015	Rap1 signaling pathway	1	ADORA2B

ADORA2B, adenosine A2B receptor; MYLK3, myosin light chain kinase 3; GP6, glycoprotein VI platelet; SEMA6B, semaphorin 6B; NCSTN, nicastrin.

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

Canine Somatic Mutations from Whole Exome Sequencing of B-Cell Lymphomas in Six Canine Breeds

Abstract