Learn R Programming
expands (version 1.5)

buildMultiSamplePhylo: Relations between inter- and intra-sample subpopulations

Description

Predicts phylogenetic relations between subpopulations from subpopulation specific copy number and point mutation profiles, while including information about sample origin of each subpopulation. This function differs from buildPhylo in that it integrates the subpoulations predicted in multiple, geographically distinct tumor-samples into one common phylogeny and in that it includes point mutations in addition to copy number variations for phylogeny reconstruction.

Usage

buildMultiSamplePhylo(samGr,out, treeAlgorithm = "bionjs", keepAmbigSeg = FALSE, plotF=1)

Arguments

samGr
List with three fields: cbs - Matrix in which each row corresponds to a copy number segment. CBS is typically the output of a circular binary segmentation algorithm. Columns in CBS must be labeled and must include chr, startpos, endpos and CN_Esti
out
Prefix of file to which multi-sample phylogeny will be saved.
treeAlgorithm
Neighbor joining algorithm used for phylogeny reconstruction (from library ape). Options: bionjs (default), njs.
keepAmbigSeg
Input parameter for called function: assignQuantityToSP.
plotF
Option for displaying a visual representation of the phylogenetic tree (0 - no display; 1 - display). Default: 1.

Value

  • An object of class "phylo" (library ape).

Details

Reconstructs phylogenetic relationships between subpopulations using the BIONJ algorithm of Gascuel. Pairwise distances between subpopulations are calculated as: $C:=(cnv_{i=j} + snv_{i=j})/(cnv_{ij}+snv_{ij})$, where $cnv_{i=j}$ is the number of copy number segments for which subpopulations i and j have the same copy number; $snv_{i=j}$ is the number of point mutations for which subpopulations i and j have the same mutation status and $cnv_{ij}, snv_{ij}$ are the total number of copy number segments and mutations respectively, for which both subpopulations have available information. Subpopulations with insufficient ploidy and point mutations information are excluded from phylogeny.

See Also

buildPhylo