Search

Abstract Background Accurate microsatellite instability (MSI) testing is essential for identifying gastric cancer (GC) patients eligible for immunotherapy. We aimed to develop and validate a CT-based radiomics signature to predict MSI and immunotherapy outcomes in GC. Methods This retrospective multicohort study included a total of 457 GC patients from two independent medical centers in China and The Cancer Imaging Archive (TCIA) databases. The primary cohort (n = 201, center 1, 2017–2022), was used for signature development via Least Absolute Shrinkage and Selection Operator (LASSO) and logistic regression analysis. Two independent immunotherapy cohorts, one from center 1 (n = 184, 2018–2021) and another from center 2 (n = 43, 2020–2021), were utilized to assess the signature’s association with immunotherapy response and survival. Diagnostic efficiency was evaluated using the area under the receiver operating characteristic curve (AUC), and survival outcomes were analyzed via the Kaplan-Meier method. The TCIA cohort (n = 29) was included to evaluate the immune infiltration landscape of the radiomics signature subgroups using both CT images and mRNA sequencing data. Results Nine radiomics features were identified for signature development, exhibiting excellent discriminative performance in both the training (AUC: 0.851, 95%CI: 0.782, 0.919) and validation cohorts (AUC: 0.816, 95%CI: 0.706, 0.926). The radscore, calculated using the signature, demonstrated strong predictive abilities for objective response in immunotherapy cohorts (AUC: 0.734, 95%CI: 0.662, 0.806; AUC: 0.724, 95%CI: 0.572, 0.877). Additionally, the radscore showed a significant association with PFS and OS, with GC patients with a low radscore experiencing a significant survival benefit from immunotherapy. Immune infiltration analysis revealed significantly higher levels of CD8 + T cells, activated CD4 + B cells, and TNFRSF18 expression in the low radscore group, while the high radscore group exhibited higher levels of T cells regulatory and HHLA2 expression. Conclusion This study developed a robust radiomics signature with the potential to serve as a non-invasive biomarker for GC’s MSI status and immunotherapy response, demonstrating notable links to post-immunotherapy PFS and OS. Additionally, distinct immune profiles were observed between low and high radscore groups, highlighting their potential clinical implications.

Abstract Background The standard evolutionary theory of ageing proposes that ageing occurs because of a trade-off between reproduction and longevity. Eusocial insect queens exhibit positive fecundity-longevity associations and so have been suggested to be counter-examples through not expressing costs of reproduction and through remodelling conserved genetic and endocrine networks regulating ageing and reproduction. If so, eusocial evolution from solitary ancestors with negative fecundity-longevity associations must have involved a stage at which costs of reproduction were suppressed and fecundity and longevity became positively associated. Using the bumblebee (Bombus terrestris), we experimentally tested whether queens in annual eusocial insects at an intermediate level of eusocial complexity experience costs of reproduction, and, using mRNA-seq, the extent to which they exhibit a remodelling of relevant genetic and endocrine networks. Specifically, we tested whether costs of reproduction are present but latent, or whether a remodelling of relevant genetic and endocrine networks has already occurred allowing queens to reproduce without costs. Results We experimentally increased queens’ costs of reproduction by removing their eggs, which caused queens to increase their egg-laying rate. Treatment queens had significantly reduced longevity relative to control queens whose egg-laying rate was not increased. Reduced longevity in treatment queens was not caused by increased worker-to-queen aggression or by increased overall activity in queens. In addition, treatment and control queens differed in age-related gene expression based on mRNA-seq in both their overall expression profiles and the expression of ageing-related genes. Remarkably, these differences appeared to occur principally with respect to relative age, not chronological age. Conclusions This study represents the first simultaneously phenotypic and transcriptomic experimental test for a longevity cost of reproduction in eusocial insect queens. The results support the occurrence of costs of reproduction in annual eusocial insects of intermediate social complexity and suggest that reproductive costs are present but latent in queens of such species, i.e. that these queens exhibit condition-dependent positive fecundity-longevity associations. They also raise the possibility that a partial remodelling of genetic and endocrine networks underpinning ageing may have occurred in intermediately eusocial species such that, in unmanipulated conditions, age-related gene expression depends more on chronological than relative age.

Abstract Background Chinese Yunnan Province, located in the Yunnan–Guizhou Plateau, is a famous tourist paradise where acute high-altitude illness common occurs among lowland people visitors due to non-acclimatization to the acute hypobaric hypoxia (AHH) conditions. Traditional Chinese medicine, such as Qi-Long-Tian (QLT) formula, has shown effectiveness and safety in the treatment of acute high-altitude diseases. The aim of this study was to clarify the therapeutic mechanisms of this traditional formula using a rat model in a simulated plateau environment. Methods Following testing, lung tissue samples were evaluated by hematoxylin–eosin staining and for biochemical characteristics. mRNA-Seq was used to compare differentially expressed genes in control rats, and in rats exposed to AHH and AHH with QLT treatment. Results Inflammation-related effectors induced following QLT treatment for AHH included MMP9 and TIMP1, and involved several phosphorylation signaling pathways implicated in AHH pathogenesis such as PI3K/AKT and MAPK signaling. Conclusion This study provides insights into the major signaling pathways induced by AHH and in the protective mechanisms involved in QLT formula activity.

Abstract Background There are challenges in generating mRNA-Seq data from whole-blood derived RNA as globin gene and rRNA are frequent contaminants. Given the abundance of erythrocytes in whole blood, globin genes comprise some 80% or more of the total RNA. Therefore, depletion of globin gene RNA and rRNA are critical steps required to have adequate coverage of reads mapping to the reference transcripts and thus reduce the total cost of sequencing. In this study, we directly compared the performance of probe hybridization (GLOBINClear Kit and Globin-Zero Gold rRNA Removal Kit) and RNAse-H enzymatic depletion (NEBNext® Globin & rRNA Depletion Kit and Ribo-Zero Plus rRNA Depletion Kit) methods from 1 μg of whole blood-derived RNA on mRNA-Seq profiling. All RNA samples were treated with DNaseI for additional cleanup before the depletion step and were processed for poly-A selection for library generation. Results Probe hybridization revealed a better overall performance than the RNAse-H enzymatic depletion method, detecting a higher number of genes and transcripts without 3′ region bias. After depletion, samples treated with probe hybridization showed globin genes at 0.5% (±0.6%) of the total mapped reads; the RNAse-H enzymatic depletion had 3.2% (±3.8%). Probe hybridization showed more junction reads and transcripts compared with RNAse-H enzymatic depletion and also had a higher correlation (R > 0.9) than RNAse-H enzymatic depletion (R > 0.85). Conclusion In this study, our results showed that 1 μg of high-quality RNA from whole blood could be routinely used for transcriptional profiling analysis studies with globin gene and rRNA depletion pre-processing. We also demonstrated that the probe hybridization depletion method is better suited to mRNA sequencing analysis with minimal effect on RNA quality during depletion procedures.

Abstract Introduction Mesenchymal stem cells (MSCs) isolated from bone marrow have different developmental origins, including neural crest. MSCs can differentiate into neural progenitor-like cells (NPCs) under the influence of bFGF and EGF. NPCs can terminally differentiate into neurons that express beta-III-tubulin and elicit action potential. The main aim of the study was to identify key genetic markers involved in differentiation of MSCs into NPCs through transcriptomic analysis. Method Total RNA was isolated from MSCs and MSCs-derived NPCs followed by cDNA library construction for transcriptomic analysis. Sample libraries that passed the quality and quantity assessments were subjected to high throughput mRNA sequencing using NextSeq®500. Differential gene expression analysis was performed using the DESeq2 R package with MSC samples being a reference group. The expression of eight differentially regulated genes was counter validated using real-time PCR. Results In total, of the 3,252 differentially regulated genes between MSCs and NPCs with two or more folds, 1,771 were upregulated genes, whereas 1,481 were downregulated in NPCs. Amongst these differential genes, 104 transcription factors were upregulated, and 45 were downregulated in NPCs. Neurogenesis related genes were upregulated in NPCs and the main non-redundant gene ontology (GO) terms enriched in NPCs were the autonomic nervous system, cell surface receptor signalling pathways), extracellular structure organisation, and programmed cell death. The main non-redundant GO terms enriched in MSCs included cytoskeleton organisation cytoskeleton structural constituent, mitotic cell cycle), and the mitotic cell cycle process Gene set enrichment analysis also confirmed cell cycle regulated pathways as well as Biocarta integrin pathway were upregulated in MSCs. Transcription factors enrichment analysis by ChEA3 revealed Foxs1 and HEYL, amongst the top five transcription factors, inhibits and enhances, respectively, the NPCs differentiation of MSCs. Conclusions The vast differences in the transcriptomic profiles between NPCs and MSCs revealed a set of markers that can identify the differentiation stage of NPCs as well as provide new targets to enhance MSCs differentiation into NPCs.

Abstract Background Biological adhesion (bioadhesion), enables organisms to attach to surfaces as well as to a range of other targets. Bioadhesion evolved numerous times independently and is ubiquitous throughout the kingdoms of life. To date, investigations have focussed on various taxa of animals, plants and bacteria, but the fundamental processes underlying bioadhesion and the degree of conservation in different biological systems remain poorly understood. This study had two aims: 1) To characterise tissue-specific gene regulation in the pedal disc of the model cnidarian Exaiptasia pallida, and 2) to elucidate putative genes involved in pedal disc adhesion. Results Five hundred and forty-seven genes were differentially expressed in the pedal disc compared to the rest of the animal. Four hundred and twenty-seven genes were significantly upregulated and 120 genes were significantly downregulated. Forty-one condensed gene ontology terms and 19 protein superfamily classifications were enriched in the pedal disc. Eight condensed gene ontology terms and 11 protein superfamily classifications were depleted. Enriched superfamilies were consistent with classifications identified previously as important for the bioadhesion of unrelated marine invertebrates. A host of genes involved in regulation of extracellular matrix generation and degradation were identified, as well as others related to development and immunity. Ab initio prediction identified 173 upregulated genes that putatively code for extracellularly secreted proteins. Conclusion The analytical workflow facilitated identification of genes putatively involved in adhesion, immunity, defence and development of the E. pallida pedal disc. When defence, immunity and development-related genes were identified, those remaining corresponded most closely to formation of the extracellular matrix (ECM), implicating ECM in the adhesion of anemones to surfaces. This study therefore provides a valuable high-throughput resource for the bioadhesion community and lays a foundation for further targeted research to elucidate bioadhesion in the Cnidaria.

Abstract Background Feline injection-site sarcoma (FISS), an aggressive iatrogenic subcutaneous malignancy, is challenging to manage clinically and little is known about the molecular basis of its pathogenesis. Tumor transcriptome profiling has proved valuable for gaining insights into the molecular basis of cancers and for identifying new therapeutic targets. Here, we report the first study of the FISS transcriptome and the first cross-species comparison of the FISS transcriptome with those of anatomically similar soft-tissue sarcomas in dogs and humans. Methods Using high-throughput short-read paired-end sequencing, we comparatively profiled FISS tumors vs. normal tissue samples as well as cultured FISS-derived cell lines vs. skin-derived fibroblasts. We analyzed the mRNA-seq data to compare cancer/normal gene expression level, identify biological processes and molecular pathways that are associated with the pathogenesis of FISS, and identify multimegabase genomic regions of potential somatic copy number alteration (SCNA) in FISS. We additionally conducted cross-species analyses to compare the transcriptome of FISS to those of soft-tissue sarcomas in dogs and humans, at the level of cancer/normal gene expression ratios. Results We found: (1) substantial differential expression biases in feline orthologs of human oncogenes and tumor suppressor genes suggesting conserved functions in FISS; (2) a genomic region with recurrent SCNA in human sarcomas that is syntenic to a feline genomic region of probable SCNA in FISS; and (3) significant overlap of the pattern of transcriptional alterations in FISS with the patterns of transcriptional alterations in soft-tissue sarcomas in humans and in dogs. We demonstrated that a protein, BarH-like homeobox 1 (BARX1), has increased expression in FISS cells at the protein level. We identified 11 drugs and four target proteins as potential new therapies for FISS, and validated that one of them (GSK-1059615) inhibits growth of FISS-derived cells in vitro. Conclusions (1) Window-based analysis of mRNA-seq data can uncover SCNAs. (2) The transcriptome of FISS-derived cells is highly consistent with that of FISS tumors. (3) FISS is highly similar to soft-tissue sarcomas in dogs and humans, at the level of gene expression. This work underscores the potential utility of comparative oncology in improving understanding and treatment of FISS.

Abstract Current single-cell RNA-seq approaches are hindered by preamplification bias, loss of strand of origin information, and the inability to observe small-RNA and mRNA dual transcriptomes. Here, we introduce a single-cell holo-transcriptome sequencing (Holo-Seq) that overcomes all three hurdles. Holo-Seq has the same quantitative accuracy and uniform coverage with a complete strand of origin information as bulk RNA-seq. Most importantly, Holo-Seq can simultaneously observe small RNAs and mRNAs in a single cell. Furthermore, we acquire small RNA and mRNA dual transcriptomes of 32 human hepatocellular carcinoma single cells, which display the genome-wide super-enhancer activity and hepatic neoplasm kinetics of these cells.

Abstract Background Oleaginous microalgae are promising production hosts for the sustainable generation of lipid-based bioproducts and as bioenergy carriers such as biodiesel. Transcriptomics of the lipid accumulation phase, triggered efficiently by nitrogen starvation, is a valuable approach for the identification of gene targets for metabolic engineering. Results An explorative analysis of the detailed transcriptional response to different stages of nitrogen availability was performed in the oleaginous green alga Monoraphidium neglectum. Transcript data were correlated with metabolic data for cellular contents of starch and of different lipid fractions. A pronounced transcriptional down-regulation of photosynthesis became apparent in response to nitrogen starvation, whereas glucose catabolism was found to be up-regulated. An in-depth reconstruction and analysis of the pathways for glycerolipid, central carbon, and starch metabolism revealed that distinct transcriptional changes were generally found only for specific steps within a metabolic pathway. In addition to pathway analyses, the transcript data were also used to refine the current genome annotation. The transcriptome data were integrated into a database and complemented with data for other microalgae which were also subjected to nitrogen starvation. It is available at https://tdbmn.cebitec.uni-bielefeld.de . Conclusions Based on the transcriptional responses to different stages of nitrogen availability, a model for triacylglycerol and lipid hyperaccumulation is proposed, which involves transcriptional induction of thioesterases, differential regulation of lipases, and a re-routing of the central carbon metabolism. Over-expression of distinct thioesterases was identified to be a potential strategy to increase the oleaginous phenotype of M. neglectum, and furthermore specific lipases were identified as potential targets for future metabolic engineering approaches.

Abstract Background The filamentous fungus Aspergillus fumigatus has become the most important airborne fungal pathogen causing life-threatening infections in immuno-compromised patients. Recently developed high-throughput transcriptome and proteome technologies, such as microarrays, RNA deep-sequencing, and LC-MS/MS of peptide mixtures, are of enormous value for systematically investigating pathogenic organisms. In the field of infection biology, one of the priorities is to collect and standardise data, in order to generate datasets that can be used to investigate and compare pathways and gene responses involved in pathogenicity. The “omics” era provides a multitude of inputs that need to be integrated and assessed. We therefore evaluated the potential of paired-end mRNA-Seq for investigating the regulatory role of the central mitogen activated protein kinase (MpkA). This kinase is involved in the cell wall integrity signalling pathway of A. fumigatus and essential for maintaining an intact cell wall in response to stress. Results The comparison of the transcriptome and proteome of an A. fumigatus wild-type strain with an mpkA null mutant strain revealed that 70.4% of the genome was found to be expressed and that MpkA plays a significant role in the regulation of many genes involved in cell wall remodelling, oxidative stress and iron starvation response, and secondary metabolite biosynthesis. Moreover, absence of the mpkA gene also strongly affects the expression of genes involved in primary metabolism. The data were further processed to evaluate the potential of the mRNA-Seq technique. We comprehensively matched up our data to published transcriptome studies and were able to show an improved data comparability of mRNA-Seq experiments independently of the technique used. Analysis of transcriptome and proteome data revealed only a weak correlation between mRNA and protein abundance. Conclusions High-throughput analysis of MpkA-dependent gene expression confirmed many previous findings that this kinase is important for regulating many genes involved in metabolic pathways. Our analysis showed more than 2000 differentially regulated genes. RNA deep-sequencing is less error-prone than established microarray-based technologies. It also provides additional information in A. fumigatus studies and as a result is more suitable for the creation of extensive datasets.

Abstract Background Dysregulation of imprinted genes, which are expressed in a parent-of-origin-specific manner, plays an important role in various human diseases, such as cancer and behavioral disorder. To date, however, fewer than 100 imprinted genes have been identified in the human genome. The recent availability of high-throughput technology makes it possible to have large-scale prediction of imprinted genes. Here we propose a Bayesian model (dsPIG) to predict imprinted genes on the basis of allelic expression observed in mRNA-Seq data of independent human tissues. Results Our model (dsPIG) was capable of identifying imprinted genes with high sensitivity and specificity and a low false discovery rate when the number of sequenced tissue samples was fairly large, according to simulations. By applying dsPIG to the mRNA-Seq data, we predicted 94 imprinted genes in 20 cerebellum samples and 57 imprinted genes in 9 diverse tissue samples with expected low false discovery rates. We also assessed dsPIG using previously validated imprinted and non-imprinted genes. With simulations, we further analyzed how imbalanced allelic expression of non-imprinted genes or different minor allele frequencies affected the predictions of dsPIG. Interestingly, we found that, among biallelically expressed genes, at least 18 genes expressed significantly more transcripts from one allele than the other among different individuals and tissues. Conclusion With the prevalence of the mRNA-Seq technology, dsPIG has become a useful tool for analysis of allelic expression and large-scale prediction of imprinted genes. For ease of use, we have set up a web service and also provided an R package for dsPIG at http://www.shoudanliang.com/dsPIG/.

Abstract Background Neurodegenerative diseases are progressive and irreversible and they can be initiated by mutations in specific genes. Spalt-like genes (Sall) encode transcription factors expressed in the central nervous system. In humans, SALL mutations are associated with hereditary syndromes characterized by mental retardation, sensorineural deafness and motoneuron problems, among others. Drosophila sall mutants exhibit severe neurodegeneration of the central nervous system at embryonic stage 16, which surprisingly reverts later in development at embryonic stage 17, suggesting a potential to recover from neurodegeneration. We hypothesize that this recovery is mediated by a reorganization of the transcriptome counteracting SALL lost. To identify genes associated to neurodegeneration and neuroprotection, we used mRNA-Seq to compare the transcriptome of Drosophila sall mutant and wild type embryos from neurodegeneration and reversal stages. Results Neurodegeneration stage is associated with transcriptional changes in 220 genes, of which only 5% were already described as relevant for neurodegeneration. Genes related to the groups of Redox, Lifespan/Aging and Mitochondrial diseases are significantly represented at this stage. By contrast, neurodegeneration reversal stage is associated with significant changes in 480 genes, including 424 not previously associated with neuroprotection. Immune response and Salt stress are the most represented groups at this stage. Conclusions We identify new genes associated to neurodegeneration and neuroprotection by using an mRNA-Seq approach. The strong homology between Drosophila and human genes raises the possibility to unveil novel genes involved in neurodegeneration and neuroprotection also in humans.

Abstract Background Neonatal cerebellar cortical degeneration is a neurodegenerative disease described in several canine breeds including the Beagle. Affected Beagles are unable to ambulate normally from the onset of walking and the main pathological findings include Purkinje cell loss with swollen dendritic processes. Previous reports suggest an autosomal recessive mode of inheritance. The development of massively parallel sequencing techniques has presented the opportunity to investigate individual clinical cases using genome-wide sequencing approaches. We used genome-wide mRNA sequencing (mRNA-seq) of cerebellum tissue from a single Beagle with neonatal cerebellar cortical degeneration as a method of candidate gene sequencing, with the aim of identifying the causal mutation. Results A four-week old Beagle dog presented with progressive signs of cerebellar ataxia and the owner elected euthanasia. Histopathology revealed findings consistent with cerebellar cortical degeneration. Genome-wide mRNA sequencing (mRNA-seq) of RNA from cerebellum tissue was used as a method of candidate gene sequencing. After analysis of the canine orthologues of human spinocerebellar ataxia associated genes, we identified a homozygous 8 bp deletion in the β-III spectrin gene, SPTBN2, associated with spinocerebellar type 5 in humans. Genotype analysis of the sire, dam, ten clinically unaffected siblings, and an affected sibling from a previous litter, showed the mutation to fully segregate with the disorder. Previous studies have shown that β-III spectrin is critical for Purkinje cell development, and the absence of this protein can lead to cell damage through excitotoxicity, consistent with the observed Purkinje cell loss, degeneration of dendritic processes and associated neurological dysfunction in this Beagle. Conclusions An 8 bp deletion in the SPTBN2 gene encoding β-III spectrin is associated with neonatal cerebellar cortical degeneration in Beagle dogs. This study shows that mRNA-seq is a feasible method of screening candidate genes for mutations associated with rare diseases when a suitable tissue resource is available.

Abstract Background Transcriptome analysis is increasingly being used to study the evolutionary origins and ecology of non-model plants. One issue for both transcriptome assembly and differential gene expression analyses is the common occurrence in plants of hybridisation and whole genome duplication (WGD) and hybridization resulting in allopolyploidy. The divergence of duplicated genes following WGD creates near identical homeologues that can be problematic for de novo assembly and also reference based assembly protocols that use short reads (35 - 100 bp). Results Here we report a successful strategy for the assembly of two transcriptomes made using 75 bp Illumina reads from Pachycladon fastigiatum and Pachycladon cheesemanii. Both are allopolyploid plant species (2n = 20) that originated in the New Zealand Alps about 0.8 million years ago. In a systematic analysis of 19 different coverage cutoffs and 20 different k-mer sizes we showed that i) none of the genes could be assembled across all of the parameter space ii) assembly of each gene required an optimal set of parameter values and iii) these parameter values could be explained in part by different gene expression levels and different degrees of similarity between genes. Conclusions To obtain optimal transcriptome assemblies for allopolyploid plants, k-mer size and k-mer coverage need to be considered simultaneously across a broad parameter space. This is important for assembling a maximum number of full length ESTs and for avoiding chimeric assemblies of homeologous and paralogous gene copies.