Reactive oxygen species (ROS) are byproducts of aerobic metabolism and may cause oxidative damage to biomolecules. Plants have a complex redox system, involving enzymatic and non-enzymatic compounds. The evolutionary origin of enzymatic antioxidant defense in plants is yet unclear. Here, we describe the redox gene network for A. thaliana and investigate the evolutionary origin of this network. We gathered from public repositories 246 A. thaliana genes directly involved with ROS metabolism and proposed an A. thaliana redox gene network. Using orthology information of 238 Eukaryotes from STRINGdb, we inferred the evolutionary root of each gene to reconstruct the evolutionary history of A. thaliana antioxidant gene network. We found two interconnected clusters: one formed by SOD-related, Thiol-redox, peroxidases, and other oxido-reductase; and the other formed entirely by class III peroxidases. Each cluster emerged in different periods of evolution: the cluster formed by SOD-related, Thiol-redox, peroxidases, and other oxido-reductase emerged before opisthokonta-plant divergence; the cluster composed by class III peroxidases emerged after opisthokonta-plant divergence and therefore contained the most recent network components. According to our results, class III peroxidases are in expansion throughout plant evolution, with new orthologs emerging in each evaluated plant clade divergence.

Ewing Sarcoma (ES) is a rare malignant tumor occurring most frequently in adolescents and young adults. The ES hallmark is a chromosomal translocation between the chromosomes 11 and 22 that results in an aberrant transcription factor (TF) through the fusion of genes from the FET and ETS families, commonly EWSR1 and FLI1. The regulatory mechanisms behind the ES transcriptional alterations remain poorly understood. Here, we reconstruct the ES regulatory network using public available transcriptional data. Seven TFs were identified as potential MRs and clustered into two groups: one composed by PAX7 and RUNX3, and another composed by ARNT2, CREB3L1, GLI3, MEF2C, and PBX3. The MRs within each cluster act as reciprocal agonists regarding the regulation of shared genes, regulon activity, and implications in clinical outcome, while the clusters counteract each other. The regulons of all the seven MRs were differentially methylated. PAX7 and RUNX3 regulon activity were associated with good prognosis while ARNT2, CREB3L1, GLI3, and PBX3 were associated with bad prognosis. PAX7 and RUNX3 appear as highly expressed in ES biopsies and ES cell lines. This work contributes to the understanding of the ES regulome, identifying candidate MRs, analyzing their methilome and pointing to potential prognostic factors.

Sepsis remains a leading cause of death in ICUs all over the world, with pediatric sepsis accounting for a high percentage of mortality in pediatric ICUs. Its complexity makes it difficult to establish a consensus on genetic biomarkers and therapeutic targets. A promising strategy is to investigate the regulatory mechanisms involved in sepsis progression, but there are few studies regarding gene regulation in sepsis. This work aimed to reconstruct the sepsis regulatory network and identify transcription factors (TFs) driving transcriptional states, which we refer to here as master regulators. We used public gene expression datasets to infer the co-expression network associated with sepsis in a retrospective study. We identified a set of 15 TFs as potential master regulators of pediatric sepsis, which were divided into two main clusters. The first cluster corresponded to TFs with decreased activity in pediatric sepsis, and GATA3 and RORA, as well as other TFs previously implicated in the context of inflammatory response. The second cluster corresponded to TFs with increased activity in pediatric sepsis and was composed of TRIM25, RFX2, and MEF2A, genes not previously described as acting in a coordinated way in pediatric sepsis. Altogether, these results show how a subset of master regulators TF can drive pathological transcriptional states, with implications for sepsis biology and treatment.

Lead poisoning effects are wide and include nervous system impairment, peculiarly during development, leading to neural damage. Lead interaction with calcium and zinc-containing metalloproteins broadly affects cellular metabolism since these proteins are related to intracellular ion balance, activation of signaling transduction cascades, and gene expression regulation. In spite of lead being recognized as a neurotoxin, there are gaps in knowledge about the global effect of lead in modulating the transcription of entire cellular systems in neural cells. In order to investigate the effects of lead poisoning in a systemic perspective, we applied the transcriptogram methodology in an RNA-seq dataset of human embryonic-derived neural progenitor cells (ES-NP cells) treated with 30 µM lead acetate for 26 days. We observed early downregulation of several cellular systems involved with cell differentiation, such as cytoskeleton organization, RNA, and protein biosynthesis. The downregulated cellular systems presented big and tightly connected networks. For long treatment times (12 to 26 days), it was possible to observe a massive impairment in cell transcription profile. Taking the enriched terms together, we observed interference in all layers of gene expression regulation, from chromatin remodeling to vesicle transport. Considering that ES-NP cells are progenitor cells that can originate other neural cell types, our results suggest that lead-induced gene expression disturbance might impair cells’ ability to differentiate, therefore influencing ES-NP cells’ fate.

The origin of nervous systems is a main theme in biology and its mechanisms are largely underlied by synaptic neurotransmission. One problem to explain synapse establishment is that synaptic orthologs are present in multiple aneural organisms. We questioned how the interactions among these elements evolved and to what extent it relates to our understanding of the nervous systems complexity. We identified the human neurotransmission gene network based on genes present in GABAergic, glutamatergic, serotonergic, dopaminergic, and cholinergic systems. The network comprises 321 human genes, 83 of which act exclusively in the nervous system. We reconstructed the evolutionary scenario of synapse emergence by looking for synaptic orthologs in 476 eukaryotes. The Human-Cnidaria common ancestor displayed a massive emergence of neuroexclusive genes, mainly ionotropic receptors, which might have been crucial to the evolution of synapses. Very few synaptic genes had their origin after the Human-Cnidaria common ancestor. We also identified a higher abundance of synaptic proteins in vertebrates, which suggests an increase in the synaptic network complexity of those organisms.