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.

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.

HLA-II immunopeptidome profiling and deep learning reveal features of antigenicity to inform antigen discovery.

CD4+ T cell responses are exquisitely antigen specific and directed toward peptide epitopes displayed by human leukocyte antigen class II (HLA-II) on antigen-presenting cells. Underrepresentation of diverse alleles in ligand databases and an incomplete understanding of factors affecting antigen presentation in vivo have limited progress in defining principles of peptide immunogenicity. Here, we employed monoallelic immunopeptidomics to identify 358,024 HLA-II binders, with a particular focus on HLA-DQ and HLA-DP. We uncovered peptide-binding patterns across a spectrum of binding affinities and enrichment of structural antigen features. These aspects underpinned the development of context-aware predictor of T cell antigens (CAPTAn), a deep learning model that predicts peptide antigens based on their affinity to HLA-II and full sequence of their source proteins. CAPTAn was instrumental in discovering prevalent T cell epitopes from bacteria in the human microbiome and a pan-variant epitope from SARS-CoV-2. Together CAPTAn and associated datasets present a resource for antigen discovery and the unraveling genetic associations of HLA alleles with immunopathologies.

The CD4+ T cell response to a commensal-derived epitope transitions from a tolerant to an inflammatory state in Crohn's disease.

Reciprocal interactions between host T helper cells and gut microbiota enforce local immunological tolerance and modulate extra-intestinal immunity. However, our understanding of antigen-specific tolerance to the microbiome is limited. Here, we developed a systematic approach to predict HLA class-II-specific epitopes using the humanized bacteria-originated T cell antigen (hBOTA) algorithm. We identified a diverse set of microbiome epitopes spanning all major taxa that are compatible with presentation by multiple HLA-II alleles. In particular, we uncovered an immunodominant epitope from the TonB-dependent receptor SusC that was universally recognized and ubiquitous among Bacteroidales. In healthy human subjects, SusC-reactive T cell responses were characterized by IL-10-dominant cytokine profiles, whereas in patients with active Crohn's disease, responses were associated with elevated IL-17A. Our results highlight the potential of targeted antigen discovery within the microbiome to reveal principles of tolerance and functional transitions during inflammation.

Cell-Type-Specific Control of Brainstem Locomotor Circuits by Basal Ganglia.

The basal ganglia (BG) are critical for adaptive motor control, but the circuit principles underlying their pathway-specific modulation of target regions are not well understood. Here, we dissect the mechanisms underlying BG direct and indirect pathway-mediated control of the mesencephalic locomotor region (MLR), a brainstem target of BG that is critical for locomotion. We optogenetically dissect the locomotor function of the three neurochemically distinct cell types within the MLR: glutamatergic, GABAergic, and cholinergic neurons. We find that the glutamatergic subpopulation encodes locomotor state and speed, is necessary and sufficient for locomotion, and is selectively innervated by BG. We further show activation and suppression, respectively, of MLR glutamatergic neurons by direct and indirect pathways, which is required for bidirectional control of locomotion by BG circuits. These findings provide a fundamental understanding of how BG can initiate or suppress a motor program through cell-type-specific regulation of neurons linked to specific actions.

Reconstruction and Simulation of Neocortical Microcircuitry.

We present a first-draft digital reconstruction of the microcircuitry of somatosensory cortex of juvenile rat. The reconstruction uses cellular and synaptic organizing principles to algorithmically reconstruct detailed anatomy and physiology from sparse experimental data. An objective anatomical method defines a neocortical volume of 0.29 ± 0.01 mm(3) containing ~31,000 neurons, and patch-clamp studies identify 55 layer-specific morphological and 207 morpho-electrical neuron subtypes. When digitally reconstructed neurons are positioned in the volume and synapse formation is restricted to biological bouton densities and numbers of synapses per connection, their overlapping arbors form ~8 million connections with ~37 million synapses. Simulations reproduce an array of in vitro and in vivo experiments without parameter tuning. Additionally, we find a spectrum of network states with a sharp transition from synchronous to asynchronous activity, modulated by physiological mechanisms. The spectrum of network states, dynamically reconfigured around this transition, supports diverse information processing strategies. PAPERCLIP: VIDEO ABSTRACT.

Major southern San Andreas earthquakes modulated by lake-filling events.

Hydrologic loads can stimulate seismicity in the Earth's crust1. However, evidence for the triggering of large earthquakes remains elusive. The southern San Andreas Fault (SSAF) in Southern California lies next to the Salton Sea2, a remnant of ancient Lake Cahuilla that periodically filled and desiccated over the past millennium3-5. Here we use new geologic and palaeoseismic data to demonstrate that the past six major earthquakes on the SSAF probably occurred during highstands of Lake Cahuilla5,6. To investigate possible causal relationships, we computed time-dependent Coulomb stress changes7,8 due to variations in the lake level. Using a fully coupled model of a poroelastic crust9-11 overlying a viscoelastic mantle12,13, we find that hydrologic loads increased Coulomb stress on the SSAF by several hundred kilopascals and fault-stressing rates by more than a factor of 2, which is probably sufficient for earthquake triggering7,8. The destabilizing effects of lake inundation are enhanced by a nonvertical fault dip14-17, the presence of a fault damage zone18,19 and lateral pore-pressure diffusion20,21. Our model may be applicable to other regions in which hydrologic loading, either natural8,22 or anthropogenic1,23, was associated with substantial seismicity.

Genetic continuity and change among the Indigenous peoples of California.

Before the colonial period, California harboured more language variation than all of Europe, and linguistic and archaeological analyses have led to many hypotheses to explain this diversity1. We report genome-wide data from 79 ancient individuals from California and 40 ancient individuals from Northern Mexico dating to 7,400-200 years before present (BP). Our analyses document long-term genetic continuity between people living on the Northern Channel Islands of California and the adjacent Santa Barbara mainland coast from 7,400 years BP to modern Chumash groups represented by individuals who lived around 200 years BP. The distinctive genetic lineages that characterize present-day and ancient people from Northwest Mexico increased in frequency in Southern and Central California by 5,200 years BP, providing evidence for northward migrations that are candidates for spreading Uto-Aztecan languages before the dispersal of maize agriculture from Mexico2-4. Individuals from Baja California share more alleles with the earliest individual from Central California in the dataset than with later individuals from Central California, potentially reflecting an earlier linguistic substrate, whose impact on local ancestry was diluted by later migrations from inland regions1,5. After 1,600 years BP, ancient individuals from the Channel Islands lived in communities with effective sizes similar to those in pre-agricultural Caribbean and Patagonia, and smaller than those on the California mainland and in sampled regions of Mexico.

Pathogenic SYNGAP1 mutations impair cognitive development by disrupting maturation of dendritic spine synapses.

Mutations that cause intellectual disability (ID) and autism spectrum disorder (ASD) are commonly found in genes that encode for synaptic proteins. However, it remains unclear how mutations that disrupt synapse function impact intellectual ability. In the SYNGAP1 mouse model of ID/ASD, we found that dendritic spine synapses develop prematurely during the early postnatal period. Premature spine maturation dramatically enhanced excitability in the developing hippocampus, which corresponded with the emergence of behavioral abnormalities. Inducing SYNGAP1 mutations after critical developmental windows closed had minimal impact on spine synapse function, whereas repairing these pathogenic mutations in adulthood did not improve behavior and cognition. These data demonstrate that SynGAP protein acts as a critical developmental repressor of neural excitability that promotes the development of life-long cognitive abilities. We propose that the pace of dendritic spine synapse maturation in early life is a critical determinant of normal intellectual development.

A Novel Macrophage Subpopulation Conveys Increased Genetic Risk of Coronary Artery Disease.

BACKGROUND: Coronary artery disease (CAD), the leading cause of death worldwide, is influenced by both environmental and genetic factors. Although over 250 genetic risk loci have been identified through genome-wide association studies, the specific causal variants and their regulatory mechanisms are still largely unknown, particularly in disease-relevant cell types such as macrophages. METHODS: We utilized single-cell RNA-seq and single-cell multiomics approaches in primary human monocyte-derived macrophages to explore the transcriptional regulatory network involved in a critical pathogenic event of coronary atherosclerosis-the formation of lipid-laden foam cells. The relative genetic contribution to CAD was assessed by partitioning disease heritability across different macrophage subpopulations. Meta-analysis of single-cell RNA-seq data sets from 38 human atherosclerotic samples was conducted to provide high-resolution cross-referencing to macrophage subpopulations in vivo. RESULTS: We identified 18 782 cis-regulatory elements by jointly profiling the gene expression and chromatin accessibility of >5000 macrophages. Integration with CAD genome-wide association study data prioritized 121 CAD-related genetic variants and 56 candidate causal genes. We showed that CAD heritability was not uniformly distributed and was particularly enriched in the gene programs of a novel CD52-hi lipid-handling macrophage subpopulation. These CD52-hi macrophages displayed significantly less lipoprotein accumulation and were also found in human atherosclerotic plaques. We investigated the cis-regulatory effect of a risk variant rs10488763 on FDX1, implicating the recruitment of AP-1 and C/EBP-ß in the causal mechanisms at this locus. CONCLUSIONS: Our results provide genetic evidence of the divergent roles of macrophage subsets in atherogenesis and highlight lipid-handling macrophages as a key subpopulation through which genetic variants operate to influence disease. These findings provide an unbiased framework for functional fine-mapping of genome-wide association study results using single-cell multiomics and offer new insights into the genotype-environment interactions underlying atherosclerotic disease.

An alphacoronavirus polymerase structure reveals conserved replication factor functions.

Coronaviruses are a diverse subfamily of viruses containing pathogens of humans and animals. This subfamily of viruses replicates their RNA genomes using a core polymerase complex composed of viral non-structural proteins: nsp7, nsp8 and nsp12. Most of our understanding of coronavirus molecular biology comes from betacoronaviruses like SARS-CoV and SARS-CoV-2, the latter of which is the causative agent of COVID-19. In contrast, members of the alphacoronavirus genus are relatively understudied despite their importance in human and animal health. Here we have used cryo-electron microscopy to determine structures of the alphacoronavirus porcine epidemic diarrhea virus (PEDV) core polymerase complex bound to RNA. One structure shows an unexpected nsp8 stoichiometry despite remaining bound to RNA. Biochemical analysis shows that the N-terminal extension of one nsp8 is not required for in vitro RNA synthesis for alpha- and betacoronaviruses. Our work demonstrates the importance of studying diverse coronaviruses in revealing aspects of coronavirus replication and identifying areas of conservation to be targeted by antiviral drugs.

All exons are not created equal-exon vulnerability determines the effect of exonic mutations on splicing.

It is now widely accepted that aberrant splicing of constitutive exons is often caused by mutations affecting cis-acting splicing regulatory elements (SREs), but there is a misconception that all exons have an equal dependency on SREs and thus a similar vulnerability to aberrant splicing. We demonstrate that some exons are more likely to be affected by exonic splicing mutations (ESMs) due to an inherent vulnerability, which is context dependent and influenced by the strength of exon definition. We have developed VulExMap, a tool which is based on empirical data that can designate whether a constitutive exon is vulnerable. Using VulExMap, we find that only 25% of all exons can be categorized as vulnerable, whereas two-thirds of 359 previously reported ESMs in 75 disease genes are located in vulnerable exons. Because VulExMap analysis is based on empirical data on splicing of exons in their endogenous context, it includes all features important in determining the vulnerability. We believe that VulExMap will be an important tool when assessing the effect of exonic mutations by pinpointing whether they are located in exons vulnerable to ESMs.

Propulsive cell entry diverts pathogens from immune degradation by remodeling the phagocytic synapse.

Phagocytosis is a critical immune function for infection control and tissue homeostasis. During phagocytosis, pathogens are internalized and degraded in phagolysosomes. For pathogens that evade immune degradation, the prevailing view is that virulence factors are required to disrupt the biogenesis of phagolysosomes. In contrast, we present here that physical forces from motile pathogens during cell entry divert them away from the canonical degradative pathway. This altered fate begins with the force-induced remodeling of the phagocytic synapse formation. We used the parasite Toxoplasma gondii as a model because live Toxoplasma actively invades host cells using gliding motility. To differentiate the effects of physical forces from virulence factors in phagocytosis, we employed magnetic forces to induce propulsive entry of inactivated Toxoplasma into macrophages. Experiments and computer simulations show that large propulsive forces hinder productive activation of receptors by preventing their spatial segregation from phosphatases at the phagocytic synapse. Consequently, the inactivated parasites are engulfed into vacuoles that fail to mature into degradative units, similar to the live motile parasite's intracellular pathway. Using yeast cells and opsonized beads, we confirmed that this mechanism is general, not specific to the parasite used. These results reveal new aspects of immune evasion by demonstrating how physical forces during active cell entry, independent of virulence factors, enable pathogens to circumvent phagolysosomal degradation.

Heat shock potentiates aminoglycosides against gram-negative bacteria by enhancing antibiotic uptake, protein aggregation, and ROS.

The potentiation of antibiotics is a promising strategy for combatting antibiotic-resistant/tolerant bacteria. Herein, we report that a 5-min sublethal heat shock enhances the bactericidal actions of aminoglycoside antibiotics by six orders of magnitude against both exponential- and stationary-phase Escherichia coli. This combined treatment also effectively kills various E. coli persisters, E. coli clinical isolates, and numerous gram-negative but not gram-positive bacteria and enables aminoglycosides at 5% of minimum inhibitory concentrations to eradicate multidrug-resistant pathogens Acinetobacter baumannii and Klebsiella pneumoniae. Mechanistically, the potentiation is achieved comprehensively by heat shock-enhanced proton motive force that thus promotes the bacterial uptake of aminoglycosides, as well as by increasing irreversible protein aggregation and reactive oxygen species that further augment the downstream lethality of aminoglycosides. Consistently, protonophores, chemical chaperones, antioxidants, and anaerobic culturing abolish heat shock-enhanced aminoglycoside lethality. We also demonstrate as a proof of concept that infrared irradiation- or photothermal nanosphere-induced thermal treatments potentiate aminoglycoside killing of Pseudomonas aeruginosa in a mouse acute skin wound model. Our study advances the understanding of the mechanism of actions of aminoglycosides and demonstrates a high potential for thermal ablation in curing bacterial infections when combined with aminoglycosides.

"Self-inactivating" rabies viruses are susceptible to loss of their intended attenuating modification.

Monosynaptic tracing using rabies virus is an important technique in neuroscience, allowing brain-wide labeling of neurons directly presynaptic to a targeted neuronal population. A 2017 article reported the development of a noncytotoxic version-a major advance-based on attenuating the rabies virus by the addition of a destabilization domain to the C terminus of a viral protein. However, this modification did not appear to hinder the ability of the virus to spread between neurons. We analyzed two viruses provided by the authors and show here that both were mutants that had lost the intended modification, explaining the paper's paradoxical results. We then made a virus that actually did have the intended modification in at least the majority of virions and found that it did not spread efficiently under the conditions described in the original paper, namely, without an exogenous protease being expressed in order to remove the destabilization domain. We found that it did spread when the protease was supplied, although this also appeared to result in the deaths of most source cells by 3 wk postinjection. We conclude that the new approach is not robust but that it could become a viable technique given further optimization and validation.

Seizure Cycles under Pharmacotherapy.

OBJECTIVE: This study was undertaken to determine the effects of antiseizure medications (ASMs) on multidien (multiday) cycles of interictal epileptiform activity (IEA) and seizures and evaluate their potential clinical significance. METHODS: We retrospectively analyzed up to 10 years of data from 88 of the 256 total adults with pharmacoresistant focal epilepsy who participated in the clinical trials of the RNS System, an intracranial device that keeps records of IEA counts. Following adjunctive ASM trials, we evaluated changes over months in (1) rates of self-reported disabling seizures and (2) multidien IEA cycle strength (spectral power for periodicity between 4 and 40 days). We used a survival analysis and the receiver operating characteristics to assess changes in IEA as a predictor of seizure control. RESULTS: Among 56 (33.3%) of the 168 adjunctive ASM trials suitable for analysis, ASM introduction was followed by an average 50 to 70% decrease in multidien IEA cycle strength and a concomitant 50 to 70% decrease in relative seizure rate for up to 12 months. Individuals with a ≥50% decrease in IEA cycle strength in the first 3 months of an ASM trial had a higher probability of remaining seizure responders (≥50% seizure rate reduction, p < 10-7) or super-responders (≥90%, p < 10-8) over the next 12 months. INTERPRETATION: In this large cohort, a decrease in multidien IEA cycle strength following initiation of an adjunctive ASM correlated with seizure control for up to 12 months, suggesting that fluctuations in IEA mirror "disease activity" in pharmacoresistant focal epilepsy and may have clinical utility as a biomarker to predict treatment response. ANN NEUROL 2024;95:743-753.

MAST: Phylogenetic Inference with Mixtures Across Sites and Trees.

Hundreds or thousands of loci are now routinely used in modern phylogenomic studies. Concatenation approaches to tree inference assume that there is a single topology for the entire dataset, but different loci may have different evolutionary histories due to incomplete lineage sorting (ILS), introgression, and/or horizontal gene transfer; even single loci may not be treelike due to recombination. To overcome this shortcoming, we introduce an implementation of a multi-tree mixture model that we call mixtures across sites and trees (MAST). This model extends a prior implementation by Boussau et al. (2009) by allowing users to estimate the weight of each of a set of pre-specified bifurcating trees in a single alignment. The MAST model allows each tree to have its own weight, topology, branch lengths, substitution model, nucleotide or amino acid frequencies, and model of rate heterogeneity across sites. We implemented the MAST model in a maximum-likelihood framework in the popular phylogenetic software, IQ-TREE. Simulations show that we can accurately recover the true model parameters, including branch lengths and tree weights for a given set of tree topologies, under a wide range of biologically realistic scenarios. We also show that we can use standard statistical inference approaches to reject a single-tree model when data are simulated under multiple trees (and vice versa). We applied the MAST model to multiple primate datasets and found that it can recover the signal of ILS in the Great Apes, as well as the asymmetry in minor trees caused by introgression among several macaque species. When applied to a dataset of 4 Platyrrhine species for which standard concatenated maximum likelihood (ML) and gene tree approaches disagree, we observe that MAST gives the highest weight (i.e., the largest proportion of sites) to the tree also supported by gene tree approaches. These results suggest that the MAST model is able to analyze a concatenated alignment using ML while avoiding some of the biases that come with assuming there is only a single tree. We discuss how the MAST model can be extended in the future.

State-dependent effects of responsive neurostimulation depend on seizure localization.

Brain-responsive neurostimulation is firmly ensconced among treatment options for drug-resistant focal epilepsy, but over a quarter of patients treated with the RNS System do not experience meaningful seizure reduction. Initial titration of RNS therapy is typically similar for all patients, raising the possibility that treatment response might be enhanced by consideration of patient-specific variables. Indeed, small, single-center studies have yielded preliminary evidence that RNS System effectiveness depends on the brain state during which stimulation is applied. The generalizability of these findings remains unclear, however, and it is unknown whether state-dependent effects of responsive neurostimulation are also stratified by location of the seizure onset zone where stimulation is delivered. We aimed to determine whether state-dependent effects of the RNS System are evident in the large, diverse, multi-center cohort of RNS System clinical trial participants and to test whether these effects differ between mesiotemporal and neocortical epilepsies. Eighty-one of 256 patients who were treated with the RNS System across 31 centers during clinical trials met criteria for inclusion in this retrospective study. Risk states were defined in relation to phases of daily and multi-day cycles of interictal epileptiform activity that are thought to determine seizure likelihood. We found that the probabilities of risk state transitions depended on the stimulation parameter being changed, the starting seizure risk state, and the stimulated brain region. Changes in two commonly adjusted stimulation parameters, charge density and stimulation frequency, produced opposite effects on risk state transitions depending on seizure localization. Greater variance in acute risk state transitions was explained by state-dependent responsive neurostimulation for bipolar stimulation for neocortical epilepsies and for monopolar stimulation for mesiotemporal epilepsies. Variability in effectiveness of RNS System therapy across individuals may relate, at least partly, to the fact that current treatment paradigms do not account fully for fluctuations in brain states or locations of simulation sites. State-dependence of electrical brain stimulation may inform development of next-generation closed-loop devices that can detect changes in brain state and deliver adaptive, localization-specific patterns of stimulation to maximize therapeutic effects.

Palaeo-Eskimo genetic ancestry and the peopling of Chukotka and North America.

Much of the American Arctic was first settled 5,000 years ago, by groups of people known as Palaeo-Eskimos. They were subsequently joined and largely displaced around 1,000 years ago by ancestors of the present-day Inuit and Yup'ik1-3. The genetic relationship between Palaeo-Eskimos and Native American, Inuit, Yup'ik and Aleut populations remains uncertain4-6. Here we present genomic data for 48 ancient individuals from Chukotka, East Siberia, the Aleutian Islands, Alaska, and the Canadian Arctic. We co-analyse these data with data from present-day Alaskan Iñupiat and West Siberian populations and published genomes. Using methods based on rare-allele and haplotype sharing, as well as established techniques4,7-9, we show that Palaeo-Eskimo-related ancestry is ubiquitous among people who speak Na-Dene and Eskimo-Aleut languages. We develop a comprehensive model for the Holocene peopling events of Chukotka and North America, and show that Na-Dene-speaking peoples, people of the Aleutian Islands, and Yup'ik and Inuit across the Arctic region all share ancestry from a single Palaeo-Eskimo-related Siberian source.