A mouse DRG genetic toolkit reveals morphological and physiological diversity of somatosensory neuron subtypes.

Dorsal root ganglia (DRG) somatosensory neurons detect mechanical, thermal, and chemical stimuli acting on the body. Achieving a holistic view of how different DRG neuron subtypes relay neural signals from the periphery to the CNS has been challenging with existing tools. Here, we develop and curate a mouse genetic toolkit that allows for interrogating the properties and functions of distinct cutaneous targeting DRG neuron subtypes. These tools have enabled a broad morphological analysis, which revealed distinct cutaneous axon arborization areas and branching patterns of the transcriptionally distinct DRG neuron subtypes. Moreover, in vivo physiological analysis revealed that each subtype has a distinct threshold and range of responses to mechanical and/or thermal stimuli. These findings support a model in which morphologically and physiologically distinct cutaneous DRG sensory neuron subtypes tile mechanical and thermal stimulus space to collectively encode a wide range of natural stimuli.

Host and environmental factors shape upper airway microbiota and respiratory health across the human lifespan.

Our understanding of the normal variation in the upper respiratory tract (URT) microbiota across the human lifespan and how these relate to host, environment, and health is limited. We studied the microbiota of 3,104 saliva (<10 year-olds)/oropharynx (≥10 year-olds) and 2,485 nasopharynx samples of 3,160 Dutch individuals 0-87 years of age, participating in a cross-sectional population-wide study (PIENTER-3) using 16S-rRNA sequencing. The microbiota composition was strongly related to age, especially in the nasopharynx, with maturation occurring throughout childhood and adolescence. Clear niche- and age-specific associations were found between the microbiota composition and host/environmental factors and health outcomes. Among others, social interaction, sex, and season were associated with the nasopharyngeal microbial community. By contrast, the oral microbiota was more related to antibiotics, tobacco, and alcohol use. We present an atlas of the URT microbiota across the lifespan in association with environment and health, establishing a baseline for future research.

Detection and analysis of complex structural variation in human genomes across populations and in brains of donors with psychiatric disorders.

Complex structural variations (cxSVs) are often overlooked in genome analyses due to detection challenges. We developed ARC-SV, a probabilistic and machine-learning-based method that enables accurate detection and reconstruction of cxSVs from standard datasets. By applying ARC-SV across 4,262 genomes representing all continental populations, we identified cxSVs as a significant source of natural human genetic variation. Rare cxSVs have a propensity to occur in neural genes and loci that underwent rapid human-specific evolution, including those regulating corticogenesis. By performing single-nucleus multiomics in postmortem brains, we discovered cxSVs associated with differential gene expression and chromatin accessibility across various brain regions and cell types. Additionally, cxSVs detected in brains of psychiatric cases are enriched for linkage with psychiatric GWAS risk alleles detected in the same brains. Furthermore, our analysis revealed significantly decreased brain-region- and cell-type-specific expression of cxSV genes, specifically for psychiatric cases, implicating cxSVs in the molecular etiology of major neuropsychiatric disorders.

Population immunity predicts evolutionary trajectories of SARS-CoV-2.

The large-scale evolution of the SARS-CoV-2 virus has been marked by rapid turnover of genetic clades. New variants show intrinsic changes, notably increased transmissibility, and antigenic changes that reduce cross-immunity induced by previous infections or vaccinations. How this functional variation shapes global evolution has remained unclear. Here, we establish a predictive fitness model for SARS-CoV-2 that integrates antigenic and intrinsic selection. The model is informed by tracking of time-resolved sequence data, epidemiological records, and cross-neutralization data of viral variants. Our inference shows that immune pressure, including contributions of vaccinations and previous infections, has become the dominant force driving the recent evolution of SARS-CoV-2. The fitness model can serve continued surveillance in two ways. First, it successfully predicts the short-term evolution of circulating strains and flags emerging variants likely to displace the previously predominant variant. Second, it predicts likely antigenic profiles of successful escape variants prior to their emergence.

The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights.

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

The genetic history of Scandinavia from the Roman Iron Age to the present.

We investigate a 2,000-year genetic transect through Scandinavia spanning the Iron Age to the present, based on 48 new and 249 published ancient genomes and genotypes from 16,638 modern individuals. We find regional variation in the timing and magnitude of gene flow from three sources: the eastern Baltic, the British-Irish Isles, and southern Europe. British-Irish ancestry was widespread in Scandinavia from the Viking period, whereas eastern Baltic ancestry is more localized to Gotland and central Sweden. In some regions, a drop in current levels of external ancestry suggests that ancient immigrants contributed proportionately less to the modern Scandinavian gene pool than indicated by the ancestry of genomes from the Viking and Medieval periods. Finally, we show that a north-south genetic cline that characterizes modern Scandinavians is mainly due to the differential levels of Uralic ancestry and that this cline existed in the Viking Age and possibly earlier.

A genetic history of the Balkans from Roman frontier to Slavic migrations.

The rise and fall of the Roman Empire was a socio-political process with enormous ramifications for human history. The Middle Danube was a crucial frontier and a crossroads for population and cultural movement. Here, we present genome-wide data from 136 Balkan individuals dated to the 1st millennium CE. Despite extensive militarization and cultural influence, we find little ancestry contribution from peoples of Italic descent. However, we trace a large-scale influx of people of Anatolian ancestry during the Imperial period. Between ∼250 and 550 CE, we detect migrants with ancestry from Central/Northern Europe and the Steppe, confirming that "barbarian" migrations were propelled by ethnically diverse confederations. Following the end of Roman control, we detect the large-scale arrival of individuals who were genetically similar to modern Eastern European Slavic-speaking populations, who contributed 30%-60% of the ancestry of Balkan people, representing one of the largest permanent demographic changes anywhere in Europe during the Migration Period.

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.

Evolving ancient DNA techniques and the future of human history.

Ancient DNA (aDNA) techniques applied to human genomics have significantly advanced in the past decade, enabling large-scale aDNA research, sometimes independent of human remains. This commentary reviews the major milestones of aDNA techniques and explores future directions to expand the scope of aDNA research and insights into present-day human health.

Genome-wide data from medieval German Jews show that the Ashkenazi founder event pre-dated the 14th century.

We report genome-wide data from 33 Ashkenazi Jews (AJ), dated to the 14th century, obtained following a salvage excavation at the medieval Jewish cemetery of Erfurt, Germany. The Erfurt individuals are genetically similar to modern AJ, but they show more variability in Eastern European-related ancestry than modern AJ. A third of the Erfurt individuals carried a mitochondrial lineage common in modern AJ and eight carried pathogenic variants known to affect AJ today. These observations, together with high levels of runs of homozygosity, suggest that the Erfurt community had already experienced the major reduction in size that affected modern AJ. The Erfurt bottleneck was more severe, implying substructure in medieval AJ. Overall, our results suggest that the AJ founder event and the acquisition of the main sources of ancestry pre-dated the 14th century and highlight late medieval genetic heterogeneity no longer present in modern AJ.

En masse organoid phenotyping informs metabolic-associated genetic susceptibility to NASH.

Genotype-phenotype associations for common diseases are often compounded by pleiotropy and metabolic state. Here, we devised a pooled human organoid-panel of steatohepatitis to investigate the impact of metabolic status on genotype-phenotype association. En masse population-based phenotypic analysis under insulin insensitive conditions predicted key non-alcoholic steatohepatitis (NASH)-genetic factors including the glucokinase regulatory protein (GCKR)-rs1260326:C>T. Analysis of NASH clinical cohorts revealed that GCKR-rs1260326-T allele elevates disease severity only under diabetic state but protects from fibrosis under non-diabetic states. Transcriptomic, metabolomic, and pharmacological analyses indicate significant mitochondrial dysfunction incurred by GCKR-rs1260326, which was not reversed with metformin. Uncoupling oxidative mechanisms mitigated mitochondrial dysfunction and permitted adaptation to increased fatty acid supply while protecting against oxidant stress, forming a basis for future therapeutic approaches for diabetic NASH. Thus, "in-a-dish" genotype-phenotype association strategies disentangle the opposing roles of metabolic-associated gene variant functions and offer a rich mechanistic, diagnostic, and therapeutic inference toolbox toward precision hepatology. VIDEO ABSTRACT.

Synthetic mammalian signaling circuits for robust cell population control.

In multicellular organisms, cells actively sense and control their own population density. Synthetic mammalian quorum-sensing circuits could provide insight into principles of population control and extend cell therapies. However, a key challenge is reducing their inherent sensitivity to "cheater" mutations that evade control. Here, we repurposed the plant hormone auxin to enable orthogonal mammalian cell-cell communication and quorum sensing. We designed a paradoxical population control circuit, termed "Paradaux," in which auxin stimulates and inhibits net cell growth at different concentrations. This circuit limited population size over extended timescales of up to 42 days of continuous culture. By contrast, when operating in a non-paradoxical regime, population control became more susceptible to mutational escape. These results establish auxin as a versatile "private" communication system and demonstrate that paradoxical circuit architectures can provide robust population control.

Ancient genomes reveal origin and rapid trans-Eurasian migration of 7th century Avar elites.

The Avars settled the Carpathian Basin in 567/68 CE, establishing an empire lasting over 200 years. Who they were and where they came from is highly debated. Contemporaries have disagreed about whether they were, as they claimed, the direct successors of the Mongolian Steppe Rouran empire that was destroyed by the Turks in ∼550 CE. Here, we analyze new genome-wide data from 66 pre-Avar and Avar-period Carpathian Basin individuals, including the 8 richest Avar-period burials and further elite sites from Avar's empire core region. Our results provide support for a rapid long-distance trans-Eurasian migration of Avar-period elites. These individuals carried Northeast Asian ancestry matching the profile of preceding Mongolian Steppe populations, particularly a genome available from the Rouran period. Some of the later elite individuals carried an additional non-local ancestry component broadly matching the steppe, which could point to a later migration or reflect greater genetic diversity within the initial migrant population.

Genomic structure predicts metabolite dynamics in microbial communities.

The metabolic activities of microbial communities play a defining role in the evolution and persistence of life on Earth, driving redox reactions that give rise to global biogeochemical cycles. Community metabolism emerges from a hierarchy of processes, including gene expression, ecological interactions, and environmental factors. In wild communities, gene content is correlated with environmental context, but predicting metabolite dynamics from genomes remains elusive. Here, we show, for the process of denitrification, that metabolite dynamics of a community are predictable from the genes each member of the community possesses. A simple linear regression reveals a sparse and generalizable mapping from gene content to metabolite dynamics for genomically diverse bacteria. A consumer-resource model correctly predicts community metabolite dynamics from single-strain phenotypes. Our results demonstrate that the conserved impacts of metabolic genes can predict community metabolite dynamics, enabling the prediction of metabolite dynamics from metagenomes, designing denitrifying communities, and discovering how genome evolution impacts metabolism.

High-coverage whole-genome sequencing of the expanded 1000 Genomes Project cohort including 602 trios.

The 1000 Genomes Project (1kGP) is the largest fully open resource of whole-genome sequencing (WGS) data consented for public distribution without access or use restrictions. The final, phase 3 release of the 1kGP included 2,504 unrelated samples from 26 populations and was based primarily on low-coverage WGS. Here, we present a high-coverage 3,202-sample WGS 1kGP resource, which now includes 602 complete trios, sequenced to a depth of 30X using Illumina. We performed single-nucleotide variant (SNV) and short insertion and deletion (INDEL) discovery and generated a comprehensive set of structural variants (SVs) by integrating multiple analytic methods through a machine learning model. We show gains in sensitivity and precision of variant calls compared to phase 3, especially among rare SNVs as well as INDELs and SVs spanning frequency spectrum. We also generated an improved reference imputation panel, making variants discovered here accessible for association studies.

The genomic origins of the world's first farmers.

The precise genetic origins of the first Neolithic farming populations in Europe and Southwest Asia, as well as the processes and the timing of their differentiation, remain largely unknown. Demogenomic modeling of high-quality ancient genomes reveals that the early farmers of Anatolia and Europe emerged from a multiphase mixing of a Southwest Asian population with a strongly bottlenecked western hunter-gatherer population after the last glacial maximum. Moreover, the ancestors of the first farmers of Europe and Anatolia went through a period of extreme genetic drift during their westward range expansion, contributing highly to their genetic distinctiveness. This modeling elucidates the demographic processes at the root of the Neolithic transition and leads to a spatial interpretation of the population history of Southwest Asia and Europe during the late Pleistocene and early Holocene.

The deep population history of northern East Asia from the Late Pleistocene to the Holocene.

Northern East Asia was inhabited by modern humans as early as 40 thousand years ago (ka), as demonstrated by the Tianyuan individual. Using genome-wide data obtained from 25 individuals dated to 33.6-3.4 ka from the Amur region, we show that Tianyuan-related ancestry was widespread in northern East Asia before the Last Glacial Maximum (LGM). At the close of the LGM stadial, the earliest northern East Asian appeared in the Amur region, and this population is basal to ancient northern East Asians. Human populations in the Amur region have maintained genetic continuity from 14 ka, and these early inhabitants represent the closest East Asian source known for Ancient Paleo-Siberians. We also observed that EDAR V370A was likely to have been elevated to high frequency after the LGM, suggesting the possible timing for its selection. This study provides a deep look into the population dynamics of northern East Asia.

High-precision coding in visual cortex.

Individual neurons in visual cortex provide the brain with unreliable estimates of visual features. It is not known whether the single-neuron variability is correlated across large neural populations, thus impairing the global encoding of stimuli. We recorded simultaneously from up to 50,000 neurons in mouse primary visual cortex (V1) and in higher order visual areas and measured stimulus discrimination thresholds of 0.35° and 0.37°, respectively, in an orientation decoding task. These neural thresholds were almost 100 times smaller than the behavioral discrimination thresholds reported in mice. This discrepancy could not be explained by stimulus properties or arousal states. Furthermore, behavioral variability during a sensory discrimination task could not be explained by neural variability in V1. Instead, behavior-related neural activity arose dynamically across a network of non-sensory brain areas. These results imply that perceptual discrimination in mice is limited by downstream decoders, not by neural noise in sensory representations.

The genomic history of the Aegean palatial civilizations.

The Cycladic, the Minoan, and the Helladic (Mycenaean) cultures define the Bronze Age (BA) of Greece. Urbanism, complex social structures, craft and agricultural specialization, and the earliest forms of writing characterize this iconic period. We sequenced six Early to Middle BA whole genomes, along with 11 mitochondrial genomes, sampled from the three BA cultures of the Aegean Sea. The Early BA (EBA) genomes are homogeneous and derive most of their ancestry from Neolithic Aegeans, contrary to earlier hypotheses that the Neolithic-EBA cultural transition was due to massive population turnover. EBA Aegeans were shaped by relatively small-scale migration from East of the Aegean, as evidenced by the Caucasus-related ancestry also detected in Anatolians. In contrast, Middle BA (MBA) individuals of northern Greece differ from EBA populations in showing ∼50% Pontic-Caspian Steppe-related ancestry, dated at ca. 2,600-2,000 BCE. Such gene flow events during the MBA contributed toward shaping present-day Greek genomes.

The genomic history of the Middle East.

The Middle East region is important to understand human evolution and migrations but is underrepresented in genomic studies. Here, we generated 137 high-coverage physically phased genome sequences from eight Middle Eastern populations using linked-read sequencing. We found no genetic traces of early expansions out-of-Africa in present-day populations but found Arabians have elevated Basal Eurasian ancestry that dilutes their Neanderthal ancestry. Population sizes within the region started diverging 15-20 kya, when Levantines expanded while Arabians maintained smaller populations that derived ancestry from local hunter-gatherers. Arabians suffered a population bottleneck around the aridification of Arabia 6 kya, while Levantines had a distinct bottleneck overlapping the 4.2 kya aridification event. We found an association between movement and admixture of populations in the region and the spread of Semitic languages. Finally, we identify variants that show evidence of selection, including polygenic selection. Our results provide detailed insights into the genomic and selective histories of the Middle East.