Proteome-scale movements and compartment connectivity during the eukaryotic cell cycle.

Cell cycle progression relies on coordinated changes in the composition and subcellular localization of the proteome. By applying two distinct convolutional neural networks on images of millions of live yeast cells, we resolved proteome-level dynamics in both concentration and localization during the cell cycle, with resolution of ∼20 subcellular localization classes. We show that a quarter of the proteome displays cell cycle periodicity, with proteins tending to be controlled either at the level of localization or concentration, but not both. Distinct levels of protein regulation are preferentially utilized for different aspects of the cell cycle, with changes in protein concentration being mostly involved in cell cycle control and changes in protein localization in the biophysical implementation of the cell cycle program. We present a resource for exploring global proteome dynamics during the cell cycle, which will aid in understanding a fundamental biological process at a systems level.

Mapping transcriptomic vector fields of single cells.

Single-cell (sc)RNA-seq, together with RNA velocity and metabolic labeling, reveals cellular states and transitions at unprecedented resolution. Fully exploiting these data, however, requires kinetic models capable of unveiling governing regulatory functions. Here, we introduce an analytical framework dynamo (https://github.com/aristoteleo/dynamo-release), which infers absolute RNA velocity, reconstructs continuous vector fields that predict cell fates, employs differential geometry to extract underlying regulations, and ultimately predicts optimal reprogramming paths and perturbation outcomes. We highlight dynamo's power to overcome fundamental limitations of conventional splicing-based RNA velocity analyses to enable accurate velocity estimations on a metabolically labeled human hematopoiesis scRNA-seq dataset. Furthermore, differential geometry analyses reveal mechanisms driving early megakaryocyte appearance and elucidate asymmetrical regulation within the PU.1-GATA1 circuit. Leveraging the least-action-path method, dynamo accurately predicts drivers of numerous hematopoietic transitions. Finally, in silico perturbations predict cell-fate diversions induced by gene perturbations. Dynamo, thus, represents an important step in advancing quantitative and predictive theories of cell-state transitions.

Budding epithelial morphogenesis driven by cell-matrix versus cell-cell adhesion.

Many embryonic organs undergo epithelial morphogenesis to form tree-like hierarchical structures. However, it remains unclear what drives the budding and branching of stratified epithelia, such as in the embryonic salivary gland and pancreas. Here, we performed live-organ imaging of mouse embryonic salivary glands at single-cell resolution to reveal that budding morphogenesis is driven by expansion and folding of a distinct epithelial surface cell sheet characterized by strong cell-matrix adhesions and weak cell-cell adhesions. Profiling of single-cell transcriptomes of this epithelium revealed spatial patterns of transcription underlying these cell adhesion differences. We then synthetically reconstituted budding morphogenesis by experimentally suppressing E-cadherin expression and inducing basement membrane formation in 3D spheroid cultures of engineered cells, which required ß1-integrin-mediated cell-matrix adhesion for successful budding. Thus, stratified epithelial budding, the key first step of branching morphogenesis, is driven by an overall combination of strong cell-matrix adhesion and weak cell-cell adhesion by peripheral epithelial cells.

Gene-Specific Control of tRNA Expression by RNA Polymerase II.

Increasing evidence suggests that tRNA levels are dynamically and specifically regulated in response to internal and external cues to modulate the cellular translational program. However, the molecular players and the mechanisms regulating the gene-specific expression of tRNAs are still unknown. Using an inducible auxin-degron system to rapidly deplete RPB1 (the largest subunit of RNA Pol II) in living cells, we identified Pol II as a direct gene-specific regulator of tRNA transcription. Our data suggest that Pol II transcription robustly interferes with Pol III function at specific tRNA genes. This activity was further found to be essential for MAF1-mediated repression of a large set of tRNA genes during serum starvation, indicating that repression of tRNA genes by Pol II is dynamically regulated. Hence, Pol II plays a direct and central role in the gene-specific regulation of tRNA expression.

Differential proliferation regulates multi-tissue morphogenesis during embryonic axial extension: integrating viscous modeling and experimental approaches.

A major challenge in biology is to understand how mechanical interactions and cellular behavior affect the shapes of tissues and embryo morphology. The extension of the neural tube and paraxial mesoderm, which form the spinal cord and musculoskeletal system, respectively, results in the elongated shape of the vertebrate embryonic body. Despite our understanding of how each of these tissues elongates independently of the others, the morphogenetic consequences of their simultaneous growth and mechanical interactions are still unclear. Our study investigates how differential growth, tissue biophysical properties and mechanical interactions affect embryonic morphogenesis during axial extension using a 2D multi-tissue continuum-based mathematical model. Our model captures the dynamics observed in vivo by time-lapse imaging of bird embryos, and reveals the underestimated influence of differential tissue proliferation rates. We confirmed this prediction in quail embryos by showing that decreasing the rate of cell proliferation in the paraxial mesoderm affects long-term tissue dynamics, and shaping of both the paraxial mesoderm and the neighboring neural tube. Overall, our work provides a new theoretical platform upon which to consider the long-term consequences of tissue differential growth and mechanical interactions on morphogenesis.

BRCA1 and ELK-1 regulate neural progenitor cell fate in the optic tectum in response to visual experience in Xenopus laevis tadpoles.

In developing Xenopus tadpoles, the optic tectum begins to receive patterned visual input while visuomotor circuits are still undergoing neurogenesis and circuit assembly. This visual input regulates neural progenitor cell fate decisions such that maintaining tadpoles in the dark increases proliferation, expanding the progenitor pool, while visual stimulation promotes neuronal differentiation. To identify regulators of activity-dependent neural progenitor cell fate, we profiled the transcriptomes of proliferating neural progenitor cells and newly differentiated neurons using RNA-Seq. We used advanced bioinformatic analysis of 1,130 differentially expressed transcripts to identify six differentially regulated transcriptional regulators, including Breast Cancer 1 (BRCA1) and the ETS-family transcription factor, ELK-1, which are predicted to regulate the majority of the other differentially expressed transcripts. BRCA1 is known for its role in cancers, but relatively little is known about its potential role in regulating neural progenitor cell fate. ELK-1 is a multifunctional transcription factor which regulates immediate early gene expression. We investigated the potential functions of BRCA1 and ELK-1 in activity-regulated neurogenesis in the tadpole visual system using in vivo time-lapse imaging to monitor the fate of GFP-expressing SOX2+ neural progenitor cells in the optic tectum. Our longitudinal in vivo imaging analysis showed that knockdown of either BRCA1 or ELK-1 altered the fates of neural progenitor cells and furthermore that the effects of visual experience on neurogenesis depend on BRCA1 and ELK-1 expression. These studies provide insight into the potential mechanisms by which neural activity affects neural progenitor cell fate.

A general approach for selection of epitope-directed binders to proteins.

Directing antibodies to a particular epitope among many possible on a target protein is a significant challenge. Here, we present a simple and general method for epitope-directed selection (EDS) using a differential phage selection strategy. This involves engineering the protein of interest (POI) with the epitope of interest (EOI) mutated using a systematic bioinformatics algorithm to guide the local design of an EOI decoy variant. Using several alternating rounds of negative selection with the EOI decoy variant followed by positive selection on the wild-type POI, we were able to identify highly specific and potent antibodies to five different EOI antigens that bind and functionally block known sites of proteolysis. Among these, we developed highly specific antibodies that target the proteolytic site on the CUB domain containing protein 1 (CDCP1) to prevent its proteolysis allowing us to study the cellular maturation of this event that triggers malignancy. We generated antibodies that recognize the junction between the pro- and catalytic domains for three different matrix metalloproteases (MMPs), MMP1, MMP3, and MMP9, that selectively block activation of each of these enzymes and impair cell migration. We targeted a proteolytic epitope on the cell surface receptor, EPH Receptor A2 (EphA2), that is known to transform it from a tumor suppressor to an oncoprotein. We believe that the EDS method greatly facilitates the generation of antibodies to specific EOIs on a wide range of proteins and enzymes for broad therapeutic and diagnostic applications.

betAS: intuitive analysis and visualization of differential alternative splicing using beta distributions.

Next-generation RNA sequencing allows alternative splicing (AS) quantification with unprecedented resolution, with the relative inclusion of an alternative sequence in transcripts being commonly quantified by the proportion of reads supporting it as percent spliced-in (PSI). However, PSI values do not incorporate information about precision, proportional to the respective AS events' read coverage. Beta distributions are suitable to quantify inclusion levels of alternative sequences, using reads supporting their inclusion and exclusion as surrogates for the two distribution shape parameters. Each such beta distribution has the PSI as its mean value and is narrower when the read coverage is higher, facilitating the interpretability of its precision when plotted. We herein introduce a computational pipeline, based on beta distributions accurately modeling PSI values and their precision, to quantitatively and visually compare AS between groups of samples. Our methodology includes a differential splicing significance metric that compromises the magnitude of intergroup differences, the estimation uncertainty in individual samples, and the intragroup variability, being therefore suitable for multiple-group comparisons. To make our approach accessible and clear to both noncomputational and computational biologists, we developed betAS, an interactive web app and user-friendly R package for visual and intuitive differential splicing analysis from read count data.

mbDecoda: a debiased approach to compositional data analysis for microbiome surveys.

Potentially pathogenic or probiotic microbes can be identified by comparing their abundance levels between healthy and diseased populations, or more broadly, by linking microbiome composition with clinical phenotypes or environmental factors. However, in microbiome studies, feature tables provide relative rather than absolute abundance of each feature in each sample, as the microbial loads of the samples and the ratios of sequencing depth to microbial load are both unknown and subject to considerable variation. Moreover, microbiome abundance data are count-valued, often over-dispersed and contain a substantial proportion of zeros. To carry out differential abundance analysis while addressing these challenges, we introduce mbDecoda, a model-based approach for debiased analysis of sparse compositions of microbiomes. mbDecoda employs a zero-inflated negative binomial model, linking mean abundance to the variable of interest through a log link function, and it accommodates the adjustment for confounding factors. To efficiently obtain maximum likelihood estimates of model parameters, an Expectation Maximization algorithm is developed. A minimum coverage interval approach is then proposed to rectify compositional bias, enabling accurate and reliable absolute abundance analysis. Through extensive simulation studies and analysis of real-world microbiome datasets, we demonstrate that mbDecoda compares favorably with state-of-the-art methods in terms of effectiveness, robustness and reproducibility.

Normalization of RNA-Seq data using adaptive trimmed mean with multi-reference.

The normalization of RNA sequencing data is a primary step for downstream analysis. The most popular method used for the normalization is the trimmed mean of M values (TMM) and DESeq. The TMM tries to trim away extreme log fold changes of the data to normalize the raw read counts based on the remaining non-deferentially expressed genes. However, the major problem with the TMM is that the values of trimming factor M are heuristic. This paper tries to estimate the adaptive value of M in TMM based on Jaeckel's Estimator, and each sample acts as a reference to find the scale factor of each sample. The presented approach is validated on SEQC, MAQC2, MAQC3, PICKRELL and two simulated datasets with two-group and three-group conditions by varying the percentage of differential expression and the number of replicates. The performance of the present approach is compared with various state-of-the-art methods, and it is better in terms of area under the receiver operating characteristic curve and differential expression.

Optimizing in silico drug discovery: simulation of connected differential expression signatures and applications to benchmarking.

BACKGROUND: We present a novel simulation method for generating connected differential expression signatures. Traditional methods have struggled with the lack of reliable benchmarking data and biases in drug-disease pair labeling, limiting the rigorous benchmarking of connectivity-based approaches. OBJECTIVE: Our aim is to develop a simulation method based on a statistical framework that allows for adjustable levels of parametrization, especially the connectivity, to generate a pair of interconnected differential signatures. This could help to address the issue of benchmarking data availability for connectivity-based drug repurposing approaches. METHODS: We first detailed the simulation process and how it reflected real biological variability and the interconnectedness of gene expression signatures. Then, we generated several datasets to enable the evaluation of different existing algorithms that compare differential expression signatures, providing insights into their performance and limitations. RESULTS: Our findings demonstrate the ability of our simulation to produce realistic data, as evidenced by correlation analyses and the log2 fold-change distribution of deregulated genes. Benchmarking reveals that methods like extreme cosine similarity and Pearson correlation outperform others in identifying connected signatures. CONCLUSION: Overall, our method provides a reliable tool for simulating differential expression signatures. The data simulated by our tool encompass a wide spectrum of possibilities to challenge and evaluate existing methods to estimate connectivity scores. This may represent a critical gap in connectivity-based drug repurposing research because reliable benchmarking data are essential for assessing and advancing in the development of new algorithms. The simulation tool is available as a R package (General Public License (GPL) license) at https://github.com/cgonzalez-gomez/cosimu.

dCCA: detecting differential covariation patterns between two types of high-throughput omics data.

MOTIVATION: The advent of multimodal omics data has provided an unprecedented opportunity to systematically investigate underlying biological mechanisms from distinct yet complementary angles. However, the joint analysis of multi-omics data remains challenging because it requires modeling interactions between multiple sets of high-throughput variables. Furthermore, these interaction patterns may vary across different clinical groups, reflecting disease-related biological processes. RESULTS: We propose a novel approach called Differential Canonical Correlation Analysis (dCCA) to capture differential covariation patterns between two multivariate vectors across clinical groups. Unlike classical Canonical Correlation Analysis, which maximizes the correlation between two multivariate vectors, dCCA aims to maximally recover differentially expressed multivariate-to-multivariate covariation patterns between groups. We have developed computational algorithms and a toolkit to sparsely select paired subsets of variables from two sets of multivariate variables while maximizing the differential covariation. Extensive simulation analyses demonstrate the superior performance of dCCA in selecting variables of interest and recovering differential correlations. We applied dCCA to the Pan-Kidney cohort from the Cancer Genome Atlas Program database and identified differentially expressed covariations between noncoding RNAs and gene expressions. AVAILABILITY AND IMPLEMENTATION: The R package that implements dCCA is available at https://github.com/hwiyoungstat/dCCA.

Quantitative Control for Stoichiometric Protein Synthesis.

Bacterial protein synthesis rates have evolved to maintain preferred stoichiometries at striking precision, from the components of protein complexes to constituents of entire pathways. Setting relative protein production rates to be well within a factor of two requires concerted tuning of transcription, RNA turnover, and translation, allowing many potential regulatory strategies to achieve the preferred output. The last decade has seen a greatly expanded capacity for precise interrogation of each step of the central dogma genome-wide. Here, we summarize how these technologies have shaped the current understanding of diverse bacterial regulatory architectures underpinning stoichiometric protein synthesis. We focus on the emerging expanded view of bacterial operons, which encode diverse primary and secondary mRNA structures for tuning protein stoichiometry. Emphasis is placed on how quantitative tuning is achieved. We discuss the challenges and open questions in the application of quantitative, genome-wide methodologies to the problem of precise protein production.

Plasticity mechanisms of genetically distinct Purkinje cells.

Despite its uniform appearance, the cerebellar cortex is highly heterogeneous in terms of structure, genetics and physiology. Purkinje cells (PCs), the principal and sole output neurons of the cerebellar cortex, can be categorized into multiple populations that differentially express molecular markers and display distinctive physiological features. Such features include action potential rate, but also their propensity for synaptic and intrinsic plasticity. However, the precise molecular and genetic factors that correlate with the differential physiological properties of PCs remain elusive. In this article, we provide a detailed overview of the cellular mechanisms that regulate PC activity and plasticity. We further perform a pathway analysis to highlight how molecular characteristics of specific PC populations may influence their physiology and plasticity mechanisms.

msqrob2PTM: Differential Abundance and Differential Usage Analysis of MS-Based Proteomics Data at the Posttranslational Modification and Peptidoform Level.

In the era of open-modification search engines, more posttranslational modifications than ever can be detected by LC-MS/MS-based proteomics. This development can switch proteomics research into a higher gear, as PTMs are key in many cellular pathways important in cell proliferation, migration, metastasis, and aging. However, despite these advances in modification identification, statistical methods for PTM-level quantification and differential analysis have yet to catch up. This absence can partly be explained by statistical challenges inherent to the data, such as the confounding of PTM intensities with its parent protein abundance. Therefore, we have developed msqrob2PTM, a new workflow in the msqrob2 universe capable of differential abundance analysis at the PTM and at the peptidoform level. The latter is important for validating PTMs found as significantly differential. Indeed, as our method can deal with multiple PTMs per peptidoform, there is a possibility that significant PTMs stem from one significant peptidoform carrying another PTM, hinting that it might be the other PTM driving the perceived differential abundance. Our workflows can flag both differential peptidoform abundance (DPA) and differential peptidoform usage (DPU). This enables a distinction between direct assessment of differential abundance of peptidoforms (DPA) and differences in the relative usage of peptidoforms corrected for corresponding protein abundances (DPU). For DPA, we directly model the log2-transformed peptidoform intensities, while for DPU, we correct for parent protein abundance by an intermediate normalization step which calculates the log2-ratio of the peptidoform intensities to their summarized parent protein intensities. We demonstrated the utility and performance of msqrob2PTM by applying it to datasets with known ground truth, as well as to biological PTM-rich datasets. Our results show that msqrob2PTM is on par with, or surpassing the performance of, the current state-of-the-art methods. Moreover, msqrob2PTM is currently unique in providing output at the peptidoform level.

Proteomics-Derived Biomarker Panel Facilitates Distinguishing Primary Lung Adenocarcinomas With Intestinal or Mucinous Differentiation From Lung Metastatic Colorectal Cancer.

The diagnosis of primary lung adenocarcinomas with intestinal or mucinous differentiation (PAIM) remains challenging due to the overlapping histomorphological, immunohistochemical (IHC), and genetic characteristics with lung metastatic colorectal cancer (lmCRC). This study aimed to explore the protein biomarkers that could distinguish between PAIM and lmCRC. To uncover differences between the two diseases, we used tandem mass tagging-based shotgun proteomics to characterize proteomes of formalin-fixed, paraffin-embedded tumor samples of PAIM (n = 22) and lmCRC (n = 17).Then three machine learning algorithms, namely support vector machine (SVM), random forest, and the Least Absolute Shrinkage and Selection Operator, were utilized to select protein features with diagnostic significance. These candidate proteins were further validated in an independent cohort (PAIM, n = 11; lmCRC, n = 19) by IHC to confirm their diagnostic performance. In total, 105 proteins out of 7871 proteins were significantly dysregulated between PAIM and lmCRC samples and well-separated two groups by Uniform Manifold Approximation and Projection. The upregulated proteins in PAIM were involved in actin cytoskeleton organization, platelet degranulation, and regulation of leukocyte chemotaxis, while downregulated ones were involved in mitochondrial transmembrane transport, vasculature development, and stem cell proliferation. A set of ten candidate proteins (high-level expression in lmCRC: CDH17, ATP1B3, GLB1, OXNAD1, LYST, FABP1; high-level expression in PAIM: CK7 (an established marker), NARR, MLPH, S100A14) was ultimately selected to distinguish PAIM from lmCRC by machine learning algorithms. We further confirmed using IHC that the five protein biomarkers including CDH17, CK7, MLPH, FABP1 and NARR were effective biomarkers for distinguishing PAIM from lmCRC. Our study depicts PAIM-specific proteomic characteristics and demonstrates the potential utility of new protein biomarkers for the differential diagnosis of PAIM and lmCRC. These findings may contribute to improving the diagnostic accuracy and guide appropriate treatments for these patients.

Quantification Quality Control Emerges as a Crucial Factor to Enhance Single-Cell Proteomics Data Analysis.

Mass spectrometry (MS)-based single-cell proteomics (SCP) provides us the opportunity to unbiasedly explore biological variability within cells without the limitation of antibody availability. This field is rapidly developed with the main focuses on instrument advancement, sample preparation refinement, and signal boosting methods; however, the optimal data processing and analysis are rarely investigated which holds an arduous challenge because of the high proportion of missing values and batch effect. Here, we introduced a quantification quality control to intensify the identification of differentially expressed proteins (DEPs) by considering both within and across SCP data. Combining quantification quality control with isobaric matching between runs (IMBR) and PSM-level normalization, an additional 12% and 19% of proteins and peptides, with more than 90% of proteins/peptides containing valid values, were quantified. Clearly, quantification quality control was able to reduce quantification variations and q-values with the more apparent cell type separations. In addition, we found that PSM-level normalization performed similar to other protein-level normalizations but kept the original data profiles without the additional requirement of data manipulation. In proof of concept of our refined pipeline, six uniquely identified DEPs exhibiting varied fold-changes and playing critical roles for melanoma and monocyte functionalities were selected for validation using immunoblotting. Five out of six validated DEPs showed an identical trend with the SCP dataset, emphasizing the feasibility of combining the IMBR, cell quality control, and PSM-level normalization in SCP analysis, which is beneficial for future SCP studies.

Confidence-ranked reconstruction of census microdata from published statistics.

A reconstruction attack on a private dataset D takes as input some publicly accessible information about the dataset and produces a list of candidate elements of D. We introduce a class of data reconstruction attacks based on randomized methods for nonconvex optimization. We empirically demonstrate that our attacks can not only reconstruct full rows of D from aggregate query statistics Q(D)∈ℝm but can do so in a way that reliably ranks reconstructed rows by their odds of appearing in the private data, providing a signature that could be used for prioritizing reconstructed rows for further actions such as identity theft or hate crime. We also design a sequence of baselines for evaluating reconstruction attacks. Our attacks significantly outperform those that are based only on access to a public distribution or population from which the private dataset D was sampled, demonstrating that they are exploiting information in the aggregate statistics Q(D) and not simply the overall structure of the distribution. In other words, the queries Q(D) are permitting reconstruction of elements of this dataset, not the distribution from which D was drawn. These findings are established both on 2010 US decennial Census data and queries and Census-derived American Community Survey datasets. Taken together, our methods and experiments illustrate the risks in releasing numerically precise aggregate statistics of a large dataset and provide further motivation for the careful application of provably private techniques such as differential privacy.

An empirical Bayes method for differential expression analysis of single cells with deep generative models.

Detecting differentially expressed genes is important for characterizing subpopulations of cells. In scRNA-seq data, however, nuisance variation due to technical factors like sequencing depth and RNA capture efficiency obscures the underlying biological signal. Deep generative models have been extensively applied to scRNA-seq data, with a special focus on embedding cells into a low-dimensional latent space and correcting for batch effects. However, little attention has been paid to the problem of utilizing the uncertainty from the deep generative model for differential expression (DE). Furthermore, the existing approaches do not allow for controlling for effect size or the false discovery rate (FDR). Here, we present lvm-DE, a generic Bayesian approach for performing DE predictions from a fitted deep generative model, while controlling the FDR. We apply the lvm-DE framework to scVI and scSphere, two deep generative models. The resulting approaches outperform state-of-the-art methods at estimating the log fold change in gene expression levels as well as detecting differentially expressed genes between subpopulations of cells.

Allostatic load in children: The cost of empathic concern.

Early-life adversity affects long-term health outcomes but there is considerable interindividual variability in susceptibility to environmental influences. We proposed that positive psychological characteristics that reflect engagement with context, such as being concerned about people or performance on tasks (i.e., empathic concern), could moderate the interindividual variation in sensitivity to the quality of the early environment. We studied 526 children of various Asian nationalities in Singapore (46.6% female, 13.4% below the poverty line) with longitudinal data on perinatal and childhood experiences, maternal report on empathic concern of the child, and a comprehensive set of physiological measures reflecting pediatric allostatic load assessed at 6 y of age. The perinatal and childhood experiences included adversities and positive experiences. We found that cumulative adverse childhood experience was positively associated with allostatic load of children at 6 y of age at higher levels of empathic concern but not significantly associated at lower levels of empathic concern. This finding reveals evidence for the importance of empathic concern as a psychological characteristic that moderates the developmental impact of environmental influences, serving as a source for vulnerability to adversities in children.