Position-dependent motif characterization using non-negative matrix factorization

Efficient discovery of frequently co-occurring mutations in a sequence database with matrix factorization

We have developed a robust method for efficiently tracking multiple co-occurring mutations in a sequence database. Evolution often hinges on the interaction of several mutations to produce significant phenotypic changes that lead to the proliferation of a variant. However, identifying numerous simultaneous mutations across a vast database of sequences poses a significant computational challenge. Our approach leverages a matrix factorization technique to automatically and efficiently pinpoint subsets of positions where co-mutations occur, appearing in a substantial number of sequences within the database. We validated our method using SARS-CoV-2 receptor-binding domains, comprising approximately seven hundred thousand sequences of the Spike protein, demonstrating superior performance compared to a reasonably exhaustive brute-force method. Furthermore, we explore the biological significance of the identified co-mutational positions (CMPs) and their potential impact on the virus's evolution and functionality, identifying key mutations in Delta and Omicron variants. This analysis underscores the significant role of identified CMPs in understanding the evolutionary trajectory. By tracking the "birth" and "death" of CMPs, we can elucidate the persistence and impact of specific groups of mutations across different viral strains, providing valuable insights into the virus' adaptability and thus, possibly aiding vaccine design strategies.

A tumor-infiltrating B lymphocytes -related index based on machine-learning predicts prognosis and immunotherapy response in lung adenocarcinoma

Introduction

Tumor-infiltrating B lymphocytes (TILBs) play a pivotal role in shaping the immune microenvironment of tumors (TIME) and in the progression of lung adenocarcinoma (LUAD). However, there remains a scarcity of research that has thoroughly and systematically delineated the characteristics of TILBs in LUAD.

Method

The research employed single-cell RNA sequencing from the GSE117570 dataset to identify markers linked to TILBs. A comprehensive machine learning approach, utilizing ten distinct algorithms, facilitated the creation of a TILB-related index (BRI) across the TCGA, GSE31210, and GSE72094 datasets. We used multiple algorithms to evaluate the relationships between BRI and TIME, as well as immune therapy-related biomarkers. Additionally, we assessed the role of BRI in predicting immune therapy response in two datasets, GSE91061 and GSE126044.

Result

BRI functioned as an independent risk determinant in LUAD, demonstrating a robust and reliable capacity to predict overall survival rates. We observed significant differences in the scores of B cells, M2 macrophages, NK cells, and regulatory T cells between the high and low BRI score groups. Notably, BRI was found to inversely correlate with cytotoxic CD8+ T-cell infiltration (r = -0.43, p < 0.001) and positively correlate with regulatory T cells (r = 0.31, p = 0.008). We also found that patients with lower BRI were more likely to respond to immunotherapy and were associated with reduced IC50 values for standard chemotherapy and targeted therapy drugs, in contrast to higher BRI. Additionally, the BRI-based survival prediction nomogram demonstrated significant promise for clinical application in predicting the 1-, 3-, and 5-year overall survival rates among LUAD patients.

Discussion

Our study developed a BRI model using ten different algorithms and 101 algorithm combinations. The BRI could be a valuable tool for risk stratification, prognosis, and selection of treatment approaches.

Altered dNTP pools accelerate tumor formation in mice

Alterations in deoxyribonucleoside triphosphate (dNTP) pools have been linked to increased mutation rates and genome instability in unicellular organisms and cell cultures. However, the role of dNTP pool changes in tumor development in mammals remains unclear. In this study, we present a mouse model with a point mutation at the allosteric specificity site of ribonucleotide reductase, RRM1-Y285A. This mutation reduced ribonucleotide reductase activity, impairing the synthesis of deoxyadenosine triphosphate (dATP) and deoxyguanosine triphosphate (dGTP). Heterozygous Rrm1+/Y285A mice exhibited distinct alterations in dNTP pools across various organs, shorter lifespans and earlier tumor onset compared with wild-type controls. Mutational spectrum analysis of tumors revealed two distinct signatures, one resembling a signature extracted from a human cancer harboring a mutation of the same amino acid residue in ribonucleotide reductase, RRM1Y285C. Our findings suggest that mutations in enzymes involved in dNTP metabolism can serve as drivers of cancer development.

The genetics of spatiotemporal variation in cortical thickness in youth

Prior studies have shown strong genetic effects on cortical thickness (CT), structural covariance, and neurodevelopmental trajectories in childhood and adolescence. However, the importance of genetic factors on the induction of spatiotemporal variation during neurodevelopment remains poorly understood. Here, we explore the genetics of maturational coupling by examining 308 MRI-derived regional CT measures in a longitudinal sample of 677 twins and family members. We find dynamic inter-regional genetic covariation in youth, with the emergence of regional subnetworks in late childhood and early adolescence. Three critical neurodevelopmental epochs in genetically-mediated maturational coupling were identified, with dramatic network strengthening near eleven years of age. These changes are associated with statistically-significant (empirical p-value <0.0001) increases in network strength as measured by average clustering coefficient and assortativity. We then identify genes from the Allen Human Brain Atlas with similar co-expression patterns to genetically-mediated structural covariation in children. This set was enriched for genes involved in potassium transport and dendrite formation. Genetically-mediated CT-CT covariance was also strongly correlated with expression patterns for genes located in cells of neuronal origin.

Mutations in yeast Pcf11, a conserved protein essential for mRNA 3' end processing and transcription termination, elicit the Environmental Stress Response

The Pcf11 protein is an essential subunit of the large complex that cleaves and polyadenylates eukaryotic mRNA precursor. It has also been functionally linked to gene-looping, termination of RNA Polymerase II (Pol II) transcripts, and mRNA export. We have examined a poorly characterized but conserved domain (amino acids 142-225) of the Saccharomyces cerevisiae  Pcf11 and found that while it is not needed for mRNA 3' end processing or termination downstream of the poly(A) sites of protein-coding genes, its presence improves the interaction with Pol II and the use of transcription terminators near gene promoters. Analysis of genome-wide Pol II occupancy in cells with Pcf11 missing this region, as well as Pcf11 mutated in the Pol II CTD Interacting Domain, indicates that systematic changes in mRNA expression are mediated primarily at the level of transcription. Global expression analysis also shows that a general stress response, involving both activation and suppression of specific gene sets known to be regulated in response to a wide variety of stresses, is induced in the two pcf11 mutants, even though cells are grown in optimal conditions. The mutants also cause an unbalanced expression of cell wall-related genes that does not activate the Cell Wall Integrity pathway but is associated with strong caffeine sensitivity. Based on these findings, we propose that Pcf11 can modulate the expression level of specific functional groups of genes in ways that do not involve its well-characterized role in mRNA 3' end processing.

Sctensor detects many-to-many cell-cell interactions from single cell RNA-sequencing data

BACKGROUND

Complex biological systems are described as a multitude of cell-cell interactions (CCIs). Recent single-cell RNA-sequencing studies focus on CCIs based on ligand-receptor (L-R) gene co-expression but the analytical methods are not appropriate to detect many-to-many CCIs.

RESULTS

In this work, we propose scTensor, a novel method for extracting representative triadic relationships (or hypergraphs), which include ligand-expression, receptor-expression, and related L-R pairs.

CONCLUSIONS

Through extensive studies with simulated and empirical datasets, we have shown that scTensor can detect some hypergraphs that cannot be detected using conventional CCI detection methods, especially when they include many-to-many relationships. scTensor is implemented as a freely available R/Bioconductor package.

Identifying promoter sequence architectures via a chunking-based algorithm using non-negative matrix factorisation

Core promoters are stretches of DNA at the beginning of genes that contain information that facilitates the binding of transcription initiation complexes. Different functional subsets of genes have core promoters with distinct architectures and characteristic motifs. Some of these motifs inform the selection of transcription start sites (TSS). By discovering motifs with fixed distances from known TSS positions, we could in principle classify promoters into different functional groups. Due to the variability and overlap of architectures, promoter classification is a difficult task that requires new approaches. In this study, we present a new method based on non-negative matrix factorisation (NMF) and the associated software called seqArchR that clusters promoter sequences based on their motifs at near-fixed distances from a reference point, such as TSS. When combined with experimental data from CAGE, seqArchR can efficiently identify TSS-directing motifs, including known ones like TATA, DPE, and nucleosome positioning signal, as well as novel lineage-specific motifs and the function of genes associated with them. By using seqArchR on developmental time courses, we reveal how relative use of promoter architectures changes over time with stage-specific expression. seqArchR is a powerful tool for initial genome-wide classification and functional characterisation of promoters. Its use cases are more general: it can also be used to discover any motifs at near-fixed distances from a reference point, even if they are present in only a small subset of sequences.

Machine Learning-Based Integration Develops a Macrophage-Related Index for Predicting Prognosis and Immunotherapy Response in Lung Adenocarcinoma

BACKGROUND

Macrophages play a critical role in tumor immune microenvironment (TIME) formation and cancer progression in lung adenocarcinoma (LUAD). However, few studies have comprehensively and systematically described the characteristics of macrophages in LUAD.

METHODS

This study identified macrophage-related markers with single-cell RNA sequencing data from the GSE189487 dataset. An integrative machine learning-based procedure based on 10 algorithms was developed to construct a macrophage-related index (MRI) in The Cancer Genome Atlas (TCGA), GSE30219, GSE31210, and GSE72094 datasets. Several algorithms were used to evaluate the associations of MRI with TIME and immunotherapy-related biomarkers. The role of MRI in predicting the immunotherapy response was evaluated with the GSE91061 dataset.

RESULTS

The optimal MRI constructed by the combination of the Lasso algorithm and plsRCox was an independent risk factor in LUAD and showed a stable and powerful performance in predicting the overall survival rate of patients with LUAD. Those with low MRI scores had a higher TIME score, a higher level of immune cells, a higher immunophenoscore, and a lower Tumor Immune Dysfunction and Exclusion (TIDE) score, indicating a better response to immunotherapy. The IC50 value of common drugs for chemotherapy and target therapy with low MRI scores was higher compared to high MRI scores. Moreover, the survival prediction nomogram, developed from MRI, had good potential for clinical application in predicting the 1-, 3-, and 5-year overall survival rate of LUAD.

CONCLUSION

Our study constructed for the first time a consensus MRI for LUAD with 10 machine learning algorithms. The MRI could be helpful for risk stratification, prognosis, and selection of treatment approach in LUAD.

Gene-regulation modules in nonalcoholic fatty liver disease revealed by single-nucleus ATAC-seq

We investigated the progression of nonalcoholic fatty liver disease from fatty liver to steatohepatitis using single-nucleus and bulk ATAC-seq on the livers of rats fed a high-fat diet (HFD). Rats fed HFD for 4 wk developed fatty liver, and those fed HFD for 8 wk further progressed to steatohepatitis. We observed an increase in the proportion of inflammatory macrophages, consistent with the pathological progression. Utilizing machine learning, we divided global gene regulation into modules, wherein transcription factors within a module could regulate genes within the same module, reaffirming known regulatory relationships between transcription factors and biological processes. We identified core genes-central to co-expression and protein-protein interaction-for the biological processes discovered. Notably, a large part of the core genes overlapped with genes previously implicated in nonalcoholic fatty liver disease. Single-nucleus ATAC-seq, combined with data-driven statistical analysis, offers insight into in vivo global gene regulation as a combination of modules and assists in identifying core genes of relevant biological processes.

DeepTPpred: A Deep Learning Approach With Matrix Factorization for Predicting Therapeutic Peptides by Integrating Length Information

The abuse of traditional antibiotics has led to increased resistance of bacteria and viruses. Efficient therapeutic peptide prediction is critical for peptide drug discovery. However, most of the existing methods only make effective predictions for one class of therapeutic peptides. It is worth noting that currently no predictive method considers sequence length information as a distinct feature of therapeutic peptides. In this article, a novel deep learning approach with matrix factorization for predicting therapeutic peptides (DeepTPpred) by integrating length information are proposed. The matrix factorization layer can learn the potential features of the encoded sequence through the mechanism of first compression and then restoration. And the length features of the sequence of therapeutic peptides are embedded with encoded amino acid sequences. To automatically learn therapeutic peptide predictions, these latent features are input into the neural networks with self-attention mechanism. On eight therapeutic peptide datasets, DeepTPpred achieved excellent prediction results. Based on these datasets, we first integrated eight datasets to obtain a full therapeutic peptide integration dataset. Then, we obtained two functional integration datasets based on the functional similarity of the peptides. Finally, we also conduct experiments on the latest versions of the ACP and CPP datasets. Overall, the experimental results show that our work is effective for the identification of therapeutic peptides.

Spontaneous and information-induced bursting activities in honeybee hives

Social entrainment is important for functioning of beehive organization. By analyzing a dataset of approximately 1000 honeybees (Apis mellifera) tracked in 5 trials, we discovered that honeybees exhibit synchronized activity (bursting behavior) in their locomotion. These bursts occurred spontaneously, potentially as a result of intrinsic bee interactions. The empirical data and simulations demonstrate that physical contact is one of the mechanisms for these bursts. We found that a subset of honeybees within a hive which become active before the peak of each burst, and we refer to these bees as "pioneer bees." Pioneer bees are not selected randomly, but rather, are linked to foraging behavior and waggle dancing, which may help spread external information in the hive. By using transfer entropy, we found that information flows from pioneer bees to non-pioneer bees, which suggest that the bursting behavior is caused by foraging behavior and spreading the information through the hive and promoting integrated group behavior among individuals.

Decision of the Optimal Rank of a Nonnegative Matrix Factorization Model for Gene Expression Data Sets Utilizing the Unit Invariant Knee Method: Development and Evaluation of the Elbow Method for Rank Selection

BACKGROUND

There is a great need to develop a computational approach to analyze and exploit the information contained in gene expression data. The recent utilization of nonnegative matrix factorization (NMF) in computational biology has demonstrated the capability to derive essential details from a high amount of data in particular gene expression microarrays. A common problem in NMF is finding the proper number rank (r) of factors of the degraded demonstration, but no agreement exists on which technique is most appropriate to utilize for this purpose. Thus, various techniques have been suggested to select the optimal value of rank factorization (r).

OBJECTIVE

In this work, a new metric for rank selection is proposed based on the elbow method, which was methodically compared against the cophenetic metric.

METHODS

To decide the optimum number rank (r), this study focused on the unit invariant knee (UIK) method of the NMF on gene expression data sets. Since the UIK method requires an extremum distance estimator that is eventually employed for inflection and identification of a knee point, the proposed method finds the first inflection point of the curvature of the residual sum of squares of the proposed algorithms using the UIK method on gene expression data sets as a target matrix.

RESULTS

Computation was conducted for the UIK task using gene expression data of acute lymphoblastic leukemia and acute myeloid leukemia samples. Consequently, the distinct results of NMF were subjected to comparison on different algorithms. The proposed UIK method is easy to perform, fast, free of a priori rank value input, and does not require initial parameters that significantly influence the model's functionality.

CONCLUSIONS

This study demonstrates that the elbow method provides a credible prediction for both gene expression data and for precisely estimating simulated mutational processes data with known dimensions. The proposed UIK method is faster than conventional methods, including metrics utilizing the consensus matrix as a criterion for rank selection, while achieving significantly better computational efficiency without visual inspection on the curvatives. Finally, the suggested rank tuning method based on the elbow method for gene expression data is arguably theoretically superior to the cophenetic measure.

Polar Nitrated Aromatic Compounds in Urban Fine Particulate Matter: A Focus on Formation via an Aqueous-Phase Radical Mechanism

Polar nitrated aromatic compounds (pNACs) are key ambient brown carbon chromophores; however, their formation mechanisms, especially in the aqueous phase, remain unclear. We developed an advanced technique for pNACs and measured 1764 compounds in atmospheric fine particulate matter sampled in urban Beijing, China. Molecular formulas were derived for 433 compounds, of which 17 were confirmed using reference standards. Potential novel species with up to four aromatic rings and a maximum of five functional groups were found. Higher concentrations were detected in the heating season, with a median of 82.6 ng m^-3 for Σ₁₇pNACs. Non-negative matrix factorization analysis indicated that primary emissions particularly coal combustion were dominant in the heating season. While in the non-heating season, aqueous-phase nitration could generate abundant pNACs with the carboxyl group, which was confirmed by their significant association with the aerosol liquid water content. Aqueous-phase formation of 3- and 5-nitrosalicylic acids instead of their isomer of 4-hydroxy-3-nitrobenzoic acid suggests the existence of an intermediate where the intramolecular hydrogen bond favors kinetics-controlled NO₂^• nitration. This study provides not only a promising technique for the pNAC measurement but also evidence for their atmospheric aqueous-phase formation, facilitating further evaluation of pNACs' climatic effects.

Post COVID-19 pandemic recovery of intracity human mobility in Wuhan: Spatiotemporal characteristic and driving mechanism

After successfully inhibiting the first wave of COVID-19 transmission through a city lockdown, Wuhan implemented a series of policies to gradually lift restrictions and restore daily activities. Existing studies mainly focus on the intercity recovery under a macroscopic view. How does the intracity mobility return to normal? Is the recovery process consistent among different subareas, and what factor affects the post-pandemic recovery? To answer these questions, we sorted out policies adopted during the Wuhan resumption, and collected the long-time mobility big data in 1105 traffic analysis zones (TAZs) to construct an observation matrix (A). We then used the nonnegative matrix factorization (NMF) method to approximate A as the product of two condensed matrices (WH). The column vectors of W matrix were visualized as five typical recovery curves to reveal the temporal change. The row vectors of H matrix were visualized to identify the spatial distribution of each recovery type, and were analyzed with variables of population, GDP, land use, and key facility to explain the recovery driving mechanisms. We found that the "staggered time" policies implemented in Wuhan effectively staggered the peak mobility of several recovery types ("staggered peak"). Besides, different TAZs had heterogeneous response intensities to these policies ("staggered area") which were closely related to land uses and key facilities. The creative policies taken by Wuhan highlight the wisdom of public health crisis management, and could provide an empirical reference for the adjustment of post-pandemic intervention measures in other cities.

Whole-genome sequencing of chronic lymphocytic leukemia identifies subgroups with distinct biological and clinical features

The value of genome-wide over targeted driver analyses for predicting clinical outcomes of cancer patients is debated. Here, we report the whole-genome sequencing of 485 chronic lymphocytic leukemia patients enrolled in clinical trials as part of the United Kingdom's 100,000 Genomes Project. We identify an extended catalog of recurrent coding and noncoding genetic mutations that represents a source for future studies and provide the most complete high-resolution map of structural variants, copy number changes and global genome features including telomere length, mutational signatures and genomic complexity. We demonstrate the relationship of these features with clinical outcome and show that integration of 186 distinct recurrent genomic alterations defines five genomic subgroups that associate with response to therapy, refining conventional outcome prediction. While requiring independent validation, our findings highlight the potential of whole-genome sequencing to inform future risk stratification in chronic lymphocytic leukemia.

Scalable and model-free detection of spatial patterns and colocalization

The expeditious growth in spatial omics technologies enables the profiling of genome-wide molecular events at molecular and single-cell resolution, highlighting a need for fast and reliable methods to characterize spatial patterns. We developed SpaGene, a model-free method to discover spatial patterns rapidly in large-scale spatial omics studies. Analyzing simulation and a variety of spatially resolved transcriptomics data showed that SpaGene is more powerful and scalable than existing methods. Spatial expression patterns identified by SpaGene reconstruct unobserved tissue structures. SpaGene also successfully discovers ligand-receptor interactions through their colocalization.

Protocol for integrative subtyping of lower-grade gliomas using the SUMO pipeline

Grouping patients into subtypes with homogeneous molecular features can guide diagnosis and therapeutic interventions. SUMO is a computational pipeline that uses nonnegative matrix factorization of patient-similarity networks to integrate continuous multi-omic datasets for molecular subtyping of a disease. Here, we present a detailed protocol to demonstrate its use in determining subtypes of lower-grade gliomas by integrating gene expression, DNA methylation, and miRNA expression data from the TCGA-LGG cohort. For complete details on the use and execution of this profile, please refer to Sienkiewicz et al. (2022).

THUNDER: A reference-free deconvolution method to infer cell type proportions from bulk Hi-C data

Hi-C data provide population averaged estimates of three-dimensional chromatin contacts across cell types and states in bulk samples. Effective analysis of Hi-C data entails controlling for the potential confounding factor of differential cell type proportions across heterogeneous bulk samples. We propose a novel unsupervised deconvolution method for inferring cell type composition from bulk Hi-C data, the Two-step Hi-c UNsupervised DEconvolution appRoach (THUNDER). We conducted extensive simulations to test THUNDER based on combining two published single-cell Hi-C (scHi-C) datasets. THUNDER more accurately estimates the underlying cell type proportions compared to reference-free methods (e.g., TOAST, and NMF) and is more robust than reference-dependent methods (e.g. MuSiC). We further demonstrate the practical utility of THUNDER to estimate cell type proportions and identify cell-type-specific interactions in Hi-C data from adult human cortex tissue samples. THUNDER will be a useful tool in adjusting for varying cell type composition in population samples, facilitating valid and more powerful downstream analysis such as differential chromatin organization studies. Additionally, THUNDER estimated contact profiles provide a useful exploratory framework to investigate cell-type-specificity of the chromatin interactome while experimental data is still rare.

Radiomics of Contrast-Enhanced Computed Tomography: A Potential Biomarker for Pretreatment Prediction of the Response to Bacillus Calmette-Guerin Immunotherapy in Non-Muscle-Invasive Bladder Cancer

Background:Bacillus Calmette-Guerin (BCG) instillation is recommended postoperatively after transurethral resection of bladder cancer (TURBT) in patients with high-risk non-muscle-invasive bladder cancer (NMIBC). An accurate prediction model for the BCG response can help identify patients with NMIBC who may benefit from alternative therapy.

Objective: To investigate the value of computed tomography (CT) radiomics features in predicting the response to BCG instillation among patients with primary high-risk NMIBC.

Methods: Patients with pathologically confirmed high-risk NMIBC were retrospectively reviewed. Patients who underwent contrast-enhanced CT examination within one to 2 weeks before TURBT and received ≥5 BCG instillation treatments in two independent hospitals were enrolled. Patients with a routine follow-up of at least 1 year at the outpatient department were included in the final cohort. Radiomics features based on CT images were extracted from the tumor and its periphery in the training cohort, and a radiomics signature was built with recursive feature elimination. Selected features further underwent an unsupervised radiomics analysis using the newly introduced method, non-negative matrix factorization (NMF), to compute factor factorization decompositions of the radiomics matrix. Finally, a robust component, which was most associated with BCG failure in 1 year, was selected. The performance of the selected component was assessed and tested in an external validation cohort.

Results: Overall, 128 patients (training cohort, n = 104; external validation cohort, n = 24) were included, including 12 BCG failures in the training cohort and 11 failures in the validation cohort each. NMF revealed five components, of which component 3 was selected for the best discrimination of BCG failure; it had an area under the curve (AUC) of .79, sensitivity of .79, and specificity of .65 in the training set. In the external validation cohort, it achieved an AUC of .68, sensitivity of .73, and specificity of .69. Survival analysis showed that patients with higher component scores had poor recurrence-free survival (RFS) in both cohorts (C-index: training cohort, .69; validation cohort, .68).

Conclusion: The study suggested that radiomics components based on NMF might be a potential biomarker to predict BCG response and RFS after BCG treatment in patients with high-risk NMIBC.

Topological data analysis gives two folding paths in HP35(nle-nle), double mutant of villin headpiece subdomain

The folding dynamics of proteins is a primary area of interest in protein science. We carried out topological data analysis (TDA) of the folding process of HP35(nle-nle), a double-mutant of the villin headpiece subdomain. Using persistent homology and non-negative matrix factorization, we reduced the dimension of protein structure and investigated the flow in the reduced space. We found this protein has two folding paths, distinguished by the pairings of inter-helix residues. Our analysis showed the excellent performance of TDA in capturing the formation of tertiary structure.

Detecting molecular subtypes from multi-omics datasets using SUMO

We present a data integration framework that uses non-negative matrix factorization of patient-similarity networks to integrate continuous multi-omics datasets for molecular subtyping. It is demonstrated to have the capability to handle missing data without using imputation and to be consistently among the best in detecting subtypes with differential prognosis and enrichment of clinical associations in a large number of cancers. When applying the approach to data from individuals with lower-grade gliomas, we identify a subtype with a significantly worse prognosis. Tumors assigned to this subtype are hypomethylated genome wide with a gain of AP-1 occupancy in demethylated distal enhancers. The tumors are also enriched for somatic chromosome 7 (chr7) gain, chr10 loss, and other molecular events that have been suggested as diagnostic markers for "IDH wild type, with molecular features of glioblastoma" by the cIMPACT-NOW consortium but have yet to be included in the World Health Organization (WHO) guidelines.

Correlated mechanochemical maps of Arabidopsis thaliana primary cell walls using atomic force microscope infrared spectroscopy

Spatial heterogeneity in composition and organisation of the primary cell wall affects the mechanics of cellular morphogenesis. However, directly correlating cell wall composition, organisation and mechanics has been challenging. To overcome this barrier, we applied atomic force microscopy coupled with infrared (AFM-IR) spectroscopy to generate spatially correlated maps of chemical and mechanical properties for paraformaldehyde-fixed, intact Arabidopsis thaliana epidermal cell walls. AFM-IR spectra were deconvoluted by non-negative matrix factorisation (NMF) into a linear combination of IR spectral factors representing sets of chemical groups comprising different cell wall components. This approach enables quantification of chemical composition from IR spectral signatures and visualisation of chemical heterogeneity at nanometer resolution. Cross-correlation analysis of the spatial distribution of NMFs and mechanical properties suggests that the carbohydrate composition of cell wall junctions correlates with increased local stiffness. Together, our work establishes new methodology to use AFM-IR for the mechanochemical analysis of intact plant primary cell walls.

Unsupervised methods in LC-MS data treatment: Application for potential chemotaxonomic markers search

The combination of Liquid Chromatography and Mass Spectrometry (LC-MS) is commonly used to determine and characterize biologically active compounds because of its high resolution and sensitivity. In this work we explore the interpretation of LC-MS data using multivariate statistical analysis algorithms to extract useful chemical information and identify clusters of similar samples. Samples of leaves from 19 plants belonging to the Apiaceae family were analyzed in unified LC conditions by high- and low-resolution mass spectrometry in a wide range scan mode. LC-MS data preprocessing was performed followed by statistical analysis using tensor decomposition in the form of Parallel Factor Analysis (PARAFAC); matrix factorization following tensor unfolding with principal component analysis (PCA), independent component analysis (ICA), non-negative matrix factorization (NMF); or unsupervised feature selection (UFS). The optimal number of components for each of these methods were found and results were compared using four different metrics: silhouette score, Davies-Bouldin index, computational time, number of noisy components. It was found that PCA, ICA and UFS give the best results across the majority of the criteria for both low- and high-resolution data. An algorithm for biomarker signal selection is suggested and 23 potential chemotaxonomic markers were tentatively identified using MS² data. Dendrograms constructed by the methods were compared to the molecular phylogenic tree by calculating pixel-wise mean square error (MSE). Therefore, the suggested approach can support chemotaxonomic studies and yield valuable chemical information for biomarker discovery.

Spatially resolved transcriptomics reveals the architecture of the tumor-microenvironment interface

During tumor progression, cancer cells come into contact with various non-tumor cell types, but it is unclear how tumors adapt to these new environments. Here, we integrate spatially resolved transcriptomics, single-cell RNA-seq, and single-nucleus RNA-seq to characterize tumor-microenvironment interactions at the tumor boundary. Using a zebrafish model of melanoma, we identify a distinct "interface" cell state where the tumor contacts neighboring tissues. This interface is composed of specialized tumor and microenvironment cells that upregulate a common set of cilia genes, and cilia proteins are enriched only where the tumor contacts the microenvironment. Cilia gene expression is regulated by ETS-family transcription factors, which normally act to suppress cilia genes outside of the interface. A cilia-enriched interface is conserved in human patient samples, suggesting it is a conserved feature of human melanoma. Our results demonstrate the power of spatially resolved transcriptomics in uncovering mechanisms that allow tumors to adapt to new environments.

MutVis: Automated framework for analysis and visualization of mutational signatures in pathogenic bacterial strains

In recent years, mutational signature analysis has become a routine practice in cancer genomics for classification and diagnosis. Characterizing mutational signatures across species or within genomes of a bacteria helps in understanding their evolution and adaptation. However, an integrated framework for analysis and visualization of mutational signatures in bacterial genome is lacking. Hence, we aim to develop an integrated, automated, open-source and user-friendly framework called MutVis to analyze mutational signatures from bacterial whole genome next generation sequencing data. The current framework integrates various publicly available packages using Snakemake workflow management software, Python and R scripting. MutVis supports variant calling, transition (Ti) and transversion (Tv) graphical representation, generation of mutational count matrix, graphical visualization of base-pair substitution spectrum (BPSs) and mutation signatures extraction. TvTi plots provide the 6 base substitution classification for both genome and gene level. Further resolution of base pair substitution classification is provided as 96-profile BPSs plot. Mutation signatures is derived based on the characteristic pattern observed in BPSs using non-negative matrix factorization. Relative contribution of signatures is given as hierarchically clustered heatmap. This provides information on active signatures in the individual given sample and classify samples according to signature contributions. We demonstrated the MutVis framework using geographically different strains of Mycobacterium tuberculosis, downloaded from PATRIC TB-ARC Antibiotic Resistance Catalog (n = 963). The current framework can be used to study mutation biases and characteristic mutational signatures in bacterial genomes and is freely available at https://github.com/AkshathaPrasanna/MutVis.

Alternative Polyadenylation in Stem Cell Self-Renewal and Differentiation

Cellular function is shaped by transcriptional and post-transcriptional mechanisms, including alternative polyadenylation (APA). By directly controlling 3'- untranslated region (UTR) length and the selection of the last exon, APA regulates up to 70% of all cellular transcripts influencing RNA stability, output, and protein isoform expression. Cell-state-dependent 3'-UTR shortening has been identified as a hallmark of cellular proliferation. Hence, quiescent/dormant stem cells are characterized by long 3'-UTRs, whereas proliferative stem/progenitor cells exhibit 3'-UTR shortening. Here, the latest studies analyzing the role of APA in regulating stem cell state, self-renewal, differentiation, and metabolism are reviewed. The new role of APA in controlling stem cell fate opens novel potential therapeutic avenues in the field of regenerative medicine.

Aptardi predicts polyadenylation sites in sample-specific transcriptomes using high-throughput RNA sequencing and DNA sequence

Annotation of polyadenylation sites from short-read RNA sequencing alone is a challenging computational task. Other algorithms rooted in DNA sequence predict potential polyadenylation sites; however, in vivo expression of a particular site varies based on a myriad of conditions. Here, we introduce aptardi (alternative polyadenylation transcriptome analysis from RNA-Seq data and DNA sequence information), which leverages both DNA sequence and RNA sequencing in a machine learning paradigm to predict expressed polyadenylation sites. Specifically, as input aptardi takes DNA nucleotide sequence, genome-aligned RNA-Seq data, and an initial transcriptome. The program evaluates these initial transcripts to identify expressed polyadenylation sites in the biological sample and refines transcript 3'-ends accordingly. The average precision of the aptardi model is twice that of a standard transcriptome assembler. In particular, the recall of the aptardi model (the proportion of true polyadenylation sites detected by the algorithm) is improved by over three-fold. Also, the model-trained using the Human Brain Reference RNA commercial standard-performs well when applied to RNA-sequencing samples from different tissues and different mammalian species. Finally, aptardi's input is simple to compile and its output is easily amenable to downstream analyses such as quantitation and differential expression.

An automated approach for determining the number of components in non-negative matrix factorization with application to mutational signature learning

Non-negative matrix factorization (NMF) is a popular method for finding a low rank approximation of a matrix, thereby revealing the latent components behind it. In genomics, NMF is widely used to interpret mutation data and derive the underlying mutational processes and their activities. A key challenge in the use of NMF is determining the number of components, or rank of the factorization. Here we propose a novel method, CV2K, to choose this number automatically from data that is based on a detailed cross validation procedure combined with a parsimony consideration. We apply our method for mutational signature analysis and demonstrate its utility on both simulated and real data sets. In comparison to previous approaches, some of which involve human assessment, CV2K leads to improved predictions across a wide range of data sets.

Non-invasive focusing and imaging in scattering media with a fluorescence-based transmission matrix

In biological microscopy, light scattering represents the main limitation to image at depth. Recently, a set of wavefront shaping techniques has been developed in order to manipulate coherent light in strongly disordered materials. The Transmission Matrix approach has shown its capability to inverse the effect of scattering and efficiently focus light. In practice, the matrix is usually measured using an invasive detector or low-resolution acoustic guide stars. Here, we introduce a non-invasive and all-optical strategy based on linear fluorescence to reconstruct the transmission matrices, to and from a fluorescent object placed inside a scattering medium. It consists in demixing the incoherent patterns emitted by the object using low-rank factorizations and phase retrieval algorithms. We experimentally demonstrate the efficiency of this method through robust and selective focusing. Additionally, from the same measurements, it is possible to exploit memory effect correlations to image and reconstruct extended objects. This approach opens up a new route towards imaging in scattering media with linear or non-linear contrast mechanisms.

Light scattering represents the main limitation to image at depth in biological microscopy. The authors present a strategy to characterize light propagation in and out of a scattering medium based on linear fluorescence feedback and from the same measurements exploit memory effect correlations to image and reconstruct extended objects.

Compact representation of temporal processes in echosounder time series via matrix decomposition

The recent explosion in the availability of echosounder data from diverse ocean platforms has created unprecedented opportunities to observe the marine ecosystems at broad scales. However, the critical lack of methods capable of automatically discovering and summarizing prominent spatio-temporal echogram structures has limited the effective and wider use of these rich datasets. To address this challenge, a data-driven methodology is developed based on matrix decomposition that builds compact representation of long-term echosounder time series using intrinsic features in the data. In a two-stage approach, noisy outliers are first removed from the data by principal component pursuit, then a temporally smooth nonnegative matrix factorization is employed to automatically discover a small number of distinct daily echogram patterns, whose time-varying linear combination (activation) reconstructs the dominant echogram structures. This low-rank representation provides biological information that is more tractable and interpretable than the original data, and is suitable for visualization and systematic analysis with other ocean variables. Unlike existing methods that rely on fixed, handcrafted rules, this unsupervised machine learning approach is well-suited for extracting information from data collected from unfamiliar or rapidly changing ecosystems. This work forms the basis for constructing robust time series analytics for large-scale, acoustics-based biological observation in the ocean.

Principal curve approaches for inferring 3D chromatin architecture

Three-dimensional (3D) genome spatial organization is critical for numerous cellular processes, including transcription, while certain conformation-driven structural alterations are frequently oncogenic. Genome architecture had been notoriously difficult to elucidate, but the advent of the suite of chromatin conformation capture assays, notably Hi-C, has transformed understanding of chromatin structure and provided downstream biological insights. Although many findings have flowed from direct analysis of the pairwise proximity data produced by these assays, there is added value in generating corresponding 3D reconstructions deriving from superposing genomic features on the reconstruction. Accordingly, many methods for inferring 3D architecture from proximity data have been advanced. However, none of these approaches exploit the fact that single chromosome solutions constitute a one-dimensional (1D) curve in 3D. Rather, this aspect has either been addressed by imposition of constraints, which is both computationally burdensome and cell type specific, or ignored with contiguity imposed after the fact. Here, we target finding a 1D curve by extending principal curve methodology to the metric scaling problem. We illustrate how this approach yields a sequence of candidate solutions, indexed by an underlying smoothness or degrees-of-freedom parameter, and propose methods for selection from this sequence. We apply the methodology to Hi-C data obtained on IMR90 cells and so are positioned to evaluate reconstruction accuracy by referencing orthogonal imaging data. The results indicate the utility and reproducibility of our principal curve approach in the face of underlying structural variation.

XPC deficiency increases risk of hematologic malignancies through mutator phenotype and characteristic mutational signature

Recent studies demonstrated a dramatically increased risk of leukemia in patients with a rare genetic disorder, Xeroderma Pigmentosum group C (XP-C), characterized by constitutive deficiency of global genome nucleotide excision repair (GG-NER). The genetic mechanisms of non-skin cancers in XP-C patients remain unexplored. In this study, we analyze a unique collection of internal XP-C tumor genomes including 6 leukemias and 2 sarcomas. We observe a specific mutational pattern and an average of 25-fold increase of mutation rates in XP-C versus sporadic leukemia which we presume leads to its elevated incidence and early appearance. We describe a strong mutational asymmetry with respect to transcription and the direction of replication in XP-C tumors suggesting association of mutagenesis with bulky purine DNA lesions of probably endogenous origin. These findings suggest existence of a balance between formation and repair of bulky DNA lesions by GG-NER in human body cells which is disrupted in XP-C patients.

Xeroderma Pigmentosum group C (XP-C) is a rare genetic disorder characterised by deficient DNA repair leading to skin and internal cancer, but the latter is not well understood molecularly. Here the authors sequence genomes of non-skin cancers from XP-C patients to unravel its mutational patterns.

Identification of three glioblastoma subtypes and a six-gene prognostic risk index based on the expression of growth factors and cytokines

Glioblastoma multiforme (GBM) is the most common and invasive tumor of the central nervous system. Growth factors and cytokines (GFCKs) play a crucial role in tumor invasion. In the present study, GFCK expression profiles from GBM patients in the Chinese Glioma Genome Atlas were used to perform sample clustering with nonnegative matrix factorization. Three GBM subtypes were identified based on differences in GFCK expression, and the subtypes differed in characteristics and prognosis. A prognostic risk index (RI) comprising six GFCKs (BMP2, CCN3, GKN1, LIF, MDK, and SEMA3G) was defined using univariate Cox hazard analysis and multivariate stepwise Cox regression. The RI was validated in two independent data sets and may be independent of some known prognostic factors. Our results suggest that GBM occurs as different subtypes expressing different patterns of GFCKs and that these expression patterns can be captured in an RI that can predict prognosis.

Protein-Folding Analysis Using Features Obtained by Persistent Homology

Understanding the protein-folding process is an outstanding issue in biophysics; recent developments in molecular dynamics simulation have provided insights into this phenomenon. However, the large freedom of atomic motion hinders the understanding of this process. In this study, we applied persistent homology, an emerging method to analyze topological features in a data set, to reveal protein-folding dynamics. We developed a new, to our knowledge, method to characterize the protein structure based on persistent homology and applied this method to molecular dynamics simulations of chignolin. Using principle component analysis or nonnegative matrix factorization, our analysis method revealed two stable states and one saddle state, corresponding to the native, misfolded, and transition states, respectively. We also identified an unfolded state with slow dynamics in the reduced space. Our method serves as a promising tool to understand the protein-folding process.

Heterogeneity of MSI-H gastric cancer identifies a subtype with worse survival

BACKGROUND

Microsatellite instability-high (MSI-H) tumour patients generally have a better prognosis than microsatellite-stable (MSS) ones due to the large number of non-synonymous mutations. However, an increasing number of studies have revealed that less than half of MSI-H patients gain survival benefits or symptom alleviation from immune checkpoint-blockade treatment. Thus, an in-depth inspection of heterogeneous MSI-H tumours is urgently required.

METHODS

Here, we used non-negative matrix factorisation (non-NMF)-based consensus clustering to define stomach adenocarcinoma (STAD) MSI-H subtypes in samples from The Cancer Genome Atlas and an Asian cohort, GSE62254.

RESULTS

MSI-H STAD samples are basically clustered into two subgroups (MSI-H1 and MSI-H2). Further examination of the immune landscape showed that immune suppression factors were enriched in the MSI-H1 subgroup, which may be associated with the poor prognosis in this subgroup.

CONCLUSIONS

Our results illustrate the genetic heterogeneity within MSI-H STADs, with important implications for cancer patient risk stratification, prognosis and treatment.

Revealing Relationships Among Cognitive Functions Using Functional Connectivity and a Large-Scale Meta-Analysis Database

To characterize each cognitive function per se and to understand the brain as an aggregate of those functions, it is vital to relate dozens of these functions to each other. Knowledge about the relationships among cognitive functions is informative not only for basic neuroscientific research but also for clinical applications and developments of brain-inspired artificial intelligence. In the present study, we propose an exhaustive data mining approach to reveal relationships among cognitive functions based on functional brain mapping and network analysis. We began our analysis with 109 pseudo-activation maps (cognitive function maps; CFM) that were reconstructed from a functional magnetic resonance imaging meta-analysis database, each of which corresponds to one of 109 cognitive functions such as ‘emotion,’ ‘attention,’ ‘episodic memory,’ etc. Based on the resting-state functional connectivity between the CFMs, we mapped the cognitive functions onto a two-dimensional space where the relevant functions were located close to each other, which provided a rough picture of the brain as an aggregate of cognitive functions. Then, we conducted so-called conceptual analysis of cognitive functions using clustering of voxels in each CFM connected to the other 108 CFMs with various strengths. As a result, a CFM for each cognitive function was subdivided into several parts, each of which is strongly associated with some CFMs for a subset of the other cognitive functions, which brought in sub-concepts (i.e., sub-functions) of the cognitive function. Moreover, we conducted network analysis for the network whose nodes were parcels derived from whole-brain parcellation based on the whole-brain voxel-to-CFM resting-state functional connectivities. Since each parcel is characterized by associations with the 109 cognitive functions, network analyses using them are expected to inform about relationships between cognitive and network characteristics. Indeed, we found that informational diversities of interaction between parcels and densities of local connectivity were dependent on the kinds of associated functions. In addition, we identified the homogeneous and inhomogeneous network communities about the associated functions. Altogether, we suggested the effectiveness of our approach in which we fused the large-scale meta-analysis of functional brain mapping with the methods of network neuroscience to investigate the relationships among cognitive functions.

Learning a Latent Space of Highly Multidimensional Cancer Data

We introduce a Unified Disentanglement Network (UFDN) trained on The Cancer Genome Atlas (TCGA), which we refer to as UFDN-TCGA. We demonstrate that UFDN-TCGA learns a biologically relevant, low-dimensional latent space of high-dimensional gene expression data by applying our network to two classification tasks of cancer status and cancer type. UFDN-TCGA performs comparably to random forest methods. The UFDN allows for continuous, partial interpolation between distinct cancer types. Furthermore, we perform an analysis of differentially expressed genes between skin cutaneous melanoma (SKCM) samples and the same samples interpolated into glioblastoma (GBM). We demonstrate that our interpolations consist of relevant metagenes that recapitulate known glioblastoma mechanisms.

Implication of DNA repair genes in Lynch-like syndrome

Many colorectal cancers (CRCs) that exhibit microsatellite instability (MSI) are not explained by MLH1 promoter methylation or germline mutations in mismatch repair (MMR) genes, which cause Lynch syndrome (LS). Instead, these Lynch-like syndrome (LLS) patients have somatic mutations in MMR genes. However, many of these patients are young and have relatives with cancer, suggesting a hereditary entity. We performed germline sequence analysis in LLS patients and determined their tumor's mutational profiles using FFPE DNA. Six hundred and fifty-four consecutive CRC patients were screened for suspected LS using MSI and absence of MLH1 methylation. Suspected LS cases were exome sequenced to identify germline and somatic mutations. Single nucleotide variants were used to characterize mutational signatures. We identified 23 suspected LS cases. Germline sequence analysis of 16 available samples identified five cases with LS mutations and 11 cases without LS mutations, LLS. Most LLS tumors had a combination of somatic MMR gene mutation and loss of heterozygosity. LLS patients were relatively young and had excess first-degree relatives with cancer. Four of the 11 LLS patients had rare likely pathogenic variants in genes that maintain genome integrity. Moreover, tumors from this group had a distinct mutational signature compared to tumors from LLS patients lacking germline mutations in these genes. In summary, more than a third of the LLS patients studied had germline mutations in genes that maintain genome integrity and their tumors had a distinct mutational signature. The possibility of hereditary factors in LLS warrants further studies so counseling can be properly informed.

Integrative construction of regulatory region networks in 127 human reference epigenomes by matrix factorization

Despite large experimental and computational efforts aiming to dissect the mechanisms underlying disease risk, mapping cis-regulatory elements to target genes remains a challenge. Here, we introduce a matrix factorization framework to integrate physical and functional interaction data of genomic segments. The framework was used to predict a regulatory network of chromatin interaction edges linking more than 20 000 promoters and 1.8 million enhancers across 127 human reference epigenomes, including edges that are present in any of the input datasets. Our network integrates functional evidence of correlated activity patterns from epigenomic data and physical evidence of chromatin interactions. An important contribution of this work is the representation of heterogeneous data with different qualities as networks. We show that the unbiased integration of independent data sources suggestive of regulatory interactions produces meaningful associations supported by existing functional and physical evidence, correlating with expected independent biological features.

A Spatiotemporal Constraint Non-Negative Matrix Factorization Model to Discover Intra-Urban Mobility Patterns from Taxi Trips

Taxi services provide an urban transport option to citizens. Massive taxi trajectories contain rich information for understanding human travel activities, which are essential to sustainable urban mobility and transportation. The origin and destination (O-D) pairs of urban taxi trips can reveal the spatiotemporal patterns of human mobility and then offer fundamental information to interpret and reform formal, functional, and perceptual regions of cities. Matrices are one of the most effective models to represent taxi trajectories and O-D trips. Among matrix representations, non-negative matrix factorization (NMF) gives meaningful interpretations of complex latent relationships. However, the independence assumption for observations is violated by spatial and temporal autocorrelation in taxi flows, which is not compensated in classical NMF models. In order to discover human intra-urban mobility patterns, a novel spatiotemporal constraint NMF (STC-NMF) model that explicitly solves spatial and temporal dependencies is proposed in this paper. It factorizes taxi flow matrices in both spatial and temporal aspects, thus revealing inherent spatiotemporal patterns. With three-month taxi trajectories harvested in Beijing, China, the STC-NMF model is employed to investigate taxi travel patterns and their spatial interaction modes. As the results, four departure patterns, three arrival patterns, and eight spatial interaction patterns during weekdays and weekends are discovered. Moreover, it is found that intensive movements within certain time windows are significantly related to region functionalities and the spatial interaction flows exhibit an obvious distance decay tendency. The outcome of the proposed model is more consistent with the inherent spatiotemporal characteristics of human intra-urban movements. The knowledge gained in this research would be useful to taxi services and transportation management for promoting sustainable urban development.

Cancer classification and pathway discovery using non-negative matrix factorization

OBJECTIVES

Extracting genetic information from a full range of sequencing data is important for understanding disease. We propose a novel method to effectively explore the landscape of genetic mutations and aggregate them to predict cancer type.

DESIGN

We applied non-smooth non-negative matrix factorization (nsNMF) and support vector machine (SVM) to utilize the full range of sequencing data, aiming to better aggregate genetic mutations and improve their power to predict disease type. More specifically, we introduce a novel classifier to distinguish cancer types using somatic mutations obtained from whole-exome sequencing data. Mutations were identified from multiple cancers and scored using SIFT, PP2, and CADD, and collapsed at the individual gene level. nsNMF was then applied to reduce dimensionality and obtain coefficient and basis matrices. A feature matrix was derived from the obtained matrices to train a classifier for cancer type classification with the SVM model.

RESULTS

We have demonstrated that the classifier was able to distinguish four cancer types with reasonable accuracy. In five-fold cross-validations using mutation counts as features, the average prediction accuracy was 80% (SEM = 0.1%), significantly outperforming baselines and outperforming models using mutation scores as features.

CONCLUSION

Using the factor matrices derived from the nsNMF, we identified multiple genes and pathways that are significantly associated with each cancer type. This study presents a generic and complete pipeline to study the associations between somatic mutations and cancers. The proposed method can be adapted to other studies for disease status classification and pathway discovery.

Subtyping of microsatellite instability-high colorectal cancer

Background

Patients with microsatellite instability-high (MSI-H) colorectal cancer (CRC) generally have a better prognosis than patients with microsatellite stable (MSS) CRC. However, some MSI-H CRC patients do not gain overall survival benefits from immune checkpoint-blockade treatment. In other words, heterogeneity within the subgroup of MSI-H tumors remains poorly understood. Thus, an in-depth molecular characterization of MSI-H tumors is urgently required.

Methods

Here, we use nonnegative matrix factorization (NMF)-based consensus clustering to define CRC MSI-H subtypes in The Cancer Genome Atlas and a French multicenter cohort GSE39582. CIBERSORT was used to calculate the proportions of 22 lymphocytes in tumor tissue in MSI-H subtypes.

Results

MSI-H CRC samples basically clustered into two subgroups (MSI-H1 and MSI-H2). MSI-H1 was characterized by a lower BRAF mutational status, higher frequency of chromosomal instability, global hypomethylation, and worse survival than MSI-H2. Further examination of the immune landscape showed that macrophages of the M2 phenotype were enriched in MSI-H1, which may be associated with poor prognosis in this subgroup.

Conclusions

Our results illustrate the genetic heterogeneity in MSI-H CRCs and macrophages may serve as good targets for anticancer therapy in MSI-H1.

Graphical abstract

Electronic supplementary material

The online version of this article (10.1186/s12964-019-0397-4) contains supplementary material, which is available to authorized users.

Computational approaches for discovery of mutational signatures in cancer

The accumulation of somatic mutations in a genome is the result of the activity of one or more mutagenic processes, each of which leaves its own imprint. The study of these DNA fingerprints, termed mutational signatures, holds important potential for furthering our understanding of the causes and evolution of cancer, and can provide insights of relevance for cancer prevention and treatment. In this review, we focus our attention on the mathematical models and computational techniques that have driven recent advances in the field.

Role of inhibitory control in modulating focal seizure spread

Focal seizure propagation is classically thought to be spatially contiguous. However, distribution of seizures through a large-scale epileptic network has been theorized. Here, we used a multielectrode array, wide field calcium imaging, and two-photon calcium imaging to study focal seizure propagation pathways in an acute rodent neocortical 4-aminopyridine model. Although ictal neuronal bursts did not propagate beyond a 2-3-mm region, they were associated with hemisphere-wide field potential fluctuations and parvalbumin-positive interneuron activity outside the seizure focus. While bicuculline surface application enhanced contiguous seizure propagation, focal bicuculline microinjection at sites distant to the 4-aminopyridine focus resulted in epileptic network formation with maximal activity at the two foci. Our study suggests that both classical and epileptic network propagation can arise from localized inhibition defects, and that the network appearance can arise in the context of normal brain structure without requirement for pathological connectivity changes between sites.

Non-negative Matrix Factorization Reveals Resting-State Cortical Alpha Network Abnormalities in the First-Episode Schizophrenia Spectrum

BACKGROUND

Little is known about neural oscillatory dynamics in first-episode psychosis. Pathophysiology of functional connectivity can be measured through network activity of alpha oscillations, reflecting long-range communication between distal brain regions.

METHODS

Resting magnetoencephalographic activity was collected from 31 individuals with first-episode schizophrenia spectrum psychosis and 22 healthy control individuals. Activity was projected to the realistic cortical surface, based on structural magnetic resonance imaging. The first principal component of activity in 40 Brodmann areas per hemisphere was Hilbert transformed within the alpha range. Non-negative matrix factorization was applied to single-trial alpha phase-locking values from all subjects to determine alpha networks. Within networks, energy and entropy were compared.

RESULTS

Four cortical alpha networks were pathological in individuals with first-episode schizophrenia spectrum psychosis. The networks involved the bilateral anterior and posterior cingulate; left auditory, medial temporal, and cingulate cortex; right inferior frontal gyrus and widespread areas; and right posterior parietal cortex and widespread areas. Energy and entropy were associated with the Positive and Negative Syndrome Scale total and thought disorder factors for the first three networks. In addition, the left posterior temporal network was associated with positive and negative factors, and the right inferior frontal network was associated with the positive factor.

CONCLUSIONS

Machine learning network analysis of resting alpha-band neural activity identified several aberrant networks in individuals with first-episode schizophrenia spectrum psychosis, including the left temporal, right inferior frontal, right posterior parietal, and bilateral cingulate cortices. Abnormal long-range alpha communication is evident at the first presentation for psychosis and may provide clues about mechanisms of dysconnectivity in psychosis and novel targets for noninvasive brain stimulation.

Modeling of Genome-Wide Polyadenylation Signals in Xenopus tropicalis

Alternative polyadenylation (APA) is an important post-transcriptional modification event to process messenger RNA (mRNA) for transcriptional termination, transport, and translation. In the present study, we characterized poly(A) signals in Xenopus tropicalis using 70,918 highly confident poly(A) sites derived from 16,511 protein-coding genes to understand their roles in the regulation of embryo development and gender difference. We examined potential factors, including the gene length, the number of introns in a gene, and the intron length, that may affect the prevalence of APA. We observed 12 prominent poly(A) signal patterns, which accounted for approximately 92% of total APA sites in Xenopus tropicalis. Among them, three patterns are specific to X. tropicalis, so they are absent in other animals such as humans or mice. We catalogued APA sites based on their genomic regions and developed a bioinformatics pipeline to identify over-represented signal patterns for each class. Then the schema of cis elements for APA sites in each genomic region was proposed. More importantly, APA usage is dramatically dynamic in embryos along five developmental stages and well-coordinated with the maternal-to-zygotic transition event. We used an entropy-based method to identify developmental stage-specific APA sites and identified significant signal patterns around specific sites and constitutive sites. We found that the APA frequency in different genomic regions varies with developmental stages and that those sites located in intron or coding sequence regions contribute most to the dynamics of gene expression during developmental stages. This study deciphers the characteristics and poly(A) signal patterns for both canonical APA sites and non-canonical APA sites across different developmental stages and gender dimorphisms in X. tropicalis, providing new insights into the dynamic regulation of distal and proximal APA.

Molecular Classification of Epithelial Ovarian Cancer Based on Methylation Profiling: Evidence for Survival Heterogeneity

PURPOSE

Ovarian cancer is a heterogeneous disease that can be divided into multiple subtypes with variable etiology, pathogenesis, and prognosis. We analyzed DNA methylation profiling data to identify biologic subgroups of ovarian cancer and study their relationship with histologic subtypes, copy number variation, RNA expression data, and outcomes.

EXPERIMENTAL DESIGN

A total of 162 paraffin-embedded ovarian epithelial tumor tissues, including the five major epithelial ovarian tumor subtypes (high- and low-grade serous, endometrioid, mucinous, and clear cell) and tumors of low malignant potential were selected from two different sources: The Polish Ovarian Cancer study, and the Surveillance, Epidemiology, and End Results Residual Tissue Repository (SEER RTR). Analyses were restricted to Caucasian women. Methylation profiling was conducted using the Illumina 450K methylation array. For 45 tumors array copy number data were available. NanoString gene expression data for 39 genes were available for 61 high-grade serous carcinomas (HGSC).

RESULTS

Consensus nonnegative matrix factorization clustering of the 1,000 most variable CpG sites showed four major clusters among all epithelial ovarian cancers. We observed statistically significant differences in survival (log-rank test, P = 9.1 × 10^-7) and genomic instability across these clusters. Within HGSC, clustering showed three subgroups with survival differences (log-rank test, P = 0.002). Comparing models with and without methylation subgroups in addition to previously identified gene expression subtypes suggested that the methylation subgroups added significant survival information (P = 0.007).

CONCLUSIONS

DNA methylation profiling of ovarian cancer identified novel molecular subgroups that had significant survival difference and provided insights into the molecular underpinnings of ovarian cancer.See related commentary by Ishak et al., p. 5729.