Logo: to the web site of Uppsala University

Early mammals were nocturnal while dinosaurs dominated the daytime. Mammalian transition to daytime activity accelerated after the Cretaceous-Paleogene extinction, but the underlying mechanisms remain unclear. We identified a conserved cell-intrinsic, thermodynamic mechanism that likely facilitated this shift. In cells from diurnal mammals, protein synthesis, phosphorylation, and circadian timing were less sensitive to temperature changes than were cells from nocturnal mammals. Comparative genomics revealed accelerated evolution within essential signaling pathways, including mechanistic target of rapamycin (mTOR), that increase the robustness of diurnal cellular clocks to thermal and osmotic perturbation. In nocturnal mice, mTOR inhibition shifted cells, tissues, and behavior toward diurnal activity. These findings uncover a genetic and biochemical basis for nocturnal-diurnal switching, emphasizing how cellular signaling networks can encode complex phenotypes such as temporal niche selection.

Early mammals were nocturnal while dinosaurs dominated the daytime. Mammalian transition to daytime activity accelerated after the Cretaceous-Paleogene extinction, but the underlying mechanisms remain unclear. We identified a conserved cell-intrinsic, thermodynamic mechanism that likely facilitated this shift. In cells from diurnal mammals, protein synthesis, phosphorylation, and circadian timing were less sensitive to temperature changes than were cells from nocturnal mammals. Comparative genomics revealed accelerated evolution within essential signaling pathways, including mechanistic target of rapamycin (mTOR), that increase the robustness of diurnal cellular clocks to thermal and osmotic perturbation. In nocturnal mice, mTOR inhibition shifted cells, tissues, and behavior toward diurnal activity. These findings uncover a genetic and biochemical basis for nocturnal-diurnal switching, emphasizing how cellular signaling networks can encode complex phenotypes such as temporal niche selection.

Mammalian genomes are biased towards GC bases at third codon positions, likely due to a GC-biased ancestral genome and the selectively neutral recombination-related process of GC-biased gene conversion. The unwanted transcript hypothesis posits that this high GC content at synonymous sites may be beneficial for protecting against spurious transcripts, particularly in species with low effective population sizes. Utilising a 240 placental mammal genome alignment and single-base resolution conservation scores, we interpret sequence conservation at mammalian four-fold degenerate sites in this context and find evidence in support of the unwanted transcript hypothesis, including a strong GC bias, high conservation at sites relating to exon splicing, less human genetic variation at conserved four-fold degenerate sites, and conservation of sites important for epigenetic regulation of developmental genes. Additionally, we show that high conservation of four-fold degenerate sites in essential developmental genes, including homeobox genes, likely relates to the low mutation rates experienced by these genes.

Chromosomal inversions, by suppressing recombination, can profoundly shape genome evolution and drive adaptation. In the common quail (Coturnix coturnix), a highly mobile bird with a vast Palearctic breeding range, we previously identified a massive inversion on chromosome 1 associated with distinct phenotypes and restricted geographic distribution. Here, using a new de novo genome assembly, we characterise this inversion and uncover additional, ancient structural variation on chromosome 2 that segregates across the species' range: either two putatively linked inversions or a single, large inversion that appears as two due to scaffolding limitations. Together, the inversions encompass a remarkable 15.6% of the quail genome (153.6 Mbp), creating highly divergent haplotypes that diverged over a million years ago. While the chromosome 1 inversion is linked to phenotypic differences, including morphology and migratory behaviour, the chromosome 2 inversion(s) show no such association. Notably, all inversion regions exhibit reduced effective population size and a relaxation of purifying selection, evidenced by elevated nonsynonymous-to-synonymous substitution ratios (N/S). This suggests that inversions, particularly the geographically restricted one on chromosome 1, may act as engines of diversification, accelerating the accumulation of functional variation and potentially contributing to local adaptation, especially within isolated island populations. Our findings demonstrate how large-scale chromosomal rearrangements can compartmentalise a genome, fostering distinct evolutionary trajectories within a single, highly mobile species.

Noncoding DNA is central to our understanding of human gene regulation and complex diseases^1,2, and measuring the evolutionary sequence constraint can establish the functional relevance of putative regulatory elements in the human genome^{3,4,5,6,7,8,9}. Identifying the genomic elements that have become constrained specifically in primates has been hampered by the faster evolution of noncoding DNA compared to protein-coding DNA¹⁰, the relatively short timescales separating primate species¹¹, and the previously limited availability of whole-genome sequences¹². Here we construct a whole-genome alignment of 239 species, representing nearly half of all extant species in the primate order. Using this resource, we identified human regulatory elements that are under selective constraint across primates and other mammals at a 5% false discovery rate. We detected 111,318 DNase I hypersensitivity sites and 267,410 transcription factor binding sites that are constrained specifically in primates but not across other placental mammals and validate their cis-regulatory effects on gene expression. These regulatory elements are enriched for human genetic variants that affect gene expression and complex traits and diseases. Our results highlight the important role of recent evolution in regulatory sequence elements differentiating primates, including humans, from other placental mammals.

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (similar to 10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Background: The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 x 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.

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

A range of different genetic architectures underpin local adaptation in nature. Honey bees (Apis mellifera) in the Eastern African Mountains harbor high frequencies of two chromosomal inversions that likely govern adaptation to this high-elevation habitat. In the Americas, honey bees are hybrids of European and African ancestries and adaptation to latitudinal variation in climate correlates with the proportion of these ancestries across the genome. It is unknown which, if either, of these forms of genetic variation governs adaptation in honey bees living at high elevations in the Americas. Here, we performed whole-genome sequencing of 29 honey bees from both high- and low-elevation populations in Colombia. Analysis of genetic ancestry indicated that both populations were predominantly of African ancestry, but the East African inversions were not detected. However, individuals in the higher elevation population had significantly higher proportions of European ancestry, likely reflecting local adaptation. Several genomic regions exhibited particularly high differentiation between highland and lowland bees, containing candidate loci for local adaptation. Genes that were highly differentiated between highland and lowland populations were enriched for functions related to reproduction and sperm competition. Furthermore, variation in levels of European ancestry across the genome was correlated between populations of honey bees in the highland population and populations at higher latitudes in South America. The results are consistent with the hypothesis that adaptation to both latitude and elevation in these hybrid honey bees are mediated by variation in ancestry at many loci across the genome.

Over the last six decades, populations of the bumblebees Bombus sylvicola and Bombus balteatus in Colorado have experienced decreases in tongue length, a trait important for plant-pollinator mutualisms. It has been hypothesized that this observation reflects selection resulting from shifts in floral composition under climate change. Here we used morphometrics and population genomics to determine whether morphological change is ongoing, investigate the genetic basis of morphological variation, and analyse population structure in these populations. We generated a genome assembly of B. balteatus. We then analysed whole-genome sequencing data and morphometric measurements of 580 samples of both species from seven high-altitude localities. Out of 281 samples originally identified as B. sylvicola, 67 formed a separate genetic cluster comprising a newly-discovered cryptic species ("incognitus"). However, an absence of genetic structure within species suggests that gene flow is common between mountains. We found a significant decrease in tongue length between bees collected between 2012-2014 and in 2017, indicating that morphological shifts are ongoing. We did not discover any genetic associations with tongue length, but a SNP related to production of a proteolytic digestive enzyme was implicated in body size variation. We identified evidence of covariance between kinship and both tongue length and body size, which is suggestive of a genetic component of these traits, although it is possible that shared environmental effects between colonies are responsible. Our results provide evidence for ongoing modification of a morphological trait important for pollination and indicate that this trait probably has a complex genetic and environmental basis.