Dire wolves (Aenocyon dirus) are extinct predators of Pleistocene North America. Although phenotypically similar to living wolves (Canis lupus), dire wolves have yet to be placed confidently in the canid family tree. We generated 3.4× and 12.8× paleogenomes from two well-preserved dire wolves dating to >13,000 and >72,000 years ago and estimated consensus species trees for these and 10 canid species. Our results revealed that ∼2/3 of dire wolf ancestry is most likely derived from a lineage sister to the clade comprising the gray wolf, coyote, and dhole and the remaining ∼1/3 from a lineage near the base of Canini diversity. Our results underscore the power of paleogenomes to resolve long-standing taxonomic questions and contribute to growing evidence of the role of post-speciation gene flow as an evolutionary force.

BACKGROUND:While liquid biopsy has potential to transform cancer diagnostics through minimally-invasive detection and monitoring of tumors, the impact of preanalytical factors such as the timing and anatomical location of blood draw is not well understood. METHODS:To address this gap, we leveraged pet dogs with spontaneous cancer as a model system for liquid biopsy of plasma cell-free DNA (cfDNA), as their compressed disease timeline facilitates rapid diagnostic benchmarking. We examined key cfDNA metrics, including DNA concentration in plasma, as well as tumor fraction and fragment size ratio derived from ultra low pass whole genome sequencing. RESULTS:We show that liquid biopsy metrics in dogs are consistent with reported human metrics. The tumor content of samples is slightly higher when blood is obtained from a central vein closer to the tumor. Metrics also differ between lymphoma and non-hematopoietic cancers, supporting cancer-type-specific interpretation. Disease status tracks with liquid biopsy findings over the course of treatment over both short (hours to days) and long (weeks to months) time frames, and trends of increased tumor fraction and other metrics are observed prior to clinical relapse in dogs with lymphoma and osteosarcoma. CONCLUSIONS:Together, these data support the utility of pet dogs with cancer as a relevant system for advancing liquid biopsy platforms.

Pet dogs share human-like environments while aging on a compressed timescale, making them a powerful translational model for aging research. Using genomic and phenotypic data from 7,627 dogs in the Dog Aging Project, including 976 profiled for 159 blood metabolites and clinical analytes, we generated the first GWAS catalog in dogs. Blood traits map to orthologous loci in dogs and humans, indicating deeply conserved pathways. Breed ancestry explains substantial variance in blood traits, and selection on visible characteristics such as fur type has pleiotropic metabolic effects. Leveraging mosaic ancestry in mixed-breed dogs and longitudinal mortality data, we identify blood traits elevated in short-lived breeds that predict individual mortality risk - including globulin and potassium - and protective traits enriched in long-lived breeds, such as ethanolamine. Although some aging-associated traits relate to growth hormone pathways, many do not, indicating that aging in dogs is multifactorial. These findings establish dogs as a translational system for identifying genetic determinants and biomarkers of aging relevant to extending healthy lifespans.

Timber rattlesnakes ( Crotalus horridus ) face escalating threats in the Northeastern Appalachians, including habitat fragmentation, human encroachment, and the fungal pathogen Ophidiomyces ophiodiicola . Using untargeted sequencing of DNA extracted from scale clips, we generated both host whole-genome and metagenomic data for 97 snakes from eight populations. Analysis of the snake genomes shows the populations surveyed exhibit relatively low levels of inbreeding and are genetically distinct, but that the degree of separation correlates only weakly with geographic distance. A genome-wide association analysis identified a locus associated with black-to-yellow color variation that contains an aldehyde dehydrogenase gene ( ALDH4A1 ) related to genes involved in hair color differences in humans. Metagenomic analysis showed that O. ophiodiicola read counts were generally higher in snakes exhibiting clinical signs of Snake Fungal Disease, but some visually asymptomatic snakes had high pathogen loads. Together, these findings highlight the dual utility of untargeted sequencing for population genetics and pathogen surveillance, providing a foundation for future studies of adaptation, disease dynamics, and conservation in this declining species. ### Competing Interest Statement The authors have declared no competing interest.

Genetic tests for behavioral and personality traits in dogs are now being marketed to pet owners, but their predictive accuracy has not been validated. To evaluate the reliability of such tests, we analyzed data from Darwin's Ark, a community science initiative that includes over 3,000 dogs with both genetic data and individual-level behavioral phenotypes. None of the candidate variants had significant associations or predictive power for behavioral traits as previously reported. However, we found strong associations with aesthetic traits that differentiate breeds, such as height, leg length, and ear shape. Our results suggest that earlier studies using breed-average phenotypes, rather than individually measured phenotypes, were confounded by population structure. Behavior in dogs is polygenic and complex, and thus cannot be accurately predicted using tests that consider only a few genetic variants. Furthermore, behavior in dogs is only moderately heritable, and environmental influences inherently limit the potential accuracy of genomic predictions. Developing meaningful, accurate genetic predictions for complex traits that can improve dog health and welfare will require very large cohorts of individually phenotyped dogs.

Cancers in pet dogs are prevalent, progress rapidly, and closely resemble human cancers, positioning them as powerful models for precision oncology. While genetic drivers of human cancer often transcend histologic boundaries, most comparative studies have focused on matched cancer types, leaving the broader scope of genomic similarity unresolved. We performed the first exome-wide, histology-agnostic comparison of canine and human cancers, analyzing 429 dog and 14,966 human tumors across 39 types. Mutational signatures and genes under selection are widely shared between species, and cancer types are as genomically similar between species as within species, with no greater similarity within dog breeds than between breeds. Machine-learning models identify genetic features shared by dog and human tumors of different histologies, mirroring cross-histology patterns in human cancer. These findings establish dog cancer as a powerful system for genomics-informed precision oncology and support pan-cancer approaches to discover translationally relevant models beyond histologic classification. ### Competing Interest Statement Drs. Karlsson and London are board of medical and veterinary advisors at FidoCure National Cancer Institute, https://ror.org/040gcmg81, R01CA255319 National Cancer Institute, F32CA247088 National Cancer Institute, https://ror.org/040gcmg81, R37CA218570

Obsessive-compulsive disorder (OCD) affects ~1% of children and adults and is partly caused by genetic factors. We conducted a genome-wide association study (GWAS) meta-analysis combining 53,660 OCD cases and 2,044,417 controls and identified 30 independent genome-wide significant loci. Gene-based approaches identified 249 potential effector genes for OCD, with 25 of these classified as the most likely causal candidates, including WDR6, DALRD3 and CTNND1 and multiple genes in the major histocompatibility complex (MHC) region. We estimated that ~11,500 genetic variants explained 90% of OCD genetic heritability. OCD genetic risk was associated with excitatory neurons in the hippocampus and the cortex, along with D1 and D2 type dopamine receptor-containing medium spiny neurons. OCD genetic risk was shared with 65 of 112 additional phenotypes, including all the psychiatric disorders we examined. In particular, OCD shared genetic risk with anxiety, depression, anorexia nervosa and Tourette syndrome and was negatively associated with inflammatory bowel diseases, educational attainment and body mass index.

Osteosarcoma (OS) is the most common bone tumor in both dogs and humans. It predominantly occurs in the appendicular skeleton, with about 25% of cases occurring within the axial skeleton. Progression of local disease is often the life-limiting factor for patients with axial OS, in contrast to appendicular OS, where local disease is addressed surgically, and metastatic disease remains the primary obstacle. While OS is a rare human cancer, limiting the availability of samples, its higher incidence in dogs provides a valuable comparative model for study. Both canines and humans share commonalities in clinical presentation, but dogs have an accelerated progression. Similarly, complex structural genetic changes define appendicular OS in both dogs and people, but it is unclear whether the genomic landscape of axial OS exhibits different alterations that may separate it from appendicular OS. We performed pilot whole genome sequencing of canine primary oral (maxillary or mandibular) OS tumors (n = 8) and matched normal tissue. We found that the genetics of canine oral OS largely parallel the genetics of canine appendicular OS, including an overall low number of recurrent point mutations affecting the same gene (TP53 and SETD2, 3/8 samples). Structural variants dominated the landscape of mutational changes, with recurrent variants in DMD (4/8) and DLG2 (3/8) found at a similar incidence to appendicular OS. This pilot suggests genomic similarities between oral and appendicular OS.

The process of domestication is commonly perceived as a human achievement, and domestic species are typically assumed to be those under human control. Domestic species have emerged from a greater diversity of interactions than this perspective allows, and none of the many definitions proposed for domestication can readily, reliably, and consistently distinguish domestic and nondomestic populations. Here, we propose that the process of domestication should instead be defined solely as evolution of a nonhuman population in response to an anthropogenic niche and that a domestic population is one that cannot sustain itself outside of an anthropogenic niche. As a result, this definition does not require comparisons with a presumed and largely unobservable ancestor. Instead, it focuses on the observable relationship between a nonhuman population and humans. It also avoids making assumptions about how domestication happens, thus enabling an exploration of the mechanisms underlying the process of adaptation to an anthropogenic niche. By applying this definition to plants, animals, and microbes, we illustrate its utility for investigating the evolution of the relationship between humans and other species and for anticipating which species are likely to survive in an increasingly human-influenced world. Domestication is simply an evolutionary process resulting from the interaction between two species, one of which is human. As we work to protect Earth's biodiversity, this definition allows us to understand why, in response to the conditions human societies create, some species survive and thrive, while others struggle and go extinct.

Infection with Lassa virus (LASV) can cause Lassa fever, a haemorrhagic illness with an estimated fatality rate of 29.7%, but causes no or mild symptoms in many individuals. Here, to investigate whether human genetic variation underlies the heterogeneity of LASV infection, we carried out genome-wide association studies (GWAS) as well as seroprevalence surveys, human leukocyte antigen typing and high-throughput variant functional characterization assays. We analysed Lassa fever susceptibility and fatal outcomes in 533 cases of Lassa fever and 1,986 population controls recruited over a 7 year period in Nigeria and Sierra Leone. We detected genome-wide significant variant associations with Lassa fever fatal outcomes near GRM7 and LIF in the Nigerian cohort. We also show that a haplotype bearing signatures of positive selection and overlapping LARGE1, a required LASV entry factor, is associated with decreased risk of Lassa fever in the Nigerian cohort but not in the Sierra Leone cohort. Overall, we identified variants and genes that may impact the risk of severe Lassa fever, demonstrating how GWAS can provide insight into viral pathogenesis. GWAS in difficult-to-recruit populations identifies variants associated with Lassa fever outcome and susceptibility at loci proximal to LIF, GRM7 and LARGE1.

Background The degree of gene and sequence preservation across species provides valuable insights into the relative necessity of genes from the perspective of natural selection. Here, we developed novel interspecies metrics across 462 mammalian species, GISMO (Gene identity score of mammalian orthologs) and GISMO-mis (GISMO-missense), to quantify gene loss traversing millions of years of evolution. GISMO is a measure of gene loss across mammals weighed by evolutionary distance relative to humans, whereas GISMO-mis quantifies the ratio of missense to synonymous variants across mammalian species for a given gene. Rationale Despite large sample sizes, current human constraint metrics are still not well calibrated for short genes. Traversing over 100 million years of evolution across hundreds of mammals can identify the most essential genes and improve gene-disease association. Beyond human genetics, these metrics provide measures of gene constraint to further enable mammalian genetics research. Results Our analyses showed that both metrics are strongly correlated with measures of human gene constraint for loss-of-function, missense, and copy number dosage derived from upwards of a million human samples, which highlight the power of interspecies constraint. Importantly, neither GISMO nor GISMO-mis are strongly correlated with coding sequence length. Therefore both metrics can identify novel constrained genes that were too small for existing human constraint metrics to capture. We also found that GISMO scores capture rare variant association signals across a range of phenotypes associated with decreased fecundity, such as schizophrenia, autism, and neurodevelopmental disorders. Moreover, common variant heritability of disease traits are highly enriched in the most constrained deciles of both metrics, further underscoring the biological relevance of these metrics in identifying functionally important genes. We further showed that both scores have the lowest duplication and deletion rate in the most constrained deciles for copy number variants in the UK Biobank, suggesting that it may be an important metric for dosage sensitivity. We additionally demonstrate that GISMO can improve prioritization of recessive disorder genes and captures homozygous selection. Conclusions Overall, we demonstrate that the most constrained genes for gene loss and missense variation capture the largest fraction of heritability, GISMO can help prioritize recessive disorder genes, and identify the most conserved genes across the mammalian tree. ![Figure][1] ### Competing Interest Statement The authors have declared no competing interest. [1]: pending:yes

Age in dogs is associated with the risk of many diseases, and canine size is a major factor in that risk. However, the size patterns are complex. While small size dogs tend to live longer, some diseases are more prevalent among small dogs. In this study we seek to quantify how the pattern of disease history varies across the spectrum of dog size, dog age, and their interaction. Utilizing owner-reported data on disease history from a substantial number of companion dogs enrolled in the Dog Aging Project, we investigate how body size, as measured by weight, associates with the lifetime prevalence of a reported condition and its pattern across age for various disease categories. We found significant positive associations between dog size and the lifetime prevalence of skin, bone/orthopedic, gastrointestinal, ear/nose/throat, cancer/tumor, brain/neurologic, endocrine, and infectious diseases. Similarly, dog size was negatively associated with lifetime prevalence of ocular, cardiac, liver/pancreas, and respiratory disease categories. Kidney/urinary disease prevalence did not vary by size. We also found that the association between age and lifetime disease prevalence varied by dog size for many conditions including ocular, cardiac, orthopedic, ear/nose/throat, and cancer. Controlling for sex, purebred vs. mixed-breed status, and geographic region made little difference in all disease categories we studied. Our results align with the reduced lifespan in larger dogs for most of the disease categories and suggest potential avenues for further examination.

Current genetic research on obsessive-compulsive disorder (OCD) supports contributions to risk specifically from common single nucleotide variants (SNVs), along with rare coding SNVs and small insertion-deletions (indels). The contribution to OCD risk from rare copy number variants (CNVs), however, has not been formally assessed at a similar scale. Here we describe an analysis of rare CNVs called from genotype array data in 2248 deeply phenotyped OCD cases and 3608 unaffected controls from Sweden and Norway. Cases carry an elevated burden of CNVs >= 30 kb in size (OR = 1.12, P = 1.77 x 10(-3)). The excess rate of these CNVs in cases versus controls was around 0.07 (95% CI 0.02-0.11, P = 2.58 x 10(-3)). This signal was largely driven by CNVs overlapping protein-coding regions (OR = 1.19, P = 3.08 x 10(-4)), particularly deletions impacting loss-of-function intolerant genes (pLI >0.995, OR = 4.12, P = 2.54 x 10(-5)). We did not identify any specific locus where CNV burden was associated with OCD case status at genome-wide significance, but we noted non-random recurrence of CNV deletions in cases (permutation P = 2.60 x 10(-3)). In cases where sufficient clinical data were available (n = 1612) we found that carriers of neurodevelopmental duplications were more likely to have comorbid autism (P < 0.001), and that carriers of deletions overlapping neurodevelopmental genes had lower treatment response (P = 0.02). The results demonstrate a contribution of rare CNVs to OCD risk, and suggest that studies of rare coding variation in OCD would have increased power to identify risk genes if this class of variation were incorporated into formal tests.

To date, four genome-wide association studies (GWAS) of obsessive-compulsive disorder (OCD) have been published, reporting a high single-nucleotide polymorphism (SNP)-heritability of 28% but finding only one significant SNP. A sub-stantial increase in sample size will likely lead to further identification of SNPs, genes, and biological pathways mediating the susceptibility to OCD. We conducted a GWAS meta-analysis with a 2-3-fold increase in case sample size (OCD cases: N = 37,015, controls: N = 948,616) compared to the last OCD GWAS, including six previously published cohorts (OCGAS, IOCDF-GC, IOCDF-GC-trio, NORDiC-nor, NORDiC-swe, and iPSYCH) and unpublished self-report data from 23andMe Inc. We explored the genetic architecture of OCD by conducting gene-based tests, tissue and celltype enrichment analyses, and estimating heritability and genetic correlations with 74 pheno-types. To examine a potential heterogeneity in our data, we conducted multivariable GWASs with MTAG. We found support for 15 independent genome-wide significant loci (14 new) and 79 protein-coding genes. Tissue enrichment analyses implicate multiple cortical regions, the amygdala, and hypothalamus, while cell type analyses yielded 12 cell types linked to OCD (all neurons). The SNP-based heritability of OCD was estimated to be 0.08. Using MTAG we found evidence for specific genetic underpinnings characteristic of different cohort-ascertainment and identified additional significant SNPs. OCD was genetically correlated with 40 disorders or traits-positively with all psychiatric disorders and negatively with BMI, age at first birth and multiple autoimmune diseases. The GWAS meta-analysis identified several biologically informative genes as important contributors to the aetiology of OCD. Overall, we have begun laying the groundwork through which the biology of OCD will be understood and described.

Dogs have played an outsized role in the field of behavioral genetics since its earliest days. Their unique evolutionary history and ubiquity in the modern world make them a potentially powerful model system for discovering how genetic changes lead to changes in behavior. Genomic technology has supercharged this potential by enabling scientists to sequence the DNA of thousands of dogs and test for correlations with behavioral traits. However, fractures in the early history of animal behavior between biological and psychological subfields may be impeding progress. In addition, canine behavioral genetics has included almost exclusively dogs from modern breeds, who represent just a small fraction of all dog diversity. By expanding the scope of dog behavior studies, and incorporating an evolutionary perspective on canine behavioral genetics, we can move beyond associations to understanding the complex interactions between genes and environment that lead to dog behavior.

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

Understanding the regulatory landscape of the human genome is a long-standing objective of modern biology. Using the reference-free alignment across 241 mammalian genomes produced by the Zoonomia Consortium, we charted evolutionary trajectories for 0.92 million human candidate cis-regulatory elements (cCREs) and 15.6 million human transcription factor binding sites (TFBSs). We identified 439,461 cCREs and 2,024,062 TFBSs under evolutionary constraint. Genes near constrained elements perform fundamental cellular processes, whereas genes near primate-specific elements are involved in environmental interaction, including odor perception and immune response. About 20% of TFBSs are transposable element-derived and exhibit intricate patterns of gains and losses during primate evolution whereas sequence variants associated with complex traits are enriched in constrained TFBSs. Our annotations illuminate the regulatory functions of the human genome.

We examined transposable element (TE) content of 248 placental mammal genome assemblies, the largest de novo TE curation effort in eukaryotes to date. We found that although mammals resemble one another in total TE content and diversity, they show substantial differences with regard to recent TE accumulation. This includes multiple recent expansion and quiescence events across the mammalian tree. Young TEs, particularly long interspersed elements, drive increases in genome size, whereas DNA transposons are associated with smaller genomes. Mammals tend to accumulate only a few types of TEs at any given time, with one TE type dominating. We also found association between dietary habit and the presence of DNA transposon invasions. These detailed annotations will serve as a benchmark for future comparative TE analyses among placental mammals.

Horizontal transfer of transposable elements (TEs) is an important mechanism contributing to genetic diversity and innovation. Bats (order Chiroptera) have repeatedly been shown to experience horizontal transfer of TEs at what appears to be a high rate compared with other mammals. We investigated the occurrence of horizontally transferred (HT) DNA transposons involving bats. We found over 200 putative HT elements within bats; 16 transposons were shared across distantly related mammalian clades, and 2 other elements were shared with a fish and two lizard species. Our results indicate that bats are a hotspot for horizontal transfer of DNA transposons. These events broadly coincide with the diversification of several bat clades, supporting the hypothesis that DNA transposon invasions have contributed to genetic diversification of bats.