Search

A single dose of psilocybin, a psychedelic that acutely causes distortions of space-time perception and ego dissolution, produces rapid and persistent therapeutic effects in human clinical trials1-4. In animal models, psilocybin induces neuroplasticity in cortex and hippocampus5-8. It remains unclear how human brain network changes relate to subjective and lasting effects of psychedelics. Here we tracked individual-specific brain changes with longitudinal precision functional mapping (roughly 18 magnetic resonance imaging visits per participant). Healthy adults were tracked before, during and for 3 weeks after high-dose psilocybin (25 mg) and methylphenidate (40 mg), and brought back for an additional psilocybin dose 6-12 months later. Psilocybin massively disrupted functional connectivity (FC) in cortex and subcortex, acutely causing more than threefold greater change than methylphenidate. These FC changes were driven by brain desynchronization across spatial scales (areal, global), which dissolved network distinctions by reducing correlations within and anticorrelations between networks. Psilocybin-driven FC changes were strongest in the default mode network, which is connected to the anterior hippocampus and is thought to create our sense of space, time and self. Individual differences in FC changes were strongly linked to the subjective psychedelic experience. Performing a perceptual task reduced psilocybin-driven FC changes. Psilocybin caused persistent decrease in FC between the anterior hippocampus and default mode network, lasting for weeks. Persistent reduction of hippocampal-default mode network connectivity may represent a neuroanatomical and mechanistic correlate of the proplasticity and therapeutic effects of psychedelics.

Microglia, the immune cells of the central nervous system, are dynamic and heterogenous cells. While single cell RNA sequencing has become the conventional methodology for evaluating microglial state, transcriptomics do not provide insight into functional changes, identifying a critical gap in the field. Here, we propose a novel organelle phenotyping approach in which we treat live human induced pluripotent stem cell-derived microglia (iMGL) with organelle dyes staining mitochondria, lipids, lysosomes and acquire data by live-cell spectral microscopy. Dimensionality reduction techniques and unbiased cluster identification allow for recognition of microglial subpopulations with single-cell resolution based on organelle function. We validated this methodology using lipopolysaccharide and IL-10 treatment to polarize iMGL to an "inflammatory" and "anti-inflammatory" state, respectively, and then applied it to identify a novel regulator of iMGL function, complement protein C1q. While C1q is traditionally known as the initiator of the complement cascade, here we use organelle phenotyping to identify a role for C1q in regulating iMGL polarization via fatty acid storage and mitochondria membrane potential. Follow up evaluation of microglia using traditional read outs of activation state confirm that C1q drives an increase in microglia pro-inflammatory gene production and migration, while suppressing microglial proliferation. These data together validate the use of a novel organelle phenotyping approach and enable better mechanistic investigation of molecular regulators of microglial state.

In this work, we examine the use of environment-sensitive fluorescent dyes in fluorescence lifetime imaging microscopy (FLIM) biosensors. We screened merocyanine dyes to find an optimal combination of environment-induced lifetime changes, photostability, and brightness at wavelengths suitable for live-cell imaging. FLIM was used to monitor a biosensor reporting conformational changes of endogenous Cdc42 in living cells. The ability to quantify activity using phasor analysis of a single fluorophore (e.g., rather than ratio imaging) eliminated potential artifacts. We leveraged these properties to determine specific concentrations of activated Cdc42 across the cell.

Fluorescence microscopy can provide valuable information about cell interior dynamics. Particularly, mean squared displacement (MSD) analysis is widely used to characterize proteins and sub-cellular structures' mobility providing the laws of molecular diffusion. The MSD curve is traditionally extracted from individual trajectories recorded by single-particle tracking-based techniques. More recently, image correlation methods like iMSD have been shown capable of providing averaged dynamic information directly from images, without the need for isolation and localization of individual particles. iMSD is a powerful technique that has been successfully applied to many different biological problems, over a wide spatial and temporal scales. The aim of this work is to review and compare these two well-established methodologies and their performance in different situations, to give an insight on how to make the most out of their unique characteristics. We show the analysis of the same datasets by the two methods. Regardless of the experimental differences in the input data for MSD or iMSD analysis, our results show that the two approaches can address equivalent questions for free diffusing systems. We focused on studying a range of diffusion coefficients between D = 0.001μm2s-1and D = 0.1μm2s-1, where we verified that the equivalence is maintained even for the case of isolated particles. This opens new opportunities for studying intracellular dynamics using equipment commonly available in any biophysical laboratory.

Fluorescence lifetime imaging microscopy (FLIM) has become a standard tool in the quantitative characterization of subcellular environments. However, quantitative FLIM analyses face several challenges. First, spatial correlations between pixels are often ignored as signal from individual pixels is analyzed independently thereby limiting spatial resolution. Second, existing methods deduce photon ratios instead of absolute lifetime maps. Next, the number of fluorophore species contributing to the signal is unknown, while excited state lifetimes with

The cell membrane has a fundamental role in the cell life cycle but there's still much to be learned about its heterogeneous structure, regulation, and protein interaction. Additionally, the protein-membrane interaction is often overlooked when studying specific protein dynamics. In this work, we present a new tool for a better understanding of protein dynamics and membrane function using live cells and fast non-invasive techniques without the need for individual particle tracking. To this end, we used the 2D-pair correlation function (2D-pCF) to study protein interactions across cellular membranes. We performed numerical simulations and confocal experiments using a GAP-mEGFP fusion construct known to interact with the plasmatic membrane. Our results demonstrate that based on a quantitative correlation analysis as the 2D pair correlation of the signal intensities, is possible to characterize protein-membrane interactions in live systems and real-time. Combining experimental and numerical results this work presents a new powerful approach to the study of the dynamic protein-membrane interaction.

Chromosomal instability (CIN) and epigenetic alterations are characteristics of advanced and metastatic cancers1-4, but whether they are mechanistically linked is unknown. Here we show that missegregation of mitotic chromosomes, their sequestration in micronuclei5,6 and subsequent rupture of the micronuclear envelope7 profoundly disrupt normal histone post-translational modifications (PTMs), a phenomenon conserved across humans and mice, as well as in cancer and non-transformed cells. Some of the changes in histone PTMs occur because of the rupture of the micronuclear envelope, whereas others are inherited from mitotic abnormalities before the micronucleus is formed. Using orthogonal approaches, we demonstrate that micronuclei exhibit extensive differences in chromatin accessibility, with a strong positional bias between promoters and distal or intergenic regions, in line with observed redistributions of histone PTMs. Inducing CIN causes widespread epigenetic dysregulation, and chromosomes that transit in micronuclei experience heritable abnormalities in their accessibility long after they have been reincorporated into the primary nucleus. Thus, as well as altering genomic copy number, CIN promotes epigenetic reprogramming and heterogeneity in cancer.

Motor cortex (M1) has been thought to form a continuous somatotopic homunculus extending down the precentral gyrus from foot to face representations1,2, despite evidence for concentric functional zones3 and maps of complex actions4. Here, using precision functional magnetic resonance imaging (fMRI) methods, we find that the classic homunculus is interrupted by regions with distinct connectivity, structure and function, alternating with effector-specific (foot, hand and mouth) areas. These inter-effector regions exhibit decreased cortical thickness and strong functional connectivity to each other, as well as to the cingulo-opercular network (CON), critical for action5 and physiological control6, arousal7, errors8 and pain9. This interdigitation of action control-linked and motor effector regions was verified in the three largest fMRI datasets. Macaque and pediatric (newborn, infant and child) precision fMRI suggested cross-species homologues and developmental precursors of the inter-effector system. A battery of motor and action fMRI tasks documented concentric effector somatotopies, separated by the CON-linked inter-effector regions. The inter-effectors lacked movement specificity and co-activated during action planning (coordination of hands and feet) and axial body movement (such as of the abdomen or eyebrows). These results, together with previous studies demonstrating stimulation-evoked complex actions4 and connectivity to internal organs10 such as the adrenal medulla, suggest that M1 is punctuated by a system for whole-body action planning, the somato-cognitive action network (SCAN). In M1, two parallel systems intertwine, forming an integrate-isolate pattern: effector-specific regions (foot, hand and mouth) for isolating fine motor control and the SCAN for integrating goals, physiology and body movement.

Fluorescence lifetime imaging captures the spatial distribution of chemical species across cellular environments employing pulsed illumination confocal setups. However, quantitative interpretation of lifetime data continues to face critical challenges. For instance, fluorescent species with known in vitro excited-state lifetimes may split into multiple species with unique lifetimes when introduced into complex living environments. What is more, mixtures of species, which may be both endogenous and introduced into the sample, may exhibit 1) very similar lifetimes as well as 2) wide ranges of lifetimes including lifetimes shorter than the instrumental response function or whose duration may be long enough to be comparable to the interpulse window. By contrast, existing methods of analysis are optimized for well-separated and intermediate lifetimes. Here, we broaden the applicability of fluorescence lifetime analysis by simultaneously treating unknown mixtures of arbitrary lifetimes-outside the intermediate, Goldilocks, zone-for data drawn from a single confocal spot leveraging the tools of Bayesian nonparametrics (BNP). We benchmark our algorithm, termed BNP lifetime analysis, using a range of synthetic and experimental data. Moreover, we show that the BNP lifetime analysis method can distinguish and deduce lifetimes using photon counts as small as 500.

Nonalcoholic steatohepatitis (NASH) is the most common chronic liver disease in the US, partly due to the increasing incidence of metabolic syndrome, obesity, and type 2 diabetes. The roles of bile acids and their receptors, such as the nuclear receptor farnesoid X receptor (FXR) and the G protein-coupled receptor TGR5, on the development of NASH are not fully clear. C57BL/6J male mice fed a Western diet (WD) develop characteristics of NASH, allowing determination of the effects of FXR and TGR5 agonists on this disease. Here we show that the FXR-TGR5 dual agonist INT-767 prevents progression of WD-induced hepatic steatosis, inflammation, and fibrosis, as determined by histological and biochemical assays and novel label-free microscopy imaging techniques, including third harmonic generation, second harmonic generation, and fluorescence lifetime imaging microscopy. Furthermore, we show INT-767 decreases liver fatty acid synthesis and fatty acid and cholesterol uptake, as well as liver inflammation. INT-767 markedly changed bile acid composition in the liver and intestine, leading to notable decreases in the hydrophobicity index of bile acids, known to limit cholesterol and lipid absorption. In addition, INT-767 upregulated expression of liver p-AMPK, SIRT1, PGC-1α, and SIRT3, which are master regulators of mitochondrial function. Finally, we found INT-767 treatment reduced WD-induced dysbiosis of gut microbiota. Interestingly, the effects of INT-767 in attenuating NASH were absent in FXR-null mice, but still present in TGR5-null mice. Our findings support treatment and prevention protocols with the dual FXR-TGR5 agonist INT-767 arrest progression of WD-induced NASH in mice mediated by FXR-dependent, TGR5-independent mechanisms.

The excited-state lifetime is an intrinsic property of fluorescent molecules that can be leveraged for multiplexed imaging. An advantage of fluorescence lifetime-based multiplexing is that signals from multiple probes can be gathered simultaneously, whereas traditional spectral fluorescence imaging typically requires multiple images at different excitation and emission wavelengths. Additionally, lifetime and spectra could both be utilized to expand the multiplexing capacity of fluorescence. However, resolving exogenous molecular probes based exclusively on the fluorescence lifetime has been limited by technical challenges in analyzing lifetime data. The phasor approach to lifetime analysis offers a simple, graphical solution that has increasingly been used to assess endogenous cellular autofluorescence to quantify metabolic factors. In this study, we employed the phasor analysis of FLIM to quantitatively resolve three exogenous, antibody-targeted fluorescent probes with similar spectral properties based on lifetime information alone. First, we demonstrated that three biomarkers that were spatially restricted to the cell membrane, cytosol, or nucleus could be accurately distinguished using FLIM and phasor analysis. Next, we successfully resolved and quantified three probes that were all targeted to cell surface biomarkers. Finally, we demonstrated that lifetime-based quantitation accuracy can be improved through intensity matching of various probe-biomarker combinations, which will expand the utility of this technique. Importantly, we reconstructed images for each individual probe, as well as an overlay of all three probes, from a single FLIM image. Our results demonstrate that FLIM and phasor analysis can be leveraged as a powerful tool for simultaneous detection of multiple biomarkers with high sensitivity and accuracy.

Aβ oligomers play a central role in the neurodegeneration observed with Alzheimer's disease. Our laboratory has developed covalently stabilized trimers derived from residues 17-36 of Aβ as model systems for studying Aβ oligomers. In the current study, we apply the emerging techniques of fluorescence lifetime imaging microscopy (FLIM) and native mass spectrometry (native MS) to better understand the assembly and interactions of the oligomer model system 2AT-L in aqueous solutions and with cells. 2AT-L and fluorescently labeled 2AT-L analogues assemble in the membrane-like environment of SDS-PAGE, showing diffuse bands of oligomers in equilibrium. Native ion mobility-mass spectrometry (native IM-MS) of 2AT-L allows for the identification of discrete oligomers in solution and shows similar patterns of oligomer formation between 2AT-L and fluorescently labeled analogues. Fluorescence microscopy with SH-SY5Y cells reveals that fluorescently labeled 2AT-L analogues colocalize within lysosomes. FLIM studies with phasor analysis further elucidate the assembly of 2AT-L within cells and establish the occurrence of FRET, indicating the presence of oligomers within cells. Collectively, these multiple complementary techniques help better understand the complex behavior of the 2AT-L model system.

Teixobactin has been the source of intensive study and interest as a promising antibiotic, because of its excellent activity against drug-resistant Gram-positive pathogens and its novel but not yet fully understood mechanism of action that precludes drug resistance. Recent studies have demonstrated that the mode of action of teixobactin is more complicated than initially thought, with supramolecular assembly of the antibiotic appearing to play a critical role in the binding process. Further studies of the interactions of teixobactin with bacteria and its molecular targets offer the promise of providing deeper insights into its novel mechanism of action and guiding the design of additional drug candidates and analogues. The current study reports the preparation and study of teixobactin analogues bearing a variety of fluorophores. Structured illumination microscopy of the fluorescent teixobactin analogues with B. subtilis enables super-resolution visualization of the interaction of teixobactin with bacterial cell walls and permits the observation of aggregated clusters of the antibiotic on the bacteria. Förster resonance energy transfer (FRET) microscopy further elucidates the supramolecular assembly by showing that fluorescent teixobactin molecules co-localize within a few nanometers on B. subtilis. Fluorescence microscopy over time with a fluorescent teixobactin analogue and propidium iodide in B. subtilis reveals a correlation between cell death and binding of the antibiotic to cellular targets, followed by lysis of cells. Collectively, these studies provide new insights into the binding of teixobactin to Gram-positive bacteria, its supramolecular mechanism of action, and the lysis of bacteria that follows.

The original version of this Article contained an error in the order of the Figures. Figure 6 was published as Figure 7 and Figure 7 was published as Figure 6. The Figure legends were correct. The original Figures 6 and 7 and accompanying legends appear below. The original Article has been corrected.

Nuclear proteins can modulate their DNA binding activity and the exploration volume available during DNA target search by self-associating into higher-order oligomers. Directly tracking this process in the nucleoplasm of a living cell is, however, a complex task. Thus, here we present a microscopy method based on radial pair correlation of molecular brightness fluctuations (radial pCOMB) that can extract the mobility of a fluorescently tagged nuclear protein as a function of its oligomeric state and spatiotemporally map the anisotropy of this parameter with respect to nuclear architecture. By simply performing a rapid frame scan acquisition, radial pCOMB has the capacity to detect, within each pixel, protein oligomer formation and the size-dependent obstruction nuclear architecture imparts on this complex's transport across sub-micrometer distances. From application of radial pCOMB to an oligomeric transcription factor and DNA repair protein, we demonstrate that homo-oligomer formation differentially regulates chromatin accessibility and interaction with the DNA template.

Fluorescence lifetime imaging microscopy (FLIM) may reveal subcellular spatial lifetime maps of key molecular species. Yet, such a quantitative picture of life necessarily demands high photon budgets at every pixel under the current analysis paradigm, thereby increasing acquisition time and photodamage to the sample. Motivated by recent developments in computational statistics, we provide a direct means to update our knowledge of the lifetime maps of species of different lifetimes from direct photon arrivals, while accounting for experimental features such as arbitrary forms of the instrument response function (IRF) and exploiting information from empty laser pulses not resulting in photon detection. Our ability to construct lifetime maps holds for arbitrary lifetimes, from short lifetimes (comparable to the IRF) to lifetimes exceeding interpulse times. As our method is highly data efficient, for the same amount of data normally used to determine lifetimes and photon ratios, working within the Bayesian paradigm, we report direct blind unmixing of lifetimes with subnanosecond resolution and subpixel spatial resolution using standard raster scan FLIM images. We demonstrate our method using a wide range of simulated and experimental data.

Macromolecular crowding is crucial for cellular homeostasis. In vivo studies of macromolecular crowding and water dynamics are needed to understand their roles in cellular physiology and fate determination. Macromolecular crowding in the lens is essential for normal optics, and an understanding of its regulation will help prevent cataract and presbyopia. Here, we combine the use of the nanoenvironmental sensor [6-acetyl-2-dimethylaminonaphthalene (ACDAN)] to visualize lens macromolecular crowding with in vivo studies of aquaporin 0 zebrafish mutants that disrupt its regulation. Spectral phasor analysis of ACDAN fluorescence reveals water dipolar relaxation and demonstrates that mutations in two zebrafish aquaporin 0s, Aqp0a and Aqp0b, alter water state and macromolecular crowding in living lenses. Our results provide in vivo evidence that Aqp0a promotes fluid influx in the deeper lens cortex, whereas Aqp0b facilitates fluid efflux. This evidence reveals previously unidentified spatial regulation of macromolecular crowding and spatially distinct roles for Aqp0 in the lens.

How do children within the autism spectrum (ASD) integrate linguistic and visual information during real-time processing of sentences and dynamic scenes? We employed a visual word paradigm with sentences containing different verb types (causatives such as 'crack'; perception/psychological verbs such as 'examine') and with different scene contexts (agent in the scene moving towards a target object or remaining neutral). Participants were a group of 10-year old ASD and typically developing children (TD). We reasoned that the two dominant theories about ASD cognitive processing style—the weak central coherence theory (WCC) and the enhanced perceptual function (EPF)—make different predictions on how and where information about linguistic and visual information may be integrated. WCC predicts that ASD children would have greater difficulty with sentence and scene semantics while EPF predicts an early, post-perceptual integration via verb-specific and object (but not full scene) representations. The results suggest that ASD children differ from TD children by computing the link between verb and object faster, suggesting an enhanced post-perceptual integration.

Multiplexed mRNA profiling in the spatial context provides new information enabling basic research and clinical applications. Unfortunately, existing spatial transcriptomics methods are limited due to either low multiplexing or complexity. Here, we introduce a spatialomics technology, termed Multi Omic Single-scan Assay with Integrated Combinatorial Analysis (MOSAICA), that integrates in situ labeling of mRNA and protein markers in cells or tissues with combinatorial fluorescence spectral and lifetime encoded probes, spectral and time-resolved fluorescence imaging, and machine learning-based decoding. We demonstrate MOSAICA's multiplexing scalability in detecting 10-plex targets in fixed colorectal cancer cells using combinatorial labeling of five fluorophores with facile error-detection and removal of autofluorescence. MOSAICA's analysis is strongly correlated with sequencing data (Pearson's r = 0.96) and was further benchmarked using RNAscopeTM and LGC StellarisTM. We further apply MOSAICA for multiplexed analysis of clinical melanoma Formalin-Fixed Paraffin-Embedded (FFPE) tissues. We finally demonstrate simultaneous co-detection of protein and mRNA in cancer cells.

UTX (also known as KDM6A) encodes a histone H3K27 demethylase and is an important tumour suppressor that is frequently mutated in human cancers1. However, as the demethylase activity of UTX is often dispensable for mediating tumour suppression and developmental regulation2-8, the underlying molecular activity of UTX remains unknown. Here we show that phase separation of UTX underlies its chromatin-regulatory activity in tumour suppression. A core intrinsically disordered region (cIDR) of UTX forms phase-separated liquid condensates, and cIDR loss caused by the most frequent cancer mutation of UTX is mainly responsible for abolishing tumour suppression. Deletion, mutagenesis and replacement assays of the intrinsically disordered region demonstrate a critical role of UTX condensation in tumour suppression and embryonic stem cell differentiation. As shown by reconstitution in vitro and engineered systems in cells, UTX recruits the histone methyltransferase MLL4 (also known as KMT2D) to the same condensates and enriches the H3K4 methylation activity of MLL4. Moreover, UTX regulates genome-wide histone modifications and high-order chromatin interactions in a condensation-dependent manner. We also found that UTY, the Y chromosome homologue of UTX with weaker tumour-suppressive activity, forms condensates with reduced molecular dynamics. These studies demonstrate a crucial biological function of liquid condensates with proper material states in enabling the tumour-suppressive activity of a chromatin regulator.

Current methods for the diagnosis of urinary tract infections with antimicrobial susceptibility testing take 2-3 days and require a clinical laboratory. The lack of a rapid, point-of-care antibiotic susceptibility test (AST) has contributed to the misuse of antibiotics when treating urinary tract infections (UTIs) and consequently to the rise of multi-drug-resistant organisms. The current clinical approach has led to reduced treatment options and increased costs of diagnosis and therapy. To address this issue, novel diagnostics are needed for the timely determination of antimicrobial susceptibility. We present a rapid, point-of-care, phenotypic AST device that can report the antibiotic susceptibility/resistance of a uropathogen to a panel of antibiotics in as few as 3 h by utilizing fluorescent-labelling chemistry and a highly sensitive particle-counting instrument. We analysed 744 urine samples from the outpatients and inpatients of two Italian hospitals. The 130 UTI-positive patient urine samples we found were measured using a panel of six common UTI antibiotics plus a growth control. By comparing our results to hospital laboratory urine cultures, we obtained an overall sensitivity = 81%, a specificity = 83%, an SPV (sensitivity predicted value) = 95%, and an RPV (resistance predicted value) = 54%. According to our preliminary data, the sensitivity predicted value for a single antibiotic agent was 95%, thus allowing (in the vast majority of cases) an early (within 3 h) recognition of an effective agent for a single patient.

The hippocampus is critically important for a diverse range of cognitive processes, such as episodic memory, prospective memory, affective processing, and spatial navigation. Using individual-specific precision functional mapping of resting-state functional MRI data, we found the anterior hippocampus (head and body) to be preferentially functionally connected to the default mode network (DMN), as expected. The hippocampal tail, however, was strongly preferentially functionally connected to the parietal memory network (PMN), which supports goal-oriented cognition and stimulus recognition. This anterior-posterior dichotomy of resting-state functional connectivity was well-matched by differences in task deactivations and anatomical segmentations of the hippocampus. Task deactivations were localized to the hippocampal head and body (DMN), relatively sparing the tail (PMN). The functional dichotomization of the hippocampus into anterior DMN-connected and posterior PMN-connected parcels suggests parallel but distinct circuits between the hippocampus and medial parietal cortex for self- versus goal-oriented processing.

In the original version, the bounds given in Eqs. (87a) and (87b) on the contribution to the early-time optical depth, (15,30), contained a numerical error in deriving the 95th percentile from the Monte Carlo samples. The corrected 95% upper bounds are: τ(15,30) < 0:018 (lowE, flat τ(15, 30), FlexKnot), (1) τ(15, 30) < 0:023 (lowE, flat knot, FlexKnot): (2) These bounds are a factor of 3 larger than the originally reported results. Consequently, the new bounds do not significantly improve upon previous results from Planck data presented in Millea & Bouchet (2018) as was stated, but are instead comparable. Equations (1) and (2) give results that are now similar to those of Heinrich & Hu (2021), who used the same Planck 2018 data to derive a 95% upper bound of 0.020 using the principal component analysis (PCA) model and uniform priors on the PCA mode amplitudes.

In the original version, the bounds given in Eqs. (87a) and (87b) on the contribution to the early-time optical depth, (15,30), contained a numerical error in deriving the 95th percentile from the Monte Carlo samples. The corrected 95% upper bounds are: τ(15,30) < 0:018 (lowE, flat τ(15, 30), FlexKnot), (1) τ(15, 30) < 0:023 (lowE, flat knot, FlexKnot): (2) These bounds are a factor of 3 larger than the originally reported results. Consequently, the new bounds do not significantly improve upon previous results from Planck data presented in Millea & Bouchet (2018) as was stated, but are instead comparable. Equations (1) and (2) give results that are now similar to those of Heinrich & Hu (2021), who used the same Planck 2018 data to derive a 95% upper bound of 0.020 using the principal component analysis (PCA) model and uniform priors on the PCA mode amplitudes.

We describe a method based on a pair of transmission filters placed in the emission path of a microscope to resolve the emission wavelength of every point in an image. The method can be applied to any type of imaging device that provides the light in the wavelength transmission range of the filters. Unique characteristics of the filter approach are that the light does not need to be collimated and the wavelength response does not depend on the scattering of the sample or tissue. The pair of filters are used to produce the spectral phasor of the transmitted light, which is sufficient to perform spectral deconvolution over a broad wavelength range. The method is sensitive enough to distinguish free and protein-bound NADH and can be used in metabolic studies.

Hyperspectral imaging is highly sought after in many fields including mineralogy and geology, environment and agriculture, astronomy and, importantly, biomedical imaging and biological fluorescence. We developed ultrafast phasor-based hyperspectral snapshot microscopy based on sine/cosine interference filters for biomedical imaging not feasible with conventional hyperspectral detection methods. Current approaches rely on slow spatial or spectral scanning limiting their application in living biological tissues, while faster snapshot methods such as image mapping spectrometry and multispectral interferometry are limited in spatial and/or spectral resolution, are computationally demanding, and imaging devices are very expensive to manufacture. Leveraging light sheet microscopy, phasor-based hyperspectral snapshot microscopy improved imaging speed 10-100 fold which, combined with minimal light exposure and high detection efficiency, enabled hyperspectral metabolic imaging of live, three-dimensional mouse tissues not feasible with other methods. As a fit-free method that does not require any a priori information often unavailable in complex and evolving biological systems, the rule of linear combinations of the phasor could spectrally resolve subtle differences between cell types in the developing zebrafish retina and spectrally separate and track multiple organelles in 3D cultured cells over time. The sine/cosine snapshot method is adaptable to any microscope or imaging device thus making hyperspectral imaging and fit-free analysis based on linear combinations broadly available to researchers and the public.

The phasor approach is a well-established method for data visualization and image analysis in spectral and lifetime fluorescence microscopy. Nevertheless, it is typically applied in a user-dependent manner by manually selecting regions of interest on the phasor space to find distinct regions in the fluorescence images. In this paper we present our work on using machine learning clustering techniques to establish an unsupervised and automatic method that can be used for identifying populations of fluorescent species in spectral and lifetime imaging. We demonstrate our method using both synthetic data, created by sampling photon arrival times and plotting the distributions on the phasor plot, and real live cells samples, by staining cellular organelles with a selection of commercial probes.

Fluorescence lifetime imaging microscopy (FLIM) and spectral imaging are two broadly applied methods for increasing dimensionality in microscopy. However, their combination is typically inefficient and slow in terms of acquisition and processing. By integrating technological and computational advances, we developed a robust and unbiased spectral FLIM (S-FLIM) system. Our method, Phasor S-FLIM, combines true parallel multichannel digital frequency domain electronics with a multidimensional phasor approach to extract detailed and precise information about the photophysics of fluorescent specimens at optical resolution. To show the flexibility of the Phasor S-FLIM technology and its applications to the biological and biomedical field, we address four common, yet challenging, problems: the blind unmixing of spectral and lifetime signatures from multiple unknown species, the unbiased bleedthrough- and background-free Förster resonance energy transfer analysis of biosensors, the photophysical characterization of environment-sensitive probes in living cells and parallel four-color FLIM imaging in tumor spheroids.

Any chemist studying the interaction of molecules with lipid assemblies will eventually be confronted by the topic of membrane bilayer heterogeneity and may ultimately encounter the heterogeneity of natural membranes. In artificial bilayers, heterogeneity is defined by phase segregation that can be in the nano- and micrometer range. In biological bilayers, heterogeneity is considered in the context of small (10-200 nm) sterol and sphingolipid-enriched heterogeneous and highly dynamic domains. Several techniques can be used to assess membrane heterogeneity in living systems. Our approach is to use a fluorescent reporter molecule immersed in the bilayer, which, by changes in its spectroscopic properties, senses physical-chemistry aspects of the membrane. This dye in combination with microscopy and fluctuation techniques can give information about membrane heterogeneity at different temporal and spatial levels: going from average fluidity to number and diffusion coefficient of nanodomains. LAURDAN (6-dodecanoyl-2-(dimethylamino) naphthalene), is a fluorescent probe designed and synthesized in 1979 by Gregorio Weber with the purpose to study the phenomenon of dipolar relaxation. The spectral displacement observed when LAURDAN is either in fluid or gel phase permitted the use of the technique in the field of membrane dynamics. The quantitation of the spectral displacement was first addressed by the generalized polarization (GP) function in the cuvette, a ratio of the difference in intensity at two wavelengths divided by their sum. In 1997, GP measurements were done for the first time in the microscope, adding to the technique the spatial resolution and allowing the visualization of lipid segregation both in liposomes and cells. A new prospective to the membrane heterogeneity was obtained when LAURDAN fluorescent lifetime measurements were done in the microscope. Two channel lifetime imaging provides information on membrane polarity and dipole relaxation (the two parameters responsible for the spectral shift of LAURDAN), and the application of phasor analysis allows pixel by pixel understanding of these two parameters in the membrane. To increase temporal resolution, LAURDAN GP was combined with fluctuation correlation spectroscopy (FCS) and the motility of nanometric highly packed structures in biological membranes was registered. Lately the application of phasor analysis to spectral images from membranes labeled with LAURDAN allows us to study the full spectra pixel by pixel in an image. All these methodologies, using LAURDAN, offer the possibility to address different properties of membranes depending on the question being asked. In this Account, we will focus on the principles, advantages, and limitations of different approaches to orient the reader to select the most appropriate technique for their research.

Flaviviruses are major human disease-causing pathogens, including dengue virus (DENV), Zika virus, yellow fever virus and others. DENV infects hundreds of millions of people per year around the world, causing a tremendous social and economic burden. DENV capsid (C) protein plays an essential role during genome encapsidation and viral particle formation. It has been previously shown that DENV C enters the nucleus in infected cells. However, whether DENV C protein exhibits nuclear export remains unclear. By spatially cross-correlating different regions of the cell, we investigated DENV C movement across the nuclear envelope during the infection cycle. We observed that transport takes place in both directions and with similar translocation times (in the ms time scale) suggesting a bidirectional movement of both C protein import and export.Furthermore, from the pair cross-correlation functions in cytoplasmic or nuclear regions we found two populations of C molecules in each compartment with fast and slow mobilities. While in the cytoplasm the correlation times were in the 2-6 and 40-110 ms range for the fast and slow mobility populations respectively, in the cell nucleus they were 1-10 and 25-140 ms range, respectively. The fast mobility of DENV C in cytoplasmic and nuclear regions agreed with the diffusion coefficients from Brownian motion previously reported from correlation analysis. These studies provide the first evidence of DENV C shuttling from and to the nucleus in infected cells, opening new venues for antiviral interventions.

Circadian gene expression driven by transcription activators CLOCK and BMAL1 is intimately associated with dynamic chromatin remodeling. However, how cellular metabolism directs circadian chromatin remodeling is virtually unexplored. We report that the S-adenosylhomocysteine (SAH) hydrolyzing enzyme adenosylhomocysteinase (AHCY) cyclically associates to CLOCK-BMAL1 at chromatin sites and promotes circadian transcriptional activity. SAH is a potent feedback inhibitor of S-adenosylmethionine (SAM)-dependent methyltransferases, and timely hydrolysis of SAH by AHCY is critical to sustain methylation reactions. We show that AHCY is essential for cyclic H3K4 trimethylation, genome-wide recruitment of BMAL1 to chromatin, and subsequent circadian transcription. Depletion or targeted pharmacological inhibition of AHCY in mammalian cells markedly decreases the amplitude of circadian gene expression. In mice, pharmacological inhibition of AHCY in the hypothalamus alters circadian locomotor activity and rhythmic transcription within the suprachiasmatic nucleus. These results reveal a previously unappreciated connection between cellular metabolism, chromatin dynamics, and circadian regulation.

We describe an extension of the most recent version of the Planck Catalogue of Compact Sources (PCCS2), produced using a new multi-band Bayesian Extraction and Estimation Package (BeeP). BeeP assumes that the compact sources present in PCCS2 at 857 GHz have a dust-like spectral energy distribution (SED), which leads to emission at both lower and higher frequencies, and adjusts the parameters of the source and its SED to fit the emission observed in Planck's three highest frequency channels at 353, 545, and 857 GHz, as well as the IRIS map at 3000 GHz. In order to reduce confusion regarding diffuse cirrus emission, BeeP's data model includes a description of the background emission surrounding each source, and it adjusts the confidence in the source parameter extraction based on the statistical properties of the spatial distribution of the background emission. BeeP produces the following three new sets of parameters for each source: (a) fits to a modified blackbody (MBB) thermal emission model of the source; (b) SED-independent source flux densities at each frequency considered; and (c) fits to an MBB model of the background in which the source is embedded. BeeP also calculates, for each source, a reliability parameter, which takes into account confusion due to the surrounding cirrus. This parameter can be used to extract sub-samples of high-frequency sources with statistically well-understood properties. We define a high-reliability subset (BeeP/base), containing 26 083 sources (54.1% of the total PCCS2 catalogue), the majority of which have no information on reliability in the PCCS2. We describe the characteristics of this specific high-quality subset of PCCS2 and its validation against other data sets, specifically for: the sub-sample of PCCS2 located in low-cirrus areas; the Planck Catalogue of Galactic Cold Clumps; the Herschel GAMA15-field catalogue; and the temperature-and spectral-index-reconstructed dust maps obtained with Planck's Generalized Needlet Internal Linear Combination method. The results of the BeeP extension of PCCS2, which are made publicly available via the Planck Legacy Archive, will enable the study of the thermal properties of well-defined samples of compact Galactic and extragalactic dusty sources.

We describe an extension of the most recent version of the Planck Catalogue of Compact Sources (PCCS2), produced using a new multi-band Bayesian Extraction and Estimation Package (BeeP). BeeP assumes that the compact sources present in PCCS2 at 857 GHz have a dust-like spectral energy distribution (SED), which leads to emission at both lower and higher frequencies, and adjusts the parameters of the source and its SED to fit the emission observed in Planck's three highest frequency channels at 353, 545, and 857 GHz, as well as the IRIS map at 3000 GHz. In order to reduce confusion regarding diffuse cirrus emission, BeeP's data model includes a description of the background emission surrounding each source, and it adjusts the confidence in the source parameter extraction based on the statistical properties of the spatial distribution of the background emission. BeeP produces the following three new sets of parameters for each source: (a) fits to a modified blackbody (MBB) thermal emission model of the source; (b) SED-independent source flux densities at each frequency considered; and (c) fits to an MBB model of the background in which the source is embedded. BeeP also calculates, for each source, a reliability parameter, which takes into account confusion due to the surrounding cirrus. This parameter can be used to extract sub-samples of high-frequency sources with statistically well-understood properties. We define a high-reliability subset (BeeP/base), containing 26 083 sources (54.1% of the total PCCS2 catalogue), the majority of which have no information on reliability in the PCCS2. We describe the characteristics of this specific high-quality subset of PCCS2 and its validation against other data sets, specifically for: the sub-sample of PCCS2 located in low-cirrus areas; the Planck Catalogue of Galactic Cold Clumps; the Herschel GAMA15-field catalogue; and the temperature-and spectral-index-reconstructed dust maps obtained with Planck's Generalized Needlet Internal Linear Combination method. The results of the BeeP extension of PCCS2, which are made publicly available via the Planck Legacy Archive, will enable the study of the thermal properties of well-defined samples of compact Galactic and extragalactic dusty sources.

The phasor approach is used in fluorescence lifetime imaging microscopy for several purposes, notably to calculate the metabolic index of single cells and tissues. An important feature of the phasor approach is that it is a fit-free method allowing immediate and easy to interpret analysis of images. In a recent paper, we showed that three or four intensity fractions of exponential components can be resolved in each pixel of an image by the phasor approach using simple algebra, provided the component phasors are known. This method only makes use of the rule of linear combination of phasors rather than fits. Without prior knowledge of the components and their single exponential decay times, resolution of components and fractions is much more challenging. Blind decomposition has been carried out only for cuvette experiments wherein the statistics in terms of the number of photons collected is very good. In this paper, we show that using the phasor approach and measurements of the decay at phasor harmonics 2 and 3, available using modern electronics, we could resolve the decay in each pixel of an image in live cells or mice liver tissues with two or more exponential components without prior knowledge of the values of the components. In this paper, blind decomposition is achieved using a graphical method for two components and a minimization method for three components. This specific use of the phasor approach to resolve multicomponents in a pixel enables applications where multiplexing species with different lifetimes and potentially different spectra can provide a different type of super-resolved image content.

Diabetic kidney disease continues to be the leading cause of chronic kidney disease, often advancing to end stage kidney disease. In addition to the well characterized glomerular alterations including mesangial expansion, podocyte injury, and glomerulosclerosis, tubulointerstitial fibrosis is also an important component of diabetic kidney injury. Similarly, tubulointerstitial fibrosis is a critical component of any chronic kidney injury. Therefore, sensitive and quantitative identification of tubulointerstitial fibrosis is critical for the assessment of long-term prognosis of kidney disease. Here, we employed phasor approach to fluorescence lifetime imaging, commonly known as FLIM, to understand tissue heterogeneity and calculate changes in the tissue autofluorescence lifetime signatures due to diabetic kidney disease. FLIM imaging was performed on cryostat sections of snap-frozen biopsy material of patients with diabetic nephropathy. There was an overall increase in phase lifetime (τphase) with increased disease severity. Multicomponent phasor analysis shows the distinctive differences between the different disease states. Thus, phasor autofluorescence lifetime imaging, which does not involve any staining, can be used to understand and evaluate the severity of kidney disease.

Serine incorporator protein 5 (SERINC5) is the host antiretroviral factor that reduces HIV-1 infectivity by incorporating into virions and inhibiting the envelope glycoprotein (Env) mediated virus fusion with target cells. We and others have shown that SERINC5 incorporation into virions alters the Env structure and sensitizes the virus to broadly neutralizing antibodies targeting cryptic Env epitopes. We have also found that SERINC5 accelerates the loss of Env function over time compared to control viruses. However, the exact mechanism by which SERINC5 inhibits HIV-1 fusion is not understood. Here, we utilized 2D and 3D super-resolution microscopy to examine the effect of SERINC5 on the distribution of Env glycoproteins on single HIV-1 particles. We find that, in agreement with a previous report, Env glycoproteins form clusters on the surface of mature virions. Importantly, incorporation of SERINC5, but not SERINC2, which lacks antiviral activity, disrupted Env clusters without affecting the overall Env content. We also show that SERINC5 and SERINC2 also form clusters on single virions. Unexpectedly, Env and SERINC molecules exhibited poor codistribution on virions, as evidenced by much greater Env-SERINC pairwise distances compared to Env-Env distances. This observation is inconsistent with the previously reported interaction between Env and SERINC5 and suggests an indirect effect of SERINC5 on Env cluster formation. Collectively, our results reveal a multifaceted mechanism of SERINC5-mediated restriction of HIV-1 fusion that, aside from the effects on individual Env trimers, involves disruption of Env clusters, which likely serve as sites of viral fusion with target cells.

High-resolution imaging techniques capable of detecting identifiable endogenous fluorophores in the eye along with genetic testing will dramatically improve diagnostic capabilities in the ophthalmology clinic and accelerate the development of new treatments for blinding diseases. Two-photon excitation (TPE)-based imaging overcomes the filtering of ultraviolet light by the lens of the human eye and thus can be utilized to discover defects in vitamin A metabolism during the regeneration of the visual pigments required for the detection of light. Combining TPE with fluorescence lifetime imaging (FLIM) and spectral analyses offers the potential of detecting diseases of the retina at earlier stages before irreversible structural damage has occurred. The main barriers to realizing the benefits of TPE for imaging the human retina arise from concerns about the high light exposure typically needed for informative TPE imaging and the requirement to correlate the ensuing data with different states of health and disease. To overcome these hurdles, we improved TPE efficiency by controlling temporal properties of the excitation light and employed phasor analyses to FLIM and spectral data in mouse models of retinal diseases. Modeling of retinal photodamage revealed that plasma-mediated effects do not play a role and that melanin-related thermal effects are mitigated by reducing pulse repetition frequency. By using noninvasive TPE imaging we identified molecular components of individual granules in the retinal pigment epithelium and present their analytical characteristics.

We present full-sky maps of the cosmic microwave background (CMB) and polarized synchrotron and thermal dust emission, derived from the third set of Planck frequency maps. These products have significantly lower contamination from instrumental systematic effects than previous versions. The methodologies used to derive these maps follow closely those described in earlier papers, adopting four methods (Commander, NILC, SEVEM, and SMICA) to extract the CMB component, as well as three methods (Commander, GNILC, and SMICA) to extract astrophysical components. Our revised CMB temperature maps agree with corresponding products in the Planck 2015 delivery, whereas the polarization maps exhibit significantly lower large-scale power, reflecting the improved data processing described in companion papers; however, the noise properties of the resulting data products are complicated, and the best available end-to-end simulations exhibit relative biases with respect to the data at the few percent level. Using these maps, we are for the first time able to fit the spectral index of thermal dust independently over 3° regions. We derive a conservative estimate of the mean spectral index of polarized thermal dust emission of βd = 1.55 ± 0.05, where the uncertainty marginalizes both over all known systematic uncertainties and different estimation techniques. For polarized synchrotron emission, we find a mean spectral index of βs = -3.1 ± 0.1, consistent with previously reported measurements. We note that the current data processing does not allow for construction of unbiased single-bolometer maps, and this limits our ability to extract CO emission and correlated components. The foreground results for intensity derived in this paper therefore do not supersede corresponding Planck 2015 products. For polarization the new results supersede the corresponding 2015 products in all respects.

We present a final description of the data-processing pipeline for the Planck Low Frequency Instrument (LFI), implemented for the 2018 data release. Several improvements have been made with respect to the previous release, especially in the calibration process and in the correction of instrumental features such as the effects of nonlinearity in the response of the analogue-to-digital converters. We provide a brief pedagogical introduction to the complete pipeline, as well as a detailed description of the important changes implemented. Self-consistency of the pipeline is demonstrated using dedicated simulations and null tests. We present the final version of the LFI full sky maps at 30, 44, and 70 GHz, both in temperature and polarization, together with a refined estimate of the solar dipole and a final assessment of the main LFI instrumental parameters.

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E-mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and optimal modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following final results: flocalNL= -0.9 ± 5.1; fequilNL= -26 ± 47; and forthoNL= -38 ± 24 (68% CL, statistical). These results include low-multipole (4 ≤ ℓ < 40) polarization data that are not included in our previous analysis. The results also pass an extensive battery of tests (with additional tests regarding foreground residuals compared to 2015), and they are stable with respect to our 2015 measurements (with small fluctuations, at the level of a fraction of a standard deviation, which is consistent with changes in data processing). Polarizationonly bispectra display a significant improvement in robustness; they can now be used independently to set primordial NG constraints with a sensitivity comparable to WMAP temperature-based results and they give excellent agreement. In addition to the analysis of the standard local, equilateral, and orthogonal bispectrum shapes, we consider a large number of additional cases, such as scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The nonprimordial lensing bispectrum is, however, detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5σ. Beyond estimates of individual shape amplitudes, we also present model-independent reconstructions and analyses of the Planck CMB bispectrum. Our final constraint on the local primordial trispectrum shape is glocalNL= (-5.8 ± 6.5) × 104(68% CL, statistical), while constraints for other trispectrum shapes are also determined. Exploiting the tight limits on various bispectrum and trispectrum shapes, we constrain the parameter space of different early-Universe scenarios that generate primordial NG, including general single-field models of inflation, multi-field models (e.g. curvaton models), models of inflation with axion fields producing parity-violation bispectra in the tensor sector, and inflationary models involving vector-like fields with directionally-dependent bispectra. Our results provide a high-precision test for structure-formation scenarios, showing complete agreement with the basic picture of the CDM cosmology regarding the statistics of the initial conditions, with cosmic structures arising from adiabatic, passive, Gaussian, and primordial seed perturbations.

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E-mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and optimal modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following final results: flocalNL= -0.9 ± 5.1; fequilNL= -26 ± 47; and forthoNL= -38 ± 24 (68% CL, statistical). These results include low-multipole (4 ≤ ℓ < 40) polarization data that are not included in our previous analysis. The results also pass an extensive battery of tests (with additional tests regarding foreground residuals compared to 2015), and they are stable with respect to our 2015 measurements (with small fluctuations, at the level of a fraction of a standard deviation, which is consistent with changes in data processing). Polarizationonly bispectra display a significant improvement in robustness; they can now be used independently to set primordial NG constraints with a sensitivity comparable to WMAP temperature-based results and they give excellent agreement. In addition to the analysis of the standard local, equilateral, and orthogonal bispectrum shapes, we consider a large number of additional cases, such as scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The nonprimordial lensing bispectrum is, however, detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5σ. Beyond estimates of individual shape amplitudes, we also present model-independent reconstructions and analyses of the Planck CMB bispectrum. Our final constraint on the local primordial trispectrum shape is glocalNL= (-5.8 ± 6.5) × 104(68% CL, statistical), while constraints for other trispectrum shapes are also determined. Exploiting the tight limits on various bispectrum and trispectrum shapes, we constrain the parameter space of different early-Universe scenarios that generate primordial NG, including general single-field models of inflation, multi-field models (e.g. curvaton models), models of inflation with axion fields producing parity-violation bispectra in the tensor sector, and inflationary models involving vector-like fields with directionally-dependent bispectra. Our results provide a high-precision test for structure-formation scenarios, showing complete agreement with the basic picture of the CDM cosmology regarding the statistics of the initial conditions, with cosmic structures arising from adiabatic, passive, Gaussian, and primordial seed perturbations.

We describe the legacy Planck cosmic microwave background (CMB) likelihoods derived from the 2018 data release. The overall approach is similar in spirit to the one retained for the 2013 and 2015 data release, with a hybrid method using different approximations at low (ℓ < 30) and high (ℓ ≥ 30) multipoles, implementing several methodological and data-analysis refinements compared to previous releases. With more realistic simulations, and better correction and modelling of systematic effects, we can now make full use of the CMB polarization observed in the High Frequency Instrument (HFI) channels. The low-multipole EE cross-spectra from the 100 GHz and 143 GHz data give a constraint on the λCDM reionization optical-depth parameter τ to better than 15% (in combination with the TT low-ℓ data and the high-ℓ temperature and polarization data), tightening constraints on all parameters with posterior distributions correlated with τ. We also update the weaker constraint on τ from the joint TEB likelihood using the Low Frequency Instrument (LFI) channels, which was used in 2015 as part of our baseline analysis. At higher multipoles, the CMB temperature spectrum and likelihood are very similar to previous releases. A better model of the temperature-to-polarization leakage and corrections for the effective calibrations of the polarization channels (i.e., the polarization efficiencies) allow us to make full use of polarization spectra, improving the λCDM constraints on the parameters θMC, ωc, ωb, and H0 by more than 30%, and ns by more than 20% compared to TT-only constraints. Extensive tests on the robustness of the modelling of the polarization data demonstrate good consistency, with some residual modelling uncertainties. At high multipoles, we are now limited mainly by the accuracy of the polarization efficiency modelling. Using our various tests, simulations, and comparison between different high-multipole likelihood implementations, we estimate the consistency of the results to be better than the 0.5σ level on the λCDM parameters, as well as classical single-parameter extensions for the joint likelihood (to be compared to the 0.3σ levels we achieved in 2015 for the temperature data alone on λCDM only). Minor curiosities already present in the previous releases remain, such as the differences between the best-fit λCDM parameters for the ℓ < 800 and ℓ > 800 ranges of the power spectrum, or the preference for more smoothing of the power-spectrum peaks than predicted in λCDM fits. These are shown to be driven by the temperature power spectrum and are not significantly modified by the inclusion of the polarization data. Overall, the legacy Planck CMB likelihoods provide a robust tool for constraining the cosmological model and represent a reference for future CMB observations.

This paper presents the High Frequency Instrument (HFI) data processing procedures for the Planck 2018 release. Major improvements in mapmaking have been achieved since the previous Planck 2015 release, many of which were used and described already in an intermediate paper dedicated to the Planck polarized data at low multipoles. These improvements enabled the first significant measurement of the reionization optical depth parameter using Planck-HFI data. This paper presents an extensive analysis of systematic effects, including the use of end-to-end simulations to facilitate their removal and characterize the residuals. The polarized data, which presented a number of known problems in the 2015 Planck release, are very significantly improved, especially the leakage from intensity to polarization. Calibration, based on the cosmic microwave background (CMB) dipole, is now extremely accurate and in the frequency range 100-353 GHz reduces intensity-to-polarization leakage caused by calibration mismatch. The Solar dipole direction has been determined in the three lowest HFI frequency channels to within one arc minute, and its amplitude has an absolute uncertainty smaller than 0.35 μK, an accuracy of order 10-4. This is a major legacy from the Planck HFI for future CMB experiments. The removal of bandpass leakage has been improved for the main high-frequency foregrounds by extracting the bandpass-mismatch coefficients for each detector as part of the mapmaking process; these values in turn improve the intensity maps. This is a major change in the philosophy of "frequency maps", which are now computed from single detector data, all adjusted to the same average bandpass response for the main foregrounds. End-to-end simulations have been shown to reproduce very well the relative gain calibration of detectors, as well as drifts within a frequency induced by the residuals of the main systematic effect (analogue-to-digital convertor non-linearity residuals). Using these simulations, we have been able to measure and correct the small frequency calibration bias induced by this systematic effect at the 10-4 level. There is no detectable sign of a residual calibration bias between the first and second acoustic peaks in the CMB channels, at the 10-3 level.

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857 GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter ΛCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (θ∗) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the ΛCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

We report on the implications for cosmic inflation of the 2018 release of the Planck cosmic microwave background (CMB) anisotropy measurements. The results are fully consistent with those reported using the data from the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be ns = 0.9649 ± 0.0042 at 68% CL. We find no evidence for a scale dependence of ns, either as a running or as a running of the running. The Universe is found to be consistent with spatial flatness with a precision of 0.4% at 95% CL by combining Planck with a compilation of baryon acoustic oscillation data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, r0.002 < 0.10, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain r0.002 < 0.056. In the framework of standard single-field inflationary models with Einstein gravity, these results imply that: (a) the predictions of slow-roll models with a concave potential, V″(φ) < 0, are increasingly favoured by the data; and (b) based on two different methods for reconstructing the inflaton potential, we find no evidence for dynamics beyond slow roll. Three different methods for the non-parametric reconstruction of the primordial power spectrum consistently confirm a pure power law in the range of comoving scales 0.005 Mpc-1 k 0.2 Mpc-1. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectra. For the case of oscillatory features that are logarithmic or linear in k, this result is further strengthened by a new combined analysis including the Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions for the cosmological fluctuations. In correlated, mixed adiabatic and isocurvature models, the non-adiabatic contribution to the observed CMB temperature variance is constrained to 1.3%, 1.7%, and 1.7% at 95% CL for cold dark matter, neutrino density, and neutrino velocity, respectively. Planck power spectra plus lensing set constraints on the amplitude of compensated cold dark matter-baryon isocurvature perturbations that are consistent with current complementary measurements. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadupolar modulation of the primordial fluctuations. However, the polarization data do not support physical models for a scale-dependent dipolar modulation. All these findings support the key predictions of the standard single-field inflationary models, which will be further tested by future cosmological observations.

We present cosmological parameter results from the final full-mission Planck measurements of the cosmic microwave background (CMB) anisotropies, combining information from the temperature and polarization maps and the lensing reconstruction. Compared to the 2015 results, improved measurements of large-scale polarization allow the reionization optical depth to be measured with higher precision, leading to significant gains in the precision of other correlated parameters. Improved modelling of the small-scale polarization leads to more robust constraints on many parameters, with residual modelling uncertainties estimated to affect them only at the 0.5σ level. We find good consistency with the standard spatially-flat 6-parameter ΛCDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base ΛCDM"in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density ωch2 = 0.120 ± 0.001, baryon density ωbh2 = 0.0224 ± 0.0001, scalar spectral index ns = 0.965 ± 0.004, and optical depth τ = 0.054 ± 0.007 (in this abstract we quote 68% confidence regions on measured parameters and 95% on upper limits). The angular acoustic scale is measured to 0.03% precision, with 100θ∗ = 1.0411 ± 0.0003. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-ΛCDM cosmology, the inferred (model-dependent) late-Universe parameters are: Hubble constant H0 = (67.4 ± 0.5) km s-1 Mpc-1; matter density parameter ωm = 0.315 ± 0.007; and matter fluctuation amplitude σ8 = 0.811 ± 0.006. We find no compelling evidence for extensions to the base-ΛCDM model. Combining with baryon acoustic oscillation (BAO) measurements (and considering single-parameter extensions) we constrain the effective extra relativistic degrees of freedom to be Neff = 2.99 ± 0.17, in agreement with the Standard Model prediction Neff = 3.046, and find that the neutrino mass is tightly constrained to mν < 0.12 eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base ΛCDM at over 2σ, which pulls some parameters that affect the lensing amplitude away from the ΛCDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. The joint constraint with BAO measurements on spatial curvature is consistent with a flat universe, ωK = 0.001 ± 0.002. Also combining with Type Ia supernovae (SNe), the dark-energy equation of state parameter is measured to be w0 = -1.03 ± 0.03, consistent with a cosmological constant. We find no evidence for deviations from a purely power-law primordial spectrum, and combining with data from BAO, BICEP2, and Keck Array data, we place a limit on the tensor-to-scalar ratio r0.002 < 0.06. Standard big-bang nucleosynthesis predictions for the helium and deuterium abundances for the base-ΛCDM cosmology are in excellent agreement with observations. The Planck base-ΛCDM results are in good agreement with BAO, SNe, and some galaxy lensing observations, but in slight tension with the Dark Energy Survey's combined-probe results including galaxy clustering (which prefers lower fluctuation amplitudes or matter density parameters), and in significant, 3.6σ, tension with local measurements of the Hubble constant (which prefer a higher value). Simple model extensions that can partially resolve these tensions are not favoured by the Planck data.

We present measurements of the cosmic microwave background (CMB) lensing potential using the final Planck 2018 temperature and polarization data. Using polarization maps filtered to account for the noise anisotropy, we increase the significance of the detection of lensing in the polarization maps from 5σ to 9σ. Combined with temperature, lensing is detected at 40σ. We present an extensive set of tests of the robustness of the lensingpotential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles 8 ≤ L ≤ 400 (extending the range to lower L compared to 2015), which we use to constrain cosmological parameters. We find good consistency between lensing constraints and the results from the Planck CMB power spectra within the σCDM model. Combined with baryon density and other weak priors, the lensing analysis alone constrains σ8Ω0.25m= 0.589 ± 0.020 (1σ errors). Also combining with baryon acoustic oscillation data, we find tight individual parameter constraints, σ8 = 0.811 ± 0.019, H0 = 67.9+1.2-1.3km s-1Mpc-1, and m = 0.303+0.016-0.018. Combining with Planck CMB power spectrum data, we measure σ8 to better than 1% precision, finding σ8 = 0.811 ± 0.006. CMB lensing reconstruction data are complementary to galaxy lensing data at lower redshift, having a different degeneracy direction in σ8 - m space; we find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using the Planck cosmic infrared background (CIB) maps as an additional tracer of high-redshift matter, we make a combined Planck-only estimate of the lensing potential over 60% of the sky with considerably more small-scale signal.We additionally demonstrate delensing of the Planck power spectra using the joint and individual lensing potential estimates, detecting a maximum removal of 40% of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance.