Search

One of the justifiable criticisms of human genetic studies is the underrepresentation of participants from diverse populations. Lack of inclusion must be addressed at-scale to identify causal disease factors and understand the genetic causes of health disparities. We present genome-wide associations for 2068 traits from 635,969 participants in the Department of Veterans Affairs Million Veteran Program, a longitudinal study of diverse United States Veterans. Systematic analysis revealed 13,672 genomic risk loci; 1608 were only significant after including non-European populations. Fine-mapping identified causal variants at 6318 signals across 613 traits. One-third (n = 2069) were identified in participants from non-European populations. This reveals a broadly similar genetic architecture across populations, highlights genetic insights gained from underrepresented groups, and presents an extensive atlas of genetic associations.

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.

Nucleotide changes in gene regulatory elements are important determinants of neuronal development and diseases. Using massively parallel reporter assays in primary human cells from mid-gestation cortex and cerebral organoids, we interrogated the cis-regulatory activity of 102,767 open chromatin regions, including thousands of sequences with cell type-specific accessibility and variants associated with brain gene regulation. In primary cells, we identified 46,802 active enhancer sequences and 164 variants that alter enhancer activity. Activity was comparable in organoids and primary cells, suggesting that organoids provide an adequate model for the developing cortex. Using deep learning we decoded the sequence basis and upstream regulators of enhancer activity. This work establishes a comprehensive catalog of functional gene regulatory elements and variants in human neuronal development.

One-carbon metabolism is an essential branch of cellular metabolism that intersects with epigenetic regulation. In this work, we show how formaldehyde (FA), a one-carbon unit derived from both endogenous sources and environmental exposure, regulates one-carbon metabolism by inhibiting the biosynthesis of S-adenosylmethionine (SAM), the major methyl donor in cells. FA reacts with privileged, hyperreactive cysteine sites in the proteome, including Cys120 in S-adenosylmethionine synthase isoform type-1 (MAT1A). FA exposure inhibited MAT1A activity and decreased SAM production with MAT-isoform specificity. A genetic mouse model of chronic FA overload showed a decrease n SAM and in methylation on selected histones and genes. Epigenetic and transcriptional regulation of Mat1a and related genes function as compensatory mechanisms for FA-dependent SAM depletion, revealing a biochemical feedback cycle between FA and SAM one-carbon units.

The metal-halide ionic octahedron is the optoelectronic unit for halide perovskites, and a crown ether-assisted supramolecular assembly approach can pack various ionic octahedra into tunable symmetries. In this work, we demonstrate near-unity photoluminescence quantum yield (PLQY) blue and green emission with the supramolecular assembly of hafnium (Hf) and zirconium (Zr) halide octahedral clusters. (18C6@K)2HfBr6 powders showed blue emission with a near-unity PLQY (96.2%), and green emission was also achieved with (18C6@K)2ZrCl4Br2 powders at a PLQY of 82.7%. These highly emissive powders feature facile low-temperature solution-based synthesis conditions and maintain high PLQY in solution-processable semiconductor inks under ambient conditions, and they were used in thin-film displays and emissive three-dimensional-printed architectures that exhibited high spatial resolution.

In nature, nonheme iron enzymes use dioxygen to generate high-spin iron(IV)=O species for a variety of oxygenation reactions. Although synthetic chemists have long sought to mimic this reactivity, the enzyme-like activation of O2 to form high-spin iron(IV) = O species remains an unrealized goal. Here, we report a metal-organic framework featuring iron(II) sites with a local structure similar to that in α-ketoglutarate-dependent dioxygenases. The framework reacts with O2 at low temperatures to form high-spin iron(IV) = O species that are characterized using in situ diffuse reflectance infrared Fourier transform, in situ and variable-field Mössbauer, Fe Kβ x-ray emission, and nuclear resonance vibrational spectroscopies. In the presence of O2, the framework is competent for catalytic oxygenation of cyclohexane and the stoichiometric conversion of ethane to ethanol.

Among mammals, post-reproductive life spans are currently documented only in humans and a few species of toothed whales. Here we show that a post-reproductive life span exists among wild chimpanzees in the Ngogo community of Kibale National Park, Uganda. Post-reproductive representation was 0.195, indicating that a female who reached adulthood could expect to live about one-fifth of her adult life in a post-reproductive state, around half as long as human hunter-gatherers. Post-reproductive females exhibited hormonal signatures of menopause, including sharply increasing gonadotropins after age 50. We discuss whether post-reproductive life spans in wild chimpanzees occur only rarely, as a short-term response to favorable ecological conditions, or instead are an evolved species-typical trait as well as the implications of these alternatives for our understanding of the evolution of post-reproductive life spans.

Recent advances in single-cell transcriptomics have illuminated the diverse neuronal and glial cell types within the human brain. However, the regulatory programs governing cell identity and function remain unclear. Using a single-nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq), we explored open chromatin landscapes across 1.1 million cells in 42 brain regions from three adults. Integrating this data unveiled 107 distinct cell types and their specific utilization of 544,735 candidate cis-regulatory DNA elements (cCREs) in the human genome. Nearly a third of the cCREs demonstrated conservation and chromatin accessibility in the mouse brain cells. We reveal strong links between specific brain cell types and neuropsychiatric disorders including schizophrenia, bipolar disorder, Alzheimers disease (AD), and major depression, and have developed deep learning models to predict the regulatory roles of noncoding risk variants in these disorders.

Delineating the gene-regulatory programs underlying complex cell types is fundamental for understanding brain function in health and disease. Here, we comprehensively examined human brain cell epigenomes by probing DNA methylation and chromatin conformation at single-cell resolution in 517 thousand cells (399 thousand neurons and 118 thousand non-neurons) from 46 regions of three adult male brains. We identified 188 cell types and characterized their molecular signatures. Integrative analyses revealed concordant changes in DNA methylation, chromatin accessibility, chromatin organization, and gene expression across cell types, cortical areas, and basal ganglia structures. We further developed single-cell methylation barcodes that reliably predict brain cell types using the methylation status of select genomic sites. This multimodal epigenomic brain cell atlas provides new insights into the complexity of cell-type-specific gene regulation in adult human brains.

The precise pattern and timing of speciation events that gave rise to all living placental mammals remain controversial. We provide a comprehensive phylogenetic analysis of genetic variation across an alignment of 241 placental mammal genome assemblies, addressing prior concerns regarding limited genomic sampling across species. We compared neutral genome-wide phylogenomic signals using concatenation and coalescent-based approaches, interrogated phylogenetic variation across chromosomes, and analyzed extensive catalogs of structural variants. Interordinal relationships exhibit relatively low rates of phylogenomic conflict across diverse datasets and analytical methods. Conversely, X-chromosome versus autosome conflicts characterize multiple independent clades that radiated during the Cenozoic. Genomic time trees reveal an accumulation of cladogenic events before and immediately after the Cretaceous-Paleogene (K-Pg) boundary, implying important roles for Cretaceous continental vicariance and the K-Pg extinction in the placental radiation.

We reconstruct the phenotype of Balto, the heroic sled dog renowned for transporting diphtheria antitoxin to Nome, Alaska, in 1925, using evolutionary constraint estimates from the Zoonomia alignment of 240 mammals and 682 genomes from dogs and wolves of the 21st century. Balto shares just part of his diverse ancestry with the eponymous Siberian husky breed. Balto's genotype predicts a combination of coat features atypical for modern sled dog breeds, and a slightly smaller stature. He had enhanced starch digestion compared with Greenland sled dogs and a compendium of derived homozygous coding variants at constrained positions in genes connected to bone and skin development. We propose that Balto's population of origin, which was less inbred and genetically healthier than that of modern breeds, was adapted to the extreme environment of 1920s Alaska.

Species persistence can be influenced by the amount, type, and distribution of diversity across the genome, suggesting a potential relationship between historical demography and resilience. In this study, we surveyed genetic variation across single genomes of 240 mammals that compose the Zoonomia alignment to evaluate how historical effective population size (Ne) affects heterozygosity and deleterious genetic load and how these factors may contribute to extinction risk. We find that species with smaller historical Ne carry a proportionally larger burden of deleterious alleles owing to long-term accumulation and fixation of genetic load and have a higher risk of extinction. This suggests that historical demography can inform contemporary resilience. Models that included genomic data were predictive of species' conservation status, suggesting that, in the absence of adequate census or ecological data, genomic information may provide an initial risk assessment.

Conserved genomic sequences disrupted in humans may underlie uniquely human phenotypic traits. We identified and characterized 10,032 human-specific conserved deletions (hCONDELs). These short (average 2.56 base pairs) deletions are enriched for human brain functions across genetic, epigenomic, and transcriptomic datasets. Using massively parallel reporter assays in six cell types, we discovered 800 hCONDELs conferring significant differences in regulatory activity, half of which enhance rather than disrupt regulatory function. We highlight several hCONDELs with putative human-specific effects on brain development, including HDAC5, CPEB4, and PPP2CA. Reverting an hCONDEL to the ancestral sequence alters the expression of LOXL2 and developmental genes involved in myelination and synaptic function. Our data provide a rich resource to investigate the evolutionary mechanisms driving new traits in humans and other species.

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

Human accelerated regions (HARs) are conserved genomic loci that evolved at an accelerated rate in the human lineage and may underlie human-specific traits. We generated HARs and chimpanzee accelerated regions with an automated pipeline and an alignment of 241 mammalian genomes. Combining deep learning with chromatin capture experiments in human and chimpanzee neural progenitor cells, we discovered a significant enrichment of HARs in topologically associating domains containing human-specific genomic variants that change three-dimensional (3D) genome organization. Differential gene expression between humans and chimpanzees at these loci suggests rewiring of regulatory interactions between HARs and neurodevelopmental genes. Thus, comparative genomics together with models of 3D genome folding revealed enhancer hijacking as an explanation for the rapid evolution of HARs.

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Molecular clocks are the basis for dating the divergence between lineages over macroevolutionary timescales (~105 to 108 years). However, classical DNA-based clocks tick too slowly to inform us about the recent past. Here, we demonstrate that stochastic DNA methylation changes at a subset of cytosines in plant genomes display a clocklike behavior. This "epimutation clock" is orders of magnitude faster than DNA-based clocks and enables phylogenetic explorations on a scale of years to centuries. We show experimentally that epimutation clocks recapitulate known topologies and branching times of intraspecies phylogenetic trees in the self-fertilizing plant Arabidopsis thaliana and the clonal seagrass Zostera marina, which represent two major modes of plant reproduction. This discovery will open new possibilities for high-resolution temporal studies of plant biodiversity.

Water is essential for survival, but one in three individuals worldwide (2.2 billion people) lacks access to safe drinking water. Water intake requirements largely reflect water turnover (WT), the water used by the body each day. We investigated the determinants of human WT in 5604 people from the ages of 8 days to 96 years from 23 countries using isotope-tracking (2H) methods. Age, body size, and composition were significantly associated with WT, as were physical activity, athletic status, pregnancy, socioeconomic status, and environmental characteristics (latitude, altitude, air temperature, and humidity). People who lived in countries with a low human development index (HDI) had higher WT than people in high-HDI countries. On the basis of this extensive dataset, we provide equations to predict human WT in relation to anthropometric, economic, and environmental factors.

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.

Loss of nuclear TDP-43 is a hallmark of neurodegeneration in TDP-43 proteinopathies, including amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). TDP-43 mislocalization results in cryptic splicing and polyadenylation of pre-messenger RNAs (pre-mRNAs) encoding stathmin-2 (also known as SCG10), a protein that is required for axonal regeneration. We found that TDP-43 binding to a GU-rich region sterically blocked recognition of the cryptic 3 splice site in STMN2 pre-mRNA. Targeting dCasRx or antisense oligonucleotides (ASOs) suppressed cryptic splicing, which restored axonal regeneration and stathmin-2-dependent lysosome trafficking in TDP-43-deficient human motor neurons. In mice that were gene-edited to contain human STMN2 cryptic splice-polyadenylation sequences, ASO injection into cerebral spinal fluid successfully corrected Stmn2 pre-mRNA misprocessing and restored stathmin-2 expression levels independently of TDP-43 binding.

Honey bees use a complex form of spatial referential communication. Their "waggle dance" communicates the direction, distance, and quality of a resource to nestmates by encoding celestial cues, retinal optic flow, and relative food value into motion and sound within the nest. We show that correct waggle dancing requires social learning. Bees without the opportunity to follow any dances before they first danced produced significantly more disordered dances with larger waggle angle divergence errors and encoded distance incorrectly. The former deficit improved with experience, but distance encoding was set for life. The first dances of bees that could follow other dancers showed neither impairment. Social learning, therefore, shapes honey bee signaling, as it does communication in human infants, birds, and multiple other vertebrate species.

Cell-cell interactions in the central nervous system play important roles in neurologic diseases. However, little is known about the specific molecular pathways involved, and methods for their systematic identification are limited. Here, we developed a forward genetic screening platform that combines CRISPR-Cas9 perturbations, cell coculture in picoliter droplets, and microfluidic-based fluorescence-activated droplet sorting to identify mechanisms of cell-cell communication. We used SPEAC-seq (systematic perturbation of encapsulated associated cells followed by sequencing), in combination with in vivo genetic perturbations, to identify microglia-produced amphiregulin as a suppressor of disease-promoting astrocyte responses in multiple sclerosis preclinical models and clinical samples. Thus, SPEAC-seq enables the high-throughput systematic identification of cell-cell communication mechanisms.

The induction of proinflammatory T cells by dendritic cell (DC) subtypes is critical for antitumor responses and effective immune checkpoint blockade (ICB) therapy. Here, we show that human CD1c+CD5+ DCs are reduced in melanoma-affected lymph nodes, with CD5 expression on DCs correlating with patient survival. Activating CD5 on DCs enhanced T cell priming and improved survival after ICB therapy. CD5+ DC numbers increased during ICB therapy, and low interleukin-6 (IL-6) concentrations promoted their de novo differentiation. Mechanistically, CD5 expression by DCs was required to generate optimally protective CD5hi T helper and CD8+ T cells; further, deletion of CD5 from T cells dampened tumor elimination in response to ICB therapy in vivo. Thus, CD5+ DCs are an essential component of optimal ICB therapy.

Molecular characterization of cell types using single-cell transcriptome sequencing is revolutionizing cell biology and enabling new insights into the physiology of human organs. We created a human reference atlas comprising nearly 500,000 cells from 24 different tissues and organs, many from the same donor. This atlas enabled molecular characterization of more than 400 cell types, their distribution across tissues, and tissue-specific variation in gene expression. Using multiple tissues from a single donor enabled identification of the clonal distribution of T cells between tissues, identification of the tissue-specific mutation rate in B cells, and analysis of the cell cycle state and proliferative potential of shared cell types across tissues. Cell type-specific RNA splicing was discovered and analyzed across tissues within an individual.

In cases of severe wildlife population decline, a key question is whether recovery efforts will be impeded by genetic factors, such as inbreeding depression. Decades of excess mortality from gillnet fishing have driven Mexico's vaquita porpoise (Phocoena sinus) to ~10 remaining individuals. We analyzed whole-genome sequences from 20 vaquitas and integrated genomic and demographic information into stochastic, individual-based simulations to quantify the species' recovery potential. Our analysis suggests that the vaquita's historical rarity has resulted in a low burden of segregating deleterious variation, reducing the risk of inbreeding depression. Similarly, genome-informed simulations suggest that the vaquita can recover if bycatch mortality is immediately halted. This study provides hope for vaquitas and other naturally rare endangered species and highlights the utility of genomics in predicting extinction risk.

The mass of the W boson, a mediator of the weak force between elementary particles, is tightly constrained by the symmetries of the standard model of particle physics. The Higgs boson was the last missing component of the model. After observation of the Higgs boson, a measurement of the W boson mass provides a stringent test of the model. We measure the W boson mass, MW, using data corresponding to 8.8 inverse femtobarns of integrated luminosity collected in proton-antiproton collisions at a 1.96 tera-electron volt center-of-mass energy with the CDF II detector at the Fermilab Tevatron collider. A sample of approximately 4 million W boson candidates is used to obtain [Formula: see text], the precision of which exceeds that of all previous measurements combined (stat, statistical uncertainty; syst, systematic uncertainty; MeV, mega-electron volts; c, speed of light in a vacuum). This measurement is in significant tension with the standard model expectation.

Compared to its predecessors, the Telomere-to-Telomere CHM13 genome adds nearly 200 million base pairs of sequence, corrects thousands of structural errors, and unlocks the most complex regions of the human genome for clinical and functional study. We show how this reference universally improves read mapping and variant calling for 3202 and 17 globally diverse samples sequenced with short and long reads, respectively. We identify hundreds of thousands of variants per sample in previously unresolved regions, showcasing the promise of the T2T-CHM13 reference for evolutionary and biomedical discovery. Simultaneously, this reference eliminates tens of thousands of spurious variants per sample, including reduction of false positives in 269 medically relevant genes by up to a factor of 12. Because of these improvements in variant discovery coupled with population and functional genomic resources, T2T-CHM13 is positioned to replace GRCh38 as the prevailing reference for human genetics.

Since its initial release in 2000, the human reference genome has covered only the euchromatic fraction of the genome, leaving important heterochromatic regions unfinished. Addressing the remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium presents a complete 3.055 billion-base pair sequence of a human genome, T2T-CHM13, that includes gapless assemblies for all chromosomes except Y, corrects errors in the prior references, and introduces nearly 200 million base pairs of sequence containing 1956 gene predictions, 99 of which are predicted to be protein coding. The completed regions include all centromeric satellite arrays, recent segmental duplications, and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies.

From the preparation of pharmaceuticals to enzymatic construction of natural products, carbocations are central to molecular synthesis. Although these reactive intermediates are engaged in stereoselective processes in nature, exerting enantiocontrol over carbocations with synthetic catalysts remains challenging. Many resonance-stabilized tricoordinated carbocations, such as iminium and oxocarbenium ions, have been applied in catalytic enantioselective reactions. However, their dicoordinated counterparts (aryl and vinyl carbocations) have not, despite their emerging utility in chemical synthesis. We report the discovery of a highly enantioselective vinyl carbocation carbon-hydrogen (C-H) insertion reaction enabled by imidodiphosphorimidate organocatalysts. Active site confinement featured in this catalyst class not only enables effective enantiocontrol but also expands the scope of vinyl cation C-H insertion chemistry, which broadens the utility of this transition metal-free C(sp3)-H functionalization platform.

Adaptations to infectious and dietary pressures shape mammalian physiology and disease risk. How such adaptations affect sex-biased diseases remains insufficiently studied. In this study, we show that sex-dependent hepatic gene programs confer a robust (~300%) survival advantage for male mice during lethal bacterial infection. The transcription factor B cell lymphoma 6 (BCL6), which masculinizes hepatic gene expression at puberty, is essential for this advantage. However, protection by BCL6 protein comes at a cost during conditions of dietary excess, which result in overt fatty liver and glucose intolerance in males. Deleting hepatic BCL6 reverses these phenotypes but markedly lowers male survival during infection, thus establishing a sex-dependent trade-off between host defense and metabolic systems. Our findings offer strong evidence that some current sex-biased diseases are rooted in ancient evolutionary trade-offs between immunity and metabolism.

We engineered an ultrasensitive reporter of p16INK4a, a biomarker of cellular senescence. Our reporter detected p16INK4a-expressing fibroblasts with certain senescent characteristics that appeared shortly after birth in the basement membrane adjacent to epithelial stem cells in the lung. Furthermore, these p16INK4a+ fibroblasts had enhanced capacity to sense tissue inflammation and respond through their increased secretory capacity to promote epithelial regeneration. In addition, p16INK4a expression was required in fibroblasts to enhance epithelial regeneration. This study highlights a role for p16INK4a+ fibroblasts as tissue-resident sentinels in the stem cell niche that monitor barrier integrity and rapidly respond to inflammation to promote tissue regeneration.

Urbanization transforms environments in ways that alter biological evolution. We examined whether urban environmental change drives parallel evolution by sampling 110,019 white clover plants from 6169 populations in 160 cities globally. Plants were assayed for a Mendelian antiherbivore defense that also affects tolerance to abiotic stressors. Urban-rural gradients were associated with the evolution of clines in defense in 47% of cities throughout the world. Variation in the strength of clines was explained by environmental changes in drought stress and vegetation cover that varied among cities. Sequencing 2074 genomes from 26 cities revealed that the evolution of urban-rural clines was best explained by adaptive evolution, but the degree of parallel adaptation varied among cities. Our results demonstrate that urbanization leads to adaptation at a global scale.

The study of quantum phase transitions that are not clearly associated with broken symmetry is a major effort in condensed matter physics, particularly in regard to the problem of high-temperature superconductivity, for which such transitions are thought to underlie the mechanism of superconductivity itself. Here we argue that the putative quantum critical point in the prototypical unconventional superconductor CeCoIn5 is characterized by the delocalization of electrons in a transition that connects two Fermi surfaces of different volumes, with no apparent broken symmetry. Drawing on established theory of f-electron metals, we discuss an interpretation for such a transition that involves the fractionalization of spin and charge, a model that effectively describes the anomalous transport behavior we measured for the Hall effect.

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.

Genetic admixture is central to primate evolution. We combined 50 years of field observations of immigration and group demography with genomic data from ~9 generations of hybrid baboons to investigate the consequences of admixture in the wild. Despite no obvious fitness costs to hybrids, we found signatures of selection against admixture similar to those described for archaic hominins. These patterns were concentrated near genes where ancestry is strongly associated with gene expression. Our analyses also show that introgression is partially predictable across the genome. This study demonstrates the value of integrating genomic and field data for revealing how "genomic signatures of selection" (e.g., reduced introgression in low-recombination regions) manifest in nature; moreover, it underscores the importance of other primates as living models for human evolution.

We analyzed 131 human brains (44 neurotypical, 19 with Tourette syndrome, 9 with schizophrenia, and 59 with autism) for somatic mutations after whole genome sequencing to a depth of more than 200×. Typically, brains had 20 to 60 detectable single-nucleotide mutations, but ~6% of brains harbored hundreds of somatic mutations. Hypermutability was associated with age and damaging mutations in genes implicated in cancers and, in some brains, reflected in vivo clonal expansions. Somatic duplications, likely arising during development, were found in ~5% of normal and diseased brains, reflecting background mutagenesis. Brains with autism were associated with mutations creating putative transcription factor binding motifs in enhancer-like regions in the developing brain. The top-ranked affected motifs corresponded to MEIS (myeloid ectopic viral integration site) transcription factors, suggesting a potential link between their involvement in gene regulation and autism.

Progenitors of the developing human neocortex reside in the ventricular and outer subventricular zones (VZ and OSVZ, respectively). However, whether cells derived from these niches have similar developmental fates is unknown. By performing fate mapping in primary human tissue, we demonstrate that astrocytes derived from these niches populate anatomically distinct layers. Cortical plate astrocytes emerge from VZ progenitors and proliferate locally, while putative white matter astrocytes are morphologically heterogeneous and emerge from both VZ and OSVZ progenitors. Furthermore, via single-cell sequencing of morphologically defined astrocyte subtypes using Patch-seq, we identify molecular distinctions between VZ-derived cortical plate astrocytes and OSVZ-derived white matter astrocytes that persist into adulthood. Together, our study highlights a complex role for cell lineage in the diversification of human neocortical astrocytes.

Superconductivity and charge density waves (CDWs) are competitive, yet coexisting, orders in cuprate superconductors. To understand their microscopic interdependence, a probe capable of discerning their interaction on its natural length and time scale is necessary. We use ultrafast resonant soft x-ray scattering to track the transient evolution of CDW correlations in YBa2Cu3O6+x after the quench of superconductivity by an infrared laser pulse. We observe a nonthermal response of the CDW order characterized by a near doubling of the correlation length within ≈1 picosecond of the superconducting quench. Our results are consistent with a model in which the interaction between superconductivity and CDWs manifests inhomogeneously through disruption of spatial coherence, with superconductivity playing the dominant role in stabilizing CDW topological defects, such as discommensurations.

In this work, we find that CD8+ T cells expressing inhibitory killer cell immunoglobulin-like receptors (KIRs) are the human equivalent of Ly49+CD8+ regulatory T cells in mice and are increased in the blood and inflamed tissues of patients with a variety of autoimmune diseases. Moreover, these CD8+ T cells efficiently eliminated pathogenic gliadin-specific CD4+ T cells from the leukocytes of celiac disease patients in vitro. We also find elevated levels of KIR+CD8+ T cells, but not CD4+ regulatory T cells, in COVID-19 patients, correlating with disease severity and vasculitis. Selective ablation of Ly49+CD8+ T cells in virus-infected mice led to autoimmunity after infection. Our results indicate that in both species, these regulatory CD8+ T cells act specifically to suppress pathogenic T cells in autoimmune and infectious diseases.

Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease. Knowledge of circulating immune cell types and states associated with SLE remains incomplete. We profiled more than 1.2 million peripheral blood mononuclear cells (162 cases, 99 controls) with multiplexed single-cell RNA sequencing (mux-seq). Cases exhibited elevated expression of type 1 interferon-stimulated genes (ISGs) in monocytes, reduction of naïve CD4+ T cells that correlated with monocyte ISG expression, and expansion of repertoire-restricted cytotoxic GZMH+ CD8+ T cells. Cell type-specific expression features predicted case-control status and stratified patients into two molecular subtypes. We integrated dense genotyping data to map cell type-specific cis-expression quantitative trait loci and to link SLE-associated variants to cell type-specific expression. These results demonstrate mux-seq as a systematic approach to characterize cellular composition, identify transcriptional signatures, and annotate genetic variants associated with SLE.

Mobile elements and repetitive genomic regions are sources of lineage-specific genomic innovation and uniquely fingerprint individual genomes. Comprehensive analyses of such repeat elements, including those found in more complex regions of the genome, require a complete, linear genome assembly. We present a de novo repeat discovery and annotation of the T2T-CHM13 human reference genome. We identified previously unknown satellite arrays, expanded the catalog of variants and families for repeats and mobile elements, characterized classes of complex composite repeats, and located retroelement transduction events. We detected nascent transcription and delineated CpG methylation profiles to define the structure of transcriptionally active retroelements in humans, including those in centromeres. These data expand our insight into the diversity, distribution, and evolution of repetitive regions that have shaped the human genome.

The completion of a telomere-to-telomere human reference genome, T2T-CHM13, has resolved complex regions of the genome, including repetitive and homologous regions. Here, we present a high-resolution epigenetic study of previously unresolved sequences, representing entire acrocentric chromosome short arms, gene family expansions, and a diverse collection of repeat classes. This resource precisely maps CpG methylation (32.28 million CpGs), DNA accessibility, and short-read datasets (166,058 previously unresolved chromatin immunoprecipitation sequencing peaks) to provide evidence of activity across previously unidentified or corrected genes and reveals clinically relevant paralog-specific regulation. Probing CpG methylation across human centromeres from six diverse individuals generated an estimate of variability in kinetochore localization. This analysis provides a framework with which to investigate the most elusive regions of the human genome, granting insights into epigenetic regulation.

Existing human genome assemblies have almost entirely excluded repetitive sequences within and near centromeres, limiting our understanding of their organization, evolution, and functions, which include facilitating proper chromosome segregation. Now, a complete, telomere-to-telomere human genome assembly (T2T-CHM13) has enabled us to comprehensively characterize pericentromeric and centromeric repeats, which constitute 6.2% of the genome (189.9 megabases). Detailed maps of these regions revealed multimegabase structural rearrangements, including in active centromeric repeat arrays. Analysis of centromere-associated sequences uncovered a strong relationship between the position of the centromere and the evolution of the surrounding DNA through layered repeat expansions. Furthermore, comparisons of chromosome X centromeres across a diverse panel of individuals illuminated high degrees of structural, epigenetic, and sequence variation in these complex and rapidly evolving regions.

The COVID-19 pandemic has revealed critical knowledge gaps in our understanding of and a need to update the traditional view of transmission pathways for respiratory viruses. The long-standing definitions of droplet and airborne transmission do not account for the mechanisms by which virus-laden respiratory droplets and aerosols travel through the air and lead to infection. In this Review, we discuss current evidence regarding the transmission of respiratory viruses by aerosols-how they are generated, transported, and deposited, as well as the factors affecting the relative contributions of droplet-spray deposition versus aerosol inhalation as modes of transmission. Improved understanding of aerosol transmission brought about by studies of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection requires a reevaluation of the major transmission pathways for other respiratory viruses, which will allow better-informed controls to reduce airborne transmission.

Total daily energy expenditure ("total expenditure") reflects daily energy needs and is a critical variable in human health and physiology, but its trajectory over the life course is poorly studied. We analyzed a large, diverse database of total expenditure measured by the doubly labeled water method for males and females aged 8 days to 95 years. Total expenditure increased with fat-free mass in a power-law manner, with four distinct life stages. Fat-free mass-adjusted expenditure accelerates rapidly in neonates to ~50% above adult values at ~1 year; declines slowly to adult levels by ~20 years; remains stable in adulthood (20 to 60 years), even during pregnancy; then declines in older adults. These changes shed light on human development and aging and should help shape nutrition and health strategies across the life span.

Yu et al. suggested calculating precisely the size ranges of the three parts of our figure 3A, adjusting the free-energy levels in figure 3B, and considering the shape effect in the first-principles calculation. The first and second suggestions raise strong concerns for misinterpretation and overinterpretation of our experiments. The original calculation is sufficient to support our claim about crystalline-to-disordered transformations.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern evades antibody-mediated immunity that comes from vaccination or infection with earlier variants due to accumulation of numerous spike mutations. To understand the Omicron antigenic shift, we determined cryo-electron microscopy and x-ray crystal structures of the spike protein and the receptor-binding domain bound to the broadly neutralizing sarbecovirus monoclonal antibody (mAb) S309 (the parent mAb of sotrovimab) and to the human ACE2 receptor. We provide a blueprint for understanding the marked reduction of binding of other therapeutic mAbs that leads to dampened neutralizing activity. Remodeling of interactions between the Omicron receptor-binding domain and human ACE2 likely explains the enhanced affinity for the host receptor relative to the ancestral virus.