#예스24 #북클럽 에서 지난 번에 읽은 모든 것을 기억하는 남자 이후에 집어든 같은 작가의 같은 시리즈입니다. 영문 제목은 The Fix.
한국어판 표지
영문판 원서 표지
모든것을 기억하는 남자와 같은 주인공인 #에이머스데커 #알렉스재미슨 이 나와요. 같이 FBI에서 미제사건 해결 특수반을 맡아서 일하게 됩니다. 데커는 출근하던 길에 FBI 빌딩 앞에서 어떤 남자가 지나가는 여자를 총으로 쏘고 자신도 이어서 자살해버리는 사건을 목격하게 됩니다. 단순한 살인 사건으로 보였지만, 범인과 희생자의 정체가 밝혀지면서 사건의 진실은 점점 미궁 속으로 빠져들어갑니다.
사실 시리즈 첫번째 책인 모든 것을 기억하는 남자와 이 책 사이에 한 권이 더 있어요 (이 책은 3편이라는 이야기). 그걸 모르고 1편->3편으로 바로 갔는데 그렇게 따로 읽어도 전혀 문제없긴 합니다. 하지만 이제 2편으로 돌아가볼 생각입니다. 북클럽에 있는 지 몰랐어요.
전작과 마찬가지로 흥미로운 사건 전개와 반전, 그리고 쉴 새 없는 액션에 시간가는 줄 모르고 술술 읽힙니다. 번역이 좀 되게 번역체로 되어 있긴 한데, 오히려 단순하게 번역한 게 더 이해가 직관적으로 되는 장점이 있는 것 같습니다.
사실, 반전과 트릭을 보고 아, 이거 어디서 봤는데? 내가 이 책을 읽었던가??? 하고 엄청나게 기시감에 사로잡혀 있었어요. 원문으로 읽었던가? 아니면 영화화나 드라마화 된 걸 본건가? 라고 계속 생각했었는데 결국 정답은 찾지 못했습니다 ㅠㅠ
에이머스데커라는 본 건 절대 까먹지 않는 거구의 형사를 주인공으로 내세워 흥미진진하게 이야기를 풀어나가긴 하지만, 주인공의 또다른 특이점인 공감각이 어떨 때는 일상생활이 불가능할 정도로 영향을 끼친다고 묘사되어 있지만, 또 대부분의 경우에는 언급조차 없어서 조금 아쉬운 부분도 있었어요.
하지만 미스터리 스릴러, 형사물 등등 좋아하시는 분들에게는 추천! 하는 시리즈입니다!!
아 좀 근데 읽고나니 마음이 어두워져서 조금 밝은 책 읽고 시리즈를 계속 봐야할 것 같아요. 달러구트의 꿈백화점 2편을 이럴려고 아껴두었죠 ㅎ
#scRNASeq 분석을 진행하던 도중에 UMAP을 보니 데이터에 #Doublet 이 제법 많이 있다는 사실을 알게 되었습니다. #UMAP 에서 보면 서로 다른 클러스터 사이에 브릿지같은 연결 다리가 보이죠? 그게 서로 다른 두 종류의 세포가 같은 방울 안에 묶여서 같이 시퀀싱 된 더블렛이라고 볼 수 있습니다.
이렇게 데이터 상에서 더블렛이 나타나는 경우에 다음 단계의 분석에 영향을 미칠 수도 있고, 클러스터링이 어려울 수가 있습니다. 그래서 #Seurat 에서는 기본적으로 feature number 혹은 gene number 등에 범위를 줘서 어느정도 doublet을 거르고 있습니다.
Seurat 기본 튜토리얼에서 너무 많거나 적은 RNA가 한 세포 안에 들어가 있는 걸 걸러주는 명령어, 출처: https://satijalab.org/seurat/articles/pbmc3k_tutorial.html
하지만 데이터의 성격이나 품질에 따라서 이러한 스크리닝으로도 더블렛이 제거가 안되는 경우가 있습니다. 그럴 때에 따로 Doublet 을 찾아내는 프로그램을 사용하게 됩니다. 이러한 프로그램에는 여러가지가 있지만 저는 그 중에 #DoubletFinder 라는 프로그램을 선택해서 사용하였습니다. 이 소프트웨어를 쓴 이유는:
사용 설명이 비교적 자세하고 친절하게 되어있음 (그렇지 않은 프로그램들이 참 많죠 ㅠㅠ)
비교적 최근까지 업데이트가 이루어짐
Seurat 파이프라인과 연계가 간편함 (대부분의 다른 프로그램들은 SingleCellExperiment 스트럭쳐로 옮긴 다음에 실행을 해야해서, 데이터 구조를 변화를 시켜서 실행시킨 다음, 결과가 나온 후에 다시 Seurat 파이프라인으로 옮기는 과정이 들어가야 합니다.)
우리 학교 사람이 만들었어요.
이 소프트웨어는 아래의 링크에서 무료로 구하실 수 있습니다. R 패키지로 설치도 가능하구요.
DoubletFinder에 관한 배경이나 설명, 알고리즘은 위의 웹페이지에 자세하게 설명이 되어 있으니 넘어가고 실질적으로 사용하는 과정만 기록하겠습니다.
1. 일단 Seurat package를 불러들입니다. 만약 샘플끼리 이미 Integration을 했다면, 더블렛을 제거하고 다시 Integrate를 하는 것을 권장합니다. 아래의 R 코드는 Cellranger에서 count를 마치고난 결과를 읽어들여서 Seurat object로 만드는 것부터 시작합니다.
#다음 두줄은 나중에 합쳐서 분석을 위해 정보를 추가했습니다 myData$patient<-"p10" myData$visit<-"v3"
#익숙한 방법으로 mt 유전자들을 마크하고 걸러냅니다. 이때, Seurat 기본 튜토리얼에서는 5% 이상은 걸러내라고 적혀있는데, 실제 나중에 나온 실험 논문을 찾아보면 이 수치가 종마다, 조직마다 다르니까 찾아보시는 게 좋습니다 myData[["percent.mt"]]<-PercentageFeatureSet(myData, pattern = "^MT-") VlnPlot(myData,features=c("nFeature_RNA","nCount_RNA","percent.mt"),ncol=3) myData<-subset(myData,subset=nFeature_RNA > 200 & nFeature_RNA<4500 & percent.mt<15)
#이 시점에서 그래프가 나옵니다. pK를 계산해주는 그래프인데, 그래프만으로는 최대값 지점을 한눈에 알기가 어려우므로 다음 줄을 실행시킵니다 bcmvn_myData
pK 그래프, y값이 최고인 지점의 x값의 pK를 찾으면 됩니다.
pK 그래프 표. 맨 오른쪽 BCmetric이 가장 큰 줄의 pK를 찾아서 기억하세요.
4. 이전 단계에서 pK 최대값 지점을 확인한 후에 다음 코드를 실행시킵니다. 여기서 붉은색 숫자는 샘플 안에 Doublet이 얼마나 들어있는지를 예측하는 퍼센트입니다. 이 수치는 시퀀싱 테크닉마다 다르고, 샘플 내의 세포 갯수에 따라 또 다릅니다. 저의 경우에는 10X를 이용했는데 10X Genomics 홈페이지에 다음과 같이 퍼센트가 나와있습니다. 자신 샘플의 세포 갯수에 따라 조절해주세요.
#위의 코드를 실행하고 나면, myData의 구조를 다음 코드로 확인해봅니다 head(myData[[]])
#그 후에 아래의 pANN_0.25_0.005_763 을 자신의 데이터에 맞는 이름으로 고쳐주면 됩니다 myData <- doubletFinder_v3(myData, PCs = 1:30, pN = 0.25, pK =0.005, nExp = nExp_poi.adj, reuse.pANN = "pANN_0.25_0.005_763", sct = FALSE) head(myData[[]])
#그리고, DimPlot을 만들어 Doublet과 Singlet을 관찰해 봅시다 물론 DF.classifications 뒤에 오는 숫자는 알맞게 고쳐줘야합니다DimPlot(myData, group.by = "DF.classifications_0.25_0.005_763", raster = FALSE)
6. 이제 실제로 Doublet이 몇개인지 Singlet이 몇개인지 전체는 몇개인지 알아봅시다. 이런 수치들은 그때그때 적어서 기록으로 남겨두는 편입니다
