CuffCompareOptions

Option set for cuffcompare

Description

A CuffCompareOptions object specifies options for the cuffcompare function, which compares assembled transcripts across several experiments [1].

Creation

Description

example

cuffcompareOpt = CuffCompareOptions creates a CuffCompareOptions object with the default property values.

CuffCompareOptions requires the Cufflinks Support Package for Bioinformatics Toolbox™. If the support package is not installed, then the function provides a download link.

Note

CuffCompareOptions is supported on the Mac and UNIX® platforms only.

cuffcompareOpt = CuffCompareOptions(Name,Value) sets the object properties using one or more name-value pair arguments. Enclose each property name in quotes. For example, cuffcompareOpt = CuffCompareOptions('SupressMapFiles',true) prevents the creation of .tmap and .refmap files.

cuffcompareOpt = CuffCompareOptions(S) specifies optional parameters using the string or character vector S.

Input Arguments

expand all

cuffcompare options, specified as a string or character vector. S must be in the original cuffcompare option syntax (prefixed by one or two dashes).

Example: '-d 100 -e 80'

Properties

expand all

Prefix for consensus transcript names in the output combined.gtf file, specified as a string or character vector. This option must be a string or character vector with a non-zero length.

Example: "consensusTs"

Data Types: char | string

Flag to ignore intron-redundant transfrags if they have the same 5' end but different 3' ends, specified as true or false.

Example: true

Data Types: logical

Flag to discard single-exon transfrags and reference transcripts, specified as true or false.

Example: true

Data Types: logical

Flag to discard single-exon reference transcripts, specified as true or false.

Example: true

Data Types: logical

Additional commands, specified as a string or character vector. The commands must be in the original syntax (prefixed by one or two dashes). Use this option to apply undocumented flags and flags without corresponding MATLAB properties. When the function converts the original flags to MATLAB properties, it stores any unrecognized flags in this option.

Example: "--library-type fr-secondstrand"

Data Types: char | string

Name of the text file containing a list of GTF files to process, specified as a string or character vector. The file must contain one GTF file path per line. You can use this option as an alternative to passing an array of file names to cuffcompare.

Example: "gtfManifestFile.txt"

Data Types: char | string

Flag to treat input GTF files as GFF files, specified as true or false. Use this option when the input GFF or GTF files are not produced by cufflinks.

Example: true

Data Types: logical

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.

Example: true

Data Types: logical

Flag to include transfrags contained by other transfrags in the same locus in the output combined.gtf, specified as true or false. By default, cuffcompare does not include these contained transfrags. If the value is true, the contained transfrags include a contained_in attribute indicating the first container transfrag found.

Example: true

Data Types: logical

Number of bases from the free ends of terminal exons to use when assessing exon accuracy, specified as a positive integer.

Example: 80

Data Types: double

Number of bases to use for grouping transcript start sites, specified as a positive integer.

Example: 90

Data Types: double

Prefix for cuffcompare output files, specified as a string or character vector. This option must be a string or character vector with a non-zero length.

Example: "cuffcompareOut"

Data Types: char | string

Name of the GTF or GFF file containing reference transcripts to compare to each sample, specified as a string or character vector. If you provide a file, the function compares each sample to the references in the file and marks isoforms as overlapping, matching, or novel. The function stores these tags in the output files .refmap and .tmap files.

Example: "references.gtf"

Data Types: char | string

Name of directory containing FASTA sequences to classify input transcripts as repeats, specified as a string or character vector. The directory must contain FASTA-format files with the underlying genomic sequences and contain one FASTA file per reference. Name each FASTA file after the chromosome with the extension .fa or .fasta.

Example: "./SequenceDirectory/"

Data Types: char | string

Flag to consider only reference transcripts that overlap with any of the input transfrags, specified as true or false. If the value is true:

  • The function ignores any reference transcripts that do not overlap with any of the input transfrags.

  • You must also specify the ReferenceGTF option.

Example: true

Data Types: logical

Flag to consider only input transcripts that overlap with any of the reference transcripts, specified as true or false. If the value is true:

  • The function ignores any input transcripts that do not overlap with any of the reference transcripts and reports no novel loci.

  • You must also specify the ReferenceGTF option.

Example: true

Data Types: logical

Flag to prevent the creation of .tmap and .refmap files, specified as true or false. Set the value to true to prevent the function from generating the files.

Example: true

Data Types: logical

This property is read-only.

Supported version of the original cufflinks software, returned as a string.

Example: "2.2.1"

Data Types: string

Object Functions

getCommandTranslate object properties to original options syntax
getOptionsTableReturn table with all properties and equivalent options in original syntax

Examples

collapse all

Create a CuffCompareOptions object with the default values.

opt = CuffCompareOptions;

Create an object using name-value pairs.

opt2 = CuffCompareOptions('GenericGFF',true,'MaxAccuracyRange',80)

Create an object using the original syntax.

opt3 = CuffCompareOptions('-d 100 -e 80')

Create a CufflinksOptions object to define cufflinks options, such as the number of parallel threads and the output directory to store the results.

cflOpt = CufflinksOptions;
cflOpt.NumThreads = 8;
cflOpt.OutputDirectory = "./cufflinksOut";

The SAM files provided for this example contain aligned reads for Mycoplasma pneumoniae from two samples with three replicates each. The reads are simulated 100bp-reads for two genes (gyrA and gyrB) located next to each other on the genome. All the reads are sorted by reference position, as required by cufflinks.

sams = ["Myco_1_1.sam","Myco_1_2.sam","Myco_1_3.sam",...
        "Myco_2_1.sam", "Myco_2_2.sam", "Myco_2_3.sam"];

Assemble the transcriptome from the aligned reads.

[gtfs,isofpkm,genes,skipped] = cufflinks(sams,cflOpt);

gtfs is a list of GTF files that contain assembled isoforms.

Compare the assembled isoforms using cuffcompare.

stats = cuffcompare(gtfs);

Merge the assembled transcripts using cuffmerge.

mergedGTF = cuffmerge(gtfs,'OutputDirectory','./cuffMergeOutput');

mergedGTF reports only one transcript. This is because the two genes of interest are located next to each other, and cuffmerge cannot distinguish two distinct genes. To guide cuffmerge, use a reference GTF (gyrAB.gtf) containing information about these two genes. If the file is not located in the same directory that you run cuffmerge from, you must also specify the file path.

gyrAB = which('gyrAB.gtf');
mergedGTF2 = cuffmerge(gtfs,'OutputDirectory','./cuffMergeOutput2',...
			'ReferenceGTF',gyrAB);

Calculate abundances (expression levels) from aligned reads for each sample.

abundances1 = cuffquant(mergedGTF2,["Myco_1_1.sam","Myco_1_2.sam","Myco_1_3.sam"],...
                        'OutputDirectory','./cuffquantOutput1');
abundances2 = cuffquant(mergedGTF2,["Myco_2_1.sam", "Myco_2_2.sam", "Myco_2_3.sam"],...
                        'OutputDirectory','./cuffquantOutput2');

Assess the significance of changes in expression for genes and transcripts between conditions by performing the differential testing using cuffdiff. The cuffdiff function operates in two distinct steps: the function first estimates abundances from aligned reads, and then performs the statistical analysis. In some cases (for example, distributing computing load across multiple workers), performing the two steps separately is desirable. After performing the first step with cuffquant, you can then use the binary CXB output file as an input to cuffdiff to perform statistical analysis. Because cuffdiff returns several files, specify the output directory is recommended.

isoformDiff = cuffdiff(mergedGTF2,[abundances1,abundances2],...
                      'OutputDirectory','./cuffdiffOutput');

Display a table containing the differential expression test results for the two genes gyrB and gyrA.

readtable(isoformDiff,'FileType','text')
ans =

  2×14 table

        test_id            gene_id        gene              locus             sample_1    sample_2    status     value_1       value_2      log2_fold_change_    test_stat    p_value    q_value    significant
    ________________    _____________    ______    _______________________    ________    ________    ______    __________    __________    _________________    _________    _______    _______    ___________

    'TCONS_00000001'    'XLOC_000001'    'gyrB'    'NC_000912.1:2868-7340'      'q1'        'q2'       'OK'     1.0913e+05    4.2228e+05          1.9522           7.8886      5e-05      5e-05        'yes'   
    'TCONS_00000002'    'XLOC_000001'    'gyrA'    'NC_000912.1:2868-7340'      'q1'        'q2'       'OK'     3.5158e+05    1.1546e+05         -1.6064          -7.3811      5e-05      5e-05        'yes'   

You can use cuffnorm to generate normalized expression tables for further analyses. cuffnorm results are useful when you have many samples and you want to cluster them or plot expression levels for genes that are important in your study. Note that you cannot perform differential expression analysis using cuffnorm.

Specify a cell array, where each element is a string vector containing file names for a single sample with replicates.

alignmentFiles = {["Myco_1_1.sam","Myco_1_2.sam","Myco_1_3.sam"],...
                  ["Myco_2_1.sam", "Myco_2_2.sam", "Myco_2_3.sam"]}
isoformNorm = cuffnorm(mergedGTF2, alignmentFiles,...
                      'OutputDirectory', './cuffnormOutput');

Display a table containing the normalized expression levels for each transcript.

readtable(isoformNorm,'FileType','text')
ans =

  2×7 table

      tracking_id          q1_0          q1_2          q1_1          q2_1          q2_0          q2_2   
    ________________    __________    __________    __________    __________    __________    __________

    'TCONS_00000001'    1.0913e+05         78628    1.2132e+05    4.3639e+05    4.2228e+05    4.2814e+05
    'TCONS_00000002'    3.5158e+05    3.7458e+05    3.4238e+05    1.0483e+05    1.1546e+05    1.1105e+05

Column names starting with q have the format: conditionX_N, indicating that the column contains values for replicate N of conditionX.

References

[1] Trapnell, C., B. Williams, G. Pertea, A. Mortazavi, G. Kwan, J. van Baren, S. Salzberg, B. Wold, and L. Pachter. 2010. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechnology. 28:511–515.

See Also

|

External Websites

Introduced in R2019a