- 1Division of Applied Life Science (BK21Four), Institute of Agriculture and Life Science, Gyeongsang National University, Jinju, South Korea
- 2Biotechnology Division, Bangladesh Livestock Research Institute, Savar, Bangladesh
- 3Dairy Science Division, National Institute of Animal Science, Cheonan, South Korea
- 4Faculty of Animal Science, Universitas Gadjah Mada, Yogyakarta, Indonesia
This study estimated the effect of essential fatty acid (FA) supplementation on fermentation indices, greenhouse gases, microbes, and FA profiles in the rumen. The treatments used pure FAs consisting of C18:2n-6 FA (LA), C18:3n-3 FA (LNA), or a mixture of these FAs at 1:1 ratio (Combo). In vitro rumen incubation was performed in 50 mL glass serum bottles containing 2 mg of pure FAs, 15 mL of rumen buffer (rumen fluid+anaerobe culture medium = 1:2), and 150 mg of synthetic diet (411 g cellulose, 411 g starch, and 178 g casein/kg dry matter) at 39°C for 8 h with five replications and three blanks. In rumen fermentation indices, LA exhibited highest (P < 0.05) ammonia-N and total gas volume after 8 h of incubation. Furthermore, LA presented lower (P < 0.05) pH with higher (P < 0.05) total volatile fatty acid (P = 0.034) than Combo, while LNA was not different compared with those in the other treatments. Additionally, Combo produced highest (P < 0.05) CO2 with lowest (P < 0.05) CH4. In the early hours of incubation, LA improved (P < 0.005) Fibrobacter succinogenes and Ruminococcus flavefaciens, while LNA improved (P < 0.005) Ruminococcus albus. After 8 h of incubation, LNA had lower (P < 0.05) methanogenic archaea than LA and Combo but had higher (P < 0.05) rumen ciliates than LA. R. albus was higher (P < 0.05) in LA than in LNA and Combo. It was observed that the rate of biohydrogenation of n-6 and n-3 FAs was comparatively lowest (P < 0.05) in Combo, characterized by higher C18:2n-6 and/or C18:3n-3 FA and polyunsaturated FA (PUFA) concentrations with lower (P < 0.05) concentrations of C18:0 and saturated FA and the ratio of saturated FAs to PUFAs. Therefore, this study concluded that dietary C18:2n-6 could improve populations of fibrolytic bacteria and rumen fermentation indices, but dietary mixture of pure C18:2n-6 and C18:3n-3 is recommended because it is effective in reducing enteric methane emissions and resisting biohydrogenation in the rumen with less effect on rumen microbes.
Introduction
Oil-rich supplements have purposes to manipulate rumen ecosystem through a number of mechanisms, including reductions in organic matter fermentation, rumen ciliate numbers, methanogenic activity, and the use of hydrogen for biohydrogenation (Beauchemin et al., 2009). In addition, the rumen fermentation pattern is significantly influenced by the type of fat (Machmüller et al., 1998). Essential fatty acids (FAs) consisting of linolenic acid (C18:3n-3) and linoleic acid (C18:2n-6) are polyunsaturated FAs (PUFAs) that are required for animals because they are not synthesized in the body. In human, the essential FAs present beneficial effects for human health, such as decreasing cardiovascular morbidity, hypertension, and diabetes mellitus. The C18:3n-3 and C18:2n-6 in linseed and sunflower oils, respectively, have been reported to depress rumen methanogenesis (Machmüller et al., 1998, 2000; Martin et al., 2006; Beauchemin et al., 2009). Methane (CH4) emissions by ruminants produce around 5.55% of total global greenhouse gas emission from human activity (FAO, 2016). Additionally, CH4 represents a significant energy loss of ruminants in the range of 2–12% of gross energy intake (Johnson and Johnson, 1995; FAO, 2016).
In ruminants, dietary essential FAs undergo biohydrogenation to saturated fatty acids (SFAs) in the rumen, which can be an alternative approach to reduce the amount of free hydrogen for methanogenesis (Machmüller et al., 1998, 2000; Martin et al., 2006; Beauchemin et al., 2009). Additionally, dietary essential FAs have been repeatedly demonstrated to have adverse effects on rumen microbes such as protozoa and methanogenic archaea (Doreau and Ferlay, 1995) and could depress methanogenesis (Fievez et al., 2003). As a negative effect of dietary unsaturated fatty acids (UFAs), Zhang et al. (2008) reported that oil-rich supplement containing C18:2n-6 and C18:3n-3 could inhibit the growth of fibrolytic fungi, Fibrobacter succinogenes, and Ruminococcus flavefaciens in the rumen but suppressed the methanogenesis. On the other study, dietary soybean oil and linseed oil potentially decrease the populations of Butyrivibrio fibrisolvens, Ruminococcus albus, and F. succinogenes, and also decrease total volatile fatty acid (VFA) production in the rumen (Yang et al., 2009). As a positive effect, Jouany et al. (2008) indicated that C18:3n-3 FA had a greater effect on reducing CH4 production than did C18:2n-6 by in vitro measurement. Otherwise, dietary UFAs usually undergo extensive biohydrogenation by rumen microbes (Beam et al., 2000), which is the major challenge for achieving the targeted effects of these FAs through dietary supplementation. The extent of biohydrogenation of long-chain PUFAs varies largely depending on the source. A disappearance of C18:2n-6 of approximately 90–98% occurred after 9 h of incubation when it was supplied in a pure form (Honkanen et al., 2012). On the other hand, another study (Carriquiry et al., 2008) observed that approximately 40–65% of the C18:2n-6 fraction could remain after 36 h of incubation when commercial mixed fat sources were used. These results mean that the extent of biohydrogenation of PUFAs in the rumen may vary depending on the form of fat, which may further affect the fermentation pattern. Most of the studies regarding FA supplementation in ruminants were conducted by oil or oil seed sources of FAs (Machmüller et al., 1998, 2000; Martin et al., 2006; Beauchemin et al., 2009). However, the effect of single C18 n-6 or n-3 FA supplementation on rumen ecosystem has been scarcely studied, which might show different effects from C18 n-6 FA and C18 n-3 FA on rumen fermentation, methanogenesis, and microbial population. Especially, the effects of single C18 n-6 or n-3 FA supplementation on population of fibrolytic bacteria were in limited information. Therefore, this study was conducted to estimate the effect of pure C18:2n-6, C18:3n-3, or their mixture on rumen microbes with its fermentation indices, greenhouse gases, and FA profiles.
Materials and Methods
In vitro Incubation
Animal care and handling in the present study was approved by the Animal Ethics Committee of Gyeongsang National University, South Korea (GNU-191011-E0050). Rumen fluid was collected from two rumen-cannulated Hanwoo heifers (Average body weight = 432.63 kg) fed rice straw (CP = 5.40%, NDF = 63.85%) and concentrate (CP = 12.51%, NDF = 47.51%) mixture at a 8:2 ratio of dry matter (DM) weight. The concentrate mainly consisted of corn meal, soybean meal, soybean hull, and cotton seed pellet. This fluid was collected before morning feeding, filtered with double cheese cloths, stirred to obtain a uniform fluid, and then mixed with anaerobic culture medium at a 1:2 ratio to make rumen buffer. The detail protocol to prepare a rumen buffer was described by Adesogan et al. (2005). Carbon dioxide (CO2) gas was flushed into rumen buffer continuously to maintain anaerobic conditions (Carriquiry et al., 2008; Jayanegara et al., 2015). A fat-free synthetic diet was prepared using as a substrate, containing 411 g cellulose (Sigma, C6413), 411 g starch (Sigma, S4180), and 178 g casein (Sigma, C3400) per kg of DM (Carriquiry et al., 2008). The treatments consisted of pure C18:2n-6 FA (Sigma, L1376; LA), pure C18:3n-3 FA (Sigma, L2376; LNA), or a mixture of these FAs (Combo) at 1:1 ratio (w/w). Five hundred milligrams of LA and LNA was dissolved with 5 mL methanol (A 4524, Fisher Chemical), respectively. After then, 2 mL from each dissolved FAs were mixed for Combo. Before ruminal incubation, 20 μL of dissolved FA (2 mg of FA) was applied into 150 mg of synthetic diet. This synthetic diet applied FA was stored at room temperature to evaporate methanol for 12 h. In the present study, the application rate of each FA treatment was 1.3% DM, which demonstrate the general C18:2n-6 FA and C18:3n-3 FA concentrations in ruminant diet (Kim et al., 2007, 2016). Ruminal incubation was performed in the 50-mL glass serum bottle containing 15 mL of rumen buffer with synthetic diet applied FA at 39°C for 0, 1, 2, 4, and 8 h. Each FA treatments used five incubation bottles as replications along with three blanks for each hour. Thus, total 90 incubation bottles were used in the present study.
Sampling
In an assigned hour (0, 1, 2, 4, and 8 h, respectively), all bottles were withdrawn from the incubator, the gas pressure was quickly measured by a pressure transducer (Fisher Scientific, TraceableTM, Friendswood, TX, United States), and then, the bottles were placed into ice to stop microbial activity (Honkanen et al., 2012). All incubated rumen samples including blank were subsampled 1 mL for microbial quantification using PCR. The remain incubated rumen samples were transferred to 50 mL conical tube to separate sample residue and supernatant of incubated rumen buffer through centrifugation at 2,568 × g for 15 min (Supra 21k, Hanil Electric Corporation, Seoul, South Korea, with rotor A50S-6C No.6). The supernatant of incubated rumen buffer from centrifugation was subsampled 10 mL and stored for further analyses of rumen fermentation indices. For FA analysis, 2 mL rumen buffer was frozen at −70°C for 2 days and freeze-dried using Cascade Console Freeze Dry System (LABCONCO, FreeZone Plus 12 Liter, MO, United States) as the procedure to sample preparation according to previous studies (Paradhipta et al., 2020).
Analyses
Greenhouse Gas Emissions and Fermentation Indices
The total gas volume was calculated from the gas pressure (psi) according to Mauricio et al. (1999). Gas samples were collected in vacutainer tubes without additives for CH4 and CO2 analysis using a multigas analyzer (Yes Plus LGA, Critical Environment Technologies, Canada Inc., Delta, BC, Canada). The results of CH4 and CO2 were expressed as milligrams per gram of fresh weight. For rumen fermentation indices, the subsampled supernatant measured the pH using an electric pH meter (SevenEasy, Mettler Toledo, Greifensee, Switzerland). After then, the supernatant was centrifuged at 21,500 × g for 15 min for the analyses of ammonia-N (NH3-N) and VFA. The NH3-N concentration was measured by distillation of the sample in a Buchi apparatus (B-342, BÜCHI, Flawil, Switzerland) followed by titration with 0.1 N H2SO4 in a burette according to Chaney and Marbach (1962). The concentrations of VFAs were determined using an HPLC system (L-2200, Hitachi, Tokyo, Japan) fitted with a UV detector (L-2400; Hitachi) and a column (Metacarb 87H; Varian, Palo Alto, CA, United States) as described by Muck and Dickerson (1988).
Fatty Acid Profiles
Two-step methylation procedure was used for the preparation of FA methyl esters (FAMEs), which was described by Jenkins et al. (2001). One milligram of internal standard (C19:0) was added to the previously freeze-dried sample to calculate the total FA concentration. The FAs were esterified by adding 2 mL of sodium methoxide (Sigma Aldrich, St. Louis, MO, United States), followed by vortexing and then incubating in a 50°C water bath for 10 min. After cooling for 5 min, 3 mL of 5% methanolic HCl was added followed by vortexing. The tubes were then incubated at 80°C for 10 min in a water bath. After incubation, the solution was allowed to cool for 7 min, and 1 mL of hexane and 7.5 mL of K2CO3 were consecutively added. The tubes were shaken and centrifuged at 9,861 × g at 4°C. The upper layer containing FAMEs was transferred to a 2-mL top crimp vial (Agilent Technology, Cranberry Township, PA, United States) with Pasteur pipettes (Hilgenberg, Strauchgraben, Malsfeld, Germany). The FAMEs were analyzed using a gas chromatograph (450-GC, Varian, Palo Alto, CA, United States) equipped with an autosampler (CP-8400, Varian, Palo Alto, CA, United States), a flame ionization detector, and a Varian capillary column (CP-Sil 88, Palo Alto, CA, United States, 100 m × 0.25 mm × 0.2 μm). Hydrogen was the carrier gas. The injector and detector were maintained at 230°C. The oven temperature was initially set at 120°C for 1 min, increased by 5°C/min up to 190°C, held at 190°C for 30 min, increased again by 2°C/min up to 220°C, and held at 220°C for 40 min. The peak of the samples was identified, and concentrations were calculated based on the retention time and peak area of known standards.
DNA Extraction and Real-Time PCR
The DNA was extracted from incubated rumen sample by physical disruption using a mini bead-beater (BioSpec Products, Bartlesville, OK, United States) followed by isolation and purification using a commercial DNA extraction kit (QIAamp DNA mini kit, Qiagen, Germantown, MD, United States). A 700-μL aliquot of homogenously incubated rumen sample was transferred into a 1.5-mL microcentrifuge tube, and 180 μL of buffer animal tissue lysis (ATL) and 20 μL of proteinase K (supplied in the QIAamp DNA Mini Kit) were added to the tube followed by vortexing. The tubes were then incubated at 56°C for 12 h for cell lysis in a heating block (HB-48, Daihan Scientific Co., Ltd., Seoul, South Korea). Next, the DNA purification protocol was followed as described in the manufacturer’s instruction manual (Qiagen, Germantown, MD, United States). The DNA concentrations were measured by using a NanoDrop Spectrophotometer (ND-1000, United States). The primer information of general bacteria, F. succinogenes, R. flavefaciens, R. albus, rumen methanogenic archaea, and ciliates (Entodinium) are given in Table 1. Species-specific real-time qPCR was performed using a Bio-Rad C1000 TouchTM Thermal Cycler Real-Time PCR Detection System (CFX96TM Real-Time System, Bio-Rad Laboratories, Inc., Hercules, CA, United States), with fluorescence detection by SYBR Green Real-time Master Mix (TOYOBO Co., Ltd., Osaka, Japan). The values of the cycle threshold (Ct) after real-time PCR were used to determine the fold change of different microbial populations relative to the respective blanks without FA treatments. The abundance of these microbes was expressed by the following equation: relative quantification = 2–ΔCT(Target)–ΔCT(Blank), where Ct represents the threshold cycle. Rumen general bacteria were used as reference genes for internal controls according to Kim et al. (2012). All quantitative PCR mixtures (final volume of 20 μL) contained forward and reverse primers, SYBR Green Master Mix, DNA template, and sterilized distilled water. A negative control without template DNA was used in every qPCR assay for each primer.
Statistical Analyses
The experiment was a completely randomized design. All data were analyzed using the general linear model procedure of Statistical Analysis System (SAS) ver. 9.3 (SAS, 2000). Its model was Yij = μ + Ti + eij, where Yij = response variable, μ = overall mean, T = effect of treatment i, and eij = error effect. In addition, rumen pH, NH3-N, acetate, propionate, CH4, CO2, rumen microbes, and major 18-carbon FAs were analyzed using the PROC MIXED procedure of SAS to test the significant levels of supplementary treatment, incubation hour, and the interaction between supplementary treatment and incubation hour. Its model was Yijk = μ + αi + βj + (αβ)ij + eijk, where Yijk = response variable, μ = overall mean, αi = the effect of supplementary treatment, βj = the effect of incubation hour, (αβ)ij = the interaction effect between supplementary treatment and incubation hour, and eijk = error effect. Mean differences were tested for significance using Tukey’s test at P ≤ 0.05.
Results
Fatty Acid Profiles Before Incubation
The total FAs before incubation from LA, LNA, and Combo were 4.11, 3.87, and 3.77 mg/mL, respectively (Table 2). The concentrations of C18:2n-6 was highest (P < 0.001; 18.4 vs. 11.02 vs. 3.07%) in LA, followed by Combo and LNA, while the concentrations of C18:3n-3 was highest (P < 0.001; 6.78 vs. 3.72 vs. 0.26%) in LNA, followed by Combo and LA. The concentrations of C14:0, C14:1, C16:0, C16:1, C18:0, SFA, monounsaturated fatty acid (MUFA), and PUFA were higher (P < 0.05) in LNA than LA. The ratio of SFA to PUFA was higher (P = 0.005; 7.65 vs. 3.74 and 4.92) in LNA than LA and Combo. In addition, the ratio of n-6 to n-3 was higher (P < 0.001; 35.76 vs. 0.48 and 2.84) in LA than LNA and Combo.
Table 2. Effects of fatty acid supplementation on the fatty acid profiles of rumen buffer just before incubation.
Fermentation Indices and Greenhouse Gas Emissions
After 8 h of incubation, the pH was lower (P = 0.034; 6.52 vs. 6.65) in LA than in Combo, while the pH in LNA was not different compared with that in the other treatments (Table 3). The NH3-N concentration was higher (P = 0.002; 33.46 vs. 30.38 and 30.38 mg/100 mL) in LA than in LNA and Combo after 8 h of incubation. The total VFA concentration was higher (P = 0.039; 85.22 vs. 78.93 mmol/L) in LA than in Combo, while the concentration in LNA did not differ compared with that in the other treatments. Among the individual VFAs, the concentrations of acetate, propionate, butyrate, isovalerate, and valerate remained unaffected. Only the concentration of isobutyrate was found to be higher (P < 0.001; 0.93 vs. 0.83 and 0.74%) in LA than in LNA and Combo. The total gas volume in LA was higher (P < 0.001; 35.43 vs. 31.13 and 30.29 mL/g) than that in LNA and Combo. The CH4 emissions were highest (P = 0.001; 55.50 vs. 46.60 vs. 40.20 mg/g) in LA, followed by LNA and then Combo. The CO2 gas concentration was higher in Combo than in LA and LNA (P = 0.042; 2.05 vs. 1.61 vs. 1.80 mg/g).
Table 3. Effects of fatty acid supplementation on the fermentation indices and greenhouse gas emissions of rumen buffer after 8 h of incubation.
The rumen pH was decreased by increasing incubation hour, but supplementary treatment had no effect from 0 to 4 h of incubation (Figure 1A). The NH3-N concentration was increased by incubation hour in all supplementary treatments, which Combo presented higher concentration than LA at 0 h (P = 0.043; 20.44 vs. 18.90 mg/100 mL) and 1 h (P = 0.032; 23.94 vs. 22.40 mg/100 mL) of incubation (Figure 1B). The acetate and propionate concentrations were affected by the supplementary treatment at 1 and 2 h of incubation, but the patterns of concentration change over hour were cubical pattern (Figures 1C,D). Furthermore, LA resulted in highest concentrations of acetate (P = 0.001; 45.41 vs. 41.01 and 41.49%) and propionate (P < 0.001; 25.14 vs. 21.86 and 21.92%) at 1 h, then had lowest (P < 0.01) concentrations of acetate (P < 0.001; 44.22 vs. 46.86 and 47.00%) and propionate (P = 0.001; 22.74 vs. 23.61 and 24.42%) at 2 h of incubation. There were interaction effects between the supplementary treatment and hour on rumen pH (P = 0.041) and NH3-N (P < 0.001), acetate (P < 0.001), and propionate (P < 0.001) concentrations, which might be a reason for the cubical patterns of these results during incubation. The concentration of CH4 was increased by increasing incubation hour in all supplementary treatments (Figure 2A). LA had highest (P < 0.01; 3.68 vs. 2.68 and 2.60 mg/g) CH4 concentration at 2 h of incubation, while LA and Combo presented higher (P < 0.05; 4.50 and 4.80 vs. 3.85 mg/g) concentration than LNA at 4 h of incubation. Similar to CH4, concentration of CO2 was also increased by increasing incubation hour, which Combo and LA presented higher (P < 0.05; 0.83 and 0.92 vs. 0.76 mg/g) concentration than LNA at 1 h or incubation (Figure 2B). Additionally, there were interaction effects between the supplementary treatment and hour on the CH4 (P < 0.001) and CO2 (P = 0.005) concentrations.
Figure 1. Changes in rumen pH (A) and concentrations of NH3-N (B), acetate (C), and propionate (D) during rumen incubation with FAs for 8 h. Diet with C18:2n-6, LA (filled circles); diet with C18:3n-3, LNA (filled triangles); and diet with the mixture of C18:2n-6 and C18:3n-3 at a ratio of 1:1 (filled squares). The respective significance levels of the supplementary treatment, incubation period, supplementary treatment × incubation period, and SEM for pH, NH3-N, acetate, and propionate are P < 0.001, P = 0.054, P = 0.041, and SEM = 0.057; P < 0.001, P = 0.115, P < 0.001, and SEM = 1.064; P < 0.001, P = 0.071, P < 0.001, and SEM = 0.715; and P < 0.001, P = 0.037, P ≤ 0.001, and SEM = 0.487. Different symbols at the same hour differ significantly. ∗P < 0.05; ∗∗P < 0.01.
Figure 2. Changes in CH4 (A) and CO2 (B) emissions during rumen incubation with FAs for 8 h. Diet with C18:2n-6, LA (filled circles); diet with C18:3n-3, LNA (filled triangles); and diet with the mixture of C18:2n-6 and C18:3n-3 at a ratio of 1:1 (filled squares). The respective significance levels of the supplementary treatment, incubation period, supplementary treatment × incubation period, and SEM for CH4 and CO2 are P < 0.001, P = 0.001, P < 0.001, and SEM = 3.021 and P < 0.001, P < 0.001, P = 0.005, and SEM = 91.69. Different symbols at the same hour differ significantly. ∗P < 0.05; ∗∗P < 0.01.
Microbial Populations
It was observed that methanogenic archaea were less abundant (P = 0.03; 0.74 vs. 0.96 and 0.96) in LNA than in LA and Combo (Table 4). Rumen ciliates were less abundant (P = 0.03; 0.20 vs. 0.41 of fold change) in LA than in LNA, respectively, while the population in Combo had no difference compared with that in the other treatments. Among fibrolytic bacteria, R. albus was observed to be more abundant (P = 0.005; 0.17 vs. 0.09 and 0.10 of fold change) in LA than in LNA and Combo. F. succinogenes and R. flavefaciens remained unaffected by the treatments at 8 h of incubation.
Table 4. Effects of fatty acid supplementation on microbial populations of rumen buffer after 8 h of incubation (fold change unit).
It was observed that the population of methanogenic archaea was not affected by supplementary treatment during 4 h of incubation (Figure 3A). The population of rumen ciliates dropped drastically from 0 to 2 h of incubation in all supplementary treatments (Figure 3B). The population of F. succinogenes seemed to decrease by increasing incubation hour (Figure 3C). Furthermore, LA presented higher F. succinogenes population than LNA and Combo at 2 h (P = 0.003; 1.31 vs. 0.88 and 0.80) and 4 h (P = 0.001; 1.01 vs. 0.72 and 0.86) of incubation. In general, R. flavefaciens population increased during 4 h of incubation in all supplementary treatments, then decreased after that (Figure 3D). The R. flavefaciens population was higher (P = 0.007; 1.08 vs. 0.77 and 0.63) in LA and presented higher population than in LNA and Combo at 2 h of incubation. The population of R. albus was higher in LNA than in LA and Combo at 1 h (P = 0.015; 0.87 vs. 0.71 and 0.57) and 4 h (P = 0.008; 0.84 vs. 0.52 and 042) of incubation (Figure 3E). Generally, all supplementary treatments seemed to decrease the population of R. albus after 8 h of incubation. An interaction between the supplementary treatment and hour was reported for rumen ciliates (P = 0.019), F. succinogenes (P = 0.003), and R. albus (P < 0.001), which caused cubical patterns of those microbes during incubation.
Figure 3. Fold change in methanogenic archaea (A), rumen ciliate (B), F. succinogenes (C), R. flavefaciens (D), and R. albus (E) compared with that of the blank during rumen incubation with FAs for 8 h. Diet with C18:2n-6, LA (filled circles); diet with C18:3n-3, LNA (filled triangles); and diet with the mixture of C18:2n-6 and C18:3n-3 at a ratio of 1:1 (filled squares). The respective significance levels of the supplementary treatment, incubation period, supplementary treatment × incubation period, and SEM for methanogenic archaea, rumen ciliates, F. succinogenes, R. flavefaciens, and R. albus are P = 0.255, P = 0.157, P = 0.492, and SEM = 2.046; P < 0.001, P = 0.458, P = 0.019, and SEM = 0.344; P < 0.001, P < 0.001, P = 0.003, and SEM = 0.124; P < 0.001, P = 0.032, P = 0.301, and SEM = 0.168; and P < 0.001; P = 0.001, P < 0.001, and SEM = 0.099. Different symbols at the same hour differ significantly. ∗P < 0.05; ∗∗P < 0.01.
Fatty Acid Profiles After Incubation
The concentrations of total FAs and MUFAs were found to be unaffected (P > 0.05) by the treatments after 8 h of incubation (Table 5). The concentrations of C14:1 (P = 0.001; 1.45 and 1.44 vs. 1.30%) and C16:1 (P = 0.001; 0.89 and 0.87 vs. 0.77%) were higher in LNA and Combo than in LA. Otherwise, the concentrations of C18:0 (P = 0.008; 56.62 vs. 55.45 and 55.28%) and C20:3n-6 (P = 0.001; 0.36 vs. 0.33 and 0.33%) were higher in LA than in LNA and Combo. On the other hand, the C18:2n-6 concentration was higher (P < 0.001; 3.05 and 3.66 vs. 1.82%) in LA and Combo than in LNA. The C18:3n-3 concentration was higher (P < 0.001; 1.88 vs. 1.26 vs. 0.20%) in LNA, followed by Combo and then LA. LA and LNA produced higher concentrations of C20:0 (P = 0.001; 0.80 and 0.78 vs. 0.76%) and C24:1 (P = 0.005; 0.46 and 0.46 vs. 0.42%) than did Combo. LA had a higher SFA concentration (P = 0.005; 86.87 vs. 85.11%) than that in Combo, while the concentration of SFA in LNA was not different compared with that in the other treatments. Combo had the highest PUFA concentration (P = 0.002; 5.58 vs. 3.96 and 4.37%), but it had the lowest ratio of SFAs to PUFAs (P = 0.002; 15.25 vs. 21.93 and 19.69) compared with those of the other treatments.
Table 5. Effects of fatty acid supplementation on the fatty acid profiles of rumen buffer after 8 h of incubation.
Generally, applications of LA and Combo decreased concentration of C18:2n-6 during incubation, while applications of LNA and Combo decreased concentration of C18:3n-3. At 1 h (P = 0.012; 7.61 vs. 5.61 vs. 2.68%), 2 h (P = 0.019; 3.33 and 3.48 vs. 1.96%), and 4 h (P = 0.001; 4.24 and 4.81 vs. 1.69%) were observed that the C18:2n-6 concentration was higher in LA and Combo than in LNA (Figure 4A). However, the C18:3n-3 concentration was higher in LNA and Combo than in LA at 1 h (P = 0.015; 2.51 and 1.98 vs. 0.27%), 2 h (P = 0.012; 2.13 and 2.49 vs. 0.23%), and 4 h (P = 0.008; 1.74 and 2.00 vs. 0.20%) (Figure 4B). All supplementary treatments seemed to increase C18:1cis-9 concentration during 2 h of incubation, and then decreased it after that (Figure 4C). The concentration of C18:1cis-9 was higher (P = 0.028; 7.64 and 747 vs. 7.06) in LNA and combo at 1 h. At 2 h, LNA had the highest (P = 0.008; 7.65 vs. 7.03 and 7.04%) concentration of C18:1cis-9. The C18:0 concentration was not affected by the supplementary treatment during incubation (Figure 4D). The concentration of C18:0 was increased by incubation hour in all supplementary treatments. Interaction effects between the supplementary treatment and hour were observed on C18:2n-6 (P < 0.001), C18:3n-3 (P < 0.001), and C18:0 (P = 0.036).
Figure 4. Changes in C18:2n-6 (A), C18:3n-3 (B), C18:1cis-9 (C), and C18:0 (D) concentrations (%, FA) during rumen incubation with FAs for 8 h. Diet with C18:2n-6, LA (filled circles); diet with C18:3n-3, LNA (filled triangles); and diet with the mixture of C18:2n-6 and C18:3n-3 at a ratio of 1:1 (filled squares). The respective significance levels of the supplementary treatment, incubation period, supplementary treatment × incubation period, and SEM for concentrations of C18:2n-6, C18:3n-3, C18:1n-9, and C18:0 are P < 0.001, P < 0.001, P < 0.001, and SEM = 1.915; P < 0.001, P < 0.001, P < 0.001, and SEM = 1.111; P < 0.001, P = 0.001, P = 0.218, and SEM = 0.306; and P < 0.001, P < 0.001, P = 0.036, and SEM = 1.741. Different symbols at the same hour differ significantly. ∗P < 0.05; ∗∗P < 0.01.
Discussion
Rumen Fermentation, Greenhouse Gases, and Microbes
The increases of NH3-N and total VFA concentrations and total gas volume with a decrease of rumen pH in LA indicated a higher relative rumen fermentation rate by LA compared with LNA or Combo. The reason for the high VFA concentration and total gas volume in LA could be supported by the higher total amount of fibrolytic bacteria in LA than in LNA and Combo (Table 4), which increased rumen digestion. It was found in Martin et al. (2006) that supplementary treatment with C18:3n-3 at either 100 or 50% could depress degradations of organic matter and fiber in the rumen. This previous study also reported that the population of total fibrinolytic bacteria was lower in LNA and Combo than in LA. In addition, LA had the highest acetate and propionate concentrations at 1 h of incubation, and then presented the lowest at 2 h of incubation. The reason of these results was unknown in the present study. However, on the other side, LA increased the CH4 emissions. In agreement with the present study, Zhang et al. (2008) observed higher total gas production and CH4 concentration by supplementary C18:2n-6 compared with C18:3n-3. Combo had the lowest CH4 emissions in the present study, even though it presented a high population of methanogenic archaea. In addition, supplementary n-9 FA such as oleic acid was reported to reduce methane production in an in vitro study (Wu et al., 2016). The reason for the low CH4
emissions by Combo was not clear. However, this observation supported the result of the highest CO2 concentration occurring in Combo followed by LNA and LA numerically, which could indicate a low conversion rate of CO2 to CH4 by methanogenic archaea (Wolin et al., 2007; Tapio et al., 2017). Combination of n-3 and n-6 FAs might effectively inhibit methanogenic archaea to utilize free CO2 during methanogenesis process, which could reduce methane gas emission without decrease its population. However, the further investigations need to be conducted.
Supplementary PUFAs are considered to depress methanogenesis (Fievez et al., 2003). Microbial DNA analysis showed that methanogenic archaea were suppressed only by LNA, while ciliates were suppressed by both LA and Combo. Hristov et al. (2013) reported that both C18:2n-6 and C18:3n-3 were effective to decrease rumen ciliate. Supporting the result of the present study, other previous study reported that C18:2n-6 was more effective than C18:3n-3 to decrease ciliates. However, effectiveness of C18:2n-6 against C18:3n-3 to inhibit ciliates could be influenced by supplementation levels on rumen fluid (Zhang et al., 2008). Combo also contained C18:2n-6 that could present a similar result of ciliates with LA in the present study. In the present study, the observed lower population of methanogenic archaea in LNA than in other treatments (Table 3) is in agreement with other findings (Fievez et al., 2003; Zhang et al., 2008). Methanogenic archaea and rumen ciliate might had different sensibility with each pure essential FA that caused the decrease of methanogenic archaea in LNA and rumen ciliate in LA and Combo. Along with methanogenic archaea, rumen ciliates are responsible for CH4 emissions in the rumen (Whitelaw et al., 1984). Rumen ciliates are associated with ecto- and endosymbiotic methanogenic archaea (Finlay et al., 1994). Therefore, the degradation of fiber and suppressive effects of FAs on methanogens and ciliates might have occurred and contributed to the CH4 emissions in the present study. Increased fiber degradation in the rumen is also responsible for an increase of enteric CH4 emission (Hristov et al., 2013). Decreased methanogenesis sometimes shifts rumen fermentation from acetogenic to propionigenic (Wettstein et al., 2000) because of increased hydrogen availability for propionate synthesis. However, this relation was only observed at 4 h of incubation, when CH4 production was lowest but propionate production was numerically highest in the LNA treatment (Figures 1D, 2A).
As shown in Figure 3, the population of methanogenic archaea in LNA tended to decrease over 8 h of incubation, while the other treatments showed similar populations before and after incubation. LNA has a higher degree of unsaturation than LA. Although most long-chain PUFAs could exert toxic effects on rumen microbes (Maia et al., 2007), greater effects were reported with a higher degree of unsaturation as well as a higher applied dose (Zhang et al., 2008). Both the direct toxic effects of the FA and decreased amount of ruminal hydrogen in LNA could negatively affect methanogens, compared with the effects in the other treatments. A higher degree of unsaturation in LNA consumed more hydrogen for biohydrogenation of the FA (Zhang et al., 2008). However, it was reported that only 1 to 2% of the metabolic hydrogen in the rumen is used for the purpose of biohydrogenation (Jenkins et al., 2008). Therefore, the CH4 emissions by LNA were reported to be lower than those by LA in the present study. On the other hand, LNA has previously been shown to be toxic to cellulolytic bacteria, particularly F. succinogenes, R. albus, and R. flavefaciens (Maia et al., 2007). This previous study supported the results of the present study, which reported lower populations of R. albus and F. succinogenes (numerically, P = 0.063) bacteria in LNA than in LA after 8 h of incubation. Martin et al. (2006) explained that the negative effects of PUFAs on fibrolytic bacteria were not only the result of direct toxicity of PUFAs themselves but also indirectly related to the inhibition of methanogens. Nevertheless, compared with LA, LNA presented higher R. albus populations at 1 and 4 h of incubation in the present study, but the reason for this result was not clear. The inhibition of methanogens leads to the accumulation of free hydrogen in the rumen, which can negatively affect the growth of cellulolytic bacteria (Wolin et al., 2007). The LNA treatment inhibited methanogenic archaea to a greater extent and consequently exerted maximum negative effects on fibrolytic bacteria overall compared with the effects of the other treatments. Generally, populations of all microbes decreased after 8 h of incubation because substrate was completely degraded and indicated the end of ruminal fermentation. In addition, it also supported the results of major C18 FAs, which also decrease after 8 h of incubation (Figure 4).
Fatty Acid Profiles in the Rumen
It is well known that unsaturated FAs undergo extensive biohydrogenation in the rumen by microbes (Beam et al., 2000; Petri et al., 2014). However, the rate of biohydrogenation might vary depending on the source or form of FA supplied (Carriquiry et al., 2008; Honkanen et al., 2012; Petri et al., 2014). In the present study, pure C18:2n-6 FA and C18:3n-3 FA were used either alone or in combination in equal ratios, and the results suggested less FA degradation in the combined treatment. At 0 h of incubation, the highest concentration of C18:2n-6 in LA and the highest concentration of C18:3n-3 in LNA were expected before considering the respective FAs were added accordingly in those treatments (Table 2). Furthermore, the application of mixed n-6 and n-3 FAs increased the concentrations of both C18:2n-6 FA and C18:3n-3 FA in the Combo treatment. Supplementary C18:2n-6 or C18:3n-3 did not have effect on total FA, but it could presented a minor change on proportion of FA profiles in the diet that could be a reason for results for C14:0, C14:1, C16:0, C16:1, and C18:0 at 0 h. Although the PUFAs were biohydrogenated after 8 h of incubation, the concentration patterns of remained C18:2n-6 FA and C18:3n-3 FA was similar (Table 5) to those in the beginning. This result indicated that the supplementation of these FAs could ensure their existence at an increased rate in the rumen, at least up to 8 h, which is in agreement with Carriquiry et al. (2008). Moreover, the use of these FAs in combination exerted inhibition of biohydrogenation to some extent, in contrast to the individual FA treatments. Lower C18:0 FA and total SFA concentrations and moderate C18:3n-3 concentrations and n-6 to n-3 ratios but higher C18:2n-6 FA and total PUFA concentrations in Combo indicated restricted biohydrogenation in this treatment (Table 5). This indication is more pronounced in Figure 4. This figure shows that from 0 to 2 h of incubation, approximately at 81.6% of linoleic acid was biohydrogenated in the LA treatment, and only 16.2% remained at the end of incubation (8 h). In Combo, 51.4% was biohydrogenated at 2 h, and 32.4% of C18:2n-6 remained at the end of incubation. In the case of C18:3n-3 FA, biohydrogenation was approximately at 64.2 vs. 48.3% in LNA vs. Combo, respectively at the first hour of incubation. And at the end, the remaining C18:3n-3 concentrations were approximately at 26.8 vs. 32.9%. In contrast, the concentration of C18:0 FA, the final end product of biohydrogenation of 18-carbon PUFAs, remained lower in Combo throughout the incubation. These results indicated the partial restriction of biohydrogenation of PUFAs when used in mixtures rather than individually.
Conclusion
The results indicated that better fermentation occurred in the LA treatment than in the LNA and Combo treatments, as expressed by the higher total gas volume and NH3-N and total VFA concentrations. However, the CH4 emissions were suppressed in Combo, followed by LNA and LA. Population of rumen ciliates was decreased by LA, but populations of methanogenic archaea and R. albus were decreased by LNA compared with the other treatments. Adding C18:2n-6 and C18:3n-3 in the diet can increase their proportion in the ruminal contents. A lower ratio of n-6 to n-3 can be achieved by dietary supplementation with n-3 FAs. The CH4 emissions and the extent of biohydrogenation were reduced when these FAs were used in combination rather than individually, as indicated by the lower SFA but higher PUFA concentrations in the Combo treatment than in the other treatments. Therefore, the present study recommended Combo treatment for ruminant because it was more effective to reduce enteric methane emissions and resist biohydrogenation in the rumen compared with LA and LNA.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The animal study was reviewed and approved by Animal Ethics Committee of Gyeongsang National University, South Korea (GNU-191011-E0050).
Author Contributions
SMA, DHK, and SCK designed the study. SMA, DHK, DHVP, and SCK wrote and edited manuscript. All authors performed the experiment and analyzed the data.
Funding
This research was supported (Project No. 315017-05-2-SB030) by IPET (Korea Institute of Planning and Evaluation for Technology in Food, Agriculture, Forestry and Fisheries) and Ministry of Agriculture, Food and Rural Affairs, South Korea.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Adesogan, A. T., Krueger, N. K., and Kim, S. C. (2005). A novel, wireless, automated system for measuring fermentation gas production kinetics of feeds and its application to feed characterization. Anim. Feed Sci. Technol. 12, 211–223. doi: 10.1016/j.anifeedsci.2005.04.058
Beam, T. M., Jenkins, T. C., Moate, P. J., Kohn, R. A., and Palmquist, D. L. (2000). Effects of amount and source of fat on the rates of lipolysis and biohydrogenation of fatty acids in ruminal contents. J. Dairy Sci. 83, 2564–2573. doi: 10.3168/jds.S0022-0302(00)75149-6
Beauchemin, K. A., McGinn, S. M., Benchaar, C., and Holthausen, L. (2009). Crushed sunflower, flax, or canola seeds in lactating dairy cow diets: effects on methane production, rumen fermentation, and milk production. J. Dairy Sci. 92, 2118–2127. doi: 10.3168/jds.2008-1903
Carriquiry, M., Weber, W. J., Baumgard, L. H., and Crooker, B. A. (2008). In vitro biohydrogenation of four dietary fats. Anim. Feed Sci. Technol. 141, 339–355. doi: 10.1016/j.anifeedsci.2007.06.028
Chaney, A. L., and Marbach, E. P. (1962). Modified reagents for determination of urea and ammonia. Clin. Chem. 8, 130–132.
Denman, S. E., and McSweeney, C. S. (2006). Development of a real-time PCR assay for monitoring anaerobic fungal and cellulolytic bacterial populations within the rumen. FEMS Microbiol. Ecol. 58, 572–582. doi: 10.1111/j.1574-6941.2006.00190.x
Denman, S. E., Tomkins, N. W., and McSweeney, C. S. (2007). Quantitation and diversity analysis of ruminal methanogenic populations in response to the antimethanogenic compound bromochloromethane. FEMS Microbiol. Ecol. 62, 313–322. doi: 10.1111/j.1574-6941.2007.00394.x
Doreau, M., and Ferlay, A. (1995). Effects of dietary lipids on nitrogen metabolism in the rumen: a review. Livest. Prod. Sci. 43, 97–110. doi: 10.1016/0301-6226(95)00041-I
FAO (2016). Reducing Enteric Methane for Improving Food Security and Livelihoods. New Zealand: The Food and Agriculture Organization of the United Nation.
Fievez, V., Dohme, F., Danneels, M., Raes, K., and Demeyer, D. (2003). Fish oil as potent rumen methane inhibitors and associated effects on rumen fermentation in vitro and in vivo. Anim. Feed Sci. Technol. 104, 41–58. doi: 10.1016/S0377-8401(02)00330-9
Finlay, B. J., Esteban, G., Clarke, K. J., Williams, A. G., Embley, T. M., and Hirt, R. P. (1994). Some rumen ciliates have endosymbiotic methanogens. FEMS Microbiol. Lett. 117, 157–162. doi: 10.1111/j.1574-6968.1994.tb06758.x
Honkanen, A. M., Griinari, J. M., Vanhatalo, A., Ahvenjärvi, S., Toivonen, V., and Shingfield, K. J. (2012). Characterization of the disappearance and formation of biohydrogenation intermediates during incubations of linoleic acid with rumen fluid in vitro. J. Dairy Sci. 95, 1376–1394. doi: 10.3168/jds.2011-4390
Hristov, A. N., Oh, J., Firkins, J. L., Dijkstra, J., Kebreab, E., and Waghorn, et al. (2013). Special topics—Mitigation of methane and nitrous oxide emissions from animal operations: I. a review of enteric methane mitigation options. J. Anim. Sci. 91, 5045–5069. doi: 10.2527/jas.2013-6583
Jayanegara, A., Goel, G., Makkar, H. P. S., and Becker. (2015). Divergence between prufied hydrolysable and condensed tannin effects on methane emission, rumen fermentation and microbial population in vitro. Anim. Feed Sci. Technol. 209, 60–68. doi: 10.1016/j.anifeedsci.2015.08.002
Jenkins, T. C., Thies, E. J., and Mosley, E. E. (2001). Direct methylation procedure for converting fatty amides to fatty acid methyl ester in feed and digesta samples. J. Agric. Food Chem. 49, 2142–2145. doi: 10.1021/jf001356x
Jenkins, T. C., Wallace, R. J., Moate, P. J., and Mosley, E. E. (2008). Board-invited review: recent advances in biohydrogenation of unsaturated fatty acids within the rumen microbial ecosystem. J. Anim. Sci. 86, 397–412. doi: 10.2527/jas.2007-0588
Johnson, K. A., and Johnson, D. E. (1995). Methane emissions from cattle. J. Anim. Sci. 73, 2483–2492. doi: 10.2527/1995.7382483x
Jouany, J. P., Papon, Y., Morgavi, D. P., and Doreau, M. (2008). “Linseed oil and a combination of sunflower oil and malic acid decrease rumen methane emissions in vitro,” in Livestock and Global Climate Change, eds P. Rowlinson, M. Steele, and A. Nefzaoui (Cambridge: Cambridge University Press), 140–143.
Kim, D. H., Lee, H. J., Amanullah, S. M., Adesogan, A. T., and Kim, S. C. (2016). Effects of dietary n-6/n-3 fatty acid ratio on nutrient digestibility and blood metabolites of Hanwoo heifers. Anim. Sci. J. 87, 46–53. doi: 10.1111/asj.12401
Kim, E. T., Kim, C. H., Min, K. S., and Lee, S. S. (2012). Effects of plant extracts on microbial population, methane emission and ruminal fermentation characteristics in in vitro. Asian-Aust. J. Anim. Sci. 25, 806–811. doi: 10.5713/ajas.2011.11447
Kim, S. C., Adesogan, A. T., Badinga, L., and Staples, C. R. (2007). Effects of dietary n-6:n-3 fatty acid ratio on feed intake, digestibility, and fatty acid profiles of the ruminal contents, liver, and muscle of growing lambs. J. Anim. Sci. 85, 706–716. doi: 10.2527/jas.2006-289
Koike, S., and Kobayashi, Y. (2001). Development and use of competitive PCR assays for the rumen cellulolytic bacteria: Fibrobacter succinogens, Ruminococcus albus, and Ruminococcus falvefaciens. FEMS Microbiol. Lett. 204, 361–366. doi: 10.1111/j.1574-6968.2001.tb10911.x
Machmüller, A., Ossowski, D. A., and Kreuzer, M. (2000). Comparative evaluation of the effects of coconut oil, oilseeds and crystalline fat on methane release, digestion and energy balance in lambs. Anim. Feed Sci. Technol. 85, 41–60. doi: 10.1016/S0377-8401(00)00126-7
Machmüller, A., Ossowski, D. A., Wanner, M., and Kreuzer, M. (1998). Potential of various fatty feeds to reduce methane release from rumen fermentation in vitro (RUSITEC). Anim. Feed Sci. Technol. 71, 117–130. doi: 10.1016/S0377-8401(97)00126-0
Maia, M. R. G., Chaudhary, L. C., Figueres, L., and Wallace, R. J. (2007). Metabolism of polyunsaturated fatty acids and their toxicity to the microflora of the rumen. Antonie Van Leeuwenhoek. 91, 303–314. doi: 10.1007/s10482-006-9118-2
Martin, C., Rouel, J., Jouany, J. P., Doreau, M., and Chilliard, Y. (2006). Methane output and diet digestibility in response to feeding dairy cows crude linseed, extruded linseed, or linseed oil. J. Anim. Sci. 86, 2642–2650. doi: 10.2527/jas.2007-0774
Mauricio, R. M., Mould, F. L., Dhanoa, M. S., Owen, E., Channa, K. S., and Theodorou, M. K. (1999). A semi-automated in vitro gas production technique for ruminant feedstuff evaluation. Anim. Feed Sci. Technol. 79, 321–330. doi: 10.1016/S0377-8401(99)00033-4
Muck, R. E., and Dickerson, J. T. (1988). Storage temperature effects on proteolysis in alfalfa silage. Trans ASAE. 31, 1005–1009. doi: 10.13031/2013.30813
Paradhipta, D. H. V., Joo, Y. H., Lee, H. J., Lee, S. S., Kwak, Y. G., Han, O. K., et al. (2020). Effects if wild or mutated inoculants on rye silage and its rumen fermentation indices. Asian-Australas. J. Anim. Sci. 33, 949–956. doi: 10.5713/ajas.19.0308
Petri, R. M., Mapiye, C., Dugan, M. E. R., and McAllister, T. A. (2014). Subcutaneous adipose fatty acid profiles and related rumen bacterial populations of steers fed red clover or grass hay diets containing flax or sunflower-seed. PLoS One 9:e104167. doi: 10.1371/journal.pone.0104167
Skillman, L. C., Toovey, A. F., Williams, A. J., and Wright, A. G. (2006). Development and validation of a real-time PCR method to quantify rumen protozoa and examination of variability between Endodinium populations in sheep offered a hay-based diet. Appl. Environ. Microb. 72, 200–206. doi: 10.1128/AEM.72.1.200-206.2006
Tapio, I., Snelling, T. J., Strozzi, F., and Wallace, R. J. (2017). The ruminal microbiome associated with methane emissions from ruminant livestock. J. Anim. Sci. Biotechnol. 8:7. doi: 10.1186/s40104-017-0141-0
Wettstein, H. R., Machmüller, A., and Kreuzer, M. (2000). Effect of raw and modified canola lecithins compared to canola oil, canola seed and soy lecithin on ruminal fermentation measured with rumen simulation technique. Anim. Feed Sci. Technol. 85, 153–169. doi: 10.1016/S0377-8401(00)00149-8
Whitelaw, F. G., Eadie, J. M., Bruce, L. A., and Shand, W. J. (1984). Methane formation in faunated and ciliate-free cattle and its relationship with rumen volatile fatty acid proportions. Br. J. Nutr. 52, 261–275. doi: 10.1079/BJN19840094
Wolin, M. J., Miller, T. L., and Stewart, C. S. (2007). “Microbe-microbe interactions,” in The Rumen Microbial Ecosystem, eds P. N. Hobson and C. S. Stewart (London: Blackie Academic and Professional Press), 467–491.
Wu, D., Xu, L., Tang, S., Guan, L., He, Z., and Guan, Y. (2016). Influence of oleic acid on rumen fermentation and fatty acid formation in vitro. Plos One 11:e0156835. doi: 10.1371/journal.pone.0156835
Yang, S. L., Bu, D. P., Wang, J. Q., Hu, Z. Y., Li, D., Wei, H. Y., et al. (2009). Soybean oil and linseed oil supplementation affect profiles of ruminal microorganisms in dairy cows. Animal 3, 156–1569. doi: 10.1017/S1751731109990462
Keywords: biohydrogenation, fatty acid, methane emission, rumen fermentation, rumen microbes
Citation: Amanullah SM, Kim DH, Paradhipta DHV, Lee HJ, Joo YH, Lee SS, Kim ET and Kim SC (2021) Effects of Essential Fatty Acid Supplementation on in vitro Fermentation Indices, Greenhouse Gas, Microbes, and Fatty Acid Profiles in the Rumen. Front. Microbiol. 12:637220. doi: 10.3389/fmicb.2021.637220
Received: 03 December 2020; Accepted: 10 February 2021;
Published: 11 March 2021.
Edited by:
Kwangcheol Casey Jeong, University of Florida, United StatesReviewed by:
Raies Mir, United States Department of Agriculture (USDA), United StatesSang-Suk Lee, Sunchon National University, South Korea
Copyright © 2021 Amanullah, Kim, Paradhipta, Lee, Joo, Lee, Kim and Kim. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sam Churl Kim, a2ltc2NAZ251LmFjLmty
†These authors have contributed equally to this work