kallisto icon indicating copy to clipboard operation
kallisto copied to clipboard

allele-specific expression

Open afadda1 opened this issue 3 years ago • 2 comments

Hi. i have provided the two different alleles in my reference sequence. however, i'm not sure if kallisto can handle this properly. also, what happens in the case of homozygosity? is the tpm reported for each allele is the tpm for the gene (transcripts from the 2 chromosomes)?

afadda1 avatar Jan 11 '22 17:01 afadda1

Yes, kallisto can do allele-specific expression; just put the two transcripts (one for each of the two alleles) in the transcriptome index. You'll get a TPM value for each transcript. If only one allele is expressed (as in the case of homozygosity), the other allele will have a much smaller TPM (its TPM may not necessarily be zero though because the abundances are assigned to the alleles based on a probablistic framework). The gene-level TPM should be obtained by summing the TPMs of all transcripts (including alleles) associated with that gene.

Yenaled avatar Jan 12 '22 02:01 Yenaled

thanks a lot for the explanation. very helpful!

afadda1 avatar Jan 12 '22 11:01 afadda1