# Gene Ontology Enrichment Analysis

# Installation

    BiocManager::install("TADCompare")

# Introduction

Using the output of `TADCompare` and `TimeCompare`, we can do a range of analyses. One common one is gene ontology enrichment analysis to determine the pathways in which genes near TAD boundaries occur in. To do this, we use rGREAT an R package for performing gene ontology enrichment analysis.

# Performing gene ontology analysis using TADCompare

In the first example, we show how to perform gene ontology enrichment using differential boundaries. Here, we perform the analysis on shifted boundaries detected in matrix 1.

    library(rGREAT)
    # Reading in data
    data("rao_chr22_prim")
    data("rao_chr22_rep")

    # Performing differential analysis
    results <- TADCompare(rao_chr22_prim, rao_chr22_rep, resolution = 50000)

    # Saving the results into its own data frame
    TAD_Frame <- results$TAD_Frame

    # Filter data to only include complex boundaries enriched in the second
    # contact matrix
    TAD_Frame <- TAD_Frame %>% dplyr::filter((Type == "Shifted") &
                                             (Enriched_In == "Matrix 2"))

    # Assign a chromosome and convert to a bed format
    TAD_Frame <- TAD_Frame %>% dplyr::select(Boundary) %>% mutate(chr = "chr22",
        start = Boundary, end = Boundary) %>% dplyr::select(chr, start, end)

    # Set up rGREAT job with default parameters
    great_shift <- submitGreatJob(TAD_Frame, request_interval = 1, version = "2.0")

    # Submit the job
    enrichment_table <- getEnrichmentTables(great_shift)

    # Subset to only include vital information
    enrichment_table <- bind_rows(enrichment_table, .id = "source") %>%
      dplyr::select(Ontology = source, Description = name,
                    `P-value` = Hyper_Raw_PValue)

    # Print head organizaed by p-values
    head(enrichment_table %>% dplyr::arrange(`P-value`))

                   Ontology                                          Description     P-value
    1 GO Molecular Function                   gamma-glutamyltransferase activity 0.001802360
    2 GO Biological Process                     glutathione biosynthetic process 0.002927596
    3 GO Cellular Component         anchored to external side of plasma membrane 0.003152529
    4 GO Cellular Component        intrinsic to external side of plasma membrane 0.004051881
    5 GO Biological Process                         peptide biosynthetic process 0.004725996
    6 GO Molecular Function transferase activity, transferring amino-acyl groups 0.004950625

The first column, “Ontology”, is simply the domain from which the corresponding ontology (“Description” column) comes from. Here, we use the default, which is the GO ontologies. For more available ontologies, see the rGREAT vignette. “Description” is the pathway itself. “P-value” is the unadjusted hypergeometric p-value, as output by `rGREAT`. `rGREAT` also provides binomial p-values (Binom_Raw_Pvalue, Binom_Adjp_BH) and adjusted hypergeometric p-values (Hyper_Adjp_BH).

Now we demonstrate how to perform the same analysis but for all boundary types simultaneously. In this case, we use time-varying data.

    # Read in time course data
    data("time_mats")
    # Identifying boundaries
    results <- TimeCompare(time_mats, resolution = 50000)

    # Pulling out the frame of TADs
    TAD_Frame <- results$TAD_Bounds

    # Getting coordinates for TAD boundaries and converting into bed format
    Bound_List <- lapply(unique(TAD_Frame$Category), function(x) {
        TAD_Frame %>% filter((Category == x)) %>% mutate(chr = "chr22") %>%
            dplyr::select(chr, Coordinate) %>%
            mutate(start = Coordinate, end = Coordinate) %>%
            dplyr::select(chr, start, end)
    })

    # Performing rGREAT analysis for each boundary Category
    TAD_Enrich <- lapply(Bound_List, function(x) {
      getEnrichmentTables(submitGreatJob(x, request_interval = 1, version = "2.0"))
    })

    # Name list of data frames to keep track of which enrichment belongs to which
    names(TAD_Enrich) <- unique(TAD_Frame$Category)

    # Bind each category of pathway and create new column for each pathway
    TAD_Enrich <- lapply(names(TAD_Enrich), function(x) {
      bind_rows(lapply(TAD_Enrich[[x]], function(y) {
        y %>% mutate(Category = x)
      }), .id = "source")
    })

    # Bind each boundary category together and pull out important variables
    enrichment_table <- bind_rows(TAD_Enrich) %>%
      dplyr::select(Ontology = source, Description = name,
                    `P-value` = Hyper_Raw_PValue, Category)

    # Get the top enriched pathways
    head(enrichment_table %>% dplyr::arrange(`P-value`))

                   Ontology                                              Description      P-value            Category
    1 GO Biological Process positive regulation of B cell receptor signaling pathway 0.0002254283         Dynamic TAD
    2 GO Molecular Function                               lipid transporter activity 0.0003760684   Highly Common TAD
    3 GO Biological Process          regulation of B cell receptor signaling pathway 0.0004508185         Dynamic TAD
    4 GO Biological Process                               Schwann cell proliferation 0.0006761897 Early Appearing TAD
    5 GO Biological Process                            lipoprotein metabolic process 0.0007139088   Highly Common TAD
    6 GO Molecular Function        peptide-methionine (R)-S-oxide reductase activity 0.0009015354   Highly Common TAD

These columns are the same as the differential analysis but with an extra column, “Category”, indicating the type of time-varying TAD boundary.

# Session Info

    sessionInfo()

    ## R version 4.5.2 (2025-10-31)
    ## Platform: x86_64-pc-linux-gnu
    ## Running under: Ubuntu 24.04.3 LTS
    ##
    ## Matrix products: default
    ## BLAS:   /usr/lib/x86_64-linux-gnu/openblas-pthread/libblas.so.3
    ## LAPACK: /usr/lib/x86_64-linux-gnu/openblas-pthread/libopenblasp-r0.3.26.so;  LAPACK version 3.12.0
    ##
    ## locale:
    ##  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
    ##  [3] LC_TIME=en_US.UTF-8        LC_COLLATE=C
    ##  [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8
    ##  [7] LC_PAPER=en_US.UTF-8       LC_NAME=C
    ##  [9] LC_ADDRESS=C               LC_TELEPHONE=C
    ## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
    ##
    ## time zone: Etc/UTC
    ## tzcode source: system (glibc)
    ##
    ## attached base packages:
    ## [1] stats     graphics  grDevices utils     datasets  methods   base
    ##
    ## other attached packages:
    ## [1] microbenchmark_1.5.0 TADCompare_1.21.0    SpectralTAD_1.27.0
    ## [4] dplyr_1.2.0          BiocStyle_2.39.0
    ##
    ## loaded via a namespace (and not attached):
    ##  [1] tidyselect_1.2.1            farver_2.1.2
    ##  [3] S7_0.2.1                    fastmap_1.2.0
    ##  [5] digest_0.6.39               lifecycle_1.0.5
    ##  [7] cluster_2.1.8.2             HiCcompare_1.33.0
    ##  [9] magrittr_2.0.4              compiler_4.5.2
    ## [11] rlang_1.1.7                 sass_0.4.10
    ## [13] tools_4.5.2                 yaml_2.3.12
    ## [15] data.table_1.18.2.1         knitr_1.51
    ## [17] ggsignif_0.6.4              S4Arrays_1.11.1
    ## [19] PRIMME_3.2-6                DelayedArray_0.37.0
    ## [21] plyr_1.8.9                  RColorBrewer_1.1-3
    ## [23] abind_1.4-8                 BiocParallel_1.45.0
    ## [25] KernSmooth_2.23-26          withr_3.0.2
    ## [27] purrr_1.2.1                 BiocGenerics_0.57.0
    ## [29] sys_3.4.3                   grid_4.5.2
    ## [31] stats4_4.5.2                ggpubr_0.6.2
    ## [33] Rhdf5lib_1.33.0             ggplot2_4.0.2
    ## [35] scales_1.4.0                gtools_3.9.5
    ## [37] SummarizedExperiment_1.41.1 cli_3.6.5
    ## [39] rmarkdown_2.30              generics_0.1.4
    ## [41] reshape2_1.4.5              cachem_1.1.0
    ## [43] rhdf5_2.55.13               stringr_1.6.0
    ## [45] splines_4.5.2               parallel_4.5.2
    ## [47] BiocManager_1.30.27         XVector_0.51.0
    ## [49] matrixStats_1.5.0           vctrs_0.7.1
    ## [51] Matrix_1.7-4                jsonlite_2.0.0
    ## [53] carData_3.0-6               car_3.1-5
    ## [55] IRanges_2.45.0              S4Vectors_0.49.0
    ## [57] rstatix_0.7.3               Formula_1.2-5
    ## [59] maketools_1.3.2             jquerylib_0.1.4
    ## [61] tidyr_1.3.2                 glue_1.8.0
    ## [63] codetools_0.2-20            cowplot_1.2.0
    ## [65] stringi_1.8.7               gtable_0.3.6
    ## [67] GenomicRanges_1.63.1        tibble_3.3.1
    ## [69] pillar_1.11.1               htmltools_0.5.9
    ## [71] Seqinfo_1.1.0               rhdf5filters_1.23.3
    ## [73] R6_2.6.1                    evaluate_1.0.5
    ## [75] lattice_0.22-9              Biobase_2.71.0
    ## [77] backports_1.5.0             pheatmap_1.0.13
    ## [79] broom_1.0.12                bslib_0.10.0
    ## [81] Rcpp_1.1.1                  InteractionSet_1.39.0
    ## [83] gridExtra_2.3               SparseArray_1.11.10
    ## [85] nlme_3.1-168                mgcv_1.9-4
    ## [87] xfun_0.56                   MatrixGenerics_1.23.0
    ## [89] buildtools_1.0.0            pkgconfig_2.0.3