Unraveling the Complex Interplay Between T Cell Metabolism and Function.

Metabolism drives function, on both an organismal and a cellular level. In T cell biology, metabolic remodeling is intrinsically linked to cellular development, activation, function, differentiation, and survival. After naive T cells are activated, increased demands for metabolic currency in the form of ATP, as well as biomass for cell growth, proliferation, and the production of effector molecules, are met by rewiring cellular metabolism. Consequently, pharmacological strategies are being developed to perturb or enhance selective metabolic processes that are skewed in immune-related pathologies. Here we review the most recent advances describing the metabolic changes that occur during the T cell lifecycle. We discuss how T cell metabolism can have profound effects on health and disease and where it might be a promising target to treat a variety of pathologies.

Human T cell generation is restored in CD3δ severe combined immunodeficiency through adenine base editing.

CD3δ SCID is a devastating inborn error of immunity caused by mutations in CD3D, encoding the invariant CD3δ chain of the CD3/TCR complex necessary for normal thymopoiesis. We demonstrate an adenine base editing (ABE) strategy to restore CD3δ in autologous hematopoietic stem and progenitor cells (HSPCs). Delivery of mRNA encoding a laboratory-evolved ABE and guide RNA into a CD3δ SCID patient's HSPCs resulted in a 71.2% ± 7.85% (n = 3) correction of the pathogenic mutation. Edited HSPCs differentiated in artificial thymic organoids produced mature T cells exhibiting diverse TCR repertoires and TCR-dependent functions. Edited human HSPCs transplanted into immunodeficient mice showed 88% reversion of the CD3D defect in human CD34+ cells isolated from mouse bone marrow after 16 weeks, indicating correction of long-term repopulating HSCs. These findings demonstrate the preclinical efficacy of ABE in HSPCs for the treatment of CD3δ SCID, providing a foundation for the development of a one-time treatment for CD3δ SCID patients.

A cross-disorder dosage sensitivity map of the human genome.

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

Phosphoinositide acyl chain saturation drives CD8+ effector T cell signaling and function.

How lipidome changes support CD8+ effector T (Teff) cell differentiation is not well understood. Here we show that, although naive T cells are rich in polyunsaturated phosphoinositides (PIPn with 3-4 double bonds), Teff cells have unique PIPn marked by saturated fatty acyl chains (0-2 double bonds). PIPn are precursors for second messengers. Polyunsaturated phosphatidylinositol bisphosphate (PIP2) exclusively supported signaling immediately upon T cell antigen receptor activation. In late Teff cells, activity of phospholipase C-γ1, the enzyme that cleaves PIP2 into downstream mediators, waned, and saturated PIPn became essential for sustained signaling. Saturated PIP was more rapidly converted to PIP2 with subsequent recruitment of phospholipase C-γ1, and loss of saturated PIPn impaired Teff cell fitness and function, even in cells with abundant polyunsaturated PIPn. Glucose was the substrate for de novo PIPn synthesis, and was rapidly utilized for saturated PIP2 generation. Thus, separate PIPn pools with distinct acyl chain compositions and metabolic dependencies drive important signaling events to initiate and then sustain effector function during CD8+ T cell differentiation.

Bacteria-to-Human Protein Networks Reveal Origins of Endogenous DNA Damage.

DNA damage provokes mutations and cancer and results from external carcinogens or endogenous cellular processes. However, the intrinsic instigators of endogenous DNA damage are poorly understood. Here, we identify proteins that promote endogenous DNA damage when overproduced: the DNA "damage-up" proteins (DDPs). We discover a large network of DDPs in Escherichia coli and deconvolute them into six function clusters, demonstrating DDP mechanisms in three: reactive oxygen increase by transmembrane transporters, chromosome loss by replisome binding, and replication stalling by transcription factors. Their 284 human homologs are over-represented among known cancer drivers, and their RNAs in tumors predict heavy mutagenesis and a poor prognosis. Half of the tested human homologs promote DNA damage and mutation when overproduced in human cells, with DNA damage-elevating mechanisms like those in E. coli. Our work identifies networks of DDPs that provoke endogenous DNA damage and may reveal DNA damage-associated functions of many human known and newly implicated cancer-promoting proteins.

Dissecting the Causal Mechanism of X-Linked Dystonia-Parkinsonism by Integrating Genome and Transcriptome Assembly.

X-linked Dystonia-Parkinsonism (XDP) is a Mendelian neurodegenerative disease that is endemic to the Philippines and is associated with a founder haplotype. We integrated multiple genome and transcriptome assembly technologies to narrow the causal mutation to the TAF1 locus, which included a SINE-VNTR-Alu (SVA) retrotransposition into intron 32 of the gene. Transcriptome analyses identified decreased expression of the canonical cTAF1 transcript among XDP probands, and de novo assembly across multiple pluripotent stem-cell-derived neuronal lineages discovered aberrant TAF1 transcription that involved alternative splicing and intron retention (IR) in proximity to the SVA that was anti-correlated with overall TAF1 expression. CRISPR/Cas9 excision of the SVA rescued this XDP-specific transcriptional signature and normalized TAF1 expression in probands. These data suggest an SVA-mediated aberrant transcriptional mechanism associated with XDP and may provide a roadmap for layered technologies and integrated assembly-based analyses for other unsolved Mendelian disorders.

Mitochondrial Priming by CD28.

T cell receptor (TCR) signaling without CD28 can elicit primary effector T cells, but memory T cells generated during this process are anergic, failing to respond to secondary antigen exposure. We show that, upon T cell activation, CD28 transiently promotes expression of carnitine palmitoyltransferase 1a (Cpt1a), an enzyme that facilitates mitochondrial fatty acid oxidation (FAO), before the first cell division, coinciding with mitochondrial elongation and enhanced spare respiratory capacity (SRC). microRNA-33 (miR33), a target of thioredoxin-interacting protein (TXNIP), attenuates Cpt1a expression in the absence of CD28, resulting in cells that thereafter are metabolically compromised during reactivation or periods of increased bioenergetic demand. Early CD28-dependent mitochondrial engagement is needed for T cells to remodel cristae, develop SRC, and rapidly produce cytokines upon restimulation-cardinal features of protective memory T cells. Our data show that initial CD28 signals during T cell activation prime mitochondria with latent metabolic capacity that is essential for future T cell responses.

Inflammatory macrophage dependence on NAD+ salvage is a consequence of reactive oxygen species-mediated DNA damage.

The adoption of Warburg metabolism is critical for the activation of macrophages in response to lipopolysaccharide. Macrophages stimulated with lipopolysaccharide increase their expression of nicotinamide phosphoribosyltransferase (NAMPT), a key enzyme in NAD+ salvage, and loss of NAMPT activity alters their inflammatory potential. However, the events that lead to the cells' becoming dependent on NAD+ salvage remain poorly defined. We found that depletion of NAD+ and increased expression of NAMPT occurred rapidly after inflammatory activation and coincided with DNA damage caused by reactive oxygen species (ROS). ROS produced by complex III of the mitochondrial electron-transport chain were required for macrophage activation. DNA damage was associated with activation of poly(ADP-ribose) polymerase, which led to consumption of NAD+. In this setting, increased NAMPT expression allowed the maintenance of NAD+ pools sufficient for glyceraldehyde-3-phosphate dehydrogenase activity and Warburg metabolism. Our findings provide an integrated explanation for the dependence of inflammatory macrophages on the NAD+ salvage pathway.

Matrix viscoelasticity promotes liver cancer progression in the pre-cirrhotic liver.

Type 2 diabetes mellitus is a major risk factor for hepatocellular carcinoma (HCC). Changes in extracellular matrix (ECM) mechanics contribute to cancer development1,2, and increased stiffness is known to promote HCC progression in cirrhotic conditions3,4. Type 2 diabetes mellitus is characterized by an accumulation of advanced glycation end-products (AGEs) in the ECM; however, how this affects HCC in non-cirrhotic conditions is unclear. Here we find that, in patients and animal models, AGEs promote changes in collagen architecture and enhance ECM viscoelasticity, with greater viscous dissipation and faster stress relaxation, but not changes in stiffness. High AGEs and viscoelasticity combined with oncogenic ß-catenin signalling promote HCC induction, whereas inhibiting AGE production, reconstituting the AGE clearance receptor AGER1 or breaking AGE-mediated collagen cross-links reduces viscoelasticity and HCC growth. Matrix analysis and computational modelling demonstrate that lower interconnectivity of AGE-bundled collagen matrix, marked by shorter fibre length and greater heterogeneity, enhances viscoelasticity. Mechanistically, animal studies and 3D cell cultures show that enhanced viscoelasticity promotes HCC cell proliferation and invasion through an integrin-ß1-tensin-1-YAP mechanotransductive pathway. These results reveal that AGE-mediated structural changes enhance ECM viscoelasticity, and that viscoelasticity can promote cancer progression in vivo, independent of stiffness.

A genomic mutational constraint map using variation in 76,156 human genomes.

The depletion of disruptive variation caused by purifying natural selection (constraint) has been widely used to investigate protein-coding genes underlying human disorders1-4, but attempts to assess constraint for non-protein-coding regions have proved more difficult. Here we aggregate, process and release a dataset of 76,156 human genomes from the Genome Aggregation Database (gnomAD)-the largest public open-access human genome allele frequency reference dataset-and use it to build a genomic constraint map for the whole genome (genomic non-coding constraint of haploinsufficient variation (Gnocchi)). We present a refined mutational model that incorporates local sequence context and regional genomic features to detect depletions of variation. As expected, the average constraint for protein-coding sequences is stronger than that for non-coding regions. Within the non-coding genome, constrained regions are enriched for known regulatory elements and variants that are implicated in complex human diseases and traits, facilitating the triangulation of biological annotation, disease association and natural selection to non-coding DNA analysis. More constrained regulatory elements tend to regulate more constrained protein-coding genes, which in turn suggests that non-coding constraint can aid the identification of constrained genes that are as yet unrecognized by current gene constraint metrics. We demonstrate that this genome-wide constraint map improves the identification and interpretation of functional human genetic variation.

Quantum textures of the many-body wavefunctions in magic-angle graphene.

Interactions among electrons create novel many-body quantum phases of matter with wavefunctions that reflect electronic correlation effects, broken symmetries and collective excitations. Many quantum phases have been discovered in magic-angle twisted bilayer graphene (MATBG), including correlated insulating1, unconventional superconducting2-5 and magnetic topological6-9 phases. The lack of microscopic information10,11 of possible broken symmetries has hampered our understanding of these phases12-17. Here we use high-resolution scanning tunnelling microscopy to study the wavefunctions of the correlated phases in MATBG. The squares of the wavefunctions of gapped phases, including those of the correlated insulating, pseudogap and superconducting phases, show distinct broken-symmetry patterns with a √3 × âˆš3 super-periodicity on the graphene atomic lattice that has a complex spatial dependence on the moiré scale. We introduce a symmetry-based analysis using a set of complex-valued local order parameters, which show intricate textures that distinguish the various correlated phases. We compare the observed quantum textures of the correlated insulators at fillings of ±2 electrons per moiré unit cell to those expected for proposed theoretical ground states. In typical MATBG devices, these textures closely match those of the proposed incommensurate Kekulé spiral order15, whereas in ultralow-strain samples, our data have local symmetries like those of a time-reversal symmetric intervalley coherent phase12. Moreover, the superconducting state of MATBG shows strong signatures of intervalley coherence, only distinguishable from those of the insulator with our phase-sensitive measurements.

Inhibiting membrane rupture with NINJ1 antibodies limits tissue injury.

Plasma membrane rupture (PMR) in dying cells undergoing pyroptosis or apoptosis requires the cell-surface protein NINJ11. PMR releases pro-inflammatory cytoplasmic molecules, collectively called damage-associated molecular patterns (DAMPs), that activate immune cells. Therefore, inhibiting NINJ1 and PMR may limit the inflammation that is associated with excessive cell death. Here we describe an anti-NINJ1 monoclonal antibody that specifically targets mouse NINJ1 and blocks oligomerization of NINJ1, preventing PMR. Electron microscopy studies showed that this antibody prevents NINJ1 from forming oligomeric filaments. In mice, inhibition of NINJ1 or Ninj1 deficiency ameliorated hepatocellular PMR induced with TNF plus D-galactosamine, concanavalin A, Jo2 anti-Fas agonist antibody or ischaemia-reperfusion injury. Accordingly, serum levels of lactate dehydrogenase, the liver enzymes alanine aminotransaminase and aspartate aminotransferase, and the DAMPs interleukin 18 and HMGB1 were reduced. Moreover, in the liver ischaemia-reperfusion injury model, there was an attendant reduction in neutrophil infiltration. These data indicate that NINJ1 mediates PMR and inflammation in diseases driven by aberrant hepatocellular death.

Entwined African and Asian genetic roots of medieval peoples of the Swahili coast.

The urban peoples of the Swahili coast traded across eastern Africa and the Indian Ocean and were among the first practitioners of Islam among sub-Saharan people1,2. The extent to which these early interactions between Africans and non-Africans were accompanied by genetic exchange remains unknown. Here we report ancient DNA data for 80 individuals from 6 medieval and early modern (AD 1250-1800) coastal towns and an inland town after AD 1650. More than half of the DNA of many of the individuals from coastal towns originates from primarily female ancestors from Africa, with a large proportion-and occasionally more than half-of the DNA coming from Asian ancestors. The Asian ancestry includes components associated with Persia and India, with 80-90% of the Asian DNA originating from Persian men. Peoples of African and Asian origins began to mix by about AD 1000, coinciding with the large-scale adoption of Islam. Before about AD 1500, the Southwest Asian ancestry was mainly Persian-related, consistent with the narrative of the Kilwa Chronicle, the oldest history told by people of the Swahili coast3. After this time, the sources of DNA became increasingly Arabian, consistent with evidence of growing interactions with southern Arabia4. Subsequent interactions with Asian and African people further changed the ancestry of present-day people of the Swahili coast in relation to the medieval individuals whose DNA we sequenced.

Mechanical integrated circuit materials.

Recent developments in autonomous engineered matter have introduced the ability for intelligent materials to process environmental stimuli and functionally adapt1-4. To formulate a foundation for such an engineered living material paradigm, researchers have introduced sensing5-11 and actuating12-16 functionalities in soft matter. Yet, information processing is the key functional element of autonomous engineered matter that has been recently explored through unconventional techniques with limited computing scalability17-20. Here we uncover a relation between Boolean mathematics and kinematically reconfigurable electrical circuits to realize all combinational logic operations in soft, conductive mechanical materials. We establish an analytical framework that minimizes the canonical functions of combinational logic by the Quine-McCluskey method, and governs the mechanical design of reconfigurable integrated circuit switching networks in soft matter. The resulting mechanical integrated circuit materials perform higher-level arithmetic, number comparison, and decode binary data to visual representations. We exemplify two methods to automate the design on the basis of canonical Boolean functions and individual gate-switching assemblies. We also increase the computational density of the materials by a monolithic layer-by-layer design approach. As the framework established here leverages mathematics and kinematics for system design, the proposed approach of mechanical integrated circuit materials can be realized on any length scale and in a wide variety of physics.

Two-Parameter Mobility Assessments Discriminate Diverse Regulatory Factor Behaviors in Chromatin.

Enzymatic probes of chromatin structure reveal accessible versus inaccessible chromatin states, while super-resolution microscopy reveals a continuum of chromatin compaction states. Characterizing histone H2B movements by single-molecule tracking (SMT), we resolved chromatin domains ranging from low to high mobility and displaying different subnuclear localizations patterns. Heterochromatin constituents correlated with the lowest mobility chromatin, whereas transcription factors varied widely with regard to their respective mobility with low- or high-mobility chromatin. Pioneer transcription factors, which bind nucleosomes, can access the low-mobility chromatin domains, whereas weak or non-nucleosome binding factors are excluded from the domains and enriched in higher mobility domains. Nonspecific DNA and nucleosome binding accounted for most of the low mobility of strong nucleosome interactor FOXA1. Our analysis shows how the parameters of the mobility of chromatin-bound factors, but not their diffusion behaviors or SMT-residence times within chromatin, distinguish functional characteristics of different chromatin-interacting proteins.

Diverse heterochromatin states restricting cell identity and reprogramming.

Heterochromatin is defined as a chromosomal domain harboring repressive H3K9me2/3 or H3K27me3 histone modifications and relevant factors that physically compact the chromatin. Heterochromatin can restrict where transcription factors bind, providing a barrier to gene activation and changes in cell identity. While heterochromatin thus helps maintain cell differentiation, it presents a barrier to overcome during efforts to reprogram cells for biomedical purposes. Recent findings have revealed complexity in the composition and regulation of heterochromatin, and shown that transiently disrupting the machinery of heterochromatin can enhance reprogramming. Here, we discuss how heterochromatin is established and maintained during development, and how our growing understanding of the mechanisms regulating H3K9me3 heterochromatin can be leveraged to improve our ability to direct changes in cell identity.

Evidence for unconventional superconductivity in twisted bilayer graphene.

The emergence of superconductivity and correlated insulators in magic-angle twisted bilayer graphene (MATBG) has raised the intriguing possibility that its pairing mechanism is distinct from that of conventional superconductors1-4, as described by the Bardeen-Cooper-Schrieffer (BCS) theory. However, recent studies have shown that superconductivity persists even when Coulomb interactions are partially screened5,6. This suggests that pairing in MATBG might be conventional in nature and a consequence of the large density of states of its flat bands. Here we combine tunnelling and Andreev reflection spectroscopy with a scanning tunnelling microscope to observe several key experimental signatures of unconventional superconductivity in MATBG. We show that the tunnelling spectra below the transition temperature Tc are inconsistent with those of a conventional s-wave superconductor, but rather resemble those of a nodal superconductor with an anisotropic pairing mechanism. We observe a large discrepancy between the tunnelling gap ΔT, which far exceeds the mean-field BCS ratio (with 2ΔT/kBTc ~ 25), and the gap ΔAR extracted from Andreev reflection spectroscopy (2ΔAR/kBTc ~ 6). The tunnelling gap persists even when superconductivity is suppressed, indicating its emergence from a pseudogap phase. Moreover, the pseudogap and superconductivity are both absent when MATBG is aligned with hexagonal boron nitride. These findings and other observations reported here provide a preponderance of evidence for a non-BCS mechanism for superconductivity in MATBG.

Genome-wide enhancer maps link risk variants to disease genes.

Genome-wide association studies (GWAS) have identified thousands of noncoding loci that are associated with human diseases and complex traits, each of which could reveal insights into the mechanisms of disease1. Many of the underlying causal variants may affect enhancers2,3, but we lack accurate maps of enhancers and their target genes to interpret such variants. We recently developed the activity-by-contact (ABC) model to predict which enhancers regulate which genes and validated the model using CRISPR perturbations in several cell types4. Here we apply this ABC model to create enhancer-gene maps in 131 human cell types and tissues, and use these maps to interpret the functions of GWAS variants. Across 72 diseases and complex traits, ABC links 5,036 GWAS signals to 2,249 unique genes, including a class of 577 genes that appear to influence multiple phenotypes through variants in enhancers that act in different cell types. In inflammatory bowel disease (IBD), causal variants are enriched in predicted enhancers by more than 20-fold in particular cell types such as dendritic cells, and ABC achieves higher precision than other regulatory methods at connecting noncoding variants to target genes. These variant-to-function maps reveal an enhancer that contains an IBD risk variant and that regulates the expression of PPIF to alter the membrane potential of mitochondria in macrophages. Our study reveals principles of genome regulation, identifies genes that affect IBD and provides a resource and generalizable strategy to connect risk variants of common diseases to their molecular and cellular functions.

Intelligent electroactive material systems with self-adaptive mechanical memory and sequential logic.

By synthesizing the requisite functionalities of intelligence in an integrated material system, it may become possible to animate otherwise inanimate matter. A significant challenge in this vision is to continually sense, process, and memorize information in a decentralized way. Here, we introduce an approach that enables all such functionalities in a soft mechanical material system. By integrating nonvolatile memory with continuous processing, we develop a sequential logic-based material design framework. Soft, conductive networks interconnect with embedded electroactive actuators to enable self-adaptive behavior that facilitates autonomous toggling and counting. The design principles are scaled in processing complexity and memory capacity to develop a model 8-bit mechanical material that can solve linear algebraic equations based on analog mechanical inputs. The resulting material system operates continually to monitor the current mechanical configuration and to autonomously search for solutions within a desired error. The methods created in this work are a foundation for future synthetic general intelligence that can empower materials to autonomously react to diverse stimuli in their environment.

Selection pressures on evolution of ribonuclease H explored with rigorous free-energy-based design.

Understanding natural protein evolution and designing novel proteins are motivating interest in development of high-throughput methods to explore large sequence spaces. In this work, we demonstrate the application of multisite λ dynamics (MSλD), a rigorous free energy simulation method, and chemical denaturation experiments to quantify evolutionary selection pressure from sequence-stability relationships and to address questions of design. This study examines a mesophilic phylogenetic clade of ribonuclease H (RNase H), furthering its extensive characterization in earlier studies, focusing on E. coli RNase H (ecRNH) and a more stable consensus sequence (AncCcons) differing at 15 positions. The stabilities of 32,768 chimeras between these two sequences were computed using the MSλD framework. The most stable and least stable chimeras were predicted and tested along with several other sequences, revealing a designed chimera with approximately the same stability increase as AncCcons, but requiring only half the mutations. Comparing the computed stabilities with experiment for 12 sequences reveals a Pearson correlation of 0.86 and root mean squared error of 1.18 kcal/mol, an unprecedented level of accuracy well beyond less rigorous computational design methods. We then quantified selection pressure using a simple evolutionary model in which sequences are selected according to the Boltzmann factor of their stability. Selection temperatures from 110 to 168 K are estimated in three ways by comparing experimental and computational results to evolutionary models. These estimates indicate selection pressure is high, which has implications for evolutionary dynamics and for the accuracy required for design, and suggests accurate high-throughput computational methods like MSλD may enable more effective protein design.