A Paleogenomic Reconstruction of the Deep Population History of the Andes.

There are many unanswered questions about the population history of the Central and South Central Andes, particularly regarding the impact of large-scale societies, such as the Moche, Wari, Tiwanaku, and Inca. We assembled genome-wide data on 89 individuals dating from ∼9,000-500 years ago (BP), with a particular focus on the period of the rise and fall of state societies. Today's genetic structure began to develop by 5,800 BP, followed by bi-directional gene flow between the North and South Highlands, and between the Highlands and Coast. We detect minimal admixture among neighboring groups between ∼2,000-500 BP, although we do detect cosmopolitanism (people of diverse ancestries living side-by-side) in the heartlands of the Tiwanaku and Inca polities. We also highlight cases of long-range mobility connecting the Andes to Argentina and the Northwest Andes to the Amazon Basin. VIDEO ABSTRACT.

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

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

Global Trends in Marine Plankton Diversity across Kingdoms of Life.

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

Fasting reshapes tissue-specific niches to improve NK cell-mediated anti-tumor immunity.

Fasting is associated with improved outcomes in cancer. Here, we investigated the impact of fasting on natural killer (NK) cell anti-tumor immunity. Cyclic fasting improved immunity against solid and metastatic tumors in an NK cell-dependent manner. During fasting, NK cells underwent redistribution from peripheral tissues to the bone marrow (BM). In humans, fasting also reduced circulating NK cell numbers. NK cells in the spleen of fasted mice were metabolically rewired by elevated concentrations of fatty acids and glucocorticoids, augmenting fatty acid metabolism via increased expression of the enzyme CPT1A, and Cpt1a deletion impaired NK cell survival and function in this setting. In parallel, redistribution of NK cells to the BM during fasting required the trafficking mediators S1PR5 and CXCR4. These cells were primed by an increased pool of interleukin (IL)-12-expressing BM myeloid cells, which improved IFN-γ production. Our findings identify a link between dietary restriction and optimized innate immune responses, with the potential to enhance immunotherapy strategies.

Somatic and Germline Mutation Periodicity Follow the Orientation of the DNA Minor Groove around Nucleosomes.

Mutation rates along the genome are highly variable and influenced by several chromatin features. Here, we addressed how nucleosomes, the most pervasive chromatin structure in eukaryotes, affect the generation of mutations. We discovered that within nucleosomes, the somatic mutation rate across several tumor cohorts exhibits a strong 10 base pair (bp) periodicity. This periodic pattern tracks the alternation of the DNA minor groove facing toward and away from the histones. The strength and phase of the mutation rate periodicity are determined by the mutational processes active in tumors. We uncovered similar periodic patterns in the genetic variation among human and Arabidopsis populations, also detectable in their divergence from close species, indicating that the same principles underlie germline and somatic mutation rates. We propose that differential DNA damage and repair processes dependent on the minor groove orientation in nucleosome-bound DNA contribute to the 10-bp periodicity in AT/CG content in eukaryotic genomes.

Reconstructing the Deep Population History of Central and South America.

We report genome-wide ancient DNA from 49 individuals forming four parallel time transects in Belize, Brazil, the Central Andes, and the Southern Cone, each dating to at least ∼9,000 years ago. The common ancestral population radiated rapidly from just one of the two early branches that contributed to Native Americans today. We document two previously unappreciated streams of gene flow between North and South America. One affected the Central Andes by ∼4,200 years ago, while the other explains an affinity between the oldest North American genome associated with the Clovis culture and the oldest Central and South Americans from Chile, Brazil, and Belize. However, this was not the primary source for later South Americans, as the other ancient individuals derive from lineages without specific affinity to the Clovis-associated genome, suggesting a population replacement that began at least 9,000 years ago and was followed by substantial population continuity in multiple regions.

Modified mRNA Vaccines Protect against Zika Virus Infection.

The emergence of ZIKV infection has prompted a global effort to develop safe and effective vaccines. We engineered a lipid nanoparticle (LNP) encapsulated modified mRNA vaccine encoding wild-type or variant ZIKV structural genes and tested immunogenicity and protection in mice. Two doses of modified mRNA LNPs encoding prM-E genes that produced virus-like particles resulted in high neutralizing antibody titers (∼1/100,000) that protected against ZIKV infection and conferred sterilizing immunity. To offset a theoretical concern of ZIKV vaccines inducing antibodies that cross-react with the related dengue virus (DENV), we designed modified prM-E RNA encoding mutations destroying the conserved fusion-loop epitope in the E protein. This variant protected against ZIKV and diminished production of antibodies enhancing DENV infection in cells or mice. A modified mRNA vaccine can prevent ZIKV disease and be adapted to reduce the risk of sensitizing individuals to subsequent exposure to DENV, should this become a clinically relevant concern.

EGFR Ligands Differentially Stabilize Receptor Dimers to Specify Signaling Kinetics.

Epidermal growth factor receptor (EGFR) regulates many crucial cellular programs, with seven different activating ligands shaping cell signaling in distinct ways. Using crystallography and other approaches, we show how the EGFR ligands epiregulin (EREG) and epigen (EPGN) stabilize different dimeric conformations of the EGFR extracellular region. As a consequence, EREG or EPGN induce less stable EGFR dimers than EGF-making them partial agonists of EGFR dimerization. Unexpectedly, this weakened dimerization elicits more sustained EGFR signaling than seen with EGF, provoking responses in breast cancer cells associated with differentiation rather than proliferation. Our results reveal how responses to different EGFR ligands are defined by receptor dimerization strength and signaling dynamics. These findings have broad implications for understanding receptor tyrosine kinase (RTK) signaling specificity. Our results also suggest parallels between partial and/or biased agonism in RTKs and G-protein-coupled receptors, as well as new therapeutic opportunities for correcting RTK signaling output.

Inflammation Improves Glucose Homeostasis through IKKß-XBP1s Interaction.

It is widely believed that inflammation associated with obesity has an important role in the development of type 2 diabetes. IκB kinase beta (IKKß) is a crucial kinase that responds to inflammatory stimuli such as tumor necrosis factor α (TNF-α) by initiating a variety of intracellular signaling cascades and is considered to be a key element in the inflammation-mediated development of insulin resistance. We show here, contrary to expectation, that IKKß-mediated inflammation is a positive regulator of hepatic glucose homeostasis. IKKß phosphorylates the spliced form of X-Box Binding Protein 1 (XBP1s) and increases the activity of XBP1s. We have used three experimental approaches to enhance the IKKß activity in the liver of obese mice and observed increased XBP1s activity, reduced ER stress, and a significant improvement in insulin sensitivity and consequently in glucose homeostasis. Our results reveal a beneficial role of IKKß-mediated hepatic inflammation in glucose homeostasis.

A Vulnerability of a Subset of Colon Cancers with Potential Clinical Utility.

BRAF(V600E) mutant colon cancers (CCs) have a characteristic gene expression signature that is also found in some tumors lacking this mutation. Collectively, they are referred to as "BRAF-like" tumors and represent some 20% of CCs. We used a shRNA-based genetic screen focused on genes upregulated in BRAF(V600E) CCs to identify vulnerabilities of this tumor subtype that might be exploited therapeutically. Here, we identify RANBP2 (also known as NUP358) as essential for survival of BRAF-like, but not for non-BRAF-like, CC cells. Suppression of RANBP2 results in mitotic defects only in BRAF-like CC cells, leading to cell death. Mechanistically, RANBP2 silencing reduces microtubule outgrowth from the kinetochores, thereby inducing spindle perturbations, providing an explanation for the observed mitotic defects. We find that BRAF-like CCs display far greater sensitivity to the microtubule poison vinorelbine both in vitro and in vivo, suggesting that vinorelbine is a potential tailored treatment for BRAF-like CCs.

Treatment of obesity with celastrol.

Despite all modern advances in medicine, an effective drug treatment of obesity has not been found yet. Discovery of leptin two decades ago created hopes for treatment of obesity. However, development of leptin resistance has been a big obstacle, mitigating a leptin-centric treatment of obesity. Here, by using in silico drug-screening methods, we discovered that Celastrol, a pentacyclic triterpene extracted from the roots of Tripterygium Wilfordi (thunder god vine) plant, is a powerful anti-obesity agent. Celastrol suppresses food intake, blocks reduction of energy expenditure, and leads to up to 45% weight loss in hyperleptinemic diet-induced obese (DIO) mice by increasing leptin sensitivity, but it is ineffective in leptin-deficient (ob/ob) and leptin receptor-deficient (db/db) mouse models. These results indicate that Celastrol is a leptin sensitizer and a promising agent for the pharmacological treatment of obesity.

Inter-cellular forces orchestrate contact inhibition of locomotion.

Contact inhibition of locomotion (CIL) is a multifaceted process that causes many cell types to repel each other upon collision. During development, this seemingly uncoordinated reaction is a critical driver of cellular dispersion within embryonic tissues. Here, we show that Drosophila hemocytes require a precisely orchestrated CIL response for their developmental dispersal. Hemocyte collision and subsequent repulsion involves a stereotyped sequence of kinematic stages that are modulated by global changes in cytoskeletal dynamics. Tracking actin retrograde flow within hemocytes in vivo reveals synchronous reorganization of colliding actin networks through engagement of an inter-cellular adhesion. This inter-cellular actin-clutch leads to a subsequent build-up in lamellar tension, triggering the development of a transient stress fiber, which orchestrates cellular repulsion. Our findings reveal that the physical coupling of the flowing actin networks during CIL acts as a mechanotransducer, allowing cells to haptically sense each other and coordinate their behaviors.

Hybrid speciation driven by multilocus introgression of ecological traits.

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

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.

Recognition and control of neutrophil extracellular trap formation by MICL.

Regulation of neutrophil activation is critical for disease control. Neutrophil extracellular traps (NETs), which are web-like structures composed of DNA and neutrophil-derived proteins, are formed following pro-inflammatory signals; however, if this process is uncontrolled, NETs contribute to disease pathogenesis, exacerbating inflammation and host tissue damage1,2. Here we show that myeloid inhibitory C-type lectin-like (MICL), an inhibitory C-type lectin receptor, directly recognizes DNA in NETs; this interaction is vital to regulate neutrophil activation. Loss or inhibition of MICL functionality leads to uncontrolled NET formation through the ROS-PAD4 pathway and the development of an auto-inflammatory feedback loop. We show that in the context of rheumatoid arthritis, such dysregulation leads to exacerbated pathology in both mouse models and in human patients, where autoantibodies to MICL inhibit key functions of this receptor. Of note, we also detect similarly inhibitory anti-MICL autoantibodies in patients with other diseases linked to aberrant NET formation, including lupus and severe COVID-19. By contrast, dysregulation of NET release is protective during systemic infection with the fungal pathogen Aspergillus fumigatus. Together, we show that the recognition of NETs by MICL represents a fundamental autoregulatory pathway that controls neutrophil activity and NET formation.

TMPRSS2 is a functional receptor for human coronavirus HKU1.

Four endemic seasonal human coronaviruses causing common colds circulate worldwide: HKU1, 229E, NL63 and OC43 (ref. 1). After binding to cellular receptors, coronavirus spike proteins are primed for fusion by transmembrane serine protease 2 (TMPRSS2) or endosomal cathepsins2-9. NL63 uses angiotensin-converting enzyme 2 as a receptor10, whereas 229E uses human aminopeptidase-N11. HKU1 and OC43 spikes bind cells through 9-O-acetylated sialic acid, but their protein receptors remain unknown12. Here we show that TMPRSS2 is a functional receptor for HKU1. TMPRSS2 triggers HKU1 spike-mediated cell-cell fusion and pseudovirus infection. Catalytically inactive TMPRSS2 mutants do not cleave HKU1 spike but allow pseudovirus infection. Furthermore, TMPRSS2 binds with high affinity to the HKU1 receptor binding domain (Kd 334 and 137 nM for HKU1A and HKU1B genotypes) but not to SARS-CoV-2. Conserved amino acids in the HKU1 receptor binding domain are essential for binding to TMPRSS2 and pseudovirus infection. Newly designed anti-TMPRSS2 nanobodies potently inhibit HKU1 spike attachment to TMPRSS2, fusion and pseudovirus infection. The nanobodies also reduce infection of primary human bronchial cells by an authentic HKU1 virus. Our findings illustrate the various evolution strategies of coronaviruses, which use TMPRSS2 to either directly bind to target cells or prime their spike for membrane fusion and entry.

A transcriptomic taxonomy of mouse brain-wide spinal projecting neurons.

The brain controls nearly all bodily functions via spinal projecting neurons (SPNs) that carry command signals from the brain to the spinal cord. However, a comprehensive molecular characterization of brain-wide SPNs is still lacking. Here we transcriptionally profiled a total of 65,002 SPNs, identified 76 region-specific SPN types, and mapped these types into a companion atlas of the whole mouse brain1. This taxonomy reveals a three-component organization of SPNs: (1) molecularly homogeneous excitatory SPNs from the cortex, red nucleus and cerebellum with somatotopic spinal terminations suitable for point-to-point communication; (2) heterogeneous populations in the reticular formation with broad spinal termination patterns, suitable for relaying commands related to the activities of the entire spinal cord; and (3) modulatory neurons expressing slow-acting neurotransmitters and/or neuropeptides in the hypothalamus, midbrain and reticular formation for 'gain setting' of brain-spinal signals. In addition, this atlas revealed a LIM homeobox transcription factor code that parcellates the reticulospinal neurons into five molecularly distinct and spatially segregated populations. Finally, we found transcriptional signatures of a subset of SPNs with large soma size and correlated these with fast-firing electrophysiological properties. Together, this study establishes a comprehensive taxonomy of brain-wide SPNs and provides insight into the functional organization of SPNs in mediating brain control of bodily functions.

Structural insights into the mechanism of the sodium/iodide symporter.

The sodium/iodide symporter (NIS) is the essential plasma membrane protein that mediates active iodide (I-) transport into the thyroid gland, the first step in the biosynthesis of the thyroid hormones-the master regulators of intermediary metabolism. NIS couples the inward translocation of I- against its electrochemical gradient to the inward transport of Na+ down its electrochemical gradient1,2. For nearly 50 years before its molecular identification3, NIS was the molecule at the centre of the single most effective internal radiation cancer therapy: radioiodide (131I-) treatment for thyroid cancer2. Mutations in NIS cause congenital hypothyroidism, which must be treated immediately after birth to prevent stunted growth and cognitive deficiency2. Here we report three structures of rat NIS, determined by single-particle cryo-electron microscopy: one with no substrates bound; one with two Na+ and one I- bound; and one with one Na+ and the oxyanion perrhenate bound. Structural analyses, functional characterization and computational studies show the substrate-binding sites and key residues for transport activity. Our results yield insights into how NIS selects, couples and translocates anions-thereby establishing a framework for understanding NIS function-and how it transports different substrates with different stoichiometries and releases substrates from its substrate-binding cavity into the cytosol.

Biosynthetic potential of the global ocean microbiome.

Natural microbial communities are phylogenetically and metabolically diverse. In addition to underexplored organismal groups1, this diversity encompasses a rich discovery potential for ecologically and biotechnologically relevant enzymes and biochemical compounds2,3. However, studying this diversity to identify genomic pathways for the synthesis of such compounds4 and assigning them to their respective hosts remains challenging. The biosynthetic potential of microorganisms in the open ocean remains largely uncharted owing to limitations in the analysis of genome-resolved data at the global scale. Here we investigated the diversity and novelty of biosynthetic gene clusters in the ocean by integrating around 10,000 microbial genomes from cultivated and single cells with more than 25,000 newly reconstructed draft genomes from more than 1,000 seawater samples. These efforts revealed approximately 40,000 putative mostly new biosynthetic gene clusters, several of which were found in previously unsuspected phylogenetic groups. Among these groups, we identified a lineage rich in biosynthetic gene clusters ('Candidatus Eudoremicrobiaceae') that belongs to an uncultivated bacterial phylum and includes some of the most biosynthetically diverse microorganisms in this environment. From these, we characterized the phospeptin and pythonamide pathways, revealing cases of unusual bioactive compound structure and enzymology, respectively. Together, this research demonstrates how microbiomics-driven strategies can enable the investigation of previously undescribed enzymes and natural products in underexplored microbial groups and environments.

IκB kinase-α coordinates BRD4 and JAK/STAT signaling to subvert DNA damage-based anticancer therapy.

Activation of the IκB kinase (IKK) complex has recurrently been linked to colorectal cancer (CRC) initiation and progression. However, identification of downstream effectors other than NF-κB has remained elusive. Here, analysis of IKK-dependent substrates in CRC cells after UV treatment revealed that phosphorylation of BRD4 by IKK-α is required for its chromatin-binding at target genes upon DNA damage. Moreover, IKK-α induces the NF-κB-dependent transcription of the cytokine LIF, leading to STAT3 activation, association with BRD4 and recruitment to specific target genes. IKK-α abrogation results in defective BRD4 and STAT3 functions and consequently irreparable DNA damage and apoptotic cell death upon different stimuli. Simultaneous inhibition of BRAF-dependent IKK-α activity, BRD4, and the JAK/STAT pathway enhanced the therapeutic potential of 5-fluorouracil combined with irinotecan in CRC cells and is curative in a chemotherapy-resistant xenograft model. Finally, coordinated expression of LIF and IKK-α is a poor prognosis marker for CRC patients. Our data uncover a functional link between IKK-α, BRD4, and JAK/STAT signaling with clinical relevance.