Gene set analysis approaches for RNA-seq data: performance evaluation and application guideline

Comparative Analysis of Immune Response Genes Induced by a Virulent or Attenuated Strain of Babesia bigemina

RNA-seq technology has been widely used for the characterization of the transcriptome profile induced by several diseases in both humans and animals. In the present study, RNA-seq was used to identify the differential expression of genes associated with the immune response in cattle infected with two different strains of Babesia bigemina, both derived from the same Mexican field isolate, which exhibit distinct phenotypic characteristics: the virulent strain, capable of producing acute clinical signs, and the attenuated strain, capable of stimulating a protective immune response when used as an immunogen with an efficacy greater than 80%. The differential gene expression analysis performed revealed a total of 620 differentially expressed genes (DEGs). However, the intersection of the edgeR and DESeq2 programs used in the bioinformatics analysis only identified 247 DEGs, of which 108 genes were enriched to be closely correlated with the bovine immune response based on gene ontology terms; most of the DEGs obtained encode proteins associated with the major histocompatibility complex, immunoglobulins, and T-cell surface receptors. The infection caused by the attenuated strain induced higher transcription of immune response genes compared to the infection caused by the virulent strain; nonetheless, in both infections, a greater down-regulation than up-regulation was observed. Different immunoglobulin-associated genes were found to be up-regulated in the group inoculated with the attenuated strain, whereas these were down-regulated in the virulent strain-inoculated group. In addition, an up-regulation of the HSPA6, CD163, and SLC11a1 genes was observed in the group inoculated with the virulent strain, previously reported in other Apicomplexan infections. The findings provide relevant information that could contribute to clarifying the immune response associated with an acute bovine babesiosis infection by B. bigemina.

Improving data interpretability with new differential sample variance gene set tests

Background

Gene set analysis methods have played a major role in generating biological interpretations from omics data such as gene expression datasets. However, most methods focus on detecting homogenous pattern changes in mean expression and methods detecting pattern changes in variance remain poorly explored. While a few studies attempted to use gene-level variance analysis, such approach remains under-utilized. When comparing two phenotypes, gene sets with distinct changes in subgroups under one phenotype are overlooked by available methods although they reflect meaningful biological differences between two phenotypes. Multivariate sample-level variance analysis methods are needed to detect such pattern changes.

Results

We use ranking schemes based on minimum spanning tree to generalize the Cramer-Von Mises and Anderson-Darling univariate statistics into multivariate gene set analysis methods to detect differential sample variance or mean. We characterize these methods in addition to two methods developed earlier using simulation results with different parameters. We apply the developed methods to microarray gene expression dataset of prednisolone-resistant and prednisolone-sensitive children diagnosed with B-lineage acute lymphoblastic leukemia and bulk RNA-sequencing gene expression dataset of benign hyperplastic polyps and potentially malignant sessile serrated adenoma/polyps. One or both of the two compared phenotypes in each of these datasets have distinct molecular subtypes that contribute to heterogeneous differences. Our results show that methods designed to detect differential sample variance are able to detect specific hallmark signaling pathways associated with the two compared phenotypes as documented in available literature.

Conclusions

The results in this study demonstrate the usefulness of methods designed to detect differential sample variance in providing biological interpretations when biologically relevant but heterogeneous changes between two phenotypes are prevalent in specific signaling pathways. Software implementation of the developed methods is available with detailed documentation from Bioconductor package GSAR. The available methods are applicable to gene expression datasets in a normalized matrix form and could be used with other omics datasets in a normalized matrix form with available collection of feature sets.

Benchmarking Algorithms for Gene Set Scoring of Single-cell ATAC-seq Data

Gene set scoring (GSS) has been routinely conducted for gene expression analysis of bulk or single-cell RNA sequencing (RNA-seq) data, which helps to decipher single-cell heterogeneity and cell type-specific variability by incorporating prior knowledge from functional gene sets. Single-cell assay for transposase accessible chromatin using sequencing (scATAC-seq) is a powerful technique for interrogating single-cell chromatin-based gene regulation, and genes or gene sets with dynamic regulatory potentials can be regarded as cell type-specific markers as if in single-cell RNA-seq (scRNA-seq). However, there are few GSS tools specifically designed for scATAC-seq, and the applicability and performance of RNA-seq GSS tools on scATAC-seq data remain to be investigated. Here, we systematically benchmarked ten GSS tools, including four bulk RNA-seq tools, five scRNA-seq tools, and one scATAC-seq method. First, using matched scATAC-seq and scRNA-seq datasets, we found that the performance of GSS tools on scATAC-seq data was comparable to that on scRNA-seq, suggesting their applicability to scATAC-seq. Then, the performance of different GSS tools was extensively evaluated using up to ten scATAC-seq datasets. Moreover, we evaluated the impact of gene activity conversion, dropout imputation, and gene set collections on the results of GSS. Results show that dropout imputation can significantly promote the performance of almost all GSS tools, while the impact of gene activity conversion methods or gene set collections on GSS performance is more dependent on GSS tools or datasets. Finally, we provided practical guidelines for choosing appropriate preprocessing methods and GSS tools in different application scenarios.

Assessment of Gene Set Enrichment Analysis using curated RNA-seq-based benchmarks

Pathway enrichment analysis is a ubiquitous computational biology method to interpret a list of genes (typically derived from the association of large-scale omics data with phenotypes of interest) in terms of higher-level, predefined gene sets that share biological function, chromosomal location, or other common features. Among many tools developed so far, Gene Set Enrichment Analysis (GSEA) stands out as one of the pioneering and most widely used methods. Although originally developed for microarray data, GSEA is nowadays extensively utilized for RNA-seq data analysis. Here, we quantitatively assessed the performance of a variety of GSEA modalities and provide guidance in the practical use of GSEA in RNA-seq experiments. We leveraged harmonized RNA-seq datasets available from The Cancer Genome Atlas (TCGA) in combination with large, curated pathway collections from the Molecular Signatures Database to obtain cancer-type-specific target pathway lists across multiple cancer types. We carried out a detailed analysis of GSEA performance using both gene-set and phenotype permutations combined with four different choices for the Kolmogorov-Smirnov enrichment statistic. Based on our benchmarks, we conclude that the classic/unweighted gene-set permutation approach offered comparable or better sensitivity-vs-specificity tradeoffs across cancer types compared with other, more complex and computationally intensive permutation methods. Finally, we analyzed other large cohorts for thyroid cancer and hepatocellular carcinoma. We utilized a new consensus metric, the Enrichment Evidence Score (EES), which showed a remarkable agreement between pathways identified in TCGA and those from other sources, despite differences in cancer etiology. This finding suggests an EES-based strategy to identify a core set of pathways that may be complemented by an expanded set of pathways for downstream exploratory analysis. This work fills the existing gap in current guidelines and benchmarks for the use of GSEA with RNA-seq data and provides a framework to enable detailed benchmarking of other RNA-seq-based pathway analysis tools.

Integrative bioinformatics and experimental validation of hub genetic markers in acne vulgaris: Toward personalized diagnostic and therapeutic strategies

BACKGROUND

Acne vulgaris is a widespread chronic inflammatory dermatological condition. The precise molecular and genetic mechanisms of its pathogenesis remain incompletely understood. This research synthesizes existing databases, targeting a comprehensive exploration of core genetic markers.

METHODS

Gene expression datasets (GSE6475, GSE108110, and GSE53795) were retrieved from the GEO. Differentially expressed genes (DEGs) were identified using the limma package. Enrichment analyses were conducted using GSVA for pathway assessment and clusterProfiler for GO and KEGG analyses. PPI networks and immune cell infiltration were analyzed using the STRING database and ssGSEA, respectively. We investigated the correlation between hub gene biomarkers and immune cell infiltration using Spearman's rank analysis. ROC curve analysis validated the hub genes' diagnostic accuracy. miRNet, TarBase v8.0, and ChEA3 identified miRNA/transcription factor-gene interactions, while DrugBank delineated drug-gene interactions. Experiments utilized HaCaT cells stimulated with Propionibacterium acnes, treated with retinoic acid and methotrexate, and evaluated using RT-qPCR, ELISA, western blot, lentiviral transduction, CCK-8, wound-healing, and transwell assays.

RESULTS

There were 104 genes with consistent differences across the three datasets of paired acne and normal skin. Functional analyses emphasized the significant enrichment of these DEGs in immune-related pathways. PPI network analysis pinpointed hub genes PTPRC, CXCL8, ITGB2, and MMP9 as central players in acne pathogenesis. Elevated levels of specific immune cell infiltration in acne lesions corroborated the inflammatory nature of the disease. ROC curve analysis identified the acne diagnostic potential of four hub genes. Key miRNAs, particularly hsa-mir-124-3p, and central transcription factors like TFEC were noted as significant regulators. In vitro validation using HaCaT cells confirmed the upregulation of hub genes following Propionibacterium acnes exposure, while CXCL8 knockdown reduced pro-inflammatory cytokines, cell proliferation, and migration. DrugBank insights led to the exploration of retinoic acid and methotrexate, both of which mitigated gene expression upsurge and inflammatory mediator secretion.

CONCLUSION

This comprehensive study elucidated pivotal genes associated with acne pathogenesis, notably PTPRC, CXCL8, ITGB2, and MMP9. The findings underscore potential biomarkers, therapeutic targets, and the therapeutic potential of agents like retinoic acid and methotrexate. The congruence between bioinformatics and experimental validations suggests promising avenues for personalized acne treatments.

Integrative analysis of single-cell and bulk transcriptome data reveal the significant role of macrophages in lupus nephritis

OBJECTIVE

We attempted to identify abnormal immune cell components and signaling pathways in lupus nephritis (LN) and to identify potential therapeutic targets.

METHODS

Differentially expressed genes (DEGs) between LN and normal kidney tissues were identified from bulk transcriptome data, and functional annotation was performed. The phenotypic changes in macrophages and aberrant intercellular signaling communications within immune cells were imputed from LN scRNA-seq data using trajectory analysis and verified using immunofluorescence staining. Finally, lentivirus-mediated overexpression of LGALS9, the gene encoding Galectin 9, in THP-1 cells was used to study the functional effect of this gene on monocytic cells.

RESULTS

From bulk transcriptome data, a significant activation of interferon (IFN) signaling was observed, and its intensity showed a significantly positive correlation with the abundance of infiltrating macrophages in LN. Analysis of scRNA-seq data revealed 17 immune cell clusters, with macrophages showing the highest enrichment of intercellular signal communication in LN. Trajectory analysis revealed macrophages in LN undergo a phenotypic change from inflammatory patrolling macrophages to phagocytic and then to antigen-presenting macrophages, and secrete various pro-inflammatory factors and complement components. LGALS9 was found significantly upregulated in macrophages in LN, which was confirmed by the immunofluorescence assay. Gene functional study showed that LGALS9 overexpression in THP-1 cells significantly elicited pro-inflammatory activation, releasing multiple immune cell chemoattractants.

CONCLUSION

Our results present an important pathophysiological role for macrophages in LN, and our preliminary results demonstrate significant pro-inflammatory effects of LGALS9 gene in LN macrophages.

Transcriptome Analysis of Compensatory Growth and Meat Quality Alteration after Varied Restricted Feeding Conditions in Beef Cattle

Compensatory growth (CG) is a physiological response that accelerates growth following a period of nutrient limitation, with the potential to improve growth efficiency and meat quality in cattle. However, the underlying molecular mechanisms remain poorly understood. In this study, 60 Huaxi cattle were divided into one ad libitum feeding (ALF) group and two restricted feeding groups (75% restricted, RF75; 50% restricted, RF50) undergoing a short-term restriction period followed by evaluation of CG. Detailed comparisons of growth performance during the experimental period, as well as carcass and meat quality traits, were conducted, complemented by a comprehensive transcriptome analysis of the longissimus dorsi muscle using differential expression analysis, gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), and weighted correlation network analysis (WGCNA). The results showed that irrespective of the restriction degree, the restricted animals exhibited CG, achieving final body weights comparable to the ALF group. Compensating animals showed differences in meat quality traits, such as pH, cooking loss, and fat content, compared to the ALF group. Transcriptomic analysis revealed 57 genes and 31 pathways differentially regulated during CG, covering immune response, acid-lipid metabolism, and protein synthesis. Notably, complement-coagulation-fibrinolytic system synergy was identified as potentially responsible for meat quality optimization in RF75. This study provides novel and valuable genetic insights into the regulatory mechanisms of CG in beef cattle.

Roastgsa: a comparison of rotation-based scores for gene set enrichment analysis

BACKGROUND

Gene-wise differential expression is usually the first major step in the statistical analysis of high-throughput data obtained from techniques such as microarrays or RNA-sequencing. The analysis at gene level is often complemented by interrogating the data in a broader biological context that considers as unit of measure groups of genes that may have a common function or biological trait. Among the vast number of publications about gene set analysis (GSA), the rotation test for gene set analysis, also referred to as roast, is a general sample randomization approach that maintains the integrity of the intra-gene set correlation structure in defining the null distribution of the test.

RESULTS

We present roastgsa, an R package that contains several enrichment score functions that feed the roast algorithm for hypothesis testing. These implemented methods are evaluated using both simulated and benchmarking data in microarray and RNA-seq datasets. We find that computationally intensive measures based on Kolmogorov-Smirnov (KS) statistics fail to improve the rates of simpler measures of GSA like mean and maxmean scores. We also show the importance of accounting for the gene linear dependence structure of the testing set, which is linked to the loss of effective signature size. Complete graphical representation of the results, including an approximation for the effective signature size, can be obtained as part of the roastgsa output.

CONCLUSIONS

We encourage the usage of the absmean (non-directional), mean (directional) and maxmean (directional) scores for roast GSA analysis as these are simple measures of enrichment that have presented dominant results in all provided analyses in comparison to the more complex KS measures.

Integrative pathway and network analysis provide insights on flooding-tolerance genes in soybean

Soybean is highly sensitive to flooding and extreme rainfall. The phenotypic variation of flooding tolerance is a complex quantitative trait controlled by many genes and their interaction with environmental factors. We previously constructed a gene-pool relevant to soybean flooding-tolerant responses from integrated multiple omics and non-omics databases, and selected 144 prioritized flooding tolerance genes (FTgenes). In this study, we proposed a comprehensive framework at the systems level, using competitive (hypergeometric test) and self-contained (sum-statistic, sum-square-statistic) pathway-based approaches to identify biologically enriched pathways through evaluating the joint effects of the FTgenes within annotated pathways. These FTgenes were significantly enriched in 36 pathways in the Gene Ontology database. These pathways were related to plant hormones, defense-related, primary metabolic process, and system development pathways, which plays key roles in soybean flooding-induced responses. We further identified nine key FTgenes from important subnetworks extracted from several gene networks of enriched pathways. The nine key FTgenes were significantly expressed in soybean root under flooding stress in a qRT-PCR analysis. We demonstrated that this systems biology framework is promising to uncover important key genes underlying the molecular mechanisms of flooding-tolerant responses in soybean. This result supplied a good foundation for gene function analysis in further work.

Comprehensive investigation of pathway enrichment methods for functional interpretation of LC-MS global metabolomics data

BACKGROUND

Global or untargeted metabolomics is widely used to comprehensively investigate metabolic profiles under various pathophysiological conditions such as inflammations, infections, responses to exposures or interactions with microbial communities. However, biological interpretation of global metabolomics data remains a daunting task. Recent years have seen growing applications of pathway enrichment analysis based on putative annotations of liquid chromatography coupled with mass spectrometry (LC-MS) peaks for functional interpretation of LC-MS-based global metabolomics data. However, due to intricate peak-metabolite and metabolite-pathway relationships, considerable variations are observed among results obtained using different approaches. There is an urgent need to benchmark these approaches to inform the best practices.

RESULTS

We have conducted a benchmark study of common peak annotation approaches and pathway enrichment methods in current metabolomics studies. Representative approaches, including three peak annotation methods and four enrichment methods, were selected and benchmarked under different scenarios. Based on the results, we have provided a set of recommendations regarding peak annotation, ranking metrics and feature selection. The overall better performance was obtained for the mummichog approach. We have observed that a ~30% annotation rate is sufficient to achieve high recall (~90% based on mummichog), and using semi-annotated data improves functional interpretation. Based on the current platforms and enrichment methods, we further propose an identifiability index to indicate the possibility of a pathway being reliably identified. Finally, we evaluated all methods using 11 COVID-19 and 8 inflammatory bowel diseases (IBD) global metabolomics datasets.

The hitchhikers' guide to RNA sequencing and functional analysis

DNA and RNA sequencing technologies have revolutionized biology and biomedical sciences, sequencing full genomes and transcriptomes at very high speeds and reasonably low costs. RNA sequencing (RNA-Seq) enables transcript identification and quantification, but once sequencing has concluded researchers can be easily overwhelmed with questions such as how to go from raw data to differential expression (DE), pathway analysis and interpretation. Several pipelines and procedures have been developed to this effect. Even though there is no unique way to perform RNA-Seq analysis, it usually follows these steps: 1) raw reads quality check, 2) alignment of reads to a reference genome, 3) aligned reads' summarization according to an annotation file, 4) DE analysis and 5) gene set analysis and/or functional enrichment analysis. Each step requires researchers to make decisions, and the wide variety of options and resulting large volumes of data often lead to interpretation challenges. There also seems to be insufficient guidance on how best to obtain relevant information and derive actionable knowledge from transcription experiments. In this paper, we explain RNA-Seq steps in detail and outline differences and similarities of different popular options, as well as advantages and disadvantages. We also discuss non-coding RNA analysis, multi-omics, meta-transcriptomics and the use of artificial intelligence methods complementing the arsenal of tools available to researchers. Lastly, we perform a complete analysis from raw reads to DE and functional enrichment analysis, visually illustrating how results are not absolute truths and how algorithmic decisions can greatly impact results and interpretation.

Synergy between oligosaccharides and probiotics: From metabolic properties to beneficial effects

Synbiotic is defined as the dietary mixture that comprises both probiotic microorganisms and prebiotic substrates. The concept has been steadily gaining attention owing to the rising recognition of probiotic, prebiotics, and gut health. Among prebiotic substances, oligosaccharides demonstrated considerable health beneficial effects in varieties of food products and their combination with probiotics have been subjected to full range of evaluations. This review delineated the landscape of studies using microbial cultures, cell lines, animal model, and human subjects to explore the functional properties and host impacts of these combinations. Overall, the results suggested that these combinations possess respective metabolic properties that could facilitate beneficial activities therefore could be employed as dietary interventions for human health improvement and therapeutic purposes. However, uncertainties, such as applicational practicalities, underutilized analytical tools, contradictory results in studies, unclear mechanisms, and legislation hurdles, still challenges the broad utilization of these combinations. Future studies to address these issues may not only advance current knowledge on probiotic-prebiotic-host interrelationship but also promote respective applications in food and nutrition.

Topological Distribution of Wound Stiffness Modulates Wound-Induced Hair Follicle Neogenesis

In the large full-thickness mouse skin regeneration model, wound-induced hair neogenesis (WIHN) occurs in the wound center. This implies a spatial regulation of hair regeneration. The role of mechanotransduction during tissue regeneration is poorly understood. Here, we created wounds with equal area but different shapes to understand if perturbing mechanical forces change the area and quantity of de novo hair regeneration. Atomic force microscopy of wound stiffness demonstrated a stiffness gradient across the wound with the wound center softer than the margin. Reducing mechanotransduction signals using FAK or myosin II inhibitors significantly increased WIHN and, conversely, enhancing these signals with an actin stabilizer reduced WIHN. Here, α-SMA was downregulated in FAK inhibitor-treated wounds and lowered wound stiffness. Wound center epithelial cells exhibited a spherical morphology relative to wound margin cells. Differential gene expression analysis of FAK inhibitor-treated wound RNAseq data showed that cytoskeleton-, integrin-, and matrix-associated genes were downregulated, while hair follicular neogenesis, cell proliferation, and cell signaling genes were upregulated. Immunohistochemistry staining showed that FAK inhibition increased pSTAT3 nuclear staining in the regenerative wound center, implying enhanced signaling for hair follicular neogenesis. These findings suggest that controlling wound stiffness modulates tissue regeneration encompassing epithelial competence, tissue patterning, and regeneration during wound healing.

enrichMiR predicts functionally relevant microRNAs based on target collections

MicroRNAs (miRNAs) are small non-coding RNAs that are among the main post-transcriptional regulators of gene expression. A number of data collections and prediction tools have gathered putative or confirmed targets of these regulators. It is often useful, for discovery and validation, to harness such collections to perform target enrichment analysis in given transcriptional signatures or gene-sets in order to predict involved miRNAs. While several methods have been proposed to this end, a flexible and user-friendly interface for such analyses using various approaches and collections is lacking. enrichMiR (https://ethz-ins.org/enrichMiR/) addresses this gap by enabling users to perform a series of enrichment tests, based on several target collections, to rank miRNAs according to their likely involvement in the control of a given transcriptional signature or gene-set. enrichMiR results can furthermore be visualised through interactive and publication-ready plots. To guide the choice of the appropriate analysis method, we benchmarked various tests across a panel of experiments involving the perturbation of known miRNAs. Finally, we showcase enrichMiR functionalities in a pair of use cases.

On the influence of several factors on pathway enrichment analysis

Pathway enrichment analysis has become a widely used knowledge-based approach for the interpretation of biomedical data. Its popularity has led to an explosion of both enrichment methods and pathway databases. While the elegance of pathway enrichment lies in its simplicity, multiple factors can impact the results of such an analysis, which may not be accounted for. Researchers may fail to give influential aspects their due, resorting instead to popular methods and gene set collections, or default settings. Despite ongoing efforts to establish set guidelines, meaningful results are still hampered by a lack of consensus or gold standards around how enrichment analysis should be conducted. Nonetheless, such concerns have prompted a series of benchmark studies specifically focused on evaluating the influence of various factors on pathway enrichment results. In this review, we organize and summarize the findings of these benchmarks to provide a comprehensive overview on the influence of these factors. Our work covers a broad spectrum of factors, spanning from methodological assumptions to those related to prior biological knowledge, such as pathway definitions and database choice. In doing so, we aim to shed light on how these aspects can lead to insignificant, uninteresting or even contradictory results. Finally, we conclude the review by proposing future benchmarks as well as solutions to overcome some of the challenges, which originate from the outlined factors.

Phytochrome B and PHYTOCHROME INTERACTING FACTOR8 modulate seasonal growth in trees

The seasonally synchronized annual growth cycle that is regulated mainly by photoperiod and temperature cues is a crucial adaptive strategy for perennial plants in boreal and temperate ecosystems. Phytochrome B (phyB), as a light and thermal sensor, has been extensively studied in Arabidopsis. However, the specific mechanisms for how the phytochrome photoreceptors control the phenology in tree species remain poorly understood. We characterized the functions of PHYB genes and their downstream PHYTOCHROME INTERACTING FACTOR (PIF) targets in the regulation of shade avoidance and seasonal growth in hybrid aspen trees. We show that while phyB1 and phyB2, as phyB in other plants, act as suppressors of shoot elongation during vegetative growth, they act as promoters of tree seasonal growth. Furthermore, while the Populus homologs of both PIF4 and PIF8 are involved in the shade avoidance syndrome (SAS), only PIF8 plays a major role as a suppressor of seasonal growth. Our data suggest that the PHYB-PIF8 regulon controls seasonal growth through the regulation of FT and CENL1 expression while a genome-wide transcriptome analysis suggests how, in Populus trees, phyB coordinately regulates SAS responses and seasonal growth cessation.

Statistical Approach of Gene Set Analysis with Quantitative Trait Loci for Crop Gene Expression Studies

Genome-wide expression study is a powerful genomic technology to quantify expression dynamics of genes in a genome. In gene expression study, gene set analysis has become the first choice to gain insights into the underlying biology of diseases or stresses in plants. It also reduces the complexity of statistical analysis and enhances the explanatory power of the obtained results from the primary downstream differential expression analysis. The gene set analysis approaches are well developed in microarrays and RNA-seq gene expression data analysis. These approaches mainly focus on analyzing the gene sets with gene ontology or pathway annotation data. However, in plant biology, such methods may not establish any formal relationship between the genotypes and the phenotypes, as most of the traits are quantitative and controlled by polygenes. The existing Quantitative Trait Loci (QTL)-based gene set analysis approaches only focus on the over-representation analysis of the selected genes while ignoring their associated gene scores. Therefore, we developed an innovative statistical approach, GSQSeq, to analyze the gene sets with trait enriched QTL data. This approach considers the associated differential expression scores of genes while analyzing the gene sets. The performance of the developed method was tested on five different crop gene expression datasets obtained from real crop gene expression studies. Our analytical results indicated that the trait-specific analysis of gene sets was more robust and successful through the proposed approach than existing techniques. Further, the developed method provides a valuable platform for integrating the gene expression data with QTL data.

From Transcriptomics to Treatment in Inherited Optic Neuropathies

Inherited optic neuropathies, including Leber Hereditary Optic Neuropathy (LHON) and Dominant Optic Atrophy (DOA), are monogenetic diseases with a final common pathway of mitochondrial dysfunction leading to retinal ganglion cell (RGC) death and ultimately loss of vision. They are, therefore, excellent models with which to investigate this ubiquitous disease process-implicated in both common polygenetic ocular diseases (e.g., Glaucoma) and late-onset central nervous system neurodegenerative diseases (e.g., Parkinson disease). In recent years, cellular and animal models of LHON and DOA have matured in parallel with techniques (such as RNA-seq) to determine and analyze the transcriptomes of affected cells. This confluence leaves us at a particularly exciting time with the potential for the identification of novel pathogenic players and therapeutic targets. Here, we present a discussion of the importance of inherited optic neuropathies and how transcriptomic techniques can be exploited in the development of novel mutation-independent, neuroprotective therapies.

Development of a clinically relevant ovarian cancer model incorporating surgical cytoreduction to evaluate treatment of micro-metastatic disease

OBJECTIVES

Mouse models of ovarian cancer commonly transfer large numbers of tumor cells into the peritoneal cavity to establish experimental metastatic disease, which may not adequately model early metastatic spread from a primary tumor site. We hypothesized we could develop an ovarian cancer model that predictably represents micro-metastatic disease.

METHODS

Murine ID8_VEGF ovarian cancer cells were transduced to express enhanced luciferase (eLuc) to enable intravital detection of microscopic disease burden and injected beneath the ovarian bursa of C57Bl/6 mice. At 6 or 10 weeks after orthotopic injection, when mice had detectable metastases, hysterectomy and bilateral salpingo-oophorectomy was performed to remove all macroscopic disease, and survival monitored. Immunohistochemistry and gene expression profiling were performed on primary and metastatic tumors.

RESULTS

eLuc-transduced ID8_VEGF cells were brighter than cells transduced with standard luciferase, enabling in vivo visualization of microscopic intra-abdominal metastases developing after orthotopic injection. Primary surgical cytoreduction removed the primary tumor mass but left minimal residual disease in all mice. Metastatic sites that developed following orthotopic injection were similar to metastatic human ovarian cancer sites. Gene expression and immune infiltration were similar between primary and metastatic mouse tumors. Surgical cytoreduction prolonged survival compared to no surgery, with earlier cytoreduction more beneficial than delayed, despite micro-metastatic disease in both settings.

CONCLUSIONS

Mice with primary ovarian tumors established through orthotopic injection develop progressively fatal metastatic ovarian cancer, and benefit from surgical cytoreduction to remove bulky disease. This model enables the analysis of therapeutic regimens designed to target and potentially eradicate established minimal residual disease.

Discovery of Selenocysteine as a Potential Nanomedicine Promotes Cartilage Regeneration With Enhanced Immune Response by Text Mining and Biomedical Databases

Background

Unlike bone tissue, little progress has been made regarding cartilage regeneration, and many challenges remain. Furthermore, the key roles of cartilage lesion caused by traumas, focal lesion, or articular overstress remain unclear. Traumatic injuries to the meniscus as well as its degeneration are important risk factors for long-term joint dysfunction, degenerative joint lesions, and knee osteoarthritis (OA) a chronic joint disease characterized by degeneration of articular cartilage and hyperosteogeny. Nearly 50% of the individuals with meniscus injuries develop OA over time. Due to the limited inherent self-repair capacity of cartilage lesion, the Biomaterial drug-nanomedicine is considered to be a promising alternative. Therefore, it is important to elucidate the gene potential regeneration mechanisms and discover novel precise medication, which are identified through this study to investigate their function and role in pathogenesis.

Methods

We downloaded the mRNA microarray statistics GSE117999, involving paired cartilage lesion tissue samples from 12 OA patients and 12 patients from a control group. First, we analyzed these statistics to recognize the differentially expressed genes (DEGs). We then exposed the gene ontology (GO) annotation and the Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway enrichment analyses for these DEGs. Protein-protein interaction (PPI) networks were then constructed, from which we attained eight significant genes after a functional interaction analysis. Finally, we identified a potential nanomedicine attained from this assay set, using a wide range of inhibitor information archived in the Search Tool for the Retrieval of Interacting Genes (STRING) database.

Results

Sixty-six DEGs were identified with our standards for meaning (adjusted P-value < 0.01, |log2 - FC| ≥1.2). Furthermore, we identified eight hub genes and one potential nanomedicine - Selenocysteine based on these integrative data.

Conclusion

We identified eight hub genes that could work as prospective biomarkers for the diagnostic and biomaterial drug treatment of cartilage lesion, involving the novel genes CAMP, DEFA3, TOLLIP, HLA-DQA2, SLC38A6, SLC3A1, FAM20A, and ANO8. Meanwhile, these genes were mainly associated with immune response, immune mediator induction, and cell chemotaxis. Significant support is provided for obtaining a series of novel gene targets, and we identify potential mechanisms for cartilage regeneration and final nanomedicine immunotherapy in regenerative medicine.

Gene Set Analysis: Challenges, Opportunities, and Future Research

Gene set analysis methods are widely used to provide insight into high-throughput gene expression data. There are many gene set analysis methods available. These methods rely on various assumptions and have different requirements, strengths and weaknesses. In this paper, we classify gene set analysis methods based on their components, describe the underlying requirements and assumptions for each class, and provide directions for future research in developing and evaluating gene set analysis methods.

Transcriptomics in Toxicogenomics, Part II: Preprocessing and Differential Expression Analysis for High Quality Data

Preprocessing of transcriptomics data plays a pivotal role in the development of toxicogenomics-driven tools for chemical toxicity assessment. The generation and exploitation of large volumes of molecular profiles, following an appropriate experimental design, allows the employment of toxicogenomics (TGx) approaches for a thorough characterisation of the mechanism of action (MOA) of different compounds. To date, a plethora of data preprocessing methodologies have been suggested. However, in most cases, building the optimal analytical workflow is not straightforward. A careful selection of the right tools must be carried out, since it will affect the downstream analyses and modelling approaches. Transcriptomics data preprocessing spans across multiple steps such as quality check, filtering, normalization, batch effect detection and correction. Currently, there is a lack of standard guidelines for data preprocessing in the TGx field. Defining the optimal tools and procedures to be employed in the transcriptomics data preprocessing will lead to the generation of homogeneous and unbiased data, allowing the development of more reliable, robust and accurate predictive models. In this review, we outline methods for the preprocessing of three main transcriptomic technologies including microarray, bulk RNA-Sequencing (RNA-Seq), and single cell RNA-Sequencing (scRNA-Seq). Moreover, we discuss the most common methods for the identification of differentially expressed genes and to perform a functional enrichment analysis. This review is the second part of a three-article series on Transcriptomics in Toxicogenomics.

Fifteen Years of Gene Set Analysis for High-Throughput Genomic Data: A Review of Statistical Approaches and Future Challenges

Over the last decade, gene set analysis has become the first choice for gaining insights into underlying complex biology of diseases through gene expression and gene association studies. It also reduces the complexity of statistical analysis and enhances the explanatory power of the obtained results. Although gene set analysis approaches are extensively used in gene expression and genome wide association data analysis, the statistical structure and steps common to these approaches have not yet been comprehensively discussed, which limits their utility. In this article, we provide a comprehensive overview, statistical structure and steps of gene set analysis approaches used for microarrays, RNA-sequencing and genome wide association data analysis. Further, we also classify the gene set analysis approaches and tools by the type of genomic study, null hypothesis, sampling model and nature of the test statistic, etc. Rather than reviewing the gene set analysis approaches individually, we provide the generation-wise evolution of such approaches for microarrays, RNA-sequencing and genome wide association studies and discuss their relative merits and limitations. Here, we identify the key biological and statistical challenges in current gene set analysis, which will be addressed by statisticians and biologists collectively in order to develop the next generation of gene set analysis approaches. Further, this study will serve as a catalog and provide guidelines to genome researchers and experimental biologists for choosing the proper gene set analysis approach based on several factors.

Identification of altered biological processes in heterogeneous RNA-sequencing data by discretization of expression profiles

Heterogeneity is a fundamental feature of complex phenotypes. So far, genomic screenings have profiled thousands of samples providing insights into the transcriptome of the cell. However, disentangling the heterogeneity of these transcriptomic Big Data to identify defective biological processes remains challenging. Here we present GSECA, a method exploiting the bimodal behavior of RNA-sequencing gene expression profiles to identify altered gene sets in heterogeneous patient cohorts. Using simulated and experimental RNA-sequencing data sets, we show that GSECA provides higher performances than other available algorithms in detecting truly altered biological processes in large cohorts. Applied to 5941 samples from 14 different cancer types, GSECA correctly identified the alteration of the PI3K/AKT signaling pathway driven by the somatic loss of PTEN and verified the emerging role of PTEN in modulating immune-related processes. In particular, we showed that, in prostate cancer, PTEN loss appears to establish an immunosuppressive tumor microenvironment through the activation of STAT3, and low PTEN expression levels have a detrimental impact on patient disease-free survival. GSECA is available at https://github.com/matteocereda/GSECA.

Housekeeping gene validation for RT-qPCR studies on synovial fibroblasts derived from healthy and osteoarthritic patients with focus on mechanical loading

Selection of appropriate housekeeping genes is essential for the validity of data normalization in reverse transcription quantitative PCR (RT-qPCR). Synovial fibroblasts (SF) play a mediating role in the development and progression of osteoarthritis (OA) pathogenesis, but there is no information on reliable housekeeping genes available. Therefore the goal of this study was to identify a set of reliable housekeeping genes suitable for studies of mechanical loading on SF from healthy and OA patients. Nine genes were evaluated towards expression stability and ranked according their relative stability determined by four different mathematical procedures (geNorm, NormFinder, BestKeeper and comparative ΔC_q). We observed that RPLP0 (ribosomal protein, large, P0) and EEF1A1 (eukaryotic translation elongation factor 1 alpha 1) turned out to be the genes with the most stable expression in SF from non-OA or OA patients treated with or without mechanical loading. According to geNorm two genes are sufficient for normalization throughout. Expression of one tested target gene varied considerably, if normalized to different candidate housekeeping genes. Our study provides a tool for accurate and valid housekeeping gene selection in gene expression experiments on SF from healthy and OA patients with and without mechanical loading in consistent with the MIQE (Minimum Information for Publication of Quantitative Real-Time PCR Experiments) guidelines and additionally demonstrates the impact of proper housekeeping gene selection on the expression of the gene of interest.

Longitudinal linear combination test for gene set analysis

Background

Although microarray studies have greatly contributed to recent genetic advances, lack of replication has been a continuing concern in this area. Complex study designs have the potential to address this concern, though they remain undervalued by investigators due to the lack of proper analysis methods. The primary challenge in the analysis of complex microarray study data is handling the correlation structure within data while also dealing with the combination of large number of genetic measurements and small number of subjects that are ubiquitous even in standard microarray studies. Motivated by the lack of available methods for analysis of repeatedly measured phenotypic or transcriptomic data, herein we develop a longitudinal linear combination test (LLCT).

Results

LLCT is a two-step method to analyze multiple longitudinal phenotypes when there is high dimensionality in response and/or explanatory variables. Alternating between calculating within-subjects and between-subjects variations in two steps, LLCT examines if the maximum possible correlation between a linear combination of the time trends and a linear combination of the predictors given by the gene expressions is statistically significant. A generalization of this method can handle family-based study designs when the subjects are not independent. This method is also applicable to time-course microarray, with the ability to identify gene sets that exhibit significantly different expression patterns over time. Based on the results from a simulation study, LLCT outperformed its alternative: pathway analysis via regression. LLCT was shown to be very powerful in the analysis of large gene sets even when the sample size is small.

Conclusions

This self-contained pathway analysis method is applicable to a wide range of longitudinal genomics, proteomics, metabolomics (OMICS) data, allows adjusting for potentially time-dependent covariates and works well with unbalanced and incomplete data. An important potential application of this method could be time-course linkage of OMICS, an attractive possibility for future genetic researchers.

Availability: R package of LLCT is available at: https://github.com/its-likeli-jeff/LLCT

Proteome-transcriptome alignment of molecular portraits achieved by self-contained gene set analysis: Consensus colon cancer subtypes case study

Gene set analysis (GSA) has become the common methodology for analyzing transcriptomics data. However, self-contained GSA techniques are rarely, if ever, used for proteomics data analysis. Here we present a self-contained proteome level GSA of four consensus molecular subtypes (CMSs) previously established by transcriptome dissection of colon carcinoma specimens. Despite notable difference in structure of proteomics and transcriptomics data, many pathway-wide characteristic features of CMSs found at the mRNA level were reproduced at the protein level. In particular, CMS1 features show heavy involvement of immune system as well as the pathways related to mismatch repair, DNA replication and functioning of proteasome, while CMS4 tumors upregulate complement pathway and proteins participating in epithelial-to-mesenchymal transition (EMT). In addition, protein level GSA yielded a set of novel observations visible at the proteome, but not at the transcriptome level, including possible involvement of major histocompatibility complex II (MHC-II) antigens in the known immunogenicity of CMS1 and a connection between cholesterol trafficking and the regulation of Integrin-linked kinase (ILK) in CMS3. Overall, this study proves utility of self-contained GSA approaches as a critical tool for analyzing proteomics data in general and dissecting protein-level molecular portraits of human tumors in particular.

Construction of subtype-specific prognostic gene signatures for early-stage non-small cell lung cancer using meta feature selection methods

Feature selection in the framework of meta-analyses (meta feature selection), combines meta-analysis with a feature selection process and thus allows meta-analysis feature selection across multiple datasets. In the present study, a meta feature selection procedure that fitted a multiple Cox regression model to estimate the effect size of a gene in individual studies and to identify the overall effect of the gene using a meta-analysis model was proposed. The method was used to identify prognostic gene signatures for lung adenocarcinoma and lung squamous cell carcinoma. Furthermore, redundant gene elimination (RGE) is of crucial importance during feature selection, and is also essential for a meta feature selection process. The current study demonstrated that the proposed meta feature selection procedure with RGE outperforms that without RGE in terms of predictive ability, model parsimony and biological interpretation.

Engineered Adoptive T-cell Therapy Prolongs Survival in a Preclinical Model of Advanced-Stage Ovarian Cancer

Adoptive T-cell therapy using high-affinity T-cell receptors (TCR) to target tumor antigens has potential for improving outcomes in high-grade serous ovarian cancer (HGSOC) patients. Ovarian tumors develop a hostile, multicomponent tumor microenvironment containing suppressive cells, inhibitory ligands, and soluble factors that facilitate evasion of antitumor immune responses. Developing and validating an immunocompetent mouse model of metastatic ovarian cancer that shares antigenic and immunosuppressive qualities of human disease would facilitate establishing effective T-cell therapies. We used deep transcriptome profiling and IHC analysis of human HGSOC tumors and disseminated mouse ID8_VEGF tumors to compare immunologic features. We then evaluated the ability of CD8 T cells engineered to express a high-affinity TCR specific for mesothelin, an ovarian cancer antigen, to infiltrate advanced ID8_VEGF murine ovarian tumors and control tumor growth. Human CD8 T cells engineered to target mesothelin were also evaluated for ability to kill HLA-A2⁺ HGSOC lines. IHC and gene-expression profiling revealed striking similarities between tumors of both species, including processing/presentation of a leading candidate target antigen, suppressive immune cell infiltration, and expression of molecules that inhibit T-cell function. Engineered T cells targeting mesothelin infiltrated mouse tumors but became progressively dysfunctional and failed to persist. Treatment with repeated doses of T cells maintained functional activity, significantly prolonging survival of mice harboring late-stage disease at treatment onset. Human CD8 T cells engineered to target mesothelin were tumoricidal for three HGSOC lines. Treatment with engineered T cells may have clinical applicability in patients with advanced-stage HGSOC.

Simultaneous Enrichment Analysis of all Possible Gene-sets: Unifying Self-Contained and Competitive Methods

Studying sets of genomic features is increasingly popular in genomics, proteomics and metabolomics since analyzing at set level not only creates a natural connection to biological knowledge but also offers more statistical power. Currently, there are two gene-set testing approaches, self-contained and competitive, both of which have their advantages and disadvantages, but neither offers the final solution. We introduce simultaneous enrichment analysis (SEA), a new approach for analysis of feature sets in genomics and other omics based on a new unified null hypothesis, which includes the self-contained and competitive null hypotheses as special cases. We employ closed testing using Simes tests to test this new hypothesis. For every feature set, the proportion of active features is estimated, and a confidence bound is provided. Also, for every unified null hypotheses, a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$P$\end{document}-value is calculated, which is adjusted for family-wise error rate. SEA does not need to assume that the features are independent. Moreover, users are allowed to choose the feature set(s) of interest after observing the data. We develop a novel pipeline and apply it on RNA-seq data of dystrophin-deficient mdx mice, showcasing the flexibility of the method. Finally, the power properties of the method are evaluated through simulation studies.

A Novel Joint Gene Set Analysis Framework Improves Identification of Enriched Pathways in Cross Disease Transcriptomic Analysis

Motivation: Gene set enrichment analysis is a widely accepted expression analysis tool which aims at detecting coordinated expression change within a pre-defined gene sets rather than individual genes. The benefit of gene set analysis over individual differentially expressed (DE) gene analysis includes more reproducible and interpretable results and detecting small but consistent change among gene set which could not be detected by DE gene analysis. There have been many successful gene set analysis applications in human diseases. However, when the sample size of a disease study is small and no other public data sets of the same disease are available, it will lead to lack of power to detect pathways of importance to the disease. Results: We have developed a novel joint gene set analysis statistical framework which aims at improving the power of identifying enriched gene sets through integrating multiple similar disease data sets. Through comprehensive simulation studies, we demonstrated that our proposed frameworks obtained much better AUC scores than single data set analysis and another meta-analysis method in identification of enriched pathways. When applied to two real data sets, the proposed framework could retain the enriched gene sets identified by single data set analysis and exclusively obtained up to 200% more disease-related gene sets demonstrating the improved identification power through information shared between similar diseases. We expect that the proposed framework would enable researchers to better explore public data sets when the sample size of their study is limited.

Comprehensive and critical evaluation of individualized pathway activity measurement tools on pan-cancer data

Motivation

: Biological pathways are extensively used for the analysis of transcriptome data to characterize biological mechanisms underlying various phenotypes. There are a number of computational tools that summarize transcriptome data at the pathway level. However, there is no comparative study on how well these tools produce useful information at the cohort level, enabling comparison of many samples or patients.

Results

: In this study, we systematically compared and evaluated 13 different pathway activity inference tools based on 5 comparison criteria using pan-cancer data set. This study has two major contributions. First, our study provides a comprehensive survey on computational techniques used by existing pathway activity inference tools. The tools use different strategies and assume different requirements on data: input transformation, use of labels, necessity of cohort-level input data, use of gene relations and scoring metric. Second, we performed extensive evaluations on the performance of these tools. Because different tools use different methods to map samples to the pathway dimension, the tools are evaluated at the pathway level using five comparison criteria. Starting from measuring how well a tool maintains the characteristics of original gene expression values, robustness was also investigated by adding noise into gene expression data. Classification tasks on three clinical variables (tumor versus normal, survival and cancer subtypes) were performed to evaluate the utility of tools for their clinical applications. In addition, the inferred activity values were compared between the tools to see how similar they are along with the scoring schemes they use.

Probabilistic prioritization of candidate pathway association with pathway score

Background

Current methods for gene-set or pathway analysis are usually designed to test the enrichment of a single gene-set. Once the analysis is carried out for each of the sets under study, a list of significant sets can be obtained. However, if one wishes to further prioritize the importance or strength of association of these sets, no such quantitative measure is available. Using the magnitude of p-value to rank the pathways may not be appropriate because p-value is not a measure for strength of significance. In addition, when testing each pathway, these analyses are often implicitly affected by the number of differentially expressed genes included in the set and/or affected by the dependence among genes.

Results

Here we propose a two-stage procedure to prioritize the pathways/gene-sets. In the first stage we develop a pathway-level measure with three properties. First, it contains all genes (differentially expressed or not) in the same set, and summarizes the collective effect of all genes per sample. Second, this pathway score accounts for the correlation between genes by synchronizing their correlation directions. Third, the score includes a rank transformation to enhance the variation among samples as well as to avoid the influence of extreme heterogeneity among genes. In the second stage, all scores are included simultaneously in a Bayesian logistic regression model which can evaluate the strength of association for each set and rank the sets based on posterior probabilities. Simulations from Gaussian distributions and human microarray data, and a breast cancer study with RNA-Seq are considered for demonstration and comparison with other existing methods.

Conclusions

The proposed summary pathway score provides for each sample an overall evaluation of gene expression in a gene-set. It demonstrates the advantages of including all genes in the set and the synchronization of correlation direction. The simultaneous utilization of all pathway-level scores in a Bayesian model not only offers a probabilistic evaluation and ranking of the pathway association but also presents good accuracy in identifying the top-ranking pathways. The resulting recommendation list of ranked pathways can be a reference for potential target therapy or for future allocation of research resources.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2411-z) contains supplementary material, which is available to authorized users.

CGPS: A machine learning-based approach integrating multiple gene set analysis tools for better prioritization of biologically relevant pathways

Gene set enrichment (GSE) analyses play an important role in the interpretation of large-scale transcriptome datasets. Multiple GSE tools can be integrated into a single method as obtaining optimal results is challenging due to the plethora of GSE tools and their discrepant performances. Several existing ensemble methods lead to different scores in sorting pathways as integrated results; furthermore, it is difficult for users to choose a single ensemble score to obtain optimal final results. Here, we develop an ensemble method using a machine learning approach called Combined Gene set analysis incorporating Prioritization and Sensitivity (CGPS) that integrates the results provided by nine prominent GSE tools into a single ensemble score (R score) to sort pathways as integrated results. Moreover, to the best of our knowledge, CGPS is the first GSE ensemble method built based on a priori knowledge of pathways and phenotypes. Compared with 10 widely used individual methods and five types of ensemble scores from two ensemble methods, we demonstrate that sorting pathways based on the R score can better prioritize relevant pathways, as established by an evaluation of 120 simulated datasets and 45 real datasets. Additionally, CGPS is applied to expression data involving the drug panobinostat, which is an anticancer treatment against multiple myeloma. The results identify cell processes associated with cancer, such as the p53 signaling pathway (hsa04115); by contrast, according to two ensemble methods (EnrichmentBrowser and EGSEA), this pathway has a rank higher than 20, which may cause users to miss the pathway in their analyses. We show that this method, which is based on a priori knowledge, can capture valuable biological information from numerous types of gene set collections, such as KEGG pathways, GO terms, Reactome, and BioCarta. CGPS is publicly available as a standalone source code at ftp://ftp.cbi.pku.edu.cn/pub/CGPS_download/cgps-1.0.0.tar.gz.

Computational systems biology approaches for Parkinson's disease

Parkinson's disease (PD) is a prime example of a complex and heterogeneous disorder, characterized by multifaceted and varied motor- and non-motor symptoms and different possible interplays of genetic and environmental risk factors. While investigations of individual PD-causing mutations and risk factors in isolation are providing important insights to improve our understanding of the molecular mechanisms behind PD, there is a growing consensus that a more complete understanding of these mechanisms will require an integrative modeling of multifactorial disease-associated perturbations in molecular networks. Identifying and interpreting the combinatorial effects of multiple PD-associated molecular changes may pave the way towards an earlier and reliable diagnosis and more effective therapeutic interventions. This review provides an overview of computational systems biology approaches developed in recent years to study multifactorial molecular alterations in complex disorders, with a focus on PD research applications. Strengths and weaknesses of different cellular pathway and network analyses, and multivariate machine learning techniques for investigating PD-related omics data are discussed, and strategies proposed to exploit the synergies of multiple biological knowledge and data sources. A final outlook provides an overview of specific challenges and possible next steps for translating systems biology findings in PD to new omics-based diagnostic tools and targeted, drug-based therapeutic approaches.

Variance component score test for time-course gene set analysis of longitudinal RNA-seq data

As gene expression measurement technology is shifting from microarrays to sequencing, the statistical tools available for their analysis must be adapted since RNA-seq data are measured as counts. It has been proposed to model RNA-seq counts as continuous variables using nonparametric regression to account for their inherent heteroscedasticity. In this vein, we propose tcgsaseq, a principled, model-free, and efficient method for detecting longitudinal changes in RNA-seq gene sets defined a priori. The method identifies those gene sets whose expression varies over time, based on an original variance component score test accounting for both covariates and heteroscedasticity without assuming any specific parametric distribution for the (transformed) counts. We demonstrate that despite the presence of a nonparametric component, our test statistic has a simple form and limiting distribution, and both may be computed quickly. A permutation version of the test is additionally proposed for very small sample sizes. Applied to both simulated data and two real datasets, tcgsaseq is shown to exhibit very good statistical properties, with an increase in stability and power when compared to state-of-the-art methods ROAST (rotation gene set testing), edgeR, and DESeq2, which can fail to control the type I error under certain realistic settings. We have made the method available for the community in the R package tcgsaseq.

Extracting the Strongest Signals from Omics Data: Differentially Expressed Pathways and Beyond

The analysis of gene sets (in a form of functionally related genes or pathways) has become the method of choice for extracting the strongest signals from omics data. The motivation behind using gene sets instead of individual genes is two-fold. First, this approach incorporates pre-existing biological knowledge into the analysis and facilitates the interpretation of experimental results. Second, it employs a statistical hypotheses testing framework. Here, we briefly review main Gene Set Analysis (GSA) approaches for testing differential expression of gene sets and several GSA approaches for testing statistical hypotheses beyond differential expression that allow extracting additional biological information from the data. We distinguish three major types of GSA approaches testing: (1) differential expression (DE), (2) differential variability (DV), and (3) differential co-expression (DC) of gene sets between two phenotypes. We also present comparative power analysis and Type I error rates for different approaches in each major type of GSA on simulated data. Our evaluation presents a concise guideline for selecting GSA approaches best performing under particular experimental settings. The value of the three major types of GSA approaches is illustrated with real data example. While being applied to the same data set, major types of GSA approaches result in complementary biological information.

A Novel Unsupervised Algorithm for Biological Process-based Analysis on Cancer

The aberrant alterations of biological functions are well known in tumorigenesis and cancer development. Hence, with advances in high-throughput sequencing technologies, capturing and quantifying the functional alterations in cancers based on expression profiles to explore cancer malignant process is highlighted as one of the important topics among cancer researches. In this article, we propose an algorithm for quantifying biological processes by using gene expression profiles over a sample population, which involves the idea of constructing principal curves to condense information of each biological process by a novel scoring scheme on an individualized manner. After applying our method on several large-scale breast cancer datasets in survival analysis, a subset of these biological processes extracted from corresponding survival model is then found to have significant associations with clinical outcomes. Further analyses of these biological processes enable the study of the interplays between biological processes and cancer phenotypes of interest, provide us valuable insights into cancer biology in biological process level and guide the precision treatment for cancer patients. And notably, prognosis predictions based on our method are consistently superior to the existing state of art methods with the same intention.

A transcriptomics data-driven gene space accurately predicts liver cytopathology and drug-induced liver injury

Predicting unanticipated harmful effects of chemicals and drug molecules is a difficult and costly task. Here we utilize a 'big data compacting and data fusion'-concept to capture diverse adverse outcomes on cellular and organismal levels. The approach generates from transcriptomics data set a 'predictive toxicogenomics space' (PTGS) tool composed of 1,331 genes distributed over 14 overlapping cytotoxicity-related gene space components. Involving ∼2.5 × 10⁸ data points and 1,300 compounds to construct and validate the PTGS, the tool serves to: explain dose-dependent cytotoxicity effects, provide a virtual cytotoxicity probability estimate intrinsic to omics data, predict chemically-induced pathological states in liver resulting from repeated dosing of rats, and furthermore, predict human drug-induced liver injury (DILI) from hepatocyte experiments. Analysing 68 DILI-annotated drugs, the PTGS tool outperforms and complements existing tests, leading to a hereto-unseen level of DILI prediction accuracy.

Combining multiple tools outperforms individual methods in gene set enrichment analyses

Motivation

Gene set enrichment (GSE) analysis allows researchers to efficiently extract biological insight from long lists of differentially expressed genes by interrogating them at a systems level. In recent years, there has been a proliferation of GSE analysis methods and hence it has become increasingly difficult for researchers to select an optimal GSE tool based on their particular dataset. Moreover, the majority of GSE analysis methods do not allow researchers to simultaneously compare gene set level results between multiple experimental conditions.

Results

The ensemble of genes set enrichment analyses (EGSEA) is a method developed for RNA-sequencing data that combines results from twelve algorithms and calculates collective gene set scores to improve the biological relevance of the highest ranked gene sets. EGSEA’s gene set database contains around 25 000 gene sets from sixteen collections. It has multiple visualization capabilities that allow researchers to view gene sets at various levels of granularity. EGSEA has been tested on simulated data and on a number of human and mouse datasets and, based on biologists’ feedback, consistently outperforms the individual tools that have been combined. Our evaluation demonstrates the superiority of the ensemble approach for GSE analysis, and its utility to effectively and efficiently extrapolate biological functions and potential involvement in disease processes from lists of differentially regulated genes.

Availability and Implementation

EGSEA is available as an R package at http://www.bioconductor.org/packages/EGSEA/. The gene sets collections are available in the R package EGSEAdata from http://www.bioconductor.org/packages/EGSEAdata/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Altered expression of histamine signaling genes in autism spectrum disorder

The histaminergic system (HS) has a critical role in cognition, sleep and other behaviors. Although not well studied in autism spectrum disorder (ASD), the HS is implicated in many neurological disorders, some of which share comorbidity with ASD, including Tourette syndrome (TS). Preliminary studies suggest that antagonism of histamine receptors 1-3 reduces symptoms and specific behaviors in ASD patients and relevant animal models. In addition, the HS mediates neuroinflammation, which may be heightened in ASD. Together, this suggests that the HS may also be altered in ASD. Using RNA sequencing (RNA-seq), we investigated genome-wide expression, as well as a focused gene set analysis of key HS genes (HDC, HNMT, HRH1, HRH2, HRH3 and HRH4) in postmortem dorsolateral prefrontal cortex (DLPFC) initially in 13 subjects with ASD and 39 matched controls. At the genome level, eight transcripts were differentially expressed (false discovery rate <0.05), six of which were small nucleolar RNAs (snoRNAs). There was no significant diagnosis effect on any of the individual HS genes but expression of the gene set of HNMT, HRH1, HRH2 and HRH3 was significantly altered. Curated HS gene sets were also significantly differentially expressed. Differential expression analysis of these gene sets in an independent RNA-seq ASD data set from DLPFC of 47 additional subjects confirmed these findings. Understanding the physiological relevance of an altered HS may suggest new therapeutic options for the treatment of ASD.

Transcriptomics as a tool to discover new antibacterial targets

The emergence of antibiotic-resistant pathogens, multiple drug-resistance, and extremely drug-resistant strains demonstrates the need for improved strategies to discover new drug-based compounds. The development of transcriptomics, proteomics, and metabolomics has provided new tools for global studies of living organisms. However, the compendium of expression profiles produced by these methods has introduced new scientific challenges into antimicrobial research. In this review, we discuss the practical value of transcriptomic techniques as well as their difficulties and pitfalls. We advocate the construction of new databases of transcriptomic data, using standardized formats in addition to standardized models of bacterial and yeast similar to those used in systems biology. The inclusion of proteomic and metabolomic data is also essential, as the resulting networks can provide a landscape to rationally predict and exploit new drug targets and to understand drug synergies.

Symmetric Directional False Discovery Rate Control

This research is motivated from the analysis of a real gene expression data that aims to identify a subset of "interesting" or "significant" genes for further studies. When we blindly applied the standard false discovery rate (FDR) methods, our biology collaborators were suspicious or confused, as the selected list of significant genes was highly unbalanced: there were ten times more under-expressed genes than the over-expressed genes. Their concerns led us to realize that the observed two-sample t-statistics were highly skewed and asymmetric, and thus the standard FDR methods might be inappropriate. To tackle this case, we propose a symmetric directional FDR control method that categorizes the genes into "over-expressed" and "under-expressed" genes, pairs "over-expressed" and "under-expressed" genes, defines the p-values for gene pairs via column permutations, and then applies the standard FDR method to select "significant" gene pairs instead of "significant" individual genes. We compare our proposed symmetric directional FDR method with the standard FDR method by applying them to simulated data and several well-known real data sets.

Resolving host-pathogen interactions by dual RNA-seq

The transcriptome is a powerful proxy for the physiological state of a cell, healthy or diseased. As a result, transcriptome analysis has become a key tool in understanding the molecular changes that accompany bacterial infections of eukaryotic cells. Until recently, such transcriptomic studies have been technically limited to analyzing mRNA expression changes in either the bacterial pathogen or the infected eukaryotic host cell. However, the increasing sensitivity of high-throughput RNA sequencing now enables "dual RNA-seq" studies, simultaneously capturing all classes of coding and noncoding transcripts in both the pathogen and the host. In the five years since the concept of dual RNA-seq was introduced, the technique has been applied to a range of infection models. This has not only led to a better understanding of the physiological changes in pathogen and host during the course of an infection but has also revealed hidden molecular phenotypes of virulence-associated small noncoding RNAs that were not visible in standard infection assays. Here, we use the knowledge gained from these recent studies to suggest experimental and computational guidelines for the design of future dual RNA-seq studies. We conclude this review by discussing prospective applications of the technique.

GSAR: Bioconductor package for Gene Set analysis in R

BACKGROUND

Gene set analysis (in a form of functionally related genes or pathways) has become the method of choice for analyzing omics data in general and gene expression data in particular. There are many statistical methods that either summarize gene-level statistics for a gene set or apply a multivariate statistic that accounts for intergene correlations. Most available methods detect complex departures from the null hypothesis but lack the ability to identify the specific alternative hypothesis that rejects the null.

RESULTS

GSAR (Gene Set Analysis in R) is an open-source R/Bioconductor software package for gene set analysis (GSA). It implements self-contained multivariate non-parametric statistical methods testing a complex null hypothesis against specific alternatives, such as differences in mean (shift), variance (scale), or net correlation structure. The package also provides a graphical visualization tool, based on the union of two minimum spanning trees, for correlation networks to examine the change in the correlation structures of a gene set between two conditions and highlight influential genes (hubs).

CONCLUSIONS

Package GSAR provides a set of multivariate non-parametric statistical methods that test a complex null hypothesis against specific alternatives. The methods in package GSAR are applicable to any type of omics data that can be represented in a matrix format. The package, with detailed instructions and examples, is freely available under the GPL (> = 2) license from the Bioconductor web site.

A Robust Analytical Pipeline for Genome-Wide Identification of the Genes Regulated by a Transcription Factor: Combinatorial Analysis Performed Using gSELEX-Seq and RNA-Seq

For identifying the genes that are regulated by a transcription factor (TF), we have established an analytical pipeline that combines genomic systematic evolution of ligands by exponential enrichment (gSELEX)-Seq and RNA-Seq. Here, SELEX was used to select DNA fragments from an Aspergillus nidulans genomic library that bound specifically to AmyR, a TF from A. nidulans. High-throughput sequencing data were obtained for the DNAs enriched through the selection, following which various in silico analyses were performed. Mapping reads to the genome revealed the binding motifs including the canonical AmyR-binding motif, CGGN8CGG, as well as the candidate promoters controlled by AmyR. In parallel, differentially expressed genes related to AmyR were identified by using RNA-Seq analysis with samples from A. nidulans WT and amyR deletant. By obtaining the intersecting set of genes detected using both gSELEX-Seq and RNA-Seq, the genes directly regulated by AmyR in A. nidulans can be identified with high reliability. This analytical pipeline is a robust platform for comprehensive genome-wide identification of the genes that are regulated by a target TF.