To infer the functional potential of 16S rRNA gene data among AH, DU, and DL, the program Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) was used in QIIME2-2021.11 [82 (link)]. Only SCTLD-susceptible corals were evaluated and only ASVs that were present in at least 100 samples were selected. The picrust2 full-pipeline was used with the hidden state set to “mp” and the placement tool to place sequences into a tree set to “epa-ng.” The outputs were predicted metagenomes for Kyoto Encyclopedia of Genes and Genomes (KEGG [83 (link)]) orthologs and MetaCyc pathway [84 (link)] abundances. To assess the differential abundance of these outputs among disease states, the R package Maaslin2 was utilized [85 ]. For both KEGG and MetaCyc tests, data were log-transformed, a random effect was set to coral species, and the data were subsequently analyzed with a linear model. In the KEGG assessment, the minimum abundance = 0.05 and the minimum prevalence = 0.1. There were no minimums set for the MetaCyc test due to the lower number of pathways found in MetaCyc. The top 10 predicted pathways were selected based on values with the lowest Padj and effect sizes < −0.5 and >0.5. The top pathways were manually annotated on KEGG and MetaCyc websites.
Free full text: Click here