Interactions between ploidy and resource availability shape clonal interference at initiation and recurrence of glioblastoma.

Glioblastoma (GBM) is the most aggressive form of primary brain tumor. Complete surgical resection of GBM is almost impossible due to the infiltrative nature of the cancer. While no evidence for recent selection events have been found after diagnosis, the selective forces that govern gliomagenesis are strong, shaping the tumor's cell composition during the initial progression to malignancy with late consequences for invasiveness and therapy response. We present a mathematical model that simulates the growth and invasion of a glioma, given its ploidy level and the nature of its brain tissue micro-environment (TME), and use it to make inferences about GBM initiation and response to standard-of-care treatment. We approximate the spatial distribution of resource access in the TME through integration of in-silico modelling, multi-omics data and image analysis of primary and recurrent GBM. In the pre-malignant setting, our in-silico results suggest that low ploidy cancer cells are more resistant to starvation-induced cell death. In the malignant setting, between first and second surgery, simulated tumors with different ploidy compositions progressed at different rates. Whether higher ploidy predicted fast recurrence, however, depended on the TME. Historical data supports this dependence on TME resources, as shown by a significant correlation between the median glucose uptake rates in human tissues and the median ploidy of cancer types that arise in the respective tissues (Spearman r = -0.70; P = 0.026). Taken together our findings suggest that availability of metabolic substrates in the TME drives different cell fate decisions for cancer cells with different ploidy and shapes GBM disease initiation and relapse characteristics.

Mathematical Modeling of Clonal Interference by Density-Dependent Selection in Heterogeneous Cancer Cell Lines.

Many cancer cell lines are aneuploid and heterogeneous, with multiple karyotypes co-existing within the same cell line. Karyotype heterogeneity has been shown to manifest phenotypically, thus affecting how cells respond to drugs or to minor differences in culture media. Knowing how to interpret karyotype heterogeneity phenotypically would give insights into cellular phenotypes before they unfold temporally. Here, we re-analyzed single cell RNA (scRNA) and scDNA sequencing data from eight stomach cancer cell lines by placing gene expression programs into a phenotypic context. Using live cell imaging, we quantified differences in the growth rate and contact inhibition between the eight cell lines and used these differences to prioritize the transcriptomic biomarkers of the growth rate and carrying capacity. Using these biomarkers, we found significant differences in the predicted growth rate or carrying capacity between multiple karyotypes detected within the same cell line. We used these predictions to simulate how the clonal composition of a cell line would change depending on density conditions during in-vitro experiments. Once validated, these models can aid in the design of experiments that steer evolution with density-dependent selection.

Intra-tumor heterogeneity, turnover rate and karyotype space shape susceptibility to missegregation-induced extinction.

The phenotypic efficacy of somatic copy number alterations (SCNAs) stems from their incidence per base pair of the genome, which is orders of magnitudes greater than that of point mutations. One mitotic event stands out in its potential to significantly change a cell's SCNA burden-a chromosome missegregation. A stochastic model of chromosome mis-segregations has been previously developed to describe the evolution of SCNAs of a single chromosome type. Building upon this work, we derive a general deterministic framework for modeling missegregations of multiple chromosome types. The framework offers flexibility to model intra-tumor heterogeneity in the SCNAs of all chromosomes, as well as in missegregation- and turnover rates. The model can be used to test how selection acts upon coexisting karyotypes over hundreds of generations. We use the model to calculate missegregation-induced population extinction (MIE) curves, that separate viable from non-viable populations as a function of their turnover- and missegregation rates. Turnover- and missegregation rates estimated from scRNA-seq data are then compared to theoretical predictions. We find convergence of theoretical and empirical results in both the location of MIE curves and the necessary conditions for MIE. When a dependency of missegregation rate on karyotype is introduced, karyotypes associated with low missegregation rates act as a stabilizing refuge, rendering MIE impossible unless turnover rates are exceedingly high. Intra-tumor heterogeneity, including heterogeneity in missegregation rates, increases as tumors progress, rendering MIE unlikely.

Backbone and side chain NMR resonance assignments for an archaeal homolog of the endonuclease Nob1 involved in ribosome biogenesis.

Eukaryotic ribosome biogenesis requires the concerted action of ~200 auxiliary protein factors on the nascent ribosome. For many of these factors structural and functional information is still lacking. The endonuclease Nob1 has been recently identified in yeast as the enzyme responsible for the final cytoplasmatic trimming step of the pre-18S rRNA during the biogenesis of the small ribosomal subunit. Here we report the NMR resonance assignments for a Nob1 homolog from the thermophilic archeon Pyrococcus horikoshii as a prerequisite for further structural studies of this class of proteins.

Structural and functional analysis of the archaeal endonuclease Nob1.

Eukaryotic ribosome biogenesis requires the concerted action of numerous ribosome assembly factors, for most of which structural and functional information is currently lacking. Nob1, which can be identified in eukaryotes and archaea, is required for the final maturation of the small subunit ribosomal RNA in yeast by catalyzing cleavage at site D after export of the preribosomal subunit into the cytoplasm. Here, we show that this also holds true for Nob1 from the archaeon Pyrococcus horikoshii, which efficiently cleaves RNA-substrates containing the D-site of the preribosomal RNA in a manganese-dependent manner. The structure of PhNob1 solved by nuclear magnetic resonance spectroscopy revealed a PIN domain common with many nucleases and a zinc ribbon domain, which are structurally connected by a flexible linker. We show that amino acid residues required for substrate binding reside in the PIN domain whereas the zinc ribbon domain alone is sufficient to bind helix 40 of the small subunit rRNA. This suggests that the zinc ribbon domain acts as an anchor point for the protein on the nascent subunit positioning it in the proximity of the cleavage site.

Identification of a specific fucoxanthin-chlorophyll protein in the light harvesting complex of photosystem I in the diatom Cyclotella meneghiniana.

Thylakoids of the diatom Cyclotella meneghiniana were separated by discontinuous gradient centrifugation into photosystem (PS) I, PSII, and fucoxanthin-chlorophyll protein (FCP) fractions. FCPs are homologue to light harvesting complexes of higher plants with similar function in e.g. brown algae and diatoms. Still, it is unclear if FCP complexes are specifically associated with either PSI or PSII, or if FCP complexes function as one antenna for both photosystems. However, a trimeric FCP complex, FCPa, and a higher FCP oligomer, FCPb, have been described for C. meneghiniana, already. In this study, biochemical and spectroscopical evidences are provided that reveal a different subset of associated Fcp polypeptides within the isolated photosystem complexes. Whereas the PSII associated Fcp antenna resembles FCPa since it contains Fcp2 and Fcp6, at least three different Fcp polypeptides are associated with PSI. By re-solubilisation and a further purification step Fcp polypeptides were partially removed from PSI and both fractions were analysed again by biochemical and spectroscopical means, as well as by HPLC. Thereby a protein related to Fcp4 and a so far undescribed 17 kDa Fcp were found to be strongly coupled to PSI, whereas presumably Fcp5, a subunit of the FCPb complex, is only loosely bound to the PSI core. Thus, an association of FCPb and PSI is assumed.

The monomeric photosystem I-complex of the diatom Phaeodactylum tricornutum binds specific fucoxanthin chlorophyll proteins (FCPs) as light-harvesting complexes.

A photosystem I (PSI)-fucoxanthin chlorophyll protein (FCP) complex with a chlorophyll a/P700 ratio of approximately 200:1 was isolated from the diatom Phaeodactylum tricornutum. Spectroscopic analysis proved that the more tightly bound FCP functions as a light-harvesting complex, actively transferring light energy from its accessory pigments chlorophyll c and fucoxanthin to the PSI core. Using an antibody against all FCP polypeptides of Cyclotella cryptica it could be shown that the polypeptides of the major FCP fraction differ from the FCPs found in the PSI fraction. Since these FCPs are tightly bound to PSI, active in energy transfer, and not found in the main FCP fraction, we suppose them to be PSI specific. Blue Native-PAGE, gel filtration and first electron microscopy studies of the PSI-FCP sample revealed a monomeric complex comparable in size and shape to the PSI-LHCI complex of green algae.