How Far Are We from the Completion of the Human Protein Interactome Reconstruction?

After more than fifteen years from the first high-throughput experiments for human protein–protein interaction (PPI) detection, we are still wondering how close the completion of the genome-scale human PPI network reconstruction is, what needs to be further explored and whether the biological insights gained from the holistic investigation of the current network are valid and useful. The unique structure of PICKLE, a meta-database of the human experimentally determined direct PPI network developed by our group, presently covering ~80% of the UniProtKB/Swiss-Prot reviewed human complete proteome, enables the evaluation of the interactome expansion by comparing the successive PICKLE releases since 2013. We observe a gradual overall increase of 39%, 182%, and 67% in protein nodes, PPIs, and supporting references, respectively. Our results indicate that, in recent years, (a) the PPI addition rate has decreased, (b) the new PPIs are largely determined by high-throughput experiments and mainly concern existing protein nodes and (c), as we had predicted earlier, most of the newly added protein nodes have a low degree. These observations, combined with a largely overlapping k-core between PICKLE releases and a network density increase, imply that an almost complete picture of a structurally defined network has been reached. The comparative unsupervised application of two clustering algorithms indicated that exploring the full interactome topology can reveal the protein neighborhoods involved in closely related biological processes as transcriptional regulation, cell signaling and multiprotein complexes such as the connexon complex associated with cancers. A well-reconstructed human protein interactome is a powerful tool in network biology and medicine research forming the basis for multi-omic and dynamic analyses.

Three Reagents for in-Solution Enrichment of Ancient Human DNA at More than a Million SNPs

In-solution enrichment for hundreds of thousands of single nucleotide polymorphisms (SNPs) has been the source of >70% of all genome-scale ancient human DNA data published to date. This approach has made it possible to generate data for one to two orders of magnitude lower cost than random shotgun sequencing, making it economical to study ancient samples with low proportions of human DNA, and increasing the rate of conversion of sampled remains into working data thereby facilitating ethical stewardship of human remains. So far, nearly all ancient DNA data obtained using in-solution enrichment has been generated using a set of bait sequences targeting about 1.24 million SNPs (the 1240k reagent). These sequences were published in 2015, but synthesis of the reagent has been cost-effective for only a few laboratories. In 2021, two companies made available reagents that target the same core set of SNPs along with supplementary content. Here, we test the properties of the three reagents on a common set of 27 ancient DNA libraries across a range of richness of DNA content and percentages of human molecules. All three reagents are highly effective at enriching many hundreds of thousands of SNPs. For all three reagents and a wide range of conditions, one round of enrichment produces data that is as useful as two rounds when tens of millions of sequences are read out as is typical for such experiments. In our testing, the Twist Ancient DNA reagent produces the highest coverages, greatest uniformity on targeted positions, and almost no bias toward enriching one allele more than another relative to shotgun sequencing. Allelic bias in 1240k enrichment has made it challenging to carry out joint analysis of these data with shotgun data, creating a situation where the ancient DNA community has been publishing two important bodes of data that cannot easily be co-analyzed by population genetic methods. To address this challenge, we introduce a subset of hundreds of thousands of SNPs for which 1240k data can be effectively co-analyzed with all other major data types.

Machine-learning from Pseudomonas putida transcriptomes reveals its transcriptional regulatory network

Bacterial gene expression is orchestrated by numerous transcription factors (TFs). Elucidating how gene expression is regulated is fundamental to understanding bacterial physiology and engineering it for practical use. In this study, a machine-learning approach was applied to uncover the genome-scale transcriptional regulatory network (TRN) in Pseudomonas putida, an important organism for bioproduction. We performed independent component analysis of a compendium of 321 high-quality gene expression profiles, which were previously published or newly generated in this study. We identified 84 groups of independently modulated genes (iModulons) that explain 75.7% of the total variance in the compendium. With these iModulons, we (i) expand our understanding of the regulatory functions of 39 iModulon associated TFs (e.g., HexR, Zur) by systematic comparison with 1,993 previously reported TF-gene interactions; (ii) outline transcriptional changes after the transition from the exponential growth to stationary phases; (iii) capture group of genes required for utilizing diverse carbon sources and increased stationary response with slower growth rates; (iv) unveil multiple evolutionary strategies of transcriptome reallocation to achieve fast growth rates; and (v) define an osmotic stimulon, which includes the Type VI secretion system, as coordination of multiple iModulon activity changes. Taken together, this study provides the first quantitative genome-scale TRN for P. putida and a basis for a comprehensive understanding of its complex transcriptome changes in a variety of physiological states.

Genome-Scale Metabolic Modelling of Lifestyle Changes in Rhizobium leguminosarum

Rhizobia are soil bacteria that induce nodule formation on plant roots and differentiate into nitrogen-fixing bacteroids. A detailed understanding of this complex symbiosis is essential for advancing ongoing efforts to engineer novel symbioses with cereal crops for sustainable agriculture.

Comparison of Different Label-Free Techniques for the Semi-Absolute Quantification of Protein Abundance

In proteomics, it is essential to quantify proteins in absolute terms if we wish to compare results among studies and integrate high-throughput biological data into genome-scale metabolic models. While labeling target peptides with stable isotopes allow protein abundance to be accurately quantified, the utility of this technique is constrained by the low number of quantifiable proteins that it yields. Recently, label-free shotgun proteomics has become the “gold standard” for carrying out global assessments of biological samples containing thousands of proteins. However, this tool must be further improved if we wish to accurately quantify absolute levels of proteins. Here, we used different label-free quantification techniques to estimate absolute protein abundance in the model yeast Saccharomyces cerevisiae. More specifically, we evaluated the performance of seven different quantification methods, based either on spectral counting (SC) or extracted-ion chromatogram (XIC), which were applied to samples from five different proteome backgrounds. We also compared the accuracy and reproducibility of two strategies for transforming relative abundance into absolute abundance: a UPS2-based strategy and the total protein approach (TPA). This study mentions technical challenges related to UPS2 use and proposes ways of addressing them, including utilizing a smaller, more highly optimized amount of UPS2. Overall, three SC-based methods (PAI, SAF, and NSAF) yielded the best results because they struck a good balance between experimental performance and protein quantification.

Whole-genome-scale identification of novel non-protein-coding RNAs controlling cell proliferation and survival through a functional forward genetics strategy

Identification of cell fate-controlling lncRNAs is essential to our understanding of molecular cell biology. Here we present a human genome-scale forward-genetics approach for the identification of lncRNAs based on gene function. This approach can identify genes that play a causal role, and immediately distinguish them from those that are differentially expressed but do not affect cell function. Our genome-scale library plus next-generation-sequencing and bioinformatic approach, radically upscales the breadth and rate of functional ncRNA discovery. Human gDNA was digested to produce a lentiviral expression library containing inserts in both sense and anti-sense orientation. The library was used to transduce human Jurkat T-leukaemic cells. Cell populations were selected using continuous culture ± anti-FAS IgM, and sequencing used to identify sequences controlling cell proliferation. This strategy resulted in the identification of thousands of new sequences based solely on their function including many ncRNAs previously identified as being able to modulate cell survival or to act as key cancer regulators such as AC084816.1*, AC097103.2, AC087473.1, CASC15*, DLEU1*, ENTPD1-AS1*, HULC*, MIRLET7BHG*, PCAT-1, SChLAP1, and TP53TG1. Independent validation confirmed 4 out of 5 sequences that were identified by this strategy, conferred a striking resistance to anti-FAS IgM-induced apoptosis.

A Computational Framework to Identify Metabolic Engineering Strategies for the Co-Production of Metabolites

Microbial production of chemicals is a more sustainable alternative to traditional chemical processes. However, the shift to bioprocess is usually accompanied by a drop in economic feasibility. Co-production of more than one chemical can improve the economy of bioprocesses, enhance carbon utilization and also ensure better exploitation of resources. While a number of tools exist for in silico metabolic engineering, there is a dearth of computational tools that can co-optimize the production of multiple metabolites. In this work, we propose co-FSEOF (co-production using Flux Scanning based on Enforced Objective Flux), an algorithm designed to identify intervention strategies to co-optimize the production of a set of metabolites. Co-FSEOF can be used to identify all pairs of products that can be co-optimized with ease using a single intervention. Beyond this, it can also identify higher-order intervention strategies for a given set of metabolites. We have employed this tool on the genome-scale metabolic models of Escherichia coli and Saccharomyces cerevisiae, and identified intervention targets that can co-optimize the production of pairs of metabolites under both aerobic and anaerobic conditions. Anaerobic conditions were found to support the co-production of a higher number of metabolites when compared to aerobic conditions in both organisms. The proposed computational framework will enhance the ease of study of metabolite co-production and thereby aid the design of better bioprocesses.

Inspecting the Solution Space of Genome-Scale Metabolic Models

Genome-scale metabolic models are frequently used in computational biology. They offer an integrative view on the metabolic network of an organism without the need to know kinetic information in detail. However, the huge solution space which comes with the analysis of genome-scale models by using, e.g., Flux Balance Analysis (FBA) poses a problem, since it is hard to thoroughly investigate and often only an arbitrarily selected individual flux distribution is discussed as an outcome of FBA. Here, we introduce a new approach to inspect the solution space and we compare it with other approaches, namely Flux Variability Analysis (FVA) and CoPE-FBA, using several different genome-scale models of lactic acid bacteria. We examine the extent to which different types of experimental data limit the solution space and how the robustness of the system increases as a result. We find that our new approach to inspect the solution space is a good complementary method that offers additional insights into the variance of biological phenotypes and can help to prevent wrong conclusions in the analysis of FBA results.