Construction of single-cell cross-species chromatin accessibility landscapes with combinatorial-hybridization-based ATAC-seq.

Despite recent advances in single-cell genomics, the lack of maps for single-cell candidate cis-regulatory elements (cCREs) in non-mammal species has limited our exploration of conserved regulatory programs across vertebrates and invertebrates. Here, we developed a combinatorial-hybridization-based method for single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) named CH-ATAC-seq, enabling the construction of single-cell accessible chromatin landscapes for zebrafish, Drosophila, and earthworms (Eisenia andrei). By integrating scATAC censuses of humans, monkeys, and mice, we systematically identified 152 distinct main cell types and around 0.8 million cell-type-specific cCREs. Our analysis provided insights into the conservation of neural, muscle, and immune lineages across species, while epithelial cells exhibited a higher organ-origin heterogeneity. Additionally, a large-scale gene regulatory network (GRN) was constructed in four vertebrates by integrating scRNA-seq censuses. Overall, our study provides a valuable resource for comparative epigenomics, identifying the evolutionary conservation and divergence of gene regulation across different species.

DUB3 is a MAGEA3 deubiquitinase and a potential therapeutic target in hepatocellular carcinoma.

Although melanoma-associated antigen A3 and A6 (MAGEA3/6)-specific tumor vaccines have shown antitumor effects in melanoma and non-small cell lung cancer (NSCLC), many cancers do not respond because MAGEA3 can promote cancer without triggering an immune response. Here, we identified DUB3 as the MAGEA3 deubiquitinase. DUB3 interacts with, deubiquitinates and stabilizes MAGEA3. Depletion of DUB3 in hepatocellular carcinoma (HCC) cells results in MAGEA3 degradation and P53-dependent growth inhibition. Moreover, DUB3 knockout attenuates HCC tumorigenesis in vivo, which can be rescued by restoration of MAGEA3. Intriguingly, pharmacological inhibition of DUB3 by palbociclib promotes degradation of MAGEA3 and inhibits tumor growth in preclinical models implanted with parental HCC cells but not with DUB3 knockout HCC cells. In patients with HCC, DUB3 is highly expressed, and its levels positively correlate with MAGEA3 levels. Taken together, DUB3 is a MAGEA3 deubiquitinase, and abrogating DUB3 enzymatic activity by palbociclib is a promising therapeutic strategy for HCC.

Pan-Cancer Single-Nucleus Total RNA Sequencing Using snHH-Seq.

Tumor heterogeneity and its drivers impair tumor progression and cancer therapy. Single-cell RNA sequencing is used to investigate the heterogeneity of tumor ecosystems. However, most methods of scRNA-seq amplify the termini of polyadenylated transcripts, making it challenging to perform total RNA analysis and somatic mutation analysis.Therefore, a high-throughput and high-sensitivity method called snHH-seq is developed, which combines random primers and a preindex strategy in the droplet microfluidic platform. This innovative method allows for the detection of total RNA in single nuclei from clinically frozen samples. A robust pipeline to facilitate the analysis of full-length RNA-seq data is also established. snHH-seq is applied to more than 730 000 single nuclei from 32 patients with various tumor types. The pan-cancer study enables it to comprehensively profile data on the tumor transcriptome, including expression levels, mutations, splicing patterns, clone dynamics, etc. New malignant cell subclusters and exploring their specific function across cancers are identified. Furthermore, the malignant status of epithelial cells is investigated among different cancer types with respect to mutation and splicing patterns. The ability to detect full-length RNA at the single-nucleus level provides a powerful tool for studying complex biological systems and has broad implications for understanding tumor pathology.

Characterization of human pluripotent stem cell differentiation by single-cell dual-omics analyses.

In vivo differentiation of human pluripotent stem cells (hPSCs) has unique advantages, such as multilineage differentiation, angiogenesis, and close cell-cell interactions. To systematically investigate multilineage differentiation mechanisms of hPSCs, we constructed the in vivo hPSC differentiation landscape containing 239,670 cells using teratoma models. We identified 43 cell types, inferred 18 cell differentiation trajectories, and characterized common and specific gene regulation patterns during hPSC differentiation at both transcriptional and epigenetic levels. Additionally, we developed the developmental single-cell Basic Local Alignment Search Tool (dscBLAST), an R-based cell identification tool, to simplify the identification processes of developmental cells. Using dscBLAST, we aligned cells in multiple differentiation models to normally developing cells to further understand their differentiation states. Overall, our study offers new insights into stem cell differentiation and human embryonic development; dscBLAST shows favorable cell identification performance, providing a powerful identification tool for developmental cells.

USP39-Mediated Non-Proteolytic Control of ETS2 Suppresses Nuclear Localization and Activity.

ETS2 is a member of the ETS family of transcription factors and has been implicated in the regulation of cell proliferation, differentiation, apoptosis, and tumorigenesis. The aberrant activation of ETS2 is associated with various human cancers, highlighting its importance as a therapeutic target. Understanding the regulatory mechanisms and interacting partners of ETS2 is crucial for elucidating its precise role in cellular processes and developing novel strategies to modulate its activity. In this study, we conducted binding assays using a human deubiquitinase (DUB) library and identified USP39 as a novel ETS2-binding DUB. USP39 interacts with ETS2 through their respective amino-terminal regions, and the zinc finger and PNT domains are not required for this binding. USP39 deubiquitinates ETS2 without affecting its protein stability. Interestingly, however, USP39 significantly suppresses the transcriptional activity of ETS2. Furthermore, we demonstrated that USP39 leads to a reduction in the nuclear localization of ETS2. Our findings provide valuable insights into the intricate regulatory mechanisms governing ETS2 function. Understanding the interplay between USP39 and ETS2 may have implications for therapeutic interventions targeting ETS2-related diseases, including cancer, where the dysregulation of ETS2 is frequently observed.

The deubiquitinase ZRANB1 is an E3 ubiquitin ligase for SLC7A11 and regulates ferroptotic resistance.

The dependency of cancer cells on iron increases their susceptibility to ferroptosis, thus providing new opportunities for patients with treatment-resistant tumors. However, we show that lipid peroxidation, a hallmark of ferroptosis, was found in various areas of patient samples, indicating the potential resistance of ferroptosis. Using whole deubiquitinases (DUBs) sgRNA screening, we found that loss of ZRANB1 confers cancer cell resistance to ferroptosis. Intriguingly, functional studies revealed that ZRANB1 ubiquitinates and represses SLC7A11 expression as an E3 ubiquitin ligase and that ZRANB1 inhibits glutathione (GSH) synthesis through SLC7A11 degradation, leading to elevated lipid peroxidation and ferroptosis. Deletion of the region (residues 463-584) abolishes the E3 activity of ZRANB1. Moreover, we show that ZRANB1 has lower expression in tumors, which is positively correlated with lipid peroxidation. Collectively, our results demonstrate the role of ZRANB1 in ferroptosis resistance and unveil mechanisms involving modulation of E3 ligase activity through an unconventional catalytic domain.

ZHX2 deficiency enriches hybrid MET cells through regulating E-cadherin expression.

Growing evidence indicates that the epithelial to mesenchymal (E/M) hybrid state plays a key role in tumorigenesis. Importantly, a hybrid mesenchymal to epithelial transition (MET) state in which individual cells express both epithelial and mesenchymal markers was recently identified in vivo, further strengthening the bonds between the hybrid EMT state and cancer progression. However, the role and the molecular mechanisms by which the hybrid MET state is maintained in triple-negative breast cancer cells (TNBC) remain elusive. Here, we find that loss of ZHX2 expression results in the hybrid MET phenotype in mesenchymal TNBC cells. Mechanistically, through directly binding to the CDH1 promoter, depletion of ZHX2 specifically reactivates expression of CDH1 encoding E-cadherin, an epithelial marker that is crucial for maintaining epithelial phenotype. Functionally, loss of ZHX2 expression enriches the hybrid MET cells and inhibits the migration and dissemination of TNBC cells or organoids, which could be reversed by restoration of E-cadherin. Moreover, depletion of ZHX2 suppresses lung metastasis in preclinical models of TNBC. In patients with TNBC, ZHX2 expression was amplified and negatively correlated with the expression of E-cadherin. These findings suggest that loss of ZHX2 promotes the hybrid MET state to impair TNBC progression.

An analytical framework for decoding cell type-specific genetic variation of gene regulation.

A deeper understanding of genetic regulation and functional mechanisms underlying genetic associations with complex traits and diseases is impeded by cellular heterogeneity and linkage disequilibrium. To address these limits, we introduce Huatuo, a framework to decode genetic variation of gene regulation at cell type and single-nucleotide resolutions by integrating deep-learning-based variant predictions with population-based association analyses. We apply Huatuo to generate a comprehensive cell type-specific genetic variation landscape across human tissues and further evaluate their potential roles in complex diseases and traits. Finally, we show that Huatuo's inferences permit prioritizations of driver cell types associated with complex traits and diseases and allow for systematic insights into the mechanisms of phenotype-causal genetic variation.

Network integration and protein structural binding analysis of neurodegeneration-related interactome.

Neurodegenerative diseases (NDs) usually connect with aggregation and molecular interactions of pathological proteins. The integration of accumulative data from clinical and biomedical research will allow for the excavation of pathological proteins and related interactors. It is also important to systematically study their interacting proteins in order to find more related proteins and potential therapeutic targets. Understanding binding regions in protein interactions will help functional proteomics and provide an alternative method for predicting novel interactions. This study integrated data from biomedical research to achieve systematic mining and analysis of pathogenic proteins and their interaction network. A workflow has been built as a solution for the collective information of proteins involved in NDs, related protein-protein interactions (PPIs) and interactive visualizations. It also included protein isoforms and mapped them in a disease-related PPI network to illuminate the impact of alternative splicing on protein binding. The interacting proteins enriched by diseases and biological processes (BPs) revealed possible regulatory modules. A high-resolution network with structural affinity information was generated. Finally, Neurodegenerative Disease Atlas (NDAtlas) was constructed with an interactive and intuitive view of protein docking with 3D molecular graphics beyond the traditional 2D network. NDAtlas is available at http://bis.zju.edu.cn/ndatlas.

Construction of a cross-species cell landscape at single-cell level.

Individual cells are basic units of life. Despite extensive efforts to characterize the cellular heterogeneity of different organisms, cross-species comparisons of landscape dynamics have not been achieved. Here, we applied single-cell RNA sequencing (scRNA-seq) to map organism-level cell landscapes at multiple life stages for mice, zebrafish and Drosophila. By integrating the comprehensive dataset of > 2.6 million single cells, we constructed a cross-species cell landscape and identified signatures and common pathways that changed throughout the life span. We identified structural inflammation and mitochondrial dysfunction as the most common hallmarks of organism aging, and found that pharmacological activation of mitochondrial metabolism alleviated aging phenotypes in mice. The cross-species cell landscape with other published datasets were stored in an integrated online portal-Cell Landscape. Our work provides a valuable resource for studying lineage development, maturation and aging.

How many cell types are there in nature? How do they change during the life cycle? These are two fundamental questions that researchers have been trying to understand in the area of biology. In this study, single-cell mRNA sequencing data were used to profile over 2.6 million individual cells from mice, zebrafish and Drosophila at different life stages, 1.3 million of which were newly collected. The comprehensive datasets allow investigators to construct a cross-species cell landscape that helps to reveal the conservation and diversity of cell taxonomies at genetic and regulatory levels. The resources in this study are assembled into a publicly available website at http://bis.zju.edu.cn/cellatlas/.

Activation of ALOX12 by a multi-organelle-orienting photosensitizer drives ACSL4-independent cell ferroptosis.

Ferroptosis is a recently-defined tumor suppression mechanism, but the sensitivity of many tumorigenic cells to ferroptosis is limited by their deficient expression of acyl-CoA synthetase long-chain family member 4 (ACSL4). Here, we report the discovery of a photosensitizer, namely TPCI, which can evoke ACSL4-independent ferroptosis of cancer cells in photodynamic therapy. Through co-localization with 12-lipoxygenase (ALOX12) in multiple subcellular organelles, TPCI activates ALOX12 to generate lipid reactive oxygen species in large quantity and trigger cell ferroptosis. Intriguingly, confining TPCI exclusively in lysosomes switches the cell death from ferroptosis to apoptosis. More strikingly, the ferroptosis mediated by TPCI-induced ALOX12 activation does not require the participation of ACSL4. Therefore, our study identifies TPCI as the first ALOX12 activator to induce ferroptosis independent of ACSL4, which renders a viable therapeutic approach on the basis of distinct ferroptosis of cancer cells, regardless their ACSL4 expressions.

Deep learning of cross-species single-cell landscapes identifies conserved regulatory programs underlying cell types.

Despite extensive efforts to generate and analyze reference genomes, genetic models to predict gene regulation and cell fate decisions are lacking for most species. Here, we generated whole-body single-cell transcriptomic landscapes of zebrafish, Drosophila and earthworm. We then integrated cell landscapes from eight representative metazoan species to study gene regulation across evolution. Using these uniformly constructed cross-species landscapes, we developed a deep-learning-based strategy, Nvwa, to predict gene expression and identify regulatory sequences at the single-cell level. We systematically compared cell-type-specific transcription factors to reveal conserved genetic regulation in vertebrates and invertebrates. Our work provides a valuable resource and offers a new strategy for studying regulatory grammar in diverse biological systems.

Screening of cell-virus, cell-cell, gene-gene crosstalk among animal kingdom at single cell resolution.

BACKGROUND: The exact animal origin of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains obscure and understanding its host range is vital for preventing interspecies transmission. METHODS: Herein, we applied single-cell sequencing to multiple tissues of 20 species (30 data sets) and integrated them with public resources (45 data sets covering 26 species) to expand the virus receptor distribution investigation. While the binding affinity between virus and receptor is essential for viral infectivity, understanding the receptor distribution could predict the permissive organs and tissues when infection occurs. RESULTS: Based on the transcriptomic data, the expression profiles of receptor or associated entry factors for viruses capable of causing respiratory, blood, and brain diseases were described in detail. Conserved cellular connectomes and regulomes were also identified, revealing fundamental cell-cell and gene-gene cross-talks from reptiles to humans. CONCLUSIONS: Overall, our study provides a resource of the single-cell atlas of the animal kingdom which could help to identify the potential host range and tissue tropism of viruses and reveal the host-virus co-evolution.

[Clinical observation on wrist-ankle acupuncture for shoulder-hand syndrome phaseâ after stroke].

OBJECTIVE: To compare the clinical efficacy between wrist-ankle acupuncture and conventional acupuncture on shoulder-hand syndrome (SHS) phaseⅠafter stroke. METHODS: A total of 64 patients with SHS phaseⅠafter stroke were randomized into a wrist-ankle acupuncture group and a conventional acupuncture group, 32 cases in each group. On the basis treatment of internal medicine and conventional rehabilitation, wrist-ankle acupuncture was applied at upper 4 area, upper 5 area and upper 6 area on the affected side in the wrist-ankle acupuncture group, while acupuncture was applied at Jianyu (LI 15), Quchi (LI 11), Shousanli (LI 10), etc. on the affected side in the conventional acupuncture group. The treatment was given 30 min each time, once a day, 5 days a week for 3 weeks in both groups. Before and after treatment, the visual analogue scale (VAS) score, degree of hand swelling, shoulder-hand syndrome scale (SHSS) score, Fugl-Meyer assessment for upper extremity (FMA-UE) score and modified Barthel index (MBI) score were observed, and the clinical therapeutic effect was evaluated in both groups. RESULTS: After treatment, the VAS scores, degree of hand swelling and SHSS scores were decreased (P<0.05), and the FMA-UE scores and MBI scores were increased (P<0.05) compared before treatment in both groups; in the wrist-ankle acupuncture group, the VAS score, degree of hand swelling and SHSS score were lower (P<0.05), and the FMA-UE score and MBI score were higher (P<0.05) than those in the conventional acupuncture group. The total effective rate was 96.9% (31/32) in the wrist-ankle acupuncture group, which was superior to 90.6% (29/32) in the conventional acupuncture group (P<0.05). CONCLUSION: Wrist-ankle acupuncture can effectively relieve pain and hand swelling, improve motor function of upper extremity and self-care ability of daily life in patients with shoulder-hand syndrome phaseⅠafter stroke, the therapeutic effect is superior to conventional acupuncture.

Cell landscape of larval and adult Xenopus laevis at single-cell resolution.

The rapid development of high-throughput single-cell RNA sequencing technology offers a good opportunity to dissect cell heterogeneity of animals. A large number of organism-wide single-cell atlases have been constructed for vertebrates such as Homo sapiens, Macaca fascicularis, Mus musculus and Danio rerio. However, an intermediate taxon that links mammals to vertebrates of more ancient origin is still lacking. Here, we construct the first Xenopus cell landscape to date, including larval and adult organs. Common cell lineage-specific transcription factors have been identified in vertebrates, including fish, amphibians and mammals. The comparison of larval and adult erythrocytes identifies stage-specific hemoglobin subtypes, as well as a common type of cluster containing both larval and adult hemoglobin, mainly at NF59. In addition, cell lineages originating from all three layers exhibits both antigen processing and presentation during metamorphosis, indicating a common regulatory mechanism during metamorphosis. Overall, our study provides a large-scale resource for research on Xenopus metamorphosis and adult organs.

Construction of the axolotl cell landscape using combinatorial hybridization sequencing at single-cell resolution.

The Mexican axolotl (Ambystoma mexicanum) is a well-established tetrapod model for regeneration and developmental studies. Remarkably, neotenic axolotls may undergo metamorphosis, a process that triggers many dramatic changes in diverse organs, accompanied by gradually decline of their regeneration capacity and lifespan. However, the molecular regulation and cellular changes in neotenic and metamorphosed axolotls are still poorly investigated. Here, we develop a single-cell sequencing method based on combinatorial hybridization to generate a tissue-based transcriptomic landscape of the neotenic and metamorphosed axolotls. We perform gene expression profiling of over 1 million single cells across 19 tissues to construct the first adult axolotl cell landscape. Comparison of single-cell transcriptomes between the tissues of neotenic and metamorphosed axolotls reveal the heterogeneity of non-immune parenchymal cells in different tissues and established their regulatory network. Furthermore, we describe dynamic gene expression patterns during limb development in neotenic axolotls. This system-level single-cell analysis of molecular characteristics in neotenic and metamorphosed axolotls, serves as a resource to explore the molecular identity of the axolotl and facilitates better understanding of metamorphosis.

ChIP-Hub provides an integrative platform for exploring plant regulome.

Plant genomes encode a complex and evolutionary diverse regulatory grammar that forms the basis for most life on earth. A wealth of regulome and epigenome data have been generated in various plant species, but no common, standardized resource is available so far for biologists. Here, we present ChIP-Hub, an integrative web-based platform in the ENCODE standards that bundles >10,000 publicly available datasets reanalyzed from >40 plant species, allowing visualization and meta-analysis. We manually curate the datasets through assessing ~540 original publications and comprehensively evaluate their data quality. As a proof of concept, we extensively survey the co-association of different regulators and construct a hierarchical regulatory network under a broad developmental context. Furthermore, we show how our annotation allows to investigate the dynamic activity of tissue-specific regulatory elements (promoters and enhancers) and their underlying sequence grammar. Finally, we analyze the function and conservation of tissue-specific promoters, enhancers and chromatin states using comparative genomics approaches. Taken together, the ChIP-Hub platform and the analysis results provide rich resources for deep exploration of plant ENCODE. ChIP-Hub is available at https://biobigdata.nju.edu.cn/ChIPHub/ .

Rice LEAFY COTYLEDON1 Hinders Embryo Greening During the Seed Development.

LEAFY COTYLEDON1 (LEC1) is the central regulator of seed development in Arabidopsis, while its function in monocots is largely elusive. We generated Oslec1 mutants using CRISPR/Cas9 technology. Oslec1 mutant seeds lost desiccation tolerance and triggered embryo greening at the early development stage. Transcriptome analysis demonstrated that Oslec1 mutation altered diverse hormonal pathways and stress response in seed maturation, and promoted a series of photosynthesis-related genes. Further, genome-wide identification of OsLEC1-binding sites demonstrated that OsLEC1 bound to genes involved in photosynthesis, photomorphogenesis, as well as abscisic acid (ABA) and gibberellin (GA) pathways, involved in seed maturation. We illustrated an OsLEC1-regulating gene network during seed development, including the interconnection between photosynthesis and ABA/GA biosynthesis/signaling. Our findings suggested that OsLEC1 acts as not only a central regulator of seed maturation but also an inhibitor of embryo greening during rice seed development. This study would provide new understanding for the OsLEC1 regulatory mechanisms on photosynthesis in the monocot seed development.

Biogenesis, Functions, Interactions, and Resources of Non-Coding RNAs in Plants.

Plant transcriptomes encompass a large number of functional non-coding RNAs (ncRNAs), only some of which have protein-coding capacity. Since their initial discovery, ncRNAs have been classified into two broad categories based on their biogenesis and mechanisms of action, housekeeping ncRNAs and regulatory ncRNAs. With advances in RNA sequencing technology and computational methods, bioinformatics resources continue to emerge and update rapidly, including workflow for in silico ncRNA analysis, up-to-date platforms, databases, and tools dedicated to ncRNA identification and functional annotation. In this review, we aim to describe the biogenesis, biological functions, and interactions with DNA, RNA, protein, and microorganism of five major regulatory ncRNAs (miRNA, siRNA, tsRNA, circRNA, lncRNA) in plants. Then, we systematically summarize tools for analysis and prediction of plant ncRNAs, as well as databases. Furthermore, we discuss the silico analysis process of these ncRNAs and present a protocol for step-by-step computational analysis of ncRNAs. In general, this review will help researchers better understand the world of ncRNAs at multiple levels.