Center for Computational and Theoretical Biology
Refine
Has Fulltext
- yes (47)
Is part of the Bibliography
- yes (47)
Year of publication
Document Type
- Journal article (40)
- Doctoral Thesis (6)
- Preprint (1)
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 (47)
- 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)
EU-Project number / Contract (GA) number
- 2020010013 (1)
- 250194-Carnivorom (1)
- 835102) (1)
Spatial biological networks are abundant on all scales of life, from single cells to ecosystems, and perform various important functions including signal transmission and nutrient transport. These biological functions depend on the architecture of the network, which emerges as the result of a dynamic, feedback-driven developmental process. While cell behavior during growth can be genetically encoded, the resulting network structure depends on spatial constraints and tissue architecture. Since network growth is often difficult to observe experimentally, computer simulations can help to understand how local cell behavior determines the resulting network architecture. We present here a computational framework based on directional statistics to model network formation in space and time under arbitrary spatial constraints. Growth is described as a biased correlated random walk where direction and branching depend on the local environmental conditions and constraints, which are presented as 3D multilayer grid. To demonstrate the application of our tool, we perform growth simulations of a dense network between cells and compare the results to experimental data from osteocyte networks in bone. Our generic framework might help to better understand how network patterns depend on spatial constraints, or to identify the biological cause of deviations from healthy network function.
Author summary
We present a novel modeling approach and computational implementation to better understand the development of spatial biological networks under the influence of external signals. Our tool allows us to study the relationship between local biological growth parameters and the emerging macroscopic network function using simulations. This computational approach can generate plausible network graphs that take local feedback into account and provide a basis for comparative studies using graph-based methods.
White Blood Cell (WBC) Leukaemia is caused by excessive production of leukocytes in the bone marrow, and image-based detection of malignant WBCs is important for its detection. Convolutional Neural Networks (CNNs) present the current state-of-the-art for this type of image classification, but their computational cost for training and deployment can be high. We here present an improved hybrid approach for efficient classification of WBC Leukemia. We first extract features from WBC images using VGGNet, a powerful CNN architecture, pre-trained on ImageNet. The extracted features are then filtered using a statistically enhanced Salp Swarm Algorithm (SESSA). This bio-inspired optimization algorithm selects the most relevant features and removes highly correlated and noisy features. We applied the proposed approach to two public WBC Leukemia reference datasets and achieve both high accuracy and reduced computational complexity. The SESSA optimization selected only 1 K out of 25 K features extracted with VGGNet, while improving accuracy at the same time. The results are among the best achieved on these datasets and outperform several convolutional network models. We expect that the combination of CNN feature extraction and SESSA feature optimization could be useful for many other image classification tasks.
In the framework of the presented doctoral thesis, the plant ubiquitous, non-selective vacuolar cation channel TPC1/SV was electrophysiologically studied in Arabidopsis thaliana mesophyll vacuoles to further enlighten its physiological role in plant stress responses. For this, the hyperactive channel version fou2 (D454N), gaining a non-functional vacuolar calcium sensor, strong retarded growth phenotype and upregulated JA signalling pathway, and eight fou2 reverting WT-like ouf mutants were used. Except of ouf4, all other seven ouf mutants carried a 2nd mutation in the TPC1 gene. Therefore, the TPC1 electrical features of all ouf mutants were electrophysiologically characterized with the patch clamp method and compared with fou2 and WT.
Due to a missense mutation, ouf1 and ouf7 mutants harboured a truncated TPC1 channel protein, resulting in an impaired protein integrity and in turn loss of TPC1 channel activity. Accordingly, ouf1 and ouf7 mimicked the tpc1-2 null mutant with a WT- rather fou2-like phenotype. The ouf2 (G583D D454N) mutant exhibited inactive TPC1 channels, probably because the G583D mutation located in luminal part of the S11 helix caused (i) a shift of the activation threshold to much more positive voltages (i.e. to more than +110 mV) (ii) or channel blockage. As a result of the TPC1 channel inactivity, the ouf2 mutant also imitates the WT-like phenotype of the tpc1-2 null mutant. In the ouf6 mutant (A669V D454N) the 2nd reverting mutation selectively influenced fou2-like SV channel features. Both, the fast activation kinetics and reduced luminal calcium sensitivity were similar in ouf6 and fou2. However, deviations in both, the relative and absolute open channel probability, resulted in strongly reduced (80 %) current density at 0 mM and channel inactivity in the voltage range between -30 mV to +40 mV compared to fou2 and WT. Furthermore, the TPC1 channels in ouf6 exhibited a higher susceptibility to inhibitory luminal Ca2+ than fou2. As a result of these different effects, the TPC1 channel activity almost vanished at high luminal Ca2+ loads, what is very likely the reason that ouf6 lost the fou2-like phenotype. The ouf4 mutation did not change the fou2 TPC1-channel features like fast channel activation, single channel conductance and voltage-dependent gating behaviour. Nevertheless, the TPC1 current density was 80% less in ouf4 than in fou2. Since the TPC1 gene was not the target of the 2nd mutation, it can be assumed that it is modulated via external, yet unknown factor. In the ouf8 mutant the TPC1 channels additionally possess M629I mutation within the selectivity filter II resulting in a 50% decrease in the TPC1 unitary conductance. However, the slightly increased relative open channel probability of the TPC1 channels in ouf8 compared to fou2 appeared to be sufficient to compensate the reduced transport capacity of individual TPC1 channels. As a result, a similar macroscopic outward current density of ouf8 and fou2 was detected in the absence of vacuolar Ca2+. Furthermore, ouf8 mutation did not drastically change the typical fou2 TPC1 channel features such as fast activation, vacuolar calcium insensitivity and voltage dependency. However, a reversible block of the cytosol-directed potassium efflux at increased vacuolar calcium concentration in ouf8 mutant was found. Further inspection of transiently expressed TPC1 channel variants (M629I, M629T) on the single channel level suggest that Met629 of AtTPC1 in the channel pore region is crucial for the unitary channel conductance.
Taken together, current membrane recordings from ouf mutants revealed one common feature: All of them lacked or showed a strongly impaired ability for TPC1-mediated potassium release from the vacuole into the cytosol. Additionally, considering the detected dependence of the vacuolar membrane voltage on TPC1 activity, it thus seems that the TPC1-triggered vacuolar membrane depolarization caused by vacuolar K+ release plays a key role in generation of the fou2-like phenotype. Accordingly, one can conclude that TPC1-dependent vacuolar membrane depolarization and initiation of jasmonate production are likely linked. This statement is supported also by the complete restoration of WT-like plant phenotype and JA signalling in the ouf mutants. Finally, as a control element of the vacuolar membrane voltage TPC1 is probably upstream located in JA signalling pathway and therefore a perfect junction for linking multiple physiological stimuli and response to them.
Im Rahmen der vorgelegten Doktorarbeit wurde der in Pflanzen ubiquitär exprimierte, nicht-selektive vakuoläre Kationenkanal TPC1/SV elektrophysiologisch in Arabidopsis thaliana Mesophyllvakuolen untersucht, um seine physiologische Rolle in der pflanzlichen Stressantwort weiter aufzuklären. Hierfür wurde die hyperaktive Kanalvariante fou2 (D454N), die einen nicht-funktionalen vakuolären Calciumsensor, ein stark verzögertes Pflanzenwachstum und einen hochregulierten Jasmonsäure-Signalweg aufweist, sowie acht ouf Mutanten mit fou2-umkehrenden Phänotyp benutzt. Mit Ausnahme von ouf4 enthalten alle anderen ouf Mutanten eine weitere Mutation im TPC1-Gen. Daher wurden die elektrischen Eigenschaften von TPC1 in allen ouf Mutanten elektrophysiologisch mittels der Patch clamp Technik charakterisiert und mit fou2 und dem Wildtyp verglichen.
Aufgrund einer Missense-Mutation beinhalten die Mutanten ouf1 und ouf7 ein verkürztes TPC1 Protein, woraus eine gestörte Proteinintegrität resultiert und daraus wiederum ein Fehlen der TCP1-Kanalaktivität. Dementsprechend ähneln ouf1 und ouf7 der tpc1-2 Nullmutante mit einem WT- oder eher fou2-artigen Phänotyp. Wahrscheinlich weist die ouf2 (G583D D454N) Mutante einen inaktiven TPC1-Kanal auf, weil die G583D Mutation, die in einem luminalen Teil der S11 Helix sitzt, eine Verschiebung der Aktivierungsschwelle hin zu einer höheren Spannung (z. B. mehr als +110 mV) oder einen Kanalblock verursacht. Als Folge der TPC1 Kanal Inaktivität, ahmt die ouf2 Mutante auch den WT-ähnlichen Phänotyp der tpc1-2 Nullmutante nach. In der ouf6 Mutante (A669V D454N) beeinflusst die zweite Mutation selektiv die fou2-ähnlichen SV-Kanaleigenschaften. Sowohl die schnelle Aktivierungskinetik als auch die verringerte luminale Calciumsensitivität waren denen von ouf6 und fou2 ähnlich. Die Abweichungen in der relativen sowie der absoluten Offenwahrscheinlichkeit resultierten jedoch in einer stark reduzierten (80 %) Stromdichte bei 0 mM luminalem Calcium verglichen mit fou2 und dem WT, sowie einer Kanalinaktivität bei Spannungen zwischen -30 mV und +40 mV. Darüber hinaus zeigten die TPC1 Kanäle in ouf6 eine höhere Anfälligkeit für inhibitorisches, luminales Calcium als die in fou2. Das Ergebnis der beiden unterschiedlichen Effekte ist, dass die TPC1 Kanalaktivität bei einer hohen luminalen Calciumkonzentration fast verschwindet, woraus zu schließen ist, dass ouf6 den fou2-ähnlichen Phänotyp verlor. Die ouf4 Mutation veränderte nicht die fou2 TPC1 Kanaleigenschaften, wie die schnelle Kanalaktivierung, die Einzelkanalleitfähigkeit und das spannungsabhängige Verhalten. Nichtsdestotrotz war die TCP1 Stromdichte in ouf4 um 80 % geringer als in fou2. Da das TPC1 Gen nicht das Ziel der zweiten Mutation war, kann angenommen werden, dass es durch äußere, bisher noch unbekannte Faktoren, reguliert wird. In der ouf8 Mutante haben die TPC1 Kanäle zusätzlich eine M629I Mutation innerhalb des zweiten Selektivitätsfilters, welche in einem 50 % Rückgang der TCP1 Einzelkanalleitfähigkeit resultiert. Jedoch scheint die leicht erhöhte Offenwahrscheinlichkeit der TCP1 Kanäle in ouf8, verglichen mit fou2, ausreichend zu sein, um die reduzierte Transportkapazität der individuellen TPC1 Kanäle zu kompensieren. Schlussfolgernd wurde eine ähnliche makroskopische auswärts gerichtete Stromdichte des ouf8 und des fou2 in Abwesenheit vakuolären Calciums entdeckt. Des Weiteren änderte eine ouf8 Mutation die fou2 TPC1 Kanaleigenschaften wie eine schnelle Aktivierung, vakuoläre Calciuminsensitivität und die Spannungsabhängigkeit nicht drastisch. Jedoch wurde ein reversibler Block des Zytosol-gerichteten Kalium Ausstroms bei erhöhten vakuolären Calcium Konzentrationen in ouf8 gefunden. Eine weitere Betrachtung transient exprimierter TPC1 Kanalvarianten (M629I, M629T) auf Einzelkanalebene weist darauf hin, dass das Met629 des AtTPC1 in der Kanalporenregion entscheidend ist für die Einzelkanalleitfähigkeit.
Zusammengefasst zeigt der über die Membran von ouf Mutanten gemessene Strom eine Gemeinsamkeit: Alle zeigten keinen oder einen stark beeinträchtigten TPC1-vermittelten Kaliumausstrom aus der Vakuole ins Zytosol. Unter Berücksichtigung der beobachteten Abhängigkeit der vakuolären Membranspannung von der TPC1 Aktivität, scheint es, als ob die durch TPC1 angeregte Depolarisation der Vakuolenmembran, welche durch die vakuoläre Kaliumfreisetzung bedingt wird, in der Ausbildung des fou2 Phänotyps eine Rolle spielt. Daraus lässt sich ableiten, dass die TPC1-abhängige Depolarisation der Vakuolenmembran und die Jasmonat Bildung vermutlich verbunden sind. Diese Behauptung wird auch gestützt durch die komplette Wiederherstellung des WT-ähnlichen Pflanzenphänotyps und des Jasmonsäure Signalwegs in den ouf Mutanten. Letztendlich ist TPC1 als kontrollierendes Element der vakuolären Membranspannung wahrscheinlich dem Jasmonsäure Signalweg vorgeschaltet und deswegen ein perfekter Knotenpunkt, der verschiedene physiologische Stimuli und ihre Antworten verbindet.
Bees need food of appropriate nutritional quality to maintain their metabolic functions. They largely obtain all required nutrients from floral resources, i.e., pollen and nectar. However, the diversity, composition and nutritional quality of floral resources varies with the surrounding environment and can be strongly altered in human-impacted habitats. We investigated whether differences in plant species richness as found in the surrounding environment correlated with variation in the floral diversity and nutritional quality of larval provisions (i.e., mixtures of pollen, nectar and salivary secretions) composed by the mass-provisioning stingless bee Tetragonula carbonaria (Apidae: Meliponini). We found that the floral diversity of larval provisions increased with increasing plant species richness. The sucrose and fat (total fatty acid) content and the proportion and concentration of the omega-6 fatty acid linoleic acid decreased, whereas the proportion of the omega-3 fatty acid linolenic acid increased with increasing plant species richness. Protein (total amino acid) content and amino acid composition did not change. The protein to fat (P:F) ratio, known to affect bee foraging, increased on average by more than 40% from plantations to forests and gardens, while the omega-6:3 ratio, known to negatively affect cognitive performance, decreased with increasing plant species richness. Our results suggest that plant species richness may support T. carbonaria colonies by providing not only a continuous resource supply (as shown in a previous study), but also floral resources of high nutritional quality.
How genomic and ecological traits shape island biodiversity - insights from individual-based models
(2020)
Life on oceanic islands provides a playground and comparably easy\-/studied basis
for the understanding of biodiversity in general. Island biota feature many
fascinating patterns: endemic species, species radiations and species with
peculiar trait syndromes. However, classic and current island biogeography
theory does not yet consider all the factors necessary to explain many of these
patterns. In response to this, there is currently a shift in island biogeography
research to systematically consider species traits and thus gain a more
functional perspective. Despite this recent development, a set of species
characteristics remains largely ignored in island biogeography, namely genomic
traits. Evidence suggests that genomic factors could explain many of the
speciation and adaptation patterns found in nature and thus may be highly
informative to explain the fascinating and iconic phenomena known for oceanic
islands, including species radiations and susceptibility to biotic invasions.
Unfortunately, the current lack of comprehensive meaningful data makes studying
these factors challenging. Even with paleontological data and space-for-time
rationales, data is bound to be incomplete due to the very environmental
processes taking place on oceanic islands, such as land slides and volcanism,
and lacks causal information due to the focus on correlative approaches. As
promising alternative, integrative mechanistic models can explicitly consider
essential underlying eco\-/evolutionary mechanisms. In fact, these models have
shown to be applicable to a variety of different systems and study questions.
In this thesis, I therefore examined present mechanistic island models to
identify how they might be used to address some of the current open questions in
island biodiversity research. Since none of the models simultaneously considered
speciation and adaptation at a genomic level, I developed a new genome- and
niche-explicit, individual-based model. I used this model to address three
different phenomena of island biodiversity: environmental variation, insular
species radiations and species invasions.
Using only a single model I could show that small-bodied species with flexible
genomes are successful under environmental variation, that a complex combination
of dispersal abilities, reproductive strategies and genomic traits affect the
occurrence of species radiations and that invasions are primarily driven by the
intensity of introductions and the trait characteristics of invasive
species. This highlights how the consideration of functional traits can promote
the understanding of some of the understudied phenomena in island biodiversity.
The results presented in this thesis exemplify the generality of integrative
models which are built on first principles. Thus, by applying such models to
various complex study questions, they are able to unveil multiple biodiversity
dynamics and patterns. The combination of several models such as the one I
developed to an eco\-/evolutionary model ensemble could further help to identify
fundamental eco\-/evolutionary principles. I conclude the thesis with an outlook
on how to use and extend my developed model to investigate geomorphological
dynamics in archipelagos and to allow dynamic genomes, which would further
increase the model's generality.
Understanding extinction debts: spatio-temporal scales, mechanisms and a roadmap for future research
(2019)
Extinction debt refers to delayed species extinctions expected as a consequence of ecosystem perturbation. Quantifying such extinctions and investigating long‐term consequences of perturbations has proven challenging, because perturbations are not isolated and occur across various spatial and temporal scales, from local habitat losses to global warming. Additionally, the relative importance of eco‐evolutionary processes varies across scales, because levels of ecological organization, i.e. individuals, (meta)populations and (meta)communities, respond hierarchically to perturbations. To summarize our current knowledge of the scales and mechanisms influencing extinction debts, we reviewed recent empirical, theoretical and methodological studies addressing either the spatio–temporal scales of extinction debts or the eco‐evolutionary mechanisms delaying extinctions. Extinction debts were detected across a range of ecosystems and taxonomic groups, with estimates ranging from 9 to 90% of current species richness. The duration over which debts have been sustained varies from 5 to 570 yr, and projections of the total period required to settle a debt can extend to 1000 yr. Reported causes of delayed extinctions are 1) life‐history traits that prolong individual survival, and 2) population and metapopulation dynamics that maintain populations under deteriorated conditions. Other potential factors that may extend survival time such as microevolutionary dynamics, or delayed extinctions of interaction partners, have rarely been analyzed. Therefore, we propose a roadmap for future research with three key avenues: 1) the microevolutionary dynamics of extinction processes, 2) the disjunctive loss of interacting species and 3) the impact of multiple regimes of perturbation on the payment of debts. For their ability to integrate processes occurring at different levels of ecological organization, we highlight mechanistic simulation models as tools to address these knowledge gaps and to deepen our understanding of extinction dynamics.
Although the concept of botanical carnivory has been known since Darwin's time, the molecular mechanisms that allow animal feeding remain unknown, primarily due to a complete lack of genomic information. Here, we show that the transcriptomic landscape of the Dionaea trap is dramatically shifted toward signal transduction and nutrient transport upon insect feeding, with touch hormone signaling and protein secretion prevailing. At the same time, a massive induction of general defense responses is accompanied by the repression of cell death-related genes/processes. We hypothesize that the carnivory syndrome of Dionaea evolved by exaptation of ancient defense pathways, replacing cell death with nutrient acquisition.
Experimental high-throughput analysis of molecular networks is a central approach to characterize the adaptation of plant metabolism to the environment. However, recent studies have demonstrated that it is hardly possible to predict in situ metabolic phenotypes from experiments under controlled conditions, such as growth chambers or greenhouses. This is particularly due to the high molecular variance of in situ samples induced by environmental fluctuations. An approach of functional metabolome interpretation of field samples would be desirable in order to be able to identify and trace back the impact of environmental changes on plant metabolism. To test the applicability of metabolomics studies for a characterization of plant populations in the field, we have identified and analyzed in situ samples of nearby grown natural populations of Arabidopsis thaliana in Austria. A. thaliana is the primary molecular biological model system in plant biology with one of the best functionally annotated genomes representing a reference system for all other plant genome projects. The genomes of these novel natural populations were sequenced and phylogenetically compared to a comprehensive genome database of A. thaliana ecotypes. Experimental results on primary and secondary metabolite profiling and genotypic variation were functionally integrated by a data mining strategy, which combines statistical output of metabolomics data with genome-derived biochemical pathway reconstruction and metabolic modeling. Correlations of biochemical model predictions and population-specific genetic variation indicated varying strategies of metabolic regulation on a population level which enabled the direct comparison, differentiation, and prediction of metabolic adaptation of the same species to different habitats. These differences were most pronounced at organic and amino acid metabolism as well as at the interface of primary and secondary metabolism and allowed for the direct classification of population-specific metabolic phenotypes within geographically contiguous sampling sites.
Die NFκB-Signalwege, Apoptose und Nekroptose sind essentielle Prozesse in der Immunantwort. Außerdem sind diese Signalwege Teil der Regulation von Zelldifferenzierung, -proliferation, -tod und Entzündungsreaktionen. Dabei wird zuerst der Rezeptor (TNFR1 oder TRAILR 1/2) aktiviert, die rekrutierten DD-Adapterproteine TRADD, FADD und RIPK1 leiten dann die entsprechende Signalkaskade weiter und bestimmen durch ihre Zusammenwirkung, ob der NFκB-Signalweg, Apoptose oder Nekroptose induziert wird.
TNFR1 und TRAILR 1/2 benötigen die DD-Adapterproteine TRADD, FADD und RIPK1 für die Zelltodinduktion, deren konkrete Bedeutung in Bezug auf Rezeptor-Spezifität, Zusammenwirken und Relevanz allerdings noch unklar ist. Um das Zusammenspiel dieser Proteine besser zu verstehen, wurden in dieser Arbeit Nekroptose-kompetente RIPK3-exprimierende HeLa-Zellen verwendet, bei denen die DD-Adapterproteine FADD, TRADD und RIPK1 einzeln oder in Kombination von zweien ausgeknockt wurden. Es stellte sich heraus, dass RIPK1 essentiell für die TNFR1- und TRAILR 1/2-vermittelte Nekroptose-Induktion ist, doch RIPK1 alleine, d.h. ohne FADD- oder TRADD-Mitbeteiligung, nur bei der TNFR1-Nekroptose-Induktion ausreicht. Wiederum inhibiert TRADD die TNFR1- und TRAILR 1/2-induzierte Nekroptose. RIPK1 und TRADD sind aber unverzichtbar für die NFκB-Aktivierung durch TNFR1 oder TRAILR 1/2 und spielen eine wichtige Rolle bei TNFR1-induzierter Apoptose. Andererseits ist FADD alleine ausreichend für die TRAILR 1/2-bezogene Caspase-8 Aktivierung. Zudem ist FADD notwendig für die TRAIL-induzierte NFκB-Signalaktivierung. In Abwesenheit von FADD und TRADD vermittelt RIPK1 die TNF-induzierte Caspase-8 Aktivierung. FADD wird für die TRAIL-induzierte Nekroptose benötigt, aber gegenläufig wirkt die TNF-induzierte Nektroptose in einer Caspase-8 abhängigen und unabhängigen Weise. Zudem sensitiviert TWEAK die TNF- und TRAIL-induzierte Nekroptose.
Zusammenfassend wurde in dieser Arbeit die Auswirkung von TNFR1 und TRAILR 1/2 auf die Aktivierung der unterschiedlichen Signalkaskaden untersucht. Des Weiteren wurde gezeigt, in welcher Weise sich das Zusammenspiel von TRADD, FADD und RIPK1 auf die Induktion von NFκB, Apoptose und Nekroptose auswirkt.
Although many genes have been identified using high throughput technologies in endometriosis (ES), only a small number of individual genes have been analyzed functionally. This is due to the complexity of the disease that has different stages and is affected by various genetic and environmental factors. Many genes are upregulated or downregulated at each stage of the disease, thus making it difficult to identify key genes. In addition, little is known about the differences between the different stages of the disease. We assumed that the study of the identified genes in ES at a system-level can help to better understand the molecular mechanism of the disease at different stages of the development. We used publicly available microarray data containing archived endometrial samples from women with minimal/mild endometriosis (MMES), mild/severe endometriosis (MSES) and without endometriosis. Using weighted gene co-expression analysis (WGCNA), functional modules were derived from normal endometrium (NEM) as the reference sample. Subsequently, we tested whether the topology or connectivity pattern of the modules was preserved in MMES and/or MSES. Common and specific hub genes were identified in non-preserved modules. Accordingly, hub genes were detected in the non-preserved modules at each stage. We identified sixteen co-expression modules. Of the 16 modules, nine were non-preserved in both MMES and MSES whereas five were preserved in NEM, MMES, and MSES. Importantly, two non-preserved modules were found in either MMES or MSES, highlighting differences between the two stages of the disease. Analyzing the hub genes in the non-preserved modules showed that they mostly lost or gained their centrality in NEM after developing the disease into MMES and MSES. The same scenario was observed, when the severeness of the disease switched from MMES to MSES. Interestingly, the expression analysis of the new selected gene candidates including CC2D2A, AEBP1, HOXB6, IER3, and STX18 as well as IGF-1, CYP11A1 and MMP-2 could validate such shifts between different stages. The overrepresented gene ontology (GO) terms were enriched in specific modules, such as genetic disposition, estrogen dependence, progesterone resistance and inflammation, which are known as endometriosis hallmarks. Some modules uncovered novel co-expressed gene clusters that were not previously discovered.