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

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these 2 fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we lay out a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., genome-wide association studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur analytically and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study, we re-examine an analysis testing for coevolution of expression levels between genes across a fungal phylogeny and show that including eigenvectors of the covariance matrix as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

Comparative Population Genomics of Arctic Sled Dogs Reveals a Deep and Complex History

Recent evidence demonstrates genomic and morphological continuity in the Arctic ancestral lineage of dogs. Here, we use the Siberian Husky to investigate the genomic legacy of the northeast Eurasian Arctic lineage and model the deep population history using genome-wide single nucleotide polymorphisms. Utilizing ancient dog-calibrated molecular clocks, we found that at least two distinct lineages of Arctic dogs existed in ancient Eurasia at the end of the Pleistocene. This pushes back the origin of sled dogs in the northeast Siberian Arctic with humans likely intentionally selecting dogs to perform different functions and keeping breeding populations that overlap in time and space relatively reproductively isolated. In modern Siberian Huskies, we found significant population structure based on how they are used by humans, recent European breed introgression in about half of the dogs that participate in races, moderate levels of inbreeding, and fewer potentially harmful variants in populations under strong selection for form and function (show, sled show, and racing populations of Siberian Huskies). As the struggle to preserve unique evolutionary lineages while maintaining genetic health intensifies across pedigreed dogs, understanding the genomic history to guide policies and best practices for breed management is crucial to sustain these ancient lineages and their unique evolutionary identity.

Parallel Evolution at the Regulatory Base-Pair Level Contributes to Mammalian Interspecific Differences in Polygenic Traits

Parallel evolution occurs when distinct lineages with similar ancestral states converge on a new phenotype. Parallel evolution has been well documented at the organ, gene pathway, and amino acid sequence level but in theory, it can also occur at individual nucleotides within noncoding regions. To examine the role of parallel evolution in shaping the biology of mammalian complex traits, we used data on single-nucleotide polymorphisms (SNPs) influencing human intraspecific variation to predict trait values in other species for 11 complex traits. We found that the alleles at SNP positions associated with human intraspecific height and red blood cell (RBC) count variation are associated with interspecific variation in the corresponding traits across mammals. These associations hold for deeper branches of mammalian evolution as well as between strains of collaborative cross mice. While variation in RBC count between primates uses both ancient and more recently evolved genomic regions, we found that only primate-specific elements were correlated with primate body size. We show that the SNP positions driving these signals are flanked by conserved sequences, maintain synteny with target genes, and overlap transcription factor binding sites. This work highlights the potential of conserved but tunable regulatory elements to be reused in parallel to facilitate evolutionary adaptation in mammals.

Dogs and their genes: what ever will they think of next?

The Edward Novitski Prize recognizes creativity and intellectual ingenuity in the solution of problems in genetics research. The prize honors scientific experimental work-either a single experimental accomplishment or a body of work. Ostrander is recognized for work developing the domestic dog as an experimental system for solving fundamental biological problems and identifying genetic sequences of relevance to human health and disease. Including work on disease and behavioral health, Ostrander has shown a dedication to creative methods for understanding canine genetics and the value of translating research organisms to human genetics.

A Nonsynonymous Substitution of Lhx3 Leads to Changes in Body Size in Dogs and Mice

Lhx3 is a LIM-homeodomain transcription factor that affects body size in mammals by regulating the secretion of pituitary hormones. Akita, Shiba Inu, and Mame Shiba Inu dogs are Japanese native dog breeds that have different body sizes. To determine whether Lhx3 plays a role in the differing body sizes of these three dog breeds, we sequenced the Lhx3 gene in the three breeds, which led to the identification of an SNP in codon 280 (S280N) associated with body size. The allele frequency at this SNP differed significantly between the large Akita and the two kinds of smaller Shiba dogs. To validate the function of this SNP on body size, we introduced this change into the Lhx3 gene of mice. Homozygous mutant mice (S279N+/+) were found to have significantly increased body lengths and weights compared to heterozygous mutant (S279N+/-) and wild-type (S279N-/-) mice several weeks after weaning. These results demonstrate that a nonsynonymous substitution in Lhx3 plays an important role in regulating body size in mammals.

Heritability within groups is uninformative about differences among groups: Cases from behavioral, evolutionary, and statistical genetics

Without the ability to control or randomize environments (or genotypes), it is difficult to determine the degree to which observed phenotypic differences between two groups of individuals are due to genetic vs. environmental differences. However, some have suggested that these concerns may be limited to pathological cases, and methods have appeared that seem to give-directly or indirectly-some support to claims that aggregate heritable variation within groups can be related to heritable variation among groups. We consider three families of approaches: the "between-group heritability" sometimes invoked in behavior genetics, the statistic [Formula: see text] used in empirical work in evolutionary quantitative genetics, and methods based on variation in ancestry in an admixed population, used in anthropological and statistical genetics. We take up these examples to show mathematically that information on within-group genetic and phenotypic information in the aggregate cannot separate among-group differences into genetic and environmental components, and we provide simulation results that support our claims. We discuss these results in terms of the long-running debate on this topic.

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these two fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we derive a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., Genome-Wide Association Studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur using analytical theory and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate this by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study of this, we re-examine an analysis testing for co-evolution of expression levels between genes across a fungal phylogeny, and show that including covariance matrix eigenvectors as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

A High-Quality Blue Whale Genome, Segmental Duplications, and Historical Demography

The blue whale, Balaenoptera musculus, is the largest animal known to have ever existed, making it an important case study in longevity and resistance to cancer. To further this and other blue whale-related research, we report a reference-quality, long-read-based genome assembly of this fascinating species. We assembled the genome from PacBio long reads and utilized Illumina/10×, optical maps, and Hi-C data for scaffolding, polishing, and manual curation. We also provided long read RNA-seq data to facilitate the annotation of the assembly by NCBI and Ensembl. Additionally, we annotated both haplotypes using TOGA and measured the genome size by flow cytometry. We then compared the blue whale genome with other cetaceans and artiodactyls, including vaquita (Phocoena sinus), the world's smallest cetacean, to investigate blue whale's unique biological traits. We found a dramatic amplification of several genes in the blue whale genome resulting from a recent burst in segmental duplications, though the possible connection between this amplification and giant body size requires further study. We also discovered sites in the insulin-like growth factor-1 gene correlated with body size in cetaceans. Finally, using our assembly to examine the heterozygosity and historical demography of Pacific and Atlantic blue whale populations, we found that the genomes of both populations are highly heterozygous and that their genetic isolation dates to the last interglacial period. Taken together, these results indicate how a high-quality, annotated blue whale genome will serve as an important resource for biology, evolution, and conservation research.

Heritability within groups is uninformative about differences among groups: cases from behavioral, evolutionary, and statistical genetics

Without the ability to control or randomize environments (or genotypes), it is difficult to determine the degree to which observed phenotypic differences between two groups of individuals are due to genetic vs. environmental differences. However, some have suggested that these concerns may be limited to pathological cases, and methods have appeared that seem to give-directly or indirectly-some support to claims that aggregate heritable variation within groups can be related to heritable variation among groups. We consider three families of approaches: the "between-group heritability" sometimes invoked in behavior genetics, the statistic P S T used in empirical work in evolutionary quantitative genetics, and methods based on variation in ancestry in an admixed population, used in anthropological and statistical genetics. We take up these examples to show mathematically that information on within-group genetic and phenotypic information in the aggregate cannot separate among-group differences into genetic and environmental components, and we provide simulation results that support our claims. We discuss these results in terms of the long-running debate on this topic.

Gene purging and the evolution of Neoave metabolism and longevity

Maintenance of the proteasome requires oxidative phosphorylation (ATP) and mitigation of oxidative damage, in an increasingly dysfunctional relationship with aging. SLC3A2 plays a role on both sides of this dichotomy as an adaptor to SLC7A5, a transporter of branched-chain amino acids (BCAA: Leu, Ile, Val), and to SLC7A11, a cystine importer supplying cysteine to the synthesis of the antioxidant glutathione. Endurance in mammalian muscle depends in part on oxidation of BCAA; however, elevated serum levels are associated with insulin resistance and shortened lifespans. Intriguingly, the evolution of modern birds (Neoaves) has entailed the purging of genes including SLC3A2, SLC7A5, -7, -8, -10, and SLC1A4, -5, largely removing BCAA exchangers and their interacting Na+/Gln symporters in pursuit of improved energetics. Additional gene purging included mitochondrial BCAA aminotransferase (BCAT2), pointing to reduced oxidation of BCAA and increased hepatic conversion to triglycerides and glucose. Fat deposits are anhydrous and highly reduced, maximizing the fuel/weight ratio for prolonged flight, but fat accumulation in muscle cells of aging humans contributes to inflammation and senescence. Duplications of the bidirectional α-ketoacid transporters SLC16A3, SLC16A7, the cystine transporters SLC7A9, SLC7A11, and N-glycan branching enzymes MGAT4B, MGAT4C in Neoaves suggests a shift to the transport of deaminated essential amino acid, and stronger mitigation of oxidative stress supported by the galectin lattice. We suggest that Alfred Lotka's theory of natural selection as a maximum power organizer (PNAS 8:151,1922) made an unusually large contribution to Neoave evolution. Further molecular analysis of Neoaves may reveal novel rewiring with applications for human health and longevity.

Patterns of enrichment and acceleration in evolutionary rates of promoters suggest a role of regulatory regions in cetacean gigantism

BACKGROUND

Cetaceans (whales, porpoises, and dolphins) are a lineage of aquatic mammals from which some species became giants. Only recently, gigantism has been investigated from the molecular point of view. Studies focused mainly on coding regions, and no data on the influence of regulatory regions on gigantism in this group was available. Accordingly, we investigated the molecular evolution of non-coding regulatory regions of genes already described in the literature for association with size in mammals, focusing mainly on the promoter regions. For this, we used Ciiider and phyloP tools. Ciiider identifies significantly enriched transcription factor binding sites, and phyloP estimates the molecular evolution rate of the promoter.

RESULTS

We found evidence of enrichment of transcription binding factors related to large body size, with distinct patterns between giant and non-giant cetaceans in the IGFBP7 and NCAPG promoters, in which repressive agents are present in small cetaceans and those that stimulate transcription, in giant cetaceans. In addition, we found evidence of acceleration in the IGF2, IGFBP2, IGFBP7, and ZFAT promoters.

CONCLUSION

Our results indicate that regulatory regions may also influence cetaceans' body size, providing candidate genes for future research to understand the molecular basis of the largest living animals.

Insulin signaling in development

Nutrient intake is obligatory for animal growth and development, but nutrients alone are not sufficient. Indeed, insulin and homologous hormones are required for normal growth even in the presence of nutrients. These hormones communicate nutrient status between organs, allowing animals to coordinate growth and metabolism with nutrient supply. Insulin and related hormones, such as insulin-like growth factors and insulin-like peptides, play important roles in development and metabolism, with defects in insulin production and signaling leading to hyperglycemia and diabetes. Here, we describe the insulin hormone family and the signal transduction pathways activated by these hormones. We highlight the roles of insulin signaling in coordinating maternal and fetal metabolism and growth during pregnancy, and we describe how secretion of insulin is regulated at different life stages. Additionally, we discuss the roles of insulin signaling in cell growth, stem cell proliferation and cell differentiation. We provide examples of the role of insulin in development across multiple model organisms: Caenorhabditis elegans, Drosophila, zebrafish, mouse and human.

The extinct Sicilian wolf shows a complex history of isolation and admixture with ancient dogs

The Sicilian wolf remained isolated in Sicily from the end of the Pleistocene until its extermination in the 1930s-1960s. Given its long-term isolation on the island and distinctive morphology, the genetic origin of the Sicilian wolf remains debated. We sequenced four nuclear genomes and five mitogenomes from the seven existing museum specimens to investigate the Sicilian wolf ancestry, relationships with extant and extinct wolves and dogs, and diversity. Our results show that the Sicilian wolf is most closely related to the Italian wolf but carries ancestry from a lineage related to European Eneolithic and Bronze Age dogs. The average nucleotide diversity of the Sicilian wolf was half of the Italian wolf, with 37-50% of its genome contained in runs of homozygosity. Overall, we show that, by the time it went extinct, the Sicilian wolf had high inbreeding and low-genetic diversity, consistent with a population in an insular environment.

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.

The Role of Socialisation in the Taming and Management of Wild Dingoes by Australian Aboriginal People

Simple Summary

The dingo (Canis dingo) is a wild-living canid endemic to mainland Australia; the descendent of an early lineage of dog introduced thousands of years ago to the continent, where it was isolated from further introductions of domestic canines until European colonisation began in 1788. Dingoes are notoriously difficult to maintain in captivity and owing to their predatory nature it is also known that they can pose a serious risk to children. Yet, written records and oral histories indicate that Aboriginal people in mainland Australia routinely practiced the rearing and keeping of dingoes in a tame state within their home communities and domestic spaces. This paper reviews historical and archaeological evidence for the management of wild and captive dingoes by Indigenous communities, revealing a substantial divide between the nature and outcomes of these interactions between historical/pre-contact Aboriginal societies and those in contemporary Australia. It is concluded that this special human-wild canid relationship has implications for the understanding of the domestication of dogs from wolves during the Late Pleistocene.

Abstract

Historical sources and Indigenous oral traditions indicate that Australian Aboriginal people commonly reared and kept the wild-caught pups of dingoes (C. dingo) as tamed companion animals. A review of the available evidence suggests Indigenous communities employed an intense socialisation process that forged close personal bonds between humans and their tame dingoes from an early age. This was complemented by oral traditions which passed down awareness of the dangers to children posed by wild or unfamiliar dingoes, and which communicated the importance of treating dingoes with respect. Together, these practices resulted in what can be interpreted as substantially altered behaviours in tamed dingoes, which, despite their naturally high prey drive, were not considered a serious threat to children and were thus able to be maintained as companion animals in the long term. This relationship is of importance for understanding the original domestication of the dog, as it demonstrates a means by which careful and deliberate socialisation by foragers could both manage risks to children’s safety posed by keeping wild canids in the domestic realm and retain them well into reproductive maturity—both issues which have been highlighted as obstacles to the domestication of dogs from wolves.

Comparative genomics uncovers the evolutionary history, demography, and molecular adaptations of South American canids

The remarkable radiation of South American (SA) canids produced 10 extant species distributed across diverse habitats, including disparate forms such as the short-legged, hypercarnivorous bush dog and the long-legged, largely frugivorous maned wolf. Despite considerable research spanning nearly two centuries, many aspects of their evolutionary history remain unknown. Here, we analyzed 31 whole genomes encompassing all extant SA canid species to assess phylogenetic relationships, interspecific hybridization, historical demography, current genetic diversity, and the molecular bases of adaptations in the bush dog and maned wolf. We found that SA canids originated from a single ancestor that colonized South America 3.9 to 3.5 Mya, followed by diversification east of the Andes and then a single colonization event and radiation of Lycalopex species west of the Andes. We detected extensive historical gene flow between recently diverged lineages and observed distinct patterns of genomic diversity and demographic history in SA canids, likely induced by past climatic cycles compounded by human-induced population declines. Genome-wide scans of selection showed that disparate limb proportions in the bush dog and maned wolf may derive from mutations in genes regulating chondrocyte proliferation and enlargement. Further, frugivory in the maned wolf may have been enabled by variants in genes associated with energy intake from short-chain fatty acids. In contrast, unique genetic variants detected in the bush dog may underlie interdigital webbing and dental adaptations for hypercarnivory. Our analyses shed light on the evolution of a unique carnivoran radiation and how it was shaped by South American topography and climate change.

The evolution of dog diet and foraging: Insights from archaeological canids in Siberia

Research on the evolution of dog foraging and diet has largely focused on scavenging during their initial domestication and genetic adaptations to starch-rich food environments following the advent of agriculture. The Siberian archaeological record evidences other critical shifts in dog foraging and diet that likely characterize Holocene dogs globally. By the Middle Holocene, body size reconstruction for Siberia dogs indicates that most were far smaller than Pleistocene wolves. This contributed to dogs' tendencies to scavenge, feed on small prey, and reduce social foraging. Stable carbon and nitrogen isotope analysis of Siberian dogs reveals that their diets were more diverse than those of Pleistocene wolves. This included habitual consumption of marine and freshwater foods by the Middle Holocene and reliance on C4 foods by the Late Holocene. Feeding on such foods and anthropogenic waste increased dogs' exposure to microbes, affected their gut microbiomes, and shaped long-term dog population history.