banner

Notizia

Aug 11, 2023

Inter

The ISME Journal (2023) Citare questo articolo

Dettagli sulle metriche

L'analisi multi-omica è un potente strumento per il rilevamento e lo studio delle interazioni tra regni, come quelle tra membri batterici e archeologici di complesse comunità microbiche produttrici di biogas. Nel presente studio, i microbiomi di tre digestori di biogas su scala industriale, ciascuno alimentato con substrati diversi, sono stati analizzati utilizzando un quadro metagenomico genomico-centrico guidato dall’apprendimento automatico integrato con dati metatrascrittomici. Questi dati ci hanno permesso di chiarire la relazione tra abbondanti comunità metanogeniche centrali e i loro partner batterici sintrofici. In totale, abbiamo rilevato 297 genomi assemblati con metagenoma (nrMAG) di alta qualità e non ridondanti. Inoltre, i profili genetici assemblati di 16 S rRNA di questi nrMAG hanno mostrato che il phylum Firmicutes possedeva il numero di copie più alto, mentre i rappresentanti del dominio archaeal avevano il numero più basso. Ulteriori indagini sulle tre comunità microbiche anaerobiche hanno mostrato alterazioni caratteristiche nel tempo, ma sono rimaste specifiche per ciascun impianto di biogas su scala industriale. L'abbondanza relativa di vari microrganismi rivelata dai dati del metagenoma era indipendente dai corrispondenti dati sull'attività del metatrascrittoma. Gli Archaea hanno mostrato un'attività considerevolmente più elevata di quanto ci si aspettasse dalla loro abbondanza. Abbiamo rilevato 51 nrMAG presenti in tutti e tre i microbiomi delle piante di biogas con abbondanze diverse. Il microbioma centrale era correlato ai principali parametri chimici della fermentazione e nessun parametro individuale è emerso come un fattore predominante nella composizione della comunità. Vari meccanismi di trasferimento di H2/elettroni tra specie sono stati assegnati ai metanogeni idrogenotrofi negli impianti di biogas che funzionavano con biomassa agricola e acque reflue. L'analisi dei dati del metatrascrittoma ha rivelato che le vie della metanogenesi erano le più attive tra tutte le principali vie metaboliche.

Nei sistemi di digestione anaerobica (AD) ingegnerizzati, la decomposizione della materia organica e la produzione di biogas si basano su un efficiente riciclo dei nutrienti [1, 2]. Questo processo coinvolge una comunità microbica diversificata, in cui ciascun microbo ha un ruolo specifico. La composizione e la funzione della comunità microbica coinvolta negli stadi dell'AD svolgono un ruolo importante nell'efficienza del processo complessivo, che è influenzato anche da numerosi fattori, come le interazioni microbo-microbo, la composizione del substrato, i parametri fisico-chimici e le condizioni operative [3,4,5,6].

Svelare le interazioni microbiche e i loro meccanismi sottostanti è un compito complesso poiché molti microbi non possono sopravvivere senza partner microbici specifici [7]. Sono stati sviluppati metodi innovativi per l’isolamento e la coltivazione microbica, alcuni dei quali si sono rivelati efficaci nel portare nuovi microrganismi in coltura [8]. Tuttavia, sono necessari ulteriori sviluppi e avanzamenti di queste tecnologie per sfruttare appieno il potenziale della diversità microbica [9, 10]. Poiché esistono forti interazioni sintrofiche tra i microbi nei consorzi microbici di AD, è necessaria una ricerca approfondita per comprenderne la diversità, il ruolo metabolico e la distribuzione. Questi studi sono stati inizialmente ostacolati dalle limitazioni inerenti ai metodi microbiologici coltura-dipendenti, che richiedono l'isolamento dei microrganismi e possono essere difficili quando le relazioni sintrofiche sono onnipresenti [11,12,13]. Gli strumenti di sequenziamento e bioinformatica ad alto rendimento consentono l'analisi di massa del materiale genomico e quindi forniscono informazioni sulla tassonomia e sulle funzioni di intere comunità microbiche [14,15,16].

L’identificazione e la caratterizzazione del genoma dei microbi comunemente presenti nei microbiomi di AD possono rivelare percorsi cruciali e caratteristiche ecologiche coinvolte nella catena alimentare microbica [17, 18]. Studi precedenti che utilizzavano il sequenziamento dell'amplicone del gene 16 S rRNA hanno mostrato che i microrganismi centrali creano una comunità stabile in grado di resistere a varie perturbazioni (ad esempio, cambiamenti dei parametri AD) [19,20,21]. Tuttavia, il sequenziamento degli ampliconi e gli approcci metagenomici basati sulla lettura potrebbero non essere in grado di identificare specie sconosciute mediante l’analisi del microbioma centrale a causa della loro dipendenza dai database di riferimento [22, 23]. Inoltre, le comunità produttrici di biogas rappresentano un contingente ampio e diversificato di microrganismi non caratterizzati che sono stati precedentemente descritti come “materia oscura microbica” [19, 24, 25]. Per affrontare questo problema, è possibile utilizzare algoritmi bioinformatici sempre più sofisticati per ricostruire i genomi di singole specie (o MAG: genomi assemblati con metagenoma) da queste comunità complesse [26,27,28].

5% contamination by CheckM2 analysis (Fig. 2A and B). Furthermore, 24% of nrMAGs (n = 70) were identified at a quality level that was higher than their representatives in the Genome Taxonomy Database (Supplimentary Table 3). These observations confirm the effectiveness of Semibin as a binning procedure for complex anaerobic biogas-producing communities [35]. In addition, the larger number of unique nrMAGs binned by Semibin may also lead to better mapping of metatranscriptome data./p>90%), low contamination (<5%) and have not been found in any available databases (Archaea: MAG_165_3; Bacteria: MAG_114_1, _18_1, 29_1, 77_1, 87_1 and 95_2). Purple symbols represent the core nrMAGs. The outer ring represents bacteria that are significantly more prevalent in the particular BP (using lefser, p < 0.05)./p>90%; n = 107) were present in the Biogas Microbiome database (ANI cutoff: 95%). It is worth noting that seven high-quality nrMAGs, which accounted for 7% of the total high-quality nrMAGs (n = 107) and 2% of all nrMAGs (n = 297), could not be associated with a nearest representative in either database. These nrMAGs were deemed novel because they did not meet the following criteria: species-level identification with 95% reference radius for GTDB and ≥95% ANI for Biogas Microbiome (Fig. 4 and Suppl. Table 3). Three high-quality nrMAGs were found to contain 16 S rRNA gene sequences as well (Supplimentary Table 3). In addition, the seven putative novel nrMAGs demonstrated high abundance and activity, accounting for 3% of the total count per million (CPM) and 5% of the total transcripts per million (TPM), respectively, in the examined microbiome. These nrMAGs may be members of the hypothesised microbial dark matter that can now be released into the realm of known participants in AD communities./p> 0.5). Based on this observation, eight clusters showing characteristic correlations with AD chemical parameters were plotted (Fig. 6B). The eight clusters can be divided into two groups, with microorganisms from clusters I–V correlating positively with the main influencing parameters and microorganisms from clusters VI–VIII correlating negatively with these parameters. We also found that top core nrMAGs belonging to the phyla Firmicutes, Spirochaetota, and Methanobacteriota were positively correlated with TAN, VOAs, and TIC, while members of the phylum Bacteroidota were also correlated with many of the measured parameters. Finally, the C/N ratio showed a more pronounced impact on Bacteroidota than Firmicutes among the top core nrMAGs (Fig. 6B)./p> 0.7 based on network. The blue lines represent positive correlations (Pearson's rho > 0.7), red lines represent negative correlations (Pearson's rho < −0.7). A phylum marked with an asterisk is not present among the correlating microorganisms. B Pearson correlations between core nrMAGs and measured AD chemical parameters. Asterisks represent the significant correlations (considered statistically significant as follows: p < 0.05 (*), p < 0.001 (**), p < 0.0001 (***), ns (no significant difference)). Using the cladogram shown on the left and data for TAN (total ammonia nitrogen), VOAs (volatile organic acids), TIC (total inorganic carbon), TS (total solids) and C/N (carbon to nitrogen ratio), eight clusters were distinguished. These are shown here separated by colour as follows: Cluster I: grey; Cluster II: brown; Cluster III: yellow; Cluster IV: green; Cluster V: purple; Cluster VI: blue; Cluster VII: orange; and Cluster VIII: dark green. Coloured boxes represent the specific phyla to which the nrMAGs belong to. Red dots represent the top hydrolysing nrMAGs (n = 10)./p> −0.7). Based on genomic data, these archaea belong to the class of methanogens capable of maintaining hydrogenotrophic methanogenesis using formate as an electron donor. Typically, formate is oxidised to CO2 by formate dehydrogenase, upon which it is further reduced to methane [85]. Here, interspecies formate transfer (IFT) was maintained by microbial partners, namely SR-FBR-E99 sp002497965 and LD21 sp012519515 (from the phylum Bacteroidota; MAG_72_3 and MAG_394_2), which were positively correlated with MAG_26_2 and MAG_103_2 (Fig. 6A). These MAGs have been previously detected in biogas digesters and were described for their versatile hydrolysing ability. Aside from carbohydrate utilisation, they are protein- and amino-acid degrading microorganisms capable of producing formate. They were also consistently detected alongside formate-utilising methanogens [86]./p> 0.6). However, there were exceptions involving microorganisms belonging to Cluster VII./p> 2.0, p < 0.05). The heat map depicts the average gene activity of each methanogen in the various BPs. The blanks represent genes that were not significantly different in the provided methanogen, or the indicated gene is absent in the given nrMAG./p> 2.0, p ≤ 0.05) (Supplimentary Table 5). KBP and SZBP were shown to contain fewer differentially expressed genes (DEGs) than these compared to MWBP (Fig. 8B and Supplimentary Fig. 2). A comprehensive investigation of DEGs found 215 genes involved in methanogenesis. The activities of alpha, beta, and gamma subunits of methyl coenzyme M reductase showed the most substantial changes (log2 FC > 10, p < 0.0001). Additionally, the expression of genes implicated in hydrogenotrophic methanogenesis, such as the ion-sulphur subunit of F420-non-reducing hydrogenase and methylenetetrahydromethanopterin dehydrogenase, displayed substantial differences (log2 FC > 5, p < 0.001). Acetyl-CoA synthetase demonstrated variances in expression across genes involved in acetotrophic methanogenesis (log2 FC > 5, p < 0.001) [99]. The majority of the different gene activities in KBP were related to Methanothrix A harundinaceae D (MAG_5_1), while Methanothrix_sp016706325 (MAG_97_2) in MWBP and Methanobacterium sp012838205, Methanoculleus sp12797575, and Methanosarcina species (MAG_669_3, _57_1, and _165_3) in SZBP (Fig. 8C)./p>90% cont. <5%) nrMAGs are deposited in the NCBI SRA under accession numbers from SAMN32989584 to SAMN32989690. The main data generated or analysed for this study are included in this published article and its supplementary information files. Workflows and R scripts are available from the first author on reasonable request./p>

CONDIVIDERE