get_regulatory_summary: Return summary of gene regulation.

In regutools: regutools: an R package for data extraction from RegulonDB

Description

This function takes the output of get_gene_regulators() with format multirow, onerow or table, or a vector with genes and retrieves information about the TFs and their regulated genes

Usage

get_regulatory_summary(regulondb, gene_regulators)

Arguments

regulondb	A regulondb() object.
gene_regulators	Result from get_gene_regulators() or vector of genes

Value

A data frame with the following columns:

• The name or gene of TF

• Regulated Genes per TF

• Percent of regulated genes per TF

• positive, negative or dual regulation

• Name(s) of regulated genes

Author(s)

Carmina Barberena Jonas, Jesús Emiliano Sotelo Fonseca, José Alquicira Hernández, Joselyn Chávez

Examples

## Connect to the RegulonDB database if necessary
if (!exists("regulondb_conn")) regulondb_conn <- connect_database()

## Build the regulon db object
e_coli_regulondb <-
    regulondb(
        database_conn = regulondb_conn,
        organism = "E.coli",
        database_version = "1",
        genome_version = "1"
    )

## Get the araC regulators
araC_regulation <-
    get_gene_regulators(
        e_coli_regulondb,
        genes = c("araC"),
        format = "multirow",
        output.type = "TF"
    )

## Summarize the araC regulation
get_regulatory_summary(e_coli_regulondb, araC_regulation)

## Retrieve summary of genes 'araC' and 'modB'
get_regulatory_summary(e_coli_regulondb,
    gene_regulators = c("araC", "modB")
)

## Obtain the summary for 'ECK120000050' and 'modB'
get_regulatory_summary(e_coli_regulondb,
    gene_regulators = c("ECK120000050", "modB")
)

regutools documentation built on Dec. 20, 2020, 2 a.m.