2.1 Harvesting and Silage
Alfalfa (Medicago sativa L.) from the Hunan Province, China, was harvested at the initial flowering stage, leaving stubble of 5 cm with a forage harvester. The harvested alfalfa was chopped to about 2.0 cm in length using a straw chopper after it was wilted until DM to 65 % on the plot. After wilting, alfalfa was silage with/without urea in a special plastic storage tank (plastic drums of 5 dm3 each, diameter 17.5 cm, and height 28.3 cm). Urea was prepared on the day of silage, applied in powder at 0 % (AL), 0.5 % (AU1) and 1 % (AU2) of the fresh weight, respectively. Before being packed into the tanks, urea was spread on the chopped forage on the floor and thoroughly mixed. Approximately 2 kg of chopped fresh alfalfa was immediately packed into a special plastic storage tank, compacted, sealed with plastic wraps (internally and externally), and stored at 20℃-25℃ until sampled. There were 36 silage tanks in each treatment.
2.2 Chemical Analysis
After silage for 0, 15, 30, and 60 d, the silage tanks were opened and samples were taken for chemical analysis with three replicates for each treatment. In total, 36 samples were collected.
Silage samples (20 g each) were diluted with 180 mL of distilled water, homogenized for 1 min in an Erlenmeyer flask to obtain silage alfalfa extracts, filtered through two layers of cheesecloth, and immediately analyzed for pH using a digital pH meter (Model PH S210-Basic, METTLER TOLEDO Solutions, Shanghai Precision & Scientific Instrument Co., LTD, China).
An aliquot of the silage alfalfa extract was centrifuged at 12,000 × g for 10 min at 4 °C (Weatherburn, 1967). Taking 15 μL supernatant added to 300 μL phenol-Nitpuna solution in a 1.5 ml centrifuge tube, mix well, then add 300μL base sodium hypochlorite solution to mix well again and incubated at 37 ℃ for 30 min. The supernatant was analyzed for ammonia-N using a microplate reader (TACAN, Switzerland ) set at 625 nm. The percentage of total nitrogen was calculated using Eq. 1:
Five hundred grams silage were dried in the air-blower-driver drying closets at 65 °C for 48 h and used for further chemical analysis. Dry matter (DM), crude protein (CP), crude fiber (CF), neutral detergent fiber (NDF), acid detergent fiber (ADF),in the silage were analyzed according to the methods of the Association of Official Analytical Chemists (AOAC) (William and George, 2007).
2.3 Total Deoxyribonucleic Acid Extraction, Polymerase Chain Reaction Amplification and 16S Ribosomal Ribonucleic Acid Sequencing
Total DNA was extracted using the E.Z.N.A. ®Stool DNA Kit (D4015, Omega, Inc., Norcross, GA, USA) according to manufacturer’s instructions. After DNA isolation, the purity and concentration of genomic DNA were measured with a spectrophotometer (Nanodrop 115 1000, Thermo Fisher Scientific, Waltham, MA, USA) and its integrity was verified by agarose (0.8 %) gel electrophoresis. After being eluted in Elution buffer, the total DNA was stored at -80 °C until measurement in the polymerase chain reaction (PCR).
Polymerase chain reaction amplification was performed in a total volume of 25 μL of reaction mixture containing 25 ng of template DNA, 12.5 μL PCR Premix, 2.5 μL of each primer, and PCR-grade water to adjust the volume. The V3–V4 region of the bacteria 16S rRNA gene was amplified with primers 338F (5'-ACTCCTACGGGAGGCAGCAG-3') and 806R (5'-GGACTACHVGGGTWTCTAAT-3') by thermocycler PCR system (GeneAmp 9700, ABI Inc., Edison, NJ, USA). The PCR conditions for the amplification of the prokaryotic 16S fragments consisted of an initial denaturation at 98 °C for 5 min; 35 cycles of 10 s for denaturation at 98 °C, 30 s for annealing at 54°C, 45 s for extension at 72 °C; and a final extension at 72 °C for 10 min. The PCR products were confirmed through 2 % agarose gel electrophoresis. Then, the PCR products were extracted from the agarose gel, purified by AMPure XT beads (Beckman Coulter Genomics, Danvers, MA, United States), and quantified by Qubit (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s protocol. The amplicon pools were used for sequencing. The size of amplicon library was assessed on the Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA), while the number was evaluated using the Library Quantification for Illumina (Kapa Biosciences, Woburn, MA, USA). Finally, the Phix Control library (V3, Illumina, Inc., San Diego, CA, USA) was unified with the amplicon library (expected at 30 %).
2.4 Sequence Processing and Bioinformatics Analysis
The amplicon sequences were conducted using an Illumina MiSeq platform according to the manufacturer's recommendations, provided by LC-Bio (LC-Bio, Hangzhou, Zhejiang, China). All paired-end reads were merged using FLASH (version 2.1.11, Magoč and Salzberg, 2011). Quality filtering on the raw tags was performed according to the FQTRIM (V 0.94) to obtain the high-quality clean tags. Low-quality sequences were filtered out using Vsearch software (v2.3.4, Rognes et al. 2016). The same operational taxonomic units (OTUs) were assigned by Vsearch (v2.3.4), with ≥ 97 % sequence similarity. The Ribosomal Database Project (RDP) classifier was applied to perform sequence-level taxonomic assignment. Alpha diversity was applied in analyzing the complexity of species diversity for a sample through five indices, including Chao1, observed species, Good’s coverage, Shannon, and Simpson’s diversity induces, and these were calculated in our samples using QIIME (Version 1.8.0). Beta diversity analysis was used to evaluate the differences in the samples in terms of species complexity. Beta diversity was calculated by principal coordinate’s analysis (PCoA) and cluster analysis using the QIIME software (Version 1.8.0).
2.5 Statistical Analysis
The experiment was a completely randomized design with three treatments, three replicates, and four-time points per treatment. The model included the treatment (T) effect, silage time (ST) effect, and the treatment × silage time (T*ST) reaction effect. The nutrient data were analyzed by the general linear model (GLM) procedure of SAS (version 8.0; SAS Institute, Inc., Cary, NC, USA), and the means were separated. Tukey–Kramer multiple comparison test was used for the different sample means. Differences among means were considered significant at P < 0.05 and extremely significant at P < 0.01, whereas 0.05 ≤ P < 0.1 was considered as a tendency of significant difference. The SAS and R (v3.0.3) software were used for analyzing the 16S rRNA gene data. Continuous variables such as number of reads, the relative abundance of bacteria, and species richness and diversity were analyzed using the GLM procedure of SAS. After statistical comparison of the taxa, Benjamini–Hochberg correction was used to control the false discovery rate by the package “p.adjust” in R; P < 0.05 was considered as statistically significant.