Philipp Spahn

Differential gene expression during metamorphosis: a promising approach for age estimation of forensically important Calliphora vicina pupae (Diptera: Calliphoridae)

Necrophagous blow fly larvae can provide accurate estimates of the minimum postmortem interval in death investigations. During larval development, predictable morphological changes occur and measurements of weight, length, and width are compared to species-specific growth curves for reliable age estimates. However, aging blow fly pupae is more challenging because morphological and anatomical changes are not visible with the naked eye. Thus, delicate preparation of the pupae or rearing to the adult stage seems unavoidable. Conversely, metamorphosis evokes a remodelling of the larval shape to adult structures, and gene expression analysis potentially serves as a molecular tool to mirror the ageing process of a pupa. The present study focuses on the differential expression of two newly described, arbitrarily named genes (15_2, 2014192) and two previously identified genes (actin, arylphorin receptor) during Calliphora vicina (Diptera: Calliphoridae) metamorphosis. Quantification through real-time PCR revealed significant up- and downregulation of these transcripts found to be temperature dependent and age specific, hence, a new possibility to age forensically important blow fly pupae.

Tracking the evolutionary history of Cortinarius species in section Calochroi, with transoceanic disjunct distributions.

BackgroundCortinarius species in section Calochroi display local, clinal and circumboreal patterns of distribution across the Northern Hemisphere where these ectomycorrhizal fungi occur with host trees throughout their geographical range within a continent, or have disjunct intercontinental distributions, the origins of which are not understood. We inferred evolutionary histories of four species, 1) C. arcuatorum, 2) C. aureofulvus, 3) C. elegantior and 4) C. napus, from populations distributed throughout the Old World, and portions of the New World (Central- and North America) based on genetic variation of 154 haplotype internal transcribed spacer (ITS) sequences from 83 population samples. By describing the population structure of these species across their geographical distribution, we attempt to identify their historical migration and patterns of diversification.ResultsModels of population structure from nested clade, demographic and coalescent-based analyses revealed genetically differentiated and geographically structured haplotypes in C. arcuatorum and C. elegantior, while C. aureofulvus showed considerably less population structure and C. napus lacked sufficient genetic differentiation to resolve any population structure. Disjunct populations within C. arcuatorum, C. aureofulvus and C. elegantior show little or no morphological differentiation, whereas in C. napus there is a high level of homoplasy and phenotypic plasticity for veil and lamellae colour. The ITS sequences of the type specimens of C. albobrunnoides and C. albobrunnoides var. violaceovelatus were identical to one another and are treated as one species with a wider range of geographic distribution under C. napus.ConclusionsOur results indicate that each of the Calochroi species has undergone a relatively independent evolutionary history, hypothesised as follows: 1) a widely distributed ancestral population of C. arcuatorum diverged into distinctive sympatric populations in the New World; 2) two divergent lineages in C. elegantior gave rise to the New World and Old World haplotypes, respectively; and 3) the low levels of genetic divergence within C. aureofulvus and C. napus may be the result of more recent demographic population expansions. The scenario of migration via the Bering Land Bridge provides the most probable explanation for contemporaneous disjunct geographic distributions of these species, but it does not offer an explanation for the low degree of genetic divergence between populations of C. aureofulvus and C. napus. Our findings are mostly consistent with the designation of New World allopatric populations as separate species from the European counterpart species C. arcuatorum and C. elegantior. We propose the synonymy of C. albobrunnoides, C. albobrunnoides var. violaceovelatus and C. subpurpureophyllus var. sulphureovelatus with C. napus. The results also reinforce previous observations that linked C. arcuatorum and C. aureofulvus displaying distributions in parts of North America and Europe. Interpretations of the population structure of these fungi suggest that host tree history has heavily influenced their modern distributions; however, the complex issues related to co-migration of these fungi with their tree hosts remain unclear at this time.

A Consensus Genome-scale Reconstruction of Chinese Hamster Ovary Cell Metabolism

Chinese hamster ovary (CHO) cells dominate biotherapeutic protein production and are widely used in mammalian cell line engineering research. To elucidate metabolic bottlenecks in protein production and to guide cell engineering and bioprocess optimization, we reconstructed the metabolic pathways in CHO and associated them with >1,700 genes in the Cricetulus griseus genome. The genome-scale metabolic model based on this reconstruction, iCHO1766, and cell-line-specific models for CHO-K1, CHO-S, and CHO-DG44 cells provide the biochemical basis of growth and recombinant protein production. The models accurately predict growth phenotypes and known auxotrophies in CHO cells. With the models, we quantify the protein synthesis capacity of CHO cells and demonstrate that common bioprocess treatments, such as histone deacetylase inhibitors, inefficiently increase product yield. However, our simulations show that the metabolic resources in CHO are more than three times more efficiently utilized for growth or recombinant protein synthesis following targeted efforts to engineer the CHO secretory pathway. This model will further accelerate CHO cell engineering and help optimize bioprocesses.

Systems glycobiology for glycoengineering

Glycosylation serves essential functions on many proteins produced in biopharmaceutical manufacturing, making it mandatory to thoroughly consider its biogenesis during the production process. Glycoengineering efforts involve the rational design of glycosylation through adjustments in culturing conditions or genetic modifications. Computational models have been developed to aid this process, aiming to offer cheaper and faster alternatives to costly screening strategies. Recently, these models have been successfully utilized to predict glycosylation of products of industrial relevance. Furthermore, systems-level analyses of glycan diversity are elucidating deeper insights into the mechanisms underlying glycosylation. As computational models of glycosylation continue to be expanded, refined, and leveraged for detailed analysis of glycomics data, they will become invaluable resources for cell line development and glycoengineering.

A Markov chain model for N-linked protein glycosylation – towards a low-parameter tool for model-driven glycoengineering

Glycosylation is a critical quality attribute of most recombinant biotherapeutics. Consequently, drug development requires careful control of glycoforms to meet bioactivity and biosafety requirements. However, glycoengineering can be extraordinarily difficult given the complex reaction networks underlying glycosylation and the vast number of different glycans that can be synthesized in a host cell. Computational modeling offers an intriguing option to rationally guide glycoengineering, but the high parametric demands of current modeling approaches pose challenges to their application. Here we present a novel low-parameter approach to describe glycosylation using flux-balance and Markov chain modeling. The model recapitulates the biological complexity of glycosylation, but does not require user-provided kinetic information. We use this method to predict and experimentally validate glycoprofiles on EPO, IgG as well as the endogenous secretome following glycosyltransferase knock-out in different Chinese hamster ovary (CHO) cell lines. Our approach offers a flexible and user-friendly platform that can serve as a basis for powerful computational engineering efforts in mammalian cell factories for biopharmaceutical production.

A Vertex Model of Drosophila Ventral Furrow Formation

Ventral furrow formation in Drosophila is an outstanding model system to study the mechanisms involved in large-scale tissue rearrangements. Ventral cells accumulate myosin at their apical sides and, while being tightly coupled to each other via apical adherens junctions, execute actomyosin contractions that lead to reduction of their apical cell surface. Thereby, a band of constricted cells along the ventral epithelium emerges which will form a tissue indentation along the ventral midline (the ventral furrow). Here we adopt a 2D vertex model to simulate ventral furrow formation in a surface view allowing easy comparison with confocal live-recordings. We show that in order to reproduce furrow morphology seen in vivo, a gradient of contractility must be assumed in the ventral epithelium which renders cells more contractile the closer they lie to the ventral midline. The model predicts previous experimental findings, such as the gain of eccentric morphology of constricting cells and an incremental fashion of apical cell area reduction. Analysis of the model suggests that this incremental area reduction is caused by the dynamical interplay of cell elasticity and stochastic contractility as well as by the opposing forces from contracting neighbour cells. We underpin results from the model through in vivo analysis of ventral furrow formation in wildtype and twi mutant embryos. Our results show that ventral furrow formation can be accomplished as a “tug-of-war” between stochastically contracting, mechanically coupled cells and may require less rigorous regulation than previously thought. Summary For the developmental biologist it is a fascinating question how cells can coordinate major tissue movements during embryonic development. The so-called ventral furrow of the Drosophila embryo is a well-studied example of such a process when cells from a ventral band, spanning nearly the entire length of the embryo, undergo dramatic shape change by contracting their tips and then fold inwards into the interior of the embryo. Although numerous genes have been identified that are critical for ventral furrow formation, it is an open question how cells work together to elicit this tissue rearrangement. We use a computational model to mimic the physical properties of cells in the ventral epithelium and simulate the formation of the furrow. We find that the ventral furrow can form through stochastic self-organisation and that previous experimental observations can be readily explained in our model by considering forces that arise when cells execute contractions while being coupled to each other in a mechanically coherent epithelium. The model highlights the importance of a physical perspective when studying tissue morphogenesis and shows that only a minimal genetic regulation may be required to drive complex processes in embryonic development.

The PDZ-GEF protein Dizzy regulates the establishment of adherens junctions required for ventral furrow formation in Drosophila

Summary The PDZ-GEF protein Dizzy (Dzy) and its downstream GTPase Rap1 have pleiotropic roles during development of the Drosophila embryo. Here, we show that maternally provided Dzy and Rap1 first function during ventral furrow formation (VFF) where they are critical to guarantee rapid apical cell constrictions. Contraction of the apical actomyosin filament system occurs independently of Dzy and Rap1, but loss of Dzy results in a delayed establishment of the apical adherens junction (AJ) belt, whereas in the absence of Rap1 only a fragmentary apical AJ belt is formed in the epithelium. The timely establishment of apical AJs appears to be essential for coupling actomyosin contractions to cell shape change and to assure completion of the ventral furrow. Immediately after VFF, the downregulation of Dzy and Rap1 is necessary to allow normal mesodermal development to continue after the epithelial-to-mesenchymal transition, as overexpression of Dzy or of constitutively active Rap1 compromises mesodermal migration and monolayer formation. We propose that Dzy and Rap1 are crucial factors regulating the dynamics of AJs during gastrulation.

The analysis of temporal gene expression to estimate the age of forensically important blow fly pupae: results from three blind studies

The estimation of the minimum postmortem interval can be implemented by age estimation of corpse-associated primary colonizers such as the blow flies (Diptera: Calliphoridae). In cases where pupae represent the oldest stages found on a body, their age must be added to the duration of prepupal development to estimate the minimum postmortem interval. Although methods to age blow fly larvae have been well established using morphology, length or weight and age estimation of pupae has proved challenging. In a previous work, we quantified the changes in mRNA levels of four differentially expressed genes during the metamorphosis of Calliphora vicina pupae, hence representing molecular markers for pupal age (i.e., time elapsed since pupariation). Here, we demonstrate how these data can be used to estimate pupal age with inverse prediction. We present three blind studies conducted under various conditions and show that age of C. vicina pupae can be well estimated based on gene expression data. As these data are quantitative and can be processed automatically, gene expression has the potential to outperform morphological analysis in age estimation of forensically relevant blow fly pupae.

Predictive glycoengineering of biosimilars using a Markov chain glycosylation model.

Biosimilar drugs must closely resemble the pharmacological attributes of innovator products to ensure safety and efficacy to obtain regulatory approval. Glycosylation is one critical quality attribute that must be matched, but it is inherently difficult to control due to the complexity of its biogenesis. This usually implies that costly and time-consuming experimentation is required for clone identification and optimization of biosimilar glycosylation. Here, a computational method that utilizes a Markov model of glycosylation to predict optimal glycoengineering strategies to obtain a specific glycosylation profile with desired properties is described. The approach uses a genetic algorithm to find the required quantities to perturb glycosylation reaction rates that lead to the best possible match with a given glycosylation profile. Furthermore, the approach can be used to identify cell lines and clones that will require minimal intervention while achieving a glycoprofile that is most similar to the desired profile. Thus, this approach can facilitate biosimilar design by providing computational glycoengineering guidelines that can be generated with a minimal time and cost.

Modulating carbohydrate-protein interactions through glycoengineering of monoclonal antibodies to impact cancer physiology.

Diverse glycans on proteins impact cell and organism physiology, along with drug activity. Since many protein-based biotherapeutics are glycosylated and these glycans have biological activity, there is a desire to engineer glycosylation for recombinant protein-based biotherapeutics. Engineered glycosylation can impact the recombinant protein efficacy and also influence many cell pathways by first changing glycan-protein interactions and consequently modulating disease physiologies. However, its complexity is enormous. Recent advances in glycoengineering now make it easier to modulate protein-glycan interactions. Here, we discuss how engineered glycans contribute to therapeutic monoclonal antibodies (mAbs) in the treatment of cancers, how these glycoengineered therapeutic mAbs affect the transformed phenotypes and downstream cell pathways. Furthermore, we suggest how systems biology can help in the next generation mAb glycoengineering process by aiding in data analysis and guiding engineering efforts to tailor mAb glycan and ultimately drug efficacy, safety and affordability.