INTRODUCTION
The rapid growth of meat and milk production in China requires development of an efficient silage industry. Maize straw is an important feed resource for ruminants. In the north of China, such as Inner Mongolia, it is not only a main base for milk, beef and mutton production, but also a main production base of maize. Many of silage maize variety are cultivated for ruminants feed.
Ensiled silage has the benefit of preventing loss in nutrition value and can keep well throughout a year or more. The essence of the silage is that lactic acid bacteria (LAB) multiply in the anaerobic environment and produce lactate which drops the pH to below 4.0. The low pH in turn prevents or restrains the growth of other infectious microbes.
Microbial succession in food fermentation or silage has been previously studied. For instance, Zheng et al [
1] examined microbial community dynamics of ensiling alfalfa
Medicago sativa and found that
Lactobacillus dominated the whole process of ensiling like the dynamics patterns of ensiling reported by McGarvey et al [
2]. Li et al [
3] reported the bacterial community succession in
Vicia faba fermentation persisted for 77 days using MiSeq analysis base 16S rRNA genes. Yang et al [
4] used 16S rRNA sequences to report the bacterial community diversity and dynamics for 15 days fermentation of douchi sauce. For maize ensiling, Lin et al [
5] reported dominant establishment of
lactobacillus in early stages of maize ensiling (1 to 7 days). Subsequent studies by Zhou et al [
6] also demonstrated the dominant position of
lactobacillus after 60 days of ensiling in whole-plant corn silage.
However, research on maize silage production has mainly used culture-dependent methods or denaturing gradient gel electrophoresis (DGGE) [
5–
7], and lack using new techniques such as high-through sequence which can give a full-scale insight into the microbial community without culture-dependent process. In this study, we used Illumina MiSeq to describe the bacterial community dynamics during ensiling of whole-plant maize; this will enlighten the understanding of the dynamic process of bacterial community during the fermentation of maize.
MATERIALS AND METHODS
Dealing of samples
Raw maize samples were collected from a local meadow named Chajintai in Kerzuozhong County (E122° 15′, N43° 38′), Inner Mongolia, China, on Oct 8, 2016. The fresh maize was harvested by combine-harvester 20 cm from ground and chopped into 1 to 2 cm in lengths and immediately transported to the laboratory for processing. Three kilograms of the chopped maize material was compressed and sealed in 18 polyethylene silo bags (Hiryu KN type, 400×600 mm, Asahikasei, Tokyo, Japan) using a vacuum sealer (BH 950, Matsushita, Tokyo, Japan) for air removal. The silo bags were stored at room temperature (25°C). Samples were taken from three bags per timing at 5, 10, 20, 40, 60, and 63 days (the bags were opened for exposure to air at day 60 and lasted for 3 days), respectively, and divided into two portions. One part was sampled in 50 mL cryogenic vials and stored in −80°C for 16S rRNA gene amplicons sequencing, and another was stored in −20°C for biochemical analysis.
Chemical analyses
The dry matter (DM) content of silage samples was determined by drying 100 g of samples collected on each sampling day at 65°C for 48 h in a dry oven [
8]. The pH value was measured as follows: 10 g of silage was mixed with 90 mL of deionized water and shaken for 5 min every other 50 min, after 3 h, the pH value of filtrate form silage was measured by electrode pH meter (S20K, Mettler Toledo, Greifensee, Switzerland). Concentrations of lactic acid and acetic acid were measured in high performance liquid chromatography (1100, Agilent, Palo Alto, CA, USA) fitted with a UV detector (210 nm) and a column (ICSep COREGEL-87H, Transgenomic, Inc., New Haven, CT, USA). The mobile phase was 0.005 M H
2SO
4 at a flow rate of 0.6 mL/min at 55°C [
9]. Fifteen grams of fresh forage or silage at each sampling day was combined with 135 g of de-ionized H
2O, blended for 30 s and strained through four layers of cheesecloth. The liquid fraction was sampled and analyzed for and ammonia N (NH
3-N) as described by Zahiroddini et al [
10].
DNA extraction of the microbial community
The fermented maize straw was sampled (20 g for each bag) and quickly frozen at −80°C for DNA extraction. Silage (20 g) was homogenized with distilled water (80 mL), then filtered through two layers of medical gauze and centrifuged for 15 min at 10,000 g at 4°C.
Triplicate samples of centrifuge production were dissolved and blended together. DNA was extracted using the FastDNA SPIN Kit for Soil (MP Biomedicals, Santa Ana, CA, USA), according to the manufacturer’s instructions [
11]. The DNA quality was checked by 1% agarose gel electrophoresis and spectrophotometer (NanoDrop2000, Thermo Scientific, Waltham, MA, USA, optical density at 260 nm/280 nm ratio).
Polymerase chain reaction amplification of target region
Primers 338F (ACTCCTACGGGAGGCAGCA) and 806R (GGACTACHVGGGTWTCTAAT) [
12] both with 12 bps barcodes were used to amplify the V3–4 hyper-variable region of the bacterial 16S rRNA gene. The sequences were assigned to their respective samples according to their barcodes of primers. The polymerase chain reaction (PCR) amplification of bacterial 16S rRNA V3–4 region was performed as described in TransStart FastPfu system (Transgen Biotech Company, Beijing, China) [
13]. The PCR cycle conditions were the following: an initial denaturation at 95°C for 3 min, 27 cycles of 95°C for 30 s, 55°C for 30 s, and 72°C for 45 s, and a final extension at 72°C for 10 min. Each 20 μL reaction mixture consisted of 4.0 μL of 5×FastPfu Buffer, 2 μL of deoxynucleoside triphosphate mix (2.5 mM each), 0.4 μL of FastPfu DNA Polymerase, 10 ng of template DNA, 0.8 μL of Forward Primer 338F (5 μM), 0.8 μL of Forward Primer 806R (5 μM), 0.2 μL of bovine serum albumin, and others of double-distilled H
2O.
Library preparation and sequencing
NEB Next Ultra DNA Library Prep Kit for Illumina (New England Biolabs Inc., Ipswich, MA, USA) was used to generate sequencing libraries according to the manufacturer’s instructions. The library quality was assessed by spectrophotometry and then sequenced on Illumina Miseq PE300 platform (Illumina Corporation, San Diego, CA, USA) at the Shanghai Majorbio Bio-Pharm Technology co., Ltd. (Shanghai, China).
Bioinformatics analysis of sequencing data
Statistical analysis
All results are reported as the mean±standard error of three replicate groups, and data were subjected to one-way analysis of variance. When there were significant differences (p<0.05), the group means were further compared with Duncan’s multiple range tests. All statistical analyses were performed using SPSS 17.0 (SPSS, Chicago, IL, USA).
DISCUSSION
Compared to extensive work with culture-dependent methods and LAB-focusing [
19–
21], advanced molecular techniques could better demonstrate the bacterial diversity of maize silage. The present study examined the bacterial diversity of whole-plant maize silage during fermentation (for 60 days) and subsequent aerobic exposure (from d 60 to d 63) by Illumina Miseq sequence platform.
As known from previous studies, 5 to 7 days is a good time to check the microbial composition and biochemical index [
5,
6]. This was verified by the pH value which reached below 4.0 and the rapid growth of
Lactobacilus in silage after 5 days of fermentation. For 60 days of fermentation, the pH remained below 4.0 stably and even after opening the package for aerobic exposure. The low pH value insured the anaerobic stability and eventually contributed to the favorable fermentation of whole-plant maize. Lactic acid and acetic acid, produced by LAB, were the reason for pH decline [
22]. In the present, the lactic acid increased significantly after 20 days of ensiling, while the significant rise of acetic acid occurred at later stage and aerobic exposure stage. Both acids have benefits for the control of spoilage yeasts or other pathogenic microorganisms [
23]. During the ensiling, the content of DM dropped slightly. The content of lactic acid and acetic acid was at a high level, and the ammoniac nitrogen content was in limit ranges. All these indexes showed the excellent fermentation characteristics of the silage.
High-through sequencing has the benefit of finding more bacterial species in silage. We found 295 species belonging to 207 genera in the silages. Except 5 genera of LAB, more bacteria belong to non-LAB. It’s worth noting that previous research mainly focused on LAB, or pathogenic bacteria [
5,
6], and hardly on others which may play an equally important role during the ensiling process.
It has been well documented that the high content of LAB, such as
Lactobacilus, is an assurance of rapid decline in pH, which is a key factor ensuring fermentation quality during the ensiling process [
23,
24]. In the study, the species of
Lactobacillus dominated the bacterial communities during the whole process of ensiling, even after aerobic exposure. The increasing of
Lactobacillus abundance (from 49.57% in d 5 to 63.20% in d 60) displayed a positive correlation with the increasing of lactic acid (from 48.37 mg/g in d 5 to 96.9 mg/g in d 60) and acetic acid (from 4.79 mg/g to 16.87 mg/g), indicating that the increase of lactic acid and acetic acid was mainly due to
Lactobacillus. Some other studies also reported the predominance of
Lactobacillus in silages [
1,
5,
25].
Lactobacillus species, such as
L. plantarum,
L. brevis, and
L. buchneri, are extensively used as inoculants in silage [
7,
21,
24,
26]. In our study, the sequences in the genus
Lactobacillus were mainly assigned to
L. plantarum,
L. brevis and
L. crustorum.
L. plantarum, a homo-fermentative bacteria, produces only lactic acid, and
L. brevis, hetero-fermentative LAB, produces both lactic and acetic acid. Their stable presence contributes to the rising of lactic and acetic acid during fermentations. Above all,
L. crustorum was also a main constituted species of
Lactobacillus during ensiling and played an important role for good fermentation quality of ensilage. Although
L. crustorum has been widely reported in food fermentation, little has been reported in silage. As reported by Arasu et al [
27],
L. crustorum was isolated from silage, and showed high growth rates and organic acid (lactic acid, acetic acid) production capability. So,
L. crustorum may be a potential probiotic for use silage.
In addition to
Lactobacillus, several other acids-producing bacteria appear in the early stages (d 5 and d 10), and decrease or disappear in subsequent periods, namely
Pediococcus,
Lactococcus,
Weissella, and
Leuconostoc. This situation was demonstrated in several studies [
5,
25].
Pediococcus,
Lactococcus,
Weissella, and
Leuconostoc, hetero-fermentative bacteria, can produce lactic and acetic acid. The low pH caused by these acids is important for initiating the fermentation of silage. These species have a role of creating a suitable environment for the further growth of
lactobacilli [
21,
24]. The LAB profiles of our study are like the DGGE analysis of maize silage from Inner Mongolia which reported the dominant of
Lactobacillus and the appearace of
Lactococcus,
Leuconostoc, and
Weissella after half a year of ensiling from epiphytic bacteria [
28], indicating the richness of epiphytic
Lactobacillus in this region.
Gluconobacter,
Acetobacter, and
Flavobacterium were found in minor abundance during the early stages of ensilage.
Gluconobacter and
Acetobacter were reported by Gharechahi et al [
25] in natural fermented maize silages.
Acetobacter, an acetic acid producing and nitrogen-fixing bacteria, can be found in many plants [
29], and may contributes to the pH decline at the early stage of silages.
Flavobacterium was found in the silages and developed for inoculants because of its ability to produce cellulase [
30]. Cellulase can degrade cellulose fiber to monosaccharides, such as hexose and pentose. From this point, water soluble carbohydrates or sugars were not completely illustrated the substrate available for LAB [
22].
The genus
Klebsiella, an epiphytic microbiota, sub-dominated the bacterial community during the entire process of ensiling. Like our findings, other studies have found that
Klebsiella is present throughout the fermentation process of corn silage [
7]. However, according to the report of Ni et al [
31],
Klebsiella was only found in pre-ensiled corn silage, not ensiling process. Most of
Klebsiella species, belonging to Enterobacteriaceae family, are found to produce 1, 3-propanediol during fermentation such as
Klebsiella pneumoniae ATCC 8724 [
32], but the propanediol and butyric acid produced by
Klebsiella was not previously considered to be benefit indicators in silage [
22]. In addition to the
Klebsiella, other bacteria, such as
Pantoea,
Rahnella,
Sphingomonas,
Serratia,
Chryseobacterium,
Sphingobacterium,
Stenotrophomonas,
Acinetobacter, and
Enterobacter, were evenly distributed throughout the fermentation process of silages in lower amounts. In fact, the above bacteria were found in various maize and grass silages [
25,
30,
31,
33], but our study covered all these species for the first time. This may be related to our study methods, which cover a greater depth and breadth than other methods, or to different epiphytic bacteria of different silages. Among them,
Pantoea,
Rahnella,
Serratia,
Enterobacter, and
Enterobacter belong to the family of Enterobacteriaceae which produces not only organic acids but also ethanol. Such as
Enterobacter and
Rahnella, they were found to produce high level of alcohol if they were prevalent on pre-ensiled crops [
34].
Aerobic exposure had a noticeable impact on the structure of bacterial communities. In the study, the remarkable changes were the declining of
Lactobacillus (from 63.20% to 45.03%) and the increasing of
Sporolactobacillus (from 5.32% to 28.15%). The slightly reduction in the concentration of lactic acid may be due to the decrease of
Lactobacillus. The increases in the acetic acid contents may be due to the high proportion of
L. crustorum and
L. brevis, both of which produce acetic acid. As a facultative anaerobic,
Sporolactobacillus, a spore-forming, gram-positive bacterium, can produce D-lactic acids by homo-fermentation of hexoses [
35]. Although the presence of
Sporolactobacillus in corn silage and its antifungal activities have been reported by Kharazian et al [
36], there is still little coverage of the genus in silage. In our study,
Sporolactobacillus may be contributed to the maintaining of low pH during the aerobic exposure stage, since the silage in this period was at risk of pH rise due to the decreasing of LAB
Lactobacillus.
Spoilage and pathogenic microorganisms may be present in silage during fermentation process, such as
Clostridium, Cyanobacteria, and
Listeria [
37].
Listeria was not found in our study, but Cyanobacteria appeared in most of the samples (d 5, d 10, d 20, d 60).
Clostridium appeared in aerobic exposure period. The species of
Clostridium may be responsible for the production of butyric acid, which can reduce the content of lactic acid, as well as the production of biogenic amines by protein fermentation [
7,
22]. It may be related to
Clostridium that biogenic amines increased, and lactic acid slightly decreased in aerobic exposure samples.
The term succession is defined as the biological changes of an ecosystem in a predictable way of species composition shifts [
38]. The predictable dynamics of bacterial diversity were shown in our present reflected by the dominance of
Lactobacillus over the whole time, indicating the good participation of this epiphytic bacterium. This agreed with other studies [
5,
25], which reported the predominance of
Lactobacillus in maize silage.