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

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

Colorectal Cancer Cells Enter a Diapause-like DTP State to Survive Chemotherapy.

Cancer cells enter a reversible drug-tolerant persister (DTP) state to evade death from chemotherapy and targeted agents. It is increasingly appreciated that DTPs are important drivers of therapy failure and tumor relapse. We combined cellular barcoding and mathematical modeling in patient-derived colorectal cancer models to identify and characterize DTPs in response to chemotherapy. Barcode analysis revealed no loss of clonal complexity of tumors that entered the DTP state and recurred following treatment cessation. Our data fit a mathematical model where all cancer cells, and not a small subpopulation, possess an equipotent capacity to become DTPs. Mechanistically, we determined that DTPs display remarkable transcriptional and functional similarities to diapause, a reversible state of suspended embryonic development triggered by unfavorable environmental conditions. Our study provides insight into how cancer cells use a developmentally conserved mechanism to drive the DTP state, pointing to novel therapeutic opportunities to target DTPs.

Key Parameters of Tumor Epitope Immunogenicity Revealed Through a Consortium Approach Improve Neoantigen Prediction.

Many approaches to identify therapeutically relevant neoantigens couple tumor sequencing with bioinformatic algorithms and inferred rules of tumor epitope immunogenicity. However, there are no reference data to compare these approaches, and the parameters governing tumor epitope immunogenicity remain unclear. Here, we assembled a global consortium wherein each participant predicted immunogenic epitopes from shared tumor sequencing data. 608 epitopes were subsequently assessed for T cell binding in patient-matched samples. By integrating peptide features associated with presentation and recognition, we developed a model of tumor epitope immunogenicity that filtered out 98% of non-immunogenic peptides with a precision above 0.70. Pipelines prioritizing model features had superior performance, and pipeline alterations leveraging them improved prediction performance. These findings were validated in an independent cohort of 310 epitopes prioritized from tumor sequencing data and assessed for T cell binding. This data resource enables identification of parameters underlying effective anti-tumor immunity and is available to the research community.

Structural transitions enable interleukin-18 maturation and signaling.

Several interleukin-1 (IL-1) family members, including IL-1ß and IL-18, require processing by inflammasome-associated caspases to unleash their activities. Here, we unveil, by cryoelectron microscopy (cryo-EM), two major conformations of the complex between caspase-1 and pro-IL-18. One conformation is similar to the complex of caspase-4 and pro-IL-18, with interactions at both the active site and an exosite (closed conformation), and the other only contains interactions at the active site (open conformation). Thus, pro-IL-18 recruitment and processing by caspase-1 is less dependent on the exosite than the active site, unlike caspase-4. Structure determination by nuclear magnetic resonance uncovers a compact fold of apo pro-IL-18, which is similar to caspase-1-bound pro-IL-18 but distinct from cleaved IL-18. Binding sites for IL-18 receptor and IL-18 binding protein are only formed upon conformational changes after pro-IL-18 cleavage. These studies show how pro-IL-18 is selected as a caspase-1 substrate, and why cleavage is necessary for its inflammatory activity.

High-Dimensional Analysis Delineates Myeloid and Lymphoid Compartment Remodeling during Successful Immune-Checkpoint Cancer Therapy.

Although current immune-checkpoint therapy (ICT) mainly targets lymphoid cells, it is associated with a broader remodeling of the tumor micro-environment. Here, using complementary forms of high-dimensional profiling, we define differences across all hematopoietic cells from syngeneic mouse tumors during unrestrained tumor growth or effective ICT. Unbiased assessment of gene expression of tumor-infiltrating cells by single-cell RNA sequencing (scRNAseq) and longitudinal assessment of cellular protein expression by mass cytometry (CyTOF) revealed significant remodeling of both the lymphoid and myeloid intratumoral compartments. Surprisingly, we observed multiple subpopulations of monocytes/macrophages, distinguishable by the markers CD206, CX3CR1, CD1d, and iNOS, that change over time during ICT in a manner partially dependent on IFNγ. Our data support the hypothesis that this macrophage polarization/activation results from effects on circulatory monocytes and early macrophages entering tumors, rather than on pre-polarized mature intratumoral macrophages.

A Periplasmic Polymer Curves Vibrio cholerae and Promotes Pathogenesis.

Pathogenic Vibrio cholerae remains a major human health concern. V. cholerae has a characteristic curved rod morphology, with a longer outer face and a shorter inner face. The mechanism and function of this curvature were previously unknown. Here, we identify and characterize CrvA, the first curvature determinant in V. cholerae. CrvA self-assembles into filaments at the inner face of cell curvature. Unlike traditional cytoskeletons, CrvA localizes to the periplasm and thus can be considered a periskeletal element. To quantify how curvature forms, we developed QuASAR (quantitative analysis of sacculus architecture remodeling), which measures subcellular peptidoglycan dynamics. QuASAR reveals that CrvA asymmetrically patterns peptidoglycan insertion rather than removal, causing more material insertions into the outer face than the inner face. Furthermore, crvA is quorum regulated, and CrvA-dependent curvature increases at high cell density. Finally, we demonstrate that CrvA promotes motility in hydrogels and confers an advantage in host colonization and pathogenesis.

Neoantigen-specific cytotoxic Tr1 CD4 T cells suppress cancer immunotherapy.

CD4+ T cells can either enhance or inhibit tumour immunity. Although regulatory T cells have long been known to impede antitumour responses1-5, other CD4+ T cells have recently been implicated in inhibiting this response6,7. Yet, the nature and function of the latter remain unclear. Here, using vaccines containing MHC class I (MHC-I) neoantigens (neoAgs) and different doses of tumour-derived MHC-II neoAgs, we discovered that whereas the inclusion of vaccines with low doses of MHC-II-restricted peptides (LDVax) promoted tumour rejection, vaccines containing high doses of the same MHC-II neoAgs (HDVax) inhibited rejection. Characterization of the inhibitory cells induced by HDVax identified them as type 1 regulatory T (Tr1) cells expressing IL-10, granzyme B, perforin, CCL5 and LILRB4. Tumour-specific Tr1 cells suppressed tumour rejection induced by anti-PD1, LDVax or adoptively transferred tumour-specific effector T cells. Mechanistically, HDVax-induced Tr1 cells selectively killed MHC-II tumour antigen-presenting type 1 conventional dendritic cells (cDC1s), leading to low numbers of cDC1s in tumours. We then documented modalities to overcome this inhibition, specifically via anti-LILRB4 blockade, using a CD8-directed IL-2 mutein, or targeted loss of cDC2/monocytes. Collectively, these data show that cytotoxic Tr1 cells, which maintain peripheral tolerance, also inhibit antitumour responses and thereby function to impede immune control of cancer.

Bipolar impact and phasing of Heinrich-type climate variability.

During the last ice age, the Laurentide Ice Sheet exhibited extreme iceberg discharge events that are recorded in North Atlantic sediments1. These Heinrich events have far-reaching climate impacts, including widespread disruptions to hydrological and biogeochemical cycles2-4. They occurred during Heinrich stadials-cold periods with strongly weakened Atlantic overturning circulation5-7. Heinrich-type variability is not distinctive in Greenland water isotope ratios, a well-dated site temperature proxy8, complicating efforts to assess their regional climate impact and phasing against Antarctic climate change. Here we show that Heinrich events have no detectable temperature impact on Greenland and cooling occurs at the onset of several Heinrich stadials, and that both types of Heinrich variability have a distinct imprint on Antarctic climate. Antarctic ice cores show accelerated warming that is synchronous with increases in methane during Heinrich events, suggesting an atmospheric teleconnection9, despite the absence of a Greenland climate signal. Greenland ice-core nitrogen stable isotope ratios, a sensitive temperature proxy, indicate an abrupt cooling of about three degrees Celsius at the onset of Heinrich Stadial 1 (17.8 thousand years before present, where present is defined as 1950). Antarctic warming lags this cooling by 133 ± 93 years, consistent with an oceanic teleconnection. Paradoxically, proximal sites are less affected by Heinrich events than remote sites, suggesting spatially complex event dynamics.

Widespread six degrees Celsius cooling on land during the Last Glacial Maximum.

The magnitude of global cooling during the Last Glacial Maximum (LGM, the coldest multimillennial interval of the last glacial period) is an important constraint for evaluating estimates of Earth's climate sensitivity1,2. Reliable LGM temperatures come from high-latitude ice cores3,4, but substantial disagreement exists between proxy records in the low latitudes1,5-8, where quantitative low-elevation records on land are scarce. Filling this data gap, noble gases in ancient groundwater record past land surface temperatures through a direct physical relationship that is rooted in their temperature-dependent solubility in water9,10. Dissolved noble gases are suitable tracers of LGM temperature because of their complete insensitivity to biological and chemical processes and the ubiquity of LGM-aged groundwater around the globe11,12. However, although several individual noble gas studies have found substantial tropical LGM cooling13-16, they have used different methodologies and provide limited spatial coverage. Here we use noble gases in groundwater to show that the low-altitude, low-to-mid-latitude land surface (45 degrees south to 35 degrees north) cooled by 5.8 ± 0.6 degrees Celsius (mean ± 95% confidence interval) during the LGM. Our analysis includes four decades of groundwater noble gas data from six continents, along with new records from the tropics, all of which were interpreted using the same physical framework. Our land-based result broadly supports a recent reconstruction based on marine proxy data assimilation1 that suggested greater climate sensitivity than previous estimates5-7.

Genome copy number predicts extreme evolutionary rate variation in plant mitochondrial DNA.

Nuclear and organellar genomes can evolve at vastly different rates despite occupying the same cell. In most bilaterian animals, mitochondrial DNA (mtDNA) evolves faster than nuclear DNA, whereas this trend is generally reversed in plants. However, in some exceptional angiosperm clades, mtDNA substitution rates have increased up to 5,000-fold compared with closely related lineages. The mechanisms responsible for this acceleration are generally unknown. Because plants rely on homologous recombination to repair mtDNA damage, we hypothesized that mtDNA copy numbers may predict evolutionary rates, as lower copy numbers may provide fewer templates for such repair mechanisms. In support of this hypothesis, we found that copy number explains 47% of the variation in synonymous substitution rates of mtDNA across 60 diverse seed plant species representing ~300 million years of evolution. Copy number was also negatively correlated with mitogenome size, which may be a cause or consequence of mutation rate variation. Both relationships were unique to mtDNA and not observed in plastid DNA. These results suggest that homologous recombinational repair plays a role in driving mtDNA substitution rates in plants and may explain variation in mtDNA evolution more broadly across eukaryotes. Our findings also contribute to broader questions about the relationships between mutation rates, genome size, selection efficiency, and the drift-barrier hypothesis.

Estimates of actual and potential lives saved in the United States from the use of COVID-19 convalescent plasma.

In the Spring of 2020, the United States of America (USA) deployed COVID-19 convalescent plasma (CCP) to treat hospitalized patients. Over 500,000 patients were treated with CCP during the first year of the pandemic. In this study, we estimated the number of actual inpatient lives saved by CCP treatment in the United States of America based on CCP weekly use, weekly national mortality data, and CCP mortality reduction data from meta-analyses of randomized controlled trials and real-world data. We also estimate the potential number of lives saved if CCP had been deployed for 100% of hospitalized patients or used in 15 to 75% of outpatients. Depending on the assumptions modeled in stratified analyses, we estimated that CCP saved between 16,476 and 66,296 lives. The CCP ideal use might have saved as many as 234,869 lives and prevented 1,136,133 hospitalizations. CCP deployment was a successful strategy for ameliorating the impact of the COVID-19 pandemic in the USA. This experience has important implications for convalescent plasma use in future infectious disease emergencies.

A randomized trial of behavioral interventions yielding sustained reductions in distracted driving.

Distracted driving is responsible for nearly 1 million crashes each year in the United States alone, and a major source of driver distraction is handheld phone use. We conducted a randomized, controlled trial to compare the effectiveness of interventions designed to create sustained reductions in handheld use while driving (NCT04587609). Participants were 1,653 consenting Progressive® Snapshot® usage-based auto insurance customers ages 18 to 77 who averaged at least 2 min/h of handheld use while driving in the month prior to study invitation. They were randomly assigned to one of five arms for a 10-wk intervention period. Arm 1 (control) got education about the risks of handheld phone use, as did the other arms. Arm 2 got a free phone mount to facilitate hands-free use. Arm 3 got the mount plus a commitment exercise and tips for hands-free use. Arm 4 got the mount, commitment, and tips plus weekly goal gamification and social competition. Arm 5 was the same as Arm 4, plus offered behaviorally designed financial incentives. Postintervention, participants were monitored until the end of their insurance rating period, 25 to 65 d more. Outcome differences were measured using fractional logistic regression. Arm 4 participants, who received gamification and competition, reduced their handheld use by 20.5% relative to control (P < 0.001); Arm 5 participants, who additionally received financial incentives, reduced their use by 27.6% (P < 0.001). Both groups sustained these reductions through the end of their insurance rating period.

cDC1 prime and are licensed by CD4+ T cells to induce anti-tumour immunity.

Conventional type 1 dendritic cells (cDC1)1 are thought to perform antigen cross-presentation, which is required to prime CD8+ T cells2,3, whereas cDC2 are specialized for priming CD4+ T cells4,5. CD4+ T cells are also considered to help CD8+ T cell responses through a variety of mechanisms6-11, including a process whereby CD4+ T cells 'license' cDC1 for CD8+ T cell priming12. However, this model has not been directly tested in vivo or in the setting of help-dependent tumour rejection. Here we generated an Xcr1Cre mouse strain to evaluate the cellular interactions that mediate tumour rejection in a model requiring CD4+ and CD8+ T cells. As expected, tumour rejection required cDC1 and CD8+ T cell priming required the expression of major histocompatibility class I molecules by cDC1. Unexpectedly, early priming of CD4+ T cells against tumour-derived antigens also required cDC1, and this was not simply because they transport antigens to lymph nodes for processing by cDC2, as selective deletion of major histocompatibility class II molecules in cDC1 also prevented early CD4+ T cell priming. Furthermore, deletion of either major histocompatibility class II or CD40 in cDC1 impaired tumour rejection, consistent with a role for cognate CD4+ T cell interactions and CD40 signalling in cDC1 licensing. Finally, CD40 signalling in cDC1 was critical not only for CD8+ T cell priming, but also for initial CD4+ T cell activation. Thus, in the setting of tumour-derived antigens, cDC1 function as an autonomous platform capable of antigen processing and priming for both CD4+ and CD8+ T cells and of the direct orchestration of their cross-talk that is required for optimal anti-tumour immunity.

Mixed, nonclassical behavior in a classic allosteric protein.

Allostery is a major driver of biological processes requiring coordination. Thus, it is one of the most fundamental and remarkable phenomena in nature, and there is motivation to understand and manipulate it to a multitude of ends. Today, it is often described in terms of two phenomenological models proposed more than a half-century ago involving only T(tense) or R(relaxed) conformations. Here, methyl-based NMR provides extensive detail on a dynamic T to R switch in the classical dimeric allosteric protein, yeast chorismate mutase (CM), that occurs in the absence of substrate, but only with the activator bound. Switching of individual subunits is uncoupled based on direct observation of mixed TR states in the dimer. This unique finding excludes both classic models and solves the paradox of a coexisting hyperbolic binding curve and highly skewed substrate-free T-R equilibrium. Surprisingly, structures of the activator-bound and effector-free forms of CM appear the same by NMR, providing another example of the need to account for dynamic ensembles. The apo enzyme, which has a sigmoidal activity profile, is shown to switch, not to R, but to a related high-energy state. Thus, the conformational repertoire of CM does not just change as a matter of degree depending on the allosteric input, be it effector and/or substrate. Rather, the allosteric model appears to completely change in different contexts, which is only consistent with modern ensemble-based frameworks.

Twist response of actin filaments.

Actin cytoskeleton force generation, sensing, and adaptation are dictated by the bending and twisting mechanics of filaments. Here, we use magnetic tweezers and microfluidics to twist and pull individual actin filaments and evaluate their response to applied loads. Twisted filaments bend and dissipate torsional strain by adopting a supercoiled plectoneme. Pulling prevents plectoneme formation, which causes twisted filaments to sever. Analysis over a range of twisting and pulling forces and direct visualization of filament and single subunit twisting fluctuations yield an actin filament torsional persistence length of ~10 µm, similar to the bending persistence length. Filament severing by cofilin is driven by local twist strain at boundaries between bare and decorated segments and is accelerated by low pN pulling forces. This work explains how contractile forces generated by myosin motors accelerate filament severing by cofilin and establishes a role for filament twisting in the regulation of actin filament stability and assembly dynamics.

Lineage-specific genes are clustered with HET-domain genes and respond to environmental and genetic manipulations regulating reproduction in Neurospora.

Lineage-specific genes (LSGs) have long been postulated to play roles in the establishment of genetic barriers to intercrossing and speciation. In the genome of Neurospora crassa, most of the 670 Neurospora LSGs that are aggregated adjacent to the telomeres are clustered with 61% of the HET-domain genes, some of which regulate self-recognition and define vegetative incompatibility groups. In contrast, the LSG-encoding proteins possess few to no domains that would help to identify potential functional roles. Possible functional roles of LSGs were further assessed by performing transcriptomic profiling in genetic mutants and in response to environmental alterations, as well as examining gene knockouts for phenotypes. Among the 342 LSGs that are dynamically expressed during both asexual and sexual phases, 64% were detectable on unusual carbon sources such as furfural, a wildfire-produced chemical that is a strong inducer of sexual development, and the structurally-related furan 5-hydroxymethyl furfural (HMF). Expression of a significant portion of the LSGs was sensitive to light and temperature, factors that also regulate the switch from asexual to sexual reproduction. Furthermore, expression of the LSGs was significantly affected in the knockouts of adv-1 and pp-1 that regulate hyphal communication, and expression of more than one quarter of the LSGs was affected by perturbation of the mating locus. These observations encouraged further investigation of the roles of clustered lineage-specific and HET-domain genes in ecology and reproduction regulation in Neurospora, especially the regulation of the switch from the asexual growth to sexual reproduction, in response to dramatic environmental conditions changes.

The effect of climate change on avian offspring production: A global meta-analysis.

Climate change affects timing of reproduction in many bird species, but few studies have investigated its influence on annual reproductive output. Here, we assess changes in the annual production of young by female breeders in 201 populations of 104 bird species (N = 745,962 clutches) covering all continents between 1970 and 2019. Overall, average offspring production has declined in recent decades, but considerable differences were found among species and populations. A total of 56.7% of populations showed a declining trend in offspring production (significant in 17.4%), whereas 43.3% exhibited an increase (significant in 10.4%). The results show that climatic changes affect offspring production through compounded effects on ecological and life history traits of species. Migratory and larger-bodied species experienced reduced offspring production with increasing temperatures during the chick-rearing period, whereas smaller-bodied, sedentary species tended to produce more offspring. Likewise, multi-brooded species showed increased breeding success with increasing temperatures, whereas rising temperatures were unrelated to reproductive success in single-brooded species. Our study suggests that rapid declines in size of bird populations reported by many studies from different parts of the world are driven only to a small degree by changes in the production of young.

Uropathogenic Escherichia coli wield enterobactin-derived catabolites as siderophores.

Uropathogenic Escherichia coli (UPEC) secrete multiple siderophore types to scavenge extracellular iron(III) ions during clinical urinary tract infections, despite the metabolic costs of biosynthesis. Here, we find the siderophore enterobactin (Ent) and its related products to be prominent components of the iron-responsive extracellular metabolome of a model UPEC strain. Using defined Ent biosynthesis and import mutants, we identify lower molecular weight dimeric exometabolites as products of incomplete siderophore catabolism, rather than prematurely released biosynthetic intermediates. In E. coli, iron acquisition from iron(III)-Ent complexes requires intracellular esterases that hydrolyze the siderophore. Although UPEC are equipped to consume the products of completely hydrolyzed Ent, we find that Ent and its derivatives may be incompletely hydrolyzed to yield products with retained siderophore activity. These results are consistent with catabolic inefficiency as means to obtain more than one iron ion per siderophore molecule. This is compatible with an evolved UPEC strategy to maximize the nutritional returns from metabolic investments in siderophore biosynthesis.

Genome-wide association identifies a BAHD acyltransferase activity that assembles an ester of glucuronosylglycerol and phenylacetic acid.

Genome-wide association studies (GWAS) are an effective approach to identify new specialized metabolites and the genes involved in their biosynthesis and regulation. In this study, GWAS of Arabidopsis thaliana soluble leaf and stem metabolites identified alleles of an uncharacterized BAHD-family acyltransferase (AT5G57840) associated with natural variation in three structurally related metabolites. These metabolites were esters of glucuronosylglycerol, with one metabolite containing phenylacetic acid as the acyl component of the ester. Knockout and overexpression of AT5G57840 in Arabidopsis and heterologous overexpression in Nicotiana benthamiana and Escherichia coli demonstrated that it is capable of utilizing phenylacetyl-CoA as an acyl donor and glucuronosylglycerol as an acyl acceptor. We, thus, named the protein Glucuronosylglycerol Ester Synthase (GGES). Additionally, phenylacetyl glucuronosylglycerol increased in Arabidopsis CYP79A2 mutants that overproduce phenylacetic acid and was lost in knockout mutants of UDP-sulfoquinovosyl: diacylglycerol sulfoquinovosyl transferase, an enzyme required for glucuronosylglycerol biosynthesis and associated with glycerolipid metabolism under phosphate-starvation stress. GGES is a member of a well-supported clade of BAHD family acyltransferases that arose by duplication and neofunctionalized during the evolution of the Brassicales within a larger clade that includes HCT as well as enzymes that synthesize other plant-specialized metabolites. Together, this work extends our understanding of the catalytic diversity of BAHD acyltransferases and uncovers a pathway that involves contributions from both phenylalanine and lipid metabolism.