Development and validation of a 23-gene expression signature for molecular subtyping of medulloblastoma in a long-term Chinese cohort

PURPOSE

Medulloblastoma is the most common childhood malignant brain tumor and is a leading cause of cancer-related death in children. Recent transcriptional studies have shown that medulloblastomas comprise at least four molecular subgroups, each with distinct demographics, genetics, and clinical outcomes. Medulloblastoma subtyping has become critical for subgroup-specific therapies. The use of gene expression assays to determine the molecular subgroup of clinical specimens is a long-awaited application of molecular biology for this pediatric cancer.

METHODS

In the current study, we established a medulloblastoma transcriptome database of 460 samples retrieved from three published datasets (GSE21140, GSE37382, and GSE37418). With this database, we identified a 23-gene signature that is significantly associated with the medulloblastoma subgroups and achieved a classification accuracy of 95.2%.

RESULTS

The 23-gene signature was further validated in a long-term cohort of 142 Chinese medulloblastoma patients. The 23-gene signature classified 21 patients as WNT (15%), 41 as SHH (29%), 16 as Group 3 (11%), and 64 as Group 4 (45%). For patients of WNT, SHH, Group 3, and Group 4, 5-year overall-survival rate reached 80%, 62%, 27%, and 47%, respectively (p < 0.0001), meanwhile 5-year progression-free survival reached 80%, 52%, 27%, and 45%, respectively (p < 0.0001). Besides, SHH/TP53-mutant tumors were associated with worse prognosis compared with SHH/TP53 wild-type tumors and other subgroups. We demonstrated that subgroup assignments by the 23-gene signature and Northcott's NanoString assay were highly comparable with a concordance rate of 96.4%.

CONCLUSIONS

In conclusion, we present a novel gene signature that is capable of accurately and reliably assigning FFPE medulloblastoma samples to their molecular subgroup, which may serve as an auxiliary tool for medulloblastoma subtyping in the clinic. Future incorporation of this gene signature into prospective clinical trials is warranted to further evaluate its clinical.

Development of Nanopore sequencing-based full-length transcriptome database toward functional genome annotation of the Pacific oyster, Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is a widely cultivated shellfish in the world, while its transcriptome diversity remains less unexplored due to the limitation of short reads. In this study, we used Oxford Nanopore sequencing to develop the full-length transcriptome database of C. gigas. We identified 77,920 full-length transcripts from 21,523 genes, and uncovered 9668 alternative splicing events and 87,468 alternative polyadenylation sites. Notably, a total of 16,721 novel transcripts were annotated in this work. Furthermore, integrative analysis of 25 publicly available RNA-seq datasets revealed the transcriptome diversity involved in post-transcriptional regulation in C. gigas. We further developed a Drupal based webserver, Cgtdb, which can be used for transcriptome visualization, sequence alignment, and functional genome annotation analyses. This work provides valuable resources and a useful tool for integrative analysis of various transcriptome datasets in C. gigas, which will serve as an essential reference for functional annotation of the oyster genome.

De novo transcriptome analysis for examination of the nutrition metabolic system related to the evolutionary process through which stick insects gain the ability of flight (Phasmatodea)

Objective

Insects are the most evolutionarily successful groups of organisms, and this success is largely due to their flight ability. Interestingly, some stick insects have lost their flight ability despite having wings. To elucidate the shift from wingless to flying forms during insect evolution, we compared the nutritional metabolism system among flight-winged, flightless-winged, and flightless-wingless stick insect groups.

Results

Here, we report RNA sequencing of midgut transcriptome of Entoria okinawaensis, a prominent Japanese flightless-wingless stick insect, and the comparative analysis of its transcriptome in publicly available midgut transcriptomes obtained from seven stick insect species. A gene enrichment analysis for differentially expressed genes, including those obtained from winged vs wingless and flight vs flightless genes comparisons, revealed that carbohydrate metabolic process-related genes were highly expressed in the winged stick insect group. We also found that the expression of the mitochondrial enolase superfamily member 1 transcript was significantly higher in the winged stick insect group than in the wingless stick insect group. Our findings could indicate that carbohydrate metabolic processes are related to the evolutionary process through which stick insects gain the ability of flight.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-021-05600-0.

Peroxisome proliferator-activated receptor α agonist-induced histidine decarboxylase gene expression in the rat and mouse liver

By analysis of the data from the Toxicogenomics Database (TG-GATEs), histidine decarboxylase gene (Hdc) was identified as largely and commonly upregulated by three fibrates, clofibrate, fenofibrate, and WY-14,643, which are known to induce hepatocellular hypertrophy and proliferation via stimulation of peroxisome proliferator-activated receptor α (PPARα) in rodents. As histamine has been reported to be involved in the proliferation of liver cells, the present study was conducted to focus on Hdc. Among other genes related to histidine and histamine, the expression of the gene of histamine ammonia lyase (Hal) was exclusively mobilized by the three fibrates. The expression of Hdc, which was usually very low in the liver, was increased with the repeated administration of fibrates, and concomitantly, the constitutive expression of Hal was suppressed. An interpretation is that the formation of urocanic acid from histidine under the normal condition switches to the formation of histamine. The mobilization of gene expression of Hdc and Hal by PPARα agonists could not be reproduced in primary cultured hepatocytes. The Hdc mRNA appeared to be translated to a protein which is processed differently from brain but similarly to gastric mucosa. Surprisingly, the fibrates caused hepatic hypertrophy but no induction of Hdc mRNA at all in mice. These results revealed that the changes in the histidine catabolism by PPARα agonists might be partially, but not directly, involved in the hepatocyte proliferation in rats, and there is a large genetic distance even between rat and mouse.

Extraction of peroxisome proliferator-activated receptor α agonist-induced lipid metabolism-related and unrelated genes in rat liver and analysis of their genomic location

Although peroxisome proliferator-activated receptor α (PPARα) agonists are obviously hepatocarcinogenic in rodents, they have been widely used for dyslipidemia and proven to be safe for clinical use without respect to the species difference. It is established that PPARα acts as a part of the transcription factor complex, but its precise mechanism is still unknown. Using the data of Toxicogenomics Database, reliable genes responsive to PPARα agonists, clofibrate, fenofibrate and WY-14,643, in rat liver, were extracted from both in vivo and in vitro data, and sorted by their fold increase. It was found that there were many genes responding to fibrates exclusively in vivo. Most of the in vivo specific genes appear to be unrelated to lipid metabolism and are not upregulated in the kidney. Fifty-seven genes directly related to cell proliferation were extracted from in vivo data, but they were not induced in vitro at all. Analysis of PPAR-responsive elements could not explain the observed difference in induction. To evaluate possible interaction between neighboring genes in gene expression, the correlation of the fold changes of neighboring genes for 22 drugs with various PPARα agonistic potencies were calculated for the genes showing more than 2.5 fold induction by 3 fibrates in vivo, and their genomic location was compared with that of the human orthologue. In the present study, many candidates of genes other than lipid metabolism were selected, and these could be good starting points to elucidate the mechanism of PPARα agonist-induced rodent-specific toxicity.