Exploring the Mechanisms of Yishen Tongluo Decoction on Repairing DNA Damage in Mouse Spermatogonia Cells Based on Whole Transcriptome Sequencing

The aim of this study was to investigate the potential mechanism through which Yishen Tongluo decoction (YSTL) repairs DNA damage caused by benzo(a)pyrene diol epoxide (BPDE) in mouse spermatocytes (GC-2). The GC-2 cells were divided randomly into the control group, BPDE group, and low-, medium-, and high-dose YSTL groups of YSTL decoction. A comet assay was used to detect the DNA fragment index (DFI) of cells in each group. Based on the DFI results, whole transcriptome sequencing was conducted, followed by trend analysis, gene ontology (GO) enrichment analysis, kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis, and ceRNA network analysis. Compared with the control group, the BPDE group reported a significant increase in the DNA fragmentation index (DFI) (p < .05). Compared with the BPDE group, the low-, high- and medium-dose YSTL groups had a significantly reduced DFI (p < .05). Whole-transcriptome sequencing revealed seven differentially expressed circRNAs, 203 differentially expressed miRNAs, and 3,662 differentially expressed mRNAs between the control group and the BPDE group. There was a total of 12 differentially expressed circRNAs, 204 miRNAs, and 2150 mRNAs between the BPDE group and the traditional Chinese medicine group. The pathways involved include DNA repair pathway, nucleotide excision repair pathway, base excision repair pathway, etc. The ceRNA network reported that Hmga2 was the core protein involved, novel_cir_000117 and mmu-miR-466c-3p were located upstream of Hmga2, and they were regulatory factors associated with Hmga2. Finally, we conclude that YSTL decoction may repair sperm DNA damage caused by BPDE through the novel_cir_000117-mmu-miR-466c-3p-Hmga2 pathway.

Signaling pathways and regulatory networks in quail skeletal muscle development: insights from whole transcriptome sequencing

Quail, as an advantageous avian model organism due to its compact size and short reproductive cycle, holds substantial potential for enhancing our understanding of skeletal muscle development. The quantity of skeletal muscle represents a vital economic trait in poultry production. Unraveling the molecular mechanisms governing quail skeletal muscle development is of paramount importance for optimizing meat and egg yield through selective breeding programs. However, a comprehensive characterization of the regulatory dynamics and molecular control underpinning quail skeletal muscle development remains elusive. In this study, through the application of HE staining on quail leg muscle sections, coupled with preceding fluorescence quantification PCR of markers indicative of skeletal muscle differentiation, we have delineated embryonic day 9 (E9) and embryonic day 14 (E14) as the start and ending points, respectively, of quail skeletal muscle differentiation. Then, we employed whole transcriptome sequencing to investigate the temporal expression profiles of leg muscles in quail embryos at the initiation of differentiation (E9) and upon completion of differentiation (E14). Our analysis revealed the expression patterns of 12,012 genes, 625 lncRNAs, 14,457 circRNAs, and 969 miRNAs in quail skeletal muscle samples. Differential expression analysis between the E14 and E9 groups uncovered 3,479 differentially expressed mRNAs, 124 lncRNAs, 292 circRNAs, and 154 miRNAs. Furthermore, enrichment analysis highlighted the heightened activity of signaling pathways related to skeletal muscle metabolism and intermuscular fat formation, such as the ECM-receptor interaction, focal adhesion, and PPAR signaling pathway during E14 skeletal muscle development. Conversely, the E9 stage exhibited a prevalence of pathways associated with myoblast proliferation, exemplified by cell cycle processes. Additionally, we constructed regulatory networks encompassing lncRNA‒mRNA, miRNA‒mRNA, lncRNA‒miRNA-mRNA, and circRNA-miRNA‒mRNA interactions, thus shedding light on their putative roles within quail skeletal muscle. Collectively, our findings illuminate the gene and non-coding RNA expression characteristics during quail skeletal muscle development, serving as a foundation for future investigations into the regulatory mechanisms governing non-coding RNA and quail skeletal muscle development in poultry production.

Case report: Molecular profiling facilitates the diagnosis of a challenging case of lung cancer with choriocarcinoma features

Accurate diagnoses are crucial in determining the most effective treatment across different cancers. In challenging cases, morphology-based traditional pathology methods have important limitations, while molecular profiling can provide valuable information to guide clinical decisions. We present a 35-year female with lung cancer with choriocarcinoma features. Her disease involved the right lower lung, brain, and thoracic lymph nodes. The pathology from brain metastasis was reported as "metastatic choriocarcinoma" (a germ cell tumor) by local pathologists. She initiated carboplatin and etoposide, a regimen for choriocarcinoma. Subsequently, her case was assessed by pathologists from an academic cancer center, who gave the diagnosis of "adenocarcinoma with aberrant expression of β-hCG" and finally pathologists at our hospital, who gave the diagnosis of "poorly differentiated carcinoma with choriocarcinoma features". Genomic profiling detected a KRAS G13R mutation and transcriptomics profiling was suggestive of lung origin. The patient was treated with carboplatin/paclitaxel/ipilimumab/nivolumab followed by consolidation radiation therapy. She had no evidence of progression to date, 16 months after the initial presentation. The molecular profiling could facilitate diagnosing of challenging cancer cases. In addition, chemoimmunotherapy and local consolidation radiation therapy may provide promising therapeutic options for patients with lung cancer exhibiting choriocarcinoma features.

Novel MED15::ATF1 fusion in a pediatric melanoma with spitzoid features and aggressive presentation

Childhood melanoma is a rare and biologically heterogeneous pediatric malignancy. The differential diagnosis of pediatric melanoma is usually broad, including a wide variety of spindle cell or epithelioid neoplasms. Different molecular alterations affecting the MAPK and PI3K/AKT/mTOR pathways, tumor suppressor genes, and telomerase reactivation have been implicated in melanoma tumorigenesis and progression. Here, we report a novel MED15::ATF1 fusion in a pediatric melanoma with spitzoid features and an aggressive clinical course.

Novel HNRNPM::LEUTX fusion resulting from chromothripsis of chromosome 19 in a pediatric undifferentiated small round cell neoplasm

Small round cell neoplasms comprise a diverse group of tumors characterized by a primitive/undifferentiated appearance. Although several entities are associated with recurrent gene fusions, many of these neoplasms have not been fully characterized, and novel molecular alterations are being discovered. Here, we report an undifferentiated small round cell neoplasm arising in the anterior mediastinum of a 17-month-old female. The tumor harbored a novel HNRNPM::LEUTX fusion resulting from chromothripsis of chromosome 19, which was identified by whole transcriptome sequencing, but not by targeted sequencing. The structural variations caused by the chromothripsis event also challenged the interpretation of the targeted sequencing findings. This report expands the spectrum of gene partners involved in LEUTX fusions and underscores the value of whole transcriptome sequencing in the diagnostic workup of undifferentiated small round cell tumors. It also highlights the interpretive challenges associated with complex genomic alterations. A careful evidence-based analysis of sequencing data along with histopathologic correlation is essential to ensure correct categorization of fusions.

Analytic validation of NeXT Dx™, a comprehensive genomic profiling assay

We describe the analytic validation of NeXT Dx, a comprehensive genomic profiling assay to aid therapy and clinical trial selection for patients diagnosed with solid tumor cancers. Proprietary methods were utilized to perform whole exome and whole transcriptome sequencing for detection of single nucleotide variants (SNVs), insertions/deletions (indels), copy number alterations (CNAs), and gene fusions, and determination of tumor mutation burden and microsatellite instability. Variant calling is enhanced by sequencing a patient-specific normal sample from, for example, a blood specimen. This provides highly accurate somatic variant calls as well as the incidental reporting of pathogenic and likely pathogenic germline alterations. Fusion detection via RNA sequencing provides more extensive and accurate fusion calling compared to DNA-based tests. NeXT Dx features the proprietary Accuracy and Content Enhanced technology, developed to optimize sequencing and provide more uniform coverage across the exome. The exome was validated at a median sequencing depth of >500x. While variants from 401 cancer-associated genes are currently reported from the assay, the exome/transcriptome assay is broadly validated to enable reporting of additional variants as they become clinically relevant. NeXT Dx demonstrated analytic sensitivities as follows: SNVs (99.4%), indels (98.2%), CNAs (98.0%), and fusions (95.8%). The overall analytic specificity was >99.0%.

Primary refractory plasmablastic lymphoma: A precision oncology approach

Introduction

Hematologic malignancies are currently underrepresented in multidisciplinary molecular-tumor-boards (MTB). This study assesses the potential of precision-oncology in primary-refractory plasmablastic-lymphoma (prPBL), a highly lethal blood cancer.

Methods

We evaluated clinicopathological and molecular-genetic data of 14 clinically annotated prPBL-patients from initial diagnosis. For this proof-of-concept study, we employed our certified institutional MTB-pipeline (University-Cancer-Center-Schleswig-Holstein, UCCSH) to annotate a comprehensive dataset within the scope of a virtual MTB-setting, ultimately recommending molecularly stratified therapies. Evidence-levels for MTB-recommendations were defined in accordance with the NCT/DKTK and ESCAT criteria.

Results

Median age in the cohort was 76.5 years (range 56-91), 78.6% of patients were male, 50% were HIV-positive and clinical outcome was dismal. Comprehensive genomic/transcriptomic analysis revealed potential recommendations of a molecularly stratified treatment option with evidence-levels according to NCT/DKTK of at least m2B/ESCAT of at least IIIA were detected for all 14 prPBL-cases. In addition, immunohistochemical-assessment (CD19/CD30/CD38/CD79B) revealed targeted treatment-recommendations in all 14 cases. Genetic alterations were classified by treatment-baskets proposed by Horak et al. Hereby, we identified tyrosine-kinases (TK; n=4), PI3K-MTOR-AKT-pathway (PAM; n=3), cell-cycle-alterations (CC; n=2), RAF-MEK-ERK-cascade (RME; n=2), immune-evasion (IE; n=2), B-cell-targets (BCT; n=25) and others (OTH; n=4) for targeted treatment-recommendations. The minimum requirement for consideration of a drug within the scope of the study was FDA-fast-track development.

Discussion

The presented proof-of-concept study demonstrates the clinical potential of precision-oncology, even in prPBL-patients. Due to the aggressive course of the disease, there is an urgent medical-need for personalized treatment approaches, and this population should be considered for MTB inclusion at the earliest time.

Identification of a Novel Long Non-coding RNA, lnc-ATMIN-4:2, and its Clinicopathological and Prognostic Significance in Advanced Gastric Cancer

BACKGROUND/AIM

Long non-coding RNAs (lncRNAs) are emerging as significant regulators of gene expression and a novel promising biomarker for cancer diagnosis and prognosis. This study identified a novel, differentially expressed lncRNA in advanced gastric cancer (AGC), Inc-ATMIN-4:2, and evaluated its clinicopathological and prognostic significance.

PATIENTS AND METHODS

Whole transcriptome sequencing was performed to identify differentially expressed lncRNAs in AGC tissue samples. We also analyzed lnc-ATMIN-4:2 expression in 317 patients with AGC using RNA in situ hybridization.

RESULTS

High (>30 dots) lnc-ATMIN-4:2 expression significantly correlated with younger age, poorly differentiated histology, diffuse type, deeper invasion depth, perineural invasion, lymph node metastasis, and higher stage group. In addition, high lnc-ATMIN-4:2 expression was significantly associated with worse overall survival in patients with AGC.

CONCLUSION

This study elucidated the significance of lncRNAs in AGC and indicated the value of lnc-ATMIN-4:2 expression as a predictive biomarker for the overall survival of patients with AGC.

Analysis of ceRNA networks during mechanical tension-induced osteogenic differentiation of periodontal ligament stem cells

The molecular mechanisms underlying osteogenic differentiation of periodontal ligament stem cells (PDLSCs) under mechanical tension remain unclear. This study aimed to identify a potential long non-coding ribonucleic acids (lncRNAs)/circular RNAs (circRNAs)-microRNAs (miRNAs)-messenger RNAs (mRNAs) network in mechanical tension-induced osteogenic differentiation of PDLSCs. PDLSCs were isolated from the healthy human periodontal ligament, identified, cultured, and exposed to tensile force. The expression of osteogenic markers was examined, and whole transcriptome sequencing was performed to identify the expression patterns of lncRNA, circRNA, miRNAs, and mRNAs. Enrichment analyses were also performed. Candidate targets of differentially expressed non-coding RNAs (ncRNAs) were predicted, and potential competitive endogenous RNA (ceRNA) networks were constructed by Cytoscape. We found that the osteogenic differentiation of PDLSCs was significantly enhanced under dynamic tension (magnitude: 12%, frequency: 0.7 Hz). Overall, 344 lncRNAs, 57 miRNAs, 41 circRNAs, and 70 mRNAs were differentially expressed in the tension group and the control group. Functional enrichment analysis showed that differentially expressed mRNAs were mainly enriched in osteogenesis-related and mechanical stress-related biological processes and signal transduction pathways (e.g., tumor necrosis factor [TNF] and Hippo signaling pathways). The lncRNA/circRNA-miRNA-mRNA networks were depicted, and potential key ceRNA networks were identified. Our findings may help to further explore the underlying regulatory mechanism of osteogenic differentiation of PDLSCs under mechanical tensile stress.

Integrated Analysis of lncRNA and circRNA Mediated ceRNA Regulatory Networks in Skin Reveals Innate Immunity Differences Between Wild-Type and Yellow Mutant Rainbow Trout (Oncorhynchus mykiss)

Fish skin is a vital immune organ that forms the first protective barrier preventing entry of external pathogens. Rainbow trout is an important aquaculture fish species that is farmed worldwide. However, our knowledge of innate immunity differences between wild-type (WR_S) and yellow mutant rainbow trout (YR_S) remains limited. In this study, we performed whole transcriptome analysis of skin from WR_S and YR_S cultured in a natural flowing water pond. A total of 2448 mRNAs, 1630 lncRNAs, 22 circRNAs and 50 miRNAs were found to be differentially expressed (DE). Among these DEmRNAs, numerous key immune-related genes, including ifih1, dhx58, trim25, atp6v1e1, tap1, tap2, cd209, hsp90a.1, nlrp3, nlrc3, and several other genes associated with metabolism (gstp1, nampt, naprt and cd38) were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEmRNAs revealed that many were significantly enriched in innate immune-related GO terms and pathways, including NAD+ADP-ribosyltransferase activity, complement binding, immune response and response to bacterium GO terms, and RIG-I-like receptor signaling, NOD-like receptor signaling and phagosome KEGG pathways. Furthermore, the immune-related competing endogenous RNA networks were constructed, from which we found that lncRNAs MSTRG.11484.2, MSTRG.32014.1 and MSTRG.29012.1 regulated at least three immune-related genes (ifih1, dhx58 and irf3) through PC-5p-43254_34, PC-3p-28352_70 and bta-miR-11987_L-1R-1_1ss8TA, and tap2 was regulated by two circRNAs (circRNA5279 and circRNA5277) by oni-mir-124a-2-p5_1ss13GA. The findings expand our understanding of the innate immune system of rainbow trout, and lay the foundation for further study of immune mechanisms and disease resistance breeding.

Transcriptional Profiling of Malignant Melanoma Reveals Novel and Potentially Targetable Gene Fusions

Simple Summary

Malignant melanoma is a complex disease that is estimated to claim over 7000 lives in the United States in 2021. Although recent advances in genomic technology have helped with the identification of driver variants, molecular studies and clinical trials have often focused on prevalent alterations, such as the BRAF-V600E mutation. With the inclusion of whole transcriptome sequencing, molecular profiling of melanomas has identified gene fusions and revealed gene expression profiles that are consistent with the activation of signaling pathways by common driver mutations. Patients harboring such fusions may benefit from currently approved targeted therapies and should be considered in the design of future clinical trials to further personalize treatments for patients with malignant melanoma.

Abstract

Invasive melanoma is the deadliest type of skin cancer, with 101,110 expected cases to be diagnosed in 2021. Recurrent BRAF and NRAS mutations are well documented in melanoma. Biologic implications of gene fusions and the efficacy of therapeutically targeting them remains unknown. Retrospective review of patient samples that underwent next-generation sequencing of the exons of 592 cancer-relevant genes and whole transcriptome sequencing for the detection of gene fusion events and gene expression profiling. Expression of PDL1 and ERK1/2 was assessed by immunohistochemistry (IHC). There were 33 (2.6%) cases with oncogenic fusions (14 novel), involving BRAF, RAF1, PRKCA, TERT, AXL, and FGFR3. MAPK pathway-associated genes were over-expressed in BRAF and RAF1 fusion-positive tumors in absence of other driver alterations. Increased expression in tumors with PRKCA and TERT fusions was concurrent with MAPK pathway alterations. For a subset of samples with available tissue, increased phosphorylation of ERK1/2 was observed in BRAF, RAF1, and PRKCA fusion-positive tumors. Oncogenic gene fusions are associated with transcriptional activation of the MAPK pathway, suggesting they could be therapeutic targets with available inhibitors. Additional analyses to fully characterize the oncogenic effects of these fusions may support biomarker driven clinical trials.

The Competitive Endogenous RNA (ceRNA) Regulation in Porcine Alveolar Macrophages (3D4/21) Infected by Swine Influenza Virus (H1N1 and H3N2)

H1N1 and H3N2 are the two most common subtypes of swine influenza virus (SIV). They not only endanger the pig industry, but are also a huge risk of zoonotic diseases. However, the molecular mechanism and regulatory network of pigs (hosts) against influenza virus infection are still unclear. In this study, porcine alveolar macrophage cell (3D4/21) models infected by swine influenza virus (H1N1 and H3N2) were constructed. The expression profiles of miRNAs, mRNAs, lncRNAs and circRNAs after H1N1 and H3N2 infected 3D4/21 cells were revealed in this study. Then, two ceRNAs (TCONS_00166432-miR10391-MAN2A1 and novel_circ_0004733-miR10391-MAN2A1) that regulated H1N1 and H3N2 infection in 3D4/21 cells were verified by the methods of bioinformatics analysis, gene overexpression, gene interference, real-time quantitative PCR (qPCR), dual luciferase activity assay and RNA immunoprecipitation (RIP). In addition, the important candidate molecules (miR-10391, TCONS_00166432, and novel_circ_0004733) were identified by qPCR and enzyme linked immunosorbent assay (ELISA). Finally, the regulatory effect and possible molecular mechanism of the target gene MAN2A1 were identified by the methods of gene interference, qPCR, Western blot and ELISA. The results of this study suggested that TCONS_00166432 and novel_circ_0004733 could competitively bind miR-10391 to target the MAN2A1 gene to regulate swine influenza virus infecting 3D4/21 cells. This study reported for the first time the ceRNA networks involved in the regulation of the swine influenza virus infecting 3D4/21 cells, which provided a new insight into the molecular mechanism of 3D4/21 cells against swine influenza virus infection.

NRG1 and NRG2 fusion positive solid tumor malignancies: a paradigm of ligand-fusion oncogenesis

Neuregulins (NRGs) are a family of six related physiological ligands all containing a receptor-binding epidermal growth factor (EGF)-like domain that mediate their binding to cellular receptors. Neuregulin-1 (NRG1) is the main physiological ligand to HER3. NRG1 fusion (NRG1+) was first reported in a breast cancer cell line and NRG2 fusions have recently been identified in solid tumors. It is postulated that NRG1 fusions, through mostly transmembrane fusion partners, result in NRG1 being concentrated in proximity to HER3, leading to its constitutive activation and oncogenesis. Recently, a monoclonal antibody that disrupts the binding of NRG1 to HER3 and HER3/HER2 heterodimerization has resulted in NRG1+ tumor shrinkage, suggesting that 'ligand-fusion' may be a novel mechanism of oncogenesis.

Integrative Studies of Human Cord Blood Derived Mononuclear Cells and Umbilical Cord Derived Mesenchyme Stem Cells in Ameliorating Bronchopulmonary Dysplasia

Bronchopulmonary dysplasia (BPD) is a common pulmonary complication observed in preterm infants that is composed of multifactorial pathogenesis. Current strategies, albeit successful in moderately reducing morbidity and mortality of BPD, failed to draw overall satisfactory conclusion. Here, using a typical mouse model mimicking hallmarks of BPD, we revealed that both cord blood-derived mononuclear cells (CB-MNCs) and umbilical cord-derived mesenchymal stem cells (UC-MSCs) are efficient in alleviating BPD. Notably, infusion of CB-MNCs has more prominent effects in preventing alveolar simplification and pulmonary vessel loss, restoring pulmonary respiratory functions and balancing inflammatory responses. To further elucidate the underlying mechanisms within the divergent therapeutic effects of UC-MSC and CB-MNC, we systematically investigated the long noncoding RNA (lncRNA)-microRNA (miRNA)-messenger RNA (mRNA) and circular RNA (circRNA)-miRNA-mRNA networks by whole-transcriptome sequencing. Importantly, pathway analysis integrating Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG)/gene set enrichment analysis (GSEA) method indicates that the competing endogenous RNA (ceRNA) network is mainly related to the regulation of GTPase activity (GO: 0043087), extracellular signal-regulated kinase 1 (ERK1) and ERK2 signal cascade (GO: 0070371), chromosome regulation (GO: 0007059), and cell cycle control (GO: 0044770). Through rigorous selection of the lncRNA/circRNA-based ceRNA network, we demonstrated that the hub genes reside in UC-MSC- and CB-MNC-infused networks directed to the function of cell adhesion, motor transportation (Cdk13, Lrrn2), immune homeostasis balance, and autophagy (Homer3, Prkcd) relatively. Our studies illustrate the first comprehensive mRNA-miRNA-lncRNA and mRNA-miRNA-circRNA networks in stem cell-infused BPD model, which will be valuable in identifying reliable biomarkers or therapeutic targets for BPD pathogenesis and shed new light in the priming and conditioning of UC-MSCs or CB-MNCs in the treatment of neonatal lung injury.

Genomic Landscape of Angiosarcoma: A Targeted and Immunotherapy Biomarker Analysis

Simple Summary

Angiosarcomas (AS) are rare, highly aggressive sarcomas with limited therapeutic options. Genomic sequencing techniques have identified recurrent genetic abnormalities. Nevertheless, the association of these findings with etiology, site of origin, prognosis, and therapeutic implications is not well understood. We analyzed Next Generation Sequencing (NGS) and Whole Transcriptome Sequencing (WTS) data in a cohort of 143 AS cases. We identified distinct genomic biology according to the AS primary site. Head and neck AS cases primarily have Immunotherapy (IO) response markers and mutations in TP53 and POT1. On the other hand, breast AS is enriched for cell cycle alterations, predominately MYC amplification. Additionally, a microenvironment with abundant immune cells is present in a minority of cases but distributed evenly among primary sites. Our findings can facilitate the design and optimization of therapeutic strategies for AS according to its biology at different primary sites.

Abstract

We performed a retrospective analysis of angiosarcoma (AS) genomic biomarkers and their associations with the site of origin in a cohort of 143 cases. Primary sites were head and neck (31%), breast (22%), extremity (11%), viscera (20%), skin at other locations (8%), and unknown (9%). All cases had Next Generation Sequencing (NGS) data with a 592 gene panel, and 53 cases had Whole Exome Sequencing (WES) data, which we used to study the microenvironment phenotype. The immunotherapy (IO) response biomarkers Tumor Mutation Burden (TMB), Microsatellite Instability (MSI), and PD-L1 status were the most frequently encountered alteration, present in 36.4% of the cohort and 65% of head and neck AS (H/N-AS) (p < 0.0001). In H/N-AS, TMB-High was seen in 63.4% of cases (p < 0.0001) and PDL-1 positivity in 33% of cases. The most common genetic alterations were TP53 (29%), MYC amplification (23%), ARID1A (17%), POT1 (16%), and ATRX (13%). H/N-AS cases had predominantly mutations in TP53 (50.0%, p = 0.0004), POT1 (40.5%, p < 0.0001), and ARID1A (33.3%, p = 0.5875). In breast AS, leading alterations were MYC amplification (63.3%, p < 0.0001), HRAS (16.1%, p = 0.0377), and PIK3CA (16.1%, p = 0.2352). At other sites, conclusions are difficult to generate due to the small number of cases. A microenvironment with a high immune signature, previously associated with IO response, was evenly distributed in 13% of the cases at different primary sites. Our findings can facilitate the design and optimization of therapeutic strategies for AS.

Comprehensive Analysis of the Expression Profiles of Hepatic lncRNAs in the Mouse Model of Alcoholic Liver Disease

Background and Aim: The worldwide prevalence of alcoholic liver disease (ALD) due to escalating alcohol consumption has presented an unprecedented pressure on human health. A few studies have determined long non-coding RNAs (lncRNAs) involved in the pathogenesis of liver diseases. However, the roles of lncRNAs in ALD development is still poorly understood.

Methods: An ALD mouse model was established and confirmed. Expression profiles of lncRNAs were obtained by whole transcriptome sequencing. The altered lncRNAs in ALD mice were further verified by qRT-PCR. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were used to enrich the functions of these lncRNAs. In combination with miRNA and mRNA profiles, we constructed concise endogenous RNA (ceRNA) networks. The function of the most up/downregulated lnRNA was further verified and investigated in both ALD model and AML-12 cells.

Results: Totally, five downregulated lncRNAs were obtained and verified in ALD mice. The GO term and KEGG pathway analyses revealed that the identified lncRNAs were associated with alcohol-induced hepatic oxidative damage, cellular inflammation, and lipid metabolism. Combination the differentially modulated miRNAs and mRNAs with ceRNA network analysis, we constructed five ceRNA networks and obtained 30 miRNAs and 25 mRNAs that may participate in ALD. Further, we verified and investigate the function of the most downregulated lnc_1700023H06Rik. Depletion lnc_1700023H06Rik reduced genes encoding for lipid metabolism, especially mRNA Acat2 (ENSMUST00000159697) and Pgrmc2 (ENSMUST00000058578) both in vivo and in vitro. Knocking down lnc_1700023H06Rik induced triglyceride accumulation and lactate dehydrogenase leakage in AML12 cells, consisting with that in alcohol-treated cells.

Conclusion: The five remarkably downregulated lncRNAs in ALD mouse model were identified as novel biomarkers, highlighting the key role of lncRNAs in the development of ALD. The effect of lnc_1700023H06Rik plays a pivotal role in lipid deposition and its pathological pathway in ALD needs further investigation.

Males and Females Have Distinct Molecular Events in the Articular Cartilage during Knee Osteoarthritis

Osteoarthritis (OA) is a major public health challenge that imposes a remarkable burden on the affected individuals and the healthcare system. Based on the clinical observation, males and females have different prevalence rates and severity levels of OA. Thus, sex-based differences may play essential roles in OA’s prognosis and treatment outcomes. To date, the comprehensive understanding of the relationship between sex and OA is still largely lacking. In the current study, we analyzed a published transcriptome dataset of knee articular cartilage (GSE114007) from 18 healthy (five females, 13 males) and 20 OA (11 females, nine males) donors to provide a slight insight into this important but complex issue. First, comparing female healthy cartilage samples with those of males revealed 36 differential expression genes (DEGs), indicating the fundamental sex-related differences at the molecular level. Meanwhile, 923 DEGs were distinguished between OA and healthy female cartilage, which can be enriched to 15 Reactome pathways. On the other hand, when comparing OA and healthy male cartilage, there are only 419 DEGs were identified, and only six pathways were enriched against the Reactome database. The different signaling response to OA in the male and female cartilage was further enforced by recognizing 50 genes with significantly different OA-responsive expression fold changes in males and females. Particularly, 14 Reactome pathways, such as “Extracellular matrix organization”, “Collagen biosynthesis and modifying enzymes”, “Dissolution of fibrin clot”, and “Platelet Aggregation (Plug formation)”, can be noted from these 50 sex-dependent OA-responsive genes. Overall, the current study explores the Sex as a Biological Variable (SABV) at the transcriptomic level in the knee articular cartilage in both healthy status and OA event, which could help predict the differential OA prognosis and treatment outcome of males and female patients.

Identification of Novel Fusion Transcripts in High Grade Serous Ovarian Cancer

Fusion genes are structural chromosomal rearrangements resulting in the exchange of DNA sequences between genes. This results in the formation of a new combined gene. They have been implicated in carcinogenesis in a number of different cancers, though they have been understudied in high grade serous ovarian cancer. This study used high throughput tools to compare the transcriptome of high grade serous ovarian cancer and normal fallopian tubes in the interest of identifying unique fusion transcripts within each group. Indeed, we found that there were significantly more fusion transcripts in the cancer samples relative to the normal fallopian tubes. Following this, the role of fusion transcripts in chemo-response and overall survival was investigated. This led to the identification of fusion transcripts significantly associated with overall survival. Validation was performed with different analytical platforms and different algorithms to find fusion transcripts.

RNA interactions in right ventricular dysfunction induced type II cardiorenal syndrome

Right ventricular (RV) dysfunction induced type II cardiorenal syndrome (CRS) has a high mortality rate, but little attention has been paid to this disease, and its unique molecular characteristics remain unclear. This study aims to investigate the transcriptomic expression profile in this disease and identify key RNA pairs that regulate related molecular signaling networks. We established an RV dysfunction-induced type II CRS mouse model by pulmonary artery constriction (PAC). PAC mice developed severe RV hypertrophy and fibrosis; renal atrophy and dysfunction with elevated creatinine were subsequently observed. Expression profiles in RV and kidney tissues were obtained by whole transcriptome sequencing, revealing a total of 741 and 86 differentially expressed (DE) mRNAs, 159 and 29 DEmiRNAs and 233 and 104 DEcircRNAs between RV and kidney tissue, respectively. Competing endogenous RNA (ceRNA) networks were established. A significant alteration in proliferative, fibrotic and metabolic pathways was found based on GO and KEGG analyses, and the network revealed key ceRNA pairs, such as novel_circ_002631/miR-181a-5p/Creb1 and novel_circ_002631/miR-33-y/Kpan6. These findings indicate that significantly dysregulated pathways in RV dysfunction induced type II CRS include Ras, PI3K/Akt, cGMP-PKG pathways, and thyroid metabolic pathways. These ceRNA pairs can be considered potential targets for the treatment of type II CRS.

PDA Indolylmaleimides Induce Anti-Tumor Effects in Prostate Carcinoma Cell Lines Through Mitotic Death

Castrate resistant prostate cancer in men shares several characteristics with canine prostate cancer (PCa). Due to current insufficient therapies, evaluating novel therapeutic agents for late-stage PCa is of considerable interest for both species. PDA indolylmaleimides showed anticancer effects in several neoplastic cell lines. Herein, a comparative characterization of PDA-66 and PDA-377 mediated effects was performed in human and canine PCa cell lines, which is also the first detailed characterization of these agents on cells derived from solid tumors in general. While PDA-377 showed only weak growth inhibition on human PCa cell lines, PDA-66 inhibited proliferation and induced apoptosis in human and canine cell lines with concentrations in the low micromolar range. Morphological characterization and whole transcriptome sequencing revealed that PDA-66 induces mitotic death through its microtubule-depolymerizing ability. PDA-66 appears to be a worthwhile anti-mitotic agent for further evaluation. The similarities in cellular and molecular response observed in the cell lines of both origins form a solid basis for the use of canine PCa in vivo models to gain valuable interchangeable data to the advantage of both species.

Cartilage Targets of Knee Osteoarthritis Shared by Both Genders

As the leading cause of disability, osteoarthritis (OA) affects people of all ages, sexes, and races. With the increasing understanding of OA, the sex differences have attracted specific attention as the burden of OA is greater in women. There is no doubt that gender-specific OA management has great potential for precision treatment. On the other hand, from the marketing aspect, a medication targeting the OA-responsive biomarker(s) shared by both genders is more favorable for drug development. Thus, in the current study, a published transcriptome dataset of knee articular cartilage was used to compare OA and healthy samples for identifying the genes with the same significantly different expression trend in both males and females. With 128 genes upregulated and 143 genes downregulated in both OA males and females, 9 KEGG pathways have been enriched based on the current knowledge, including 'renal cell carcinoma,' 'ECM-receptor interaction,' 'HIF-1 signaling pathway,' 'MicroRNAs in cancer,' 'focal adhesion,' 'Relaxin signaling pathway,' 'breast cancer,' 'PI3K-Akt signaling pathway,' and 'human papillomavirus infection.' Here, we explore the potential impacts of these clusters in OA. We also analyze the identified 'cell plasma membrane related genes' in-depth to identify the potential chondrocyte cell surface target(s) of OA management.

Proper Read Filtering Method to Adequately Analyze Whole-Transcriptome Sequencing and RNA Based Immune Repertoire Sequencing Data for Tumor Milieu Research

Simple Summary

The recent advancement in high-throughput sequencing has become indispensable for immune-genomics and profiling the T- and B-cell receptor repertoires. Immune repertoire sequencing (IR-seq) and whole transcriptome sequencing (WTS) can be implemented to investigate and quantitatively characterize the complex pattern of the CDR3 region. We conducted T-cell diversity analysis result comparisons of these sequencing methods and suggest an intuitive approach to discriminate reliable TCR sequences and clonotype patterns from capturing errors. Although bulk-RNA sequencing is commonly used for cancer analysis, we confirmed capturing highly enriched TCR transcripts with IR-seq is more reliable for accurate immune repertoire discovery, and singleton read filtering criteria should be applied to capture true clonotypes from error-prone sequencing data. The use of such well-established data and analytical methodologies can broaden understanding of antigen specificity in immunity and enabling efficient therapeutic antibody finding.

Abstract

Analysis of the T-cell receptor (TCR) repertoire is essential to characterize the extensive collections of T-cell populations with recognizing antigens in cancer research, and whole transcriptome sequencing (WTS) and immune repertoire sequencing (IR-seq) are commonly used for this measure. To date, no standard read filtering method for IR measurement has been presented. We assessed the diversity of the TCR repertoire results from the paired WTS and IR-seq data of 31 multiple myeloma (MM) patients. To invent an adequate read filtering strategy for IR analysis, we conducted comparisons with WTS results. First, our analyses for determining an optimal threshold for selecting clonotypes showed that the clonotypes supported by a single read largely affected the shared clonotypes and manifested distinct patterns of mapping qualities, unlike clonotypes with multiple reads. Second, although IR-seq could reflect a wider TCR region with a higher capture rate than WTS, an adequate comparison with the removal of unwanted bias from potential sequencing errors was possible only after applying our read filtering strategy. As a result, we suggest that TCR repertoire analysis be carried out through IR-seq to produce reliable and accurate results, along with the removal of single-read clonotypes, to conduct immune research in cancer using high-throughput sequencing.

Changing Expression Profiles of Messenger RNA, MicroRNA, Long Non-coding RNA, and Circular RNA Reveal the Key Regulators and Interaction Networks of Competing Endogenous RNA in Pulmonary Fibrosis

Pulmonary fibrosis is a kind of interstitial lung disease with architectural remodeling of tissues and excessive matrix deposition. Apart from messenger RNA (mRNA), microRNA (miRNA), long non-coding RNA (lncRNA), and circular RNA (circRNA) could also play important roles in the regulatory processes of occurrence and progression of pulmonary fibrosis. In the present study, the pulmonary fibrosis model was administered with bleomycin. Whole transcriptome sequencing analysis was applied to investigate the expression profiles of mRNAs, lncRNAs, circRNAs, and miRNAs. After comparing bleomycin-induced pulmonary fibrosis model lung samples and controls, 286 lncRNAs, 192 mRNAs, 605 circRNAs, and 32 miRNAs were found to be differentially expressed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to investigate the potential functions of these differentially expressed (DE) mRNAs and non-coding RNAs (ncRNAs). The terms related to inflammatory response and tumor necrosis factor (TNF) signaling pathway were enriched, implying potential roles in regulatory process. In addition, two co-expression networks were also constructed to understand the internal regulating relationships of these mRNAs and ncRNAs. Our study provides a systematic perspective on the potential functions of these DE mRNAs and ncRNAs during PF process and could help pave the way for effective therapeutics for this devastating and complex disease.

Identification of Novel Targets of Knee Osteoarthritis Shared by Cartilage and Synovial Tissue

Arthritis is the leading cause of disability among adults, while osteoarthritis (OA) is the most common form of arthritis that results in cartilage loss. However, accumulating evidence suggests that the protective hyaline cartilage should not be the sole focus of OA treatment. Particularly, synovium also plays essential roles in OA’s initiation and progression and warrants serious consideration when battling against OA. Thus, biomarkers with similar OA-responsive expressions in cartilage and synovium should be the potential targets for OA treatment. On the other hand, molecules with a distinguished response during OA in cartilage and synovium should be ruled out as OA therapeutic(s) to avoid controversial effects in different tissues. Here, to pave the path for developing a new generation of OA therapeutics, two published transcriptome datasets of knee articular cartilage and synovium were analyzed in-depth. Genes with statistically significantly different expression in OA and healthy cartilage were compared with those in the synovium. Thirty-five genes with similar OA-responsive expression in both tissues were identified while recognizing three genes with opposite OA-responsive alteration trends in cartilage and synovium. These genes were clustered based on the currently available knowledge, and the potential impacts of these clusters in OA were explored.

Whole transcriptome analysis on blue light-induced eye damage

AIM

To analyze abnormal gene expressions of mice eyes exposed to blue light using RNA-seq and analyze the related signaling pathways.

METHODS

Kunming mice were divided into an experimental group that was exposed to blue light and a control group that was exposed to natural light. After 14d, the mice were euthanized and their eyeballs were collected. Whole transcriptome analysis was attempted to analyze the gene expression of the eyeballs using RNA-seq to reconstruct genetic networks. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were used to reveal the related signaling pathways.

RESULTS

The 737 differentially expressed genes were identified, including 430 up and 307 down regulated genes, by calculating the gene FPKM in each sample and conducting differential gene analysis. GO and KEGG pathway enrichment analysis showed that blue light damage may associated with the visual perception, sensory perception of light stimulus, phototransduction, and JAK-STAT signaling pathways. Differential lncRNA, circRNA and miRNA analysis showed that blue light exposure affected pathways for retinal cone cell development and phototransduction, among others.

CONCLUSION

Exposure to blue light can cause a certain degree of abnormal gene expression and modulate signaling pathways in the eye.

Integrative Molecular Characterization of Resistance to Neoadjuvant Chemoradiation in Rectal Cancer

PURPOSE

Molecular properties associated with complete response or acquired resistance to concurrent chemotherapy and radiotherapy (CRT) are incompletely characterized.Experimental Design: We performed integrated whole-exome/transcriptome sequencing and immune infiltrate analysis on rectal adenocarcinoma tumors prior to neoadjuvant CRT (pre-CRT) and at time of resection (post-CRT) in 17 patients [8 complete/partial responders, 9 nonresponders (NR)].

RESULTS

CRT was not associated with increased tumor mutational burden or neoantigen load and did not alter the distribution of established somatic tumor mutations in rectal cancer. Concurrent KRAS/TP53 mutations (KP) associated with NR tumors and were enriched for an epithelial-mesenchymal transition transcriptional program. Furthermore, NR was associated with reduced CD4/CD8 T-cell infiltrates and a post-CRT M2 macrophage phenotype. Absence of any local tumor recurrences, KP/NR status predicted worse progression-free survival, suggesting that local immune escape during or after CRT with specific genomic features contributes to distant progression.

CONCLUSIONS

Overall, while CRT did not impact genomic profiles, CRT impacted the tumor immune microenvironment, particularly in resistant cases.

A Possible Mechanism: Genistein Improves Metabolism and Induces White Fat Browning Through Modulating Hypothalamic Expression of Ucn3, Depp, and Stc1

Bioactive food components have gained growing attention in recent years. Multiple studies demonstrated that genistein had beneficial effects on metabolism. However, the exact mechanism by which genistein improves metabolism remains unclear, especially the central regulation. This study was designed to evaluate whether addition of genistein to the high-fat diet could counter metabolic disorders and whether these alterations were associated with gene expression in hypothalamus. C57BL/6 mice were fed either a high-fat diet (HF), high-fat diet with genistein (0.25 g/kg diet) (HFG) or a normal control diet (CON) for 8 weeks. Body weight was assessed during the study. After 8-week intervention, content of inguinal subcutaneous adipose tissue (SAT), perirenal visceral adipose tissue (VAT) and brown adipose tissue (BAT) were weighed. Glucose tolerance test, the serum levels of insulin and lipid were assessed. The mRNA of browning marker was detected in the white fat. The hypothalamus was collected for whole transcriptome sequencing and reverse transcription quantitative PCR validation. The results demonstrated that mice fed HFG diet had lower body weight and SAT mass, decrease levels of low-density lipoprotein cholesterol and free fatty acids, higher browning marker of Ucp1 and Cidea in WAT and an improvement in glucose tolerance and insulin sensitivity compared with those in HF group. Transcriptome sequencing showed that there were three differentially expressed genes in hypothalamus among the three groups, including Ucn3, Depp, and Stc1, which were significantly correlated with the browning markers in WAT and insulin sensitivity. Thus, regulating gene expressions in hypothalamus is a potential mechanism for genistein improving metabolism and inducing WAT browning, which may provide a novel target for the precaution and treatment of T2DM.

Whole-transcriptome sequencing identified gene expression signatures associated with aggressive clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most prevalent subtype of kidney cancer, yet molecular biomarkers have not been used for the prognosis of ccRCC to aide clinical decision making. This study aimed to identify genes associated with ccRCC aggressiveness and overall survival (OS). Samples of ccRCC tumor tissue were obtained from 33 patients who underwent nephrectomy. Gene expression was determined using whole-transcriptome sequencing. The Cancer Genome Atlas Kidney Renal Clear Cell Carcinoma (TCGA-KIRC) RNA-seq data was used to test association with OS. 290 genes were differentially expressed between tumors with high and low stage, size, grade, and necrosis (SSIGN) score (≥7 vs. ≤3) with P_ADJ<0.05. Four genes, G6PD, APLP1, GCNT3, and PLPP2, were also over-expressed in advanced stage (III and IV) and high grade (3 and 4) ccRCC and tumor with necrosis (P_ADJ<0.05). Investigation stratifying by stage found that APLP1 and PLPP2 overexpression were significantly associated with poorer OS in the early stage (Quartile 1 vs. Quartile 4, HR = 3.87, 95% CI:1.25-11.97, P = 0.02 and HR = 4.77, 95% CI:1.37-16.57, P = 0.04 respectively). These genes are potential biomarkers of ccRCC aggressiveness and prognosis that direct clinical and surgical management.

Utility of DNA, RNA, Protein, and Functional Approaches to Solve Cryptic Immunodeficiencies

PURPOSE

We report a female infant identified by newborn screening for severe combined immunodeficiencies (NBS SCID) with T cell lymphopenia (TCL). The patient had persistently elevated alpha-fetoprotein (AFP) with IgA deficiency, and elevated IgM. Gene sequencing for a SCID panel was uninformative. We sought to determine the cause of the immunodeficiency in this infant.

METHODS

We performed whole-exome sequencing (WES) on the patient and parents to identify a genetic diagnosis. Based on the WES result, we developed a novel flow cytometric panel for rapid assessment of DNA repair defects using blood samples. We also performed whole transcriptome sequencing (WTS) on fibroblast RNA from the patient and father for abnormal transcript analysis.

RESULTS

WES revealed a pathogenic paternally inherited indel in ATM. We used the flow panel to assess several proteins in the DNA repair pathway in lymphocyte subsets. The patient had absent phosphorylation of ATM, resulting in absent or aberrant phosphorylation of downstream proteins, including γH2AX. However, ataxia-telangiectasia (AT) is an autosomal recessive condition, and the abnormal functional data did not correspond with a single ATM variant. WTS revealed in-frame reciprocal fusion transcripts involving ATM and SLC35F2 indicating a chromosome 11 inversion within 11q22.3, of maternal origin. Inversion breakpoints were identified within ATM intron 16 and SLC35F2 intron 7.

CONCLUSIONS

We identified a novel ATM-breaking chromosome 11 inversion in trans with a pathogenic indel (compound heterozygote) resulting in non-functional ATM protein, consistent with a diagnosis of AT. Utilization of several molecular and functional assays allowed successful resolution of this case.

Identification and comparison of novel circular RNAs with associated co-expression and competing endogenous RNA networks in pulmonary tuberculosis

Pulmonary tuberculosis (PTB) is caused by Mycobacterium tuberculosis and is one of the most serious diseases worldwide. Circular RNAs (circRNAs) are a large class of non-coding RNAs that were identified with potential regulatory roles in disease pathogenesis and progression. In this study, we used whole transcriptome sequencing to identify circRNAs from 3 PTB patients and 3 healthy individuals to determine the expression pattern of circRNAs in blood and the circRNA molecular regulatory networks in PTB pathogenesis. One hundred and seventy differentially expressed (≥ 2-fold change) circRNAs were dysregulated in PTB, compared with in healthy individuals. Quantitative real-time polymerase chain reaction was used to validate the RNA sequencing analysis from 20 PTB patients, and the results were consistent with the sequencing data. Gene Ontology annotation and Kyoto Encyclopedia of Genes and Genomes pathway analysis were applied to explore the potential circRNA functions of the significantly deregulated genes. Several immunity pathways, including endocytosis pathways in cancer, mitogen-activated protein kinase signaling pathway, human T-lymphotropic virus type 1 infection, and ubiquitin-mediated proteolysis, were involved in PTB pathogenesis. Competing endogenous RNAs (ceRNA) were constructed and inferred that aberrant expression of circRNA-associated ceRNA resulted in extensive variation in gene expression by miRNA-mediated circRNA-mRNA crosstalk interactions. Our study revealed that the circRNA-miRNA-mRNA network may shed light on the biological functions of circRNAs in PTB and provide useful information for exploring potential roles of circRNA in PTB.

Identification of novel follicular dendritic cell sarcoma markers, FDCSP and SRGN, by whole transcriptome sequencing

Follicular dendritic cell (FDC)-sarcoma is a rare neoplasm with morphologic and phenotypic features of FDCs. It shows an extremely heterogeneous morphology, therefore, its diagnosis relys on the phenotype of tumor cells. Aim of the present study was the identification of new specific markers for FDC-sarcoma by whole transcriptome sequencing (WTS). Candidate markers were selected based on gene expression level and biological function. Immunohistochemistry was performed on reactive tonsils, on 22 cases of FDC-sarcomas and 214 control cases including 114 carcinomas, 87 soft tissue tumors, 5 melanomas, 5 thymomas and 3 interdigitating dendritic cell sarcomas. FDC secreted protein (FDCSP) and Serglycin (SRGN) proved to be specific markers of FDC and related tumor. They showed better specificity and sensitivity values than some well known markers used in FDC sarcoma diagnosis (specificity: 98.6%, and 100%, respectively; sensitivity: 72.73% and 68.18%, respectively). In our cohorts CXCL13, CD21, CD35, FDCSP and SRGN were the best markers for FDC-sarcoma diagnosis and could discriminate 21/22 FDC sarcomas from other mesenchymal tumors by linear discriminant analysis. In summary, by WTS we identified two novel FDC markers and by the analysis of a wide cohort of cases and controls we propose an efficient marker panel for the diagnosis of this rare and enigmatic tumor.

Classification based on mutations of TERT promoter and IDH characterizes subtypes in grade II/III gliomas

BACKGROUND

Grade II and III gliomas have variable clinical behaviors, showing the distinct molecular genetic alterations from glioblastoma (GBM), many of which eventually transform into more aggressive tumors. Since the classifications of grade II/III gliomas based on the genetic alterations have been recently emerging, it is now a trend to include molecular data into the standard diagnostic algorithm of glioma.

METHODS

Here we sequenced TERT promoter mutational status (TERTp-mut) in the DNA of 377 grade II/III gliomas and analyzed the clinical factors, molecular aberrations, and transcriptome profiles.

RESULTS

We found that TERTp-mut occurred in 145 of 377 grade II and III gliomas (38.5%), mutually exclusive with a TP53 mutation (TP53-mut; P < .001) and coincident with a 1p/19q co-deletion (P = .002). TERTp-mut was an independent predictive factor of a good prognosis in all patients (P = .048). It has been an independent factor associated with a good outcome in the IDH mutation (IDH-mut) subgroup (P = .018), but it has also been associated with a poor outcome in the IDH wild-type (IDH-wt) subgroup (P = .049). Combining TERTp-mut and IDH-mut allowed the grade II/III malignancies to be reclassified into IDH-mut/TERTp-mut, IDH-mut only, TERTp-mut only, and IDH-wt/TERTp-wt. 1p/19q co-deletion, TP53-muts, Ki-67 expression differences, and p-MET expression differences characterized IDH-mut/TERTp-mut, IDH-mut only, TERTp-mut only, and IDH-wt/TERTp-wt subtypes, respectively.

CONCLUSIONS

Our results showed that TERTp-mut combined with IDH-mut allowed simple classification of grade II/III gliomas for stratifying patients and clarifying diagnostic accuracy by supplementing standard histopathological criteria.

Genetic and clinical characteristics of primary and secondary glioblastoma is associated with differential molecular subtype distribution

Glioblastoma multiforme (GBM) is classified into primary (pGBM) or secondary (sGBM) based on clinical progression. However, there are some limits to this classification for insight into genetically and clinically distinction between pGBM and sGBM. The aim of this study is to characterize pGBM and sGBM associating with differential molecular subtype distribution. Whole transcriptome sequencing data was used to assess the distribution of molecular subtypes and genetic alterations in 88 pGBM and 34 sGBM in a Chinese population-based cohort, and the biological progression and prognostic impact were analyzed by combining clinical information. Forty-one percentage of pGBM were designated as Mesenchymal subtype, while only 15% were the Proneural subtype. However, sGBM displayed the opposite ratio of Mesenchymal (15%) and Proneural (44%) subtypes. Mutations in isocitrate dehydrogenase-1 (IDH1) were found to be highly concentrated in the Proneural subtypes. In addition, patients with sGBM were 10 years younger on average than those with pGBM, and exhibited clinical features of shorter overall survival and frontal lobe tumor location tendency. Furthermore, in sGBM, gene sets related to malignant progression were found to be enriched. Overall, these results reveal the intrinsic distinction between pGBM and sGBM, and provide insight into the genetic and clinical attributes of GBM.

Whole transcriptome sequencing identifies BCOR internal tandem duplication as a common feature of clear cell sarcoma of the kidney

PURPOSE

Clear cell sarcoma of the kidney (CCSK) is a rare pediatric renal tumor that is frequently difficult to distinguish among other childhood renal tumors due to its histological heterogeneity. This work evaluates genetic abnormalities carried by a series of CCSK samples by whole transcriptome sequencing (WTS), to identify molecular biomarkers that could improve the diagnostic process.

METHODS

WTS was performed on tumor RNA from 8 patients with CCSK. Bioinformatic analysis, with implementation of a pipeline for detection of intragenic rearrangements, was executed. Sanger sequencing and gene expression were evaluated to validate BCOR internal tandem duplication (ITD).

RESULTS

WTS did not identify any shared SNVs, Ins/Del or fusion event. Conversely, analysis of intragenic rearrangements enabled the detection of a breakpoint within BCOR transcript recurrent in all samples. Three different in-frame ITD in exon15 of BCOR, were detected. The presence of the ITD was confirmed on tumor DNA and cDNA, and resulted in overexpression of BCOR.

CONCLUSIONS

WTS coupled with specific bioinformatic analysis is able to detect rare genetic events, as intragenic rearrangements. ITD in the last exon of BCOR is recurrent in all CCSK samples analyzed, representing a valuable molecular marker to improve diagnosis of this rare childhood renal tumor.

Instruction of haematopoietic lineage choices, evolution of transcriptional landscapes and cancer stem cell hierarchies derived from an AML1-ETO mouse model

The t(8;21) chromosomal translocation activates aberrant expression of the AML1-ETO (AE) fusion protein and is commonly associated with core binding factor acute myeloid leukaemia (CBF AML). Combining a conditional mouse model that closely resembles the slow evolution and the mosaic AE expression pattern of human t(8;21) CBF AML with global transcriptome sequencing, we find that disease progression was characterized by two principal pathogenic mechanisms. Initially, AE expression modified the lineage potential of haematopoietic stem cells (HSCs), resulting in the selective expansion of the myeloid compartment at the expense of normal erythro- and lymphopoiesis. This lineage skewing was followed by a second substantial rewiring of transcriptional networks occurring in the trajectory to manifest leukaemia. We also find that both HSC and lineage-restricted granulocyte macrophage progenitors (GMPs) acquired leukaemic stem cell (LSC) potential being capable of initiating and maintaining the disease. Finally, our data demonstrate that long-term expression of AE induces an indolent myeloproliferative disease (MPD)-like myeloid leukaemia phenotype with complete penetrance and that acute inactivation of AE function is a potential novel therapeutic option.