Scalable identification of lineage-specific gene regulatory networks from metacells with NetID.

The identification of gene regulatory networks (GRNs) is crucial for understanding cellular differentiation. Single-cell RNA sequencing data encode gene-level covariations at high resolution, yet data sparsity and high dimensionality hamper accurate and scalable GRN reconstruction. To overcome these challenges, we introduce NetID leveraging homogenous metacells while avoiding spurious gene-gene correlations. Benchmarking demonstrates superior performance of NetID compared to imputation-based methods. By incorporating cell fate probability information, NetID facilitates the prediction of lineage-specific GRNs and recovers known network motifs governing bone marrow hematopoiesis, making it a powerful toolkit for deciphering gene regulatory control of cellular differentiation from large-scale single-cell transcriptome data.

The landscape of RNA binding proteins in mammalian spermatogenesis.

Despite continuous expansion of the RNA binding protein (RBP) world, there is a lack of systematic understanding of RBPs in the mammalian testis, which harbors one of the most complex tissue transcriptomes. We adapted RNA interactome capture to mouse male germ cells, building an RBP atlas characterized by multiple layers of dynamics along spermatogenesis. Trapping of RNA-cross-linked peptides showed that the glutamic acid-arginine (ER) patch, a residue-coevolved polyampholytic element present in coiled coils, enhances RNA binding of its host RBPs. Deletion of this element in NONO (non-POU domain-containing octamer-binding protein) led to a defective mitosis-to-meiosis transition due to compromised NONO-RNA interactions. Whole-exome sequencing of over 1000 infertile men revealed a prominent role of RBPs in the human genetic architecture of male infertility and identified risk ER patch variants.

Design, synthesis, and biological evaluation of novel benzimidazole derivatives as anti-cervical cancer agents through PI3K/Akt/mTOR pathway and tubulin inhibition.

Phosphatidylinositol 3-kinase (PI3K) is one of the most attractive therapeutic targets for cervical cancer treatment. In this study, we designed and synthesized a series of benzimidazole derivatives and evaluated their anti-cervical cancer activity. Compound 4r exhibited strong antiproliferative activity in different cervical cancer cell lines HeLa, SiHa and Ca Ski, and relative lower cytotoxicity to normal hepatic and renal cell lines LO2 and HEK-293t (IC50 values were at 21.08 µM and 23.96 µM respectively). Its IC50 value was at 3.38 µM to the SiHa cells. Further mechanistic studies revealed that 4r induced apoptosis, arrested cell cycle in G2/M phase, suppressed PI3K/Akt/mTOR pathway and inhibit the polymerization of tubulin. Molecular docking study suggested that 4r formed key H-bonds action with PI3Kα (PDB ID:8EXU) and tubulin (PDB ID:1SA0). Zebrafish acute toxicity experiments showed that high concentrations of 4r did not cause death or malformation of zebrafish embryos. All these results demonstrated that 4r would be a promising lead candidate for further development of novel PI3K and tubulin dual inhibitors in cervical cancer treatment.

Comparative assessment of orthodontic clear aligner versus fixed appliance for anterior retraction: a finite element study.

BACKGROUND: The aim of this study is to conduct a comparative evaluation of different designs of clear aligners and examine the disparities between clear aligners and fixed appliances. METHODS: 3D digital models were created, consisting of a maxillary dentition without first premolars, maxilla, periodontal ligaments, attachments, micro-implant, 3D printed lingual retractor, brackets, archwire and clear aligner. The study involved the creation of five design models for clear aligner maxillary anterior internal retraction and one design model for fixed appliance maxillary anterior internal retraction, which were subsequently subjected to finite element analysis. These design models included: (1) Model C0 Control, (2) Model C1 Posterior Micro-implant, (3) Model C2 Anterior Micro-implant, (4) Model C3 Palatal Plate, (5) Model C4 Lingual Retractor, and (6) Model F0 Fixed Appliance. RESULTS: In the clear aligner models, a consistent pattern of tooth movement was observed. Notably, among all tested models, the modified clear aligner Model C3 exhibited the smallest differences in sagittal displacement of the crown-root of the central incisor, vertical displacement of the central incisor, sagittal displacement of the second premolar and second molar, as well as vertical displacement of posterior teeth. However, distinct variations in tooth movement trends were observed between the clear aligner models and the fixed appliance model. Furthermore, compared to the fixed appliance model, significant increases in tooth displacement were achieved with the use of clear aligner models. CONCLUSIONS: In the clear aligner models, the movement trend of the teeth remained consistent, but there were variations in the amount of tooth displacement. Overall, the Model C3 exhibited better torque control and provided greater protection for posterior anchorage teeth compared to the other four clear aligner models. On the other hand, the fixed appliance model provides superior anterior torque control and better protection of the posterior anchorage teeth compared to clear aligner models. The clear aligner approach and the fixed appliance approach still exhibit a disparity; nevertheless, this study offers a developmental direction and establishes a theoretical foundation for future non-invasive, aesthetically pleasing, comfortable, and efficient modalities of clear aligner treatment.

Approximate estimation of cell-type resolution transcriptome in bulk tissue through matrix completion.

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for uncovering cellular heterogeneity. However, the high costs associated with this technique have rendered it impractical for studying large patient cohorts. We introduce ENIGMA (Deconvolution based on Regularized Matrix Completion), a method that addresses this limitation through accurately deconvoluting bulk tissue RNA-seq data into a readout with cell-type resolution by leveraging information from scRNA-seq data. By employing a matrix completion strategy, ENIGMA minimizes the distance between the mixture transcriptome obtained with bulk sequencing and a weighted combination of cell-type-specific expression. This allows the quantification of cell-type proportions and reconstruction of cell-type-specific transcriptomes. To validate its performance, ENIGMA was tested on both simulated and real datasets, including disease-related tissues, demonstrating its ability in uncovering novel biological insights.

Tox4 regulates transcriptional elongation and reinitiation during murine T cell development.

HMG protein Tox4 is a regulator of PP1 phosphatases with unknown function in development. Here we show that Tox4 conditional knockout in mice reduces thymic cellularity, partially blocks T cell development, and decreases ratio of CD8 to CD4 through decreasing proliferation and increasing apoptosis of CD8 cells. In addition, single-cell RNA-seq discovered that Tox4 loss also impairs proliferation of the fast-proliferating double positive (DP) blast population within DP cells in part due to downregulation of genes critical for proliferation, notably Cdk1. Moreover, genes with high and low expression level are more dependent on Tox4 than genes with medium expression level. Mechanistically, Tox4 may facilitate transcriptional reinitiation and restrict elongation in a dephosphorylation-dependent manner, a mechanism that is conserved between mouse and human. These results provide insights into the role of TOX4 in development and establish it as an evolutionarily conserved regulator of transcriptional elongation and reinitiation.

Single-cell RNA-seq uncovers dynamic processes orchestrated by RNA-binding protein DDX43 in chromatin remodeling during spermiogenesis.

Mammalian spermatogenesis shows prominent chromatin and transcriptomic switches in germ cells, but it is unclear how such dynamics are controlled. Here we identify RNA helicase DDX43 as an essential regulator of the chromatin remodeling process during spermiogenesis. Testis-specific Ddx43 knockout mice show male infertility with defective histone-to-protamine replacement and post-meiotic chromatin condensation defects. The loss of its ATP hydrolysis activity by a missense mutation replicates the infertility phenotype in global Ddx43 knockout mice. Single-cell RNA sequencing analyses of germ cells depleted of Ddx43 or expressing the Ddx43 ATPase-dead mutant reveals that DDX43 regulates dynamic RNA regulatory processes that underlie spermatid chromatin remodeling and differentiation. Transcriptomic profiling focusing on early-stage spermatids combined with enhanced crosslinking immunoprecipitation and sequencing further identifies Elfn2 as DDX43-targeted hub gene. These findings illustrate an essential role for DDX43 in spermiogenesis and highlight the single-cell-based strategy to dissect cell-state-specific regulation of male germline development.

Assessment of labially impacted canines traction mode with clear aligners vs. fixed appliance: A comparative study based on 3D finite element analysis.

Purpose: The objective of this study was to evaluate and compare the biomechanical differences between clear aligner and fixed appliance in the traction of labially impacted canines based on 3D finite element analysis. Methods: A series of patient-oriented finite element models were constructed, including a maxillary dentition with a right labially canine, maxilla, periodontal ligaments, traction attachments, and clear aligners. The two most common clinical scenarios were investigated: Scenario A: impacted canine (distal) and Scenario B: impacted canine (mesial). For each clinical scenario, three traction models with clear aligners and one fixed appliance model were established. Results: In all four models, the impacted canines exhibited similar initial displacement tendencies of mesially rotated in Scenario A and distally rotated in Scenario B, and with small differences in periodontal ligament stress magnitude. However, the sum of the periodontal ligament stresses of the anchorage teeth in the clear aligner mode was in the range of 56.28-76.21 kPa and in the fixed appliance mode was in the range of 6.61-7.22 kPa. The maximum value of initial displacement of the anchorage teeth in the clear aligner mode was in the range of 13.71-19.72 µm, while in the fixed appliance mode was 3.10-3.92 µm. Conclusion: For impacted canines, clear aligner mode and fixed appliance mode have little difference in biomechanical effect. However, the anchorage teeth in the clear aligner mode endure higher stress and show a more pronounced displacement tendency. In addition, the biomechanical effects of different clear aligner traction models are various but not obvious.

Cancer-associated dynamics and potential regulators of intronic polyadenylation revealed by IPAFinder using standard RNA-seq data.

Intronic polyadenylation (IpA) usually leads to changes in the coding region of an mRNA, and its implication in diseases has been recognized, although at its very beginning status. Conveniently and accurately identifying IpA is of great importance for further evaluating its biological significance. Here, we developed IPAFinder, a bioinformatic method for the de novo identification of intronic poly(A) sites and their dynamic changes from standard RNA-seq data. Applying IPAFinder to 256 pan-cancer tumor/normal pairs across six tumor types, we discovered 490 recurrent dynamically changed IpA events, some of which are novel and derived from cancer-associated genes such as TSC1, SPERD2, and CCND2 Furthermore, IPAFinder revealed that IpA could be regulated by factors related to splicing and m6A modification. In summary, IPAFinder enables the global discovery and characterization of biologically regulated IpA with standard RNA-seq data and should reveal the biological significance of IpA in various processes.

Comprehensive characterization of somatic variants associated with intronic polyadenylation in human cancers.

Somatic single nucleotide variants (SNVs) in cancer genome affect gene expression through various mechanisms depending on their genomic location. While somatic SNVs near canonical splice sites have been reported to cause abnormal splicing of cancer-related genes, whether these SNVs can affect gene expression through other mechanisms remains an open question. Here, we analyzed RNA sequencing and exome data from 4,998 cancer patients covering ten cancer types and identified 152 somatic SNVs near splice sites that were associated with abnormal intronic polyadenylation (IPA). IPA-associated somatic variants favored the localization near the donor splice sites compared to the acceptor splice sites. A proportion of SNV-associated IPA events overlapped with premature cleavage and polyadenylation events triggered by U1 small nuclear ribonucleoproteins (snRNP) inhibition. GC content, intron length and polyadenylation signal were three genomic features that differentiated between SNV-associated IPA and intron retention. Notably, IPA-associated SNVs were enriched in tumor suppressor genes (TSGs), including the well-known TSGs such as PTEN and CDH1 with recurrent SNV-associated IPA events. Minigene assay confirmed that SNVs from PTEN, CDH1, VEGFA, GRHL2, CUL3 and WWC2 could lead to IPA. This work reveals that IPA acts as a novel mechanism explaining the functional consequence of somatic SNVs in human cancer.

Independent component analysis based gene co-expression network inference (ICAnet) to decipher functional modules for better single-cell clustering and batch integration.

With the tremendous increase of publicly available single-cell RNA-sequencing (scRNA-seq) datasets, bioinformatics methods based on gene co-expression network are becoming efficient tools for analyzing scRNA-seq data, improving cell type prediction accuracy and in turn facilitating biological discovery. However, the current methods are mainly based on overall co-expression correlation and overlook co-expression that exists in only a subset of cells, thus fail to discover certain rare cell types and sensitive to batch effect. Here, we developed independent component analysis-based gene co-expression network inference (ICAnet) that decomposed scRNA-seq data into a series of independent gene expression components and inferred co-expression modules, which improved cell clustering and rare cell-type discovery. ICAnet showed efficient performance for cell clustering and batch integration using scRNA-seq datasets spanning multiple cells/tissues/donors/library types. It works stably on datasets produced by different library construction strategies and with different sequencing depths and cell numbers. We demonstrated the capability of ICAnet to discover rare cell types in multiple independent scRNA-seq datasets from different sources. Importantly, the identified modules activated in acute myeloid leukemia scRNA-seq datasets have the potential to serve as new diagnostic markers. Thus, ICAnet is a competitive tool for cell clustering and biological interpretations of single-cell RNA-seq data analysis.

Genome-wide study of C2H2 zinc finger gene family in Medicago truncatula.

BACKGROUND: C2H2 zinc finger proteins (C2H2 ZFPs) play vital roles in shaping many aspects of plant growth and adaptation to the environment. Plant genomes harbor hundreds of C2H2 ZFPs, which compose one of the most important and largest transcription factor families in higher plants. Although the C2H2 ZFP gene family has been reported in several plant species, it has not been described in the model leguminous species Medicago truncatula. RESULTS: In this study, we identified 218 C2H2 type ZFPs with 337 individual C2H2 motifs in M. truncatula. We showed that the high rate of local gene duplication has significantly contributed to the expansion of the C2H2 gene family in M. truncatula. The identified ZFPs exhibit high variation in motif arrangement and expression pattern, suggesting that the short C2H2 zinc finger motif has been adopted as a scaffold by numerous transcription factors with different functions to recognize cis-elements. By analyzing the public expression datasets and quantitative RT-PCR (qRT-PCR), we identified several C2H2 ZFPs that are specifically expressed in certain tissues, such as the nodule, seed, and flower. CONCLUSION: Our genome-wide work revealed an expanded C2H2 ZFP gene family in an important legume M. truncatula, and provides new insights into the diversification and expansion of C2H2 ZFPs in higher plants.

Pan-cancer analysis identifies telomerase-associated signatures and cancer subtypes.

BACKGROUND: Cancer cells become immortalized through telomere maintenance mechanisms, such as telomerase reverse transcriptase (TERT) activation. In addition to maintaining telomere length, TERT activates manifold cell survival signaling pathways. However, telomerase-associated gene signatures in cancer remain elusive. METHODS: We performed a systematic analysis of TERT high (TERThigh) and low (TERTlow) cancers using multidimensional data from The Cancer Genome Atlas (TCGA). Multidimensional data were analyzed by propensity score matching weight algorithm. Coexpression networks were constructed by weight gene coexpression network analysis (WGCNA). Random forest classifiers were generated to identify cancer subtypes. RESULTS: The TERThigh-specific mRNA expression signature is associated with cell cycle-related coexpression modules across cancer types. Experimental screening of hub genes in the cell cycle module suggested TPX2 and EXO1 as potential regulators of telomerase activity and cell survival. MiRNA analysis revealed that the TERThigh-specific miR-17-92 cluster can target biological processes enriched in TERTlow cancer and that its expression is negatively correlated with the tumor/normal telomere length ratio. Intriguingly, TERThigh cancers tend to have mutations in extracellular matrix organization genes and amplify MAPK signaling. By mining the clinical actionable gene database, we uncovered a number of TERThigh-specific somatic mutations, amplifications and high expression genes containing therapeutic targets. Finally, a random forest classifier integrating telomerase-associated multi-omics signatures identifies two cancer subtypes showed profound differences in telomerase activity and patient survival. CONCLUSIONS: In summary, our results depict a telomerase-associated molecular landscape in cancers and provide therapeutic opportunities for cancer treatment.

High autophagic flux guards ESC identity through coordinating autophagy machinery gene program by FOXO1.

Although much is known about transcriptional networks that control embryonic stem cell (ESC) self-renewal and differentiation, the metabolic regulation of ESC is less clear. Autophagy is a catabolic process that is activated under both stress and normal conditions to degrade damaged organelles and aggregated proteins, and thus plays pivotal roles in somatic and adult stem cell function. However, if and how ESCs harness autophagy to regulate stemness remains largely unknown. Recently, we have defined that autophagy is essential for mitochondrial homeostasis regulation in pluripotency acquirement and maintenance. Here we identified high autophagic flux as an essential mechanism to maintain ESC identity. We show that mouse ESCs exhibit a high autophagic flux that is maintained by coordinating expression of autophagy core molecular machinery genes through FOXO1, a forkhead family transcription factor. Tapering autophagic flux by manipulating either Atg3 or Foxo1 expression compromised ESC self-renewal, pluripotency, and differentiation that could be restored by gain of wild-type but not function-deficient Atg3 or Foxo1 mutants, respectively. Our results define a newly recognized role of autophagic flux in mouse ESC identity maintenance that links cellular catabolism to ESC fate regulation.

ATG3-dependent autophagy mediates mitochondrial homeostasis in pluripotency acquirement and maintenance.

Pluripotent stem cells, including induced pluripotent and embryonic stem cells (ESCs), have less developed mitochondria than somatic cells and, therefore, rely more heavily on glycolysis for energy production. 1-3 However, how mitochondrial homeostasis matches the demands of nuclear reprogramming and regulates pluripotency in ESCs is largely unknown. Here, we identified ATG3-dependent autophagy as an executor for both mitochondrial remodeling during somatic cell reprogramming and mitochondrial homeostasis regulation in ESCs. Dysfunctional autophagy by Atg3 deletion inhibited mitochondrial removal during pluripotency induction, resulting in decreased reprogramming efficiency and accumulation of abnormal mitochondria in established iPSCs. In Atg3 null mouse ESCs, accumulation of aberrant mitochondria was accompanied by enhanced ROS generation, defective ATP production and attenuated pluripotency gene expression, leading to abnormal self-renewal and differentiation. These defects were rescued by reacquisition of wild-type but not lipidation-deficient Atg3 expression. Taken together, our findings highlight a critical role of ATG3-dependent autophagy for mitochondrial homeostasis regulation in both pluripotency acquirement and maintenance.

Treatment of multiple sclerosis by transplantation of neural stem cells derived from induced pluripotent stem cells.

Multiple sclerosis (MS) is an autoimmune disease of the central nervous system (CNS), with focal T lymphocytic infiltration and damage of myelin and axons. The underlying mechanism of pathogenesis remains unclear and there are currently no effective treatments. The development of neural stem cell (NSC) transplantation provides a promising strategy to treat neurodegenerative disease. However, the limited availability of NSCs prevents their application in neural disease therapy. In this study, we generated NSCs from induced pluripotent stem cells (iPSCs) and transplanted these cells into mice with experimental autoimmune encephalomyelitis (EAE), a model of MS. The results showed that transplantation of iPSC-derived NSCs dramatically reduced T cell infiltration and ameliorated white matter damage in the treated EAE mice. Correspondingly, the disease symptom score was greatly decreased, and motor ability was dramatically rescued in the iPSC-NSC-treated EAE mice, indicating the effectiveness of using iPSC-NSCs to treat MS. Our study provides pre-clinical evidence to support the feasibility of treating MS by transplantation of iPSC-derived NSCs.

Immunogenicity and functional evaluation of iPSC-derived organs for transplantation.

Whether physiologically induced pluripotent stem cell (iPSC)-derived organs are immunogenic and can be used for transplantation is unclear. Here, we generated iPSC-derived skin, islet, and heart representing three germ layers of the body through 4n complementation and evaluated their immunogenicity and therapeutic efficacy. Upon transplantation into recipient mice, iPSC-derived skin successfully survived and repaired local tissue wounds. In diabetic mouse models, explanted iPSC-derived islets effectively produced insulin and lowered blood glucose to basal levels. iPSC-derived heart grafts maintained normal beating for more than 3 months in syngeneic recipients. Importantly, no obvious immune rejection responses against iPSC-derived organs were detected long after transplantation. Our study not only demonstrates the fundamental immunogenicity and function of iPSC derivatives, but also provides preclinical evidence to support the feasibility of using iPSC-derived skin, islet, and heart for therapeutic use.