getGenes

Return table of unique genes in GTFAnnotation object

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 one or more genes that belong to one or more references specified by R.

genes = getGenes(AnnotObj,"Gene",G) returns one or more genes specified by G.

genes = getGenes(AnnotObj,"Transcript",T) returns one or more genes that contains one or more transcripts specified by T.

Examples

Retrieve Genes from a GTF-formatted File

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

Input Arguments

AnnotObj — GTF annotation

GTFAnnotation object

GTF annotation, specified as a GTFAnnotation object.

R — Names of reference sequences

character vector | string | string vector | cell array of character vectors | categorical array

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.

Data Types: char | string | cell | categorical

G — Names of genes

character vector | string | string vector | cell array of character vectors | categorical array

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.

Data Types: char | string | cell | categorical

T — Names of transcripts

character vector | string | string vector | cell array of character vectors | categorical array

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.

Data Types: char | string | cell | categorical

Output Arguments

genes — Genes referenced by exons in AnnotObj

table

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

Variable Name	Description
GeneID	Cell array of character vectors containing gene IDs as listed in AnnotObj, obtained from the Gene field of AnnotObj.
GeneName	Cell 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.
Reference	Categorical array representing the names of reference sequences to which the genes belong, obtained from the Reference field of AnnotObj.
Start	Start location of the first exon in each gene.
Stop	Stop location of the last exon in each gene.
Strand	Categorical array containing the strand of each gene.
NumTranscripts	Integer array listing the number of transcripts in each gene.

Version History

Introduced in R2014b