Este estudo foi revisado pelo Comitê de Ética do Primeiro Hospital Afiliado da Universidade de Medicina Chinesa de Anhui (Número da revisão ética: 2025AH-101, Data de aprovação ética: 2025-07-02).
Temas de pesquisa
Pacientes com DP e controles saudáveis do Primeiro Hospital Afiliado da Universidade de Medicina Chinesa de Anhui foram inscritos. Amostras frescas de fezes foram coletadas para análise de acompanhamento. Amostras fecais frescas eram coletadas, alicotadas em tubos EP estéreis, imediatamente congeladas em nitrogênio líquido e então armazenadas a −80 °C para preservação a longo prazo. Os critérios de inclusão dos pacientes com DP foram os seguintes: (1) diagnóstico clínico de DP; (2) idade >70 anos; (3) para o grupo Piribedil: iniciou a monoterapia com piribedil após o registro e não teve exposição prévia a este medicamento; para o grupo de DP: não receber tratamento com medicamentos anti-Parkinson; (4) não há outras doenças sistêmicas graves, como doenças cardíacas graves, doença hepática, doença renal ou câncer. Antes da inscrição, foi obtido consentimento informado por escrito de cada participante. Os critérios de exclusão foram os seguintes: (1) complicações graves da DP, como disfunção motora severa; (2) recebeu tratamento antibiótico recentemente; (3) presença de outras condições de saúde que possam afetar os resultados do estudo, como doença inflamatória intestinal, síndrome do intestino irritável, etc.; (4) alérgico ao piribedil. Os participantes saudáveis do controle (grupo Blank) eram indivíduos com idade pareada (idade >70 anos) sem histórico de doenças neurológicas ou gastrointestinais graves.
Todos os participantes foram obrigados a manter seus hábitos alimentares regulares antes da inscrição. Nenhum recebeu probióticos, prebióticos, suplementos de fibras alimentares ou outros suplementos nutricionais. Além disso, todos os participantes não haviam recebido qualquer tratamento antibiótico por pelo menos três meses antes da coleta da amostra. Nenhuma intervenção dietética padronizada foi realizada durante o estudo, mas todos os participantes foram instruídos a manter um estilo de vida e uma dieta estáveis para minimizar possíveis efeitos de confusão sobre a microbiota intestinal.
Agrupamento de pesquisa
Essa coorte do estudo foi dividida em três grupos: grupo DP (n = 10): amostras fecais coletadas de pacientes com doença de Parkinson que não estavam tomando medicação anti-Parkinson. Grupo Piribedil (n=10): amostras fecais coletadas de pacientes com doença de Parkinson após 21 dias consecutivos de tratamento com piribedil (50 mg, duas vezes ao dia). Grupo controle em branco (n=5): amostras fecais coletadas de indivíduos saudáveis.
Sequenciamento de RNA ribossomal 16S
O sequenciamento foi realizado em uma plataforma de sequenciamento, conforme relatado em outrosestudos 12,13. Resumidamente, DNA microbiano total foi isolado de amostras fecais usando o método do dodecilsulfato de sódio (SDS) 14,15,16. A integridade do DNA extraído foi verificada por eletroforese em gel de agarosa a 1%, e a concentração e pureza do DNA foram determinadas com um espectrofotômetro. A região hipervariável V3-V4 do gene bacteriano 16S rRNA foi amplificada a partir do DNA genômico usando primers degenerados 341F (5'-CCTAYGGGRBGCASCAG-3') e 806R (5'-GGACTACNNGGGTATCTAAT-3') com termociclador de Reação em Cadeia da Polimerase (PCR). A amplificação do gene bacteriano 16S rRNA foi realizada sob as seguintes condições térmicas: uma etapa inicial de pré-desnaturação a 94 °C por 2 minutos, seguida por 28 ciclos consistindo em desnaturação a 94 °C por 30 s, recozimento a 55 °C por 30 s e elongação a 72 °C por 30 s. Uma extensão final foi realizada a 72 °C por 10 minutos, e a reação foi então mantida em 4 °C. As misturas de PCR contêm 2× mistura mestra PCR de 25 μL, primer direto (10 μM), 2 μL, primer reverso (10 μM) 2 μL, DNA molde de 10 ng e, finalmente, ddH2O até 50 μL. As reações PCR foram realizadas em triplicada. As misturas de PCR contêm 2× mistura mestre de PCR de 25 μL, primer direto (10 μM) 2 μL, primer reverso (10 μM) 2 μL, DNA molde de 10 ng e, finalmente, ddH2O até 50 μL. As reações PCR foram realizadas em triplicada. O produto PCR foi extraído de gel de agarose a 2% e purificado usando um kit de extração em gel. O Kit de Preparação da Biblioteca Universal Plus DNA para MGI V2 foi usado para construir a biblioteca: (1) Ligação Splinter; (2) Utilizar triagem magnética de esferas para remover os segmentos autocontíguos da junta; (3) Enriquecimento de modelos de biblioteca por amplificação por PCR; (4) Esferas magnéticas para recuperar produtos de PCR; (5) Ciclização de produtos PCR. Após a validação da biblioteca, o sequenciamento era realizado.
Pontos de verificação de qualidade: Checkpoint 1 (Coleta de amostras): As amostras foram congeladas imediatamente após a coleta e armazenadas a −80°C. Qualquer amostra com sinais visíveis de descongelamento ou contaminação foi descartada. Ponto de controle 2 (extração de DNA): Pureza aceitável: A260/A280 = 1,8–2,0, A260/A230 ≥1,8; Integridade verificada por uma única faixa de alto peso molecular em gel de agarose a 1%. Checkpoint 3 (amplificação por PCR): Uma única faixa a ~460bp em gel de agarose a 2% indica amplificação bem-sucedida (primers 341F/806R, região V3-V4); Se ausente ou inespecífico, ajuste a temperatura de recozimento ou a concentração do molde. A concentração de recuperação ≥0,5 ng/μl foi necessária para a preparação da biblioteca. Checkpoint 4 (Preparação e sequenciamento de biblioteca): Realizado pela instalação de sequenciamento seguindo seus protocolos validados padrão; nenhum ponto de verificação adicional interno (por exemplo, dimensionamento do bioanalisador) foi aplicado.
Segurança e descarte de resíduos: Luvas e jalecos eram usados durante todo o manuseio das amostras. As superfícies de trabalho foram limpas com 70% de etanol. Todos os consumíveis contaminados (ponteiras, tubos, géis) foram autoclavados (121°C, 20 min) antes de serem descartados como resíduos bio-perigosos.
Bioinformática
As leituras de sequenciamento bruto do gene 16S rRNA foram primeiro desmultiplexadas, submetidas a filtragem de qualidade usando fastp (versão 0.20.0) e montadas via FLASH (versão 1.2.11) sob critérios padronizados de triagem. Resumidamente: (i) Sequências brutas de 300 pb foram reduzidas quando a média de qualidade caiu para abaixo de 20 dentro de uma janela deslizante de 50 pb. Leituras com comprimentos truncados menores que 50 bp ou contendo bases ambíguas foram removidas; (ii) A montagem da sequência era realizada apenas para sobreposições acima de 10 pb, com a razão máxima de desajuste da região sobreposta definida em 0,2. Leituras não pareadas que não puderam ser montadas foram eliminadas; (iii) Amostras individuais foram particionadas com base em sequências de código de barras e espoletas para ajuste subsequente da orientação. A correspondência de códigos de barras foi realizada com alinhamento exato, enquanto a correspondência de primers permitiu no máximo 2 desajustes de nucleotídeos. As leituras de rRNA 16S foram analisadas usando o pipeline Quantitative Insights into Microbial Ecology (QIIME) (versão 1.8.0)17. Sequências de tags foram adicionadas a leituras de extremidade pareada usando o FLASH (versão 1.2.11)18. O corte foi realizado usando UPARSE (versão 7.0.1090) com comprimento mínimo de 200 bp e qualidade média de 25, seguido por agrupamento com um corte de 97% para gerar unidades taxonômicas operacionais (OTUs)19. As leituras representativas da OTU foram classificadas taxonomicamente via RDP Classifier 2.2, com limite de confiança ≥0,6. As abundâncias de etiquetas e leituras foram determinadas por comparação usando o USEARCH global20. A identificação de sequências quiméricas foi baseada em comparações com o banco de dados Gold usando UCHIME (v4.2.40)21. A análise de coordenadas principais (PCoA) foi posteriormente realizada no conjunto de dados esparso baseado na dissimilaridade de Bray–Curtis, com o limiar de similaridade de sequência estabelecido em 97% para agrupamento. OTUs representativas e perfis taxonômicos correspondentes foram analisados mais detalhadamente usando o software R (versão 3.1.1). O algoritmo de tamanho do efeito da análise discriminante linear (LEfSe), baseado no tamanho do efeito da análise discriminante linear (LDA), foi aplicado para triagem de biomarcadores com abundância diferencial significativa no nível da OTU. Esse método enfatiza tanto a significância estatística quanto a relevância biológica, identificando primeiro características com diferenças consistentes entre os grupos (usando o teste de Kruskal-Wallis para comparação multigrupo) e depois estimando o tamanho do efeito de cada característica diferencialmente abundante usando LDA. Uma pontuação LDA >2,0 é definida como um aumento significativo daabundância 22.
Métodos Estatísticos
Análises estatísticas foram realizadas usando SPSS. Os dados normalmente distribuídos foram analisados usando a análise unidirecional da variância (ANOVA). Dados que não seguiram a distribuição normal foram analisados usando o teste de soma de ranks de Wilcoxon. Todos os dados quantitativos foram expressos como média ± erro padrão (SE) para parâmetros que seguem uma distribuição normal, enquanto dados não distribuídos normalmente foram apresentados como mediana com intervalo interquartil. Um valor-p < 0,05 foi considerado estatisticamente significativo.