The 50 most-cited articles in eLife, 2022

Background: Measures to quantify changes in the pace of biological aging in response to intervention are needed to evaluate geroprotective interventions for humans. Previously, we showed that quantification of the pace of biological aging from a DNA-methylation blood test was possible (Belsky et al., 2020). Here, we report a next-generation DNA-methylation biomarker of Pace of Aging, DunedinPACE (for Pace of Aging Calculated from the Epigenome). Methods: We used data from the Dunedin Study 1972-1973 birth cohort tracking within-individual decline in 19 indicators of organ-system integrity across four time points spanning two decades to model Pace of Aging. We distilled this two-decade Pace of Aging into a single-time-point DNA-methylation blood-test using elastic-net regression and a DNA-methylation dataset restricted to exclude probes with low test-retest reliability. We evaluated the resulting measure, named DunedinPACE, in five additional datasets. Results: DunedinPACE showed high test-retest reliability, was associated with morbidity, disability, and mortality, and indicated…

ABSTRACT Equilibrium fluctuations and triggered conformational changes often underlie the functional cycles of membrane proteins. For example, transporters mediate the passage of molecules across cell membranes by alternating between inward-facing (IF) and outward-facing (OF) states, while receptors undergo intracellular structural rearrangements that initiate signaling cascades. Although the conformational plasticity of these proteins has historically posed a challenge for traditional de novo protein structure prediction pipelines, the recent success of AlphaFold2 (AF2) in CASP14 culminated in the modeling of a transporter in multiple conformations to high accuracy. Given that AF2 was designed to predict static structures of proteins, it remains unclear if this result represents an underexplored capability to accurately predict multiple conformations and/or structural heterogeneity. Here, we present an approach to drive AF2 to sample alternative conformations of topologically diverse transporters and G-protein coupled receptors (GPCRs) that are absent from the AF2 training set.…

Contemporary bioinformatic and chemoinformatic capabilities hold promise to reshape knowledge management, analysis and interpretation of data in natural products research. Currently, reliance on a disparate set of non-standardized, insular, and specialized databases presents a series of challenges for data access, both within the discipline and for integration and interoperability between related fields. The fundamental elements of exchange are referenced structure-organism pairs that establish relationships between distinct molecular structures and the living organisms from which they were identified. Consolidating and sharing such information via an open platform has strong transformative potential for natural products research and beyond. This is the ultimate goal of the newly established LOTUS initiative, which has now completed the first steps toward the harmonization, curation, validation and open dissemination of 750,000+ referenced structure-organism pairs. LOTUS data is hosted on Wikidata and regularly mirrored on https://lotus.naturalproducts.net. Data sharing within the…

Brain organoids have been used to recapitulate the processes of brain development and related diseases. However, the lack of vasculatures, which regulate neurogenesis and brain disorders, limits the utility of brain organoids. In this study, we induced vessel and brain organoids, respectively, and then fused two types of organoids together to obtain vascularized brain organoids. The fused brain organoids were engrafted with robust vascular network-like structures and exhibited increased number of neural progenitors, in line with the possibility that vessels regulate neural development. Fusion organoids also contained functional blood-brain barrier-like structures, as well as microglial cells, a specific population of immune cells in the brain. The incorporated microglia responded actively to immune stimuli to the fused brain organoids and showed ability of engulfing synapses. Thus, the fusion organoids established in this study allow modeling interactions between the neuronal and non-neuronal components in vitro, particularly the vasculature and microglia niche.

The development of multiple chromosome-scale reference genome sequences in many taxonomic groups has yielded a high-resolution view of the patterns and processes of molecular evolution. Nonetheless, leveraging information across multiple genomes remains a significant challenge in nearly all eukaryotic systems. These challenges range from studying the evolution of chromosome structure, to finding candidate genes for quantitative trait loci, to testing hypotheses about speciation and adaptation. Here, we present GENESPACE, which addresses these challenges by integrating conserved gene order and orthology to define the expected physical position of all genes across multiple genomes. We demonstrate this utility by dissecting presence-absence, copy-number, and structural variation at three levels of biological organization: spanning 300 million years of vertebrate sex chromosome evolution, across the diversity of the Poaceae (grass) plant family, and among 26 maize cultivars. The methods to build and visualize syntenic orthology in the GENESPACE R package offer a significant addition to…

Abstract The peptidoglycan cell wall is a predominant structure of bacteria, determining cell shape and supporting survival in diverse conditions. Peptidoglycan is dynamic and requires regulated synthesis of new material, remodeling, and turnover – or autolysis – of old material. Despite exploitation of peptidoglycan synthesis as an antibiotic target, we lack a fundamental understanding of how peptidoglycan synthesis and autolysis intersect to maintain the cell wall. Here, we uncover a critical physiological role for a widely misunderstood class of autolytic enzymes, lytic transglycosylases (LTGs). We demonstrate that LTG activity is essential to survival by contributing to periplasmic processes upstream and independent of peptidoglycan recycling. Defects accumulate in Vibrio cholerae LTG mutants due to generally inadequate LTG activity, rather than absence of specific enzymes, and essential LTG activities are likely independent of protein-protein interactions, as heterologous expression of a non-native LTG rescues growth of a conditionally LTG-null mutant. Lastly, we demonstrate…

We present a new approach for macromolecular structure determination from multiple particles in electron cryo-tomography (cryo-ET) data sets. Whereas existing subtomogram averaging approaches are based on 3D data models, we propose to optimise a regularised likelihood target that approximates a function of the 2D experimental images. In addition, analogous to Bayesian polishing and contrast transfer function (CTF) refinement in single-particle analysis, we describe the approaches that exploit the increased signal-to-noise ratio in the averaged structure to optimise tilt-series alignments, beam-induced motions of the particles throughout the tilt-series acquisition, defoci of the individual particles, as well as higher-order optical aberrations of the microscope. Implementation of our approaches in the open-source software package RELION aims to facilitate their general use, particularly for those researchers who are already familiar with its single-particle analysis tools. We illustrate for three applications that our approaches allow structure determination from cryo-ET data to…

Abundant filamentous inclusions of tau are characteristic of more than 20 neurodegenerative diseases that are collectively termed tauopathies. Electron cryo-microscopy (cryo-EM) structures of tau amyloid filaments from human brain revealed that distinct tau folds characterise many different diseases. A lack of laboratory-based model systems to generate these structures has hampered efforts to uncover the molecular mechanisms that underlie tauopathies. Here, we report in vitro assembly conditions with recombinant tau that replicate the structures of filaments from both Alzheimer's disease (AD) and chronic traumatic encephalopathy (CTE), as determined by cryo-EM. Our results suggest that post-translational modifications of tau modulate filament assembly, and that previously observed additional densities in AD and CTE filaments may arise from the presence of inorganic salts, like phosphates and sodium chloride. In vitro assembly of tau into disease-relevant filaments will facilitate studies to determine their roles in different diseases, as well as the development of compounds that…

Abstract The recognition that individual GPCRs can activate multiple signaling pathways has raised the possibility of developing drugs selectively targeting therapeutically relevant ones. This requires tools to determine which G proteins and βarrestins are activated by a given receptor. Here, we present a set of BRET sensors monitoring the activation of the 12 G protein subtypes based on the translocation of their effectors to the plasma membrane (EMTA). Unlike most of the existing detection systems, EMTA does not require modification of receptors or G proteins (except for G s ). EMTA was found to be suitable for the detection of constitutive activity, inverse agonism, biased signaling and polypharmacology. Profiling of 100 therapeutically relevant human GPCRs resulted in 1,500 pathway-specific concentration-response curves and revealed a great diversity of coupling profiles ranging from exquisite selectivity to broad promiscuity. Overall, this work describes unique resources for studying the complexities underlying GPCR signaling and pharmacology.

Abstract Background Type 2 diabetes mellitus is known to be associated with neurobiological and cognitive deficits; however, their extent, overlap with aging effects, and the effectiveness of existing treatments in the context of the brain are currently unknown. Methods We characterized neurocognitive effects independently associated with T2DM and age in a large cohort of human subjects from the UK Biobank with cross-sectional neuroimaging and cognitive data. We then proceeded to evaluate the extent of overlap between the effects related to T2DM and age by applying correlation measures to the separately characterized neurocognitive changes. Our findings were complemented by meta-analyses of published reports with cognitive or neuroimaging measures for T2DM and healthy controls (HC). We also evaluated in a cohort of T2DM diagnosed individuals using UK Biobank how disease chronicity and metformin treatment interact with the identified neurocognitive effects. Results The UK Biobank dataset included cognitive and neuroimaging data (N=20,314) including 1,012 T2DM and 19,302 HC, aged…

Defining reference models for population variation, and the ability to study individual deviations is essential for understanding inter-individual variability and its relation to the onset and progression of medical conditions. In this work, we assembled a reference cohort of neuroimaging data from 82 sites (N=58,836; ages 2-100) and used normative modeling to characterize lifespan trajectories of cortical thickness and subcortical volume. Models are validated against a manually quality checked subset (N=24,354) and we provide an interface for transferring to new data sources. We showcase the clinical value by applying the models to a transdiagnostic psychiatric sample (N=1985), showing they can be used to quantify variability underlying multiple disorders whilst also refining case-control inferences. These models will be augmented with additional samples and imaging modalities as they become available. This provides a common reference platform to bind results from different studies and ultimately paves the way for personalized clinical decision-making.

Expansion microscopy (ExM) is a powerful technique to overcome the diffraction limit of light microscopy that can be applied in both tissues and cells. In ExM, samples are embedded in a swellable polymer gel to physically expand the sample and isotropically increase resolution in x, y, and z. The maximum resolution increase is limited by the expansion factor of the gel, which is four-fold for the original ExM protocol. Variations on the original ExM method have been reported that allow for greater expansion factors but at the cost of ease of adoption or versatility. Here, we systematically explore the ExM recipe space and present a novel method termed Ten-fold Robust Expansion Microscopy (TREx) that, like the original ExM method, requires no specialized equipment or procedures. We demonstrate that TREx gels expand 10-fold, can be handled easily, and can be applied to both thick mouse brain tissue sections and cultured human cells enabling high-resolution subcellular imaging with a single expansion step. Furthermore, we show that TREx can provide ultrastructural context to…

Abstract Ageing is the gradual decline in organismal fitness that occurs over time leading to tissue dysfunction and disease. At the cellular level, ageing is associated with reduced function, altered gene expression and a perturbed epigenome. Somatic cell reprogramming, the process of converting somatic cells to induced pluripotent stem cells (iPSCs), can reverse these age-associated changes. However, during iPSC reprogramming somatic cell identity is lost, and can be difficult to reacquire as re-differentiated iPSCs often resemble foetal rather than mature adult cells. Recent work has demonstrated that the epigenome is already rejuvenated by the maturation phase of reprogramming, which suggests full iPSC reprogramming is not required to reverse ageing of somatic cells. Here we have developed the first “maturation phase transient reprogramming” (MPTR) method, where reprogramming factors are expressed until this rejuvenation point followed by withdrawal of their induction. Using dermal fibroblasts from middle age donors, we found that cells reacquire their fibroblast identity…

The neurophysiology of cells and tissues are monitored electrophysiologically and optically in diverse experiments and species, ranging from flies to humans. Understanding the brain requires integration of data across this diversity, and thus these data must be findable, accessible, interoperable, and reusable (FAIR). This requires a standard language for data and metadata that can coevolve with neuroscience. We describe design and implementation principles for a language for neurophysiology data. Our open-source software (Neurodata Without Borders, NWB) defines and modularizes the interdependent, yet separable, components of a data language. We demonstrate NWB's impact through unified description of neurophysiology data across diverse modalities and species. NWB exists in an ecosystem, which includes data management, analysis, visualization, and archive tools. Thus, the NWB data language enables reproduction, interchange, and reuse of diverse neurophysiology data. More broadly, the design principles of NWB are generally applicable to enhance discovery across biology through data…

Throughout development, the brain transits from early highly synchronous activity patterns to a mature state with sparse and decorrelated neural activity, yet the mechanisms underlying this process are poorly understood. The developmental transition has important functional consequences, as the latter state is thought to allow for more efficient storage, retrieval, and processing of information. Here, we show that, in the mouse medial prefrontal cortex (mPFC), neural activity during the first two postnatal weeks decorrelates following specific spatial patterns. This process is accompanied by a concomitant tilting of excitation-inhibition (E-I) ratio toward inhibition. Using optogenetic manipulations and neural network modeling, we show that the two phenomena are mechanistically linked, and that a relative increase of inhibition drives the decorrelation of neural activity. Accordingly, in mice mimicking the etiology of neurodevelopmental disorders, subtle alterations in E-I ratio are associated with specific impairments in the correlational structure of spike trains. Finally,…

ABSTRACT Advances in modern artificial intelligence (AI) have inspired a paradigm shift in human neuroscience, yielding large-scale functional magnetic resonance imaging (fMRI) datasets that provide high-resolution brain responses to tens of thousands of naturalistic visual stimuli. Because such experiments necessarily involve brief stimulus durations and few repetitions of each stimulus, achieving sufficient signal-to-noise ratio can be a major challenge. We address this challenge by introducing GLMsingle , a scalable, user-friendly toolbox available in MATLAB and Python that enables accurate estimation of single-trial fMRI responses ( glmsingle.org ). Requiring only fMRI time-series data and a design matrix as inputs, GLMsingle integrates three techniques for improving the accuracy of trial-wise general linear model (GLM) beta estimates. First, for each voxel, a custom hemodynamic response function (HRF) is identified from a library of candidate functions. Second, cross-validation is used to derive a set of noise regressors from voxels unrelated to the experimental paradigm.…

Background: The combined impact of immunity and SARS-CoV-2 variants on viral kinetics during infections has been unclear. Methods: We characterized 1,280 infections from the National Basketball Association occupational health cohort identified between June 2020 and January 2022 using serial RT-qPCR testing. Logistic regression and semi-mechanistic viral RNA kinetics models were used to quantify the effect of age, variant, symptom status, infection history, vaccination status and antibody titer to the founder SARS-CoV-2 strain on the duration of potential infectiousness and overall viral kinetics. The frequency of viral rebounds was quantified under multiple cycle threshold (Ct) value-based definitions. Results: Among individuals detected partway through their infection, 51.0% (95% credible interval [CrI]: 48.3-53.6%) remained potentially infectious (Ct <30) 5 days post detection, with small differences across variants and vaccination status. Only seven viral rebounds (0.7%; N=999) were observed, with rebound defined as 3+days with Ct <30 following an initial clearance of 3+days with…

The parabrachial nucleus (PBN) is a major hub that receives sensory information from both internal and external environments. Specific populations of PBN neurons are involved in behaviors including food and water intake, nociceptive responses, breathing regulation, as well as learning and responding appropriately to threatening stimuli. However, it is unclear how many PBN neuron populations exist and how different behaviors may be encoded by unique signaling molecules or receptors. Here we provide a repository of data on the molecular identity, spatial location, and projection patterns of dozens of PBN neuron subclusters. Using single-cell RNA sequencing, we identified 21 subclusters of neurons in the PBN and neighboring regions. Multiplexed in situ hybridization showed many of these subclusters are enriched within specific PBN subregions with scattered cells in several other regions. We also provide detailed visualization of the axonal projections from 21 Cre-driver lines of mice. These results are all publicly available for download and provide a foundation for further…

Abstract Pharmacologic perturbation projects, such as Connectivity Map (CMap) and Library of Integrated Network-based Cellular Signatures (LINCS), have produced many perturbed expression data, providing enormous opportunities for computational therapeutic discovery. However, currently there is no consensus on which methodologies and parameters are the most optimal to conduct such analysis. Aiming to fill this gap, we developed new benchmarking standards for quantitatively estimating drug retrieval performance. Investigations of potential factors influencing drug retrieval were conducted based on these standards. As a result, we determined an optimal strategy for LINCS data-based therapeutic discovery. With this approach, we further identified new therapeutics for liver cancer of which the current treatment modalities remain imperfect. Both computational and experimental results demonstrated homoharringtonine (HHT) could be a promising anti-liver cancer agent. In summary, our findings will not only impact the future applications of LINCS data but also offer new opportunities for…

-linked genes. However, despite the availability of gene expression and epigenomic datasets, few variant-to-gene links have emerged. It is unclear whether these sparse results are due to limitations in available data and methods, or to deficiencies in the underlying assumed model. To better distinguish between these possibilities, we identified 220 gene-trait pairs in which protein-coding variants influence a complex trait or its Mendelian cognate. Despite the presence of expression quantitative trait loci near most GWAS associations, by applying a gene-based approach we found limited evidence that the baseline expression of trait-related genes explains GWAS associations, whether using colocalization methods (8% of genes implicated), transcription-wide association (2% of genes implicated), or a combination of regulatory annotations and distance (4% of genes implicated). These results contradict the hypothesis that most complex trait-associated variants coincide with homeostatic expression QTLs, suggesting that better models are needed. The field must confront this deficit and pursue…

Concerns about systemic racism at academic and research institutions have increased over the past decade. Here, we investigate data from the National Science Foundation (NSF), a major funder of research in the United States, and find evidence for pervasive racial disparities. In particular, white principal investigators (PIs) are consistently funded at higher rates than most non-white PIs. Funding rates for white PIs have also been increasing relative to annual overall rates with time. Moreover, disparities occur across all disciplinary directorates within the NSF and are greater for research proposals. The distributions of average external review scores also exhibit systematic offsets based on PI race. Similar patterns have been described in other research funding bodies, suggesting that racial disparities are widespread. The prevalence and persistence of these racial disparities in funding have cascading impacts that perpetuate a cumulative advantage to white PIs across all of science, technology, engineering, and mathematics.

CRISPR interference (CRISPRi) enables programmable, reversible, and titratable repression of gene expression (knockdown) in mammalian cells. Initial CRISPRi-mediated genetic screens have showcased the potential to address basic questions in cell biology, genetics, and biotechnology, but wider deployment of CRISPRi screening has been constrained by the large size of single guide RNA (sgRNA) libraries and challenges in generating cell models with consistent CRISPRi-mediated knockdown. Here, we present next-generation CRISPRi sgRNA libraries and effector expression constructs that enable strong and consistent knockdown across mammalian cell models. First, we combine empirical sgRNA selection with a dual-sgRNA library design to generate an ultra-compact (1-3 elements per gene), highly active CRISPRi sgRNA library. Next, we compare CRISPRi effectors to show that the recently published Zim3-dCas9 provides an excellent balance between strong on-target knockdown and minimal non-specific effects on cell growth or the transcriptome. Finally, we engineer a suite of cell lines with stable…

Abstract Monitoring autophagic flux is necessary for most autophagy studies. The autophagic flux assays currently available for mammalian cells are generally complicated and do not yield highly quantitative results. Yeast autophagic flux is routinely monitored with the GFP-based processing assay, whereby the amount of GFP proteolytically released from GFP-containing reporters (e.g., GFP-Atg8), detected by immunoblotting, reflects autophagic flux. However, this simple and effective assay is typically inapplicable to mammalian cells because GFP is efficiently degraded in lysosomes while the more proteolytically resistant RFP accumulates in lysosomes under basal conditions. Here, we report a HaloTag (Halo)-based reporter processing assay to monitor mammalian autophagic flux. We found that Halo is sensitive to lysosomal proteolysis but becomes resistant upon ligand binding. When delivered into lysosomes by autophagy, pulse-labeled Halo-based reporters (e.g., Halo-LC3 and Halo-GFP) are proteolytically processed to generate Halo ligand when delivered into lysosomes by autophagy. Hence,…

Abstract Recent studies have shown that gene flow or admixture has been pervasive throughout human history. While several methods exist for dating admixture in contemporary populations, they are not suitable for sparse, low coverage data available from ancient specimens. To overcome this limitation, we developed DATES that leverages ancestry covariance patterns across the genome of a single individual to infer the timing of admixture. By performing simulations, we show that DATES provides reliable results under a range of demographic scenarios and outperforms available methods for ancient DNA applications. We apply DATES to ~1,100 ancient genomes to reconstruct gene flow events during the European Holocene. Present-day Europeans derive ancestry from three distinct groups, local Mesolithic hunter-gatherers, Anatolian farmers, and Yamnaya Steppe pastoralists. These ancestral groups were themselves admixed. By studying the formation of Anatolian farmers, we infer that the gene flow related to Iranian Neolithic farmers occurred no later than 9,600 BCE, predating agriculture in Anatolia.…

The pupillary light response is an important automatic physiological response which optimizes light reaching the retina. Recent work has shown that the pupil also adjusts in response to illusory brightness and a range of cognitive functions, however, it remains unclear what exactly drives these endogenous changes. Here, we show that the imagery pupillary light response correlates with objective measures of sensory imagery strength. Further, the trial-by-trial phenomenological vividness of visual imagery is tracked by the imagery pupillary light response. We also demonstrated that a group of individuals without visual imagery (aphantasia) do not show any significant evidence of an imagery pupillary light response, however they do show perceptual pupil light responses and pupil dilation with larger cognitive load. Our results provide evidence that the pupillary light response indexes the sensory strength of visual imagery. This work also provides the first physiological validation of aphantasia.

Drosophila melanogaster olfactory neurons have long been thought to express only one chemosensory receptor gene family. There are two main olfactory receptor gene families in Drosophila , the odorant receptors (ORs) and the ionotropic receptors (IRs). The dozens of odorant-binding receptors in each family require at least one co-receptor gene in order to function: Orco for ORs, and Ir25a , Ir8a , and Ir76b for IRs. Using a new genetic knock-in strategy, we targeted the four co-receptors representing the main chemosensory families in D. melanogaster ( Orco, Ir8a, Ir76b, Ir25a ). Co-receptor knock-in expression patterns were verified as accurate representations of endogenous expression. We find extensive overlap in expression among the different co-receptors. As defined by innervation into antennal lobe glomeruli, Ir25a is broadly expressed in 88% of all olfactory sensory neuron classes and is co-expressed in 82% of Orco+ neuron classes, including all neuron classes in the maxillary palp. Orco , Ir8a , and Ir76b expression patterns are also more expansive than previously assumed.…

Two-thirds of human hormones and one-third of clinical drugs act on membrane receptors that couple to G proteins to achieve appropriate functional responses. While G protein transducers from literature are annotated in the Guide to Pharmacology database, two recent large-scale datasets now expand the receptor-G protein 'couplome'. However, these three datasets differ in scope and reported G protein couplings giving different coverage and conclusions on G protein-coupled receptor (GPCR)-G protein signaling. Here, we report a common coupling map uncovering novel couplings supported by both large-scale studies, the selectivity/promiscuity of GPCRs and G proteins, and how the co-coupling and co-expression of G proteins compare to the families from phylogenetic relationships. The coupling map and insights on GPCR-G protein selectivity will catalyze advances in receptor research and cellular signaling toward the exploitation of G protein signaling bias in design of safer drugs.

Neural activity underlying working memory is not a local phenomenon but distributed across multiple brain regions. To elucidate the circuit mechanism of such distributed activity, we developed an anatomically constrained computational model of large-scale macaque cortex. We found that mnemonic internal states may emerge from inter-areal reverberation, even in a regime where none of the isolated areas is capable of generating self-sustained activity. The mnemonic activity pattern along the cortical hierarchy indicates a transition in space, separating areas engaged in working memory and those which do not. A host of spatially distinct attractor states is found, potentially subserving various internal processes. The model yields testable predictions, including the idea of counterstream inhibitory bias, the role of prefrontal areas in controlling distributed attractors, and the resilience of distributed activity to lesions or inactivation. This work provides a theoretical framework for identifying large-scale brain mechanisms and computational principles of distributed cognitive…

Spatial partitioning is a propensity of biological systems orchestrating cell activities in space and time. The dynamic regulation of plasma membrane nano-environments has recently emerged as a key fundamental aspect of plant signaling, but the molecular components governing it are still mostly unclear. The receptor kinase FERONIA (FER) controls ligand-induced complex formation of the immune receptor kinase FLAGELLIN SENSING 2 (FLS2) with its co-receptor BRASSINOSTEROID-INSENSITIVE 1-ASSOCIATED KINASE 1 (BAK1), and perception of the endogenous peptide hormone RAPID ALKALANIZATION FACTOR 23 (RALF23) by FER inhibits immunity. Here, we show that FER regulates the plasma membrane nanoscale organization of FLS2 and BAK1. Our study demonstrates that akin to FER, leucine-rich repeat (LRR) extensin proteins (LRXs) contribute to RALF23 responsiveness and regulate BAK1 nanoscale organization and immune signaling. Furthermore, RALF23 perception leads to rapid modification of FLS2 and BAK1 nanoscale organization, and its inhibitory activity on immune signaling relies on FER kinase activity. Our…

Fixing cells with paraformaldehyde (PFA) is an essential step in numerous biological techniques as it is thought to preserve a snapshot of biomolecular transactions in living cells. Fixed-cell imaging techniques such as immunofluorescence have been widely used to detect liquid-liquid phase separation (LLPS) in vivo. Here, we compared images, before and after fixation, of cells expressing intrinsically disordered proteins that are able to undergo LLPS. Surprisingly, we found that PFA fixation can both enhance and diminish putative LLPS behaviors. For specific proteins, fixation can even cause their droplet-like puncta to artificially appear in cells that do not have any detectable puncta in the live condition. Fixing cells in the presence of glycine, a molecule that modulates fixation rates, can reverse the fixation effect from enhancing to diminishing LLPS appearance. We further established a kinetic model of fixation in the context of dynamic protein-protein interactions. Simulations based on the model suggest that protein localization in fixed cells depends on an intricate balance…

Abstract Zebrafish are an established research organism that has made many contributions to our understanding of vertebrate tissue and organ development, yet there are still significant gaps in our understanding of the genes that regulate gonad development, sex, and reproduction. Unlike the development of many organs, such as the brain and heart that form during the first few days of development, zebrafish gonads do not begin to form until the larval stage (≥5 dpf). Thus, forward genetic screens have identified very few genes required for gonad development. In addition, bulk RNA sequencing studies which identify genes expressed in the gonads do not have the resolution necessary to define minor cell populations that may play significant roles in development and function of these organs. To overcome these limitations, we have used single-cell RNA sequencing to determine the transcriptomes of cells isolated from juvenile zebrafish ovaries. This resulted in the profiles of 10,658 germ cells and 14,431 somatic cells. Our germ cell data represents all developmental stages from germline…

Genes of unknown function are among the biggest challenges in molecular biology, especially in microbial systems, where 40–60% of the predicted genes are unknown. Despite previous attempts, systematic approaches to include the unknown fraction into analytical workflows are still lacking. Here, we present a conceptual framework, its translation into the computational workflow AGNOSTOS and a demonstration on how we can bridge the known-unknown gap in genomes and metagenomes. By analyzing 415,971,742 genes predicted from 1749 metagenomes and 28,941 bacterial and archaeal genomes, we quantify the extent of the unknown fraction, its diversity, and its relevance across multiple organisms and environments. The unknown sequence space is exceptionally diverse, phylogenetically more conserved than the known fraction and predominantly taxonomically restricted at the species level. From the 71 M genes identified to be of unknown function, we compiled a collection of 283,874 lineage-specific genes of unknown function for Cand . Patescibacteria (also known as Candidate Phyla Radiation, CPR),…

Infections at the maternal-fetal interface can directly harm the fetus and induce complications that adversely impact pregnancy outcomes. Innate immune signaling by both fetal-derived placental trophoblasts and the maternal decidua must provide antimicrobial defenses at this critical interface without compromising its integrity. Here, we developed matched trophoblast (TO) and decidua organoids (DO) from human placentas to define the relative contributions of these cells to antiviral defenses at the maternal-fetal interface. We demonstrate that TO and DO basally secrete distinct immunomodulatory factors, including the constitutive release of the antiviral type III interferon IFN-λ2 from TOs, and differentially respond to viral infections through the induction of organoid-specific factors. Finally, we define the differential susceptibility and innate immune signaling of TO and DO to human cytomegalovirus (HCMV) and develop a co-culture model of TO and DO which showed that trophoblast-derived factors protect decidual cells from HCMV infection. Our findings establish matched TO and DO…

Abstract Cohesin and CTCF are major drivers of 3D genome organization, but their role in neurons is still emerging. Here we show a prominent role for cohesin in the expression of genes that facilitate neuronal maturation and homeostasis. Unexpectedly, we observed two major classes of activity-regulated genes with distinct reliance on cohesin in primary cortical neurons. Immediate early genes remained fully inducible by KCl and BDNF, and short-range enhancer-promoter contacts at the Immediate early gene Fos formed robustly in the absence of cohesin. In contrast, cohesin was required for full expression of a subset of secondary response genes characterised by long-range chromatin contacts. Cohesin-dependence of constitutive neuronal genes with key functions in synaptic transmission and neurotransmitter signaling also scaled with chromatin loop length. Our data demonstrate that key genes required for the maturation and activation of primary cortical neurons depend cohesin for their full expression, and that the degree to which these genes rely on cohesin scales with the genomic…

Abstract Learning from experience depends at least in part on changes in neuronal connections. We present the largest map of connectivity to date between cortical neurons of a defined type (L2/3 pyramidal cells), which was enabled by automated analysis of serial section electron microscopy images with improved handling of image defects. We used the map to identify constraints on the learning algorithms employed by the cortex. Previous cortical studies modeled a continuum of synapse sizes (Arellano et al. 2007) by a log-normal distribution (Loewenstein, Kuras, and Rumpel 2011; de Vivo et al. 2017; Santuy et al. 2018). A continuum is consistent with most neural network models of learning, in which synaptic strength is a continuously graded analog variable. Here we show that synapse size, when restricted to synapses between L2/3 pyramidal cells, is well-modeled by the sum of a binary variable and an analog variable drawn from a log-normal distribution. Two synapses sharing the same presynaptic and postsynaptic cells are known to be correlated in size (Sorra and Harris 1993; Koester and…

Abstract Background We present a multidisciplinary approach to clinical trial design and analysis in a pilot study ( NCT02415621 ) in which evolution-based mathematical models guide patient-specific dosing for Abiraterone treatment in men with castrate resistant metastatic prostate cancer. Methods Abiraterone plus prednisone were administered intermittently based on an evolutionary mathematical model. Outcomes are compared to historical controls and a matched contemporaneous cohort who met trial eligibility but received SOC dosing. Longitudinal cohort data allowed modification of pre-trial model parameter estimates. Model simulations of each patient using updated parameters critically evaluated trial design. Results Trial patients, on average, received no abiraterone during 59% of time on treatment. Median Time to Radiographic Progression (TTP) was 30.4 months compared to 14.3 months in the contemporaneous SOC group (p<0.001). All patients in the SOC group have progressed but 4 in the adaptive cohort remain on treatment at >1800 days. Longitudinal trial data found the competition…

Abstract Mitochondrial DNA copy number (mtDNA-CN) is an accessible blood-based measurement believed to capture underlying mitochondrial function. The specific biological processes underpinning its regulation, and whether those processes are causative for disease, is an area of active investigation. We developed a novel method for array-based mtDNA-CN estimation suitable for biobank-scale studies, called “AutoMitoC”. We applied AutoMitoC to 395,781 UKBiobank study participants and performed genome and exome-wide association studies, identifying novel common and rare genetic determinants. Overall, genetic analyses identified 73 loci for mtDNA-CN, which implicated several genes involved in rare mtDNA depletion disorders, dNTP metabolism, and the mitochondrial central dogma. Rare variant analysis identified SAMHD1 mutation carriers as having higher mtDNA-CN (beta=0.23 SDs; 95% CI, 0.18-0.29; P=2.6×10 −19 ), a potential therapeutic target for patients with mtDNA depletion disorders, but at increased risk of breast cancer (OR=1.91; 95% CI, 1.52-2.40; P=2.7×10 −8 ). Finally, Mendelian…

First recognized more than 30 years ago, glycine protects cells against rupture from diverse types of injury. This robust and widely observed effect has been speculated to target a late downstream process common to multiple modes of tissue injury. The molecular target of glycine that mediates cytoprotection, however, remains elusive. Here, we show that glycine works at the level of NINJ1, a newly identified executioner of plasma membrane rupture in pyroptosis, necrosis, and post-apoptosis lysis. NINJ1 is thought to cluster within the plasma membrane to cause cell rupture. We demonstrate that the execution of pyroptotic cell rupture is similar for human and mouse NINJ1 and that NINJ1 knockout functionally and morphologically phenocopies glycine cytoprotection in macrophages undergoing lytic cell death. Next, we show that glycine prevents NINJ1 clustering by either direct or indirect mechanisms. In pyroptosis, glycine preserves cellular integrity but does not affect upstream inflammasome activities or accompanying energetic cell death. By positioning NINJ1 clustering as a glycine…

Positional information is a central concept in developmental biology. In developing organs, positional information can be idealized as a local coordinate system that arises from morphogen gradients controlled by organizers at key locations. This offers a plausible mechanism for the integration of the molecular networks operating in individual cells into the spatially coordinated multicellular responses necessary for the organization of emergent forms. Understanding how positional cues guide morphogenesis requires the quantification of gene expression and growth dynamics in the context of their underlying coordinate systems. Here, we present recent advances in the MorphoGraphX software (Barbier de Reuille et al., 2015⁠) that implement a generalized framework to annotate developing organs with local coordinate systems. These coordinate systems introduce an organ-centric spatial context to microscopy data, allowing gene expression and growth to be quantified and compared in the context of the positional information thought to control them.

To fight the COVID-19 pandemic caused by the RNA virus SARS-CoV-2, a global vaccination campaign is in progress to achieve the immunization of billions of people mainly with adenoviral vector- or mRNA-based vaccines, all of which encode the SARS-CoV-2 Spike protein. In some rare cases, cerebral venous sinus thromboses (CVST) have been reported as a severe side effect occurring 4-14 days after the first vaccination and were often accompanied by thrombocytopenia. Besides CVST, splanchnic vein thromboses (SVT) and other thromboembolic events have been observed. These events only occurred following vaccination with adenoviral vector-based vaccines but not following vaccination with mRNA-based vaccines. Meanwhile, scientists have proposed an immune-based pathomechanism and the condition has been coined vaccine-induced immune thrombotic thrombocytopenia (VITT). Here, we describe an unexpected mechanism that could explain thromboembolic events occurring with DNA-based but not with RNA-based vaccines. We show that DNA-encoded mRNA coding for Spike protein can be spliced in a way that the…

Multiple learning processes contribute to successful goal-directed actions in the face of changing physiological states, biomechanical constraints, and environmental contexts. Amongst these processes, implicit sensorimotor adaptation is of primary importance, ensuring that movements remain well-calibrated and accurate. A large body of work on reaching movements has emphasized how adaptation centers on an iterative process designed to minimize visual errors. The role of proprioception has been largely neglected, thought to play a passive role in which proprioception is affected by the visual error but does not directly contribute to adaptation. Here, we present an alternative to this visuo-centric framework, outlining a model in which implicit adaptation acts to minimize a proprioceptive error, the distance between the perceived hand position and its intended goal. This proprioceptive re-alignment model (PReMo) is consistent with many phenomena that have previously been interpreted in terms of learning from visual errors, and offers a parsimonious account of numerous unexplained…

Abstract With the continual evolution of new strains of SARS-CoV-2 that are more virulent, transmissible, and able to evade current vaccines, there is an urgent need for effective anti-viral drugs. SARS-CoV-2 main protease (M pro ) is a leading target for drug design due to its conserved and indispensable role in the viral life cycle. Drugs targeting M pro appear promising but will elicit selection pressure for resistance. To understand resistance potential in M pro , we performed a comprehensive mutational scan of the protease that analyzed the function of all possible single amino acid changes. We developed three separate high-throughput assays of M pro function in yeast, based on either the ability of M pro variants to cleave at a defined cut-site or on the toxicity of their expression to yeast. We used deep sequencing to quantify the functional effects of each variant in each screen. The protein fitness landscapes from all three screens were strongly correlated, indicating that they captured the biophysical properties critical to M pro function. The fitness landscapes revealed a…

Background: Neuronal- and circuit-level abnormalities of excitation and inhibition are shown to be associated with tau and amyloid-beta (Aβ) in preclinical models of Alzheimer's disease (AD). These relationships remain poorly understood in patients with AD. Methods: Using empirical spectra from magnetoencephalography and computational modeling (neural mass model), we examined excitatory and inhibitory parameters of neuronal subpopulations and investigated their specific associations to regional tau and Aβ, measured by positron emission tomography, in patients with AD. Results: Patients with AD showed abnormal excitatory and inhibitory time-constants and neural gains compared to age-matched controls. Increased excitatory time-constants distinctly correlated with higher tau depositions while increased inhibitory time-constants distinctly correlated with higher Aβ depositions. Conclusions: Our results provide critical insights about potential mechanistic links between abnormal neural oscillations and cellular correlates of impaired excitatory and inhibitory synaptic functions…

Protein biomarkers have been identified across many age-related morbidities. However, characterising epigenetic influences could further inform disease predictions. Here, we leverage epigenome-wide data to study links between the DNA methylation (DNAm) signatures of the circulating proteome and incident diseases. Using data from four cohorts, we trained and tested epigenetic scores (EpiScores) for 953 plasma proteins, identifying 109 scores that explained between 1% and 58% of the variance in protein levels after adjusting for known protein quantitative trait loci (pQTL) genetic effects. By projecting these EpiScores into an independent sample (Generation Scotland; n = 9537) and relating them to incident morbidities over a follow-up of 14 years, we uncovered 137 EpiScore-disease associations. These associations were largely independent of immune cell proportions, common lifestyle and health factors, and biological aging. Notably, we found that our diabetes-associated EpiScores highlighted previous top biomarker associations from proteome-wide assessments of diabetes. These EpiScores…

Like neocortical structures, the archicortical hippocampus differs in its folding patterns across individuals. Here, we present an automated and robust BIDS-App, HippUnfold, for defining and indexing individual-specific hippocampal folding in MRI, analogous to popular tools used in neocortical reconstruction. Such tailoring is critical for inter-individual alignment, with topology serving as the basis for homology. This topological framework enables qualitatively new analyses of morphological and laminar structure in the hippocampus or its subfields. It is critical for refining current neuroimaging analyses at a meso- as well as micro-scale. HippUnfold uses state-of-the-art deep learning combined with previously developed topological constraints to generate uniquely folded surfaces to fit a given subject's hippocampal conformation. It is designed to work with commonly employed sub-millimetric MRI acquisitions, with possible extension to microscopic resolution. In this paper, we describe the power of HippUnfold in feature extraction, and highlight its unique value compared to several…

The outcome of an encounter with Mycobacterium tuberculosis ( Mtb ) depends on the pathogen’s ability to adapt to the variable immune pressures exerted by the host. Understanding this interplay has proven difficult, largely because experimentally tractable animal models do not recapitulate the heterogeneity of tuberculosis disease. We leveraged the genetically diverse Collaborative Cross (CC) mouse panel in conjunction with a library of Mtb mutants to create a resource for associating bacterial genetic requirements with host genetics and immunity. We report that CC strains vary dramatically in their susceptibility to infection and produce qualitatively distinct immune states. Global analysis of Mtb transposon mutant fitness (TnSeq) across the CC panel revealed that many virulence pathways are only required in specific host microenvironments, identifying a large fraction of the pathogen’s genome that has been maintained to ensure fitness in a diverse population. Both immunological and bacterial traits can be associated with genetic variants distributed across the mouse genome, making…

Polygenic adaptation is thought to be ubiquitous, yet remains poorly understood. Here, we model this process analytically, in the plausible setting of a highly polygenic, quantitative trait that experiences a sudden shift in the fitness optimum. We show how the mean phenotype changes over time, depending on the effect sizes of loci that contribute to variance in the trait, and characterize the allele dynamics at these loci. Notably, we describe the two phases of the allele dynamics: The first is a rapid phase, in which directional selection introduces small frequency differences between alleles whose effects are aligned with or opposed to the shift, ultimately leading to small differences in their probability of fixation during a second, longer phase, governed by stabilizing selection. As we discuss, key results should hold in more general settings and have important implications for efforts to identify the genetic basis of adaptation in humans and other species.

Abstract Reinforcement Learning (RL) models have revolutionized the cognitive and brain sciences, promising to explain behavior from simple conditioning to complex problem solving, to shed light on developmental and individual differences, and to anchor cognitive processes in specific brain mechanisms. However, the RL literature increasingly reveals contradictory results, which might cast doubt on these claims. We hypothesized that many contradictions arise from two commonly-held assumptions about computational model parameters that are actually often invalid: That parameters generalize between contexts (e.g., tasks, models) and that they capture interpretable (i.e., unique, distinctive) neurocognitive processes. To test this, we asked 291 participants aged 8-30 years to complete three learning tasks in one experimental session, and fitted RL models to each. We found that some parameters (exploration / decision noise) showed significant generalization: they followed similar developmental trajectories, and were reciprocally predictive between tasks. However, generalization was…

Mosquitoes transmit numerous pathogens, but large gaps remain in our understanding of their physiology. To facilitate explorations of mosquito biology, we have created Aegypti-Atlas (http://aegyptiatlas.buchonlab.com/), an online resource hosting RNAseq profiles of Ae. aegypti body parts (head, thorax, abdomen, gut, Malpighian tubules, ovaries), gut regions (crop, proventriculus, anterior and posterior midgut, hindgut), and a gut time course of blood meal digestion. Using Aegypti-Atlas, we provide insights into regionalization of gut function, blood feeding response, and immune defenses. We find that the anterior and posterior midgut possess digestive specializations which are preserved in the blood-fed state. Blood feeding initiates the sequential induction and repression/depletion of multiple cohorts of peptidases. With respect to defense, immune signaling components, but not recognition or effector molecules, show enrichment in ovaries. Basal expression of antimicrobial peptides is dominated by holotricin and gambicin, which are expressed in carcass and digestive tissues,…

Abstract Dendritic cells (DCs) play an important role in anti-tumor immunity by inducing T cell differentiation. Herein, we found that the mechanical sensor Piezo1 expressed by DCs integrates innate inflammatory stimuli and stiffness signals and directs the reciprocal differentiation of T H 1 and regulatory T (T reg ) cells in cancer. Genetic deletion of Piezo1 in DCs inhibited the generation of T H 1 cells while driving the development of T reg cells in promoting cancer growth. Mechanistically, Piezo1-deficient DCs regulated the secretion of the polarizing cytokines TGFβ1 and IL-12, leading to increased TGFβR2-p-Smad3 activity and decreased IL-12Rβ2-p-STAT4 activity while inducing the reciprocal differentiation of T reg and T H 1 cells. In addition, Piezo1 integrated the SIRT1-hypoxia-inducible factor-1 alpha (HIF1α)-dependent metabolic pathway and calcium-calcineurin-NFAT signaling pathway to orchestrate reciprocal T H 1 and T reg lineage commitment through DC-derived IL-12 and TGFβ1. Our studies provide critical insight for understanding the role of the DC-based mechanical…