Emx2 underlies the development and evolution of marsupial gliding membranes.

Phenotypic variation among species is a product of evolutionary changes to developmental programs1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium-an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.

Signalling by senescent melanocytes hyperactivates hair growth.

Niche signals maintain stem cells in a prolonged quiescence or transiently activate them for proper regeneration1. Altering balanced niche signalling can lead to regenerative disorders. Melanocytic skin nevi in human often display excessive hair growth, suggesting hair stem cell hyperactivity. Here, using genetic mouse models of nevi2,3, we show that dermal clusters of senescent melanocytes drive epithelial hair stem cells to exit quiescence and change their transcriptome and composition, potently enhancing hair renewal. Nevus melanocytes activate a distinct secretome, enriched for signalling factors. Osteopontin, the leading nevus signalling factor, is both necessary and sufficient to induce hair growth. Injection of osteopontin or its genetic overexpression is sufficient to induce robust hair growth in mice, whereas germline and conditional deletions of either osteopontin or CD44, its cognate receptor on epithelial hair cells, rescue enhanced hair growth induced by dermal nevus melanocytes. Osteopontin is overexpressed in human hairy nevi, and it stimulates new growth of human hair follicles. Although broad accumulation of senescent cells, such as upon ageing or genotoxic stress, is detrimental for the regenerative capacity of tissue4, we show that signalling by senescent cell clusters can potently enhance the activity of adjacent intact stem cells and stimulate tissue renewal. This finding identifies senescent cells and their secretome as an attractive therapeutic target in regenerative disorders.

Spatial transition tensor of single cells.

Spatial transcriptomics and messenger RNA splicing encode extensive spatiotemporal information for cell states and transitions. The current lineage-inference methods either lack spatial dynamics for state transition or cannot capture different dynamics associated with multiple cell states and transition paths. Here we present spatial transition tensor (STT), a method that uses messenger RNA splicing and spatial transcriptomes through a multiscale dynamical model to characterize multistability in space. By learning a four-dimensional transition tensor and spatial-constrained random walk, STT reconstructs cell-state-specific dynamics and spatial state transitions via both short-time local tensor streamlines between cells and long-time transition paths among attractors. Benchmarking and applications of STT on several transcriptome datasets via multiple technologies on epithelial-mesenchymal transitions, blood development, spatially resolved mouse brain and chicken heart development, indicate STT's capability in recovering cell-state-specific dynamics and their associated genes not seen using existing methods. Overall, STT provides a consistent multiscale description of single-cell transcriptome data across multiple spatiotemporal scales.

Screening cell-cell communication in spatial transcriptomics via collective optimal transport.

Spatial transcriptomic technologies and spatially annotated single-cell RNA sequencing datasets provide unprecedented opportunities to dissect cell-cell communication (CCC). However, incorporation of the spatial information and complex biochemical processes required in the reconstruction of CCC remains a major challenge. Here, we present COMMOT (COMMunication analysis by Optimal Transport) to infer CCC in spatial transcriptomics, which accounts for the competition between different ligand and receptor species as well as spatial distances between cells. A collective optimal transport method is developed to handle complex molecular interactions and spatial constraints. Furthermore, we introduce downstream analysis tools to infer spatial signaling directionality and genes regulated by signaling using machine learning models. We apply COMMOT to simulation data and eight spatial datasets acquired with five different technologies to show its effectiveness and robustness in identifying spatial CCC in data with varying spatial resolutions and gene coverages. Finally, COMMOT identifies new CCCs during skin morphogenesis in a case study of human epidermal development.

exFINDER: identify external communication signals using single-cell transcriptomics data.

Cells make decisions through their communication with other cells and receiving signals from their environment. Using single-cell transcriptomics, computational tools have been developed to infer cell-cell communication through ligands and receptors. However, the existing methods only deal with signals sent by the measured cells in the data, the received signals from the external system are missing in the inference. Here, we present exFINDER, a method that identifies such external signals received by the cells in the single-cell transcriptomics datasets by utilizing the prior knowledge of signaling pathways. In particular, exFINDER can uncover external signals that activate the given target genes, infer the external signal-target signaling network (exSigNet), and perform quantitative analysis on exSigNets. The applications of exFINDER to scRNA-seq datasets from different species demonstrate the accuracy and robustness of identifying external signals, revealing critical transition-related signaling activities, inferring critical external signals and targets, clustering signal-target paths, and evaluating relevant biological events. Overall, exFINDER can be applied to scRNA-seq data to reveal the external signal-associated activities and maybe novel cells that send such signals.

The promising application of cell-cell interaction analysis in cancer from single-cell and spatial transcriptomics.

Cell-cell interactions instruct cell fate and function. These interactions are hijacked to promote cancer development. Single-cell transcriptomics and spatial transcriptomics have become powerful new tools for researchers to profile the transcriptional landscape of cancer at unparalleled genetic depth. In this review, we discuss the rapidly growing array of computational tools to infer cell-cell interactions from non-spatial single-cell RNA-sequencing and the limited but growing number of methods for spatial transcriptomics data. Downstream analyses of these computational tools and applications to cancer studies are highlighted. We finish by suggesting several directions for further extensions that anticipate the increasing availability of multi-omics cancer data.

Tipping points in epithelial-mesenchymal lineages from single-cell transcriptomics data.

Understanding cell fate decision-making during complex biological processes is an open challenge that is now aided by high-resolution single-cell sequencing technologies. Specifically, it remains challenging to identify and characterize transition states corresponding to "tipping points" whereby cells commit to new cell states. Here, we present a computational method that takes advantage of single-cell transcriptomics data to infer the stability and gene regulatory networks (GRNs) along cell lineages. Our method uses the unspliced and spliced counts from single-cell RNA sequencing data and cell ordering along lineage trajectories to train an RNA splicing multivariate model, from which cell-state stability along the lineage is inferred based on spectral analysis of the model's Jacobian matrix. Moreover, the model infers the RNA cross-species interactions resulting in GRNs and their variation along the cell lineage. When applied to epithelial-mesenchymal transition in ovarian and lung cancer-derived cell lines, our model predicts a saddle-node transition between the epithelial and mesenchymal states passing through an unstable, intermediate cell state. Furthermore, we show that the underlying GRN controlling epithelial-mesenchymal transition rearranges during the transition, resulting in denser and less modular networks in the intermediate state. Overall, our method represents a flexible tool to study cell lineages with a combination of theory-driven modeling and single-cell transcriptomics data.

DURIAN: an integrative deconvolution and imputation method for robust signaling analysis of single-cell transcriptomics data.

Single-cell RNA sequencing trades read-depth for dimensionality, often leading to loss of critical signaling gene information that is typically present in bulk data sets. We introduce DURIAN (Deconvolution and mUltitask-Regression-based ImputAtioN), an integrative method for recovery of gene expression in single-cell data. Through systematic benchmarking, we demonstrate the accuracy, robustness and empirical convergence of DURIAN using both synthetic and published data sets. We show that use of DURIAN improves single-cell clustering, low-dimensional embedding, and recovery of intercellular signaling networks. Our study resolves several inconsistent results of cell-cell communication analysis using single-cell or bulk data independently. The method has broad application in biomarker discovery and cell signaling analysis using single-cell transcriptomics data sets.

UFold: fast and accurate RNA secondary structure prediction with deep learning.

For many RNA molecules, the secondary structure is essential for the correct function of the RNA. Predicting RNA secondary structure from nucleotide sequences is a long-standing problem in genomics, but the prediction performance has reached a plateau over time. Traditional RNA secondary structure prediction algorithms are primarily based on thermodynamic models through free energy minimization, which imposes strong prior assumptions and is slow to run. Here, we propose a deep learning-based method, called UFold, for RNA secondary structure prediction, trained directly on annotated data and base-pairing rules. UFold proposes a novel image-like representation of RNA sequences, which can be efficiently processed by Fully Convolutional Networks (FCNs). We benchmark the performance of UFold on both within- and cross-family RNA datasets. It significantly outperforms previous methods on within-family datasets, while achieving a similar performance as the traditional methods when trained and tested on distinct RNA families. UFold is also able to predict pseudoknots accurately. Its prediction is fast with an inference time of about 160 ms per sequence up to 1500 bp in length. An online web server running UFold is available at https://ufold.ics.uci.edu. Code is available at https://github.com/uci-cbcl/UFold.

Integrated analysis of multimodal single-cell data with structural similarity.

Multimodal single-cell sequencing technologies provide unprecedented information on cellular heterogeneity from multiple layers of genomic readouts. However, joint analysis of two modalities without properly handling the noise often leads to overfitting of one modality by the other and worse clustering results than vanilla single-modality analysis. How to efficiently utilize the extra information from single cell multi-omics to delineate cell states and identify meaningful signal remains as a significant computational challenge. In this work, we propose a deep learning framework, named SAILERX, for efficient, robust, and flexible analysis of multi-modal single-cell data. SAILERX consists of a variational autoencoder with invariant representation learning to correct technical noises from sequencing process, and a multimodal data alignment mechanism to integrate information from different modalities. Instead of performing hard alignment by projecting both modalities to a shared latent space, SAILERX encourages the local structures of two modalities measured by pairwise similarities to be similar. This strategy is more robust against overfitting of noises, which facilitates various downstream analysis such as clustering, imputation, and marker gene detection. Furthermore, the invariant representation learning part enables SAILERX to perform integrative analysis on both multi- and single-modal datasets, making it an applicable and scalable tool for more general scenarios.

High dietary wheat starch negatively regulated growth performance, glucose and lipid metabolisms, liver and intestinal health of juvenile largemouth bass, Micropterus salmoides.

Largemouth bass (Micropterus salmoides) were fed with three diets containing 6%, 12%, and 18% wheat starch for 70 days to examine their impacts on growth performance, glucose and lipid metabolisms, and liver and intestinal health. The results suggested that the 18% starch group inhibited the growth, and improved the hepatic glycogen content compared with the 6% and 12% starch groups (P < 0.05). High starch significantly improved the activities of glycolysis-related enzymes, hexokinase (HK), glucokinase (GK), phosphofructokinase (PFK), and pyruvate kinase (PK) (P < 0.05); promoted the mRNA expression of glycolysis-related phosphofructokinase (pfk); decreased the activities of gluconeogenesis-related enzymes, pyruvate carboxylase (PC), and phosphoenolpyruvate carboxykinase (PEPCK); and reduced the mRNA expression of gluconeogenesis-related fructose-1,6-bisphosphatase-1(fbp1) (P < 0.05). High starch reduced the hepatic mRNA expressions of bile acid metabolism-related cholesterol hydroxylase (cyp7a1) and small heterodimer partner (shp) (P < 0.05), increased the activity of hepatic fatty acid synthase (FAS) (P < 0.05), and reduced the hepatic mRNA expressions of lipid metabolism-related peroxisome proliferator-activated receptor α (ppar-α) and carnitine palmitoyltransferase 1α (cpt-1α) (P < 0.05). High starch promoted inflammation; significantly reduced the mRNA expressions of anti-inflammatory cytokines transforming growth factor-ß1 (tgf-ß1), interleukin-10 (il-10), and interleukin-11ß (il-11ß); and increased the mRNA expressions of pro-inflammatory cytokine tumor necrosis factor-α (tnf-α), interleukin-1ß (il-1ß), and interleukin-8 (il-8) in the liver and intestinal tract (P < 0.05). Additionally, high starch negatively influenced the intestinal microbiota, with the reduced relative abundance of Trichotes and Actinobacteria and the increased relative abundance of Firmicutes and Proteobacteria. In conclusion, low dietary wheat starch level (6%) was more profitable to the growth and health of M. salmoides, while high dietary starch level (12% and 18%) could regulate the glucose and lipid metabolisms, impair the liver and intestinal health, and thus decrease the growth performance of M. salmoides.

Theoretical and computational tools to model multistable gene regulatory networks.

The last decade has witnessed a surge of theoretical and computational models to describe the dynamics of complex gene regulatory networks, and how these interactions can give rise to multistable and heterogeneous cell populations. As the use of theoretical modeling to describe genetic and biochemical circuits becomes more widespread, theoreticians with mathematical and physical backgrounds routinely apply concepts from statistical physics, non-linear dynamics, and network theory to biological systems. This review aims at providing a clear overview of the most important methodologies applied in the field while highlighting current and future challenges. It also includes hands-on tutorials to solve and simulate some of the archetypical biological system models used in the field. Furthermore, we provide concrete examples from the existing literature for theoreticians that wish to explore this fast-developing field. Whenever possible, we highlight the similarities and differences between biochemical and regulatory networks and 'classical' systems typically studied in non-equilibrium statistical and quantum mechanics.

Identifying age-specific gene signatures of the human cerebral cortex with joint analysis of transcriptomes and functional connectomes.

The human cerebral cortex undergoes profound structural and functional dynamic variations across the lifespan, whereas the underlying molecular mechanisms remain unclear. Here, with a novel method transcriptome-connectome correlation analysis (TCA), which integrates the brain functional magnetic resonance images and region-specific transcriptomes, we identify age-specific cortex (ASC) gene signatures for adolescence, early adulthood and late adulthood. The ASC gene signatures are significantly correlated with the cortical thickness (P-value <2.00e-3) and myelination (P-value <1.00e-3), two key brain structural features that vary in accordance with brain development. In addition to the molecular underpinning of age-related brain functions, the ASC gene signatures allow delineation of the molecular mechanisms of neuropsychiatric disorders, such as the regulation between ARNT2 and its target gene ETF1 involved in Schizophrenia. We further validate the ASC gene signatures with published gene sets associated with the adult cortex, and confirm the robustness of TCA on other brain image datasets. Availability: All scripts are written in R. Scripts for the TCA method and related statistics result can be freely accessed at https://github.com/Soulnature/TCA. Additional data related to this paper may be requested from the authors.

spliceJAC: transition genes and state-specific gene regulation from single-cell transcriptome data.

Extracting dynamical information from single-cell transcriptomics is a novel task with the promise to advance our understanding of cell state transition and interactions between genes. Yet, theory-oriented, bottom-up approaches that consider differences among cell states are largely lacking. Here, we present spliceJAC, a method to quantify the multivariate mRNA splicing from single-cell RNA sequencing (scRNA-seq). spliceJAC utilizes the unspliced and spliced mRNA count matrices to constructs cell state-specific gene-gene regulatory interactions and applies stability analysis to predict putative driver genes critical to the transitions between cell states. By applying spliceJAC to biological systems including pancreas endothelium development and epithelial-mesenchymal transition (EMT) in A549 lung cancer cells, we predict genes that serve specific signaling roles in different cell states, recover important differentially expressed genes in agreement with pre-existing analysis, and predict new transition genes that are either exclusive or shared between different cell state transitions.

Circadian control of interferon-sensitive gene expression in murine skin.

The circadian clock coordinates a variety of immune responses with signals from the external environment to promote survival. We investigated the potential reciprocal relationship between the circadian clock and skin inflammation. We treated mice topically with the Toll-like receptor 7 (TLR7) agonist imiquimod (IMQ) to activate IFN-sensitive gene (ISG) pathways and induce psoriasiform inflammation. IMQ transiently altered core clock gene expression, an effect mirrored in human patient psoriatic lesions. In mouse skin 1 d after IMQ treatment, ISGs, including the key ISG transcription factor IFN regulatory factor 7 (Irf7), were more highly induced after treatment during the day than the night. Nuclear localization of phosphorylated-IRF7 was most prominently time-of-day dependent in epidermal leukocytes, suggesting that these cell types play an important role in the diurnal ISG response to IMQ. Mice lacking Bmal1 systemically had exacerbated and arrhythmic ISG/Irf7 expression after IMQ. Furthermore, daytime-restricted feeding, which affects the phase of the skin circadian clock, reverses the diurnal rhythm of IMQ-induced ISG expression in the skin. These results suggest a role for the circadian clock, driven by BMAL1, as a negative regulator of the ISG response, and highlight the finding that feeding time can modulate the skin immune response. Since the IFN response is essential for the antiviral and antitumor effects of TLR activation, these findings are consistent with the time-of-day-dependent variability in the ability to fight microbial pathogens and tumor initiation and offer support for the use of chronotherapy for their treatment.

SAILER: scalable and accurate invariant representation learning for single-cell ATAC-seq processing and integration.

MOTIVATION: Single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) provides new opportunities to dissect epigenomic heterogeneity and elucidate transcriptional regulatory mechanisms. However, computational modeling of scATAC-seq data is challenging due to its high dimension, extreme sparsity, complex dependencies and high sensitivity to confounding factors from various sources. RESULTS: Here, we propose a new deep generative model framework, named SAILER, for analyzing scATAC-seq data. SAILER aims to learn a low-dimensional nonlinear latent representation of each cell that defines its intrinsic chromatin state, invariant to extrinsic confounding factors like read depth and batch effects. SAILER adopts the conventional encoder-decoder framework to learn the latent representation but imposes additional constraints to ensure the independence of the learned representations from the confounding factors. Experimental results on both simulated and real scATAC-seq datasets demonstrate that SAILER learns better and biologically more meaningful representations of cells than other methods. Its noise-free cell embeddings bring in significant benefits in downstream analyses: clustering and imputation based on SAILER result in 6.9% and 18.5% improvements over existing methods, respectively. Moreover, because no matrix factorization is involved, SAILER can easily scale to process millions of cells. We implemented SAILER into a software package, freely available to all for large-scale scATAC-seq data analysis. AVAILABILITY AND IMPLEMENTATION: The software is publicly available at https://github.com/uci-cbcl/SAILER. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Induction and application of ferroptosis in cancer therapy.

At present, more than one cell death pathways have been found, one of which is ferroptosis. Ferroptosis was discovered in 2012 and described as an iron-dependent and lipid peroxidation-driven regulated cell death pathway. In the past few years, ferroptosis has been shown to induce tumor cell death, providing new ideas for tumor treatment. In this article, we summarize the latest advances in ferroptosis-induced tumor therapy at the intersection of tumor biology, molecular biology, redox biology, and materials chemistry. First, we state the characteristics of ferroptosis in cells, then introduce the key molecular mechanism of ferroptosis, and describes the relationship between ferroptosis and oxidative stress signaling pathways. Finally, we focused on several types of ferroptosis inducers discovered by scholars, and the application of ferroptosis in systemic chemotherapy, radiotherapy, immunotherapy and nanomedicine, in the hope that ferroptosis can exert its potential in the treatment of tumors.

Improved urine DNA methylation panel for early bladder cancer detection.

BACKGROUND: Bladder cancer is one of the most common malignancies but the corresponding diagnostic methods are either invasive or limited in specificity and/or sensitivity. This study aimed to develop a urine-based methylation panel for bladder cancer detection by improving published panels and validate performance of the new panel with clinical samples. METHODS: Related researches were reviewed and 19 potential panels were selected. RRBS was performed on a cohort with 45 samples to reassess these panels and a new panel inherited best markers was developed. The new panel was applied with qMSP platform to 33 samples from the RRBS cohort and the results were compared to those of RRBS. Lastly, another larger cohort with 207 samples was used to validate new panel performance with qMSP. RESULTS: Three biomarkers (PCDH17, POU4F2 and PENK) were selected to construct a new panel P3. P3 panel achieved 100% specificity and 71% sensitivity with RRBS in corresponding cohort and then showed a better performance of 100% specificity and 84% sensitivity with qMSP platforms in a balanced cohort. When validated with 207-sample cohort, P3 with qMSP showed a performance of 97% specificity and 87% sensitivity which was modestly improved compared to the panels it derided from. CONCLUSIONS: Overall, the P3 panel achieved relatively high sensitivity and accuracy in bladder cancer detection.

Urinary metabolomics for discovering metabolic biomarkers of bladder cancer by UPLC-MS.

Bladder cancer (BC) is one of the most frequent cancer in the world, and its incidence is rising worldwide, especially in developed countries. Urine metabolomics is a powerful approach to discover potential biomarkers for cancer diagnosis. In this study, we applied an ultra-performance liquid chromatography coupled to mass spectrometry (UPLC-MS) method to profile the metabolites in urine from 29 bladder cancer patients and 15 healthy controls. The differential metabolites were extracted and analyzed by univariate and multivariate analysis methods. Together, 19 metabolites were discovered as differently expressed biomarkers in the two groups, which mainly related to the pathways of phenylacetate metabolism, propanoate metabolism, fatty acid metabolism, pyruvate metabolism, arginine and proline metabolism, glycine and serine metabolism, and bile acid biosynthesis. In addition, a subset of 11 metabolites of those 19 ones were further filtered as potential biomarkers for BC diagnosis by using logistic regression model. The results revealed that the area under the curve (AUC) value, sensitivity and specificity of receiving operator characteristic (ROC) curve were 0.983, 95.3% and 100%, respectively, indicating an excellent discrimination power for BC patients from healthy controls. It was the first time to reveal the potential diagnostic markers of BC by metabolomics, and this will provide a new sight for exploring the biomarkers of the other disease in the future work.

Inferring latent temporal progression and regulatory networks from cross-sectional transcriptomic data of cancer samples.

Unraveling molecular regulatory networks underlying disease progression is critically important for understanding disease mechanisms and identifying drug targets. The existing methods for inferring gene regulatory networks (GRNs) rely mainly on time-course gene expression data. However, most available omics data from cross-sectional studies of cancer patients often lack sufficient temporal information, leading to a key challenge for GRN inference. Through quantifying the latent progression using random walks-based manifold distance, we propose a latent-temporal progression-based Bayesian method, PROB, for inferring GRNs from the cross-sectional transcriptomic data of tumor samples. The robustness of PROB to the measurement variabilities in the data is mathematically proved and numerically verified. Performance evaluation on real data indicates that PROB outperforms other methods in both pseudotime inference and GRN inference. Applications to bladder cancer and breast cancer demonstrate that our method is effective to identify key regulators of cancer progression or drug targets. The identified ACSS1 is experimentally validated to promote epithelial-to-mesenchymal transition of bladder cancer cells, and the predicted FOXM1-targets interactions are verified and are predictive of relapse in breast cancer. Our study suggests new effective ways to clinical transcriptomic data modeling for characterizing cancer progression and facilitates the translation of regulatory network-based approaches into precision medicine.