Transcriptomic profiling of Dip2a in the neural differentiation of mouse embryonic stem cells.

Introduction: The disconnected-interacting protein 2 homolog A (DIP2A), a member of disconnected-interacting 2 protein family, has been shown to be involved in human nervous system-related mental illness. This protein is highly expressed in the nervous system of mouse. Mutation of mouse DIP2A causes defects in spine morphology and synaptic transmission, autism-like behaviors, and defective social novelty [5], [27], indicating that DIP2A is critical to the maintenance of neural development. However, the role of DIP2A in neural differentiation has yet to be investigated. Objective: To determine the role of DIP2A in neural differentiation, a neural differentiation model was established using mouse embryonic stem cells (mESCs) and studied by using gene-knockout technology and RNA-sequencing-based transcriptome analysis. Results: We found that DIP2A is not required for mESCs pluripotency maintenance, but loss of DIP2A causes the neural differentiation abnormalities in both N2B27 and KSR medium. Functional knockout of Dip2a gene also decreased proliferation of mESCs by perturbation of the cell cycle and profoundly inhibited the expression of a large number of neural development-associated genes which mainly enriched in spinal cord development and postsynapse assembly. Conclusions: The results of this report demonstrate that DIP2A plays an essential role in regulating differentiation of mESCs towards the neural fate.

CAPG interference induces apoptosis and ferroptosis in colorectal cancer cells through the P53 pathway.

PURPOSE: Given the high incidence and mortality rates of colorectal cancer (CRC) and the inadequacy of existing treatments for many patients, this study aimed to explore the potential of Capping Actin Protein (CAPG), a protein involved in actin-related movements, as a novel therapeutic target for CRC. METHODS: Bioinformatic analysis of gene expression was conducted using the UALCAN website. Cell proliferation was measured using the CCK-8 kit. Cell cycle, apoptosis, and ferroptosis were analyzed using flow cytometry. Tumorigenesis was evaluated by the subcutaneous inoculation of CRC cells into BALB/c nude female mice. Differentially expressed genes and signaling pathways were identified using RNA sequencing. RESULTS: CAPG was significantly overexpressed in human CRC tissues and its upregulation was correlated with poor overall survival. CAPG knockdown led to notable inhibition of CRC cells in vitro and in vivo. Interference with CAPG blocked the cell cycle at the G1 phase and triggered apoptosis and ferroptosis by upregulating the P53 pathway in CRC cells. CONCLUSION: CRC patients with higher CAPG levels have a poorer prognosis. CAPG inhibits apoptosis and ferroptosis, while promoting CRC cell proliferation by repressing the P53 pathway. Our study suggests that CAPG may be a potential therapeutic target for CRC prognosis and treatment.

Virulence and Adhesion of the Treponema pallidum Nichols Strain Simultaneously Decrease in a Continuous-Infection New Zealand White Rabbit Model.

Syphilis is a sexually transmitted disease caused by T. pallidum, and the T. pallidum Nichols strain is widely used with the New Zealand white rabbit model for evaluating drug and vaccine protection. However, changes in the virulence of T. pallidum during transmission are still unknown. Herein, we explored the virulence of T. pallidum in the rabbit model of continuous infection through phenotype observation and further investigated the relationship between virulence and adhesion. During the construction of the syphilis rabbit model, the optimal dose of 104/site of T. pallidum was determined to effectively observe the depiction of syphilis lesions and immune responses for further virulence evaluation. Its virulence was gradually weakened during the interaction with host cells or the testicular passage, which was also proven using the pathological phenotype of the syphilis rabbit model. In addition, the adhesive ability of T. pallidum was reduced with increasing generation, which was verified via the co-incubation of the pathogen with Sf1Ep cells. This study provides insight into the relationship by which the virulence and adhesion of T. pallidum were decreased in a New Zealand white rabbit model of continuous infection and contributes to our knowledge regarding the development of syphilis.

Loss-of-function mutations in CST6 cause dry skin, desquamation and abnormal keratosis without hypotrichosis.

Cystatin M/E (encoded by the CST6 gene) is a cysteine protease inhibitor, that exerts regulatory and protective effects against uncontrolled proteolysis mainly by directly regulating cathepsin V, cathepsin L, and legumain activities. Previous studies have suggested that CST6 may exert a regulatory role in epidermal differentiation and hair follicle formation by inhibiting the activity of respective cognate target proteases. However, until recently, studies have revealed that loss- or gain-of-function of the CST6 gene causes dry skin with hypotrichosis in humans. Here, we reported two siblings of Chinese origin with dry skin, desquamation and abnormal keratosis without hypotrichosis. By applying whole-exome sequencing, we identified homozygous loss-of-function mutation c.251G > A (p.Gly84Asp) in the CST6 gene as the underlying genetic cause. Further fluorimetric enzyme assays demonstrated the mutant cystatin M/E protein lost its inhibitory function on the protease activity of cathepsins. Moreover, the corresponding mutation in mice resulted in excessive cornification, desquamation, impaired skin barrier function, and abnormal proliferation and differentiation of keratinocytes. In conclusion, the homozygous missense mutation c.251G > A in CST6 gene resulted in dry skin, desquamation, as well as abnormal keratosis of the skin, promoting our understanding of the role of protease-antiprotease balance in human skin disorders.

NanoDeep: a deep learning framework for nanopore adaptive sampling on microbial sequencing.

Nanopore sequencers can enrich or deplete the targeted DNA molecules in a library by reversing the voltage across individual nanopores. However, it requires substantial computational resources to achieve rapid operations in parallel at read-time sequencing. We present a deep learning framework, NanoDeep, to overcome these limitations by incorporating convolutional neural network and squeeze and excitation. We first showed that the raw squiggle derived from native DNA sequences determines the origin of microbial and human genomes. Then, we demonstrated that NanoDeep successfully classified bacterial reads from the pooled library with human sequence and showed enrichment for bacterial sequence compared with routine nanopore sequencing setting. Further, we showed that NanoDeep improves the sequencing efficiency and preserves the fidelity of bacterial genomes in the mock sample. In addition, NanoDeep performs well in the enrichment of metagenome sequences of gut samples, showing its potential applications in the enrichment of unknown microbiota. Our toolkit is available at https://github.com/lysovosyl/NanoDeep.

Reciprocal Regulation between lncRNA ANRIL and p15 in Steroid-Induced Glaucoma.

Steroid-induced glaucoma (SIG) is the most common adverse steroid-related effect on the eyes. SIG patients can suffer from trabecular meshwork (TM) dysfunction, intraocular pressure (IOP) elevation, and irreversible vision loss. Previous studies have mainly focused on the role of extracellular matrix turnover in TM dysfunction; however, whether the cellular effects of TM cells are involved in the pathogenesis of SIG remains unclear. Here, we found that the induction of cellular senescence was associated with TM dysfunction, causing SIG in cultured cells and mouse models. Especially, we established the transcriptome landscape in the TM tissue of SIG mice via microarray screening and identified ANRIL as the most differentially expressed long non-coding RNA, with a 5.4-fold change. The expression level of ANRIL was closely related to ocular manifestations (IOP elevation, cup/disc ratio, and retinal nerve fiber layer thickness). Furthermore, p15, the molecular target of ANRIL, was significantly upregulated in SIG and was correlated with ocular manifestations in an opposite direction to ANRIL. The reciprocal regulation between ANRIL and p15 was validated using luciferase reporter assay. Through depletion in cultured cells and a mouse model, ANRIL/p15 signaling was confirmed in cellular senescence via cyclin-dependent kinase activity and, subsequently, by phosphorylation of the retinoblastoma protein. ANRIL depletion imitated the SIG phenotype, most importantly IOP elevation. ANRIL depletion-induced IOP elevation in mice can be effectively suppressed by p15 depletion. Analyses of the single-cell atlas and transcriptome dynamics of human TM tissue showed that ANRIL/p15 expression is spatially enriched in human TM cells and is correlated with TM dysfunction. Moreover, ANRIL is colocalized with a GWAS risk variant (rs944800) of glaucoma, suggesting its potential role underlying genetic susceptibility of glaucoma. Together, our findings suggested that steroid treatment promoted cellular senescence, which caused TM dysfunction, IOP elevation, and irreversible vision loss. Molecular therapy targeting the ANRIL/p15 signal exerted a protective effect against steroid treatment and shed new light on glaucoma management.

Engineered Cas12a-Plus nuclease enables gene editing with enhanced activity and specificity.

BACKGROUND: The CRISPR-Cas12a (formerly Cpf1) system is a versatile gene-editing tool with properties distinct from the broadly used Cas9 system. Features such as recognition of T-rich protospacer-adjacent motif (PAM) and generation of sticky breaks, as well as amenability for multiplex editing in a single crRNA and lower off-target nuclease activity, broaden the targeting scope of available tools and enable more accurate genome editing. However, the widespread use of the nuclease for gene editing, especially in clinical applications, is hindered by insufficient activity and specificity despite previous efforts to improve the system. Currently reported Cas12a variants achieve high activity with a compromise of specificity. Here, we used structure-guided protein engineering to improve both editing efficiency and targeting accuracy of Acidaminococcus sp. Cas12a (AsCas12a) and Lachnospiraceae bacterium Cas12a (LbCas12a). RESULTS: We created new AsCas12a variant termed "AsCas12a-Plus" with increased activity (1.5~2.0-fold improvement) and specificity (reducing off-targets from 29 to 23 and specificity index increased from 92% to 94% with 33 sgRNAs), and this property was retained in multiplex editing and transcriptional activation. When used to disrupt the oncogenic BRAFV600E mutant, AsCas12a-Plus showed less off-target activity while maintaining comparable editing efficiency and BRAFV600E cancer cell killing. By introducing the corresponding substitutions into LbCas12a, we also generated LbCas12a-Plus (activity improved ~1.1-fold and off-targets decreased from 20 to 12 while specificity index increased from 78% to 89% with 15 sgRNAs), suggesting this strategy may be generally applicable across Cas12a orthologs. We compared Cas12a-Plus, other variants described in this study, and the reported enCas12a-HF, enCas12a, and Cas12a-ultra, and found that Cas12a-Plus outperformed other variants with a good balance for enhanced activity and improved specificity. CONCLUSIONS: Our discoveries provide alternative AsCas12a and LbCas12a variants with high specificity and activity, which expand the gene-editing toolbox and can be more suitable for clinical applications.

Immune Cell Infiltration Analysis Demonstrates Excessive Mast Cell Activation in Psoriasis.

Psoriasis represents multiple inflammatory processes and exaggerated physiological responses to epithelial damage by innate and adaptive immune components, thus it is critical to compare the immune cell niche in disease and healthy skin. Here, we inferred the proportions of different immune cell types in psoriatic and healthy skin using the CIBERSORT algorithm with expression profiles as input. As a result, we observed a dramatic change of immune cell profiles in psoriatic skin compared with healthy skin. Interestingly, the resting mast cells is almost eliminated in psoriatic skin. In contrast, the activated mast cells are enriched in psoriatic skin, indicating that mast cells activation may play an important role in psoriasis pathogenesis. In addition, we found that the proportion of the resting mast cells gradually come back to the normal level in lesioned skin upon etanercept treatment, suggesting that mast cells play a critical role in immune cell niche maintenance. Further experiments validated a significant decrease in mast cell population and an excessive mast cell activation in psoriatic skin compared with healthy skin. In conclusion, our integrative analyses of the immune cell profiles and the corresponding marker genes expression provide a better understanding of the inflammation response in psoriasis and important clues for clinical applications.

Cabazitaxel suppresses colorectal cancer cell growth via enhancing the p53 antitumor pathway.

There were approximately 1.93 million new cases and 940 000 deaths from colorectal cancer in 2020. The first-line chemotherapeutic drugs for colorectal cancer are mainly based on 5-fluorouracil, although the use of these drugs is limited by the development of drug resistance. Consequently, there is a need for novel chemotherapeutic drugs for the efficient treatment of colorectal cancer patients. In the present study, we screened 160 drugs approved by the Food and Drug Administration and identified that cabazitaxel (CBT), a microtube inhibitor, can suppress colony formation and cell migration of colorectal cancer cells in vitro. CBT also induces G2/M phase arrest and apoptosis of colorectal cancer cells. Most importantly, it inhibits the growth of colorectal cancer cell xenograft tumors in vivo. Transcriptome analysis by RNA-sequencing revealed that Tub family genes are abnormally expressed in CBT-treated colorectal cancer cells. The expression of several p53 downstream genes that are associated with cell cycle arrest, apoptosis, and inhibition of angiogenesis and metastasis is induced by CBT in colorectal cancer cells. Overall, our results suggests that CBT suppresses colorectal cancer by upregulating the p53 pathway, and thus CBT may have potential as an alternative chemotherapeutic drug for colorectal cancer.

Tag-seq: a convenient and scalable method for genome-wide specificity assessment of CRISPR/Cas nucleases.

Genome-wide identification of DNA double-strand breaks (DSBs) induced by CRISPR-associated protein (Cas) systems is vital for profiling the off-target events of Cas nucleases. However, current methods for off-target discovery are tedious and costly, restricting their widespread applications. Here we present an easy alternative method for CRISPR off-target detection by tracing the integrated oligonucleotide Tag using next-generation-sequencing (CRISPR-Tag-seq, or Tag-seq). Tag-seq enables rapid and convenient profiling of nuclease-induced DSBs by incorporating the optimized double-stranded oligodeoxynucleotide sequence (termed Tag), adapters, and PCR primers. Moreover, we employ a one-step procedure for library preparation in Tag-seq, which can be applied in the routine workflow of a molecular biology laboratory. We further show that Tag-seq successfully determines the cleavage specificity of SpCas9 variants and Cas12a/Cpf1 in a large-scale manner, and discover the integration sites of exogenous genes introduced by the Sleeping Beauty transposon. Our results demonstrate that Tag-seq is an efficient and scalable approach to genome-wide identification of Cas-nuclease-induced off-targets.

Large scale RNA-binding proteins/LncRNAs interaction analysis to uncover lncRNA nuclear localization mechanisms.

Long non-coding RNAs (lncRNAs) are key regulators of major biological processes and their functional modes are dictated by their subcellular localization. Relative nuclear enrichment of lncRNAs compared to mRNAs is a prevalent phenomenon but the molecular mechanisms governing their nuclear retention in cells remain largely unknown. Here in this study, we harness the recently released eCLIP data for a large number of RNA-binding proteins (RBPs) in K562 and HepG2 cells and utilize multiple bioinformatics methods to comprehensively survey the roles of RBPs in lncRNA nuclear retention. We identify an array of splicing RBPs that bind to nuclear-enriched lincRNAs (large intergenic non-coding RNAs) thus may act as trans-factors regulating their nuclear retention. Further analyses reveal that these RBPs may bind with distinct core motifs, flanking sequence compositions, or secondary structures to drive lincRNA nuclear retention. Moreover, network analyses uncover potential co-regulatory RBP clusters and the physical interaction between HNRNPU and SAFB2 proteins in K562 cells is further experimentally verified. Altogether, our analyses reveal previously unknown factors and mechanisms that govern lincRNA nuclear localization in cells.

PARMAP: A Pan-Genome-Based Computational Framework for Predicting Antimicrobial Resistance.

Antimicrobial resistance (AMR) has emerged as one of the most urgent global threats to public health. Accurate detection of AMR phenotypes is critical for reducing the spread of AMR strains. Here, we developed PARMAP (Prediction of Antimicrobial Resistance by MAPping genetic alterations in pan-genome) to predict AMR phenotypes and to identify AMR-associated genetic alterations based on the pan-genome of bacteria by utilizing machine learning algorithms. When we applied PARMAP to 1,597 Neisseria gonorrhoeae strains, it successfully predicted their AMR phenotypes based on a pan-genome analysis. Furthermore, it identified 328 genetic alterations in 23 known AMR genes and discovered many new AMR-associated genetic alterations in ciprofloxacin-resistant N. gonorrhoeae, and it clearly indicated the genetic heterogeneity of AMR genes in different subtypes of resistant N. gonorrhoeae. Additionally, PARMAP performed well in predicting the AMR phenotypes of Mycobacterium tuberculosis and Escherichia coli, indicating the robustness of the PARMAP framework. In conclusion, PARMAP not only precisely predicts the AMR of a population of strains of a given species but also uses whole-genome sequencing data to prioritize candidate AMR-associated genetic alterations based on their likelihood of contributing to AMR. Thus, we believe that PARMAP will accelerate investigations into AMR mechanisms in other human pathogens.

ß-Catenin safeguards the ground state of mousepluripotency by strengthening the robustness of the transcriptional apparatus.

Mouse embryonic stem cells cultured with MEK (mitogen-activated protein kinase kinase) and GSK3 (glycogen synthase kinase 3) inhibitors (2i) more closely resemble the inner cell mass of preimplantation blastocysts than those cultured with SL [serum/leukemia inhibitory factor (LIF)]. The transcriptional mechanisms governing this pluripotent ground state are unresolved. Release of promoter-proximal paused RNA polymerase II (Pol2) is a multistep process necessary for pluripotency and cell cycle gene transcription in SL. We show that ß-catenin, stabilized by GSK3 inhibition in medium with 2i, supplies transcriptional coregulators at pluripotency loci. This selectively strengthens pluripotency loci and renders them addicted to transcription initiation for productive gene body elongation in detriment to Pol2 pause release. By contrast, cell cycle genes are not bound by ß-catenin, and proliferation/self-renewal remains tightly controlled by Pol2 pause release under 2i conditions. Our findings explain how pluripotency is reinforced in the ground state and also provide a general model for transcriptional resilience/adaptation upon network perturbation in other contexts.

Hippo-YAP signaling controls lineage differentiation of mouse embryonic stem cells through modulating the formation of super-enhancers.

Hippo-YAP signaling pathway functions in early lineage differentiation of pluripotent stem cells, but the detailed mechanisms remain elusive. We found that knockout (KO) of Mst1 and Mst2, two key components of the Hippo signaling in mouse embryonic stem cells (ESCs), resulted in a disruption of differentiation into mesendoderm lineage. To further uncover the underlying regulatory mechanisms, we performed a series of ChIP-seq experiments with antibodies against YAP, ESC master transcription factors and some characterized histone modification markers as well as RNA-seq assays using wild type and Mst KO samples at ES and day 4 embryoid body stage respectively. We demonstrate that YAP is preferentially co-localized with super-enhancer (SE) markers such as Nanog, Sox2, Oct4 and H3K27ac in ESCs. The hyper-activation of nuclear YAP in Mst KO ESCs facilitates the binding of Nanog, Sox2 and Oct4 as well as H3K27ac modification at the loci where YAP binds. Moreover, Mst depletion results in novel SE formation and enhanced liquid-liquid phase-separated Med1 condensates on lineage associated genes, leading to the upregulation of these genes and the distortion of ESC differentiation. Our study reveals a novel mechanism on how Hippo-YAP signaling pathway dictates ESC lineage differentiation.

MyoD induced enhancer RNA interacts with hnRNPL to activate target gene transcription during myogenic differentiation.

Emerging evidence supports roles of enhancer RNAs (eRNAs) in regulating target gene. Here, we study eRNA regulation and function during skeletal myoblast differentiation. We provide a panoramic view of enhancer transcription and categorization of eRNAs. Master transcription factor MyoD is crucial in activating eRNA production. Super enhancer (se) generated seRNA-1 and -2 promote myogenic differentiation in vitro and in vivo. seRNA-1 regulates expression levels of two nearby genes, myoglobin (Mb) and apolipoprotein L6 (Apol6), by binding to heterogeneous nuclear ribonucleoprotein L (hnRNPL). A CAAA tract on seRNA-1 is essential in mediating seRNA-1/hnRNPL binding and function. Disruption of seRNA-1-hnRNPL interaction attenuates Pol II and H3K36me3 deposition at the Mb locus, in coincidence with the reduction of its transcription. Furthermore, analyses of hnRNPL binding transcriptome-wide reveal its association with eRNAs is a general phenomenon in multiple cells. Collectively, we propose that eRNA-hnRNPL interaction represents a mechanism contributing to target mRNA activation.

Elevated H3K27ac in aged skeletal muscle leads to increase in extracellular matrix and fibrogenic conversion of muscle satellite cells.

Epigenetic alterations occur in various cells and tissues during aging, but it is not known if such alterations are also associated with aging in skeletal muscle. Here, we examined the changes of a panel of histone modifications and found H3K27ac (an active enhancer mark) is markedly increased in aged human skeletal muscle tissues. Further analyses uncovered that the H3K27ac increase and enhancer activation are associated with the up-regulation of extracellular matrix (ECM) genes; this may result in alteration of the niche environment for skeletal muscle stem cells, also called satellite cells (SCs), which causes decreased myogenic potential and fibrogenic conversion of SCs. In mice, treatment of aging muscles with JQ1, an inhibitor of enhancer activation, inhibited the ECM up-regulation and fibrogenic conversion of SCs and restored their myogenic differentiation potential. Altogether, our findings not only uncovered a novel aspect of skeletal muscle aging that is associated with enhancer remodeling but also implicated JQ1 as a potential treatment approach for restoring SC function in aging muscle.

YY1 regulates skeletal muscle regeneration through controlling metabolic reprogramming of satellite cells.

Skeletal muscle satellite cells (SCs) are adult muscle stem cells responsible for muscle regeneration after acute or chronic injuries. The lineage progression of quiescent SC toward activation, proliferation, and differentiation during the regeneration is orchestrated by cascades of transcription factors (TFs). Here, we elucidate the function of TF Yin Yang1 (YY1) in muscle regeneration. Muscle-specific deletion of YY1 in embryonic muscle progenitors leads to severe deformity of diaphragm muscle formation, thus neonatal death. Inducible deletion of YY1 in SC almost completely blocks the acute damage-induced muscle repair and exacerbates the chronic injury-induced dystrophic phenotype. Examination of SC revealed that YY1 loss results in cell-autonomous defect in activation and proliferation. Mechanistic search revealed that YY1 binds and represses mitochondrial gene expression. Simultaneously, it also stabilizes Hif1α protein and activates Hif1α-mediated glycolytic genes to facilitate a metabolic reprogramming toward glycolysis which is needed for SC proliferation. Altogether, our findings have identified YY1 as a key regulator of SC metabolic reprogramming through its dual roles in modulating both mitochondrial and glycolytic pathways.

SKmDB: an integrated database of next generation sequencing information in skeletal muscle.

MOTIVATION: Skeletal muscles have indispensable functions and also possess prominent regenerative ability. The rapid emergence of Next Generation Sequencing (NGS) data in recent years offers us an unprecedented perspective to understand gene regulatory networks governing skeletal muscle development and regeneration. However, the data from public NGS database are often in raw data format or processed with different procedures, causing obstacles to make full use of them. RESULTS: We provide SKmDB, an integrated database of NGS information in skeletal muscle. SKmDB not only includes all NGS datasets available in the human and mouse skeletal muscle tissues and cells, but also provide preliminary data analyses including gene/isoform expression levels, gene co-expression subnetworks, as well as assembly of putative lincRNAs, typical and super enhancers and transcription factor hotspots. Users can efficiently search, browse and visualize the information with the well-designed user interface and server side. SKmDB thus will offer wet lab biologists useful information to study gene regulatory mechanisms in the field of skeletal muscle development and regeneration. AVAILABILITY AND IMPLEMENTATION: Freely available on the web at http://sunlab.cpy.cuhk.edu.hk/SKmDB. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

lncFunTK: a toolkit for functional annotation of long noncoding RNAs.

Motivation: Thousands of long noncoding RNAs (lncRNAs) were newly identified from high throughput RNA-seq data. Functional annotation and prioritization of these lncRNAs for further experimental validation as well as the functional investigation is the bottleneck step for many noncoding RNA studies. Results: Here we describe lncFunTK that can run either as standard application or webserver for this purpose. It integrates high throughput sequencing data (i.e. ChIP-seq, CLIP-seq and RNA-seq) to construct the regulatory network associated with lncRNAs. Through the network, it calculates the Functional Information Score (FIS) of each individual lncRNA for prioritizing and inferring its functions through Gene Ontology (GO) terms of neighboring genes. In addition, it also provides utility scripts to support the input data preprocessing and the parameter optimizing. We further demonstrate that lncFunTK can be widely used in various biological systems for lncRNA prioritization and functional annotation. Availability and implementation: The lncFunTK standalone version is an open source package and freely available at http://sunlab.cpy.cuhk.edu.hk/lncfuntk under the MIT license. A webserver implementation is also available at http://sunlab.cpy.cuhk.edu.hk/lncfuntk/runlncfuntk.html. Supplementary information: Supplementary data are available at Bioinformatics online.