Center for Computational and Theoretical Biology
Refine
Has Fulltext
- yes (45)
Is part of the Bibliography
- yes (45)
Year of publication
Document Type
- Journal article (39)
- Doctoral Thesis (5)
- Preprint (1)
Language
- English (45) (remove)
Keywords
- machine learning (5)
- active zone (4)
- Caenorhabditis elegans (2)
- bee decline (2)
- dSTORM (2)
- deep learning (2)
- evolution (2)
- foraging (2)
- genetic architecture (2)
- juvenile hormone (2)
Institute
- Center for Computational and Theoretical Biology (45)
- Theodor-Boveri-Institut für Biowissenschaften (20)
- Pathologisches Institut (4)
- Physiologisches Institut (4)
- Julius-von-Sachs-Institut für Biowissenschaften (3)
- Klinik und Poliklinik für Unfall-, Hand-, Plastische und Wiederherstellungschirurgie (Chirurgische Klinik II) (3)
- Neurochirurgische Klinik und Poliklinik (3)
- Urologische Klinik und Poliklinik (3)
- Comprehensive Cancer Center Mainfranken (2)
- Klinik und Poliklinik für Dermatologie, Venerologie und Allergologie (2)
- Klinik und Poliklinik für Psychiatrie, Psychosomatik und Psychotherapie (2)
- Deutsches Zentrum für Herzinsuffizienz (DZHI) (1)
- Graduate School of Life Sciences (1)
- Institut für Experimentelle Biomedizin (1)
- Institut für Molekulare Infektionsbiologie (1)
- Institut für Theoretische Physik und Astrophysik (1)
- Institut für diagnostische und interventionelle Radiologie (Institut für Röntgendiagnostik) (1)
- Klinik und Poliklinik für Nuklearmedizin (1)
- Medizinische Klinik und Poliklinik I (1)
- Medizinische Klinik und Poliklinik II (1)
EU-Project number / Contract (GA) number
- 2020010013 (1)
- 250194-Carnivorom (1)
- 835102) (1)
Tropical peat swamp forests sequester globally significant stores of carbon in deep layers of waterlogged, anoxic, acidic and nutrient-depleted peat. The roles of microbes in supporting these forests through the formation of peat, carbon sequestration and nutrient cycling are virtually unknown. This study investigated physicochemical peat properties and microbial diversity between three dominant tree species: Shorea uliginosa (Dipterocarpaceae), Koompassia malaccensis (legumes associated with nitrogen-fixing bacteria), Eleiodoxa conferta (palm) and depths (surface, 45 and 90 cm) using microbial 16S rRNA gene amplicon sequencing. Water pH, oxygen, nitrogen, phosphorus, total phenolic contents and C/N ratio differed significantly between depths, but not tree species. Depth also strongly influenced microbial diversity and composition, while both depth and tree species exhibited significant impact on the archaeal communities. Microbial diversity was highest at the surface, where fresh leaf litter accumulates, and nutrient supply is guaranteed. Nitrogen was the core parameter correlating to microbial communities, but the interactive effects from various environmental variables displayed significant correlation to relative abundance of major microbial groups. Proteobacteria was the dominant phylum and the most abundant genus, Rhodoplanes, might be involved in nitrogen fixation. The most abundant methanogens and methanotrophs affiliated, respectively, to families Methanomassiliicoccaceae and Methylocystaceae. Our results demonstrated diverse microbial communities and provide valuable insights on microbial ecology in these extreme ecosystems.
Young grapevines (Vitis vinifera) suffer and eventually can die from the crown gall disease caused by the plant pathogen Allorhizobium vitis (Rhizobiaceae). Virulent members of A. vitis harbor a tumor-inducing plasmid and induce formation of crown galls due to the oncogenes encoded on the transfer DNA. The expression of oncogenes in transformed host cells induces unregulated cell proliferation and metabolic and physiological changes. The crown gall produces opines uncommon to plants, which provide an important nutrient source for A. vitis harboring opine catabolism enzymes. Crown galls host a distinct bacterial community, and the mechanisms establishing a crown gall–specific bacterial community are currently unknown. Thus, we were interested in whether genes homologous to those of the tumor-inducing plasmid coexist in the genomes of the microbial species coexisting in crown galls. We isolated 8 bacterial strains from grapevine crown galls, sequenced their genomes, and tested their virulence and opine utilization ability in bioassays. In addition, the 8 genome sequences were compared with 34 published bacterial genomes, including closely related plant-associated bacteria not from crown galls. Homologous genes for virulence and opine anabolism were only present in the virulent Rhizobiaceae. In contrast, homologs of the opine catabolism genes were present in all strains including the nonvirulent members of the Rhizobiaceae and non-Rhizobiaceae. Gene neighborhood and sequence identity of the opine degradation cluster of virulent and nonvirulent strains together with the results of the opine utilization assay support the important role of opine utilization for cocolonization in crown galls, thereby shaping the crown gall community.
Abstract
Cell lineage decisions occur in three-dimensional spatial patterns that are difficult to identify by eye. There is an ongoing effort to replicate such patterns using mathematical modeling. One approach uses long ranging cell-cell communication to replicate common spatial arrangements like checkerboard and engulfing patterns. In this model, the cell-cell communication has been implemented as a signal that disperses throughout the tissue. On the other hand, machine learning models have been developed for pattern recognition and pattern reconstruction tasks. We combined synthetic data generated by the mathematical model with spatial summary statistics and deep learning algorithms to recognize and reconstruct cell fate patterns in organoids of mouse embryonic stem cells. Application of Moran’s index and pair correlation functions for in vitro and synthetic data from the model showed local clustering and radial segregation. To assess the patterns as a whole, a graph neural network was developed and trained on synthetic data from the model. Application to in vitro data predicted a low signal dispersion value. To test this result, we implemented a multilayer perceptron for the prediction of a given cell fate based on the fates of the neighboring cells. The results show a 70% accuracy of cell fate imputation based on the nine nearest neighbors of a cell. Overall, our approach combines deep learning with mathematical modeling to link cell fate patterns with potential underlying mechanisms.
Author summary
Mammalian embryo development relies on organized differentiation of stem cells into different lineages. Particularly at the early stages of embryogenesis, cells of different fates form three-dimensional spatial patterns that are difficult to identify by eye. Pattern quantification and mathematical modeling have produced first insights into potential mechanisms for the cell fate arrangements. However, these approaches have relied on classifications of the patterns such as inside-out or random, or used summary statistics such as pair correlation functions or cluster radii. Deep neural networks allow characterizing patterns directly. Since the tissue context can be readily reproduced by a graph, we implemented a graph neural network to characterize the patterns of embryonic stem cell organoids as a whole. In addition, we implemented a multilayer perceptron model to reconstruct the fate of a given cell based on its neighbors. To train and test the models, we used synthetic data generated by our mathematical model for cell-cell communication. This interplay of deep learning and mathematical modeling in combination with summary statistics allowed us to identify a potential mechanism for cell fate determination in mouse embryonic stem cells. Our results agree with a mechanism with a dispersion of the intercellular signal that links a cell’s fate to those of the local neighborhood.
Summary
Embryos develop in a concerted sequence of spatiotemporal arrangements of cells. In the preimplantation mouse embryo, the distribution of the cells in the inner cell mass evolves from a salt-and-pepper pattern to spatial segregation of two distinct cell types. The exact properties of the salt-and-pepper pattern have not been analyzed so far. We investigate the spatiotemporal distribution of NANOG- and GATA6-expressing cells in the ICM of the mouse blastocysts with quantitative three-dimensional single-cell-based neighborhood analyses. A combination of spatial statistics and agent-based modeling reveals that the cell fate distribution follows a local clustering pattern. Using ordinary differential equations modeling, we show that this pattern can be established by a distance-based signaling mechanism enabling cells to integrate information from the whole inner cell mass into their cell fate decision. Our work highlights the importance of longer-range signaling to ensure coordinated decisions in groups of cells to successfully build embryos.
Highlights
• The local cell neighborhood and global ICM population composition correlate
• ICM cells show characteristics of local clustering in early and mid mouse blastocysts
• ICM patterning requires integration of signals from cells beyond the first neighbors
Integrative, three-dimensional \(in\) \(silico\) modeling of gas exchange in the human alveolus
(2024)
The lung plays a vital role by exchanging respiratory gases. At the core of this gas exchange is a simple yet crucial passive diffusion process occurring within the alveoli. These balloon-like structures, connected to the peripheral airways, are surrounded by a dense network
of small capillaries. Here, inhaled air comes into close proximity with deoxygenated blood coming from the heart, enabling the exchange of oxygen and carbon dioxide across their concentration gradients.
The efficiency of gas exchange can be measured through indicators such as the diffusion capacity of the lung for oxygen and the reaction half-time. A notable discrepancy exists in humans between physiological estimates of diffusion capacity and the theoretical maximum capacity under optimal structural conditions (morphological estimate). This discrepancy is influenced by a range of interrelated factors, including structural elements like the surface area and thickness of the diffusion barrier, as well as physiological factors such as blood flow dynamics. To unravel the different roles of these factors, we investigated how morphological and physiological properties of the human alveolar micro-environment collectively and individually influence the process of gas exchange. To this end, we developed an integrative in silico approach combining 3D morphological modeling and simulation of blood flow and of oxygen transport.
At the core of our approach lies the simulation software Alvin, serving as an interactive platform for the underlying mathematical model of oxygen transport within the alveolus. Developed by integrating and expanding existing mathematical models, our spatio-temporal model produces results in agreement with experimental data. Alvin allows for real-time parameter adjustments and the execution of multiple simultaneous simulation instances and provides detailed quantitative feedback, offering an immersive exploration of the simulated gas exchange process. The morphological and physiological parameters at play were further investigated with a focus on the microvasculature. By compiling a stereological database from the literature and 3D geometric modeling, we created a sheet-flow model as a realistic representation of the morphology of the human alveolar capillary network. Blood flow was simulated using computational fluid dynamics. Our findings were in line with previous estimations and highlighted the crucial role of viscosity models in predicting pressure drop across the microvasculature. Furthermore, we showcased how our approach can be harnessed to explore structural details, such as the connectivity of the alveolar capillary network with the vascular tree, using blood flow indices. It is important to emphasize that
so far we have relied on different data sources and that experimental validation is needed to move forward.
Integration of our findings into Alvin allowed quantification of the simulated gas exchange process through the diffusion capacity for oxygen and reaction half-time. In addition to evaluating the collective influences of the morphological and physiological properties, our interactive software facilitates the assessment of individual parameter value changes. Exploring blood volume and surface area available for gas exchange revealed linear correlations with diffusion capacity. The blood flow velocity had a positive, non-linear effect on diffusion capacity. The reaction half-time confirmed that under normal conditions, the gas exchange process is not diffusion-limited. Collectively, our alveolar model yielded a diffusion capacity value that fell in the middle of previous physiological and morphological estimates, implying that alveolar-level phenomena contribute to 50% of the diffusion capacity limitations that occur in vivo.
In summary, our integrative in silico approach disentangles various structural and functional influences on alveolar gas exchange, complementing traditional investigations in respiratory
research. We further showcase its utility in teaching and the interpretation of published data. To advance our understanding, future work should prioritize obtaining a cohesive experimental data set and identifying an appropriate viscosity model for blood flow simulations.
(1) Background: C-X-C Motif Chemokine Receptor 4 (CXCR4) and Fibroblast Activation Protein Alpha (FAP) are promising theranostic targets. However, it is unclear whether CXCR4 and FAP positivity mark distinct microenvironments, especially in solid tumors. (2) Methods: Using Random Forest (RF) analysis, we searched for entity-independent mRNA and microRNA signatures related to CXCR4 and FAP overexpression in our pan-cancer cohort from The Cancer Genome Atlas (TCGA) database — representing n = 9242 specimens from 29 tumor entities. CXCR4- and FAP-positive samples were assessed via StringDB cluster analysis, EnrichR, Metascape, and Gene Set Enrichment Analysis (GSEA). Findings were validated via correlation analyses in n = 1541 tumor samples. TIMER2.0 analyzed the association of CXCR4 / FAP expression and infiltration levels of immune-related cells. (3) Results: We identified entity-independent CXCR4 and FAP gene signatures representative for the majority of solid cancers. While CXCR4 positivity marked an immune-related microenvironment, FAP overexpression highlighted an angiogenesis-associated niche. TIMER2.0 analysis confirmed characteristic infiltration levels of CD8+ cells for CXCR4-positive tumors and endothelial cells for FAP-positive tumors. (4) Conclusions: CXCR4- and FAP-directed PET imaging could provide a non-invasive decision aid for entity-agnostic treatment of microenvironment in solid malignancies. Moreover, this machine learning workflow can easily be transferred towards other theranostic targets.
Purpose
To fully automatically derive quantitative parameters from late gadolinium enhancement (LGE) cardiac MR (CMR) in patients with myocardial infarction and to investigate if phase sensitive or magnitude reconstructions or a combination of both results in best segmentation accuracy.
Methods
In this retrospective single center study, a convolutional neural network with a U-Net architecture with a self-configuring framework (“nnU-net”) was trained for segmentation of left ventricular myocardium and infarct zone in LGE-CMR. A database of 170 examinations from 78 patients with history of myocardial infarction was assembled. Separate fitting of the model was performed, using phase sensitive inversion recovery, the magnitude reconstruction or both contrasts as input channels.
Manual labelling served as ground truth. In a subset of 10 patients, the performance of the trained models was evaluated and quantitatively compared by determination of the Sørensen-Dice similarity coefficient (DSC) and volumes of the infarct zone compared with the manual ground truth using Pearson’s r correlation and Bland-Altman analysis.
Results
The model achieved high similarity coefficients for myocardium and scar tissue. No significant difference was observed between using PSIR, magnitude reconstruction or both contrasts as input (PSIR and MAG; mean DSC: 0.83 ± 0.03 for myocardium and 0.72 ± 0.08 for scars). A strong correlation for volumes of infarct zone was observed between manual and model-based approach (r = 0.96), with a significant underestimation of the volumes obtained from the neural network.
Conclusion
The self-configuring nnU-net achieves predictions with strong agreement compared to manual segmentation, proving the potential as a promising tool to provide fully automatic quantitative evaluation of LGE-CMR.
Osmotic stress can be detrimental to plants, whose survival relies heavily on proteomic plasticity. Protein ubiquitination is a central post-translational modification in osmotic-mediated stress. In this study, we used the K-Ɛ-GG antibody enrichment method integrated with high-resolution mass spectrometry to compile a list of 719 ubiquitinated lysine (K-Ub) residues from 450 Arabidopsis root membrane proteins (58% of which are transmembrane proteins), thereby adding to the database of ubiquitinated substrates in plants. Although no ubiquitin (Ub) motifs could be identified, the presence of acidic residues close to K-Ub was revealed. Our ubiquitinome analysis pointed to a broad role of ubiquitination in the internalization and sorting of cargo proteins. Moreover, the simultaneous proteome and ubiquitinome quantification showed that ubiquitination is mostly not involved in membrane protein degradation in response to short osmotic treatment but that it is putatively involved in protein internalization, as described for the aquaporin PIP2;1. Our in silico analysis of ubiquitinated proteins shows that two E2 Ub-conjugating enzymes, UBC32 and UBC34, putatively target membrane proteins under osmotic stress. Finally, we revealed a positive role for UBC32 and UBC34 in primary root growth under osmotic stress.
Automatic image reconstruction is critical to cope with steadily increasing data from advanced microscopy. We describe here the Fiji macro 3D ART VeSElecT which we developed to study synaptic vesicles in electron tomograms. We apply this tool to quantify vesicle properties (i) in embryonic Danio rerio 4 and 8 days past fertilization (dpf) and (ii) to compare Caenorhabditis elegans N2 neuromuscular junctions (NMJ) wild-type and its septin mutant (unc-59(e261)). We demonstrate development-specific and mutant-specific changes in synaptic vesicle pools in both models. We confirm the functionality of our macro by applying our 3D ART VeSElecT on zebrafish NMJ showing smaller vesicles in 8 dpf embryos then 4 dpf, which was validated by manual reconstruction of the vesicle pool. Furthermore, we analyze the impact of C. elegans septin mutant unc-59(e261) on vesicle pool formation and vesicle size. Automated vesicle registration and characterization was implemented in Fiji as two macros (registration and measurement). This flexible arrangement allows in particular reducing false positives by an optional manual revision step. Preprocessing and contrast enhancement work on image-stacks of 1nm/pixel in x and y direction. Semi-automated cell selection was integrated. 3D ART VeSElecT removes interfering components, detects vesicles by 3D segmentation and calculates vesicle volume and diameter (spherical approximation, inner/outer diameter). Results are collected in color using the RoiManager plugin including the possibility of manual removal of non-matching confounder vesicles. Detailed evaluation considered performance (detected vesicles) and specificity (true vesicles) as well as precision and recall. We furthermore show gain in segmentation and morphological filtering compared to learning based methods and a large time gain compared to manual segmentation. 3D ART VeSElecT shows small error rates and its speed gain can be up to 68 times faster in comparison to manual annotation. Both automatic and semi-automatic modes are explained including a tutorial.
Ultrastructural analysis of wild-type and RIM1α knockout active zones in a large cortical synapse
(2022)
Rab3A-interacting molecule (RIM) is crucial for fast Ca\(^{2+}\)-triggered synaptic vesicle (SV) release in presynaptic active zones (AZs). We investigated hippocampal giant mossy fiber bouton (MFB) AZ architecture in 3D using electron tomography of rapid cryo-immobilized acute brain slices in RIM1α\(^{−/−}\) and wild-type mice. In RIM1α\(^{−/−}\), AZs are larger with increased synaptic cleft widths and a 3-fold reduced number of tightly docked SVs (0–2 nm). The distance of tightly docked SVs to the AZ center is increased from 110 to 195 nm, and the width of their electron-dense material between outer SV membrane and AZ membrane is reduced. Furthermore, the SV pool in RIM1α\(^{−/−}\) is more heterogeneous. Thus, RIM1α, besides its role in tight SV docking, is crucial for synaptic architecture and vesicle pool organization in MFBs.