Whole-genome sequencing reveals a complex African population demographic history and signatures of local adaptation.

We conduct high coverage (>30×) whole-genome sequencing of 180 individuals from 12 indigenous African populations. We identify millions of unreported variants, many predicted to be functionally important. We observe that the ancestors of southern African San and central African rainforest hunter-gatherers (RHG) diverged from other populations >200 kya and maintained a large effective population size. We observe evidence for ancient population structure in Africa and for multiple introgression events from "ghost" populations with highly diverged genetic lineages. Although currently geographically isolated, we observe evidence for gene flow between eastern and southern Khoesan-speaking hunter-gatherer populations lasting until ∼12 kya. We identify signatures of local adaptation for traits related to skin color, immune response, height, and metabolic processes. We identify a positively selected variant in the lightly pigmented San that influences pigmentation in vitro by regulating the enhancer activity and gene expression of PDPK1.

The genomic history of the Middle East.

The Middle East region is important to understand human evolution and migrations but is underrepresented in genomic studies. Here, we generated 137 high-coverage physically phased genome sequences from eight Middle Eastern populations using linked-read sequencing. We found no genetic traces of early expansions out-of-Africa in present-day populations but found Arabians have elevated Basal Eurasian ancestry that dilutes their Neanderthal ancestry. Population sizes within the region started diverging 15-20 kya, when Levantines expanded while Arabians maintained smaller populations that derived ancestry from local hunter-gatherers. Arabians suffered a population bottleneck around the aridification of Arabia 6 kya, while Levantines had a distinct bottleneck overlapping the 4.2 kya aridification event. We found an association between movement and admixture of populations in the region and the spread of Semitic languages. Finally, we identify variants that show evidence of selection, including polygenic selection. Our results provide detailed insights into the genomic and selective histories of the Middle East.

Predator mass mortality events restructure food webs through trophic decoupling.

Predators have a key role in structuring ecosystems1-4. However, predator loss is accelerating globally4-6, and predator mass-mortality events7 (MMEs)-rapid large-scale die-offs-are now emblematic of the Anthropocene epoch6. Owing to their rare and unpredictable nature7, we lack an understanding of how MMEs immediately impact ecosystems. Past predator-removal studies2,3 may be insufficient to understand the ecological consequences of MMEs because, in nature, dead predators decompose in situ and generate a resource pulse8, which could alter ensuing ecosystem dynamics by temporarily enhancing productivity. Here we experimentally induce MMEs in tritrophic, freshwater lake food webs and report ecological dynamics that are distinct from predator losses2,3 or resource pulses9 alone, but that can be predicted from theory8. MMEs led to the proliferation of diverse consumer and producer communities resulting from weakened top-down predator control1-3 and stronger bottom-up effects through predator decomposition8. In contrast to predator removals alone, enhanced primary production after MMEs dampened the consumer community response. As a consequence, MMEs generated biomass dynamics that were most similar to those of undisturbed systems, indicating that they may be cryptic disturbances in nature. These biomass dynamics led to trophic decoupling, whereby the indirect beneficial effects of predators on primary producers are lost and later materialize as direct bottom-up effects that stimulate primary production amid intensified herbivory. These results reveal ecological signatures of MMEs and demonstrate the feasibility of forecasting novel ecological dynamics arising with intensifying global change.

Disproportionate declines of formerly abundant species underlie insect loss.

Studies have reported widespread declines in terrestrial insect abundances in recent years1-4, but trends in other biodiversity metrics are less clear-cut5-7. Here we examined long-term trends in 923 terrestrial insect assemblages monitored in 106 studies, and found concomitant declines in abundance and species richness. For studies that were resolved to species level (551 sites in 57 studies), we observed a decline in the number of initially abundant species through time, but not in the number of very rare species. At the population level, we found that species that were most abundant at the start of the time series showed the strongest average declines (corrected for regression-to-the-mean effects). Rarer species were, on average, also declining, but these were offset by increases of other species. Our results suggest that the observed decreases in total insect abundance2 can mostly be explained by widespread declines of formerly abundant species. This counters the common narrative that biodiversity loss is mostly characterized by declines of rare species8,9. Although our results suggest that fundamental changes are occurring in insect assemblages, it is important to recognize that they represent only trends from those locations for which sufficient long-term data are available. Nevertheless, given the importance of abundant species in ecosystems10, their general declines are likely to have broad repercussions for food webs and ecosystem functioning.

Complexity of avian evolution revealed by family-level genomes.

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

Unequal climate impacts on global values of natural capital.

Ecosystems generate a wide range of benefits for humans, including some market goods as well as other benefits that are not directly reflected in market activity1. Climate change will alter the distribution of ecosystems around the world and change the flow of these benefits2,3. However, the specific implications of ecosystem changes for human welfare remain unclear, as they depend on the nature of these changes, the value of the affected benefits and the extent to which communities rely on natural systems for their well-being4. Here we estimate country-level changes in economic production and the value of non-market ecosystem benefits resulting from climate-change-induced shifts in terrestrial vegetation cover, as projected by dynamic global vegetation models (DGVMs) driven by general circulation climate models. Our results show that the annual population-weighted mean global flow of non-market ecosystem benefits valued in the wealth accounts of the World Bank will be reduced by 9.2% in 2100 under the Shared Socioeconomic Pathway SSP2-6.0 with respect to the baseline no climate change scenario and that the global population-weighted average change in gross domestic product (GDP) by 2100 is -1.3% of the baseline GDP. Because lower-income countries are more reliant on natural capital, these GDP effects are regressive. Approximately 90% of these damages are borne by the poorest 50% of countries and regions, whereas the wealthiest 10% experience only 2% of these losses.

Elevated genetic risk for multiple sclerosis emerged in steppe pastoralist populations.

Multiple sclerosis (MS) is a neuro-inflammatory and neurodegenerative disease that is most prevalent in Northern Europe. Although it is known that inherited risk for MS is located within or in close proximity to immune-related genes, it is unknown when, where and how this genetic risk originated1. Here, by using a large ancient genome dataset from the Mesolithic period to the Bronze Age2, along with new Medieval and post-Medieval genomes, we show that the genetic risk for MS rose among pastoralists from the Pontic steppe and was brought into Europe by the Yamnaya-related migration approximately 5,000 years ago. We further show that these MS-associated immunogenetic variants underwent positive selection both within the steppe population and later in Europe, probably driven by pathogenic challenges coinciding with changes in diet, lifestyle and population density. This study highlights the critical importance of the Neolithic period and Bronze Age as determinants of modern immune responses and their subsequent effect on the risk of developing MS in a changing environment.

Indigenous Australian genomes show deep structure and rich novel variation.

The Indigenous peoples of Australia have a rich linguistic and cultural history. How this relates to genetic diversity remains largely unknown because of their limited engagement with genomic studies. Here we analyse the genomes of 159 individuals from four remote Indigenous communities, including people who speak a language (Tiwi) not from the most widespread family (Pama-Nyungan). This large collection of Indigenous Australian genomes was made possible by careful community engagement and consultation. We observe exceptionally strong population structure across Australia, driven by divergence times between communities of 26,000-35,000 years ago and long-term low but stable effective population sizes. This demographic history, including early divergence from Papua New Guinean (47,000 years ago) and Eurasian groups1, has generated the highest proportion of previously undescribed genetic variation seen outside Africa and the most extended homozygosity compared with global samples. A substantial proportion of this variation is not observed in global reference panels or clinical datasets, and variation with predicted functional consequence is more likely to be homozygous than in other populations, with consequent implications for medical genomics2. Our results show that Indigenous Australians are not a single homogeneous genetic group and their genetic relationship with the peoples of New Guinea is not uniform. These patterns imply that the full breadth of Indigenous Australian genetic diversity remains uncharacterized, potentially limiting genomic medicine and equitable healthcare for Indigenous Australians.

The global wildland-urban interface.

The wildland-urban interface (WUI) is where buildings and wildland vegetation meet or intermingle1,2. It is where human-environmental conflicts and risks can be concentrated, including the loss of houses and lives to wildfire, habitat loss and fragmentation and the spread of zoonotic diseases3. However, a global analysis of the WUI has been lacking. Here, we present a global map of the 2020 WUI at 10 m resolution using a globally consistent and validated approach based on remote sensing-derived datasets of building area4 and wildland vegetation5. We show that the WUI is a global phenomenon, identify many previously undocumented WUI hotspots and highlight the wide range of population density, land cover types and biomass levels in different parts of the global WUI. The WUI covers only 4.7% of the land surface but is home to nearly half its population (3.5 billion). The WUI is especially widespread in Europe (15% of the land area) and the temperate broadleaf and mixed forests biome (18%). Of all people living near 2003-2020 wildfires (0.4 billion), two thirds have their home in the WUI, most of them in Africa (150 million). Given that wildfire activity is predicted to increase because of climate change in many regions6, there is a need to understand housing growth and vegetation patterns as drivers of WUI change.

Continent-wide declines in shallow reef life over a decade of ocean warming.

Human society is dependent on nature1,2, but whether our ecological foundations are at risk remains unknown in the absence of systematic monitoring of species' populations3. Knowledge of species fluctuations is particularly inadequate in the marine realm4. Here we assess the population trends of 1,057 common shallow reef species from multiple phyla at 1,636 sites around Australia over the past decade. Most populations decreased over this period, including many tropical fishes, temperate invertebrates (particularly echinoderms) and southwestern Australian macroalgae, whereas coral populations remained relatively stable. Population declines typically followed heatwave years, when local water temperatures were more than 0.5 °C above temperatures in 2008. Following heatwaves5,6, species abundances generally tended to decline near warm range edges, and increase near cool range edges. More than 30% of shallow invertebrate species in cool latitudes exhibited high extinction risk, with rapidly declining populations trapped by deep ocean barriers, preventing poleward retreat as temperatures rise. Greater conservation effort is needed to safeguard temperate marine ecosystems, which are disproportionately threatened and include species with deep evolutionary roots. Fundamental among such efforts, and broader societal needs to efficiently adapt to interacting anthropogenic and natural pressures, is greatly expanded monitoring of species' population trends7,8.

Evolutionary Genomics of High Fecundity.

Natural highly fecund populations abound. These range from viruses to gadids. Many highly fecund populations are economically important. Highly fecund populations provide an important contrast to the low-fecundity organisms that have traditionally been applied in evolutionary studies. A key question regarding high fecundity is whether large numbers of offspring are produced on a regular basis, by few individuals each time, in a sweepstakes mode of reproduction. Such reproduction characteristics are not incorporated into the classical Wright-Fisher model, the standard reference model of population genetics, or similar types of models, in which each individual can produce only small numbers of offspring relative to the population size. The expected genomic footprints of population genetic models of sweepstakes reproduction are very different from those of the Wright-Fisher model. A key, immediate issue involves identifying the footprints of sweepstakes reproduction in genomic data. Whole-genome sequencing data can be used to distinguish the patterns made by sweepstakes reproduction from the patterns made by population growth in a population evolving according to the Wright-Fisher model (or similar models). If the hypothesis of sweepstakes reproduction cannot be rejected, then models of sweepstakes reproduction and associated multiple-merger coalescents will become at least as relevant as the Wright-Fisher model (or similar models) and the Kingman coalescent, the cornerstones of mathematical population genetics, in further discussions of evolutionary genomics of highly fecund populations.

Lidar reveals pre-Hispanic low-density urbanism in the Bolivian Amazon.

Archaeological remains of agrarian-based, low-density urbananism1-3 have been reported to exist beneath the tropical forests of Southeast Asia, Sri Lanka and Central America4-6. However, beyond some large interconnected settlements in southern Amazonia7-9, there has been no such evidence for pre-Hispanic Amazonia. Here we present lidar data of sites belonging to the Casarabe culture (around AD 500 to AD 1400)10-13 in the Llanos de Mojos savannah-forest mosaic, southwest Amazonia, revealing the presence of two remarkably large sites (147 ha and 315 ha) in a dense four-tiered settlement system. The Casarabe culture area, as far as known today, spans approximately 4,500 km2, with one of the large settlement sites controlling an area of approximately 500 km2. The civic-ceremonial architecture of these large settlement sites includes stepped platforms, on top of which lie U-shaped structures, rectangular platform mounds and conical pyramids (which are up to 22 m tall). The large settlement sites are surrounded by ranked concentric polygonal banks and represent central nodes that are connected to lower-ranked sites by straight, raised causeways that stretch over several kilometres. Massive water-management infrastructure, composed of canals and reservoirs, complete the settlement system in an anthropogenically modified landscape. Our results indicate that the Casarabe-culture settlement pattern represents a type of tropical low-density urbanism that has not previously been described in Amazonia.

Africa's oldest dinosaurs reveal early suppression of dinosaur distribution.

The vertebrate lineages that would shape Mesozoic and Cenozoic terrestrial ecosystems originated across Triassic Pangaea1-11. By the Late Triassic (Carnian stage, ~235 million years ago), cosmopolitan 'disaster faunas' (refs. 12-14) had given way to highly endemic assemblages12,13 on the supercontinent. Testing the tempo and mode of the establishment of this endemism is challenging-there were few geographic barriers to dispersal across Pangaea during the Late Triassic. Instead, palaeolatitudinal climate belts, and not continental boundaries, are proposed to have controlled distribution15-18. During this time of high endemism, dinosaurs began to disperse and thus offer an opportunity to test the timing and drivers of this biogeographic pattern. Increased sampling can test this prediction: if dinosaurs initially dispersed under palaeolatitudinal-driven endemism, then an assemblage similar to those of South America4,19-21 and India19,22-including the earliest dinosaurs-should be present in Carnian deposits in south-central Africa. Here we report a new Carnian assemblage from Zimbabwe that includes Africa's oldest definitive dinosaurs, including a nearly complete skeleton of the sauropodomorph Mbiresaurus raathi gen. et sp. nov. This assemblage resembles other dinosaur-bearing Carnian assemblages, suggesting that a similar vertebrate fauna ranged high-latitude austral Pangaea. The distribution of the first dinosaurs is correlated with palaeolatitude-linked climatic barriers, and dinosaurian dispersal to the rest of the supercontinent was delayed until these barriers relaxed, suggesting that climatic controls influenced the initial composition of the terrestrial faunas that persist to this day.

Tropical tree mortality has increased with rising atmospheric water stress.

Evidence exists that tree mortality is accelerating in some regions of the tropics1,2, with profound consequences for the future of the tropical carbon sink and the global anthropogenic carbon budget left to limit peak global warming below 2 °C. However, the mechanisms that may be driving such mortality changes and whether particular species are especially vulnerable remain unclear3-8. Here we analyse a 49-year record of tree dynamics from 24 old-growth forest plots encompassing a broad climatic gradient across the Australian moist tropics and find that annual tree mortality risk has, on average, doubled across all plots and species over the last 35 years, indicating a potential halving in life expectancy and carbon residence time. Associated losses in biomass were not offset by gains from growth and recruitment. Plots in less moist local climates presented higher average mortality risk, but local mean climate did not predict the pace of temporal increase in mortality risk. Species varied in the trajectories of their mortality risk, with the highest average risk found nearer to the upper end of the atmospheric vapour pressure deficit niches of species. A long-term increase in vapour pressure deficit was evident across the region, suggesting that thresholds involving atmospheric water stress, driven by global warming, may be a primary cause of increasing tree mortality in moist tropical forests.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

A genetic history of the pre-contact Caribbean.

Humans settled the Caribbean about 6,000 years ago, and ceramic use and intensified agriculture mark a shift from the Archaic to the Ceramic Age at around 2,500 years ago1-3. Here we report genome-wide data from 174 ancient individuals from The Bahamas, Haiti and the Dominican Republic (collectively, Hispaniola), Puerto Rico, Curaçao and Venezuela, which we co-analysed with 89 previously published ancient individuals. Stone-tool-using Caribbean people, who first entered the Caribbean during the Archaic Age, derive from a deeply divergent population that is closest to Central and northern South American individuals; contrary to previous work4, we find no support for ancestry contributed by a population related to North American individuals. Archaic-related lineages were >98% replaced by a genetically homogeneous ceramic-using population related to speakers of languages in the Arawak family from northeast South America; these people moved through the Lesser Antilles and into the Greater Antilles at least 1,700 years ago, introducing ancestry that is still present. Ancient Caribbean people avoided close kin unions despite limited mate pools that reflect small effective population sizes, which we estimate to be a minimum of 500-1,500 and a maximum of 1,530-8,150 individuals on the combined islands of Puerto Rico and Hispaniola in the dozens of generations before the individuals who we analysed lived. Census sizes are unlikely to be more than tenfold larger than effective population sizes, so previous pan-Caribbean estimates of hundreds of thousands of people are too large5,6. Confirming a small and interconnected Ceramic Age population7, we detect 19 pairs of cross-island cousins, close relatives buried around 75 km apart in Hispaniola and low genetic differentiation across islands. Genetic continuity across transitions in pottery styles reveals that cultural changes during the Ceramic Age were not driven by migration of genetically differentiated groups from the mainland, but instead reflected interactions within an interconnected Caribbean world1,8.

Cryptic transmission of SARS-CoV-2 and the first COVID-19 wave.

Considerable uncertainty surrounds the timeline of introductions and onsets of local transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) globally1-7. Although a limited number of SARS-CoV-2 introductions were reported in January and February 2020 (refs.8,9), the narrowness of the initial testing criteria, combined with a slow growth in testing capacity and porous travel screening10, left many countries vulnerable to unmitigated, cryptic transmission. Here we use a global metapopulation epidemic model to provide a mechanistic understanding of the early dispersal of infections and the temporal windows of the introduction of SARS-CoV-2 and onset of local transmission in Europe and the USA. We find that community transmission of SARS-CoV-2 was likely to have been present in several areas of Europe and the USA by January 2020, and estimate that by early March, only 1 to 4 in 100 SARS-CoV-2 infections were detected by surveillance systems. The modelling results highlight international travel as the key driver of the introduction of SARS-CoV-2, with possible introductions and transmission events as early as December 2019 to January 2020. We find a heterogeneous geographic distribution of cumulative infection attack rates by 4 July 2020, ranging from 0.78% to 15.2% across US states and 0.19% to 13.2% in European countries. Our approach complements phylogenetic analyses and other surveillance approaches and provides insights that can be used to design innovative, model-driven surveillance systems that guide enhanced testing and response strategies.

A dynamical limit to evolutionary adaptation.

Natural selection makes evolutionary adaptation possible even if the overwhelming majority of new mutations are deleterious. However, in rapidly evolving populations where numerous linked mutations occur and segregate simultaneously, clonal interference and genetic hitchhiking can limit the efficiency of selection, allowing deleterious mutations to accumulate over time. This can in principle overwhelm the fitness increases provided by beneficial mutations, leading to an overall fitness decline. Here, we analyze the conditions under which evolution will tend to drive populations to higher versus lower fitness. Our analysis focuses on quantifying the boundary between these two regimes, as a function of parameters such as population size, mutation rates, and selection pressures. This boundary represents a state in which adaptation is precisely balanced by Muller's ratchet, and we show that it can be characterized by rapid molecular evolution without any net fitness change. Finally, we consider the implications of global fitness-mediated epistasis and find that under some circumstances, this can drive populations toward the boundary state, which can thus represent a long-term evolutionary attractor.

High germline mutation rates, but not extreme population outbreaks, influence genetic diversity in a keystone coral predator.

Lewontin's paradox, the observation that levels of genetic diversity (π) do not scale linearly with census population size (Nc) variation, is an evolutionary conundrum. The most extreme mismatches between π and Nc are found for highly abundant marine invertebrates. Yet, the influences of new mutations on π relative to extrinsic processes such as Nc fluctuations are unknown. Here, we provide the first germline mutation rate (µ) estimate for a marine invertebrate in corallivorous crown-of-thorns sea stars (Acanthaster cf. solaris). We use high-coverage whole-genome sequencing of 14 parent-offspring trios alongside empirical estimates of Nc in Australia's Great Barrier Reef to jointly examine the determinants of π in populations undergoing extreme Nc fluctuations. The A. cf. solaris mean µ was 9.13 x 10-09 mutations per-site per-generation (95% CI: 6.51 x 10-09 to 1.18 x 10-08), exceeding estimates for other invertebrates and showing greater concordance with vertebrate mutation rates. Lower-than-expected Ne (~70,000-180,000) and low Ne/Nc values (0.0047-0.048) indicated weak influences of population outbreaks on long-term π. Our findings are consistent with elevated µ evolving in response to reduced Ne and generation time length, with important implications for explaining high mutational loads and the determinants of genetic diversity in marine invertebrate taxa.

The long-term expansion and recession of human populations.

Over the last 12,000 y, human populations have expanded and transformed critical earth systems. Yet, a key unresolved question in the environmental and social sciences remains: Why did human populations grow and, sometimes, decline in the first place? Our research builds on 20 y of archaeological research studying the deep time dynamics of human populations to propose an explanation for the long-term growth and stability of human populations. Innovations in the productive capacity of populations fuels exponential-like growth over thousands of years; however, innovations saturate over time and, often, may leave populations vulnerable to large recessions in their well-being and population density. Empirically, we find a trade-off between changes in land use that increase the production and consumption of carbohydrates, driving repeated waves of population growth over thousands of years, and the susceptibility of populations to large recessions due to a lag in the impact of humans on resources. These results shed light on the long-term drivers of human population growth and decline.