Digital Gene Expression Profiling Research Articles

Heterogeneity of hepatocellular carcinoma that responds differently to combination therapy with TACE and Sorafenib as determined by digital spatial gene expression profiling.

The combination of Sorafenib and transcatheter arterial chemoembolization (TACE) exhibits limited efficacy in the treatment of certain advanced hepatocellular carcinomas (HCC), and the molecular mechanisms underlying resistance to this combination remain unclear. This study aims to underscore the distinctive contribution of GeoMx DSP technology in elucidating the molecular intricacies of HCC resistance to the Sorafenib and TACE combination. Patients with advanced HCC during the waiting period before liver transplantation were classified into sensitive and resistant groups based on their response to Sorafenib and TACE combination therapy. Employing GeoMx DSP technology for comprehensive gene expression profiling, we identified pivotal molecular targets linked to resistance against combination therapy. The investigation scrutinized intra-tumoral and inter-individual variances, unveiling a spectrum of crucial molecular targets, such as PLG, PLVAP, immunoglobulin genes, ORM1, and NR4A1, among others. Additionally, we explored signaling pathways associated with treatment responsiveness, including the PPAR signaling pathway. Notably, we emphasized the significance of the immune microenvironment characterized by heightened SPP1 expression in HCC resistance to combination therapy. In the resistant group, SPP1+ tumor-associated macrophage (TAM) infiltration was notably pronounced (p = 0.037), while T-cell depletion showed a mitigated presence (p = 0.013). The study reveals intra- and inter-individual heterogeneity in HCC that is differentially responsive to the combination of Sorafenib and TACE, highlighting multiple key molecular targets associated with treatment resistance. The immune microenvironment is important, and in particular, SPP1+ TAM infiltration may play a key role. Meanwhile, the introduction of immunotherapy in patients resistant to combination therapy may lead to positive results.

Abstract PO3-14-01: Heterogeneity of tumor immune microenvironment to predict everolimus efficacy in premenopausal women with hormone receptor-positive/human epidermal growth factor receptor 2-negative adavanced breast cancer

Abstract Objective Everolimus (EVE), a mammalian target of rapamycin (mTOR) inhibitor, has been shown to prolong survival of patients with hormone receptor-positive/human epidermal growth factor receptor 2-negative (HR+/HER2-), advanced breast cancer (ABC) with endocrine resistance. Preclinical studies have suggested that mTOR is closely associated with the tumor immune microenvironment (TME). We aimed to explore the tumor heterogeneity of TME between patients with different response to everolimus (EVE), and identify potential markers to predict efficacy of EVE in premenopausal patients with HR+/HER2−, advanced breast cancer (ABC) suffering endocrine resistance. Methods Patients were divided into two groups according to the best response to EVE: the responsive subgroup (RS) and the non-responsive subgroup (NRS). Postoperative tumor tissue from patients who received EVE was collected for whole exome sequencing. Transcript abundance of 126 distinct tumor regions of interest (ROIs) were quantitated using digital spatial profiling, including 42 tumor cell zones (TCZs), 42 immune cell zones (ICZs), and 42 stroma cell zones (SCZs). Based on investigation in multiple discrete areas, differences in cell components, function and pathway, gene expression and immune cell infiltration between RS and NRS were investigated. Results Totally, 93 differentially expressed genes (DEGs) were identified in the ICZs, 174 DEGs in the TCZs, and 197 DEGs in the SCZs, using p< 0.05 and |log2FC >1| as the screening criterion. Among all the subgroup analyses on different ROIs, the statistical significance of PIP gene was the highest (ICZs: p=2.6E-12, TCZs: p=1.25E-20, SCZs: p=2.0E-18) in the comparison of RS and NRS. Moreover, we found that PIP was highly expressed in NRS and patients with a progression-free survival < 1 year. The receptor-legend analysis showed that, the intra-and inter-cellular communication in RS was mainly mediated by interaction between CCL-CCR family, while NRS mainly by COL family complex. Comparing all cell components extracted from each ROI, no statistical differences were found between RS and NRS, only a statistical trend in the SCZs (p=0.098). Further cell classification analysis showed fibroblasts highest in TCZs and SCZs, macrophages in ICZs. Compared to NRS, in TCZs, fibroblasts were significantly increased (p=0.01), and plasma cells (p=0.02) and neutrophils (p=0.03) were decreased in RS. In SCZs, the downregulated fibroblasts (p=0.04) and enhanced plasma cells (p=0.01) were observed in NRS. In ICZs, pDCs (p=0.048) and neutrophils (p=0.006) were significantly elevated in NRS compared with RS. Conclusion The utility of digital spatial gene expression profiling might provide some references to predict the efficacy of EVE in premenopausal patients with HR+/HER2- ABC with endocrine resistance. Citation Format: Yujing Tan, Fei Ma, Jiayu Wang, Pin Zhang, Binghe Xu, Liyan Xue, Ying Fan. Heterogeneity of tumor immune microenvironment to predict everolimus efficacy in premenopausal women with hormone receptor-positive/human epidermal growth factor receptor 2-negative adavanced breast cancer [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO3-14-01.

Transcriptome response of diatom Skeletonema marinoi to lower temperature

Temperature is an important environmental factor for phytoplankton. Phytoplankton growth, metabolism, biodiversity, productivity, and distribution are influenced by temperature-driven nutrient stratification and mixing, as well as species’ optimal growth temperatures. There have been a number of studies focused on physiological and biochemical mechanisms of environmental–biological interactions in diatoms, yet the underlying transcriptional regulators remain limited. Here, we performed an RNA-seq-based gene expression analysis to explore the Skeletonema marinoi (isolated from Jiaozhou Bay of Qingdao, 36.13°N, 120.16°E on July 5th, 2013) cellular responses induced by low temperature (12 °C). Digital gene expression profiling of S. marinoi generated 20,319 unigenes, of which 573 differentially expressed genes appeared in the low-temperature treatment group. According to GO and KEGG enrichment analysis, different genes were involved in ten metabolic and biosynthesis pathways: ribosome, lipid, porphyrin, and chlorophyll metabolism showed strong transcriptional cold tolerances. The regulation of genes related to translation processes (e.g., pentatricopeptide repeat), fatty acid metabolism (e.g., acyl-CoA synthetase), and photosynthesis (e.g., porphyrin enzymes) provides new molecular-level insight into cold stress responses in eukaryotic marine phytoplankton. Our study suggests that this Skeletonema species could be a potential candidate for understanding the fate of thermo-sensitive diatom communities and oceanic ecosystems facing climate change.

Roles of the APETALA3-3 ortholog in the petal identity specification and morphological differentiation in Delphinium anthriscifolium flowers.

The genus Delphinium (Ranunculaceae) with its unique and highly complex floral structure is an ideal system to address some key questions in terms of morphological and evolutionary studies in flowers. In D. anthriscifolium, for example, the original eight petal primordia differentiate into three types at maturity (i.e., two dorsal spurred, two lateral flat, and four ventral reduced petals). The mechanisms underlying their identity determination and morphological differentiation remain unclear. Here, through a comprehensive approach combining digital gene expression (DGE) profiles, in situ hybridization, and virus-induced gene silencing (VIGS), we explore the role of the APETALLATA3-3 (AP3-3) ortholog in D. anthriscifolium. Our findings reveal that the DeanAP3-3 not only functions as a traditionally known petal identity gene but also plays a critical role in petal morphological differentiation. The DeanAP3-3 gene is expressed in all the petal primordia before their morphological differentiation at earlier stages, but shows a gradient expression level difference along the dorsventral floral axis, with higher expression level in the dorsal spurred petals, intermediate level in the lateral flat petals and lower level in the ventral reduced petals. VIGS experiments revealed that flowers with strong phenotypic changes showed a complete transformation of all the three types of petals into non-spurred sepals. However, in the flowers with moderate phenotypic changes, the transformation of spurred petals into flat petals is associated with moderate silencing of the DeanAP3-3 gene, suggesting a significant impact of expression level on petal morphological differentiation. This research also shed some insights into the role of changes in gene expression levels on morphological differentiation in plants.

Refining Diagnostic Subtypes of Peripheral T-Cell Lymphoma Using a Multiparameter Approach

Peripheral T-cell lymphoma (PTCL) encompasses a diverse group of post-thymic lymphomas, with ~40% of cases not further classifiable and designated as PTCL-not otherwise specified (PTCL-NOS). Two molecular prognostic subtypes of PTCL-NOS, PTCL-TBX21 and PTCL-GATA3, were identified through gene expression profiling (GEP). A subset of PTCL-NOS with T-follicular helper (TFH) differentiation was subsequently categorized as nodal PTCL with TFH phenotype (PTCL-TFH), also known in current classifications as nodal PTCL-TFH, NOS, or TFH lymphoma, NOS. However, the boundary between PTCL-NOS and PTCL-TFH remains poorly defined. To refine these subtypes, 101 PTCL-NOS with centralized pathology review and extensive immunophenotyping by immunohistochemistry (IHC) using 22 antibodies, digital GEP (nCounter, NanoString Inc) (n= 70), and RNA sequencing (n= 73) were included in this study. IHC was used to identify PTCL-TFH cases (those with strong expression of more than two TFH markers); the remaining patients were classified as PTCL-GATA3 and PTCL-TBX21 using pre-defined GEP signatures. Cases were reclassified into PTCL-NOS (n= 63) and PTCL-TFH (n= 38), with PTCL-NOS subclassified into PTCL-GATA3 (n= 22; 34%) and PTCL-TBX21 (n= 41; 66%) and showing significant differences in OS (p<0.001). PTCL-GATA3 was characterized by medium to large transformed cells (average= 90%, 70-100%) and had minimal tumor microenvironment (TME) (TME-poor: 100%). In contrast, PTCL-TBX21 was more heterogeneous, associated with pleomorphic cells in a polymorphous background (TME-rich: 78%), including Lennert lymphoma-like cases (22%). mRNA and CIBERSORT analysis substantiated the findings and identified a subset enriched in cytotoxic features in the PTCL-TBX21 subtype. IHC indicated that PTCL-GATA3 cases were CD4+/CD8- (83%) or CD4-/CD8- (17%) and lacked expression of CD8 or cytotoxic markers compared to PTCL-TBX21 (p<0.01). PTCL-GATA3 showed significantly higher expression of LEF1 (average= 80%), Ki67 (80%), and MYC (25%) than PTCL-TBX21 (25%, 30%, <5%, respectively, p<0.01). mRNA and protein expression of these biomarkers showed a significant positive correlation (r=0.5, p<0.001), and expectedly higher mRNA expression of LEF1 and MYC was observed in the PTCL-GATA3 versus PTCL-TBX21 (p<0.05). Strong expression of CD30 (>50% of cells) was only seen in PTCL-GATA3 cases. EBER positivity, found only in rare background cells, was not significantly different (18% of PTCL-GATA3 vs. 12% of PTCL-TBX21). Within PTCL-TBX21, we identified cytotoxic and non-cytotoxic subsets with divergent morphological, phenotypic, and clinical findings. The cytotoxic PTCL-TBX21 exhibited an activated cytotoxic phenotype (23/25, 96%), denoted by TIA1 and granzyme-B and/or perforin expression. Extranodal involvement and single-cell apoptosis were observed in 17% and 45% of the cytotoxic cases, respectively, but absent in non-cytotoxic PTCL-TBX21 cases. A trend towards higher Ki67 expression (average= 40% vs. 20%, p= 0.08) was seen in the cytotoxic subgroup. In contrast, the non-cytotoxic PTCL-TBX21 was associated with a CD4+/CD8- phenotype and higher ICOS expression (average= 30%) and CCR4 (60%) compared to the cytotoxic PTCL-TBX21 (p= 0.001). PTCL-TFH cases showed a CD4+/CD8- phenotype (90%) and rarely a CD4-/CD8- phenotype (10%). While a subset of PTCL-TFH cases had AITL-like features such as numerous clear cells and/or prominent vasculature (31% of cases), which were not seen in the other subgroups, the remaining cases exhibited morphology that was indistinguishable from other PTLC-NOS cases including sheets of transformed cells or pleomorphic cells in a polymorphous background. Morphologically, PTCL-TFH cases, like PTCL-TBX21 were associated with a rich TME (TME-rich: 75%). CIBERSORT analysis showed enrichment of plasma cells (p<0.01) in PTCL-TFH, compared to PTCL-TBX21, and validated by morphology (30% of cases vs. 5%). Expression of one TFH marker was frequent in some PTCL-NOS cases (PTCL-GATA3= 41% of cases, PTCL-TBX21-non cytotoxic= 47%, PTCL-TBX21-cytotoxic= 5%). Conclusion: Our comprehensive evaluation underscores the importance of integrating morphology, immunophenotyping, and GEP in achieving an accurate diagnosis, potentially leading to more tailored treatment strategies for PTCL. Correlation with pending whole-exome sequencing studies will be provided at the time of presentation.

Genomic Stability, DNA Repair and the SKP1-CUL1-F-Box Complex in HIV Associated Diffuse Large B-Cell Lymphoma

Introduction: Previously we demonstrated that HIV associated diffuse large B-cell lymphomas (DLBCL) of Germinal Center B-cell subtype (GCB) are more proliferative, with enhanced genomic stability and increased expression of DNA repair genes, particularly Fanconi anemia (FA) and homologous recombination (HR) genes compared to their HIV negative counterparts (Maguire et al, 2019). It is unknown if this signature is unique to the GCB subtype or if it is also observed in Activated B-cell (ABC) and Unclassified (UNC) DLBCL. Herein we validate our HIV associated GCB DLBCL findings in a larger independent cohort and report for the first time the aberrational profile of molecularly profiled HIV associated ABC DLBCL compared to their HIV negative counterparts respectively. Methods: An independent cohort of 353 DLBCL FFPE cases was assembled, including 125 HIV-DLBCL samples. All samples underwent hematopathological review to confirm diagnosis and determine tumor content. Samples with <60% tumor content were macro-dissected before they were RNA/DNA extracted. Cell of origin (COO) profiles were assessed using the Lymph3Cx molecular assay which can distinguish primary mediastinal B-cell lymphoma (PMBL) like from GCB, ABC and UNC DLBCL (Mottok et al 2018). Copy number (CN) analyses were performed on samples with sufficient DNA (41 HIV neg, 32 HIV pos) using the OncoScan FFPE assay kit and the resulting data analyzed using Chromosome Analysis Suite (ChAS) v4.1 and the multisample viewer module of Reproductive Health Research Analysis Software (RHAS) v1.1. Additional digital gene expression profiling was performed using the Pan-Cancer panel (NanoString, Seattle, WA) customized to include 9 BCL2 family genes for a total of 739 genes and analyzed using the NanoString nSolver Advanced Analysis module. Results: GCB, ABC, UNC and PMBL-like proportions were 50%, 30%, 14% and 6% in HIV unassociated (HIV neg) and 50%, 23%, 25% and 2% in HIV pos tumors respectively, confirming both cohorts were GCB-enriched by molecular profiling. PMBL-like cases were excluded from further analyses. OncoScan CN analyses showed that HIV pos DLBCL were significantly less aberrant than HIV neg DLBCL (median CN per sample 0.56 vs 1.01, p<0.001). HIV pos GCB tumors were significantly less aberrant than both HIV neg GCB tumors (median CN per sample 0.50 vs 1.11, p<0.001), and HIV pos ABC tumors (median CN per sample 0.50 vs 0.81, p=0.009). HIV pos ABC tumors were also less aberrant than HIV neg ABC tumors, however, the difference was not significant overall (median CN per sample 0.81 vs 1.03, p=0.162) or for gain/amplifications (median CN per sample 0.38 vs 0.42, p=0.814), but was borderline significant for loss/deletions (median CN per sample 0.31 vs 0.47, p=0.06). GEP analyses identified Cell Cycle and DNA repair as the top two upregulated pathways in HIV pos GCB tumors compared to HIV neg GCB tumors. In contrast, Cell Cycle and DNA repair were downregulated in HIV pos ABC tumors compared to HIV neg ABC tumors. However, one gene ( RBX1) was identified as significantly upregulated in both pathways in HIV pos GCB (adjusted p-value 0.0215) and HIV pos ABC (adjusted p-value 0.0046) tumors compared to their HIV neg counterparts. A second gene ( SKP1) was also significantly upregulated in HIV pos GCB (adjusted p-value 3.94E-5) and HIV pos ABC (adjusted p-value 0.0031) tumors. Together RBX1 and SKP1 encode the two principal functional subunits of the SKP1, CUL1, F-box protein (SCF) complex which is part of the ubiquitin proteasome system (UPS) and has been reported as essential for maintaining genome and chromosome stability (Thompson et al 2021). Conclusions: These results demonstrate that HIV pos DLBCL tumors are less aberrant than their HIV neg counterparts irrespective of COO subtype. However, the degree of genomic stability differs by COO, with HIV pos GCB tumors being more genomically stable than HIV pos ABC tumors. Although both have enhanced expression of UPS related SCF genes ( RBX1, SKP1), only HIV pos GCB tumors have enhanced expression of key DNA repair genes as well including Fanconi Anemia ( BRCA2/FANCD1, FANCL) and mismatch repair ( MSH2, MSH6) genes (adjusted p-values of ≤ 0.0263). Given the reported role of SCF in genomic stability maintenance, the HIV pos COO CN differences could be due to the contributions of both DNA repair and SCF mechanisms in HIV pos GCB and mainly SCF mechanisms in HIV pos ABC tumors.

Transcriptome and digital gene expression analysis reveal immune responses of mantle and visceral mass pearl culturing in Hyriopsis cumingii

Biomineralization is a widespread phenomenon in marine mollusks and is responsible for the production of shells and pearls. However, the regulatory mechanisms governing the adaptive immune responses in the mollusk mantle and visceral mass during mineralization remain unclear. In this work, we examined the mantle and visceral mass immune responses of Hyriopsis cumingii during pearl culture using high-throughput sequencing techniques. A mantle transcriptome database was established using transcriptome sequencing technology and reference to the major databases. Digital gene expression profiling was used to identify the differentially expressed genes of mantle and visceral mass at different insertion periods. Moreover, quantitative real-time PCR was used to verify the expression of five immune-related genes. Transcriptome sequencing results showed 257,457 unigenes were identified. Digital gene expression profiles showed 1389, 3572, 1888, and 2613 differentially expressed genes (DEGs) in the mantle and visceral mass at 5, 20, 50, and 90 d after insertion, respectively, with the highest number at 20 d and the lowest at 5 d after insertion (q < 0.05). A cluster analysis of the DEGs showed similar clustering and expression features in the mantle to the control group, and at 5, 50 and 90 d, after mantle insertion. The DEGs in the visceral mass showed similar clustering and expression features to the control group and at 5, 20 and 50 d after insertion. We also screened 22 immune-related DEGs in the mantle and visceral mass during the same pearl culture period, including serine/threonine-protein kinase NLK, C-type lectin, and galectin. The greatest number of DEGs was found 90 d after insertion. Compared with the mantle, more immune-related DEGs were down-regulated than up-regulated in the visceral mass during pearl culture, indicating that the immune regulatory mechanisms in the visceral mass and the mantle differ during pearl culture, and that the visceral mass is liable to higher infection and mortality rates. Quantitative real-time PCR results showed that the expression of five immune-related genes was consistent with DGE results. Our findings will further knowledge of the immune systems that are present in the mantle and visceral mass during pearl culture.

Long non-coding RNA mitophagy and ALK-negative anaplastic lymphoma-associated transcript: a novel regulator of mitophagy in T-cell lymphoma.

Long non-coding RNA (lncRNA) are emerging as powerful and versatile regulators of transcriptional programs and distinctive biomarkers of progression of T-cell lymphoma. Their role in the aggressive anaplastic lymphoma kinase-negative (ALK-) subtype of anaplastic large cell lymphoma (ALCL) has been elucidated only in part. Starting from our previously identified ALCL-associated lncRNA signature and performing digital gene expression profiling of a retrospective cohort of ALCL, we defined an 11 lncRNA signature able to discriminate among ALCL subtypes. We selected a not previously characterized lncRNA, MTAAT, with preferential expression in ALK- ALCL, for molecular and functional studies. We demonstrated that lncRNA MTAAT contributes to an aberrant mitochondrial turnover restraining mitophagy and promoting cellular proliferation. Functionally, lncRNA MTAAT acts as a repressor of a set of genes related to mitochondrial quality control via chromatin reorganization. Collectively, our work demonstrates the transcriptional role of lncRNA MTAAT in orchestrating a complex transcriptional program sustaining the progression of ALK- ALCL.

Intelectin enhances the phagocytosis of macrophages via CDC42-WASF2-ARPC2 signaling axis in Megalobrama amblycephala

Intelectin has been identified in various vertebrates and plays an important role in the host immune system. In our previous studies, recombinant Megalobrama amblycephala intelectin (rMaINTL) protein with excellent bacterial binding and agglutination activities enhances the phagocytic and killing activities of macrophages in M. amblycephala; however, the underlying regulatory mechanisms remain unclear. The present study showed that treatment with Aeromonas hydrophila and LPS induced the expression of rMaINTL in macrophages, and its level and distribution in macrophages or kidney tissue markedly increased after incubation or injection with rMaINTL. The cellular structure of macrophages was significantly affected after incubation with rMaINTL, resulting in an increased surface area and pseudopodia extension, which might contribute to enhancing the phagocytic ability of macrophages. Then, digital gene expression profiling analysis of the kidneys from rMaINTL-treated juvenile M. amblycephala identified some phagocytosis-related signaling factors that were enriched in pathways involved in the regulation of the actin cytoskeleton. In addition, qRT-PCR and western blotting verified that rMaINTL upregulated the expression of CDC42, WASF2, and ARPC2 in vitro and in vivo; however, the expression of these proteins was inhibited by a CDC42 inhibitor in macrophages. Moreover, CDC42 mediated the promotion of rMaINTL on actin polymerization by increasing the F-actin/G-actin ratio, which led to the extension of pseudopodia and remodeling of the macrophage cytoskeleton. Furthermore, the enhancement of macrophage phagocytosis by rMaINTL was blocked by the CDC42 inhibitor. These results suggested that rMaINTL induced the expression of CDC42 as well as the downstream signaling molecules WASF2 and ARPC2, thereby facilitating actin polymerization to promote cytoskeletal remodeling and phagocytosis. Overall, MaINTL enhanced the phagocytosis activity of macrophages in M. amblycephala via activation of the CDC42-WASF2-ARPC2 signaling axis.

Transcriptome analysis of Ganoderma lingzhi (Agaricomycetes) response to Trichoderma hengshanicum infection.

Green mold caused by Trichoderma spp. has become one of the most serious diseases which threatening the production of Ganoderma lingzhi. To understand the possible resistance mechanism of the G. lingzhi response to T. hengshanicum infection, we examined the G. lingzhi transcript accumulation at 0, 12, and 24 h after T. hengshanicum inoculation. The gene expression analysis was conducted on the interaction between G. lingzhi and T. hengshanicum using RNA-seq and digital gene expression (DGE) profiling methods. Transcriptome sequencing indicated that there were 162 differentially expressed genes (DEGs) at three infection time points, containing 15 up-regulated DEGs and 147 down-regulated DEGs. Resistance-related genes thaumatin-like proteins (TLPs) (PR-5s), phenylalanine ammonia-lyase, and Beta-1,3-glucan binding protein were significantly up-regulated. At the three time points of infection, the heat shock proteins (HSPs) genes of G. lingzhi were down-regulated. The down-regulation of HSPs genes led to the inhibition of HSP function, which may compromise the HSP-mediated defense signaling transduction pathway, leading to G. lingzhi susceptibility. Pathway enrichment analyses showed that the main enriched pathways by G. lingzhi after infection were sphingolipid metabolism, ether lipid metabolism, and valine, leucine and isoleucine degradation pathway. Overall, the results described here improve fundamental knowledge of molecular responses to G. lingzhi defense and contribute to the design of strategies against Trichoderma spp.

Genomic and proteomic characterization of Philadelphia-like B-lineage acute lymphoblastic leukemia: A report of Indian patients.

The gold standard for the identification of Philadelphia (Ph)-like acute lymphoblastic leukemia (ALL) patients is gene expression profiling. Because of its diverse nature, its identification is extremely difficult and expensive. On the genomic and proteomic landscape of Ph-like ALL patients, there is a paucity of published literature from developing countries. The authors used digital barcoded nCounter NanoString gene expression profiling for its detection, followed by molecular and proteomic characterization using fluorescence in situ hybridization and liquid chromatography-tandem mass spectrometry (LC-MS/MS). The authors found 32.05% Ph-like ALL patients and their median age at presentation was considerably higher than Ph-negative ALL cases (p=.0306). Furthermore, we identified 20% CRLF2 overexpressed cases having 8.33% CRLF2-IGH translocation with concomitant R683S mutation and 8.33% CRLF2-P2RY8 translocation. In 80% of CRLF2 downregulated cases, we identified 10% as having JAK2 rearrangement. Minimal residual disease-positivity was more common in Ph-like ALL cases (55.55% vs. 25% in Ph-negative ALL cases). Immunoglobulin J chain (Jchain), small nuclear ribonucleoprotein SmD1 (SNRPD1), immunoglobulin κ constant (IGKC), NADH dehydrogenase (ubiquinone) 1 α subcomplex subunit 2 (NDUFA2), histone H2AX (H2AFX), charged multivesicular body protein 4b (CHMP4B), and carbonyl reductase (NADPH) (CBR1) proteins were identified to be substantially expressed in Ph-like ALL patients, using LC-MS/MS. Gene enrichment analysis indicated that involvement of spliceosomal mediated messenger RNA splicing pathway and four microRNAs was statistically significant in Ph-like ALL patients. For the first time, we have described incidence, molecular, and proteomic characterization of Ph-like ALL, in developing nations. In developing countries, detecting Philadelphia (Ph)-like B-lineage acute lymphoblastic leukemia is complicated and challenging due to its diverse genetic landscape. There is no well-defined and cost-effective methodology for its detection. The incidence of this high-risk subtype is very high in adult cases, and there is an urgent need for its accurate detection. We have developed an online PHi-RACE classifier for its rapid detection, followed by delineating the genomic and proteomic landscape of Ph-like acute lymphoblastic leukemias for the first time in Indian patients.

Comparative Expression Profiling Reveals the Regulatory Effects of Dietary Mannan Oligosaccharides on the Intestinal Immune Response of Juvenile Megalobrama amblycephala against Aeromonas hydrophila Infection.

Mannan oligosaccharides (MOS) are functional oligosaccharides with beneficial effects on the non-specific immunity of Megalobrama amblycephala, but systematic studies on the immunomodulatory mechanisms of MOS are still lacking. To investigate the protective mechanisms of three different levels of dietary MOS supplementation on the intestinal immunity of juvenile M. amblycephala, comparative digital gene expression (DGE) profiling was performed. In this study, 622 differentially expressed genes (DEGs) were identified, while the similar expression tendency of 34 genes by qRT-PCR validated the accuracy of the DGE analyses. Gene Ontology (GO) enrichment revealed that the DEGs were mainly enriched in two functional categories of biological process and molecular function. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the DEGs were mainly related to complement and coagulation cascades, coagulation cascades, platelet activation, natural killer cell mediated cytotoxicity, Fc gamma R-mediated phagocytosis and antigen processing and presentation. In addition, the pro-inflammatory, apoptosis and tight junction-related genes were more significantly up-regulated upon infection in the dietary MOS groups to enhance host immune functions and maintain the stability of the intestinal barrier. These results will be helpful to clarify the regulatory mechanism of MOS on the intestinal immunity of M. amblycephala and lay the theoretical foundation for the prevention and protection of fish bacterial diseases.

Comparative transcriptome, digital gene expression and proteome profiling analyses provide insights into the brachyurization from the megalopa to the first juvenile in Eriocheir sinensis

Eriocheir sinensis larva normally experiences 11 stages. The reduced abdomen folded beneath the thorax is the most prominent characteristic of morphological alteration from megalopa to juvenile crab. Up to date, the molecular mechanisms of brachyurization remain a mystery. Here, transcriptome library, digital gene expression (DGE) libraries and proteome libraries at two developmental stages [the megalopa stage of E. sinensis (stage M) and the first stage of juvenile crab (stage J1)] of the Chinese mitten crab larva were constructed for RNA sequencing and iTRAQ approaches followed by bioinformatics analysis, respectively. In total, 1106 genes and 871 proteins were differentially expressed between the stage M and stage J1. Moreover, several important pathways were identified, including biosynthesis of secondary metabolites, metabolic pathways, focal adhesion, and some disease pathways. Besides, muscle contraction, oxidative phosphorylation, calcium signaling, PI3K-Akt, DNA replication pathway, and integrin signaling pathway also had important functions in brachyurization process. Furthermore, the components, actin, actin-related protein, collagens, filamin-A/B, laminin, integrins, paxillin, and fibronectin had up-regulated expression levels in M stage compared to J1 stage.

Assessment of 2-Year Storage Conditions on Protein, RNA, and DNA in Unstained Human Tissue Sections, Including a Novel Multiplex Digital Gene Expression Profiling Method with Implications for Biobanking.

Background: Formalin-fixed, paraffin-embedded (FFPE) tissues are a valuable resource for clinical and basic science research. Paraffin blocks and the resulting unstained sections (USS) are often stored for years before being used. Previous studies have evaluated the effects of time, temperature, humidity, and inert gases on preservation of USS; however, no study has examined all four variables together. Methods: In the current work, we prospectively and blindly assessed time points from 0 to 24 months, room versus refrigerated temperature, and presence of a desiccant and/or nitrogen atmosphere on a variety of benign and malignant tissues from North America and Africa. End points included immunohistochemistry (IHC), in situ hybridization (ISH), extracted RNA and DNA quantity and quality, and messenger RNA performance in a novel, multiplexed digital gene expression profiling assay of both housekeeping and tumor-specific genes. Results: We found that using current methods of antigen retrieval, staining, and extraction, the end points of IHC, ISH, RNA, and DNA were well preserved under the various conditions tested, with implications that pre-embedding factors contribute to variability in subsequent tissue integrity. We also document that spectrophotometric estimations of nucleic acid concentrations were in general estimated to be higher than with fluorimetric methods, which may be pertinent to end assay development. We further describe a new multiplex assay, the PlexSet digital gene expression assay, suitable for evaluating RNA quality in FFPE tissues. Conclusion: Altogether, these results may provide helpful guidance with regard to approaches for long-term storage conditions for USS.

Characterization and Clinical Impact of the Tumor Microenvironment of Post-Transplant Large B-Cell Lymphoma

Post-transplant lymphoproliferative disorders (PTLDs) arise after solid organ or allogeneic hematopoietic stem cell transplantation. They are characterized by abnormal proliferation of lymphoid cells associated with immunosuppression and comprise three histological categories (non-destructive, polymorphic, monomorphic) entailing a heterogeneous clinical behavior. The majority of monomorphic PTLDs fulfill the criteria of aggressive large B-cell lymphomas and are often associated with Epstein-Barr virus (EBV) positivity. The composition and function of the tumor microenvironment (TME) has a major impact on lymphoma development and survival, but to date, only few studies investigating the TME in PTLD have been reported. We used digital gene expression profiling and in silico immunophenotyping to characterize the TME of 75 post-transplant monomorphic aggressive B-cell lymphomas from three Nordic countries and correlated the findings with patient demographics (Table 1) and survival. Unsupervised clustering demonstrated a high degree of heterogeneity and variation in the immune response profile among the PTLDs. EBV+ PTLDs clustered closely together and were enriched for type I interferon signaling and antiviral response genes. In addition, a signature for cytotoxicity associated with EBV-positivity and better overall survival (OS; HR 0.61, 95%CI 0.40-0.92, P=0.019). In silico immunophenotyping utilizing CIBERSORTx suggested that the PTLDs can be divided into two subgroups with distinct immune cell compositions (Figure 1). The patients in the inflamed subgroup with high proportions of T cells and NK cells, as well as high T cell/macrophage ratio, had significantly better OS and PFS compared to the patients in the non-inflamed subgroup (median OS >200.0 vs. 15.2 months, P=0.006; median PFS 89.0 vs. 14.4 months, P=0.049). In multivariable analysis with the EBV-status, International Prognostic Index (IPI) and rituximab-containing treatment, the inflamed TME remained as an independent predictor for favorable OS and PFS. The inflamed subgroup was associated with EBV-positivity (P=0.008) and earlier onset of the disease (P=0.004). We also compared the TME composition between post-transplant and de novo diffuse large B-cell lymphomas (DLBCLs, n=82) and discovered that the PTLD TME was enriched in dendritic cells, NK cells, and activated mast cells, whereas the proportion of T cells, memory B cells and resting mast cells was lower compared to the de novo DLBCL. In conclusion, we provide a comprehensive phenotypic characterization of PTLD, highlighting the importance of TME biology, and in particular T cell infiltration, in the clinical behavior and prognosis of PTLD. Figure 1View largeDownload PPTFigure 1View largeDownload PPT Close modal

Tumor Microenvironment and It's Prognostic Value Is Different in Testicular and Systemic Diffuse Large B-Cell Lymphoma

Testicular diffuse large B-cell lymphoma (T-DLBCL) is an aggressive lymphoid cancer with distinct clinical and biological features. Majority of the cases represent non-GCB phenotype, but the genetic landscape with frequent aberrations in MYD88, NFKBIZ, PDCD1LG2, and CD274 differ from systemic non-GCB DLBCL and highlight the significance of immune-escape and tumor microenvironment (TME). We have previously shown that high T-cell gene expression signature and high tumor-infiltrating lymphocyte (TIL) content translate to better survival in testicular but not in systemic DLBCL. We have also described a subpopulation of T-bet+FoxP3+ Tregs that have a significant adverse impact on survival in T-DLBCL but not in in systemic DLBCL. In order to further characterize the differences in the TME between T- and systemic DLBCL, we used NanoString digital gene expression profiling to assess the expression of 770 immune response genes in 60 T-DLBCL and 84 systemic DLBCL samples. In addition, we used multiplex immunohistochemistry (mIHC) to characterize TILs, tumor-associated macrophages (TAMs), NK-cells, checkpoint molecules, cell proportions, and cell interactions in 80 T-DLBCLs and 192 systemic DLBCLs. Unsupervised hierarchical clustering of the gene expression data revealed clusters of genes with differential expression between T- and systemic DLBCLs. In T-DLBCL, 11 HLA genes were downregulated and genes related to B-cell receptor (BCR) signaling upregulated compared to systemic DLBCL. All T-DLBCL patients that clustered among the patients with systemic DLBCL had high T-cell signature expression, and systemic DLBCL patients that clustered among patients with T-DLBCL had significantly higher Ki-67 expression (p=0.00085). Based on the mIHC data, we observed not only a significantly lower proportion of FoxP3+ Tregs altogether but also a higher proportion of T-bet+FoxP3+ double positive Tregs in T-DLBCL compared to systemic DLBCL, confirming our previous findings on the significant role of T-bet+FoxP3+ double positive Tregs in the T-DLBCL TME. We also observed less TAMs in T-DLBCL TME. Unsupervised hierarchical clustering of the cell interaction data verified that distinct interaction clusters are found and separate the patients into different subgroups. Interactions of cytotoxic T-cells were more commonly seen in T-DLBCL whereas interactions of PD-1 expressing cells and interactions between TILs and TAMs were more frequent in systemic DLBCL. Taken together, our results reveal differences in the TME between T- and systemic DLBCL, and further emphasize the impact of immune escape and distinct immune cells in T-DLBCL. Future studies are ongoing in order to further characterize the cell-to-cell interactions, their association with survival, and TME differences between T- and systemic DLBCL, in particular non-GCB DLBCL.

Newly Diagnosed, Untreated, Multiple Myeloma (MM) Patient Samples Already Harbor Cereblon (CRBN) Exon 10 Deletions Associated with Drug Resistance

Background: Acquisition of drug resistance, particularly immunomodulatory drugs (IMiDs) in MM is one of the most important challenges in patient treatment today. Early on in the disease, the plasma cells in the MM clone are often sensitive to chemotherapy. However, over time, after exposure to chemotherapy, the cells acquire increasing genomic complexity and mutations leading to resistance to therapy and ultimately patient death. Cereblon (CRBN) is a binding protein, which is required for the action of the widely used IMiDs such as lenalidomide and pomolidomide. CRBN is also essential to the newer CRBN E3 ligase modulator drugs and proteolysis targeting chimeras. We recently described an algorithm based on digital gene expression profiling of purified CD138(+) MM cells from patients at different clinical stages that we term the "MM-IP-7", which correlated with clinical drug resistance, poor survival, and disease progression following combination treatment containing IMIDs and/or PIs (XY Zhu et al, eHaem 2022). Given the recent description of genetic changes in the CRBN gene that are associated with IMiD resistance particularly the exon 10-spliced CRBN transcripts that increase in incidence in parallel with IMiD-refractory states (S. Gooding et al, Blood 2021), we sought to build an integrated, single-platform test to measure both gene expression levels and genetic deletions. Methods: We designed a novel CRBN splicing isoform probe with homology to both exon 9 and exon 11. This probe will thus capture CRBN mRNA sequences in which exon 10 is deleted and portions of exons 9 and 11 and juxtaposed. We also included probes to CRBN overall (Exons 5 and 6). These probes were combined with the MM-IP-7 probes to create a dual gene expression and genetic assay that was used to analyze 30 cases of newly diagnosed MM, 30 cases of end stage MM, and 3 MM cell lines (MM1.S, OPM2, XG1). To define increased exon 10 deleted allele frequency, we calculated the ratio of exon 10 deletion/total exon 10. Cases with apparent exon 10 deletion by expression profiling, were further evaluated with RT-PCR. Results: 40% MM samples in the end stage group versus 23% of samples in the newly diagnosed group were identified with more than 20% exon 10 deleted mRNA as compared to total CRBN mRNA levels. Three cases were shown to have more than 45% exon 10 deleted mRNA, including two newly diagnosed cases and one end stage case. RT-PCR confirmed exon 10 deletion in a newly diagnosed patient sample with sufficient remaining RNA that was tested (cloning and sequencing for confirmation are underway). Cell line results showed very low exon 10 splicing variant, which correlated with our mRNAseq data for these lines. Conclusions: Our findings demonstrate that inclusion of probes bridging CRBN exon 10 deletion with homology to exons 9 and 11, is feasible approach to use on the gene expression platform, which allows for combined evaluation of expression and genetic changes. CRBN exon 10 deletion is a known treatment resistance feature in MM, which may already be present at diagnosis and selected by subsequent therapy with possible therapeutic planning implications. Additional studies to confirm the results and expand the study set are underway.

Effects of DNA methylation and gene expression on rats with protein malnutrition in early life

Although the mechanisms underlying how malnutrition in early life affects the susceptibility to diseases later in life remain unclear, there is considerable interest in the potential role of DNA methylation in the intrauterine programming of diseases. In this study, digital gene expression profiles were used to analyze changes in gene expression of protein-restricted rats early in life, while intergenerational rat models were used to explore differences in whole blood genomic DNA methylation in an environment stimulated by maternal protein intervention. Nine rats were randomly divided into early-life low-protein group (LPE), fetal low-protein group (LPF), and normal control group (CON). The LPE group was fed a low-protein diet on the 1st day of pregnancy until the end of lactation. The LPF group was given low-protein diet during pregnancy. The CON group was given a 20% protein diet from the 1st day of pregnancy. Total mRNA was extracted from the sacrificed rats at the 48th week. The number of differentially expressed genes for LPE versus CON, LPF versus CON, and LPF versus LPE was 178, 223, and 302, respectively. Comparing LPE versus CON and LPF versus CON, the upregulated genes common to the two groups were Gimap-9, Serinc-4, Dnah-2, Sf3b-5, and Sat-2, and the downregulated genes were Ppp1r-3. Comparing LPF versus CON and LPF versus LPE, the upregulated genes were Mgat2 and Cars, and the downregulated genes were Ddx28 and Slc12a9. The differentially expressed genes were mainly related to cell metabolism, immune response, signaling pathway, endocrine metabolism, stress response, ATP binding, and other functions. Early-life protein malnutrition affects gene expression of rat offspring and involves multiple aspects of growth and development, with different stages of early-life malnutrition leading to altered DNA methylation expression of corresponding genes, mainly in mitochondrial genes.

Digital gene expression profiling of the transcriptional response to Sclerotinia sclerotiorum and its antagonistic bacterium Bacillus amyloliquefaciens in soybean.

Soybean Sclerotinia stem rot caused by Sclerotinia sclerotiorum is a common disease in soybean, and effective biological control is urgently needed. We have previously confirmed that Bacillus amyloliquefaciens can effectively antagonize S. sclerotiorum in a plate competition experiment and a soybean seedling inoculation experiment. In this study, the mechanisms underlying plant death caused by S. sclerotiorum and soybean resistance to S. sclerotiorum induced by B. amyloliquefaciens were evaluated. The stems of potted soybean seedlings were inoculated with S. sclerotiorum (Gm-Ss), B. amyloliquefaciens (Gm-Ba), and their combination (Gm-Ba-Ss), using scratch treatments as a control, followed by dual RNA sequencing and bioinformatics analyses. Global gene expression levels in the Gm-Ss treatment were much lower than those in the Gm-Ba, Gm-Ba-Ss, and Gm groups, suggesting that S. sclerotiorum strongly inhibited global gene expression in soybean. In a pairwise comparison of Gm-Ss vs. Gm, 19983 differentially expressed genes (DEGs) were identified. Down-regulated DEGs were involved in various KEGG pathways, including ko01110 (biosynthesis of secondary metabolites), ko01100 (metabolic pathways), ko01120 (microbial metabolism in diverse environments), ko00500 (starch and sucrose metabolism), and ko04075 (plant hormone signal transmission), suggesting that S. sclerotiorum inoculation had a serious negative effect on soybean metabolism. In Gm-Ba vs. Gm, 13091 DEGs were identified, and these DEGs were significantly enriched in ko03010 (ribosome) and ko03008 (ribosome biogenesis in eucaryotes). Our results suggest that B. amyloliquefaciens increases the expression of genes encoding the ribosomal subunit, promotes cell wall biogenesis, and induces systemic resistance. S. sclerotiorum strongly inhibited metabolism in soybean, inhibited the synthesis of the cytoskeleton, and induced the up-regulation of programmed death and senescence-related genes via an ethylene signal transduction pathway. These results improve our understanding of S. sclerotiorum-induced plant death and soybean resistance to S. sclerotiorum induced by B. amyloliquefaciens and may contribute to the improvement of strategies to avoid yield losses.

Identification of castration-dependent and -independent driver genes and pathways in castration-resistant prostate cancer (CRPC)

BackgroundProstate cancer (PCa) is one of the most diagnosed cancers in the world. PCa inevitably progresses to castration-resistant prostate cancer (CRPC) after androgen deprivation therapy treatment, and castration-resistant state means a shorter survival time than other causes. Here we aimed to define castration-dependent and -independent diver genes and molecular pathways in CRPC which are responsible for such lethal metastatic events.MethodsBy employing digital gene expression (DGE) profiling, the alterations of the epididymal gene expression profile in the mature and bilateral castrated rat were explored. Then we detect and characterize the castration-dependent and -independent genes and pathways with two data set of CPRC-associated gene expression profiles publicly available on the NCBI.ResultsWe identified 1,632 up-regulated and 816 down-regulated genes in rat’s epididymis after bilateral castration. Differential expression analysis of CRPC samples compared with the primary PCa samples was also done. In contrast to castration, we identified 97 up-regulated genes and 128 down-regulated genes that changed in both GEO dataset and DGE profile, and 120 up-regulated genes and 136 down-regulated genes changed only in CRPC, considered as CRPC-specific genes independent of castration. CRPC-specific DEGs were mainly enriched in cell proliferation, while CRPC-castration genes were associated with prostate gland development. NUSAP1 and NCAPG were identified as key genes, which might be promising biomarkers of the diagnosis and prognosis of CRPC.ConclusionOur study will provide insights into gene regulation of CRPC dependent or independent of castration and will improve understandings of CRPC development and progression.