今天我們繼續(xù)掃盲，學(xué)習(xí)一些基礎(chǔ)的知識和概念甫煞。
Gene enrichment and covariation analysis
其實我們在做TCR分析的時候菇曲，應(yīng)該也是實驗組 + 對照組進行分析，其中做重要的就是我們要尋找實驗組在接受病原刺激后TCR重排選擇基因的偏好性抚吠。Gene usage preferences were quantified by calculating a normalized Jensen–Shannon divergence (JSD) between the observed gene segment frequencies for each repertoire and background gene frequencies calculated from large-scale repertoire profiling studies常潮，這里其實就是相對于正常的樣本，疾病樣本在TCR重排基因選擇的偏好性埃跷，當然蕊玷，這里用到的是JS散度邮利，大家可以參考文章KL散度、JS散度垃帅、Wasserstein距離延届，JSD 是 Kullback-Leibler 散度的對稱版本，further normalize the JSD values by dividing them by the mean Shannon entropy（香農(nóng)熵贸诚，又叫信息熵方庭，大家參考我之前的文章10X單細胞（10X空間轉(zhuǎn)錄組）基礎(chǔ)算法之KL散度） of the two distributions being compared, which helps to correct for variation in total gene number across segments。To set lower significance thresholds for the JSD heat maps(that is, the values below which the mapped colour is a uniform dark blue)酱固。
圖片.png
we compared the 2–4 different background repertoire datasets(這里就設(shè)置成我們的對照樣本) for each chain/organism to one another and took the largest observed JSD value across all comparisons.
Covariation（協(xié)變械念，協(xié)方差） between gene usage in different segments was quantified using the adjusted mutual information,a variant of the mutual information metric that corrects for the numbers and frequencies of the observed genes (mutual information between pairs of distributions tends to increase with the number of observation classes)。當然运悲，這個在單細胞數(shù)據(jù)中其實應(yīng)該用到的不多龄减。
CDR3 motif discovery.
used a simple, depth-first search procedure to identify over-represented sequence patterns in the CDR3 amino sequences of each repertoire.Motifs were represented as fixed-length patterns consisting of fully-specified amino acid positions, wild card positions, and amino acid group positions,The score of a motif was calculated using a chi-squared formalism:
where ‘observed’ represents the number of times the motif was observed in the repertoire sequences and ‘expected’ represents an estimate of the expected number of observations based on a background set of TCR sequences with V and J gene compositions that match the observed repertoire（這里的背景我們設(shè)置為單細胞的對照樣本）。（這一部分才是最為關(guān)鍵的地方）班眯。
Starting with two-position motifs scoring above a seed threshold, each motif was iteratively extended by adding new specified positions (that is, replacing an internal wild card or lengthening the motif at either end) that increased the motif score.The set of identified motifs were sorted by motif score and filtered for redundancy希停。Finally, motifs scoring above a threshold were extended to include near-neighbour TCRs using a stringent distance threshold; this allowed us to capture additional pattern instances that were not captured by our limited set of amino acid groupings. The final set of motifs for each repertoire were visualized using the TCR logo representation。（看來這才是TCR分析正確的打開方式）署隘。
TCRdiv 多樣性的衡量(也很重要)
為了衡量多樣性宠能，generalizes Simpson’s diversity index by accounting for TCR similarity as well as exact identity（關(guān)于辛普森多樣性指數(shù)，大家可以百度百科一下）磁餐。辛普森多樣性可以被認為是衡量從混合總體中抽取兩個獨立樣本中相同物種或類別的項目的概率违崇，或者換句話說，如果樣本是返回 1 的兩個抽取樣本的函數(shù)的期望值 相同诊霹，否則為 0 羞延。We instead estimate the expected value of a Gaussian function（高斯函數(shù)，確實需要很多的數(shù)學(xué)知識） of the inter-sample distance that returns 1 if the two samples are identical and exp(? (TCRdist(a,b) / s.d.)2) otherwise, where the s.d. was taken to be 18.45 for single-chain distances and twice that for paired analyses based on empirical assessments of receptor distance distributions for multiple epitopes畅哑。Taking the inverse of this estimate gives a diversity measure (TCRdiv) that can be interpreted as an effective population size for similarity-weighted sharing.(這部分有點難以理解肴楷，大家需要多一些耐心和學(xué)習(xí)了).
這部分的代碼在tcr-dist，作者已經(jīng)都封裝好了荠呐，我們用一下就可以赛蔫，感興趣大家可以多多學(xué)習(xí)一下。
到目前為止泥张，算是把基礎(chǔ)說完了呵恢，接下來的分析，就要更上一層樓了媚创。
生活很好渗钉，有你更好
