Idiopathic Epilepsy Risk Allele Trends in Belgian Tervuren: A Longitudinal Genetic Analysis

Idiopathic epilepsy (IE) has been known to be inherited in the Belgian Tervuren for many decades. Risk genotypes for IE in this breed have recently been identified on Canis familiaris chromosomes (CFA) 14 and 37. In the current study, the allele frequencies of these loci were analyzed to determine whether dog breeders had employed a purposeful selection against IE, leading to a reduction in risk-associated allele frequency within the breed over time. The allele frequencies of two generational groupings of Belgian Tervuren with and without IE were compared. Allele frequencies for risk-associated alleles on CFA14 were unchanged between 1985 and 2015, whereas those on CFA37 increased during that time in the control population (p < 0.05). In contrast, dogs with IE showed a decrease (p < 0.05) in the IE risk-associated allele frequency at the CFA37 locus. Seizure prevalence in the Belgian Tervuren appears to be increasing. These results suggest that, despite awareness that IE is inherited, selection against IE has not been successful.

The history of Coast Salish "woolly dogs" revealed by ancient genomics and Indigenous Knowledge

Ancestral Coast Salish societies in the Pacific Northwest kept long-haired "woolly dogs" that were bred and cared for over millennia. However, the dog wool-weaving tradition declined during the 19th century, and the population was lost. In this study, we analyzed genomic and isotopic data from a preserved woolly dog pelt from "Mutton," collected in 1859. Mutton is the only known example of an Indigenous North American dog with dominant precolonial ancestry postdating the onset of settler colonialism. We identified candidate genetic variants potentially linked with their distinct woolly phenotype. We integrated these data with interviews from Coast Salish Elders, Knowledge Keepers, and weavers about shared traditional knowledge and memories surrounding woolly dogs, their importance within Coast Salish societies, and how colonial policies led directly to their disappearance.

Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

Centuries of genome instability and evolution in soft-shell clam, Mya arenaria, bivalve transmissible neoplasia

Transmissible cancers are infectious parasitic clones that metastasize to new hosts, living past the death of the founder animal in which the cancer initiated. We investigated the evolutionary history of a cancer lineage that has spread though the soft-shell clam (Mya arenaria) population by assembling a chromosome-scale soft-shell clam reference genome and characterizing somatic mutations in transmissible cancer. We observe high mutation density, widespread copy-number gain, structural rearrangement, loss of heterozygosity, variable telomere lengths, mitochondrial genome expansion and transposable element activity, all indicative of an unstable cancer genome. We also discover a previously unreported mutational signature associated with overexpression of an error-prone polymerase and use this to estimate the lineage to be >200 years old. Our study reveals the ability for an invertebrate cancer lineage to survive for centuries while its genome continues to structurally mutate, likely contributing to the evolution of this lineage as a parasitic cancer.

Universal DNA methylation age across mammalian tissues

Aging, often considered a result of random cellular damage, can be accurately estimated using DNA methylation profiles, the foundation of pan-tissue epigenetic clocks. Here, we demonstrate the development of universal pan-mammalian clocks, using 11,754 methylation arrays from our Mammalian Methylation Consortium, which encompass 59 tissue types across 185 mammalian species. These predictive models estimate mammalian tissue age with high accuracy (r > 0.96). Age deviations correlate with human mortality risk, mouse somatotropic axis mutations and caloric restriction. We identified specific cytosines with methylation levels that change with age across numerous species. These sites, highly enriched in polycomb repressive complex 2-binding locations, are near genes implicated in mammalian development, cancer, obesity and longevity. Our findings offer new evidence suggesting that aging is evolutionarily conserved and intertwined with developmental processes across all mammals.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

DNA methylation networks underlying mammalian traits

Using DNA methylation profiles (n = 15,456) from 348 mammalian species, we constructed phyloepigenetic trees that bear marked similarities to traditional phylogenetic ones. Using unsupervised clustering across all samples, we identified 55 distinct cytosine modules, of which 30 are related to traits such as maximum life span, adult weight, age, sex, and human mortality risk. Maximum life span is associated with methylation levels in HOXL subclass homeobox genes and developmental processes and is potentially regulated by pluripotency transcription factors. The methylation state of some modules responds to perturbations such as caloric restriction, ablation of growth hormone receptors, consumption of high-fat diets, and expression of Yamanaka factors. This study reveals an intertwined evolution of the genome and epigenome that mediates the biological characteristics and traits of different mammalian species.

Genome-wide variant analyses reveal new patterns of admixture and population structure in Australian dingoes

Admixture between species is a cause for concern in wildlife management. Canids are particularly vulnerable to interspecific hybridisation, and genetic admixture has shaped their evolutionary history. Microsatellite DNA testing, relying on a small number of genetic markers and geographically restricted reference populations, has identified extensive domestic dog admixture in Australian dingoes and driven conservation management policy. But there exists a concern that geographic variation in dingo genotypes could confound ancestry analyses that use a small number of genetic markers. Here, we apply genome-wide single-nucleotide polymorphism (SNP) genotyping to a set of 402 wild and captive dingoes collected from across Australia and then carry out comparisons to domestic dogs. We then perform ancestry modelling and biogeographic analyses to characterise population structure in dingoes and investigate the extent of admixture between dingoes and dogs in different regions of the continent. We show that there are at least five distinct dingo populations across Australia. We observed limited evidence of dog admixture in wild dingoes. Our work challenges previous reports regarding the occurrence and extent of dog admixture in dingoes, as our ancestry analyses show that previous assessments severely overestimate the degree of domestic dog admixture in dingo populations, particularly in south-eastern Australia. These findings strongly support the use of genome-wide SNP genotyping as a refined method for wildlife managers and policymakers to assess and inform dingo management policy and legislation moving forwards.

The dogs of Chernobyl: Demographic insights into populations inhabiting the nuclear exclusion zone

The 1986 Chernobyl nuclear disaster initiated a series of catastrophic events resulting in long-term and widespread environmental contamination. We characterize the genetic structure of 302 dogs representing three free-roaming dog populations living within the power plant itself, as well as those 15 to 45 kilometers from the disaster site. Genome-wide profiles from Chernobyl, purebred and free-breeding dogs, worldwide reveal that the individuals from the power plant and Chernobyl City are genetically distinct, with the former displaying increased intrapopulation genetic similarity and differentiation. Analysis of shared ancestral genome segments highlights differences in the extent and timing of western breed introgression. Kinship analysis reveals 15 families, with the largest spanning all collection sites within the radioactive exclusion zone, reflecting migration of dogs between the power plant and Chernobyl City. This study presents the first characterization of a domestic species in Chernobyl, establishing their importance for genetic studies into the effects of exposure to long-term, low-dose ionizing radiation.

Abstract IA021: Studies of aging, morphology, and cancer in the domestic dog

Each of the nearly 450 domestic dog breeds in existence today has a unique history and genetic profile. Early in the domestication process, dogs were bred to perform specific tasks such as herding, guarding and hunting. Later, fanciers developed hundreds of breeds, often displaying variation in body morphology, including shape and size, in addition to coat color, skull shape, leg length, etc. We are interested in understanding the genetic underpinnings of that variation, and applying those results to studies of growth regulation, normal and abnormal behavior and disease susceptibility. Dogs are a naturally-occurring model for many complex human disorders, including cancer. One in two dogs that live to be ten or older will develop cancer in their lifetimes, with most common human cancers occur spontaneously in dogs. This offers unique opportunities to study cancer susceptibility, progression and treatment response. The high rate of single types of cancer in some breeds, such as histiocytic sarcoma in Bernese Mountain Dogs or invasive bladder cancer in Scottish Terriers offers specific advantages for genetic studies that do not exist in human population studies. As with cancer, domestic dogs present interesting opportunities for studies of aging. Surprisingly, small breeds typically live longer than larger breeds, offering the chance to study the interaction between morphology and lifespan. To take advantage of these opportunities, we have assembled the largest and most diverse dataset of dog breeds and other canines to date, with canines originating from six continents, reflecting the extensive phenotypic variation and heritage that coalesce in modern dogs, and identified over 90 million single nucleotide variants. In this presentation we explore the genetics of cancer and aging in dogs, drawing parallels to the human condition.

Citation Format: Elaine A. Ostrander. Studies of aging, morphology, and cancer in the domestic dog [abstract]. In: Proceedings of the AACR Special Conference: Aging and Cancer; 2022 Nov 17-20; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2022;83(2 Suppl_1):Abstract nr IA021.

Domestic dog lineages reveal genetic drivers of behavioral diversification

Selective breeding of domestic dogs has generated diverse breeds often optimized for performing specialized tasks. Despite the heritability of breed-typical behavioral traits, identification of causal loci has proven challenging due to the complexity of canine population structure. We overcome longstanding difficulties in identifying genetic drivers of canine behavior by developing a framework for understanding relationships between breeds and the behaviors that define them, utilizing genetic data for over 4,000 domestic, semi-feral, and wild canids and behavioral survey data for over 46,000 dogs. We identify ten major canine genetic lineages and their behavioral correlates and show that breed diversification is predominantly driven by non-coding regulatory variation. We determine that lineage-associated genes converge in neurodevelopmental co-expression networks, identifying a sheepdog-associated enrichment for interrelated axon guidance functions. This work presents a scaffold for canine diversification that positions the domestic dog as an unparalleled system for revealing the genetic origins of behavioral diversity.

Identification of a naturally-occurring canine model for early detection and intervention research in high grade urothelial carcinoma

BACKGROUND

Early detection and intervention research is expected to improve the outcomes for patients with high grade muscle invasive urothelial carcinoma (InvUC). With limited patients in suitable high-risk study cohorts, relevant animal model research is critical. Experimental animal models often fail to adequately represent human cancer. The purpose of this study was to determine the suitability of dogs with high breed-associated risk for naturally-occurring InvUC to serve as relevant models for early detection and intervention research. The feasibility of screening and early intervention, and similarities and differences between canine and human tumors, and early and later canine tumors were determined.

METHODS

STs (n=120) ≥ 6 years old with no outward evidence of urinary disease were screened at 6-month intervals for 3 years with physical exam, ultrasonography, and urinalysis with sediment exam. Cystoscopic biopsy was performed in dogs with positive screening tests. The pathological, clinical, and molecular characteristics of the "early" cancer detected by screening were determined. Transcriptomic signatures were compared between the early tumors and published findings in human InvUC, and to more advanced "later" canine tumors from STs who had the typical presentation of hematuria and urinary dysfunction. An early intervention trial of an oral cyclooxygenase inhibitor, deracoxib, was conducted in dogs with cancer detected through screening.

RESULTS

Biopsy-confirmed bladder cancer was detected in 32 (27%) of 120 STs including InvUC (n=29, three starting as dysplasia), grade 1 noninvasive cancer (n=2), and carcinoma in situ (n=1). Transcriptomic signatures including druggable targets such as EGFR and the PI3K-AKT-mTOR pathway, were very similar between canine and human InvUC, especially within luminal and basal molecular subtypes. Marked transcriptomic differences were noted between early and later canine tumors, particularly within luminal subtype tumors. The deracoxib remission rate (42% CR+PR) compared very favorably to that with single-agent cyclooxygenase inhibitors in more advanced canine InvUC (17-25%), supporting the value of early intervention.

CONCLUSIONS

The study defined a novel naturally-occurring animal model to complement experimental models for early detection and intervention research in InvUC. Research incorporating the canine model is expected to lead to improved outcomes for humans, as well as pet dogs, facing bladder cancer.

DNA methylation clocks for dogs and humans

DNA methylation profiles have been used to develop biomarkers of aging known as epigenetic clocks, which predict chronological age with remarkable accuracy and show promise for inferring health status as an indicator of biological age. Epigenetic clocks were first built to monitor human aging, but their underlying principles appear to be evolutionarily conserved, as they have now been successfully developed for many mammalian species. Here, we describe reliable and highly accurate epigenetic clocks shown to apply to 93 domestic dog breeds. The methylation profiles were generated using the mammalian methylation array, which utilizes DNA sequences that are conserved across all mammalian species. Canine epigenetic clocks were constructed to estimate age and also average time to death. We also present two highly accurate human–dog dual species epigenetic clocks (R = 0.97), which may facilitate the ready translation from canine to human use (or vice versa) of antiaging treatments being developed for longevity and preventive medicine. Finally, epigenome-wide association studies here reveal individual methylation sites that may underlie the inverse relationship between breed weight and lifespan. Overall, we describe robust biomarkers to measure aging and, potentially, health status in canines.

Genetic Origins of the Two Canis lupus familiaris (Dog) Freight Dog Populations

Despite periodic drops in popularity, Arctic sled dogs continue to play a vital role in northern societies, providing both freight transit and recreational race activities. In this study, we selected the Mackenzie River Husky, a freight dog of complex history, and the Chinook, an American Kennel Club recognized freight dog breed whose heritage reportedly overlaps that of the MKRH, for detailed population analysis. We tested each to determine their component breeds and used admixture analysis to ascertain their population structure. We utilized haplotype analysis to identify genomic regions shared between each population and their founding breeds. Our data show that the Alaskan Malamutes and modern Greenland sled dog contributed to both populations, but there are also unexpected contributions from the German Shepherd dog and Collie. We used haplotype analysis to identify genomic regions nearing fixation in population type and identify provocative genes in each region. Finally, in response to recent reports regarding the importance of dietary lipid genes in Arctic dogs, we analyzed 8 such genes in a targeted analysis observing signatures of selection in both populations at the MLXIPL gene loci. These data highlight the genetic routes that breeds of similar function have taken toward their occupation as successful sled dogs.

The effects of age, sex, weight, and breed on canid methylomes

Unlike genomes, which are static throughout the lifespan of an organism, DNA methylomes are dynamic. To study these dynamics, we developed quantitative models that measure the effect of multiple factors on DNA methylomes including, age, sex, weight, and genetics. We conducted our study in canids, which prove to be an ideal species to assess epigenetic moderators due to their extreme variability in size and well-characterized genetic structure. We collected buccal swabs from 217 canids (207 domestic dogs and 10 grey wolves) and used targeted bisulphite sequencing to measure methylomes. We also measured genotypes at over one thousand single nucleotide polymorphisms (SNPs). As expected, we found that DNA methylomes are strongly associated with age, enabling the construction of epigenetic clocks. However, we also identify novel associations between methylomes and sex, weight, and sterilization status, leading to accurate models that predict these factors. Methylomes are also affected by genetics, and we observe multiple associations between SNP loci and methylated CpGs. Finally, we show that several factors moderate the relationship between epigenetic ages and real ages, such as body weight, which increases epigenetic ageing. In conclusion, we demonstrate that the plasticity of DNA methylomes is impacted by myriad genetics and physiological factors, and that DNA methylation biomarkers are accurate predictors of age, sex and sterilization status.

Patagonian sheepdog: Genomic analyses trace the footprints of extinct UK herding dogs to South America

Most modern dog breeds were developed within the last two hundred years, following strong and recent human selection based predominantly on aesthetics, with few modern breeds constructed solely to maximize their work potential. In many cases, these working breeds represent the last remnants of now lost populations. The Patagonian sheepdog (PGOD), a rare herding breed, is a remarkable example of such a population. Maintained as an isolated population for over 130 years, the PGOD offers a unique opportunity to understand the genetic relationship amongst modern herding breeds, determine key genomic structure of the founder PGOD populations, and investigate how canine genomic data can mirror human migration patterns. We thus analyzed the population structure of 159 PGOD, comparing them with 1514 dogs representing 175 established breeds. Using 150,069 SNPs from a high-density SNP genotyping array, we establish the genomic composition, ancestry, and genetic diversity of the population, complementing genomic data with the PGOD's migratory history to South America. Our phylogenetic analysis reveals that PGODs are most closely related to modern herding breeds hailing from the United Kingdom. Admixture models illustrate a greater degree of diversity and genetic heterogeneity within the very small PGOD population than in Western European herding breeds, suggesting the PGOD predates the 200-year-old construction of most pure breeds known today. We thus propose that PGODs originated from the foundational herding dogs of the UK, prior to the Victorian explosion of breeds, and that they are the closest link to a now-extinct population of herding dogs from which modern herding breeds descended.

Darwinian genomics and diversity in the tree of life

Genomics encompasses the entire tree of life, both extinct and extant, and the evolutionary processes that shape this diversity. To date, genomic research has focused on humans, a small number of agricultural species, and established laboratory models. Fewer than 18,000 of ∼2,000,000 eukaryotic species (<1%) have a representative genome sequence in GenBank, and only a fraction of these have ancillary information on genome structure, genetic variation, gene expression, epigenetic modifications, and population diversity. This imbalance reflects a perception that human studies are paramount in disease research. Yet understanding how genomes work, and how genetic variation shapes phenotypes, requires a broad view that embraces the vast diversity of life. We have the technology to collect massive and exquisitely detailed datasets about the world, but expertise is siloed into distinct fields. A new approach, integrating comparative genomics with cell and evolutionary biology, ecology, archaeology, anthropology, and conservation biology, is essential for understanding and protecting ourselves and our world. Here, we describe potential for scientific discovery when comparative genomics works in close collaboration with a broad range of fields as well as the technical, scientific, and social constraints that must be addressed.

A Rare Germline HOXB13 Variant Contributes to Risk of Prostate Cancer in Men of African Ancestry

A rare African ancestry-specific germline deletion variant in HOXB13 (X285K, rs77179853) was recently reported in Martinican men with early-onset prostate cancer. Given the role of HOXB13 germline variation in prostate cancer, we investigated the association between HOXB13 X285K and prostate cancer risk in a large sample of 22 361 African ancestry men, including 11 688 prostate cancer cases. The risk allele was present only in men of West African ancestry, with an allele frequency in men that ranged from 0.40% in Ghana and 0.31% in Nigeria to 0% in Uganda and South Africa, with a range of frequencies in men with admixed African ancestry from North America and Europe (0-0.26%). HOXB13 X285K was associated with 2.4-fold increased odds of prostate cancer (95% confidence interval [CI] = 1.5-3.9, p = 2 × 10^-4), with greater risk observed for more aggressive and advanced disease (Gleason ≥8: odds ratio [OR] = 4.7, 95% CI = 2.3-9.5, p = 2 × 10^-5; stage T3/T4: OR = 4.5, 95% CI = 2.0-10.0, p = 2 × 10^-4; metastatic disease: OR = 5.1, 95% CI = 1.9-13.7, p = 0.001). We estimated that the allele arose in West Africa 1500-4600 yr ago. Further analysis is needed to understand how the HOXB13 X285K variant impacts the HOXB13 protein and function in the prostate. Understanding who carries this mutation may inform prostate cancer screening in men of West African ancestry. PATIENT SUMMARY: A rare African ancestry-specific germline deletion in HOXB13, found only in men of West African ancestry, was reported to be associated with an increased risk of overall and advanced prostate cancer. Understanding who carries this mutation may help inform screening for prostate cancer in men of West African ancestry.

Human-modified canids in human-modified landscapes: The evolutionary consequences of hybridization for grey wolves and free-ranging domestic dogs

Introgressive hybridization between domestic animals and their wild relatives is an indirect form of human-induced evolution, altering gene pools and phenotypic traits of wild and domestic populations. Although this process is well documented in many taxa, its evolutionary consequences are poorly understood. In this study, we assess introgression patterns in admixed populations of Eurasian wolves and free-ranging domestic dogs (FRDs), identifying chromosomal regions with significantly overrepresented hybrid ancestry and assessing whether genes located within these regions show signatures of selection. Although the dog admixture proportion in West Eurasian wolves (2.7%) was greater than the wolf admixture proportion in FRDs (0.75%), the number and average length of chromosomal blocks showing significant overrepresentation of hybrid ancestry were smaller in wolves than FRDs. In wolves, 6% of genes located within these blocks showed signatures of positive selection compared to 23% in FRDs. We found that introgression from wolves may provide a considerable adaptive advantage to FRDs, counterbalancing some of the negative effects of domestication, which can include reduced genetic diversity and excessive tameness. In wolves, introgression from FRDs is mostly driven by drift, with a small number of positively selected genes associated with brain function and behaviour. The predominance of drift may be the consequence of small effective size of wolf populations, which reduces efficiency of selection for weakly advantageous or against weakly disadvantageous introgressed variants. Small wolf population sizes result largely from human-induced habitat loss and hunting, thus linking introgression rates to anthropogenic processes. Our results imply that maintenance of large population sizes should be an important element of wolf management strategies aimed at reducing introgression rates of dog-derived variants.

Marital status and prostate cancer incidence: a pooled analysis of 12 case-control studies from the PRACTICAL consortium

While being in a committed relationship is associated with a better prostate cancer prognosis, little is known about how marital status relates to its incidence. Social support provided by marriage/relationship could promote a healthy lifestyle and an increased healthcare seeking behavior. We investigated the association between marital status and prostate cancer risk using data from the PRACTICAL Consortium. Pooled analyses were conducted combining 12 case-control studies based on histologically-confirmed incident prostate cancers and controls with information on marital status prior to diagnosis/interview. Marital status was categorized as married/partner, separated/divorced, single, or widowed. Tumours with Gleason scores ≥ 8 defined high-grade cancers, and low-grade otherwise. NCI-SEER's summary stages (local, regional, distant) indicated the extent of the cancer. Logistic regression was used to derive odds ratios (ORs) and 95% confidence intervals (CI) for the association between marital status and prostate cancer risk, adjusting for potential confounders. Overall, 14,760 cases and 12,019 controls contributed to analyses. Compared to men who were married/with a partner, widowed men had an OR of 1.19 (95% CI 1.03-1.35) of prostate cancer, with little difference between low- and high-grade tumours. Risk estimates among widowers were 1.14 (95% CI 0.97-1.34) for local, 1.53 (95% CI 1.22-1.92) for regional, and 1.56 (95% CI 1.05-2.32) for distant stage tumours. Single men had elevated risks of high-grade cancers. Our findings highlight elevated risks of incident prostate cancer among widowers, more often characterized by tumours that had spread beyond the prostate at the time of diagnosis. Social support interventions and closer medical follow-up in this sub-population are warranted.

Basal and Luminal Molecular Subtypes in Naturally-Occurring Canine Urothelial Carcinoma are Associated with Tumor Immune Signatures and Dog Breed

BACKGROUND: Improved therapies are needed for patients with invasive urothelial carcinoma (InvUC). Tailoring treatment to molecular subtypes holds promise, but requires further study, including studies in pre-clinical animal models. Naturally-occurring canine InvUC harbors luminal and basal subtypes, mimicking those observed in humans, and could offer a relevant model for the disease in people. OBJECTIVE: To further validate the canine InvUC model, clinical and tumor characteristics associated with luminal and basal subtypes in dogs were determined, with comparison to findings from humans. METHODS: RNA sequencing (RNA-seq) analyses were performed on 56 canine InvUC tissues and bladder mucosa from four normal dogs. Data were aligned to CanFam 3.1, and differentially expressed genes identified. Data were interrogated with panels of genes defining luminal and basal subtypes, immune signatures, and other tumor features. Subject and tumor characteristics, and outcome data were obtained from medical records. RESULTS: Twenty-nine tumors were classified as luminal and 27 tumors as basal subtype. Basal tumors were strongly associated with immune infiltration (OR 52.22, 95%CI 4.68–582.38, P = 0.001) and cancer progression signatures in RNA-seq analyses, more advanced clinical stage, and earlier onset of distant metastases in exploratory analyses (P = 0.0113). Luminal tumors were strongly associated with breeds at high risk for InvUC (OR 0.06, 95%CI 0.01 –0.37, P = 0.002), non-immune infiltrative signatures, and less advanced clinical stage. CONCLUSIONS: Dogs with InvUC could provide a valuable model for testing new treatment strategies in the context of molecular subtype and immune status, and the search for germline variants impacting InvUC onset and subtype.

Abstract 829: Genetic risk factors for histiocytic sarcoma in a canine disease model

Histiocytic sarcoma is a rare, aggressive hematopoietic malignancy of mature tissue histiocytes with a poorly understood etiology in humans. A histologically and clinically similar counterpart in dogs affects two breeds at unusually high frequency, with 20-25% of Flatcoated Retrievers and Bernese Mountain Dogs developing histiocytic sarcoma, while the disease is extremely rare in other breeds. In this study, we leveraged the canine model system to investigate genetic risk factors for histiocytic sarcoma in the Flatcoated Retriever. Population bottlenecks and artificial selection during breed formation within last 200 years have created the unique breed structure in dogs, characterized by long linkage disequilibrium and low within-breed genetic diversity. This allows for the identification of disease-associated traits with relatively few markers and individuals. We conducted multiple genome-wide association studies in Flatcoated Retrievers in which we identified two loci conferring risk for histiocytic sarcoma (177 cases and 132 controls, total). The first locus, on canine chromosome (CFA) 5, also colocalizes with previously identified risk for hemangiosarcoma and B-cell lymphoma, two hematopoietic malignancies in the closely related golden retriever breed. The second locus is unique to the Flatcoated Retriever and increases risk for histiocytic sarcoma in combination with CFA5. There have been multiple reports of human histiocytic sarcoma diagnosis concurrent with other hematological cancers, including lymphoma, suggested to be the result of transdifferentiation or a common neoplastic progenitor. Here we link germline predisposition for histiocytic sarcoma and B-cell lymphoma to the same locus in retrievers. Subsequent whole genome, transcriptome, and ChIP sequencing, highlighted candidate genes at each locus with roles in cell migration, as well as candidate regulatory variants with predicted impact on transcription factor binding. This comprehensive, multi-omics approach has been previously largely unexplored in the canine model. Here we advance our knowledge of the biology of histiocytic sarcoma, which will have direct impact on canine health and elucidate the genetic mechanisms underlying a rare human cancer.

Citation Format: Jacquelyn M. Evans, Heidi G. Parker, Jocelyn Plassais, Gerard R. Rutteman, Elaine A. Ostrander. Genetic risk factors for histiocytic sarcoma in a canine disease model [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 829.

Pleistocene origins, western ghost lineages, and the emerging phylogeographic history of the red wolf and coyote

The red wolf (Canis rufus) of the eastern US was driven to near-extinction by colonial-era persecution and habitat conversion, which facilitated coyote (C. latrans) range expansion and widespread hybridization with red wolves. The observation of some grey wolf (C. lupus) ancestry within red wolves sparked controversy over whether it was historically a subspecies of grey wolf with its predominant "coyote-like" ancestry obtained from post-colonial coyote hybridization (2-species hypothesis) versus a distinct species closely related to the coyote that hybridized with grey wolf (3-species hypothesis). We analysed mitogenomes sourced from before the 20th century bottleneck and coyote invasion, along with hundreds of modern amplicons, which led us to reject the 2-species model and to investigate a broader phylogeographic 3-species model suggested by the fossil record. Our findings broadly support this model, in which red wolves ranged the width of the American continent prior to arrival of the grey wolf to the mid-continent 60-80 ka; red wolves subsequently disappeared from the mid-continent, relegated to California and the eastern forests, which ushered in emergence of the coyote in their place (50-30 ka); by the early Holocene (12-10 ka), coyotes had expanded into California, where they admixed with and phenotypically replaced western red wolves in a process analogous to the 20th century coyote invasion of the eastern forests. Findings indicate that the red wolf pre-dated not only European colonization but human, and possibly coyote, presence in North America. These findings highlight the urgency of expanding conservation efforts for the red wolf.

Biology’s best friend: Bridging disciplinary gaps to advance canine science

The rich history and global abundance of domestic dogs (Canis familiaris) present a unique opportunity and an ideal model for interdisciplinary research. Canine evolutionary history demonstrates unprecedented changes across all levels of biological organization. These include diversification from highly social, pack-dwelling wild carnivores (extant gray wolves, C. lupus), to increased dependence on humans (domestication), to modern in-home colonization featuring close physical proximity to humans (interspecies bonding). The young, emerging field of “canine science” comprises diverse biological disciplines including evolution, genetics, cognition, behavior, physiology, comparative medicine, and ecology, drawing on studies of both natural and experimental systems and scaling across all levels of biological organization, from genomes to ecosystems. However, limited connections bridge the various fields associated with canine science, although in every branch it is recognized that this species is one of the most phenotypically variable mammals. Yet there has been growing interest in integrating the insights from genomic evolution with those from ecophysiology and ecology, thus facilitating a more biologically comprehensive perspective of dogs. In particular, integrative, mechanistic, and/or ecological studies have been generally underrepresented. To address these emerging interests, we have collected the most compelling questions in the field of canine biology and present avenues of current and future research. This paper serves to both orient the reader to this special issue, as well as offer a forward-looking perspective from diverse biological sub-disciplines to highlight current and future goals in canine research.

Multi-omics approach identifies germline regulatory variants associated with hematopoietic malignancies in retriever dog breeds

Histiocytic sarcoma is an aggressive hematopoietic malignancy of mature tissue histiocytes with a poorly understood etiology in humans. A histologically and clinically similar counterpart affects flat-coated retrievers (FCRs) at unusually high frequency, with 20% developing the lethal disease. The similar clinical presentation combined with the closed population structure of dogs, leading to high genetic homogeneity, makes dogs an excellent model for genetic studies of cancer susceptibility. To determine the genetic risk factors underlying histiocytic sarcoma in FCRs, we conducted multiple genome-wide association studies (GWASs), identifying two loci that confer significant risk on canine chromosomes (CFA) 5 (P_wald = 4.83x10^-9) and 19 (P_wald = 2.25x10^-7). We subsequently undertook a multi-omics approach that has been largely unexplored in the canine model to interrogate these regions, generating whole genome, transcriptome, and chromatin immunoprecipitation sequencing. These data highlight the PI3K pathway gene PIK3R6 on CFA5, and proximal candidate regulatory variants that are strongly associated with histiocytic sarcoma and predicted to impact transcription factor binding. The CFA5 association colocalizes with susceptibility loci for two hematopoietic malignancies, hemangiosarcoma and B-cell lymphoma, in the closely related golden retriever breed, revealing the risk contribution this single locus makes to multiple hematological cancers. By comparison, the CFA19 locus is unique to the FCR and harbors risk alleles associated with upregulation of TNFAIP6, which itself affects cell migration and metastasis. Together, these loci explain ~35% of disease risk, an exceptionally high value that demonstrates the advantages of domestic dogs for complex trait mapping and genetic studies of cancer susceptibility.

We have identified two regions of the canine genome that explain a striking 35% of risk for developing histiocytic sarcoma in FCRs. The disease is uniformly lethal, affects 20% of FCRs, and parallels a cancer of the same name in humans. Both regions harbor genes involved in cell migration and cancer-related pathways. The first includes variants in regulatory regions at the tumor suppressor PIK3R6 locus that are strongly associated with histiocytic sarcoma and likely confer risk for other hematopoietic cancers. FCRs with risk alleles at the second locus demonstrate increased expression of TNFAIP6, which correlates with poor prognosis in multiple human cancers. In identifying genomic differences between affected and unaffected dogs, we advance our understanding of both canine and human health biology and set the stage for the development of diagnostic and therapeutic strategies.

The effects of ionizing radiation on domestic dogs: a review of the atomic bomb testing era

Dogs were frequently employed as laboratory subjects during the era of atomic bomb testing (1950–1980), particularly in studies used to generate predictive data regarding the expected effects of accidental human occupational exposure to radiation. The bulk of these studies were only partly reported in the primary literature, despite providing vital information regarding the effects of radiation exposure on a model mammalian species. Herein we review this literature and summarize the biological effects in relation to the isotopes used and the method of radionuclide exposure. Overall, these studies demonstrate the wide range of developmental and physiological effects of exposure to radiation and radionuclides in a mid‐sized mammal.

Copy number variation underlies complex phenotypes in domestic dog breeds and other canids

Extreme phenotypic diversity, a history of artificial selection, and socioeconomic value make domestic dog breeds a compelling subject for genomic research. Copy number variation (CNV) is known to account for a significant part of inter-individual genomic diversity in other systems. However, a comprehensive genome-wide study of structural variation as it relates to breed-specific phenotypes is lacking. We have generated whole genome CNV maps for more than 300 canids. Our data set extends the canine structural variation landscape to more than 100 dog breeds, including novel variants that cannot be assessed using microarray technologies. We have taken advantage of this data set to perform the first CNV-based genome-wide association study (GWAS) in canids. We identify 96 loci that display copy number differences across breeds, which are statistically associated with a previously compiled set of breed-specific morphometrics and disease susceptibilities. Among these, we highlight the discovery of a long-range interaction involving a CNV near MED13L and TBX3, which could influence breed standard height. Integration of the CNVs with chromatin interactions, long noncoding RNA expression, and single nucleotide variation highlights a subset of specific loci and genes with potential functional relevance and the prospect to explain trait variation between dog breeds.

Rare Germline Variants in ATM Predispose to Prostate Cancer: A PRACTICAL Consortium Study

BACKGROUND

Germline ATM mutations are suggested to contribute to predisposition to prostate cancer (PrCa). Previous studies have had inadequate power to estimate variant effect sizes.

OBJECTIVE

To precisely estimate the contribution of germline ATM mutations to PrCa risk.

DESIGN, SETTING, AND PARTICIPANTS

We analysed next-generation sequencing data from 13 PRACTICAL study groups comprising 5560 cases and 3353 controls of European ancestry.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Variant Call Format files were harmonised, annotated for rare ATM variants, and classified as tier 1 (likely pathogenic) or tier 2 (potentially deleterious). Associations with overall PrCa risk and clinical subtypes were estimated.

RESULTS AND LIMITATIONS

PrCa risk was higher in carriers of a tier 1 germline ATM variant, with an overall odds ratio (OR) of 4.4 (95% confidence interval [CI]: 2.0-9.5). There was also evidence that PrCa cases with younger age at diagnosis (<65 yr) had elevated tier 1 variant frequencies (p_difference = 0.04). Tier 2 variants were also associated with PrCa risk, with an OR of 1.4 (95% CI: 1.1-1.7).

CONCLUSIONS

Carriers of pathogenic ATM variants have an elevated risk of developing PrCa and are at an increased risk for earlier-onset disease presentation. These results provide information for counselling of men and their families.

PATIENT SUMMARY

In this study, we estimated that men who inherit a likely pathogenic mutation in the ATM gene had an approximately a fourfold risk of developing prostate cancer. In addition, they are likely to develop the disease earlier.

Trans-ancestry genome-wide association meta-analysis of prostate cancer identifies new susceptibility loci and informs genetic risk prediction

Prostate cancer is a highly heritable disease with large disparities in incidence rates across ancestry populations. We conducted a multiancestry meta-analysis of prostate cancer genome-wide association studies (107,247 cases and 127,006 controls) and identified 86 new genetic risk variants independently associated with prostate cancer risk, bringing the total to 269 known risk variants. The top genetic risk score (GRS) decile was associated with odds ratios that ranged from 5.06 (95% confidence interval (CI), 4.84-5.29) for men of European ancestry to 3.74 (95% CI, 3.36-4.17) for men of African ancestry. Men of African ancestry were estimated to have a mean GRS that was 2.18-times higher (95% CI, 2.14-2.22), and men of East Asian ancestry 0.73-times lower (95% CI, 0.71-0.76), than men of European ancestry. These findings support the role of germline variation contributing to population differences in prostate cancer risk, with the GRS offering an approach for personalized risk prediction.

Genetic analysis of the modern Australian labradoodle dog breed reveals an excess of the poodle genome

The genomic diversity of the domestic dog is an invaluable resource for advancing understanding of mammalian biology, evolutionary biology, morphologic variation, and behavior. There are approximately 350 recognized breeds in the world today, many established through hybridization and selection followed by intense breeding programs aimed at retaining or enhancing specific traits. As a result, many breeds suffer from an excess of particular diseases, one of many factors leading to the recent trend of "designer breed" development, i.e. crossing purebred dogs from existing breeds in the hope that offspring will be enriched for desired traits and characteristics of the parental breeds. We used a dense panel of 150,106 SNPs to analyze the population structure of the Australian labradoodle (ALBD), to understand how such breeds are developed. Haplotype and admixture analyses show that breeds other than the poodle (POOD) and Labrador retriever (LAB) contributed to ALBD formation, but that the breed is, at the genetic level, predominantly POOD, with all small and large varieties contributing to its construction. Allele frequency analysis reveals that the breed is enhanced for variants associated with a poodle-like coat, which is perceived by breeders to have an association with hypoallergenicity. We observed little enhancement for LAB-specific alleles. This study provides a blueprint for understanding how dog breeds are formed, highlighting the limited scope of desired traits in defining new breeds.

A Germline Variant at 8q24 Contributes to Familial Clustering of Prostate Cancer in Men of African Ancestry

Although men of African ancestry have a high risk of prostate cancer (PCa), no genes or mutations have been identified that contribute to familial clustering of PCa in this population. We investigated whether the African ancestry-specific PCa risk variant at 8q24, rs72725854, is enriched in men with a PCa family history in 9052 cases, 143 cases from high-risk families, and 8595 controls of African ancestry. We found the risk allele to be significantly associated with earlier age at diagnosis, more aggressive disease, and enriched in men with a PCa family history (32% of high-risk familial cases carried the variant vs 23% of cases without a family history and 12% of controls). For cases with two or more first-degree relatives with PCa who had at least one family member diagnosed at age <60 yr, the odds ratios for TA heterozygotes and TT homozygotes were 3.92 (95% confidence interval [CI] = 2.13-7.22) and 33.41 (95% CI = 10.86-102.84), respectively. Among men with a PCa family history, the absolute risk by age 60 yr reached 21% (95% CI = 17-25%) for TA heterozygotes and 38% (95% CI = 13-65%) for TT homozygotes. We estimate that in men of African ancestry, rs72725854 accounts for 32% of the total familial risk explained by all known PCa risk variants. PATIENT SUMMARY: We found that rs72725854, an African ancestry-specific risk variant, is more common in men with a family history of prostate cancer and in those diagnosed with prostate cancer at younger ages. Men of African ancestry may benefit from the knowledge of their carrier status for this genetic risk variant to guide decisions about prostate cancer screening.

Abstract 3517: A germline variant at 8q24 contributes to familial clustering of prostate cancer in men of African ancestry

Germline variation at 8q24 is the strongest risk factor for prostate cancer (PCa) across all racial and ethnic populations. While most 8q24 associations have been observed across populations, rs72725854 [T risk allele frequency ~6%] is only found in men of African ancestry and is the strongest known genome-wide association signal for PCa in this population. We investigated whether the T allele of rs72725854 is associated with PCa family history and age at diagnosis, characteristics known to have a strong genetic component. Analyses were performed using a sample of 9,052 cases and 8,595 controls from the African Ancestry Prostate Cancer (AAPC) GWAS Consortium and the ELLIPSE/PRACTICAL OncoArray Consortium. Participants were unselected for PCa family history. Among cases, 23.7% carried at least one copy of the T allele versus 11.6% of controls. The OR was 2.29 (95% CI=2.10–2.49) for TA heterozygotes and 5.04 (95% CI=3.36–7.55) for TT homozygotes. The percentage of cases carrying the T allele was significantly greater for men with a PCa family history (27.4% vs. 22.7% without a family history, p=0.002) and for men diagnosed &lt;60 (28.2% vs. 21.6% if over ≥60, p=0.002). The mean age at diagnosis for men with the TT genotype was 61.1 years (sd=8.7), compared to 62.7 (sd=9.1) for TA heterozygotes and 64.3 (sd=8.9) for AA homozygotes (p=5.7E-14). Carrier frequency was highest among men with both a positive family history and an early diagnosis (30.8%). The T allele was also over-represented in cases with more advanced PCa, with carrier frequencies ranging from 26% for lethal PCa (metastatic disease, PSA&gt;100 ng/ml or death from PCa), 25.4% for high-risk disease (stage T3/T4, Gleason 8-10, or PSA=20-100 ng/ml), 24.6% for intermediate-risk disease (Gleason=7, stage T1/T2, and PSA=10-20 ng/ml), and 21.4% for low-risk disease (Gleason&lt;7, stage T1/T2, and PSA&lt;10 ng/ml) (p=0.026). We also examined whether the risk allele is over-represented in 144 men from PCa families, with multiple first- and/or second-degree relatives with PCa or men diagnosed with PCa ≤55 years old. Among affected probands, 32.7% carried the risk allele, with 3.5% being homozygous carriers. The OR for TA heterozygotes and TT homozygotes was 3.41 (95% CI=2.33–4.98) and 11.06 (95% CI=3.92–31.18), respectively. Among men without a family history, the absolute risk for PCa by age 60 for non-risk allele carriers was 4.3%, compared to 9.0% and 15.6% for TA heterozygotes and TT homozygotes, respectively. Absolute risks by age 60 were higher among men with a family history of prostate cancer, reaching 9.0% for non-risk allele carriers, compared to 20.8% and 37.7% for TA heterozygotes and TT homozygotes, respectively. Given the high PCa risk conveyed by rs72725854 and the greater frequency of the allele in men with more aggressive and lethal disease, carriers of the risk allele would benefit from earlier and more regular PSA screening.

Citation Format: Burcu F. Darst, Jeannette T. Bensen, Sue A. Ingles, Benjamin A. Rybicki, Barbara Nemesure, Esther M. John, Jay H. Fowke, Victoria L. Stevens, Sonja I. Berndt, Chad D. Huff, Jong Y. Park, Wei Zheng, Elaine A. Ostrander, Shiv Srivastava, John Carpten, Thomas A. Sellers, Maureen Sanderson, Dana C. Crawford, Olivier Cussenot, Jennifer Cullen, Rick A. Kittles, Jianfeng Xu, Zsofia Kote-Jarai, Luc Multigner, Marie-Elise Parent, Florence Menegaux, Geraldine Cancel-Tassin, Adam S. Kibel, Eric A. Klein, Phyllis J. Goodman, Jennifer J. Hu, Graham Casey, Anselm J. Hennis, Ian M. Thompson, Robin Leach, James L. Mohler, Elizabeth T. Fontham, Gary J. Smith, Jack A. Taylor, Rosalind A. Eeles, Laurent Brureau, Stephen J. Chanock, Stephen Watya, Janet L. Stanford, Diptasri Mandal, William B. Isaacs, Kathleen A. Cooney, William J. Blot, David V. Conti, Christopher A. Haiman. A germline variant at 8q24 contributes to familial clustering of prostate cancer in men of African ancestry [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 3517.

Two-stage Study of Familial Prostate Cancer by Whole-exome Sequencing and Custom Capture Identifies 10 Novel Genes Associated with the Risk of Prostate Cancer

BACKGROUND

Family history of prostate cancer (PCa) is a well-known risk factor, and both common and rare genetic variants are associated with the disease.

OBJECTIVE

To detect new genetic variants associated with PCa, capitalizing on the role of family history and more aggressive PCa.

DESIGN, SETTING, AND PARTICIPANTS

A two-stage design was used. In stage one, whole-exome sequencing was used to identify potential risk alleles among affected men with a strong family history of disease or with more aggressive disease (491 cases and 429 controls). Aggressive disease was based on a sum of scores for Gleason score, node status, metastasis, tumor stage, prostate-specific antigen at diagnosis, systemic recurrence, and time to PCa death. Genes identified in stage one were screened in stage two using a custom-capture design in an independent set of 2917 cases and 1899 controls.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Frequencies of genetic variants (singly or jointly in a gene) were compared between cases and controls.

RESULTS AND LIMITATIONS

Eleven genes previously reported to be associated with PCa were detected (ATM, BRCA2, HOXB13, FAM111A, EMSY, HNF1B, KLK3, MSMB, PCAT1, PRSS3, and TERT), as well as an additional 10 novel genes (PABPC1, QK1, FAM114A1, MUC6, MYCBP2, RAPGEF4, RNASEH2B, ULK4, XPO7, and THAP3). Of these 10 novel genes, all but PABPC1 and ULK4 were primarily associated with the risk of aggressive PCa.

CONCLUSIONS

Our approach demonstrates the advantage of gene sequencing in the search for genetic variants associated with PCa and the benefits of sampling patients with a strong family history of disease or an aggressive form of disease.

PATIENT SUMMARY

Multiple genes are associated with prostate cancer (PCa) among men with a strong family history of this disease or among men with an aggressive form of PCa.

Quantitative Translation of Dog-to-Human Aging by Conserved Remodeling of the DNA Methylome

All mammals progress through similar physiological stages throughout life, from early development to puberty, aging, and death. Yet, the extent to which this conserved physiology reflects underlying genomic events is unclear. Here, we map the common methylation changes experienced by mammalian genomes as they age, focusing on comparison of humans with dogs, an emerging model of aging. Using oligo-capture sequencing, we characterize methylomes of 104 Labrador retrievers spanning a 16-year age range, achieving >150× coverage within mammalian syntenic blocks. Comparison with human methylomes reveals a nonlinear relationship that translates dog-to-human years and aligns the timing of major physiological milestones between the two species, with extension to mice. Conserved changes center on developmental gene networks, which are sufficient to translate age and the effects of anti-aging interventions across multiple mammals. These results establish methylation not only as a diagnostic age readout but also as a cross-species translator of physiological aging milestones.

RNAseq expression patterns of canine invasive urothelial carcinoma reveal two distinct tumor clusters and shared regions of dysregulation with human bladder tumors

BACKGROUND

Invasive urothelial carcinoma (iUC) is highly similar between dogs and humans in terms of pathologic presentation, molecular subtypes, response to treatment and age at onset. Thus, the dog is an established and relevant model for testing and development of targeted drugs benefiting both canine and human patients. We sought to identify gene expression patterns associated with two primary types of canine iUC tumors: those that express a common somatic mutation in the BRAF gene, and those that do not.

METHODS

We performed RNAseq on tumor and normal tissues from pet dogs. Analysis of differential expression and clustering, and positional and individual expression was used to develop gene set enrichment profiles distinguishing iUC tumors with and without BRAFV595E mutations, as well as genomic regions harboring excessive numbers of dysregulated genes.

RESULTS

We identified two expression clusters that are defined by the presence/absence of a BRAFV595E (BRAFV600E in humans) somatic mutation. BRAFV595E tumors shared significantly more dysregulated genes than BRAF wild-type tumors, and vice versa, with 398 genes differentiating the two clusters. Key genes fall into clades of limited function: tissue development, cell cycle regulation, immune response, and membrane transport. The genomic site with highest number of dysregulated genes overall lies in a locus corresponding to human chromosome 8q24, a region frequently amplified in human urothelial cancers.

CONCLUSIONS

These data identify critical sets of genes that are differently regulated in association with an activating mutation in the MAPK/ERK pathway in canine iUC tumors. The experiments also highlight the value of the canine system in identifying expression patterns associated with a common, shared cancer.

DNA methylation and cis-regulation of gene expression by prostate cancer risk SNPs

Genome-wide association studies have identified more than 100 SNPs that increase the risk of prostate cancer (PrCa). We identify and compare expression quantitative trait loci (eQTLs) and CpG methylation quantitative trait loci (meQTLs) among 147 established PrCa risk SNPs in primary prostate tumors (n = 355 from a Seattle-based study and n = 495 from The Cancer Genome Atlas, TCGA) and tumor-adjacent, histologically benign samples (n = 471 from a Mayo Clinic study). The role of DNA methylation in eQTL regulation of gene expression was investigated by data triangulation using several causal inference approaches, including a proposed adaptation of the Causal Inference Test (CIT) for causal direction. Comparing eQTLs between tumors and benign samples, we show that 98 of the 147 risk SNPs were identified as eQTLs in the tumor-adjacent benign samples, and almost all 34 eQTL identified in tumor sets were also eQTLs in the benign samples. Three lines of results support the causal role of DNA methylation. First, nearly 100 of the 147 risk SNPs were identified as meQTLs in one tumor set, and almost all eQTLs in tumors were meQTLs. Second, the loss of eQTLs in tumors relative to benign samples was associated with altered DNA methylation. Third, among risk SNPs identified as both eQTLs and meQTLs, mediation analyses suggest that over two-thirds have evidence of a causal role for DNA methylation, mostly mediating genetic influence on gene expression. In summary, we provide a comprehensive catalog of eQTLs, meQTLs and putative cancer genes for known PrCa risk SNPs. We observe that a substantial portion of germline eQTL regulatory mechanisms are maintained in the tumor development, despite somatic alterations in tumor genome. Finally, our mediation analyses illuminate the likely intermediary role of CpG methylation in eQTL regulation of gene expression.

A four-gene transcript score to predict metastatic-lethal progression in men treated for localized prostate cancer: Development and validation studies

BACKGROUND

Molecular studies have tried to address the unmet need for prognostic biomarkers in prostate cancer (PCa). Some gene expression tests improve upon clinical factors for prediction of outcomes, but additional tools for accurate prediction of tumor aggressiveness are needed.

METHODS

Based on a previously published panel of 23 gene transcripts that distinguished patients with metastatic progression, we constructed a prediction model using independent training and testing datasets. Using the validated messenger RNAs and Gleason score (GS), we performed model selection in the training set to define a final locked model to classify patients who developed metastatic-lethal events from those who remained recurrence-free. In an independent testing dataset, we compared our locked model to established clinical prognostic factors and utilized Kaplan-Meier curves and receiver operating characteristic analyses to evaluate the model's performance.

RESULTS

Thirteen of 23 previously identified gene transcripts that stratified patients with aggressive PCa were validated in the training dataset. These biomarkers plus GS were used to develop a four-gene (CST2, FBLN1, TNFRSF19, and ZNF704) transcript (4GT) score that was significantly higher in patients who progressed to metastatic-lethal events compared to those without recurrence in the testing dataset (P = 5.7 × 10^-11 ). The 4GT score provided higher prediction accuracy (area under the ROC curve [AUC] = 0.76; 95% confidence interval [CI] = 0.69-0.83; partial area under the ROC curve [pAUC] = 0.008) than GS alone (AUC = 0.63; 95% CI = 0.56-0.70; pAUC = 0.002), and it improved risk stratification in subgroups defined by a combination of clinicopathological features (ie, Cancer of the Prostate Risk Assessment-Surgery).

CONCLUSION

Our validated 4GT score has prognostic value for metastatic-lethal progression in men treated for localized PCa and warrants further evaluation for its clinical utility.

Dog10K: an international sequencing effort to advance studies of canine domestication, phenotypes and health

Dogs are the most phenotypically diverse mammalian species, and they possess more known heritable disorders than any other non-human mammal. Efforts to catalog and characterize genetic variation across well-chosen populations of canines are necessary to advance our understanding of their evolutionary history and genetic architecture. To date, no organized effort has been undertaken to sequence the world's canid populations. The Dog10K Consortium (http://www.dog10kgenomes.org) is an international collaboration of researchers from across the globe who will generate 20× whole genomes from 10 000 canids in 5 years. This effort will capture the genetic diversity that underlies the phenotypic and geographical variability of modern canids worldwide. Breeds, village dogs, niche populations and extended pedigrees are currently being sequenced, and de novo assemblies of multiple canids are being constructed. This unprecedented dataset will address the genetic underpinnings of domestication, breed formation, aging, behavior and morphological variation. More generally, this effort will advance our understanding of human and canine health.

Homozygosity for Mobile Element Insertions Associated with WBSCR17 Could Predict Success in Assistance Dog Training Programs

Assistance dog training programs can see as many as 60% of their trainees dismissed. Many training programs utilize behavioral assays prior to admittance to identify likely successful candidates, yet such assays can be insconsistent. Recently, four canine retrotransposon mobile element insertions (MEIs) in or near genes WBSCR17 (Cfa6.6 and Cfa6.7), GTF2I (Cfa6.66) and POM121 (Cfa6.83) were identified in domestic dogs and gray wolves. Variations in these MEIs were significantly associated with a heightened propensity to initiate prolonged social contact or hypersociability. Using our dataset of 837 dogs, 228 of which had paired survey-based behavioral data, we discovered that one of the insertions in WBSCR17 is the most important predictor of dog sociable behaviors related to human proximity, measured by the Canine Behavioral Assessment Research Questionnaire (C-BARQ©). We found a positive correlation between insertions at Cfa6.6 and dog separation distress in the form of restlessness when about to be left alone by the owner. Lastly, assistance dogs showed significant heterozygosity deficiency at locus Cfa6.6 and higher frequency of insertions at Cfa6.6 and Cfa6.7. We suggest that training programs could utilize this genetic survey to screen for MEIs at WBSCR17 to identify dogs with sociable traits compatible with successful assistance dog performance.

Targeted Resequencing of the Coding Sequence of 38 Genes Near Breast Cancer GWAS Loci in a Large Case-Control Study

BACKGROUND

Genes regulated by breast cancer risk alleles identified through genome-wide association studies (GWAS) may harbor rare coding risk alleles.

METHODS

We sequenced the coding regions for 38 genes within 500 kb of 38 lead GWAS SNPs in 13,538 breast cancer cases and 5,518 controls.

RESULTS

Truncating variants in these genes were rare, and were not associated with breast cancer risk. Burden testing of rare missense variants highlighted 5 genes with some suggestion of an association with breast cancer, although none met the multiple testing thresholds: MKL1, FTO, NEK10, MDM4, and COX11. Six common alleles in COX11, MAP3K1 (two), and NEK10 (three) were associated at the P < 0.0001 significance level, but these likely reflect linkage disequilibrium with causal regulatory variants.

CONCLUSIONS

There was no evidence that rare coding variants in these genes confer substantial breast cancer risks. However, more modest effect sizes could not be ruled out.

IMPACT

We tested the hypothesis that rare variants in 38 genes near breast cancer GWAS loci may mediate risk. These variants do not appear to play a major role in breast cancer heritability.

Natural Selection and Origin of a Melanistic Allele in North American Gray Wolves

Pigmentation is often used to understand how natural selection affects genetic variation in wild populations since it can have a simple genetic basis, and can affect a variety of fitness-related traits (e.g., camouflage, thermoregulation, and sexual display). In gray wolves, the K locus, a β-defensin gene, causes black coat color via a dominantly inherited KB allele. The allele is derived from dog-wolf hybridization and is at high frequency in North American wolf populations. We designed a DNA capture array to probe the geographic origin, age, and number of introgression events of the KB allele in a panel of 331 wolves and 20 dogs. We found low diversity in KB, but not ancestral ky, wolf haplotypes consistent with a selective sweep of the black haplotype across North America. Further, North American wolf KB haplotypes are monophyletic, suggesting that a single adaptive introgression from dogs to wolves most likely occurred in the Northwest Territories or Yukon. We use a new analytical approach to date the origin of the KB allele in Yukon wolves to between 1,598 and 7,248 years ago, suggesting that introgression with early Native American dogs was the source. Using population genetic simulations, we show that the K locus is undergoing natural selection in four wolf populations. We find evidence for balancing selection, specifically in Yellowstone wolves, which could be a result of selection for enhanced immunity in response to distemper. With these data, we demonstrate how the spread of an adaptive variant may have occurred across a species' geographic range.

Author Correction: Association analyses of more than 140,000 men identify 63 new prostate cancer susceptibility loci

In the version of this article initially published, the name of author Manuela Gago-Dominguez was misspelled as Manuela Gago Dominguez. The error has been corrected in the HTML and PDF version of the article.

DNA methylation profiles in African American prostate cancer patients in relation to disease progression

This study examined whether differential DNA methylation is associated with clinical features of more aggressive disease at diagnosis and prostate cancer recurrence in African American men, who are more likely to die from prostate cancer than other populations. Tumor tissues from 76 African Americans diagnosed with prostate cancer who had radical prostatectomy as their primary treatment were profiled for epigenome-wide DNA methylation levels. Long-term follow-up identified 19 patients with prostate cancer recurrence. Twenty-three CpGs were differentially methylated (FDR q≤0.25, mean methylation difference≥0.10) in patients with vs. without recurrence, including CpGs in GCK, CDKL2, PRDM13, and ZFR2. Methylation differences were also observed between men with metastatic-lethal prostate cancer vs. no recurrence (five CpGs), regional vs. local pathological stage (two CpGs), and higher vs. lower tumor aggressiveness (one CpG). These results indicate that differentially methylated CpG sites identified in tumor tissues of African American men may contribute to prostate cancer aggressiveness.

Canine Cancer Genomics: Lessons for Canine and Human Health

Dogs are second only to humans in medical surveillance and preventative health care, leading to a recent perception of increased cancer incidence. Scientific priorities in veterinary oncology have thus shifted, with a demand for cancer genetic screens, better diagnostics, and more effective therapies. Most dog breeds came into existence within the last 300 years, and many are derived from small numbers of founders. Each has undergone strong artificial selection, in which dog fanciers selected for many traits, including body size, fur type, color, skull shape, and behavior, to create novel breeds. The adoption of the breed barrier rule-no dog may become a registered member of a breed unless both its dam and its sire are registered members-ensures a relatively closed genetic pool within each breed. As a result, there is strong phenotypic homogeneity within breeds but extraordinary phenotypic variation between breeds. One consequence of this is the high level of breed-associated genetic disease. We and others have taken advantage of this to identify genes for a large number of canine maladies for which mouse models do not exist, particularly with regard to cancer.