Divergent Learning-Related Transcriptional States of Cortical Glutamatergic Neurons.

Experience-dependent gene expression reshapes neural circuits, permitting the learning of knowledge and skills. Most learning involves repetitive experiences during which neurons undergo multiple stages of functional and structural plasticity. Currently, the diversity of transcriptional responses underlying dynamic plasticity during repetition-based learning is poorly understood. To close this gap, we analyzed single-nucleus transcriptomes of L2/3 glutamatergic neurons of the primary motor cortex after 3 d motor skill training or home cage control in water-restricted male mice. "Train" and "control" neurons could be discriminated with high accuracy based on expression patterns of many genes, indicating that recent experience leaves a widespread transcriptional signature across L2/3 neurons. These discriminating genes exhibited divergent modes of coregulation, differentiating neurons into discrete clusters of transcriptional states. Several states showed gene expressions associated with activity-dependent plasticity. Some of these states were also prominent in the previously published reference, suggesting that they represent both spontaneous and task-related plasticity events. Markedly, however, two states were unique to our dataset. The first state, further enriched by motor training, showed gene expression suggestive of late-stage plasticity with repeated activation, which is suitable for expected emergent neuronal ensembles that stably retain motor learning. The second state, equally found in both train and control mice, showed elevated levels of metabolic pathways and norepinephrine sensitivity, suggesting a response to common experiences specific to our experimental conditions, such as water restriction or circadian rhythm. Together, we uncovered divergent transcriptional responses across L2/3 neurons, each potentially linked with distinct features of repetition-based motor learning such as plasticity, memory, and motivation.

Global genome decompaction leads to stochastic activation of gene expression as a first step toward fate commitment in human hematopoietic cells.

When human cord blood-derived CD34+ cells are induced to differentiate, they undergo rapid and dynamic morphological and molecular transformations that are critical for fate commitment. In particular, the cells pass through a transitory phase known as "multilineage-primed" state. These cells are characterized by a mixed gene expression profile, different in each cell, with the coexpression of many genes characteristic for concurrent cell lineages. The aim of our study is to understand the mechanisms of the establishment and the exit from this transitory state. We investigated this issue using single-cell RNA sequencing and ATAC-seq. Two phases were detected. The first phase is a rapid and global chromatin decompaction that makes most of the gene promoters in the genome accessible for transcription. It results 24 h later in enhanced and pervasive transcription of the genome leading to the concomitant increase in the cell-to-cell variability of transcriptional profiles. The second phase is the exit from the multilineage-primed phase marked by a slow chromatin closure and a subsequent overall down-regulation of gene transcription. This process is selective and results in the emergence of coherent expression profiles corresponding to distinct cell subpopulations. The typical time scale of these events spans 48 to 72 h. These observations suggest that the nonspecificity of genome decompaction is the condition for the generation of a highly variable multilineage expression profile. The nonspecific phase is followed by specific regulatory actions that stabilize and maintain the activity of key genes, while the rest of the genome becomes repressed again by the chromatin recompaction. Thus, the initiation of differentiation is reminiscent of a constrained optimization process that associates the spontaneous generation of gene expression diversity to subsequent regulatory actions that maintain the activity of some genes, while the rest of the genome sinks back to the repressive closed chromatin state.

GlycoEnzOnto: a GlycoEnzyme pathway and molecular function ontology.

MOTIVATION: The 'glycoEnzymes' include a set of proteins having related enzymatic, metabolic, transport, structural and cofactor functions. Currently, there is no established ontology to describe glycoEnzyme properties and to relate them to glycan biosynthesis pathways. RESULTS: We present GlycoEnzOnto, an ontology describing 403 human glycoEnzymes curated along 139 glycosylation pathways, 134 molecular functions and 22 cellular compartments. The pathways described regulate nucleotide-sugar metabolism, glycosyl-substrate/donor transport, glycan biosynthesis and degradation. The role of each enzyme in the glycosylation initiation, elongation/branching and capping/termination phases is described. IUPAC linear strings present systematic human/machine-readable descriptions of individual reaction steps and enable automated knowledge-based curation of biochemical networks. All GlycoEnzOnto knowledge is integrated with the Gene Ontology biological processes. GlycoEnzOnto enables improved transcript overrepresentation analyses and glycosylation pathway identification compared to other available schema, e.g. KEGG and Reactome. Overall, GlycoEnzOnto represents a holistic glycoinformatics resource for systems-level analyses. AVAILABILITY AND IMPLEMENTATION: https://github.com/neel-lab/GlycoEnzOnto. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

ΔFBA-Predicting metabolic flux alterations using genome-scale metabolic models and differential transcriptomic data.

Genome-scale metabolic models (GEMs) provide a powerful framework for simulating the entire set of biochemical reactions in a cell using a constraint-based modeling strategy called flux balance analysis (FBA). FBA relies on an assumed metabolic objective for generating metabolic fluxes using GEMs. But, the most appropriate metabolic objective is not always obvious for a given condition and is likely context-specific, which often complicate the estimation of metabolic flux alterations between conditions. Here, we propose a new method, called ΔFBA (deltaFBA), that integrates differential gene expression data to evaluate directly metabolic flux differences between two conditions. Notably, ΔFBA does not require specifying the cellular objective. Rather, ΔFBA seeks to maximize the consistency and minimize inconsistency between the predicted flux differences and differential gene expression. We showcased the performance of ΔFBA through several case studies involving the prediction of metabolic alterations caused by genetic and environmental perturbations in Escherichia coli and caused by Type-2 diabetes in human muscle. Importantly, in comparison to existing methods, ΔFBA gives a more accurate prediction of flux differences.

DeltaNeTS+: elucidating the mechanism of drugs and diseases using gene expression and transcriptional regulatory networks.

BACKGROUND: Knowledge on the molecular targets of diseases and drugs is crucial for elucidating disease pathogenesis and mechanism of action of drugs, and for driving drug discovery and treatment formulation. In this regard, high-throughput gene transcriptional profiling has become a leading technology, generating whole-genome data on the transcriptional alterations caused by diseases or drug compounds. However, identifying direct gene targets, especially in the background of indirect (downstream) effects, based on differential gene expressions is difficult due to the complexity of gene regulatory network governing the gene transcriptional processes. RESULTS: In this work, we developed a network analysis method, called DeltaNeTS+, for inferring direct gene targets of drugs and diseases from gene transcriptional profiles. DeltaNeTS+ uses a gene regulatory network model to identify direct perturbations to the transcription of genes using gene expression data. Importantly, DeltaNeTS+ is able to combine both steady-state and time-course expression profiles, as well as leverage information on the gene network structure. We demonstrated the power of DeltaNeTS+ in predicting gene targets using gene expression data in complex organisms, including Caenorhabditis elegans and human cell lines (T-cell and Calu-3). More specifically, in an application to time-course gene expression profiles of influenza A H1N1 (swine flu) and H5N1 (avian flu) infection, DeltaNeTS+ shed light on the key differences of dynamic cellular perturbations caused by the two influenza strains. CONCLUSION: DeltaNeTS+ is a powerful network analysis tool for inferring gene targets from gene expression profiles. As demonstrated in the case studies, by incorporating available information on gene network structure, DeltaNeTS+ produces accurate predictions of direct gene targets from a small sample size (~ 10 s). Integrating static and dynamic expression data with transcriptional network structure extracted from genomic information, as enabled by DeltaNeTS+, is crucial toward personalized medicine, where treatments can be tailored to individual patients. DeltaNeTS+ can be freely downloaded from http://www.github.com/cabsel/deltanetsplus .

Stochastic and Age-Dependent Proteostasis Decline Underlies Heterogeneity in Heat-Shock Response Dynamics.

Significant non-genetic stochastic factors affect aging, causing lifespan differences among individuals, even those sharing the same genetic and environmental background. In Caenorhabditis elegans, differences in heat-shock response (HSR) are predictive of lifespan. However, factors contributing to the heterogeneity of HSR are still not fully elucidated. Here, the authors characterized HSR dynamics in isogenic C. elegans expressing GFP reporter for hsp-16.2 for identifying the key contributors of HSR heterogeneity. Specifically, microfluidic devices that enable cross-sectional and longitudinal measurements of HSR dynamics in C. elegans at different scales are developed: in populations, within individuals, and in embryos. The authors adapted a mathematical model of HSR to single C. elegans and identified model parameters associated with proteostasis-maintenance of protein homeostasis-more specifically, protein turnover, as the major drivers of heterogeneity in HSR dynamics. It is verified that individuals with enhanced proteostasis fidelity in early adulthood live longer. The model-based comparative analysis of protein turnover in day-1 and day-2 adult C. elegans revealed a stochastic-onset of age-related proteostasis decline that increases the heterogeneity of HSR capacity. Finally, the analysis of C. elegans embryos showed higher HSR and proteostasis capacity than young adults and established transgenerational contribution to HSR heterogeneity that depends on maternal age.

A systems-based framework to computationally describe putative transcription factors and signaling pathways regulating glycan biosynthesis.

Glycosylation is a common posttranslational modification, and glycan biosynthesis is regulated by a set of glycogenes. The role of transcription factors (TFs) in regulating the glycogenes and related glycosylation pathways is largely unknown. In this work, we performed data mining of TF-glycogene relationships from the Cistrome Cancer database (DB), which integrates chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA-Seq data to constitute regulatory relationships. In total, we observed 22,654 potentially significant TF-glycogene relationships, which include interactions involving 526 unique TFs and 341 glycogenes that span 29 the Cancer Genome Atlas (TCGA) cancer types. Here, TF-glycogene interactions appeared in clusters or so-called communities, suggesting that changes in single TF expression during both health and disease may affect multiple carbohydrate structures. Upon applying the Fisher's exact test along with glycogene pathway classification, we identified TFs that may specifically regulate the biosynthesis of individual glycan types. Integration with Reactome DB knowledge provided an avenue to relate cell-signaling pathways to TFs and cellular glycosylation state. Whereas analysis results are presented for all 29 cancer types, specific focus is placed on human luminal and basal breast cancer disease progression. Overall, the article presents a computational approach to describe TF-glycogene relationships, the starting point for experimental system-wide validation.

Network perturbation analysis of gene transcriptional profiles reveals protein targets and mechanism of action of drugs and influenza A viral infection.

Genome-wide transcriptional profiling provides a global view of cellular state and how this state changes under different treatments (e.g. drugs) or conditions (e.g. healthy and diseased). Here, we present ProTINA (Protein Target Inference by Network Analysis), a network perturbation analysis method for inferring protein targets of compounds from gene transcriptional profiles. ProTINA uses a dynamic model of the cell-type specific protein-gene transcriptional regulation to infer network perturbations from steady state and time-series differential gene expression profiles. A candidate protein target is scored based on the gene network's dysregulation, including enhancement and attenuation of transcriptional regulatory activity of the protein on its downstream genes, caused by drug treatments. For benchmark datasets from three drug treatment studies, ProTINA was able to provide highly accurate protein target predictions and to reveal the mechanism of action of compounds with high sensitivity and specificity. Further, an application of ProTINA to gene expression profiles of influenza A viral infection led to new insights of the early events in the infection.

SINCERITIES: inferring gene regulatory networks from time-stamped single cell transcriptional expression profiles.

MOTIVATION: Single cell transcriptional profiling opens up a new avenue in studying the functional role of cell-to-cell variability in physiological processes. The analysis of single cell expression profiles creates new challenges due to the distributive nature of the data and the stochastic dynamics of gene transcription process. The reconstruction of gene regulatory networks (GRNs) using single cell transcriptional profiles is particularly challenging, especially when directed gene-gene relationships are desired. RESULTS: We developed SINCERITIES (SINgle CEll Regularized Inference using TIme-stamped Expression profileS) for the inference of GRNs from single cell transcriptional profiles. We focused on time-stamped cross-sectional expression data, commonly generated from transcriptional profiling of single cells collected at multiple time points after cell stimulation. SINCERITIES recovers directed regulatory relationships among genes by employing regularized linear regression (ridge regression), using temporal changes in the distributions of gene expressions. Meanwhile, the modes of the gene regulations (activation and repression) come from partial correlation analyses between pairs of genes. We demonstrated the efficacy of SINCERITIES in inferring GRNs using in silico time-stamped single cell expression data and single cell transcriptional profiles of THP-1 monocytic human leukemia cells. The case studies showed that SINCERITIES could provide accurate GRN predictions, significantly better than other GRN inference algorithms such as TSNI, GENIE3 and JUMP3. Moreover, SINCERITIES has a low computational complexity and is amenable to problems of extremely large dimensionality. Finally, an application of SINCERITIES to single cell expression data of T2EC chicken erythrocytes pointed to BATF as a candidate novel regulator of erythroid development. AVAILABILITY AND IMPLEMENTATION: MATLAB and R version of SINCERITIES are freely available from the following websites: http://www.cabsel.ethz.ch/tools/sincerities.html and https://github.com/CABSEL/SINCERITIES. The single cell THP-1 and T2EC transcriptional profiles are available from the original publications (Kouno et al., 2013; Richard et al., 2016). The in silico single cell data are available on SINCERITIES websites. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Single-Cell-Based Analysis Highlights a Surge in Cell-to-Cell Molecular Variability Preceding Irreversible Commitment in a Differentiation Process.

In some recent studies, a view emerged that stochastic dynamics governing the switching of cells from one differentiation state to another could be characterized by a peak in gene expression variability at the point of fate commitment. We have tested this hypothesis at the single-cell level by analyzing primary chicken erythroid progenitors through their differentiation process and measuring the expression of selected genes at six sequential time-points after induction of differentiation. In contrast to population-based expression data, single-cell gene expression data revealed a high cell-to-cell variability, which was masked by averaging. We were able to show that the correlation network was a very dynamical entity and that a subgroup of genes tend to follow the predictions from the dynamical network biomarker (DNB) theory. In addition, we also identified a small group of functionally related genes encoding proteins involved in sterol synthesis that could act as the initial drivers of the differentiation. In order to assess quantitatively the cell-to-cell variability in gene expression and its evolution in time, we used Shannon entropy as a measure of the heterogeneity. Entropy values showed a significant increase in the first 8 h of the differentiation process, reaching a peak between 8 and 24 h, before decreasing to significantly lower values. Moreover, we observed that the previous point of maximum entropy precedes two paramount key points: an irreversible commitment to differentiation between 24 and 48 h followed by a significant increase in cell size variability at 48 h. In conclusion, when analyzed at the single cell level, the differentiation process looks very different from its classical population average view. New observables (like entropy) can be computed, the behavior of which is fully compatible with the idea that differentiation is not a "simple" program that all cells execute identically but results from the dynamical behavior of the underlying molecular network.

Acoustic Compressibility of Caenorhabditis elegans.

The acoustic compressibility of Caenorhabditis elegans is a necessary parameter for further understanding the underlying physics of acoustic manipulation techniques of this widely used model organism in biological sciences. In this work, numerical simulations were combined with experimental trajectory velocimetry of L1 C. elegans larvae to estimate the acoustic compressibility of C. elegans. A method based on bulk acoustic wave acoustophoresis was used for trajectory velocimetry experiments in a microfluidic channel. The model-based data analysis took into account the different sizes and shapes of L1 C. elegans larvae (255 ± 26 µm in length and 15 ± 2 µm in diameter). Moreover, the top and bottom walls of the microfluidic channel were considered in the hydrodynamic drag coefficient calculations, for both the C. elegans and the calibration particles. The hydrodynamic interaction between the specimen and the channel walls was further minimized by acoustically levitating the C. elegans and the particles to the middle of the measurement channel. Our data suggest an acoustic compressibility κCe of 430 TPa-1 with an uncertainty range of ±20 TPa-1 for C. elegans, a much lower value than what was previously reported for adult C. elegans using static methods. Our estimated compressibility is consistent with the relative volume fraction of lipids and proteins that would mainly make up for the body of C. elegans. This work is a departing point for practical engineering and design criteria for integrated acoustofluidic devices for biological applications.

Inferring gene targets of drugs and chemical compounds from gene expression profiles.

MOTIVATION: Finding genes which are directly perturbed or targeted by drugs is of great interest and importance in drug discovery. Several network filtering methods have been created to predict the gene targets of drugs from gene expression data based on an ordinary differential equation model of the gene regulatory network (GRN). A critical step in these methods involves inferring the GRN from the expression data, which is a very challenging problem on its own. In addition, existing network filtering methods require computationally intensive parameter tuning or expression data from experiments with known genetic perturbations or both. RESULTS: We developed a method called DeltaNet for the identification of drug targets from gene expression data. Here, the gene target predictions were directly inferred from the data without a separate step of GRN inference. DeltaNet formulation led to solving an underdetermined linear regression problem, for which we employed least angle regression (DeltaNet-LAR) or LASSO regularization (DeltaNet-LASSO). The predictions using DeltaNet for expression data of Escherichia coli, yeast, fruit fly and human were significantly more accurate than those using network filtering methods, namely mode of action by network identification (MNI) and sparse simultaneous equation model (SSEM). Furthermore, DeltaNet using LAR did not require any parameter tuning and could provide computational speed-up over existing methods. CONCLUSION: DeltaNet is a robust and numerically efficient tool for identifying gene perturbations from gene expression data. Importantly, the method requires little to no expert supervision, while providing accurate gene target predictions. AVAILABILITY AND IMPLEMENTATION: DeltaNet is available on http://www.cabsel.ethz.ch/tools/DeltaNet CONTACT: rudi.gunawan@chem.ethz.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Optimal design of gene knockout experiments for gene regulatory network inference.

MOTIVATION: We addressed the problem of inferring gene regulatory network (GRN) from gene expression data of knockout (KO) experiments. This inference is known to be underdetermined and the GRN is not identifiable from data. Past studies have shown that suboptimal design of experiments (DOE) contributes significantly to the identifiability issue of biological networks, including GRNs. However, optimizing DOE has received much less attention than developing methods for GRN inference. RESULTS: We developed REDuction of UnCertain Edges (REDUCE) algorithm for finding the optimal gene KO experiment for inferring directed graphs (digraphs) of GRNs. REDUCE employed ensemble inference to define uncertain gene interactions that could not be verified by prior data. The optimal experiment corresponds to the maximum number of uncertain interactions that could be verified by the resulting data. For this purpose, we introduced the concept of edge separatoid which gave a list of nodes (genes) that upon their removal would allow the verification of a particular gene interaction. Finally, we proposed a procedure that iterates over performing KO experiments, ensemble update and optimal DOE. The case studies including the inference of Escherichia coli GRN and DREAM 4 100-gene GRNs, demonstrated the efficacy of the iterative GRN inference. In comparison to systematic KOs, REDUCE could provide much higher information return per gene KO experiment and consequently more accurate GRN estimates. CONCLUSIONS: REDUCE represents an enabling tool for tackling the underdetermined GRN inference. Along with advances in gene deletion and automation technology, the iterative procedure brings an efficient and fully automated GRN inference closer to reality. AVAILABILITY AND IMPLEMENTATION: MATLAB and Python scripts of REDUCE are available on www.cabsel.ethz.ch/tools/REDUCE CONTACT: rudi.gunawan@chem.ethz.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Glycosylation flux analysis reveals dynamic changes of intracellular glycosylation flux distribution in Chinese hamster ovary fed-batch cultures.

N-linked glycosylation of proteins has both functional and structural significance. Importantly, the glycan structure of a therapeutic protein influences its efficacy, pharmacokinetics, pharmacodynamics and immunogenicity. In this work, we developed glycosylation flux analysis (GFA) for predicting intracellular production and consumption rates (fluxes) of glycoforms, and applied this analysis to CHO fed-batch immunoglobulin G (IgG) production using two different media compositions, with and without additional manganese feeding. The GFA is based on a constraint-based modeling of the glycosylation network, employing a pseudo steady state assumption. While the glycosylation fluxes in the network are balanced at each time point, the GFA allows the fluxes to vary with time by way of two scaling factors: (1) an enzyme-specific factor that captures the temporal changes among glycosylation reactions catalysed by the same enzyme, and (2) the cell specific productivity factor that accounts for the dynamic changes in the IgG production rate. The GFA of the CHO fed-batch cultivations showed that regardless of the media composition, galactosylation fluxes decreased with the cultivation time more significantly than the other glycosylation reactions. Furthermore, the GFA showed that the addition of Mn, a cofactor of galactosyltransferase, has the effect of increasing the galactosylation fluxes but only during the beginning of the cultivation period. The results thus demonstrated the power of the GFA in delineating the dynamic alterations of the glycosylation fluxes by local (enzyme-specific) and global (cell specific productivity) factors.

Are mutagenic non D-loop direct repeat motifs in mitochondrial DNA under a negative selection pressure?

Non D-loop direct repeats (DRs) in mitochondrial DNA (mtDNA) have been commonly implicated in the mutagenesis of mtDNA deletions associated with neuromuscular disease and ageing. Further, these DRs have been hypothesized to put a constraint on the lifespan of mammals and are under a negative selection pressure. Using a compendium of 294 mammalian mtDNA, we re-examined the relationship between species lifespan and the mutagenicity of such DRs. Contradicting the prevailing hypotheses, we found no significant evidence that long-lived mammals possess fewer mutagenic DRs than short-lived mammals. By comparing DR counts in human mtDNA with those in selectively randomized sequences, we also showed that the number of DRs in human mtDNA is primarily determined by global mtDNA properties, such as the bias in synonymous codon usage (SCU) and nucleotide composition. We found that SCU bias in mtDNA positively correlates with DR counts, where repeated usage of a subset of codons leads to more frequent DR occurrences. While bias in SCU and nucleotide composition has been attributed to nucleotide mutational bias, mammalian mtDNA still exhibit higher SCU bias and DR counts than expected from such mutational bias, suggesting a lack of negative selection against non D-loop DRs.

TRaCE+: Ensemble inference of gene regulatory networks from transcriptional expression profiles of gene knock-out experiments.

BACKGROUND: The inference of gene regulatory networks (GRNs) from transcriptional expression profiles is challenging, predominantly due to its underdetermined nature. One important consequence of underdetermination is the existence of many possible solutions to this inference. Our previously proposed ensemble inference algorithm TRaCE addressed this issue by inferring an ensemble of network directed graphs (digraphs) using differential gene expressions from gene knock-out (KO) experiments. However, TRaCE could not deal with the mode of the transcriptional regulations (activation or repression), an important feature of GRNs. RESULTS: In this work, we developed a new algorithm called TRaCE+ for the inference of an ensemble of signed GRN digraphs from transcriptional expression data of gene KO experiments. The sign of the edges indicates whether the regulation is an activation (positive) or a repression (negative). TRaCE+ generates the upper and lower bounds of the ensemble, which define uncertain regulatory interactions that could not be verified by the data. As demonstrated in the case studies using Escherichia coli GRN and 100-gene gold-standard GRNs from DREAM 4 network inference challenge, by accounting for regulatory signs, TRaCE+ could extract more information from the KO data than TRaCE, leading to fewer uncertain edges. Importantly, iterating TRaCE+ with an optimal design of gene KOs could resolve the underdetermined issue of GRN inference in much fewer KO experiments than using TRaCE. CONCLUSIONS: TRaCE+ expands the applications of ensemble GRN inference strategy by accounting for the mode of the gene regulatory interactions. In comparison to TRaCE, TRaCE+ enables a better utilization of gene KO data, thereby reducing the cost of tackling underdetermined GRN inference. TRaCE+ subroutines for MATLAB are freely available at the following website: http://www.cabsel.ethz.ch/tools/trace.html .

REDEMPTION: reduced dimension ensemble modeling and parameter estimation.

UNLABELLED: Here, we present REDEMPTION ( RE: duced D: imension E: nsemble M: odeling and P: arameter estima TION: ), a toolbox for parameter estimation and ensemble modeling of ordinary differential equations (ODEs) using time-series data. For models with more reactions than measured species, a common scenario in biological modeling, the parameter estimation is formulated as a nested optimization problem based on incremental parameter estimation strategy. REDEMPTION also includes a tool for the identification of an ensemble of parameter combinations that provide satisfactory goodness-of-fit to the data. The functionalities of REDEMPTION are accessible through a MATLAB user interface (UI), as well as through programming script. For computational speed-up, REDEMPTION provides a numerical parallelization option using MATLAB Parallel Computing toolbox. AVAILABILITY AND IMPLEMENTATION: REDEMPTION can be downloaded from http://www.cabsel.ethz.ch/tools/redemption. CONTACT: rudi.gunawan@chem.ethz.ch.

Context-Dependent Role of Mitochondrial Fusion-Fission in Clonal Expansion of mtDNA Mutations.

The accumulation of mutant mitochondrial DNA (mtDNA) molecules in aged cells has been associated with mitochondrial dysfunction, age-related diseases and the ageing process itself. This accumulation has been shown to often occur clonally, where mutant mtDNA grow in number and overpopulate the wild-type mtDNA. However, the cell possesses quality control (QC) mechanisms that maintain mitochondrial function, in which dysfunctional mitochondria are isolated and removed by selective fusion and mitochondrial autophagy (mitophagy), respectively. The aim of this study is to elucidate the circumstances related to mitochondrial QC that allow the expansion of mutant mtDNA molecules. For the purpose of the study, we have developed a mathematical model of mitochondrial QC process by extending our previous validated model of mitochondrial turnover and fusion-fission. A global sensitivity analysis of the model suggested that the selectivity of mitophagy and fusion is the most critical QC parameter for clearing de novo mutant mtDNA molecules. We further simulated several scenarios involving perturbations of key QC parameters to gain a better understanding of their dynamic and synergistic interactions. Our model simulations showed that a higher frequency of mitochondrial fusion-fission can provide a faster clearance of mutant mtDNA, but only when mutant-rich mitochondria that are transiently created are efficiently prevented from re-fusing with other mitochondria and selectively removed. Otherwise, faster fusion-fission quickens the accumulation of mutant mtDNA. Finally, we used the insights gained from model simulations and analysis to propose a possible circumstance involving deterioration of mitochondrial QC that permits mutant mtDNA to expand with age.

Maximizing signal-to-noise ratio in the random mutation capture assay.

The 'Random Mutation Capture' assay allows for the sensitive quantitation of DNA mutations at extremely low mutation frequencies. This method is based on PCR detection of mutations that render the mutated target sequence resistant to restriction enzyme digestion. The original protocol prescribes an end-point dilution to about 0.1 mutant DNA molecules per PCR well, such that the mutation burden can be simply calculated by counting the number of amplified PCR wells. However, the statistical aspects associated with the single molecular nature of this protocol and several other molecular approaches relying on binary (on/off) output can significantly affect the quantification accuracy, and this issue has so far been ignored. The present work proposes a design of experiment (DoE) using statistical modeling and Monte Carlo simulations to obtain a statistically optimal sampling protocol, one that minimizes the coefficient of variance in the measurement estimates. Here, the DoE prescribed a dilution factor at about 1.6 mutant molecules per well. Theoretical results and experimental validation revealed an up to 10-fold improvement in the information obtained per PCR well, i.e. the optimal protocol achieves the same coefficient of variation using one-tenth the number of wells used in the original assay. Additionally, this optimization equally applies to any method that relies on binary detection of a small number of templates.

Cell and tissue-specific glycosylation pathways informed by single-cell transcriptomics.

While single cell studies have made significant impacts in various subfields of biology, they lag in the Glycosciences. To address this gap, we analyzed single-cell glycogene expressions in the Tabula Sapiens dataset of human tissues and cell types using a recent glycosylation-specific gene ontology (GlycoEnzOnto). At the median sequencing (count) depth, ~40-50 out of 400 glycogenes were detected in individual cells. Upon increasing the sequencing depth, the number of detectable glycogenes saturates at ~200 glycogenes, suggesting that the average human cell expresses about half of the glycogene repertoire. Hierarchies in glycogene and glycopathway expressions emerged from our analysis: nucleotide-sugar synthesis and transport exhibited the highest gene expressions, followed by genes for core enzymes, glycan modification and extensions, and finally terminal modifications. Interestingly, the same cell types showed variable glycopathway expressions based on their organ or tissue origin, suggesting nuanced cell- and tissue-specific glycosylation patterns. Probing deeper into the transcription factors (TFs) of glycogenes, we identified distinct groupings of TFs controlling different aspects of glycosylation: core biosynthesis, terminal modifications, etc. We present webtools to explore the interconnections across glycogenes, glycopathways, and TFs regulating glycosylation in human cell/tissue types. Overall, the study presents an overview of glycosylation across multiple human organ systems.