Genome-resolved functional characterisation of the broiler caecal microbiome following commercial gut health interventions

Resumo Técnico: Caracterização funcional em nível de genoma do microbioma cecal de frangos de corte após intervenções comerciais de saúde intestinal

Definição do Problema
O microbioma cecal de frangos de corte é uma comunidade metabolicamente ativa, crítica para a utilização de nutrientes e a homeostase do hospedeiro. Embora as Intervenções de Saúde Intestinal (GHI — Gut Health Interventions) — incluindo probióticos, prebióticos, óleos essenciais, vacinas contra coccidiose e dietas com redução de proteína bruta — sejam amplamente utilizadas para modular este ecossistema, o seu impacto é frequentemente avaliado através de mudanças taxonômicas. No entanto, as mudanças taxonômicas podem ser sutis ou inconsistentes em ecossistemas microbianos maduros. Além disso, as abordagens metagenômicas existentes, centradas em genes, não conseguem resolver como os traços funcionais, tais como resistência antimicrobiana (ARG), virulência (VAG) e adaptação ao estresse (SAG), estão distribuídos entre populações microbianas específicas. Há uma carência de dados de resolução genômica que caracterizem como as combinações comerciais de GHI moldam o repertório de genes funcionais do intestino de frangos de corte, particularmente no que diz respeito ao ciclagem biogeoquímica e à adaptação ecológica de táxons específicos a estas intervenções.

Metodologia
Este estudo empregou uma abordagem metagenômica de resolução genômica utilizando 307 Metagenome-Assembled Genomes (MAGs) de alta qualidade, recuperados de 118 amostras de conteúdo cecal através de dois ciclos de produção (C1 e C2). As amostras originaram-se de um ensaio comercial envolvendo seis grupos de tratamento (T1–T6) que receberam várias combinações de GHI, incluindo ionóforos, vacinas vivas atenuadas contra coccidiose, probióticos à base de Bacillus (A e B), óleos essenciais, prebióticos e dietas com redução de proteína bruta.

• Reconstrução Genômica: Os MAGs foram montados a partir de dados de metagenômica shotgun da Illumina, filtrados para ≥50% de completude e ≤10% de contaminação.
• Anotação Funcional: Os genes foram anotados utilizando o pipeline METABOLIC, AMRFinderPlus e ABRicate contra múltiplas bases de dados (NDARO, CARD, MEGARes, VFDB, etc.) para identificar ARGs, SAGs e VAGs. Foram aplicados limiares de ≥70% de identidade de nucleotídeos e ≥50% de cobertura.
• Análise Filogenética: Uma árvore filogenética foi construída utilizando 175 MAGs contendo genes de cópia única (SCGs) suficientes para avaliar a novidade e a distribuição taxonômica.
• Modelagem Estatística: Métricas de diversidade (alfa e beta) foram calculadas utilizando o pacote vegan. Associações entre componentes de GHI e abundância de genes funcionais foram modeladas usando Modelos de Variáveis Latentes de Função Linear Generalizada (GLLVM). A análise de abundância diferencial foi realizada com DESeq2, utilizando o grupo apenas com ionóforo (T1) como referência.
• Reconstrução Metabólica: Vias para o ciclo do carbono, nitrogênio e enxofre foram identificadas para avaliar a potencial produção de gases de efeito estufa (metano, óxido nitroso) e outras sínteses de metabólitos (amônio, sulfeto de hidrogênio).

Principais Resultados

• Panorama de Genes Funcionais: Um total de 208 genes funcionais putativos foram anotados. Genes associados à virulência (VAGs) foram os mais prevalentes (51,92%), seguidos por genes de resistência antimicrobiana (ARGs, 43,75%) e genes associados ao estresse (SAGs, 4,33%).
  • Distribuição Taxonômica: Firmicutes A carregou as classes de ARG mais diversas. A resistência a glicopeptídeos foi o ARG mais amplamente distribuído, encontrado em 65,08% dos MAGs portadores de ARG. Notavelmente, nenhum determinante de resistência a ionóforos foi detectado em qualquer MAG.
  • Concentração de Genes: Os genes funcionais não estavam distribuídos aleatoriamente; MAGs específicos, como Escherichia flexneri (bin.48), abrigavam um elevado número de ARGs, SAGs e VAGs, sugerindo que o risco microbiano está ligado ao repertório genético de membros chave da comunidade, e não apenas à abundância.
• Impacto das Combinações de GHI:
  • Redução da Diversidade: No Ciclo 1, a combinação de probióticos à base de Bacillus A, óleos essenciais e vacinação contra coccidiose (T5) reduziu significativamente a alfa diversidade do resistoma e do viruloma em comparação com outros tratamentos.
  • Abundância de Genes: O T5 foi associado a uma abundância significativamente menor de genes de resistência específicos a macrolídeos (mef(En2)) e lincosamidas (lnuA, lnu(AN2)), bem como resistência ao arsênico (arsB). Estes genes estavam restritos a Phocaeicola vulgatus (bin.96).
  • Intervenções Individuais: O Probiótico A foi unicamente associado à redução da abundância de SAG e VAG sem alterar significativamente a abundância de ARG, indicando um perfil funcionalmente neutro-a-benéfico. Em contraste, a vacinação contra coccidiose mostrou associações positivas amplas com a abundância de ARG.
• Ciclagem Biogeoquímica:
  • Gases de Efeito Estufa: Nenhum MAG possui a capacidade de metanogênese (produção de metano) ou produção/consumo de óxido nitroso. No entanto, 19 MAGs (principalmente Firmicutes A) codificam vias para o consumo de metano (metanotrofia).
  • Ciclagem de Nutrientes: Previu-se que os MAGs participam da produção de amônio via fixação de nitrogênio e amonificação de nitrito, e síntese de sulfeto de hidrogênio via redução de sulfito. Nenhum MAG foi previsto para consumir amônio ou sulfeto de hidrogênio.
  • Efeitos das Intervenções: Probiótico A, Probiótico B, ionóforos de fase de acabamento e prebióticos não alteraram as funções biogeoquímicas centrais (metanotrofia, síntese de amônio ou sulfeto de hidrogênio). A vacinação contra coccidiose foi associada a mudanças em todas as 14 etapas de reação biogeoquímica elegíveis.

Principais Contribuições

1. Linha de Base de Resolução Genômica: Este estudo fornece a primeira caracterização sistemática, em nível de genoma, de como as combinações comerciais de GHI moldam o panorama de genes funcionais do microbioma cecal de frangos de corte, indo além das mudanças taxonômicas para distribuições de genes específicos.
2. Diferenciação dos Efeitos das Intervenções: Demonstra que combinações específicas de GHI, particularmente aquelas contendo probióticos, podem reduzir a diversidade e a abundância de genes de resistência e virulência, enquanto outras intervenções (como a vacinação isolada) podem aumentar a abundância de genes funcionais.
3. Adaptação Ecológica vs. Pressão Seletiva: A distribuição ampla da resistência a glicopeptídeos e a ausência de resistência a ionóforos sugerem que o resistoma observado reflete a adaptação ecológica ao ambiente intestinal, em vez de pressão seletiva dos antimicrobianos específicos utilizados no ensaio.
4. Potencial Biogeoquímico: Estabelece que o microbioma cecal de frangos de corte tem capacidade intrínseca limitada para produção de gases de efeito estufa, mas possui potencial para o consumo de metano e ciclagem de nutrientes (síntese de amônio e sulfeto de hidrogênio), que permanece amplamente estável através da maioria das estratégias de GHI.

Significância e Alegações
Os autores alegam que esta investigação avança a base de evidências para a seleção racional de GHI na produção avícola. Ao demonstrar que a seleção estratégica de GHI — especificamente combinações envolvendo probióticos — pode modular o panorama de genes funcionais para reduzir a diversidade do resistoma e do viruloma, o estudo apoia iniciativas de One Health visando mitigar riscos microbianos. As descobertas sugerem que intervenções contendo probióticos oferecem uma abordagem favorável para manter um ecossistema intestinal estável sem interromper as funções biogeoquímicas centrais. O estudo enfatiza que estes resultados representam o potencial genômico e servem como um passo exploratório para compreender como as intervenções comerciais moldam o microbioma funcional ao nível da população, em vez de confirmar a resistência fenotípica ou a dinâmica de transmissão.