findNeighbors {annotate} | R Documentation

## A function to locate neighboring genes within a defined range around a target gene represented by a LocusLink id

### Description

Give a data package with mappings between LocsuLink ids and their locations on chromosomes, this function locates genes that are within a defined range on a given chromosome. If a LocusLink id is passed as one of the arguments, genes located will be neighbors to the gene represented by the LocusLink id within a defined range on the chromosome the target gene resides

### Usage

    findNeighbors(chrLoc, llID, chromosome, upBase, downBase, mergeOrNot = TRUE)
    checkArgs(llID, chromosome, upBase, downBase)
    findChr4LL(llID, chrEnv, organism)
    getValidChr(organism)
    getBoundary(loc, base, lower = TRUE)
    weightByConfi(foundLLs)

### Arguments

`chrLoc` | `chrLoc` a character string for the name of the data package that contains mappings between LocusLink ids and their locations on chromosomes. For each chromosome, there assumed to be mappings for the start and end locations of genes represented by LocusLink ids. The data package needs to be built using chrLocPkgBuilder of AnnBuilder
`llID` | `llID` a character string for the LocusLink id representing a gene whose nighbors are sought. llID can be missing
`chromosome` | `chromosome` a character string for the number of the chromosome of interest. chromosome is only required for locating genes within a range on the chromosome

`upBase`

`upBase` a numeric or character string for the number of base pairs that defines the upper limit of the range to locate genes. If neighbors of a given gene is sourght, the value will be the distance in number of base pairs from the target gene upstream, to which search for genes will be conducted. Otherwise, the value will be the upper limit in number of base pairs from the p arm, to which search for genes will be conducted

`downBase`

`downBase` a numeric or character string for the number of base pairs that defines the lower limit of the range to locat gene. If neighbors of a given gene is sourght, the value will be the distance in number of base pairs from the target gene downstream, to which search for genes will be conducted. Otherwise, the value will be the lower limit in number of base pairs from the p arm, to which search for genes will be conducted

`organism` | `organism` a character string for the name of the organism of interest
`chrEnv` | `chrEnv` an environment object with keys for LocusLink ids and values for the chromosomes where genes reside
`loc` | `loc` a numeric of character string for the chromosomal location of gene of interest
`base` | `base` either a `downBase` or `upBase`
`lower` | `lower` a boolean indicating whether the lower or upper boundary of search limit is sought
`mergeOrNot` | `mergeOrNot` a boolean to indicate whether gene found up and down streams will be merged (TRUE)
`foundLLs` | `foundLLs` a vector of character strings for LocusLink ids

### Details

A chrLoc data package can be created using function chrLocPkgBuilder of AnnBuilder, in which locusLink ids are mapped to location data on individual chromosomes.

Genes are considered to be neighbors to a given target gene or within a given range when the transcription of genes start and end within the given range.

findNeighbors, checkArgs, findChr4LL, getValidChr, and getBoundary are accessory functions called by findNeighbors and may not have real values outside.

### Value

The function returns a list of named vectors. The length of the list is one when genes in a given rwgion are sought but varies depending on whether a given gene can be mapped to one or more chromosomes when neighoring genes of a target gene are sought. Names of vector can be "Confident" when a gene can be confidently placed on a chromosome or "Unconfident" when a gene can be placed on a chromosome but its exact location can not be determined with great confidence.

### Note

This function is part of the Bioconductor project at Dana-Farber Cancer Institute to provide bioinformatics functionalities through R

### Author(s)

Jianhua Zhang

### References

http://www.genome.ucsc.edu/goldenPath/

### Examples

    if(require("humanCHRLOC")){
       findNeighbors("humanCHRLOC", "10228", upBase = 600000, downBase = 600000)
    }else{
       print("Can not find neighbors without the required data package")
    }

[Package annotate version 1.8.0 Index]