The tuatara genome reveals ancient features of amniote evolution.

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

Publisher Correction: The tuatara genome reveals ancient features of amniote evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Accumulation and ineffective silencing of transposable elements on an avian W Chromosome.

One of the defining features of transposable elements (TEs) is their ability to move to new locations in the host genome. To minimize the potentially deleterious effects of de novo TE insertions, hosts have evolved several mechanisms to control TE activity, including recombination-mediated removal and epigenetic silencing; however, increasing evidence suggests that silencing of TEs is often incomplete. The crow family experienced a recent radiation of LTR retrotransposons (LTRs), offering an opportunity to gain insight into the regulatory control of young, potentially still active TEs. We quantified the abundance of TE-derived transcripts across several tissues in 15 Eurasian crows (Corvus (corone) spp.) raised under common garden conditions and find evidence for ineffective TE suppression on the female-specific W Chromosome. Using RNA-seq data, we show that ∼9.5% of all transcribed TEs had considerably greater (average, 16-fold) transcript abundance in female crows and that >85% of these female-biased TEs originated on the W Chromosome. After accounting for differences in TE density among chromosomal classes, W-linked TEs were significantly more highly expressed than TEs residing on other chromosomes, consistent with ineffective silencing on the former. Together, our results suggest that the crow W Chromosome acts as a source of transcriptionally active TEs, with possible negative fitness consequences for female birds analogous to Drosophila (an X/Y system), in which overexpression of Y-linked TEs is associated with male-specific aging and fitness loss ("toxic Y").

Major population splits coincide with episodes of rapid climate change in a forest-dependent bird.

Climate change influences population demography by altering patterns of gene flow and reproductive isolation. Direct mutation rates offer the possibility for accurate dating on the within-species level but are currently only available for a handful of vertebrate species. Here, we use the first directly estimated mutation rate in birds to study the evolutionary history of pied flycatchers (Ficedula hypoleuca). Using a combination of demographic inference and species distribution modelling, we show that all major population splits in this forest-dependent system occurred during periods of increased climate instability and rapid global temperature change. We show that the divergent Spanish subspecies originated during the Eemian-Weichselian transition 115-104 thousand years ago (kya), and not during the last glacial maximum (26.5-19 kya), as previously suggested. The magnitude and rates of climate change during the glacial-interglacial transitions that preceded population splits in pied flycatchers were similar to, or exceeded, those predicted to occur in the course of the current, human-induced climate crisis. As such, our results provide a timely reminder of the strong impact that episodes of climate instability and rapid temperature changes can have on species' evolutionary trajectories, with important implications for the natural world in the Anthropocene.

The genome of a Late Pleistocene human from a Clovis burial site in western Montana.

Clovis, with its distinctive biface, blade and osseous technologies, is the oldest widespread archaeological complex defined in North America, dating from 11,100 to 10,700 (14)C years before present (bp) (13,000 to 12,600 calendar years bp). Nearly 50 years of archaeological research point to the Clovis complex as having developed south of the North American ice sheets from an ancestral technology. However, both the origins and the genetic legacy of the people who manufactured Clovis tools remain under debate. It is generally believed that these people ultimately derived from Asia and were directly related to contemporary Native Americans. An alternative, Solutrean, hypothesis posits that the Clovis predecessors emigrated from southwestern Europe during the Last Glacial Maximum. Here we report the genome sequence of a male infant (Anzick-1) recovered from the Anzick burial site in western Montana. The human bones date to 10,705 ± 35 (14)C years bp (approximately 12,707-12,556 calendar years bp) and were directly associated with Clovis tools. We sequenced the genome to an average depth of 14.4× and show that the gene flow from the Siberian Upper Palaeolithic Mal'ta population into Native American ancestors is also shared by the Anzick-1 individual and thus happened before 12,600 years bp. We also show that the Anzick-1 individual is more closely related to all indigenous American populations than to any other group. Our data are compatible with the hypothesis that Anzick-1 belonged to a population directly ancestral to many contemporary Native Americans. Finally, we find evidence of a deep divergence in Native American populations that predates the Anzick-1 individual.

Genome-wide evidence supports mitochondrial relationships and pervasive parallel phenotypic evolution in open-habitat chats.

In wheatears and related species ('open-habitat chats'), molecular phylogenetics has led to a comprehensively revised understanding of species relationships and species diversity. Phylogenetic analyses have suggested that, in many cases, phenotypic similarities do not reflect species' relationships, revealing traditionally defined genera as non-monophyletic. This led to the suggestion of pervasive parallel evolution of open-habitat chats' plumage coloration and ecological phenotypes. However, to date, the molecular evidence for the phylogenetic relationships among open-habitat chats is mainly limited to mitochondrial DNA. Here, we assessed whether the mitochondrial relationships are supported by genome-wide data. To this end, we reconstructed the species tree among 14 open-habitat chat taxa using multi-species coalescent analyses based on ~1'300 SNPs. Our results confirm previous ones based chiefly on mitochondrial DNA; notably the paraphyly of the Oenanthe lugens complex and the clustering of individual species formerly placed in the genera Cercomela and Myrmecocichla within Oenanthe. Since several variable morphological and ecological characteristics occur in multiple places across the open-habitat chat phylogeny, our study consolidates the evidence for pervasive parallel evolution in the plumage coloration and ecology of open-habitat chats.

Parallel plumage colour evolution and introgressive hybridization in wheatears.

Genetic and phenotypic mosaics, in which various phenotypes and different genomic regions show discordant patterns of species or population divergence, offer unique opportunities to study the role of ancestral and introgressed genetic variation in phenotypic evolution. Here, we investigated the evolution of discordant phenotypic and genetic divergence in a monophyletic clade of four songbird taxa-pied wheatear (O. pleschanka), Cyprus wheatear (Oenanthe cypriaca), and western and eastern subspecies of black-eared wheatear (O. h. hispanica and O. h. melanoleuca). Phenotypically, black back and neck sides distinguish pied and Cyprus wheatears from the white-backed/necked black-eared wheatears. Meanwhile, mitochondrial variation only distinguishes western black-eared wheatear. In the absence of nuclear genetic data, and given frequent hybridization among eastern black-eared and pied wheatear, it remains unclear whether introgression is responsible for discordance between mitochondrial divergence patterns and phenotypic similarities, or whether plumage coloration evolved in parallel. Multispecies coalescent analyses of about 20,000 SNPs obtained from RAD data mapped to a draft genome assembly resolve the species tree, provide evidence for the parallel evolution of colour phenotypes and establish western and eastern black-eared wheatears as independent taxa that should be recognized as full species. The presence of the entire admixture spectrum in the Iranian hybrid zone and the detection of footprints of introgression from pied into eastern black-eared wheatear beyond the hybrid zone despite strong geographic structure of ancestry proportions furthermore suggest a potential role for introgression in parallel plumage colour evolution. Our results support the importance of standing heterospecific and/or ancestral variation in phenotypic evolution.

Improved calibration of the human mitochondrial clock using ancient genomes.

Reliable estimates of the rate at which DNA accumulates mutations (the substitution rate) are crucial for our understanding of the evolution and past demography of virtually any species. In humans, there are considerable uncertainties around these rates, with substantial variation among recent published estimates. Substitution rates have traditionally been estimated by associating dated events to the root (e.g., the divergence between humans and chimpanzees) or to internal nodes in a phylogenetic tree (e.g., first entry into the Americas). The recent availability of ancient mitochondrial DNA sequences allows for a more direct calibration by assigning the age of the sequenced samples to the tips within the human phylogenetic tree. But studies also vary greatly in the methodology employed and in the sequence panels analyzed, making it difficult to tease apart the causes for the differences between previous estimates. To clarify this issue, we compiled a comprehensive data set of 350 ancient and modern human complete mitochondrial DNA genomes, among which 146 were generated for the purpose of this study and estimated substitution rates using calibrations based both on dated nodes and tips. Our results demonstrate that, for the same data set, estimates based on individual dated tips are far more consistent with each other than those based on nodes and should thus be considered as more reliable.

Reconstructing the origin and spread of horse domestication in the Eurasian steppe.

Despite decades of research across multiple disciplines, the early history of horse domestication remains poorly understood. On the basis of current evidence from archaeology, mitochondrial DNA, and Y-chromosomal sequencing, a number of different domestication scenarios have been proposed, ranging from the spread of domestic horses out of a restricted primary area of domestication to the domestication of numerous distinct wild horse populations. In this paper, we reconstruct both the population genetic structure of the extinct wild progenitor of domestic horses, Equus ferus, and the origin and spread of horse domestication in the Eurasian steppes by fitting a spatially explicit stepping-stone model to genotype data from >300 horses sampled across northern Eurasia. We find strong evidence for an expansion of E. ferus out of eastern Eurasia about 160 kya, likely reflecting the colonization of Eurasia by this species. Our best-fitting scenario further suggests that horse domestication originated in the western part of the Eurasian steppe and that domestic herds were repeatedly restocked with local wild horses as they spread out of this area. By showing that horse domestication was initiated in the western Eurasian steppe and that the spread of domestic herds across Eurasia involved extensive introgression from the wild, the scenario of horse domestication proposed here unites evidence from archaeology, mitochondrial DNA, and Y-chromosomal DNA.

Ancient trade routes shaped the genetic structure of horses in eastern Eurasia.

Animal exchange networks have been shown to play an important role in determining gene flow among domestic animal populations. The Silk Road is one of the oldest continuous exchange networks in human history, yet its effectiveness in facilitating animal exchange across large geographical distances and topographically challenging landscapes has never been explicitly studied. Horses are known to have been traded along the Silk Roads; however, extensive movement of horses in connection with other human activities may have obscured the genetic signature of the Silk Roads. To investigate the role of the Silk Roads in shaping the genetic structure of horses in eastern Eurasia, we analysed microsatellite genotyping data from 455 village horses sampled from 17 locations. Using least-cost path methods, we compared the performance of models containing the Silk Roads as corridors for gene flow with models containing single landscape features. We also determined whether the recent isolation of former Soviet Union countries from the rest of Eurasia has affected the genetic structure of our samples. The overall level of genetic differentiation was low, consistent with historically high levels of gene flow across the study region. The spatial genetic structure was characterized by a significant, albeit weak, pattern of isolation by distance across the continent with no evidence for the presence of distinct genetic clusters. Incorporating landscape features considerably improved the fit of the data; however, when we controlled for geographical distance, only the correlation between genetic differentiation and the Silk Roads remained significant, supporting the effectiveness of this ancient trade network in facilitating gene flow across large geographical distances in a topographically complex landscape.

Autosomal genetic diversity in non-breed horses from eastern Eurasia provides insights into historical population movements.

Many events in the history of eastern Eurasia, including the process of domestication itself, the initial spread of domestic horses and subsequent movements, are believed to have affected the genetic structure of domestic horse populations in this area. We investigated levels of within- and between-population genetic diversity in 'non-breed horses' (working horses sampled in remote areas) from 17 locations in Asia and parts of Eastern Europe, using 26 autosomal microsatellite loci. Non-breed horses have not been subject to the same intensity of artificial selection and closed breeding as have most breed animals and are thus expected to better reflect the population history of domestic horses. Despite geographic distances of between 300 and 7000 km between sampling locations, pairwise F (ST) was very low (range: <0.001 to -0.033), suggesting historically high levels of gene flow. Our analyses of non-breed horses revealed a pattern of isolation by distance and a significant decline in genetic diversity (expected heterozygosity and allelic richness) from east to west, consistent with a westward expansion of horses out of East Asia. Although the timing of this putative expansion is unclear, our results highlight the benefit of studying animals that do not belong to particular breeds when investigating aspects of a population's history.

Human disturbance increases coronavirus prevalence in bats.

Human land modification is a known driver of animal-to-human transmission of infectious agents (zoonotic spillover). Infection prevalence in the reservoir is a key predictor of spillover, but landscape-level associations between the intensity of land modification and infection rates in wildlife remain largely untested. Bat-borne coronaviruses have caused three major disease outbreaks in humans: severe acute respiratory syndrome (SARS), Middle East respiratory syndrome, and coronavirus disease 2019 (COVID-19). We statistically link high-resolution land modification data with bat coronavirus surveillance records and show that coronavirus prevalence significantly increases with the intensity of human impact across all climates and levels of background biodiversity. The most significant contributors to the overall human impact are agriculture, deforestation, and mining. Regions of high predicted bat coronavirus prevalence coincide with global disease hotspots, suggesting that infection prevalence in wildlife may be an important factor underlying links between human land modification and zoonotic disease emergence.

RUNX2 tandem repeats and the evolution of facial length in placental mammals.

BACKGROUND: When simple sequence repeats are integrated into functional genes, they can potentially act as evolutionary 'tuning knobs', supplying abundant genetic variation with minimal risk of pleiotropic deleterious effects. The genetic basis of variation in facial shape and length represents a possible example of this phenomenon. Runt-related transcription factor 2 (RUNX2), which is involved in osteoblast differentiation, contains a functionally-important tandem repeat of glutamine and alanine amino acids. The ratio of glutamines to alanines (the QA ratio) in this protein seemingly influences the regulation of bone development. Notably, in domestic breeds of dog, and in carnivorans in general, the ratio of glutamines to alanines is strongly correlated with facial length. RESULTS: In this study we examine whether this correlation holds true across placental mammals, particularly those mammals for which facial length is highly variable and related to adaptive behavior and lifestyle (e.g., primates, afrotherians, xenarthrans). We obtained relative facial length measurements and RUNX2 sequences for 41 mammalian species representing 12 orders. Using both a phylogenetic generalized least squares model and a recently-developed Bayesian comparative method, we tested for a correlation between genetic and morphometric data while controlling for phylogeny, evolutionary rates, and divergence times. Non-carnivoran taxa generally had substantially lower glutamine-alanine ratios than carnivorans (primates and xenarthrans with means of 1.34 and 1.25, respectively, compared to a mean of 3.1 for carnivorans), and we found no correlation between RUNX2 sequence and face length across placental mammals. CONCLUSIONS: Results of our diverse comparative phylogenetic analyses indicate that QA ratio does not consistently correlate with face length across the 41 mammalian taxa considered. Thus, although RUNX2 might function as a 'tuning knob' modifying face length in carnivorans, this relationship is not conserved across mammals in general.

Discovery and population genomics of structural variation in a songbird genus.

Structural variation (SV) constitutes an important type of genetic mutations providing the raw material for evolution. Here, we uncover the genome-wide spectrum of intra- and interspecific SV segregating in natural populations of seven songbird species in the genus Corvus. Combining short-read (N = 127) and long-read re-sequencing (N = 31), as well as optical mapping (N = 16), we apply both assembly- and read mapping approaches to detect SV and characterize a total of 220,452 insertions, deletions and inversions. We exploit sampling across wide phylogenetic timescales to validate SV genotypes and assess the contribution of SV to evolutionary processes in an avian model of incipient speciation. We reveal an evolutionary young (~530,000 years) cis-acting 2.25-kb LTR retrotransposon insertion reducing expression of the NDP gene with consequences for premating isolation. Our results attest to the wealth and evolutionary significance of SV segregating in natural populations and highlight the need for reliable SV genotyping.

Genotype-free estimation of allele frequencies reduces bias and improves demographic inference from RADSeq data.

Restriction-site associated DNA sequencing (RADSeq) facilitates rapid generation of thousands of genetic markers at relatively low cost; however, several sources of error specific to RADSeq methods often lead to biased estimates of allele frequencies and thereby to erroneous population genetic inference. Estimating the distribution of sample allele frequencies without calling genotypes was shown to improve population inference from whole genome sequencing data, but the ability of this approach to account for RADSeq-specific biases remains unexplored. Here we assess in how far genotype-free methods of allele frequency estimation affect demographic inference from empirical RADSeq data. Using the well-studied pied flycatcher (Ficedula hypoleuca) as a study system, we compare allele frequency estimation and demographic inference from whole genome sequencing data with that from RADSeq data matched for samples using both genotype-based and genotype free methods. The demographic history of pied flycatchers as inferred from RADSeq data was highly congruent with that inferred from whole genome resequencing (WGS) data when allele frequencies were estimated directly from the read data. In contrast, when allele frequencies were derived from called genotypes, RADSeq-based estimates of most model parameters fell outside the 95% confidence interval of estimates derived from WGS data. Notably, more stringent filtering of the genotype calls tended to increase the discrepancy between parameter estimates from WGS and RADSeq data, respectively. The results from this study demonstrate the ability of genotype-free methods to improve allele frequency spectrum- (AFS-) based demographic inference from empirical RADSeq data and highlight the need to account for uncertainty in NGS data regardless of sequencing method.

Evolutionary analysis of the female-specific avian W chromosome.

The typically repetitive nature of the sex-limited chromosome means that it is often excluded from or poorly covered in genome assemblies, hindering studies of evolutionary and population genomic processes in non-recombining chromosomes. Here, we present a draft assembly of the non-recombining region of the collared flycatcher W chromosome, containing 46 genes without evidence of female-specific functional differentiation. Survival of genes during W chromosome degeneration has been highly non-random and expression data suggest that this can be attributed to selection for maintaining gene dose and ancestral expression levels of essential genes. Re-sequencing of large population samples revealed dramatically reduced levels of within-species diversity and elevated rates of between-species differentiation (lineage sorting), consistent with low effective population size. Concordance between W chromosome and mitochondrial DNA phylogenetic trees demonstrates evolutionary stable matrilineal inheritance of this nuclear-cytonuclear pair of chromosomes. Our results show both commonalities and differences between W chromosome and Y chromosome evolution.

POPULATION GENETICS. Genomic evidence for the Pleistocene and recent population history of Native Americans.

How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative "Paleoamerican" relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.

European domestic horses originated in two holocene refugia.

The role of European wild horses in horse domestication is poorly understood. While the fossil record for wild horses in Europe prior to horse domestication is scarce, there have been suggestions that wild populations from various European regions might have contributed to the gene pool of domestic horses. To distinguish between regions where domestic populations are mainly descended from local wild stock and those where horses were largely imported, we investigated patterns of genetic diversity in 24 European horse breeds typed at 12 microsatellite loci. The distribution of high levels of genetic diversity in Europe coincides with the distribution of predominantly open landscapes prior to domestication, as suggested by simulation-based vegetation reconstructions, with breeds from Iberia and the Caspian Sea region having significantly higher genetic diversity than breeds from central Europe and the UK, which were largely forested at the time the first domestic horses appear there. Our results suggest that not only the Eastern steppes, but also the Iberian Peninsula provided refugia for wild horses in the Holocene, and that the genetic contribution of these wild populations to local domestic stock may have been considerable. In contrast, the consistently low levels of diversity in central Europe and the UK suggest that domestic horses in these regions largely derive from horses that were imported from the Eastern refugium, the Iberian refugium, or both.