Option set for cuffgffread
A CuffGFFReadOptions
object contains options for the
cuffgffread
function, which filters and converts GFF and GTF files
[1].
creates
a cuffgffreadOpt
= CuffGFFReadOptionsCuffGFFReadOptions
object with the default property values.
CuffGFFReadOptions
requires the Cufflinks Support Package for the Bioinformatics Toolbox™. If the support package is not installed, then the function provides a download
link. For details, see Bioinformatics Toolbox Software Support Packages.
Note
CuffGFFReadOptions
is supported on the Mac and UNIX® platforms only.
sets the object properties using
one or more name-value pair arguments. Enclose each property name in quotes. For example,
cuffgffreadOpt
= CuffGFFReadOptions(Name,Value)cuffgffreadOpt = CuffGFFReadOptions('DiscardSingleExon',true)
discards transcripts spanning a single exon.
specifies optional parameters using the string or character vector
cuffgffreadOpt
= CuffGFFReadOptions(S
)S
.
S
— cuffgffread
optionscuffgffread
options, specified as a string or character vector.
S
must be in the original gffread
option
syntax (prefixed by one or two dashes).
Example: '-U'
AppendDescription
— Flag to add file descriptions to descr
attributefalse
(default) | true
Flag to add file descriptions from sequence files to the
descr
attribute of the output GFF record, specified as
true
or false
. Specify the sequence files using the
SequenceInfo
option.
Example:
true
Data Types: logical
CheckOppositeStrand
— Flag to check opposite strand when checking for in-frame stop codonsfalse
(default) | true
Flag to check opposite strand when checking for in-frame stop codons, specified as true
or false
.
Example:
true
Data Types: logical
CheckPhase
— Flag to adjust coding sequence phasefalse
(default) | true
Flag to adjust coding sequence phase when checking for in-frame stop codons, specified as true
or false
.
Example:
true
Data Types: logical
Cluster
— Flag to cluster input transcripts into locitrue
(default) | false
Flag to cluster the input transcripts into loci, specified as
true
or false
. This option is the same as the
Merge
property, except that it does not collapse fully contained
transcripts with identical introns.
Example:
false
Data Types: logical
CodingOnly
— Flag to discard transcripts with no coding sequencefalse
(default) | true
Flag to discard transcripts with no coding sequence feature (CDS), specified as true
or false
.
Example:
true
Data Types: logical
CollapseContainer
— Flag to collapse fully-contained transcriptsfalse
(default) | true
Flag to collapse fully contained transcripts that are shorter
with fewer introns than the container, specified as true
or
false
. This property applies only when you set Merge
to true
.
Example:
true
Data Types: logical
CollapseFull
— Flag to collapse shorter transcripts overlapping at least 80% with another exonfalse
(default) | true
Flag to collapse shorter transcripts overlapping at least 80%
with another single exon transcript, specified as true
or
false
. This property applies only when you set Merge
to true
.
Example:
true
Data Types: logical
CoordinateRange
— Genomic range to filter transcriptsGenomic range to filter transcripts, specified as a string or character vector. The format must be "[[<strand>]<chr>:]<start>..<end>"
, where start
and end
are genomic positions, chr
is an optional chromosome or contig name, and an optional strand
('+'
or '-'
).
Example:
“+NC_000912.1:4821..7340”
Data Types: char
| string
DiscardInvalidCDS
— Flag to ignore mRNA transcripts either lacking start or stop codon or having in-frame stop codonfalse
(default) | true
Flag to ignore mRNA transcripts either lacking a start or stop codon or having an in-frame stop codon, specified as true
or false
.
Example:
true
Data Types: logical
DiscardNonCanonicalSplice
— Flag to ignore multiexon mRNA transcripts that have intron with noncanonical splice sequencefalse
(default) | true
Flag to ignore multiexon mRNA transcripts that have an intron
with a noncanonical splice sequence, specified as true
or
false
. A noncanonical splice sequence is any splice sequence other than
"GT-AG"
, "CG-AG"
, or
"AT-AC"
.
Example:
true
Data Types: logical
DiscardSingleExon
— Flag to ignore transcripts spanning single exonfalse
(default) | true
Flag to ignore transcripts spanning a single exon, specified as true
or false
.
Example:
true
Data Types: logical
DiscardTerminatedCDS
— Flag to ignore transcripts with in-frame stop codonfalse
(default) | true
Flag to ignore transcripts with an in-frame stop codon, specified as true
or false
.
Example:
true
Data Types: logical
ExtraCommand
— Additional commands""
(default) | character vector | stringThe commands must be in the native syntax (prefixed by one or two dashes). Use this option to apply undocumented flags and flags without corresponding MATLAB properties.
When the software converts the original flags to MATLAB properties, it stores any unrecognized flags in this property.
Example: "-E"
Data Types: char
| string
FastaCDSFile
— Name of file to save spliced coding sequencesName of a file to save the spliced coding sequences in the FASTA format, specified as a string or character vector.
Example:
"splicedCoding.FASTA"
Data Types: char
| string
FastaExonsFile
— Name of file to save spliced exonsName of a file to save the spliced exons in the FASTA format, specified as a string or character vector.
Example:
"splicedExon.FASTA"
Data Types: char
| string
FastaProteinFile
— Name of file to save protein translation of coding sequencesName of a file to save the protein translation of coding sequences in the FASTA format, specified as a string or character vector.
Example:
"translated.FASTA"
Data Types: char
| string
FirstExonOnly
— Flag to parse additional attributes only from first exonfalse
(default) | true
Flag to parse additional attributes only from the first exon, specified as true
or false
.
Example: true
Data Types: logical
ForceExons
— Flag to list lowest-level GFF features as exon featuresfalse
(default) | true
Flag to list the lowest-level GFF features as exon features in
the output file, specified as true
or false
.
Example:
true
Data Types: logical
FullyContained
— Flag to discard transcripts not contained fullyfalse
(default) | true
Flag to discard transcripts not contained fully within the
range, specified as true
or false
. Specify the range using
the CoordinateRange
option.
Example:
true
Data Types: logical
GTFOutput
— Flag to output GTF-format transcript filesfalse
(default) | true
Flag to output GTF-format transcript files, specified as
true
or false
.
Example:
true
Data Types: logical
IncludeAll
— Flag to use all object propertiesfalse
(default) | true
Flag to include all the object properties with the
corresponding default values when converting to the original options syntax, specified as
true
or false
. You can convert the properties to the
original syntax prefixed by one or two dashes (such as '-d 100 -e 80'
) by
using getCommand
. The
default value false
means that when you call
getCommand(optionsObject)
, it converts only the specified properties.
If the value is true
, getCommand
converts all available
properties, with default values for unspecified properties, to the original syntax.
Note
If you set IncludeAll
to true
, the software
translates all available properties, with default values for unspecified properties. The
only exception is that when the default value of a property is NaN
,
Inf
, []
, ''
, or
""
, then the software does not translate the corresponding
property.
Example: true
Data Types: logical
MaxIntronLength
— Maximum intron length for transcript to be included in outputInf
(default) | positive integerMaximum intron length for a transcript to include in the output
file, specified as a positive integer. Inf
, the default value, sets no limit
on the intron length.
Example:
500
Data Types: double
Merge
— Flag to merge transcripts to locifalse
(default) | true
Flag to merge transcripts into loci by collapsing transcripts with identical introns, specified as true
or false
.
Example:
true
Data Types: logical
MergeCloseExons
— Flag to merge exons into single exonfalse
(default) | true
Flag to merge exons into a single exon when separated by fewer than 4 base-pair introns, specified as true
or false
.
Example:
true
Data Types: logical
MergeInfoFile
— Name of file to save information on duplicates when mergingName of a file to save information on duplicates when merging,
specified as a string or character vector. This property applies only when you set
Merge
to true
.
Example:
"duplicates.txt"
Data Types: char
| string
PreserveAttributes
— Flag to retain all attributes in outputfalse
(default) | true
Flag to retain all attributes in the output file, specified as true
or false
.
Example:
true
Data Types: logical
Pseudo
— Flag to filter out records containing "pseudo"true
(default) | false
Flag to filter out records containing the word "pseudo,"
specified as true
or false
.
Example:
false
Data Types: logical
ReplacementTable
— Name of file containing replacement tableName of a file containing a replacement table, specified as a string or character vector. The table must have two columns, where the first column contains the original transcript IDs and the second column contains the new transcript IDs. An example table follows.
origTranscript1 | newTranscript1 |
origTranscript2 | newTranscript2 |
origTranscript3 | newTranscript3 |
If you provide a replacement table, the function replaces the transcript IDs found in the first column with the new transcripts IDs from the second column and filters out those transcripts not found.
Example:
"replaceTbl.txt"
Data Types: char
| string
SequenceFile
— Name of FASTA-format file containing genomic sequencesName of a FASTA-format file containing genomic sequences for all input mappings, specified as a string or character vector.
Example:
"seqs.fasta"
Data Types: char
| string
SequenceInfo
— Name of tab-delimited file with additional information on input sequenceName of a tab-delimited file with additional information on
each input sequence, specified as a string or character vector. This file must have three
columns: a sequence name column, a sequence length column, and a sequence description column. If
AppendDescription
is true
, the sequence description
is included as an attribute in the output GFF file.
Example:
"seqinfo.txt"
Data Types: char
| string
UrlDecode
— Flag to decode URL-encoded characters in attribute namesfalse
(default) | true
Flag to decode url-encoded characters in attribute names,
specified as true
or false
. For instance,
"transcript%20description" is decoded to "transcript description".
Example:
true
Data Types: logical
UseEnsemblConversion
— Flag to use GTF-to-GFF3 conversion method from Ensemblfalse
(default) | true
Flag to use the GTF-to-GFF3 conversion method from Ensembl, specified as true
or false
.
Example:
true
Data Types: logical
UseNonTranscript
— Flag to include nontranscript GFF records in output filefalse
(default) | true
Flag to include nontranscript GFF records in the output file, specified as true
or false
.
Example:
true
Data Types: logical
UseTrackName
— Flag to use track name in second column of GFF output linefalse
(default) | true
Flag to use the track name in the second column of the GFF output line, specified as true
or false
.
Example:
true
Data Types: logical
Version
— Supported versionThis property is read-only.
Supported version of the original cufflinks software, returned as a string.
Example: "2.2.1"
Data Types: string
WriteCoordinates
— Flag to write exon coordinates projected onto spliced sequencefalse
(default) | true
Flag to write the exon coordinates projected onto the spliced
sequence, specified as true
or false
. This property
applies only when FastaExonsFile
or FastaCDSFile
is
specified.
Example:
true
Data Types: logical
getCommand | Translate object properties to original options syntax |
getOptionsTable | Return table with all properties and equivalent options in original syntax |
Create a CuffGFFReadOptions
object with the default values.
opt = CuffGFFReadOptions;
Create an object using name-value pairs.
opt2 = CuffGFFReadOptions('DiscardSingleExon',true,'FastaExonsFile','exons.fa');
Create an object by using the original syntax.
opt3 = CuffGFFReadOptions('-U -w exons.fa')
Convert a GTF file to a GFF file while retaining all attributes.
cuffgffread('gyrAB.gtf','gyrABOut.gff','PreserveAttributes',true)
You can also set the options using an object. For instance, specify the output to be in the GTF format.
opt = CuffGFFReadOptions; opt.GTFOutput = true; opt.PreserveAttributes = true; cuffgffread('gyrAB.gtf','gyrABOut.gtf',opt);
Once you have the options object, you can retrieve the equivalent original options for all object properties using getOptionsTable
.
getOptionsTable(opt)
ans = 33×3 table PropertyName FlagName FlagShortName ___________________________ ________________ _____________ AppendDescription 'AppendDescription' '-A' '' CheckOppositeStrand 'CheckOppositeStrand' '-B' '' CheckPhase 'CheckPhase' '-H' '' Cluster 'Cluster' '--cluster-only' '' CodingOnly 'CodingOnly' '-C' '' CollapseContainer 'CollapseContainer' '-K' '' CollapseFull 'CollapseFull' '-Q' '' CoordinateRange 'CoordinateRange' '-r' '' DiscardInvalidCDS 'DiscardInvalidCDS' '-J' '' DiscardNonCanonicalSplice 'DiscardNonCanonicalSplice' '-N' '' DiscardSingleExon 'DiscardSingleExon' '-U' '' DiscardTerminatedCDS 'DiscardTerminatedCDS' '-V' '' FastaCDSFile 'FastaCDSFile' '-x' '' FastaExonsFile 'FastaExonsFile' '-w' '' FastaProteinFile 'FastaProteinFile' '-y' '' FirstExonOnly 'FirstExonOnly' '-G' '' ForceExons 'ForceExons' '--force-exons' '' FullyContained 'FullyContained' '-R' '' GTFOutput 'GTFOutput' '-T' '' MaxIntronLength 'MaxIntronLength' '-i' '' Merge 'Merge' '--merge' '-M' MergeCloseExons 'MergeCloseExons' '-Z' '' MergeInfoFile 'MergeInfoFile' '-d' '' PreserveAttributes 'PreserveAttributes' '-F' '' Pseudo 'Pseudo' '--no-pseudo' '' ReplacementTable 'ReplacementTable' '-m' '' SequenceFile 'SequenceFile' '-g' '' SequenceInfo 'SequenceInfo' '-s' '' UrlDecode 'UrlDecode' '-D' '' UseEnsemblConversion 'UseEnsemblConversion' '-L' '' UseNonTranscript 'UseNonTranscript' '-O' '' UseTrackName 'UseTrackName' '-t' '' WriteCoordinates 'WriteCoordinates' '-W' ''
[1] Trapnell, Cole, Brian A Williams, Geo Pertea, Ali Mortazavi, Gordon Kwan, Marijke J van Baren, Steven L Salzberg, Barbara J Wold, and Lior Pachter. “Transcript Assembly and Quantification by RNA-Seq Reveals Unannotated Transcripts and Isoform Switching during Cell Differentiation.” Nature Biotechnology 28, no. 5 (May 2010): 511–15.