Genome-wide analyses reveals an association between invasive urothelial carcinoma in the Shetland sheepdog and NIPAL1.

Naturally occurring canine invasive urinary carcinoma (iUC) closely resembles human muscle invasive bladder cancer in terms of histopathology, metastases, response to therapy, and low survival rate. The heterogeneous nature of the disease has led to the association of large numbers of risk loci in humans, however most are of small effect. There exists a need for new and accurate animal models of invasive bladder cancer. In dogs, distinct breeds show markedly different rates of iUC, thus presenting an opportunity to identify additional risk factors and overcome the locus heterogeneity encountered in human mapping studies. In the association study presented here, inclusive of 100 Shetland sheepdogs and 58 dogs of other breeds, we identify a homozygous protein altering point mutation within the NIPAL1 gene which increases risk by eight-fold (OR = 8.42, CI = 3.12-22.71), accounting for nearly 30% of iUC risk in the Shetland sheepdog. Inclusion of six additional loci accounts for most of the disease risk in the breed and explains nearly 75% of the phenotypes in this study. When combined with sequence data from tumors, we show that variation in the MAPK signaling pathway is an overarching cause of iUC susceptibility in dogs.

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.

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.

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.

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.

Genomic Diversity and Runs of Homozygosity in Bernese Mountain Dogs.

Bernese mountain dogs are a large dog breed formed in the early 1900s in Switzerland. While originally farm dogs that were used for pulling carts, guarding, and driving cattle, today they are considered multi-purpose companion and family dogs. The breed is predisposed to several complex diseases, such as histiocytic sarcoma, degenerative myelopathy, or hip dysplasia. Using whole-genome sequencing (WGS) data, we assessed the genomic architecture of 33 unrelated dogs from four countries: France, Sweden, Switzerland, and the United States. Analysis of runs of homozygosity (ROH) identified 12,643 ROH with an average length of 2.29 Mb and an average inbreeding coefficient of 0.395. Multidimensional scaling analysis of the genetic relatedness revealed limited clustering of European versus USA dogs, suggesting exchanges of breeding stock between continents. Furthermore, only two mtDNA haplotypes were detected in the 33 studied dogs, both of which are widespread throughout multiple dog breeds. WGS-based ROH analyses revealed several fixed or nearly fixed regions harboring discreet morphological trait-associated as well as disease-associated genetic variants. Several genes involved in the regulation of immune cells were found in the ROH shared by all dogs, which is notable in the context of the breed's strong predisposition to hematopoietic cancers. High levels of inbreeding and relatedness, strongly exaggerated in the last 30 years, have likely led to the high prevalence of specific genetic disorders in this breed.

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.

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.

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.

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.

Natural and human-driven selection of a single non-coding body size variant in ancient and modern canids.

Domestic dogs (Canis lupus familiaris) are the most variable-sized mammalian species on Earth, displaying a 40-fold size difference between breeds.1 Although dogs of variable size are found in the archeological record,2-4 the most dramatic shifts in body size are the result of selection over the last two centuries, as dog breeders selected and propagated phenotypic extremes within closed breeding populations.5 Analyses of over 200 domestic breeds have identified approximately 20 body size genes regulating insulin processing, fatty acid metabolism, TGFß signaling, and skeletal formation.6-10 Of these, insulin-like growth factor 1 (IGF1) predominates, controlling approximately 15% of body size variation between breeds.8 The identification of a functional mutation associated with IGF1 has thus far proven elusive.6,10,11 Here, to identify and elucidate the role of an ancestral IGF1 allele in the propagation of modern canids, we analyzed 1,431 genome sequences from 13 species, including both ancient and modern canids, thus allowing us to define the evolutionary history of both ancestral and derived alleles at this locus. We identified a single variant in an antisense long non-coding RNA (IGF1-AS) that interacts with the IGF1 gene, creating a duplex. While the derived mutation predominates in both modern gray wolves and large domestic breeds, the ancestral allele, which predisposes to small size, was common in small-sized breeds and smaller wild canids. Our analyses demonstrate that this major regulator of canid body size nearly vanished in Pleistocene wolves, before its recent resurgence resulting from human-imposed selection for small-sized breed dogs.

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.