$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Visão geral da estratégia de visualização, quantificação e mapeamento de populações celulares de interesse no TME
Para quantificar as populações celulares de interesse (COIs) em diferentes compartimentos teciduais (TCs) e caracterizar sua organização espacial, projetamos um fluxo de trabalho que integra técnicas acessíveis e fáceis de usar e maximiza as informações posicionais que podem ser obtidas a partir de preciosas amostras clínicas de FFPE(Figura 1). Em primeiro lugar, as seções de FFPE de tecido inteiro serial foram manchadas para visualização de COIs (por exemplo, células imunes) e TCs (por exemplo, estroma versus parenchyma) (Figura 1, passo 1). O número de seções consecutivas a serem manchadas deve ser mantido ao mínimo que permita a visualização das células de interesse ou características teciduais necessárias para abordar a questão da pesquisa. Quanto menor o número de seções seriais, maior a semelhança da arquitetura tecidual e a concordância entre seções contíguas. Além disso, a capacidade de multiplexação pode ser expandida através da reutilização de seções manchadas fluorescentes através de técnicas de descascamento e ressarcimento19.
Uma vez que as etapas de coloração foram feitas, um scanner de slides inteiro foi usado para digitalizar as imagens. As imagens adquiridas a partir de seções seriais foram alinhadas e consolidadas em um slide multiplex virtual de forma automatizada(Figura 1, seção 2). Em seguida, um ROI para o tecido foi delineado com um protocolo definido pelo usuário que identificou pixels associados ao tecido (TAPs) (Figura 1, passo 3). Posteriormente, o tecido de ROI foi segmentado em TCs definidos como ROIs adicionais. (Figura 1, passo 4). Em seguida, os protocolos definidos pelo usuário detectaram e quantificaram cois em diferentes TCs(Figura 1, etapa 5). Finalmente, foram gerados mapas térmicos de tecidos de COIs com base em suas densidades e suas coordenadas teciduais (Figura 1, passo 6).

Figura 1: Representação esquemática da estratégia de visualização, quantificação e mapeamento de células imunes no TME. (1) As seções de tecido sumidérdia foram manchadas para rotular COIs e TCs. As seções de tecido inteiros manchadas foram digitalizadas usando um scanner slide inteiro. (2) As imagens adquiridas a partir de seções seriais foram vinculadas, alinhadas e co-registradas de forma automatizada, utilizando um módulo de análise Tissuealign. Uma imagem composta foi gerada a partir do alinhamento de alta precisão de imagens individuais. (3) Foi utilizado um protocolo definido pelo usuário para detecção automatizada de pixels associados ao tecido (TAPs) na imagem composta. (4) O tecido foi segmentado em TCs (por exemplo, estroma e parenchyma) definidos como ROIs. (5) Foram utilizados protocolos definidos pelo usuário para a detecção e quantificação automatizada de COIs em diferentes TCs. (6) Foram gerados mapas térmicos de COIs de tecido. Clique aqui para ver uma versão maior desta figura.
COIs e TCs de imagem
Três seções de tecido inteiro de FFPE serial de tumor ressecado de um sujeito com carcinoma hepatocelular associado ao HBV foram manchadas em uma ou mais rodadas de coloração como na Figura 2A. A seção I foi manchada com H&E para mostrar a arquitetura tecidual, a morfologia celular e determinar parâmetros clinicamente relevantes, como tipo de malignidade, grau de tumor e avaliação geral da infiltração imunológica(Figura 2C). Na seção Contígua II, foram utilizadas duas rodadas de IFM foram utilizadas para rotular células parênquias hepáticas e não parênquias(Figura 2A). Na primeira rodada, os vasos normais e tumorais foram visualizados utilizando-se a coloração de CD34 de células endoteliais. Além disso, as células epitelial (hepatócitos e colangiocitos) foram identificadas por meio de citoqueratina 8/18, e as células hepáticas hepáticas ativadas fibrogênicas foram identificadas como células alfa lisas do músculo (αSMA+)(Figura 2C). Após a aquisição da imagem, as seções teciduais foram despojadas e rescadas com anticorpos contra macrófagos (CD68) e miofibroblastos (desmin). Para melhor caracterizar o infiltrado imunológico tumoral, a seção serial adjacente III foi manchada usando duas rodadas de IFM para os marcadores celulares CD3, CD4, CD8, bifurcação P3 (FoxP3) e mieloperoxidase (MPO). Em todos os casos, o DAPI foi usado como contra-mancha nuclear. Finalmente, a seção III foi manchada com mancha PSR e contrariada com verde rápido para visualizar o colágeno fibrilar e segmentar o tecido em estroma e parenchyma(Figura 2C).
Todo um scanner slide equipado com uma lente objetiva de 20X foi usado para digitalizar seções manchadas e criar slides virtuais. Seis imagens foram adquiridas a partir das três seções seriais (Figura 2B) e os slides virtuais posteriormente analisados utilizando o software VIS de acordo com a representação esquemmática na Figura 1.
Análise de Imagens
A análise da imagem compreendeu cinco etapas: 1) alinhamento tecidual; 2) detecção de tecido; 3) segmentação de tecidos; 4) quantificação automatizada de COIs; e 5) mapeamento térmico do tecido. Todos os protocolos para análise de imagem foram desenvolvidos utilizando o módulo Autor do software de análise de imagem e são referidos no texto como APP.
Alinhamento tecidual
Seis slides virtuais de três seções seriais, abrangendo 11 marcadores mais manchas de H&E e PSR, foram carregados no módulo Tissualign do software de análise de imagens. Em seguida, as imagens foram vinculadas, alinhadas e co-registradas de forma automatizada, gerando uma imagem composta virtual de 11 plex mais H&E e PSR, contendo todas as camadas das imagens individuais(Figuras 2A-C). O alinhamento foi preciso no caso de imagens originárias de seções seriais adjacentes, mostrando estruturas teciduais correspondentes posicionadas e dispostas de forma homóloga no alinhamento(Figura 2C e Figura S1A). Além disso, o alinhamento foi preciso no nível individual da célula para imagens originárias da mesma seção(Figura S1B). O tempo para alinhamento automático depende do número, tamanho, complexidade e semelhança das imagens a serem alinhadas. O alinhamento dos seis slides virtuais acima mencionados levou 15 minutos em nossa estação VIS.

Figura 2: Coloração de seções de tecido serial e alinhamento de imagem. (A) Resumo das manchas feitas em três seções seriais para visualização de COIs e TCs. Os números entre parênteses indicam a designação da imagem. Para as seções II e III, os tecidos foram despidos e ressonados com um segundo coquetel de anticorpos. (B) Visão geral de seis imagens individuais de tecido inteiro antes e depois do alinhamento do tecido (esquerda e direita, respectivamente). Barra de escala = 3.500 μm.(C) Visão ampliada de imagens alinhadas. Barra de escala = 80 μm. Clique aqui para ver uma versão maior desta figura.
Detecção de tecidos
Uma vez que as imagens foram vinculadas e alinhadas, procuramos identificar os TAPs(Figura 3A). Para projetar um APP para a detecção automatizada de TAPs (APP 1, Tabela 1),aproveitamos duas propriedades que diferenciam os TAPs de pixels não associados ao tecido. Primeiro, o sinal DAPI (banda azul) é restrito aos núcleos, que estão localizados exclusivamente no tecido, o que significa que todos os pixels DAPI+ são um subconjunto de TAPs. Em segundo lugar, os TAPs têm maior sinal de autofluorescência nas faixas verde e amarela em comparação com pixels não associados ao tecido. Consequentemente, desenvolvemos o APP 1 para detecção de tecidos(Tabela 1),que detecta os TAPs com base no sinal de linha de base nesses canais usando técnicas simples de limiar. Limiares para as faixas azul, verde e amarela foram definidos de modo que os TAPs tinham valores de intensidade de fundo acima dos limiares, enquanto os pixels não associados ao tecido tinham valores abaixo. O APP 1 para detecção de tecido foi aplicado à imagem IIA, que contém camadas nos canais azul, verde e amarelo(Figura 3A). Como saídas do APP 1, uma máscara verde brilhante foi colocada no topo dos TAPs, e um ROI chamado "Tissue" foi delineado (saída, Figura 3A). Além disso, a área do tecido foi determinada como variável de saída quantitativa. Como o APP 1 não incorpora os pixels não associados ao tecido no tecido DO ROI, eles foram excluídos da análise subseqüente com base neste ROI(Figura 3A). A precisão do APP 1 na identificação de TAPs é mostrada na Figura 3A.
Segmentação de tecidos e delineamento de ROIs para TCs
Em seguida, começamos a definir diferentes compartimentos dentro do tecido de ROI segmentando o tecido em estroma versus parenchyma. Utilizou-se a imagem manchada de PSR (IIIC, Figura 2C), onde o estroma pode ser definido como a área associada à deposição de colagens fibrilais (banda vermelha), a parenchyma como a área onde faltam as colagens fibrilas, e prevalece o rápido tine de contracoloração verde (faixa verde)(Figura 3B). Criamos o APP 2(Tabela 1)para delimitar digitalmente os TCs Stroma e Parenchyma. Este APP funciona no tecido DE ROI predefinido (saída, Figura 3A) e utiliza áreas representativas de estroma e parenchyma para treinar a ferramenta Classifier integrada no módulo de Análise de Imagens. O Classificador treinado atribui os pixels a um estroma ou a um rótulo parenchyma (salmão e verde, respectivamente, Figura 3B). Após a classificação dos pixels, o APP 2 executou operações morfológicas com o objetivo de definir os ROIs Stroma e Parenchyma(Figura 3B e Tabela 1). O desempenho do APP 2 na classificação de pixels e na geração dos respectivos ROIs é mostrado na Figura 3B. Além disso, o APP 2 quantifica a área do estroma e do parenchyma. Finalmente, mesmo que a segmentação seja feita usando a seção manchada de PSR, as regiões de stroma e parenchyma delineadas podem ser transferidas para qualquer imagem alinhada à imagem PSR.

Figura 3: Detecção/segmentação automatizada de tecidos e geração de respectivos ROIs. (A) A Imagem IIA foi utilizada para identificar os TAPs (imagem esquerda, barra de escala = 6.000 μm). Uma máscara verde brilhante foi atribuída aos TAPs usando app 1 (Tabela 1) gerando um ROI chamado Tissue (saída 1). Certo, inset mostra visualização ampliada demonstrando a precisão do APP 1 na detecção de TAPs. Barra de escala = 350 μm.(B) O tecido de ROI (saída 1) é segmentado em estroma e parenchyma usando o APP 2. A imagem à esquerda mostra uma visão do tecido roi segmentado em estromo roi (salmão) e parenchyma ROI (verde). Barra de escala = 4.500 μm. À direita, vistas ampliadas do inset para o ROI Tissue, a coloração psr original (imagem IIIC), e o estromo rois e parenchyma. Barra de escala = 250 μm. Clique aqui para ver uma versão maior desta figura.
Quantificação automatizada de COIs
Em seguida, passamos a identificar, localizar e quantificar COIs nos ROIs Stroma e Parenchyma. ApPs 3 a 8(Tabela 1) foram criados para localizar e contar os seguintes COIs: CD4+FoxP3+, CD8+, CD68+, MPO+, αSMA+, e células CD34+, respectivamente. O APP 3 foi projetado para localizar e contar células CD4+FoxP3+ (imagem IIIA, Figura 2C) como marcadores substitutos de células T regulatórias (Tregs). Este protocolo detecta a colocalização do sinal a partir do fator de transcrição nuclear FoxP3 (banda vermelha) e do dna de rotulagem de dna DAPI (banda azul). Dado que as células T recentemente ativadas regulam o FoxP3, para enriquecer para Tregs, estabelecemos limites para a pré-seleção apenas de células FoxP3+ brilhantes (FoxP3hi). Em seguida, de todas as células deoi DAPI+FoxP3 pré-selecionadas, apenas aquelas que estavam cercadas por sinais CD4 em forma de anel brilhante (banda verde) foram rotuladas e contadas como células FoxP3hiCD4+ (etiqueta rosa, Figura 4A). A densidade das células FoxP3hiCD4+ nos ROIs Stroma e Parenchyma foram determinadas como variáveis de saída quantitativa satisfatoida de APP 3 (Figura 4A).
Da mesma forma, as APPs 4 a 6 foram projetadas para a detecção de células CD8+, CD68+e MPO+. Esses APPs compartilham o mesmo desenho de linha de base para detectar e quantificar COIs. Especificamente, os COIs são identificados com base na intensidade do sinal do biomarcador específico da população celular e, em seguida, várias etapas morfológicas pós-processamento são executadas para delinear células individuais(Tabela 1). As células individuais ou COIs são rotuladas, contadas e suas coordenadas teciduais registradas. As APPs 4 a 6 também determinam a densidade dos COIs nos ROIs Stroma e Parenchyma (Figura 4B-D).
A qualidade da nossa coloração DAPI não foi boa o suficiente para integrar a segmentação de núcleos em APPs 3 a 6, por isso não podemos garantir que todos os objetos rotulados individualmente sejam células individuais. Por essa razão, expressamos a densidade das células em contagens de objetos rotulados/mm2 (Figura 4). No entanto, os agregados celulares foram separados com sucesso em células individuais nas etapas de pós-processamento incorporadas em APPs 3 a 6, e uma extensa inspeção visual mostrou que a maioria dos objetos rotulados correspondia a células únicas.
Para detectar a área αSMA+ e CD34+, desenvolvemos APPs 7 e 8, respectivamente (Tabela 1). Ambas as APPs detectam o sinal específico com base em limiares e determinam a porcentagem de área positiva nos ROIs Stroma e Parenchyma (Figura 4E-F).
Uma das possibilidades mais interessantes de gerar slides multiplex virtuais é a análise da expressão de colocalização. Geramos o APP 10 para detectar a colocalização entre αSMA e desmin, dois marcadores co-expressos por miofibroblastos no fígado. O APP 10 usa limiares para encontrar pixels positivos para αSMA, desmin e αSMA plus desmin(Tabela 1). Como variáveis de saída quantitativas, o APP 10 determina a área αSMA+, a área de desmin+ e a área de expressão colocalizada desses dois marcadores(Figura S3).

Figura 4: Identificação e quantificação de COIs no estroma e parenchyma. (A-F) Detecção e quantificação automatizadas de CoIs CD4+FoxP3+, CD8+, CD68+, MPO+, αSMA+, e CD34+ COIs nos ROIs Stroma e Parenchyma utilizando os protocolos 3, 4, 5, 6, 7 e 8, respectivamente(Tabela 1). Mostradas à esquerda as imagens originais, no meio as imagens processadas, e à direita as quantificações. Para figuras 4A-D, barra de escala = 40 μm. Para as figuras 4E e F, barra de escala = 350 μm. Clique aqui para ver uma versão maior desta figura.
Como alternativa para quantificar os COIs nos TCs Stroma e Parenchyma, determinamos a densidade de células imunes nos diferentes nódulos malignos nomeados 1 a 4(Figura 5A, He I). O ROI para cada nódulo foi delineado manualmente conforme indicado na Figura 5A. As assinaturas imunológicas teciduais distintas caracterizaram cada nódulo, revelando ainda mais a heterogeneidade intrínseca do TME.
Mapas térmicos de tecidos
Como mencionado acima, as APPs 3 a 8 armazenam as coordenadas teciduais de cada objeto rotulado individualmente. Este recurso permite a geração automatizada de mapas teciduais onde regiões de alta densidade de uma determinada população celular são exibidas como pontos quentes (vermelho), e regiões com densidade relativamente baixa como pontos frios (azul escuro). Os valores de densidade intermediária são atribuídos cores de acordo com a escala de cor mostrada na Figura 5. Mapas térmicos de tecidos foram gerados por APPs que dividiram as imagens em círculos de 50 μm de diâmetro e atribuíram uma cor de acordo com a densidade relativa de um dado COI dentro do círculo. Como mostrado na Figura 5B-G,os padrões de posicionamento e a distribuição de intensidade dos diferentes COIs no TME foram bastante variados. Além disso, ao nível dos nódulos individuais, o arranjo de diferentes populações na área tecidual foi único (Figura S2A-C). Para dar um exemplo do poder dessa técnica e visualizar a organização espacial de pontos quentes de diferentes populações no mesmo nódulo, os pontos quentes de tipos de células individuais foram extraídos manualmente e mapeados juntos no contorno do nódulo 2 (Figura S2, Figura De Figura E).

Figura 5: Mapas de calor de tecido de COIs no TME. (A)Coloração vermelha picrosirius mostrando localização dos nódulos 1, 2, 3 e 4. (B-G) Mapas de calor de tecido para COIs CD4+FoxP3+, CD8+, CD68+, MPO+, CD34+e αSMA+, respectivamente. Azul escuro indica densidade relativa baixa, e vermelho indica relativa alta densidade. Os valores de densidade intermediária são atribuídos cores de acordo com a escala de cor mostrada. (H e I) Quantificação de COIs nos nódulos 1, 2 e 3 + 4 organizados por tipo celular e por nódulo, respectivamente. Clique aqui para ver uma versão maior desta figura.

Figura Suplementar S1: Validação do alinhamento tecidual. (A) A coloração CD34 (em vermelho) feita na seção II (entrada 1) é usada para gerar uma máscara CD34 em verde (saída 1). A máscara verde (saída 1) é sobreposta na imagem H&E da seção serial alinhada I (entrada 2). A imagem de fusão mostra correspondência perfeita de estruturas vasculares. Barra de escala = 50 μm. (B) Imagem IIIA mostrando a fusão de DAPI, CD4 e FoxP3 (entrada 1) foi usada para gerar uma etiqueta para células CD4+FoxP3+ (saída 1 em magenta). O rótulo de saída 1 foi transferido para a imagem alinhada IIIB (entrada 2) e mostra correspondência perfeita entre os pares FoxP3/DAPI e CD4/CD3 na imagem de mesclagem. Barra de escala = 15 μm. Clique aqui para ver uma versão maior desta figura.

Figura suplementar S2: Visão ampliada dos mapas de calor do tecido. (A–C) Mapas de calor tecidual para células CD4+FoxP3+, CD8+, CD68+e MPO+ em nódulos 1-4. As barras de escala nos nódulos 1, 2 e 3 + 4 representam 1.500 μm, 700 μm e 500 μm, respectivamente. (D) Contorno do nódulo 2 com linha sólida preta. (E) Os pontos quentes para as células CD4+FoxP3+, CD8+, CD68+e MPO+ no nódulo 2 foram extraídos e mapeados juntos no contorno do nódulo 2 definido em D. Clique aqui para ver uma versão maior desta figura.

Figura Suplementar S3: Análise de Colocalização. (A) À esquerda e no meio estão imagens da etiqueta αSMA em verde e etiqueta de smin em vermelho, respectivamente. À direita está uma área dupla positiva αSMA/desmin em amarelo. (B) Quantificação da área αSMA+, área de desmin + e αSMA/desmin área dupla positiva. Barra de escala = 150 μm. Clique aqui para ver uma versão maior desta figura.
| App | Propósito | Classificação | Classificação | Etapas pós-processamento | Variáveis de saída |
| Método | Características |
| (valor de pixel) |
| 1 | Detecção de tecido | Limite | Canal DAPI (150) | o Rotular objetos com valores acima do limiar colocalizados para os 3 canais | o TECIDO ROI |
| Canal FITC/A488 (120) | o Fechar objeto positivo 5 pixels | o Área de Tecido |
| Canal TRITC/A568 (40) | o Criar tecido de ROI | |
| 2 | Segmentação de tecidos | Floresta de Decisão | Mediana RGB-R | o Preencher buracos | o ROI Stroma |
| Mediana RGB-G | o Criar ROI Stroma | o Área de Stroma |
| Mediana RGB-B | o Criar ROI Parenchyma | o ROI Parenchyma |
| Mediana do IHS-S | | o Parenchyma Área |
| H&E Eosin mediana | | |
| 3 | Para localizar e quantificar células CD4+ FoxP3+ | Limite | Canal DAPI (>600) | o Etiquetar objetos com colocalização de DAPI e Cy5/A647, cercados pelo sinal FITC/A488 | o Contagem e densidade de células CD4+FoxP3+ em ROIs Stroma e Parenchyma |
| Canal FITC/A488 poli suavização (>850) | o Limpar objetos menores que 7 μm2
| o Coordenadas de células CD4+FoxP3+ individuais |
| Canal Cy5/A647(>800) | | |
| 4 | Para localizar e quantificar células CD8+ | Limite | Canal DAPI (<1200) | o Limpar objetos positivos menores que 15 μm2
| o Contagem e densidade de células CD8+ em ROIs Stroma e Parenchyma |
| Canal Cy5/A647 mediana (>80) | o Fechar objetos positivos 2 pixels | o Coordenadas de células individuais |
| o Objetos separados | |
| 5 | Para localizar e quantificar células CD68+ | Limite | Canal FITC/A488 (>200) | o Limpar objetos positivos menores que 20 μm2
| o Contagem e densidade de células CD68+ em ROIs Stroma e Parenchyma |
| o Dilate objetos positivos 3 pixels | o Coordenadas de células CD68+ individuais |
| o Objetos separados | |
| 6 | Para localizar e quantificar células MPO+ | Limite | Canal DAPI (>400) | o Limpar objetos menores que 5 μm2
| o Contagem e densidade de células MPO+ em ROIs Stroma e Parenchyma. |
| Canal TRITC/A568 (900-4000) | o Dilado 3 pixels objetos positivos | o Coordenadas de células MPO+ individuais. |
| o Objetos separados | |
| 7 | Para localizar e quantificar a área αSMA+ | Limite | Canal TRITC/CF568 (>1050) | o Limpar objetos positivos menores que 25 μm2
| o Contagem e densidade da área αSMA+ em ROIs Stroma e Parenchyma |
| o Dilado 3 pixels objetos positivos | o Coordenadas de pixels αSMA+ |
| 8 | Para localizar e quantificar a área CD34+ | Limite | Canal DAPI (<5000) | o Limpar objetos positivos menores que 25 μm2
| o Contagem e densidade da área CD34+ em ROIs Stroma e Parenchyma |
| Canal Cy5/A647 mediana (>120) | o Dilado 3 pixels objetos positivos | o Coordenadas de pixels CD34+ |
| 9 | Crie mapas de calor de tecido para uma determinada população celular | Mapa de calor do objeto | Mapa de calor do objeto | | o Mapa de Calor |
| Raio de desenho 50 μm | --- |
| 10 | Quantificar a colocalização entre αSMA e Desmin | Limite | Canal TRITC (CF568) (>1050) | o Rotular objetos com valores de limite acima para TRITC (CF568) | o Quantificar expressão co-localizada de αSMA e Desmin |
| Canal Cy5 (A647) (>1000) | o Rotular objetos com valores de limite acima para Cy5 (A647) |
| o Rotular objetos com colocalização de valores limiares acima para TRITC (CF568) e Cy5 (A647) |
| o Limpar objetos positivos menores que 25 μm2
|
Tabela 1: Parâmetros gerais utilizados para o desenho de APPs empregados para análise de imagem. Os parâmetros especificados nesta tabela são ajustados às características únicas das imagens utilizadas nesta análise (por exemplo, fundo, artefatos, etc.) e podem não ser aplicáveis a outras imagens. Como as etapas de pós-processamento mencionadas foram definidas para as imagens específicas analisadas neste estudo, elas não são intencionalmente detalhadas. O usuário deve personalizar as APPs para as imagens a serem analisadas.
| Seção/Coloração | Anticorpo Primário | Anticorpo Secundário |
| Seção II/1st Coloração | Mouse IgG2a anti-humano αSMA Mouse IgG1 anti-humano CD34 Coelho anti-humano Cytokeratin 8/18 | Cabra anti-rato IgG2a CF568 Rato anti-rato IgG1 A647 Burro anti-coelho A488 |
| Seção II/2de Coloração | Coelho anti-humano Desmin CD68 anti-humano do rato | Burro anti-coelho A647 Burro anti-mouse DyLight 755 |
| Seção III/1st Coloração | CD anti-humano do rato Coelho anti-humano FoxP3 MPO anti-humano de cabra | Burro anti-mouse A488 Burro anti-coelho A647 Burro anti-cabra A568 |
| Seção III/2de Coloração | Coelho anti-humano CD3 CD anti-humano do rato | Burro anti-mouse DyLight 755 Burro anti-coelho A647 |
Tabela 2: Pares de anticorpos secundários primários para if