chimeraSeqs: A function to generate the nucleotide sequences of a fusion...

Description Usage Arguments Value Author(s) See Also Examples

View source: R/chimeraSeqs.R

Description

A function generating the nucleotide sequences of a chimera.

Usage

1
chimeraSeqs(fset, extend=1000, type="transcripts")

Arguments

fset

A fSet object.

extend

number of nucleotides used to extend a genomic region that is not an annotated gene. Default is 1000 nts

type

Chimera can be build at transcript level, i.e. transcript level will generate multiple fusion transcripts depending of the number of transcripts associated to each of the two genes involved in the fusion.

Value

A DNAStringSet encompassing the fusions generated using all the isoforms for each gene involved in the fusion.The name of each element of the DNAStringSet has the following format: gene1-lengthOfGeneFragment:gene2-lengthOfGeneFragment. In case the fusion junction is located in an intronic sequence, a warning is provided. The presence of a partial intron in the fusion is an indication that the fusion does not generate an active chimeric peptide.

Author(s)

Raffaele A Calogero

See Also

fusionName, chimeraSeqSet

Examples

1
2
3
4
5
6
	tmp <- importFusionData("fusionmap", paste(find.package(package="chimera"),"/examples/mcf7.FMFusionReport", sep=""), org="hg19")
	fusion.names <- fusionName(tmp)
	fusion.names
	myset <- tmp[[13]]
	tmp.seq <- chimeraSeqs(myset, type="transcripts")
	#writeXStringSet(tmp.seq, paste(sub(":","_",fusion.names[[13]]),".fa",sep=""), format="fasta")

chimera documentation built on Nov. 8, 2020, 5:16 p.m.