O estudo foi conduzido de acordo com a Declaração de Helsinque, e o protocolo foi aprovado pelo Comitê de Ética do Segundo Hospital Afiliado da Universidade Médica de Nanjing (número de aprovação: 2021-KY-100-01). O consentimento informado por escrito foi obtido de todos os participantes.
Este estudo observacional retrospectivo de centro único foi conduzido no Departamento de Neurologia, Hospital Afiliado à Segunda Universidade Médica de Nanjing. De janeiro de 2022 a dezembro de 2023, foram incluídos um total de 32 pacientes com DP primária que atenderam aos Critérios Diagnósticos para a Doença de Parkinson (edição de 2020). No final, seis pacientes que atenderam a critérios rigorosos de seleção e completaram todos os procedimentos do estudo foram incluídos na análise subsequente. Os detalhes do processo de triagem, exclusão e acompanhamento são apresentados na Figura Suplementar 1.
Coleta de amostras fecais e análise da comunidade microbiana
O microbiota fecal foi obtido de voluntários saudáveis. O processo de triagem dos doadores foi rigorosamente conduzido de acordo com o Consenso Chinês de Especialistas sobre a Aplicação Clínica do Transplante de Microbiota Fecal. Critérios de inclusão dos doadores: Idade entre 18 e 40 anos, sem histórico de doenças crônicas; ausência de histórico de doenças gastrointestinais, distúrbios metabólicos, doenças autoimunes ou transtornos psiquiátricos; ausência de uso de antibióticos, probióticos, imunossupressores ou inibidores da bomba de prótons nos últimos 3 meses; ausência de viagens a áreas epidêmicas, transfusões sanguíneas ou histórico cirúrgico nos últimos 6 meses. A triagem sanguínea e fecal para exclusão de doenças infecciosas e patógenos inclui: vírus da imunodeficiência humana (HIV), vírus da hepatite B (HBV), vírus da hepatite C (HCV), Treponema pallidum, Clostridium difficile, Salmonella, Shigella, Campylobacter, patogénico Escherichia coli, vírus entéricos e parasitas. O volume de uma única amostra de fezes válida deve ser de pelo menos 50 g. Utilizando um sistema automatizado de purificação de microbiota fecal (consulte Tabela de Materiais), foram adicionados 500 mL de soro fisiológico estéril a cada 100 g de fezes frescas, seguidos por três ciclos de centrifugação e lavagem padronizados (700 × g, 3 min), e imediatamente inoculados após o preparo de uma suspensão com solução salina estéril em uma proporção volumétrica de 1:2. Os indivíduos receberam uma única injeção de 150 mL da suspensão contendo aproximadamente 5 × 10^13 bactérias viáveis13.
Todos os receptores receberam TFM por colonoscopia, com o instrumento avançado até o ceco e o cólon ascendente; a suspensão foi infundida lentamente e de forma uniforme através do canal de biópsia para promover a colonização preferencial do cólon direito. Após a infusão, os pacientes permaneceram na posição decúbito lateral direito por pelo menos 2 h. Antes do procedimento, foi seguida uma dieta de baixo resíduo durante 3 dias, seguida por uma dieta líquida completa nas 24 h anteriores ao procedimento, com jejum de sólidos por 6 h e de líquidos por 2 h antes da infusão. Na noite anterior ao transplante, foi realizada uma limpeza completa do intestino com solução eletrolítica de polietilenoglicol até obtenção de efluente claro e livre de fezes. As restrições medicamentosas incluíram a suspensão de antibióticos sistêmicos por pelo menos 72 h antes da TFM e a sua evitação durante 1 semana após o procedimento. Inibidores da bomba de prótons (IBPs) são descontinuados entre 24 e 48 h antes do procedimento, e laxantes e probióticos são suspensos no dia do transplante.
Durante o procedimento, a frequência cardíaca, a pressão arterial e a saturação de oxigênio foram monitoradas continuamente, e quaisquer sintomas, como distensão abdominal, dor ou náusea, foram observados; a infusão foi interrompida imediatamente caso ocorresse dor abdominal grave ou hipotensão. Durante os primeiros 7 dias após o transplante, foi realizado acompanhamento diário para registrar a frequência das evacuações intestinais, a consistência das fezes e todos os eventos adversos, incluindo dor abdominal, diarreia, febre e hematocázia. O acompanhamento de longo prazo estendeu-se até 5 meses, com coleta programada de amostras fecais para análise da microbiota intestinal e avaliações clínicas utilizando a Escala Unificada de Avaliação da Doença de Parkinson (UPDRS) e o Índice de Constipação de Wexner, juntamente com o registro de alterações nos sintomas e do uso de medicações concomitantes, a fim de avaliar de forma abrangente tanto a eficácia quanto a segurança.
As amostras fecais foram coletadas em seis momentos: antes do TFM e aos 1, 2, 3, 4 e 5 meses após o TFM. Todos os participantes receberam instruções padronizadas para a coleta das amostras de fezes. Aproximadamente 10–20 g da porção mediana das fezes foram coletados utilizando pinças estéreis e colocados em recipientes estéreis para coleta. As amostras foram transportadas para o laboratório em até 2–3 h após a coleta e fracionadas em tubos estéreis de microcentrífuga de 1,5 mL (ver Tabela de Materiais) por pesquisadores treinados em condições assépticas. Todas as alíquotas foram imediatamente armazenadas a −80 °C até a análise posterior.
Após a coleta das amostras fecais, todos os espécimes foram transportados em gelo seco para um fornecedor comercial de sequenciamento (consulte a Tabela de Materiais) para o sequenciamento do gene do RNA ribossômico 16S. A preparação das bibliotecas, o sequenciamento e o controle inicial de qualidade foram realizados de acordo com o protocolo padrão do fornecedor. O DNA genômico foi extraído de cada amostra fecal, e as regiões hipervariáveis V3–V4 do gene bacteriano do RNA ribossômico 16S foram amplificadas por reação em cadeia da polimerase (PCR) utilizando os iniciadores 341F e 805R. O sequenciamento foi realizado em uma plataforma de sequenciamento de leitura curta de bancada (consulte a Tabela de Materiais) utilizando química de leitura pareada 2 × 300 pb para caracterizar os perfis microbianos intestinais dos pacientes antes do TFM e em diferentes pontos temporais de acompanhamento pós-TFM. Após o sequenciamento, as leituras brutas foram processadas para obter sequências de alta qualidade para análises microbianas posteriores. A classificação taxonômica foi realizada utilizando um pacote comercial de atribuição taxonômica 16S e sua base de dados associada de identificação microbiana (consulte a Tabela de Materiais). Os perfis taxonômicos resultantes foram utilizados para análises subsequentes de diversidade microbiana, composição comunitária, abundância diferencial e função microbiana predita.
Avaliações clínicas e rótulos analíticos
Os pacientes com DP foram estratificados em dois grupos de gravidade de acordo com a duração da doença: um grupo de alta gravidade (>6 anos, designado como “_H”) e um grupo de baixa gravidade (≤6 anos, designado como “_L”). As amostras fecais coletadas antes do TFM foram definidas como amostras basais e rotuladas como FB_H ou FB_L de acordo com o grupo de gravidade correspondente. As amostras coletadas aos 1, 2, 3, 4 e 5 meses após o TFM foram rotuladas como FP1M_H, FP2M_H, FP3M_H, FP4M_H e FP5M_H para o grupo de alta gravidade, e como FP1M_L, FP2M_L, FP3M_L, FP4M_L e FP5M_L para o grupo de baixa gravidade. Os parâmetros clínicos avaliados incluíram o escore de constipação de Wexner, Escala de Sintomas Não Motores (NMSS), Questionário de Sintomas Não Motores (NMSQ), Inventário Neuropsiquiátrico (NPI), estágio de Hoehn–Yahr, Escala de Atividades de Vida Diária (ADL), Índice de Qualidade do Sono de Pittsburgh (PSQI), a escala original (UPDRS I–VI) e o escore total da UPDRS. Uma combinação de métodos estatísticos descritivos e inferenciais foi utilizada para analisar as alterações nos indicadores clínicos dos pacientes com DP no grupo de alta gravidade (grupo _H) e no grupo de baixa gravidade (grupo _L) em seis pontos temporais. Para explorar preliminarmente as diferenças entre os grupos em todos os pontos temporais, as distribuições dos indicadores entre os grupos _H e _L foram comparadas em cada ponto temporal utilizando o teste de soma de postos de Wilcoxon, e aplicou-se um modelo linear misto (LMM) mais rigoroso com interceptos aleatórios específicos por paciente. Devido ao tamanho limitado da amostra, os resultados do modelo de efeitos mistos são fornecidos apenas como referências exploratórias. O desfecho primário deste estudo exploratório foi a alteração na diversidade e composição microbiana intestinal do ponto basal até os pontos temporais pós-transplante de microbiota fecal, sendo as alterações clínicas avaliadas como desfechos secundários.
Agrupamento de OTUs, análise de diversidade e abundância taxonômica
Unidades taxonômicas operacionais (OTUs) foram utilizadas para resumir a composição da comunidade microbiana, e sequências com ≥97% de similaridade foram agrupadas na mesma OTU. A profundidade de sequenciamento foi avaliada por meio de curvas de rarefação baseadas no índice de Shannon, cujos platôs indicam, em geral, cobertura suficiente para análises de diversidade subsequentes.
A diversidade alfa foi avaliada utilizando os índices Chao1 e estimador de cobertura baseado na abundância (ACE) para riqueza microbiana, e os índices de Shannon e Simpson para diversidade geral e equitatividade. A diversidade beta foi avaliada utilizando distâncias UniFrac ponderadas e não ponderadas, seguidas por análise de coordenadas principais (PCoA) para visualizar as diferenças na composição da comunidade microbiana entre grupos. As diferenças entre grupos na diversidade alfa e beta foram avaliadas usando testes não paramétricos e Análise de Variância Multivariada Permutacional (PERMANOVA), respectivamente, conforme descrito na seção de análise estatística.
A composição taxonômica foi resumida nos níveis de filo, classe, ordem, família e gênero com base em dados de sequenciamento do gene do RNA ribossômico 16S. Utilizando um pacote de análise de ecologia comunitária para R (veja Tabela de Materiais), identifique e visualize os 10 táxons mais abundantes em cada nível taxonômico. Os perfis taxonômicos foram exibidos usando gráficos de barras, e táxons selecionados foram comparados entre grupos utilizando gráficos de caixa.
Análise da abundância diferencial de táxons microbianos
Os táxons diferencialmente abundantes entre grupos pré-definidos foram identificados utilizando a análise linear do tamanho do efeito (LEfSe). O LEfSe integra o teste de Kruskal–Wallis, o teste de Wilcoxon e a análise discriminante linear (LDA) para identificar táxons com significância estatística e valor discriminatório. Considerou-se haver diferença na abundância para táxons com P < 0,05 e escore LDA (log10) > 3. As abundâncias relativas dos táxons discriminatórios foram visualizadas utilizando gráficos de barras taxonômicos.
Previsão funcional e análise de correlação clínica
Os perfis funcionais previstos da microbiota intestinal foram inferidos a partir de dados de sequenciamento do gene 16S rRNA utilizando o Tax4Fun2 e anotados de acordo com a base de dados Kyoto Encyclopedia of Genes and Genomes (KEGG). As vias previstas diferencialmente enriquecidas foram identificadas usando o LEfSe, com P < 0,05 e um escore LDA (log10) > 3 como limiares. Como esses perfis funcionais foram inferidos em vez de diretamente medidos por abordagens metagenômicas ou metabolômicas, foram interpretados como diferenças funcionais previstas. Foi realizada análise de correlação de postos de Spearman para explorar associações entre características microbianas e indicadores clínicos.
Análise estatística
As análises estatísticas foram realizadas utilizando o software R (consulte a Tabela de Materiais). As variáveis clínicas foram resumidas como medianas e intervalos (mínimo–máximo). Os índices de diversidade alfa foram comparados entre grupos utilizando o teste de soma de postos de Wilcoxon, enquanto a diversidade beta foi avaliada por meio de PERMANOVA com 999 permutações. Os táxons e vias funcionais preditos com abundância diferencial foram identificados utilizando o LEfSe, com P < 0,05 e escore LDA (log10) > 3 como limiares de significância. Para a análise longitudinal de abundância diferencial, foi aplicado um modelo de efeitos mistos com interceptos aleatórios específicos por paciente, seguido de comparações pareadas pós-hoc com correção da taxa de falsa descoberta de Benjamini-Hochberg. A análise de correlação de postos de Spearman foi utilizada para avaliar as associações entre características microbianas intestinais e indicadores clínicos, com aplicação da correção da taxa de falsa descoberta de Benjamini–Hochberg aos valores de P de correlação. Um valor de P bicaudal < 0,05 foi considerado estatisticamente significativo. A correção para múltiplas comparações foi aplicada às análises de correlação e de abundância taxonômica. Dado o tamanho reduzido da amostra e o delineamento com amostragem longitudinal repetida, as análises estatísticas foram consideradas exploratórias, e os resultados foram interpretados com cautela.