## 3.0 contains all regulons that are considered for AUCell scoring
load("./int/3.0_regulons_forAUCell.RData")

## Define the gene you would like to find
gene_of_interest <- ""


## Iterates over all regulon names in regulons and searches the gene of interest inside the regulon.
for(reg in names(regulons)){
  if(gene_of_interest %in% regulons[[reg]]){
  ## Prints regulon name if gene of interest was found in regulon
    print(reg)
  }
}