Cell Line L1 Analysis
Nasir Mahmood Abbasi
2024-07-08
1. load libraries
2. Load Seurat Object
#Load Seurat Object L1
load("../Documents/1-SS-STeps/4-Analysis_and_Robj_Marie/analyse juillet 2023/ObjetsR/L1.Robj")
L13. QC
# Set identity classes to an existing column in meta data
Idents(object = L1) <- "cell_line"
L1[["percent.rb"]] <- PercentageFeatureSet(L1, pattern = "^RP[SL]")
VlnPlot(L1, features = c("nFeature_RNA",
"nCount_RNA",
"percent.mito",
"percent.rb"),
ncol = 4, pt.size = 0.1) &
theme(plot.title = element_text(size=10))
FeatureScatter(L1, feature1 = "percent.mito",
feature2 = "percent.rb")
VlnPlot(L1, features = c("nFeature_RNA",
"nCount_RNA",
"percent.mito"),
ncol = 3)
FeatureScatter(L1,
feature1 = "percent.mito",
feature2 = "percent.rb") +
geom_smooth(method = 'lm')
FeatureScatter(L1,
feature1 = "nCount_RNA",
feature2 = "nFeature_RNA") +
geom_smooth(method = 'lm')
##FeatureScatter is typically used to visualize feature-feature relationships ##for anything calculated by the object, ##i.e. columns in object metadata, PC scores etc.
FeatureScatter(L1,
feature1 = "nCount_RNA",
feature2 = "percent.mito")+
geom_smooth(method = 'lm')
FeatureScatter(L1,
feature1 = "nCount_RNA",
feature2 = "nFeature_RNA")+
geom_smooth(method = 'lm')
Assign Cell-Cycle Scores
Running SCTransform on assay: RNA
vst.flavor='v2' set. Using model with fixed slope and excluding poisson genes.
Calculating cell attributes from input UMI matrix: log_umi
Variance stabilizing transformation of count matrix of size 17646 by 5825
Model formula is y ~ log_umi
Get Negative Binomial regression parameters per gene
Using 2000 genes, 5000 cells
Found 99 outliers - those will be ignored in fitting/regularization step
Second step: Get residuals using fitted parameters for 17646 genes
Computing corrected count matrix for 17646 genes
Calculating gene attributes
Wall clock passed: Time difference of 19.01677 secs
Determine variable features
|
| | 0%
|
|============================== | 25%
|
|============================================================= | 50%
|
|============================================================================================ | 75%
|
|==========================================================================================================================| 100%
Place corrected count matrix in counts slot
Set default assay to SCT
Warning: The following features are not present in the object: MLF1IP, not searching for symbol synonymsWarning: The following features are not present in the object: FAM64A, HN1, not searching for symbol synonyms4. Normalize data
# Apply SCTransform
L1 <- SCTransform(L1, vars.to.regress = c("percent.rb","percent.mito", "CC.Difference"),
do.scale=TRUE,
do.center=TRUE,
verbose = TRUE)Running SCTransform on assay: RNA
vst.flavor='v2' set. Using model with fixed slope and excluding poisson genes.
Calculating cell attributes from input UMI matrix: log_umi
Variance stabilizing transformation of count matrix of size 17646 by 5825
Model formula is y ~ log_umi
Get Negative Binomial regression parameters per gene
Using 2000 genes, 5000 cells
Found 99 outliers - those will be ignored in fitting/regularization step
Second step: Get residuals using fitted parameters for 17646 genes
Computing corrected count matrix for 17646 genes
Calculating gene attributes
Wall clock passed: Time difference of 18.97784 secs
Determine variable features
Regressing out percent.rb, percent.mito, CC.Difference
|
| | 0%
|
|= | 0%
|
|= | 1%
|
|== | 1%
|
|== | 2%
|
|=== | 2%
|
|=== | 3%
|
|==== | 3%
|
|==== | 4%
|
|===== | 4%
|
|====== | 5%
|
|======= | 5%
|
|======= | 6%
|
|======== | 6%
|
|======== | 7%
|
|========= | 7%
|
|========= | 8%
|
|========== | 8%
|
|========== | 9%
|
|=========== | 9%
|
|============ | 9%
|
|============ | 10%
|
|============= | 10%
|
|============= | 11%
|
|============== | 11%
|
|============== | 12%
|
|=============== | 12%
|
|=============== | 13%
|
|================ | 13%
|
|================ | 14%
|
|================= | 14%
|
|================== | 14%
|
|================== | 15%
|
|=================== | 15%
|
|=================== | 16%
|
|==================== | 16%
|
|==================== | 17%
|
|===================== | 17%
|
|===================== | 18%
|
|====================== | 18%
|
|======================= | 18%
|
|======================= | 19%
|
|======================== | 19%
|
|======================== | 20%
|
|========================= | 20%
|
|========================= | 21%
|
|========================== | 21%
|
|========================== | 22%
|
|=========================== | 22%
|
|=========================== | 23%
|
|============================ | 23%
|
|============================= | 23%
|
|============================= | 24%
|
|============================== | 24%
|
|============================== | 25%
|
|=============================== | 25%
|
|=============================== | 26%
|
|================================ | 26%
|
|================================ | 27%
|
|================================= | 27%
|
|================================== | 27%
|
|================================== | 28%
|
|=================================== | 28%
|
|=================================== | 29%
|
|==================================== | 29%
|
|==================================== | 30%
|
|===================================== | 30%
|
|===================================== | 31%
|
|====================================== | 31%
|
|====================================== | 32%
|
|======================================= | 32%
|
|======================================== | 32%
|
|======================================== | 33%
|
|========================================= | 33%
|
|========================================= | 34%
|
|========================================== | 34%
|
|========================================== | 35%
|
|=========================================== | 35%
|
|=========================================== | 36%
|
|============================================ | 36%
|
|============================================= | 36%
|
|============================================= | 37%
|
|============================================== | 37%
|
|============================================== | 38%
|
|=============================================== | 38%
|
|=============================================== | 39%
|
|================================================ | 39%
|
|================================================ | 40%
|
|================================================= | 40%
|
|================================================= | 41%
|
|================================================== | 41%
|
|=================================================== | 41%
|
|=================================================== | 42%
|
|==================================================== | 42%
|
|==================================================== | 43%
|
|===================================================== | 43%
|
|===================================================== | 44%
|
|====================================================== | 44%
|
|====================================================== | 45%
|
|======================================================= | 45%
|
|======================================================== | 46%
|
|========================================================= | 46%
|
|========================================================= | 47%
|
|========================================================== | 47%
|
|========================================================== | 48%
|
|=========================================================== | 48%
|
|=========================================================== | 49%
|
|============================================================ | 49%
|
|============================================================ | 50%
|
|============================================================= | 50%
|
|============================================================== | 50%
|
|============================================================== | 51%
|
|=============================================================== | 51%
|
|=============================================================== | 52%
|
|================================================================ | 52%
|
|================================================================ | 53%
|
|================================================================= | 53%
|
|================================================================= | 54%
|
|================================================================== | 54%
|
|================================================================== | 55%
|
|=================================================================== | 55%
|
|==================================================================== | 55%
|
|==================================================================== | 56%
|
|===================================================================== | 56%
|
|===================================================================== | 57%
|
|====================================================================== | 57%
|
|====================================================================== | 58%
|
|======================================================================= | 58%
|
|======================================================================= | 59%
|
|======================================================================== | 59%
|
|========================================================================= | 59%
|
|========================================================================= | 60%
|
|========================================================================== | 60%
|
|========================================================================== | 61%
|
|=========================================================================== | 61%
|
|=========================================================================== | 62%
|
|============================================================================ | 62%
|
|============================================================================ | 63%
|
|============================================================================= | 63%
|
|============================================================================= | 64%
|
|============================================================================== | 64%
|
|=============================================================================== | 64%
|
|=============================================================================== | 65%
|
|================================================================================ | 65%
|
|================================================================================ | 66%
|
|================================================================================= | 66%
|
|================================================================================= | 67%
|
|================================================================================== | 67%
|
|================================================================================== | 68%
|
|=================================================================================== | 68%
|
|==================================================================================== | 68%
|
|==================================================================================== | 69%
|
|===================================================================================== | 69%
|
|===================================================================================== | 70%
|
|====================================================================================== | 70%
|
|====================================================================================== | 71%
|
|======================================================================================= | 71%
|
|======================================================================================= | 72%
|
|======================================================================================== | 72%
|
|======================================================================================== | 73%
|
|========================================================================================= | 73%
|
|========================================================================================== | 73%
|
|========================================================================================== | 74%
|
|=========================================================================================== | 74%
|
|=========================================================================================== | 75%
|
|============================================================================================ | 75%
|
|============================================================================================ | 76%
|
|============================================================================================= | 76%
|
|============================================================================================= | 77%
|
|============================================================================================== | 77%
|
|=============================================================================================== | 77%
|
|=============================================================================================== | 78%
|
|================================================================================================ | 78%
|
|================================================================================================ | 79%
|
|================================================================================================= | 79%
|
|================================================================================================= | 80%
|
|================================================================================================== | 80%
|
|================================================================================================== | 81%
|
|=================================================================================================== | 81%
|
|=================================================================================================== | 82%
|
|==================================================================================================== | 82%
|
|===================================================================================================== | 82%
|
|===================================================================================================== | 83%
|
|====================================================================================================== | 83%
|
|====================================================================================================== | 84%
|
|======================================================================================================= | 84%
|
|======================================================================================================= | 85%
|
|======================================================================================================== | 85%
|
|======================================================================================================== | 86%
|
|========================================================================================================= | 86%
|
|========================================================================================================== | 86%
|
|========================================================================================================== | 87%
|
|=========================================================================================================== | 87%
|
|=========================================================================================================== | 88%
|
|============================================================================================================ | 88%
|
|============================================================================================================ | 89%
|
|============================================================================================================= | 89%
|
|============================================================================================================= | 90%
|
|============================================================================================================== | 90%
|
|============================================================================================================== | 91%
|
|=============================================================================================================== | 91%
|
|================================================================================================================ | 91%
|
|================================================================================================================ | 92%
|
|================================================================================================================= | 92%
|
|================================================================================================================= | 93%
|
|================================================================================================================== | 93%
|
|================================================================================================================== | 94%
|
|=================================================================================================================== | 94%
|
|=================================================================================================================== | 95%
|
|==================================================================================================================== | 95%
|
|===================================================================================================================== | 96%
|
|====================================================================================================================== | 96%
|
|====================================================================================================================== | 97%
|
|======================================================================================================================= | 97%
|
|======================================================================================================================= | 98%
|
|======================================================================================================================== | 98%
|
|======================================================================================================================== | 99%
|
|========================================================================================================================= | 99%
|
|========================================================================================================================= | 100%
|
|==========================================================================================================================| 100%
Centering and scaling data matrix
|
| | 0%
|
|============================== | 25%
|
|============================================================= | 50%
|
|============================================================================================ | 75%
|
|==========================================================================================================================| 100%
Place corrected count matrix in counts slot
Set default assay to SCT 5. Perform PCA
Variables_genes <- L1@assays$SCT@var.features
# Exclude genes starting with "HLA-" AND "Xist" AND "TRBV, TRAV"
Variables_genes_after_exclusion <- Variables_genes[!grepl("^HLA-|^XIST|^TRBV|^TRAV", Variables_genes)]
# These are now standard steps in the Seurat workflow for visualization and clustering
L1 <- RunPCA(L1,
features = Variables_genes_after_exclusion,
do.print = TRUE,
pcs.print = 1:5,
genes.print = 15,
npcs = 50)PC_ 1
Positive: TUBA1B, H2AFZ, HMGB1, STMN1, DUT, TUBB, HMGB2, MCM7, HMGN2, HIST1H4C
UBE2C, TYMS, PCLAF, PTMA, TUBB4B, TK1, DEK, RRM2, MCM3, CDT1
CKS2, PKMYT1, SIVA1, TOP2A, DNMT1, CENPM, UBE2S, ALYREF, RANBP1, HIST1H1A
Negative: CCND2, LTB, RNF213, BTG2, CD44, SLFN5, CD48, NFATC3, CTSC, DDIT4
GAS5, CSF1, INPP4B, SEMA4A, CD74, CALD1, IKZF3, TGFBR3, BTG1, ZYX
KDM5B, EGLN3, DLEU1, PNRC1, PTGIS, PTPRC, VIM, KDM7A, IL4, RALGAPA2
PC_ 2
Positive: ACTB, TMSB10, PFN1, NME1, MIF, SRRT, CLIC1, ATP5MC3, EIF4A1, PRDX1
MRPL4, CFL1, PRELID1, HSP90AB1, PPIB, ATP5MF, POMP, TIMM13, LDHA, SSBP1
NDUFS5, MRPL16, PARK7, BCAP31, PSMB3, CORO1A, PSMA7, PSME2, ACTG1, GADD45GIP1
Negative: MALAT1, ASPM, CENPF, MKI67, TOP2A, HIST1H1E, HIST1H1B, TUBB, PCLAF, NUSAP1
CAMK4, RRM2, TYMS, ATAD2, MBNL1, DHFR, SYNE2, RCSD1, HIST1H1C, STMN1
RABGAP1L, PDE3B, KIF11, IKZF2, PTPRC, MAL, NEIL3, NEK7, NCAPG2, BRCA1
PC_ 3
Positive: RPS28, RPL32, RPL29, RPL35A, RPS15, RPL36, RPL14, GZMA, CEBPD, HSP90AB1
LSR, SLCO3A1, ADGRE5, OAZ1, FXYD5, EEF2, AC002069.2, BSG, RPSA, MLLT3
FAM107B, AC069410.1, PDGFD, LEF1, P2RY14, FOXP1, SPINT2, H1FX, PRDX2, TRBC2
Negative: GAPDH, HPGDS, PKM, NKG7, KRT1, FABP5, LY6E, S100A4, NPDC1, EEF1A1
VIM, CD48, NQO1, S100P, C12orf75, CFH, FSCN1, SLC25A5, LMNA, SH3BGRL3
SH2D2A, CYP1B1, PRDX1, SIX3, ALOX5AP, ACTG1, P4HA2, TMSB4X, BACE2, TPST2
PC_ 4
Positive: RPL10, MRPL16, ID3, SEPTIN11, HTATSF1, LAGE3, TCN1, GIMAP7, BCAP31, RIPOR2
PON2, ARHGEF6, FYB1, NRXN3, KLF2, LINC00892, TMSB10, AIFM1, STK26, DKC1
TCF7, DSC1, ZP3, CD81, MKRN1, TAFA2, IDH3G, COA1, COBL, IKZF2
Negative: RPL14, RPL32, RPL29, IL32, RPS27, CYBA, CTSC, CISH, RPSA, CDKN1A
RHOC, RPL35A, CACYBP, ACTB, GAS5, IFITM1, ATP1B1, SOCS1, FTH1, FXYD5
CXCR3, CEBPD, LAT, S100A11, H3F3A, LTB, GZMM, ITGB7, PIM1, B2M
PC_ 5
Positive: KLF2, HNRNPA1, S100A10, IMPDH2, GIMAP7, NUCKS1, SH3BP5, MAL, MT-CYB, TLE5
VIM, ETS1, S1PR4, SELPLG, KLF3, CD55, IL7R, PTMA, RPSA, HDGF
ITGB1, CCR7, ATP5MC2, NOSIP, HSD17B11, CDKN2B, CACYBP, FGFBP2, HMGB2, RPL28
Negative: HIST1H1A, SKAP1, HIST1H1B, TUBA4A, HIST1H1E, HIST1H2AH, HIST1H4C, HIST1H1D, HIST1H1C, SRRT
CTSC, HIST1H3B, MRPL16, TUBB, SSR4, EIF4G3, ARHGAP4, RPL10, HIST1H3D, CD96
SRM, MRPL41, HIST2H2AC, PTPN7, HIST1H3C, HIST1H4F, HIST1H2AL, SEPTIN11, TMC8, RALGAPA2 # determine dimensionality of the data
ElbowPlot(L1, ndims =50)
NA
NAPerform PCA TEST
library(ggplot2)
library(RColorBrewer)
# Assuming you have 10 different cell lines, generating a color palette with 10 colors
cell_line_colors <- brewer.pal(10, "Set3")
# Assuming L1$cell_line is a factor or character vector containing cell line names
data <- as.data.frame(table(L1$cell_line))
colnames(data) <- c("cell_line", "nUMI") # Change column name to nUMI
ncells <- ggplot(data, aes(x = cell_line, y = nUMI, fill = cell_line)) +
geom_col() +
theme_classic() +
geom_text(aes(label = nUMI),
position = position_dodge(width = 0.9),
vjust = -0.25) +
scale_fill_manual(values = cell_line_colors) +
theme(axis.text.x = element_text(angle = 45, hjust = 1),
plot.title = element_text(hjust = 0.5)) + # Adjust the title position
ggtitle("Filtered cells per sample") +
xlab("Cell lines") + # Adjust x-axis label
ylab("Frequency") # Adjust y-axis label
print(ncells)
# TEST-1
# given that the output of RunPCA is "pca"
# replace "so" by the name of your seurat object
pct <- L1[["pca"]]@stdev / sum(L1[["pca"]]@stdev) * 100
cumu <- cumsum(pct) # Calculate cumulative percents for each PC
# Determine the difference between variation of PC and subsequent PC
co2 <- sort(which((pct[-length(pct)] - pct[-1]) > 0.1), decreasing = T)[1] + 1
# last point where change of % of variation is more than 0.1%. -> co2
co2[1] 15# TEST-2
# get significant PCs
stdv <- L1[["pca"]]@stdev
sum.stdv <- sum(L1[["pca"]]@stdev)
percent.stdv <- (stdv / sum.stdv) * 100
cumulative <- cumsum(percent.stdv)
co1 <- which(cumulative > 90 & percent.stdv < 5)[1]
co2 <- sort(which((percent.stdv[1:length(percent.stdv) - 1] -
percent.stdv[2:length(percent.stdv)]) > 0.1),
decreasing = T)[1] + 1
min.pc <- min(co1, co2)
min.pc[1] 15# Create a dataframe with values
plot_df <- data.frame(pct = percent.stdv,
cumu = cumulative,
rank = 1:length(percent.stdv))
# Elbow plot to visualize
ggplot(plot_df, aes(cumulative, percent.stdv, label = rank, color = rank > min.pc)) +
geom_text() +
geom_vline(xintercept = 90, color = "grey") +
geom_hline(yintercept = min(percent.stdv[percent.stdv > 5]), color = "grey") +
theme_bw()
NA
NA
NA6. Clustering
L1 <- FindNeighbors(L1,
dims = 1:min.pc,
verbose = FALSE)
# understanding resolution
L1 <- FindClusters(L1,
resolution = c(0.1, 0.2, 0.3, 0.4, 0.5, 0.6,
0.7,0.8, 0.9, 1, 1.1, 1.2))Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 5825
Number of edges: 200767
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.9215
Number of communities: 3
Elapsed time: 0 seconds
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 5825
Number of edges: 200767
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.8916
Number of communities: 6
Elapsed time: 0 seconds
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 5825
Number of edges: 200767
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.8701
Number of communities: 7
Elapsed time: 0 seconds
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 5825
Number of edges: 200767
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.8494
Number of communities: 7
Elapsed time: 0 seconds
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 5825
Number of edges: 200767
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.8287
Number of communities: 7
Elapsed time: 0 seconds
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 5825
Number of edges: 200767
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.8146
Number of communities: 9
Elapsed time: 0 seconds
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 5825
Number of edges: 200767
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.8022
Number of communities: 10
Elapsed time: 0 seconds
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 5825
Number of edges: 200767
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.7907
Number of communities: 10
Elapsed time: 0 seconds
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 5825
Number of edges: 200767
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.7793
Number of communities: 11
Elapsed time: 0 seconds
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 5825
Number of edges: 200767
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.7681
Number of communities: 12
Elapsed time: 0 seconds
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 5825
Number of edges: 200767
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.7579
Number of communities: 13
Elapsed time: 0 seconds
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 5825
Number of edges: 200767
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.7486
Number of communities: 14
Elapsed time: 0 seconds
# non-linear dimensionality reduction --------------
L1 <- RunUMAP(L1,
dims = 1:min.pc,
verbose = FALSE)Warning: The default method for RunUMAP has changed from calling Python UMAP via reticulate to the R-native UWOT using the cosine metric
To use Python UMAP via reticulate, set umap.method to 'umap-learn' and metric to 'correlation'
This message will be shown once per session
# note that you can set `label = TRUE` or use the Label Clusters function to help label
# individual clusters
DimPlot(L1,group.by = "cell_line",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1,
group.by = "SCT_snn_res.0.1",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1,
group.by = "SCT_snn_res.0.2",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1,
group.by = "SCT_snn_res.0.3",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1,
group.by = "SCT_snn_res.0.4",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1,
group.by = "SCT_snn_res.0.5",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1,
group.by = "SCT_snn_res.0.6",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1,
group.by = "SCT_snn_res.0.7",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1,
group.by = "SCT_snn_res.0.8",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1,
group.by = "SCT_snn_res.0.9",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1,
group.by = "SCT_snn_res.1",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1,
group.by = "SCT_snn_res.1.1",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1,
group.by = "SCT_snn_res.1.2",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
7. Azimuth Annotation
InstallData("pbmcref")Warning: The following packages are already installed and will not be reinstalled: pbmcref# The RunAzimuth function can take a Seurat object as input
L1 <- RunAzimuth(L1, reference = "pbmcref")Warning: Overwriting miscellanous data for modelWarning: Adding a dimensional reduction (refUMAP) without the associated assay being presentWarning: Adding a dimensional reduction (refUMAP) without the associated assay being presentdetected inputs from HUMAN with id type Gene.name
reference rownames detected HUMAN with id type Gene.name
Normalizing query using reference SCT model
Warning: 113 features of the features specified were not present in both the reference query assays.
Continuing with remaining 4887 features.Projecting cell embeddings
Finding query neighbors
Finding neighborhoods
Finding anchors
Found 489 anchors
Finding integration vectors
Finding integration vector weights
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Predicting cell labels
Predicting cell labels
Warning: Feature names cannot have underscores ('_'), replacing with dashes ('-')Predicting cell labels
Warning: Feature names cannot have underscores ('_'), replacing with dashes ('-')
| | 0 % ~calculating
Integrating dataset 2 with reference dataset
Finding integration vectors
Integrating data
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=01s Warning: Keys should be one or more alphanumeric characters followed by an underscore, setting key from integrated_dr_ to integrateddr_Computing nearest neighbors
Running UMAP projection
Warning: Number of neighbors between query and reference is not equal to the number of neighbors within reference12:44:57 Read 5825 rows
12:44:57 Processing block 1 of 1
12:44:57 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 20
12:44:57 Initializing by weighted average of neighbor coordinates using 1 thread
12:44:57 Commencing optimization for 67 epochs, with 116500 positive edges
Using method 'umap'
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
12:44:58 Finished
Warning: No assay specified, setting assay as RNA by default.Projecting reference PCA onto query
Finding integration vector weights
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Projecting back the query cells into original PCA space
Finding integration vector weights
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Computing scores:
Finding neighbors of original query cells
Finding neighbors of transformed query cells
Computing query SNN
Determining bandwidth and computing transition probabilities
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Total elapsed time: 2.2985680103302UMAPPlot(L1, group.by = "predicted.celltype.l2",
reduction = "umap",
label.size = 3,
repel = T,
label = F)
UMAPPlot(L1, group.by = "predicted.celltype.l2",
reduction = "umap",
label.size = 3,
repel = T,
label = T)
8. Cell type annotation using ProjectTils
#Load reference atlas and query data
ref <- readRDS(file = "CD4T_human_ref_v1.rds")
#Run Projection algorithm
query.projected <- Run.ProjecTILs(L1, ref = ref)
|
| | 0%[1] "Using assay SCT for query"Pre-filtering cells with scGate...
### Detected a total of 5396 pure 'Target' cells (92.64% of total)[1] "429 out of 5825 ( 7% ) non-pure cells removed. Use filter.cells=FALSE to avoid pre-filtering"Performing log-normalization
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|[1] "Aligning query to reference map for batch-correction..."Warning: Layer counts isn't present in the assay object[[assay]]; returning NULLPreparing PCA embeddings for objects...
Warning: Number of dimensions changing from 50 to 20
| | 0 % ~calculating
|+++++++++++++++++++++++++ | 50% ~05s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=06s
| | 0 % ~calculating
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=13s
Projecting corrected query onto Reference PCA space
Projecting corrected query onto Reference UMAP spaceWarning: Not all features provided are in this Assay object, removing the following feature(s): CXCL13, GNLY, TNFRSF4, CD177, CCL22, LAG3, CGA, G0S2, IL17A, CCL4L2, TRDC, IL21, TNFRSF18, FOXP3, HLA-DQA1, KRT86, CAV1, IL10, CCL3, IL1R2, ACTG2, CCR8, EBI3, DIRAS3, KLRD1, THBS1, ZBED2, CAVIN3, GNG8, MS4A6A, TYROBP, PROK2, BEX3, BASP1, GLUL, PVALB, HES4, MRC1, CTLA4, TASL, IL1RN, LHFPL6, GBP5, MYADM, H1-4, KLF4, CTSH, H2AZ1, TENT5A, CD9, MSC, RCAN2, GPR25, CCR6, TNS3, ELAPOR1, NR4A3, IGFBP3, CCL3L3, PRRG4, POLR1F, SDC4, FCRL3, GPR55, PLK2, CRLF2, AHSP, FAIM2, CSF2RB, CDCA7L, HTRA1, H1-2, H1-0, ETV7, FLT1, LMCD1, WARS1, H1-3, ASCL2, CLNK, H2BC11, ADRB2, GPX1, STAC, H2AC6, IRAG2, GNA15, CD80, PLAAT3, H1-10, MYO7A, DAPK2, H3C10, NAP1L2, IL22, CHGB, KRT81, SLC28A3, PDLIM4, ZNF683, ECEL1, HES1, HSD11B1, CPE, CD40, CA2, NT5E, CD86, H4C3, RORC, CXXC5, NEBL, GSTM2
|
|==========================================================================================================================| 100%Creating slots functional.cluster and functional.cluster.conf in query object#reference atlas
DimPlot(ref, label = T)
#Visualize projection
plot.projection(ref, query.projected, linesize = 0.5, pointsize = 0.5)
#Plot the predicted composition of the query in terms of reference T cell subtypes
plot.statepred.composition(ref, query.projected, metric = "Percent")
L1 <- ProjecTILs.classifier(query = L1, ref = ref)
|
| | 0%[1] "Using assay SCT for query"Pre-filtering cells with scGate...
### Detected a total of 5396 pure 'Target' cells (92.64% of total)[1] "429 out of 5825 ( 7% ) non-pure cells removed. Use filter.cells=FALSE to avoid pre-filtering"Performing log-normalization
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|[1] "Aligning query to reference map for batch-correction..."Warning: Layer counts isn't present in the assay object[[assay]]; returning NULLPreparing PCA embeddings for objects...
Warning: Number of dimensions changing from 50 to 20
| | 0 % ~calculating
|+++++++++++++++++++++++++ | 50% ~05s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=06s
| | 0 % ~calculating
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=13s
Projecting corrected query onto Reference PCA space
Projecting corrected query onto Reference UMAP spaceWarning: Not all features provided are in this Assay object, removing the following feature(s): CXCL13, GNLY, TNFRSF4, CD177, CCL22, LAG3, CGA, G0S2, IL17A, CCL4L2, TRDC, IL21, TNFRSF18, FOXP3, HLA-DQA1, KRT86, CAV1, IL10, CCL3, IL1R2, ACTG2, CCR8, EBI3, DIRAS3, KLRD1, THBS1, ZBED2, CAVIN3, GNG8, MS4A6A, TYROBP, PROK2, BEX3, BASP1, GLUL, PVALB, HES4, MRC1, CTLA4, TASL, IL1RN, LHFPL6, GBP5, MYADM, H1-4, KLF4, CTSH, H2AZ1, TENT5A, CD9, MSC, RCAN2, GPR25, CCR6, TNS3, ELAPOR1, NR4A3, IGFBP3, CCL3L3, PRRG4, POLR1F, SDC4, FCRL3, GPR55, PLK2, CRLF2, AHSP, FAIM2, CSF2RB, CDCA7L, HTRA1, H1-2, H1-0, ETV7, FLT1, LMCD1, WARS1, H1-3, ASCL2, CLNK, H2BC11, ADRB2, GPX1, STAC, H2AC6, IRAG2, GNA15, CD80, PLAAT3, H1-10, MYO7A, DAPK2, H3C10, NAP1L2, IL22, CHGB, KRT81, SLC28A3, PDLIM4, ZNF683, ECEL1, HES1, HSD11B1, CPE, CD40, CA2, NT5E, CD86, H4C3, RORC, CXXC5, NEBL, GSTM2
|
|===================================================================================================================================| 100%Creating slots functional.cluster and functional.cluster.conf in query objectUMAPPlot(L1, group.by = "functional.cluster",
reduction = "umap",
label.size = 3,
repel = T,
label = T)
NA
NA9.Cell type annotation using SingleR
# Monaco Annotations
DimPlot(L1, group.by = "monaco.main",
reduction = "umap",
label.size = 3,
repel = T,
label = F)
DimPlot(L1, group.by = "monaco.main",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1, group.by = "monaco.fine",
reduction = "umap",
label.size = 3,
repel = T,
label = F)
DimPlot(L1, group.by = "monaco.fine",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
# HPCA Annotations
DimPlot(L1, group.by = "hpca.main",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1, group.by = "hpca.fine",
reduction = "umap",
label.size = 3,
repel = T,
label = F)
DimPlot(L1, group.by = "hpca.fine",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
# DICE Annotations
DimPlot(L1, group.by = "dice.main",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1, group.by = "dice.fine",
reduction = "umap",
label.size = 3,
repel = T,
label = F)
DimPlot(L1, group.by = "dice.fine",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
# BPE Annotations
DimPlot(L1, group.by = "bpe.main",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
DimPlot(L1, group.by = "bpe.fine",
reduction = "umap",
label.size = 3,
repel = T,
label = F)
DimPlot(L1, group.by = "bpe.fine",
reduction = "umap",
label.size = 3,
repel = T,
label = T, label.box = T)
NA
NA10. clusTree
library(clustree)
clustree(L1, prefix = "SCT_snn_res.")
11.Save the Seurat object as an Robj file
LS0tCnRpdGxlOiAiQ2VsbCBMaW5lIEwxIEFuYWx5c2lzIgphdXRob3I6IE5hc2lyIE1haG1vb2QgQWJiYXNpCmRhdGU6ICJgciBTeXMuRGF0ZSgpYCIKb3V0cHV0OgogICNybWRmb3JtYXRzOjpyZWFkdGhlZG93bgogIGh0bWxfbm90ZWJvb2s6CiAgICB0b2M6IHRydWUKICAgIHRvY19mbG9hdDogdHJ1ZQogICAgdG9jX2NvbGxhcHNlZDogdHJ1ZQoKLS0tCgojIDEuIGxvYWQgbGlicmFyaWVzCmBgYHtyIHNldHVwLCBpbmNsdWRlPUZBTFNFfQoKbGlicmFyeShTZXVyYXQpCmxpYnJhcnkoU2V1cmF0T2JqZWN0KQpsaWJyYXJ5KFNldXJhdERhdGEpCmxpYnJhcnkocGF0Y2h3b3JrKQoKbGlicmFyeShkcGx5cikKbGlicmFyeSh0aWR5dmVyc2UpCmxpYnJhcnkoZ2dwbG90MikKbGlicmFyeShSQ29sb3JCcmV3ZXIpCmxpYnJhcnkobWFncml0dHIpCmxpYnJhcnkoZGJwbHlyKQpsaWJyYXJ5KHJtYXJrZG93bikKbGlicmFyeShrbml0cikKbGlicmFyeSh0aW55dGV4KQojQXppbXV0aCBBbm5vdGF0aW9uIGxpYnJhcmllcwpsaWJyYXJ5KEF6aW11dGgpCiNQcm9qZWNUaWxzIEFubm90YXRpb24gbGlicmFyaWVzCmxpYnJhcnkoU1RBQ0FTKQpsaWJyYXJ5KFByb2plY1RJTHMpCiNzaW5nbGVSIEFubm90YXRpb24gbGlicmFyaWVzCmxpYnJhcnkoU2luZ2xlUikKbGlicmFyeShjZWxsZGV4KQpsaWJyYXJ5KFNpbmdsZUNlbGxFeHBlcmltZW50KQoKYGBgCgoKIyAyLiBMb2FkIFNldXJhdCBPYmplY3QgCmBgYHtyIGxvYWRfc2V1cmF0fQoKI0xvYWQgU2V1cmF0IE9iamVjdCBMMQpsb2FkKCIuLi9Eb2N1bWVudHMvMS1TUy1TVGVwcy80LUFuYWx5c2lzX2FuZF9Sb2JqX01hcmllL2FuYWx5c2UganVpbGxldCAyMDIzL09iamV0c1IvTDEuUm9iaiIpCgpMMQpgYGAKCgojIDMuIFFDCmBgYHtyIFFDLCBmaWcuaGVpZ2h0PTYsIGZpZy53aWR0aD0xMH0KIyBTZXQgaWRlbnRpdHkgY2xhc3NlcyB0byBhbiBleGlzdGluZyBjb2x1bW4gaW4gbWV0YSBkYXRhCklkZW50cyhvYmplY3QgPSBMMSkgPC0gImNlbGxfbGluZSIKCkwxW1sicGVyY2VudC5yYiJdXSA8LSBQZXJjZW50YWdlRmVhdHVyZVNldChMMSwgcGF0dGVybiA9ICJeUlBbU0xdIikKClZsblBsb3QoTDEsIGZlYXR1cmVzID0gYygibkZlYXR1cmVfUk5BIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIm5Db3VudF9STkEiLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAicGVyY2VudC5taXRvIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAicGVyY2VudC5yYiIpLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgIG5jb2wgPSA0LCBwdC5zaXplID0gMC4xKSAmIAogICAgICAgICAgICAgIHRoZW1lKHBsb3QudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZT0xMCkpCgpGZWF0dXJlU2NhdHRlcihMMSwgZmVhdHVyZTEgPSAicGVyY2VudC5taXRvIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBmZWF0dXJlMiA9ICJwZXJjZW50LnJiIikKClZsblBsb3QoTDEsIGZlYXR1cmVzID0gYygibkZlYXR1cmVfUk5BIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJuQ291bnRfUk5BIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJwZXJjZW50Lm1pdG8iKSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbmNvbCA9IDMpCgpGZWF0dXJlU2NhdHRlcihMMSwgCiAgICAgICAgICAgICAgIGZlYXR1cmUxID0gInBlcmNlbnQubWl0byIsIAogICAgICAgICAgICAgICBmZWF0dXJlMiA9ICJwZXJjZW50LnJiIikgKwogICAgICAgIGdlb21fc21vb3RoKG1ldGhvZCA9ICdsbScpCgpGZWF0dXJlU2NhdHRlcihMMSwgCiAgICAgICAgICAgICAgIGZlYXR1cmUxID0gIm5Db3VudF9STkEiLCAKICAgICAgICAgICAgICAgZmVhdHVyZTIgPSAibkZlYXR1cmVfUk5BIikgKwogICAgICAgIGdlb21fc21vb3RoKG1ldGhvZCA9ICdsbScpCgpgYGAKIyNGZWF0dXJlU2NhdHRlciBpcyB0eXBpY2FsbHkgdXNlZCB0byB2aXN1YWxpemUgZmVhdHVyZS1mZWF0dXJlIHJlbGF0aW9uc2hpcHMKIyNmb3IgYW55dGhpbmcgY2FsY3VsYXRlZCBieSB0aGUgb2JqZWN0LCAKIyNpLmUuIGNvbHVtbnMgaW4gb2JqZWN0IG1ldGFkYXRhLCBQQyBzY29yZXMgZXRjLgoKYGBge3IgRkMsIGZpZy5oZWlnaHQ9NiwgZmlnLndpZHRoPTEwfQoKRmVhdHVyZVNjYXR0ZXIoTDEsIAogICAgICAgICAgICAgICBmZWF0dXJlMSA9ICJuQ291bnRfUk5BIiwgCiAgICAgICAgICAgICAgIGZlYXR1cmUyID0gInBlcmNlbnQubWl0byIpKwogIGdlb21fc21vb3RoKG1ldGhvZCA9ICdsbScpCgpGZWF0dXJlU2NhdHRlcihMMSwgCiAgICAgICAgICAgICAgIGZlYXR1cmUxID0gIm5Db3VudF9STkEiLCAKICAgICAgICAgICAgICAgZmVhdHVyZTIgPSAibkZlYXR1cmVfUk5BIikrCiAgZ2VvbV9zbW9vdGgobWV0aG9kID0gJ2xtJykKCmBgYAojIyBBc3NpZ24gQ2VsbC1DeWNsZSBTY29yZXMKYGBge3IgUmVncmVzcywgZWNobz1GQUxTRSwgZmlnLmhlaWdodD02LCBmaWcud2lkdGg9MTB9CgpMMSA8LSBTQ1RyYW5zZm9ybShMMSwgZG8uc2NhbGU9RkFMU0UsIGRvLmNlbnRlcj1GQUxTRSkKCgojIEEgbGlzdCBvZiBjZWxsIGN5Y2xlIG1hcmtlcnMsIGZyb20gVGlyb3NoIGV0IGFsLCAyMDE1LCBpcyBsb2FkZWQgd2l0aCBTZXVyYXQuICBXZSBjYW4KIyBzZWdyZWdhdGUgdGhpcyBsaXN0IGludG8gbWFya2VycyBvZiBHMi9NIHBoYXNlIGFuZCBtYXJrZXJzIG9mIFMgcGhhc2UKcy5nZW5lcyA8LSBjYy5nZW5lcyRzLmdlbmVzCmcybS5nZW5lcyA8LSBjYy5nZW5lcyRnMm0uZ2VuZXMKCgpMMSA8LSBDZWxsQ3ljbGVTY29yaW5nKEwxLCBzLmZlYXR1cmVzID0gcy5nZW5lcywgCiAgICAgICAgICAgICAgICAgICAgICAgICAgIGcybS5mZWF0dXJlcyA9IGcybS5nZW5lcywgCiAgICAgICAgICAgICAgICAgICAgICAgICAgIHNldC5pZGVudCA9IFRSVUUpCgpEZWZhdWx0QXNzYXkoTDEpIDwtICJSTkEiCkwxJENDLkRpZmZlcmVuY2UgPC0gTDEkUy5TY29yZSAtIEwxJEcyTS5TY29yZQoKYGBgCgojIDQuIE5vcm1hbGl6ZSBkYXRhCmBgYHtyIE5vcm1hbGl6ZSwgZmlnLmhlaWdodD00LCBmaWcud2lkdGg9Nn0KCgojIEFwcGx5IFNDVHJhbnNmb3JtCkwxIDwtIFNDVHJhbnNmb3JtKEwxLCB2YXJzLnRvLnJlZ3Jlc3MgPSBjKCJwZXJjZW50LnJiIiwicGVyY2VudC5taXRvIiwgIkNDLkRpZmZlcmVuY2UiKSwgCiAgICAgICAgICAgICAgICAgIGRvLnNjYWxlPVRSVUUsIAogICAgICAgICAgICAgICAgICBkby5jZW50ZXI9VFJVRSwgCiAgICAgICAgICAgICAgICAgIHZlcmJvc2UgPSBUUlVFKQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIApgYGAKCgojIDUuIFBlcmZvcm0gUENBCmBgYHtyIFBDQSwgZmlnLmhlaWdodD00LCBmaWcud2lkdGg9Nn0KClZhcmlhYmxlc19nZW5lcyA8LSBMMUBhc3NheXMkU0NUQHZhci5mZWF0dXJlcwoKIyBFeGNsdWRlIGdlbmVzIHN0YXJ0aW5nIHdpdGggIkhMQS0iIEFORCAiWGlzdCIgQU5EICJUUkJWLCBUUkFWIgpWYXJpYWJsZXNfZ2VuZXNfYWZ0ZXJfZXhjbHVzaW9uIDwtIFZhcmlhYmxlc19nZW5lc1shZ3JlcGwoIl5ITEEtfF5YSVNUfF5UUkJWfF5UUkFWIiwgVmFyaWFibGVzX2dlbmVzKV0KCgojIFRoZXNlIGFyZSBub3cgc3RhbmRhcmQgc3RlcHMgaW4gdGhlIFNldXJhdCB3b3JrZmxvdyBmb3IgdmlzdWFsaXphdGlvbiBhbmQgY2x1c3RlcmluZwpMMSA8LSBSdW5QQ0EoTDEsCiAgICAgICAgICAgICBmZWF0dXJlcyA9IFZhcmlhYmxlc19nZW5lc19hZnRlcl9leGNsdXNpb24sCiAgICAgICAgICAgICBkby5wcmludCA9IFRSVUUsIAogICAgICAgICAgICAgcGNzLnByaW50ID0gMTo1LCAKICAgICAgICAgICAgIGdlbmVzLnByaW50ID0gMTUsCiAgICAgICAgICAgICBucGNzID0gNTApCgojIGRldGVybWluZSBkaW1lbnNpb25hbGl0eSBvZiB0aGUgZGF0YQpFbGJvd1Bsb3QoTDEsIG5kaW1zID01MCkKCgpgYGAKIyBQZXJmb3JtIFBDQSBURVNUCmBgYHtyIFBDQS1URVNULCBmaWcuaGVpZ2h0PTYsIGZpZy53aWR0aD0xMH0KCgpsaWJyYXJ5KGdncGxvdDIpCmxpYnJhcnkoUkNvbG9yQnJld2VyKSAgCgojIEFzc3VtaW5nIHlvdSBoYXZlIDEwIGRpZmZlcmVudCBjZWxsIGxpbmVzLCBnZW5lcmF0aW5nIGEgY29sb3IgcGFsZXR0ZSB3aXRoIDEwIGNvbG9ycwpjZWxsX2xpbmVfY29sb3JzIDwtIGJyZXdlci5wYWwoMTAsICJTZXQzIikKCiMgQXNzdW1pbmcgTDEkY2VsbF9saW5lIGlzIGEgZmFjdG9yIG9yIGNoYXJhY3RlciB2ZWN0b3IgY29udGFpbmluZyBjZWxsIGxpbmUgbmFtZXMKZGF0YSA8LSBhcy5kYXRhLmZyYW1lKHRhYmxlKEwxJGNlbGxfbGluZSkpCmNvbG5hbWVzKGRhdGEpIDwtIGMoImNlbGxfbGluZSIsICJuVU1JIikgICMgQ2hhbmdlIGNvbHVtbiBuYW1lIHRvIG5VTUkKCm5jZWxscyA8LSBnZ3Bsb3QoZGF0YSwgYWVzKHggPSBjZWxsX2xpbmUsIHkgPSBuVU1JLCBmaWxsID0gY2VsbF9saW5lKSkgKyAKICBnZW9tX2NvbCgpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIGdlb21fdGV4dChhZXMobGFiZWwgPSBuVU1JKSwgCiAgICAgICAgICAgIHBvc2l0aW9uID0gcG9zaXRpb25fZG9kZ2Uod2lkdGggPSAwLjkpLCAKICAgICAgICAgICAgdmp1c3QgPSAtMC4yNSkgKwogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IGNlbGxfbGluZV9jb2xvcnMpICsgCiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA0NSwgaGp1c3QgPSAxKSwKICAgICAgICBwbG90LnRpdGxlID0gZWxlbWVudF90ZXh0KGhqdXN0ID0gMC41KSkgKyAgIyBBZGp1c3QgdGhlIHRpdGxlIHBvc2l0aW9uCiAgZ2d0aXRsZSgiRmlsdGVyZWQgY2VsbHMgcGVyIHNhbXBsZSIpICsKICB4bGFiKCJDZWxsIGxpbmVzIikgKyAgIyBBZGp1c3QgeC1heGlzIGxhYmVsCiAgeWxhYigiRnJlcXVlbmN5IikgICAgIyBBZGp1c3QgeS1heGlzIGxhYmVsCgpwcmludChuY2VsbHMpCgoKCiMgVEVTVC0xCiMgZ2l2ZW4gdGhhdCB0aGUgb3V0cHV0IG9mIFJ1blBDQSBpcyAicGNhIgojIHJlcGxhY2UgInNvIiBieSB0aGUgbmFtZSBvZiB5b3VyIHNldXJhdCBvYmplY3QKCnBjdCA8LSBMMVtbInBjYSJdXUBzdGRldiAvIHN1bShMMVtbInBjYSJdXUBzdGRldikgKiAxMDAKY3VtdSA8LSBjdW1zdW0ocGN0KSAjIENhbGN1bGF0ZSBjdW11bGF0aXZlIHBlcmNlbnRzIGZvciBlYWNoIFBDCiMgRGV0ZXJtaW5lIHRoZSBkaWZmZXJlbmNlIGJldHdlZW4gdmFyaWF0aW9uIG9mIFBDIGFuZCBzdWJzZXF1ZW50IFBDCmNvMiA8LSBzb3J0KHdoaWNoKChwY3RbLWxlbmd0aChwY3QpXSAtIHBjdFstMV0pID4gMC4xKSwgZGVjcmVhc2luZyA9IFQpWzFdICsgMQojIGxhc3QgcG9pbnQgd2hlcmUgY2hhbmdlIG9mICUgb2YgdmFyaWF0aW9uIGlzIG1vcmUgdGhhbiAwLjElLiAtPiBjbzIKY28yCgojIFRFU1QtMgojIGdldCBzaWduaWZpY2FudCBQQ3MKc3RkdiA8LSBMMVtbInBjYSJdXUBzdGRldgpzdW0uc3RkdiA8LSBzdW0oTDFbWyJwY2EiXV1Ac3RkZXYpCnBlcmNlbnQuc3RkdiA8LSAoc3RkdiAvIHN1bS5zdGR2KSAqIDEwMApjdW11bGF0aXZlIDwtIGN1bXN1bShwZXJjZW50LnN0ZHYpCmNvMSA8LSB3aGljaChjdW11bGF0aXZlID4gOTAgJiBwZXJjZW50LnN0ZHYgPCA1KVsxXQpjbzIgPC0gc29ydCh3aGljaCgocGVyY2VudC5zdGR2WzE6bGVuZ3RoKHBlcmNlbnQuc3RkdikgLSAxXSAtIAogICAgICAgICAgICAgICAgICAgICAgIHBlcmNlbnQuc3RkdlsyOmxlbmd0aChwZXJjZW50LnN0ZHYpXSkgPiAwLjEpLCAKICAgICAgICAgICAgICBkZWNyZWFzaW5nID0gVClbMV0gKyAxCm1pbi5wYyA8LSBtaW4oY28xLCBjbzIpCm1pbi5wYwoKIyBDcmVhdGUgYSBkYXRhZnJhbWUgd2l0aCB2YWx1ZXMKcGxvdF9kZiA8LSBkYXRhLmZyYW1lKHBjdCA9IHBlcmNlbnQuc3RkdiwgCiAgICAgICAgICAgY3VtdSA9IGN1bXVsYXRpdmUsIAogICAgICAgICAgIHJhbmsgPSAxOmxlbmd0aChwZXJjZW50LnN0ZHYpKQoKIyBFbGJvdyBwbG90IHRvIHZpc3VhbGl6ZSAKICBnZ3Bsb3QocGxvdF9kZiwgYWVzKGN1bXVsYXRpdmUsIHBlcmNlbnQuc3RkdiwgbGFiZWwgPSByYW5rLCBjb2xvciA9IHJhbmsgPiBtaW4ucGMpKSArIAogIGdlb21fdGV4dCgpICsgCiAgZ2VvbV92bGluZSh4aW50ZXJjZXB0ID0gOTAsIGNvbG9yID0gImdyZXkiKSArIAogIGdlb21faGxpbmUoeWludGVyY2VwdCA9IG1pbihwZXJjZW50LnN0ZHZbcGVyY2VudC5zdGR2ID4gNV0pLCBjb2xvciA9ICJncmV5IikgKwogIHRoZW1lX2J3KCkKCiAgCgpgYGAKCiMgNi4gQ2x1c3RlcmluZwpgYGB7ciBDMSwgZmlnLmhlaWdodD00LCBmaWcud2lkdGg9Nn0KTDEgPC0gRmluZE5laWdoYm9ycyhMMSwgCiAgICAgICAgICAgICAgICAgICAgZGltcyA9IDE6bWluLnBjLCAKICAgICAgICAgICAgICAgICAgICB2ZXJib3NlID0gRkFMU0UpCgojIHVuZGVyc3RhbmRpbmcgcmVzb2x1dGlvbgpMMSA8LSBGaW5kQ2x1c3RlcnMoTDEsIAogICAgICAgICAgICAgICAgICByZXNvbHV0aW9uID0gYygwLjEsIDAuMiwgMC4zLCAwLjQsIDAuNSwgMC42LCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAwLjcsMC44LCAwLjksIDEsIDEuMSwgMS4yKSkKCgpgYGAKCgpgYGB7ciBDMiwgZmlnLmhlaWdodD00LCBmaWcud2lkdGg9Nn0KCiMgbm9uLWxpbmVhciBkaW1lbnNpb25hbGl0eSByZWR1Y3Rpb24gLS0tLS0tLS0tLS0tLS0KTDEgPC0gUnVuVU1BUChMMSwgCiAgICAgICAgICAgICAgZGltcyA9IDE6bWluLnBjLAogICAgICAgICAgICAgIHZlcmJvc2UgPSBGQUxTRSkKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAoKIyBub3RlIHRoYXQgeW91IGNhbiBzZXQgYGxhYmVsID0gVFJVRWAgb3IgdXNlIHRoZSBMYWJlbCBDbHVzdGVycyBmdW5jdGlvbiB0byBoZWxwIGxhYmVsCiMgaW5kaXZpZHVhbCBjbHVzdGVycwpEaW1QbG90KEwxLGdyb3VwLmJ5ID0gImNlbGxfbGluZSIsIAogICAgICAgIHJlZHVjdGlvbiA9ICJ1bWFwIiwKICAgICAgICBsYWJlbC5zaXplID0gMywKICAgICAgICByZXBlbCA9IFQsCiAgICAgICAgbGFiZWwgPSBULCBsYWJlbC5ib3ggPSBUKQoKCkRpbVBsb3QoTDEsCiAgICAgICAgZ3JvdXAuYnkgPSAiU0NUX3Nubl9yZXMuMC4xIiwgCiAgICAgICAgcmVkdWN0aW9uID0gInVtYXAiLAogICAgICAgIGxhYmVsLnNpemUgPSAzLAogICAgICAgIHJlcGVsID0gVCwKICAgICAgICBsYWJlbCA9IFQsIGxhYmVsLmJveCA9IFQpCgpEaW1QbG90KEwxLAogICAgICAgIGdyb3VwLmJ5ID0gIlNDVF9zbm5fcmVzLjAuMiIsIAogICAgICAgIHJlZHVjdGlvbiA9ICJ1bWFwIiwKICAgICAgICBsYWJlbC5zaXplID0gMywKICAgICAgICByZXBlbCA9IFQsCiAgICAgICAgbGFiZWwgPSBULCBsYWJlbC5ib3ggPSBUKQoKRGltUGxvdChMMSwKICAgICAgICBncm91cC5ieSA9ICJTQ1Rfc25uX3Jlcy4wLjMiLCAKICAgICAgICByZWR1Y3Rpb24gPSAidW1hcCIsCiAgICAgICAgbGFiZWwuc2l6ZSA9IDMsCiAgICAgICAgcmVwZWwgPSBULAogICAgICAgIGxhYmVsID0gVCwgbGFiZWwuYm94ID0gVCkKCkRpbVBsb3QoTDEsCiAgICAgICAgZ3JvdXAuYnkgPSAiU0NUX3Nubl9yZXMuMC40IiwgCiAgICAgICAgcmVkdWN0aW9uID0gInVtYXAiLAogICAgICAgIGxhYmVsLnNpemUgPSAzLAogICAgICAgIHJlcGVsID0gVCwKICAgICAgICBsYWJlbCA9IFQsIGxhYmVsLmJveCA9IFQpCgoKRGltUGxvdChMMSwKICAgICAgICBncm91cC5ieSA9ICJTQ1Rfc25uX3Jlcy4wLjUiLCAKICAgICAgICByZWR1Y3Rpb24gPSAidW1hcCIsCiAgICAgICAgbGFiZWwuc2l6ZSA9IDMsCiAgICAgICAgcmVwZWwgPSBULAogICAgICAgIGxhYmVsID0gVCwgbGFiZWwuYm94ID0gVCkKCkRpbVBsb3QoTDEsCiAgICAgICAgZ3JvdXAuYnkgPSAiU0NUX3Nubl9yZXMuMC42IiwgCiAgICAgICAgcmVkdWN0aW9uID0gInVtYXAiLAogICAgICAgIGxhYmVsLnNpemUgPSAzLAogICAgICAgIHJlcGVsID0gVCwKICAgICAgICBsYWJlbCA9IFQsIGxhYmVsLmJveCA9IFQpCgpEaW1QbG90KEwxLAogICAgICAgIGdyb3VwLmJ5ID0gIlNDVF9zbm5fcmVzLjAuNyIsIAogICAgICAgIHJlZHVjdGlvbiA9ICJ1bWFwIiwKICAgICAgICBsYWJlbC5zaXplID0gMywKICAgICAgICByZXBlbCA9IFQsCiAgICAgICAgbGFiZWwgPSBULCBsYWJlbC5ib3ggPSBUKQoKRGltUGxvdChMMSwKICAgICAgICBncm91cC5ieSA9ICJTQ1Rfc25uX3Jlcy4wLjgiLCAKICAgICAgICByZWR1Y3Rpb24gPSAidW1hcCIsCiAgICAgICAgbGFiZWwuc2l6ZSA9IDMsCiAgICAgICAgcmVwZWwgPSBULAogICAgICAgIGxhYmVsID0gVCwgbGFiZWwuYm94ID0gVCkKRGltUGxvdChMMSwKICAgICAgICBncm91cC5ieSA9ICJTQ1Rfc25uX3Jlcy4wLjkiLCAKICAgICAgICByZWR1Y3Rpb24gPSAidW1hcCIsCiAgICAgICAgbGFiZWwuc2l6ZSA9IDMsCiAgICAgICAgcmVwZWwgPSBULAogICAgICAgIGxhYmVsID0gVCwgbGFiZWwuYm94ID0gVCkKRGltUGxvdChMMSwKICAgICAgICBncm91cC5ieSA9ICJTQ1Rfc25uX3Jlcy4xIiwgCiAgICAgICAgcmVkdWN0aW9uID0gInVtYXAiLAogICAgICAgIGxhYmVsLnNpemUgPSAzLAogICAgICAgIHJlcGVsID0gVCwKICAgICAgICBsYWJlbCA9IFQsIGxhYmVsLmJveCA9IFQpCkRpbVBsb3QoTDEsCiAgICAgICAgZ3JvdXAuYnkgPSAiU0NUX3Nubl9yZXMuMS4xIiwgCiAgICAgICAgcmVkdWN0aW9uID0gInVtYXAiLAogICAgICAgIGxhYmVsLnNpemUgPSAzLAogICAgICAgIHJlcGVsID0gVCwKICAgICAgICBsYWJlbCA9IFQsIGxhYmVsLmJveCA9IFQpCkRpbVBsb3QoTDEsCiAgICAgICAgZ3JvdXAuYnkgPSAiU0NUX3Nubl9yZXMuMS4yIiwgCiAgICAgICAgcmVkdWN0aW9uID0gInVtYXAiLAogICAgICAgIGxhYmVsLnNpemUgPSAzLAogICAgICAgIHJlcGVsID0gVCwKICAgICAgICBsYWJlbCA9IFQsIGxhYmVsLmJveCA9IFQpCgpgYGAKCiMgNy4gQXppbXV0aCBBbm5vdGF0aW9uCmBgYHtyIGF6aW11dGhfQW5ub3RhdGlvbiwgZmlnLmhlaWdodD00LCBmaWcud2lkdGg9Nn0KSW5zdGFsbERhdGEoInBibWNyZWYiKQoKIyBUaGUgUnVuQXppbXV0aCBmdW5jdGlvbiBjYW4gdGFrZSBhIFNldXJhdCBvYmplY3QgYXMgaW5wdXQKTDEgPC0gUnVuQXppbXV0aChMMSwgcmVmZXJlbmNlID0gInBibWNyZWYiKQoKVU1BUFBsb3QoTDEsIGdyb3VwLmJ5ID0gInByZWRpY3RlZC5jZWxsdHlwZS5sMiIsIAogICAgICAgIHJlZHVjdGlvbiA9ICJ1bWFwIiwKICAgICAgICBsYWJlbC5zaXplID0gMywKICAgICAgICByZXBlbCA9IFQsCiAgICAgICAgbGFiZWwgPSBGKQoKClVNQVBQbG90KEwxLCBncm91cC5ieSA9ICJwcmVkaWN0ZWQuY2VsbHR5cGUubDIiLCAKICAgICAgICByZWR1Y3Rpb24gPSAidW1hcCIsCiAgICAgICAgbGFiZWwuc2l6ZSA9IDMsCiAgICAgICAgcmVwZWwgPSBULAogICAgICAgIGxhYmVsID0gVCkKYGBgCgojIDguIENlbGwgdHlwZSBhbm5vdGF0aW9uIHVzaW5nIFByb2plY3RUaWxzCmBgYHtyIFByb2plY1RpbHMsIGZpZy5oZWlnaHQ9NCwgZmlnLndpZHRoPTZ9CgoKI0xvYWQgcmVmZXJlbmNlIGF0bGFzIGFuZCBxdWVyeSBkYXRhCnJlZiA8LSByZWFkUkRTKGZpbGUgPSAiQ0Q0VF9odW1hbl9yZWZfdjEucmRzIikKCiNSdW4gUHJvamVjdGlvbiBhbGdvcml0aG0KcXVlcnkucHJvamVjdGVkIDwtIFJ1bi5Qcm9qZWNUSUxzKEwxLCByZWYgPSByZWYpCgoKCiNyZWZlcmVuY2UgYXRsYXMKRGltUGxvdChyZWYsIGxhYmVsID0gVCkKCiNWaXN1YWxpemUgcHJvamVjdGlvbgpwbG90LnByb2plY3Rpb24ocmVmLCBxdWVyeS5wcm9qZWN0ZWQsIGxpbmVzaXplID0gMC41LCBwb2ludHNpemUgPSAwLjUpCgojUGxvdCB0aGUgcHJlZGljdGVkIGNvbXBvc2l0aW9uIG9mIHRoZSBxdWVyeSBpbiB0ZXJtcyBvZiByZWZlcmVuY2UgVCBjZWxsIHN1YnR5cGVzCnBsb3Quc3RhdGVwcmVkLmNvbXBvc2l0aW9uKHJlZiwgcXVlcnkucHJvamVjdGVkLCBtZXRyaWMgPSAiUGVyY2VudCIpCgoKTDEgPC0gUHJvamVjVElMcy5jbGFzc2lmaWVyKHF1ZXJ5ID0gTDEsIHJlZiA9IHJlZikKVU1BUFBsb3QoTDEsIGdyb3VwLmJ5ID0gImZ1bmN0aW9uYWwuY2x1c3RlciIsIAogICAgICAgIHJlZHVjdGlvbiA9ICJ1bWFwIiwKICAgICAgICBsYWJlbC5zaXplID0gMywKICAgICAgICByZXBlbCA9IFQsCiAgICAgICAgbGFiZWwgPSBUKQoKCmBgYAoKIyA5LkNlbGwgdHlwZSBhbm5vdGF0aW9uIHVzaW5nIFNpbmdsZVIKYGBge3Igc2luZ2xlUiwgZmlnLmhlaWdodD00LCBmaWcud2lkdGg9Nn0KI2dldCByZWZlcmVuY2UgZGF0YXNldHMgZnJvbSBjZWxsZGV4IHBhY2thZ2UKCm1vbmFjby5yZWYgPC0gY2VsbGRleDo6TW9uYWNvSW1tdW5lRGF0YSgpCmhwY2EucmVmIDwtIGNlbGxkZXg6Okh1bWFuUHJpbWFyeUNlbGxBdGxhc0RhdGEoKQpkaWNlLnJlZiA8LSBjZWxsZGV4OjpEYXRhYmFzZUltbXVuZUNlbGxFeHByZXNzaW9uRGF0YSgpCmJwZS5yZWYgPC0gY2VsbGRleDo6Qmx1ZXByaW50RW5jb2RlRGF0YSgpCgojY29udmVydCBvdXIgU2V1cmF0IG9iamVjdCB0byBzaW5nbGUgY2VsbCBleHBlcmltZW50IChTQ0UpCnNjZSA8LSBhcy5TaW5nbGVDZWxsRXhwZXJpbWVudChEaWV0U2V1cmF0KEwxKSkKCiN1c2luZyBTaW5nbGVSCm1vbmFjby5tYWluIDwtIFNpbmdsZVIodGVzdCA9IHNjZSxhc3NheS50eXBlLnRlc3QgPSAxLHJlZiA9IG1vbmFjby5yZWYsbGFiZWxzID0gbW9uYWNvLnJlZiRsYWJlbC5tYWluKQptb25hY28uZmluZSA8LSBTaW5nbGVSKHRlc3QgPSBzY2UsYXNzYXkudHlwZS50ZXN0ID0gMSxyZWYgPSBtb25hY28ucmVmLGxhYmVscyA9IG1vbmFjby5yZWYkbGFiZWwuZmluZSkKaHBjYS5tYWluIDwtIFNpbmdsZVIodGVzdCA9IHNjZSxhc3NheS50eXBlLnRlc3QgPSAxLHJlZiA9IGhwY2EucmVmLGxhYmVscyA9IGhwY2EucmVmJGxhYmVsLm1haW4pCmhwY2EuZmluZSA8LSBTaW5nbGVSKHRlc3QgPSBzY2UsYXNzYXkudHlwZS50ZXN0ID0gMSxyZWYgPSBocGNhLnJlZixsYWJlbHMgPSBocGNhLnJlZiRsYWJlbC5maW5lKQpkaWNlLm1haW4gPC0gU2luZ2xlUih0ZXN0ID0gc2NlLGFzc2F5LnR5cGUudGVzdCA9IDEscmVmID0gZGljZS5yZWYsbGFiZWxzID0gZGljZS5yZWYkbGFiZWwubWFpbikKZGljZS5maW5lIDwtIFNpbmdsZVIodGVzdCA9IHNjZSxhc3NheS50eXBlLnRlc3QgPSAxLHJlZiA9IGRpY2UucmVmLGxhYmVscyA9IGRpY2UucmVmJGxhYmVsLmZpbmUpCmJwZS5tYWluIDwtIFNpbmdsZVIodGVzdCA9IHNjZSxhc3NheS50eXBlLnRlc3QgPSAxLHJlZiA9IGJwZS5yZWYsbGFiZWxzID0gYnBlLnJlZiRsYWJlbC5tYWluKQpicGUuZmluZSA8LSBTaW5nbGVSKHRlc3QgPSBzY2UsYXNzYXkudHlwZS50ZXN0ID0gMSxyZWYgPSBicGUucmVmLGxhYmVscyA9IGJwZS5yZWYkbGFiZWwuZmluZSkKCgojc3VtbWFyeSBvZiBnZW5lcmFsIGNlbGwgdHlwZSBhbm5vdGF0aW9ucwoKI3RhYmxlKG1vbmFjby5tYWluJHBydW5lZC5sYWJlbHMpCiN0YWJsZShocGNhLm1haW4kcHJ1bmVkLmxhYmVscykKI3RhYmxlKGRpY2UubWFpbiRwcnVuZWQubGFiZWxzKQojdGFibGUoYnBlLm1haW4kcHJ1bmVkLmxhYmVscykKCiNUaGUgZmluZXIgY2VsbCB0eXBlcyBhbm5vdGF0aW9ucyBhcmUgeW91IGFmdGVyLCB0aGUgaGFyZGVyIHRoZXkgYXJlIHRvIGdldCByZWxpYWJseS4gCiNUaGlzIGlzIHdoZXJlIGNvbXBhcmluZyBtYW55IGRhdGFiYXNlcywgYXMgd2VsbCBhcyB1c2luZyBpbmRpdmlkdWFsIG1hcmtlcnMgZnJvbSBsaXRlcmF0dXJlLCAKI3dvdWxkIGFsbCBiZSB2ZXJ5IHZhbHVhYmxlLgoKI3RhYmxlKG1vbmFjby5maW5lJHBydW5lZC5sYWJlbHMpCiN0YWJsZShocGNhLmZpbmUkcHJ1bmVkLmxhYmVscykKI3RhYmxlKGRpY2UuZmluZSRwcnVuZWQubGFiZWxzKQojdGFibGUoYnBlLmZpbmUkcHJ1bmVkLmxhYmVscykKCgoKI2FkZCB0aGUgYW5ub3RhdGlvbnMgdG8gdGhlIFNldXJhdCBvYmplY3QgbWV0YWRhdGEKTDFAbWV0YS5kYXRhJG1vbmFjby5tYWluIDwtIG1vbmFjby5tYWluJHBydW5lZC5sYWJlbHMKTDFAbWV0YS5kYXRhJG1vbmFjby5maW5lIDwtIG1vbmFjby5maW5lJHBydW5lZC5sYWJlbHMKIwpMMUBtZXRhLmRhdGEkaHBjYS5tYWluICAgPC0gaHBjYS5tYWluJHBydW5lZC5sYWJlbHMKTDFAbWV0YS5kYXRhJGhwY2EuZmluZSAgIDwtIGhwY2EuZmluZSRwcnVuZWQubGFiZWxzCiMgIApMMUBtZXRhLmRhdGEkZGljZS5tYWluICAgPC0gZGljZS5tYWluJHBydW5lZC5sYWJlbHMKTDFAbWV0YS5kYXRhJGRpY2UuZmluZSAgIDwtIGRpY2UuZmluZSRwcnVuZWQubGFiZWxzCiMgCkwxQG1ldGEuZGF0YSRicGUubWFpbiAgIDwtIGJwZS5tYWluJHBydW5lZC5sYWJlbHMKTDFAbWV0YS5kYXRhJGJwZS5maW5lICAgPC0gYnBlLmZpbmUkcHJ1bmVkLmxhYmVscwoKIyBNb25hY28gQW5ub3RhdGlvbnMKRGltUGxvdChMMSwgZ3JvdXAuYnkgPSAibW9uYWNvLm1haW4iLCAKICAgICAgICByZWR1Y3Rpb24gPSAidW1hcCIsCiAgICAgICAgbGFiZWwuc2l6ZSA9IDMsCiAgICAgICAgcmVwZWwgPSBULAogICAgICAgIGxhYmVsID0gRikKCkRpbVBsb3QoTDEsIGdyb3VwLmJ5ID0gIm1vbmFjby5tYWluIiwgCiAgICAgICAgcmVkdWN0aW9uID0gInVtYXAiLAogICAgICAgIGxhYmVsLnNpemUgPSAzLAogICAgICAgIHJlcGVsID0gVCwKICAgICAgICBsYWJlbCA9IFQsIGxhYmVsLmJveCA9IFQpCgpEaW1QbG90KEwxLCBncm91cC5ieSA9ICJtb25hY28uZmluZSIsIAogICAgICAgIHJlZHVjdGlvbiA9ICJ1bWFwIiwKICAgICAgICBsYWJlbC5zaXplID0gMywKICAgICAgICByZXBlbCA9IFQsCiAgICAgICAgbGFiZWwgPSBGKQoKRGltUGxvdChMMSwgZ3JvdXAuYnkgPSAibW9uYWNvLmZpbmUiLCAKICAgICAgICByZWR1Y3Rpb24gPSAidW1hcCIsCiAgICAgICAgbGFiZWwuc2l6ZSA9IDMsCiAgICAgICAgcmVwZWwgPSBULAogICAgICAgIGxhYmVsID0gVCwgbGFiZWwuYm94ID0gVCkKCiMgSFBDQSBBbm5vdGF0aW9ucwpEaW1QbG90KEwxLCBncm91cC5ieSA9ICJocGNhLm1haW4iLCAKICAgICAgICByZWR1Y3Rpb24gPSAidW1hcCIsCiAgICAgICAgbGFiZWwuc2l6ZSA9IDMsCiAgICAgICAgcmVwZWwgPSBULAogICAgICAgIGxhYmVsID0gVCwgbGFiZWwuYm94ID0gVCkKCkRpbVBsb3QoTDEsIGdyb3VwLmJ5ID0gImhwY2EuZmluZSIsIAogICAgICAgIHJlZHVjdGlvbiA9ICJ1bWFwIiwKICAgICAgICBsYWJlbC5zaXplID0gMywKICAgICAgICByZXBlbCA9IFQsCiAgICAgICAgbGFiZWwgPSBGKQoKRGltUGxvdChMMSwgZ3JvdXAuYnkgPSAiaHBjYS5maW5lIiwgCiAgICAgICAgcmVkdWN0aW9uID0gInVtYXAiLAogICAgICAgIGxhYmVsLnNpemUgPSAzLAogICAgICAgIHJlcGVsID0gVCwKICAgICAgICBsYWJlbCA9IFQsIGxhYmVsLmJveCA9IFQpCiMgRElDRSBBbm5vdGF0aW9ucwpEaW1QbG90KEwxLCBncm91cC5ieSA9ICJkaWNlLm1haW4iLCAKICAgICAgICByZWR1Y3Rpb24gPSAidW1hcCIsCiAgICAgICAgbGFiZWwuc2l6ZSA9IDMsCiAgICAgICAgcmVwZWwgPSBULAogICAgICAgIGxhYmVsID0gVCwgbGFiZWwuYm94ID0gVCkKCkRpbVBsb3QoTDEsIGdyb3VwLmJ5ID0gImRpY2UuZmluZSIsIAogICAgICAgIHJlZHVjdGlvbiA9ICJ1bWFwIiwKICAgICAgICBsYWJlbC5zaXplID0gMywKICAgICAgICByZXBlbCA9IFQsCiAgICAgICAgbGFiZWwgPSBGKQoKRGltUGxvdChMMSwgZ3JvdXAuYnkgPSAiZGljZS5maW5lIiwgCiAgICAgICAgcmVkdWN0aW9uID0gInVtYXAiLAogICAgICAgIGxhYmVsLnNpemUgPSAzLAogICAgICAgIHJlcGVsID0gVCwKICAgICAgICBsYWJlbCA9IFQsIGxhYmVsLmJveCA9IFQpCgojIEJQRSBBbm5vdGF0aW9ucwpEaW1QbG90KEwxLCBncm91cC5ieSA9ICJicGUubWFpbiIsIAogICAgICAgIHJlZHVjdGlvbiA9ICJ1bWFwIiwKICAgICAgICBsYWJlbC5zaXplID0gMywKICAgICAgICByZXBlbCA9IFQsCiAgICAgICAgbGFiZWwgPSBULCBsYWJlbC5ib3ggPSBUKQoKRGltUGxvdChMMSwgZ3JvdXAuYnkgPSAiYnBlLmZpbmUiLCAKICAgICAgICByZWR1Y3Rpb24gPSAidW1hcCIsCiAgICAgICAgbGFiZWwuc2l6ZSA9IDMsCiAgICAgICAgcmVwZWwgPSBULAogICAgICAgIGxhYmVsID0gRikKCkRpbVBsb3QoTDEsIGdyb3VwLmJ5ID0gImJwZS5maW5lIiwgCiAgICAgICAgcmVkdWN0aW9uID0gInVtYXAiLAogICAgICAgIGxhYmVsLnNpemUgPSAzLAogICAgICAgIHJlcGVsID0gVCwKICAgICAgICBsYWJlbCA9IFQsIGxhYmVsLmJveCA9IFQpCgoKCmBgYAoKIyAxMC4gY2x1c1RyZWUKYGBge3IgY2x1c1RyZWUsIGZpZy5oZWlnaHQ9MTIsIGZpZy53aWR0aD0xMH0KbGlicmFyeShjbHVzdHJlZSkKY2x1c3RyZWUoTDEsIHByZWZpeCA9ICJTQ1Rfc25uX3Jlcy4iKQpgYGAKCiMgMTEuU2F2ZSB0aGUgU2V1cmF0IG9iamVjdCBhcyBhbiBSb2JqIGZpbGUKYGBge3Igc2F2ZVJPQkosIGVjaG89RkFMU0V9CgpzYXZlKEwxLCBmaWxlID0gIkwxX0FuYWx5c2lzLlJvYmoiKQoKCmBgYAoKCgoKCgoK