Massive expansion of human gut bacteriophage diversity.

Bacteriophages drive evolutionary change in bacterial communities by creating gene flow networks that fuel ecological adaptions. However, the extent of viral diversity and its prevalence in the human gut remains largely unknown. Here, we introduce the Gut Phage Database, a collection of ∼142,000 non-redundant viral genomes (>10 kb) obtained by mining a dataset of 28,060 globally distributed human gut metagenomes and 2,898 reference genomes of cultured gut bacteria. Host assignment revealed that viral diversity is highest in the Firmicutes phyla and that ∼36% of viral clusters (VCs) are not restricted to a single species, creating gene flow networks across phylogenetically distinct bacterial species. Epidemiological analysis uncovered 280 globally distributed VCs found in at least 5 continents and a highly prevalent phage clade with features reminiscent of p-crAssphage. This high-quality, large-scale catalog of phage genomes will improve future virome studies and enable ecological and evolutionary analysis of human gut bacteriophages.

SnapShot: Microbial Extremophiles.

Extremophiles are remarkable examples of life's resilience, thriving in hot springs at boiling temperatures, in brine lakes saturated with salt, and in the driest deserts. We review the biogeography, currently known limits of life, and molecular adaptations to extremes. See the online interactive map for additional detail on biogeography, environmental microbiology, and exemplary species. To view this SnapShot, open or download the PDF.

Charting the Complexity of the Marine Microbiome through Single-Cell Genomics.

Marine bacteria and archaea play key roles in global biogeochemistry. To improve our understanding of this complex microbiome, we employed single-cell genomics and a randomized, hypothesis-agnostic cell selection strategy to recover 12,715 partial genomes from the tropical and subtropical euphotic ocean. A substantial fraction of known prokaryoplankton coding potential was recovered from a single, 0.4 mL ocean sample, which indicates that genomic information disperses effectively across the globe. Yet, we found each genome to be unique, implying limited clonality within prokaryoplankton populations. Light harvesting and secondary metabolite biosynthetic pathways were numerous across lineages, highlighting the value of single-cell genomics to advance the identification of ecological roles and biotechnology potential of uncultured microbial groups. This genome collection enabled functional annotation and genus-level taxonomic assignments for >80% of individual metagenome reads from the tropical and subtropical surface ocean, thus offering a model to improve reference genome databases for complex microbiomes.

Relaxed Selection Limits Lifespan by Increasing Mutation Load.

To uncover the selective forces shaping life-history trait evolution across species, we investigate the genomic basis underlying adaptations to seasonal habitat desiccation in African killifishes, identifying the genetic variants associated with positive and relaxed purifying selection in 45 killifish species and 231 wild individuals distributed throughout sub-Saharan Africa. In annual species, genetic drift led to the expansion of nuclear and mitochondrial genomes and caused the accumulation of deleterious genetic variants in key life-history modulating genes such as mtor, insr, ampk, foxo3, and polg. Relaxation of purifying selection is also significantly associated with mitochondrial function and aging in human populations. We find that relaxation of purifying selection prominently shapes genomes and is a prime candidate force molding the evolution of lifespan and the distribution of genetic variants associated with late-onset diseases in different species. VIDEO ABSTRACT.

From "A"IV to "Z"IKV: Attacks from Emerging and Re-emerging Pathogens.

100 years after the infamous "Spanish flu" pandemic, the 2017-2018 flu season has been severe, with numerous infections worldwide. In between, there have been continuous, relentless attacks from (re-)emerging viruses. To fully understand viral pathogenesis and develop effective medical countermeasures, we must strengthen current surveillance and basic research efforts.

Mutualisms weaken the latitudinal diversity gradient among oceanic islands.

The latitudinal diversity gradient (LDG) dominates global patterns of diversity1,2, but the factors that underlie the LDG remain elusive. Here we use a unique global dataset3 to show that vascular plants on oceanic islands exhibit a weakened LDG and explore potential mechanisms for this effect. Our results show that traditional physical drivers of island biogeography4-namely area and isolation-contribute to the difference between island and mainland diversity at a given latitude (that is, the island species deficit), as smaller and more distant islands experience reduced colonization. However, plant species with mutualists are underrepresented on islands, and we find that this plant mutualism filter explains more variation in the island species deficit than abiotic factors. In particular, plant species that require animal pollinators or microbial mutualists such as arbuscular mycorrhizal fungi contribute disproportionately to the island species deficit near the Equator, with contributions decreasing with distance from the Equator. Plant mutualist filters on species richness are particularly strong at low absolute latitudes where mainland richness is highest, weakening the LDG of oceanic islands. These results provide empirical evidence that mutualisms, habitat heterogeneity and dispersal are key to the maintenance of high tropical plant diversity and mediate the biogeographic patterns of plant diversity on Earth.

Biogeographic response of marine plankton to Cenozoic environmental changes.

In palaeontological studies, groups with consistent ecological and morphological traits across a clade's history (functional groups)1 afford different perspectives on biodiversity dynamics than do species and genera2,3, which are evolutionarily ephemeral. Here we analyse Triton, a global dataset of Cenozoic macroperforate planktonic foraminiferal occurrences4, to contextualize changes in latitudinal equitability gradients1, functional diversity, palaeolatitudinal specialization and community equitability. We identify: global morphological communities becoming less specialized preceding the richness increase after the Cretaceous-Palaeogene extinction; ecological specialization during the Early Eocene Climatic Optimum, suggesting inhibitive equatorial temperatures during the peak of the Cenozoic hothouse; increased specialization due to circulation changes across the Eocene-Oligocene transition, preceding the loss of morphological diversity; changes in morphological specialization and richness about 19 million years ago, coeval with pelagic shark extinctions5; delayed onset of changing functional group richness and specialization between hemispheres during the mid-Miocene plankton diversification. The detailed nature of the Triton dataset permits a unique spatiotemporal view of Cenozoic pelagic macroevolution, in which global biogeographic responses of functional communities and richness are decoupled during Cenozoic climate events. The global response of functional groups to similar abiotic selection pressures may depend on the background climatic state (greenhouse or icehouse) to which a group is adapted.

Amazon forest biogeography predicts resilience and vulnerability to drought.

Amazonia contains the most extensive tropical forests on Earth, but Amazon carbon sinks of atmospheric CO2 are declining, as deforestation and climate-change-associated droughts1-4 threaten to push these forests past a tipping point towards collapse5-8. Forests exhibit complex drought responses, indicating both resilience (photosynthetic greening) and vulnerability (browning and tree mortality), that are difficult to explain by climate variation alone9-17. Here we combine remotely sensed photosynthetic indices with ground-measured tree demography to identify mechanisms underlying drought resilience/vulnerability in different intact forest ecotopes18,19 (defined by water-table depth, soil fertility and texture, and vegetation characteristics). In higher-fertility southern Amazonia, drought response was structured by water-table depth, with resilient greening in shallow-water-table forests (where greater water availability heightened response to excess sunlight), contrasting with vulnerability (browning and excess tree mortality) over deeper water tables. Notably, the resilience of shallow-water-table forest weakened as drought lengthened. By contrast, lower-fertility northern Amazonia, with slower-growing but hardier trees (or, alternatively, tall forests, with deep-rooted water access), supported more-drought-resilient forests independent of water-table depth. This functional biogeography of drought response provides a framework for conservation decisions and improved predictions of heterogeneous forest responses to future climate changes, warning that Amazonia's most productive forests are also at greatest risk, and that longer/more frequent droughts are undermining multiple ecohydrological strategies and capacities for Amazon forest resilience.

Revising the global biogeography of annual and perennial plants.

There are two main life cycles in plants-annual and perennial1,2. These life cycles are associated with different traits that determine ecosystem function3,4. Although life cycles are textbook examples of plant adaptation to different environments, we lack comprehensive knowledge regarding their global distributional patterns. Here we assembled an extensive database of plant life cycle assignments of 235,000 plant species coupled with millions of georeferenced datapoints to map the worldwide biogeography of these plant species. We found that annual plants are half as common as initially thought5-8, accounting for only 6% of plant species. Our analyses indicate that annuals are favoured in hot and dry regions. However, a more accurate model shows that the prevalence of annual species is driven by temperature and precipitation in the driest quarter (rather than yearly means), explaining, for example, why some Mediterranean systems have more annuals than desert systems. Furthermore, this pattern remains consistent among different families, indicating convergent evolution. Finally, we demonstrate that increasing climate variability and anthropogenic disturbance increase annual favourability. Considering future climate change, we predict an increase in annual prevalence for 69% of the world's ecoregions by 2060. Overall, our analyses raise concerns for ecosystem services provided by perennial plants, as ongoing changes are leading to a higher proportion of annual plants globally.

Late Cenozoic cooling restructured global marine plankton communities.

The geographic ranges of marine organisms, including planktonic foraminifera1, diatoms, dinoflagellates2, copepods3 and fish4, are shifting polewards owing to anthropogenic climate change5. However, the extent to which species will move and whether these poleward range shifts represent precursor signals that lead to extinction is unclear6. Understanding the development of marine biodiversity patterns over geological time and the factors that influence them are key to contextualizing these current trends. The fossil record of the macroperforate planktonic foraminifera provides a rich and phylogenetically resolved dataset that provides unique opportunities for understanding marine biogeography dynamics and how species distributions have responded to ancient climate changes. Here we apply a bipartite network approach to quantify group diversity, latitudinal specialization and latitudinal equitability for planktonic foraminifera over the past eight million years using Triton, a recently developed high-resolution global dataset of planktonic foraminiferal occurrences7. The results depict a global, clade-wide shift towards the Equator in ecological and morphological community equitability over the past eight million years in response to temperature changes during the late Cenozoic bipolar ice sheet formation. Collectively, the Triton data indicate the presence of a latitudinal equitability gradient among planktonic foraminiferal functional groups which is coupled to the latitudinal biodiversity gradient only through the geologically recent past (the past two million years). Before this time, latitudinal equitability gradients indicate that higher latitudes promoted community equitability across ecological and morphological groups. Observed range shifts among marine planktonic microorganisms1,2,8 in the recent and geological past suggest substantial poleward expansion of marine communities even under the most conservative future global warming scenarios.

Origination of the modern-style diversity gradient 15 million years ago.

The latitudinal diversity gradient (LDG) is a prevalent feature of modern ecosystems across diverse clades1-4. Recognized for well over a century, the causal mechanisms for LDGs remain disputed, in part because numerous putative drivers simultaneously covary with latitude1,3,5. The past provides the opportunity to disentangle LDG mechanisms because the relationships among biodiversity, latitude and possible causal factors have varied over time6-9. Here we quantify the emergence of the LDG in planktonic foraminifera at high spatiotemporal resolution over the past 40 million years, finding that a modern-style gradient arose only 15 million years ago. Spatial and temporal models suggest that LDGs for planktonic foraminifera may be controlled by the physical structure of the water column. Steepening of the latitudinal temperature gradient over 15 million years ago, associated with an increased vertical temperature gradient at low latitudes, may have enhanced niche partitioning and provided more opportunities for speciation at low latitudes. Supporting this hypothesis, we find that higher rates of low-latitude speciation steepened the diversity gradient, consistent with spatiotemporal patterns of depth partitioning by planktonic foraminifera. Extirpation of species from high latitudes also strengthened the LDG, but this effect tended to be weaker than speciation. Our results provide a step change in understanding the evolution of marine LDGs over long timescales.

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.

Climate change increases cross-species viral transmission risk.

At least 10,000 virus species have the ability to infect humans but, at present, the vast majority are circulating silently in wild mammals1,2. However, changes in climate and land use will lead to opportunities for viral sharing among previously geographically isolated species of wildlife3,4. In some cases, this will facilitate zoonotic spillover-a mechanistic link between global environmental change and disease emergence. Here we simulate potential hotspots of future viral sharing, using a phylogeographical model of the mammal-virus network, and projections of geographical range shifts for 3,139 mammal species under climate-change and land-use scenarios for the year 2070. We predict that species will aggregate in new combinations at high elevations, in biodiversity hotspots, and in areas of high human population density in Asia and Africa, causing the cross-species transmission of their associated viruses an estimated 4,000 times. Owing to their unique dispersal ability, bats account for the majority of novel viral sharing and are likely to share viruses along evolutionary pathways that will facilitate future emergence in humans. Notably, we find that this ecological transition may already be underway, and holding warming under 2 °C within the twenty-first century will not reduce future viral sharing. Our findings highlight an urgent need to pair viral surveillance and discovery efforts with biodiversity surveys tracking the range shifts of species, especially in tropical regions that contain the most zoonoses and are experiencing rapid warming.

Opportunities and challenges of macrogenetic studies.

The rapidly emerging field of macrogenetics focuses on analysing publicly accessible genetic datasets from thousands of species to explore large-scale patterns and predictors of intraspecific genetic variation. Facilitated by advances in evolutionary biology, technology, data infrastructure, statistics and open science, macrogenetics addresses core evolutionary hypotheses (such as disentangling environmental and life-history effects on genetic variation) with a global focus. Yet, there are important, often overlooked, limitations to this approach and best practices need to be considered and adopted if macrogenetics is to continue its exciting trajectory and reach its full potential in fields such as biodiversity monitoring and conservation. Here, we review the history of this rapidly growing field, highlight knowledge gaps and future directions, and provide guidelines for further research.

A Triassic stem lepidosaur illuminates the origin of lizard-like reptiles.

The early evolution of diapsid reptiles is marked by a deep contrast between our knowledge of the origin and early evolution of archosauromorphs (crocodiles, avian and non-avian dinosaurs) to that of lepidosauromorphs (squamates (lizards, snakes) and sphenodontians (tuataras)). Whereas the former include hundreds of fossil species across various lineages during the Triassic period1, the latter are represented by an extremely patchy early fossil record comprising only a handful of fragmentary fossils, most of which have uncertain phylogenetic affinities and are confined to Europe1-3. Here we report the discovery of a three-dimensionally preserved reptile skull, assigned as Taytalura alcoberi gen. et sp. nov., from the Late Triassic epoch of Argentina that is robustly inferred phylogenetically as the earliest evolving lepidosauromorph, using various data types and optimality criteria. Micro-computed tomography scans of this skull reveal details about the origin of the lepidosaurian skull from early diapsids, suggesting that several traits traditionally associated with sphenodontians in fact originated much earlier in lepidosauromorph evolution. Taytalura suggests that the strongly evolutionarily conserved skull architecture of sphenodontians represents the plesiomorphic condition for all lepidosaurs, that stem and crown lepidosaurs were contemporaries for at least ten million years during the Triassic, and that early lepidosauromorphs had a much broader geographical distribution than has previously been thought.

The genomic origins of the Bronze Age Tarim Basin mummies.

The identity of the earliest inhabitants of Xinjiang, in the heart of Inner Asia, and the languages that they spoke have long been debated and remain contentious1. Here we present genomic data from 5 individuals dating to around 3000-2800 BC from the Dzungarian Basin and 13 individuals dating to around 2100-1700 BC from the Tarim Basin, representing the earliest yet discovered human remains from North and South Xinjiang, respectively. We find that the Early Bronze Age Dzungarian individuals exhibit a predominantly Afanasievo ancestry with an additional local contribution, and the Early-Middle Bronze Age Tarim individuals contain only a local ancestry. The Tarim individuals from the site of Xiaohe further exhibit strong evidence of milk proteins in their dental calculus, indicating a reliance on dairy pastoralism at the site since its founding. Our results do not support previous hypotheses for the origin of the Tarim mummies, who were argued to be Proto-Tocharian-speaking pastoralists descended from the Afanasievo1,2 or to have originated among the Bactria-Margiana Archaeological Complex3 or Inner Asian Mountain Corridor cultures4. Instead, although Tocharian may have been plausibly introduced to the Dzungarian Basin by Afanasievo migrants during the Early Bronze Age, we find that the earliest Tarim Basin cultures appear to have arisen from a genetically isolated local population that adopted neighbouring pastoralist and agriculturalist practices, which allowed them to settle and thrive along the shifting riverine oases of the Taklamakan Desert.

Emergence and expansion of SARS-CoV-2 B.1.526 after identification in New York.

SARS-CoV-2 infections have surged across the globe in recent months, concomitant with considerable viral evolution1-3. Extensive mutations in the spike protein may threaten the efficacy of vaccines and therapeutic monoclonal antibodies4. Two signature spike mutations of concern are E484K, which has a crucial role in the loss of neutralizing activity of antibodies, and N501Y, a driver of rapid worldwide transmission of the B.1.1.7 lineage. Here we report the emergence of the variant lineage B.1.526 (also known as the Iota variant5), which contains E484K, and its rise to dominance in New York City in early 2021. This variant is partially or completely resistant to two therapeutic monoclonal antibodies that are in clinical use and is less susceptible to neutralization by plasma from individuals who had recovered from SARS-CoV-2 infection or serum from vaccinated individuals, posing a modest antigenic challenge. The presence of the B.1.526 lineage has now been reported in all 50 states in the United States and in many other countries. B.1.526 rapidly replaced earlier lineages in New York, with an estimated transmission advantage of 35%. These transmission dynamics, together with the relative antibody resistance of its E484K sub-lineage, are likely to have contributed to the sharp rise and rapid spread of B.1.526. Although SARS-CoV-2 B.1.526 initially outpaced B.1.1.7 in the region, its growth subsequently slowed concurrently with the rise of B.1.1.7 and ensuing variants.

Detection of a SARS-CoV-2 variant of concern in South Africa.

Continued uncontrolled transmission of SARS-CoV-2 in many parts of the world is creating conditions for substantial evolutionary changes to the virus1,2. Here we describe a newly arisen lineage of SARS-CoV-2 (designated 501Y.V2; also known as B.1.351 or 20H) that is defined by eight mutations in the spike protein, including three substitutions (K417N, E484K and N501Y) at residues in its receptor-binding domain that may have functional importance3-5. This lineage was identified in South Africa after the first wave of the epidemic in a severely affected metropolitan area (Nelson Mandela Bay) that is located on the coast of the Eastern Cape province. This lineage spread rapidly, and became dominant in Eastern Cape, Western Cape and KwaZulu-Natal provinces within weeks. Although the full import of the mutations is yet to be determined, the genomic data-which show rapid expansion and displacement of other lineages in several regions-suggest that this lineage is associated with a selection advantage that most plausibly results from increased transmissibility or immune escape6-8.

Untangling introductions and persistence in COVID-19 resurgence in Europe.

After the first wave of SARS-CoV-2 infections in spring 2020, Europe experienced a resurgence of the virus starting in late summer 2020 that was deadlier and more difficult to contain1. Relaxed intervention measures and summer travel have been implicated as drivers of the second wave2. Here we build a phylogeographical model to evaluate how newly introduced lineages, as opposed to the rekindling of persistent lineages, contributed to the resurgence of COVID-19 in Europe. We inform this model using genomic, mobility and epidemiological data from 10 European countries and estimate that in many countries more than half of the lineages circulating in late summer resulted from new introductions since 15 June 2020. The success in onward transmission of newly introduced lineages was negatively associated with the local incidence of COVID-19 during this period. The pervasive spread of variants in summer 2020 highlights the threat of viral dissemination when restrictions are lifted, and this needs to be carefully considered in strategies to control the current spread of variants that are more transmissible and/or evade immunity. Our findings indicate that more effective and coordinated measures are required to contain the spread through cross-border travel even as vaccination is reducing disease burden.

West Nile virus spread in Europe: Phylogeographic pattern analysis and key drivers.

BACKGROUND: West Nile virus (WNV) outbreaks in birds, humans, and livestock have occurred in multiple areas in Europe and have had a significant impact on animal and human health. The patterns of emergence and spread of WNV in Europe are very different from those in the US and understanding these are important for guiding preparedness activities. METHODS: We mapped the evolution and spread history of WNV in Europe by incorporating viral genome sequences and epidemiological data into phylodynamic models. Spatially explicit phylogeographic models were developed to explore the possible contribution of different drivers to viral dispersal direction and velocity. A "skygrid-GLM" approach was used to identify how changes in environments would predict viral genetic diversity variations over time. FINDINGS: Among the six lineages found in Europe, WNV-2a (a sub-lineage of WNV-2) has been predominant (accounting for 73% of all sequences obtained in Europe that have been shared in the public domain) and has spread to at least 14 countries. In the past two decades, WNV-2a has evolved into two major co-circulating clusters, both originating from Central Europe, but with distinct dynamic history and transmission patterns. WNV-2a spreads at a high dispersal velocity (88km/yr-215 km/yr) which is correlated to bird movements. Notably, amongst multiple drivers that could affect the spread of WNV, factors related to land use were found to strongly influence the spread of WNV. Specifically, the intensity of agricultural activities (defined by factors related to crops and livestock production, such as coverage of cropland, pasture, cultivated and managed vegetation, livestock density) were positively associated with both spread direction and velocity. In addition, WNV spread direction was associated with high coverage of wetlands and migratory bird flyways. CONCLUSION: Our results suggest that-in addition to ecological conditions favouring bird- and mosquito- presence-agricultural land use may be a significant driver of WNV emergence and spread. Our study also identified significant gaps in data and the need to strengthen virological surveillance in countries of Central Europe from where WNV outbreaks are likely seeded. Enhanced monitoring for early detection of further dispersal could be targeted to areas with high agricultural activities and habitats of migratory birds.