Genome-wide Transcriptome Profiling Research Articles

Testing the Significance of Ranked Gene Sets in Genome-wide Transcriptome Profiling Data Using Weighted Rank Correlation Statistics

Background: Popular gene set enrichment analysis approaches assumed that genes in the gene set contributed to the statistics equally. However, the genes in the transcription factors (TFs) derived gene sets, or gene sets constructed by TF targets identified by the ChIP-Seq experiment, have a rank attribute, as each of these genes have been assigned with a p-value which indicates the true or false possibilities of the ownerships of the genes belong to the gene sets. Objective: Popular gene set enrichment analysis approaches assumed that genes in the gene set contributed to the statistics equally. However, the genes in the transcription factors (TFs) derived gene sets, or gene sets constructed by TF targets identified by the ChIP-Seq experiment, have a rank attribute, as each of these genes have been assigned with a p-value which indicates the true or false possibilities of the ownerships of the genes belong to the gene sets. Methods: A method was proposed by first creating ranked gene sets and gene lists and then applying weighted Kendall's tau rank correlation statistics to the test. After introducing top-down weights to the genes in the gene set, a new software called "Flaver" was developed. Results: Theoretical properties of the proposed method were established, and its differences over the GSEA approach were demonstrated when analyzing the transcriptome profiling data across 55 human tissues and 176 human cell-lines. The results indicated that the TFs identified by our method have higher tendency to be differentially expressed across the tissues analyzed than its competitors. It significantly outperforms the well-known gene set enrichment analyzing tools, GOStats (9%) and GSEA (17%), in analyzing well-documented human RNA transcriptome datasets. Conclusions: The method is outstanding in detecting gene sets of which the gene ranks were correlated with the expression levels of the genes in the transcriptome data.

Cotton SNARE complex component GhSYP121 regulates salicylic acid signaling during defense against Verticillium dahliae.

Syntaxin of plant (SYP) plays a crucial role in SNARE-mediated membrane trafficking during endocytic and secretory pathways, contributing to the regulation and execution of plant immunity against pathogens. Verticillium wilt is among the most destructive fungal diseases affecting cotton worldwide. However, information regarding SYP family genes in cotton is scarce. Through genome-wide identification and transcriptome profiling, we identified GhSYP121, a Qa SNARE gene in Gossypium hirsutum. GhSYP121 is notably induced by Verticillium dahliae, the causal agent of Verticillium wilt in cotton, and acts as a negative regulator of defense against V. dahliae. This is evidenced by the reduced resistance of GhSYP121-deficient cotton and the increased susceptibility of GhSYP121-overexpressing lines. Furthermore, the activation of the salicylic acid (SA) pathway by V. dahliae is inversely correlated with the expression level of GhSYP121. GhSYP121 interacts with its cognate SNARE component, GhSNAP33, which is required for the penetration resistance against V. dahliae in cotton. Collectively, GhSYP121, as a member of the cotton SNARE complex, is involved in regulating the SA pathway during plant defense against V. dahliae. This finding enhances our understanding of the potential role of GhSYP121 in these distinct pathways that contribute to plant defense against V. dahliae infection.

Genome-wide DNA methylome and transcriptome profiling reveals key genes involved in the dysregulation of adipose-stem cells in Crohn's disease.

Crohn's disease (CD) is characterised by the expansion of mesenteric adipose tissue (MAT), named creeping fat (CF), which seems to be directly related to disease activity. Adipose-stem cells (ASCs) isolated from the CF of patients with CD are extremely pro-inflammatory, which persists during disease remission. We hypothesised that the dysfunctional ASCs in CD accumulate epigenetic modifications triggered by the inflammatory environment that could serve as molecular markers. Genome-wide DNA methylome and transcriptome profiling were performed in ASCs isolated from MAT adipose-tissue biopsies of patients with active and inactive disease and from non-Crohn's disease patients (non-CD). A validation cohort was used to test the main candidate genes via qPCR in other fat depots and immune cells. We found differences in DNA-methylation and gene expression between ASCs isolated from patients with CD and from non-CD subjects, but we found no differences related to disease activity. Pathway enrichment analysis revealed that oxidative stress and immune response were significantly enriched in active CD and integration analysis identified MAB21L2, a cell fate-determining gene, as the most affected gene in CD. Validation analysis confirmed the elevated gene expression of MAB21L2 in MAT and in adipose tissue macrophages in active CD. We also found a strong association between expression of the calcium channel subunit gene CACNA1H and disease remission, as CACNA1H expression was higher in ASCs and MAT from patients with inactive CD, and correlates negatively with C-reactive protein in peripheral blood mononuclear cells. We identified a potential gene signature of CD in ASCs obtained from MAT. Integration analysis highlighted two novel genes demonstrating a negative correlation between promoter DNA methylation and transcription: one linked to ASCs in CD (MAB21L2) and the other (CACNA1H) related to disease remission.

TabHLH27 orchestrates root growth and drought tolerance to enhance water use efficiency in wheat.

Cultivating high-yield wheat under limited water resources is crucial for sustainable agriculture in semiarid regions. Amid water scarcity, plants activate drought response signaling, yet the delicate balance between drought tolerance and development remains unclear. Through genome-wide association studies and transcriptome profiling, we identified a wheat atypical basic helix-loop-helix (bHLH) transcription factor (TF), TabHLH27-A1, as a promising quantitative trait locus candidate for both relative root dry weight and spikelet number per spike in wheat. TabHLH27-A1/B1/D1 knock-out reduced wheat drought tolerance, yield, and water use efficiency (WUE). TabHLH27-A1 exhibited rapid induction with polyethylene glycol (PEG) treatment, gradually declining over days. It activated stress response genes such as TaCBL8-B1 and TaCPI2-A1 while inhibiting root growth genes like TaSH15-B1 and TaWRKY70-B1 under short-term PEG stimulus. The distinct transcriptional regulation of TabHLH27-A1 involved diverse interacting factors such as TaABI3-D1 and TabZIP62-D1. Natural variations of TabHLH27-A1 influence its transcriptional responses to drought stress, with TabHLH27-A1Hap-II associated with stronger drought tolerance, larger root system, more spikelets, and higher WUE in wheat. Significantly, the excellent TabHLH27-A1Hap-II was selected during the breeding process in China, and introgression of TabHLH27-A1Hap-II allele improved drought tolerance and grain yield, especially under water-limited conditions. Our study highlights TabHLH27-A1's role in balancing root growth and drought tolerance, providing a genetic manipulation locus for enhancing WUE in wheat.

Genome-wide characterization, transcriptome profiling, and functional analysis of the ALMT gene family in Medicago for aluminum resistance

Aluminum (Al) is the major limiting factor affecting plant productivity in acidic soils. Al3+ ions exhibit increased solubility at a pH below 5, leading to plant root tip toxicity. Alternatively, plants can perceive very low concentrations of Al3+, and Al triggers downstream signaling even at pH 5.7 without causing Al toxicity. The ALUMINUM-ACTIVATED-MALATE-TRANSPORTER (ALMT) family members act as anion channels, with some regulating the secretion of malate from root apices to chelate Al, which is a crucial mechanism for plant Al resistance. To date, the role of the ALMT gene family within the legume Medicago species has not been fully characterized. In this study, we investigated the ALMT gene family in M. sativa and M. truncatula and identified 68 MsALMTs and 18 MtALMTs, respectively. Phylogenetic analysis classified these genes into five clades, and synteny analysis uncovered genuine paralogs and orthologs. The real-time quantitative reverse transcription PCR (qRT-PCR) analysis revealed that MtALMT8, MtALMT9, and MtALMT15 in clade 2-2b are expressed in both roots and root nodules, and MtALMT8 and MtALMT9 are significantly upregulated by Al in root tips. We also observed that MtALMT8 and MtALMT9 can partially restore the Al sensitivity of Atalmt1 in Arabidopsis. Moreover, transcriptome analysis examined the expression patterns of these genes in M. sativa in response to Al at both pH 5.7 and pH 4.6, as well as to protons, and found that Al and protons can independently induce some Al-resistance genes. Overall, our findings indicate that MtALMT8 and MtALMT9 may play a role in Al resistance, and highlight the resemblance between the ALMT genes in Medicago species and those in Arabidopsis.

Abstract 7192: Teriflunomide overcomes gemcitabine resistance by inhibiting the GSK3β/mTOR pathway in pancreatic cancer

Abstract Background: Pancreatic cancer (PC) is a lethal malignancy that is partly attributable to the acquisition of chemoresistance to conventional chemotherapies. Hence, developing strategies to overcome such chemoresistance is essential. Teriflunomide, an active metabolite of leflunomide, has been widely utilized in managing autoimmune disorders, and emerging evidence suggests that it can augment the efficacy of chemotherapy in cancer. We investigated the therapeutic effect of teriflunomide in overcoming resistance to Gemcitabine-based therapy in PC. Methods: We performed docking and genomewide transcriptomic profiling analyses to identify therapeutic targets of Teriflunomide associated with Gemcitabine resistance (Gem-R) in PC. Subsequently, we analyzed the TCGA dataset to determine the prognostic significance of candidate genes in PC patients. A series of experiments were performed in PC cells, and the key findings were validated in patient-derived 3D organoids and a Gem-R Mia-PaCa2 xenograft model. Results: Docking and transcriptomic profiling studies identified GSK3β, MAPK8, and MAPK9 as potential targets of Teriflunomide, which were significantly upregulated in Gem-R cell lines (P < 0.05). GSK3β was independently associated with poor overall and disease-free survival in PC patients. Combined Teriflunomide and Gemcitabine treatment exhibited superior synergistic anti-tumor potential, as evidenced by reducing PC cell proliferation, clonogenicity, invasion, and migration. The combined treatment with Teriflunomide and Gemcitabine primarily inhibited GSK3β phosphorylation, resulting in the silencing of the downstream AKT/mTOR pathway. Lentiviral-based silencing of GSK3β inhibited the AKT/mTOR pathway and reduced the IC50s of Gemcitabine in PC cells. Furthermore, co-administrating both drugs to Gem-R Mia-PaCa-2 xenograft mice significantly reduced tumor growth rates and volumes. These findings were recapitulated in patient-derived 3D-organoids, where the combined treatment resulted in fewer and smaller organoids vs. each drug individually (P < 0.05). Conclusions: We provide novel evidence that Teriflunomide helps overcome resistance to Gemcitabine by inhibiting the GSK3β/AKT/mTOR pathway, which could improve treatment outcomes in PC. Citation Format: Caiming Xu, Yuan Li, Silei Sui, Meiling Chen, Laleh Melstorm, Haiyong Han, Steven T. Rosen, Ajay Goel. Teriflunomide overcomes gemcitabine resistance by inhibiting the GSK3β/mTOR pathway in pancreatic cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7192.

Abstract 523: Aronia berry and oligomeric proanthocyanidins synergistically exhibit anticancer effects in colorectal cancer by inhibiting LMNB1

Abstract Background: Colorectal cancer (CRC) is one of the most common and recurrent diseases globally, ranking as the second leading cause of mortality. The high mortality rates associated with this disease are largely due to the ineffectiveness of current treatments and the development of chemotherapeutic resistance to conventional therapies. In recent years, several natural agents, such as Aronia berry extract (ABE) and oligomeric proanthocyanidins (OPCs), have shown promise in providing a safe, cost-effective, and synergistically multi-targeted treatment approach across various cancers. Herein, we hypothesized that the combined treatment with ABE and OPCs may synergistically regulate multiple oncogenic pathways, thereby presenting a more potent anti-cancer activity, particularly targeting CRC. Methods: We performed a series of cell culture experiments to evaluate the synergistic anticancer effect of the combination of ABE and OPCs in CRC cell lines. By performing genome-wide transcriptomic profiling, we identified specific growth-signaling pathways targeted by the combined treatment and subsequently validated these findings in cell lines and patient-derived 3D organoids. Results: The combined treatment with ABE and OPCs exhibited superior synergistic anti-tumor potential, affecting cell viability, proliferation, migration, and invasion in CRC cells. The transcriptomic analysis identified that LMNB1, associated with the cell apoptosis pathway Fold enrichment [FC]: 4.29, P =0.005), plays a crucial role in the anti-tumorigenic effects of these two natural products. Moreover, the combination of ABE and OPCs significantly altered the expression of key proteins involved in cell apoptosis and suppressed the expression level of LMNB1 in HCT116 (Combination vs. ABE, FC= 0.28; Combination vs. OPCs, FC= 0.24) and SW480 (Combination vs. ABE, FC= 0.65; Combination vs. OPCs, FC= 0.47), thereby inhibiting downstream Akt phosphorylation. Finally, we successfully demonstrated that the combined treatment displayed a superior anti-cancer activity in 3D organoids, significantly reducing both their number and size (P < 0.05). Conclusion: We present novel evidence for the synergistic anticancer effects of ABE and OPCs in CRC cells, which are in part achieved by regulating cellular apoptosis and the oncogene LMNB1 within the Akt signaling pathway. Citation Format: Yuan Li, Caiming Xu, Ajay Goel. Aronia berry and oligomeric proanthocyanidins synergistically exhibit anticancer effects in colorectal cancer by inhibiting LMNB1 [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 523.

Genome-wide transcriptome profiling and development of age prediction models in the human brain.

Aging-related transcriptome changes in various regions of the healthy human brain have been explored in previous works, however, a study to develop prediction models for age based on the expression levels of specific panels of transcripts is lacking. Moreover, studies that have assessed sexually dimorphic gene activities in the aging brain have reported discrepant results, suggesting that additional studies would be advantageous. The prefrontal cortex (PFC) region was previously shown to have a particularly large number of significant transcriptome alterations during healthy aging in a study that compared different regions in the human brain. We harmonized neuropathologically normal PFC transcriptome datasets obtained from the Gene Expression Omnibus (GEO) repository, ranging in age from 21 to 105 years, and found a large number of differentially regulated transcripts in the old and elderly, compared to young samples overall, and compared female and male-specific expression alterations. We assessed the genes that were associated with age by employing ontology, pathway, and network analyses. Furthermore, we applied various established (least absolute shrinkage and selection operator (Lasso) and Elastic Net (EN)) and recent (eXtreme Gradient Boosting (XGBoost) and Light Gradient Boosting Machine (LightGBM)) machine learning algorithms to develop accurate prediction models for chronological age and validated them. Studies to further validate these models in other large populations and molecular studies to elucidate the potential mechanisms by which the transcripts identified may be related to aging phenotypes would be advantageous.

Genome-wide identification and transcriptome profiling expression analysis of the U-box E3 ubiquitin ligase gene family related to abiotic stress in maize (Zea mays L.)

BackgroundThe U-box gene family encodes E3 ubiquitin ligases involved in plant hormone signaling pathways and abiotic stress responses. However, there has yet to be a comprehensive analysis of the U-box gene family in maize (Zea mays L.) and its responses to abiotic stress.ResultsIn this study, 85 U-box family proteins were identified in maize and were classified into four subfamilies based on phylogenetic analysis. In addition to the conserved U-box domain, we identified additional functional domains, including Pkinase, ARM, KAP and Tyr domains, by analyzing the conserved motifs and gene structures. Chromosomal localization and collinearity analysis revealed that gene duplications may have contributed to the expansion and evolution of the U-box gene family. GO annotation and KEGG pathway enrichment analysis identified a total of 105 GO terms and 21 KEGG pathways that were notably enriched, including ubiquitin-protein transferase activity, ubiquitin conjugating enzyme activity and ubiquitin-mediated proteolysis pathway. Tissue expression analysis showed that some ZmPUB genes were specifically expressed in certain tissues and that this could be due to their functions. In addition, RNA-seq data for maize seedlings under salt stress revealed 16 stress-inducible plant U-box genes, of which 10 genes were upregulated and 6 genes were downregulated. The qRT-PCR results for genes responding to abiotic stress were consistent with the transcriptome analysis. Among them, ZmPUB13, ZmPUB18, ZmPUB19 and ZmPUB68 were upregulated under all three abiotic stress conditions. Subcellular localization analysis showed that ZmPUB19 and ZmPUB59 were located in the nucleus.ConclusionsOverall, our study provides a comprehensive analysis of the U-box gene family in maize and its responses to abiotic stress, suggesting that U-box genes play an important role in the stress response and providing insights into the regulatory mechanisms underlying the response to abiotic stress in maize.

Genome-wide transcriptome profiling reveals molecular response pathways of Trichoderma harzianum in response to salt stress.

Trichoderma harzianum exhibits a strong biological control effect on many important plant pathogens, such as Fusarium oxysporum, Botrytis cinerea, and Meloidogyne. However, its biocontrol effectiveness is weakened or reduced under salt stress. The aim of this study was to investigate the molecular response of T. harzianum to salt stress at the whole-genome level. Here, we present a 44.47 Mb near-complete genome assembly of the T. harzianum qt40003 strain for the first time, which was assembled de novo with 7.59 Gb Nanopore sequencing long reads (~170-fold) and 5.2 Gb Illumina short reads (~116-fold). The assembled qt40003 genome contains 12 contigs, with a contig N50 of 4.81 Mb, in which four of the 12 contigs were entirely reconstructed in a single chromosome from telomere to telomere. The qt40003 genome contains 4.27 Mb of repeat sequences and 12,238 protein-coding genes with a BUSCO completeness of 97.5%, indicating the high accuracy and completeness of our gene annotations. Genome-wide transcriptomic analysis was used to investigate gene expression changes related to salt stress in qt40003 at 0, 2% (T2), and 4% (T4) sodium chloride concentrations. A total of 2,937 and 3,527 differentially expressed genes (DEGs) were obtained under T2 and T4 conditions, respectively. GO enrichment analysis showed that the T2-treatment DEGs were highly enriched in detoxification (p < 0.001), while the T4 DEGs were mainly enriched in cell components, mostly in cellular detoxification, cell surface, and cell wall. KEGG metabolic pathway analysis showed that 91 and 173 DEGs were significantly enriched in the T2 and T4 treatments, respectively (p < 0.01), mainly in the glutathione metabolism pathway. We further experimentally analyzed the differentially expressed glutathione transferase genes in the glutathione metabolic pathway, most of which were downregulated (13/15). In addition, we screened 13 genes related to active oxygen clearance, including six upregulated and seven downregulated genes, alongside five fungal hydrophobic proteins, of which two genes were highly expressed. Our study provides high-quality genome information for the use of T. harzianum for biological control and offers significant insights into the molecular responses of T. harzianum under salt-stress conditions.

The SlGRAS9-SlZHD17 transcriptional cascade regulates chlorophyll and carbohydrate metabolism contributing to fruit quality traits in tomato.

Some essential components of fleshy fruits are dependent on photosynthetic activity and carbohydrate metabolism. Nevertheless, the regulatory mechanisms linking chlorophyll and carbohydrate metabolism remain partially understood. Here, we uncovered the role of SlGRAS9 and SlZHD17 transcription factors in controlling chlorophyll and carbohydrate accumulation in tomato fruit. Knockout or knockdown of SlGRAS9 or SlZHD17 resulted in marked increase in chlorophyll content, reprogrammed chloroplast biogenesis and enhanced accumulation of starch and soluble sugars. Combined genome-wide transcriptomic profiling and promoter-binding experiments unveiled a complex mechanism in which the SlGRAS9/SlZHD17 regulatory module modulates the expression of chloroplast and sugar metabolism either via a sequential transcriptional cascade or through binding of both TFs to the same gene promoters, or, alternatively, via parallel pathways where each of the TFs act on different target genes. For instance, the regulation of SlAGPaseS1 and SlSUS1 is mediated by SlZHD17 whereas that of SlVI and SlGLK1 occurs only through SlGRAS9 without the intervention of SlZHD17. Both SlGRAS9 and SlZHD17 can also directly bind the promoter of SlPOR-B to regulate its expression. Taken together, our findings uncover two important regulators acting synergistically to manipulate chlorophyll and carbohydrate accumulation and provide new potential breeding targets for improving fruit quality in fleshy fruits.

Epigenetic Regulation of Cardiomyocyte Maturation by Arginine Methyltransferase CARM1.

During the neonatal stage, the cardiomyocyte undergoes a constellation of molecular, cytoarchitectural, and functional changes known collectively as cardiomyocyte maturation to increase myocardial contractility and cardiac output. Despite the importance of cardiomyocyte maturation, the molecular mechanisms governing this critical process remain largely unexplored. We leveraged an in vivo mosaic knockout system to characterize the role of Carm1, the founding member of protein arginine methyltransferase, in cardiomyocyte maturation. Using a battery of assays, including immunohistochemistry, immuno-electron microscopy imaging, and action potential recording, we assessed the effect of loss of Carm1 function on cardiomyocyte cell growth, myofibril expansion, T-tubule formation, and electrophysiological maturation. Genome-wide transcriptome profiling, H3R17me2a chromatin immunoprecipitation followed by sequencing, and assay for transposase-accessible chromatin with high-throughput sequencing were used to investigate the mechanisms by which CARM1 (coactivator-associated arginine methyltransferase 1) regulates cardiomyocyte maturation. Finally, we interrogated the human syntenic region to the H3R17me2a chromatin immunoprecipitation followed by sequencing peaks for single-nucleotide polymorphisms associated with human heart diseases. We report that mosaic ablation of Carm1 disrupts multiple aspects of cardiomyocyte maturation cell autonomously, leading to reduced cardiomyocyte size and sarcomere thickness, severe loss and disorganization of T tubules, and compromised electrophysiological maturation. Genomics study demonstrates that CARM1 directly activates genes that underlie cardiomyocyte cytoarchitectural and electrophysiological maturation. Moreover, our study reveals significant enrichment of human heart disease-associated single-nucleotide polymorphisms in the human genomic region syntenic to the H3R17me2a chromatin immunoprecipitation followed by sequencing peaks. This study establishes a critical and multifaceted role for CARM1 in regulating cardiomyocyte maturation and demonstrates that deregulation of CARM1-dependent cardiomyocyte maturation gene expression may contribute to human heart diseases.

Transcriptomic insights into the molecular mechanism for response of wild emmer wheat to stripe rust fungus.

Continuous identification and application of novel resistance genes against stripe rust are of great importance for wheat breeding. Wild emmer wheat, Triticum dicoccoides, has adapted to a broad range of environments and is a valuable genetic resource that harbors important beneficial traits, including resistance to stripe rust caused by Puccinia striiformis f. sp. tritici (Pst). However, there has been a lack of systematic exploration of genes against Pst races in wild emmer wheat. Genome-wide transcriptome profiles were conducted on two wild emmer wheat genotypes with different levels of resistance to (Pst (DR3 exhibiting moderate (Pst resistance, and D7 displaying high (Pst resistance). qRT-PCR was performed to verify findings by RNA-seq. A higher number of DEGs were identified in the moderately (Pst-resistant genotype, while the highly (Pst-resistant genotype exhibited a greater enrichment of pathways. Nonetheless, there were consistent patterns in the enrichment of pathways between the two genotypes at the same time of inoculation. At 24 hpi, a majority of pathways such as the biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, phenylalanine metabolism, and alpha-Linolenic acid metabolism exhibited significant enrichment in both genotypes. At 72 hpi, the biosynthesis of secondary metabolites and circadian rhythm-plant pathways were notably and consistently enriched in both genotypes. The majority of (WRKY, MADs , and AP2-ERF families were found to be involved in the initial stage of response to Pst invasion (24 hpi), while the MYB, NAC, TCP, and b-ZIP families played a role in defense during the later stage of Pst infection (72 hpi). In this present study, we identified numerous crucial genes, transcription factors, and pathways associated with the response and regulation of wild emmer wheat to Pst infection. Our findings offer valuable information for understanding the function of crucial Pst-responsive genes, and will deepen the understanding of the complex resistance mechanisms against Pst in wheat.

Genome-wide transcriptome profiling of mulberry (Morus alba) response to boron deficiency and toxicity reveal candidate genes associated with boron tolerance in leaves

Mulberry (Morus alba) is an essential plant with countless economic benefits; however, its growth and metabolic processes are hampered by boron (B) stresses. Very little research has been performed to elucidate boron tolerance and detoxification mechanisms in this species. The M. alba cultivar, Yu-711, was exposed to five different concentrations of boric acid (H3BO3), including deficient (T1; 0 mM) moderate B deficiency (T2; 0.02 mM), sufficient (CK; 0.1 mM) and toxic (T3 and T4; 0.5 and 1 mM) levels for 18 days of growth in pots experiment. Transcriptome analysis of B deficiency and toxicity treatments was performed on mulberry leaves. The transcriptome data reveal that a total of 6114 genes were differentially expressed (DEGs), of which 3830 were up-regulated and 2284 were down-regulated. A comparative analysis between treatment groups CK-vs-T1 (deficiency) and CK-vs-T4 (toxicity) indicates that 590 and 1383 genes were down-regulated in both deficiency and B toxicity, respectively. The results show that 206 genes were differentially expressed in all treatments. B deficiency and toxicity significantly altered the expression of the key aquaporins (PIP2-1, PIP2-7, PIP2-4 and NIP3-1) and high-affinity boron transporter genes (BOR1 and BOR7). In addition, boron stress also altered the expression of antioxidants and photosynthesis-related genes. B stresses were found to alter several transcription factors including ERF1B, which is associated with the regulation of boron uptake and the synthesis and signaling of phytohormones. Unravelling the mechanisms of B tolerance and detoxification is important and would give us further insight into how B stresses affect mulberry plants.

Genome-wide transcriptomic and biochemical profiling of major depressive disorder: Unravelling association with susceptibility, severity, and antidepressant response

Identifying biomarkers for diagnosing Major Depressive Disorder (MDD), assessing its severity, and guiding treatment is crucial. We conducted whole genome transcriptomic study in North Indian population, and analyzed biochemical parameters. Our longitudinal study investigated gene-expression profiles from 72 drug-free MDD patients and 50 healthy controls(HCs) at baseline and 24 patients after 12-weeks of treatment. Gene expression analyses identified differentially expressed genes(DEGs) associated with MDD susceptibility, symptom severity and treatment response, independently validated by qPCR. Hierarchical clustering revealed distinct expression patterns between MDD and HCs, also between mild and severe cases. Enrichment analyses of significant DEGs revealed inflammatory, apoptosis, and immune-related pathways in MDD susceptibility, severity, and treatment response. Simultaneously, we assessed thirty biochemical parameters in the same cohort, showed significant differences between MDD and HCs in 13 parameters with monocytes, eosinophils, creatinine, SGPT, and total protein remained independent predictors of MDD in a multivariate-regression model. Our study supports the role of altered immune/inflammatory signaling in MDD pathophysiology, offering clinically relevant biochemical parameters and insights into transcriptomic gene regulation in MDD pathogenesis and treatment response.

Genome-Wide Methylome and Transcriptome Profiling of Acute Myeloid Leukemia Derived Bone Marrow Mesenchymal Cells Identify Age Group Specific Biological Pathway Dysregulation in the Bone Marrow Microenvironment of AML

Adult patients with acute myeloid leukemia (AML) have a proportionate increase in mortality from their disease and therapy as age increases. In great part, this mortality is due to disease relapse that occurs at a higher rate irrespective of genetic assigned risk. While intensification of therapy in patients with favorable genomics such as core-binding-factor AML (CBF-AML) has yielded a high proportion of long-term remission and potential cure, this benefit has been less frequently observed in patients with AML age 60 and older despite receipt of similar therapy. Moreover, even less benefit has been observed in patients with unfavorable genomics such when bearing a complex karyotype (karyotype with ≥ 3 chromosome abnormalities). Reasons for this unfavorable outcome among older patients is unknown. We hypothesized that tumor-induced aging in cells of bone marrow microenvironment, specifically, in bone marrow mesenchymal stromal cells (BMSCs) plays a significant role in AML survival and progression in older compared to younger patients and that the difference in prognosis may depend on specific aging induced changes in the BMSCs. We expanded BMSCs under physiologic hypoxic conditions (91% nitrogen, 4% oxygen and 5% CO 2) to mimic the bone marrow microenvironment, isolated DNA and RNA from early passaged cells and performed methylation analysis using 850k Illumina DNAmehtylation arrays and gene expression analysis using mRNAseq. We used the Horvath Skin and Blood clock to calculate epigenetic age (DNA methylation age). We determined differentially expressed genes with increasing epigenetic or chronological age in the different AML genetic risk groups (n=9 CBF-AML, n=11 CK-AML and n=7 Normal donor BMSC) as well as in older (&amp;gt;60) and younger (&amp;lt;60) AML-BMSCs, and performed Gene Set Enrichment Analysis (GSEA) for the Hallmarks as well as aging and senescence related genesets (n=33 genesets). We found that BMSCs derived from AML bone marrow cells are more age accelerated than normal donor BMSCs. Our GSEA analysis revealed enrichment for age-associated gene networks that reflected hallmarks of aging, including mitochondrial dysfunction, inflammation, and cellular senescence. Specifically, we observed a positive enrichment for the SenMayo gene signature ( Saul et al. 2022, Nature Communications) and a negative enrichment of proliferation and mitotic gene sets with increasing epigenetic and chronological age (padj&amp;lt;0.05). Gene sets that are related to inhibition of adipogenesis were positively enriched in the CBF-AML BMSCs compared to the CK-AML BMSCs. Our integrated analysis revealed novel observations regarding the cellular changes that occur in the AML stromal microenvironment with increasing age. Most interestingly, our findings suggest that AML derived BMSCs have increased epigenetic age compared to normal donor BMSCs. Additional experiments are ongoing to functionally validate the findings predicted by the enrichment analysis. Our study is the first to suggest an increased epigenetic age in the stromal bone marrow microenvironment of older AML patients.

Polycystin-2-dependent transcriptome reveals early response of autosomal dominant polycystic kidney disease.

Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in polycystin genes, Pkd1 and Pkd2, but the underlying pathogenic mechanisms are poorly understood. To identify genes and pathways that operate downstream of polycystin-2 (PC2), a comprehensive gene expression database was created, cataloging changes in the transcriptome immediately following PC2 protein depletion. To explore cyst initiation processes, an immortalized mouse inner medullary collecting duct line was developed with the ability to knock out the Pkd2 gene conditionally. Genome-wide transcriptome profiling was performed using RNA sequencing in the cells immediately after PC2 was depleted and compared with isogenic control cells. Differentially expressed genes were identified, and a bioinformatic analysis pipeline was implemented. Altered expression of candidate cystogenic genes was validated in Pkd2 knockout mice. The expression of nearly 900 genes changed upon PC2 depletion. Differentially expressed genes were enriched for genes encoding components of the primary cilia, the canonical Wnt pathway, and MAPK signaling. Among the PC2-dependent ciliary genes, the transcription factor Glis3 was significantly downregulated. MAPK signaling formed a key node at the epicenter of PC2-dependent signaling networks. Activation of Wnt and MAPK signaling, concomitant with the downregulation of Glis3, was corroborated in Pkd2 knockout mice. The data identify a PC2 cilia-to-nucleus signaling axis and dysregulation of the Gli-similar subfamily of transcription factors as a potential initiator of cyst formation in ADPKD. The catalog of PC2-regulated genes should provide a valuable resource for future ADPKD research and new opportunities for drug development.NEW & NOTEWORTHY Autosomal dominant polycystic kidney disease (ADPKD) is the most common inherited kidney disease. Mutations in polycystin genes cause the disease, but the underlying mechanisms of cystogenesis are unknown. To help fill this knowledge gap, we created an inducible cell model of ADPKD and assembled a catalog of genes that respond in immediate proximity to polycystin-2 depletion using transcriptomic profiling. The catalog unveils a ciliary signaling-to-nucleus axis proximal to polycystin-2 dysfunction, highlighting Glis, Wnt, and MAPK signaling.

Comparing P53 expression and genome-wide transcriptome profiling to Comet assay in lymphocytes from melanoma patients and healthy controls.

This study compared the expression of TP53 in lymphocytes from malignant melanoma (MM) patients with positive sentinel nodes to healthy controls (HCs) following exposure to various doses of UVA radiation. The Lymphocyte Genome Sensitivity (LGS) assay indicated significant differences in DNA damage in lymphocytes between MM patients and HCs. qPCR data demonstrated an overall 3.4-fold increase in TP53 expression in lymphocytes from MM patients compared to healthy controls, following treatment with 0.5 mW/cm2 UVA radiation. Western blotting confirmed that p53 expression was increased in MM lymphocytes following UVA exposure compared to healthy individuals. Genome transcriptome profiling data displayed differences in gene expression between UVA-treated lymphocytes from MM patients and HCs. Peripheral lymphocytes from MM patients are more susceptible to the genotoxic effects of UVA compared to healthy individuals. Our previous studies showed that UVA exposure of various intensities caused significant differences in the levels of DNA damage between lymphocytes from cancer patients compared to HCs through the LGS assay. The present study's results provide further credibility to the LGS assay as a screening test for cancer detection. Peripheral lymphocytes could be a promising blood biopsy biomarker for staging of carcinomas and prevention of carcinoma progression at early stages.

Methamphetamine exposure drives cell cycle exit and aberrant differentiation in rat hippocampal-derived neurospheres.

Introduction: Methamphetamine (METH) abuse by pregnant drug addicts causes toxic effects on fetal neurodevelopment; however, the mechanism underlying such effect of METH is poorly understood. Methods: In the present study, we applied three-dimensional (3D) neurospheres derived from the embryonic rat hippocampal tissue to investigate the effect of METH on neurodevelopment. Through the combination of whole genome transcriptional analyses, the involved cell signalings were identified and investigated. Results: We found that METH treatment for 24h significantly and concentration-dependently reduced the size of neurospheres. Analyses of genome-wide transcriptomic profiles found that those down-regulated differentially expressed genes (DEGs) upon METH exposure were remarkably enriched in the cell cycle progression. By measuring the cell cycle and the expression of cell cycle-related checkpoint proteins, we found that METH exposure significantly elevated the percentage of G0/G1 phase and decreased the levels of the proteins involved in the G1/S transition, indicating G0/G1 cell cycle arrest. Furthermore, during the early neurodevelopment stage of neurospheres, METH caused aberrant cell differentiation both in the neurons and astrocytes, and attenuated migration ability of neurospheres accompanied by increased oxidative stress and apoptosis. Conclusion: Our findings reveal that METH induces an aberrant cell cycle arrest and neuronal differentiation, impairing the coordination of migration and differentiation of neurospheres.

LncRNA expression signature identified using genome-wide transcriptomic profiling to predict lymph node metastasis in patients with stage T1 and T2 gastric cancer

BackgroundLymph node (LN) status is vital to evaluate the curative potential of relatively early gastric cancer (GC; T1–T2) treatment (endoscopic or surgery). Currently, there is a lack of robust and convenient methods to identify LN metastasis before therapeutic decision-making.MethodsGenome-wide expression profiles of long noncoding RNA (lncRNA) in primary T1 gastric cancer data from The Cancer Genome Atlas (TCGA) was used to identify lncRNA expression signature capable of detecting LN metastasis of GC and establish a 10-lncRNA risk-prediction model based on deep learning. The performance of the lncRNA panel in diagnosing LN metastasis was evaluated both in silico and clinical validation methods. In silico validation was conducted using TCGA and Asian Cancer Research Group (ACRG) datasets. Clinical validation was performed on T1 and T2 patients, and the panel’s efficacy was compared with that of traditional tumor markers and computed tomography (CT) scans.ResultsProfiling of genome-wide RNA expression identified a panel of lncRNA to predict LN metastasis in T1 stage gastric cancer (AUC = 0.961). A 10-lncRNA risk-prediction model was then constructed, which was validated successfully in T1 and T2 datasets (TCGA, AUC = 0.852; ACRG, AUC = 0.834). Thereafter, the clinical performance of the lncRNA panel was validated in clinical cohorts (T1, AUC = 0.812; T2, AUC = 0.805; T1 + T2, AUC = 0.764). Notably, the panel demonstrated significantly better performance compared with CT and traditional tumor markers. ConclusionsThe novel 10-lncRNA could diagnose LN metastasis robustly in relatively early gastric cancer (T1–T2), with promising clinical potential.