Time resolution in cryo-EM using a PDMS-based microfluidic chip assembly and its application to the study of HflX-mediated ribosome recycling.

The rapid kinetics of biological processes and associated short-lived conformational changes pose a significant challenge in attempts to structurally visualize biomolecules during a reaction in real time. Conventionally, on-pathway intermediates have been trapped using chemical modifications or reduced temperature, giving limited insights. Here, we introduce a time-resolved cryo-EM method using a reusable PDMS-based microfluidic chip assembly with high reactant mixing efficiency. Coating of PDMS walls with SiO2 virtually eliminates non-specific sample adsorption and ensures maintenance of the stoichiometry of the reaction, rendering it highly reproducible. In an operating range from 10 to 1,000 ms, the device allows us to follow in vitro reactions of biological molecules at resolution levels in the range of 3 Å. By employing this method, we show the mechanism of progressive HflX-mediated splitting of the 70S E. coli ribosome in the presence of the GTP via capture of three high-resolution reaction intermediates within 140 ms.

Time-series reconstruction of the molecular architecture of human centriole assembly.

Centriole biogenesis, as in most organelle assemblies, involves the sequential recruitment of sub-structural elements that will support its function. To uncover this process, we correlated the spatial location of 24 centriolar proteins with structural features using expansion microscopy. A time-series reconstruction of protein distributions throughout human procentriole assembly unveiled the molecular architecture of the centriole biogenesis steps. We found that the process initiates with the formation of a naked cartwheel devoid of microtubules. Next, the bloom phase progresses with microtubule blade assembly, concomitantly with radial separation and rapid cartwheel growth. In the subsequent elongation phase, the tubulin backbone grows linearly with the recruitment of the A-C linker, followed by proteins of the inner scaffold (IS). By following six structural modules, we modeled 4D assembly of the human centriole. Collectively, this work provides a framework to investigate the spatial and temporal assembly of large macromolecules.

Time-resolved fate mapping identifies the intestinal upper crypt zone as an origin of Lgr5+ crypt base columnar cells.

In the prevailing model, Lgr5+ cells are the only intestinal stem cells (ISCs) that sustain homeostatic epithelial regeneration by upward migration of progeny through elusive upper crypt transit-amplifying (TA) intermediates. Here, we identify a proliferative upper crypt population marked by Fgfbp1, in the location of putative TA cells, that is transcriptionally distinct from Lgr5+ cells. Using a kinetic reporter for time-resolved fate mapping and Fgfbp1-CreERT2 lineage tracing, we establish that Fgfbp1+ cells are multi-potent and give rise to Lgr5+ cells, consistent with their ISC function. Fgfbp1+ cells also sustain epithelial regeneration following Lgr5+ cell depletion. We demonstrate that FGFBP1, produced by the upper crypt cells, is an essential factor for crypt proliferation and epithelial homeostasis. Our findings support a model in which tissue regeneration originates from upper crypt Fgfbp1+ cells that generate progeny propagating bi-directionally along the crypt-villus axis and serve as a source of Lgr5+ cells in the crypt base.

Managing the Steady State Chromatin Landscape by Nucleosome Dynamics.

Gene regulation arises out of dynamic competition between nucleosomes, transcription factors, and other chromatin proteins for the opportunity to bind genomic DNA. The timescales of nucleosome assembly and binding of factors to DNA determine the outcomes of this competition at any given locus. Here, we review how these properties of chromatin proteins and the interplay between the dynamics of different factors are critical for gene regulation. We discuss how molecular structures of large chromatin-associated complexes, kinetic measurements, and high resolution mapping of protein-DNA complexes in vivo set the boundary conditions for chromatin dynamics, leading to models of how the steady state behaviors of regulatory elements arise.

Fundamental behaviors emerge from simulations of a living minimal cell.

We present a whole-cell fully dynamical kinetic model (WCM) of JCVI-syn3A, a minimal cell with a reduced genome of 493 genes that has retained few regulatory proteins or small RNAs. Cryo-electron tomograms provide the cell geometry and ribosome distributions. Time-dependent behaviors of concentrations and reaction fluxes from stochastic-deterministic simulations over a cell cycle reveal how the cell balances demands of its metabolism, genetic information processes, and growth, and offer insight into the principles of life for this minimal cell. The energy economy of each process including active transport of amino acids, nucleosides, and ions is analyzed. WCM reveals how emergent imbalances lead to slowdowns in the rates of transcription and translation. Integration of experimental data is critical in building a kinetic model from which emerges a genome-wide distribution of mRNA half-lives, multiple DNA replication events that can be compared to qPCR results, and the experimentally observed doubling behavior.

Quantitative fate mapping: A general framework for analyzing progenitor state dynamics via retrospective lineage barcoding.

Natural and induced somatic mutations that accumulate in the genome during development record the phylogenetic relationships of cells; whether these lineage barcodes capture the complex dynamics of progenitor states remains unclear. We introduce quantitative fate mapping, an approach to reconstruct the hierarchy, commitment times, population sizes, and commitment biases of intermediate progenitor states during development based on a time-scaled phylogeny of their descendants. To reconstruct time-scaled phylogenies from lineage barcodes, we introduce Phylotime, a scalable maximum likelihood clustering approach based on a general barcoding mutagenesis model. We validate these approaches using realistic in silico and in vitro barcoding experiments. We further establish criteria for the number of cells that must be analyzed for robust quantitative fate mapping and a progenitor state coverage statistic to assess the robustness. This work demonstrates how lineage barcodes, natural or synthetic, enable analyzing progenitor fate and dynamics long after embryonic development in any organism.

Recording of elapsed time and temporal information about biological events using Cas9.

DNA has not been utilized to record temporal information, although DNA has been used to record biological information and to compute mathematical problems. Here, we found that indel generation by Cas9 and guide RNA can occur at steady rates, in contrast to typical dynamic biological reactions, and the accumulated indel frequency can be a function of time. By measuring indel frequencies, we developed a method for recording and measuring absolute time periods over hours to weeks in mammalian cells. These time-recordings were conducted in several cell types, with different promoters and delivery vectors for Cas9, and in both cultured cells and cells of living mice. As applications, we recorded the duration of chemical exposure and the lengths of elapsed time since the onset of biological events (e.g., heat exposure and inflammation). We propose that our systems could serve as synthetic "DNA clocks."

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

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

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.

Genomic Epidemiology of SARS-CoV-2 in Guangdong Province, China.

Coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2 infection and was first reported in central China in December 2019. Extensive molecular surveillance in Guangdong, China's most populous province, during early 2020 resulted in 1,388 reported RNA-positive cases from 1.6 million tests. In order to understand the molecular epidemiology and genetic diversity of SARS-CoV-2 in China, we generated 53 genomes from infected individuals in Guangdong using a combination of metagenomic sequencing and tiling amplicon approaches. Combined epidemiological and phylogenetic analyses indicate multiple independent introductions to Guangdong, although phylogenetic clustering is uncertain because of low virus genetic variation early in the pandemic. Our results illustrate how the timing, size, and duration of putative local transmission chains were constrained by national travel restrictions and by the province's large-scale intensive surveillance and intervention measures. Despite these successes, COVID-19 surveillance in Guangdong is still required, because the number of cases imported from other countries has increased.

Molecular Choreography of Acute Exercise.

Acute physical activity leads to several changes in metabolic, cardiovascular, and immune pathways. Although studies have examined selected changes in these pathways, the system-wide molecular response to an acute bout of exercise has not been fully characterized. We performed longitudinal multi-omic profiling of plasma and peripheral blood mononuclear cells including metabolome, lipidome, immunome, proteome, and transcriptome from 36 well-characterized volunteers, before and after a controlled bout of symptom-limited exercise. Time-series analysis revealed thousands of molecular changes and an orchestrated choreography of biological processes involving energy metabolism, oxidative stress, inflammation, tissue repair, and growth factor response, as well as regulatory pathways. Most of these processes were dampened and some were reversed in insulin-resistant participants. Finally, we discovered biological pathways involved in cardiopulmonary exercise response and developed prediction models revealing potential resting blood-based biomarkers of peak oxygen consumption.

Bacteriorhodopsin: Structural Insights Revealed Using X-Ray Lasers and Synchrotron Radiation.

Directional transport of protons across an energy transducing membrane-proton pumping-is ubiquitous in biology. Bacteriorhodopsin (bR) is a light-driven proton pump that is activated by a buried all-trans retinal chromophore being photoisomerized to a 13-cis conformation. The mechanism by which photoisomerization initiates directional proton transport against a proton concentration gradient has been studied by a myriad of biochemical, biophysical, and structural techniques. X-ray free electron lasers (XFELs) have created new opportunities to probe the structural dynamics of bR at room temperature on timescales from femtoseconds to milliseconds using time-resolved serial femtosecond crystallography (TR-SFX). Wereview these recent developments and highlight where XFEL studies reveal new details concerning the structural mechanism of retinal photoisomerization and proton pumping. We also discuss the extent to which these insights were anticipated by earlier intermediate trapping studies using synchrotron radiation. TR-SFX will open up the field for dynamical studies of other proteins that are not naturally light-sensitive.

X-Ray Free-Electron Lasers for the Structure and Dynamics of Macromolecules.

X-ray free-electron lasers provide femtosecond-duration pulses of hard X-rays with a peak brightness approximately one billion times greater than is available at synchrotron radiation facilities. One motivation for the development of such X-ray sources was the proposal to obtain structures of macromolecules, macromolecular complexes, and virus particles, without the need for crystallization, through diffraction measurements of single noncrystalline objects. Initial explorations of this idea and of outrunning radiation damage with femtosecond pulses led to the development of serial crystallography and the ability to obtain high-resolution structures of small crystals without the need for cryogenic cooling. This technique allows the understanding of conformational dynamics and enzymatics and the resolution of intermediate states in reactions over timescales of 100 fs to minutes. The promise of more photons per atom recorded in a diffraction pattern than electrons per atom contributing to an electron micrograph may enable diffraction measurements of single molecules, although challenges remain.

Biophysical Techniques in Structural Biology.

Over the past six decades, steadily increasing progress in the application of the principles and techniques of the physical sciences to the study of biological systems has led to remarkable insights into the molecular basis of life. Of particular significance has been the way in which the determination of the structures and dynamical properties of proteins and nucleic acids has so often led directly to a profound understanding of the nature and mechanism of their functional roles. The increasing number and power of experimental and theoretical techniques that can be applied successfully to living systems is now ushering in a new era of structural biology that is leading to fundamentally new information about the maintenance of health, the origins of disease, and the development of effective strategies for therapeutic intervention. This article provides a brief overview of some of the most powerful biophysical methods in use today, along with references that provide more detailed information about recent applications of each of them. In addition, this article acts as an introduction to four authoritative reviews in this volume. The first shows the ways that a multiplicity of biophysical methods can be combined with computational techniques to define the architectures of complex biological systems, such as those involving weak interactions within ensembles of molecular components. The second illustrates one aspect of this general approach by describing how recent advances in mass spectrometry, particularly in combination with other techniques, can generate fundamentally new insights into the properties of membrane proteins and their functional interactions with lipid molecules. The third reviewdemonstrates the increasing power of rapidly evolving diffraction techniques, employing the very short bursts of X-rays of extremely high intensity that are now accessible as a result of the construction of free-electron lasers, in particular to carry out time-resolved studies of biochemical reactions. The fourth describes in detail the application of such approaches to probe the mechanism of the light-induced changes associated with bacteriorhodopsin's ability to convert light energy into chemical energy.

Characterizing the Major Structural Variant Alleles of the Human Genome.

In order to provide a comprehensive resource for human structural variants (SVs), we generated long-read sequence data and analyzed SVs for fifteen human genomes. We sequence resolved 99,604 insertions, deletions, and inversions including 2,238 (1.6 Mbp) that are shared among all discovery genomes with an additional 13,053 (6.9 Mbp) present in the majority, indicating minor alleles or errors in the reference. Genotyping in 440 additional genomes confirms the most common SVs in unique euchromatin are now sequence resolved. We report a ninefold SV bias toward the last 5 Mbp of human chromosomes with nearly 55% of all VNTRs (variable number of tandem repeats) mapping to this portion of the genome. We identify SVs affecting coding and noncoding regulatory loci improving annotation and interpretation of functional variation. These data provide the framework to construct a canonical human reference and a resource for developing advanced representations capable of capturing allelic diversity.

Processing Temporal Growth Factor Patterns by an Epidermal Growth Factor Receptor Network Dynamically Established in Space.

The proto-oncogenic epidermal growth factor (EGF) receptor (EGFR) is a tyrosine kinase whose sensitivity and response to growth factor signals that vary over time and space determine cellular behavior within a developing tissue. The molecular reorganization of the receptors on the plasma membrane and the enzyme-kinetic mechanisms of phosphorylation are key determinants that couple growth factor binding to EGFR signaling. To enable signal initiation and termination while simultaneously accounting for suppression of aberrant signaling, a coordinated coupling of EGFR kinase and protein tyrosine phosphatase activity is established through space by vesicular dynamics. The dynamical operation mode of this network enables not only time-varying growth factor sensing but also adaptation of the response depending on cellular context. By connecting spatially coupled enzymatic kinase/phosphatase processes and the corresponding dynamical systems description of the EGFR network, we elaborate on the general principles necessary for processing complex growth factor signals.

Reversible, tunable epigenetic silencing of TCF1 generates flexibility in the T cell memory decision.

The immune system encodes information about the severity of a pathogenic threat in the quantity and type of memory cells it forms. This encoding emerges from lymphocyte decisions to maintain or lose self-renewal and memory potential during a challenge. By tracking CD8+ T cells at the single-cell and clonal lineage level using time-resolved transcriptomics, quantitative live imaging, and an acute infection model, we find that T cells will maintain or lose memory potential early after antigen recognition. However, following pathogen clearance, T cells may regain memory potential if initially lost. Mechanistically, this flexibility is implemented by a stochastic cis-epigenetic switch that tunably and reversibly silences the memory regulator, TCF1, in response to stimulation. Mathematical modeling shows how this flexibility allows memory T cell numbers to scale robustly with pathogen virulence and immune response magnitudes. We propose that flexibility and stochasticity in cellular decisions ensure optimal immune responses against diverse threats.

An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics.

For a decade, The Cancer Genome Atlas (TCGA) program collected clinicopathologic annotation data along with multi-platform molecular profiles of more than 11,000 human tumors across 33 different cancer types. TCGA clinical data contain key features representing the democratized nature of the data collection process. To ensure proper use of this large clinical dataset associated with genomic features, we developed a standardized dataset named the TCGA Pan-Cancer Clinical Data Resource (TCGA-CDR), which includes four major clinical outcome endpoints. In addition to detailing major challenges and statistical limitations encountered during the effort of integrating the acquired clinical data, we present a summary that includes endpoint usage recommendations for each cancer type. These TCGA-CDR findings appear to be consistent with cancer genomics studies independent of the TCGA effort and provide opportunities for investigating cancer biology using clinical correlates at an unprecedented scale.

The Clockwork Embryo: Mechanisms Regulating Developmental Rate.

Organismal development requires the reproducible unfolding of an ordered sequence of discrete steps (cell fate determination, migration, tissue folding, etc.) in both time and space. Here, we review the mechanisms that grant temporal specificity to developmental steps, including molecular clocks and timers. Individual timing mechanisms must be coordinated with each other to maintain the overall developmental sequence. However, phenotypic novelties can also arise through the modification of temporal patterns over the course of evolution. Two main types of variation in temporal patterning characterize interspecies differences in developmental time: allochrony, where the overall developmental sequence is either accelerated or slowed down while maintaining the relative duration of individual steps, and heterochrony, where the duration of specific developmental steps is altered relative to the rest. New advances in in vitro modeling of mammalian development using stem cells have recently enabled the revival of mechanistic studies of allochrony and heterochrony. In both cases, differences in the rate of basic cellular functions such as splicing, translation, protein degradation, and metabolism seem to underlie differences in developmental time. In the coming years, these studies should identify the genetic differences that drive divergence in developmental time between species.

A self-sustaining layer of early-life-origin B cells drives steady-state IgA responses in the adult gut.

The adult immune system consists of cells that emerged at various times during ontogeny. We aimed to define the relationship between developmental origin and composition of the adult B cell pool during unperturbed hematopoiesis. Lineage tracing stratified murine adult B cells based on the timing of output, revealing that a substantial portion originated within a restricted neonatal window. In addition to B-1a cells, early-life time-stamped B cells included clonally interrelated IgA plasma cells in the gut and bone marrow. These were actively maintained by B cell memory within gut chronic germinal centers and contained commensal microbiota reactivity. Neonatal rotavirus infection recruited recurrent IgA clones that were distinct from those arising by infection with the same antigen in adults. Finally, gut IgA plasma cells arose from the same hematopoietic progenitors as B-1a cells during ontogeny. Thus, a complex layer of neonatally imprinted B cells confer unique antibody responses later in life.