Biomarker potential of competing endogenous RNA networks in Polycystic Ovary Syndrome (PCOS)

Polycystic ovary syndrome (PCOS) is the most common condition affecting women of reproductive age globally. PCOS continues to be the largest contributing factor to female infertility despite significant progress in our knowledge of the molecular underpinnings and treatment of the condition. The fact that PCOS is a very diverse condition makes it one of the key reasons why we haven't been able to overcome it. Non-coding RNAs (ncRNAs) are implicated in the development of PCOS, according to growing evidence. However, it is unclear how the complex regulatory relationships between the many ncRNA types contribute to the growth of this malignancy. Competing endogenous RNA (ceRNA), a recently identified mechanism in the RNA world, suggests regulatory interactions between various RNAs, including long non-coding RNAs (lncRNAs), microRNAs (miRNAs), transcribed pseudogenes, and circular RNAs (circRNAs). Recent studies on PCOS have shown that dysregulation of multiple ceRNA networks (ceRNETs) between these ncRNAs plays crucial roles in developing the defining characteristics of PCOS development. And it is believed that such a finding may open a new door for a deeper comprehension of PCOS's unexplored facets. In addition, it may be able to provide fresh biomarkers and effective therapy targets for PCOS. This review will go over the body of information that exists about the primary roles of ceRNETs before highlighting the developing involvement of several newly found ceRNETs in a number of PCOS characteristics.

sylvestris

Transcriptional regulatory networks are pivotal components of plant's response to salt stress. However, plant adaptation strategies varied as a function of stress intensity, which is mainly modulated by climate change. Here, we determined the gene regulatory networks based on transcription factor (TF) TF_gene co-expression, using two transcriptomic data sets generated from the salt-tolerant "Tebaba" roots either treated with 50 mM NaCl (mild stress) or 150 mM NaCl (severe stress). The analysis of these regulatory networks identified specific TFs as key regulatory hubs as evidenced by their multiple interactions with different target genes related to stress response. Indeed, under mild stress, NAC and bHLH TFs were identified as central hubs regulating nitrogen storage process. Moreover, HSF TFs were revealed as a regulatory hub regulating various aspects of cellular metabolism including flavonoid biosynthesis, protein processing, phenylpropanoid metabolism, galactose metabolism, and heat shock proteins. These processes are essentially linked to short-term acclimatization under mild salt stress. This was further consolidated by the protein-protein interaction (PPI) network analysis showing structural and plant growth adjustment. Conversely, under severe salt stress, dramatic metabolic changes were observed leading to novel TF members including MYB family as regulatory hubs controlling isoflavonoid biosynthesis, oxidative stress response, abscisic acid signaling pathway, and proteolysis. The PPI network analysis also revealed deeper stress defense changes aiming to restore plant metabolic homeostasis when facing severe salt stress. Overall, both the gene co-expression and PPI network provided valuable insights on key transcription factor hubs that can be employed as candidates for future genetic crop engineering programs.

HSFA3 functions as a positive regulator of HSFA2a to enhance thermotolerance in perennial ryegrass

Perennial ryegrass (Lolium perenne) is a widely used cool season turfgrass with outstanding turf quality and grazing tolerance. High temperature is the key factor restricting the distribution of perennial ryegrass in temperate and sub-tropic regions. In this study, we found that one HEAT SHCOK TRANSCRIPTION FACOTR (HSF) class A gene from perennial ryegrass, LpHSFA3, was highly induced by heat stress. LpHSFA3 is localized in nucleus and functions as a transcription factor. Ectopic overexpression of LpHSFA3 in Arabidopsis improved thermotolerance and rescued heat sensitive deficiency of athsfa3 mutant. Overexpression of LpHSFA3 in perennial ryegrass enhanced heat tolerance and increased survival rate in summer season as evidenced by decreased EL and MDA, increased number of green leaves and total chlorophyll content. LpHSFA3 binds to the HSE region in LpHSFA2a promoter to constitutively activate the expression of LpHSFA2a and downstream heat stress responsive genes. Ectopic overexpression of LpHSFA2a consequently rescued thermal sensitivity of athsfa3 mutant and enhanced thermotolerance of athsfa2 mutant. Perennial ryegrass protoplasts with overexpression of LpHSFA3 and LpHSFA2a exhibited induction of similar subsets of heat responsive genes. These results indicated that transcription factor LpHSFA3 functions as positive regulator of LpHSFA2a to improve thermotolerance of perennial ryegrass, providing further evidence to understand the regulatory networks of plant heat stress response.

Harnessing regulatory networks in Actinobacteria for natural product discovery

Microbes typically live in complex habitats where they need to rapidly adapt to continuously changing growth conditions. To do so, they produce an astonishing array of natural products with diverse structures and functions. Actinobacteria stand out for their prolific production of bioactive molecules, including antibiotics, anticancer agents, antifungals, and immunosuppressants. Attention has been directed especially towards the identification of the compounds they produce and the mining of the large diversity of biosynthetic gene clusters (BGCs) in their genomes. However, the current return on investment in random screening for bioactive compounds is low, while it is hard to predict which of the millions of BGCs should be prioritized. Moreover, many of the BGCs for yet undiscovered natural products are silent or cryptic under laboratory growth conditions. To identify ways to prioritize and activate these BGCs, knowledge regarding the way their expression is controlled is crucial. Intricate regulatory networks control global gene expression in Actinobacteria, governed by a staggering number of up to 1000 transcription factors per strain. This review highlights recent advances in experimental and computational methods for characterizing and predicting transcription factor binding sites and their applications to guide natural product discovery. We propose that regulation-guided genome mining approaches will open new avenues toward eliciting the expression of BGCs, as well as prioritizing subsets of BGCs for expression using synthetic biology approaches.

ONE-SENTENCE SUMMARY

This review provides insights into advances in experimental and computational methods aimed at predicting transcription factor binding sites and their applications to guide natural product discovery.

Genome-wide identification of potassium channels in maize showed evolutionary patterns and variable functional responses to abiotic stresses

Potassium (K) channels are essential components of plant biology, mediating not only K ion (K+) homeostasis but also regulating several physiological processes and stress tolerance. In the current investigation, we identified 27 K+ channels in maize and deciphered the evolution and divergence pattern with four monocots and five dicot species. Chromosomal localization and expansion of K+ channel genes showed uneven distribution and were independent of genome size. The dispersed duplication is the major force in expanding K+ channels in the target genomes. The mean Ka/Ks ratio of <0.5 in paralogs and orthologs indicates horizontal and vertical expansions of K+ channel genes under strong purifying selection. The one-to-one K+ channel orthologs were prominent among the closely related species, with higher synteny between maize and the rest of the monocots. Comprehensive K+ channels promoter analysis revealed various cis-regulatory elements mediating stress tolerance with the predominance of MYB and STRE binding sites. The regulatory network showed AP2-EREBP TFs, miR164 and miR399 are prominent regulatory elements of K+ channels. The qRT-PCR analysis of K+ channels and regulatory miRNAs showed significant expressions in response to drought and waterlogging stresses. The present study expanded the knowledge on K+ channels in maize and will serve as a basis for an in-depth functional analysis.

Hsa_circITGA4/ miR-1468/EGFR/ PTEN a Master Regulators Axis in Glioblastoma Development and Progression

In the fight against glioblastoma, circular RNA is emerging as a functional molecule. However, how circular RNA (circRNA) is regulated and what role it plays is still a mystery. In this research, different bioinformatics approaches were used to evaluate glioblastoma circRNA sequencing and array data, with the goal of developing a putative molecular sponge mechanism control network. The circRNAs were obtained from the Gene Expression Omnibus datasets. MicroRNA-circRNA interactions were predicted using CircInteractome. The microRNAs' expression and survival trends were screened using the TCGA database. MicroRNA gene targets were predicted using the MiRnet database. Sponge network gene candidates were screened using data from the GEPIA. The roles of the targeted genes were to be explained by analyzing data from Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes. To build the network and display the outcomes, we utilized python program, and enrichment online Bioinformatics databases. The circRNAs hsa_circITGA4_002, hsa_circITGA4_001, hsa_circITGA4_003, hsa_circ_0030855, hsa_circ_0030857 were chosen from among GBM patients and control group. Upregulation of hsa-miR-1468, hsa-miR-3683, hsa-miR-1273c, and hsa-miR-4665-3p were associated with a poor prognosis in GBM. MicroRNA targets such as ITGA4, LAMA2, EGFR, PTEN, COL1A4, and NCAM2 were analyzed using expression and survival data. The Apoptosis, cell adhesion molecules, PI3K/AKT and P53 signaling pathways were the most abundant functional categories among gene targets. The circRNA molecular sponge regulatory network includes hsa-miR-1468 and hsa-miR-4665-3p. In this network, hs hsa_circITGA4_002, hsa_circITGA4_001, hsa_circ_0030857, EGFR, PTEN, and ITGA4 may represent GBM therapeutic targets. Their role in GBM needs additional study.

Regulation and signaling pathways in cancer stem cells: implications for targeted therapy for cancer

Cancer stem cells (CSCs), initially identified in leukemia in 1994, constitute a distinct subset of tumor cells characterized by surface markers such as CD133, CD44, and ALDH. Their behavior is regulated through a complex interplay of networks, including transcriptional, post-transcriptional, epigenetic, tumor microenvironment (TME), and epithelial-mesenchymal transition (EMT) factors. Numerous signaling pathways were found to be involved in the regulatory network of CSCs. The maintenance of CSC characteristics plays a pivotal role in driving CSC-associated tumor metastasis and conferring resistance to therapy. Consequently, CSCs have emerged as promising targets in cancer treatment. To date, researchers have developed several anticancer agents tailored to specifically target CSCs, with some of these treatment strategies currently undergoing preclinical or clinical trials. In this review, we outline the origin and biological characteristics of CSCs, explore the regulatory networks governing CSCs, discuss the signaling pathways implicated in these networks, and investigate the influential factors contributing to therapy resistance in CSCs. Finally, we offer insights into preclinical and clinical agents designed to eliminate CSCs.

The mechanism of ferroptosis and its related diseases

Ferroptosis, a regulated form of cellular death characterized by the iron-mediated accumulation of lipid peroxides, provides a novel avenue for delving into the intersection of cellular metabolism, oxidative stress, and disease pathology. We have witnessed a mounting fascination with ferroptosis, attributed to its pivotal roles across diverse physiological and pathological conditions including developmental processes, metabolic dynamics, oncogenic pathways, neurodegenerative cascades, and traumatic tissue injuries. By unraveling the intricate underpinnings of the molecular machinery, pivotal contributors, intricate signaling conduits, and regulatory networks governing ferroptosis, researchers aim to bridge the gap between the intricacies of this unique mode of cellular death and its multifaceted implications for health and disease. In light of the rapidly advancing landscape of ferroptosis research, we present a comprehensive review aiming at the extensive implications of ferroptosis in the origins and progress of human diseases. This review concludes with a careful analysis of potential treatment approaches carefully designed to either inhibit or promote ferroptosis. Additionally, we have succinctly summarized the potential therapeutic targets and compounds that hold promise in targeting ferroptosis within various diseases. This pivotal facet underscores the burgeoning possibilities for manipulating ferroptosis as a therapeutic strategy. In summary, this review enriched the insights of both investigators and practitioners, while fostering an elevated comprehension of ferroptosis and its latent translational utilities. By revealing the basic processes and investigating treatment possibilities, this review provides a crucial resource for scientists and medical practitioners, aiding in a deep understanding of ferroptosis and its effects in various disease situations.

Integrative single-nucleus multi-omics analysis prioritizes candidate cis and trans regulatory networks and their target genes in Alzheimer's disease brains

BACKGROUND

The genetic underpinnings of late-onset Alzheimer's disease (LOAD) are yet to be fully elucidated. Although numerous LOAD-associated loci have been discovered, the causal variants and their target genes remain largely unknown. Since the brain is composed of heterogenous cell subtypes, it is imperative to study the brain on a cell subtype specific level to explore the biological processes underlying LOAD.

METHODS

Here, we present the largest parallel single-nucleus (sn) multi-omics study to simultaneously profile gene expression (snRNA-seq) and chromatin accessibility (snATAC-seq) to date, using nuclei from 12 normal and 12 LOAD brains. We identified cell subtype clusters based on gene expression and chromatin accessibility profiles and characterized cell subtype-specific LOAD-associated differentially expressed genes (DEGs), differentially accessible peaks (DAPs) and cis co-accessibility networks (CCANs).

RESULTS

Integrative analysis defined disease-relevant CCANs in multiple cell subtypes and discovered LOAD-associated cell subtype-specific candidate cis regulatory elements (cCREs), their candidate target genes, and trans-interacting transcription factors (TFs), some of which, including ELK1, JUN, and SMAD4 in excitatory neurons, were also LOAD-DEGs. Finally, we focused on a subset of cell subtype-specific CCANs that overlap known LOAD-GWAS regions and catalogued putative functional SNPs changing the affinities of TF motifs within LOAD-cCREs linked to LOAD-DEGs, including APOE and MYO1E in a specific subtype of microglia and BIN1 in a subpopulation of oligodendrocytes.

CONCLUSIONS

To our knowledge, this study represents the most comprehensive systematic interrogation to date of regulatory networks and the impact of genetic variants on gene dysregulation in LOAD at a cell subtype resolution. Our findings reveal crosstalk between epigenetic, genomic, and transcriptomic determinants of LOAD pathogenesis and define catalogues of candidate genes, cCREs, and variants involved in LOAD genetic etiology and the cell subtypes in which they act to exert their pathogenic effects. Overall, these results suggest that cell subtype-specific cis-trans interactions between regulatory elements and TFs, and the genes dysregulated by these networks contribute to the development of LOAD.

Differentially expressed non-coding RNAs and their regulatory networks in liver cancer

The vast majority of human transcriptome is represented by various types of small RNAs with little or no protein-coding capability referred to as non-coding RNAs (ncRNAs). Functional ncRNAs include microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs), which are expressed at very low, but stable and reproducible levels in a variety of cell types. ncRNAs regulate gene expression due to miRNA capability of complementary base pairing with mRNAs, whereas lncRNAs and circRNAs can sponge miRNAs off their target mRNAs to act as competitive endogenous RNAs (ceRNAs). Each miRNA can target multiple mRNAs and a single mRNA can interact with several miRNAs, thereby creating miRNA-mRNA, lncRNA-miRNA-mRNA, and circRNA-miRNA-mRNA regulatory networks. Over the past few years, a variety of differentially expressed miRNAs, lncRNAs, and circRNAs (DEMs, DELs, and DECs, respectively) have been linked to cancer pathogenesis. They can exert both oncogenic and tumor suppressor roles. In this review, we discuss the recent advancements in uncovering the roles of DEMs, DELs, and DECs and their networks in aberrant cell signaling, cell cycle, transcription, angiogenesis, and apoptosis, as well as tumor microenvironment remodeling and metabolic reprogramming during hepatocarcinogenesis. We highlight the potential and challenges in the use of differentially expressed ncRNAs as biomarkers for liver cancer diagnosis and prognosis.

Phenotype Control techniques for Boolean gene regulatory networks

Modeling cell signal transduction pathways via Boolean networks (BNs) has become an established method for analyzing intracellular communications over the last few decades. What's more, BNs provide a course-grained approach, not only to understanding molecular communications, but also for targeting pathway components that alter the long-term outcomes of the system. This has come to be known as phenotype control theory. In this review we study the interplay of various approaches for controlling gene regulatory networks such as: algebraic methods, control kernel, feedback vertex set, and stable motifs. The study will also include comparative discussion between the methods, using an established cancer model of T-Cell Large Granular Lymphocyte Leukemia. Further, we explore possible options for making the control search more efficient using reduction and modularity. Finally, we will include challenges presented such as the complexity and the availability of software for implementing each of these control techniques.

The regulation of expression and splicing of transcription factors are related to the muscle damage caused by radiation in tree shrews

Radiotherapy-induced muscle injury (RIMI) is a major complication of radiotherapy for nasopharyngeal carcinoma. Transcription factor (TF) expression and alternative splicing are crucial events in transcriptional and posttranscriptional regulation, respectively, and are known to be involved in key signaling pathways contributing to a variety of human disorders, including radiation injury. To investigate the TFs and alternative splicing events involved in RIMI, we constructed a tree shrew model as described previously in which the RIMI group received 20 Gy of irradiation on the tensor veli palatini (TVP) muscles. The irradiated muscles were evaluated by RNA sequencing (RNA-seq) 6 months later, and the results compared with those for normal TVP muscles. The alt5p and alt3p events were the two main types of differentially regulated alternative splicing events (RASEs) identified via the Splice sites Usage Variation Analysis (SUVA) software, and these RASEs were highly conserved in RIMI. According to functional enrichment analysis, the differentially RASEs were primarily enriched in pathways related to transcriptional regulation. Furthermore, we identified 16 alternative splicing TFs (ASTFs) in ASTF-differentially expressed gene (DEG) networks based on co-expression analysis, and the regulatory networks were chiefly enriched in pathways linked to cell proliferation and differentiation. This study revealed that RASEs and ASTF-DEG networks may both play important regulatory roles in gene expression network alteration in RIMI. Future studies on the targeting mechanisms and early interventions directed at RASEs and ASTF-DEG networks may aid in the treatment of RIMI.

MicroRNAs as promising therapeutic agents: A perspective from acupuncture

MicroRNAs (miRNAs) are gaining recognition as potential therapeutic agents due to their small size, ability to target a wide range of genes, and significant role in disease progression. However, despite their promising potential, nearly half of the miRNA drugs developed for therapeutic purposes have been discontinued or put on hold, and none have advanced to phase III clinical trials. The development of miRNA therapeutics has faced obstacles such as difficulties in validating miRNA targets, conflicting evidence regarding competition and saturation effects, challenges in miRNA delivery, and determining appropriate dosages. These hurdles primarily arise from the intricate functional complexity of miRNAs. Acupuncture, a distinct, complementary therapy, offers a promising avenue to overcome these barriers, particularly by addressing the fundamental issue of preserving functional complexity through acupuncture regulatory networks. The acupuncture regulatory network consists of three main components: the acupoint network, the neuro-endocrine-immune (NEI) network, and the disease network. These networks represent the processes of information transformation, amplification, and conduction that occur during acupuncture. Notably, miRNAs serve as essential mediators and shared biological language within these interconnected networks. Harnessing the therapeutic potential of acupuncture-derived miRNAs can help reduce the time and economic resources required for miRNA drug development and alleviate the current developmental challenges miRNA therapeutics face. This review provides an interdisciplinary perspective by summarizing the interactions between miRNAs, their targets, and the three acupuncture regulatory networks mentioned earlier. The aim is to illuminate the challenges and opportunities in developing miRNA therapeutics. This review paper presents a comprehensive overview of miRNAs, their interactions with acupuncture regulatory networks, and their potential as therapeutic agents. By bridging the miRNA research and acupuncture fields, we aim to offer valuable insights into the obstacles and prospects of developing miRNA therapeutics.

Core conserved transcriptional regulatory networks define the invasive trophoblast cell lineage

The invasive trophoblast cell lineages in rat and human share crucial responsibilities in establishing the uterine-placental interface of the hemochorial placenta. These observations have led to the rat becoming an especially useful animal model for studying hemochorial placentation. However, our understanding of similarities or differences between regulatory mechanisms governing rat and human invasive trophoblast cell populations is limited. In this study, we generated single-nucleus ATAC-seq data from gestation day 15.5 and 19.5 rat uterine-placental interface tissues, and integrated the data with single-cell RNA-seq data generated at the same stages. We determined the chromatin accessibility profiles of invasive trophoblast, natural killer, macrophage, endothelial and smooth muscle cells, and compared invasive trophoblast chromatin accessibility with extravillous trophoblast cell accessibility. In comparing chromatin accessibility profiles between species, we found similarities in patterns of gene regulation and groups of motifs enriched in accessible regions. Finally, we identified a conserved gene regulatory network in invasive trophoblast cells. Our data, findings and analysis will facilitate future studies investigating regulatory mechanisms essential for the invasive trophoblast cell lineage.

Non-coding RNAs identification and regulatory networks in pathogen-host interaction in the microsporidia congenital infection

BACKGROUND

The interaction networks between coding and non-coding RNAs (ncRNAs) including long non-coding RNA (lncRNA), covalently closed circular RNA (circRNA) and miRNA are significant to elucidate molecular processes of biological activities and interactions between host and pathogen. Congenital infection caused by vertical transmission of microsporidia N. bombycis can result in severe economic losses in the silkworm-feeding industry. However, little is known about ncRNAs that take place in the microsporidia congenital infection. Here we conducted whole-transcriptome RNA-Seq analyses to identify ncRNAs and regulatory networks for both N. bombycis and host including silkworm embryos and larvae during the microsporidia congenital infection.

RESULTS

A total of 4,171 mRNAs, 403 lncRNA, 62 circRNAs, and 284 miRNAs encoded by N. bombycis were identified, among which some differentially expressed genes formed cross-talk and are involved in N. bombycis proliferation and infection. For instance, a lncRNA/circRNA competing endogenous RNA (ceRNA) network including 18 lncRNAs, one circRNA, and 20 miRNAs was constructed to describe 14 key parasites genes regulation, such as polar tube protein 3 (PTP3), ricin-B-lectin, spore wall protein 4 (SWP4), and heat shock protein 90 (HSP90). Regarding host silkworm upon N. bombycis congenital infection, a total of 14,889 mRNAs, 3,038 lncRNAs, 19,039 circRNAs, and 3,413 miRNAs were predicted based on silkworm genome with many differentially expressed coding and non-coding genes during distinct developmental stages. Different species of RNAs form interacting network to modulate silkworm biological processes, such as growth, metamorphosis and immune responses. Furthermore, a lncRNA/circRNA ceRNA network consisting of 140 lncRNAs, five circRNA, and seven miRNAs are constructed hypothetically to describe eight key host genes regulation, such as Toll-6, Serpin-6, inducible nitric oxide synthase (iNOS) and Caspase-8. Notably, cross-species analyses indicate that parasite and host miRNAs play a vital role in pathogen-host interaction in the microsporidia congenital infection.

CONCLUSION

This is the first comprehensive pan-transcriptome study inclusive of both N. bombycis and its host silkworm with a specific focus on the microsporidia congenital infection, and show that ncRNA-mediated regulation plays a vital role in the microsporidia congenital infection, which provides a new insight into understanding the basic biology of microsporidia and pathogen-host interaction.

Single and combined effects of polystyrene nanoplastics and Cd on submerged plants Ceratophyllum demersum L

Nanoplastics (NPs) and heavy metals are widely distributed in aquatic ecosystem, posing a potential threat to ecosystem function. Submerged macrophytes play an important role in water purification and maintaining ecological functions. However, the coupled effects of NPs and cadmium (Cd) on submerged macrophytes physiology and the mechanisms involved are still unclear. Here, the potential effects of single and co-Cd/PSNPs exposure on Ceratophyllum demersum L. (C. demersum) were explored. Our results showed that NPs aggravated the inhibition of Cd on plant growth ate (a decrease of 35.54 %), reduced chlorophyll synthesis (a decrease of 15.84 %), and disrupted the antioxidant enzyme system (a decrease of 25.07 % on SOD activity) of C. demersum. Massive PSNPs adhered to the surface of C. demersum when exposed to co-Cd/PSNPs while they did not adhere when exposed to single-NPs. The metabolic analysis further demonstrated that co-exposure down-regulated plant cuticle synthesis and that Cd exacerbated the physical damage and shadowing effects of NPs. In addition, co-exposure upregulated pentose phosphate metabolism, leading to the accumulation of starch grains. Furthermore, PSNPs reduced Cd enrichment capacity of C. demersum. Our results unraveled distinct regulatory networks for submerged macrophytes exposed to single and composite of Cd and PSNPs, providing a new theoretical basis for assessing the risks of heavy metals and NPs in the freshwater environment.

Identification of hub genes and potential molecular mechanisms related to radiotherapy sensitivity in rectal cancer based on multiple datasets

BACKGROUND

Radiotherapy resistance is the main cause of low tumor regression for locally advanced rectum adenocarcinoma (READ). The biomarkers correlated to radiotherapy sensitivity and potential molecular mechanisms have not been completely elucidated.

METHODS

A mRNA expression profile and a gene expression dataset of READ (GSE35452) were acquired from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Differentially expressed genes (DEGs) between radiotherapy responder and non-responder of READ were screened out. Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis for DEGs were performed. Random survival forest analysis was used to identified hub genes by randomForestSRC package. Based on CIBERSORT algorithm, Genomics of Drug Sensitivity in Cancer (GDSC) database, Gene set variation analysis (GSVA), enrichment analysis (GSEA), nomogram, motif enrichment and non-coding RNA network analyses, the associations between hub genes and immune cell infiltration, drug sensitivity, specific signaling pathways, prognosis prediction and TF - miRNA regulatory and ceRNA network were investigated. The expressions of hub genes in clinical samples were displayed with the online Human Protein Atlas (HPA).

RESULTS

In total, 544 up-regulated and 575 down-regulated DEGs in READ were enrolled. Among that, three hubs including PLAGL2, ZNF337 and ALG10 were identified. These three hub genes were significantly associated with tumor immune infiltration, different immune-related genes and sensitivity of chemotherapeutic drugs. Also, they were correlated with the expression of various disease-related genes. In addition, GSVA and GSEA analysis revealed that different expression levels of PLAGL2, ZNF337 and ALG10 affected various signaling pathways related to disease progression. A nomogram and calibration curves based on three hub genes showed excellent prognosis predictive performance. And then, a regulatory network of transcription factor (ZBTB6) - mRNA (PLAGL2) and a ceRNA network of miRNA (has-miR-133b) - lncRNA were established. Finally, the results from HPA online database demonstrated the protein expression levels of PLAGL2, ZNF337 and ALG10 varied widely in READ patients.

CONCLUSION

These findings indicated that up-regulation of PLAGL2, ZNF337 and ALG10 in READ associated with radiotherapy response and involved in multiple process of cellular biology in tumor. They might be potential predictive biomarkers for radiotherapy sensitivity and prognosis for READ.

Mathematical modeling of regulatory networks of intracellular processes - Aims and selected methods

Regulatory networks structure and signaling pathways dynamics are uncovered in time- and resource consuming experimental work. However, it is increasingly supported by modeling, analytical and computational techniques as well as discrete mathematics and artificial intelligence applied to to extract knowledge from existing databases. This review is focused on mathematical modeling used to analyze dynamics and robustness of these networks. This paper presents a review of selected modeling methods that facilitate advances in molecular biology.

FoxO3 and oxidative stress: a multifaceted role in cellular adaptation

Oxidative stress is a major cause of morbidity and mortality in human health and disease. In this review, we focus on the Forkhead Box (Fox) subclass O3 (FoxO3), an extensively studied transcription factor that plays a pleiotropic role in a wide range of physiological and pathological processes by regulating multiple gene regulatory networks involved in the modulation of numerous aspects of cellular metabolism, including fuel metabolism, cell death, and stress resistance. This review will also focus on regulatory mechanisms of FoxO3 expression and activity, such as crucial post-translational modifications and non-coding RNAs. Moreover, this work discusses and evidences some pathways to how this transcription factor and reactive oxygen species regulate each other, which may lead to the pathogenesis of various types of diseases. Therefore, in addition to being a promising therapeutic target, the FoxO3-regulated signaling pathways can also be used as reliable diagnostic and prognostic biomarkers and indicators for drug responsiveness.

Halotolerance, stress mechanisms, and circadian clock of salt-tolerant cyanobacteria

Cyanobacteria harbor a high level of physiological flexibility, which enables them to reside in virtually all available environmental niches, including extreme environments. In this review, we summarize the recent advancements in stress mechanisms of salt-tolerant (a.k.a. halotolerant) cyanobacteria. Omics approaches have been extensively employed in recent years to decipher mechanisms of halotolerance and to understand the relevance of halotolerance-associated gene regulatory networks. The vast knowledge from genome mining disclosed that halotolerant cyanobacteria possess extended gene families and/or clusters, encoding enzymes that synthesize unique osmoprotectants, including glycine betaine (GB), betaine derivatives, and mycosporine-like amino acids (MAAs). Comprehensive transcriptomic analyses were conducted using Halothece sp. PCC7418 (hereafter referred to as Halothece), a cyanobacterium that exhibits remarkable halotolerance. These studies revealed a specific transcriptional response when Halothece was subjected to salt stress, whereas salt and osmotic stresses were found to share a common transcriptomic response. Transcriptome and metabolite analyses of Halothece illustrated a complex dynamic relationship between the biosyntheses of osmoprotectants, as well as corresponding and ancillary pathways. Lastly, novel insights highlight the relationship between the molecular regulation of the circadian rhythm and salt stress tolerance. Since the circadian rhythm of gene expression was distorted under salt stress, halotolerant cyanobacteria may prioritize the adaptation to salt stress by attenuation of circadian rhythmicity. KEY POINTS: • Recent advancements in the understanding of stress mechanisms in halotolerant cyanobacteria are described based on omics analyses. • Transcriptome and metabolite analyses of Halothece illustrated a complex dynamic relationship between the biosyntheses of osmoprotectants, as well as corresponding and ancillary pathways. • Since salt stress affects the molecular regulation among clock-related proteins, salt stress may attenuate circadian rhythmicity.

Integrated metabolome and transcriptome analysis unveils novel pathway involved in the fruit coloration of Nitraria tangutorum Bobr

BACKGROUND

The desert shrub Nitraria tangutorum Bobr. is important for its resistance to salt and alkali in Northwest China. It is an ecologically important species in this region and provides edible and medicinal berries. This study showed a mutant of N. tangutorum (named Jincan, JC) that has a strong yellow pericarp vs red in a wild type (represented by NT).

RESULTS

In this study, the secondary metabolic and molecular mechanisms responsible for Nitraria fruit coloration were investigated using LC-MS-based widely targeted metabolomics and transcriptomics data. As a result of our study, 122 and 104 flavonoid metabolites were differentially expressed throughout the mature and transition stages between JC and NT, respectively. Furthermore, two cyanidin derivatives (cyanidin 3-O-glucoside and cyanidin-3-O-(2''-O-glucosyl) glucoside) and one pelargonidin derivative (pelargonidin-3-O-glucoside) were identified only in the NT phenotype. The functional genes F3H (flavanone 3-hydroxylase), F3'H (flavonoid 3'-hydroxylase) and UFGT (flavonoid 3-O-glucosyltransferase) and the transcription factors MYB, bHLH, NAC and bZIP were significantly downregulated in JC. Meanwhile, the activity of UFGT was extremely low in both periods of JC, with a five-fold higher enzymatic activity of UFGT in RT than in YT. In summary, due to the lack of catalysis of UGFT, yellow fruit of JC could not accumulate sufficient cyanidin and pelargonidin derivatives during fruit ripening.

CONCLUSION

Taken together, our data provide insights into the mechanism for the regulation of anthocyanin synthesis and N. tangutorum fruit coloration and provide a theoretical basis to develop new strategies for developing bioactive compounds from N. tangutorum fruits.

Multi-omics inference of differential breast cancer-related transcriptional regulatory network gene hubs between young Black and White patients

OBJECTIVE

Breast cancers (BrCA) are a leading cause of illness and mortality worldwide. Black women have a higher incidence rate relative to white women prior to age 40 years, and a lower incidence rate after 50 years. The objective of this study is to identify -omics differences between the two breast cancer cohorts to better understand the disparities observed in patient outcomes.

MATERIALS AND METHODS

Using Standard SQL, we queried ISB-CGC hosted Google BigQuery tables storing TCGA BrCA gene expression, methylation, and somatic mutation data and analyzed the combined multi-omics results using a variety of methods.

RESULTS

Among Stage II patients 50 years or younger, genes PIK3CA and CDH1 are more frequently mutated in White (W50) than in Black or African American patients (BAA50), while HUWE1, HYDIN, and FBXW7 mutations are more frequent in BAA50. Over-representation analysis (ORA) and Gene Set Enrichment Analysis (GSEA) results indicate that, among others, the Reactome Signaling by ROBO Receptors gene set is enriched in BAA50. Using the Virtual Inference of Protein-activity by Enriched Regulon analysis (VIPER) algorithm, putative top 20 master regulators identified include NUPR1, NFKBIL1, ZBTB17, TEAD1, EP300, TRAF6, CACTIN, and MID2. CACTIN and MID2 are of prognostic value. We identified driver genes, such as OTUB1, with suppressed expression whose DNA methylation status were inversely correlated with gene expression. Networks capturing microRNA and gene expression correlations identified notable microRNA hubs, such as miR-93 and miR-92a-2, expressed at higher levels in BAA50 than in W50.

DISCUSSION/CONCLUSION

The results point to several driver genes as being involved in the observed differences between the cohorts. The findings here form the basis for further mechanistic exploration.

Protein-Protein Interaction Network Analysis Using NetworkX

In recent years, extracting information from biological data has become a particularly valuable way of gaining knowledge. Molecular interaction networks provide a framework for visualizing cellular processes, but their complexity frequently makes their interpretation difficult. Proteins are one of the primary determinants of biological function. Indeed, most biological activities in the living cells are functionally regulated by protein-protein interactions (PPIs). Thus, studying protein interactions is critical for understanding their roles within the cell. Exploring the PPI networks can open new avenues for future experimental studies and offer interspecies predictions for effective interaction mapping. In this chapter we will demonstrate how to construct, visualize, and analyze a protein-protein interaction network using NetworkX.

Developing and validating a survival prediction model based on blood exosomal ceRNA network in patients with PAAD

BACKGROUND

Among the most lethal cancers, pancreatic adenocarcinoma (PAAD) is an essential component of digestive system malignancies that still lacks effective diagnosis and treatment methods. As exosomes and competing endogenous RNA (ceRNA) regulatory networks in tumors go deeper, we expect to construct a ceRNA regulatory network derived from blood exosomes of PAAD patients by bioinformatics methods and develop a survival prediction model based on it.

METHODS

Blood exosome sequencing data of PAAD patients and normal controls were downloaded from the exoRbase database, and the expression profiles of exosomal mRNA, lncRNA, and circRNA were differentially analyzed by R. The related mRNA, circRNA, lncRNA, and their corresponding miRNA prediction data were imported into Cytoscape software to visualize the ceRNA network. Then, we conducted GO and KEGG enrichment analysis of mRNA in the ceRNA network. Genes that express differently in pancreatic cancer tissues compared with normal tissues and associate with survival (P < 0.05) were determined as Hub genes by GEPIA. We identified optimal prognosis-related differentially expressed mRNAs (DEmRNAs) and generated a risk score model by performing univariate and multivariate Cox regression analyses.

RESULTS

205 DEmRNAs, 118 differentially expressed lncRNAs (DElncRNAs), and 98 differentially expressed circRNAs (DEcircRNAs) were screened out. We constructed the ceRNA network, and a total of 26 mRNA nodes, 7 lncRNA nodes, 6 circRNA nodes, and 16 miRNA nodes were identified. KEGG enrichment analysis showed that the DEmRNAs in the regulatory network were mainly enriched in Human papillomavirus infection, PI3K-Akt signaling pathway, Osteoclast differentiation, and ECM-receptor interaction. Next, six hub genes (S100A14, KRT8, KRT19, MAL2, MYO5B, PSCA) were determined through GEPIA. They all showed significantly increased expression in cancer tissues compared with control groups, and their high expression pointed to adverse survival. Two optimal prognostic-related DEmRNAs, MYO5B (HR = 1.41, P < 0.05) and PSCA (HR = 1.10, P < 0.05) were included to construct the survival prediction model.

CONCLUSION

In this study, we successfully constructed a ceRNA regulatory network in blood exosomes from PAAD patients and developed a two-gene survival prediction model that provided new targets which shall aid in diagnosing and treating PAAD.

Engineering of global transcription factors in Bacillus, a genetic tool for increasing product yields: a bioprocess overview

Transcriptional factors are well studied in bacteria for their global interactions and the effects they produce at the phenotypic level. Particularly, Bacillus subtilis has been widely employed as a model Gram-positive microorganism used to characterize these network interactions. Bacillus species are currently used as efficient commercial microbial platforms to produce diverse metabolites such as extracellular enzymes, antibiotics, surfactants, industrial chemicals, heterologous proteins, among others. However, the pleiotropic effects caused by the genetic modification of specific genes that codify for global regulators (transcription factors) have not been implicated commonly from a bioprocess point of view. Recently, these strategies have attracted the attention in Bacillus species because they can have an application to increase production efficiency of certain commercial interest metabolites. In this review, we update the recent advances that involve this trend in the use of genetic engineering (mutations, deletion, or overexpression) performed to global regulators such as Spo0A, CcpA, CodY and AbrB, which can provide an advantage for the development or improvement of bioprocesses that involve Bacillus species as production platforms. Genetic networks, regulation pathways and their relationship to the development of growth stages are also discussed to correlate the interactions that occur between these regulators, which are important to consider for application in the improvement of commercial-interest metabolites. Reported yields from these products currently produced mostly under laboratory conditions and, in a lesser extent at bioreactor level, are also discussed to give valuable perspectives about their potential use and developmental level directed to process optimization at large-scale.

Discovering Common Pathogenic Mechanisms of COVID-19 and Parkinson Disease: An Integrated Bioinformatics Analysis

Coronavirus disease 2019 (COVID-19) has emerged since December 2019 and was later characterized as a pandemic by WHO, imposing a major public health threat globally. Our study aimed to identify common signatures from different biological levels to enlighten the current unclear association between COVID-19 and Parkinson's disease (PD) as a number of possible links, and hypotheses were reported in the literature. We have analyzed transcriptome data from peripheral blood mononuclear cells (PBMCs) of both COVID-19 and PD patients, resulting in a total of 81 common differentially expressed genes (DEGs). The functional enrichment analysis of common DEGs are mostly involved in the complement system, type II interferon gamma (IFNG) signaling pathway, oxidative damage, microglia pathogen phagocytosis pathway, and GABAergic synapse. The protein-protein interaction network (PPIN) construction was carried out followed by hub detection, revealing 10 hub genes (MX1, IFI27, C1QC, C1QA, IFI6, NFIX, C1S, XAF1, IFI35, and ELANE). Some of the hub genes were associated with molecular mechanisms such as Lewy bodies-induced inflammation, microglia activation, and cytokine storm. We investigated regulatory elements of hub genes at transcription factor and miRNA levels. The major transcription factors regulating hub genes are SOX2, XAF1, RUNX1, MITF, and SPI1. We propose that these events may have important roles in the onset or progression of PD. To sum up, our analysis describes possible mechanisms linking COVID-19 and PD, elucidating some unknown clues in between.

Using empirical biological knowledge to infer regulatory networks from multi-omics data

Background

Integration of multi-omics data can provide a more complex view of the biological system consisting of different interconnected molecular components, the crucial aspect for developing novel personalised therapeutic strategies for complex diseases. Various tools have been developed to integrate multi-omics data. However, an efficient multi-omics framework for regulatory network inference at the genome level that incorporates prior knowledge is still to emerge.

Results

We present IntOMICS, an efficient integrative framework based on Bayesian networks. IntOMICS systematically analyses gene expression, DNA methylation, copy number variation and biological prior knowledge to infer regulatory networks. IntOMICS complements the missing biological prior knowledge by so-called empirical biological knowledge, estimated from the available experimental data. Regulatory networks derived from IntOMICS provide deeper insights into the complex flow of genetic information on top of the increasing accuracy trend compared to a published algorithm designed exclusively for gene expression data. The ability to capture relevant crosstalks between multi-omics modalities is verified using known associations in microsatellite stable/instable colon cancer samples. Additionally, IntOMICS performance is compared with two algorithms for multi-omics regulatory network inference that can also incorporate prior knowledge in the inference framework. IntOMICS is also applied to detect potential predictive biomarkers in microsatellite stable stage III colon cancer samples.

Conclusions

We provide IntOMICS, a framework for multi-omics data integration using a novel approach to biological knowledge discovery. IntOMICS is a powerful resource for exploratory systems biology and can provide valuable insights into the complex mechanisms of biological processes that have a vital role in personalised medicine.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04891-9.

Identification of miRNA-mRNA-TF regulatory networks in peripheral blood mononuclear cells of type 1 diabetes

BACKGROUND

Type 1 diabetes (T1D) is a T lymphocyte-mediated and B lymphocyte-assisted autoimmune disease. We aimed to identify abnormally expressed genes in peripheral blood mononuclear cells (PBMCs) of T1D and explore their possible molecular regulatory network.

METHODS

Expression datasets were downloaded from the Gene Expression Omnibus (GEO) database. Then, the differentially expressed genes (DEGs) and differentially expressed miRNAs (DEmiRNAs) were identified, and functional enrichment and immune cell infiltration analysis were performed. The starBase, miRTarBase, TarBase, JASPAR, ENCODE, and TRRUST databases constructed the miRNA-mRNA-TF regulatory network. The ROC curves were plotted to evaluate the sensitivity and specificity of miRNAs and mRNAs.

RESULT

A total of 216 DEGs directly or indirectly related to type I diabetes mellitus, natural killer cell-mediated cytotoxicity, Th1, and Th2 cell differentiation, and the IL-17 and TNF signaling pathways were obtained. The miRNA-mRNA-TF network indicates that miR-320a and SOX5 are the only miRNAs and TFs that both target ADM and RRAGD. The ROC curves showed that ADM (0.9375), RRAGD (0.8958), and hsa-mir-320a (0.9417) had high accuracy in T1D diagnosis.

CONCLUSION

The constructed regulatory networks, including miR-320a/ADM/SOX5 and miR-320a/RRAGD/SOX5, may provide new insight into the mechanisms of development and progression in T1D.

Phytohormone profiles and related gene expressions after endodormancy release in developing Pinus tabuliformis male strobili

Development after endo-dormancy release ensures perennial plants, such as forest trees, proper response to environmental changes and enhances their adaptability. In northern hemisphere, megasporophore and microsporophore of conifers undergo dormancy to complete their development. Here combined with transcriptome data, we used high-performance liquid chromatography/electrospray ionization tandem mass spectrometry (ESI-HPLC-MS/MS) to quantitatively analyse the various hormones (Abscisic Acid (ABA), 3-Indoleacetic acid (IAA), Gibberellins (GAs), Cytokinin (CTK), Jasmonic acid (JA) and Salicylic acid (SA)) of Chinese pine (Pinus tabuliformis Carr.) male strobili after endo-dormancy release. More specifically, we analysed endogenous hormones and their related-genes and verified the important role of ABA in plants growth and development. We observed rapid decrease in ABA content after dormancy release, resulting in reducing the inhibitory effect on male strobili growth. Similarly, rapid drop in ABA/GA ratio was observed and was associated with the start of male strobili growth and development. Combined with transcriptome data, we found that HAB2-SnRK2.10 played a central role in the ABA pathway in the entire network of hormones regulating male strobili development. Due to external environment warming, the differentially expressed HAB2-SnRK gene led to ABA content rapid decline, thus initiating male strobili growth. We constructed a network of hormone-regulated development to understand the interactions between hormones after male strobili dormancy release of male strobili. This study provided essential foundations for studying megasporophore and microsporophore growth mechanism after endo-dormancy and offered new ideas for flower development in gymnosperms and angiosperms.

Structure and behavior in Boolean monotonic model pools

One of the main theoretical questions in the field of discrete regulatory networks is the question what aspects of the dynamics - the structure of the state transition graph - are already imposed by structural descriptions of the network such as the interaction graph. For Boolean networks, prior work has concentrated on different versions of the Thomas conjectures that link feedback cycles in the network structure to attractor properties. Other approaches check algorithmically whether certain properties hold true for all models sharing specific structural constraints, e.g. by using model checking techniques. In this work we investigate the behavior of the pool of Boolean networks in agreement with a given interaction graph using a different approach. Grouping together states that are updated consistently across the pool we derive an equivalence relation and analyze a corresponding quotient graph on the state space. By construction this graph yields information about the dynamics of all functions in the pool. Our main result is that this graph can be computed efficiently without enumerating and analyzing all individual functions. This opens up new possibilities for applications, where such model pools arise when modeling under uncertainty.

lycopersicum

BACKGROUND

Transposable elements (TE) make up a large portion of many plant genomes and are playing innovative roles in genome evolution. Several TEs can contribute to gene regulation by influencing expression of nearby genes as stress-responsive regulatory motifs. To delineate TE-mediated plant stress regulatory networks, we took a 2-step computational approach consisting of identifying TEs in the proximity of stress-responsive genes, followed by searching for cis-regulatory motifs in these TE sequences and linking them to known regulatory factors. Through a systematic meta-analysis of RNA-seq expression profiles and genome annotations, we investigated the relation between the presence of TE superfamilies upstream, downstream or within introns of nearby genes and the differential expression of these genes in various stress conditions in the TE-poor Arabidopsis thaliana and the TE-rich Solanum lycopersicum.

RESULTS

We found that stress conditions frequently expressed genes having members of various TE superfamilies in their genomic proximity, such as SINE upon proteotoxic stress and Copia and Gypsy upon heat stress in A. thaliana, and EPRV and hAT upon infection, and Harbinger, LINE and Retrotransposon upon light stress in S. lycopersicum. These stress-specific gene-proximal TEs were mostly located within introns and more detected near upregulated than downregulated genes. Similar stress conditions were often related to the same TE superfamily. Additionally, we detected both novel and known motifs in the sequences of those TEs pointing to regulatory cooption of these TEs upon stress. Next, we constructed the regulatory network of TFs that act through binding these TEs to their target genes upon stress and discovered TE-mediated regulons targeted by TFs such as BRB/BPC, HD, HSF, GATA, NAC, DREB/CBF and MYB factors in Arabidopsis and AP2/ERF/B3, NAC, NF-Y, MYB, CXC and HD factors in tomato.

CONCLUSIONS

Overall, we map TE-mediated plant stress regulatory networks using numerous stress expression profile studies for two contrasting plant species to study the regulatory role TEs play in the response to stress. As TE-mediated gene regulation allows plants to adapt more rapidly to new environmental conditions, this study contributes to the future development of climate-resilient plants.

Enabling Studies of Genome-Scale Regulatory Network Evolution in Large Phylogenies with MRTLE

Transcriptional regulatory networks specify context-specific patterns of genes and play a central role in how species evolve and adapt. Inferring genome-scale regulatory networks in non-model species is the first step for examining patterns of conservation and divergence of regulatory networks. Transcriptomic data obtained under varying environmental stimuli in multiple species are becoming increasingly available, which can be used to infer regulatory networks. However, inference and analysis of multiple gene regulatory networks in a phylogenetic setting remains challenging. We developed an algorithm, Multi-species Regulatory neTwork LEarning (MRTLE), to facilitate such studies of regulatory network evolution. MRTLE is a probabilistic graphical model-based algorithm that uses phylogenetic structure, transcriptomic data for multiple species, and sequence-specific motifs in each species to simultaneously infer genome-scale regulatory networks across multiple species. We applied MRTLE to study regulatory network evolution across six ascomycete yeasts using transcriptomic measurements collected across different stress conditions. MRTLE networks recapitulated experimentally derived interactions in the model organism S. cerevisiae as well as non-model species, and it was more beneficial for network inference than methods that do not use phylogenetic information. We examined the regulatory networks across species and found that regulators associated with significant expression and network changes are involved in stress-related processes. MTRLE and its associated downstream analysis provide a scalable and principled framework to examine evolutionary dynamics of transcriptional regulatory networks across multiple species in a large phylogeny.

Unraveling molecular mechanism underlying biomaterial and stem cells interaction during cell fate commitment using high throughput data analysis

Stem cell differentiation towards various somatic cells and body organs has proven to be an effective technique in the understanding and progression of regenerative medicine. Despite the advances made, concerns regarding the low efficiency of differentiation and the remaining differences between stem cell products and their in vivo counterparts must be addressed. Biomaterials that mimic endogenous growth conditions represent one recent method used to improve the quality and efficiency of stem cell differentiation, though the mechanisms of this improvement remain to be completely understood. The effectiveness of various biomaterials can be analyzed through a multidisciplinary approach involving bioinformatics and systems biology tools. Here, we aim to use bioinformatics to accomplish two aims: 1) determine the effect of different biomaterials on stem cell growth and differentiation, and 2) understand the effect of cell of origin on the differentiation potential of multipotent stem cells. First, we demonstrate that the dimensionality (2D versus 3D) and the degradability of biomaterials affects the way that the cells are able to grow and differentiate at the transcriptional level. Additionally, according to transcriptional state of the cells, the particular cell of origin is an important factor in determining the response of stem cells to same biomaterial. Our data demonstrates the ability of bioinformatics to understand novel molecular mechanisms and context by which stem cells are most efficiently able to differentiate. These results and strategies can be used to suggest proper combinations of biomaterials and stem cells to achieve high differentiation efficiency and functionality of desired cell types.

Bifurcation Analysis of a Heat-Sensitive Epigenetic Regulatory Network

An epigenetic regulatory network that influences transgenerational inheritance of a heat-altered phenotype was recently discovered in Arabidopsis. Our analysis shows that transgenerational inheritance of the heat-altered phenotype operates in a switch-like manner and can be turned on or off as a function of heat. We also show that trans-acting small interfering RNAs act as an "inverse amplifier" of HTT5, the protein that controls the heat-altered phenotype by a currently unknown mechanism. Our analysis uses the resultant to find an analytic expression for a cusp curve in parameter space and to find a parameter bound on switch-like behavior.

An emergent Wnt5a/YAP/TAZ regulatory circuit and its possible role in cancer

Wnt5a is a ligand that plays several roles in development, homeostasis, and disease. A growing body of evidence indicates that Wnt5a is involved in cancer progression. Despite extensive research in this field, our knowledge about how Wnt5a is precisely involved in cancer is still incomplete. It is usually thought that certain combinations of Frizzled receptors and co-receptors might explain the observed effects of Wnt5a either as a tumor suppressor or by promoting migration and invasion. While accepting this 'receptor context' model, this review proposes that Wnt5a is integrated within a larger regulatory circuit involving β-catenin, YAP/TAZ, and LATS1/2. Remarkably, WNT5A and YAP1 are transcriptionally regulated by the Hippo and Wnt pathways, respectively, and might form a regulatory circuit acting through LATS kinases and secreted Wnt/β-catenin inhibitors, including Wnt5a itself. Therefore, understanding the precise role of Wnt5a and YAP in cancer requires a systems biology perspective.

Functional and regulatory aspects of oxidative stress response in X monosomy

The partial/complete loss of one X chromosome in a human female leads to Turner syndrome (TS). TS individuals display a range of phenotypes including short stature, osteoporosis, ovarian malfunction, diabetes, and thyroid dysfunction. Epigenetic factors and regulatory networks are distinctly different in X monosomy (45, X). In a lifetime, an individual is exposed to a variety of stress conditions. To study whether X monosomy cells display a differential response upon exposure to mild stress as compared to normal 46, XX cells and whether this may contribute to various co-morbidities in aneuploid individuals, we have carried out a transcriptomic analysis of human fibroblasts 45, X and 46, XX after exposure to mild oxidative stress. Under these conditions, over 350 transcripts were seen to be differentially expressed in 45, X and 46, XX cells. Pathways associated with oxidative stress were differentially expressed highlighting the differential regulation of genes and associated phenotypes. It could be seen that X monosomy cells are more susceptible to oxidative stress as compared to normal cells and have altered molecular pathways both in normal conditions and also upon exposure to mild oxidative stress. To explore this aspect in detail, we have mapped the expressions of transcription factors (TFs) in 45, X and 46, XX cells. The network of transcription activating factors is differentially regulated in 45, X and 46, XX cells under stress exposure. It is tempting to speculate that the altered ability of 45, X (Turner) cells to respond to stress may play a significant role in the physiological function and altered phenotypes in Turner syndrome.

Dynamic Host Immune and Transcriptomic Responses to Respiratory Syncytial Virus Infection in a Vaccination-Challenge Mouse Model

Respiratory syncytial virus (RSV) is the major cause of lower respiratory tract infections in children. Inactivated RSV vaccine was developed in the late 1960's, but the vaccine-enhanced disease (VED) occurred to vaccinated infants upon subsequent natural RSV infection. The excessive inflammatory immunopathology in the lungs might be involved in the VED, but the underlying mechanisms remain not fully understood. In this study, we utilized UV-inactivated RSV in the prime/boost approach followed by RSV challenge in BALB/c mice to mimic RSV VED. The dynamic virus load, cytokines, histology and transcriptome profiles in lung tissues of mice were investigated from day 1 to day 6 post-infection. Compared to PBS-treated mice, UV-RSV vaccination leads to a Th2 type inflammatory response characterized by enhanced histopathology, reduced Treg cells and increased IL4⁺CD4 T cells in the lung. Enhanced production of several Th2 type cytokines (IL-4, IL-5, IL-10) and TGF-β,  reduction of IL-6 and IL-17 were observed in UV-RSV vaccinated mice. A total of 5582 differentially expressed (DE) genes between PBS-treated or vaccinated mice and naïve mice were identified by RNA-Seq. Eleven conserved high-influential modules (HMs) were recognized, majorly grouped into regulatory networks related to cell cycle and cell metabolism, signal transduction, immune and inflammatory responses. At an early time post-infection, the vaccinated mice showed obvious decreased expression patterns of DE genes in 11 HMs compared to PBS-treated mice. The extracellular matrix (HM5) and immune responses (HM8) revealed tremendous differences in expression and regulation characteristics of transcripts between PBS-treated and vaccinated mice at both early and late time points. The highly connected genes in HM5 and HM8 networks were further validated by RT-qPCR. These findings reveal the relationship between RSV VED and immune responses, which could benefit the development of novel RSV vaccines.

Very long intergenic non-coding (vlinc) RNAs directly regulate multiple genes in cis and trans

BACKGROUND

The majority of the human genome is transcribed in the form of long non-coding (lnc) RNAs. While these transcripts have attracted considerable interest, their molecular mechanisms of function and biological significance remain controversial. One of the main reasons behind this lies in the significant challenges posed by lncRNAs requiring the development of novel methods and concepts to unravel their functionality. Existing methods often lack cross-validation and independent confirmation by different methodologies and therefore leave significant ambiguity as to the authenticity of the outcomes. Nonetheless, despite all the caveats, it appears that lncRNAs may function, at least in part, by regulating other genes via chromatin interactions. Therefore, the function of a lncRNA could be inferred from the function of genes it regulates. In this work, we present a genome-wide functional annotation strategy for lncRNAs based on identification of their regulatory networks via the integration of three distinct types of approaches: co-expression analysis, mapping of lncRNA-chromatin interactions, and assaying molecular effects of lncRNA knockdowns obtained using an inducible and highly specific CRISPR/Cas13 system.

RESULTS

We applied the strategy to annotate 407 very long intergenic non-coding (vlinc) RNAs belonging to a novel widespread subclass of lncRNAs. We show that vlincRNAs indeed appear to regulate multiple genes encoding proteins predominantly involved in RNA- and development-related functions, cell cycle, and cellular adhesion via a mechanism involving proximity between vlincRNAs and their targets in the nucleus. A typical vlincRNAs can be both a positive and negative regulator and regulate multiple genes both in trans and cis. Finally, we show vlincRNAs and their regulatory networks potentially represent novel components of DNA damage response and are functionally important for the ability of cancer cells to survive genotoxic stress.

CONCLUSIONS

This study provides strong evidence for the regulatory role of the vlincRNA class of lncRNAs and a potentially important role played by these transcripts in the hidden layer of RNA-based regulation in complex biological systems.

FGF signalling plays similar roles in development and regeneration of the skeleton in the brittle star Amphiura filiformis

Regeneration as an adult developmental process is in many aspects similar to embryonic development. Although many studies point out similarities and differences, no large-scale, direct and functional comparative analyses between development and regeneration of a specific cell type or structure in one animal exist. Here, we use the brittle star Amphiura filiformis to characterise the role of the FGF signalling pathway during skeletal development in embryos and arm regeneration. In both processes, we find ligands expressed in ectodermal cells that flank underlying skeletal mesenchymal cells, which express the receptors. Perturbation of FGF signalling showed inhibited skeleton formation in both embryogenesis and regeneration, without affecting other key developmental processes. Differential transcriptome analysis finds mostly differentiation genes rather than transcription factors to be downregulated in both contexts. Moreover, comparative gene analysis allowed us to discover brittle star-specific differentiation genes. In conclusion, our results show that the FGF pathway is crucial for skeletogenesis in the brittle star, as in other deuterostomes, and provide evidence for the re-deployment of a developmental gene regulatory module during regeneration.

The impaired biosynthetic networks in defective tapetum lead to male sterility in watermelon

Heterosis has been widely applied in watermelon breeding, because of the higher resistance and yield of hybrid. As the basis of heterosis utilization, genic male sterility (GMS) is an important tool for facilitating hybrid seed production, while the detailed mechanism in watermelon is still largely unknown. Here, we report a spontaneous mutant Se18 exhibited complete male sterility due to the uniquely multilayered tapetum and the un-meiotic pollen mother cells during pollen development. Using TMT based quantitative proteomic analyses, a total of 348 differentially abundant proteins (DAPs) were detected with the overwhelming majority down-regulated in mutant Se18. By analyzing the putative orthologs/homologs of Arabidopsis GMS related genes, the biosynthesis and transport of sporopollenin and tryphine precursors were predictably altered in mutant compared to its sibling wild type. Moreover, the general phenylpropanoid pathway as well as its related metabolisms was also expectably impaired in mutant, coincident with the pale yellow petals. Notably, some key transcriptional factors regulating tapetum development, together with their down-regulated targets, offered potentially valuable candidates regarding of male sterility. Collectively, the disrupted regulatory networks underlying male sterility of watermelon was proposed, which provide novel insights into genetic mechanism of male reproductive process and rich gene resources for future research. SIGNIFICANCE: Watermelon is an importantly economical cucurbit crop worldwide, with high nutritional value. Although several male sterile mutants have been identified in watermelon, the underlying molecular mechanism is poorly elucidated. Comparative cytological analysis revealed that the defective development of tapetum was responsible for male sterility in mutant Se18. Combined with the morphological comparison, male floral buds at 2.0-2.5 mm in diameter were confirmed with no obvious phenotypic differences but distinct cytological defects, which were in turn sampled for TMT based proteomic analyses. Referring to functionally characterized GMS related genes, the genetic pathway DYT1-TDF1-AMS-MS188-MS1 regulating tapetum development, together with some downstream targets, were considerably altered in mutant Se18. Moreover, enrichment analyses illustrated the general phenylpropanoid related metabolisms, as well as the biosynthesis and transport of sporopollenin and tryphine precursors, were significantly disrupted in defective anther development. Collectively, the proposed regulatory networks in watermelon not only contribute to a better understanding of molecular mechanisms underlying male sterility, but also provide valuable GMS related candidates for future researches.

Constitutive expression of Ribosomal Protein L6 modulates salt tolerance in rice transgenic plants

We have functionally characterized the RPL6, a Ribosomal Protein Large subunit gene for salt stress tolerance in rice. The overexpression of RPL6 resulted in tolerance to moderate (150 mM) to high (200 mM) levels of salt (NaCl). The transgenic rice plants expressing RPL6 constitutively showed better phenotypic and physiological responses with high quantum efficiency, accumulation of higher chlorophyll and proline contents, and an overall increase in seed yield compared with the wild type in salt stress treatments. An iTRAQ-based comparative proteomic analysis revealed the high expression of about 333 proteins among the 4378 DAPs in a selected overexpression line of RPL6 treated with 200 mM of NaCl. The functional analysis showed that these highly accumulated proteins (HAPs) are involved in photosynthesis, ribosome and chloroplast biogenesis, ion transportation, transcription and translation regulation, phytohormone and secondary metabolite signal transduction. An in silico network analysis of HAPs predicted that RPL6 binds with translation-related proteins and helicases, which coordinately affect the activities of a comprehensive signaling network, thereby inducing tolerance and promoting growth and productivity in response to salt stress. Our overall findings identified a novel candidate, RPL6, whose characterization contributed to the existing knowledge on the complexity of salt tolerance mechanism in plants.

A systemic approach to explore the mechanisms of drug resistance and altered signaling cascades in extensively drug-resistant tuberculosis

BACKGROUND AND AIM

The persistence of extensively drug-resistant (XDR) strains of Mycobacterium tuberculosis (MTB) continue to pose a significant challenge to the treatment and control of tuberculosis infections worldwide. XDR-MTB strains exhibit resistance against first-line anti-TB drugs, fluoroquinolones, and second-line injectable drugs. The mechanisms of drug resistance of MTB remains poorly understood. Our study aims at identifying the differentially expressed genes (DEGs), associated gene networks, and signaling cascades involved in rendering this pathogen resistant to multiple drugs, namely, isoniazid, rifampicin, and capreomycin.

METHODS

We used the microarray dataset GSE53843. The GEO2R tool was used to prioritize the most significant DEGs (top 250) of each drug exposure sample between XDR strains and non-resistant strains. The validation of the 250 DEGs was performed using volcano plots. Protein-protein interaction networks of the DEGs were created using STRING and Cytoscape tools, which helped decipher the relationship between these genes. The significant DEGs were functionally annotated using DAVID and ClueGO. The concomitant biological processes (BP) and molecular functions (MF) were represented as dot plots.

RESULTS AND CONCLUSION

We identified relevant molecular pathways and biological processes, such as cell wall biogenesis, lipid metabolic process, ion transport, phosphopantetheine binding, and triglyceride lipase activity. These processes indicated the involvement of multiple interconnected mechanisms in drug resistance. Our study highlighted the impact of cell wall permeability, with the dysregulation of the mur family of proteins, as essential factors in the inference of resistance. Additionally, upregulation of genes responsible for ion transport such as ctpF, arsC, and nark3, emphasizes the importance of transport channels and efflux pumps in potentially driving out stress-inducing compounds. This study investigated the upregulation of the Lip family of proteins, which play a crucial role in triglyceride lipase activity. Thereby illuminating the potential role of drug-induced dormancy and subsequent resistance in the mycobacterial strains. Multiple mechanisms such as carboxylic acid metabolic process, NAD biosynthetic process, triglyceride lipase activity, phosphopantetheine binding, organic acid biosynthetic process, and growth of symbiont in host cell were observed to partake in resistance of XDR-MTB. This study ultimately provides a platform for important mapping targets for potential therapeutics against XDR-MTB.

circFOXO3: Going around the mechanistic networks in cancer by interfering with miRNAs regulatory networks

Circular RNAs (circRNA) have gained recent interest due to their functional versatility due to their interactions with other RNA species and proteins, all of which underline complex regulatory networks involved in pathogenic mechanisms. As a result, recent insights in circRNA biology are investigating their biomarker and therapeutic potential. One such circRNA is CircFOXO3, which consists of the circularized second exon of the FOXO3 mRNA, a member of the forkhead box transcription factor family involved in the regulation of developmental programs. Recent research focused on the role of circFOXO3 in the context of cancer has highlighted several implications in key tumorigenesis mechanisms, thus consolidating its relevance among other identified circRNAs. In this paper, we will focus on the currently identified case-specific implications of circFOXO3 in cancer, with a focus on the circFOXO3-miRNA-mRNA regulatory networks, its interactions with different proteins, and their cumulated biological effects upon tumor development. Therefore, we aim to provide an integrated perspective of the mechanistic implications of circFOXO3 in different cancers while also highlighting its biomarker or therapeutic potential based on the current evidence.

Cyclin-dependent kinases and rare developmental disorders

Extensive studies in the past 30 years have established that cyclin-dependent kinases (CDKs) exert many diverse, important functions in a number of molecular and cellular processes that are at play during development. Not surprisingly, mutations affecting CDKs or their activating cyclin subunits have been involved in a variety of rare human developmental disorders. These recent findings are reviewed herein, giving a particular attention to the discovered mutations and their demonstrated or hypothesized functional consequences, which can account for pathological human phenotypes. The review highlights novel, important CDK or cyclin functions that were unveiled by their association with human disorders, and it discusses the shortcomings of mouse models to reveal some of these functions. It explains how human genetics can be used in combination with proteome-scale interaction databases to loom regulatory networks around CDKs and cyclins. Finally, it advocates the use of these networks to profile pathogenic CDK or cyclin variants, in order to gain knowledge on protein function and on pathogenic mechanisms.

Genome-wide identification and characterization of cucumber bHLH family genes and the functional characterization of CsbHLH041 in NaCl and ABA tolerance in Arabidopsis and cucumber

BACKGROUND

The basic/helix-loop-helix (bHLH) transcription factor family exists in all three eukaryotic kingdoms as important participants in biological growth and development. To date, the comprehensive genomic and functional analyses of bHLH genes has not been reported in cucumber (Cucumis sativus L.).

RESULTS

Here, a total of 142 bHLH genes were identified and classified into 32 subfamilies according to the conserved motifs, phylogenetic analysis and gene structures in cucumber. The sequences of CsbHLH proteins were highly conserved based on the results of multiple sequence alignment analyses. The chromosomal distribution, synteny analysis, and gene duplications of these 142 CsbHLHs were further analysed. Many elements related to stress responsiveness and plant hormones were present in the promoter regions of CsbHLH genes based on a cis-element analysis. By comparing the phylogeny of cucumber and Arabidopsis bHLH proteins, we found that cucumber bHLH proteins were clustered into different functional clades of Arabidopsis bHLH proteins. The expression analysis of selected CsbHLHs under abiotic stresses (NaCl, ABA and low-temperature treatments) identified five CsbHLH genes that could simultaneously respond to the three abiotic stresses. Tissue-specific expression profiles of these five genes were also analysed. In addition, 35S:CsbHLH041 enhanced the tolerance to salt and ABA in transgenic Arabidopsis and in cucumber seedlings, suggesting CsbHLH041 is an important regulator in response to abiotic stresses. Lastly, the functional interoperability network among the CsbHLH proteins was analysed.

CONCLUSION

This study provided a good foundation for further research into the functions and regulatory mechanisms of CsbHLH proteins and identified candidate genes for stress resistance in cucumber.

The Host Cell ViroCheckpoint: Identification and Pharmacologic Targeting of Novel Mechanistic Determinants of Coronavirus-Mediated Hijacked Cell States

Most antiviral agents are designed to target virus-specific proteins and mechanisms rather than the host cell proteins that are critically dysregulated following virus-mediated reprogramming of the host cell transcriptional state. To overcome these limitations, we propose that elucidation and pharmacologic targeting of host cell Master Regulator proteins-whose aberrant activities govern the reprogramed state of coronavirus-infected cells-presents unique opportunities to develop novel mechanism-based therapeutic approaches to antiviral therapy, either as monotherapy or as a complement to established treatments. Specifically, we propose that a small module of host cell Master Regulator proteins (ViroCheckpoint) is hijacked by the virus to support its efficient replication and release. Conventional methodologies are not well suited to elucidate these potentially targetable proteins. By using the VIPER network-based algorithm, we successfully interrogated 12h, 24h, and 48h signatures from Calu-3 lung adenocarcinoma cells infected with SARS-CoV, to elucidate the time-dependent reprogramming of host cells and associated Master Regulator proteins. We used the NYS CLIA-certified Darwin OncoTreat algorithm, with an existing database of RNASeq profiles following cell perturbation with 133 FDA-approved and 195 late-stage experimental compounds, to identify drugs capable of virtually abrogating the virus-induced Master Regulator signature. This approach to drug prioritization and repurposing can be trivially extended to other viral pathogens, including SARS-CoV-2, as soon as the relevant infection signature becomes available.

Abasy Atlas v2.2: The most comprehensive and up-to-date inventory of meta-curated, historical, bacterial regulatory networks, their completeness and system-level characterization

Some organism-specific databases about regulation in bacteria have become larger, accelerated by high-throughput methodologies, while others are no longer updated or accessible. Each database homogenize its datasets, giving rise to heterogeneity across databases. Such heterogeneity mainly encompasses different names for a gene and different network representations, generating duplicated interactions that could bias network analyses. Abasy (Across-bacteria systems) Atlas consolidates information from different sources into meta-curated regulatory networks in bacteria. The high-quality networks in Abasy Atlas enable cross-organisms analyses, such as benchmarking studies where gold standards are required. Nevertheless, network incompleteness still casts doubts on the conclusions of network analyses, and available sampling methods cannot reflect the curation process. To tackle this problem, the updated version of Abasy Atlas presented in this work provides historical snapshots of regulatory networks. Thus, network analyses can be performed at different completeness levels, making possible to identify potential bias and to predict future results. We leverage the recently found constraint in the complexity of regulatory networks to develop a novel model to quantify the total number of regulatory interactions as a function of the genome size. This completeness estimation is a valuable insight that may aid in the daunting task of network curation, prediction, and validation. The new version of Abasy Atlas provides 76 networks (204,282 regulatory interactions) covering 42 bacteria (64% Gram-positive and 36% Gram-negative) distributed in 9 species (Mycobacterium tuberculosis, Bacillus subtilis, Escherichia coli, Corynebacterium glutamicum, Staphylococcus aureus, Pseudomonas aeruginosa, Streptococcus pyogenes, Streptococcus pneumoniae, and Streptomyces coelicolor), containing 8459 regulons and 4335 modules.

Database URL: https://abasy.ccg.unam.mx/.

Developing a virus-microRNA interactome using cytoscape

It is currently difficult to determine the effect of oncogenic viruses on the global function and regulation of pathways within mammalian cells. A thorough understanding of the molecular pathways and individual genes altered by oncogenic viruses is needed for the identification of targets that can be utilised for early diagnosis, prevention, and treatment methods. We detail a logical step-by-step guide to uncover viral-protein-miRNA interactions using publically available datasets and the network building program, Cytoscape. This method may be applied to identify specific pathways that are altered in viral infection, and contribute to the oncogenic transformation of cells. To demonstrate this, we constructed a gene regulatory interactome encompassing Human Papillomavirus Type 16 (HPV16) and its control of specific miRNAs. This approach can be broadly applied to understand and map the regulatory functions of other oncogenic viruses, and determine their role in altering the cellular environment in cancer.

Availability and Implementation Cytoscape (Shannon et al. (2003), Smoot et al. (2010)) is freely available at https://cytoscape.org/.

•

This method allows for the analysis and visualization of large datasets to generate an interactome that integrates key players of molecular biology

•

This approach may be applied to any oncogenic virus to map its regulatory functions, and its secondary impact on gene regulation via microRNAs.

Transcriptome sequencing of Saccharina japonica sporophytes during whole developmental periods reveals regulatory networks underlying alginate and mannitol biosynthesis

BACKGROUND

Alginate is an important cell wall component and mannitol is a soluble storage carbon substance in the brown seaweed Saccharina japonica. Their contents vary with kelp developmental periods and harvesting time. Alginate and mannitol regulatory networks and molecular mechanisms are largely unknown.

RESULTS

With WGCNA and trend analysis of 20,940 known genes and 4264 new genes produced from transcriptome sequencing of 30 kelp samples from different stages and tissues, we deduced that ribosomal proteins, light harvesting complex proteins and "imm upregulated 3" gene family are closely associated with the meristematic growth and kelp maturity. Moreover, 134 and 6 genes directly involved in the alginate and mannitol metabolism were identified, respectively. Mannose-6-phosphate isomerase (MPI2), phosphomannomutase (PMM1), GDP-mannose 6-dehydrogenase (GMD3) and mannuronate C5-epimerase (MC5E70 and MC5E122) are closely related with the high content of alginate in the distal blade. Mannitol accumulation in the basal blade might be ascribed to high expression of mannitol-1-phosphate dehydrogenase (M1PDH1) and mannitol-1-phosphatase (M1Pase) (in biosynthesis direction) and low expression of mannitol-2-dehydrogenase (M2DH) and Fructokinase (FK) (in degradation direction). Oxidative phosphorylation and photosynthesis provide ATP and NADH for mannitol metabolism whereas glycosylated cycle and tricarboxylic acid (TCA) cycle produce GTP for alginate biosynthesis. RNA/protein synthesis and transportation might affect alginate complex polymerization and secretion processes. Cryptochrome (CRY-DASH), xanthophyll cycle, photosynthesis and carbon fixation influence the production of intermediate metabolite of fructose-6-phosphate, contributing to high content of mannitol in the basal blade.

CONCLUSIONS

The network of co-responsive DNA synthesis, repair and proteolysis are presumed to be involved in alginate polymerization and secretion, while upstream light-responsive reactions are important for mannitol accumulation in meristem of kelp. Our transcriptome analysis provides new insights into the transcriptional regulatory networks underlying the biosynthesis of alginate and mannitol during S. japonica developments.

Quantitative analysis reveals reciprocal regulations underlying recovery dynamics of thymocytes and thymic environment in mice

Thymic crosstalk, a set of reciprocal regulations between thymocytes and the thymic environment, is relevant for orchestrating appropriate thymocyte development as well as thymic recovery from various exogenous insults. In this work, interactions shaping thymic crosstalk and the resultant dynamics of thymocytes and thymic epithelial cells are inferred based on quantitative analysis and modeling of the recovery dynamics induced by irradiation. The analysis identifies regulatory interactions consistent with known molecular evidence and reveals their dynamic roles in the recovery process. Moreover, the analysis also predicts, and a subsequent experiment verifies, a previously unrecognized regulation of CD4+CD8+ double positive thymocytes which temporarily increases their proliferation rate upon the decrease in their population size. Our model establishes a pivotal step towards the dynamic understanding of thymic crosstalk as a regulatory network system.