Modeling T Cell Fate.

Many of the pathways that underlie the diversification of naive T cells into effector and memory subsets, and the maintenance of these populations, remain controversial. In recent years a variety of experimental tools have been developed that allow us to follow the fates of cells and their descendants. In this review we describe how mathematical models provide a natural language for describing the growth, loss, and differentiation of cell populations. By encoding mechanistic descriptions of cell behavior, models can help us interpret these new datasets and reveal the rules underpinning T cell fate decisions, both at steady state and during immune responses.

Understanding the cell: Future views of structural biology.

Determining the structure and mechanisms of all individual functional modules of cells at high molecular detail has often been seen as equal to understanding how cells work. Recent technical advances have led to a flush of high-resolution structures of various macromolecular machines, but despite this wealth of detailed information, our understanding of cellular function remains incomplete. Here, we discuss present-day limitations of structural biology and highlight novel technologies that may enable us to analyze molecular functions directly inside cells. We predict that the progression toward structural cell biology will involve a shift toward conceptualizing a 4D virtual reality of cells using digital twins. These will capture cellular segments in a highly enriched molecular detail, include dynamic changes, and facilitate simulations of molecular processes, leading to novel and experimentally testable predictions. Transferring biological questions into algorithms that learn from the existing wealth of data and explore novel solutions may ultimately unveil how cells work.

Systematic mapping of organism-scale gene-regulatory networks in aging using population asynchrony.

In aging, physiologic networks decline in function at rates that differ between individuals, producing a wide distribution of lifespan. Though 70% of human lifespan variance remains unexplained by heritable factors, little is known about the intrinsic sources of physiologic heterogeneity in aging. To understand how complex physiologic networks generate lifespan variation, new methods are needed. Here, we present Asynch-seq, an approach that uses gene-expression heterogeneity within isogenic populations to study the processes generating lifespan variation. By collecting thousands of single-individual transcriptomes, we capture the Caenorhabditis elegans "pan-transcriptome"-a highly resolved atlas of non-genetic variation. We use our atlas to guide a large-scale perturbation screen that identifies the decoupling of total mRNA content between germline and soma as the largest source of physiologic heterogeneity in aging, driven by pleiotropic genes whose knockdown dramatically reduces lifespan variance. Our work demonstrates how systematic mapping of physiologic heterogeneity can be applied to reduce inter-individual disparities in aging.

Human fetal brain self-organizes into long-term expanding organoids.

Human brain development involves an orchestrated, massive neural progenitor expansion while a multi-cellular tissue architecture is established. Continuously expanding organoids can be grown directly from multiple somatic tissues, yet to date, brain organoids can solely be established from pluripotent stem cells. Here, we show that healthy human fetal brain in vitro self-organizes into organoids (FeBOs), phenocopying aspects of in vivo cellular heterogeneity and complex organization. FeBOs can be expanded over long time periods. FeBO growth requires maintenance of tissue integrity, which ensures production of a tissue-like extracellular matrix (ECM) niche, ultimately endowing FeBO expansion. FeBO lines derived from different areas of the central nervous system (CNS), including dorsal and ventral forebrain, preserve their regional identity and allow to probe aspects of positional identity. Using CRISPR-Cas9, we showcase the generation of syngeneic mutant FeBO lines for the study of brain cancer. Taken together, FeBOs constitute a complementary CNS organoid platform.

Metabolic trade-offs constrain the cell size ratio in a nitrogen-fixing symbiosis.

Biological dinitrogen (N2) fixation is a key metabolic process exclusively performed by prokaryotes, some of which are symbiotic with eukaryotes. Species of the marine haptophyte algae Braarudosphaera bigelowii harbor the N2-fixing endosymbiotic cyanobacteria UCYN-A, which might be evolving organelle-like characteristics. We found that the size ratio between UCYN-A and their hosts is strikingly conserved across sublineages/species, which is consistent with the size relationships of organelles in this symbiosis and other species. Metabolic modeling showed that this size relationship maximizes the coordinated growth rate based on trade-offs between resource acquisition and exchange. Our findings show that the size relationships of N2-fixing endosymbionts and organelles in unicellular eukaryotes are constrained by predictable metabolic underpinnings and that UCYN-A is, in many regards, functioning like a hypothetical N2-fixing organelle (or nitroplast).

Fate Mapping and Quantitation of Hematopoiesis In Vivo.

Hematopoietic stem cells (HSCs) and downstream progenitors have long been studied based on phenotype, cell purification, proliferation, and transplantation into myeloablated recipients. These experiments, complemented by data on expression profiles, mouse mutants, and humans with hematopoietic defects, are the foundation for the current hematopoietic differentiation tree. However, there are fundamental gaps in our knowledge of the quantitative and qualitative operation of the HSC/progenitor system under physiological and pathological conditions in vivo. The hallmarks of HSCs, self-renewal and multipotency, are observed in in vitro assays and cell transplantation experiments; however, the extent to which these features occur naturally in HSCs and progenitors remains uncertain. We focus here on work that strives to address these unresolved questions, with emphasis on fate mapping and modeling of the hematopoietic flow from stem cells toward myeloid and lymphoid lineages during development and adult life.

High-resolution analyses of associations between medications, microbiome, and mortality in cancer patients.

Discerning the effect of pharmacological exposures on intestinal bacterial communities in cancer patients is challenging. Here, we deconvoluted the relationship between drug exposures and changes in microbial composition by developing and applying a new computational method, PARADIGM (parameters associated with dynamics of gut microbiota), to a large set of longitudinal fecal microbiome profiles with detailed medication-administration records from patients undergoing allogeneic hematopoietic cell transplantation. We observed that several non-antibiotic drugs, including laxatives, antiemetics, and opioids, are associated with increased Enterococcus relative abundance and decreased alpha diversity. Shotgun metagenomic sequencing further demonstrated subspecies competition, leading to increased dominant-strain genetic convergence during allo-HCT that is significantly associated with antibiotic exposures. We integrated drug-microbiome associations to predict clinical outcomes in two validation cohorts on the basis of drug exposures alone, suggesting that this approach can generate biologically and clinically relevant insights into how pharmacological exposures can perturb or preserve microbiota composition. The application of a computational method called PARADIGM to a large dataset of cancer patients' longitudinal fecal specimens and detailed daily medication records reveals associations between drug exposures and the intestinal microbiota that recapitulate in vitro findings and are also predictive of clinical outcomes.

Cooperative assembly confers regulatory specificity and long-term genetic circuit stability.

A ubiquitous feature of eukaryotic transcriptional regulation is cooperative self-assembly between transcription factors (TFs) and DNA cis-regulatory motifs. It is thought that this strategy enables specific regulatory connections to be formed in gene networks between otherwise weakly interacting, low-specificity molecular components. Here, using synthetic gene circuits constructed in yeast, we find that high regulatory specificity can emerge from cooperative, multivalent interactions among artificial zinc-finger-based TFs. We show that circuits "wired" using the strategy of cooperative TF assembly are effectively insulated from aberrant misregulation of the host cell genome. As we demonstrate in experiments and mathematical models, this mechanism is sufficient to rescue circuit-driven fitness defects, resulting in genetic and functional stability of circuits in long-term continuous culture. Our naturally inspired approach offers a simple, generalizable means for building high-fidelity, evolutionarily robust gene circuits that can be scaled to a wide range of host organisms and applications.

De novo protein identification in mammalian sperm using in situ cryoelectron tomography and AlphaFold2 docking.

To understand the molecular mechanisms of cellular pathways, contemporary workflows typically require multiple techniques to identify proteins, track their localization, and determine their structures in vitro. Here, we combined cellular cryoelectron tomography (cryo-ET) and AlphaFold2 modeling to address these questions and understand how mammalian sperm are built in situ. Our cellular cryo-ET and subtomogram averaging provided 6.0-Å reconstructions of axonemal microtubule structures. The well-resolved tertiary structures allowed us to unbiasedly match sperm-specific densities with 21,615 AlphaFold2-predicted protein models of the mouse proteome. We identified Tektin 5, CCDC105, and SPACA9 as novel microtubule-associated proteins. These proteins form an extensive interaction network crosslinking the lumen of axonemal doublet microtubules, suggesting their roles in modulating the mechanical properties of the filaments. Indeed, Tekt5 -/- sperm possess more deformed flagella with 180° bends. Together, our studies presented a cellular visual proteomics workflow and shed light on the in vivo functions of Tektin 5.

A transcription factor atlas of directed differentiation.

Transcription factors (TFs) regulate gene programs, thereby controlling diverse cellular processes and cell states. To comprehensively understand TFs and the programs they control, we created a barcoded library of all annotated human TF splice isoforms (>3,500) and applied it to build a TF Atlas charting expression profiles of human embryonic stem cells (hESCs) overexpressing each TF at single-cell resolution. We mapped TF-induced expression profiles to reference cell types and validated candidate TFs for generation of diverse cell types, spanning all three germ layers and trophoblasts. Targeted screens with subsets of the library allowed us to create a tailored cellular disease model and integrate mRNA expression and chromatin accessibility data to identify downstream regulators. Finally, we characterized the effects of combinatorial TF overexpression by developing and validating a strategy for predicting combinations of TFs that produce target expression profiles matching reference cell types to accelerate cellular engineering efforts.

Cholesterol and matrisome pathways dysregulated in astrocytes and microglia.

The impact of apolipoprotein E ε4 (APOE4), the strongest genetic risk factor for Alzheimer's disease (AD), on human brain cellular function remains unclear. Here, we investigated the effects of APOE4 on brain cell types derived from population and isogenic human induced pluripotent stem cells, post-mortem brain, and APOE targeted replacement mice. Population and isogenic models demonstrate that APOE4 local haplotype, rather than a single risk allele, contributes to risk. Global transcriptomic analyses reveal human-specific, APOE4-driven lipid metabolic dysregulation in astrocytes and microglia. APOE4 enhances de novo cholesterol synthesis despite elevated intracellular cholesterol due to lysosomal cholesterol sequestration in astrocytes. Further, matrisome dysregulation is associated with upregulated chemotaxis, glial activation, and lipid biosynthesis in astrocytes co-cultured with neurons, which recapitulates altered astrocyte matrisome signaling in human brain. Thus, APOE4 initiates glia-specific cell and non-cell autonomous dysregulation that may contribute to increased AD risk.

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.

Biosensors based on peptide exposure show single molecule conformations in live cells.

We describe an approach to study the conformation of individual proteins during single particle tracking (SPT) in living cells. "Binder/tag" is based on incorporation of a 7-mer peptide (the tag) into a protein where its solvent exposure is controlled by protein conformation. Only upon exposure can the peptide specifically interact with a reporter protein (the binder). Thus, simple fluorescence localization reflects protein conformation. Through direct excitation of bright dyes, the trajectory and conformation of individual proteins can be followed. Simple protein engineering provides highly specific biosensors suitable for SPT and FRET. We describe tagSrc, tagFyn, tagSyk, tagFAK, and an orthogonal binder/tag pair. SPT showed slowly diffusing islands of activated Src within Src clusters and dynamics of activation in adhesions. Quantitative analysis and stochastic modeling revealed in vivo Src kinetics. The simplicity of binder/tag can provide access to diverse proteins.

Meiotic recombination mirrors patterns of germline replication in mice and humans.

Genetic recombination generates novel trait combinations, and understanding how recombination is distributed across the genome is key to modern genetics. The PRDM9 protein defines recombination hotspots; however, megabase-scale recombination patterning is independent of PRDM9. The single round of DNA replication, which precedes recombination in meiosis, may establish these patterns; therefore, we devised an approach to study meiotic replication that includes robust and sensitive mapping of replication origins. We find that meiotic DNA replication is distinct; reduced origin firing slows replication in meiosis, and a distinctive replication pattern in human males underlies the subtelomeric increase in recombination. We detected a robust correlation between replication and both contemporary and historical recombination and found that replication origin density coupled with chromosome size determines the recombination potential of individual chromosomes. Our findings and methods have implications for understanding the mechanisms underlying DNA replication, genetic recombination, and the landscape of mammalian germline variation.

Time-resolved systems immunology reveals a late juncture linked to fatal COVID-19.

COVID-19 exhibits extensive patient-to-patient heterogeneity. To link immune response variation to disease severity and outcome over time, we longitudinally assessed circulating proteins as well as 188 surface protein markers, transcriptome, and T cell receptor sequence simultaneously in single peripheral immune cells from COVID-19 patients. Conditional-independence network analysis revealed primary correlates of disease severity, including gene expression signatures of apoptosis in plasmacytoid dendritic cells and attenuated inflammation but increased fatty acid metabolism in CD56dimCD16hi NK cells linked positively to circulating interleukin (IL)-15. CD8+ T cell activation was apparent without signs of exhaustion. Although cellular inflammation was depressed in severe patients early after hospitalization, it became elevated by days 17-23 post symptom onset, suggestive of a late wave of inflammatory responses. Furthermore, circulating protein trajectories at this time were divergent between and predictive of recovery versus fatal outcomes. Our findings stress the importance of timing in the analysis, clinical monitoring, and therapeutic intervention of COVID-19.

Changes in genome architecture and transcriptional dynamics progress independently of sensory experience during post-natal brain development.

Both transcription and three-dimensional (3D) architecture of the mammalian genome play critical roles in neurodevelopment and its disorders. However, 3D genome structures of single brain cells have not been solved; little is known about the dynamics of single-cell transcriptome and 3D genome after birth. Here, we generated a transcriptome (3,517 cells) and 3D genome (3,646 cells) atlas of the developing mouse cortex and hippocampus by using our high-resolution multiple annealing and looping-based amplification cycles for digital transcriptomics (MALBAC-DT) and diploid chromatin conformation capture (Dip-C) methods and developing multi-omic analysis pipelines. In adults, 3D genome "structure types" delineate all major cell types, with high correlation between chromatin A/B compartments and gene expression. During development, both transcriptome and 3D genome are extensively transformed in the first post-natal month. In neurons, 3D genome is rewired across scales, correlated with gene expression modules, and independent of sensory experience. Finally, we examine allele-specific structure of imprinted genes, revealing local and chromosome (chr)-wide differences. These findings uncover an unknown dimension of neurodevelopment.

Metabolic modeling of single Th17 cells reveals regulators of autoimmunity.

Metabolism is a major regulator of immune cell function, but it remains difficult to study the metabolic status of individual cells. Here, we present Compass, an algorithm to characterize cellular metabolic states based on single-cell RNA sequencing and flux balance analysis. We applied Compass to associate metabolic states with T helper 17 (Th17) functional variability (pathogenic potential) and recovered a metabolic switch between glycolysis and fatty acid oxidation, akin to known Th17/regulatory T cell (Treg) differences, which we validated by metabolic assays. Compass also predicted that Th17 pathogenicity was associated with arginine and downstream polyamine metabolism. Indeed, polyamine-related enzyme expression was enhanced in pathogenic Th17 and suppressed in Treg cells. Chemical and genetic perturbation of polyamine metabolism inhibited Th17 cytokines, promoted Foxp3 expression, and remodeled the transcriptome and epigenome of Th17 cells toward a Treg-like state. In vivo perturbations of the polyamine pathway altered the phenotype of encephalitogenic T cells and attenuated tissue inflammation in CNS autoimmunity.

A local regulatory T cell feedback circuit maintains immune homeostasis by pruning self-activated T cells.

A fraction of mature T cells can be activated by peripheral self-antigens, potentially eliciting host autoimmunity. We investigated homeostatic control of self-activated T cells within unperturbed tissue environments by combining high-resolution multiplexed and volumetric imaging with computational modeling. In lymph nodes, self-activated T cells produced interleukin (IL)-2, which enhanced local regulatory T cell (Treg) proliferation and inhibitory functionality. The resulting micro-domains reciprocally constrained inputs required for damaging effector responses, including CD28 co-stimulation and IL-2 signaling, constituting a negative feedback circuit. Due to these local constraints, self-activated T cells underwent transient clonal expansion, followed by rapid death ("pruning"). Computational simulations and experimental manipulations revealed the feedback machinery's quantitative limits: modest reductions in Treg micro-domain density or functionality produced non-linear breakdowns in control, enabling self-activated T cells to subvert pruning. This fine-tuned, paracrine feedback process not only enforces immune homeostasis but also establishes a sharp boundary between autoimmune and host-protective T cell responses.

A Patient-Derived Glioblastoma Organoid Model and Biobank Recapitulates Inter- and Intra-tumoral Heterogeneity.

Glioblastomas exhibit vast inter- and intra-tumoral heterogeneity, complicating the development of effective therapeutic strategies. Current in vitro models are limited in preserving the cellular and mutational diversity of parental tumors and require a prolonged generation time. Here, we report methods for generating and biobanking patient-derived glioblastoma organoids (GBOs) that recapitulate the histological features, cellular diversity, gene expression, and mutational profiles of their corresponding parental tumors. GBOs can be generated quickly with high reliability and exhibit rapid, aggressive infiltration when transplanted into adult rodent brains. We further demonstrate the utility of GBOs to test personalized therapies by correlating GBO mutational profiles with responses to specific drugs and by modeling chimeric antigen receptor T cell immunotherapy. Our studies show that GBOs maintain many key features of glioblastomas and can be rapidly deployed to investigate patient-specific treatment strategies. Additionally, our live biobank establishes a rich resource for basic and translational glioblastoma research.

The In Situ Structure of Parkinson's Disease-Linked LRRK2.

Mutations in leucine-rich repeat kinase 2 (LRRK2) are the most frequent cause of familial Parkinson's disease. LRRK2 is a multi-domain protein containing a kinase and GTPase. Using correlative light and electron microscopy, in situ cryo-electron tomography, and subtomogram analysis, we reveal a 14-Å structure of LRRK2 bearing a pathogenic mutation that oligomerizes as a right-handed double helix around microtubules, which are left-handed. Using integrative modeling, we determine the architecture of LRRK2, showing that the GTPase and kinase are in close proximity, with the GTPase closer to the microtubule surface, whereas the kinase is exposed to the cytoplasm. We identify two oligomerization interfaces mediated by non-catalytic domains. Mutation of one of these abolishes LRRK2 microtubule-association. Our work demonstrates the power of cryo-electron tomography to generate models of previously unsolved structures in their cellular environment.