iATMEcell: identification of abnormal tumor microenvironment cells to predict the clinical outcomes in cancer based on cell-cell crosstalk network.

Interactions between Tumor microenvironment (TME) cells shape the unique growth environment, sustaining tumor growth and causing the immune escape of tumor cells. Nonetheless, no studies have reported a systematic analysis of cellular interactions in the identification of cancer-related TME cells. Here, we proposed a novel network-based computational method, named as iATMEcell, to identify the abnormal TME cells associated with the biological outcome of interest based on a cell-cell crosstalk network. In the method, iATMEcell first manually collected TME cell types from multiple published studies and obtained their corresponding gene signatures. Then, a weighted cell-cell crosstalk network was constructed in the context of a specific cancer bulk tissue transcriptome data, where the weight between cells reflects both their biological function similarity and the transcriptional dysregulated activities of gene signatures shared by them. Finally, it used a network propagation algorithm to identify significantly dysregulated TME cells. Using the cancer genome atlas (TCGA) Bladder Urothelial Carcinoma training set and two independent validation sets, we illustrated that iATMEcell could identify significant abnormal cells associated with patient survival and immunotherapy response. iATMEcell was further applied to a pan-cancer analysis, which revealed that four common abnormal immune cells play important roles in the patient prognosis across multiple cancer types. Collectively, we demonstrated that iATMEcell could identify potentially abnormal TME cells based on a cell-cell crosstalk network, which provided a new insight into understanding the effect of TME cells in cancer. iATMEcell is developed as an R package, which is freely available on GitHub (https://github.com/hanjunwei-lab/iATMEcell).

Prediction of Transcription Factor Binding Sites on Cell-Free DNA Based on Deep Learning.

Transcription factors (TFs) are important regulatory elements for vital cellular activities, and the identification of transcription factor binding sites (TFBS) can help to explore gene regulatory mechanisms. Research studies have proved that cfDNA (cell-free DNA) shows relatively higher coverage at TFBS due to the protection by TF from degradation by nucleases and short fragments of cfDNA are enriched in TFBS. However, there are still great difficulties in the noninvasive identification of TFBSs from experimental techniques. In this study, we propose a deep learning-based approach that can noninvasively predict TFBSs of cfDNA by learning sequence information from known TFBSs through convolutional neural networks. Under the addition of long short-term memory, our model achieved an area under the curve of 84%. Based on this model to predict cfDNA, we found consistent motifs in cfDNA fragments and lower coverage occurred upstream and downstream of these cfDNA fragments, which is consistent with a previous study. We also found that the binding sites of the same TF differ in different cell lines. TF-specific target genes were detected from cfDNA and were enriched in cancer-related pathways. In summary, our method of locating TFBSs from plasma has the potential to reflect the intrinsic regulatory mechanism from a noninvasive perspective and provide technical guidance for dynamic monitoring of disease in clinical practice.

CNA2Subpathway: identification of dysregulated subpathway driven by copy number alterations in cancer.

Biological pathways reflect the key cellular mechanisms that dictate disease states, drug response and altered cellular function. The local areas of pathways are defined as subpathways (SPs), whose dysfunction has been reported to be associated with the occurrence and development of cancer. With the development of high-throughput sequencing technology, identifying dysfunctional SPs by using multi-omics data has become possible. Moreover, the SPs are not isolated in the biological system but interact with each other. Here, we propose a network-based calculated method, CNA2Subpathway, to identify dysfunctional SPs is driven by somatic copy number alterations (CNAs) in cancer through integrating pathway topology information, multi-omics data and SP crosstalk. This provides a novel way of SP analysis by using the SP interactions in the system biological level. Using data sets from breast cancer and head and neck cancer, we validate the effectiveness of CNA2Subpathway in identifying cancer-relevant SPs driven by the somatic CNAs, which are also shown to be associated with cancer immune and prognosis of patients. We further compare our results with five pathway or SP analysis methods based on CNA and gene expression data without considering SP crosstalk. With these analyses, we show that CNA2Subpathway could help to uncover dysfunctional SPs underlying cancer via the use of SP crosstalk. CNA2Subpathway is developed as an R-based tool, which is freely available on GitHub (https://github.com/hanjunwei-lab/CNA2Subpathway).

A comprehensive characterization of cell-free RNA in spent blastocyst medium and quality prediction for blastocyst.

The rate of pregnancy can be affected by many factors in assisted reproductive technology (ART), and one of which is the quality of embryos. Therefore, selecting the embryos with high potential is crucial for the outcome. Fifteen spent blastocyst medium (SBM) samples were collected from 14 patients who received in vitro fertilization (IVF) or intracytoplasmic sperm injection (ICSI), seven from high-grade embryos and eight from low-grade embryos. Cell-free RNA (cf-RNA) profile of SBM samples were analyzed by RNA sequencing in the present study. It was found that a large amount of cf-RNA were released into SBM, including protein-coding genes (68.9%) and long noncoding RNAs (lncRNAs) (17.26%). Furthermore, a high correlation was observed between blastocyst genes and SBM genes. And the cf-mRNAs of SBM were highly fragmented, and coding sequence (CDS) and untranslated (UTR) regions were released equally. Two hundred and thirty-two differentially expressed genes were identified in high-grade SBM (hSBM) and low-grade SBM (lSBM), which could be potential biomarker in distinguishing the embryos with different quality as an alternative or supplementary approach for subjective morphology criteria. Hence, cf-RNAs sequencing revealed the characterization of circulating transcriptomes of embryos with different quality. Based on the results, the genes related to blastocyst quality were screened, including the genes closely related to translation, immune-signaling pathway, and amino acid metabolism. Overall, the present study showed the types of SBM cf-RNAs, and the integrated analysis of cf-RNAs profiling with morphology grading displayed its potential in predicting blastocyst quality. The present study provided valuable scientific basis for noninvasive embryo selection in ART by RNA-profiling analysis.

Spatial Transcriptomic Analysis Reveals Regional Transcript Changes in Early and Late Stages of rd1 Model Mice with Retinitis Pigmentosa.

Retinitis pigmentosa (RP) is the leading cause of inherited blindness with a genetically heterogeneous disorder. Currently, there is no effective treatment that can protect vision for those with RP. In recent decades, the rd1 mouse has been used to study the pathological mechanisms of RP. Molecular biological studies using rd1 mice have clarified the mechanism of the apoptosis of photoreceptor cells in the early stage of RP. However, the pathological changes in RP over time remain unclear. The unknown pathology mechanism of RP over time and the difficulty of clinical treatment make it urgent to perform more refined and spatially informed molecular biology studies of RP. In this study, spatial transcriptomic analysis is used to study the changes in different retinal layers of rd1 mice at different ages. The results demonstrate the pattern of photoreceptor apoptosis between rd1 mice and the control group. Not only was oxidative stress enhanced in the late stage of RP, but it was accompanied by an up-regulation of the VEGF pathway. Analysis of temporal kinetic trends has further identified patterns of changes in the key pathways of the early and late stages, to help understand the important pathogenesis of RP. Overall, the application of spatial transcriptomics to rd1 mice can help to elucidate the important pathogenesis of RP involving photoreceptor apoptosis and retinal remodeling.

Spatial Transcriptome Profiling of Mouse Hippocampal Single Cell Microzone in Parkinson's Disease.

The hippocampus is an important part of the limbic system in the human brain that has essential roles in spatial navigation and cognitive functions. It is still unknown how gene expression changes in single-cell in different spatial locations of the hippocampus of Parkinson's disease. The purpose of this study was to analyze the gene expression features of single cells in different spatial locations of mouse hippocampus, and to explore the effects of gene expression regulation on learning and memory mechanisms. Here, we obtained 74 single-cell samples from different spatial locations in a mouse hippocampus through microdissection technology, and used single-cell RNA-sequencing and spatial transcriptome sequencing to visualize and quantify the single-cell transcriptome features of tissue sections. The results of differential expression analysis showed that the expression of Sv2b, Neurod6, Grp and Stk32b genes in a hippocampus single cell at different locations was significantly different, and the marker genes of CA1, CA3 and DG subregions were identified. The results of gene function enrichment analysis showed that the up-regulated differentially expressed genes Tubb2a, Eno1, Atp2b1, Plk2, Map4, Pex5l, Fibcd1 and Pdzd2 were mainly involved in neuron to neuron synapse, vesicle-mediated transport in synapse, calcium signaling pathway and neurodegenerative disease pathways, thus affecting learning and memory function. It revealed the transcriptome profile and heterogeneity of spatially located cells in the hippocampus of PD for the first time, and demonstrated that the impaired learning and memory ability of PD was affected by the synergistic effect of CA1 and CA3 subregions neuron genes. These results are crucial for understanding the pathological mechanism of the Parkinson's disease and making precise treatment plans.

Prediction of Time-Series Transcriptomic Gene Expression Based on Long Short-Term Memory with Empirical Mode Decomposition.

RNA degradation can significantly affect the results of gene expression profiling, with subsequent analysis failing to faithfully represent the initial gene expression level. It is urgent to have an artificial intelligence approach to better utilize the limited data to obtain meaningful and reliable analysis results in the case of data with missing destination time. In this study, we propose a method based on the signal decomposition technique and deep learning, named Multi-LSTM. It is divided into two main modules: One decomposes the collected gene expression data by an empirical mode decomposition (EMD) algorithm to obtain a series of sub-modules with different frequencies to improve data stability and reduce modeling complexity. The other is based on long short-term memory (LSTM) as the core predictor, aiming to deeply explore the temporal nonlinear relationships embedded in the sub-modules. Finally, the prediction results of sub-modules are reconstructed to obtain the final prediction results of time-series transcriptomic gene expression. The results show that EMD can efficiently reduce the nonlinearity of the original data, which provides reliable theoretical support to reduce the complexity and improve the robustness of LSTM models. Overall, the decomposition-combination prediction framework can effectively predict gene expression levels at unknown time points.

The Impact of Blood Sample Processing on Ribonucleic Acid (RNA) Sequencing.

In gene quantification and expression analysis, issues with sample selection and processing can be serious, as they can easily introduce irrelevant variables and lead to ambiguous results. This study aims to investigate the extent and mechanism of the impact of sample selection and processing on ribonucleic acid (RNA) sequencing. RNA from PBMCs and blood samples was investigated in this study. The integrity of this RNA was measured under different storage times. All the samples underwent high-throughput sequencing for comprehensive evaluation. The differentially expressed genes and their potential functions were analyzed after the samples were placed at room temperature for 0h, 4h and 8h, and different feature changes in these samples were also revealed. The sequencing results showed that the differences in gene expression were higher with an increased storage time, while the total number of genes detected did not change significantly. There were five genes showing gradient patterns over different storage times, all of which were protein-coding genes that had not been mentioned in previous studies. The effect of different storage times on seemingly the same samples was analyzed in this present study. This research, therefore, provides a theoretical basis for the long-term consideration of whether sample processing should be adequately addressed.

CsPbBr3 perovskite quantum dots-based Janus membrane with multifunction of luminescence, magnetism and aeolotropic electroconductivity.

Lead halide perovskite quantum dots (QDs) are promising semiconductors for next-generation photoelectric devices. However, the development of perovskite QDs-based multifunctional materials still needs to be addressed in order to further advance the application of perovskite QDs. Herein, a successful synthesis of Janus microfibers array Janus membrane (JMAJM) with up-down structure and multifunction of luminescence, magnetism and electroconductivity is firstly achieved based on CsPbBr3 QDs through a parallel electrospinning. JMAJM comprises up-down two layers tightly bonded together. The up-layer of JMAJM is luminescence/magnetism Janus microfibers array (L/M-JMAJM) constructed by [CsPbBr3/polymethyl methacrylate (PMMA)]//[CoFe2O4/PMMA] Janus microfibers as building elements. The down-layer of JMAJM is luminescence/electroconductivity Janus microfibers array (L/E-JMAJM) fabricated by [CsPbBr3/PMMA]//[polyaniline (PANI)/PMMA] Janus microfibers as building elements. Two independent microcosmic regions are designed and realized in a Janus microfiber, confining luminescence with magnetic or conductive substances into their respective regions, thus minimizing adverse effects of other dark-colored functional substances on the fluorescence of CsPbBr3 QDs. This peculiar Janus microfiber enables the effective separation and high integration of CsPbBr3 QDs with other functional substances. The up-down structure of JMAJM ensures a high integration of luminescence, magnetism and conductivity. Meanwhile, JMAJM addresses the environmental instability of CsPbBr3 QDs while simultaneously endows perovskite QDs-based materials with additional functions to realize multifunction. Under ultraviolet excitation, fluorescence characteristics of the CsPbBr3 QDs in JMAJM are maintained, exhibiting a vibrant green emission at 517 nm. Meanwhile, JMAJM achieves a maximum saturation magnetization of 20.32 emu·g-1, high conductance of 10-2 S and aeolotropic electroconductivity degree of 107. The combination of micro-partition with macro-partition in JMAJM receives superior concurrent luminescence-magnetic-conductive multifunction. This work provides a novel idea and strategy for advancing perovskite QDs-based multifunctional materials.

Agarose amplification based sequencing characterization cell-free RNA in preimplantation spent embryo medium.

BACKGROUND: The cell-free RNA (cf-RNA) of spent embryo medium (SEM) has aroused a concern of academic and clinical researchers for its potential use in non-invasive embryo screening. However, comprehensive characterization of cf-RNA from SEM still presents significant technical challenges, primarily due to the limited volume of SEM. Hence, there is urgently need to a small input liquid volume and ultralow amount of cf-RNA library preparation method to unbiased cf-RNA sequencing from SEM. (75) RESULT: Here, we report a high sensitivity agarose amplification-based cf-RNA sequencing method (SEM-Acf) for human preimplantation SEM cf-RNA analysis. It is a cf-RNA sequencing library preparation method by adding agarose amplification. The agarose amplification sensitivity (0.005 pg) and efficiency (105.35 %) were increased than that of without agarose addition (0.45 pg and 96.06 %) by âˆ¼ 90 fold and 9.29 %, respectively. Compared with SMART sequencing (SMART-seq), the correlation of gene expression was stronger in different SEM samples by using SEM-Acf. The cf-RNA number of detected and coverage uniformity of 3' end were significantly increased. The proportion of 5' end adenine, alternative splicing events and short fragments (<400 bp) were increased. It is also found that 4-mer end motifs of cf-RNA fragments was significantly differences between different embryonic stage by day3 spent cleavage medium and day5/6 spent blastocyst medium. (141) SIGNIFICANCE: This study established an efficient SEM amplification and library preparation method. Additionally, we successfully described the characterizations of SEM cf-RNA in preimplantation embryo using SEM-Acf, including expression features and fragment lengths. SEM-Acf facilitates the exploration of cf-RNA as a noninvasive embryo screening biomarker, and opens up potential clinical utilities of small input liquid volume and ultralow amount cf-RNA sequencing. (59).

scRNA-seq data analysis method to improve analysis performance.

With the development of single-cell RNA sequencing technology (scRNA-seq), we have the ability to study biological questions at the level of the individual cell transcriptome. Nowadays, many analysis tools, specifically suitable for single-cell RNA sequencing data, have been developed. In this review, the currently commonly used scRNA-seq protocols are discussed. The upstream processing flow pipeline of scRNA-seq data, including goals and popular tools for reads mapping and expression quantification, quality control, normalization, imputation, and batch effect removal is also introduced. Finally, methods to evaluate these tools in both cellular and genetic dimensions, clustering and differential expression analysis are presented.

Pan-cancer analysis reveals sex-specific signatures in the tumor microenvironment.

The processes of cancer initiation, progression, and response to therapy are affected by the sex of cancer patients. Immunotherapy responses largely depend on the tumor microenvironment (TME), but how sex may shape some TME features, remains unknown. Here, we analyzed immune infiltration signatures across 19 cancer types from 1771 male and 1137 female patients in The Cancer Genome Atlas to evaluate how sex may affect the tumor mutational burden (TMB), immune scores, stromal scores, tumor purity, immune cells, immune checkpoint genes, and functional pathways in the TME. Pan-cancer analyses showed higher TMB and tumor purity scores, as well as lower immune and stromal scores in male patients as compared to female patients. Lung adenocarcinoma, lung squamous carcinoma, kidney papillary carcinoma, and head and neck squamous carcinoma showed the most significant sex biases in terms of infiltrating immune cells, immune checkpoint gene expression, and functional pathways. We further focused on lung adenocarcinoma samples in order to identify and validate sex-specific immune cell biomarkers with prognostic potential. Overall, sex may affect the tumor microenvironment, and sex-specific TME biomarkers may help tailor cancer immunotherapy in certain cancer types.

Extracellular cell-free RNA profile in human large follicles and small follicles.

Background: Previous studies have shown that a large number of valuable and functional cell-free RNAs (cfRNAs) were found in follicular fluid. However, the species and characteristics of follicular fluid cfRNAs have not been reported. Furthermore, their implications are still barely understood in the evaluation of follicular fluid from follicles of different sizes, which warrants further studies. Objective: This study investigated the landscape and characteristics of follicular fluid cfRNAs, the source of organization, and the potential for distinguishing between follicles of different sizes. Methods: Twenty-four follicular fluid samples were collected from 20 patients who received in vitro fertilization (n = 9) or ICSI (n = 11), including 16 large follicular fluid and 8 small follicular fluid samples. Also, the cfRNA profile of follicular fluid samples was analyzed by RNA sequencing. Results: This result indicated that the concentration of follicular fluid cfRNAs ranged from 0.78 to 8.76 ng/ml, and fragment length was 20-200 nucleotides. The concentration and fragment length of large follicular fluid and small follicular fluid samples were not significantly different (p > 0.05). The technical replica correlation of follicular fluid samples ranged from 0.3 to 0.9, and the correlation of small follicular fluid samples was remarkably (p < 0.001) lower than that of large follicular fluid samples. Moreover, this study found that cfRNAs of the follicular fluid could be divided into 37 Ensembl RNA biotypes, and a large number of mRNAs, circRNAs, and lncRNAs were observed in the follicular fluid. The number of cfRNAs in large follicular fluid was remarkably (p < 0.05) higher than that of small follicular fluid. Furthermore, the follicular fluid contained a large amount of intact mRNA and splice junctions and a large number of tissue-derived RNAs, which are at a balanced state of supply and elimination in the follicular fluid. KEGG pathway analysis showed that differentially expressed cfRNAs were enriched in several pathways, including thyroid hormone synthesis, the cGMP-PKG signaling pathway, and inflammatory mediator regulation of TRP channels. In addition, we further showed that four cfRNAs (TK2, AHDC1, PHF21A, and TTYH1) serve as a potential indicator to distinguish the follicles of different sizes. The ROC curve shows great potential to predict follicular fluid from follicles of different sizes [area under the curve (AUC) > 0.88]. Conclusion: Overall, our study revealed that a large number of cfRNAs could be detected in follicular fluid and could serve as a potential non-invasive biomarker in distinguishing between follicles of different sizes. These results may inform the study of the utility and implementation of cfRNAs in clinical practice.

Sensitive and Low-Bias Transcriptome Sequencing Using Agarose PCR.

Transcriptome sequencing has emerged as an important research tool for exploring the mysteries of life at the single-cell level. However, its wide application is limited by the bias associated with the amplification reactions which is essential for library building of trace RNA. In this study, low-melting-point agarose was added to the amplification reactions to take advantage of its molecular crowding effect and polymer cross-linked structure to improve the sensitivity of the reactions and reduce bias. To further evaluate the performance of the method, it was applied to transcriptome sequencing of microregion samples from brain tissue sections of mice with Parkinson's disease at the single cell level. The results showed that agarose PCR had better performance than in-tube PCR. Further application of agarose PCR to transcriptome library sequencing could obtain data closer to that of unamplified. With the addition of low melting point agarose, the sensitivity of the amplification reaction was significantly increased, while homogeneity was increased by approximately 2-fold. Not only that, but this work also provides 11% sensitivity improvement for spatial transcriptomic study on Parkinson's disease-associated gene detection. The agarose PCR provides a new tool for efficient and homogeneous amplification of trace samples and can be widely used for spatial transcriptome library sequencing and studies.

Identification of Cancer Dysfunctional Subpathways by Integrating DNA Methylation, Copy Number Variation, and Gene-Expression Data.

A subpathway is defined as the local region of a biological pathway with specific biological functions. With the generation of large-scale sequencing data, there are more opportunities to study the molecular mechanisms of cancer development. It is necessary to investigate the potential impact of DNA methylation, copy number variation (CNV), and gene-expression changes in the molecular states of oncogenic dysfunctional subpathways. We propose a novel method, Identification of Cancer Dysfunctional Subpathways (ICDS), by integrating multi-omics data and pathway topological information to identify dysfunctional subpathways. We first calculated gene-risk scores by integrating the three following types of data: DNA methylation, CNV, and gene expression. Second, we performed a greedy search algorithm to identify the key dysfunctional subpathways within pathways for which the discriminative scores were locally maximal. Finally, a permutation test was used to calculate the statistical significance level for these key dysfunctional subpathways. We validated the effectiveness of ICDS in identifying dysregulated subpathways using datasets from liver hepatocellular carcinoma (LIHC), head-neck squamous cell carcinoma (HNSC), cervical squamous cell carcinoma, and endocervical adenocarcinoma. We further compared ICDS with methods that performed the same subpathway identification algorithm but only considered DNA methylation, CNV, or gene expression (defined as ICDS_M, ICDS_CNV, or ICDS_G, respectively). With these analyses, we confirmed that ICDS better identified cancer-associated subpathways than the three other methods, which only considered one type of data. Our ICDS method has been implemented as a freely available R-based tool (https://cran.r-project.org/web/packages/ICDS).

Prioritization of candidate cancer drugs based on a drug functional similarity network constructed by integrating pathway activities and drug activities.

Due to the speed, efficiency, relative risk, and lower costs compared to traditional drug discovery, the prioritization of candidate drugs for repurposing against cancers of interest has attracted the attention of experts in recent years. Herein, we present a powerful computational approach, termed prioritization of candidate drugs (PriorCD), for the prioritization of candidate cancer drugs based on a global network propagation algorithm and a drug-drug functional similarity network constructed by integrating pathway activity profiles and drug activity profiles. This provides a new approach to drug repurposing by first considering the drug functional similarities at the pathway level. The performance of PriorCD in drug repurposing was evaluated by using drug datasets of breast cancer and ovarian cancer. Cross-validation tests on the drugs approved for the treatment of these cancers indicated that our approach can achieve area under receiver-operating characteristic curve (AUROC) values greater than 0.82. Furthermore, literature searches validated our results, and comparison with other classical gene-based repurposing methods indicated that our pathway-level PriorCD is comparatively more effective at prioritizing candidate drugs with similar therapeutic effects. We hope that our study will be of benefit to the field of drug discovery. In order to expand the usage of PriorCD, a freely available R-based package, PriorCD, has been developed to prioritize candidate anticancer drugs for drug repurposing.