Melanopsin retinal ganglion cells mediate light-promoted brain development.

During development, melanopsin-expressing intrinsically photosensitive retinal ganglion cells (ipRGCs) become light sensitive much earlier than rods and cones. IpRGCs project to many subcortical areas, whereas physiological functions of these projections are yet to be fully elucidated. Here, we found that ipRGC-mediated light sensation promotes synaptogenesis of pyramidal neurons in various cortices and the hippocampus. This phenomenon depends on activation of ipRGCs and is mediated by the release of oxytocin from the supraoptic nucleus (SON) and the paraventricular nucleus (PVN) into cerebral-spinal fluid. We further characterized a direct connection between ipRGCs and oxytocin neurons in the SON and mutual projections between oxytocin neurons in the SON and PVN. Moreover, we showed that the lack of ipRGC-mediated, light-promoted early cortical synaptogenesis compromised learning ability in adult mice. Our results highlight the importance of light sensation early in life on the development of learning ability and therefore call attention to suitable light environment for infant care.

Liver-Resident NK Cells Control Antiviral Activity of Hepatic T Cells via the PD-1-PD-L1 Axis.

The tolerogenic microenvironment of the liver is associated with impaired hepatic T cell function. Here, we examined the contribution of liver-resident natural killer (LrNK) cells, a prominent hepatic NK cell compartment, to T cell antiviral responses in the liver. The number of virus-specific T cells increased in LrNK-cell-deficient mice during both acute and chronic lymphocytic choriomeningitis virus infection. Upon infection with adenovirus, hepatic T cells from these mice produced more cytokines, which was accompanied by reduced viral loads. Transfer of LrNK cells into LrNK-cell-deficient or wild-type mice inhibited hepatic T cell function, resulting in impaired viral clearance, whereas transfer of conventional NK cells promoted T cell antiviral responses. LrNK-cell-mediated inhibition of T cell function was dependent on the PD-1-PD-L1 axis. Our findings reveal a role for LrNK cells in the regulation of T cell immunity and provide insight into the mechanisms of immune tolerance in the liver.

Hierarchical mechanisms for direct reprogramming of fibroblasts to neurons.

Direct lineage reprogramming is a promising approach for human disease modeling and regenerative medicine, with poorly understood mechanisms. Here, we reveal a hierarchical mechanism in the direct conversion of fibroblasts into induced neuronal (iN) cells mediated by the transcription factors Ascl1, Brn2, and Myt1l. Ascl1 acts as an "on-target" pioneer factor by immediately occupying most cognate genomic sites in fibroblasts. In contrast, Brn2 and Myt1l do not access fibroblast chromatin productively on their own; instead, Ascl1 recruits Brn2 to Ascl1 sites genome wide. A unique trivalent chromatin signature in the host cells predicts the permissiveness for Ascl1 pioneering activity among different cell types. Finally, we identified Zfp238 as a key Ascl1 target gene that can partially substitute for Ascl1 during iN cell reprogramming. Thus, a precise match between pioneer factors and the chromatin context at key target genes is determinative for transdifferentiation to neurons and likely other cell types.

eccDNA-pipe: an integrated pipeline for identification, analysis and visualization of extrachromosomal circular DNA from high-throughput sequencing data.

Extrachromosomal circular DNA (eccDNA) is currently attracting considerable attention from researchers due to its significant impact on tumor biogenesis. High-throughput sequencing (HTS) methods for eccDNA identification are continually evolving. However, an efficient pipeline for the integrative and comprehensive analysis of eccDNA obtained from HTS data is still lacking. Here, we introduce eccDNA-pipe, an accessible software package that offers a user-friendly pipeline for conducting eccDNA analysis starting from raw sequencing data. This dataset includes data from various sequencing techniques such as whole-genome sequencing (WGS), Circle-seq and Circulome-seq, obtained through short-read sequencing or long-read sequencing. eccDNA-pipe presents a comprehensive solution for both upstream and downstream analysis, encompassing quality control and eccDNA identification in upstream analysis and downstream tasks such as eccDNA length distribution analysis, differential analysis of genes enriched with eccDNA and visualization of eccDNA structures. Notably, eccDNA-pipe automatically generates high-quality publication-ready plots. In summary, eccDNA-pipe provides a comprehensive and user-friendly pipeline for customized analysis of eccDNA research.

The DNA Methylcytosine Dioxygenase Tet2 Sustains Immunosuppressive Function of Tumor-Infiltrating Myeloid Cells to Promote Melanoma Progression.

Ten-Eleven-Translocation-2 (Tet2) is a DNA methylcytosine dioxygenase that functions as a tumor suppressor in hematopoietic malignancies. We examined the role of Tet2 in tumor-tissue myeloid cells and found that Tet2 sustains the immunosuppressive function of these cells. We found that Tet2 expression is increased in intratumoral myeloid cells both in mouse models of melanoma and in melanoma patients and that this increased expression is dependent on an IL-1R-MyD88 pathway. Ablation of Tet2 in myeloid cells suppressed melanoma growth in vivo and shifted the immunosuppressive gene expression program in tumor-associated macrophages to a proinflammatory one, with a concomitant reduction of the immunosuppressive function. This resulted in increased numbers of effector T cells in the tumor, and T cell depletion abolished the reduced tumor growth observed upon myeloid-specific deletion of Tet2. Our findings reveal a non-cell-intrinsic, tumor-promoting function for Tet2 and suggest that Tet2 may present a therapeutic target for the treatment of non-hematologic malignancies.

The native structure of the assembled matrix protein 1 of influenza A virus.

Influenza A virus causes millions of severe cases of disease during annual epidemics. The most abundant protein in influenza virions is matrix protein 1 (M1), which mediates virus assembly by forming an endoskeleton beneath the virus membrane1. The structure of full-length M1, and how it oligomerizes to mediate the assembly of virions, is unknown. Here we determine the complete structure of assembled M1 within intact virus particles, as well as the structure of M1 oligomers reconstituted in vitro. We find that the C-terminal domain of M1 is disordered in solution but can fold and bind in trans to the N-terminal domain of another M1 monomer, thus polymerizing M1 into linear strands that coat the interior surface of the membrane of the assembling virion. In the M1 polymer, five histidine residues-contributed by three different monomers of M1-form a cluster that can serve as the pH-sensitive disassembly switch after entry into a target cell. These structures therefore reveal mechanisms of influenza virus assembly and disassembly.

Structures and distributions of SARS-CoV-2 spike proteins on intact virions.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virions are surrounded by a lipid bilayer from which spike (S) protein trimers protrude1. Heavily glycosylated S trimers bind to the angiotensin-converting enzyme 2 receptor and mediate entry of virions into target cells2-6. S exhibits extensive conformational flexibility: it modulates exposure of its receptor-binding site and subsequently undergoes complete structural rearrangement to drive fusion of viral and cellular membranes2,7,8. The structures and conformations of soluble, overexpressed, purified S proteins have been studied in detail using cryo-electron microscopy2,7,9-12, but the structure and distribution of S on the virion surface remain unknown. Here we applied cryo-electron microscopy and tomography to image intact SARS-CoV-2 virions and determine the high-resolution structure, conformational flexibility and distribution of S trimers in situ on the virion surface. These results reveal the conformations of S on the virion, and provide a basis from which to understand interactions between S and neutralizing antibodies during infection or vaccination.

Cryo-EM Structures of MDA5-dsRNA Filaments at Different Stages of ATP Hydrolysis.

Double-stranded RNA (dsRNA) is a potent proinflammatory signature of viral infection. Long cytosolic dsRNA is recognized by MDA5. The cooperative assembly of MDA5 into helical filaments on dsRNA nucleates the assembly of a multiprotein type I interferon signaling platform. Here, we determined cryoelectron microscopy (cryo-EM) structures of MDA5-dsRNA filaments with different helical twists and bound nucleotide analogs at resolutions sufficient to build and refine atomic models. The structures identify the filament-forming interfaces, which encode the dsRNA binding cooperativity and length specificity of MDA5. The predominantly hydrophobic interface contacts confer flexibility, reflected in the variable helical twist within filaments. Mutation of filament-forming residues can result in loss or gain of signaling activity. Each MDA5 molecule spans 14 or 15 RNA base pairs, depending on the twist. Variations in twist also correlate with variations in the occupancy and type of nucleotide in the active site, providing insights on how ATP hydrolysis contributes to MDA5-dsRNA recognition.

Benchmarking spatial and single-cell transcriptomics integration methods for transcript distribution prediction and cell type deconvolution.

Spatial transcriptomics approaches have substantially advanced our capacity to detect the spatial distribution of RNA transcripts in tissues, yet it remains challenging to characterize whole-transcriptome-level data for single cells in space. Addressing this need, researchers have developed integration methods to combine spatial transcriptomic data with single-cell RNA-seq data to predict the spatial distribution of undetected transcripts and/or perform cell type deconvolution of spots in histological sections. However, to date, no independent studies have comparatively analyzed these integration methods to benchmark their performance. Here we present benchmarking of 16 integration methods using 45 paired datasets (comprising both spatial transcriptomics and scRNA-seq data) and 32 simulated datasets. We found that Tangram, gimVI, and SpaGE outperformed other integration methods for predicting the spatial distribution of RNA transcripts, whereas Cell2location, SpatialDWLS, and RCTD are the top-performing methods for the cell type deconvolution of spots. We provide a benchmark pipeline to help researchers select optimal integration methods to process their datasets.

Molecular dynamics simulations of HIV-1 matrix-membrane interactions at different stages of viral maturation.

Although the structural rearrangement of the membrane-bound matrix (MA) protein trimers upon HIV-1 maturation has been reported, the consequences of MA maturation on the MA-lipid interactions are not well understood. Long-timescale molecular dynamics simulations of the MA multimeric assemblies of immature and mature virus particles with our realistic asymmetric membrane model have explored MA-lipid interactions and lateral organization of lipids around MA complexes. The number of stable MA-phosphatidylserine and MA-phosphatidylinositol 4,5-bisphosphate (PIP2) interactions at the trimeric interface of the mature MA complex is observed to be greater compared to that of the immature MA complex. Our simulations identified an alternative PIP2-binding site in the immature MA complex where the multivalent headgroup of a PIP2 lipid with a greater negative charge binds to multiple basic amino acid residues such as ARG3 residues of both the MA monomers at the trimeric interface and highly basic region (HBR) residues (LYS29, LYS31) of one of the MA monomers. Our enhanced sampling simulations have explored the conformational space of phospholipids at different binding sites of the trimer-trimer interface of MA complexes that are not accessible by conventional unbiased molecular dynamics. Unlike the immature MA complex, the 2' acyl tail of two PIP2 lipids at the trimeric interface of the mature MA complex is observed to sample stable binding pockets of MA consisting of helix-4 residues. Together, our results provide molecular-level insights into the interactions of MA trimeric complexes with membrane and different lipid conformations at the specific binding sites of MA protein before and after viral maturation.

Engineered disulfide reveals structural dynamics of locked SARS-CoV-2 spike.

The spike (S) protein of SARS-CoV-2 has been observed in three distinct pre-fusion conformations: locked, closed and open. Of these, the function of the locked conformation remains poorly understood. Here we engineered a SARS-CoV-2 S protein construct "S-R/x3" to arrest SARS-CoV-2 spikes in the locked conformation by a disulfide bond. Using this construct we determined high-resolution structures confirming that the x3 disulfide bond has the ability to stabilize the otherwise transient locked conformations. Structural analyses reveal that wild-type SARS-CoV-2 spike can adopt two distinct locked-1 and locked-2 conformations. For the D614G spike, based on which all variants of concern were evolved, only the locked-2 conformation was observed. Analysis of the structures suggests that rigidified domain D in the locked conformations interacts with the hinge to domain C and thereby restrains RBD movement. Structural change in domain D correlates with spike conformational change. We propose that the locked-1 and locked-2 conformations of S are present in the acidic high-lipid cellular compartments during virus assembly and egress. In this model, release of the virion into the neutral pH extracellular space would favour transition to the closed or open conformations. The dynamics of this transition can be altered by mutations that modulate domain D structure, as is the case for the D614G mutation, leading to changes in viral fitness. The S-R/x3 construct provides a tool for the further structural and functional characterization of the locked conformations of S, as well as how sequence changes might alter S assembly and regulation of receptor binding domain dynamics.

A genome-wide long noncoding RNA CRISPRi screen identifies PRANCR as a novel regulator of epidermal homeostasis.

Genome-wide association studies indicate that many disease susceptibility regions reside in non-protein-coding regions of the genome. Long noncoding RNAs (lncRNAs) are a major component of the noncoding genome, but their biological impacts are not fully understood. Here, we performed a CRISPR interference (CRISPRi) screen on 2263 epidermis-expressed lncRNAs and identified nine novel candidate lncRNAs regulating keratinocyte proliferation. We further characterized a top hit from the screen, progenitor renewal associated non-coding RNA (PRANCR), using RNA interference-mediated knockdown and phenotypic analysis in organotypic human tissue. PRANCR regulates keratinocyte proliferation, cell cycle progression, and clonogenicity. PRANCR-deficient epidermis displayed impaired stratification with reduced expression of differentiation genes that are altered in human skin diseases, including keratins 1 and 10, filaggrin, and loricrin. Transcriptome analysis showed that PRANCR controls the expression of 1136 genes, with strong enrichment for late cell cycle genes containing a CHR promoter element. In addition, PRANCR depletion led to increased levels of both total and nuclear CDKN1A (also known as p21), which is known to govern both keratinocyte proliferation and differentiation. Collectively, these data show that PRANCR is a novel lncRNA regulating epidermal homeostasis and identify other lncRNA candidates that may have roles in this process as well.

WEDGE: imputation of gene expression values from single-cell RNA-seq datasets using biased matrix decomposition.

The low capture rate of expressed RNAs from single-cell sequencing technology is one of the major obstacles to downstream functional genomics analyses. Recently, a number of imputation methods have emerged for single-cell transcriptome data, however, recovering missing values in very sparse expression matrices remains a substantial challenge. Here, we propose a new algorithm, WEDGE (WEighted Decomposition of Gene Expression), to impute gene expression matrices by using a biased low-rank matrix decomposition method. WEDGE successfully recovered expression matrices, reproduced the cell-wise and gene-wise correlations and improved the clustering of cells, performing impressively for applications with sparse datasets. Overall, this study shows a potent approach for imputing sparse expression matrix data, and our WEDGE algorithm should help many researchers to more profitably explore the biological meanings embedded in their single-cell RNA sequencing datasets. The source code of WEDGE has been released at https://github.com/QuKunLab/WEDGE.

CALML5 is a ZNF750- and TINCR-induced protein that binds stratifin to regulate epidermal differentiation.

Outward migration of epidermal progenitors occurs with induction of hundreds of differentiation genes, but the identities of all regulators required for this process are unknown. We used laser capture microdissection followed by RNA sequencing to identify calmodulin-like 5 (CALML5) as the most enriched gene in differentiating outer epidermis. CALML5 mRNA was up-regulated by the ZNF750 transcription factor and then stabilized by the long noncoding RNA TINCR. CALML5 knockout impaired differentiation, abolished keratohyalin granules, and disrupted epidermal barrier function. Mass spectrometry identified SFN (stratifin/14-3-3σ) as a CALML5-binding protein. CALML5 interacts with SFN in suprabasal epidermis, cocontrols 13% of late differentiation genes, and modulates interaction of SFN to some of its binding partners. A ZNF750-TINCR-CALML5-SFN network is thus essential for epidermal differentiation.

Impacts of synoptic circulation on surface ozone pollution in a coastal eco-city in Southeastern China during 2014-2019.

The coastal eco-city of Fuzhou in Southeastern China has experienced severe ozone (O3) episodes at times in recent years. In this study, three typical synoptic circulations types (CTs) that influenced more than 80% of O3 polluted days in Fuzhou during 2014-2019 were identified using a subjective approach. The characteristics of meteorological conditions linked to photochemical formation and transport of O3 under the three CTs were summarized. Comprehensive Air Quality Model with extensions was applied to simulate O3 episodes and to quantify O3 sources from different regions in Fuzhou. When Fuzhou was located to the west of a high-pressure system (classified as "East-ridge"), more warm southwesterly currents flowed to Fuzhou, and the effects of cross-regional transport from Guangdong province and high local production promoted the occurrence of O3 episodes. Under a uniform pressure field with a low-pressure system occurring to the east of Fuzhou (defined as "East-low"), stagnant weather conditions caused the strongest local production of O3 in the atmospheric boundary layer. Controlled by high-pressure systems over the mainland (categorized as "Inland-high"), northerly airflows enhanced the contribution of cross-regional transport to O3 in Fuzhou. The abnormal increases of the "East-ridge" and "Inland-high" were closely related to O3 pollution in Fuzhou in April and May 2018, resulting in the annual maximum number of O3 polluted days during recent years. Furthermore, the rising number of autumn O3 episodes in 2017-2019 was mainly related to the "Inland-high", indicating the aggravation of cross-regional transport and highlighting the necessity of enhanced regional collaboration and efforts in combating O3 pollution.

Single-Cell Heterogeneity of the Liver-Infiltrating Lymphocytes in Individuals with Chronic Echinococcus multilocularis Infection.

Human alveolar echinococcosis (AE) is a tumor-like disease predominantly located in the liver. The cellular composition and heterogeneity of the lesion-infiltrating lymphocytes which produce an "immunosuppressive" microenvironment are poorly understood. Here, we profiled 83,921 CD45+ lymphocytes isolated from the peripheral blood (PB), perilesion (PL), and adjacent normal (AN) liver tissue of four advanced-stage AE patients using single-cell RNA and T-cell receptor (TCR) sequencing technology. We identified 23 large clusters, and the distributions and transcriptomes of these cell clusters in the liver and periphery were different. The cellular proportions of exhausted CD8+ T cells and group 2 innate lymphoid cells (ILC2s) were notably higher in PL tissue, and the expression features of these cell subsets were related to neoplasm metastasis and immune response suppression. In the 5 CD8+ T-cell populations, only CD8+ mucosa-associated invariant T (MAIT) cells were enriched in PL samples and the TRAV1-2_TRAJ33_TRAC TCR was clonally expanded. In the 11 subsets of CD4+ T cells, Th17 cells and induced regulatory T cells (iTregs) were preferentially enriched in PL samples, and their highly expressed genes were related to cell invasion, tumor metastasis, and inhibition of the inflammatory immune response. Exhaustion-specific genes (TIGIT, PD-1, and CTLA4) were upregulated in Tregs. Interestingly, there was a close contact between CD8+ T cells and iTregs or Th17 cells, especially for genes related to immunosuppression, such as PDCD1-FAM3C, which were highly expressed in PL tissue. This transcriptional data set provides valuable insights and a rich resource for deeply understanding the immune microenvironment in AE, which could provide potential target signatures for AE diagnosis and immunotherapies.

Chromatin accessibility landscapes of immune cells in rheumatoid arthritis nominate monocytes in disease pathogenesis.

BACKGROUND: Rheumatoid arthritis (RA) is a chronic, systemic autoimmune disease that involves a variety of cell types. However, how the epigenetic dysregulations of peripheral immune cells contribute to the pathogenesis of RA still remains largely unclear. RESULTS: Here, we analysed the genome-wide active DNA regulatory elements of four major immune cells, namely monocytes, B cells, CD4+ T cells and CD8+ T cells, in peripheral blood of RA patients, osteoarthritis (OA) patients and healthy donors using Assay of Transposase Accessible Chromatin with sequencing (ATAC-seq). We found a strong RA-associated chromatin dysregulation signature in monocytes, but no other examined cell types. Moreover, we found that serum C-reactive protein (CRP) can induce the RA-associated chromatin dysregulation in monocytes via in vitro experiments. And the extent of this dysregulation was regulated through the transcription factor FRA2. CONCLUSIONS: Together, our study revealed a CRP-induced pathogenic chromatin dysregulation signature in monocytes from RA patients and predicted the responsible signalling pathway as potential therapeutic targets for the disease.

Single-Cell RNA Sequencing Reveals the Heterogeneity of Infiltrating Immune Cell Profiles in the Hepatic Cystic Echinococcosis Microenvironment.

Human cystic echinococcosis, caused by the larval stage of Echinococcus granulosus sensu lato, has been reported a near-cosmopolitan zoonotic disease. Various infiltrating immune cells gather around the lesion and produce a lesion microenvironment; however, cellular composition and heterogeneity in hepatic cystic echinococcosis lesion microenvironments are incompletely understood. Here, 81,865 immune cells isolated from peripheral blood, perilesion liver tissue, and adjacent normal liver tissue from four cystic echinococcosis patients were profiled using single-cell RNA sequencing. We identified 23 discrete cell populations and found distinct differences in infiltrating immune cells between tissue environments. Despite the significant similarity between perilesion and adjacent normal liver tissue-resident immune cells, the cellular proportions of type 2 innate lymphoid cells (ILC2s) and plasmacytoid dendritic cells (pDCs) were higher in perilesion liver tissue. Interestingly, the immunosuppressive gene NFKBIA was upregulated in these cells. Seven subsets of CD4+ T cell populations were found, and there were more regulatory-CD4+ T cells (Treg-CD4+) and Th2-CD4+ T cells in perilesion tissue than in adjacent normal tissue. There was close contact between CD4+ T cells and ILC2s and pDCs, which caused upregulation of genes related to positive immune activity in adjacent normal liver tissue. However, expression of genes related to immunosuppression, especially the immune inhibitory checkpoint gene NKG2A/HLA-E, was obviously higher in perilesion tissue, suggesting that cellular interaction resulted in an inhibitory microenvironment in the cystic echinococcosis (CE) lesion. This work offers new insights into the transcriptional heterogeneity of infiltrating immune cells in hepatic cystic echinococcosis lesion microenvironments at a single-cell level and provides potential target signatures for diagnosis and immunotherapies.

ATAC-pipe: general analysis of genome-wide chromatin accessibility.

Assay of Transposase-Accessible Chromatin by deep sequencing (ATAC-seq) has been widely used to profile the chromatin accessibility genome-wide. For the absence of an integrated scheme for deep data mining of specific biological issues, here we present ATAC-pipe, an efficient pipeline for general analysis of chromatin accessibility data obtained from ATAC-seq experiments. ATAC-pipe captures information includes not only the quality of original data and genome-wide chromatin accessibility but also signatures of significant differential peaks, transcription factor (TF) occupancy and nucleosome positions around regulatory sites. In addition, ATAC-pipe automatically converts statistic results into intuitive plots at publication quality, such as the read length distribution, heatmaps of sample clustering and cell-type-specific regulatory elements, enriched TF occupancy with motifs footprints and TF-driven regulatory networks. ATAC-pipe provides convenient workflow for researchers to study chromatin accessibility and gene regulation. Availability https://github.com/QuKunLab/ATAC-pipe.

Structure and architecture of immature and mature murine leukemia virus capsids.

Retroviruses assemble and bud from infected cells in an immature form and require proteolytic maturation for infectivity. The CA (capsid) domains of the Gag polyproteins assemble a protein lattice as a truncated sphere in the immature virion. Proteolytic cleavage of Gag induces dramatic structural rearrangements; a subset of cleaved CA subsequently assembles into the mature core, whose architecture varies among retroviruses. Murine leukemia virus (MLV) is the prototypical γ-retrovirus and serves as the basis of retroviral vectors, but the structure of the MLV CA layer is unknown. Here we have combined X-ray crystallography with cryoelectron tomography to determine the structures of immature and mature MLV CA layers within authentic viral particles. This reveals the structural changes associated with maturation, and, by comparison with HIV-1, uncovers conserved and variable features. In contrast to HIV-1, most MLV CA is used for assembly of the mature core, which adopts variable, multilayered morphologies and does not form a closed structure. Unlike in HIV-1, there is similarity between protein-protein interfaces in the immature MLV CA layer and those in the mature CA layer, and structural maturation of MLV could be achieved through domain rotations that largely maintain hexameric interactions. Nevertheless, the dramatic architectural change on maturation indicates that extensive disassembly and reassembly are required for mature core growth. The core morphology suggests that wrapping of the genome in CA sheets may be sufficient to protect the MLV ribonucleoprotein during cell entry.