Features of R. oryzae and macroscopic appearance of adlay millet grits before and after fermentation
The colony of R. oryzae was loose on PDA medium and formed sporangia in later stage (Fig. 1a). The sporangia was black brown and spherical with erect sporangiophores (Fig. 1b). There were obvious stolons and rhizoids without diaphragm in the mycelium (Fig. 1c). The sporangiospores were oval and the funicular of sporangia was wedge-shaped (Fig. 1d). The adlay millet grits before fermentation was compact and well defined with each other (Fig. 1e,h). After 24 h of fermentation, there was clear liquid seeped out with pleasant sweet scent and some bubbles, and the adlay millet grits seemed relatively intact (Fig. 1f,i). After 54 h of fermentation, there was more glutinous liquid seeped out with clear wine aromas and more bubbles, and a membrane was formed by R. oryzae mycelium on the surface of the loosed adlay millet grits (Fig. 1g,j).
Metabolomic analysis of CAG and FAG samples
To investigate the global variations of metabolites in adlay millet seeds before and after R. oryzae fermentation, the metabolomic analysis was conducted for their accurate identification and annotation. Differences in metabolomics data from the CAG, FAG, and QC samples were visualized by PCA (Supplementary Fig. 1a). PC1 and PC2 explained 60.5% and 10.6% (total 71.1%) of the variation between CAG and FAG samples, respectively; all samples were within the 95% confidence interval, i.e., there were no outliers. The CAG samples were more closely grouped than the FAG samples, but the two groups were widely separated from each other and the QC samples, indicating a significant difference in metabolite profiles resulting from the R. oryzae fermentation.
In the OPLS-DA model, the FAG samples were more closely grouped, and the groups were still well-separated (Supplementary Fig. 1b). Permutation testing of the OPLS-DA model was performed (Supplementary Fig. 2); the abscissa represented the permutation retention, the point where the permutation retention equals 1 corresponds to the R2Y (goodness of fit) and Q2 (goodness of prediction) values of the original model, which were both close to 1, validating the authenticity, robustness and stability of the OPLS-DA model. Specifically, 718 metabolites were annotated in CAG and FAG samples from 18 compound classes (Fig. 2), of which the top 10 most abundant metabolites were flavonoids (113 metabolites), terpenoids (101 metabolites), alkaloids (78 metabolites), phenols (63 metabolites), amino acids and derivatives (53 metabolites), lipids and lipid-like molecules (44 metabolites), coumarins and lignans (37 metabolites), steroids and derivatives (28 metabolites), organic acids and derivatives (23 metabolites), and carbohydrates and conjugates (23 metabolites).
Identification and classification of differential metabolites
To analyze metabolic profile differences between FAG and CAG samples at the individual compound level, a hierarchical cluster analysis of the DMs was performed based on a Euclidean distance matrix and the complete linkage method (Supplementary Fig. 3). The abundance of most DMs was much higher in FAG than in CAG samples (Supplementary Fig. 3), indicating that the microbial fermentation had a great influence on the metabolic profile of adlay millet seed. A total of 203 differential metabolites were identified, of which 184 were more abundant and 19 were less abundant in FAG (Fig. 3a). The top ten DM compound classes were amino acids and derivatives (17.24%), flavonoids (11.82%), terpenoids (11.82%), alkaloids (8.37%), carbohydrates and conjugates (6.40%), lipids and lipid-like molecules (6.40%), organic acids and derivatives (5.42%), phenols (6.40%), nucleotides and derivatives (2.96%), and steroids (2.46%) (Fig. 3b). The abundances of amino acids and derivatives, nucleosides, nucleotides and derivatives, and terpenoids all increased (Fig. 3c), of which the first two classes included the three metabolites with the largest increases in abundance, namely, S-adenosylmethionine (SAMe), uridine diphosphate glucose (UDP-glucose), and β-Nicotinamide mononucleotide (β-NMN) (Fig. 3d). SAMe has health-beneficial biological activities, including prevention and treatment of arthritis, liver disease, depression, and jaundice22. β-NMN can ameliorate health conditions associated with nicotinamide
adenine dinucleotide (NAD+) deficiency, such as aging, cardiac, and cerebral ischemia, Alzheimer’s disease and obesity23. In addition, most compounds in the classes of flavonoids, phenols, organic acids, lipids, and lipid-like molecules were more abundant.
Metabolism of nutritional components
Amino acids and derivatives
Adlay millet seed has been identified as a good source of dietary protein because of its relatively high protein/carbohydrate ratio compared with other cereals1. In this study, the abundance of amino acids and derivatives notably increased (Fig. 3c, Supplementary Table 1), indicating that the R. oryzae fermentation induced positive effect on amino acids contents. In addition, the abundance of SAMe increased the most markedly (11,285-fold) of all the DMs (Table 1, Fig. 3d). SAMe is a key pleiotropic molecule involved in three types of cellular reactions: transmethylation, transsulfuration, and aminopropylation. It has been demonstrated to increase the antiviral effect of interferon, and
SAMe supplementation could restore hepatic glutathione (GSH) deposits and attenuate liver injury, treat arthritis, and depression, and reduce jaundice22. Gama-Aminobutyric acid (GABA) increased by 16.72-fold in FAG and was considered to be a potent bioactive compound because of its many physiological functions, including lowering blood pressure, enhancing memory, anti-inflammation, obesity prevention, enhancing immunity, and reducing insomnia24. Similarly, fermentation of coix/adlay bran by B. subtilis resulted in an almost 14-fold increase in GABA8 and yeast fermentation increased the GABA content 6-fold9. In general, the R. oryzae fermentation improved the trophism and abundance of important amino acids.
Table 1
Representative DMs detected in samples of FAG compared with CAG.
Compounds | Fold Change | Type | Compounds | Fold Change | Type |
Amino acids and derivatives | Nucleosides, nucleotides and derivatives |
S-adenosyl methionine (SAMe) | 11285.42 | up | Uridine 5'-diphospho-D- glucose (UDP-glucose) | 5789.74 | up |
L-alanine | 24.80 | up | β-Nicotinamide mononucleotide (β-NMN) | 1189.04 | up |
γ-aminobutyric acid (GABA) | 16.72 | up | Vitamins |
L-Glutamate | 10.17 | up | Nicotinic acid | 10.39 | up |
L-Aspartate | 7.98 | up | Flavonoids | | |
Carbohydrates | Gossypetin | 8.46 | up |
β-Cyclodextrin | 210.97 | up | Liquiritin | 8.05 | up |
Kojibiose | 100.31 | up | (-)-Epiafzelechin | 6.28 | up |
β-D-Fructose-2-phosphate | 52.14 | up | Terpenoids | | |
α-D-Glucose | 12.40 | up | Valechlorine | 79.65 | up |
Sucrose | 0.37 | down | Germanicol | 13.38 | up |
Lactulose | 0.28 | down | Simiarenol | 10.46 | up |
Stachyose | 0.25 | down | Trans-caryophyllene | 6.96 | up |
Lipids and lipid-like molecules | Alkaloids | | |
Docosahexaenoic acid (DHA) | 37.60 | up | Rhynchophylline | 21.10 | up |
1, 5-pentane dicarboxylic acid (pimelic acid) | 17.80 | up | Lupanine | 6.81 | up |
Organic acids and derivatives | Phenols | | |
(S)-2-Acetolactate | 37.18 | up | Maltol | 16.06 | up |
cis-Aconitic acid | 21.59 | up | 5-Heptadecylresorci–nol | 12.21 | up |
2-Isopropyl-3-oxosuccinate | 18.43 | up | 10-Gingerol | 5.80 | up |
Carbohydrates and carbohydrate conjugates
Carbohydrate metabolism is one of the most important cellular metabolic processes and provides material and energy for growth and proliferation25. A total of 12 carbohydrate and carbohydrate conjugate DMs were identified (Fig. 3c, Supplementary Table 1). Stachyose decreased significantly (Fig. 3d), followed by lactulose and sucrose (Table 1), whereas the Kojibiose, a good prebiotic, was increased by 100-fold. Kojibiose was found in sake and koji extracts, and was reported to strongly inhibit α-glucosidase I activity and showed antitoxic competence26.
Lipids and lipid-like molecules
Adlay millet seed is richer in lipids than most cereal grains3. 44 lipids and lipid-like molecules were annotated in FAG and CAG, including various unsaturated fatty acids, such as Docosahexaenoic acid (DHA), 1, 5-pentane dicarboxylic acid (Pimelic acid), and conjugated linoleic acid (Bovinic acid) (Supplementary Table 1). Among them, 14 kinds of lipids were found as DMs (Fig. 3b, c). The abundance of DHA was 38-fold (the biggest increase of the lipids) more abundant in FAG than CAG (Table 1). DHA has been linked to many positive health effects, e.g. promoting the growth of nerve cells, and alleviating cardiovascular and inflammatory diseases27. The abundance of pimelic acid increased 18-fold after R. oryzae fermentation. Pimelic acid is the initial precursor for biotin synthesis, suggesting that co-culture of R. oryzae with other biotin-producing microorganisms, such as B. subtilis, would increase the biotin yield from fermentation28.
Nucleosides, nucleotides and derivatives
Nucleosides and nucleotides are not only components of genetic material, but also key components of cellular metabolic pathways. 12 compounds in this class were annotated in FAG and CAG, 6 of which were DMs, UDP-glucose and β-NMN were increased most (Table 1). UDP-glucose increased more than 5000-fold in FAG compared with CAG (Fig. 3d, Table 1). As one of the most important glycosyl donors for biosynthesis of glycosides, oligosaccharides, polysaccharides, and glycoproteins, UDP-glucose was closely associated with the biosynthesis of nucleotide sugars, galactose metabolism, amino sugar, and nucleotide sugar metabolism. It was notably the β-NMN was upregulated about 1200-fold in FAG samples (Fig. 3d, Table 1). It is an endogenous bioactive nucleotide and a key precursor for NAD+ biosynthesis29. Beta-MNM has beneficial pharmacological effects on health problems arising from NAD+ deficiency, such as aging, cardiac, and cerebral ischemia, Alzheimer’s disease and obesity23, so NMN-rich fermented adlay millet seed merits further study to evaluate its potential health benefits.
Vitamins
Vitamins are important nutritional components of adlay millet seed; vitamin E and β-carotene have been previously identified in adlay millet seed2. In FAG and CAG, 10 vitamins were annotated, including vitamin E, β-carotene, ascorbic acid, nicotinic acid, vitamin A, pyridoxine, and thiamine, the last four of which were DMs (Table 1). The abundance of nicotinic acid increased 10.4-fold in FAG; it is a component of coenzyme I (NAD+) and coenzyme II (NADP+). At pharmacological level, its supplementation was beneficial for human health to combat cardiovascular disease, neurological problems, diabetes, and other skin disease30.
Metabolism of functional components
Flavonoids
Flavonoids form a large group of bioactive plant secondary metabolites that is abundant in fruits and seeds, with human health benefits, such as anti-inflammatory, anti-aging, cardio-protective and immunomodulatory31. A total of 24 flavonoid compounds was annotated as DMs in FAG compared with CAG, namely 16 flavones (abyssinone V, gossypetin, liquiritin, (-)-epiafzelechin, kaempferol 3-O-rhamnoside, cyanidin, homoferreirin, naringenin chalcone, methyl hesperidin, quercitrin, 3-O-acetylpinobanksin, sciadopitysin, baicalein, morin, hinokiflavone, and 3,8'-biapigenin), 6 flavanones (spinosin, neohesperidin, diosmin, hesperidin, luteolin, and tangeretin), and 2 isoflavones (cajanol and 5,7,3',4'-tetrahydroxyisoflavone). Of these, 21 were more abundant in FAG and three were less abundant (Fig. 3c, Supplementary Table 1). Gossypetin (8-fold more abundant in FAG) is a hydroxylated flavone with strong antioxidant, antimicrobial, anti-inflammatory and anti-atherosclerotic activities32. Liquiritin and Spinosin were benifical for treatment of Alzheimer’s disease33. Other flavonoids, such as Cyanidin, Hesperidin, and Quercetrin have anti-inflammatory and anti-cancer bioactivity31.
Terpenoids
Terpenoids constitute one of the largest classes of natural compounds present in many plants. They have many biological activities, such as anti-bacterial, anti-parasitic, anti-inflammatory, anti-tumor, and skin protection34. In FAG, 24 terpenoid compounds were assigned as DMs, i.e., 3 monoterpenoids, 2 iridoids, 6 sesquiterpenes, 1 diterpene, 11 triterpenes, and 1 tetraterpenoid (Fig. 3b, Supplementary Table 1). All the terpenoid DMs were more abundant in FAG (Fig. 3c); the biggest increase was for Valechlorine (79-fold increase), the eighth biggest increase of all the DMs (Fig. 3d). Valechlorine alleviates hepatic steatosis by enhancing autophagy, and reduces oleic acid-induced lipid accumulation32. The triterpene germanicol (13-fold increase) has anti-fungal activity35. R. oryzae fermentation effectively increased the contents of terpenoids in adlay millet seed, which should improve its biological functions and potential health value.
Alkaloids
Alkaloids have been isolated from many medicinal herbs and represent a rich source of bioactive molecules36. A total of 17 alkaloid DMs were identified in FAG samples, of which 13 were more abundant (Supplementary Table 1). The abundance of rhynchophylline, (21-fold increase) increased the most after fermentation; this compound is a well-established treatment for central nervous system diseases37. Lupanine is a positive modulator of insulin release, which helps to maintain glucose homeostasis38.
Phenols
Phenols are one of the most widely distributed classes of plant secondary metabolites and have been widely applied in the food and pharmaceutical industries as flavorings, antioxidants and antibacterial agents39. A total of 63 phenol compounds were annotated in FAG and CAG; of these, 13 were DMs and 12 were more abundant in FAG (Fig. 3b, c). The abundance of Maltol increased 16-fold, the largest increase of the phenolic compounds (Table 1). Maltol has various biological activities, such as neuroprotection, liver protection against acute carbon tetrachloride toxicity, reduction of drug-induced hepatotoxicity, and inhibition of osteoarthritis progression40. The abundance of 5-heptadecylresorcinol increased 12-fold (Table 1); it is a member of the 5-n-alkylresorcinol (AR) family, found in many cereals. 5-Heptadecylresorcinol has a stronger inhibitory effect on cancer-cell proliferation than most other ARs41. The compounds of 10-Gingerol and curcumin, which have antitumor, antioxidant and anti-inflammatory effects42, were both more abundant in FAG. The results were similar to previous studies observed for increased phenolics from defatted rice bran by R. oryzae fermentation12. Phenolic compounds in the cereals are predominantly bound to cell wall components, but the R. oryzae fermentation liberates them into the free form, increasing the apparent content and bioaccessibilities12.
KEGG annotation, related pathway enrichment and network analysis of differential metabolites
The metabolism of various substances in the biological body is not carried out in isolation, but through the common intermediate substances to connect and transform to each other, forming a complete metabolic network to regulate the metabolism and perform cellular biological functions. To aid understanding of the metabolic activities of R. oryzae fermentation on adlay millet seed, the 203 DMs were mapped to KEGG pathways. A total of 120 DMs was mapped to the KEGG data sets and they were found to be enriched in 58 metabolic pathways (Fig. 4); most DMs were related to metabolic pathways and biosynthesis of secondary metabolites. The pathways of amino acids biosynthesis and D-amino acid metabolism showed significant enrichment of DMs (Fig. 4). In detail, the alanine, aspartate, and glutamate metabolism pathway enriched abundant DMs and had high enrichment factor (Fig. 4), indicating it played important role in the fermentation of adlay millet seed by R. oryzae.
By analyzing the several specific amino acid metabolic pathways in which the DMs were mainly enriched, the metabolic network among them were shown in Fig. 5. All DMs in the network were more abundant in the FAG samples compared with CAG samples (Supplementary Table 1). The pathways of alanine, aspartate and glutamate metabolism (Fig. 5, purple box), cysteine and methionine metabolism (Fig. 5, blue box), and glycine, serine and threonine metabolism (Fig. 5, green box) were interconnected with each other by a common DM, i.e., aspartate. It was obvious that the aspartate was transformed from the intermediate metabolite oxaloacetate in the tricarboxylic acid cycle (TCA) cycle by the pathway of alanine, aspartate and glutamate metabolism (Fig. 5, purple box). The aspartate could be transformed to homoserine and finally to generate SAMe through the cysteine and methionine metabolism pathway (Fig. 5, blue box). Besides, abundant aspartate could also be participated in the nicotinate and nicotinamide metabolism pathway (Fig. 5, Supplementary Fig. 4), leading to increased biosynthesis of
nicotinic acid (Table 1). The nicotinic acid could increase NAD+ content through Preiss-Handler pathway and should decrease the β-NMN consumption, so resulting in β-NMN accumulation41 (Supplementary Fig. 4, Table 1). Aspartate was also enriched in the glycine, serine and threonine metabolism pathway (Fig. 5, green box), the more abundant glycine, serine and threonine should be closely related to generation of terpenoids and steroids through methylglutarate pathway44. Other intermediate metabolite, i.e. 2-oxoglutarate in the TCA cycle, could be converted to glutanate, and
then to GABA and glutamine ghrough alanine, aspartate and glutamate metabolism pathway (Fig. 5, purple box). Glutamate was also enriched in the pathways of arginine and proline metabolism and arginine biosynthesis. The increased content of arginine and proline can be used as precursors to synthesize alkaloids45. Overall, amino acid metabolism played a crucial role in the transformation of primary metabolites into secondary metabolites by R. oryzae fermentation. Specifically, L-aspartate appeared to play a vital role in the metabolic network for biosynthesis of SAMe, β-NMN and other functional compounds.
Quantitative verification of typical DMs and important amino acids
Three metabolites with the most significant increase in abundance, i.e. SAMe, UDP-glucose, and β-NMN, and the important amino acids related to their metabolism were further quantitative verified. As shown in Fig. 6, SAMe, UDP-glucose, and β-NMN all exhibited extremely significant increase in FAG compared with CAG. The content of SAMe showed a sharp increase from 0.56 to 370.26 µg/g after R. oryzae fermentation; and the content of UDP-glucose and β-NMN increased significantly from 0.64 µg/g to 16.14 µg/g and 0.55 µg/g to 8.32 µg/g, respectively. The aspartate which was closely related to both SAMe and β-NMN biosynthesis showed a significant increase in content from 0.35–0.52%. Other correlative amino acids were also obviously increased after R. oryzae fermentation, including glutamate (from 1.31–1.70%), proline (from 0.48–0.62%) and arginine (from 0.24–0.31%), glycine (from 0.13–0.19%), serine (from 0.23–0.32%), and threonine (0.16–0.25%). Previous studies have reported R. oryzae fermentation effectively increased the content of free amino acids in defatted adlay bran and other cereals13,46. The quantitative results were consistent with that of the metabolomics, confirming that the amino acid metabolism were important for R. oryzae fermentation of adlay millet seed to produce functional compounds.