Genome-wide quantification of copy-number aberration impact on gene expression in ovarian high-grade serous carcinoma.

Copy-number alterations (CNAs) are a hallmark of cancer and can regulate cancer cell states via altered gene expression values. Herein, we have developed a copy-number impact (CNI) analysis method that quantifies the degree to which a gene expression value is impacted by CNAs and leveraged this analysis at the pathway level. Our results show that a high CNA is not necessarily reflected at the gene expression level, and our method is capable of detecting genes and pathways whose activity is strongly influenced by CNAs. Furthermore, the CNI analysis enables unbiased categorization of CNA categories, such as deletions and amplifications. We identified six CNI-driven pathways associated with poor treatment response in ovarian high-grade serous carcinoma (HGSC), which we found to be the most CNA-driven cancer across 14 cancer types. The key driver in most of these pathways was amplified wild-type KRAS, which we validated functionally using CRISPR modulation. Our results suggest that wild-type KRAS amplification is a driver of chemotherapy resistance in HGSC and may serve as a potential treatment target.

The transcription factor network of E. coli steers global responses to shifts in RNAP concentration.

The robustness and sensitivity of gene networks to environmental changes is critical for cell survival. How gene networks produce specific, chronologically ordered responses to genome-wide perturbations, while robustly maintaining homeostasis, remains an open question. We analysed if short- and mid-term genome-wide responses to shifts in RNA polymerase (RNAP) concentration are influenced by the known topology and logic of the transcription factor network (TFN) of Escherichia coli. We found that, at the gene cohort level, the magnitude of the single-gene, mid-term transcriptional responses to changes in RNAP concentration can be explained by the absolute difference between the gene's numbers of activating and repressing input transcription factors (TFs). Interestingly, this difference is strongly positively correlated with the number of input TFs of the gene. Meanwhile, short-term responses showed only weak influence from the TFN. Our results suggest that the global topological traits of the TFN of E. coli shape which gene cohorts respond to genome-wide stresses.

Network-guided identification of cancer-selective combinatorial therapies in ovarian cancer.

Each patient's cancer consists of multiple cell subpopulations that are inherently heterogeneous and may develop differing phenotypes such as drug sensitivity or resistance. A personalized treatment regimen should therefore target multiple oncoproteins in the cancer cell populations that are driving the treatment resistance or disease progression in a given patient to provide maximal therapeutic effect, while avoiding severe co-inhibition of non-malignant cells that would lead to toxic side effects. To address the intra- and inter-tumoral heterogeneity when designing combinatorial treatment regimens for cancer patients, we have implemented a machine learning-based platform to guide identification of safe and effective combinatorial treatments that selectively inhibit cancer-related dysfunctions or resistance mechanisms in individual patients. In this case study, we show how the platform enables prediction of cancer-selective drug combinations for patients with high-grade serous ovarian cancer using single-cell imaging cytometry drug response assay, combined with genome-wide transcriptomic and genetic profiles. The platform makes use of drug-target interaction networks to prioritize those combinations that warrant further preclinical testing in scarce patient-derived primary cells. During the case study in ovarian cancer patients, we investigated (i) the relative performance of various ensemble learning algorithms for drug response prediction, (ii) the use of matched single-cell RNA-sequencing data to deconvolute cell population-specific transcriptome profiles from bulk RNA-seq data, (iii) and whether multi-patient or patient-specific predictive models lead to better predictive accuracy. The general platform and the comparison results are expected to become useful for future studies that use similar predictive approaches also in other cancer types.

POIBM: batch correction of heterogeneous RNA-seq datasets through latent sample matching.

MOTIVATION: RNA sequencing and other high-throughput technologies are essential in understanding complex diseases, such as cancers, but are susceptible to technical factors manifesting as patterns in the measurements. These batch patterns hinder the discovery of biologically relevant patterns. Unbiased batch effect correction in heterogeneous populations currently requires special experimental designs or phenotypic labels, which are not readily available for patient samples in existing datasets. RESULTS: We present POIBM, an RNA-seq batch correction method, which learns virtual reference samples directly from the data. We use a breast cancer cell line dataset to show that POIBM exceeds or matches the performance of previous methods, while being blind to the phenotypes. Further, we analyze The Cancer Genome Atlas RNA-seq data to show that batch effects plague many cancer types; POIBM effectively discovers the true replicates in stomach adenocarcinoma; and integrating the corrected data in endometrial carcinoma improves cancer subtyping. AVAILABILITY AND IMPLEMENTATION: https://bitbucket.org/anthakki/poibm/ (archived at https://doi.org/10.5281/zenodo.6122436). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

QuantISH: RNA in situ hybridization image analysis framework for quantifying cell type-specific target RNA expression and variability.

RNA in situ hybridization (RNA-ISH) is a powerful spatial transcriptomics technology to characterize target RNA abundance and localization in individual cells. This allows analysis of tumor heterogeneity and expression localization, which are not readily obtainable through transcriptomic data analysis. RNA-ISH experiments produce large amounts of data and there is a need for automated analysis methods. Here we present QuantISH, a comprehensive open-source RNA-ISH image analysis pipeline that quantifies marker expressions in individual carcinoma, immune, and stromal cells on chromogenic or fluorescent in situ hybridization images. QuantISH is designed to be modular and can be adapted to various image and sample types and staining protocols. We show that in chromogenic RNA in situ hybridization images of high-grade serous carcinoma (HGSC) QuantISH cancer cell classification has high precision, and signal expression quantification is in line with visual assessment. We further demonstrate the power of QuantISH by showing that CCNE1 average expression and DDIT3 expression variability, as captured by the variability factor developed herein, act as candidate biomarkers in HGSC. Altogether, our results demonstrate that QuantISH can quantify RNA expression levels and their variability in carcinoma cells, and thus paves the way to utilize RNA-ISH technology.

PRISM: recovering cell-type-specific expression profiles from individual composite RNA-seq samples.

MOTIVATION: A major challenge in analyzing cancer patient transcriptomes is that the tumors are inherently heterogeneous and evolving. We analyzed 214 bulk RNA samples of a longitudinal, prospective ovarian cancer cohort and found that the sample composition changes systematically due to chemotherapy and between the anatomical sites, preventing direct comparison of treatment-naive and treated samples. RESULTS: To overcome this, we developed PRISM, a latent statistical framework to simultaneously extract the sample composition and cell-type-specific whole-transcriptome profiles adapted to each individual sample. Our results indicate that the PRISM-derived composition-free transcriptomic profiles and signatures derived from them predict the patient response better than the composite raw bulk data. We validated our findings in independent ovarian cancer and melanoma cohorts, and verified that PRISM accurately estimates the composition and cell-type-specific expression through whole-genome sequencing and RNA in situ hybridization experiments. AVAILABILITYAND IMPLEMENTATION: https://bitbucket.org/anthakki/prism. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Agile workflow for interactive analysis of mass cytometry data.

MOTIVATION: Single-cell proteomics technologies, such as mass cytometry, have enabled characterization of cell-to-cell variation and cell populations at a single-cell resolution. These large amounts of data, require dedicated, interactive tools for translating the data into knowledge. RESULTS: We present a comprehensive, interactive method called Cyto to streamline analysis of large-scale cytometry data. Cyto is a workflow-based open-source solution that automates the use of state-of-the-art single-cell analysis methods with interactive visualization. We show the utility of Cyto by applying it to mass cytometry data from peripheral blood and high-grade serous ovarian cancer (HGSOC) samples. Our results show that Cyto is able to reliably capture the immune cell sub-populations from peripheral blood and cellular compositions of unique immune- and cancer cell subpopulations in HGSOC tumor and ascites samples. AVAILABILITYAND IMPLEMENTATION: The method is available as a Docker container at https://hub.docker.com/r/anduril/cyto and the user guide and source code are available at https://bitbucket.org/anduril-dev/cyto. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Cardiovascular disease in people born to unmarried mothers in two historical periods: The Helsinki Birth Cohort Study 1934-1944.

Aims:Socio-economic conditions in early life are important contributors to cardiovascular disease - the leading cause of mortality globally - in later life. We studied coronary heart disease (CHD) and stroke in adulthood among people born out of wedlock in two historical periods: before and during World War II in Finland. Methods: We compared offspring born out of wedlock before (1934-1939) and during (1940-1944) World War II with the offspring of married mothers in the Helsinki Birth Cohort Study. The war affected the position of unmarried mothers in society. We followed the study subjects from 1971 to 2014 and identified deaths and hospital admissions from CHD and stroke. Data were analysed using a Cox regression, adjusting for other childhood and adulthood socio-economic circumstances. Results: The rate of out-of-wedlock births was 240/4052 (5.9%) before World War II and 397/9197 (4.3%) during World War II. Among those born before World War II, out-of-wedlock birth was associated with an increased risk of stroke (hazard ratio (HR)=1.44; 95% confidence interval (CI) 1.00-2.07) and CHD (HR=1.37; 95% CI 1.02-1.86). Among those born out of wedlock during World War II, the risks of stroke (HR=0.89; 95% CI 0.58-1.36) and CHD (HR=0.70; 95% CI 0.48=1.03) were similar to those observed for the offspring of married mothers. The p-values for interaction of unmarried×World War II were (p=0.015) for stroke and (p=0.003) for CHD. Conclusions: In a society in which marriage is normative, being born out of wedlock is an important predictor of lifelong health disadvantage. However, this may change rapidly when societal circumstances change, such as during a war.

qSNE: quadratic rate t-SNE optimizer with automatic parameter tuning for large datasets.

MOTIVATION: Non-parametric dimensionality reduction techniques, such as t-distributed stochastic neighbor embedding (t-SNE), are the most frequently used methods in the exploratory analysis of single-cell datasets. Current implementations scale poorly to massive datasets and often require downsampling or interpolative approximations, which can leave less-frequent populations undiscovered and much information unexploited. RESULTS: We implemented a fast t-SNE package, qSNE, which uses a quasi-Newton optimizer, allowing quadratic convergence rate and automatic perplexity (level of detail) optimizer. Our results show that these improvements make qSNE significantly faster than regular t-SNE packages and enables full analysis of large datasets, such as mass cytometry data, without downsampling. AVAILABILITY AND IMPLEMENTATION: Source code and documentation are openly available at https://bitbucket.org/anthakki/qsne/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Anduril 2: upgraded large-scale data integration framework.

SUMMARY: Anduril is an analysis and integration framework that facilitates the design, use, parallelization and reproducibility of bioinformatics workflows. Anduril has been upgraded to use Scala for pipeline construction, which simplifies software maintenance, and facilitates design of complex pipelines. Additionally, Anduril's bioinformatics repository has been expanded with multiple components, and tutorial pipelines, for next-generation sequencing data analysis. AVAILABILITYAND IMPLEMENTATION: Freely available at http://anduril.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SCIP: a single-cell image processor toolbox.

Summary: Each cell is a phenotypically unique individual that is influenced by internal and external processes, operating in parallel. To characterize the dynamics of cellular processes one needs to observe many individual cells from multiple points of view and over time, so as to identify commonalities and variability. With this aim, we engineered a software, 'SCIP', to analyze multi-modal, multi-process, time-lapse microscopy morphological and functional images. SCIP is capable of automatic and/or manually corrected segmentation of cells and lineages, automatic alignment of different microscopy channels, as well as detect, count and characterize fluorescent spots (such as RNA tagged by MS2-GFP), nucleoids, Z rings, Min system, inclusion bodies, undefined structures, etc. The results can be exported into *mat files and all results can be jointly analyzed, to allow studying not only each feature and process individually, but also find potential relationships. While we exemplify its use on Escherichia coli, many of its functionalities are expected to be of use in analyzing other prokaryotes and eukaryotic cells as well. We expect SCIP to facilitate the finding of relationships between cellular processes, from small-scale (e.g. gene expression) to large-scale (e.g. cell division), in single cells and cell lineages. Availability and implementation: http://www.ca3-uninova.org/project_scip. Supplementary information: Supplementary data are available at Bioinformatics online.

Identifying differentially methylated sites in samples with varying tumor purity.

Motivation: DNA methylation aberrations are common in many cancer types. A major challenge hindering comparison of patient-derived samples is that they comprise of heterogeneous collection of cancer and microenvironment cells. We present a computational method that allows comparing cancer methylomes in two or more heterogeneous tumor samples featuring differing, unknown fraction of cancer cells. The method is unique in that it allows comparison also in the absence of normal cell control samples and without prior tumor purity estimates, as these are often unavailable or unreliable in clinical samples. Results: We use simulations and next-generation methylome, RNA and whole-genome sequencing data from two cancer types to demonstrate that the method is accurate and outperforms alternatives. The results show that our method adapts well to various cancer types and to a wide range of tumor content, and works robustly without a control or with controls derived from various sources. Availability and implementation: The method is freely available at https://bitbucket.org/anthakki/dmml. Supplementary information: Supplementary data are available at Bioinformatics online.

Extraction of hazardous metals from green liquor dregs by ethylenediaminetetraacetic acid.

Green liquor dregs are the major inorganic solid side stream of kraft pulp mills which contain environmentally hazardous metals. The presence of hazardous metals in this industrial residue brings statutory limits for its landfilling, although they are not easily mobilized from the solid phase. In this study, the chelating agent ethylenediaminetetraacetic acid (EDTA) is utilized to extract hazardous metals such as Cd, Pb and Zn from green liquor dregs. Furthermore, the influence of EDTA on the removal of Ca as the main mineral nutrient present in the green liquor dregs is studied. The effect of parameters such as EDTA dosage, L/S ratio and contact time on the removal rate of the elements is investigated. In addition, the experimental data are fitted to the Elovich model and the pseudo-first-order model to describe the desorption kinetics. The results show that 59 wt% of Cd, 13 wt% of Co, 62 wt% of Cu, 3 wt% of Mn, 12 wt% of Ni, 43 wt% of Pb, 16 wt% of Zn, and less than 1 wt% of Ca were extracted from green liquor dregs with EDTA dosage of 0.035 gEDTA salt/gdregs and the L/S ratio of 6.25 ml/g. The current study opens up new possibilities to use the green liquor dregs for improving the soil fertility instead of landfilling.

Increased cytoplasm viscosity hampers aggregate polar segregation in Escherichia coli.

In Escherichia coli, under optimal conditions, protein aggregates associated with cellular aging are excluded from midcell by the nucleoid. We study the functionality of this process under sub-optimal temperatures from population and time lapse images of individual cells and aggregates and nucleoids within. We show that, as temperature decreases, aggregates become homogeneously distributed and uncorrelated with nucleoid size and location. We present evidence that this is due to increased cytoplasm viscosity, which weakens the anisotropy in aggregate displacements at the nucleoid borders that is responsible for their preference for polar localisation. Next, we show that in plasmolysed cells, which have increased cytoplasm viscosity, aggregates are also not preferentially located at the poles. Finally, we show that the inability of cells with increased viscosity to exclude aggregates from midcell results in enhanced aggregate concentration in between the nucleoids in cells close to dividing. This weakens the asymmetries in aggregate numbers between sister cells of subsequent generations required for rejuvenating cell lineages. We conclude that the process of exclusion of protein aggregates from midcell is not immune to stress conditions affecting the cytoplasm viscosity. The findings contribute to our understanding of E. coli's internal organisation and functioning, and its fragility to stressful conditions.

Characterizing rate limiting steps in transcription from RNA production times in live cells.

MOTIVATION: Single-molecule measurements of live Escherichia coli transcription dynamics suggest that this process ranges from sub- to super-Poissonian, depending on the conditions and on the promoter. For its accurate quantification, we propose a model that accommodates all these settings, and statistical methods to estimate the model parameters and to select the relevant components. RESULTS: The new methodology has improved accuracy and avoids overestimating the transcription rate due to finite measurement time, by exploiting unobserved data and by accounting for the effects of discrete sampling. First, we use Monte Carlo simulations of models based on measurements to show that the methods are reliable and offer substantial improvements over previous methods. Next, we apply the methods on measurements of transcription intervals of different promoters in live E. coli, and show that they produce significantly different results, both in low- and high-noise settings, and that, in the latter case, they even lead to qualitatively different results. Finally, we demonstrate that the methods can be generalized for other similar purposes, such as for estimating gene activation kinetics. In this case, the new methods allow quantifying the inducer uptake dynamics as opposed to just comparing them between cases, which was not previously possible. We expect this new methodology to be a valuable tool for functional analysis of cellular processes using single-molecule or single-event microscopy measurements in live cells. AVAILABILITY AND IMPLEMENTATION: Source code is available under Mozilla Public License at http://www.cs.tut.fi/%7Ehakkin22/censored/ CONTACT: andre.ribeiro@tut.fi or andre.sanchesribeiro@tut.fi SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Temperature-Dependent Model of Multi-step Transcription Initiation in Escherichia coli Based on Live Single-Cell Measurements.

Transcription kinetics is limited by its initiation steps, which differ between promoters and with intra- and extracellular conditions. Regulation of these steps allows tuning both the rate and stochasticity of RNA production. We used time-lapse, single-RNA microscopy measurements in live Escherichia coli to study how the rate-limiting steps in initiation of the Plac/ara-1 promoter change with temperature and induction scheme. For this, we compared detailed stochastic models fit to the empirical data in maximum likelihood sense using statistical methods. Using this analysis, we found that temperature affects the rate limiting steps unequally, as nonlinear changes in the closed complex formation suffice to explain the differences in transcription dynamics between conditions. Meanwhile, a similar analysis of the PtetA promoter revealed that it has a different rate limiting step configuration, with temperature regulating different steps. Finally, we used the derived models to explore a possible cause for why the identified steps are preferred as the main cause for behavior modifications with temperature: we find that transcription dynamics is either insensitive or responds reciprocally to changes in the other steps. Our results suggests that different promoters employ different rate limiting step patterns that control not only their rate and variability, but also their sensitivity to environmental changes.

Estimation of GFP-tagged RNA numbers from temporal fluorescence intensity data.

MOTIVATION: MS2-GFP-tagging of RNA is currently the only method to measure intervals between consecutive transcription events in live cells. For this, new transcripts must be accurately detected from intensity time traces. RESULTS: We present a novel method for automatically estimating RNA numbers and production intervals from temporal data of cell fluorescence intensities that reduces uncertainty by exploiting temporal information. We also derive a robust variant, more resistant to outliers caused e.g. by RNAs moving out of focus. Using Monte Carlo simulations, we show that the quantification of RNA numbers and production intervals is generally improved compared with previous methods. Finally, we analyze data from live Escherichia coli and show statistically significant differences to previous methods. The new methods can be used to quantify numbers and production intervals of any fluorescent probes, which are present in low copy numbers, are brighter than the cell background and degrade slowly. AVAILABILITY: Source code is available under Mozilla Public License at http://www.cs.tut.fi/%7ehakkin22/jumpdet/.

The lifelong socioeconomic disadvantage of single-mother background - the Helsinki Birth Cohort study 1934-1944.

BACKGROUND: Growing up with one parent is associated with economic hardship and health disadvantages, but there is limited evidence of its lifetime consequences. We examined whether being born to an unmarried mother is associated with socioeconomic position and marital history over the lifespan. METHODS: We analysed data from the Helsinki Birth Cohort Study including birth, child welfare clinic and school healthcare records from people born in Helsinki, Finland, between 1934 and 1944. Using a unique personal identification number, we linked these data to information on adult socioeconomic position from census data at 5-year intervals between 1970 and 2000, obtained from Statistics Finland. RESULTS: Compared to children of married mothers, children of unmarried mothers were more likely to have lower educational attainment and occupational status (odds ratio for basic vs. tertiary education 3.40; 95 % confidence interval 2.17 to 5.20; for lowest vs. highest occupational category 2.75; 1.92 to 3.95). They were also less likely to reach the highest income third in adulthood and more likely to stay unmarried themselves. The associations were also present when adjusted for childhood socioeconomic position. CONCLUSION: Being born to an unmarried mother, in a society where marriage is the norm, is associated with socioeconomic disadvantage throughout life, over and above the disadvantage associated with childhood family occupational status. This disadvantage may in part mediate the association between low childhood socioeconomic position and health in later life.

Estimation of fluorescence-tagged RNA numbers from spot intensities.

MOTIVATION: Present research on gene expression using live cell imaging and fluorescent proteins or tagged RNA requires accurate automated methods of quantification of these molecules from the images. Here, we propose a novel automated method for classifying pixel intensities of fluorescent spots to RNA numbers. RESULTS: The method relies on a new model of intensity distributions of tagged RNAs, for which we estimated parameter values in maximum likelihood sense from measurement data, and constructed a maximum a posteriori classifier to estimate RNA numbers in fluorescent RNA spots. We applied the method to estimate the number of tagged RNAs in individual live Escherichia coli cells containing a gene coding for an RNA with MS2-GFP binding sites. We tested the method using two constructs, coding for either 96 or 48 binding sites, and obtained similar distributions of RNA numbers, showing that the method is adaptive. We further show that the results agree with a method that uses time series data and with quantitative polymerase chain reaction measurements. Lastly, using simulated data, we show that the method is accurate in realistic parameter ranges. This method should, in general, be applicable to live single-cell measurements of low-copy number fluorescence-tagged molecules. AVAILABILITY AND IMPLEMENTATION: MATLAB extensions written in C for parameter estimation and finding decision boundaries are available under Mozilla public license at http://www.cs.tut.fi/%7ehakkin22/estrna/ CONTACT: andre.ribeiro@tut.fi.

Enabling safe dry cake disposal of bauxite residue by deliquoring and washing with a membrane filter press.

Dry cake disposal is the preferred technique for the disposal of bauxite residue, when considering environmental issues together with possible future utilisation of the solids. In order to perform dry cake disposal in an economical way, the deliquoring of the residue must be carried out efficiently, and it is also important to wash the obtained solids well to minimise the amount of soluble soda within the solids. The study presented in this article aims at detecting the most important variables influencing the deliquoring and washing of bauxite residue, performed with a horizontal membrane filter press and by determining the optimal washing conditions. The results obtained from pilot-scale experiments are evaluated by considering the properties of the solids, for instance, the residual alkali and aluminium content, as well as the consumption of wash liquid. Two different cake washing techniques, namely classic washing and channel washing, are also used and their performances compared. The results show that cake washing can be performed successfully in a horizontal membrane filter press, and significant improvements in the recovery of alkali and aluminium can be achieved compared with pressure filtration carried out without washing, or especially compared with the more traditionally used vacuum filtration.