Todos los procedimientos que involucraron a participantes humanos fueron revisados y aprobados por el Comité de Ética del Hospital Chongqing Chenjiaqiao (aprobación nº 20240131). Se obtuvo el consentimiento informado por escrito de todos los participantes antes de la recogida de la muestra. Todas las muestras de heces y la información de los participantes se manejaron utilizando códigos de estudio desidentificados, según se indica en la Tabla de Materiales.
Selección de participantes e inscripción
Los pacientes recién diagnosticados con DM T2 que visitaron el Hospital Chongqing Chenjiaqiao entre junio de 2023 y diciembre de 2023 fueron evaluados para comprobar su elegibilidad. Los participantes se inscribían en el grupo DM si tenían entre 20 y 65 años, tenían un nuevo diagnóstico de DM2, aceptaban proporcionar una muestra fresca de heces y completaban la evaluación clínica y los exámenes clínicos requeridos. Los voluntarios sanos reclutados durante el mismo periodo fueron incluidos como el grupo emparejado con no diabéticos (NM) y se les emparejó con el grupo DM lo más estrechamente posible en cuanto a edad, sexo y características demográficas generales. Se excluyó a los participantes si presentaban trastornos hematológicos, enfermedades del sistema nervioso central, enfermedad reumática activa, enfermedad autoinmune, infección gastrointestinal aguda o crónica, diarrea crónica, estreñimiento, úlcera gastrointestinal activa o en proceso de curación, enfermedad inflamatoria intestinal, síndrome del intestino irritable, tuberculosis intestinal, tumores gastrointestinales, otras malignidades, insuficiencia cardíaca grave, insuficiencia hepática grave, insuficiencia renal severa, otros metabólicos graves enfermedades, desnutrición, inmunodeficiencia, trastornos metabólicos congénitos, enfermedades psiquiátricas, uso de sedantes-hipnóticos o abuso de drogas. También se excluyó a los participantes si habían recibido antibióticos, agentes supresores de ácido, fármacos para la motilidad gastrointestinal, probióticos, glucocorticoides o inmunosupresores dentro del mes previo a la extracción de heces. Además, se excluyeron las personas que habían experimentado diarrea, se habían sometido a cirugía gastrointestinal o endoscopia, o reportado cambios bruscos en el entorno vital o hábitos alimenticios en el plazo de 1 mes previo a la recogida de heces. Se asignó un código único de estudio a cada participante elegible antes de la recogida de la muestra. De la cohorte elegible, se seleccionaron 12 participantes recién diagnosticados con DT2 y 12 NM para secuenciación posterior de ADNr 16S y análisis del microbioma basados en la disponibilidad de muestras y los requisitos de calidad de secuenciación.
Recogida y almacenamiento de muestras de heces
A cada participante se le proporcionó un recipiente estéril para la recogida de heces etiquetado únicamente con el código asignado para asegurar su desidentificación. Se indicó a los participantes que vaciaran la vejiga antes de defecar para minimizar la contaminación de orina y que recogieran heces frescas directamente en el recipiente estéril sin contacto con el agua del inodoro, orina, desinfectantes u otros posibles contaminantes. Inmediatamente después de la recogida, aproximadamente 1–2 g de heces se transfirieron a un tubo criogénico estéril utilizando una cuchara de muestreo desechable estéril. El tubo se sellaba herméticamente y se colocaba inmediatamente sobre hielo seco o en un recipiente de transferencia preenfriado. Todas las muestras se transportaron al laboratorio en un plazo de 2 horas tras la recogida. Al llegar al laboratorio, se verificó el código del estudio y se inspeccionó cada tubo en busca de fugas o contaminación visible. Se registraron los tiempos de recogida y almacenamiento de todas las muestras. Las muestras de heces se almacenaban posteriormente a −80 °C hasta la extracción de ADN genómico. Las muestras se consideraron aceptables para análisis aguas abajo solo si permanecían completamente congeladas durante el almacenamiento y el transporte y no mostraban evidencia de fugas de tubo ni contaminación externa.
Extracción de ADN genómico
Las muestras de heces se extraían de almacenamiento a −80 °C y se colocaban sobre hielo antes del procesamiento. Cada muestra se descongeló solo hasta que fue posible homogeneizar para minimizar la degradación asociada a ciclos repetidos de congelación-descongelación. Aproximadamente 200 mg de heces se transfirieron a un tubo estéril de microcentrífuga, y el buffer de lisis proporcionado en el kit de extracción de ADN de heces se añadió según las instrucciones del fabricante que figuran en la Tabla de Materiales. Las muestras se vórtizaron vigorosamente durante 30 segundos para lograr una homogeneización completa y luego se incubaron a temperatura ambiente durante 5–10 minutos para facilitar la lisis celular. Tras la lisis, las muestras se centrifugaron a 12.000 × g. durante 10 minutos a 4 °C. El sobrenadante se transfirió cuidadosamente a un nuevo tubo estéril de microcentrífuga sin alterar la cápsula. Posteriormente, el ADN genómico fue eluído utilizando el tampón de elución suministrado con el kit o agua libre de nucleasa según el protocolo del fabricante. El ADN extraído se almacenaba a −20 °C para uso a corto plazo o a −80 °C para su preservación a largo plazo. Las muestras de ADN consideradas aptas para análisis aguas abajo parecían claras y libres de partículas visibles.
Evaluación de la calidad del ADN
La calidad del ADN se evaluó mediante electroforesis en gel de agarosa y cuantificación basada en fluorescencia. Se preparó un gel de agarosa al 1% en un tampón de electroforesis, y las muestras de ADN extraídas, mezcladas con tampón de carga, se cargaron en los pozos de gel. Se realizó electroforesis hasta que las bandas de ADN se separaron adecuadamente, tras lo cual el gel fue examinado bajo iluminación ultravioleta o de luz azul. La integridad del ADN se consideraba aceptable cuando se observaba una banda genómica intacta sin degradación marcada ni difuminación excesiva. La concentración de ADN fue posteriormente cuantificada mediante un sistema de cuantificación de ADN basado en fluorescencia. Las muestras se diluyeron entonces hasta la concentración requerida para la amplificación de PCR posterior. Las muestras de ADN de alta calidad solían mostrar una banda distinta en electroforesis en gel de agarosa y estaban suficientemente concentradas para los procedimientos posteriores de amplificación.
Amplificación por PCR de la región 16S rDNA V3–V4
La región hipervariable V3–V4 del ADN bacteriano 16S se amplificaba utilizando cebadores con código de barras. Las secuencias de cebadores fueron las siguientes: cebador frontal 5′-ACTCCTACGGGGGCAG-3′ y cebador inverso 5′-GGACTACHVGGGTWTCTAAT-3′. Las reacciones PCR se prepararon sobre hielo en un volumen final de 20 μL que contenía 4 μL de tampón PCR, 2 μL de mezcla de nucleótidos (2,5 mmol/L), 0,8 μL de cebadores directos y inversos (5 μmol/L), 0,4 μL de ADN polimerasa de alta fidelidad, 10 ng de ADN plantilla y agua libre de nucleasa a volumen. Las mezclas de reacción se mezclaban suavemente por pipeteo y se centrifugaban brevemente para recoger el líquido en el fondo del tubo. La amplificación por PCR se realizó utilizando las siguientes condiciones de ciclo: desnaturalización inicial a 95 °C durante 5 minutos, seguida de 25 ciclos de desnaturalización a 95 °C durante 30 s, recocido a 55 °C durante 30 s y extensión a 72 °C durante 30 s, con un paso final de extensión a 72 °C durante 10 min. Para minimizar el sesgo de amplificación, cada muestra se amplificaba en reacciones triplicadas, y los productos de PCR resultantes de la misma muestra se agrupaban posteriormente en un solo tubo. La amplificación exitosa se confirmó mediante la presencia de una banda clara de amplicón del tamaño esperado en la electroforesis en gel de agarosa.
Verificación y purificación de productos por PCR
Los productos PCR fueron verificados mediante electroforesis en gel de agarosa al 2%. Los productos combinados de PCR de cada muestra se cargaron sobre el gel y electroforaron hasta que las bandas del amplicón objetivo quedaron claramente separadas. Las bandas correspondientes al tamaño esperado del amplicón se extirparon usando una hoja limpia y estéril y se purificaron con un kit de extracción en gel según las instrucciones del fabricante que figuran en la Tabla de Materiales. Los productos purificados de PCR se eluían posteriormente en agua sin elución o sin nucleasas. La concentración de cada producto purificado de PCR se cuantificó mediante un sistema de cuantificación de ADN basado en fluorescencia, y los valores de concentración se registraron para la preparación posterior de la biblioteca. La amplificación y purificación exitosas se confirmaron mediante una sola banda clara al tamaño esperado del amplicón en electroforesis en gel de agarosa.
Preparación y agrupación de bibliotecas
Las bibliotecas de secuenciación se prepararon utilizando un kit de preparación de biblioteca de ADN según las instrucciones del fabricante que figuran en la Tabla de Materiales. Las secuencias adaptadoras de secuenciación se ligaron a los amplicones purificados siguiendo el flujo de trabajo estándar de preparación de bibliotecas. Los productos ligados con adaptador se purificaban posteriormente mediante extracción en gel o un método de purificación basado en perlas, dependiendo del protocolo de preparación de la biblioteca seleccionado. Los productos purificados de la biblioteca se evaluaron utilizando electroforesis en gel de agarosa al 2% para confirmar la distribución adecuada de fragmentos. Las concentraciones en biblioteca se cuantificaron mediante un sistema de cuantificación de ADN basado en fluorescencia. En función de la concentración de ADN y el tamaño del fragmento, cada biblioteca se normalizó a la misma concentración molar para asegurar una profundidad de secuenciación equilibrada entre muestras. A continuación, se agruparon cantidades molares iguales de bibliotecas individuales en una proporción 1:1 para generar el pool final de la biblioteca de secuenciación. Antes de la secuenciación, la biblioteca agrupada se desnaturalizaba según el protocolo de la plataforma de secuenciación antes de la carga de la muestra. Las bibliotecas consideradas adecuadas para la secuenciación mostraron una distribución clara de fragmentos dentro del rango de tamaño esperado y una concentración suficiente para el análisis de secuenciación posterior.
Secuenciación de extremos emparejados
La biblioteca final agrupada se cargaba en una plataforma de secuenciación de extremos emparejados según los procedimientos operativos estándar del fabricante. La secuenciación de extremos emparejados se realizó utilizando una configuración de 2 × 250 pb, o 2 × 300 pb, para asegurar una profundidad de secuenciación adecuada y una caracterización exhaustiva de la diversidad microbiana en cada muestra. La calidad de la secuenciación se monitorizó durante toda la partida utilizando métricas de calidad generadas por la plataforma. Tras completar la secuenciación, se exportaron archivos FASTQ en bruto de extremo emparejado para cada muestra para análisis bioinformático posterior. Se consideraron aceptables los conjuntos de datos de secuenciación cuando demostraron suficiente profundidad de lectura por muestra, puntuaciones de calidad adecuadas en el Q30 y asignación exitosa de códigos de barras.
Procesamiento de secuencias en bruto
Los archivos FASTQ en bruto de extremo emparejado se importaban en una cadena de análisis bioinformático del microbioma para su procesamiento posterior. Las secuencias se desmultiplexaron en función de los códigos de barras asignados a cada muestra. Posteriormente se realizó filtrado de calidad para eliminar lecturas que contenían puntuaciones de baja calidad, bases ambiguas o artefactos de secuenciación. Las bases de baja calidad situadas en los extremos de las lecturas se recortaron en función de la distribución de puntuaciones de calidad para mejorar la fiabilidad general de la secuencia. Tras el control de calidad, las lecturas de extremos emparejados se fusionaron a través de regiones superpuestas para reconstruir secuencias de longitud completa. Solo se conservaron secuencias fusionadas de alta calidad para análisis posteriores del microbioma.
Generación de variantes de secuencia de amplicones (ASV) y anotación taxonómica
La reducción de ruido de secuencias se realizó utilizando un algoritmo validado de desruido, incluyendo DADA2 o Deblur, para corregir errores de secuenciación y mejorar la precisión de la secuencia. Las secuencias quiméricas se identificaron y eliminaron durante el proceso de eliminación de ruido. Posteriormente se generaron secuencias representativas de ASV y tablas de abundancia correspondientes para cada muestra. La asignación taxonómica de secuencias ASV se realizó utilizando un clasificador Naive Bayes, y la anotación taxonómica se realizó con la base de datos de referencia SILVA 138. Se exportaron tablas de abundancia taxonómica a nivel de filo, familia y género para análisis aguas abajo. El conjunto final de datos de ASV incluyó secuencias representativas, información de abundancia y anotaciones taxonómicas correspondientes para cada muestra.
Análisis de diversidad alfa y beta
Se calcularon índices de diversidad alfa, incluyendo los índices de Chao, ACE, Shannon, Simpson y Coverage, para evaluar la diversidad y riqueza microbiana dentro de la muestra. Los índices de diversidad alfa se expresaron como media ± desviación estándar. La normalidad de los datos se evaluó utilizando la prueba de Shapiro–Wilk. Las variables normalmente distribuidas se compararon entre los grupos DM y NM usando las pruebas t de Student, mientras que las variables no distribuidas normalmente se analizaron mediante pruebas de suma de rangos de Wilcoxon. Posteriormente, se calcularon distancias de diversidad beta utilizando una matriz de distancias adecuada para evaluar diferencias en la composición de comunidades microbianas entre muestras. Se realizó un análisis de coordenadas principales (PCoA) para visualizar diferencias en la estructura de la comunidad microbiana entre grupos. Además, se realizó un análisis de similitudes (ANOSIM) para determinar si las diferencias entre grupos superaban la variación intragrupal, y se reportaron los valores R y p correspondientes de ANOSIM. El análisis de diversidad alfa reflejó la diversidad microbiana dentro de la muestra, mientras que el análisis de diversidad beta evaluó diferencias en la composición de las comunidades microbianas entre grupos.
Análisis taxonómico diferencial
La composición de la comunidad microbiana se resumió a nivel de filo, familia y género. Se generaron gráficos de barras apiladas para visualizar la abundancia relativa de taxones dominantes en muestras individuales y grupos de estudio. Se construyeron además curvas Pan/Core y diagramas de Venn para comparar los ASVs compartidos y específicos de grupo entre los grupos DM y NM. Las diferencias en la abundancia de taxones entre grupos se evaluaron mediante pruebas de suma de rangos de Wilcoxon. Posteriormente se realizó un análisis de LEfSe para identificar taxones microbianos que contribuyeron más fuertemente a la discriminación entre grupos. El análisis taxonómico diferencial permitió identificar taxones microbianos enriquecidos tanto en el grupo DM como en el NM.
Predicción funcional
Se utilizó PICRUSt2 o una cadena de predicción funcional validada equivalente para inferir vías funcionales microbianas a partir de perfiles de secuenciación de ADNrDNA 16S. Los datos de abundancia de ASV se normalizaron según los requisitos de la cadena analítica seleccionada, y las funciones predichas se mapearon a anotaciones de vías KEGG o MetaCyc. Posteriormente, se compararon las abundancias predichas de las vías entre los grupos DM y NM. Se visualizaron rutas significativamente diferentes utilizando mapas de calor u otros enfoques gráficos apropiados. Las funciones predichas se interpretaron como potencial metabólico microbiano inferido en lugar de una abundancia metabolita medida directamente. El análisis de predicción funcional permitió identificar vías microbianas candidatas que diferían entre grupos, incluyendo vías asociadas al metabolismo de carbohidratos y aminoácidos.
Análisis de aleatorización mendeliana
Las estadísticas resumen del estudio de asociación genómica (GWAS) para la microbiota intestinal se obtuvieron de conjuntos de datos públicos10. Los datos de asociación significativa a nivel genómico para taxones microbianos se obtuvieron del Catálogo GWAS (https://www.ebi.ac.uk/gwas/) del NHGRI-EBI utilizando números de acceso que iban de GCST90032172 a GCST90032644. Se accedieron a datos metagenómicos adicionales de la cohorte FINRISK 2002 a través del Archivo Europeo de Genomas-Fenómenos (ID de investigación: EGAS00001005020). Las estadísticas resumen GWAS para la diabetes mellitus tipo 2 también se obtuvieron del Catálogo NHGRI-EBI GWAS (número de acceso: ebi-a-GCST006867). Se seleccionaron variantes genéticas asociadas a taxones microbianos como variables instrumentales según umbrales estadísticos predefinidos, y se excluyeron variantes en el desequilibrio de ligamento. Los conjuntos de datos de exposición y resultados se armonizaron para asegurar una orientación alelada consistente. Posteriormente, se realizó un análisis de aleatorización mendeliana utilizando el método ponderado por la varianza inversa (IVW) como enfoque analítico principal. Se realizaron análisis de sensibilidad utilizando los métodos de MR-Egger y mediana ponderada. La intensidad del instrumento se evaluó utilizando estadísticas F, mientras que la heterogeneidad entre variables instrumentales se evaluó mediante la prueba Q de Cochran. La pleiotropía horizontal se examinó mediante la prueba de intercepción MR-Egger. Se aplicó la corrección de Bonferroni para tener en cuenta las comparaciones múltiples. Los hallazgos de la aleatorización mendeliana se interpretaron como asociaciones predichas genéticamente en lugar de evidencia definitiva de causalidad.
Salida de datos y punto final
Los resultados analíticos finales incluyeron tablas de abundancia de ASV, gráficos de composición taxonómica, métricas de diversidad alfa, análisis de diversidad beta, resultados taxonómicos diferenciales, perfiles de vías funcionales predichos y estimaciones de aleatorización mendeliana. Se verificaron todos los identificadores de muestra en los conjuntos de datos de secuenciación para garantizar la coherencia con los códigos de estudio desidentificados correspondientes. Se archivaron archivos de secuenciación en bruto, tablas ASV procesadas, salidas estadísticas y archivos fuente de figuras para análisis posteriores y gestión de datos. El protocolo se consideró completo una vez que se generaron y verificaron con éxito datos de secuenciación de alta calidad, perfiles taxonómicos, métricas de diversidad, análisis de vías funcionales predichas y resultados de aleatorización mendeliana.