Transcriptional Co-expression Research Articles

Mapping of flumioxazin tolerance in a snap bean diversity panel leads to the discovery of a master genomic region controlling multiple stress resistance genes.

Effective weed management tools are crucial for maintaining the profitable production of snap bean (Phaseolus vulgaris L.). Preemergence herbicides help the crop to gain a size advantage over the weeds, but the few preemergence herbicides registered in snap bean have poor waterhemp (Amaranthus tuberculatus) control, a major pest in snap bean production. Waterhemp and other difficult-to-control weeds can be managed by flumioxazin, an herbicide that inhibits protoporphyrinogen oxidase (PPO). However, there is limited knowledge about crop tolerance to this herbicide. We aimed to quantify the degree of snap bean tolerance to flumioxazin and explore the underlying mechanisms. We investigated the genetic basis of herbicide tolerance using genome-wide association mapping approach utilizing field-collected data from a snap bean diversity panel, combined with gene expression data of cultivars with contrasting response. The response to a preemergence application of flumioxazin was measured by assessing plant population density and shoot biomass variables. Snap bean tolerance to flumioxazin is associated with a single genomic location in chromosome 02. Tolerance is influenced by several factors, including those that are indirectly affected by seed size/weight and those that directly impact the herbicide's metabolism and protect the cell from reactive oxygen species-induced damage. Transcriptional profiling and co-expression network analysis identified biological pathways likely involved in flumioxazin tolerance, including oxidoreductase processes and programmed cell death. Transcriptional regulation of genes involved in those processes is possibly orchestrated by a transcription factor located in the region identified in the GWAS analysis. Several entries belonging to the Romano class, including Bush Romano 350, Roma II, and Romano Purpiat presented high levels of tolerance in this study. The alleles identified in the diversity panel that condition snap bean tolerance to flumioxazin shed light on a novel mechanism of herbicide tolerance and can be used in crop improvement.

Co-expression analysis reveals distinct alliances around two carbon fixation pathways in hydrothermal vent symbionts

Most autotrophic organisms possess a single carbon fixation pathway. The chemoautotrophic symbionts of the hydrothermal vent tubeworm Riftia pachyptila, however, possess two functional pathways: the Calvin–Benson–Bassham (CBB) and the reductive tricarboxylic acid (rTCA) cycles. How these two pathways are coordinated is unknown. Here we measured net carbon fixation rates, transcriptional/metabolic responses and transcriptional co-expression patterns of Riftia pachyptila endosymbionts by incubating tubeworms collected from the East Pacific Rise at environmental pressures, temperature and geochemistry. Results showed that rTCA and CBB transcriptional patterns varied in response to different geochemical regimes and that each pathway is allied to specific metabolic processes; the rTCA is allied to hydrogenases and dissimilatory nitrate reduction, whereas the CBB is allied to sulfide oxidation and assimilatory nitrate reduction, suggesting distinctive yet complementary roles in metabolic function. Furthermore, our network analysis implicates the rTCA and a group 1e hydrogenase as key players in the physiological response to limitation of sulfide and oxygen. Net carbon fixation rates were also exemplary, and accordingly, we propose that co-activity of CBB and rTCA may be an adaptation for maintaining high carbon fixation rates, conferring a fitness advantage in dynamic vent environments.

Antioxidant network-based signatures cluster glioblastoma into distinct redox-resistant phenotypes.

Aberrant reactive oxygen species (ROS) production is one of the hallmarks of cancer. During their growth and dissemination, cancer cells control redox signaling to support protumorigenic pathways. As a consequence, cancer cells become reliant on major antioxidant systems to maintain a balanced redox tone, while avoiding excessive oxidative stress and cell death. This concept appears especially relevant in the context of glioblastoma multiforme (GBM), the most aggressive form of brain tumor characterized by significant heterogeneity, which contributes to treatment resistance and tumor recurrence. From this viewpoint, this study aims to investigate whether gene regulatory networks can effectively capture the diverse redox states associated with the primary phenotypes of GBM. In this study, we utilized publicly available GBM datasets along with proprietary bulk sequencing data. Employing computational analysis and bioinformatics tools, we stratified GBM based on their antioxidant capacities and evaluated the distinctive functionalities and prognostic values of distinct transcriptional networks in silico. We established three distinct transcriptional co-expression networks and signatures (termed clusters C1, C2, and C3) with distinct antioxidant potential in GBM cancer cells. Functional analysis of each cluster revealed that C1 exhibits strong antioxidant properties, C2 is marked with a discrepant inflammatory trait and C3 was identified as the cluster with the weakest antioxidant capacity. Intriguingly, C2 exhibited a strong correlation with the highly aggressive mesenchymal subtype of GBM. Furthermore, this cluster holds substantial prognostic importance: patients with higher gene set variation analysis (GSVA) scores of the C2 signature exhibited adverse outcomes in overall and progression-free survival. In summary, we provide a set of transcriptional signatures that unveil the antioxidant potential of GBM, offering a promising prognostic application and a guide for therapeutic strategies in GBM therapy.

Prognostic Implications of Small Cell Lung Cancer Transcriptional Subtyping for CNS Metastases.

Small cell lung cancer (SCLC) often metastasizes to the brain and has poor prognosis. SCLC subtypes distinguished by expressing transcriptional factors ASCL1 or NEUROD1 have been identified. This study investigates the impact of transcription factor-defined SCLC subtype on incidence and outcomes of brain metastases (BMs). Patients with SCLC with ASCL1 (A) and NEUROD1 (N) immunohistochemical expression status were identified and classified: (1) A+/N-, (2) A+/N+, (3) A-/N+, and (4) A-/N-. Cumulative incidence competing risk analyses were used to assess incidence of CNS progression. Cox proportional hazards models were used for multivariable analyses of overall survival (OS) and CNS progression-free survival (CNS-PFS). Of 164 patients, most were either A+/N- or A+/N+ (n = 62, n = 63, respectively). BMs were present at diagnosis in 24 patients (15%). Among them, the 12-month cumulative incidence of subsequent CNS progression was numerically highest for A+/N- (50% [95% CI, 10.5 to 74.7]; P = .47). Among those BM-free at diagnosis, the 12-month cumulative incidence of CNS progression was numerically the highest for A+/N- (16% [95% CI, 7.5 to 27.9]) and A-/N+ (9.1% [95% CI, 0.0 to 34.8]; P = .20). Both subtypes, A+/N- and A-/N+, had worse OS compared with A+/N+ (A+/N-: hazard ratio [HR], 1.62 [95% CI, 1.01 to 2.51]; P < .05; A-/N+: HR, 3.02 [95% CI, 1.35 to 6.76]; P = .007). Excellent response rates (28, 65% CR/PR) across subtypes were seen in patients who had CNS-directed radiotherapy versus systemic therapy alone (9, 36% CR/PR). To our knowledge, this report is the first to investigate CNS-specific outcomes based on transcription factor subtypes in patients with SCLC. BM-free patients at diagnosis with A+/N- or A-/N+ subtypes had worse outcomes compared with those with transcriptional factor coexpression. Further investigation into the mechanisms and implications of SCLC subtyping on CNS-specific outcomes is warranted to ultimately guide personalized care.

Breaking down barriers: comprehensive functional analysis of the Aspergillus niger chitin synthase repertoire

BackgroundMembers of the fungal kingdom are heterotrophic eukaryotes encased in a chitin containing cell wall. This polymer is vital for cell wall stiffness and, ultimately, cell shape. Most fungal genomes contain numerous putative chitin synthase encoding genes. However, systematic functional analysis of the full chitin synthase catalogue in a given species is rare. This greatly limits fundamental understanding and potential applications of manipulating chitin synthesis across the fungal kingdom.ResultsIn this study, we conducted in silico profiling and subsequently deleted all predicted chitin synthase encoding genes in the multipurpose cell factory Aspergillus niger. Phylogenetic analysis suggested nine chitin synthases evolved as three distinct groups. Transcript profiling and co-expression network construction revealed remarkably independent expression, strongly supporting specific role(s) for the respective chitin synthases. Deletion mutants confirmed all genes were dispensable for germination, yet impacted colony spore titres, chitin content at hyphal septa, and internal architecture of submerged fungal pellets. We were also able to assign specific roles to individual chitin synthases, including those impacting colony radial growth rates (ChsE, ChsF), lateral cell wall chitin content (CsmA), chemical genetic interactions with a secreted antifungal protein (CsmA, CsmB, ChsE, ChsF), resistance to therapeutics (ChsE), and those that modulated pellet diameter in liquid culture (ChsA, ChsB). From an applied perspective, we show chsF deletion increases total protein in culture supernatant over threefold compared to the control strain, indicating engineering filamentous fungal chitin content is a high priority yet underexplored strategy for strain optimization.ConclusionThis study has conducted extensive analysis for the full chitin synthase encoding gene repertoire of A. niger. For the first time we reveal both redundant and non-redundant functional roles of chitin synthases in this fungus. Our data shed light on the complex, multifaceted, and dynamic role of chitin in fungal growth, morphology, survival, and secretion, thus improving fundamental understanding and opening new avenues for biotechnological applications in fungi.

P176 TYK2 inhibition attenuates the stimulation of sensory neurons with hyper-IL-6

Abstract Background Abdominal pain and diarrhoea are a cause of significant morbidity for people with colitis. Recent single cell RNA sequencing of sensory neurons innervating the gut highlights a likely role for sensory nerves in the pathophysiology of these symptoms based on the marked expression of cytokine receptors from the IL-6 receptor family in gut projecting sensory neurons. In keeping with the reported downstream signalling pathway for these receptors marked co-expression of TYK2 transcripts were also found along with IL-6 signal transducer (GP130) in colonic nociceptors characterised by expression of the capsaicin receptor TRPV1. These observations suggest that TYK2 inhibitors may ameliorate symptoms of pain and diarrhoea in people with colitis by preventing neurogenically mediated gut contractility and nociceptor activation. Methods TYK2 mediated sensory neuron activation: changes in the intracellular [Ca2+]i within thoracolumbar (T12-L5 spinal segments) dorsal root ganglia sensory neurons were examined in cells cultured overnight and loaded with Fluo-4-AM (30min) prior to imaging. Responses were measured to the hyper-IL-6 fusion protein of the soluble IL-6 receptor and IL-6, which transactivates the GP130 receptor subunit to mimic IL-6 family receptor signalling following pre-treatment with vehicle, the TYK2 inhibitor Deucravacitinib (100nM) or the pan-JAK inhibitor Pyridone-6 (2 μM). Afterwards cells were exposed to capsaicin (1 µM) followed by 50 mM KCl, to confirm the presence of nociceptors and neurons respectively by positive response. All experiments were performed using tissue from male C57/B6 mice euthanised by rising concentration of CO2followed by exsanguination in accordance with Schedule 1 of the UK Animals Scientific Procedures act (1986). Results Application of hyper-IL-6 produced a robust increase in [Ca2+]i levels in DRG sensory neurons which was significantly attenuated by pre-treatment with the TYK2 specific inhibitor Deucravacitinib (p&amp;lt;0.01) or the pan JAK inhibitor Pyridone-6 (p = &amp;lt;0.05). Conclusion Findings demonstrate that inhibition of TYK2 inhibitor prevents sensory neuron activation in response to IL-6 trans-signalling. This data suggests that TYK2 inhibitors may have utility for the treatment of abdominal pain in people with colitis.

Molecular mechanisms of microbiome modulation by the eukaryotic secondary metabolite azelaic acid.

Photosynthetic eukaryotes, such as microalgae and plants, foster fundamentally important relationships with their microbiome based on the reciprocal exchange of chemical currencies. Among these, the dicarboxylate metabolite azelaic acid (Aze) appears to play an important, but heterogeneous, role in modulating these microbiomes, as it is used as a carbon source for some heterotrophs but is toxic to others. However, the ability of Aze to promote or inhibit growth, as well as its uptake and assimilation mechanisms into bacterial cells are mostly unknown. Here, we use transcriptomics, transcriptional factor coexpression networks, uptake experiments, and metabolomics to unravel the uptake, catabolism, and toxicity of Aze on two microalgal-associated bacteria, Phycobacter and Alteromonas, whose growth is promoted or inhibited by Aze, respectively. We identify the first putative Aze transporter in bacteria, a 'C4-TRAP transporter', and show that Aze is assimilated through fatty acid degradation, with further catabolism occurring through the glyoxylate and butanoate metabolism pathways when used as a carbon source. Phycobacter took up Aze at an initial uptake rate of 3.8×10-9 nmol/cell/hr and utilized it as a carbon source in concentrations ranging from 10 μM to 1 mM, suggesting a broad range of acclimation to Aze availability. For growth-impeded bacteria, we infer that Aze inhibits the ribosome and/or protein synthesis and that a suite of efflux pumps is utilized to shuttle Aze outside the cytoplasm. We demonstrate that seawater amended with Aze becomes enriched in bacterial families that can catabolize Aze, which appears to be a different mechanism from that in soil, where modulation by the host plant is required. This study enhances our understanding of carbon cycling in the oceans and how microscale chemical interactions can structure marine microbial populations. In addition, our findings unravel the role of a key chemical currency in the modulation of eukaryote-microbiome interactions across diverse ecosystems.

Molecular mechanisms of microbiome modulation by the eukaryotic secondary metabolite azelaic acid

Photosynthetic eukaryotes, such as microalgae and plants, foster fundamentally important relationships with their microbiome based on the reciprocal exchange of chemical currencies. Among these, the dicarboxylate metabolite azelaic acid (Aze) appears to play an important, but heterogeneous, role in modulating these microbiomes, as it is used as a carbon source for some heterotrophs but is toxic to others. However, the ability of Aze to promote or inhibit growth, as well as its uptake and assimilation mechanisms into bacterial cells are mostly unknown. Here, we use transcriptomics, transcriptional factor coexpression networks, uptake experiments, and metabolomics to unravel the uptake, catabolism, and toxicity of Aze on two microalgal-associated bacteria, Phycobacter and Alteromonas, whose growth is promoted or inhibited by Aze, respectively. We identify the first putative Aze transporter in bacteria, a ‘C4-TRAP transporter’, and show that Aze is assimilated through fatty acid degradation, with further catabolism occurring through the glyoxylate and butanoate metabolism pathways when used as a carbon source. Phycobacter took up Aze at an initial uptake rate of 3.8×10–9 nmol/cell/hr and utilized it as a carbon source in concentrations ranging from 10 μM to 1 mM, suggesting a broad range of acclimation to Aze availability. For growth-impeded bacteria, we infer that Aze inhibits the ribosome and/or protein synthesis and that a suite of efflux pumps is utilized to shuttle Aze outside the cytoplasm. We demonstrate that seawater amended with Aze becomes enriched in bacterial families that can catabolize Aze, which appears to be a different mechanism from that in soil, where modulation by the host plant is required. This study enhances our understanding of carbon cycling in the oceans and how microscale chemical interactions can structure marine microbial populations. In addition, our findings unravel the role of a key chemical currency in the modulation of eukaryote-microbiome interactions across diverse ecosystems.

Megakaryocytes Have Immune Characteristics in Murine Immune Bone Marrow Failure

Introduction Megakaryocytes (MKs) have familiar roles in the bone marrow (BM) microenvironment, including platelet production, maintenance of hematopoietic stem and progenitor cells (HSPCs) function, and more recent evidence suggests MK engagement in immune responses. Immune aplastic anemia (AA) is characterized by oligoclonal expansion of T-cells and death of HSPCs leading to bone marrow failure (BMF). As MKs demonstrate niche maintenance and immune effector function, we explored the potential roles of MKs in the development of AA in a murine model of immune BMF. Methods BMF was induced in CByB6F1 mice by intravenous infusion of lymph node (LN) cells from major histocompatibility complex mismatched C57BL/B6 (B6) donor mice at a dose of 30-40x10 6 cells/recipient. No irradiation was given to recipients to avoid extra damage to MKs and CByB6F1 recipient mice were sacrificed at D7, D10, D14, and D17 following LN cell infusion (Chen J, et al. Blood. 2004). Peripheral blood (PB) samples, tibias, femurs, and sterna were harvested at each timepoint. Immunophenotyping of MKs and T cells was performed by flow cytometry on flushed BM cells and two-photon excitation microscopy (2P) imaging was performed on fixed sternal samples, using fluorescence-labelled CD41 antibody to visualize the physical location of MKs in BM. To enrich for MKs, BM cells were filtered through a 10 µm pluriStrainer™ through which most lymphocytes, red blood cells, and other small cells were eluted while MKs and other large size cells remained on the strainer. Large cells and residual small cells were subsequently assessed by single cell RNA sequencing (sc-RNAseq) on the Parse Bio Evercode™ platform based on combinatorial barcoding technology, which was able to accommodate the large size of MKs, as no partitioning step is needed in the analysis workflow. Results BMF was successfully induced as measured by decreasing PB counts from point of lymphocyte infusion accompanied by increasing expansion of CD4 and CD8 T-cells and marked reduction of total BM cells. MKs persisted throughout the course of BMF as measured by both flow cytometry and 2P (Figure 1A) despite PB thrombocytopenia and progressive marrow hypocellularity. Concomitantly, MKs identified by large size on forward scatter versus side scatter and CD41+ expression had progressive upregulation of immune effector markers IA-IE and CD53 over the course of BMF. Following sc-RNAseq and bioinformatics analysis of enriched BM cells at D14, cell clustering identified a distinct cell cluster representing MKs (Figure 1B). MK identity was supported by co-expression of transcripts encoding platelet factor 4, TUBB1, Glycoprotein 5, and Itga2b (CD41). This MK cell population was distinct from macrophages which highly expressed Fcgr4, Cd5l, Ccr5, Cd74 and Gbp5 genes. Gene expression profiling demonstrated that MKs in BMF mice showed upregulation of immune pathways such as interferon gamma and interferon alpha response and downregulation of pathways related to platelet activation and platelet homeostasis, as compared with normal mice. Conclusions In a murine model of BMF, using multiple modalities, we demonstrate that MKs cease platelet production, manifested by reduction in peripheral blood platelet count and downregulation of gene expression pathways related to platelet production, while concurrently acquiring characteristics of immune effector cells. This suggests that MKs are active participants in immune BMF.

Identification and characterization of the cupin_1 domain-containing proteins in ma bamboo (Dendrocalamus latiflorus) and their potential role in rhizome sprouting.

Cupin_1 domain-containing protein (CDP) family, which is a member of the cupin superfamily with the most diverse functions in plants, has been found to be involved in hormone pathways that are closely related to rhizome sprouting (RS), a vital form of asexual reproduction in plants. Ma bamboo is a typical clumping bamboo, which mainly reproduces by RS. In this study, we identified and characterized 53 Dendrocalamus latiflorus CDP genes and divided them into seven subfamilies. Comparing the genetic structures among subfamilies showed a relatively conserved gene structure within each subfamily, and the number of cupin_1 domains affected the conservation among D. latiflorus CDP genes. Gene collinearity results showed that segmental duplication and tandem duplication both contributed to the expansion of D. latiflorus CDP genes, and lineage-specific gene duplication was an important factor influencing the evolution of CDP genes. Expression patterns showed that CDP genes generally had higher expression levels in germinating underground buds, indicating that they might play important roles in promoting shoot sprouting. Transcription factor binding site prediction and co-expression network analysis indicated that D. latiflorus CDPs were regulated by a large number of transcription factors, and collectively participated in rhizome buds and shoot development. This study significantly provided new insights into the evolutionary patterns and molecular functions of CDP genes, and laid a foundation for further studying the regulatory mechanisms of plant rhizome sprouting.

Stress combinations and their interactions in plants database: a one-stop resource on combined stress responses in plants.

We have developed a compendium and interactive platform, named Stress Combinations and their Interactions in Plants Database (SCIPDb; http://www.nipgr.ac.in/scipdb.php), which offers information on morpho-physio-biochemical (phenome) and molecular (transcriptome and metabolome) responses of plants to different stress combinations. SCIPDb is a plant stress informatics hub for data mining on phenome, transcriptome, trait-gene ontology, and data-driven research for advancing mechanistic understanding of combined stress biology. We analyzed global phenome data from 939 studies to delineate the effects of various stress combinations on yield in major crops and found that yield was substantially affected under abiotic-abiotic stresses. Transcriptome datasets from 36 studies hosted in SCIPDb identified novel genes, whose roles have not been earlier established in combined stress. Integretome analysis under combined drought-heat stress pinpointed carbohydrate, amino acid, and energy metabolism pathways as the crucial metabolic, proteomic, and transcriptional components in plant tolerance to combined stress. These examples illustrate the application of SCIPDb in identifying novel genes and pathways involved in combined stress tolerance. Further, we showed the application of this database in identifying novel candidate genes and pathways for combined drought and pathogen stress tolerance. To our knowledge, SCIPDb is the only publicly available platform offering combined stress-specific omics big data visualization tools, such as an interactive scrollbar, stress matrix, radial tree, global distribution map, meta-phenome analysis, search, BLAST, transcript expression pattern table, Manhattan plot, and co-expression network. These tools facilitate a better understanding of the mechanisms underlying plant responses to combined stresses.

Higher incidence of co-expression of BCR-ABL fusion transcripts in an Eastern Indian population

BackgroundChronic myeloid leukaemia (CML) is a haematopoietic stem cell disorder, caused by a balanced reciprocal translocation (t(9;22) (q34;q11)) that leads to the formation of BCR (Breakpoint Cluster Region)-ABL (Abelson) fusion transcripts known as Philadelphia (Ph) chromosome. The prevalence of BCR-ABL fusion transcripts in Indian CML population is poorly understood, and few studies have been reported from India. The aim of the present study was to determine the frequencies as well as prognostic effects of the three fusion transcripts, i.e. b2a2, b3a2 and e1a2 in an Indian population.MethodsRNA was isolated from total 123 samples, 27 bone marrow (BM) samples and 96 peripheral blood (PB) samples, of CML patient followed by cDNA synthesis. Real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed using TaqMan® assay (ABI, CA, USA) to monitor BCR-ABL transcript.ResultsPh' chromosome was observed in 103 patients whereas it was not detected in 20 cases. qRT-PCR revealed that the b3a2 fusion transcripts were the most common transcript in CML patients (63.41%) while b2a2 fusion transcript was present in 16.26% cases. Co-expression of b3a2 + b2a2 fusion transcript was observed in 0.81% cases whereas co-expression of b3a2 + e1a2 fusion transcript was found in 1.63% cases. There was no correlation observed between b3a2 fusion transcript and platelet count. The fusion transcript b2a2 was observed in relatively younger patients compared to b3a2 fusion transcript. Although this correlation was not statistically significant.ConclusionThe co-expression of BCR-ABL fusion transcripts was higher (63.41% aggregate of b3a2) in the present population in contrast with other populations reported. This finding was consistent with the frequency data reported from Sudan.

Genetic alterations in the BCR-ABL1 fusion gene related to imatinib resistance in chronic myeloid leukemia

Use of the potent tyrosine kinase inhibitor imatinib as the first-line treatment in chronic myeloid leukemia (CML) has decreased mortality from 20% to 2%. Approximately 30% of CML patients experience imatinib resistance, however, largely because of point mutations in the kinase domain of the BCR-ABL1 fusion gene. The aim of this study was to use next-generation sequencing (NGS) to identify mutations related to imatinib resistance. The study included 22 patients diagnosed with CML and experiencing no clinical response to imatinib. Total RNA was used for cDNA synthesis, with amplification of a fragment encompassing the BCR-ABL1 kinase domain using a nested-PCR approach. Sanger and NGS were applied to detect genetic alterations. HaplotypeCaller was used for variant calling, and STAR-Fusion software was applied for fusion breakpoint identification. After sequencing analysis, F311I, F317L, and E450K mutations were detected respectively in three different participants, and in another two patients, single nucleotide variants in BCR (rs9608100, rs140506, rs16802) and ABL1 (rs35011138) were detected. Eleven patients carried e14a2 transcripts, nine had e13a2 transcripts, and both transcripts were identified in one patient. One patient had co-expression of e14a2 and e14a8 transcripts. The results identify candidate single nucleotide variants and co-expressed BCR-ABL1 transcripts in cellular resistance to imatinib.

WMDS.net: a network control framework for identifying key players in transcriptome programs.

Mammalian cells can be transcriptionally reprogramed to other cellular phenotypes. Controllability of such complex transitions in transcriptional networks underlying cellular phenotypes is an inherent biological characteristic. This network controllability can be interpreted by operating a few key regulators to guide the transcriptional program from one state to another. Finding the key regulators in the transcriptional program can provide key insights into the network state transition underlying cellular phenotypes. To address this challenge, here, we proposed to identify the key regulators in the transcriptional co-expression network as a minimum dominating set (MDS) of driver nodes that can fully control the network state transition. Based on the theory of structural controllability, we developed a weighted MDS network model (WMDS.net) to find the driver nodes of differential gene co-expression networks. The weight of WMDS.net integrates the degree of nodes in the network and the significance of gene co-expression difference between two physiological states into the measurement of node controllability of the transcriptional network. To confirm its validity, we applied WMDS.net to the discovery of cancer driver genes in RNA-seq datasets from The Cancer Genome Atlas. WMDS.net is powerful among various cancer datasets and outperformed the other top-tier tools with a better balance between precision and recall. https://github.com/chaofen123/WMDS.net. Supplementary data are available at Bioinformatics online.

Transcriptome analysis of sugarcane reveals rapid defense response of SES208 to Xanthomonas albilineans in early infection

BackgroundDiseases are the major factor affecting the quality and yield of sugarcane during its growth and development. However, our knowledge about the factors regulating disease responses remain limited. The present study focuses on identifying genes regulating transcriptional mechanisms responsible for resistance to leaf scald caused by Xanthomonas albilineans in S. spontaneum and S. officinarum.ResultsAfter inoculation of the two sugarcane varieties SES208 (S. spontaneum) and LA Purple (S. officinarum) with Xanthomonas albilineans, SES208 exhibited significantly greater resistance to leaf scald caused by X. albilineans than did LA Purple. Using transcriptome analysis, we identified a total of 4323 and 1755 differentially expressed genes (DEGs) in inoculated samples of SES208 and LA Purple, respectively. Significantly, 262 DEGs were specifically identified in SES208 that were enriched for KEGG pathway terms such as plant-pathogen interaction, MAPK signaling pathway, and plant hormone signal transduction. Furthermore, we built a transcriptional regulatory co-expression network that specifically identified 16 and 25 hub genes in SES208 that were enriched for putative functions in plant-pathogen interactions, MAPK signaling, and plant hormone signal transduction. All of these essential genes might be significantly involved in resistance-regulating responses in SES208 after X. albilineans inoculation. In addition, we found allele-specific expression in SES208 that was associated with the resistance phenotype of SES208 when infected by X. albilineans. After infection with X. albilineans, a great number of DEGs associated with the KEGG pathways ‘phenylpropanoid biosynthesis’ and ‘flavonoid biosynthesis’ exhibited significant expression changes in SES208 compared to LA Purple that might contribute to superior leaf scald resistance in SES208.ConclusionsWe provided the first systematical transcriptome map that the higher resistance of SES208 is associated with and elicited by the rapid activation of multiple clusters of defense response genes after infection by X. albilineans and not merely due to changes in the expression of genes generically associated with stress resistance. These results will serve as the foundation for further understanding of the molecular mechanisms of resistance against X. albilineans in S. spontaneum.

The Pattern of Chromosome Aberrations and Molecular Markers in the Population of Hematological Patients Diagnosed at the University Clinical Center Tuzla

Introduction: Besides cardiovascular, malignant diseases are one of the leading causes of death in Bosnia and Herzegovina. At the top of this list are hematological diseases. This research aimed to identify cytogenetic and molecular biomarkers in patients treated for different types of hematological neoplasms. Methods: The retrospective study included 1600 samples of patients with different hematological diseases in the period from January 2006 to May 2022. The Polymerase Chain Reaction (RT-PCR) method was used to determine the presence of genetic rearrangements and to confirm the findings of conventional cytogenetic analysis. Results: Chromosomal aberrations were found in 739 (46,18%) patients. Using the RT-PCR technique, positive cases were increased by 1,5%. The BCR-ABL fusion gene was present in e14-a2 transcript form in 73% of samples, e13-a2 isoform in 21%, e1-a2 in 2%, while e14-a2/e1-a2 transcript coexpression was present in a percentage of 4% of the samples. The PML-RARa fusion gene was found in the form of bcr 1 transcripts in 21%, bcr2 32% and bcr3 59% of the samples. In twelve cases A type of the CBFB-MYH11 fusion transcript was detected. The MLL-AF4 fusion was found in only one case. Conclusion: The obtained percentages of frequency of individual molecular gene isoforms are in accordance with the results of most other researchers. This refers to the Balkan population and the Caucasian ethnic group.

Probing Folate-Responsive and Stage-Sensitive Metabolomics and Transcriptional Co-Expression Network Markers to Predict Prognosis of Non-Small Cell Lung Cancer Patients

Tumour metabolomics and transcriptomics co-expression network as related to biological folate alteration and cancer malignancy remains unexplored in human non-small cell lung cancers (NSCLC). To probe the diagnostic biomarkers, tumour and pair lung tissue samples (n = 56) from 97 NSCLC patients were profiled for ultra-performance liquid chromatography tandem mass spectrometry (UPLC/MS/MS)-analysed metabolomics, targeted transcriptionomics, and clinical folate traits. Weighted Gene Co-expression Network Analysis (WGCNA) was performed. Tumour lactate was identified as the top VIP marker to predict advance NSCLC (AUC = 0.765, Sig = 0.017, CI 0.58-0.95). Low folate (LF)-tumours vs. adjacent lungs displayed higher glycolytic index of lactate and glutamine-associated amino acids in enriched biological pathways of amino sugar and glutathione metabolism specific to advance NSCLCs. WGCNA classified the green module for hub serine-navigated glutamine metabolites inversely associated with tumour and RBC folate, which module metabolites co-expressed with a predominant up-regulation of LF-responsive metabolic genes in glucose transport (GLUT1), de no serine synthesis (PHGDH, PSPH, and PSAT1), folate cycle (SHMT1/2 and PCFR), and down-regulation in glutaminolysis (SLC1A5, SLC7A5, GLS, and GLUD1). The LF-responsive WGCNA markers predicted poor survival rates in lung cancer patients, which could aid in optimizing folate intervention for better prognosis of NSCLCs susceptible to folate malnutrition.

GeneFriends: gene co-expression databases and tools for humans and model organisms.

Gene co-expression analysis has emerged as a powerful method to provide insights into gene function and regulation. The rapid growth of publicly available RNA-sequencing (RNA-seq) data has created opportunities for researchers to employ this abundant data to help decipher the complexity and biology of genomes. Co-expression networks have proven effective for inferring the relationship between the genes, for gene prioritization and for assigning function to poorly annotated genes based on their co-expressed partners. To facilitate such analyses we created previously an online co-expression tool for humans and mice entitled GeneFriends. To continue providing a valuable tool to the scientific community, we have now updated the GeneFriends database and website. Here, we present the new version of GeneFriends, which includes gene and transcript co-expression networks based on RNA-seq data from 46 475 human and 34 322 mouse samples. The new database also encompasses tissue-specific gene co-expression networks for 20 human and 21 mouse tissues, dataset-specific gene co-expression maps based on TCGA and GTEx projects and gene co-expression networks for additional seven model organisms (fruit fly, zebrafish, worm, rat, yeast, cow and chicken). GeneFriends is freely available at http://www.genefriends.org/.

Transcriptome and Phenotype Integrated Analysis Identifies Genes Controlling Ginsenoside Rb1 Biosynthesis and Reveals Their Interactions in the Process in Panax ginseng.

Genes are the keys to deciphering the molecular mechanism underlying a biological trait and designing approaches desirable for plant genetic improvement. Ginseng is an important medicinal herb in which ginsenosides have been shown to be the major bioactive component; however, only a few genes involved in ginsenoside biosynthesis have been cloned through orthologue analysis. Here, we report the identification of 21 genes controlling Rb1 biosynthesis by stepwise ginseng transcriptome and Rb1 content integrated analysis. We first identified the candidate genes for Rb1 biosynthesis by integrated analysis of genes with the trait from four aspects, including gene transcript differential expression between highest- and lowest-Rb1 content cultivars, gene transcript expression-Rb1 content correlation, and biological impacts of gene mutations on Rb1 content, followed by the gene transcript co-expression network. Twenty-two candidate genes were identified, of which 21 were functionally validated for Rb1 biosynthesis by gene regulation, genetic transformation, and mutation analysis. These genes were strongly correlated in expression with the previously cloned genes encoding key enzymes for Rb1 biosynthesis. Based on the correlations, a pathway for Rb1 biosynthesis was deduced to indicate the roles of the genes in Rb1 biosynthesis. Moreover, the genes formed a strong co-expression network with the previously cloned Rb1 biosynthesis genes, and the variation in the network was associated with the variation in the Rb1 content. These results indicate that Rb1 biosynthesis is a process of correlative interactions among Rb1 biosynthesis genes. Therefore, this study provides new knowledge, 21 new genes, and 96 biomarkers for Rb1 biosynthesis useful for enhanced research and breeding in ginseng.

Rational Approach to Finding Genes Encoding Molecular Biomarkers: Focus on Breast Cancer.

Early detection of cancer facilitates treatment and improves patient survival. We hypothesized that molecular biomarkers of cancer could be rationally predicted based on even partial knowledge of transcriptional regulation, functional pathways and gene co-expression networks. To test our data mining approach, we focused on breast cancer, as one of the best-studied models of this disease. We were particularly interested to check whether such a ‘guilt by association’ approach would lead to pan-cancer markers generally known in the field or whether molecular subtype-specific ‘seed’ markers will yield subtype-specific extended sets of breast cancer markers. The key challenge of this investigation was to utilize a small number of well-characterized, largely intracellular, breast cancer-related proteins to uncover similarly regulated and functionally related genes and proteins with the view to predicting a much-expanded range of disease markers, especially that of extracellular molecular markers, potentially suitable for the early non-invasive detection of the disease. We selected 23 previously characterized proteins specific to three major molecular subtypes of breast cancer and analyzed their established transcription factor networks, their known metabolic and functional pathways and the existing experimentally derived protein co-expression data. Having started with largely intracellular and transmembrane marker ‘seeds’ we predicted the existence of as many as 150 novel biomarker genes to be associated with the selected three major molecular sub-types of breast cancer all coding for extracellularly targeted or secreted proteins and therefore being potentially most suitable for molecular diagnosis of the disease. Of the 150 such predicted protein markers, 114 were predicted to be linked through the combination of regulatory networks to basal breast cancer, 48 to luminal and 7 to Her2-positive breast cancer. The reported approach to mining molecular markers is not limited to breast cancer and therefore offers a widely applicable strategy of biomarker mining.