Next Article in Journal
Effect of Foliar Fertigation of Chitosan Nanoparticles on Cadmium Accumulation and Toxicity in Solanum lycopersicum
Next Article in Special Issue
Survival-Based Biomarker Module Identification Associated with Oral Squamous Cell Carcinoma (OSCC)
Previous Article in Journal
Downregulation of Ribosomal Contents and Kinase Activities Is Associated with the Inhibitive Effect on the Growth of Group B Streptococcus Induced by Placental Extracellular Vesicles
Previous Article in Special Issue
Genome-Wide Atlas of Promoter Expression Reveals Contribution of Transcribed Regulatory Elements to Genetic Control of Disuse-Mediated Atrophy of Skeletal Muscle
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 10
Issue 7
10.3390/biology10070665
biology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

FLAME: A Web Tool for Functional and Literature Enrichment Analysis of Multiple Gene Lists

by
Foteini Thanati
1,†,
Evangelos Karatzas
1,†,
Fotis A. Baltoumas
1,
Dimitrios J. Stravopodis
2,
Aristides G. Eliopoulos
3,4,5 and
Georgios A. Pavlopoulos
1,4,*
1
Institute for Fundamental Biomedical Research, BSRC “Alexander Fleming”, 16672 Vari, Greece
2
Department of Biology, School of Sciences, National and Kapodistrian University of Athens, 15701 Athens, Greece
3
Department of Biology, School of Medicine, National and Kapodistrian University of Athens, 11527 Athens, Greece
4
Center for New Biotechnologies and Precision Medicine, School of Medicine, National and Kapodistrian University of Athens, 11527 Athens, Greece
5
Center of Basic Research, Biomedical Research Foundation of the Academy of Athens, 11527 Athens, Greece
*
Author to whom correspondence should be addressed.
†
Equally contributing authors.
Biology 2021, 10(7), 665; https://doi.org/10.3390/biology10070665
Submission received: 15 June 2021 / Revised: 8 July 2021 / Accepted: 9 July 2021 / Published: 14 July 2021
(This article belongs to the Special Issue Differential Gene Expression and Coexpression)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Manipulation of multiple gene lists when going through a functional enrichment analysis is often a hefty task. In this article, we present FLAME, a web application that facilitates such multiple-list manipulations, enabling the construction of intersections and unions among multiple gene lists for further enrichment analysis. Reported results are presented via multiple interactive viewers such as heatmaps, barcharts, Manhattan plots, and networks.

Abstract

Functional enrichment is a widely used method for interpreting experimental results by identifying classes of proteins/genes associated with certain biological functions, pathways, diseases, or phenotypes. Despite the variety of existing tools, most of them can process a single list per time, thus making a more combinatorial analysis more complicated and prone to errors. In this article, we present FLAME, a web tool for combining multiple lists prior to enrichment analysis. Users can upload several lists and use interactive UpSet plots, as an alternative to Venn diagrams, to handle unions or intersections among the given input files. Functional and literature enrichment, along with gene conversions, are offered by g:Profiler and aGOtool applications for 197 organisms. FLAME can analyze genes/proteins for related articles, Gene Ontologies, pathways, annotations, regulatory motifs, domains, diseases, and phenotypes, and can also generate protein–protein interactions derived from STRING. We have validated FLAME by interrogating gene expression data associated with the sensitivity of the distal part of the large intestine to experimental colitis-propelled colon cancer. FLAME comes with an interactive user-friendly interface for easy list manipulation and exploration, while results can be visualized as interactive and parameterizable heatmaps, barcharts, Manhattan plots, networks, and tables.

1. Introduction

Functional enrichment analysis is a method to identify classes of bioentities in which genes or proteins have been found to be over-represented. This type of analysis can aid researchers to reveal biological insights from various omics experiments and interpret gene lists of interest in a biologically meaningful way. For this purpose, several applications have been proposed [1,2]. Established tools include: g:Profiler [3], Panther [4], DAVID [5], WebGestalt [6], EnrichR [7], AmiGO [8], GeneSCF [9], AllEnricher [10], aGOtool [11], ClueGO [12], Metascape [13], NeVOmics [14], GSEA [15], GOrilla [16], Fuento [17], and NASQAR [18]. Most of them are offered as web applications, and among other functionalities, they associate overrepresented genes with GO terms [19], or pathways [20,21,22,23]. Nevertheless, (a) they often differ in the organisms, the identifiers, and the backend database they support; (b) they come with their own way of reporting results (mostly lists or static representations); and (c) they are often unable to handle and compare multiple lists for a more combinatorial analysis.
To address such critical issues, in this article, we present FLAME, a web application that allows for the combination of multiple input gene lists, as well as their parallel exploration and analysis with the use of interactive UpSet plots. FLAME integrates g:Profiler and aGOtool for an accurate and always up-to-date functional and literature enrichment analysis, and also utilizes STRING’s API [24] to generate interactive protein–protein interaction (PPI) networks. A major advantage is that FLAME follows a visual analytics approach to allow users to adjust and parameterize the reported results with the engagement of heatmaps, barcharts, Manhattan plots, networks, and tables, thus making information easier to be absorbed, comprehended, and interpreted, and knowledge easier to be extracted and exploited.

2. Materials and Methods

2.1. Input

FLAME allows for the uploading of multiple gene lists as separate files, or as pasted text. In the online version, FLAME can accommodate up to 10 active gene lists with a file size smaller than 1 MB each, a limitation that can be bypassed by downloading the application from GitHub and running it locally after editing the value of the corresponding variable (FILE_LIMIT configuration variable in global.R and shiny.maxRequestSize option in ui.R). The input text can be imported in comma-, tab-, or line-separated format. In this version, FLAME supports 197 organisms and several identifiers, such as gene IDs (proteins, transcripts, microarray IDs, etc.), SNP IDs, chromosomal intervals, and term IDs, in accordance with the gconvert function of the gprofiler2 library. Different lists are not allowed to have the same name, with options for renaming and deleting them being suitably provided. Once uploaded, list contents can be seen in interactive searchable tables.

2.2. List Manipulation with the Use of UpSet Plots

After uploading the gene lists of interest and prior to analysis, UpSet plots can be generated to show possible intersections, unions, and distinct elements between the selected lists. The UpSet plot is a sophisticated alternative of a Venn diagram and is preferred for many sets (>5) where a Venn diagram becomes incomprehensive. In Figure 1, we demonstrate a simple example with three lists, each containing 100 random genes, and visualize the various UpSet plot options in relation to a Venn diagram (Figure 1A–D). Similarly, in Figure 1E, we describe the distinct intersections among seven gene lists, a task that cannot be drawn as a Venn diagram.
An UpSet plot consists of two axes and a connected-dot matrix. The vertical rectangles represent the number of elements participating in each list combination. The connected-dots matrix indicates which combination of lists corresponds to which vertical rectangle. Finally, the horizontal bars (Set Size) denote the participation of hovered objects (from the vertical rectangles) in the respective lists. In its current version, FLAME supports four UpSet plot modes: (I) intersections, (II) distinct intersections, (III) distinct elements per file, and (IV) unions. The intersections mode creates all file combinations as long as they share at least one element, allowing an element to participate in more than one combination. The distinct intersections mode creates file combinations only for distinct elements that do not participate in other lists. The distinct elements per file mode shows the distinct elements of each input list. The unions mode constructs all available file combinations, showing their total unique elements. Upon mouse-hover over each UpSet plot element combination (vertical rectangles), FLAME presents the respective genes in a table, whereas on a mouse click, the user can append the selected list of files, with a new file containing the selected genes and processing it separately.

2.3. Functional Enrichment

Once an input file or an UpSet plot column (intersection/union of sets) has been selected, FLAME takes advantage of the g:Profiler library [3] and aGOtool API [11] to offer functional enrichment analysis for a list of 197 organisms.
In detail, g:Profiler is used for the identification of enriched functional terms from Gene Ontology [19]; pathways from KEGG [20], Reactome [22], and WikiPathways [23]; protein complexes from CORUM [25]; expression data from the Human Protein Atlas [26]; regulatory motifs from TRANSFAC [27] and miRTarBase [28]; and phenotypes from the Human Phenotype Ontology [29]. Similarly, aGOtool is used for the identification of enriched terms from the UniProt keyword classification system [30], protein families and domains from Pfam [31] and InterPro [32], and human diseases from the DISEASES database [33].
g:Profiler and aGOtool test for statistically significant enrichment to compare the user’s input gene list to a background set from organism-specific genes annotated in the Ensembl database [34] and UniProt Reference Proteomes [30], respectively. In the case of g:Profiler, the resulting p-values are corrected for multiple testing using g:SCS, Bonferroni correction, or Benjamini–Hochberg false discovery rate (FDR), whereas in the case of aGOtool, p-values are corrected using Bonferroni correction or FDR. Notably, the reported lists with the enrichment results can be shrunk or expanded using the aforementioned parameters as thresholds.
Enrichment analysis is essentially performed using ENSEMBL identifiers, while, on the basis of the user’s choice, results can be reported in the input format or as Entrez, UniProt [30], EMBL [35], ENSEMBL [34], ChEMBL [36], WikiGene [37], and RefSeq [38] identifiers. Conversion to ENSEMBL identifiers is done internally at the backend, regardless of the identifier type imported by the user. Input elements that are not identified in the enriched terms are presented to the user, along with elements that could not be translated into the requested output format.

2.4. Literature Enrichment

In addition to functional enrichment, FLAME enables literature enrichment analysis for a selected gene list, via the aGOtool API. To this end, FLAME allows users to retrieve scientific articles that are tightly connected to the genes/proteins provided in the uploaded input files. The literature enrichment analysis concept is very similar to the one of the enrichment analyses and aims to aid users to identify scientific publications of relevance to a given gene/protein list.
The publication enrichment analysis is based on the aGOtool, which uses a text corpus of all PubMed abstracts and full-text open access articles from PubMed Central. These documents are processed by OnTheFly’s [39] or EXTRACT’s [40] underlying Named Entity Recognition (NER) tagger [41] on a weekly basis for the identification of biomedical entities (genes/proteins, chemical compounds, organisms, tissues, environments, diseases, phenotypes, and Gene Ontology terms). As a result, all documents are being automatically annotated with the genes mentioned within them, thus, as in functional enrichment, turning every document into a ‘gene set’.

2.5. Protein–Protein Interaction Analysis

FLAME offers the capability to construct and visualize interactive PPIs for a set of 197 organisms using the STRING API [24]. Users may submit their gene list and visualize the results as networks, with the interacting entities being presented as nodes and their interactions as edges. In the online version, FLAME allows for a maximum of 500 proteins per request, a limitation that can be bypassed by downloading the application from the GitHub repository and running it locally after editing the value of the corresponding variable (STRING_LIMIT configuration variable in global.R). We note that for each input gene, we only keep one converted Ensembl protein identifier.
STRING supports both physical and functional interactions. Physical interactions refer to proteins that are part of the same biomolecular complex, whereas functional interactions refer to proteins that are involved in the same pathway or biological process. Through FLAME, users can select whether to visualize the full set of interactions (both physical and functional), or just the physical sub-network. Users can also adjust the interaction score and apply a cutoff on the edges. In addition, users can choose between the evidence or the confidence mode. In the first case, a multi-edged graph is drawn, in which each edge shows the evidence channel (e.g., fusion event, co-expression, and text mining), and the information comes from [42], whereas in the second case, the thickness of the edge reflects the interaction score. While the resulting network is presented in a separate Network Viewer panel, one can export a network as an image or as a tab-delimited file to be visualized by external viewers [43,44,45] or get redirected in STRING’s original source for further analysis.

2.6. Visual Analytics and Interactive Visualization

FLAME offers various interactive plots and visualization options for reporting results. Functional enrichment results produced by g:Profiler and aGOtool, as well as literature enrichment results produced by aGOtool, can be shown as tables, scatter plots, barcharts, heatmaps, and networks.
Resulting lists of functional terms are initially reported as interactive searchable tables displaying details about each functional term. One can expand each row of the table to see which of the identified genes/proteins were found to be associated with the functional term. For example, in the case of a KEGG pathway, one can see how many proteins or genes were found to be related to it and get redirected to the KEGG repository to see the actual schema of the pathway in a static form with all of the detected genes/proteins highlighted.
In the case of g:Profiler only, an ‘adjusted to the selected data sources’ interactive Manhattan plot is offered for a clearer overview. In this plot, functional terms are organized along the x-axis and colored by their data source, whereas the y-axis shows the significance (p-value) of each term. Hovering over a data point generates a popup window with key information about the functional term. When a set of points are selected using a lasso or a rectangle, the Manhattan plot will be redrawn, showing information about the selected items only. Upon selection, the corresponding tables will be automatically updated.
In all of the enrichment-type analyses offered by FLAME, for every source, the most significant functional terms will be shown in a bar chart or a scatter plot, in which the user can further customize to adjust the desired number of terms analyzed. Bar charts are sorted according to the enrichment score or the p-value. Finally, on mouse hovering, a tooltip with key information about the functional term will be shown.
Heatmap and network visualization options complement each other and come with two different execution modes. In the first case, genes are plotted against functional terms, while in the second case, functional terms are plotted against themselves (square matrix). In the case of heatmaps, a color-scheme is used to capture the enrichment or the statistical score for a particular functional term. In the second mode (functional terms vs. functional terms), a cell value depicts a similarity capturing the number of common genes. This is calculated as the summation of the unique common genes between a pair of functional terms divided by the number of total unique genes found to be associated with the functional terms. All of the heatmaps are fully interactive, and one can zoom in and isolate an area of interest, swap the x- and y-axes, and adjust the number of elements that will be shown. Notably, all heatmaps are clustered after applying a hierarchical clustering method, and export options are also supported.
In addition to heatmaps, a three-mode network visualization is offered. In the first case, nodes represent genes and/or functional terms (colored differently), whereas edges represent the similarity scores between them, as explained previously (weighted networks). In the second case, nodes represent the functional terms, while edges reflect the number of common genes between them. Finally, in the third mode, nodes represent genes that are connected on the basis of the number of common functions or processes they are involved in.
While a heatmap is a preferable option for observing all possible pairwise similarities, in the second and third network mode users can apply an edge-cutoff, based on the similarity score or the common number of functions, respectively. This can help reduce the networks’ density and make visualization more appealing and more informative. Networks are fully interactive, and one can zoom in/out, adjust the layout accordingly, or export the network in various file formats. Examples of all of the aforementioned visualizations are shown in Figure 2.

2.7. Gene ID Conversion and Orthology Search

Frequently, different databases and tools accept varying gene/protein identifiers as their input. Thereby, the ID mapping creates a bottleneck for the non-expert end-users. To aid researchers in overcoming this problem, FLAME utilizes g:Profiler’s converters to allow (i) cross-database ID conversions and (ii) cross-species ID conversions (orthologs). In detail, FLAME allows ID conversions between well-known name spaces, such as Entrez Gene, Uniprot, ChEMBL, ENSEMBL, and RefSeq, as well as among the 197 different organisms supported by FLAME.

2.8. Implementation

FLAME is mainly written in R and JavaScript. The R/Shiny package was used for the GUI implementation and the interoperability between R and Javascript. UpSet plots have been implemented with the use of R/upsetjs library. Functional enrichment analyses are offered by the R/gprofiler2 library [46] and aGOtool API. The latter is also used for literature enrichment analysis. Networks are stored as igraph objects [47] and visualized with the visNetwork library. Scatter and Manhattan plots are generated with the help of the R/plotly library [48], bar plots with R/ggplot2 [49], and heatmaps via the R/heatmaply library [50]. Interactive tables are generated through the R/DT library. Finally, network analysis is performed with the employment of the STRING API.

2.9. Integration with Other Applications

FLAME can be called from other applications with a simple GET request. Gene names must be encoded in the URL (url_genes variable) and be comma-separated (,), whereas multiple gene lists must be separated by the semicolon symbol (;). Lists will appear as uploaded files. Simple URL examples encoding three lists are http://flame.pavlopouloslab.info/?url_genes=MCL1,TTR;APOE,ACE2;TLR4,HMOX1,TP73 (accessed on 12 July 2021) or http://bib.fleming.gr:3838/Flame/?url_genes=MCL1,TTR;APOE,ACE2;TLR4,HMOX1,TP73 (accessed on 12 July 2021).

3. Results and Discussion

We tested the capacity of FLAME for knowledge discovery by re-analyzing gene expression data associated with colitis-propelled carcinogenesis (CAC) in mice [51]. In this model, tumorigenesis develops following a single application of the carcinogen azoxymethane (AOM) combined with four cycles of dextran sodium sulfate (DSS) administration that causes chronic colitis. However, inflammation, tissue damage, dysplasia, and cancer are manifested in the distal but not the proximal part of the large intestine. To gain insight into the biological basis of this intriguing phenomenon, we previously interrogated the transcriptome of proximal and distal colon and reported intrinsic differences in gene expression, which are augmented during CAC. These analyses also provided evidence that lipid metabolic pathways operating in the proximal part of the large intestine mediate resistance to experimental colitis and CAC [51].
Herein, we focused on exploring gene ontologies and pathways that may mediate sensitivity to CAC. We reasoned that, transcripts that are upregulated at both the ‘early’ (i.e., 2 DSS cycles) and ‘late’ (i.e., 4 DSS cycles) stages of AOM/DSS-induced carcinogenesis in the distal part of the colon and are not detected in the proximal region must play a prominent role in CAC. We identified 165 transcripts belonging to this intersection by using the ‘UpSet Plot’ option of FLAME (Figure 3A and Table S1), and termed this gene set ‘the susceptibility-associated gene signature’ (SAS).
To interpret SAS transcripts as biological functions and pathways, we combined KEGG and REACTOME pathway enrichment analysis using FLAME. Reassuringly, we found ‘inflammatory bowel disease’ to be among the most significantly enriched pathways (Figure 3B and Table S2). Other pathways found to be significantly enriched in SAS were related to immune function, including T helper 1 (Th1)/Th2 differentiation, Th17 cell differentiation, natural killer (NK) cell-mediated cytotoxicity, and cytokine/chemokine signaling, which align with the pivotal role of unresolved inflammation in promoting intestinal cancer [52]. Indeed, Th17 lymphocytes give rise to pathogenic Th1 cells implicated in colitis [53], and when activated by TGFβ1 (a SAS gene set cytokine; Table S1), Th17 cells promote CAC [54]. Moreover, chronic experimental colitis is associated with dual Th1 and Th2 cytokine profile [55], and there is evidence to suggest that Th2-driven colonic inflammation enhances the formation of colorectal tumors [56]. We also note the significant enrichment of SAS for the pro-inflammatory TNF pathway, which represents a major therapeutic target for inflammatory bowel diseases.
Gene ontology (GO) enrichment was also performed on the FLAME platform. SAS was found to be enriched for biological processes that were mostly related to T cell-mediated immunity and cytokine synthesis (Figure S1A and Table S3). The GO enrichment for molecular functions (Figure S1B and Table S3) identified several terms related to cytokine and chemokine activity, including TNF, and to death receptor signaling, which has been implicated in CAC [57]. A gene–gene network was also constructed using FLAME to uncover clusters of SAS genes with common molecular functions (Figure S1C). Analysis of putative protein–protein interactions (physical and functional associations) through the STRING option of FLAME uncovered a cluster of interacting components of T cell-associated immunity and a cluster of interferon beta (IFNβ) response proteins implicated in epithelial regeneration upon DSS-induced tissue damage [58] (Figure 3C).
To gain putative mechanistic insights into the aforementioned biological functions and pathway data, we used the FLAME platform to perform transcription factor motif enrichment and analyze transcription factor networks versus SAS transcripts. The results showed a cluster of interferon regulatory factor (IRF) family members (Figure S1D and Table S4), which may be linked to both the enrichment of SAS in IFNβ components and T helper cell differentiation. Another transcription factor cluster was associated with SMADs, which have been implicated in CAC, in part downstream of TGFβ1 [59,60]. Overall, the aforementioned analyses underscore the practical utility of FLAME to rapidly process, analyze, and visualize gene expression data, and hence assist knowledge discovery.

4. Conclusions

FLAME is a web application aimed to effectively complement existing enrichment tools and make gene list analysis, information extraction, knowledge integration, and visualization much easier and more comprehensive, accurate, and reliable. It combines functionalities from g:Profiler, aGOtool, and STRING to cover a wide spectrum of analyses, such as functional, literature enrichment, and network analysis. FLAME is designed for non-experts and via the engagement of UpSet plots; it allows users to manipulate, compare, and handle multiple gene lists prior to any type of analysis. Moreover, it follows a visual analytics approach, thus offering users various options for visualizing their results and adjusting the corresponding views upon parameterization to generate publication-accepted figures. In a future version, we expect to expand FLAME’s capabilities by supporting more organisms and integrating a greater plethora of sources. We expect FLAME to reach out to many users varying from non-experts to highly specialized bioinformaticians.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/biology10070665/s1, Figure S1: (A) Gene Ontology (GO)-enriched biological processes (BP) of ‘susceptibility-associated gene signature’ (SAS). (B) GO-enriched molecular functions (MF) of SAS. (C) An MF-enriched gene–gene network, where an edge width is proportional to the number of common molecular functions between two genes. (D) An unweighted network of transcription factors (functions vs. genes) predicted to be involved in the regulation of SAS genes. Table S1: Gene names comprising the ‘susceptibility-associated gene signature’ (SAS). Table S2: List of KEGG and REACTOME pathways enriched in SAS. Table S3: List of the most enriched GO biological processes (BP) and molecular functions (MF) in SAS. Table S4: List of transcription factors predicted to regulate SAS genes.

Author Contributions

Conceptualization and supervision, G.A.P.; methodology, E.K. and F.A.B.; software, F.T., E.K. and F.A.B.; validation, D.J.S., A.G.E. and G.A.P.; case study, A.G.E. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financially supported by the Operational Program Competitiveness, Entrepreneurship and Innovation, NSRF 2014-2020, action code: MIS 5002562, co-financed by Greece and the European Union (European Regional Development Fund). The study was also supported by the Hellenic Foundation for Research and Innovation (H.F.R.I.) under the ‘First Call for H.F.R.I. Research Projects to support faculty members and researchers and the procurement of high-cost research equipment grant’, grant ID: 1855-BOLOGNA. We also acknowledge support of this work by the project ‘The Greek Research Infrastructure for Personalised Medicine (pMedGR)’ (MIS 5002802), which is implemented under the Action ‘Reinforcement of the Research and Innovation Infrastructure’, funded by the Operational Program ‘Competitiveness, Entrepreneurship and Innovation’ (NSRF 2014-2020) and co-financed by Greece and the European Union (European Regional Development Fund).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

FLAME is available as a web service through http://bib.fleming.gr:3838/Flame/ (accessed on 12 July 2021) or http://flame.pavlopouloslab.info (accessed on 12 July 2021). The source code is publicly available for download in https://github.com/PavlopoulosLab/Flame (accessed on 12 July 2021).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maleki, F.; Ovens, K.; Hogan, D.J.; Kusalik, A.J. Gene Set Analysis: Challenges, Opportunities, and Future Research. Front. Genet. 2020, 11, 654. [Google Scholar] [CrossRef]
  2. Shi-Jing, L.; Fathiah-Muzaffar-Shah, F.; Saberi-Mohamad, M.; Moorthy, K.; Deris, S.; Zakaria, Z.; Napis, S. A Review on Bioinformatics Enrichment Analysis Tools Towards Functional Analysis of High Throughput Gene Set Data. Curr. Proteom. 2015, 12, 14–27. [Google Scholar] [CrossRef]
  3. Raudvere, U.; Kolberg, L.; Kuzmin, I.; Arak, T.; Adler, P.; Peterson, H.; Vilo, J. G:Profiler: A Web Server for Functional Enrichment Analysis and Conversions of Gene Lists (2019 Update). Nucleic Acids Res. 2019, 47, W191–W198. [Google Scholar] [CrossRef] [Green Version]
  4. Mi, H.; Muruganujan, A.; Casagrande, J.T.; Thomas, P.D. Large-Scale Gene Function Analysis with the PANTHER Classification System. Nat. Protoc. 2013, 8, 1551–1566. [Google Scholar] [CrossRef] [PubMed]
  5. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and Integrative Analysis of Large Gene Lists Using DAVID Bioinformatics Resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef] [PubMed]
  6. Wang, J.; Vasaikar, S.; Shi, Z.; Greer, M.; Zhang, B. WebGestalt 2017: A More Comprehensive, Powerful, Flexible and Interactive Gene Set Enrichment Analysis Toolkit. Nucleic Acids Res. 2017, 45, W130–W137. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, E.Y.; Tan, C.M.; Kou, Y.; Duan, Q.; Wang, Z.; Meirelles, G.V.; Clark, N.R.; Ma’ayan, A. Enrichr: Interactive and Collaborative HTML5 Gene List Enrichment Analysis Tool. BMC Bioinform. 2013, 14, 128. [Google Scholar] [CrossRef] [Green Version]
  8. Carbon, S.; Ireland, A.; Mungall, C.J.; Shu, S.; Marshall, B.; Lewis, S. AmiGO Hub, Web Presence Working Group AmiGO: Online Access to Ontology and Annotation Data. Bioinformatics 2009, 25, 288–289. [Google Scholar] [CrossRef]
  9. Subhash, S.; Kanduri, C. GeneSCF: A Real-Time Based Functional Enrichment Tool with Support for Multiple Organisms. BMC Bioinform. 2016, 17, 365. [Google Scholar] [CrossRef] [Green Version]
  10. Zhang, D.; Hu, Q.; Liu, X.; Zou, K.; Sarkodie, E.K.; Liu, X.; Gao, F. AllEnricher: A Comprehensive Gene Set Function Enrichment Tool for Both Model and Non-Model Species. BMC Bioinform. 2020, 21, 106. [Google Scholar] [CrossRef] [Green Version]
  11. Schölz, C.; Lyon, D.; Refsgaard, J.C.; Jensen, L.J.; Choudhary, C.; Weinert, B.T. Avoiding Abundance Bias in the Functional Annotation of Post-Translationally Modified Proteins. Nat. Methods 2015, 12, 1003–1004. [Google Scholar] [CrossRef]
  12. Bindea, G.; Mlecnik, B.; Hackl, H.; Charoentong, P.; Tosolini, M.; Kirilovsky, A.; Fridman, W.-H.; Pagès, F.; Trajanoski, Z.; Galon, J. ClueGO: A Cytoscape Plug-in to Decipher Functionally Grouped Gene Ontology and Pathway Annotation Networks. Bioinformatics 2009, 25, 1091–1093. [Google Scholar] [CrossRef] [Green Version]
  13. Zhou, Y.; Zhou, B.; Pache, L.; Chang, M.; Khodabakhshi, A.H.; Tanaseichuk, O.; Benner, C.; Chanda, S.K. Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat. Commun. 2019, 10, 1523. [Google Scholar] [CrossRef]
  14. Zúñiga-León, E.; Carrasco-Navarro, U.; Fierro, F. NeVOmics: An Enrichment Tool for Gene Ontology and Functional Network Analysis and Visualization of Data from OMICs Technologies. Genes 2018, 9, 569. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene Set Enrichment Analysis: A Knowledge-Based Approach for Interpreting Genome-Wide Expression Profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Eden, E.; Navon, R.; Steinfeld, I.; Lipson, D.; Yakhini, Z. GOrilla: A Tool for Discovery and Visualization of Enriched GO Terms in Ranked Gene Lists. BMC Bioinform. 2009, 10, 48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Weichselbaum, D.; Zagrovic, B.; Polyansky, A.A. Fuento: Functional Enrichment for Bioinformatics. Bioinformatics 2017, 33, 2604–2606. [Google Scholar] [CrossRef]
  18. Yousif, A.; Drou, N.; Rowe, J.; Khalfan, M.; Gunsalus, K.C. NASQAR: A Web-Based Platform for High-Throughput Sequencing Data Analysis and Visualization. BMC Bioinform. 2020, 21, 267. [Google Scholar] [CrossRef] [PubMed]
  19. Gene Ontology Consortium. The Gene Ontology (GO) Database and Informatics Resource. Nucleic Acids Res. 2004, 32, D258–D261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  21. Okuda, S.; Yamada, T.; Hamajima, M.; Itoh, M.; Katayama, T.; Bork, P.; Goto, S.; Kanehisa, M. KEGG Atlas Mapping for Global Analysis of Metabolic Pathways. Nucleic Acids Res. 2008, 36, W423–W426. [Google Scholar] [CrossRef]
  22. Fabregat, A.; Jupe, S.; Matthews, L.; Sidiropoulos, K.; Gillespie, M.; Garapati, P.; Haw, R.; Jassal, B.; Korninger, F.; May, B.; et al. The Reactome Pathway Knowledgebase. Nucleic Acids Res. 2018, 46, D649–D655. [Google Scholar] [CrossRef] [PubMed]
  23. Slenter, D.N.; Kutmon, M.; Hanspers, K.; Riutta, A.; Windsor, J.; Nunes, N.; Mélius, J.; Cirillo, E.; Coort, S.L.; Digles, D.; et al. WikiPathways: A Multifaceted Pathway Database Bridging Metabolomics to Other Omics Research. Nucleic Acids Res. 2018, 46, D661–D667. [Google Scholar] [CrossRef]
  24. Szklarczyk, D.; Franceschini, A.; Wyder, S.; Forslund, K.; Heller, D.; Huerta-Cepas, J.; Simonovic, M.; Roth, A.; Santos, A.; Tsafou, K.P.; et al. STRING V10: Protein-Protein Interaction Networks, Integrated over the Tree of Life. Nucleic Acids Res. 2015, 43, D447–D452. [Google Scholar] [CrossRef] [PubMed]
  25. Giurgiu, M.; Reinhard, J.; Brauner, B.; Dunger-Kaltenbach, I.; Fobo, G.; Frishman, G.; Montrone, C.; Ruepp, A. CORUM: The Comprehensive Resource of Mammalian Protein Complexes-2019. Nucleic Acids Res. 2019, 47, D559–D563. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Uhlén, M.; Fagerberg, L.; Hallström, B.M.; Lindskog, C.; Oksvold, P.; Mardinoglu, A.; Sivertsson, Å.; Kampf, C.; Sjöstedt, E.; Asplund, A.; et al. Proteomics. Tissue-Based Map of the Human Proteome. Science 2015, 347, 1260419. [Google Scholar] [CrossRef]
  27. Matys, V.; Fricke, E.; Geffers, R.; Gössling, E.; Haubrock, M.; Hehl, R.; Hornischer, K.; Karas, D.; Kel, A.E.; Kel-Margoulis, O.V.; et al. TRANSFAC: Transcriptional Regulation, from Patterns to Profiles. Nucleic Acids Res. 2003, 31, 374–378. [Google Scholar] [CrossRef]
  28. Huang, H.-Y.; Lin, Y.-C.-D.; Li, J.; Huang, K.-Y.; Shrestha, S.; Hong, H.-C.; Tang, Y.; Chen, Y.-G.; Jin, C.-N.; Yu, Y.; et al. MiRTarBase 2020: Updates to the Experimentally Validated MicroRNA–Target Interaction Database. Nucleic Acids Res. 2019, gkz896. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Robinson, P.N.; Köhler, S.; Bauer, S.; Seelow, D.; Horn, D.; Mundlos, S. The Human Phenotype Ontology: A Tool for Annotating and Analyzing Human Hereditary Disease. Am. J. Hum. Genet. 2008, 83, 610–615. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. UniProt Consortium, T. UniProt: The Universal Protein Knowledgebase. Nucleic Acids Res. 2018, 46, 2699. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Finn, R.D.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Mistry, J.; Mitchell, A.L.; Potter, S.C.; Punta, M.; Qureshi, M.; Sangrador-Vegas, A.; et al. The Pfam Protein Families Database: Towards a More Sustainable Future. Nucleic Acids Res. 2016, 44, D279–D285. [Google Scholar] [CrossRef] [PubMed]
  32. Blum, M.; Chang, H.-Y.; Chuguransky, S.; Grego, T.; Kandasaamy, S.; Mitchell, A.; Nuka, G.; Paysan-Lafosse, T.; Qureshi, M.; Raj, S.; et al. The InterPro Protein Families and Domains Database: 20 Years On. Nucleic Acids Res. 2021, 49, D344–D354. [Google Scholar] [CrossRef] [PubMed]
  33. Pletscher-Frankild, S.; Pallejà, A.; Tsafou, K.; Binder, J.X.; Jensen, L.J. DISEASES: Text Mining and Data Integration of Disease-Gene Associations. Methods 2015, 74, 83–89. [Google Scholar] [CrossRef] [PubMed]
  34. Zerbino, D.R.; Achuthan, P.; Akanni, W.; Amode, M.R.; Barrell, D.; Bhai, J.; Billis, K.; Cummins, C.; Gall, A.; Girón, C.G.; et al. Ensembl 2018. Nucleic Acids Res. 2018, 46, D754–D761. [Google Scholar] [CrossRef] [PubMed]
  35. Baker, W. The EMBL Nucleotide Sequence Database. Nucleic Acids Res. 2000, 28, 19–23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Gaulton, A.; Bellis, L.J.; Bento, A.P.; Chambers, J.; Davies, M.; Hersey, A.; Light, Y.; McGlinchey, S.; Michalovich, D.; Al-Lazikani, B.; et al. ChEMBL: A Large-Scale Bioactivity Database for Drug Discovery. Nucleic Acids Res. 2012, 40, D1100–D1107. [Google Scholar] [CrossRef] [Green Version]
  37. Hoffmann, R. A Wiki for the Life Sciences Where Authorship Matters. Nat. Genet. 2008, 40, 1047–1051. [Google Scholar] [CrossRef]
  38. O’Leary, N.A.; Wright, M.W.; Brister, J.R.; Ciufo, S.; Haddad, D.; McVeigh, R.; Rajput, B.; Robbertse, B.; Smith-White, B.; Ako-Adjei, D.; et al. Reference Sequence (RefSeq) Database at NCBI: Current Status, Taxonomic Expansion, and Functional Annotation. Nucleic Acids Res. 2016, 44, D733–D745. [Google Scholar] [CrossRef] [Green Version]
  39. Baltoumas, F.A.; Zafeiropoulou, S.; Karatzas, E.; Paragkamian, S.; Thanati, F.; Iliopoulos, I.; Eliopoulos, A.G.; Schneider, R.; Jensen, L.J.; Pafilis, E.; et al. OnTheFly 2.0: A Text.-Mining Web Application for Automated Biomedical Entity Recognition, Document Annotation, Network and Functional Enrichment Analysis. bioRxiv 2021. [Google Scholar] [CrossRef]
  40. Pafilis, E.; Buttigieg, P.L.; Ferrell, B.; Pereira, E.; Schnetzer, J.; Arvanitidis, C.; Jensen, L.J. EXTRACT: Interactive Extraction of Environment Metadata and Term Suggestion for Metagenomic Sample Annotation. Database 2016, 2016. [Google Scholar] [CrossRef] [Green Version]
  41. Pafilis, E.; Jensen, L.J. Real-Time Tagging of Biomedical Entities. bioRxiv 2016, 078469. [Google Scholar] [CrossRef]
  42. Koutrouli, M.; Hatzis, P.; Pavlopoulos, G.A. Exploring Networks in the STRING and Reactome Database. In Systems Medicine; Wolkenhauer, O., Ed.; Academic Press: Oxford, UK, 2021; pp. 507–520. ISBN 978-0-12-816078-7. [Google Scholar]
  43. Koutrouli, M.; Karatzas, E.; Paez-Espino, D.; Pavlopoulos, G.A. A Guide to Conquer the Biological Network Era Using Graph Theory. Front. Bioeng. Biotechnol. 2020, 8, 34. [Google Scholar] [CrossRef]
  44. Pavlopoulos, G.A.; Wegener, A.-L.; Schneider, R. A Survey of Visualization Tools for Biological Network Analysis. BioData Min. 2008, 1, 12. [Google Scholar] [CrossRef] [Green Version]
  45. O’Donoghue, S.I.; Gavin, A.-C.; Gehlenborg, N.; Goodsell, D.S.; Hériché, J.-K.; Nielsen, C.B.; North, C.; Olson, A.J.; Procter, J.B.; Shattuck, D.W.; et al. Visualizing Biological Data-Now and in the Future. Nat. Methods 2010, 7, S2–S4. [Google Scholar] [CrossRef]
  46. Kolberg, L.; Raudvere, U.; Kuzmin, I.; Vilo, J.; Peterson, H. Gprofiler2—An R Package for Gene List Functional Enrichment Analysis and Namespace Conversion Toolset g:Profiler. F1000Res 2020, 9. [Google Scholar] [CrossRef]
  47. Csardi, T.; Nepusz, T. The Igraph Software Package for Complex Network Research. InterJournal Complex. Syst. 2006, 1695, 1–9. [Google Scholar]
  48. Sievert, C. Interactive Web-Based Data Visualization with R, Plotly, and Shiny; CRC Press/Taylor and Francis Group: Boca Raton, FL, USA, 2020; ISBN 978-1-138-33149-5. [Google Scholar]
  49. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis, Use R! 2nd ed.; Springer: Cham, Switzerland, 2016; ISBN 978-3-319-24277-4. [Google Scholar]
  50. Galili, T.; O’Callaghan, A.; Sidi, J.; Sievert, C. Heatmaply: An R Package for Creating Interactive Cluster Heatmaps for Online Publishing. Bioinformatics 2018, 34, 1600–1602. [Google Scholar] [CrossRef]
  51. Gkouskou, K.K.; Ioannou, M.; Pavlopoulos, G.A.; Georgila, K.; Siganou, A.; Nikolaidis, G.; Kanellis, D.C.; Moore, S.; Papadakis, K.A.; Kardassis, D.; et al. Apolipoprotein A-I Inhibits Experimental Colitis and Colitis-Propelled Carcinogenesis. Oncogene 2016, 35, 2496–2505. [Google Scholar] [CrossRef]
  52. Grivennikov, S.I. Inflammation and Colorectal Cancer: Colitis-Associated Neoplasia. Semin. Immunopathol. 2013, 35, 229–244. [Google Scholar] [CrossRef]
  53. Harbour, S.N.; Maynard, C.L.; Zindl, C.L.; Schoeb, T.R.; Weaver, C.T. Th17 Cells Give Rise to Th1 Cells That Are Required for the Pathogenesis of Colitis. Proc. Natl. Acad. Sci. USA 2015, 112, 7061–7066. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Perez, L.G.; Kempski, J.; McGee, H.M.; Pelzcar, P.; Agalioti, T.; Giannou, A.; Konczalla, L.; Brockmann, L.; Wahib, R.; Xu, H.; et al. TGF-β Signaling in Th17 Cells Promotes IL-22 Production and Colitis-Associated Colon Cancer. Nat. Commun. 2020, 11, 2608. [Google Scholar] [CrossRef] [PubMed]
  55. Dieleman, L.A.; Palmen, M.J.; Akol, H.; Bloemena, E.; Peña, A.S.; Meuwissen, S.G.; Van Rees, E.P. Chronic Experimental Colitis Induced by Dextran Sulphate Sodium (DSS) Is Characterized by Th1 and Th2 Cytokines. Clin. Exp. Immunol. 1998, 114, 385–391. [Google Scholar] [CrossRef]
  56. Osawa, E.; Nakajima, A.; Fujisawa, T.; Kawamura, Y.I.; Toyama-Sorimachi, N.; Nakagama, H.; Dohi, T. Predominant T Helper Type 2-Inflammatory Responses Promote Murine Colon Cancers. Int. J. Cancer 2006, 118, 2232–2236. [Google Scholar] [CrossRef] [PubMed]
  57. Vlantis, K.; Wullaert, A.; Polykratis, A.; Kondylis, V.; Dannappel, M.; Schwarzer, R.; Welz, P.; Corona, T.; Walczak, H.; Weih, F.; et al. NEMO Prevents RIP Kinase 1-Mediated Epithelial Cell Death and Chronic Intestinal Inflammation by NF-ΚB-Dependent and -Independent Functions. Immunity 2016, 44, 553–567. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. McElrath, C.; Espinosa, V.; Lin, J.-D.; Peng, J.; Sridhar, R.; Dutta, O.; Tseng, H.-C.; Smirnov, S.V.; Risman, H.; Sandoval, M.J.; et al. Critical Role of Interferons in Gastrointestinal Injury Repair. Nat. Commun. 2021, 12, 2624. [Google Scholar] [CrossRef] [PubMed]
  59. Means, A.L.; Freeman, T.J.; Zhu, J.; Woodbury, L.G.; Marincola-Smith, P.; Wu, C.; Meyer, A.R.; Weaver, C.J.; Padmanabhan, C.; An, H.; et al. Epithelial Smad4 Deletion Up-Regulates Inflammation and Promotes Inflammation-Associated Cancer. Cell Mol. Gastroenterol. Hepatol. 2018, 6, 257–276. [Google Scholar] [CrossRef] [Green Version]
  60. Troncone, E.; Marafini, I.; Stolfi, C.; Monteleone, G. Transforming Growth Factor-Β1/Smad7 in Intestinal Immunity, Inflammation, and Cancer. Front. Immunol. 2018, 9, 1407. [Google Scholar] [CrossRef] [PubMed]
Biology 10 00665 g001 550
Figure 1. UpSet plot vs. Venn diagrams. (A) Intersection of the three gene lists (100 genes each) shown in a Venn diagram. (B) The UpSet plot Intersection option visualizes the total number of common elements among the selected sets, even though they may also participate in other sets. For example, lists 2 and 3 contain 29 common genes (green rectangle), with 25 being shared only between them, and 4 also shared with List 1, as it can be seen in (A). (C) The distinct intersections option visualizes the common number of genes, among the chosen sets, which do not exist in any other set. This option is the closest to a Venn diagram. For example, lists 2 and 3 share 25 distinct genes (red rectangle). (D) The union option appends the unique elements among the chosen sets and creates all possible combinations. For example, the combination of lists 2 and 3 results in 171 total unique genes (blue rectangle). (E) A ‘distinct intersections’ UpSet plot example with 7 lists, which cannot be visualized as a Venn diagram.
Figure 1. UpSet plot vs. Venn diagrams. (A) Intersection of the three gene lists (100 genes each) shown in a Venn diagram. (B) The UpSet plot Intersection option visualizes the total number of common elements among the selected sets, even though they may also participate in other sets. For example, lists 2 and 3 contain 29 common genes (green rectangle), with 25 being shared only between them, and 4 also shared with List 1, as it can be seen in (A). (C) The distinct intersections option visualizes the common number of genes, among the chosen sets, which do not exist in any other set. This option is the closest to a Venn diagram. For example, lists 2 and 3 share 25 distinct genes (red rectangle). (D) The union option appends the unique elements among the chosen sets and creates all possible combinations. For example, the combination of lists 2 and 3 results in 171 total unique genes (blue rectangle). (E) A ‘distinct intersections’ UpSet plot example with 7 lists, which cannot be visualized as a Venn diagram.
Biology 10 00665 g001
Biology 10 00665 g002 550
Figure 2. FLAME’s interactive visualization options. (A) g:Profiler’s result overview in a Manhattan plot. (B) A barchart, incorporating the top 17 results after combining different sources such as GO Biological Process, GO Cellular Component, GO Molecular Function, phenotypes, protein complexes, pathways, and miRNAs. (C) A hierarchically clustered heatmap with pairwise similarities between reactome terms. The color intensity reflects the common genes between each pair. (D) A heatmap showing genes (x-axis) associated with Disease Ontology terms (y-axis). (E) A scatter plot of enriched GO Molecular Function terms. The x-axis represents the enrichment score, while the y-axis the statistical significance. (F) An unweighted network consisting of GO Molecular Function functional terms and their related genes. (G) A weighted network of pairwise UniProt keywords returned by aGOtool. Edge weights are proportional to the number of common genes between each pair. (H) A gene–gene network based on common GOs, where an edge width is proportional to the number of common biological processes (BP) between two genes. (I) An interactive table with enrichment analysis results reported in different tabs (source of the active tab shown: Gene Ontology Cellular Component).
Figure 2. FLAME’s interactive visualization options. (A) g:Profiler’s result overview in a Manhattan plot. (B) A barchart, incorporating the top 17 results after combining different sources such as GO Biological Process, GO Cellular Component, GO Molecular Function, phenotypes, protein complexes, pathways, and miRNAs. (C) A hierarchically clustered heatmap with pairwise similarities between reactome terms. The color intensity reflects the common genes between each pair. (D) A heatmap showing genes (x-axis) associated with Disease Ontology terms (y-axis). (E) A scatter plot of enriched GO Molecular Function terms. The x-axis represents the enrichment score, while the y-axis the statistical significance. (F) An unweighted network consisting of GO Molecular Function functional terms and their related genes. (G) A weighted network of pairwise UniProt keywords returned by aGOtool. Edge weights are proportional to the number of common genes between each pair. (H) A gene–gene network based on common GOs, where an edge width is proportional to the number of common biological processes (BP) between two genes. (I) An interactive table with enrichment analysis results reported in different tabs (source of the active tab shown: Gene Ontology Cellular Component).
Biology 10 00665 g002
Biology 10 00665 g003 550
Figure 3. (A) UpSet plot of distinct intersections of transcripts up-regulated at the ‘early’ (i.e., 2 DSS cycles) and ‘late’ (i.e., 4 DSS cycles) stages of AOM/DSS-induced carcinogenesis in the disease-resistant proximal and disease-sensitive distal part of the colon. The ‘susceptibility-associated gene signature’ (SAS), comprising 165 transcripts, was used for gene ontology and pathway enrichment. (B) KEGG and Reactome pathway enrichment of SAS transcripts. (C) Analysis of putative protein–protein interactions through the STRING option of FLAME. A cluster of interacting components of interferon beta (IFNβ) response and of T cell-associated immunity, including Th1/Th2 differentiation, Th17 differentiation, and cytokine–cytokine receptor interaction, each represented by a different color, is shown.
Figure 3. (A) UpSet plot of distinct intersections of transcripts up-regulated at the ‘early’ (i.e., 2 DSS cycles) and ‘late’ (i.e., 4 DSS cycles) stages of AOM/DSS-induced carcinogenesis in the disease-resistant proximal and disease-sensitive distal part of the colon. The ‘susceptibility-associated gene signature’ (SAS), comprising 165 transcripts, was used for gene ontology and pathway enrichment. (B) KEGG and Reactome pathway enrichment of SAS transcripts. (C) Analysis of putative protein–protein interactions through the STRING option of FLAME. A cluster of interacting components of interferon beta (IFNβ) response and of T cell-associated immunity, including Th1/Th2 differentiation, Th17 differentiation, and cytokine–cytokine receptor interaction, each represented by a different color, is shown.
Biology 10 00665 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Thanati, F.; Karatzas, E.; Baltoumas, F.A.; Stravopodis, D.J.; Eliopoulos, A.G.; Pavlopoulos, G.A. FLAME: A Web Tool for Functional and Literature Enrichment Analysis of Multiple Gene Lists. Biology 2021, 10, 665. https://doi.org/10.3390/biology10070665

AMA Style

Thanati F, Karatzas E, Baltoumas FA, Stravopodis DJ, Eliopoulos AG, Pavlopoulos GA. FLAME: A Web Tool for Functional and Literature Enrichment Analysis of Multiple Gene Lists. Biology. 2021; 10(7):665. https://doi.org/10.3390/biology10070665

Chicago/Turabian Style

Thanati, Foteini, Evangelos Karatzas, Fotis A. Baltoumas, Dimitrios J. Stravopodis, Aristides G. Eliopoulos, and Georgios A. Pavlopoulos. 2021. "FLAME: A Web Tool for Functional and Literature Enrichment Analysis of Multiple Gene Lists" Biology 10, no. 7: 665. https://doi.org/10.3390/biology10070665

APA Style

Thanati, F., Karatzas, E., Baltoumas, F. A., Stravopodis, D. J., Eliopoulos, A. G., & Pavlopoulos, G. A. (2021). FLAME: A Web Tool for Functional and Literature Enrichment Analysis of Multiple Gene Lists. Biology, 10(7), 665. https://doi.org/10.3390/biology10070665

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop