Next Article in Journal
Severe COVID-19 and Sepsis: Immune Pathogenesis and Laboratory Markers
Next Article in Special Issue
Interactions between Anaerobic Fungi and Methanogens in the Rumen and Their Biotechnological Potential in Biogas Production from Lignocellulosic Materials
Previous Article in Journal
Automation of Infectious Focus Assay for Determination of Filovirus Titers and Direct Comparison to Plaque and TCID50 Assays
Previous Article in Special Issue
Cut-Lengths of Perennial Ryegrass Leaf-Blades Influences In Vitro Fermentation by the Anaerobic Fungus Neocallimastix frontalis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of a High-Grain Diet on the Rumen Microbiome of Goats with a Special Focus on Anaerobic Fungi

by
Katerina O. Fliegerova
1,*,†,
Sabine M. Podmirseg
2,†,
Julia Vinzelj
2,†,
Diego J. Grilli
3,†,
Simona Kvasnová
1,
Dagmar Schierová
1,
Hana Sechovcová
1,
Jakub Mrázek
1,
Giuliana Siddi
4,
Graciela N. Arenas
3 and
Giuseppe Moniello
4
1
Laboratory of Anaerobic Microbiology, Institute of Animal Physiology and Genetics, Czech Academy of Sciences, 14220 Prague, Czech Republic
2
Institute of Microbiology, University of Innsbruck, A-6020 Innsbruck, Austria
3
Área de Microbiología, Facultad de Ciencias Médicas, Universidad Nacional de Cuyo, Mendoza M5500, Argentina
4
Department of Veterinary Medicine, University of Sassari, 07100 Sassari, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Microorganisms 2021, 9(1), 157; https://doi.org/10.3390/microorganisms9010157
Submission received: 17 December 2020 / Revised: 11 January 2021 / Accepted: 11 January 2021 / Published: 12 January 2021
(This article belongs to the Special Issue Unleashing the Hidden Potential of Anaerobic Fungi)

Abstract

:
This work investigated the changes of the rumen microbiome of goats switched from a forage to a concentrate diet with special attention to anaerobic fungi (AF). Female goats were fed an alfalfa hay (AH) diet (0% grain; n = 4) for 20 days and were then abruptly shifted to a high-grain (HG) diet (40% corn grain, 60% AH; n = 4) and treated for another 10 days. Rumen content samples were collected from the cannulated animals at the end of each diet period (day 20 and 30). The microbiome structure was studied using high-throughput sequencing for bacteria, archaea (16S rRNA gene) and fungi (ITS2), accompanied by qPCR for each group. To further elucidate unclassified AF, clone library analyses were performed on the ITS1 spacer region. Rumen pH was significantly lower in HG diet fed goats, but did not induce subacute ruminal acidosis. HG diet altered prokaryotic communities, with a significant increase of Bacteroidetes and a decrease of Firmicutes. On the genus level Prevotella 1 was significantly boosted. Methanobrevibacter and Methanosphaera were the most abundant archaea regardless of the diet and HG induced a significant augmentation of unclassified Thermoplasmatales. For anaerobic fungi, HG triggered a considerable rise in Feramyces observed with both ITS markers, while a decline of Tahromyces was detected by ITS2 and decrease of Joblinomyces by ITS1 only. The uncultured BlackRhino group revealed by ITS1 and further elucidated in one sample by LSU analysis, formed a considerable part of the AF community of goats fed both diets. Results strongly indicate that the rumen ecosystem still acts as a source for novel microorganisms and unexplored microbial interactions and that initial rumen microbiota of the host animal considerably influences the reaction pattern upon diet change.

1. Introduction

The worldwide number of goats (Capra aegagrus f. hircus) is constantly increasing, together with the role of this livestock in food production [1]. According to American Goat Federation (americangoatfederation.org) there are approximately 450 million goats around the whole world. In some regions, especially the arid zones, goats are outnumbering cattle and are frequently more numerous than sheep, thus having a strong influence on the socio-economic life of human populations. The increase in demand for milk and meat production has induced intensive production systems, where goats are often fed relatively high-grain (HG) diets to achieve the maximum performance. This feeding pattern is able to saturate energy and protein requirements of the animal for maintenance and growth, but on the other hand is known to induce changes in the rumen microbiota [2,3,4,5] and often causes metabolic disorders [6,7,8].
The effect of HG diets on the health and diversity of the rumen microbiome of goats has received considerable attention in the past decade. Concentrate-rich diets induce an increase of the lactic acid and/or butyrate production/accumulation [9,10,11] and a decrease of the rumen pH that in turn results in higher proportion of Firmicutes [11,12,13]. At the genus level the reduction of Fibrobacter [14,15] and fibrolytic bacteria sensitive to low pH [9] has been described, while the genera Lactobacillus, Enterococcus, and Succiniclasticum were increased in goats fed a HG diet [10,11,14]. On the other hand, the abundance of some species, especially Prevotella, but also Butyrivibrio, is highly variable, giving inconsistent results [10,12,14,16,17,18], indicating the genetic diversity, diet specificity, and functional diversity of some rumen microbial groups [11].
All these studies show considerable influence of HG diet on the bacterial population, while information on eukaryotic members of the goat rumen ecosystem is very scarce. Metzler-Zebeli et al. [14] microscopically evaluated ciliate protozoa in the rumen fluid of goats fed different proportion of barley grain (0%, 30%, and 60%) showing that the increasing dietary grain level linearly increased the number of entodiniomorphids, whereas that of holotrich ciliates did not change. Similarly, Liu et al. [10] observed a higher protozoa number with the increase of grain (50–80%) in goat diets. However, at the onset of sub-acute ruminal acidosis (SARA), the protozoan population declined dramatically. Although it is known that anaerobic fungi (AF) (phylum Neocallimastigomycota, class Neocallimastigomycetes, order Neocallimastigales, family Neocallimastigaceae [19]), represent a major monophyletic group of symbiotic microorganisms, occupying the forestomach of ruminants and the hindgut of nonruminant herbivores, their diversity in relation to the animals’ diet in goats has remained unexplored. Together with bacteria and protozoa, AF contribute to the hydrolysis of dietary fibers resulting in the production of end fermentation products (volatile fatty acids), which are utilized by the host animals as nutrient source [20]. An impressive array of enzymes necessary for plant cell-wall degradation enables AF to break down plant cell walls, access fermentable substrates, degrade and weaken plant tissue, and reduce the size of the dietary particles [21,22,23]. Moreover, the fungal rhizoidal systems, which can be highly branched and invasive, are able to penetrate deeply into the plant stems and can aid to mechanically decompose the plant substrates [24]. Due to the capability of AF to colonize the sturdy plant structures and recalcitrant components, such as the sclerenchym and vascular system, it can be expected that Neocallimastigomycota play an especially important role in goats, which are used to live in harsh arid environments and are consuming low-quality feed rich in fiber, lignin, and tannins [25,26].
Although the diversity of goat rumen fungi has not been comprehensively studied, the recently described genus Liebetanzomyces [27] and the latest described genera Agriosomyces, Capellomyces, and Joblinomyces [28] were isolated from the feces of goats, pointing at a yet untapped source of novel fungal species associated with this host animal. The number of cultured AF has been enlarged up to 19 genera including as of yet Neocallimastix, Piromyces, Orpinomyces, Anaeromyces, Caecomyces, Cyllamyces, Buwchfawromyces, Oontomyces, Liebetanzomyces, Pecoramyces, Feramyces, Ghazallomyces, Aklioshbomyces, Agriosomyces, Capellomyces, Joblinomyces, Khoyollomyces, Tahromyces, and Aestipascuomyces, the last ten of which were described only very recently [28,29,30,31]. Among these known genera, the species of Piromyces spiralis [32], Neocallimastix californiae [33], Anaeromyces contortus [34], and the novel species of Agriosomyces longus, Capellomyces foraminis, Capellomyces elongates, and Joblinomyces apicalis [28] have been isolated from the rumen fluid and feces of domesticated and/or wild goats.
Further literature on goat AF is quite limited, involving the measurement of the hydrolytic activities of Neocallimastix sp. from rumen of native Indian and Korean goats [35,36,37] and the isolation of AF and associated methanogens of Neocallimastix frontalis, and Caecomyces communis strains from feces of the Alpine ibex (Capra ibex) [38]. However, a pyrosequencing study by Liggenstoffer et al. [39], did not detect any sequences of Neocallimastix and Caecomyces in the feces of domesticated goat, but determined Anaeromyces (48%), Piromyces (33%) and an uncultured group AL5 (20%) as prevalent genera. The cultured representative of the group AL5 was isolated recently from the feces of the domesticated goat and sheep and described as Joblinomyces [28,40]. A study by Kok et al. [41] evaluated the ability of goat rumen fungal populations to adapt to the dietary high concentrations of condensed tannins by real-time PCR. Finally, the sequences of xylan degrading enzymes of GH family 11 related to xylanases from Neocallimastix and Orpinomyces were detected in the rumen fluid of grazing goats [42].
To the best of our knowledge, there is no report on the influence of high grain diet compared to high-fiber diet on the diversity and abundance of the prokaryotic, and also the anaerobic fungal rumen microbiome of goat. A higher proportion of concentrate is one of the feeding management strategies used in intensive caprine production systems, however rational use of grain as the main energy resource has to be consistent not only with productivity but also with animal health. This study was, therefore, undertaken to evaluate both feeding regimes and their effect on the bacterial, archaeal, and fungal population in rumen fluid of fistulated goats. In addition, the role of a host-animal-specific AF population was evaluated. A culture-independent, high-throughput sequencing (HTS) approach together with qPCR analyses were chosen to analyse these microbial groups within rumen fluid of Creole goats during an alfalfa hay (AH) and high grain (HG) feeding phase. Furthermore, AF diversity was more deeply investigated by clone libraries targeting the ITS1 gene sequences.

2. Materials and Methods

2.1. Animals and Diet

This study was conducted at the Scientific and Technologic Center located in Mendoza City, Argentina. The experiment was permitted by The Committee for the Care and Use of Animals in Research of the University of Juan Agustín Maza, Mendoza, Argentina (Approval No. 102/2013) and performed in accordance with the Guide for Care and Use of Agricultural Animals in Research and Teaching of the Federal Animal Science Society [43]. Four Creole goats (meat breed, 5-year-old, BW: 40 ± 5 kg, never having an offspring) fitted with a rumen fistula were used in this study. The animals housed in pens were fed on the AH diet for a period of 20 days and were then abruptly shifted to a mixed HG diet consisting of 60% alfalfa hay and 40% corn grain for another 10 days. The diets, formulated to meet the animals’ nutrient requirements, as defined by the National Research Council [44], were composed of 1.96 or 2.33 Mcal of metabolizable energy (ME) kg−1 of dry matter with forage-to-concentrate ratios of 100:0 (AH diet) and 60:40 (HG diet), respectively. The goats were fed simultaneously once daily in the morning, had free access to water, and had not been co-housed with any other animals.

2.2. Nutritional Composition of AH and HG Diets

Samples of alfalfa hay (Medicago sativa) and corn grains were dried at 70 °C for 72 h, reduced to a particle size of less than 1 mm, and subjected to chemical analyses (Table 1): dry matter (DM), crude protein (CP), starch content, neutral detergent fiber (NDF), acid detergent fiber (ADF) [45]; acid detergent lignin (ADL) [46]; hemicellulose [47]; and cellulose [48]. Additionally, the ME content was estimated according to Aguilera [49].

2.3. Sample Collection and DNA Extraction

The rumen content (200 mL) was collected 5 h after feeding at the end of the AH diet period (day 20) and ten days after the change of the diet to HG (day 30). Simultaneously, the ruminal pH was measured with a glass electrode and clinical criteria for acute or subacute ruminal acidosis, such as decreased intake, weakness, and abdominal pain were used to evaluate the impact of the abrupt change of the diet on the animal health. A representative aliquot of each sample was freeze-dried and shipped to the Institute of Animal Physiology and Genetics, Czech Academy of Sciences (Prague, Czech Republic) for the DNA isolation and subsequent analyses. The genomic DNA was extracted using the method of Yu and Morrison [50]. The concentration and purity of the extracted nucleic acids were checked using a NanoDrop 2000c UV-Vis spectrophotometer (Thermo Scientific, Waltham, MA, USA). DNA extracts were stored at −20 °C until use.

2.4. Characterization of Goat Rumen Microbiome by High-Throughput Sequencing

The DNA extracts were subjected to an amplicon sequencing approach using the Illumina MiSeq chemistry v3 (300 bp paired-end) and targeting specific marker genes of the archaeal, bacterial and fungal population within the rumen liquid under AH and HG diet (n = 4). The sequencing was performed by Microsynth AG (Balgach, Switzerland) using the following primers: 515f (GTGCCAGCMGCCGCGGTAA) and 806r (GGACTACHVGGGTWTCTAA) targeting the V4 region of the archaeal and bacterial 16S rRNA gene [51] and ITS3 (GCATCGATGAAGAACGCAGC) and ITS4 (TCCTCCGCTTATTGATATGC) targeting the eukaryotic ITS2 region [52]. Both PCR steps (first stage of library preparation and barcoding) were performed with the KAPA Hifi Hotstart polymerase system. The conditions for the library generation for both V4 and ITS2 were: initial denaturation at 95 °C for 3 min, followed by 30 cycles (V4) or 35 cycles (ITS2) of denaturation at 98 °C for 20 s, annealing at 56 °C for 30 s and elongation at 72 °C for 30 s. The runs were completed with a final elongation step at 72 °C for 5 min. Both libraries were generated in separate reaction tubes and equimolarly pooled for subsequent barcoding.
The CoMA pipeline [53] was used for the data processing of obtained raw V4 reads. In brief, paired-end reads were merged, barcodes and primers trimmed and sequences with an average quality score of ≥25 and a median length of 253 bp of the 16S rRNA gene (V4) were kept for further analyses. For archaea and bacteria, the alignment and taxonomic affiliation of sequences were undertaken at a 97% similarity level using the blast algorithm and the SILVA SSU (release 128) and Greengenes (release 13_5) as primary and backup databases, respectively. Datasets were subsampled according to the sample with the lowest read number, singletons and doubletons removed and only OTUs (operational taxonomic units) retained that could be detected in at least two samples. For a general overview at phylum or genus-level, only OTUs with a minimum abundance of ≥1% are shown.
For analysis of the fungal community, an amplicon sequence variance (ASV)-based strategy was chosen over an OTU-based one since ASVs are more suitable to recover fungal community composition [54]. Fungal sequences (ITS2) were analyzed using the DADA2 algorithm [55], basically following the DADA2 ITS Pipeline Workflow (1.8) available on the official DADA2 homepage. In short, primers have been removed from the demultiplexed reads using cutadapt [56], then N-filtering and trimming has been done using the filterAndTrim function of the DADA2 package with the following parameters enforced: maxN = 0, maxEE = c(2,2), truncQ = 2, truncLen = c(245,200), rm.phix = TRUE, compress = TRUE). For taxonomic assignment the assignTaxonomy tool of DADA2 and the UNITE database v. 04.02.2020 [57] enriched with the ITS2 sequences of the novel AF genera [28] was used. Further analysis was done in R, mostly following the example workflow for DADA2 downstream analysis posted by Happy Belly Bioinformatics [58]. For better visualization of the diversity patterns, all ASVs that were not present at least three times in at least two samples were discarded. Changes in community composition at phylum and, if applicable, down to species level were normalized and analyzed using unfiltered ASV count data with DESeq2 [59] adapted for microbial count data [60].
SRA files of the 16S rRNA (V4) and ITS2 HTS approach were deposited to the NCBI database under BioProject PRJNA602992 and BioSamples SAMN13915736-SAMN13915751 with their respective SRA numbers SRR10961333-SRR10961348.

2.5. Characterization of Anaerobic Fungi Population Using Clone Libraries

The amplification of the anaerobic fungal ITS1 region from each sample was carried out with the combination of the fungal universal ITS1 forward primer (CTTGGTCATTTAGAGGAAGTA) [61] and the Neocallimastigomycetes-specific 5.8S rRNA gene reverse primer (GTGCAATATGCGTTCGAAGATT) [62] as described previously by Mura et al. [63]. The PCR reaction (50 µL) was performed with a PPP Master Mix kit (Top-Bio, Vestec, Czech Republic) and 0.3 μM of each primer using the following thermal cycling conditions: 33 cycles of denaturation at 94 °C for 1 min, annealing at 58 °C for 30 s and extension at 72 °C for 45 s with an initial cycle of 94 °C for 4 min and final cycle of 72 °C for 2 min. Each sample was PCR amplified in quadruplicate, pooled, and verified by agarose gel electrophoresis. The pooled PCR amplicons of the correct length (approximately 350 bp) were excised from an agarose gel with a sterile scalpel blade, purified, and concentrated using a QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany). For sample G3AH the ribosomal large subunit 28S rRNA (LSU) was also amplified in quadruplicate using Neocallimastigomycetes-specific primer set GGNL1F (5′-CATAGAGGGTGAGAATCCCGTA-3′) and GGNL4R (5′-TCAACATCCTAAGCGTAGGTA-3′) as described previously by Nagler et al. [64]. The rationale for this analysis is explained in the results section.
The TOPO® TA Cloning® Kit for Sequencing (Life Technologies, USA) was used for the preparation of the ITS1 (and, for sample G3AH, LSU) clone libraries for each sample as described by Mura et al. [63]. Plasmid DNA was isolated from 413 randomly selected ITS1 clones (in average 51 clones per library) and 35 LSU clones of E. coli using a GenElute™ HP Plasmid Miniprep Kit (Sigma-Aldrich, St. Louis, MI, USA) and sent for Sanger sequencing using the M13F priming sites within the vector pCR4 (SEQme, Dobris, Czech Republic). The quality of reads was manually verified, the vector contamination was removed and chimeras were identified by examining the multiple sequence alignment using bioinformatics software Geneious 10.3.1. (https://www.geneious.com). Unambiguous sequence data were processed in QIIME version 1.8 [65] by BLAST against ITS1 reference database of AF [66] updated with sequences of novel species. The 413 cloned sequences with taxonomic assignment from the eight libraries were concatenated into a single fasta file, and a sequence map (seqs_otus.txt) was manually built for input into the script “make_otu_table.py”. The resulting relative abundance table was summarized at the clade level. The analysis of LSU sequence data was carried out using BLAST against GenBank database. The ITS1 nucleotide sequences generated from rumen samples of goats by clone library approach have been deposited to the NCBI database under the following accession numbers: MH038788—MH039011 (AH diet) and MH039012—MH039235 (HG diet). The LSU sequences of the sample G3AH have been deposited to the NCBI database under the accession numbers MW364027—MW364061.

2.6. qPCR Analysis

The quantification of ITS1 gene copy numbers was performed using the MX 3005P QPCR System (Stratagene) and EliZyme Green MIX AddROX (Elizabeth Pharmacon). The same ITS1 primer pair and cycling conditions as described above for clone library preparation were used. The melting curve analysis with an interval from 58 to 95 °C was conducted to verify the amplification specificity. The standard curve (slope −3.258, correlation coefficient (R2) 0.991) was prepared from a 10-fold dilution series of a pCR4-TOPO plasmid containing the ITS1 region of Anaeromyces mucronatus BF2 (NCBI Accession No. MK418978). The quantitation was linear over seven orders of magnitude (102 to 108 gene copies μL−1). All samples and negative controls without a DNA template were analyzed in three replicates. Total bacteria and total archaea were quantified on a Rotor-Gene Q real-time cycler using the SensiFast SYBR® Hi-Rox chemistry (Bioline, London, UK). Freshly prepared ten-fold dilutions (107 to 102) of Nitrosomonas europea (DSMZ 21879) and Methanobrevibacter formicium (DSMZ 1535) stock DNA served as reference together with 1:10 dilutions of sample extracts, all analyzed in technical duplicates. For archaea the primer pair Parch519f/Arch915r [67,68] at a final concentration of 400 nM and for bacteria 515f/806r [51] at a concentration of 250 nM were used. Cycling consisted of 95 °C for 10 min, followed by 40 cycles of 95 °C for 20 s, 57 °C or 50°C for 20 s (archaea or bacteria, respectively), and 72 °C for 30 s each. A final melt curve (65 to 90 °C) with an increment of 0.25 °C was added after each run. R2-values of standard curves were above 0.999. The calculation of 16S rRNA copy numbers was performed with the Rotor-Gene Series Software (2.3.1) and for all investigated groups copy numbers were calculated per g of fresh matter (FM) rumen sample.

2.7. Statistical Analysis

The differences between the two diets were tested for significance using a paired t-test. In case of violation of normal distribution (Shapiro–Wilk test), the Wilcoxon singed rank test was used instead. Pearson correlation analysis was used to evaluate the (inter) dependencies of species data or correlation with surveyed diet characteristics. Significant differences were declared at the p < 0.05 level. Pearson correlation analyses and their visualizations were done in R using the packages corrplot, and lares. A quantitative profile of bacteria, archaea, and anaerobic fungi was calculated for each sample, using the product of obtained gene copy numbers and relative abundance of each clade as given by high-throughput sequencing.

3. Results

3.1. Influence of Diet on Ruminal pH

The corn grain supplementation significantly increased the starch content (p < 0.05) in HG diet (Table 1), which induced acidification of the rumen environment. The pH of the rumen fluids of HG-fed goats (6.49 ± 0.10) was significantly lower (p < 0.05) compared to AH-fed goats (6.82 ± 0.12). However, ruminal pH values of animals fed the HG diet still ranged in physiological levels. None of the goats showed any clinical or sub-clinical signs of acidosis indicating that the abrupt transition of animals from AH to HG diet did not cause either acute or subacute ruminal acidosis.

3.2. High-Throughput Sequencing of the 16S rRNA Gene V4 Region and the ITS2 Region

For V4 amplicons, a total of 34,939 reads for all eight goat rumen samples passed quality filtering of the original 71,533 reads obtained by high-throughput sequencing (HTS). While the average number of reads per sample was 4367 ± 1800, subsampling was performed at 1646 reads per sample. Based on the rarefaction curves, this allowed for a satisfactory sequencing depth (Supplementary Figure S1) to determine the core microbiome. Analysis of the V4 region revealed a bacterial and archaeal sequence distribution of 96.3 to 3.7%. Altogether 605 and 8 OTUs were detected for bacteria and archaea, respectively. OTUs that were not present at least three times in at least two samples were discarded. The OTU abundance was evenly distributed amongst all samples in the V4 dataset (273 ± 22) and the Shannon index (average 6.81 ± 0.2) indicated that all samples showed equally high diversity (Supplementary Figure S2).
After filtering, denoizing, merging, and chimera removal within the DADA2 pipeline, 322,326 ITS2 reads (60.8% of input reads) remained for taxonomic assignment. Then, 350 ASVs were identified, of which 132 ASVs remained after filtering of low abundance groups. On average across all samples, most of the filtered reads were assigned to Neocallimastigomycota (67.78 ± 3.53%), followed by Ascomycota (19.64 ± 2.15%), Basidiomycota (8.17 ± 0.98%), unidentified Fungi (4.40 ± 0.33%), and Mucoromycota (0.01 ± 0.00%).

3.2.1. Bacterial Community Composition

The bulk bacterial community consisted of the phyla Bacteroidetes (43.4 ± 7.24% AH; 55.7 ± 6.94% HG) and Firmicutes (43.8 ± 8.74% AH; 33.7 ± 6.09% HG) and some minor populations of Synergistetes, Spirochaetae, Chloroflexi, Tenericutes, Fibrobacteres, Saccharibacteria, and Verrucomicrobia in descending order (Figure 1a). All these groups except for Bacteroidetes and Spirochaetae were decreasing in relative abundance following the diet switch of AH to HG and the relative increase of the key group of Bacteroidetes was significant. At genus level a significant increase of Prevotella 1 could be observed (Supplementary Table S1). Other major genera found were Prevotella, members of the Rikenellaceae RG9 gut cluster, an uncultured rumen bacterium cluster, Christensenellaceae R-7, Butyrivibrio and Succiniclasticum to name only a few. Looking at the bacterial abundance a slight overall decrease of 16S rRNA copies from 2.08 × 1010 ± 5.36 × 109 (AH) to 1.8 × 1010 ± 9.29 × 109 (HG) was observed after diet change (Figure 1b). However, this went along with very individual responses of the bacterial abundance in each goat (Supplementary Table S2).

3.2.2. Archaeal Community Composition

After 20 days of AH diet, the archaeal rumen population was basically composed of Methanobrevibacter (two different OTUs; 69.3 ± 5.86%), Methanosphaera (18.9 ± 10.3%) and a yet unclassified Thermoplasmatales archaeon (Figure 2; Supplementary Table S1). Only in two animals (G1AH, G2AH) Candidatus Methanomethylophilus could be detected. All these genera were distributed at varying abundance, whereas G1AH and G2AH seemed to have a more similar archaeal rumen microbiome (Figure 2a). Following the diet switch to HG, a disappearance of Candidatus Methanomethylophilus was observed in the two latter mentioned animals, this genus emerged, however, at little abundance (3.3%) in G4HG (Figure 2a). There was no significant change in relative abundance for the two genera Methanobrevibacter and Methanosphaera upon diet change, however, it induced a significant increase in the unclassified Thermoplasmatales group from 10.6 ± 11% to 23.1 ± 8.11% (p < 0.01). Regarding archaeal abundance, a slight increase of 16S rRNA copies from 4.66 × 108 ± 2.1 × 108 (AH) to 4.93 × 108 ± 2.84 × 108 (HG) was observed after diet change. Archaeal rumen population rose in three goats and decreased in one animal (Figure 2b, Supplementary Table S2).

3.2.3. Fungal Community Composition

The fungal rumen microbiome after 20 days of AH diet was comprised to the largest extent of Neocallimastigomycota (47.7 ± 18.0%), Ascomycota (34.1 ± 17.7%), and Basidiomycota (14.9 ± 8.8%) (Figure 3). The community composition shifted more towards Neocallimastigomycota (88.3 ± 18.9%) after the change to HG diet. The increase in Neocallimastigomycota on HG diet was accompanied by a drop in Ascomycota (4.9 ± 3.4%), and Basidiomycota (1.3 ± 0.8%). Under both diet regimes, a small community of unclassified fungi could also be detected (3.3 ± 1.2% and 5.5 ± 3.7%, for AH and HG diet, respectively). Additionally, under AH diet a minor population of Mucoromycota was detected in two animals (Supplementary Table S1). While the individual community composition differed from animal to animal, the trend of community changes at phylum level upon diet switch were the same for all animals. DESeq2 analysis identified the drop in Ascomycota and Basidiomycota, as well as the increase in Neocallimastigomycota and unclassified Fungi upon diet change to be significant (p < 0.05; Supplementary Figure S3). It further identified 22 ASVs out of 350 ASVs to exhibit significant changes upon diet switch.
Among Ascomycota, the genus of Aspergillus formed a major part of the community in goats on AH diet (15% of filtered reads). However, the steep drop in relative abundance upon diet switch (0.27% of filtered reads in goats on HG diet) was not identified as significant by DESeq2 analysis. Instead, six other Ascomycota genera (e.g., Myrothecium, Lectera) significantly dropped, and only two Ascomycota genera (Ascochyta, and Neosetophoma) increased in relative abundance upon diet change.

3.3. Community Composition of Neocallimastigomycota Based on HTS of the ITS2

At the genus level, 12 ASVs across 4 genera of Neocallimastigomycota were identified by the DADA2 algorithm; one additional ASV could not be identified beyond belonging to Neocallimastigaceae. While individual differences between animals were less prominent at phylum level, at genus/species level those differences became clearer. For example, the Neocallimastigomycota community of goats G1 and G2 on AH diet mainly consisted of Tahromyces, and Orpinomyces, while goat G3 was dominated by unclassified Neocallimastigaceae and goat G4 was dominated by Feramyces (Figure 4c; Supplementary Table S1). Generally, it seemed that the microbiome of goats G1 and G2 was fairly similar compared to goats G3 and G4. The influence of the individual microbiome composition was especially obvious in the qPCR results: goats G1 and G2 reacted to the diet switch with a decrease by the factor of ten for total Neocallimastigomycota counts (from 3.66 × 108 to 3.1 × 107), while in goats G3 and G4 the numbers trippled from 1.29 × 108 up to 4.13 × 108. The overall trend was a drop in Neocallimastigomycota numbers by 10.1% upon diet change (Figure 4d). Furthermore, the increase of starch content in the diet led to a Neocallimastigomycota community conversion across all goats (Figure 4c). There was a significant increase in Feramyces (acc. DESeq2 analysis) after the diet switch. The rise in Feramyces in goats G1 and G2 (from 0 to 69.4% and 54.3%, respectively) was accompanied by a drop in Orpinomyces and Tahromyces and an increase in Piromyces and unclassified Neocallimastigacaea. For goat G4, the further increase of Feramyces (from 46.7 to 66.8%) led to a drop in all other identified genera (Figure 4c, Supplementary Table S1).
Broken down to ASV level, 7 ASVs assigned to Neocallimastigaceae showed significant changes upon diet switch (Supplementary Figure S3). Interestingly, while all three ASVs identified as Feramyces showed a significant increase, the other four ASVs (two Tahromyces and two Piromyces) behaved inconsistently (one of each increased; the other decreased).

3.4. Community Composition of Neocallimastigomycota Based on ITS1 Clone Libraries

The relatively high proportion of unclassified anaerobic fungi obtained by ITS2 sequencing, especially in the samples G3AH, G1HG, and G2HG, induced the further analysis of AF diversity using ITS1 barcode marker via a clone library approach. All of the 413 screened ITS1 fragments from eight clone libraries were taxonomically assigned within the phylum Neocallimastigomycota (Figure 4a). Overall, already cultured genera represented the 67.1% majority of sequences and were taxonomically assigned to Feramyces (34.4%), Joblinomyces (31.5%), and Tahromyces (1.2%). A portion of 27.3% of sequences were represented by genus-level groups called BlackRhino (23.2%) and Piromyces 4 (4.1%) (as defined by Koetschan et al. by [66]) which are currently without cultured representatives. A portion of 5.6% of sequences was assigned to JF423626, which is the accession number to a cultured strain of Neocallimastix frontalis NZRold-05 [69] (Figure 4a). In the ITS1-sequence, structure-based analysis of AF by Koetschan et al. [66]; however, this sequence clustered into the genus Orpinomyces. Thus, in the ITS1 reference database this specific AF strain is listed with the accession number only.
The shift from AH to HG diet promoted a change in relative abundances of fungal taxa (Figure 4a, Supplementary Table S3). The mean results per treatment showed that Joblinomyces was the dominant genus in AH-fed goats (44.9 ± 5.5%), followed by the BlackRhino group (25.5 ± 18.7%), Feramyces (12 ± 19.5%), sp. JF463626 (9.3 ± 5.5%), Piromyces 4 (6.6 ± 5.1%), and Tahromyces (1.7 ± 1.9%). The relative abundance of Feramyces quintupled in the HG-fed goats up to 61.1 ± 5.7%, while all the other fungal groups were decreased, in particular Joblinomyces dropped by a factor of three (14.5 ± 8.5%) (Figure 4a, Supplementary Table S3). The individual response of animals to the diet shift is evident, especially when ITS1 qPCR data are considered (Figure 4b). Goats G1 and G2 with initially no detectable sequences of Feramyces, but with high relative abundance of Joblinomyces (Figure 4a), exhibited a suppression of AF concentration (Figure 4b), decreased abundance of Joblinomyces and appearance of Feramyces under HG diet. In contrast, the HG diet induced higher AF concentration in the goats G3 and G4 (Figure 4b), characterized mainly by the enhancement of relative abundance of Feramyces (Figure 4a). The overall ruminal population of Neocallimastigomycota reacted with a slight drop from 2.47 × 108 ± 1.38 × 108 (AH) to 2.22 × 108 ± 2.45 × 108 (HG) upon diet change (Figure 4b).

3.5. Community Composition of Neocallimastigomycota in Sample G3AH Based on LSU Clone Library

In the rumen of goat G3 on AH diet the considerable 52.7% of the ITS1 sequences belonged to the monophyletic BlackRhino clade (Figure 4a). In order to get more information about taxonomic position of this uncultured clade, this sample was also analyzed using an LSU based clone library. Tree topology showed that the LSU sequences obtained (n = 35) formed four distinct clusters (Figure 5). The clade of Piromyces sp. included two clusters. Clones representing cluster I (n = 16) had the highest sequence identity (99.3–99.6%) to a sequence of a cultured Piromyces sp. strain SFH682 (MK775335) isolated from sheep feces. Clones of the sister cluster II (n = 3) had the highest sequence identity (99.5–99.8%) to a cultured Piromyces sp. strain CN12 (KY368606) isolated from goat rumen. Fourteen other clones had the highest sequence identity (98.3–99.8%) to Neocallimastix cf. cameroonii strain G3 (MG992493) isolated from sheep feces. Only two clones clustered with Feramyces sp. showing 99.5% sequence identity to Feramyces austinii strain F2a (MG605676) isolated from Barbary sheep feces. Interestingly, NCBI BLAST searching of region of local similarity between BlackRhino ITS1 (both sequences generated in this study and sequences from reference ITS1 database) and cultured species resulted in high sequence identity (99.1%) of 332 long ITS1 spacer to Piromyces sp. strain SFH682 (MK775329), which is the same strain identified here by LSU sequences analysis.

3.6. Correlation Analyses

A total of 71 parameters (diet characteristics, HTS results for archaea, bacteria, total fungi and Neocallimastigomycota, and clone library results) were subjected to a correlation analysis. Out of the 317 pairs that showed significant correlation (p < 0.05) at high levels (i.e., −0.8 > Pearson’s r > 0.8), the top parameters with the most connected partners were: Fibrobacteres (18 connections), Treponema 2 (15), Verrucomicrobia (14), Basidiomycota (13), Lachnospiraceae Genus 10 (13), Methanobrevibacter 1 (12), Lachnospira (11), and Ascomycota (11). For more details please see the correlation table and related figures (Figure 6, Supplementary Table S4). Here we focus on the strongest correlation pairs and interactions of Neocallimastigomycota with other microbial groups or physico-chemical parameters.
Out of the dietary metadata, only starch content showed correlations with the phylogenetic data. It positively correlated with Feramyces austinii (ITS1 clone library data; r +0.88), Feramyces (ITS2 HTS data; r +0.86), and with the bacterial groups of Ruminococcaceae NK4A214 (V4 HTS data; r +0.85) and Sphaerochaeta (V4 HTS data; r +0.96). Neocallimastigomycota, in general, negatively correlated with Ascomycota (r −0.99) and Basidiomycota (r −0.93). At higher phylogenetic resolution, however, Tahromyces (ITS2), Orpinomyces (ITS2), Piromyces 4 (ITS1), and Joblinomyces apicalis (ITS1) showed a high positive correlation with Ascomycota. Within archaea, Methanobrevibacter 1 was negatively correlated with Feramyces (both ITS1, and ITS2 data), but positively correlated to Tahromyces (ITS1 and ITS2), Orpinomyces (ITS2), and Joblinomyces apicalis (ITS1). In general, there were three different clusters of Neocallimastigomycota correlations: the first cluster was dominated by Feramyces (ITS1 and ITS2 data) and positively correlated with high starch content, the bacterial phylum of Bacteroidetes (V4 HTS data; r +0.72), and the archaeal group of unclassified Thermoplasmatales (r +0.76). A second cluster with a more diverse community that consisted of Orpinomyces (ITS2), Joblinomyces (ITS1), Piromyces 4 (ITS1), and Tahromyces (ITS1 and ITS2) and positively correlated with the bacterial phyla Fibrobacteres, and Verrucomicrobia and with the archaeal genera of Methanosphaera and Methanobrevibacter 1. Furthermore, the third cluster consisted of unclassified Neocallimastigomycota genera (ITS2), Piromyces (ITS2), and BlackRhino (ITS1). This last cluster generally shared little connections to the rest; only the BlackRhino group correlated positively with Methanosphaera and the bacterial genus of Lachnospira NK3A20, and negatively with Rikenellaceae RC9 Gut Group and a cluster of uncultured rumen bacteria.

4. Discussion

The animals’ diet is the major determinant of ruminal microbial composition [70,71] and it depends highly on the geographical location and the type of production system. Creole goats in Argentina are mostly kept extensively with continuous grazing on rangelands and little or no feed supplementation [72,73]. Economic aspects of meat and milk production, however, necessitate alternative feeding systems and it is a common strategy to add substantial portions of grains to diets to saturate the animals’ nutritional requirement. It is well known that feeding excessive amounts of non-structural carbohydrates results in subacute ruminal acidosis (SARA) [7]. The development of SARA seems to be dependent on the ratio of forage and grain with no signs of this metabolic disorder in goats fed up to 30% of barley grain [74]. In our study, 40% of corn grain in the diet has not induced sub-clinical signs of acidosis, despite pH-values being significantly lower in HG fed goats. However, in any of studied goats ruminal pH did not drop below 5.6, which is a value established by Khafipour et al. [75] to consider a SARA in individual animals. A study by Lee et al. [76] indicates, that goats can be more resistant to ruminal acidosis compared to cattle. This could be related to some goat-specific evolutionary adaptations of the digestive system such as large salivary glands, the large absorption area of rumen epithelium, the ability to rapidly change the rumen volume in response to environmental alterations, efficient recycling of urea, efficient use of water, and the ability to reduce their metabolism rate [25].
The impact of diet changes on the microbiome of goats, however, still remains unclear. The number of studies dealing with this topic has been increased recently [77,78], mainly with emphasis on bacterial diversity, while the response of anaerobic fungi to dietary treatment is less studied and in the goat model essentially unexplored.

4.1. Prokaryotic Response on Diet Shift

The HTS data of this study revealed the equal proportion of the two main bacterial phyla Firmicutes and Bacteroidetes in AH fed goats, while the transition from feeding on hay to corn-containing diet (40%) resulted in a significant increase of Bacteroidetes. This result is not consistent with previous reports describing that the high, readily fermented carbohydrates in HG diets lead to a higher proportion of Firmicutes [11,12,13,18,79]. In fact, in our study, Firmicutes even decreased by approximately 25% in relative abundance. The content of grain in the mentioned studies was, however, at least 65%, which might explain the detected differences. Interestingly, when comparing our results on the mature rumen microbiome of goats with the study of Wang et al. [78], Bacteroidetes were in the same range; however, Proteobacteria were clearly less abundant in our study (always below 1% compared to ≥10%). Firmicutes, in contrast, remained at higher levels in our study, even after the abundance drop that accompanied the change to HG diet. These discrepancies might be strongly influenced by the specific diet composition: the diet in the study of Wang et al. [78] was composed of a more diverse mixture (alfalfa 34.8%, rice straw 35.2%, corn 11.5%, soybean meal 5%, wheat bran 12%, salt, and mineral nutrients).
The methanogenic community composition revealed in this study corresponds well with published data [80,81,82] showing Methanobrevibacter and Methanosphaera as the two most prevalent genera in adult goats. The high incidence of Candidatus Methanomethylophilus described by Wang et al. [81], however, could not be confirmed in our study, as we detected this species in low amounts in three of the tested animals only. Furthermore, we report here for the first time a high abundance of an unclassified Thermoplasmatales group (two OTUs) in goat rumen (especially in HG fed goats), indicating that the rumen ecosystem still acts as a source for novel microorganisms and unexplored microbial interactions. Indeed, correlation analyses on all investigated parameters revealed that this group of unclassified Thermoplasmatales is positively connected with Feramyces (ITS1 and ITS2 data), and negatively connected to Ascomycota. Still, it remains unclear if this indicates a direct interaction of Feramyces and the unclassified Thermoplasmatales group, or if this co-increase under HG diet is triggered by other environmental conditions or other microbial guilds. Generally, groups such as Methanobrevibacter sp., Methanobacterium sp., and Methanosphaera have already been identified as the main methanogens associated with anaerobic fungi [38,83,84]. Looking at the key prokaryotic genera, similar members were detected for goats as described by Kumar et al. [83], mentioning Prevotella, Lachnospiraceae, Methanobrevibacter, and Methanosphaera amongst others in the rumen of dairy cows. In addition to the before mentioned groups, Belanche and coworkers [85] further detected Treponema and Ruminococcus, which were also present in our study. The study dealt with the microbial response upon diet shift from concentrate to high-fiber diet in sheep. However, the dominant AF in that study were Caecomyces sp. and Neocallimastix frontalis, which were not detected in our survey. A major bacterial group in our study, Prevotella, is a common rumen bacterium that exhibits a very versatile metabolism, able to thrive on peptides, proteins, monosaccharides and plant polysaccharides [86] and thus is often unaffected by diet changes. The increase with concentrate-rich diets, as determined in our study, was also observed by others [15,70].

4.2. Fungal Response on Diet Shift

The diversity of total fungi observed in our work was comparable to a similar study on Holstein heifers that dealt with the rumen microbiome under varying concentrate ratios in feed [18]. This study determined Neocallimastigomycota and Ascomycota to be the dominant fungal groups, but also detected a cluster of uncultured Fungi as well as a minor population of Basidiomycota. In another multi-kingdom study investigating the sheep rumen microbiome of grazing and non-grazing (concentrate added to the feed) animals, researchers could assign up to 80% of total fungal reads to the phylum Neocallimastigomycota [85]. In that study, as well as in our work, the phylum Neocallimastigomycota was especially prevalent in goats on HG diet. However, another study on dairy cows showed opposite trends [87]. Nevertheless, all three studies found Neocallimastigomycota to be inversely related to Ascomycota. High abundance of Ascomycota described previously in AH fed animals was associated with plant pathogens that can attack grass before the harvest [88], as well as with contamination by Ascomycota mold (especially Aspergillus and Fusarium) that can occur during storage processes [89]. In the present study, however, the high abundance of Ascomycota cannot solely be explained by these phenomenons, since the change in relative abundance of Aspergillus upon diet switch was not significant. Whether other fungal groups, except for Neocallimastigomycota, truly play a role in the fermentative rumen system beyond being ingested together with the feed, remains to be elucidated.
For Neocallimastigomycota in particular, it is already well established that altering the substrate availability through dietary modification affects their population structure within the digestive tract. Denman et al. [90] reported a greater diversity of Neocallimastigomycota in cattle fed high-fiber diets compared to those fed high grain-based diets (70%). This trend was also found in dairy cows [83,91], steer [92], sheep [85], and now in goats (this study). Interestingly, while Denman et al. [90] described the disappearing of the genus Neocallimastix in grain-fed cattle, higher levels of Neocallimastix and Piromyces were observed in another cattle study [91] as well as in sheep [85]. These inconsistencies on genus level might solely be due to differences in taxonomic databases/methods used, as the phylogenetic identification based on ITS 1 or ITS 2 barcode sequences alone, still remains challenging for the phylum of Neocallimastigomycota [93]. Furthermore, some of the genera, especially Piromyces, are known to be polyphyletic [94]. However, it might also be partially attributed to individual differences among host animals or even differences among individuals of the same host species (as shown in our study).
Feramyces, found as the prevalent genus in HG fed goats regardless of the used detection method, is an anaerobic fungal genus that has only been described recently. As an uncultured environmental sequence, this clade was first described by Liggenstoffer et al. [39], later denoted as AL6 by Kittelmann et al. [69], and finally renamed by Hanafy et al. [30] to Feramyces based on its first cultured representatives. This fungus is characterized by a wide substrate utilization pattern, being able to metabolize substrates often reported as not supportive for growth of other AF isolates, such as galactose, fucose, arabinose, and glucuronic acids. Moreover, some Feramyces strains displayed faster growth and higher extent on corn stover when compared with the majority of known cultivable AF strains [28]. All these features support our observation of increased abundance of this species in HG fed animals. Joblinomyces and Tahromyces are also newly described AF; however, there are no substrate preferences described yet [28]. While Tahromyces was detected in low amounts in only three animals, Joblinomyces formed a clear majority in AH fed goats (based on ITS1 cloning results), indicating that this genus might prefer the utilization of structural polysaccharides rather than storage carbohydrates. The sp. JF423626, here prevalent in AH fed goats, is a cultured, but not yet properly identified, strain, straddled between Neocallimastix and Orpinomyces (based on the molecular biological data). Both genera belong to the best studied AF and are known for their ability to utilize a wide spectrum of substrates, such as cellulose, xylose, glucose, starch, grass, and straw [95]. The features of the potentially new genus JF423626, however, have yet to be revealed.
Of the uncultured clades identified by ITS1 analysis, the BlackRhino group also represented a considerable part of sequences in all goats regardless of the diet and can thus be regarded as core member of the goat rumen. Further investigation of this group indicated, that AF of the BlackRhino clade, considered up to now as uncultured, could be identical to the already cultured strain Piromyces sp. SFH682 (MK775335) isolated from sheep feces. The fact that we could not clearly identify this group highlights the necessity of a need for another phylogenetic marker system for Neocallimastigomycota. Combining our results with qPCR data illustrated a decrease of Joblinomyces and increase of Feramyces after the addition of grain to the diet and supports the findings on strain-specific fermentation of certain types of carbohydrates [96,97]. This study thus encourages an opinion of Edwards et al. [95] about the need to re-evaluate the understanding on anaerobic fungi as simple “fiber degraders” and to reconsider their supposedly limited role in the rumen of animals fed cereal-supplemented diets.
Summarizing, after comparison of existing, corroborating, but also contradicting literature, the question remains if and to which extent these disparities arise from the host animal (e.g., cattle vs. sheep vs. goat) or varying feed compositions. Nevertheless, our results suggest that the initial, animal-specific community composition has a great influence on the microbiome shift upon the diet change. This fact was especially true for the anaerobic fungi diversity during the experiment (Figure 4). Our results highlight that 40% of concentrate proportion in the feed is not inducing SARA, but it leads to a very specific rumen microbiome that seems adapted to higher starch contents.

4.3. Comparison of the ITS1 and ITS2 as Taxonomic Marker for AF

According to literature the internal transcribed spacer (ITS) region is recommended as the universal DNA barcode marker for fungi [40,98]. Additionally, the 28S ribosomal large subunit (28S rRNA gene (LSU); D1/D2 region) might render higher coverage of early diverging lineages, such as Neocallimastigomycota, Mucoromycota, and Chytridiomycota [98]. Comparison of ITS1 and ITS2 has caused debates on which of these two regions was more appropriate for fungal phylogenetic analyses. Some authors could not detect significant superiority of either region for the classification of fungi, but rather an under- or over-representation of some fungal groups depending on the barcoding region [99]. Other authors clearly stated an advantage of using ITS1 over ITS2, in part due to shorter sequence lengths and lower GC content [100].
Here we focused on the taxonomic affiliation of anaerobic fungi by covering the whole range of the ITS region, and thus split the analyses into two approaches ITS2 (HTS) and ITS1 (cloning library), respectively. In this way we were able to combine the interpretation and compare the suitability of either fragment. The higher availability of ITS1 as compared to ITS2 fragments in the databases is also reflected in the slightly higher diversity obtained by ITS1 cloning than ITS2 HTS. Most of the uncultured AF have only been detected as environmental nucleic acid sequences (ENAS) by ITS1 primers, which can explain the missing BlackRhino group and Piromyces 4 in HTS results. However, the ITS1 region has also some disadvantages, as recently shown by Edwards et al. [101] and Yang et al. [102]. A community composition analysis can result in false higher diversity due to the heterogeneity of the ITS1 region, resulting in ambiguity in the annotation of sequences. Thus, integrating secondary structure predictions has been suggested to better classify AF sequences using ITS1 [66]. According to Edwards et al. [101] the accuracy of the method, and consequently of the results, is influenced by sample community composition, and there is an evident need to develop an alternative taxonomic marker for anaerobic fungi [93].
Nevertheless, despite differences between the results from ITS1 and ITS2 analysis and the aforementioned differences within individual animals, the trend of Feramyces dominance after the diet switch to HG was consistent.

5. Conclusions

Studies on the influence of a high grain diet on the rumen microbiome are certainly of great importance in terms of finding an optimal carbohydrate-to-fiber ratio of grain containing a diet with higher energetic and nutritional quality, but still not associated with disturbances of the ruminal microbiota leading to metabolic disorders. Strategic grain supplementation and a rational use of this nutritional resource have practical agricultural relevance for animal health and productivity.
Our results attest feasibility of abrupt diet switches to high-grain in goats, even with 40% concentrate. It seems to be a very robust system, able to quickly adapt to new nutrient conditions and still maintain a healthy microbial community. The analyses of microbial diversity and interrelationships among microbial groups suggest goat rumen to be a promising source for novel archaea (uncultured Thermoplasmatales cluster) but also anaerobic fungi isolates (BlackRhino group). Furthermore, studying this habitat more closely could reveal deeper knowledge on symbiotic interactions of AF, bacteria, and archaea. Despite a very drastic change in the composition of the anaerobic fungi population, overall absolute abundance was not significantly changed after the diet switch from AH to HG.
We could observe that, in addition to the diet, the animal-specific microbiome structure also heavily influences the response pattern upon the diet change, especially among Neocallimastigomycota.
The comparison of ITS regions, ITS1 and ITS2, for anaerobic fungi phylogenetic analysis revealed no clear advantage of one over the other. Instead, our results again highlight the need to develop a robust barcoding marker for deep-branching fungal lineages such as Neocallimastigomycota.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2607/9/1/157/s1, Table S1: Taxonomic composition of goat rumen liquid under different diets (AH at d20 and HG at d30) based on DNA amplicon sequencing targeting the ITS2 region (Fungi and Neocallimastigomycota) and the 16S rRNA gene V4 region (Archaea and Bacteria). Read abundance (for V4 at ≥ 1%) is given at phylum- and genus-level. Additionally, for each treatment mean values ± SD are shown (n = 4). Anaerobic fungi results highlighted in grey. Table S2: Abundance of anaerobic fungi, bacteria and archaea in rumen liquid samples of goats fed different diets (AH at d20 and HG at d30) determined by qPCR based on ITS1 and 16S rRNA copy number per g fresh matter (FM). Mean values of each sample (technical replicates n = 3 for anaerobic fungi and n = 2 for bacteria and archaea, respectively) and of treatments (n = 4) together with respective SD are shown. Table S3: Relative abundances of anaerobic fungal clades in the ITS1 clone libraries constructed from rumen samples of goats fed different diets (AH at d20 and HG at d30). Based on the study of Koetschan et al. (2014), the term clade is defined as a known species or an uncultivated subgroup (§) within a monophyletic lineage that has been identified using secondary structure-informed analysis of ITS1 region sequence data. Table S4: Part A lists all significant (p < 0.05) and positive (Pearson’s r > 0.8) correlation pairs of the investigated parameters for diet characteristics, NGS analysis of ITS2 and V4 and ITS2 clone library results. Only the lower triangle of the correlation matrix is shown. Part B lists all significant (p < 0.05) and negative (Pearson’s r < −0.8) correlation pairs of the investigated parameters for diet characteristics, NGS analysis of ITS2 and V4 and ITS2 clone library results. Only the lower triangle of the correlation matrix is shown. Below the total sum of correlation pairs for each parameter is exhibited. Figure S1: Rarefaction curve analysis of 16S rRNA gene amplicon sequences (V4 region; OTU-based) and the ITS2 region (ASV-based). Figure S2: Diversity analyses of 16S rRNA gene amplicon sequences (V4 region) on top and the ITS2 region (bottom) based on OUT/ASV-richness (i.e., species richness) and Shannon–Wiener index. Figure S3: Significant log2fold changes (DESeq2) of major fungal phyla (top) and Neocallimastigomycota ASVs (bottom) upon diet switch from forage (AH) to high grain (HG) diet. Figure S4: Phylogenetic relationships of anaerobic fungal LSU gene sequences generated from the sample G3AH inferred using the UPGMA method with bootstrap values from 1000 replications. The analysis involved a total of 45 nucleotide sequences (35 sequences generated from the G3AH sample, 9 sequences of the cultured genera of anaerobic fungi (Piromyces sp. MK775335, KY368606, and KY368602, N. cameroonii NG060329, KR920745, and MG992493, F. austinii MG584215, MG584217, and MG605676) plus 1 outgroup sequence.

Author Contributions

Conceptualization, K.O.F., D.J.G. and S.M.P.; methodology, K.O.F., S.M.P., J.V.; formal analysis and investigation, D.J.G.; G.N.A., S.K., D.S., H.S., G.S., J.M. and G.M.; resources, K.O.F., D.J.G. and G.M.; data curation, S.M.P., D.S. and J.V.; writing—original draft preparation, K.O.F. and S.M.P.; writing—review and editing, K.O.F., S.M.P., J.V., D.J.G. and G.M.; supervision, K.O.F., S.M.P. and G.M.; funding acquisition, S.M.P., G.M. and J.M. All authors have read and agreed to the published version of the manuscript.

Funding

S.P. acknowledges funding from the Austrian Science Fund (FWF): I 3808. G.M. acknowledges funding from the Visiting Professor Programme 2017 and 2018 of the University of Sassari and Regione Autonoma della Sardegna.

Institutional Review Board Statement

The experiment was permitted by The Committee for the Care and Use of Animals in Research of the University of Juan Agustín Maza, Mendoza, Argentina (Approval No. 102/2013).

Informed Consent Statement

Not applicable.

Data Availability Statement

The experiment was performed in accordance with the Guide for Care and Use of Agricultural Animals in Research and Teaching of the Federal Animal Science Society.

Acknowledgments

We would like to thank Sebastian Hupfauf for helpful advices during the HTS data analyses and Sophia FA Strobl for her valuable help during qPCR quantification of archaea and bacteria.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Skapetas, B.; Bampidis, V. Goat production in the world: Present situation and trends. Livest. Res. Rural Dev. 2016, 28, 200. [Google Scholar]
  2. Bi, Y.; Zeng, S.; Zhang, R.; Diao, Q.; Tu, Y. Effects of dietary energy levels on rumen bacterial community composition in Holstein heifers under the same forage to concentrate ratio condition. BMC Microbiol. 2018, 18, 1–11. [Google Scholar] [CrossRef] [PubMed]
  3. Mao, S.Y.; Zhang, R.Y.; Wang, D.S.; Zhu, W.Y. Impact of subacute ruminal acidosis (SARA) adaptation on rumen microbiota in dairy cattle using pyrosequencing. Anaerobe 2013, 24, 12–19. [Google Scholar] [CrossRef] [PubMed]
  4. Nagata, R.; Kim, Y.H.; Ohkubo, A.; Kushibiki, S.; Ichijo, T.; Sato, S. Effects of repeated subacute ruminal acidosis challenges on the adaptation of the rumen bacterial community in Holstein bulls. J. Dairy Sci. 2018, 101, 4424–4436. [Google Scholar] [CrossRef]
  5. Petri, R.M.; Schwaiger, T.; Penner, G.B.; Beauchemin, K.A.; Forster, R.J.; McKinnon, J.J.; McAllister, T.A. Changes in the rumen epimural bacterial diversity of beef cattle as affected by diet and induced ruminal acidosis. Appl. Environ. Microbiol. 2013, 79, 3744–3755. [Google Scholar] [CrossRef] [Green Version]
  6. González, L.A.; Manteca, X.; Calsamiglia, S.; Schwartzkopf-Genswein, K.S.; Ferret, A. Ruminal acidosis in feedlot cattle: Interplay between feed ingredients, rumen function and feeding behavior (a review). Anim. Feed Sci. Technol. 2012, 172, 66–79. [Google Scholar] [CrossRef]
  7. Khiaosa-ard, R.; Zebeli, Q. Diet-induced inflammation: From gut to metabolic organs and the consequences for the health and longevity of ruminants. Res. Vet. Sci. 2018, 120, 17–27. [Google Scholar] [CrossRef]
  8. Zebeli, Q.; Metzler-Zebeli, B.U. Interplay between rumen digestive disorders and diet-induced inflammation in dairy cattle. Res. Vet. Sci. 2012, 93, 1099–1108. [Google Scholar] [CrossRef]
  9. Sun, Y.Z.; Mao, S.Y.; Zhu, W.Y. Rumen chemical and bacterial changes during stepwise adaptation to a high-concentrate diet in goats. Animal 2010, 4, 210–217. [Google Scholar] [CrossRef] [Green Version]
  10. Liu, D.C.; Zhou, X.L.; Zhao, P.T.; Gao, M.; Han, H.; Hu, H.L. Effects of Increasing Non-Fiber Carbohydrate to Neutral Detergent Fiber Ratio on Rumen Fermentation and Microbiota in Goats. J. Integr. Agric. 2013, 12, 319–326. [Google Scholar] [CrossRef]
  11. Huo, W.; Zhu, W.; Mao, S. Impact of subacute ruminal acidosis on the diversity of liquid and solid-associated bacteria in the rumen of goats. World J. Microbiol. Biotechnol. 2014, 30, 669–680. [Google Scholar] [CrossRef] [PubMed]
  12. Wetzels, S.U.; Mann, E.; Metzler-Zebeli, B.U.; Wagner, M.; Klevenhusen, F.; Zebeli, Q.; Schmitz-Esser, S. Pyrosequencing reveals shifts in the bacterial epimural community relative to dietary concentrate amount in goats. J. Dairy Sci. 2015, 98, 5572–5587. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Liu, J.H.; Bian, G.R.; Zhu, W.Y.; Mao, S.Y. High-grain feeding causes strong shifts in ruminal epithelial bacterial community and expression of Toll-like receptor genes in goats. Front. Microbiol. 2015, 6, 167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Metzler-Zebeli, B.U.; Schmitz-Esser, S.; Klevenhusen, F.; Podstatzky-Lichtenstein, L.; Wagner, M.; Zebeli, Q. Grain-rich diets differently alter ruminal and colonic abundance of microbial populations and lipopolysaccharide in goats. Anaerobe 2013, 20, 65–73. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, R.; Ye, H.; Liu, J.; Mao, S. High-grain diets altered rumen fermentation and epithelial bacterial community and resulted in rumen epithelial injuries of goats. Appl. Microbiol. Biotechnol. 2017, 101, 6981–6992. [Google Scholar] [CrossRef]
  16. Grilli, D.J.; Fliegerová, K.; Kopečný, J.; Lama, S.P.; Egea, V.; Sohaefer, N.; Pereyra, C.; Ruiz, M.S.; Sosa, M.A.; Arenas, G.N.; et al. Analysis of the rumen bacterial diversity of goats during shift from forage to concentrate diet. Anaerobe 2016, 42, 17–26. [Google Scholar] [CrossRef]
  17. Liu, K.; Xu, Q.; Wang, L.; Wang, J.; Guo, W.; Zhou, M. The impact of diet on the composition and relative abundance of rumen microbes in goat. Asian-Australas. J. Anim. Sci. 2017, 30, 531–537. [Google Scholar] [CrossRef]
  18. Zhang, J.; Shi, H.; Wang, Y.; Li, S.; Cao, Z.; Ji, S.; He, Y.; Zhang, H. Effect of dietary forage to concentrate ratios on dynamic profile changes and interactions of ruminal microbiota and metabolites in holstein heifers. Front. Microbiol. 2017, 8, 2206. [Google Scholar] [CrossRef] [Green Version]
  19. Tedersoo, L.; Sánchez-Ramírez, S.; Kõljalg, U.; Bahram, M.; Döring, M.; Schigel, D.; May, T.; Ryberg, M.; Abarenkov, K. High-level classification of the Fungi and a tool for evolutionary ecological analyses. Fungal Divers. 2018, 90, 135–159. [Google Scholar] [CrossRef] [Green Version]
  20. Engelking, L.R. Overview of carbohydrate metabolism. In Textbook of Veterinary Physiological Chemistry; Academic Press: Amsterdam, The Netherlands, 2015; pp. 120–124. ISBN 978-0-12-384852-9. [Google Scholar]
  21. Ljungdahl, L.G. The cellulase/hemicellulase system of the anaerobic fungus Orpinomyces PC-2 and aspects of its applied use. Ann. N. Y. Acad. Sci. 2008, 1125, 308–321. [Google Scholar] [CrossRef]
  22. Couger, M.B.; Youssef, N.H.; Struchtemeyer, C.G.; Liggenstoffer, A.S.; Elshahed, M.S. Transcriptomic analysis of lignocellulosic biomass degradation by the anaerobic fungal isolate Orpinomyces sp. strain C1A. Biotechnol. Biofuels 2015, 8, 1–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Solomon, K.V.; Haitjema, C.H.; Henske, J.K.; Gilmore, S.P.; Borges-Rivera, D.; Lipzen, A.; Brewer, H.M.; Purvine, S.O.; Wright, A.T.; Theodorou, M.K.; et al. Early-branching gut fungi possess large, comprehensive array of biomass-degrading enzymes. Science 2016, 351, 1192–1195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Orpin, C.G. Anaerobic fungi: Taxonomy, biology and distribution in nature. In Anaerobic Fungi: Biology, Ecology and Function; Mountford, D.O., Orpin, C.G., Eds.; Marcel Dekker: New York, NY, USA, 1994; pp. 1–46. [Google Scholar]
  25. Silanikove, N. The physiological basis of adaptation in goats to harsh environments. Small Rumin. Res. 2000, 35, 181–193. [Google Scholar] [CrossRef]
  26. Li, D.B.; Hou, X.Z. Effect of fungal elimination on bacteria and protozoa populations and degradation of straw dry matter in the rumen of sheep and goats. Asian-Australas. J. Anim. Sci. 2007, 20, 70–74. [Google Scholar] [CrossRef]
  27. Joshi, A.; Lanjekar, V.B.; Dhakephalkar, P.K.; Callaghan, T.M.; Griffith, G.W.; Dagar, S.S. Liebetanzomyces polymorphus gen. et sp. nov., a new anaerobic fungus (Neocallimastigomycota) isolated from the rumen of a goat. MycoKeys 2018, 40, 89–110. [Google Scholar] [CrossRef]
  28. Hanafy, R.A.; Lanjekar, V.B.; Dhakephalkar, P.K.; Callaghan, T.M.; Dagar, S.S.; Griffith, G.W.; Elshahed, M.S.; Youssef, N.H. Seven new Neocallimastigomycota genera from wild, zoo-housed, and domesticated herbivores greatly expand the taxonomic diversity of the phylum. Mycologia 2020, 112, 1212–1239. [Google Scholar] [CrossRef]
  29. Hanafy, R.A.; Elshahed, M.S.; Liggenstoffer, A.S.; Griffith, G.W.; Youssef, N.H. Pecoramyces ruminantium, gen. nov., sp. nov., an anaerobic gut fungus from the feces of cattle and sheep. Mycologia 2017, 109, 231–243. [Google Scholar] [CrossRef]
  30. Hanafy, R.A.; Elshahed, M.S.; Youssef, N.H. Feramyces austinii, gen. nov., sp. nov., an anaerobic gut fungus from rumen and fecal samples of wild barbary sheep and fallow deer. Mycologia 2018, 110, 513–525. [Google Scholar] [CrossRef]
  31. Stabel, M.; Hanafy, R.A.; Schweitzer, T.; Greif, M.; Aliyu, H.; Flad, V.; Young, D.; Lebuhn, M.; Elshahed, M.S.; Ochsenreither, K.; et al. Aestipascuomyces dupliciliberans gen. nov, sp. nov., the first cultured representative of the uncultured SK4 clade from aoudad sheep and alpaca. Microorganisms 2020, 8, 1734. [Google Scholar] [CrossRef]
  32. Ho, Y.; Barr, D.J.S.; Abdullah, N.; Jalaludin, S.; Kudo, H. Piromyces spiralis, a new species of anaerobic fungus from the rumen of goat. Mycotaxon 1993, 48, 59–68. [Google Scholar]
  33. Li, G.J.; Hyde, K.D.; Zhao, R.L.; Hongsanan, S.; Abdel-Aziz, F.A.; Abdel-Wahab, M.A.; Alvarado, P.; Alves-Silva, G.; Ammirati, J.F.; Ariyawansa, H.A.; et al. Fungal diversity notes 253–366: Taxonomic and phylogenetic contributions to fungal taxa. Fungal Divers. 2016, 78, 1–237. [Google Scholar] [CrossRef]
  34. Hanafy, R.A.; Johnson, B.; Elshahed, M.S.; Youssef, N.H. Anaeromyces contortus, sp. nov., a new anaerobic gut fungal species (Neocallimastigomycota) isolated from the feces of cow and goat. Mycologia 2018, 110, 502–512. [Google Scholar] [CrossRef] [PubMed]
  35. Ha, J.K.; Lee, S.S.; Gao, Z.; Kim, C.H.; Kim, S.W.; Ko, J.Y.; Cheng, K.J. The effect of saturates fatty acid on cellulose digestion by the rumen anaerobic fungus, Neocallimastix frontalis C5-1. Asian-Australas. J. Anim. Sci. J. Anim. Sci 2001, 14, 941–946. [Google Scholar] [CrossRef]
  36. Nagpal, R.; Puniya, A.K.; Sehgal, J.P.; Singh, K. In vitro fibrolytic potential of anaerobic rumen fungi from ruminants and non-ruminant herbivores. Mycoscience 2011, 52, 31–38. [Google Scholar] [CrossRef]
  37. Thareja, A.; Puniya, A.K.; Goel, G.; Nagpal, R.; Sehgal, J.P.; Singh, P.K.; Singh, K. In vitro degradation of wheat straw by anaerobic fungi from small ruminants. Arch. Anim. Nutr. 2006, 60, 412–417. [Google Scholar] [CrossRef]
  38. Leis, S.; Dresch, P.; Peintner, U.; Fliegerová, K.; Sandbichler, A.M.; Insam, H.; Podmirseg, S.M. Finding a robust strain for biomethanation: Anaerobic fungi (Neocallimastigomycota) from the Alpine ibex (Capra ibex) and their associated methanogens. Anaerobe 2014, 29, 34–43. [Google Scholar] [CrossRef]
  39. Liggenstoffer, A.S.; Youssef, N.H.; Couger, M.B.; Elshahed, M.S. Phylogenetic diversity and community structure of anaerobic gut fungi (phylum Neocallimastigomycota) in ruminant and non-ruminant herbivores. ISME J. 2010, 4, 1225–1235. [Google Scholar] [CrossRef]
  40. Vinzelj, J.; Joshi, A.; Insam, H.; Podmirseg, S.M. Employing anaerobic fungi in biogas production: Challenges & opportunities. Bioresour. Technol. 2020, 300, 122687. [Google Scholar] [CrossRef]
  41. Kok, C.M.; Sieo, C.C.; Tan, H.Y.; Saad, W.Z.; Liang, J.B.; Ho, Y.W. Anaerobic cellulolytic rumen fungal populations in goats fed with and without Leucaena leucocephala hybrid, as determined by real-time PCR. J. Microbiol. 2013, 51, 700–703. [Google Scholar] [CrossRef] [Green Version]
  42. Wang, G.; Luo, H.; Meng, K.; Wang, Y.; Huang, H.; Shi, P.; Pan, X.; Yang, P.; Diao, Q.; Zhang, H.; et al. High genetic diversity and different distributions of glycosyl hydrolase family 10 and 11 xylanases in the goat rumen. PLoS ONE 2011, 6, e16731. [Google Scholar] [CrossRef] [Green Version]
  43. FASS. Guide for the Care and Use of Agricultural Animals in Research and Teaching, 3rd ed.; Swanson, J.C., McGlone, J.J., Eds.; Federation of Animal Science Societies Inc.: Champaign, IL, USA, 2010. [Google Scholar]
  44. National Research Council (U.S.). Nutrient Requirements of Small Ruminants: Sheep, Goats, Cervids, and New World Camelids; National Academies Press: Washington, DC, USA, 2007.
  45. AOAC International. Official Methods of AOAC International, 17th ed.; Horwitz, W., Latimer, G.W., Eds.; AOAC International: Gaithersburg, MA, USA, 2005; ISBN 0-935584-77-3. [Google Scholar]
  46. Goering, H.K.; van Soest, P.J. Forage Fiber Analysis (Apparatus Reagents, Procedures and Some Applications); Agriculture Handbook; United States Department of Agriculture, Government Printing Office: Washington, DC, USA, 1970. [Google Scholar]
  47. Van Soest, P.J.; Robertson, J.B.; Lewis, B.A. Methods for Dietary Fiber, Neutral Detergent Fiber, and Nonstarch Polysaccharides in Relation to Animal Nutrition. J. Dairy Sci. 1991, 74, 3583–3597. [Google Scholar] [CrossRef]
  48. Crampton, E.W.; Maynard, L.A. The relation of cellulose and lignin content to the nutritive value of animal feeds. J. Nutr. 1938, 15, 383–395. [Google Scholar] [CrossRef]
  49. Aguilera, J.F. Contributions to the knowledge of the energy nutrition of small ruminants, with particular reference to the goat. Arch. Zootec. 2001, 50, 565–596. [Google Scholar]
  50. Yu, Z.; Morrison, M. Improved extraction of PCR-quality community DNA from digesta and fecal samples. Biotechniques 2004, 36, 808–812. [Google Scholar] [CrossRef] [PubMed]
  51. Caporaso, J.G.; Lauber, C.L.; Walters, W.A.; Berg-Lyons, D.; Lozupone, C.A.; Turnbaugh, P.J.; Fierer, N.; Knight, R. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. USA 2011, 108, 4516–4522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. White, T.J.; Bruns, T.D.; Lee, S.B.; Taylor, J.W. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In PCR Protocols: A Guide to Methods and Applications; Innis, M.A., Gelfand, D.H., Sninsky, J.J., White, T.J., Eds.; Academic Press: London, UK, 1990; pp. 315–322. [Google Scholar]
  53. Hupfauf, S.; Etemadi, M.; Fernández-Delgado Juárez, M.; Gómez-Brandón, M.; Insam, H.; Podmirseg, S.M. CoMA—An intuitive and user-friendly pipeline for amplicon-sequencing data analysis. PLoS ONE 2020, 15, e0243241. [Google Scholar] [CrossRef] [PubMed]
  54. Pauvert, C.; Buée, M.; Laval, V.; Edel-Hermann, V.; Fauchery, L.; Gautier, A.; Lesur, I.; Vallance, J.; Vacher, C. Bioinformatics matters: The accuracy of plant and soil fungal community data is highly dependent on the metabarcoding pipeline. Fungal Ecol. 2019, 41, 23–33. [Google Scholar] [CrossRef]
  55. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [Green Version]
  56. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  57. Abarenkov, K.; Zirk, A.; Piirmann, T.; Pöhönen, R.; Ivanov, F.; Nilsson, R.; Henrik, R.H.; Kõljalg, U. UNITE General FASTA Release for Fungi 2; Version 04.02.2020; UNITE Community: London, UK, 2020. [Google Scholar]
  58. Lee, M. Happy Belly Bioinformatics: An open-source resource dedicated to helping biologists utilize bioinformatics. J. Open Source Educ. 2019, 2, 53. [Google Scholar] [CrossRef]
  59. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 1–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. McMurdie, P.J.; Holmes, S. Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Gardes, M.; Bruns, T.D. ITS primers with enhanced specificity for Basidiomycetes—Application to the identification of mycorrhizas and rusts. Mol. Ecol. 1993, 2, 113–118. [Google Scholar] [CrossRef] [PubMed]
  62. Edwards, J.E.; Kingston-Smith, A.H.; Jimenez, H.R.; Huws, S.A.; Skøt, K.P.; Griffith, G.W.; McEwan, N.R.; Theodorou, M.K. Dynamics of initial colonization of nonconserved perennial ryegrass by anaerobic fungi in the bovine rumen. FEMS Microbiol. Ecol. 2008, 66, 537–545. [Google Scholar] [CrossRef] [PubMed]
  63. Mura, E.; Edwards, J.; Kittelmann, S.; Kaerger, K.; Voigt, K.; Mrázek, J.; Moniello, G.; Fliegerova, K. Anaerobic fungal communities differ along the horse digestive tract. Fungal Biol. 2019, 123, 240–246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Nagler, M.; Podmirseg, S.M.; Griffith, G.W.; Insam, H.; Ascher-Jenull, J. The use of extracellular DNA as a proxy for specific microbial activity. Appl. Microbiol. Biotechnol. 2018, 102, 2885–2898. [Google Scholar] [CrossRef] [Green Version]
  65. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Peña, A.G.; Goodrich, K.; Gordon, J.I.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  66. Koetschan, C.; Kittelmann, S.; Lu, J.; Al-Halbouni, D.; Jarvis, G.N.; Müller, T.; Wolf, M.; Janssen, P.H. Internal transcribed spacer 1 secondary structure analysis reveals a common core throughout the anaerobic fungi (Neocallimastigomycota). PLoS ONE 2014, 9, e91928. [Google Scholar] [CrossRef] [Green Version]
  67. Øvreås, L.; Forney, L.; Daae, F.L.; Torsvik, V. Distribution of bacterioplankton in meromictic lake Saelenvannet, as determined by denaturing gradient gel electrophoresis of PCR-amplified gene fragments coding for 16S rRNA. Appl. Environ. Microbiol. 1997, 63, 3367–3373. [Google Scholar] [CrossRef] [Green Version]
  68. Stahl, D.A.; Amann, R. Development and application of nucleic acid probes in bacterial systematics. In Nucleic Acid Techniques in Bacterial Systematics; Stackebrandt, E., Goodfellow, M., Eds.; Wiley & Sons Ltd.: Chichester, UK, 1991; pp. 205–248. [Google Scholar]
  69. Kittelmann, S.; Naylor, G.E.; Koolaard, J.P.; Janssen, P.H. A proposed taxonomy of anaerobic fungi (class neocallimastigomycetes) suitable for large-scale sequence-based community structure analysis. PLoS ONE 2012, 7, e36866. [Google Scholar] [CrossRef]
  70. Henderson, G.; Cox, F.; Ganesh, S.; Jonker, A.; Young, W.; Janssen, P.H.; Abecia, L.; Angarita, E.; Aravena, P.; Arenas, G.N.; et al. Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci. Rep. 2015, 5, 14567. [Google Scholar] [CrossRef] [PubMed]
  71. Gruninger, R.J.; Ribeiro, G.O.; Cameron, A.; McAllister, T.A. Invited review: Application of meta-omics to understand the dynamic nature of the rumen microbiome and how it responds to diet in ruminants. Animal 2019, 13, 1843–1854. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Guevara, J.C.; Grünwaldt, E.G.; Estevez, O.R.; Bisigato, A.J.; Blanco, L.J.; Biurrun, F.N.; Ferrando, C.A.; Chirino, C.C.; Morici, E.; Fernández, B.; et al. Range and livestock production in the Monte Desert, Argentina. J. Arid Environ. 2009, 73, 228–237. [Google Scholar] [CrossRef]
  73. Allegretti, L.; Sartor, C.; Paez Lama, S.; Egea, V.; Fucili, M.; Passera, C. Effect of the physiological state of Criollo goats on the botanical composition of their diet in NE Mendoza, Argentina. Small Rumin. Res. 2012, 103, 152–157. [Google Scholar] [CrossRef]
  74. Klevenhusen, F.; Hollmann, M.; Podstatzky-Lichtenstein, L.; Krametter-Frötscher, R.; Aschenbach, J.R.; Zebeli, Q. Feeding barley grain-rich diets altered electrophysiological properties and permeability of the ruminal wall in a goat model. J. Dairy Sci. 2013, 96, 2293–2302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Khafipour, E.; Li, S.; Plaizier, J.C.; Krause, D.O. Rumen microbiome composition determined using two nutritional models of subacute ruminal acidosis. Appl. Environ. Microbiol. 2009, 75, 7115–7124. [Google Scholar] [CrossRef] [Green Version]
  76. Lee, M.; Jeong, S.; Seo, J.; Seo, S. Changes in the ruminal fermentation and bacterial community structure by a sudden change to a high-concentrate diet in Korean domestic ruminants. Asian-Australas. J. Anim. Sci. 2019, 32, 92–102. [Google Scholar] [CrossRef] [Green Version]
  77. Han, X.; Yang, Y.; Yan, H.; Wang, X.; Qu, L.; Chen, Y. Rumen bacterial diversity of 80 to 110-day-old goats using 16S rRNA sequencing. PLoS ONE 2015, 10, e0117811. [Google Scholar] [CrossRef]
  78. Wang, L.; Xu, Q.; Kong, F.; Yang, Y.; Wu, D.; Mishra, S.; Li, Y. Exploring the goat rumen microbiome from seven days to two years. PLoS ONE 2016, 11, e0154354. [Google Scholar] [CrossRef] [Green Version]
  79. Hua, C.; Tian, J.; Tian, P.; Cong, R.; Luo, Y.; Geng, Y.; Tao, S.; Ni, Y.; Zhao, R. Feeding a high concentration diet induces unhealthy alterations in the composition and metabolism of ruminal microbiota and host response in a goat model. Front. Microbiol. 2017, 8, 138. [Google Scholar] [CrossRef] [Green Version]
  80. Cunha, I.S.; Barreto, C.C.; Costa, O.Y.A.; Bomfim, M.A.; Castro, A.P.; Kruger, R.H.; Quirino, B.F. Bacteria and Archaea community structure in the rumen microbiome of goats (Capra hircus) from the semiarid region of Brazil. Anaerobe 2011, 17, 118–124. [Google Scholar] [CrossRef] [PubMed]
  81. Wang, Z.; Elekwachi, C.O.; Jiao, J.; Wang, M.; Tang, S.; Zhou, C.; Tan, Z.; Forster, R.J. Investigation and manipulation of metabolically active methanogen community composition during rumen development in black goats. Sci. Rep. 2017, 7, 422. [Google Scholar] [CrossRef] [PubMed]
  82. Ebrahimi, S.H.; Valizade, R.; Heidarian Miri, V. Rumen microbial community of Saanen goats adapted to a high-fiber diet in the Northeast of Iran. Iran. J. Appl. Anim. Sci. 2018, 8, 271–279. [Google Scholar]
  83. Kumar, S.; Indugu, N.; Vecchiarelli, B.; Pitta, D.W. Associative patterns among anaerobic fungi, methanogenic archaea, and bacterial communities in response to changes in diet and age in the rumen of dairy cows. Front. Microbiol. 2015, 6, 781. [Google Scholar] [CrossRef] [Green Version]
  84. Jin, W.; Cheng, Y.F.; Mao, S.Y.; Zhu, W.Y. Isolation of natural cultures of anaerobic fungi and indigenously associated methanogens from herbivores and their bioconversion of lignocellulosic materials to methane. Bioresour. Technol. 2011, 102, 7925–7931. [Google Scholar] [CrossRef]
  85. Belanche, A.; Kingston-Smith, A.H.; Griffith, G.W.; Newbold, C.J. A multi-kingdom study reveals the plasticity of the rumen microbiota in response to a shift from non-grazing to grazing diets in sheep. Front. Microbiol. 2019, 10. [Google Scholar] [CrossRef] [Green Version]
  86. Deusch, S.; Camarinha-Silva, A.; Conrad, J.; Beifuss, U.; Rodehutscord, M.; Seifert, J. A structural and functional elucidation of the rumen microbiome influenced by various diets and microenvironments. Front. Microbiol. 2017, 8, 1605. [Google Scholar] [CrossRef] [Green Version]
  87. Ishaq, S.L.; AlZahal, O.; Walker, N.; McBride, B. An investigation into rumen fungal and protozoal diversity in three rumen fractions, during high-fiber or grain-induced sub-acute ruminal acidosis conditions, with or without active dry yeast supplementation. Front. Microbiol. 2017, 8, 1943. [Google Scholar] [CrossRef] [Green Version]
  88. Alexander, H.M. Disease in natural plant populations, communities, and ecosystems: Insights into ecological and evolutionary processes. Plant Dis. 2010, 94, 492–503. [Google Scholar] [CrossRef] [Green Version]
  89. Gallo, A.; Giuberti, G.; Frisvad, J.C.; Bertuzzi, T.; Nielsen, K.F. Review on mycotoxin issues in ruminants: Occurrence in forages, effects of mycotoxin ingestion on health status and animal performance and practical strategies to counteract their negative effects. Toxins 2015, 7, 3057–3111. [Google Scholar] [CrossRef]
  90. Denman, S.E.; Nicholson, M.J.; Brookman, J.L.; Theodorou, M.K.; McSweeney, C.S. Detection and monitoring of anaerobic rumen fungi using an ARISA method. Lett. Appl. Microbiol. 2008, 47, 492–499. [Google Scholar] [CrossRef] [PubMed]
  91. Tapio, I.; Fischer, D.; Blasco, L.; Tapio, M.; Wallace, R.J.; Bayat, A.R.; Ventto, L.; Kahala, M.; Negussie, E.; Shingfield, K.J.; et al. Taxon abundance, diversity, co-occurrence and network analysis of the ruminal microbiota in response to dietary changes in dairy cows. PLoS ONE 2017, 12, 1–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Boots, B.; Lillis, L.; Clipson, N.; Petrie, K.; Kenny, D.A.; Boland, T.M.; Doyle, E. Responses of anaerobic rumen fungal diversity (phylum Neocallimastigomycota) to changes in bovine diet. J. Appl. Microbiol. 2013, 114, 626–635. [Google Scholar] [CrossRef] [PubMed]
  93. Hanafy, R.A.; Johnson, B.; Youssef, N.H.; Elshahed, M.S. Assessing anaerobic gut fungal diversity in herbivores using D1/D2 large ribosomal subunit sequencing and multi-year isolation. Environ. Microbiol. 2020, 22, 3883–3908. [Google Scholar] [CrossRef]
  94. Paul, S.S.; Bu, D.; Xu, J.; Hyde, K.D.; Yu, Z. A phylogenetic census of global diversity of gut anaerobic fungi and a new taxonomic framework. Fungal Divers. 2018, 89, 253–266. [Google Scholar] [CrossRef]
  95. Edwards, J.E.; Forster, R.J.; Callaghan, T.M.; Dollhofer, V.; Dagar, S.S.; Cheng, Y.; Chang, J.; Kittelmann, S.; Fliegerova, K.; Puniya, A.K.; et al. PCR and omics based techniques to study the diversity, ecology and biology of anaerobic fungi: Insights, challenges and opportunities. Front. Microbiol. 2017, 8, 1657. [Google Scholar] [CrossRef] [Green Version]
  96. Trinci, A.P.J.; Davies, D.R.; Gull, K.; Lawrence, M.I.; Bonde Nielsen, B.; Rickers, A.; Theodorou, M.K. Anaerobic fungi in herbivorous animals. Mycol. Res. 1994, 98, 129–152. [Google Scholar] [CrossRef]
  97. Orpin, C.G.; Joblin, K.N. The rumen anaerobic fungi. In The Rumen Microbial Ecosystem; Hobson, P.N., Stewart, C.S., Eds.; Springer: Dordrecht, Germany, 1997; pp. 140–195. ISBN 978-94-009-1453-7. [Google Scholar]
  98. Schoch, C.L.; Seifert, K.A.; Huhndorf, S.; Robert, V.; Spouge, J.L.; Levesque, C.A.; Chen, W.; Bolchacova, E.; Voigt, K.; Crous, P.W.; et al. Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. Proc. Natl. Acad. Sci. USA 2012, 109, 6241–6246. [Google Scholar] [CrossRef] [Green Version]
  99. Blaalid, R.; Kumar, S.; Nilsson, R.H.; Abarenkov, K.; Kirk, P.M.; Kauserud, H. ITS1 versus ITS2 as DNA metabarcodes for fungi. Mol. Ecol. Resour. 2013, 13, 218–224. [Google Scholar] [CrossRef] [Green Version]
  100. Wang, X.C.; Liu, C.; Huang, L.; Bengtsson-Palme, J.; Chen, H.; Zhang, J.H.; Cai, D.; Li, J.Q. ITS1: A DNA barcode better than ITS2 in eukaryotes? Mol. Ecol. Resour. 2015, 15, 573–586. [Google Scholar] [CrossRef]
  101. Edwards, J.E.; Hermes, G.D.A.; Kittelmann, S.; Nijsse, B.; Smidt, H. Assessment of the Accuracy of High-Throughput Sequencing of the ITS1 Region of Neocallimastigomycota for Community Composition Analysis. Front. Microbiol. 2019, 10, 2317. [Google Scholar] [CrossRef] [PubMed]
  102. Yang, R.H.; Su, J.H.; Shang, J.J.; Wu, Y.Y.; Li, Y.; Bao, D.P.; Yao, Y.J. Evaluation of the ribosomal DNA internal transcribed spacer (ITS), specifically ITS1 and ITS2, for the analysis of fungal diversity by deep sequencing. PLoS ONE 2018, 13, e0206428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (a) Relative ruminal bacterial community composition of goats fed an AH diet (left) or HG diet (middle) and mean values of both treatments (right), respectively (phylum-level). Results are based on amplicon sequencing of the V4 region of the 16S rRNA gene and phyla ≥ 1% abundances are depicted; (b) Absolute bacterial community composition based on relative abundance and total bacterial number as retrieved via qPCR. Significant changes between treatments are indicated by an asterisk.
Figure 1. (a) Relative ruminal bacterial community composition of goats fed an AH diet (left) or HG diet (middle) and mean values of both treatments (right), respectively (phylum-level). Results are based on amplicon sequencing of the V4 region of the 16S rRNA gene and phyla ≥ 1% abundances are depicted; (b) Absolute bacterial community composition based on relative abundance and total bacterial number as retrieved via qPCR. Significant changes between treatments are indicated by an asterisk.
Microorganisms 09 00157 g001
Figure 2. (a) Relative ruminal archaeal community composition of goats fed an AH diet (left) or HG diet (middle) and mean values of both treatments (right), respectively. Results are based on amplicon sequencing of the V4 region of the 16S rRNA gene and genera ≥ 1% abundance are depicted; (b) absolute archaeal community composition based on relative abundance and total archaeal number as retrieved via qPCR. Significant changes between treatments are indicated by an asterisk.
Figure 2. (a) Relative ruminal archaeal community composition of goats fed an AH diet (left) or HG diet (middle) and mean values of both treatments (right), respectively. Results are based on amplicon sequencing of the V4 region of the 16S rRNA gene and genera ≥ 1% abundance are depicted; (b) absolute archaeal community composition based on relative abundance and total archaeal number as retrieved via qPCR. Significant changes between treatments are indicated by an asterisk.
Microorganisms 09 00157 g002
Figure 3. Relative ruminal, fungal community composition at phylum-level of goats fed an AH diet (left) or HG diet (middle) and mean values of both treatments (right), respectively. Results are based on amplicon sequencing of the ITS2 region and phyla ≥ 1% abundance are depicted. Significant changes between treatments are indicated by an asterisk (t-test on relative genus-abundance data) or hash (based on DESeq2 analysis of ASV raw data).
Figure 3. Relative ruminal, fungal community composition at phylum-level of goats fed an AH diet (left) or HG diet (middle) and mean values of both treatments (right), respectively. Results are based on amplicon sequencing of the ITS2 region and phyla ≥ 1% abundance are depicted. Significant changes between treatments are indicated by an asterisk (t-test on relative genus-abundance data) or hash (based on DESeq2 analysis of ASV raw data).
Microorganisms 09 00157 g003
Figure 4. Neocallimastigomycota community composition of rumen samples of goats fed either AH or HG and mean values of both treatments: (a) relative and (b) absolute abundance of anaerobic fungi as detected by ITS1 clone library approach. The number of analyzed clones is given in brackets. (c) Relative and (d) absolute abundance of anaerobic fungi as detected by amplicon sequencing of the ITS2 region. Absolute abundance is based on qPCR results. Significant changes between treatments are indicated by an asterisk (t-test on relative genus-abundance data) or hash (based on DESeq2 analysis of ASV raw data).
Figure 4. Neocallimastigomycota community composition of rumen samples of goats fed either AH or HG and mean values of both treatments: (a) relative and (b) absolute abundance of anaerobic fungi as detected by ITS1 clone library approach. The number of analyzed clones is given in brackets. (c) Relative and (d) absolute abundance of anaerobic fungi as detected by amplicon sequencing of the ITS2 region. Absolute abundance is based on qPCR results. Significant changes between treatments are indicated by an asterisk (t-test on relative genus-abundance data) or hash (based on DESeq2 analysis of ASV raw data).
Microorganisms 09 00157 g004
Figure 5. Phylogenetic relationships of anaerobic fungal LSU gene sequences generated from the sample G3AH inferred using the UPGMA method with bootstrap values from 1000 replications. The analysis involved a total of 45 nucleotide sequences (35 sequences generated from the G3AH sample, 9 sequences of the cultured genera of anaerobic fungi (Piromyces sp. MK775335, KY368606, and KY368602, N. camerooni NG060329, KR920745, and MG992493, F. austinii MG584215, MG584217, and MG605676) plus 1 outgroup sequence. An uncompressed phylogenetic tree is shown in (Supplementary Material Figure S4).
Figure 5. Phylogenetic relationships of anaerobic fungal LSU gene sequences generated from the sample G3AH inferred using the UPGMA method with bootstrap values from 1000 replications. The analysis involved a total of 45 nucleotide sequences (35 sequences generated from the G3AH sample, 9 sequences of the cultured genera of anaerobic fungi (Piromyces sp. MK775335, KY368606, and KY368602, N. camerooni NG060329, KR920745, and MG992493, F. austinii MG584215, MG584217, and MG605676) plus 1 outgroup sequence. An uncompressed phylogenetic tree is shown in (Supplementary Material Figure S4).
Microorganisms 09 00157 g005
Figure 6. Correlogram of Pearson’s correlation pairs. Based on relative abundance data from the V4 HTS (OTU based), ITS2 HTS (ASV-based), and ITS1 clone library approaches. Only pairs with p < 0.05 and −0.7 < r < 0.7 are shown.
Figure 6. Correlogram of Pearson’s correlation pairs. Based on relative abundance data from the V4 HTS (OTU based), ITS2 HTS (ASV-based), and ITS1 clone library approaches. Only pairs with p < 0.05 and −0.7 < r < 0.7 are shown.
Microorganisms 09 00157 g006
Table 1. Total daily intake and daily intake of food, nutrients and metabolizable energy (ME) of Creole goats included in the experimental period. Mean values ± SD.
Table 1. Total daily intake and daily intake of food, nutrients and metabolizable energy (ME) of Creole goats included in the experimental period. Mean values ± SD.
Diets
Alfalfa HayHigh Grain
Total daily intake (g DM)1061 ± 173 a1252 ± 193 a
Daily food intake (g DM):
Alfalfa Hay1061 ± 173 a900 ± 193 a
Corn-352
Nutrient intake (g DM):
CP168 ± 27 a171 ± 30 a
NDF361 ± 59 a335 ± 65 a
ADF228 ± 37 a204 ± 41 a
ADL42 ± 7 a39 ± 7 a
Hemicellulose133 ± 22 a132 ± 24 a
Cellulose183 ± 30 a163 ± 33 a
Starch27 ± 4 a244 ± 5 b
Daily intake of ME (MCal kg DM−1)2.08 ± 0.34 a2.78 ± 0.38 b
a,b Values in the same row with different superscripts differ at p < 0.05. DM, dry matter; CP, crude protein; NDF, neutral detergent fiber; ADF, acid detergent fiber; ADL, acid detergent lignin.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fliegerova, K.O.; Podmirseg, S.M.; Vinzelj, J.; Grilli, D.J.; Kvasnová, S.; Schierová, D.; Sechovcová, H.; Mrázek, J.; Siddi, G.; Arenas, G.N.; et al. The Effect of a High-Grain Diet on the Rumen Microbiome of Goats with a Special Focus on Anaerobic Fungi. Microorganisms 2021, 9, 157. https://doi.org/10.3390/microorganisms9010157

AMA Style

Fliegerova KO, Podmirseg SM, Vinzelj J, Grilli DJ, Kvasnová S, Schierová D, Sechovcová H, Mrázek J, Siddi G, Arenas GN, et al. The Effect of a High-Grain Diet on the Rumen Microbiome of Goats with a Special Focus on Anaerobic Fungi. Microorganisms. 2021; 9(1):157. https://doi.org/10.3390/microorganisms9010157

Chicago/Turabian Style

Fliegerova, Katerina O., Sabine M. Podmirseg, Julia Vinzelj, Diego J. Grilli, Simona Kvasnová, Dagmar Schierová, Hana Sechovcová, Jakub Mrázek, Giuliana Siddi, Graciela N. Arenas, and et al. 2021. "The Effect of a High-Grain Diet on the Rumen Microbiome of Goats with a Special Focus on Anaerobic Fungi" Microorganisms 9, no. 1: 157. https://doi.org/10.3390/microorganisms9010157

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop