Single-nuclei transcriptome analysis of the shoot apex vascular system differentiation in Populus.

Differentiation of stem cells in the plant apex gives rise to aerial tissues and organs. Presently, we lack a lineage map of the shoot apex cells in woody perennials - a crucial gap considering their role in determining primary and secondary growth. Here, we used single-nuclei RNA-sequencing to determine cell type-specific transcriptomes of the Populus vegetative shoot apex. We identified highly heterogeneous cell populations clustered into seven broad groups represented by 18 transcriptionally distinct cell clusters. Next, we established the developmental trajectories of the epidermis, leaf mesophyll and vascular tissue. Motivated by the high similarities between Populus and Arabidopsis cell population in the vegetative apex, we applied a pipeline for interspecific single-cell gene expression data integration. We contrasted the developmental trajectories of primary phloem and xylem formation in both species, establishing the first comparison of vascular development between a model annual herbaceous and a woody perennial plant species. Our results offer a valuable resource for investigating the principles underlying cell division and differentiation conserved between herbaceous and perennial species while also allowing us to examine species-specific differences at single-cell resolution.

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom.

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

Xylem systems genetics analysis reveals a key regulator of lignin biosynthesis in Populus deltoides.

Despite the growing resources and tools for high-throughput characterization and analysis of genomic information, the discovery of the genetic elements that regulate complex traits remains a challenge. Systems genetics is an emerging field that aims to understand the flow of biological information that underlies complex traits from genotype to phenotype. In this study, we used a systems genetics approach to identify and evaluate regulators of the lignin biosynthesis pathway in Populus deltoides by combining genome, transcriptome, and phenotype data from a population of 268 unrelated individuals of P. deltoides The discovery of lignin regulators began with the quantitative genetic analysis of the xylem transcriptome and resulted in the detection of 6706 and 4628 significant local- and distant-eQTL associations, respectively. Among the locally regulated genes, we identified the R2R3-MYB transcription factor MYB125 (Potri.003G114100) as a putative trans-regulator of the majority of genes in the lignin biosynthesis pathway. The expression of MYB125 in a diverse population positively correlated with lignin content. Furthermore, overexpression of MYB125 in transgenic poplar resulted in increased lignin content, as well as altered expression of genes in the lignin biosynthesis pathway. Altogether, our findings indicate that MYB125 is involved in the control of a transcriptional coexpression network of lignin biosynthesis genes during secondary cell wall formation in P. deltoides.

DNA hypomethylation of the host tree impairs interaction with mutualistic ectomycorrhizal fungus.

Ectomycorrhizas are an intrinsic component of tree nutrition and responses to environmental variations. How epigenetic mechanisms might regulate these mutualistic interactions is unknown. By manipulating the level of expression of the chromatin remodeler DECREASE IN DNA METHYLATION 1 (DDM1) and two demethylases DEMETER-LIKE (DML) in Populus tremula × Populus alba lines, we examined how host DNA methylation modulates multiple parameters of the responses to root colonization with the mutualistic fungus Laccaria bicolor. We compared the ectomycorrhizas formed between transgenic and wild-type (WT) trees and analyzed their methylomes and transcriptomes. The poplar lines displaying lower mycorrhiza formation rate corresponded to hypomethylated overexpressing DML or RNAi-ddm1 lines. We found 86 genes and 288 transposable elements (TEs) differentially methylated between WT and hypomethylated lines (common to both OX-dml and RNAi-ddm1) and 120 genes/1441 TEs in the fungal genome suggesting a host-induced remodeling of the fungal methylome. Hypomethylated poplar lines displayed 205 differentially expressed genes (cis and trans effects) in common with 17 being differentially methylated (cis). Our findings suggest a central role of host and fungal DNA methylation in the ability to form ectomycorrhizas including not only poplar genes involved in root initiation, ethylene and jasmonate-mediated pathways, and immune response but also terpenoid metabolism.

Spatiotemporal cytokinin response imaging and ISOPENTENYLTRANSFERASE 3 function in Medicago nodule development.

Most legumes can establish a symbiotic association with soil rhizobia that trigger the development of root nodules. These nodules host the rhizobia and allow them to fix nitrogen efficiently. The perception of bacterial lipo-chitooligosaccharides (LCOs) in the epidermis initiates a signaling cascade that allows rhizobial intracellular infection in the root and de-differentiation and activation of cell division that gives rise to the nodule. Thus, nodule organogenesis and rhizobial infection need to be coupled in space and time for successful nodulation. The plant hormone cytokinin (CK) contributes to the coordination of this process, acting as an essential positive regulator of nodule organogenesis. However, the temporal regulation of tissue-specific CK signaling and biosynthesis in response to LCOs or Sinorhizobium meliloti inoculation in Medicago truncatula remains poorly understood. In this study, using a fluorescence-based CK sensor (pTCSn::nls:tGFP), we performed a high-resolution tissue-specific temporal characterization of the sequential activation of CK response during root infection and nodule development in M. truncatula after inoculation with S. meliloti. Loss-of-function mutants of the CK-biosynthetic gene ISOPENTENYLTRANSFERASE 3 (IPT3) showed impairment of nodulation, suggesting that IPT3 is required for nodule development in M. truncatula. Simultaneous live imaging of pIPT3::nls:tdTOMATO and the CK sensor showed that IPT3 induction in the pericycle at the base of nodule primordium contributes to CK biosynthesis, which in turn promotes expression of positive regulators of nodule organogenesis in M. truncatula.

An LCO-responsive homolog of NODULE INCEPTION positively regulates lateral root formation in Populus sp.

The transcription factor NODULE INCEPTION (NIN) has been studied extensively for its multiple roles in root nodule symbiosis within plants of the nitrogen-fixing clade (NFC) that associate with soil bacteria, such as rhizobia and Frankia. However, NIN homologs are present in plants outside the NFC, suggesting a role in other developmental processes. Here, we show that the biofuel crop Populus sp., which is not part of the NFC, contains eight copies of NIN with diversified protein sequence and expression patterns. Lipo-chitooligosaccharides (LCOs) are produced by rhizobia and a wide range of fungi, including mycorrhizal ones, and act as symbiotic signals that promote lateral root formation. RNAseq analysis of Populus sp. treated with purified LCO showed induction of the PtNIN2 subfamily. Moreover, the expression of PtNIN2b correlated with the formation of lateral roots and was suppressed by cytokinin treatment. Constitutive expression of PtNIN2b overcame the inhibition of lateral root development by cytokinin under high nitrate conditions. Lateral root induction in response to LCOs likely represents an ancestral function of NIN retained and repurposed in nodulating plants, as we demonstrate that the role of NIN in LCO-induced root branching is conserved in both Populus sp. and legumes. We further established a visual marker of LCO perception in Populus sp. roots, the putative sulfotransferase PtSS1 that can be used to study symbiotic interactions with the bacterial and fungal symbionts of Populus sp.

The uncharacterized gene EVE contributes to vessel element dimensions in Populus.

The radiation of angiosperms led to the emergence of the vast majority of today's plant species and all our major food crops. Their extraordinary diversification occurred in conjunction with the evolution of a more efficient vascular system for the transport of water, composed of vessel elements. The physical dimensions of these water-conducting specialized cells have played a critical role in angiosperm evolution; they determine resistance to water flow, influence photosynthesis rate, and contribute to plant stature. However, the genetic factors that determine their dimensions are unclear. Here we show that a previously uncharacterized gene, ENLARGED VESSEL ELEMENT (EVE), contributes to the dimensions of vessel elements in Populus, impacting hydraulic conductivity. Our data suggest that EVE is localized in the plasma membrane and is involved in potassium uptake of differentiating xylem cells during vessel development. In plants, EVE first emerged in streptophyte algae, but expanded dramatically among vessel-containing angiosperms. The phylogeny, structure and composition of EVE indicates that it may have been involved in an ancient horizontal gene-transfer event.

Temporal change in chromatin accessibility predicts regulators of nodulation in Medicago truncatula.

BACKGROUND: Symbiotic associations between bacteria and leguminous plants lead to the formation of root nodules that fix nitrogen needed for sustainable agricultural systems. Symbiosis triggers extensive genome and transcriptome remodeling in the plant, yet an integrated understanding of the extent of chromatin changes and transcriptional networks that functionally regulate gene expression associated with symbiosis remains poorly understood. In particular, analyses of early temporal events driving this symbiosis have only captured correlative relationships between regulators and targets at mRNA level. Here, we characterize changes in transcriptome and chromatin accessibility in the model legume Medicago truncatula, in response to rhizobial signals that trigger the formation of root nodules. RESULTS: We profiled the temporal chromatin accessibility (ATAC-seq) and transcriptome (RNA-seq) dynamics of M. truncatula roots treated with bacterial small molecules called lipo-chitooligosaccharides that trigger host symbiotic pathways of nodule development. Using a novel approach, dynamic regulatory module networks, we integrated ATAC-seq and RNA-seq time courses to predict cis-regulatory elements and transcription factors that most significantly contribute to transcriptomic changes associated with symbiosis. Regulators involved in auxin (IAA4-5, SHY2), ethylene (EIN3, ERF1), and abscisic acid (ABI5) hormone response, as well as histone and DNA methylation (IBM1), emerged among those most predictive of transcriptome dynamics. RNAi-based knockdown of EIN3 and ERF1 reduced nodule number in M. truncatula validating the role of these predicted regulators in symbiosis between legumes and rhizobia. CONCLUSIONS: Our transcriptomic and chromatin accessibility datasets provide a valuable resource to understand the gene regulatory programs controlling the early stages of the dynamic process of symbiosis. The regulators identified provide potential targets for future experimental validation, and the engineering of nodulation in species is unable to establish that symbiosis naturally.

Consumption of dietary fat causes loss of olfactory sensory neurons and associated circuitry that is not mitigated by voluntary exercise in mice.

Excess nutrition causes loss of olfactory sensory neurons (OSNs) and reduces odour discrimination and odour perception in mice. To separate diet-induced obesity from the consumption of dietary fat, we designed pair-feeding experiments whereby mice were maintained on isocaloric diets for 5 months, which prevented increased fat storage. To test our hypothesis that adiposity was not a prerequisite for loss of OSNs and bulbar projections, we used male and female mice with an odorant receptor-linked genetic reporter (M72tauLacZ; Olfr160) to visualize neural circuitry changes resulting from elevated fat in the diet. Simultaneously we monitored glucose clearance (diagnostic for prediabetes), body fat deposition, ingestive behaviours, select inflammatory markers and energy metabolism. Axonal projections to defined olfactory glomeruli were visualized in whole-mount brains, and the number of OSNs was manually counted across whole olfactory epithelia. After being pair fed a moderately high-fat (MHF) diet, mice of both sexes had body weight, adipose deposits, energy expenditure, respiratory exchange ratios and locomotor activity that were unchanged from control-fed mice. Despite this, they were still found to lose OSNs and associated bulbar projections. Even with unchanged adipocyte storage, pair-fed animals had an elevation in TNF cytokines and an intermediate ability for glucose clearance. Albeit improving health metrics, access to voluntary running while consuming an ad libitum fatty diet still precipitated a loss of OSNs and associated axonal projections for male mice. Our results support that long-term macronutrient imbalance can drive anatomical loss in the olfactory system regardless of total energy expenditure. KEY POINTS: Obesity can disrupt the structure and function of organ systems, including the olfactory system that is important for food selection and satiety. We designed dietary treatments in mice such that mice received fat, but the total calories provided were the same as in control diets so that they would not gain weight or increase adipose tissue. Mice that were not obese but consumed isocaloric fatty diets still lost olfactory neuronal circuits, had fewer numbers of olfactory neurons, had an elevation in inflammatory signals and had an intermediate ability to clear glucose (prediabetes). Mice were allowed access to running wheels while consuming fatty diets, yet still lost olfactory structures. We conclude that a long-term imbalance in nutrition that favours fat in the diet disrupts the olfactory system of mice in the absence of obesity.

Functional and comparative genomics reveals conserved noncoding sequences in the nitrogen-fixing clade.

Nitrogen is one of the most inaccessible plant nutrients, but certain species have overcome this limitation by establishing symbiotic interactions with nitrogen-fixing bacteria in the root nodule. This root-nodule symbiosis (RNS) is restricted to species within a single clade of angiosperms, suggesting a critical, but undetermined, evolutionary event at the base of this clade. To identify putative regulatory sequences implicated in the evolution of RNS, we evaluated the genomes of 25 species capable of nodulation and identified 3091 conserved noncoding sequences (CNS) in the nitrogen-fixing clade (NFC). We show that the chromatin accessibility of 452 CNS correlates significantly with the regulation of genes responding to lipochitooligosaccharides in Medicago truncatula. These included 38 CNS in proximity to 19 known genes involved in RNS. Five such regions are upstream of MtCRE1, Cytokinin Response Element 1, required to activate a suite of downstream transcription factors necessary for nodulation in M. truncatula. Genetic complementation of an Mtcre1 mutant showed a significant decrease of nodulation in the absence of the five CNS, when they are driving the expression of a functional copy of MtCRE1. CNS identified in the NFC may harbor elements required for the regulation of genes controlling RNS in M. truncatula.

Molecular Dynamics Simulations of Transmembrane Cyclic Peptide Nanotubes Using Classical Force Fields, Hydrogen Mass Repartitioning, and Hydrogen Isotope Exchange Methods: A Critical Comparison.

Self-assembled cyclic peptide nanotubes with alternating D- and L-amino acid residues in the sequence of each subunit have attracted a great deal of attention due to their potential for new nanotechnology and biomedical applications, mainly in the field of antimicrobial peptides. Molecular dynamics simulations can be used to characterize these systems with atomic resolution at different time scales, providing information that is difficult to obtain via wet lab experiments. However, the performance of classical force fields typically employed in the simulation of biomolecules has not yet been extensively tested with this kind of highly constrained peptide. Four different classical force fields (AMBER, CHARMM, OPLS, and GROMOS), using a nanotube formed by eight D,L-α-cyclic peptides inserted into a lipid bilayer as a model system, were employed here to fill this gap. Significant differences in the pseudo-cylindrical cavities formed by the nanotubes were observed, the most important being the diameter of the nanopores, the number and location of confined water molecules, and the density distribution of the solvent molecules. Furthermore, several modifications were performed on GROMOS54a7, aiming to explore acceleration strategies of the MD simulations. The hydrogen mass repartitioning (HMR) and hydrogen isotope exchange (HIE) methods were tested to slow down the fastest degrees of freedom. These approaches allowed a significant increase in the time step employed in the equation of the motion integration algorithm, from 2 fs up to 5-7 fs, with no serious changes in the structural and dynamical properties of the nanopores. Subtle differences with respect to the simulations with the unmodified force fields were observed in the concerted movements of the cyclic peptides, as well as in the lifetime of several H-bonds. All together, these results are expected to contribute to better understanding of the behavior of self-assembled cyclic peptide nanotubes, as well as to support the methods tested to speed up general MD simulations; additionally, they do provide a number of quantitative descriptors that are expected to be used as a reference to design new experiments intended to validate and complement computational studies of antimicrobial cyclic peptides.

Chilling-responsive DEMETER-LIKE DNA demethylase mediates in poplar bud break.

Annual dormancy-growth cycle is a developmental and physiological process essential for the survival of deciduous trees in temperate and boreal forests. Seasonal control of shoot growth in woody perennials requires specific genetic programmes responding to environmental signals. The environmental-controlled mechanisms that regulate the shift between winter dormancy and the growth-promoting genetic programmes are still unknown. Here, we show that dynamics in genomic DNA methylation levels are involved in the regulation of dormancy-growth cycle in poplar. The reactivation of growth in the apical shoot during bud break process in spring is preceded by a progressive reduction of genomic DNA methylation in apex tissue. The induction in apex tissue of a chilling-dependent poplar DEMETER-LIKE 10 (PtaDML10) DNA demethylase precedes shoot growth reactivation. Transgenic poplars showing downregulation of PtaDML8/10 caused delayed bud break. Genome-wide transcriptome and methylome analysis and data mining revealed that the gene targets of DEMETER-LIKE-dependent DNA demethylation are genetically associated with bud break. These data point to a chilling-dependent DEMETER-like DNA demethylase mechanisms being involved in the shift from winter dormancy to a condition that precedes shoot apical vegetative growth in poplar.

Overexpression of DEMETER, a DNA demethylase, promotes early apical bud maturation in poplar.

The transition from active growth to dormancy is critical for the survival of perennial plants. We identified a DEMETER-like (CsDML) cDNA from a winter-enriched cDNA subtractive library in chestnut (Castanea sativa Mill.), an economically and ecologically important species. Next, we characterized this DNA demethylase and its putative ortholog in the more experimentally tractable hybrid poplar (Populus tremula × alba), under the signals that trigger bud dormancy in trees. We performed phylogenetic and protein sequence analysis, gene expression profiling, and 5-methyl-cytosine methylation immunodetection studies to evaluate the role of CsDML and its homolog in poplar, PtaDML6. Transgenic hybrid poplars overexpressing CsDML were produced and analysed. Short days and cold temperatures induced CsDML and PtaDML6. Overexpression of CsDML accelerated short-day-induced bud formation, specifically from Stages 1 to 0. Buds acquired a red-brown coloration earlier than wild-type plants, alongside with the up-regulation of flavonoid biosynthesis enzymes and accumulation of flavonoids in the shoot apical meristem and bud scales. Our data show that the CsDML gene induces bud formation needed for the survival of the apical meristem under the harsh conditions of winter.

Shifts in evolutionary lability underlie independent gains and losses of root-nodule symbiosis in a single clade of plants.

Root nodule symbiosis (RNS) is a complex trait that enables plants to access atmospheric nitrogen converted into usable forms through a mutualistic relationship with soil bacteria. Pinpointing the evolutionary origins of RNS is critical for understanding its genetic basis, but building this evolutionary context is complicated by data limitations and the intermittent presence of RNS in a single clade of ca. 30,000 species of flowering plants, i.e., the nitrogen-fixing clade (NFC). We developed the most extensive de novo phylogeny for the NFC and an RNS trait database to reconstruct the evolution of RNS. Our analysis identifies evolutionary rate heterogeneity associated with a two-step process: An ancestral precursor state transitioned to a more labile state from which RNS was rapidly gained at multiple points in the NFC. We illustrate how a two-step process could explain multiple independent gains and losses of RNS, contrary to recent hypotheses suggesting one gain and numerous losses, and suggest a broader phylogenetic and genetic scope may be required for genome-phenome mapping.

The single-cell transcriptome program of nodule development cellular lineages in Medicago truncatula.

Legumes establish a symbiotic relationship with nitrogen-fixing rhizobia by developing nodules. Nodules are modified lateral roots that undergo changes in their cellular development in response to bacteria, but the transcriptional reprogramming that occurs in these root cells remains largely uncharacterized. Here, we describe the cell-type-specific transcriptome response of Medicago truncatula roots to rhizobia during early nodule development in the wild-type genotype Jemalong A17, complemented with a hypernodulating mutant (sunn-4) to expand the cell population responding to infection and subsequent biological inferences. The analysis identifies epidermal root hair and stele sub-cell types associated with a symbiotic response to infection and regulation of nodule proliferation. Trajectory inference shows cortex-derived cell lineages differentiating to form the nodule primordia and, posteriorly, its meristem, while modulating the regulation of phytohormone-related genes. Gene regulatory analysis of the cell transcriptomes identifies new regulators of nodulation, including STYLISH 4, for which the function is validated.

Digital entrepreneurship from cellular data: How omics afford the emergence of a new wave of digital ventures in health.

Data has become an indispensable input, throughput, and output for the healthcare industry. In recent years, omics technologies such as genomics and proteomics have generated vast amounts of new data at the cellular level including molecular, structural, and functional levels. Cellular data holds the potential to innovate therapeutics, vaccines, diagnostics, consumer products, or even ancestry services. However, data at the cellular level is generated with rapidly evolving omics technologies. These technologies use scientific knowledge from resource-rich environments. This raises the question of how new ventures can use cellular-level data from omics technologies to create new products and scale their business. We report on a series of interviews and a focus group discussion with entrepreneurs, investors, and data providers. By conceptualizing omics technologies as external enablers, we show how characteristics of cellular-level data negatively affect the combination mechanisms that drive venture creation and growth. We illustrate how data characteristics set boundary conditions for innovation and entrepreneurship and highlight how ventures seek to mitigate their impact. Supplementary Information: The online version contains supplementary material available at 10.1007/s12525-023-00669-w.

Land suitability assessment in the catchment area of four Southwestern Atlantic coastal lagoons: multicriteria and optimization modeling.

In the present study, a land suitability assessment was conducted in the basin of four Uruguayan coastal lagoons (Southwestern Atlantic) to analyze the productive development while minimizing eutrophication, biodiversity loss and conflicts among different land uses. Suitable land for agriculture, forest, livestock ranching, tourism and conservation sectors were initially established based on a multi-attribute model developed using a geographic information system. Experts were consulted to determine the requirements for each land use sector and the incompatibilities among land use types. The current and potential conflicts among incompatible land use sectors were analyzed by overlapping land suitability maps. We subsequently applied a multi-objective model where land (pixels) with similar suitability was clustered into "land suitability groups", using a two-phase cluster analysis and the Akaike Information Criterion. Finally, a linear programming optimization procedure was applied to allocate land use sectors into land suitable groups, maximizing total suitability and minimizing interference among sectors. Results indicated that current land use overlapped by 4.7 % with suitable land of other incompatible sectors. However, the suitable land of incompatible sectors overlapped in 20.3 % of the study area, indicating a high potential for the occurrence of future conflict. The highest competition was between agriculture and conservation, followed by forest and agriculture. We explored scenarios where livestock ranching and tourism intensified, and found that interference with conservation and agriculture notably increased. This methodology allowed us to analyze current and potential land use conflicts and to contribute to the strategic planning of the study area.

Metastatic extra-abdominal presentation of gastrointestinal stromal tumors in a young Hispanic male: a case presentation and literature review.

This article presents the medical history and management of a 44-year-old Hispanic male with metastatic extra-abdominal gastrointestinal stromal tumor including a literature review on this rare clinical presentation.

Decoding exceptional plant traits by comparative single-cell genomics.

Some plants have acquired traits of remarkable adaptive value to thrive under stress. Transferring these traits to crops could improve agriculture, but uncovering the toolkit required has remained largely elusive. We propose that single-cell genomics offers a framework to compare species with contrasting developmental traits and to identify the regulators of evolutionary innovations.

Mobility in pandemic times: Exploring changes and long-term effects of COVID-19 on urban mobility behavior.

The COVID-19 pandemic marked a global disruption of unprecedented scale which was closely associated with human mobility. Since mobility acts as a facilitator for spreading the virus, individuals were forced to reconsider their respective behaviors. Despite numerous studies having detected behavioral changes during the first lockdown period (spring 2020), there is a lack of longitudinal perspectives that can provide insights into the intra-pandemic dynamics and potential long-term effects. This article investigates COVID-19-induced mobility-behavioral transformations by analyzing travel patterns of Berlin residents during a 20-month pandemic period and comparing them to the pre-pandemic situation. Based on quantitative analysis of almost 800,000 recorded trips, our longitudinal examination revealed individuals having reduced average monthly travel distances by ∼20%, trip frequencies by ∼11%, and having switched to individual modes. Public transportation has suffered a continual regression, with trip frequencies experiencing a relative long-term reduction of ∼50%, and a respective decrease of traveled distances by ∼43%. In contrast, the bicycle (rather than the car) was the central beneficiary, indicated by bicycle-related trip frequencies experiencing a relative long-term increase of ∼53%, and travel distances increasing by ∼117%. Comparing behavioral responses to three pandemic waves, our analysis revealed each wave to have created unique response patterns, which show a gradual softening of individuals' mobility related self-restrictions. Our findings contribute to retracing and quantifying individuals' changing mobility behaviors induced by the pandemic, and to detecting possible long-term effects that may constitute a "new normal" of an entirely altered urban mobility landscape.