Search

Understanding how severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in 2019 is critical to preventing future zoonotic outbreaks before they become the next pandemic. The Huanan Seafood Wholesale Market in Wuhan, China, was identified as a likely source of cases in early reports, but later this conclusion became controversial. We show here that the earliest known COVID-19 cases from December 2019, including those without reported direct links, were geographically centered on this market. We report that live SARS-CoV-2-susceptible mammals were sold at the market in late 2019 and that within the market, SARS-CoV-2-positive environmental samples were spatially associated with vendors selling live mammals. Although there is insufficient evidence to define upstream events, and exact circumstances remain obscure, our analyses indicate that the emergence of SARS-CoV-2 occurred through the live wildlife trade in China and show that the Huanan market was the epicenter of the COVID-19 pandemic.

Understanding the circumstances that lead to pandemics is important for their prevention. We analyzed the genomic diversity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) early in the coronavirus disease 2019 (COVID-19) pandemic. We show that SARS-CoV-2 genomic diversity before February 2020 likely comprised only two distinct viral lineages, denoted "A" and "B." Phylodynamic rooting methods, coupled with epidemic simulations, reveal that these lineages were the result of at least two separate cross-species transmission events into humans. The first zoonotic transmission likely involved lineage B viruses around 18 November 2019 (23 October to 8 December), and the separate introduction of lineage A likely occurred within weeks of this event. These findings indicate that it is unlikely that SARS-CoV-2 circulated widely in humans before November 2019 and define the narrow window between when SARS-CoV-2 first jumped into humans and when the first cases of COVID-19 were reported. As with other coronaviruses, SARS-CoV-2 emergence likely resulted from multiple zoonotic events.

Cerebrovascular diseases are a leading cause of death and neurologic disability. Further understanding of disease mechanisms and therapeutic strategies requires a deeper knowledge of cerebrovascular cells in humans. We profiled transcriptomes of 181,388 cells to define a cell atlas of the adult human cerebrovasculature, including endothelial cell molecular signatures with arteriovenous segmentation and expanded perivascular cell diversity. By leveraging this reference, we investigated cellular and molecular perturbations in brain arteriovenous malformations, which are a leading cause of stroke in young people, and identified pathologic endothelial transformations with abnormal vascular patterning and the ontology of vascularly derived inflammation. We illustrate the interplay between vascular and immune cells that contributes to brain hemorrhage and catalog opportunities for targeting angiogenic and inflammatory programs in vascular malformations.

Water-stable isotopes in polar ice cores are a widely used temperature proxy in paleoclimate reconstruction, yet calibration remains challenging in East Antarctica. Here, we reconstruct the magnitude and spatial pattern of Last Glacial Maximum surface cooling in Antarctica using borehole thermometry and firn properties in seven ice cores. West Antarctic sites cooled ~10°C relative to the preindustrial period. East Antarctic sites show a range from ~4° to ~7°C cooling, which is consistent with the results of global climate models when the effects of topographic changes indicated with ice core air-content data are included, but less than those indicated with the use of water-stable isotopes calibrated against modern spatial gradients. An altered Antarctic temperature inversion during the glacial reconciles our estimates with water-isotope observations.

The evolutionarily conserved splicing regulator neuro-oncological ventral antigen 1 (NOVA1) plays a key role in neural development and function. NOVA1 also includes a protein-coding difference between the modern human genome and Neanderthal and Denisovan genomes. To investigate the functional importance of an amino acid change in humans, we reintroduced the archaic allele into human induced pluripotent cells using genome editing and then followed their neural development through cortical organoids. This modification promoted slower development and higher surface complexity in cortical organoids with the archaic version of NOVA1 Moreover, levels of synaptic markers and synaptic protein coassociations correlated with altered electrophysiological properties in organoids expressing the archaic variant. Our results suggest that the human-specific substitution in NOVA1, which is exclusive to modern humans since divergence from Neanderthals, may have had functional consequences for our species' evolution.

Severe and persistent 21st-century drought in southwestern North America (SWNA) motivates comparisons to medieval megadroughts and questions about the role of anthropogenic climate change. We use hydrological modeling and new 1200-year tree-ring reconstructions of summer soil moisture to demonstrate that the 2000-2018 SWNA drought was the second driest 19-year period since 800 CE, exceeded only by a late-1500s megadrought. The megadrought-like trajectory of 2000-2018 soil moisture was driven by natural variability superimposed on drying due to anthropogenic warming. Anthropogenic trends in temperature, relative humidity, and precipitation estimated from 31 climate models account for 47% (model interquartiles of 35 to 105%) of the 2000-2018 drought severity, pushing an otherwise moderate drought onto a trajectory comparable to the worst SWNA megadroughts since 800 CE.

Inflammasomes are multiprotein platforms that initiate innate immunity by recruitment and activation of caspase-1. The NLRP1B inflammasome is activated upon direct cleavage by the anthrax lethal toxin protease. However, the mechanism by which cleavage results in NLRP1B activation is unknown. In this study, we find that cleavage results in proteasome-mediated degradation of the amino-terminal domains of NLRP1B, liberating a carboxyl-terminal fragment that is a potent caspase-1 activator. Proteasome-mediated degradation of NLRP1B is both necessary and sufficient for NLRP1B activation. Consistent with our functional degradation model, we identify IpaH7.8, a Shigella flexneri ubiquitin ligase secreted effector, as an enzyme that induces NLRP1B degradation and activation. Our results provide a unified mechanism for NLRP1B activation by diverse pathogen-encoded enzymatic activities.

To discover leads for next-generation chemoprotective antimalarial drugs, we tested more than 500,000 compounds for their ability to inhibit liver-stage development of luciferase-expressing Plasmodium spp. parasites (681 compounds showed a half-maximal inhibitory concentration of less than 1 micromolar). Cluster analysis identified potent and previously unreported scaffold families as well as other series previously associated with chemoprophylaxis. Further testing through multiple phenotypic assays that predict stage-specific and multispecies antimalarial activity distinguished compound classes that are likely to provide symptomatic relief by reducing asexual blood-stage parasitemia from those which are likely to only prevent malaria. Target identification by using functional assays, in vitro evolution, or metabolic profiling revealed 58 mitochondrial inhibitors but also many chemotypes possibly with previously unidentified mechanisms of action.

Previous detections of individual astrophysical sources of neutrinos are limited to the Sun and the supernova 1987A, whereas the origins of the diffuse flux of high-energy cosmic neutrinos remain unidentified. On 22 September 2017, we detected a high-energy neutrino, IceCube-170922A, with an energy of ~290 tera-electron volts. Its arrival direction was consistent with the location of a known γ-ray blazar, TXS 0506+056, observed to be in a flaring state. An extensive multiwavelength campaign followed, ranging from radio frequencies to γ-rays. These observations characterize the variability and energetics of the blazar and include the detection of TXS 0506+056 in very-high-energy γ-rays. This observation of a neutrino in spatial coincidence with a γ-ray-emitting blazar during an active phase suggests that blazars may be a source of high-energy neutrinos.

Chemogenetic characterization through in vitro evolution combined with whole-genome analysis can identify antimalarial drug targets and drug-resistance genes. We performed a genome analysis of 262 Plasmodium falciparum parasites resistant to 37 diverse compounds. We found 159 gene amplifications and 148 nonsynonymous changes in 83 genes associated with drug-resistance acquisition, where gene amplifications contributed to one-third of resistance acquisition events. Beyond confirming previously identified multidrug-resistance mechanisms, we discovered hitherto unrecognized drug target-inhibitor pairs, including thymidylate synthase and a benzoquinazolinone, farnesyltransferase and a pyrimidinedione, and a dipeptidylpeptidase and an arylurea. This exploration of the P. falciparum resistome and druggable genome will likely guide drug discovery and structural biology efforts, while also advancing our understanding of resistance mechanisms available to the malaria parasite.

To reveal impacts of sexual mode on genome content, we compared chromosome-scale assemblies of the outcrossing nematode Caenorhabditis nigoni to its self-fertile sibling species, C. briggsaeC. nigoni's genome resembles that of outcrossing relatives but encodes 31% more protein-coding genes than C. briggsaeC. nigoni genes lacking C. briggsae orthologs were disproportionately small and male-biased in expression. These include the male secreted short (mss) gene family, which encodes sperm surface glycoproteins conserved only in outcrossing species. Sperm from mss-null males of outcrossing C. remanei failed to compete with wild-type sperm, despite normal fertility in noncompetitive mating. Restoring mss to C. briggsae males was sufficient to enhance sperm competitiveness. Thus, sex has a pervasive influence on genome content that can be used to identify sperm competition factors.

Majorana fermion is a hypothetical particle that is its own antiparticle. We report transport measurements that suggest the existence of one-dimensional chiral Majorana fermion modes in the hybrid system of a quantum anomalous Hall insulator thin film coupled with a superconductor. As the external magnetic field is swept, half-integer quantized conductance plateaus are observed at the locations of magnetization reversals, giving a distinct signature of the Majorana fermion modes. This transport signature is reproducible over many magnetic field sweeps and appears at different temperatures. This finding may open up an avenue to control Majorana fermions for implementing robust topological quantum computing.

Conservation of species and ecosystems is increasingly difficult because anthropogenic impacts are pervasive and accelerating. Under this rapid global change, maximizing conservation success requires a paradigm shift from maintaining ecosystems in idealized past states toward facilitating their adaptive and functional capacities, even as species ebb and flow individually. Developing effective strategies under this new paradigm will require deeper understanding of the long-term dynamics that govern ecosystem persistence and reconciliation of conflicts among approaches to conserving historical versus novel ecosystems. Integrating emerging information from conservation biology, paleobiology, and the Earth sciences is an important step forward on the path to success. Maintaining nature in all its aspects will also entail immediately addressing the overarching threats of growing human population, overconsumption, pollution, and climate change.

Telomerase helps maintain telomeres by processive synthesis of telomere repeat DNA at their 3'-ends, using an integral telomerase RNA (TER) and telomerase reverse transcriptase (TERT). We report the cryo-electron microscopy structure of Tetrahymena telomerase at ~9 angstrom resolution. In addition to seven known holoenzyme proteins, we identify two additional proteins that form a complex (TEB) with single-stranded telomere DNA-binding protein Teb1, paralogous to heterotrimeric replication protein A (RPA). The p75-p45-p19 subcomplex is identified as another RPA-related complex, CST (CTC1-STN1-TEN1). This study reveals the paths of TER in the TERT-TER-p65 catalytic core and single-stranded DNA exit; extensive subunit interactions of the TERT essential N-terminal domain, p50, and TEB; and other subunit identities and structures, including p19 and p45C crystal structures. Our findings provide structural and mechanistic insights into telomerase holoenzyme function.

Bone morphogenetic proteins (BMPs) act in dose-dependent fashion to regulate cell fate choices in a myriad of developmental contexts. In early vertebrate and invertebrate embryos, BMPs and their antagonists establish epidermal versus central nervous system domains. In this highly conserved system, BMP antagonists mediate the neural-inductive activities proposed by Hans Spemann and Hilde Mangold nearly a century ago. BMPs distributed in gradients subsequently function as morphogens to subdivide the three germ layers into distinct territories and act to organize body axes, regulate growth, maintain stem cell niches, or signal inductively across germ layers. In this Review, we summarize the variety of mechanisms that contribute to generating reliable developmental responses to BMP gradients and other morphogen systems.

The causal mechanisms responsible for the abrupt climate changes of the Last Glacial Period remain unclear. One major difficulty is dating ice-rafted debris deposits associated with Heinrich events: Extensive iceberg influxes into the North Atlantic Ocean linked to global impacts on climate and biogeochemistry. In a new ice core record of atmospheric methane with ultrahigh temporal resolution, we find abrupt methane increases within Heinrich stadials 1, 2, 4, and 5 that, uniquely, have no counterparts in Greenland temperature proxies. Using a heuristic model of tropical rainfall distribution, we propose that Hudson Strait Heinrich events caused rainfall intensification over Southern Hemisphere land areas, thereby producing excess methane in tropical wetlands. Our findings suggest that the climatic impacts of Heinrich events persisted for 740 to 1520 years.

Evolutionary studies necessary to dissect diverse biological processes have been limited by the lack of reverse genetic approaches in most organisms with sequenced genomes. We established a broadly applicable strategy using zinc finger nucleases (ZFNs) and transcription activator-like effector nucleases (TALENs) for targeted disruption of endogenous genes and cis-acting regulatory elements in diverged nematode species.

To provide context for the diversification of archosaurs-the group that includes crocodilians, dinosaurs, and birds-we generated draft genomes of three crocodilians: Alligator mississippiensis (the American alligator), Crocodylus porosus (the saltwater crocodile), and Gavialis gangeticus (the Indian gharial). We observed an exceptionally slow rate of genome evolution within crocodilians at all levels, including nucleotide substitutions, indels, transposable element content and movement, gene family evolution, and chromosomal synteny. When placed within the context of related taxa including birds and turtles, this suggests that the common ancestor of all of these taxa also exhibited slow genome evolution and that the comparatively rapid evolution is derived in birds. The data also provided the opportunity to analyze heterozygosity in crocodilians, which indicates a likely reduction in population size for all three taxa through the Pleistocene. Finally, these data combined with newly published bird genomes allowed us to reconstruct the partial genome of the common ancestor of archosaurs, thereby providing a tool to investigate the genetic starting material of crocodilians, birds, and dinosaurs.

Birds are the most species-rich class of tetrapod vertebrates and have wide relevance across many research fields. We explored bird macroevolution using full genomes from 48 avian species representing all major extant clades. The avian genome is principally characterized by its constrained size, which predominantly arose because of lineage-specific erosion of repetitive elements, large segmental deletions, and gene loss. Avian genomes furthermore show a remarkably high degree of evolutionary stasis at the levels of nucleotide sequence, gene synteny, and chromosomal structure. Despite this pattern of conservation, we detected many non-neutral evolutionary changes in protein-coding genes and noncoding regions. These analyses reveal that pan-avian genomic diversity covaries with adaptations to different lifestyles and convergent evolution of traits.

Greenland ice core water isotopic composition (δ(18)O) provides detailed evidence for abrupt climate changes but is by itself insufficient for quantitative reconstruction of past temperatures and their spatial patterns. We investigate Greenland temperature evolution during the last deglaciation using independent reconstructions from three ice cores and simulations with a coupled ocean-atmosphere climate model. Contrary to the traditional δ(18)O interpretation, the Younger Dryas period was 4.5° ± 2°C warmer than the Oldest Dryas, due to increased carbon dioxide forcing and summer insolation. The magnitude of abrupt temperature changes is larger in central Greenland (9° to 14°C) than in the northwest (5° to 9°C), fingerprinting a North Atlantic origin. Simulated changes in temperature seasonality closely track changes in the Atlantic overturning strength and support the hypothesis that abrupt climate change is mostly a winter phenomenon.

During speech perception, linguistic elements such as consonants and vowels are extracted from a complex acoustic speech signal. The superior temporal gyrus (STG) participates in high-order auditory processing of speech, but how it encodes phonetic information is poorly understood. We used high-density direct cortical surface recordings in humans while they listened to natural, continuous speech to reveal the STG representation of the entire English phonetic inventory. At single electrodes, we found response selectivity to distinct phonetic features. Encoding of acoustic properties was mediated by a distributed population response. Phonetic features could be directly related to tuning for spectrotemporal acoustic cues, some of which were encoded in a nonlinear fashion or by integration of multiple cues. These findings demonstrate the acoustic-phonetic representation of speech in human STG.

The shape of the human face and skull is largely genetically determined. However, the genomic basis of craniofacial morphology is incompletely understood and hypothesized to involve protein-coding genes, as well as gene regulatory sequences. We used a combination of epigenomic profiling, in vivo characterization of candidate enhancer sequences in transgenic mice, and targeted deletion experiments to examine the role of distant-acting enhancers in craniofacial development. We identified complex regulatory landscapes consisting of enhancers that drive spatially complex developmental expression patterns. Analysis of mouse lines in which individual craniofacial enhancers had been deleted revealed significant alterations of craniofacial shape, demonstrating the functional importance of enhancers in defining face and skull morphology. These results demonstrate that enhancers are involved in craniofacial development and suggest that enhancer sequence variation contributes to the diversity of human facial morphology.

A rapidly growing body of research examines whether human conflict can be affected by climatic changes. Drawing from archaeology, criminology, economics, geography, history, political science, and psychology, we assemble and analyze the 60 most rigorous quantitative studies and document, for the first time, a striking convergence of results. We find strong causal evidence linking climatic events to human conflict across a range of spatial and temporal scales and across all major regions of the world. The magnitude of climate's influence is substantial: for each one standard deviation (1σ) change in climate toward warmer temperatures or more extreme rainfall, median estimates indicate that the frequency of interpersonal violence rises 4% and the frequency of intergroup conflict rises 14%. Because locations throughout the inhabited world are expected to warm 2σ to 4σ by 2050, amplified rates of human conflict could represent a large and critical impact of anthropogenic climate change.