getGenes

Class: GTFAnnotation

Return table of unique genes in GTFAnnotation object

expand all in page

Syntax

genes = getGenes(AnnotObj)
genes = getGenes(AnnotObj,'Reference',R)
genes = getGenes(AnnotObj,'Gene',G)
genes = getGenes(AnnotObj,'Transcript',T)

Description

genes = getGenes(AnnotObj) returns genes, a table of genes referenced by exons in AnnotObj.

genes = getGenes(AnnotObj,'Reference',R) returns the gene(s) that belong to the reference(s) specified by R.

genes = getGenes(AnnotObj,'Gene',G) returns the gene(s) specified by G.

genes = getGenes(AnnotObj,'Transcript',T) returns the gene(s) that contains the transcript(s) specified by T.

Input Arguments

expand all

GTF annotation, specified as a GTFAnnotation object.

Names of reference sequences, specified as a character vector, string, string vector, cell array of character vectors, or categorical array.

The names must come from the Reference field of AnnotObj. If a name does not exist, the function provides a warning and ignores it.

Names of genes, specified as a character vector, string, string vector, cell array of character vectors, or categorical array.

The names must come from the Gene field of AnnotObj. If a name does not exist, the function provides a warning and ignores the name.

Names of transcripts, specified as a character vector, string, string vector, cell array of character vectors, or categorical array.

The names must come from the Transcript field of AnnotObj. If a name does not exist, the function gives a warning and ignores the name.

Output Arguments

expand all

Genes referenced by exons in AnnotObj, returned as a table. The table contains the following variables for each gene.

Variable NameDescription
GeneIDCell array of character vectors containing gene IDs as listed in AnnotObj, obtained from the Gene field of AnnotObj.
GeneNameCell array of character vectors containing gene names, obtained from the Attributes field of AnnotObj. This cell array can contain empty character vectors if the corresponding gene names are not found in Attributes.
ReferenceCategorical array representing the names of reference sequences to which the genes belong, obtained from the Reference field of AnnotObj.
StartStart location of the first exon in each gene.
StopStop location of the last exon in each gene.
StrandCategorical array containing the strand of each gene.
NumTranscriptsInteger array listing the number of transcripts in each gene.

Examples

expand all

Open Live Script

Create a GTFAnnotation object from a GTF-formatted file.

obj = GTFAnnotation('hum37_2_1M.gtf');

Retrieve unique reference names. In this case, there is only one reference sequence, which is chromosome 2 (chr2).

ref = getReferenceNames(obj)
ref = 1x1 cell array
    {'chr2'}

Get a table of all genes which belong to chr2.

genes = getGenes(obj,'Reference',ref)
genes=28×7 table
        GeneID         GeneName     Reference    Start      Stop     Strand    NumTranscripts
    ______________    __________    _________    ______    ______    ______    ______________

    {'uc010yim.1'}    {0x0 char}      chr2        41609     46385      -             1       
    {'uc002qvu.2'}    {0x0 char}      chr2       218138    249852      -             1       
    {'uc002qvv.2'}    {0x0 char}      chr2       218138    256690      -             1       
    {'uc002qvw.2'}    {0x0 char}      chr2       218138    260702      -             1       
    {'uc002qvx.2'}    {0x0 char}      chr2       218138    264068      -             1       
    {'uc002qvy.2'}    {0x0 char}      chr2       218138    264068      -             1       
    {'uc002qvz.2'}    {0x0 char}      chr2       218138    264392      -             1       
    {'uc002qwa.2'}    {0x0 char}      chr2       218138    264743      -             1       
    {'uc010ewe.2'}    {0x0 char}      chr2       218138    264810      -             1       
    {'uc002qwb.2'}    {0x0 char}      chr2       239563    242178      -             1       
    {'uc002qwc.1'}    {0x0 char}      chr2       243503    262786      -             1       
    {'uc002qwd.2'}    {0x0 char}      chr2       264869    272481      +             1       
    {'uc002qwe.3'}    {0x0 char}      chr2       264869    273148      +             1       
    {'uc002qwg.2'}    {0x0 char}      chr2       264869    278280      +             1       
    {'uc002qwh.2'}    {0x0 char}      chr2       264869    278280      +             1       
    {'uc002qwf.2'}    {0x0 char}      chr2       264869    278280      +             1       
      ⋮

See Also

| | | | | | | | | | |

Topics

External Websites

Bioinformatics Toolbox Documentation

Support