Author Correction: Tissue patrol by resident memory CD8+ T cells in human skin.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Tissue patrol by resident memory CD8+ T cells in human skin.

Emerging data show that tissue-resident memory T (TRM) cells play an important protective role at murine and human barrier sites. TRM cells in the epidermis of mouse skin patrol their surroundings and rapidly respond when antigens are encountered. However, whether a similar migratory behavior is performed by human TRM cells is unclear, as technology to longitudinally follow them in situ has been lacking. To address this issue, we developed an ex vivo culture system to label and track T cells in fresh skin samples. We validated this system by comparing in vivo and ex vivo properties of murine TRM cells. Using nanobody labeling, we subsequently demonstrated in human ex vivo skin that CD8+ TRM cells migrated through the papillary dermis and the epidermis, below sessile Langerhans cells. Collectively, this work allows the dynamic study of resident immune cells in human skin and provides evidence of tissue patrol by human CD8+ TRM cells.

Quorum Regulation via Nested Antagonistic Feedback Circuits Mediated by the Receptors CD28 and CTLA-4 Confers Robustness to T Cell Population Dynamics.

T cell responses upon infection display a remarkably reproducible pattern of expansion, contraction, and memory formation. If the robustness of this pattern builds entirely on signals derived from other cell types or if activated T cells themselves contribute to the orchestration of these population dynamics-akin to bacterial quorum regulation-is unclear. Here, we examined this question using time-lapse microscopy, genetic perturbation, bioinformatic predictions, and mathematical modeling. We found that ICAM-1-mediated cell clustering enabled CD8+ T cells to collectively regulate the balance between proliferation and apoptosis. Mechanistically, T cell expressed CD80 and CD86 interacted with the receptors CD28 and CTLA-4 on neighboring T cells; these interactions fed two nested antagonistic feedback circuits that regulated interleukin 2 production in a manner dependent on T cell density as confirmed by in vivo modulation of this network. Thus, CD8+ T cell-population-intrinsic mechanisms regulate cellular behavior, thereby promoting robustness of population dynamics.

Model-based translation of DNA damage signaling dynamics across cell types.

Interindividual variability in DNA damage response (DDR) dynamics may evoke differences in susceptibility to cancer. However, pathway dynamics are often studied in cell lines as alternative to primary cells, disregarding variability. To compare DDR dynamics in the cell line HepG2 with primary human hepatocytes (PHHs), we developed a HepG2-based computational model that describes the dynamics of DDR regulator p53 and targets MDM2, p21 and BTG2. We used this model to generate simulations of virtual PHHs and compared the results to those for PHH donor samples. Correlations between baseline p53 and p21 or BTG2 mRNA expression in the absence and presence of DNA damage for HepG2-derived virtual samples matched the moderately positive correlations observed for 50 PHH donor samples, but not the negative correlations between p53 and its inhibitor MDM2. Model parameter manipulation that affected p53 or MDM2 dynamics was not sufficient to accurately explain the negative correlation between these genes. Thus, extrapolation from HepG2 to PHH can be done for some DDR elements, yet our analysis also reveals a knowledge gap within p53 pathway regulation, which makes such extrapolation inaccurate for the regulator MDM2. This illustrates the relevance of studying pathway dynamics in addition to gene expression comparisons to allow reliable translation of cellular responses from cell lines to primary cells. Overall, with our approach we show that dynamical modeling can be used to improve our understanding of the sources of interindividual variability of pathway dynamics.

High-content high-throughput imaging reveals distinct connections between mitochondrial morphology and functionality for OXPHOS complex I, III, and V inhibitors.

Cells can adjust their mitochondrial morphology by altering the balance between mitochondrial fission and fusion to adapt to stressful conditions. The connection between a chemical perturbation, changes in mitochondrial function, and altered mitochondrial morphology is not well understood. Here, we made use of high-throughput high-content confocal microscopy to assess the effects of distinct classes of oxidative phosphorylation (OXPHOS) complex inhibitors on mitochondrial parameters in a concentration and time resolved manner. Mitochondrial morphology phenotypes were clustered based on machine learning algorithms and mitochondrial integrity patterns were mapped. In parallel, changes in mitochondrial membrane potential (MMP), mitochondrial and cellular ATP levels, and viability were microscopically assessed. We found that inhibition of MMP, mitochondrial ATP production, and oxygen consumption rate (OCR) using sublethal concentrations of complex I and III inhibitors did not trigger mitochondrial fragmentation. Instead, complex V inhibitors that suppressed ATP and OCR but increased MMP provoked a more fragmented mitochondrial morphology. In agreement, complex V but not complex I or III inhibitors triggered proteolytic cleavage of the mitochondrial fusion protein, OPA1. The relation between increased MMP and fragmentation did not extend beyond OXPHOS complex inhibitors: increasing MMP by blocking the mPTP pore did not lead to OPA1 cleavage or mitochondrial fragmentation and the OXPHOS uncoupler FCCP was associated with OPA1 cleavage and MMP reduction. Altogether, our findings connect vital mitochondrial functions and phenotypes in a high-throughput high-content confocal microscopy approach that help understanding of chemical-induced toxicity caused by OXPHOS complex perturbing chemicals.

Application of AOPs to assist regulatory assessment of chemical risks - Case studies, needs and recommendations.

While human regulatory risk assessment (RA) still largely relies on animal studies, new approach methodologies (NAMs) based on in vitro, in silico or non-mammalian alternative models are increasingly used to evaluate chemical hazards. Moreover, human epidemiological studies with biomarkers of effect (BoE) also play an invaluable role in identifying health effects associated with chemical exposures. To move towards the next generation risk assessment (NGRA), it is therefore crucial to establish bridges between NAMs and standard approaches, and to establish processes for increasing mechanistically-based biological plausibility in human studies. The Adverse Outcome Pathway (AOP) framework constitutes an important tool to address these needs but, despite a significant increase in knowledge and awareness, the use of AOPs in chemical RA remains limited. The objective of this paper is to address issues related to using AOPs in a regulatory context from various perspectives as it was discussed in a workshop organized within the European Union partnerships HBM4EU and PARC in spring 2022. The paper presents examples where the AOP framework has been proven useful for the human RA process, particularly in hazard prioritization and characterization, in integrated approaches to testing and assessment (IATA), and in the identification and validation of BoE in epidemiological studies. Nevertheless, several limitations were identified that hinder the optimal usability and acceptance of AOPs by the regulatory community including the lack of quantitative information on response-response relationships and of efficient ways to map chemical data (exposure and toxicity) onto AOPs. The paper summarizes suggestions, ongoing initiatives and third-party tools that may help to overcome these obstacles and thus assure better implementation of AOPs in the NGRA.

Stimulation of de novo glutathione synthesis by nitrofurantoin for enhanced resilience of hepatocytes.

Toxicity is not only a function of damage mechanisms, but is also determined by cellular resilience factors. Glutathione has been reported as essential element to counteract negative influences. The present work hence pursued the question how intracellular glutathione can be elevated transiently to render cells more resistant toward harmful conditions. The antibiotic nitrofurantoin (NFT) was identified to stimulate de novo synthesis of glutathione in the human hepatoma cell line, HepG2, and in primary human hepatocytes. In intact cells, activation of NFT yielded a radical anion, which subsequently initiated nuclear-factor-erythroid 2-related-factor-2 (Nrf2)-dependent induction of glutamate cysteine ligase (GCL). Application of siRNA-based intervention approaches confirmed the involvement of the Nrf2-GCL axis in the observed elevation of intracellular glutathione levels. Quantitative activation of Nrf2 by NFT, and the subsequent rise in glutathione, were similar as observed with the potent experimental Nrf2 activator diethyl maleate. The elevation of glutathione levels, observed even 48 h after withdrawal of NFT, rendered cells resistant to different stressors such as the mitochondrial inhibitor rotenone, the redox cycler paraquat, the proteasome inhibitors MG-132 or bortezomib, or high concentrations of NFT. Repurpose of the antibiotic NFT as activator of Nrf2 could thus be a promising strategy for a transient and targeted activation of the endogenous antioxidant machinery. Graphical abstract.

Mapping the cellular response to electron transport chain inhibitors reveals selective signaling networks triggered by mitochondrial perturbation.

Mitochondrial perturbation is a key event in chemical-induced organ toxicities that is incompletely understood. Here, we studied how electron transport chain (ETC) complex I, II, or III (CI, CII and CIII) inhibitors affect mitochondrial functionality, stress response activation, and cell viability using a combination of high-content imaging and TempO-Seq in HepG2 hepatocyte cells. CI and CIII inhibitors perturbed mitochondrial membrane potential (MMP) and mitochondrial and cellular ATP levels in a concentration- and time-dependent fashion and, under conditions preventing a switch to glycolysis attenuated cell viability, whereas CII inhibitors had no effect. TempO-Seq analysis of changes in mRNA expression pointed to a shared cellular response to CI and CIII inhibition. First, to define specific ETC inhibition responses, a gene set responsive toward ETC inhibition (and not to genotoxic, oxidative, or endoplasmic reticulum stress) was identified using targeted TempO-Seq in HepG2. Silencing of one of these genes, NOS3, exacerbated the impact of CI and CIII inhibitors on cell viability, indicating its functional implication in cellular responses to mitochondrial stress. Then by monitoring dynamic responses to ETC inhibition using a HepG2 GFP reporter panel for different classes of stress response pathways and applying pathway and gene network analysis to TempO-Seq data, we looked for downstream cellular events of ETC inhibition and identified the amino acid response (AAR) as being triggered in HepG2 by ETC inhibition. Through in silico approaches we provide evidence indicating that a similar AAR is associated with exposure to mitochondrial toxicants in primary human hepatocytes. Altogether, we (i) unravel quantitative, time- and concentration-resolved cellular responses to mitochondrial perturbation, (ii) identify a gene set associated with adaptation to exposure to active ETC inhibitors, and (iii) show that ER stress and an AAR accompany ETC inhibition in HepG2 and primary hepatocytes.

Heterogeneous, delayed-onset killing by multiple-hitting T cells: Stochastic simulations to assess methods for analysis of imaging data.

Although quantitative insights into the killing behaviour of Cytotoxic T Lymphocytes (CTLs) are necessary for the rational design of immune-based therapies, CTL killing function remains insufficiently characterised. One established model of CTL killing treats CTL cytotoxicity as a Poisson process, based on the assumption that CTLs serially kill antigen-presenting target cells via delivery of lethal hits, each lethal hit corresponding to a single injection of cytotoxic proteins into the target cell cytoplasm. Contradicting this model, a recent in vitro study of individual CTLs killing targets over a 12-hour period found significantly greater heterogeneity in CTL killing performance than predicted by Poisson-based killing. The observed killing process was dynamic and varied between CTLs, with the best performing CTLs exhibiting a marked increase in killing during the final hours of the experiments, along with a "burst killing" kinetic. Despite a search for potential differences between CTLs, no mechanistic explanation for the heterogeneous killing kinetics was found. Here we have used stochastic simulations to assess whether target cells might require multiple hits from CTLs before undergoing apoptosis, in order to verify whether multiple-hitting could explain the late onset, burst killing dynamics observed in vitro. We found that multiple-hitting from CTLs was entirely consistent with the observed killing kinetics. Moreover, the number of available targets and the spatiotemporal kinetics of CTL:target interactions influenced the realised CTL killing rate. We subsequently used realistic, spatial simulations to assess methods for estimating the hitting rate and the number of hits required for target death, to be applied to microscopy data of individual CTLs killing targets. We found that measuring the cumulative duration of individual contacts that targets have with CTLs would substantially improve accuracy when estimating the killing kinetics of CTLs.

Heritable tumor cell division rate heterogeneity induces clonal dominance.

Tumors consist of a hierarchical population of cells that differ in their phenotype and genotype. This hierarchical organization of cells means that a few clones (i.e., cells and several generations of offspring) are abundant while most are rare, which is called clonal dominance. Such dominance also occurred in published in vitro iterated growth and passage experiments with tumor cells in which genetic barcodes were used for lineage tracing. A potential source for such heterogeneity is that dominant clones derive from cancer stem cells with an unlimited self-renewal capacity. Furthermore, ongoing evolution and selection within the growing population may also induce clonal dominance. To understand how clonal dominance developed in the iterated growth and passage experiments, we built a computational model that accurately simulates these experiments. The model simulations reproduced the clonal dominance that developed in in vitro iterated growth and passage experiments when the division rates vary between cells, due to a combination of initial variation and of ongoing mutational processes. In contrast, the experimental results can neither be reproduced with a model that considers random growth and passage, nor with a model based on cancer stem cells. Altogether, our model suggests that in vitro clonal dominance develops due to selection of fast-dividing clones.

A systematic analysis of Nrf2 pathway activation dynamics during repeated xenobiotic exposure.

Oxidative stress leads to the activation of the Nuclear factor-erythroid-2-related factor 2 (Nrf2) pathway. While most studies have focused on the activation of the Nrf2 pathway after single chemical treatment, little is known about the dynamic regulation of the Nrf2 pathway in the context of repeated exposure scenarios. Here we employed single cell live imaging to quantitatively monitor the dynamics of the Nrf2 pathway during repeated exposure, making advantage of two HepG2 fluorescent protein reporter cell lines, expressing GFP tagged Nrf2 or sulfiredoxin 1 (Srxn1), a direct downstream target of Nrf2. High throughput live confocal imaging was used to measure the temporal dynamics of these two components of the Nrf2 pathway after repeated exposure to an extensive concentration range of diethyl maleate (DEM) and tert-butylhydroquinone (tBHQ). Single treatment with DEM or tBHQ induced Nrf2 and Srxn1 over time in a concentration-dependent manner. The Nrf2 response to a second treatment was lower than the response to the first exposure with the same concentration, indicating that the response is adaptive. Moreover, a limited fraction of individual cells committed themselves into the Nrf2 response during the second treatment. Despite the suppression of the Nrf2 pathway, the second treatment resulted in a three-fold higher Srxn1-GFP response compared to the first treatment, with all cells participating in the response. While after the first treatment Srxn1-GFP response was linearly related to Nrf2-GFP nuclear translocation, such a linear relationship was less clear for the second exposure. siRNA-mediated knockdown demonstrated that the second response is dependent on the activity of Nrf2. Several other, clinically relevant, compounds (i.e., sulphorophane, nitrofurantoin and CDDO-Me) also enhanced the induction of Srxn1-GFP upon two consecutive repeated exposure. Together the data indicate that adaptation towards pro-oxidants lowers the Nrf2 activation capacity, but simultaneously primes cells for the enhancement of an antioxidant response which depends on factors other than just Nrf2. These data provide further insight in the overall dynamics of stress pathway activation after repeated exposure and underscore the complexity of responses that may govern repeated dose toxicity.

A Sigmoid Functional Response Emerges When Cytotoxic T Lymphocytes Start Killing Fresh Target Cells.

Cytotoxic T lymphocyte (CTL)-mediated killing involves the formation of a synapse with a target cell, followed by delivery of perforin and granzymes. Previously, we derived a general functional response for CTL killing while considering that CTLs form stable synapses (i.e., single-stage) and that the number of conjugates remains at steady state. However, the killing of target cells sometimes requires multiple engagements (i.e., multistage). To study how multistage killing and a lack of steady state influence the functional response, we here analyze a set of differential equations as well as simulations employing the cellular Potts model, in both cases describing CTLs that kill target cells. We find that at steady state the total killing rate (i.e., the number of target cells killed by all CTLs) is well described by the previously derived double saturation function. Compared to single-stage killing, the total killing rate during multistage killing saturates at higher CTL and target cell densities. Importantly, when the killing is measured before the steady state is approached, a qualitatively different functional response emerges for two reasons: First, the killing signal of each CTL gets diluted over several targets and because this dilution effect is strongest at high target cell densities; this can result in a peak in the dependence of the total killing rate on the target cell density. Second, the total killing rate exhibits a sigmoid dependence on the CTL density when killing is a multistage process, because it takes typically more than one CTL to kill a target. In conclusion, a sigmoid dependence of the killing rate on the CTLs during initial phases of killing may be indicative of a multistage killing process. Observation of a sigmoid functional response may thus arise from a dilution effect and is not necessarily due to cooperative behavior of the CTLs.

Subtle CXCR3-Dependent Chemotaxis of CTLs within Infected Tissue Allows Efficient Target Localization.

It is well established how effector T cells exit the vasculature to enter the peripheral tissues in which an infection is ongoing. However, less is known regarding how CTLs migrate toward infected cells after entry into peripheral organs. Recently, it was shown that the chemokine receptor CXCR3 on T cells has an important role in their ability to localize infected cells and to control vaccinia virus infection. However, the search strategy of T cells for virus-infected targets has not been investigated in detail and could involve chemotaxis toward infected cells, chemokinesis (i.e., increased motility) combined with CTL arrest when targets are detected, or both. In this study, we describe and analyze the migration of CTLs within HSV-1-infected epidermis in vivo. We demonstrate that activated T cells display a subtle distance-dependent chemotaxis toward clusters of infected cells and confirm that this is mediated by CXCR3 and its ligands. Although the chemotactic migration is weak, computer simulations based on short-term experimental data, combined with subsequent long-term imaging indicate that this behavior is crucial for efficient target localization and T cell accumulation at effector sites. Thus, chemotactic migration of effector T cells within peripheral tissue forms an important factor in the speed with which T cells are able to arrive at sites of infection.

Reproducibility of Illumina platform deep sequencing errors allows accurate determination of DNA barcodes in cells.

BACKGROUND: Next generation sequencing (NGS) of amplified DNA is a powerful tool to describe genetic heterogeneity within cell populations that can both be used to investigate the clonal structure of cell populations and to perform genetic lineage tracing. For applications in which both abundant and rare sequences are biologically relevant, the relatively high error rate of NGS techniques complicates data analysis, as it is difficult to distinguish rare true sequences from spurious sequences that are generated by PCR or sequencing errors. This issue, for instance, applies to cellular barcoding strategies that aim to follow the amount and type of offspring of single cells, by supplying these with unique heritable DNA tags. RESULTS: Here, we use genetic barcoding data from the Illumina HiSeq platform to show that straightforward read threshold-based filtering of data is typically insufficient to filter out spurious barcodes. Importantly, we demonstrate that specific sequencing errors occur at an approximately constant rate across different samples that are sequenced in parallel. We exploit this observation by developing a novel approach to filter out spurious sequences. CONCLUSIONS: Application of our new method demonstrates its value in the identification of true sequences amongst spurious sequences in biological data sets.

Random migration and signal integration promote rapid and robust T cell recruitment.

To fight infections, rare T cells must quickly home to appropriate lymph nodes (LNs), and reliably localize the antigen (Ag) within them. The first challenge calls for rapid trafficking between LNs, whereas the second may require extensive search within each LN. Here we combine simulations and experimental data to investigate which features of random T cell migration within and between LNs allow meeting these two conflicting demands. Our model indicates that integrating signals from multiple random encounters with Ag-presenting cells permits reliable detection of even low-dose Ag, and predicts a kinetic feature of cognate T cell arrest in LNs that we confirm using intravital two-photon data. Furthermore, we obtain the most reliable retention if T cells transit through LNs stochastically, which may explain the long and widely distributed LN dwell times observed in vivo. Finally, we demonstrate that random migration, both between and within LNs, allows recruiting the majority of cognate precursors within a few days for various realistic infection scenarios. Thus, the combination of two-scale stochastic migration and signal integration is an efficient and robust strategy for T cell immune surveillance.

Tissue-resident memory CD8+ T cells continuously patrol skin epithelia to quickly recognize local antigen.

Recent work has demonstrated that following the clearance of infection a stable population of memory T cells remains present in peripheral organs and contributes to the control of secondary infections. However, little is known about how tissue-resident memory T cells behave in situ and how they encounter newly infected target cells. Here we demonstrate that antigen-specific CD8(+) T cells that remain in skin following herpes simplex virus infection show a steady-state crawling behavior in between keratinocytes. Spatially explicit simulations of the migration of these tissue-resident memory T cells indicate that the migratory dendritic behavior of these cells allows the detection of antigen-expressing target cells in physiologically relevant time frames of minutes to hours. Furthermore, we provide direct evidence for the identification of rare antigen-expressing epithelial cells by skin-patrolling memory T cells in vivo. These data demonstrate the existence of skin patrol by memory T cells and reveal the value of this patrol in the rapid detection of renewed infections at a previously infected site.

A general functional response of cytotoxic T lymphocyte-mediated killing of target cells.

Cytotoxic T lymphocytes (CTLs) kill virus-infected cells and tumor cells, and play a critical role in immune protection. Our knowledge of how the CTL killing efficiency varies with CTL and target cell numbers is limited. Here, we simulate a region of lymphoid tissue using a cellular Potts model to characterize the functional response of CTL killing of target cells, and find that the total killing rate saturates both with the CTL and the target cell densities. The relative saturation in CTL and target cell densities is determined by whether a CTL can kill multiple target cells at the same time, and whether a target cell can be killed by many CTLs together. We find that all the studied regimes can be well described by a double-saturation (DS) function with two different saturation constants. We show that this DS model can be mechanistically derived for the cases where target cells are killed by a single CTL. For the other cases, a biological interpretation of the parameters is still possible. Our results imply that this DS function can be used as a tool to predict the cellular interactions in cytotoxicity data.

B cells within germinal centers migrate preferentially from dark to light zone.

One of the main questions in the field of imaging immune cell migration in living tissues is whether cells fulfill their functionality via random or nonrandom migration processes. For some applications, this issue has remained controversial even after publication of various imaging studies. A prime example is B-cell migration in germinal centers (GCs) where somatic hypermutation and clonal selection of B cells are thought to occur within morphologically distinct regions termed dark zone (DZ) and light zone (LZ). Here, we reanalyze a previously published dataset on GC B-cell migration, applying a sensitive analysis technique to detect directed migration and using a procedure to correct for a number of artifacts that frequently occur in time-lapse imaging experiments. Although B cells roughly perform a persistent random walk, we present evidence that they have a small preference (of on average about 0.2-0.3 µm min(-1)) to migrate from DZ to LZ, which is consistent with classical views of the GC reaction. This preference is most pronounced among a large subset of almost half of the B-cell population migrating along relatively straight tracks. Using a computational model to generate long-lasting B-cell tracks based on the experimental motility data (including the small directional preference), we predict a time course to travel from DZ to LZ of a few hours. This is consistent with experimental observations, and we show that at the observed cellular motility such a time course cannot be explained without the small preferential migration from DZ to LZ.

Lymph node topology dictates T cell migration behavior.

Adaptive immunity is initiated by T cell recognition of foreign peptides presented on dendritic cells (DCs) by major histocompatibility molecules. These interactions take place in secondary lymphoid tissues, such as lymph nodes (LNs) and spleen, and hence the anatomical structure of these tissues plays a crucial role in the development of immune responses. Two-photon microscopy (2PM) imaging in LNs suggests that T cells walk in a consistent direction for several minutes, pause briefly with a regular period, and then take off in a new, random direction. Here, we construct a spatially explicit model of T cell and DC migration in LNs and show that all dynamical properties of T cells could be a consequence of the densely packed LN environment. By means of 2PM experiments, we confirm that the large velocity fluctuations of T cells are indeed environmentally determined rather than resulting from an intrinsic motility program. Our simulations further predict that T cells self-organize into microscopically small, highly dynamic streams. We present experimental evidence for the presence of such turbulent streams in LNs. Finally, the model allows us to estimate the scanning rates of DCs (2,000 different T cells per hour) and T cells (100 different DCs per hour).

Chemotactic migration of T cells towards dendritic cells promotes the detection of rare antigens.

In many immunological processes chemoattraction is thought to play a role in guiding cells to their sites of action. However, based on in vivo two-photon microscopy experiments in the absence of cognate antigen, T cell migration in lymph nodes (LNs) has been roughly described as a random walk. Although it has been shown that dendritic cells (DCs) carrying cognate antigen in some circumstances attract T cells chemotactically, it is currently still unclear whether chemoattraction of T cells towards DCs helps or hampers scanning. Chemoattraction towards DCs could on the one hand help T cells to rapidly find DCs. On the other hand, it could be deleterious if DCs become shielded by a multitude of attracted yet non-specific T cells. Results from a recent simulation study suggested that the deleterious effect dominates. We re-addressed the question whether T cell chemoattraction towards DCs is expected to promote or hamper the detection of rare antigens using the Cellular Potts Model, a formalism that allows for dynamic, flexible cellular shapes and cell migration. Our simulations show that chemoattraction of T cells enhances the DC scanning efficiency, leading to an increased probability that rare antigen-specific T cells find DCs carrying cognate antigen. Desensitization of T cells after contact with a DC further improves the scanning efficiency, yielding an almost threefold enhancement compared to random migration. Moreover, the chemotaxis-driven migration still roughly appears as a random walk, hence fine-tuned analysis of cell tracks will be required to detect chemotaxis within microscopy data.