Doze ratos machos específicos livres de patógenos (FPS) Sprague-Dawley (SD), com idades entre 6 e 7 semanas e pesando 220 g ± 10 g, foram adquiridos da Beijing Vital River Laboratory Animal Technology Co., Ltd. (Número da Licença do Animal: SCXK (Jing) 2025-0008). Os ratos foram alojados no Centro de Animais de Laboratório do Hospital Dongzhimen, Universidade de Medicina Chinesa de Pequim, em conformidade com a norma nacional "Requisitos de Animais de Laboratório para o Ambiente e Instalações de Habitação" (GB 14925-2010, emitida pelo Comitê Técnico Nacional de Padronização de Animais de Laboratório da China). Os procedimentos experimentais foram aprovados pelo Comitê de Bem-Estar Animal e Ética do Hospital Dongzhimen, Universidade de Medicina Chinesa de Pequim (Número de Aprovação: 24-54). Os reagentes e os equipamentos utilizados estão listados na Tabela de Materiais.
1. Estabelecimento de modelos e coleta de amostras
NOTA: Antes do experimento, submeta todos os ratos a um período de aclimatação de 1 semana para minimizar os efeitos induzidos pelo estresse.
- Dissolva lipopolissacarídeo de Escherichia coli (LPS) em soro fisiológico a 0,9% para obter uma solução de LPS de 0,4 mg/mL.
- Divida aleatoriamente 12 ratos Sprague-Dawley (SD) machos do SPF em dois grupos: um grupo controle (n = 6) e um grupo modelo ARDS (n = 6).
- Administrar soro fisiológico a 0,9 % por injeção intratraqueal, em uma dose de 0,5 mL/100g de peso corporal, ao grupo controle.
- Administrar a solução de LPS de 0,4 mg/mL via injeção intratraqueal no mesmo volume (0,5 mL/100g) no grupo modelo ARDS e, após 16 horas, anestesiar os ratos usando isoflurano (seguindo protocolos institucionalmente aprovados).
- Colete sangue pela aorta abdominal, depois centrifuge (4 °C, 1000 × g, 30 min) para obter soro.
- Abra a cavidade torácica e extraia o tecido pulmonar.
- Realize lavagem broncoalveolar através do coto brônquico usando uma seringa para obter Fluido de Lavagem Broncoalveolar (BALF)24,25.
- Colha o lobo do pulmão direito e enxágue-o com PBS para exame patológico posterior.
- Colhe o lobo do pulmão esquerdo e enxágue com PBS para subsequentes experimentos de sequenciamento unicelular ebiologia molecular 26,27.
2. Avaliação histopatológica do tecido pulmonar
- Imerga o tecido do pulmão direito colhido de ratos em paraformaldeído a 4% para fixação por 12 horas.
- Após a fixação, desidrate o tecido por meio de uma série graduada de etanol (70%, 80%, 90%, 95%, 100% e 100%), permitindo 2 horas em cada concentração.
- Trate o tecido desidratado com xileno até que fique transparente.
- Infiltre o tecido transparente com parafina derretida em um forno a 60 °C por 2 horas.
- Coloque o bloco de tecido infiltrado em um molde de embustionamento, despeje parafina fresca e deixe que esfrie e solidifique em um bloco de parafina.
- Divida o bloco de parafina em fatias de 4 μm de espessura usando um micrótomo, flutue as seções em um banho-maria a 45 °C para espalhá-las e depois as recolha com lâminas revestidas de poli-L-lisina.
- Asse os lâminas em um forno a 60 °C para garantir uma adesão firme das seções aos lâminas.
- Desparafinize as seções por imersão em xileno e reidrate-as através de um gradiente descendente de etanol (100%, 95%, 90%, 80% e 70%).
- Enxágue as seções com PBS, realize tingimento de hematoxilina e eosina (HE) e monte as coberturas usando resina neutra.
- Avalie a lesão pulmonar em seções escaneadas por lâmina inteira usando um sistema padronizado de pontuação baseado na infiltração de neutrófilos, inflamação intersticial, edema e congestão.
NOTA: Cada parâmetro foi pontuado de 0 a 4 com base no grau de envolvimento patológico: 0 (nenhum); 1 (≤25%); 2 (26-50%); 3 (51-75%); 4 (>75%)28,29.
3. Ensaio imunossorvente ligado a enzimas
NOTA: As amostras podem ser mantidas em -80 °C se a análise não puder ser realizada imediatamente após a aquisição.
- Dissolva os padrões e deixe-os em temperatura ambiente por 15 minutos, depois realize a diluição serial conforme as instruções do fabricante.
- Prepare a solução de trabalho com anticorpos biotinilados, a solução de trabalho conjugada enzimática e o tampão de lavagem conforme instruído no manual.
- Adicione 100 μL de cada amostra ou 100 μL de padrões serialmente diluídos aos poços revestidos por anticorpos e incube a 37 °C por 90 minutos.
- Lave a placa 4 vezes usando uma arruela de placa, adicione 100 μL de solução de trabalho de anticorpos biotinilados e incube a 37 °C por 30 minutos.
- Lave a placa 4 vezes, adicione 100 μL de solução de trabalho conjugada enzimática e incube a 37 °C por 30 minutos.
- Adicione 100 μL de solução de substrato por poço, incube protegido da luz a 37 °C por 10 minutos.
- Adicione 100 μL de solução de parada por poço, misture suavemente e meça a densidade óptica (OD) a 450 nm usando um leitor de microplaca.
- Gere uma curva padrão plotando os valores de OD em relação às concentrações padrão usando o software projetado.
- Determine as concentrações amostrais com base nos valores de OD medidos e na curva padrão.
4. Preparação de suspensões de célula única para sequenciamento
NOTA: Após a obtenção do tecido pulmonar, ele deve ser colocado no gelo e processado rapidamente para manter a viabilidade celular.
- Lave o tecido pulmonar 3 vezes com PBS pré-resfriado a 4 °C.
- Adicione 2 mL de PBS pré-resfriado e depois triture em fragmentos menores usando tesoura estéreis, até aproximadamente 0,5mm 3 polegadas.
- Use uma pipeta Pasteur para transferir o tecido picado junto com PBS para um tubo centrífugo de 15 mL.
- Adicione mais 1 mL de PBS pré-resfriado a 4 °C na placa de cultura, enxágue e ressuspenda os fragmentos de tecido restantes, e transfira a suspensão para o mesmo tubo.
- Ajuste o volume para 6,5 mL com PBS pré-resfriado a 4 °C, digera fragmentos de tecido em solução enzimática contendo Colagênase I (100 U/mL) e Dispase II (1 U/mL), e incube com sacudidas suaves em temperatura ambiente por 30 minutos.
- Quando os fragmentos de tecido se tornam translúcidos e não restam pedaços sólidos distintos, filtre-se o tecido digerido através de peneiras celulares de 70 μm e 40 μm, depois centrifuge (4° C, 500 × g, 5 min) e descarte o sobrenadante.
- Imediatamente resuspenda o pellet em 3 mL de tampão de lise de glóbulos vermelhos (RBC), misture suavemente pipetando e lise no gelo por 10 minutos.
- Quando o sobrenadante aparecer amarelo claro ou quase incolor, adicione 7 mL de PBS para terminar a reação, centrifuge (4 °C, 500 × g, 5 min) e descarte o sobrenadante.
- Adicione 3 mL de PBS para ressuspender as células e determinar a concentração e viabilidade celular usando coloração com azul de tripán.
NOTA: Se a viabilidade celular estiver abaixo de 70%, realize um procedimento de remoção de células mortas.
- Lave as células duas vezes com um tampão contendo 0,04% de albumina sérica bovina (BSA) preparada em DPBS.
- Suspenda as células para contagem e ajuste a concentração para 1.000 células/μL para carregamento e detecção subsequentes.
5. Controle de qualidade para dados de sequenciamento de célula única
NOTA: Bibliotecas de RNA de célula única foram preparadas a partir de suspensões celulares usando um método de construção de biblioteca de célula única baseado em gotas. As bibliotecas eram sequenciadas em uma plataforma de alta capacidade usando química de sequenciamento baseada em nanobolas de DNA. Aproximadamente 21.000 células foram carregadas por biblioteca, visando uma taxa de recuperação de células de ~60%. O sequenciamento gerou um total de 5,0 × 108 leituras por biblioteca, com uma profundidade média de 30.000 a 40.000 leituras por célula. Os dados brutos do sequenciamento foram processados usando o pipeline correspondente de análise de RNA de célula única para produzir a matriz de expressão gênica para cada amostra.
- Construa a biblioteca Oligo por PCR e rotulagem de códigos de barras, garantindo concentração >10 ng/μL e tamanho máximo de 180 pb ±10 pb.
- Construa a biblioteca de cDNA por fragmentação, reparo final, ligadura de adaptadores e PCR, garantindo concentração >10 ng/μL e tamanho máximo de 350–550 pb.
- Gerar um objeto Seurat para cada amostra a partir da matriz de expressão gênica usando a função Read10X no pacote Seurat (v4.3.0) no software R (v4.2.0).
- Remova células de baixa qualidade filtrando aquelas com proporção do gene mitocondrial >10% ou proporção do gene da hemoglobina >5%.
- Normalize os dados usando o método "LogNormalize" com um fator de escala de 10.000.
- Identifique características altamente variáveis usando o método "vst" e mantenha os 2.000 primeiros genes para análise a jusante.
- Realize a pontuação do ciclo celular usando os genes marcadores do ciclo celular incorporados ao Seurat e regresse os efeitos do ciclo celular (G2M. Score e S.Score) durante a escalonação de dados.
- Detectar e remover potenciais doublets usando o DoubletFinder (v2.0.3) com pN = 0,25, pK = 0,09 e uma taxa esperada de doublet de 0,06, e exclua esses de processamento posterior.
- Efeitos em lote corretos entre amostras usando o pacote Harmony (v0.1.1) com theta = 2 e max.iter.harmony = 20.
- Adicione colunas de informações para nomes de amostras individuais e detalhes de agrupamento do grupo de modelos ARDS e grupo de controle aos metadados no objeto Seurat.
6. Redução de dimensionalidade, agrupamento e anotação de tipo de célula
- Realize análise de componentes principais (PCA) em genes altamente variáveis usando a função RunPCA em Seurat.
- Determine o número ótimo de componentes principais por meio da visualização ElbowPlot.
- Identifique os principais clusters de células aplicando o algoritmo de Leiden via as funções FindClusters e FindNeighbors de Seurat com resolução de 0,1.
- Realize a redução de dimensionalidade do UMAP usando os componentes principais selecionados.
- Anote tipos celulares identificando genes marcadores com a função FindAllMarkers.
- Visualize a expressão gênica dos marcadores entre tipos celulares usando a função DotPlot.
- Adicione uma nova coluna de informações com os resultados das anotações de células aos metadados e visualize a anotação do tipo de célula com o Umap.
- Calcule a proporção de cada grupo experimental dentro de cada tipo celular identificado usando funções dplyr para contar células e calcular percentuais.
- Gerar gráficos de barras empilhadas com ggplot2 para ilustrar a composição de grupos por tipo de célula, incluindo rótulos percentuais e comparações estatísticas entre grupos.
7. Análise de subpopulação
- Extraia a subpopulação de células-alvo usando a função de subconjunto e crie um objeto Seurat correspondente para análise subpopulacional a jusante.
- Renormalize e reescale os dados da subpopulação enquanto regrede genes mitocondriais, genes de hemoglobina, contagens de UMI e escores do ciclo celular.
- Realize uma análise secundária de PCA e aplique correção em lote de Harmony especificamente para a subpopulação.
- Realizar reagrupamento usando funções FindNeighbors e FindClusters para identificar subtipos distintos dentro das subpopulações celulares-alvo.
- Gerar novos embeddings UMAP e visualizar subagrupamentos de subpopulações com anotações de tipo celular seguindo a mesma abordagem descrita anteriormente.
8. Análise diferencial de expressão gênica
- Identifique genes diferencialmente expressos (DEGs) entre o grupo modelo ARDS e o grupo controle para cada subcluster usando o método estatístico MAST.
- Incorpore regressão de variáveis latentes para controlar os efeitos de profundidade do sequenciamento durante testes de expressão diferencial.
- Filtre DEGs com base em limiares pré-definidos de mudança logaritarítmica e valores p ajustados usando a função FindMarkers.
- Gerar gráficos de vulcão com a função ggplot para cada subaglomerado para visualizar padrões diferenciais de expressão gênica com rotulagem de significância.
- Exportar os resultados do DEG e criar tabelas resumidas abrangentes para facilitar o enriquecimento funcional e outras análises a jusante.
9. Análise de enriquecimento
- Mapear símbolos genéticos para IDs ENTREZ usando a org. Rn.eg.db banco de dados de anotações para compatibilidade na análise de vias.
- Combine DEGs de todos os subclusters e filtre genes significativos com base nos limiaritários de mudança logaritarítmica e nos valores p ajustados (teste de soma de rank de Wilcoxon, logfc.threshold = 0,25, p_val < 0,05).
- Realize análise de enriquecimento por Ontologia Gênica (GO) usando a função de enriquecimento GO com categorias de ontologia de Processo Biológico, Função Molecular, Componentes Celulares e método de correção de Benjamini-Hochberg.
- Realizar a análise de enriquecimento de vias da Enciclopédia de Genes e Genomas de Kyoto (KEGG) usando a função de enriquecimento KEGG com o banco de dados de organismos de ratos e aplicar limiares de significância estatística.
- Gerar dot plots para os resultados de enriquecimento GO e KEGG usando a função dotplot, exibindo os 15 termos significativamente enriquecidos com razões gênicas e valores p ajustados.
- Exporte os resultados do enriquecimento como arquivos CSV e salve gráficos de visualização para uma interpretação abrangente da análise de caminhos.