SNVtest2: Test for tumors from 2 different sites

SNVtest2R Documentation

Test for tumors from 2 different sites

Description

This functions performs clonality testing of 2 tumors from the same patient that come from different types or organs. The test uses profiles of somatic mutations. The null hypothesis is that two tumors that come from different sites are independent. There are two alternative hypotheses: that they are clonal and from site 1, and clonal from site 2. This function calculates conditional maximum likelihood ratio relying only on loci where at least one of the tumors has a mutation, and p-value is calculated under the reference distribution under the hypothesis of independence.

Usage

SNVtest2(tumor1, tumor2, pfreq, nrep = 999)

Arguments

tumor1

Vector of the binary mutation calls from tumor 1, where 0 denotes no mutation, 1 denotes a mutation. Mutations should be in the same order as frequencies in pfreq.

tumor2

Vector of the binary mutation calls from tumor 2, where 0 denotes no mutation, 1 denotes a mutation. Mutations should be in the same order as frequencies in pfreq.

pfreq

Two column matrix of marginal frequencies of mutations in two sites known a priori. These can be obtained from TCGA or similar databases or calculated using function "get.mutation.frequencies" in this package.

nrep

Number of simulations used for generating the reference distribution under the hypothesis of independence.

Details

Test is related to the one described by "SNVtest" but it assumes that under null hypothesis two tumors come from 2 different sites with different marginal proabilities, thus the reference distribution is generated from 2 sets of frequencies. The p-value is significant when null hypothesis is rejected. It might be rejected sometimes in the absence of matches if the observed mutational profiles are unlikely to come from the declared tumor sites.

Value

The output is a vector with 6 values: "n.match", "n.site1only", "n.site2only", "xi.site1", "xi.site2", "p.value"

n.match

Number of matches between two tumors

n.site1only

Number of mutations in the first tumor only

n.site2only

Number of mutations in the second tumor only

xi.site1

Maximum likelhood estimate of Ksi, parameter of the likelihood representing clonality strength, for the alternative hypothesis that two clonal tumors come from site 1. Value close to 0 indicates independence, value close to 1 indicates perfect concordance in mutational profiles.

xi.site2

Maximum likelhood estimate of Ksi, parameter of the likelihood representing clonality strength, for the alternative hypothesis that two clonal tumors come from site 2. Value close to 0 indicates independence, value close to 1 indicates perfect concordance in mutational profiles.

p.value

p-value calculated using the null distribution generated using two independet tumors generated from two different sites

Author(s)

Irina Ostrovnaya ostrovni@mskcc.org

References

Ostrovnaya I, Mauguen A, Seshan V, Begg CB "Testing Tumors from Different Anatomic Sites for Clonal Relatedness Using Somatic Mutation Data", submitted

Examples



#Example of artificially generated independent tumors  from marginal mutation frequencies p1 and p2
n<-100
p1<-runif(n)/10
p2<-runif(n)/10
x1<-as.numeric(runif(n)<=p1)
x2<-as.numeric(runif(n)<=p2)
SNVtest2(x1,x2,cbind(p1,p2))



IOstrovnaya/Clonality documentation built on July 22, 2023, 4:16 a.m.