Domestication shapes the pig gut microbiome and immune traits from the scale of lineage to population
Data files
Sep 13, 2023 version files 23.44 MB
-
16S_analyses.Rmd
-
align_uniprot.txt
-
ASV_counts_decontam.tsv
-
ASV_taxonomy_decontam.tsv
-
card.txt
-
cutadapt.txt
-
dada2.txt
-
decontam.txt
-
derep.txt
-
dom_cat_metagenomic_metadata.csv
-
eggnog.txt
-
fastp.txt
-
fecal_metadata_Apr2022.csv
-
map_to_hostgenome.txt
-
megahit.txt
-
pig_metagenomics_clean.Rmd
-
PIG2.gds
-
plink.vcf
-
prodigal.txt
-
qiime_tree.txt
-
README.md
-
salmon.txt
-
SNP_metadata.csv
-
SNP.R
-
tables.txt
-
tree.nwk
-
vfdb.txt
Abstract
Animal ecology and evolution have long been known to shape host physiology, but more recently, the gut microbiome has been identified as a mediator between animal ecology and evolution and health. The gut microbiome has been shown to differ between wild and domestic animals, but the role of these differences for domestic animal evolution remains unknown. Gut microbiome responses to new animal genotypes and local environmental change during domestication may promote specific host phenotypes that are adaptive (or not) to the domestic environment. Because the gut microbiome supports host immune function, understanding the effects of animal ecology and evolution on the gut microbiome and immune phenotypes is critical. We investigated how domestication affects the gut microbiome and host immune state in multiple pig populations across five domestication contexts representing domestication status and current living conditions: free-ranging wild, captive wild, free-ranging domestic, captive domestic in research or industrial settings. We observed that domestication context explained much of the variation in gut microbiome composition, pathogen abundances, and immune markers, yet the main differences in the repertoire of metabolic genes found in the gut microbiome were between the wild and domestic genetic lineages. We also documented population-level effects within domestication contexts, demonstrating that fine scale environmental variation also shaped host and microbe features. Our findings highlight that understanding which gut microbiome and immune traits respond to host genetic lineage and/or scales of local ecology could inform targeted interventions that manipulate the gut microbiome to achieve beneficial health outcomes.
README: Domestication shapes the pig gut microbiome and immune traits from the scale of lineage to population.
https://doi.org/10.5061/dryad.9zw3r22mb
Description: We investigated how domestication affects the gut microbiome and host immune state in multiple pig populations across five domestication contexts representing domestication status and current living conditions: free-ranging wild, captive wild, free-ranging domestic, captive domestic in research or industrial settings. We observed that domestication context explained much of the variation in gut microbiome composition, pathogen abundances, and immune markers, yet the main differences in the repertoire of metabolic genes found in the gut microbiome were between the wild and domestic genetic lineages. We also documented population-level effects within domestication contexts, demonstrating that fine scale environmental variation also shaped host and microbe features. Our findings highlight that understanding which gut microbiome and immune traits respond to host genetic lineage and/or scales of local ecology could inform targeted interventions that manipulate the gut microbiome to achieve beneficial health outcomes. These files will allow you to replicate the bioinformatic and statistical analyses we used to assess the impact of domestication on the pig gut microbiome and immune state.
Description of the data and file structure
For 16S rRNA gene sequencing analyses, there is an R markdown file (16S_analyses.Rmd), a metadata file for each sample with covariates such as location where sample was collected and domestication group (fecal_metadata_Apr2022.xlsx), and files generated from read processing: the number and type of ASVs assigned to each sample (ASV_counts_decontam.tsv), a table that details the taxonomic lineage of each ASV from Kingdom to Species (ASV_taxonomy_decontam.tsv), and a phylogenetic tree (tree.nwk). Some samples in the metadata file will have N/A values for some covariates as we do not have data for that particular covariate. For cytokine measurements, some samples will be marked as "below_detection" if their values were below the detection limit.
For shotgun metagenomic analyses, there is an R markdown file (pig_metagenomics_clean.Rmd) and a metadata file (dom_cat_metagenomic_metadata.xlsx). The R markdown files will allow you to replicate any statistical analyses and visualizations using output files from the steps detailed in Code information below. The SNP.R file along with the plink.vcf, PIG2.gds, and SNP_metadata.xlsx files can be used to replicate all SNP analyses.
Sharing/Access information
All sequencing data are available in NCBI SRA (accession numbers PRJNA926635 for 16S and PRJNA926638 for metagenomics) with associated metadata files.
Code information
The following R files should be used in order to process raw 16S rRNA gene sequencing data:
- Cutadapt to remove primers (cutadapt.txt)
- Dada2 to learn errors, dereplicate, infer ASVs, merge amplicons, make sequence tables, and assign taxonomy (dada2.txt)
- Decontam to remove contaminants, mitochondria, and chloroplasts (decontam.txt)
- Tables to create final counts and taxonomy tables (tables.txt)
- Tree to create a phylogenetic tree (qiime_tree.txt)
The output will be an ASV decontaminated counts table, an ASV taxonomy table, and a phylogenetic tree.
The following files should be used in order to process raw shotgun metagenomic data:
- Quality filtering and trimming using Fastp (fastp.txt)
- Map sequences against reference genome (map_to_hostgenome.txt)
- Assemble to contigs using megahit (megahit.txt)
- Gene prediction with prodigal (prodigal.txt)
- Dereplication with cd-hit (derep.txt)
- Align to UniProt and taxonomic classification (align_uniprot.txt)
- Annotation with EggNog, CARD, and VFDB (eggnog.txt, card.txt, vfdb.txt)
- Salmon gene quantification (salmon.txt)
Outputs will be annotation files and quantification files.