Obsessive-Compulsive Behavior Isn't Necessarily a Bad Thing.

It is difficult to believe that in about 1960 practically nothing was known about the thymus and some of its products, T cells bearing αß receptors for antigen. Thus I was lucky to join the field of T cell biology almost at its beginning, when knowledge about the cells was just getting off the ground and there was so much to discover. This article describes findings about these cells made by others and myself that led us all from ignorance, via complete confusion, to our current state of knowledge. I believe I was fortunate to practice science in very supportive institutions and with very collaborative colleagues in two countries that both encourage independent research by independent scientists, while simultaneously ignoring or somehow being able to avoid some of the difficulties of being a woman in what was, at the time, a male-dominated profession.

Natural and Engineered Guide RNA-Directed Transposition with CRISPR-Associated Tn7-Like Transposons.

CRISPR-Cas (clustered regularly interspaced short palindromic repeats-CRISPR-associated nuclease) defense systems have been naturally coopted for guide RNA-directed transposition on multiple occasions. In all cases, cooption occurred with diverse elements related to the bacterial transposon Tn7. Tn7 tightly controls transposition; the transposase is activated only when special targets are recognized by dedicated target-site selection proteins. Tn7 and the Tn7-like elements that coopted CRISPR-Cas systems evolved complementary targeting pathways: one that recognizes a highly conserved site in the chromosome and a second pathway that targets mobile plasmids capable of cell-to-cell transfer. Tn7 and Tn7-like elements deliver a single integration into the site they recognize and also control the orientation of the integration event, providing future potential for use as programmable gene-integration tools. Early work has shown that guide RNA-directed transposition systems can be adapted to diverse hosts, even within microbial communities, suggesting great potential for engineering these systems as powerful gene-editing tools.

The rate and nature of mitochondrial DNA mutations in human pedigrees.

We examined the rate and nature of mitochondrial DNA (mtDNA) mutations in humans using sequence data from 64,806 contemporary Icelanders from 2,548 matrilines. Based on 116,663 mother-child transmissions, 8,199 mutations were detected, providing robust rate estimates by nucleotide type, functional impact, position, and different alleles at the same position. We thoroughly document the true extent of hypermutability in mtDNA, mainly affecting the control region but also some coding-region variants. The results reveal the impact of negative selection on viable deleterious mutations, including rapidly mutating disease-associated 3243A>G and 1555A>G and pre-natal selection that most likely occurs during the development of oocytes. Finally, we show that the fate of new mutations is determined by a drastic germline bottleneck, amounting to an average of 3 mtDNA units effectively transmitted from mother to child.

The WDR11 complex is a receptor for acidic-cluster-containing cargo proteins.

Vesicle trafficking is a fundamental process that allows for the sorting and transport of specific proteins (i.e., "cargoes") to different compartments of eukaryotic cells. Cargo recognition primarily occurs through coats and the associated proteins at the donor membrane. However, it remains unclear whether cargoes can also be selected at other stages of vesicle trafficking to further enhance the fidelity of the process. The WDR11-FAM91A1 complex functions downstream of the clathrin-associated AP-1 complex to facilitate protein transport from endosomes to the TGN. Here, we report the cryo-EM structure of human WDR11-FAM91A1 complex. WDR11 directly and specifically recognizes a subset of acidic clusters, which we term super acidic clusters (SACs). WDR11 complex assembly and its binding to SAC-containing proteins are indispensable for the trafficking of SAC-containing proteins and proper neuronal development in zebrafish. Our studies thus uncover that cargo proteins could be recognized in a sequence-specific manner downstream of a protein coat.

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.

The evolution and diversification of oakleaf butterflies.

Oakleaf butterflies in the genus Kallima have a polymorphic wing phenotype, enabling these insects to masquerade as dead leaves. This iconic example of protective resemblance provides an interesting evolutionary paradigm that can be employed to study biodiversity. We integrated multi-omic data analyses and functional validation to infer the evolutionary history of Kallima species and investigate the genetic basis of their variable leaf wing patterns. We find that Kallima butterflies diversified in the eastern Himalayas and dispersed to East and Southeast Asia. Moreover, we find that leaf wing polymorphism is controlled by the wing patterning gene cortex, which has been maintained in Kallima by long-term balancing selection. Our results provide macroevolutionary and microevolutionary insights into a model species originating from a mountain ecosystem.

Domestic dog lineages reveal genetic drivers of behavioral diversification.

Selective breeding of domestic dogs has generated diverse breeds often optimized for performing specialized tasks. Despite the heritability of breed-typical behavioral traits, identification of causal loci has proven challenging due to the complexity of canine population structure. We overcome longstanding difficulties in identifying genetic drivers of canine behavior by developing a framework for understanding relationships between breeds and the behaviors that define them, utilizing genetic data for over 4,000 domestic, semi-feral, and wild canids and behavioral survey data for over 46,000 dogs. We identify ten major canine genetic lineages and their behavioral correlates and show that breed diversification is predominantly driven by non-coding regulatory variation. We determine that lineage-associated genes converge in neurodevelopmental co-expression networks, identifying a sheepdog-associated enrichment for interrelated axon guidance functions. This work presents a scaffold for canine diversification that positions the domestic dog as an unparalleled system for revealing the genetic origins of behavioral diversity.

Primate hemorrhagic fever-causing arteriviruses are poised for spillover to humans.

Simian arteriviruses are endemic in some African primates and can cause fatal hemorrhagic fevers when they cross into primate hosts of new species. We find that CD163 acts as an intracellular receptor for simian hemorrhagic fever virus (SHFV; a simian arterivirus), a rare mode of virus entry that is shared with other hemorrhagic fever-causing viruses (e.g., Ebola and Lassa viruses). Further, SHFV enters and replicates in human monocytes, indicating full functionality of all of the human cellular proteins required for viral replication. Thus, simian arteriviruses in nature may not require major adaptations to the human host. Given that at least three distinct simian arteriviruses have caused fatal infections in captive macaques after host-switching, and that humans are immunologically naive to this family of viruses, development of serology tests for human surveillance should be a priority.

The emergence and ongoing convergent evolution of the SARS-CoV-2 N501Y lineages.

The independent emergence late in 2020 of the B.1.1.7, B.1.351, and P.1 lineages of SARS-CoV-2 prompted renewed concerns about the evolutionary capacity of this virus to overcome public health interventions and rising population immunity. Here, by examining patterns of synonymous and non-synonymous mutations that have accumulated in SARS-CoV-2 genomes since the pandemic began, we find that the emergence of these three "501Y lineages" coincided with a major global shift in the selective forces acting on various SARS-CoV-2 genes. Following their emergence, the adaptive evolution of 501Y lineage viruses has involved repeated selectively favored convergent mutations at 35 genome sites, mutations we refer to as the 501Y meta-signature. The ongoing convergence of viruses in many other lineages on this meta-signature suggests that it includes multiple mutation combinations capable of promoting the persistence of diverse SARS-CoV-2 lineages in the face of mounting host immune recognition.

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.

Cellular bases of olfactory circuit assembly revealed by systematic time-lapse imaging.

Neural circuit assembly features simultaneous targeting of numerous neuronal processes from constituent neuron types, yet the dynamics is poorly understood. Here, we use the Drosophila olfactory circuit to investigate dynamic cellular processes by which olfactory receptor neurons (ORNs) target axons precisely to specific glomeruli in the ipsi- and contralateral antennal lobes. Time-lapse imaging of individual axons from 30 ORN types revealed a rich diversity in extension speed, innervation timing, and ipsilateral branch locations and identified that ipsilateral targeting occurs via stabilization of transient interstitial branches. Fast imaging using adaptive optics-corrected lattice light-sheet microscopy showed that upon approaching target, many ORN types exhibiting "exploring branches" consisted of parallel microtubule-based terminal branches emanating from an F-actin-rich hub. Antennal nerve ablations uncovered essential roles for bilateral axons in contralateral target selection and for ORN axons to facilitate dendritic refinement of postsynaptic partner neurons. Altogether, these observations provide cellular bases for wiring specificity establishment.

Jpx RNA regulates CTCF anchor site selection and formation of chromosome loops.

Chromosome loops shift dynamically during development, homeostasis, and disease. CCCTC-binding factor (CTCF) is known to anchor loops and construct 3D genomes, but how anchor sites are selected is not yet understood. Here, we unveil Jpx RNA as a determinant of anchor selectivity. Jpx RNA targets thousands of genomic sites, preferentially binding promoters of active genes. Depleting Jpx RNA causes ectopic CTCF binding, massive shifts in chromosome looping, and downregulation of >700 Jpx target genes. Without Jpx, thousands of lost loops are replaced by de novo loops anchored by ectopic CTCF sites. Although Jpx controls CTCF binding on a genome-wide basis, it acts selectively at the subset of developmentally sensitive CTCF sites. Specifically, Jpx targets low-affinity CTCF motifs and displaces CTCF protein through competitive inhibition. We conclude that Jpx acts as a CTCF release factor and shapes the 3D genome by regulating anchor site usage.

Evolutionary Selection and Constraint on Human Knee Chondrocyte Regulation Impacts Osteoarthritis Risk.

During human evolution, the knee adapted to the biomechanical demands of bipedalism by altering chondrocyte developmental programs. This adaptive process was likely not without deleterious consequences to health. Today, osteoarthritis occurs in 250 million people, with risk variants enriched in non-coding sequences near chondrocyte genes, loci that likely became optimized during knee evolution. We explore this relationship by epigenetically profiling joint chondrocytes, revealing ancient selection and recent constraint and drift on knee regulatory elements, which also overlap osteoarthritis variants that contribute to disease heritability by tending to modify constrained functional sequence. We propose a model whereby genetic violations to regulatory constraint, tolerated during knee development, lead to adult pathology. In support, we discover a causal enhancer variant (rs6060369) present in billions of people at a risk locus (GDF5-UQCC1), showing how it impacts mouse knee-shape and osteoarthritis. Overall, our methods link an evolutionarily novel aspect of human anatomy to its pathogenesis.

Cerebellar Neurodynamics Predict Decision Timing and Outcome on the Single-Trial Level.

Goal-directed behavior requires the interaction of multiple brain regions. How these regions and their interactions with brain-wide activity drive action selection is less understood. We have investigated this question by combining whole-brain volumetric calcium imaging using light-field microscopy and an operant-conditioning task in larval zebrafish. We find global, recurring dynamics of brain states to exhibit pre-motor bifurcations toward mutually exclusive decision outcomes. These dynamics arise from a distributed network displaying trial-by-trial functional connectivity changes, especially between cerebellum and habenula, which correlate with decision outcome. Within this network the cerebellum shows particularly strong and predictive pre-motor activity (>10 s before movement initiation), mainly within the granule cells. Turn directions are determined by the difference neuroactivity between the ipsilateral and contralateral hemispheres, while the rate of bi-hemispheric population ramping quantitatively predicts decision time on the trial-by-trial level. Our results highlight a cognitive role of the cerebellum and its importance in motor planning.

VEGAS as a Platform for Facile Directed Evolution in Mammalian Cells.

Directed evolution, artificial selection toward designed objectives, is routinely used to develop new molecular tools and therapeutics. Successful directed molecular evolution campaigns repeatedly test diverse sequences with a designed selective pressure. Unicellular organisms and their viral pathogens are exceptional for this purpose and have been used for decades. However, many desirable targets of directed evolution perform poorly or unnaturally in unicellular backgrounds. Here, we present a system for facile directed evolution in mammalian cells. Using the RNA alphavirus Sindbis as a vector for heredity and diversity, we achieved 24-h selection cycles surpassing 10-3 mutations per base. Selection is achieved through genetically actuated sequences internal to the host cell, thus the system's name: viral evolution of genetically actuating sequences, or "VEGAS." Using VEGAS, we evolve transcription factors, GPCRs, and allosteric nanobodies toward functional signaling endpoints each in less than 1 weeks' time.

Tracking Five Millennia of Horse Management with Extensive Ancient Genome Time Series.

Horse domestication revolutionized warfare and accelerated travel, trade, and the geographic expansion of languages. Here, we present the largest DNA time series for a non-human organism to date, including genome-scale data from 149 ancient animals and 129 ancient genomes (≥1-fold coverage), 87 of which are new. This extensive dataset allows us to assess the modern legacy of past equestrian civilizations. We find that two extinct horse lineages existed during early domestication, one at the far western (Iberia) and the other at the far eastern range (Siberia) of Eurasia. None of these contributed significantly to modern diversity. We show that the influence of Persian-related horse lineages increased following the Islamic conquests in Europe and Asia. Multiple alleles associated with elite-racing, including at the MSTN "speed gene," only rose in popularity within the last millennium. Finally, the development of modern breeding impacted genetic diversity more dramatically than the previous millennia of human management.

Screening Human Embryos for Polygenic Traits Has Limited Utility.

The increasing proportion of variance in human complex traits explained by polygenic scores, along with progress in preimplantation genetic diagnosis, suggests the possibility of screening embryos for traits such as height or cognitive ability. However, the expected outcomes of embryo screening are unclear, which undermines discussion of associated ethical concerns. Here, we use theory, simulations, and real data to evaluate the potential gain of embryo screening, defined as the difference in trait value between the top-scoring embryo and the average embryo. The gain increases very slowly with the number of embryos but more rapidly with the variance explained by the score. Given current technology, the average gain due to screening would be ≈2.5 cm for height and ≈2.5 IQ points for cognitive ability. These mean values are accompanied by wide prediction intervals, and indeed, in large nuclear families, the majority of children top-scoring for height are not the tallest.

Exon-Mediated Activation of Transcription Starts.

The processing of RNA transcripts from mammalian genes occurs in proximity to their transcription. Here, we describe a phenomenon affecting thousands of genes that we call exon-mediated activation of transcription starts (EMATS), in which the splicing of internal exons impacts promoter choice and the expression level of the gene. We observed that evolutionary gain of internal exons is associated with gain of new transcription start sites (TSSs) nearby and increased gene expression. Inhibiting exon splicing reduced transcription from nearby promoters, and creation of new spliced exons activated transcription from cryptic promoters. The strongest effects occurred for weak promoters located proximal and upstream of efficiently spliced exons. Together, our findings support a model in which splicing recruits transcription machinery locally to influence TSS choice and identify exon gain, loss, and regulatory change as major contributors to the evolution of alternative promoters and gene expression in mammals.

Large-Scale Whole-Genome Sequencing of Three Diverse Asian Populations in Singapore.

Underrepresentation of Asian genomes has hindered population and medical genetics research on Asians, leading to population disparities in precision medicine. By whole-genome sequencing of 4,810 Singapore Chinese, Malays, and Indians, we found 98.3 million SNPs and small insertions or deletions, over half of which are novel. Population structure analysis demonstrated great representation of Asian genetic diversity by three ethnicities in Singapore and revealed a Malay-related novel ancestry component. Furthermore, demographic inference suggested that Malays split from Chinese ∼24,800 years ago and experienced significant admixture with East Asians ∼1,700 years ago, coinciding with the Austronesian expansion. Additionally, we identified 20 candidate loci for natural selection, 14 of which harbored robust associations with complex traits and diseases. Finally, we show that our data can substantially improve genotype imputation in diverse Asian and Oceanian populations. These results highlight the value of our data as a resource to empower human genetics discovery across broad geographic regions.

Angiotensin Analogs with Divergent Bias Stabilize Distinct Receptor Conformations.

"Biased" G protein-coupled receptor (GPCR) agonists preferentially activate pathways mediated by G proteins or ß-arrestins. Here, we use double electron-electron resonance spectroscopy to probe the changes that ligands induce in the conformational distribution of the angiotensin II type I receptor. Monitoring distances between 10 pairs of nitroxide labels distributed across the intracellular regions enabled mapping of four underlying sets of conformations. Ligands from different functional classes have distinct, characteristic effects on the conformational heterogeneity of the receptor. Compared to angiotensin II, the endogenous agonist, agonists with enhanced Gq coupling more strongly stabilize an "open" conformation with an accessible transducer-binding site. ß-arrestin-biased agonists deficient in Gq coupling do not stabilize this open conformation but instead favor two more occluded conformations. These data suggest a structural mechanism for biased ligand action at the angiotensin receptor that can be exploited to rationally design GPCR-targeting drugs with greater specificity of action.