Doce ratas machos Sprague-Dawley (SD) libres de patógenos específicos (SPF), de entre 6 y 7 semanas y con un peso de 220 g ± 10 g, fueron adquiridos a Beijing Vital River Laboratory Animal Technology Co., Ltd. (Número de licencia animal: SCXK (Jing) 2025-0008). Las ratas fueron alojadas en el Centro de Animales de Laboratorio del Hospital Dongzhimen, Universidad de Medicina China de Pekín, en cumplimiento con la norma nacional "Requisitos de Animales de Laboratorio para el Medio Ambiente y las Instalaciones de Vivienda" (GB 14925-2010, emitida por el Comité Técnico Nacional de Normalización de Animales de Laboratorio de China). Los procedimientos experimentales fueron aprobados por el Comité de Bienestar Animal y Ética del Hospital Dongzhimen, Universidad de Medicina China de Pekín (Número de aprobación: 24-54). Los reactivos y el equipo utilizado se enumeran en la Tabla de Materiales.
1. Establecimiento de modelos y recogida de muestras
NOTA: Antes del experimento, someta a todas las ratas a un periodo de aclimatación de 1 semana para minimizar los efectos inducidos por el estrés.
- Disolver lipopolisacáridos de Escherichia coli (LPS) en solución al 0,9% de suero para obtener una solución de LPS de 0,4 mg/mL.
- Dividir aleatoriamente 12 ratas Sprague-Dawley (SD) machos SPF en dos grupos: un grupo de control (n = 6) y un grupo modelo ARDS (n = 6).
- Administrar solución salina al 0,9% por inyección intratraqueal a una dosis de 0,5 mL/100g de peso corporal al grupo de control.
- Administrar la solución de LPS de 0,4 mg/mL mediante inyección intratraqueal al mismo volumen (0,5 mL/100g) al grupo modelo ARDS y, tras 16 horas, anestesiar a las ratas usando isoflurano (siguiendo protocolos institucionalmente aprobados).
- Recoge sangre a través de la aorta abdominal y luego centrifuga (4 °C, 1000 × g, 30 min) para obtener suero.
- Abre la cavidad torácica y extrae tejido pulmonar.
- Realizar lavado broncoalveolar mediante el muñón bronquial usando una jeringuilla para obtener Líquido de Lavado Broncoalveolar (BALF)24,25.
- Extrae el lóbulo del pulmón derecho y enjuágalo con PBS para un examen patológico posterior.
- Extraer el lóbulo pulmonar izquierdo y enjuágarlo con PBS para posteriores experimentos de secuenciación unicelular ybiología molecular 26,27.
2. Evaluación histopatológica del tejido pulmonar
- Sumergir el tejido del pulmón derecho extraído de las ratas en paraformaldehído al 4% para fijarlo durante 12 horas.
- Tras la fijación, deshidrata el tejido mediante una serie graduada de etanol (70%, 80%, 90%, 95%, 100% y 100%), permitiendo 2 horas en cada concentración.
- Trata el tejido deshidratado con xileno hasta que se vuelva transparente.
- Infiltra el tejido transparente con parafina fundida en un horno a 60 °C durante 2 horas.
- Coloca el bloque de tejido infiltrado en un molde de incrustación, vierte parafina fresca y deja que se enfríe y solidifique en un bloque de parafina.
- Secciona el bloque de parafina en rodajas de 4 μm de grosor usando un microtomo, haz flotar las secciones en un baño maria a 45 °C para extenderlas y luego recógelas con portaobjetos recubiertos de poli-L-lisina.
- Hornea las láminas en un horno a 60 °C para asegurar una adherencia firme de las secciones a las láminas.
- Desparafina las secciones por inmersión en xileno y rehidratalas mediante un gradiente descendente de etanol (100%, 95%, 90%, 80% y 70%).
- Enjuaga las secciones con PBS, realiza tinciones de hematoxilina y eosina (HE), y monta las cubiertas usando resina neutra.
- Evalúa la lesión pulmonar en secciones escaneadas por portaobjetos completos utilizando un sistema de puntuación estandarizado basado en la infiltración de neutrófilos, inflamación intersticial, edema y congestión.
NOTA: Cada parámetro fue puntuado de 0 a 4 según el grado de implicación patológica: 0 (ninguno); 1 (≤25%); 2 (26-50%); 3 (51-75%); 4 (>75%)28,29.
3. Ensayo inmunoenzimático ligado a enzimas
NOTA: Las muestras pueden mantenerse a -80 °C si no se puede realizar el análisis inmediatamente después de la adquisición.
- Disuelve los estándares y déjalos reposar a temperatura ambiente durante 15 minutos, luego realiza la dilución en serie según las instrucciones del fabricante.
- Prepara la solución de trabajo de anticuerpos biotinilados, la solución de trabajo conjugada enzimática y el tampón de lavado según las instrucciones del manual.
- Añadir 100 μL de cada muestra o 100 μL de enormes diluidos en serie a los pozos recubiertos de anticuerpos, e incubar a 37 °C durante 90 minutos.
- Lava la placa 4 veces con una arandela de placas, añade 100 μL de solución activa de anticuerpos biotinilados y incuba a 37 °C durante 30 minutos.
- Lava la placa 4 veces, añade 100 μL de solución de trabajo conjugado enzimático e incuba a 37 °C durante 30 minutos.
- Añadir 100 μL de solución de sustrato por pozo, incubar protegido de la luz a 37 °C durante 10 minutos.
- Añade 100 μL de solución de parada por pozo, mezcla suavemente y mide la densidad óptica (OD) a 450 nm usando un lector de microplacas.
- Genera una curva estándar representando los valores de OD frente a concentraciones estándar usando el software diseñado.
- Determinar las concentraciones de la muestra basándose en los valores de OD medidos y la curva estándar.
4. Preparación de suspensiones de celda única para secuenciación
NOTA: Tras obtener tejido pulmonar, debe colocarse en hielo y procesarse rápidamente para mantener la viabilidad celular.
- Lava el tejido pulmonar 3 veces con PBS pre-enfriado a 4 °C.
- Añade 2 mL de PBS preenfriado y pícalos aún más en fragmentos más pequeños usando tijeras estériles hasta que midan aproximadamente 0,5mm 3 pulgadas.
- Utiliza una pipeta Pasteur para transferir el tejido picado junto con PBS a un tubo centrífugo de 15 mL.
- Añadir 1 mL adicional de PBS preenfriado a 4 °C en la placa de cultivo, enjuagar y volver a suspender los fragmentos de tejido restantes, y transfiere la suspensión al mismo tubo.
- Ajusta el volumen a 6,5 mL con PBS preenfriado a 4 °C, digiera fragmentos de tejido en una solución enzimática que contenga colagénasa I (100 U/mL) y Dispase II (1 U/mL), e incube con suaves agitaciones a temperatura ambiente durante 30 minutos.
- Cuando los fragmentos de tejido se vuelven translúcidos y no quedan piezas sólidas distintas, filtrar el tejido digerido a través de coladores celulares de 70 μm y 40 μm, luego centrifugar (4° C, 500 × g, 5 min) y descartar el sobrenadante.
- Inmediatamente resuspende el pellet en 3 mL de tampón de lisis de glóbulos rojos (RBC), mezcla suavemente pipeteando y lisa sobre hielo durante 10 minutos.
- Cuando el sobrenadante aparece de un amarillo pálido o casi incoloro, se añaden 7 mL de PBS para terminar la reacción, centrifugar (4 °C, 500 × g, 5 minutos) y descartar el sobrenadante.
- Añade 3 mL de PBS para resuspender las células y determinar la concentración y viabilidad celular usando tinción con azul tripán.
NOTA: Si la viabilidad celular es inferior al 70%, realiza un procedimiento de extracción de células muertas.
- Lava las células dos veces con un tampón de lavado que contiene albúmina sérica bovina (BSA) al 0,04% preparada en DPBS.
- Suspende las células para contar y ajusta la concentración a 1.000 celdas/μL para la carga y detección posteriores.
5. Control de calidad para datos de secuenciación de una sola célula
NOTA: Las bibliotecas de ARN unicelular se prepararon a partir de suspensiones celulares utilizando un método de construcción de bibliotecas unicelulares basado en gotas. Las bibliotecas se secuenciaban en una plataforma de alto rendimiento utilizando la química de secuenciación basada en nanobolas de ADN. Se cargaron aproximadamente 21.000 celdas por biblioteca, con un objetivo de una tasa de recuperación de celdas del ~60%. La secuenciación generó un total de 5,0 ×10 8 lecturas por biblioteca, con una profundidad media de 30.000-40.000 lecturas por celda. Los datos de secuenciación en bruto se procesaron utilizando la correspondiente cadena de análisis de ARN unicelular para producir la matriz de expresión génica de cada muestra.
- Construir la biblioteca Oligo mediante PCR y etiquetado por código de barras, asegurando una concentración de >10 ng/μL y un tamaño máximo de 180 pb ±10 pb.
- Construir la biblioteca de ADNc mediante fragmentación, reparación de extremos, ligadura de adaptadores y PCR, asegurando una concentración >10 ng/μL y un tamaño máximo de 350–550 pb.
- Genera un objeto Seurat para cada muestra a partir de la matriz de expresión génica usando la función Read10X en el paquete Seurat (v4.3.0) en el software R (v4.2.0).
- Eliminar las células de baja calidad filtrando aquellas con proporción de genes mitocondriales >10% o proporción de genes de hemoglobina >5%.
- Normaliza los datos usando el método "LogNormalize" con un factor de escala de 10.000.
- Identificar características altamente variables utilizando el método "vst" y conservar los 2.000 genes superiores para análisis posteriores.
- Realizar la puntuación del ciclo celular utilizando los genes marcadores del ciclo celular incorporados de Seurat y regresar los efectos del ciclo celular (G2M. Score y S.Score) durante la escalada de datos.
- Detectar y eliminar dobletes potenciales usando DoubletFinder (v2.0.3) con pN = 0,25, pK = 0,09 y una tasa esperada de doblete de 0,06, y excluirlos del procesamiento posterior.
- Correcta los efectos por lotes entre muestras usando el paquete Harmony (v0.1.1) con theta = 2 y max.iter.harmony = 20.
- Añade columnas de información para nombres de muestras individuales y detalles de agrupación del grupo de modelos ARDS y grupo de control a los metadatos del objeto Seurat.
6. Reducción de dimensionalidad, agrupamiento y anotación de tipos de celda
- Realizar análisis de componentes principales (PCA) en genes altamente variables utilizando la función RunPCA en Seurat.
- Determinar el número óptimo de componentes principales mediante la visualización ElbowPlot.
- Identifica los principales clústeres celulares aplicando el algoritmo de Leiden mediante la función FindClusters y FindNeighbors de Seurat con una resolución de 0,1.
- Realizar la reducción de dimensionalidad de UMAP utilizando los componentes principales seleccionados.
- Anota los tipos celulares identificando genes marcadores con la función FindAllMarkers.
- Visualiza la expresión génica de marcadores entre tipos celulares usando la función DotPlot.
- Añade una nueva columna de información de los resultados de las anotaciones de celdas a los metadatos y visualiza la anotación por tipo de celda con Umap.
- Calcular la proporción de cada grupo experimental dentro de cada tipo celular identificado usando funciones dplyr para contar celdas y calcular porcentajes.
- Genera gráficos de barras apiladas con ggplot2 para ilustrar la composición de grupos por tipo de celda, incluyendo etiquetas porcentuales y comparaciones estadísticas entre grupos.
7. Análisis de subpoblaciones
- Extrae la subpoblación de la célula objetivo usando la función de subconjunto y crea un objeto Seurat correspondiente para el análisis a nivel de subpoblación aguas abajo.
- Renormalizar y reescalar los datos de la subpoblación mientras se regresan los genes mitocondriales, genes de hemoglobina, recuentos de UMI y puntuaciones del ciclo celular.
- Realizar un análisis secundario de PCA y aplicar la corrección por lotes de Harmony específicamente para la subpoblación.
- Realizar reagrupamiento utilizando funciones FindNeighbors y FindClusters para identificar subtipos distintos dentro de subpoblaciones celulares objetivo.
- Generar nuevas incrustaciones UMAP y visualizar los grupos de subpoblaciones con anotaciones de tipo celular siguiendo el mismo enfoque descrito anteriormente.
8. Análisis diferencial de expresión génica
- Identificar genes expresados diferencialmente (DEGs) entre el grupo modelo ARDS y el grupo de control para cada subgrupo utilizando el método estadístico MAST.
- Incorporar regresión de variables latentes para controlar los efectos de profundidad de secuenciación durante las pruebas de expresión diferencial.
- Filtrar los DEGs basándose en umbrales de cambio de pliegue logarítmico predefinidos y valores p ajustados usando la función FindMarkers.
- Genera gráficos volcánicos con la función ggplot para cada subcúmulo para visualizar patrones diferenciales de expresión génica con etiquetado de significación.
- Exportar los resultados de la DEG y crear tablas resumidas completas para facilitar el enriquecimiento funcional y otros análisis posteriores.
9. Análisis de enriquecimiento
- Mapea símbolos genéticos a IDs ENTREZ usando la organización. Rn.eg.db base de datos de anotaciones para compatibilidad con el análisis de vías.
- Combinar los DEGs de todos los subclústeres y filtrar genes significativos basándose en los umbrales de cambio de pliegue logarítmico y valores p ajustados (prueba de suma de rangos de Wilcoxon, logfc.threshold = 0,25, p_val < 0,05).
- Realizar análisis de enriquecimiento por Ontología Génica (GO) utilizando la función GO enriquecida con categorías de ontologías de Proceso Biológico, Función Molecular, Componentes Celulares y método de corrección de Benjamini-Hochberg.
- Realizar análisis de enriquecimiento de vías de la Enciclopedia de Genes y Genomas de Kyoto (KEGG) utilizando la función de enriquecer KEGG con la base de datos de organismos de ratas y aplicar umbrales de significación estadística.
- Genera gráficos de puntos para los resultados de enriquecimiento GO y KEGG usando la función de dotplot, mostrando los 15 términos más enriquecidos significativamente con proporciones génicas y valores p ajustados.
- Exporta los resultados del enriquecimiento como archivos CSV y guarda gráficos de visualización para una interpretación completa del análisis de vías.