Genomic surveillance reveals dynamic shifts in the connectivity of COVID-19 epidemics.

The maturation of genomic surveillance in the past decade has enabled tracking of the emergence and spread of epidemics at an unprecedented level. During the COVID-19 pandemic, for example, genomic data revealed that local epidemics varied considerably in the frequency of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage importation and persistence, likely due to a combination of COVID-19 restrictions and changing connectivity. Here, we show that local COVID-19 epidemics are driven by regional transmission, including across international boundaries, but can become increasingly connected to distant locations following the relaxation of public health interventions. By integrating genomic, mobility, and epidemiological data, we find abundant transmission occurring between both adjacent and distant locations, supported by dynamic mobility patterns. We find that changing connectivity significantly influences local COVID-19 incidence. Our findings demonstrate a complex meaning of "local" when investigating connected epidemics and emphasize the importance of collaborative interventions for pandemic prevention and mitigation.

Human population history at the crossroads of East and Southeast Asia since 11,000 years ago.

Past human genetic diversity and migration between southern China and Southeast Asia have not been well characterized, in part due to poor preservation of ancient DNA in hot and humid regions. We sequenced 31 ancient genomes from southern China (Guangxi and Fujian), including two ∼12,000- to 10,000-year-old individuals representing the oldest humans sequenced from southern China. We discovered a deeply diverged East Asian ancestry in the Guangxi region that persisted until at least 6,000 years ago. We found that ∼9,000- to 6,000-year-old Guangxi populations were a mixture of local ancestry, southern ancestry previously sampled in Fujian, and deep Asian ancestry related to Southeast Asian Hòabìnhian hunter-gatherers, showing broad admixture in the region predating the appearance of farming. Historical Guangxi populations dating to ∼1,500 to 500 years ago are closely related to Tai-Kadai and Hmong-Mien speakers. Our results show heavy interactions among three distinct ancestries at the crossroads of East and Southeast Asia.

The deep population history of northern East Asia from the Late Pleistocene to the Holocene.

Northern East Asia was inhabited by modern humans as early as 40 thousand years ago (ka), as demonstrated by the Tianyuan individual. Using genome-wide data obtained from 25 individuals dated to 33.6-3.4 ka from the Amur region, we show that Tianyuan-related ancestry was widespread in northern East Asia before the Last Glacial Maximum (LGM). At the close of the LGM stadial, the earliest northern East Asian appeared in the Amur region, and this population is basal to ancient northern East Asians. Human populations in the Amur region have maintained genetic continuity from 14 ka, and these early inhabitants represent the closest East Asian source known for Ancient Paleo-Siberians. We also observed that EDAR V370A was likely to have been elevated to high frequency after the LGM, suggesting the possible timing for its selection. This study provides a deep look into the population dynamics of northern East Asia.

Ancient and modern genomes unravel the evolutionary history of the rhinoceros family.

Only five species of the once-diverse Rhinocerotidae remain, making the reconstruction of their evolutionary history a challenge to biologists since Darwin. We sequenced genomes from five rhinoceros species (three extinct and two living), which we compared to existing data from the remaining three living species and a range of outgroups. We identify an early divergence between extant African and Eurasian lineages, resolving a key debate regarding the phylogeny of extant rhinoceroses. This early Miocene (∼16 million years ago [mya]) split post-dates the land bridge formation between the Afro-Arabian and Eurasian landmasses. Our analyses also show that while rhinoceros genomes in general exhibit low levels of genome-wide diversity, heterozygosity is lowest and inbreeding is highest in the modern species. These results suggest that while low genetic diversity is a long-term feature of the family, it has been particularly exacerbated recently, likely reflecting recent anthropogenic-driven population declines.

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.

A Dynamic 6,000-Year Genetic History of Eurasia's Eastern Steppe.

The Eastern Eurasian Steppe was home to historic empires of nomadic pastoralists, including the Xiongnu and the Mongols. However, little is known about the region's population history. Here, we reveal its dynamic genetic history by analyzing new genome-wide data for 214 ancient individuals spanning 6,000 years. We identify a pastoralist expansion into Mongolia ca. 3000 BCE, and by the Late Bronze Age, Mongolian populations were biogeographically structured into three distinct groups, all practicing dairy pastoralism regardless of ancestry. The Xiongnu emerged from the mixing of these populations and those from surrounding regions. By comparison, the Mongols exhibit much higher eastern Eurasian ancestry, resembling present-day Mongolic-speaking populations. Our results illuminate the complex interplay between genetic, sociopolitical, and cultural changes on the Eastern Steppe.

Pan-Genome of Wild and Cultivated Soybeans.

Soybean is one of the most important vegetable oil and protein feed crops. To capture the entire genomic diversity, it is needed to construct a complete high-quality pan-genome from diverse soybean accessions. In this study, we performed individual de novo genome assemblies for 26 representative soybeans that were selected from 2,898 deeply sequenced accessions. Using these assembled genomes together with three previously reported genomes, we constructed a graph-based genome and performed pan-genome analysis, which identified numerous genetic variations that cannot be detected by direct mapping of short sequence reads onto a single reference genome. The structural variations from the 2,898 accessions that were genotyped based on the graph-based genome and the RNA sequencing (RNA-seq) data from the representative 26 accessions helped to link genetic variations to candidate genes that are responsible for important traits. This pan-genome resource will promote evolutionary and functional genomics studies in soybean.

Global Trends in Marine Plankton Diversity across Kingdoms of Life.

The ocean is home to myriad small planktonic organisms that underpin the functioning of marine ecosystems. However, their spatial patterns of diversity and the underlying drivers remain poorly known, precluding projections of their responses to global changes. Here we investigate the latitudinal gradients and global predictors of plankton diversity across archaea, bacteria, eukaryotes, and major virus clades using both molecular and imaging data from Tara Oceans. We show a decline of diversity for most planktonic groups toward the poles, mainly driven by decreasing ocean temperatures. Projections into the future suggest that severe warming of the surface ocean by the end of the 21st century could lead to tropicalization of the diversity of most planktonic groups in temperate and polar regions. These changes may have multiple consequences for marine ecosystem functioning and services and are expected to be particularly significant in key areas for carbon sequestration, fisheries, and marine conservation. VIDEO ABSTRACT.

Extensive Unexplored Human Microbiome Diversity Revealed by Over 150,000 Genomes from Metagenomes Spanning Age, Geography, and Lifestyle.

The body-wide human microbiome plays a role in health, but its full diversity remains uncharacterized, particularly outside of the gut and in international populations. We leveraged 9,428 metagenomes to reconstruct 154,723 microbial genomes (45% of high quality) spanning body sites, ages, countries, and lifestyles. We recapitulated 4,930 species-level genome bins (SGBs), 77% without genomes in public repositories (unknown SGBs [uSGBs]). uSGBs are prevalent (in 93% of well-assembled samples), expand underrepresented phyla, and are enriched in non-Westernized populations (40% of the total SGBs). We annotated 2.85 M genes in SGBs, many associated with conditions including infant development (94,000) or Westernization (106,000). SGBs and uSGBs permit deeper microbiome analyses and increase the average mappability of metagenomic reads from 67.76% to 87.51% in the gut (median 94.26%) and 65.14% to 82.34% in the mouth. We thus identify thousands of microbial genomes from yet-to-be-named species, expand the pangenomes of human-associated microbes, and allow better exploitation of metagenomic technologies.

Gene Expression Changes and Community Turnover Differentially Shape the Global Ocean Metatranscriptome.

Ocean microbial communities strongly influence the biogeochemistry, food webs, and climate of our planet. Despite recent advances in understanding their taxonomic and genomic compositions, little is known about how their transcriptomes vary globally. Here, we present a dataset of 187 metatranscriptomes and 370 metagenomes from 126 globally distributed sampling stations and establish a resource of 47 million genes to study community-level transcriptomes across depth layers from pole-to-pole. We examine gene expression changes and community turnover as the underlying mechanisms shaping community transcriptomes along these axes of environmental variation and show how their individual contributions differ for multiple biogeochemically relevant processes. Furthermore, we find the relative contribution of gene expression changes to be significantly lower in polar than in non-polar waters and hypothesize that in polar regions, alterations in community activity in response to ocean warming will be driven more strongly by changes in organismal composition than by gene regulatory mechanisms. VIDEO ABSTRACT.

Genomic Insights into Zika Virus Emergence and Spread.

The emergence and spread of Zika virus in the Americas continues to challenge our disease surveillance systems. Virus genome sequencing during the epidemic uncovered the timescale of Zika virus transmission and spread. Yet, we are only beginning to explore how genomics can enhance our responses to emerging viruses.

Geographic variation of mutagenic exposures in kidney cancer genomes.

International differences in the incidence of many cancer types indicate the existence of carcinogen exposures that have not yet been identified by conventional epidemiology make a substantial contribution to cancer burden1. In clear cell renal cell carcinoma, obesity, hypertension and tobacco smoking are risk factors, but they do not explain the geographical variation in its incidence2. Underlying causes can be inferred by sequencing the genomes of cancers from populations with different incidence rates and detecting differences in patterns of somatic mutations. Here we sequenced 962 clear cell renal cell carcinomas from 11 countries with varying incidence. The somatic mutation profiles differed between countries. In Romania, Serbia and Thailand, mutational signatures characteristic of aristolochic acid compounds were present in most cases, but these were rare elsewhere. In Japan, a mutational signature of unknown cause was found in more than 70% of cases but in less than 2% elsewhere. A further mutational signature of unknown cause was ubiquitous but exhibited higher mutation loads in countries with higher incidence rates of kidney cancer. Known signatures of tobacco smoking correlated with tobacco consumption, but no signature was associated with obesity or hypertension, suggesting that non-mutagenic mechanisms of action underlie these risk factors. The results of this study indicate the existence of multiple, geographically variable, mutagenic exposures that potentially affect tens of millions of people and illustrate the opportunities for new insights into cancer causation through large-scale global cancer genomics.

The geography of climate and the global patterns of species diversity.

Climate's effect on global biodiversity is typically viewed through the lens of temperature, humidity and resulting ecosystem productivity1-6. However, it is not known whether biodiversity depends solely on these climate conditions, or whether the size and fragmentation of these climates are also crucial. Here we shift the common perspective in global biodiversity studies, transitioning from geographic space to a climate-defined multidimensional space. Our findings suggest that larger and more isolated climate conditions tend to harbour higher diversity and species turnover among terrestrial tetrapods, encompassing more than 30,000 species. By considering both the characteristics of climate itself and its geographic attributes, we can explain almost 90% of the variation in global species richness. Half of the explanatory power (45%) may be attributed either to climate itself or to the geography of climate, suggesting a nuanced interplay between them. Our work evolves the conventional idea that larger climate regions, such as the tropics, host more species primarily because of their size7,8. Instead, we underscore the integral roles of both the geographic extent and degree of isolation of climates. This refined understanding presents a more intricate picture of biodiversity distribution, which can guide our approach to biodiversity conservation in an ever-changing world.

The giant diploid faba genome unlocks variation in a global protein crop.

Increasing the proportion of locally produced plant protein in currently meat-rich diets could substantially reduce greenhouse gas emissions and loss of biodiversity1. However, plant protein production is hampered by the lack of a cool-season legume equivalent to soybean in agronomic value2. Faba bean (Vicia faba L.) has a high yield potential and is well suited for cultivation in temperate regions, but genomic resources are scarce. Here, we report a high-quality chromosome-scale assembly of the faba bean genome and show that it has expanded to a massive 13 Gb in size through an imbalance between the rates of amplification and elimination of retrotransposons and satellite repeats. Genes and recombination events are evenly dispersed across chromosomes and the gene space is remarkably compact considering the genome size, although with substantial copy number variation driven by tandem duplication. Demonstrating practical application of the genome sequence, we develop a targeted genotyping assay and use high-resolution genome-wide association analysis to dissect the genetic basis of seed size and hilum colour. The resources presented constitute a genomics-based breeding platform for faba bean, enabling breeders and geneticists to accelerate the improvement of sustainable protein production across the Mediterranean, subtropical and northern temperate agroecological zones.

The universal visitation law of human mobility.

Human mobility impacts many aspects of a city, from its spatial structure1-3 to its response to an epidemic4-7. It is also ultimately key to social interactions8, innovation9,10 and productivity11. However, our quantitative understanding of the aggregate movements of individuals remains incomplete. Existing models-such as the gravity law12,13 or the radiation model14-concentrate on the purely spatial dependence of mobility flows and do not capture the varying frequencies of recurrent visits to the same locations. Here we reveal a simple and robust scaling law that captures the temporal and spatial spectrum of population movement on the basis of large-scale mobility data from diverse cities around the globe. According to this law, the number of visitors to any location decreases as the inverse square of the product of their visiting frequency and travel distance. We further show that the spatio-temporal flows to different locations give rise to prominent spatial clusters with an area distribution that follows Zipf's law15. Finally, we build an individual mobility model based on exploration and preferential return to provide a mechanistic explanation for the discovered scaling law and the emerging spatial structure. Our findings corroborate long-standing conjectures in human geography (such as central place theory16 and Weber's theory of emergent optimality10) and allow for predictions of recurrent flows, providing a basis for applications in urban planning, traffic engineering and the mitigation of epidemic diseases.

The importance of geography in forecasting future fire patterns under climate change.

An increasing amount of California's landscape has burned in wildfires in recent decades, in conjunction with increasing temperatures and vapor pressure deficit due to climate change. As the wildland-urban interface expands, more people are exposed to and harmed by these extensive wildfires, which are also eroding the resilience of terrestrial ecosystems. With future wildfire activity expected to increase, there is an urgent demand for solutions that sustain healthy ecosystems and wildfire-resilient human communities. Those who manage disaster response, landscapes, and biodiversity rely on mapped projections of how fire activity may respond to climate change and other human factors. California wildfire is complex, however, and climate-fire relationships vary across the state. Given known geographical variability in drivers of fire activity, we asked whether the geographical extent of fire models used to create these projections may alter the interpretation of predictions. We compared models of fire occurrence spanning the entire state of California to models developed for individual ecoregions and then projected end-of-century future fire patterns under climate change scenarios. We trained a Maximum Entropy model with fire records and hydroclimatological variables from recent decades (1981 to 2010) as well as topographic and human infrastructure predictors. Results showed substantial variation in predictors of fire probability and mapped future projections of fire depending upon geographical extents of model boundaries. Only the ecoregion models, accounting for the unique patterns of vegetation, climate, and human infrastructure, projected an increase in fire in most forested regions of the state, congruent with predictions from other studies.

Relational geographies of urban unsustainability: The entanglement of California's housing crisis with WUI growth and climate change.

One of California's most pressing social and environmental challenges is the rapid expansion of the wildlands-urban interface (WUI). Multiple issues associated with WUI growth compared to more dense and compact urban form are of concern-including greatly increased fire risk, greenhouse gas emissions, and fragmentation of habitat. However, little is understood about the factors driving this growth in the first place and, specifically, its relationship to urban-regional housing dynamics. This paper connects work in urban social science, urban and regional planning, and natural sciences to highlight the potential role of housing crises in driving displacement from the urban core to relatively more affordable exurbs, and with this, WUI growth. We analyze this relationship in California, which leads the nation in lack of affordable housing, scale of WUI growth, and many associated WUI hazards, including wildfire. We offer three related arguments: first, that California's affordable housing crisis, with its effect of driving migration to exurban areas, should be recognized as a significant urban form-related sustainability challenge; second, that to understand this challenge scholars must expand the spatial scale and analytic toolkit of both urban and WUI analysis through relational, mixed methods research; and third, that political and programmatic efforts to address California's housing crisis should undergird efforts to address WUI growth and climate change. Ultimately, we argue that expanding access to affordable urban housing can produce a more sustainable and just urban form that mitigates WUI-related climate and environmental impacts and reduces the vulnerability of growing numbers of WUI residents living in harm's way.

Metapopulation dynamics of SARS-CoV-2 transmission in a small-scale Amazonian society.

The severity of infectious disease outbreaks is governed by patterns of human contact, which vary by geography, social organization, mobility, access to technology and healthcare, economic development, and culture. Whereas globalized societies and urban centers exhibit characteristics that can heighten vulnerability to pandemics, small-scale subsistence societies occupying remote, rural areas may be buffered. Accordingly, voluntary collective isolation has been proposed as one strategy to mitigate the impacts of COVID-19 and other pandemics on small-scale Indigenous populations with minimal access to healthcare infrastructure. To assess the vulnerability of such populations and the viability of interventions such as voluntary collective isolation, we simulate and analyze the dynamics of SARS-CoV-2 infection among Amazonian forager-horticulturalists in Bolivia using a stochastic network metapopulation model parameterized with high-resolution empirical data on population structure, mobility, and contact networks. Our model suggests that relative isolation offers little protection at the population level (expected approximately 80% cumulative incidence), and more remote communities are not conferred protection via greater distance from outside sources of infection, due to common features of small-scale societies that promote rapid disease transmission such as high rates of travel and dense social networks. Neighborhood density, central household location in villages, and household size greatly increase the individual risk of infection. Simulated interventions further demonstrate that without implausibly high levels of centralized control, collective isolation is unlikely to be effective, especially if it is difficult to restrict visitation between communities as well as travel to outside areas. Finally, comparison of model results to empirical COVID-19 outcomes measured via seroassay suggest that our theoretical model is successful at predicting outbreak severity at both the population and community levels. Taken together, these findings suggest that the social organization and relative isolation from urban centers of many rural Indigenous communities offer little protection from pandemics and that standard control measures, including vaccination, are required to counteract effects of tight-knit social structures characteristic of small-scale populations.

NMPFamsDB: a database of novel protein families from microbial metagenomes and metatranscriptomes.

The Novel Metagenome Protein Families Database (NMPFamsDB) is a database of metagenome- and metatranscriptome-derived protein families, whose members have no hits to proteins of reference genomes or Pfam domains. Each protein family is accompanied by multiple sequence alignments, Hidden Markov Models, taxonomic information, ecosystem and geolocation metadata, sequence and structure predictions, as well as 3D structure models predicted with AlphaFold2. In its current version, NMPFamsDB hosts over 100 000 protein families, each with at least 100 members. The reported protein families significantly expand (more than double) the number of known protein sequence clusters from reference genomes and reveal new insights into their habitat distribution, origins, functions and taxonomy. We expect NMPFamsDB to be a valuable resource for microbial proteome-wide analyses and for further discovery and characterization of novel functions. NMPFamsDB is publicly available in http://www.nmpfamsdb.org/ or https://bib.fleming.gr/NMPFamsDB.

Quantification and modeling of turnover dynamics of de novo transcripts in Drosophila melanogaster.

Most of the transcribed eukaryotic genomes are composed of non-coding transcripts. Among these transcripts, some are newly transcribed when compared to outgroups and are referred to as de novo transcripts. De novo transcripts have been shown to play a major role in genomic innovations. However, little is known about the rates at which de novo transcripts are gained and lost in individuals of the same species. Here, we address this gap and estimate the de novo transcript turnover rate with an evolutionary model. We use DNA long reads and RNA short reads from seven geographically remote samples of inbred individuals of Drosophila melanogaster to detect de novo transcripts that are gained on a short evolutionary time scale. Overall, each sampled individual contains around 2500 unspliced de novo transcripts, with most of them being sample specific. We estimate that around 0.15 transcripts are gained per year, and that each gained transcript is lost at a rate around 5× 10-5 per year. This high turnover of transcripts suggests frequent exploration of new genomic sequences within species. These rate estimates are essential to comprehend the process and timescale of de novo gene birth.