Search

Genome-wide association studies (GWAS) in long-lived human populations have led to identification of variants associated with Alzheimer's disease and cardiovascular disease, the latter being the most common cause of mortality in people worldwide. In contrast, naturally occurring cancer represents the leading cause of death in pet dogs, and specific breeds like the Golden Retriever (GR) carry up to a 65% cancer-related death rate. We hypothesized that GWAS of long-lived GRs might lead to the identification of genetic variants capable of modifying longevity within this cancer-predisposed breed. A GWAS was performed comparing GR dogs ≥ 14 years to dogs dying prior to age 12 which revealed a significant association to ERBB4, the only member of the epidermal growth factor receptor family capable of serving as both a tumor suppressor gene and an oncogene. No coding variants were identified, however, distinct haplotypes in the 5'UTR were associated with reduced lifespan in two separate populations of GR dogs. When all GR dogs were analyzed together (n = 304), the presence of haplotype 3 was associated with shorter survival (11.8 years vs. 12.8 years, p = 0.024). GRs homozygous for haplotype 3 had the shortest survival, and GRs homozygous for haplotype 1 had the longest survival (11.6 years vs. 13.5 years, p = 0.0008). Sub-analyses revealed that the difference in lifespan for GRs carrying at least 1 copy of haplotype 3 was specific to female dogs (p = 0.009), whereas survival remained significantly different in both male and female GRs homozygous for haplotype 1 or haplotype 3 (p = 0.026 and p = 0.009, respectively). Taken together, these findings implicate a potential role for ERBB4 in GR longevity and provide evidence that within-breed canine lifespan studies could serve as a mechanism to identify favorable or disease-modifying variants important to the axis of aging and cancer.

Gene retrocopies arise from the reverse transcription and insertion into the genome of processed mRNA transcripts. Although many retrocopies have acquired mutations that render them functionally inactive, most mammals retain active LINE-1 sequences capable of producing new retrocopies. New retrocopies, referred to as retro copy number variants (retroCNVs), may not be identified by standard variant calling techniques in high-throughput sequencing data. Although multiple functional FGF4 retroCNVs have been associated with skeletal dysplasias in dogs, the full landscape of canid retroCNVs has not been characterized. Here, retroCNV discovery was performed on a whole-genome sequencing data set of 293 canids from 76 breeds. We identified retroCNV parent genes via the presence of mRNA-specific 30-mers, and then identified retroCNV insertion sites through discordant read analysis. In total, we resolved insertion sites for 1911 retroCNVs from 1179 parent genes, 1236 of which appeared identical to their parent genes. Dogs had on average 54.1 total retroCNVs and 1.4 private retroCNVs. We found evidence of expression in testes for 12% (14/113) of the retroCNVs identified in six Golden Retrievers, including four chimeric transcripts, and 97 retroCNVs also had significantly elevated F ST across dog breeds, possibly indicating selection. We applied our approach to a subset of human genomes and detected an average of 4.2 retroCNVs per sample, highlighting a 13-fold relative increase of retroCNV frequency in dogs. Particularly in canids, retroCNVs are a largely unexplored source of genetic variation that can contribute to genome plasticity and that should be considered when investigating traits and diseases.

PurposeAlthough recombinant human interleukin-15 (rhIL-15) has generated much excitement as an immunotherapeutic agent for cancer, activity in human clinical trials has been modest to date, in part due to the risks of toxicity with significant dose escalation. Since pulmonary metastases are a major site of distant failure in human and dog cancers, we sought to investigate inhaled rhIL-15 in dogs with naturally occurring lung metastases from osteosarcoma (OSA) or melanoma. We hypothesized a favorable benefit/risk profile given the concentrated delivery to the lungs with decreased systemic exposure.Experimental designWe performed a phase I trial of inhaled rhIL-15 in dogs with gross pulmonary metastases using a traditional 3+3 cohort design. A starting dose of 10 µg twice daily × 14 days was used based on human, non-human primate, and murine studies. Safety, dose-limiting toxicities (DLT), and maximum tolerated dose (MTD) were the primary objectives, while response rates, progression-free and overall survival (OS), and pharmacokinetic and immune correlative analyses were secondary.ResultsFrom October 2018 to December 2020, we enrolled 21 dogs with 18 dogs reaching the 28-day response assessment to be evaluable. At dose level 5 (70 μg), we observed two DLTs, thereby establishing 50 µg twice daily × 14 days as the MTD and recommended phase 2 dose. Among 18 evaluable dogs, we observed one complete response >1 year, one partial response with resolution of multiple target lesions, and five stable disease for an overall clinical benefit rate of 39%. Plasma rhIL-15 quantitation revealed detectable and sustained rhIL-15 concentrations between 1-hour and 6 hour postnebulization. Decreased pretreatment lymphocyte counts were significantly associated with clinical benefit. Cytotoxicity assays of banked peripheral blood mononuclear cells revealed significant increases in peak cytotoxicity against canine melanoma and OSA targets that correlated with OS.ConclusionsIn this first-in-dog clinical trial of inhaled rhIL-15 in dogs with advanced metastatic disease, we observed promising clinical activity when administered as a monotherapy for only 14 days. These data have significant clinical and biological implications for both dogs and humans with refractory lung metastases and support exploration of combinatorial therapies using inhaled rhIL-15.

The authors wish to make the following corrections to this paper [...].

Natural killer (NK) cells are key effectors of the innate immune system, but major differences between human and murine NK cells have impeded translation. Outbred dogs offer an important link for studies of NK biology and immunotherapy. We analyzed gene expression of putative NK populations from healthy dogs and dogs with naturally-occurring cancers examining differential gene expression across multiple conditions, including steady-state, in vitro activation with cytokines and co-culture, and in vivo activation with inhaled IL-15 in dogs receiving IL-15 immunotherapy. We also compared dog, mouse and human CD3-NKp46+ NK cells using a novel orthologous transcriptome. Distinct transcriptional profiles between NK populations exist between conditions and in vitro versus in vivo treatments. In cross-species analysis, canine NK cells were globally more similar to human NK cells than mice. These data define canine NK cell gene expression under multiple conditions and across species, filling an important gap in translational NK studies.

Dogs exhibit a wide variety of coat color types, and many genes have been identified that control pigment production, appearance, and distribution. Some breeds, such as the Nova Scotia Duck Tolling Retriever (NSDTR), exhibit variation in pheomelanin pigment intensity that is not explained by known genetic variants. A genome-wide association study comparing light red to dark red in the NSDTR identified a significantly associated region on canine chromosome 15 (CFA 15:23 Mb-38 Mb). Coverage analysis of whole genome sequence data from eight dogs identified a 6 kb copy number variant (CNV) 152 kb upstream of KITLG. Genotyping with digital droplet PCR (ddPCR) confirmed a significant association between an increased copy number with the dark-red coat color in NSDTR (p = 6.1 × 10-7). The copy number of the CNV was also significantly associated with coat color variation in both eumelanin and pheomelanin-based Poodles (p = 1.5 × 10-8, 4.0 × 10-9) and across other breeds. Moreover, the copy number correlated with pigment intensity along the hair shaft in both pheomelanin and eumelanin coats. KITLG plays an important role in melanogenesis, and variants upstream of KITLG have been associated with coat color variation in mice as well as hair color in humans consistent with its role in the domestic dog.

Our lack of understanding of the immune microenvironment in canine osteosarcoma (cOSA) has limited the identification of potential immunotherapeutic targets. In particular, our ability to utilize readily available tissue from a dog's primary tumour to predict the type and extent of immune response in their pulmonary metastatic lesions is unknown. We, therefore, collected 21 matched pairs of primary tumours and pulmonary metastatic lesions from dogs with OSA and performed immunohistochemistry to quantify T-lymphocyte (CD3), FOXP3+ cell, B-lymphocyte (Pax-5), and CD204+ macrophage infiltration. We found that T-lymphocytes and FOXP3+ infiltrates in primary tumours positively correlated with that of metastatic lesions (ρ = 0.512, P = 0.038 and ρ = 0.698, P = 0.007, respectively), while a strong trend existed for CD204+ infiltrates (ρ = 0.404, P = 0.087). We also observed T- and B-lymphocytes, and CD204+ macrophages to be significantly higher in a dog's pulmonary metastasis compared to their primary tumour (P = 0.018, P = 0.018, P = 0.016, respectively), while FOXP3+ cells were only significantly higher in metastases when all primary tumour and metastasis lesions were compared without pairing (P = 0.036). Together, these findings suggest that the metastatic immune microenvironment may be influenced by that of the primary cOSA, and that primary tumour immune biomarkers could potentially be applied to predict immunotherapeutic responses in gross metastatic disease. We, therefore, provide a rationale for the treatment of cOSA pulmonary metastases with immunotherapeutics that enhance the anti-tumour activity of these immune cells, particularly in dogs with moderate to high immune cell infiltration in their primary tumours.

Immunotherapeutic strategies have shown promise for the treatment of canine osteosarcoma (cOSA). Very little is known about the immune microenvironment within cOSA, however, limiting our ability to identify potential immune targets and biomarkers of therapeutic response. We therefore prospectively assessed the disease-free interval (DFI) and overall survival time (ST) of 30 dogs with cOSA treated with amputation and six doses of adjuvant carboplatin. We then quantified lymphocytic (CD3+, FOXP3+) and macrophage (CD204+) infiltrates within the primary tumours of this cohort using immunohistochemistry, and evaluated their association with outcome. Overall, the median DFI and ST were 392 and 455 days, respectively. The median number of CD3+ and FOXP3+ infiltrates were 45.8 cells/mm2 (4.6-607.6 cells/mm2 ) and 8.5 mm2 (0-163.1 cells/mm2 ), respectively. The median area of CD204+ macrophages was 4.7% (1.3%-23.3%), and dogs with tumours containing greater than 4.7% CD204+ macrophages experienced a significantly longer DFI (P = 0.016). Interestingly, a significantly lower percentage of CD204+ macrophages was detected in cOSA arising from the proximal humerus compared to other appendicular bone locations (P = 0.016). Lymphocytic infiltrates did not appear to correlate with outcome in cOSA. Overall, our findings suggest that macrophages may play a role in inhibiting cOSA progression, as has been suggested in human osteosarcoma.

Domestic dog breeds exhibit remarkable morphological variations that result from centuries of artificial selection and breeding. Identifying the genetic changes that contribute to these variations could provide critical insights into the molecular basis of tissue and organismal morphogenesis. Bulldogs, French Bulldogs and Boston Terriers share many morphological and disease-predisposition traits, including brachycephalic skull morphology, widely set eyes and short stature. Unlike other brachycephalic dogs, these breeds also exhibit vertebral malformations that result in a truncated, kinked tail (screw tail). Whole genome sequencing of 100 dogs from 21 breeds identified 12.4 million bi-allelic variants that met inclusion criteria. Whole Genome Association of these variants with the breed defining phenotype of screw tail was performed using 10 cases and 84 controls and identified a frameshift mutation in the WNT pathway gene DISHEVELLED 2 (DVL2) (Chr5: 32195043_32195044del, p = 4.37 X 10-37) as the most strongly associated variant in the canine genome. This DVL2 variant was fixed in Bulldogs and French Bulldogs and had a high allele frequency (0.94) in Boston Terriers. The DVL2 variant segregated with thoracic and caudal vertebral column malformations in a recessive manner with incomplete and variable penetrance for thoracic vertebral malformations between different breeds. Importantly, analogous frameshift mutations in the human DVL1 and DVL3 genes cause Robinow syndrome, a congenital disorder characterized by similar craniofacial, limb and vertebral malformations. Analysis of the canine DVL2 variant protein showed that its ability to undergo WNT-induced phosphorylation is reduced, suggesting that altered WNT signaling may contribute to the Robinow-like syndrome in the screwtail breeds.

Curcumin has well-established anti-cancer properties in vitro, however, its therapeutic potential has been hindered by its poor bioavailability. Lipocurc is a proprietary liposome-encapsulated curcumin formulation that enables intravenous delivery and has been shown to reach its highest concentration within lung tissue. The goal of this study was to characterize the anti-cancer and anti-angiogenic activity of Lipocurc in vitro, in addition to evaluating Lipocurc infusions in dogs with naturally occurring cancer. We therefore evaluated the effect of Lipocurc, relative to free curcumin, on the viability of canine osteosarcoma, melanoma and mammary carcinoma cell lines, as well as the ability of Lipocurc to inhibit endothelial cell viability, migration and tube formation. We also undertook a pilot clinical trial consisting of four weekly 8-hour Lipocurc infusions in 10 cancer-bearing dogs. Tumour cell proliferation was inhibited by curcumin at concentrations exceeding those achievable in the lung tissue of dogs. Similarly, equivalent high concentrations of Lipocurc and curcumin also inhibited endothelial cell viability, migration and tube formation. Four out of six dogs completing planned infusions of Lipocurc experienced stable disease; however, no radiographic responses were detected.

Gliomas are the most common form of malignant primary brain tumors in humans and second most common in dogs, occurring with similar frequencies in both species. Dogs are valuable spontaneous models of human complex diseases including cancers and may provide insight into disease susceptibility and oncogenesis. Several brachycephalic breeds such as Boxer, Bulldog and Boston Terrier have an elevated risk of developing glioma, but others, including Pug and Pekingese, are not at higher risk. To identify glioma-associated genetic susceptibility factors, an across-breed genome-wide association study (GWAS) was performed on 39 dog glioma cases and 141 controls from 25 dog breeds, identifying a genome-wide significant locus on canine chromosome (CFA) 26 (p = 2.8 x 10-8). Targeted re-sequencing of the 3.4 Mb candidate region was performed, followed by genotyping of the 56 SNVs that best fit the association pattern between the re-sequenced cases and controls. We identified three candidate genes that were highly associated with glioma susceptibility: CAMKK2, P2RX7 and DENR. CAMKK2 showed reduced expression in both canine and human brain tumors, and a non-synonymous variant in P2RX7, previously demonstrated to have a 50% decrease in receptor function, was also associated with disease. Thus, one or more of these genes appear to affect glioma susceptibility.

Over the last 20-80 million years the mammalian placenta has taken on a variety of morphologies through both divergent and convergent evolution. Recently we have shown that the human placenta genome has a unique epigenetic pattern of large partially methylated domains (PMDs) and highly methylated domains (HMDs) with gene body DNA methylation positively correlating with level of gene expression. In order to determine the evolutionary conservation of DNA methylation patterns and transcriptional regulatory programs in the placenta, we performed a genome-wide methylome (MethylC-seq) analysis of human, rhesus macaque, squirrel monkey, mouse, dog, horse, and cow placentas as well as opossum extraembryonic membrane. We found that, similar to human placenta, mammalian placentas and opossum extraembryonic membrane have globally lower levels of methylation compared to somatic tissues. Higher relative gene body methylation was the conserved feature across all mammalian placentas, despite differences in PMD/HMDs and absolute methylation levels. Specifically, higher methylation over the bodies of genes involved in mitosis, vesicle-mediated transport, protein phosphorylation, and chromatin modification was observed compared with the rest of the genome. As in human placenta, higher methylation is associated with higher gene expression and is predictive of genic location across species. Analysis of DNA methylation in oocytes and preimplantation embryos shows a conserved pattern of gene body methylation similar to the placenta. Intriguingly, mouse and cow oocytes and mouse early embryos have PMD/HMDs but their placentas do not, suggesting that PMD/HMDs are a feature of early preimplantation methylation patterns that become lost during placental development in some species and following implantation of the embryo.

BMI1, a stem cell factor and member of the polycomb group of genes, has been shown to contribute to growth and chemoresistance of several human malignancies including primary osteosarcoma (OSA). Naturally occurring OSA in the dog represents a large animal model of human OSA, however the potential role of BMI1 in canine primary and metastatic OSA has not been examined. Immunohistochemical staining of canine primary and metastatic OSA tumors revealed strong nuclear expression of BMI1. An identical staining pattern was found in both primary and metastatic human OSA tissues. Canine OSA cell lines (Abrams, Moresco, and D17) expressed high levels of BMI1 compared with canine osteoblasts and knockdown or inhibition of BMI1 by siRNA or by small molecule BMI1-inhibitor PTC-209 demonstrated a role for BMI1 in canine OSA cell growth and resistance to carboplatin and doxorubicin chemotherapy. These findings suggest that inhibition of BMI1 in primary or metastatic OSA may improve response to chemotherapy and that the dog may serve as a large animal model to evaluate such therapy.