Genetic, transcriptional, and regulatory landscape of monolignol biosynthesis pathway in Miscanthus × giganteus
Biotechnology for Biofuels volume 13, Article number: 179 (2020)
Miscanthus × giganteus is widely recognized as a promising lignocellulosic biomass crop due to its advantages of high biomass production, low environmental impacts, and the potential to be cultivated on marginal land. However, the high costs of bioethanol production still limit the current commercialization of lignocellulosic bioethanol. The lignin in the cell wall and its by-products released in the pretreatment step is the main component inhibiting the enzymatic reactions in the saccharification and fermentation processes. Hence, genetic modification of the genes involved in lignin biosynthesis could be a feasible strategy to overcome this barrier by manipulating the lignin content and composition of M. × giganteus. For this purpose, the essential knowledge of these genes and understanding the underlying regulatory mechanisms in M. × giganteus is required.
In this study, MgPAL1, MgPAL5, Mg4CL1, Mg4CL3, MgHCT1, MgHCT2, MgC3′H1, MgCCoAOMT1, MgCCoAOMT3, MgCCR1, MgCCR2, MgF5H, MgCOMT, and MgCAD were identified as the major monolignol biosynthetic genes in M. × giganteus based on genetic and transcriptional evidence. Among them, 12 genes were cloned and sequenced. By combining transcription factor binding site prediction and expression correlation analysis, MYB46, MYB61, MYB63, WRKY24, WRKY35, WRKY12, ERF021, ERF058, and ERF017 were inferred to regulate the expression of these genes directly. On the basis of these results, an integrated model was summarized to depict the monolignol biosynthesis pathway and the underlying regulatory mechanism in M. × giganteus.
This study provides a list of potential gene targets for genetic improvement of lignocellulosic biomass quality of M. × giganteus, and reveals the genetic, transcriptional, and regulatory landscape of the monolignol biosynthesis pathway in M. × giganteus.
Miscanthus × giganteus is a triploid perennial rhizomatous C4 grass that originated from the natural hybridization between diploid Miscanthus sinensis and tetraploid Miscanthus sacchariflorus . Owing to its outstanding features, such as high biomass production , low environmental impacts , and the potential to be cultivated on marginal land , M. × giganteus is widely recognized as a promising lignocellulosic biomass crop for bioethanol production. However, the recalcitrant nature of lignocellulosic feedstocks leads to the high costs of pretreatment, saccharification, and fermentation processes, limiting the current commercialization of lignocellulosic bioethanol [5, 6].
Among the biopolymers in lignocellulosic biomass, the lignin enriched in the secondary cell wall is one of the main factors that account for recalcitrance . Besides, its by-products released in the pretreatment step are the primary inhibitors of enzymatic reactions in the saccharification and fermentation processes [8, 9]. To overcome this barrier, researchers have focused on manipulating the lignin content and composition of the lignocellulosic feedstocks via genetic engineering approaches. In maize [10, 11] and switchgrass [12,13,14], the reduction of lignin content and optimization of lignin composition significantly promoted the saccharification efficiency and ethanol productivity. A similar correlation was observed in the natural Miscanthus accessions [15, 16]. Additionally, the sterile nature of triploid M. × giganteus enhances the environmental safety of genetic engineering. These results suggest the potential utilization of genetic manipulation of lignin biosynthesis in M. × giganteus. To this end, basic knowledge of the genes involved in lignin biosynthesis and how these genes are regulated in M. × giganteus is essential.
The lignin in the cell wall is mainly composed of p-hydroxyphenyl (H), guaiacyl (G), and syringyl (S) units, which are polymerized from the corresponding monolignols, p-coumaryl, coniferyl, and sinapyl alcohols, respectively. In flowering plants, these monolignols are synthesized from the general phenylpropanoid pathway and the following monolignol-specific pathway, as shown in Fig. 1a . Phenylalanine ammonia-lyase (PAL) and cinnamic acid 4-hydroxylase (C4H) are the first two enzymes in the general phenylpropanoid pathway catalyzing the synthesis of p-coumaric acid from phenylalanine. Recently, a bi-functional cytosolic ascorbate peroxidase (APX) was reported to function as 4-coumarate 3-hydroxylase (C3H) synthesizing caffeic acid through the 3-hydroxylation of p-coumaric acid . The subsequent 3-O-methylation from caffeic acid to ferulic acid is catalyzed by caffeic acid/5-hydroxyconiferaldehyde 3/5-O-methyltransferase (COMT). Then, these hydroxycinnamates are converted to the corresponding CoA esters by 4-hydroxycinnamate: CoA ligase (4CL). A more generally accepted route for caffeoyl CoA synthesis is derived from p-coumaroyl CoA via the catalyzing of hydroxycinnamoyl CoA: shikimate hydroxycinnamoyl transferase (HCT) and coumaroyl shikimate 3′-hydroxylase (C3′H). The resulting caffeoyl CoA is 3-O-methylated by caffeoyl CoA 3-O-methyltransferase (CCoAOMT). In the monolignol-specific pathway, these CoA esters are converted to corresponding aldehydes and alcohols by cinnamoyl CoA reductase (CCR) and cinnamyl alcohol dehydrogenase (CAD), respectively. Besides, some enzymes in the monolignol biosynthesis pathway are involved in the flux from G to S monolignol at the levels of aldehydes and alcohols, including COMT and ferulic acid/coniferaldehyde 5-hydroxylase (F5H).
However, the monolignol biosynthesis pathway is variable across different species. In some monocots, PAL also exhibits the capability to catalyze the non-oxidative elimination of ammonia from tyrosine to p-coumaric acid directly, known as the tyrosine shortcut pathway [19,20,21,22]. Another example is caffeoyl shikimate esterase (CSE). It converts caffeoyl CoA to caffeic acid in some species, while its orthologs are absent in the genomes of Brachypodium distachyon, Zea mays, and Sorghum bicolor . This diversity highlights the need to understand the monolignol biosynthesis pathway in every lignocellulosic crop.
Indeed, some studies have focused on the genetic and regulatory basis of the monolignol biosynthesis pathway in some bioenergy grasses [24,25,26,27,28]. However, in Miscanthus species, comprehensive research is still lacking, and no integrated model has been built yet. To fill this gap, we identified the major monolignol biosynthetic genes in M. × giganteus and predicted the probable transcription factors (TFs) that directly regulate these genes. Based on these results, an integrated model was summarized to depict the monolignol biosynthesis pathway and the underlying regulatory mechanism (Fig. 1a, b). This study provides a list of potential gene targets for genetic improvement of lignocellulosic biomass quality of M. × giganteus. Also, it reveals the genetic, transcriptional, and regulatory landscape of the monolignol biosynthesis pathway in M. × giganteus.
Evolutionary history of monolignol biosynthetic genes in angiosperms
The cDNAs of 20 monolignol biosynthetic genes were cloned and sequenced from M. × giganteus. The basic information and GenBank accession number of each cDNA are listed in Table 1.
The orthologous relationships of genes could provide evidence for inferring the gene evolutionary history and functions . For this reason, we performed the genome synteny and phylogenetic analyses on the monolignol biosynthetic genes from the basal angiosperm Amborella trichopoda, dicot Arabidopsis thaliana and monocots Z. mays, S. bicolor, M. sinensis and M. × giganteus (Fig. 2a, b and Additional file 1: Fig. S2). It can be observed that the isozyme genes in angiosperms shared the same origin. Compared to Amborella trichopoda, PAL, 4CL, CCoAOMT, and CCR genes in both dicots and monocots were remarkably expanded, which could be explained by a series of whole-genome duplication (WGD) and small-scale duplication events (e.g., tandem duplication). Besides, the genes in monocots were divided into several clades that are independent of the dicot clades. This result is opposite to our previous research on cellulose synthase genes, which formed six clades posterior to the divergence between dicots and monocots .
Here we take PAL as an example. There is only one PAL gene in Amborella trichopoda (AMTR_s00148p00088930), but in Arabidopsis thaliana, Z. mays, S. bicolor, and M. sinensis, the gene numbers increase to 4, 11, 8, and 13, respectively (Fig. 2b). Monocot PAL genes were grouped into five clades, parallel to the clade of Arabidopsis thaliana. In the chromosome 7 and 8 of M. sinensis genome, Misin07G42300 and Misin08G206300, Misin07G412400 and Misin08G206400 are two pairs of PAL genes derived from a recent genus-specific WGD event (red and blue ribbons in Fig. 2c, respectively). After that, Misin07G412400 underwent additional tandem duplication events, forming Misin07G412500 and Misin07G412600 (blue arrows in Fig. 2c).
These results suggest that the monolignol biosynthetic genes were expanded and independently evolved in monocots and dicots, implying the more complex nature of organization and regulation of the pathway than in the basal angiosperm Amborella trichopoda.
Expression analysis reveals major monolignol biosynthetic genes
The relative expression levels of the monolignol biosynthetic genes from all monocot clades were determined in leaves, sheaths, roots, rhizome buds, nodes, and internodes (Fig. 3a). The isozyme genes exhibited similar or different expression patterns to each other. For instance, the relative expression pattern of MgPAL1 was quite identical to MgPAL5, whereas distinct from MgPAL4. It indicates that these duplicated genes have partially specialized at the expression level.
The internode is the primary site of lignin biosynthesis. In Arabidopsis thaliana  and Z. mays , the major monolignol biosynthetic genes showed remarkably higher expression in internodes than in leaves (also see Additional file 2: Table S5). Additionally, most of them were highly expressed in roots, which are also rich in lignified vascular tissues. Based on this conserved expression pattern, the major monolignol biosynthetic genes in M. × giganteus could be estimated. For a more intuitive comparison, the relative expression levels of all the monolignol biosynthetic genes in leaves, roots, and internodes were clustered by gene and visualized in a heatmap (Fig. 3b). The result shows that the 14 genes, MgPAL1, MgPAL5, Mg4CL1, Mg4CL3, MgHCT1, MgHCT2, MgC3′H1, MgCCoAOMT1, MgCCoAOMT3, MgCCR1, MgCCR2, MgF5H, MgCOMT, and MgCAD, were clustered together, sharing consistent expression patterns with the major monolignol biosynthetic genes in Arabidopsis thaliana and Z. mays.
However, the expression levels between two genes cannot be directly compared using relative qPCR due to different amplicon lengths, amplification efficiencies, and fluorescence thresholds. Therefore, we performed transcriptome analysis in M. × giganteus to determine the absolute expression levels of isozyme genes and narrow down the number of major monolignol biosynthetic gene candidates. The overall read mapping rate and E90N50 of the transcriptome assembly were 93.58% and 1,778 bp, respectively, suggesting the high completeness and continuity of the transcripts. As we expected, the two reference genes used in this study, eEF-1a and UBQ, were steadily expressed in all vegetative organs (Additional file 1: Fig. S3). Furthermore, Mg4CL1, MgHCT1, MgHCT2, MgCCoAOMT1, MgCCoAOMT3, MgCCR1, and MgCCR2 were highly or moderately expressed in the vegetative and reproductive organs above the ground. In contrast, Mg4CL3 was rarely expressed in any samples, indicating that Mg4CL3 is not likely to act as the major monolignol biosynthetic gene (Fig. 3c).
In conclusion, MgPAL1, MgPAL5, Mg4CL1, Mg4CL3, MgHCT1, MgHCT2, MgC3′H1, MgCCoAOMT1, MgCCoAOMT3, MgCCR1, MgCCR2, MgF5H, MgCOMT, and MgCAD are most likely to be the major monolignol biosynthetic genes. It is worth mentioning that neither MgC4H1 nor MgC4H2 was relatively highly expressed in internodes, indicating tyrosine shortcut pathway may be the primary route for 4-coumaric acid biosynthesis. Therefore, MgPAL1 or MgPAL5 should be able to utilize tyrosine as the substrate in theory. Consistently, the two orthologs of MgPAL1 and MgPAL5 in S. bicolor, SORBI_3004G220300 and SORBI_3006G148800 (Fig. 2b), were examined to have such substrate affinity . The key residue, histidine in the 4-methylidene-imidazole-5-one (MIO) domain conferring the function, could also be found in MgPAL1 and MgPAL5 (Additional file 1: Fig. S4).
Asymmetric evolution between major and non-major monolignol biosynthetic genes
According to the neofunctionalization model and plenty of studies on duplicated genes, the genes that preserved the original functions evolve slower than other copies, which is referred to as “asymmetric evolution” . We wondered whether asymmetric evolution could also be observed in the monolignol biosynthetic isozyme genes of M. × giganteus. As expected, the coding sequences of major monolignol biosynthetic genes exhibited significantly higher percent identities than the non-major genes compared to the corresponding orthologs in S. bicolor (Wilcoxon-rank sum test, p value = 0.01262) (Fig. 4a). Furthermore, the Ka/Ks ratios showed a contrary tendency (Fig. 4b). These results demonstrate that the major monolignol biosynthetic genes have higher sequence conservation and underwent stronger purifying selection. In comparison, the rest genes evolved faster at both the transcription level and the sequence level. This finding is also consistent with the critical role of the monolignol biosynthesis in plant survival and development.
MgCCR1 and MgCCR2 are a pair of genes formed in the recent genus-specific WGD event. Interestingly, compared to MgCCR2, MgCCR1 has rapidly accumulated mutations to the extent of MgCCR3 and MgCCR4 in the short-term independent evolutionary history. It suggests that the asymmetric evolution of monolignol biosynthetic genes might be accelerated at the early stage after WGD and declined in the later period. The inference agrees with the observation in the yeast WGD .
Co-regulation of genes involved in monolignol biosynthesis and closely related pathways
The transcriptome data of various M. × giganteus organs make it possible to explore the underlying gene transcription regulatory mechanisms on a border range of genes and independent samples. Firstly, we paid our attention to the functional relationship of co-expressed genes. The genes that showed significantly positively correlated expression patterns (Spearman’s correlation coefficients ≥ 0.4 and p value < 0.05) were regarded as the co-expressed genes. For each major monolignol biosynthetic gene, the co-expressed genes account for 3.19% to 18.42% of total expressed genes. The GO and KEGG enrichment analyses showed that these genes were significantly overrepresented in the GO terms and KEGG pathways that are related to secondary cell wall formation or share common intermediates with monolignol biosynthesis (Fisher’s exact test, Benjamini–Hochberg multiple testing corrected p value < 0.05) (for detailed results, see “Availability of data and materials”). For example, the co-expressed genes of MgHCT1 were significantly enriched in GO terms of “lignin biosynthetic process”, “phenylpropanoid biosynthetic process”, “phenylpropanoid metabolic process”, “plant-type secondary cell wall biogenesis” (Fig. 5a), and KEGG pathways of “phenylalanine metabolism”, “flavonoid biosynthesis”, “flavone and flavonol biosynthesis”, “cutin suberine and wax biosynthesis” (Fig. 5b). The consistency between gene expression patterns and gene functions indicates that the major monolignol biosynthetic genes and those genes in closely related pathways are under co-regulation in M. × giganteus.
Transcription factors of the major monolignol biosynthetic genes
TFs regulate the gene expression by specifically binding to the gene promoter regions. Based on this mechanism, the TF binding sites (TFBSs) could be predicted using the promoter sequences. To reduce false positives, we performed the expression correlation analysis between the TF genes and their target genes. The detailed results are listed in Additional file 2: Tables S6–S8.
MYB and secondary wall NAC (SWN) are two dominant TF families involved in lignin biosynthesis and other secondary cell wall formation-related pathways. By the TFBS prediction, possible MYB binding sites were significantly enriched in the promoters of major monolignol biosynthetic genes (Fisher’s exact test, p value = 0.0359, Additional file 1: Table S4). Furthermore, the expression levels of most major genes, including Mg4CL1, MgHCT1, MgHCT2, MgC3′H1, MgCCoAOMT1, MgCCoAOMT3, MgCCR1, MgCCR2, MgF5H, and MgCOMT, were significantly correlated with at least one corresponding MYB gene (Spearman’s correlation coefficients ≥ 0.4 or ≤ -0.4 and p value < 0.05, Additional file 2: Table S8). Among these, MYB61, MYB63, and MYB46 were predicted as the top three MYBs capable to directly activate the expression of multiple major monolignol biosynthetic genes in M. × giganteus (Additional file 2: Table S9). Similar to our result (Fig. 6a, b), the expression of the MYB63 gene in M. sinensis, MsSCM4, exhibited a strong positive correlation with MsHCT, and the heterologous expression of MsSCM4 in Nicotiana benthamiana mesophyll cells promoted lignin deposition .
In contrast to MYB, NAC binding sites were not significantly enriched in the promoters of major monolignol biosynthetic genes (p value = 0.1, Additional file 1: Table S4). Additionally, only the expression of Mg4CL1 and MgHCT2 showed correlation with a NAC gene (Additional file 2: Table S8). The result indicates that NACs are not likely to regulate major monolignol biosynthetic genes directly. These findings are consistent with the NAC-MYB-based gene regulatory network (NAC-MYB-GRN) model demonstrated in vascular plants. In this model, NACs function as the master switches that regulate the expression of MYB genes, e.g., MYB46/83 . Then, these MYBs promote lignin biosynthesis by activating downstream MYB genes like MYB58/63 and MYB103 . The difference is that the MYB46 in M. × giganteus was predicted to function through directly activating monolignol biosynthetic genes in our study.
WRKY is another TF family considered to be involved in the regulation of secondary cell wall formation. The enrichment of WRKY binding sites (p value = 0.00132, Additional file 1: Table S4) and significant expression correlation with Mg4CL1, MgCCoAOMT3, MgC3′H1, and MgCAD were also observed in our study (Additional file 2: Table S8). In dicots and some grasses, WRKY12 represses lignin biosynthesis via SWNs [36,37,38]. In contrast, WRKY12 in M. × giganteus might have the capability to reduce the expression of MgCAD by directly binding to the promoter based on our analysis (Fig. 6c).
In recent years, ERFs were reported to activate the lignin biosynthesis in dicots [39,40,41,42]. However, few studies have revealed the function of ERFs in the secondary cell wall formation of monocots . Surprisingly, we noticed that ERF binding sites were also significantly enriched in the promoters (p value = 2.06E-04, Additional file 1: Table S4), and some ERF genes were highly correlated with the expression of MgPAL1, MgPAL5, MgCCoAOMT1, MgCCoAOMT3, MgHCT2, MgCOMT, MgF5H and MgCAD (Additional file 2: Table S8). Our study indicates that ERF may also be another important TF family involved in monolignol biosynthesis in M. × giganteus. Therefore, ERFs are promising candidates for lignin content and composition manipulation by genetic engineering approaches.
Integrated model of monolignol biosynthesis pathway and gene regulation in M. × giganteus
In this study, MgPAL1, MgPAL5, Mg4CL1, Mg4CL3, MgHCT1, MgHCT2, MgC3′H1, MgCCoAOMT1, MgCCoAOMT3, MgCCR1, MgCCR2, MgF5H, MgCOMT, and MgCAD were inferred to be the most probable major monolignol biosynthetic genes in M. × giganteus. The evidence from phylogenetic relationships, expression patterns, and expression levels was combined. Besides, the result was strongly supported by significant sequence conservation. Concordant results could also be observed in other monocots. Most maize monolignol biosynthetic genes have similar expression patterns to the major genes of M. × giganteus in the same clades (Additional file 2: Table S5). In recent years, CSE was reported to catalyze the reaction from caffeoyl CoA to caffeic acid; however, this enzyme is not always present in plants . By aligning the switchgrass and rice CSE genes to the genomes of M. sinensis and M. sacchariflorus, as well as the transcriptome assembly of M. × giganteus, no CSE gene was found.
The involvement of isozymes is common in the monolignol biosynthesis pathway. Modification of a single gene may have little effect on the lignin content and composition. Furthermore, the presence of recent genus-specific WGD events makes things more complicated. Elucidation of the TFs that control the expression of monolignol biosynthetic genes thus becomes another important topic due to its ability to regulate multiple genes. Through TFBS prediction and expression correlation analysis, the TFs from MYB, WRKY, and ERF families were estimated to function by directly binding to the promoters of monolignol biosynthetic genes in M. × giganteus (Additional file 2: Table S9). Among these TFs, MYB61 and MYB63, which were estimated to be the dominant MYBs involved in the direct regulation of monolignol synthetic genes in M. × giganteus, were also reported in other monocots and dicots [44, 45]. Based on these results, an integrated model of the monolignol biosynthesis pathway and gene regulation in M. × giganteus was summarized (Fig. 1a, b).
However, sequence-based and expression-based approaches are the primary evidence in this work. The catalytic efficiency and substrate affinity should be further determined to figure out the actual contribution and preferred substrate(s) of each monolignol biosynthetic enzyme in M. × giganteus. The direct regulatory relationship between TFs and its target genes also needs confirmation using more straightforward evidence such as knockdown or knockout of the TF genes and ChIP-seq.
Independent evolutionary history accounts for the functional variations of monolignol biosynthetic genes in higher plants
Although the monolignol biosynthetic genes between monocots and dicots share the same origins, their functional variations have accumulated during the approximate 160 million-year independent evolutionary history after divergence . PAL is such a typical example. In dicots, PAL catalyzes the non-oxidative elimination of ammonia from phenylalanine to p-cinnamic acid. However, PAL from some monocots also exhibits tyrosine ammonia-lyase (TAL) activity, which directly catalyzes the reaction from tyrosine to 4-coumaric acid bypassing C4H [19,20,21,22].
Besides, the expression patterns of monolignol biosynthetic genes could also have changed even in monocots. In maize, the C4H gene Zm00001d009858 was relatively highly expressed in internodes (Additional file 2: Table S5), whereas in M. × giganteus, neither MgC4H1 nor MgC4H2 showed this trend. Based on these results, the 4-coumaric acid in M. × giganteus might be mainly synthesized from tyrosine, rather than phenylalanine. This finding implies the TF function estimated by heterologous expression should be interpreted with caution, owing to the probable differentiated regulatory mechanisms between the two species.
The WGD events in angiosperms could accelerate the evolution of monolignol biosynthetic genes. MgCCR1 and MgCCR2 are a pair of genes formed in the Miscanthus-specific WGD event. Although the expression patterns of the two genes are still similar (Fig. 3c), MgCCR1 has rapidly accumulated variations in the short term after WGD. This observation may also be appropriate in other monolignol biosynthetic genes of M. × giganteus.
Reference genome facilitates the gene expressional and functional studies
Although the expression levels of MgCCR1 and MgCCR2 were successfully distinguished using the gene-specific SNPs, the method is cloning-dependent. By taking advantage of the M. sinensis genome, this analysis could be simplified. In addition, the genome facilitated the evolutionary analysis, selective pressure analysis, and TFBS prediction of the monolignol biosynthetic genes.
However, the M. × giganteus genome is actually more complex than the M. sinensis genome. M. × giganteus is originated from the hybridization between M. sinensis and M. sacchariflorus, resulting in an allotriploid genome. That is to say, there should be approximately six similar copies in the M. × giganteus genome corresponding to a pair of alleles in S. bicolor. These genes cannot be distinguished using the M. sinensis genome only. Molecular cloning and qPCR experiment are also challenging . For further gene expressional and functional studies, a high-quality and haplotype-phased genome of M. × giganteus is necessary.
In this study, 14 genes were inferred as the major monolignol biosynthetic genes in M. × giganteus, including MgPAL1, MgPAL5, Mg4CL1, Mg4CL3, MgHCT1, MgHCT2, MgC3′H1, MgCCoAOMT1, MgCCoAOMT3, MgCCR1, MgCCR2, MgF5H, MgCOMT, and MgCAD. Furthermore, the TFs from MYB, WRKY, and ERF families were predicted to directly regulate the expression of these major monolignol biosynthetic genes by binding to their promoters. Based on these results, an integrated model of the monolignol biosynthesis pathway and the underlying regulatory mechanism was summarized. This study provides essential information for understanding the genetic, transcriptional, and regulatory landscape of the monolignol biosynthesis pathway in M. × giganteus. Moreover, a list of potential gene candidates was identified for genetic improvement of lignocellulosic biomass quality by manipulating the lignin content and composition.
Plant materials and sampling
M. × giganteus rhizomes were collected from the Miscanthus Resources Garden of Wuhan University at Ezhou, China (30°21′07′’N, 114°42′55′’E) and transplanted to a greenhouse at Wuhan University. When the plants were grown to the eight- to ten-leaf stage, various vegetative organ samples were collected for the molecular cloning and quantification of monolignol biosynthetic genes, including the first, third and fifth fully expanded leaves from the top to the bottom of the plants (namely L1, L3, and L5, respectively), sheaths (S), roots (R), rhizome buds (B), nodes (N), and the first to second, third to fourth, fifth to sixth internodes from the bottom to the top of the plants (IN12, IN34, IN56, respectively). After removal from the plants, the samples were washed and frozen in liquid nitrogen immediately for RNA extraction.
RNA extraction and cDNA synthesis
The samples were ground to fine powders in liquid nitrogen using chilled mortars. Total RNA extraction and genomic DNA (gDNA) removal were performed with an RNAprep Pure Plant kit (DP432, TIANGEN Biotech, Beijing, China) following the manufacturer’s instruction. RNA integrity was assessed by 1.2% agarose gel electrophoresis and a NanoDrop 2000/2000c spectrophotometer (Thermo Scientific, Waltham, USA).
For molecular cloning experiments, the cDNAs were synthesized using M-MLV Reverse Transcriptase (M1701, Promega, Madison, USA). In each reaction, 10 μl of RNA and 2 μl of Oligo(dT)15 primer (C1101, Promega) were mixed, denatured at 70 °C for 5 min and cooled on ice immediately to open the secondary structure of RNA. The mixture was added with 5 μl of 5× M-MLV buffer, 2 μl of RNase-free ddH2O, 1 μl of M-MLV Reverse Transcriptase and then incubated on a LifePro Thermal Cycler (BIOER, Hangzhou, China) at 42 °C for 1.5 h. For quantitative PCR (qPCR) experiments, the cDNAs were synthesized using a FastQuant RT Kit (KR106, TIANGEN Biotech). To minimize gDNA contamination, we treated the cDNAs with gDNase at 42 °C for 3 min once again, then mixed it with 2 μl of 10× Fast RT Buffer, 1 μl of RT Enzyme Mix, 2 μl of FQ-RT Primer Mix, and 5 μl of RNase-free ddH2O. The reverse transcription reactions were performed at 42 °C for 15 min on the thermal cycler. All products were denatured at 95 °C for 3 min to inactivate the reverse transcriptase before storage at − 20 °C.
Molecular cloning of monolignol biosynthetic genes
The cDNAs from various samples were mixed and diluted with nine volumes of ddH2O as the PCR template. Primers were designed based on the sequences obtained by rapid amplification of cDNA ends (RACE) or the transcriptome assembly of five Miscanthus species we published in the previous study , and orthologs in closely related species (Additional file 1: Table S1). To avoid mutations introduced by PCR, we used a high-fidelity DNA polymerase KOD-Plus-Neo (KOD-401, TOYOBO, Osaka, Japan) for amplification. For each reaction, the mixture consisted of 5 μl of 10 × PCR Buffer, 5 μl of dNTPs, 3 μl of MgSO4, 1.5 μl of each primer (10 μM), 2 μl of cDNA template, 2 μl of dimethylsulphoxide (DMSO), 29 μl of ddH2O and 1 μl of KOD-Plus-Neo in a total volume of 50 μl. PCRs were carried out on the thermal cycler using the two-step or three-step method based on the criteria described in our previous study . For two-step PCR, the program was set as follows: initial denaturation at 94 °C for 2 min, 36 cycles of denaturation at 98 °C for 10 s, annealing and extension at 68 °C for 4 min. The annealing and extension of three-step PCR were modified as: 30 s at the minimum melting temperature of the primer pair and 3.5 min at 68 °C. After that, deoxyadenosine residues were added to the blunt 3′-end of the amplicons by mixing 1 μl of Taq DNA polymerase (EP0405, Thermo Scientific) to the PCR products and incubating at 72 °C for 30 min. The PCR fragments were purified from 2% agarose gel using AxyPrep DNA Gel Extraction Kit (AP-GX-250, Axygen, CA, USA). The purified fragments were ligated to pGEM-T vectors (A3600, Promega) with T4 ligase at 16 °C for 12 h and transformed into Trans5α Chemically Competent Cells (CD201-02, TransGen Biotech, Beijing, China). The positive clones that harboring the recombinant plasmids were identified by blue-white screening and colony PCR with corresponding primer pairs. The insert fragments were sequenced by Sanger sequencing.
Phylogenetic and genome synteny analyses
The genomes and gene annotations of Amborella trichopoda (AMTR1.0), Arabidopsis thaliana (Araport11), Z. mays (NCBI B73_RefGen_v4), S. bicolor (Sorghum_bicolor_NCBIv3), and M. sinensis (v7.1 DOE-JGI, https://phytozome.jgi.doe.gov/) were downloaded for phylogenetic and genome synteny analyses. The possible monolignol biosynthetic genes in these species were identified by BLASTP (NCBI BLAST + version 2.7.1) . To explore whether CSE genes are present in the M. × giganteus, we aligned the CSE genes from Panicum virgatum (v1.1 DOE-JGI, https://phytozome.jgi.doe.gov/) and Oryza sativa to the genomes of M. sinensis and M. sacchariflorus (NCBI Msac_v3) using BLASTN (NCBI BLAST + version 2.7.1). The protein sequences of each enzyme were aligned together using MAFFT (version 7.453)  with the method “--localpair” and the maximum iterative refinement of 1000 (--maxiterate 1000). After alignment, the phylogenetic trees were constructed by IQ-TREE (version 2.0-rc2)  with the parameters “-B 1000 --bnni” and illustrated using FigTree (version 1.4.3, https://tree.bio.ed.ac.uk/software/figtree/).
The python version of MCscan in the JCVI package (version 1.0.5+3.g843d2f9)  was utilized to intuitively visualize the duplication events of monolignol biosynthetic genes and the orthology relationships across different species. Genome synteny blocks were identified between Amborella trichopoda versus Arabidopsis thaliana, and Amborella trichopoda versus S. bicolor using the default parameters “--cscore = 0.7, --dist = 20, --min_size = 4”. The orthologous relationships of the monolignol biosynthetic genes between species were highlighted in the macro- and micro-synteny plots.
Primers were designed using Oligo Primer Analysis Software (version 7.60)  based on the cDNA sequences of cloned monolignol biosynthetic genes. While for those genes failed to be cloned, the transcriptome assembly was used for primer design. According to our previous study, eEF-1a and UBQ were selected as the reference gene combination for inter-sample normalization . The sequences of qPCR primers are listed in Additional file 1: Table S2. To ensure the reliability of qPCR experiments, we accessed the amplification efficiency and specificity of each primer pair by standard curve analysis (Additional file 1: Table S3) and 2% agarose gel electrophoresis, respectively (Additional file 1: Fig. S1). Reaction mixtures were prepared with the SuperReal PreMix Plus Kit (FP205-02, TIAGEN Biotech), containing 10 μl of 2× SuperReal PreMix Plus (with SYBR Green I), 2 μl of 50 × ROX Reference Dye for fluorescence signal normalization, 2 μl of three- to ten-fold diluted cDNA template, 0.6 μl of each primer (10 μM) and 4.8 μl of RNase-free ddH2O. qPCR experiments were performed on a StepOne Real-Time PCR System (Applied Biosystems, Waltham, USA) with the program: initial denaturation at 95 °C for 15 min, 40 cycles of denaturation at 95 °C for 15 s, followed with annealing and extension at 60 °C for 1 min. Additional melting curve analysis was conducted for each reaction to assess the amplification specificity (Additional file 1: Fig. S1).
Before quantification, the fluorescence thresholds of each gene across plates were adjusted to the same value manually. Then, the relative expression levels to L1 were calculated using the efficiency-corrected − ΔΔCt method . The results were illustrated in bar charts with the R package ggplot2 (version 3.3.0) .
Relative expression pattern clustering and transcriptome analysis
Relative expression pattern clustering and transcriptome analysis were combined to identify the major monolignol biosynthetic genes in M. × giganteus. The relative expression levels in leaves, roots, and internodes were log2 transformed and clustered using the R package pheatmap (version 1.0.12)  with the default parameters.
Raw RNA-seq data of M. × giganteus were collected from NCBI BioProject (PRJNA183625, 17 samples)  and our previous study (NCBI SRA accession number: SRR1734721, 1 sample) . After quality filtering and adaptor trimming with fastp (version 0.20.0) , the clean reads from different samples were concatenated together and assembled using Trinity (version 2.8.6) . The completeness and continuity of the assembly were assessed by the overall mapping rate using bowtie2 (version 22.214.171.124)  and the “contig N50 of the most highly expressed genes that represent 90% of the total normalized expression” (E90N50). To annotate corresponding genes in the assembly, we aligned the protein sequences of monolignol biosynthetic genes and the coding sequences of transcription factor genes to the longest transcript of the assembly using BLASTP and BLASTN, respectively. The best hits were then selected.
In the quantification analysis, only the 17 samples in the PRJNA183625 were selected to minimize the batch effect. The gene-level expression across these samples was determined and normalized using the Perl script “abundance_estimates_to_matrix.pl” in Trinity with the parameters “--est_method RSEM, --cross_sample_norm TMM”. MgCCR1 and MgCCR2 are two genes formed in the genus-specific whole-genome duplication (WGD) event. Although they were assembled into one gene in the transcriptome assembly due to high similarity, their expression could be distinguished by the sequencing depths of gene-specific SNPs. Firstly, 15 gene-specific SNPs were identified by multiple sequence alignment and written to a Variant Call Format (VCF) file. Then the sequencing depth of each SNP was counted using ASEReadCounter in GATK (version 126.96.36.199) . Finally, the expression of each gene was calculated by multiplying the total expression level and the average depth proportion of gene-specific SNP. The absolute expression levels of isozyme genes and qPCR reference genes were illustrated with bar charts using ggplot2.
The possible functions of assembled genes were annotated using the eggNOG-mapper (version 2.0.1-14-gbf04860) [62, 63]. Only the genes assigned to Viridiplantae were kept for the downstream functional enrichment analyses. The expression correlation between genes was determined using Spearman’s correlation coefficient. Then Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses were performed and visualized using clusterProfiler (version 3.10.1) .
Calculation of percent identities and Ka/Ks ratios
The monolignol biosynthetic gene orthologs between M. × giganteus and S. bicolor were compared to determine whether the major genes have higher sequence conservation and undergo stronger purification selection. The genes that failed to be cloned were replaced by the orthologs in M. sinensis. The protein sequences were aligned using MAFFT as described above. The results were used for guiding the codon alignments of coding sequences with PAL2NAL (version 14) . The percent identities were calculated with a custom Python script and compared between the major genes and non-major genes using single-tailed Wilcoxon’s rank-sum test. The selective pressures were measured by nonsynonymous to synonymous substitution rate (Ka/Ks) ratios using the KaKs_Calculator (version 2.0) .
Transcription factor analysis
The 500-bp upstream of each gene was regarded as the promoter. The DNA sequence of each region was extracted from the M. sinensis genome with a custom Python script. The TFBS were predicted by PlantRegMap  using maize transcription factors as the targets. The motifs on both positive and negative strands were taken into consideration. The enrichment of transcription factors was determined by single-tailed Fisher’s exact test.
The Spearman’s correlation tests, Wilcoxon’s rank-sum test, and Fisher’s exact tests were performed in R (version 3.5.3) using the functions cor.test, wilcox.test, and phyper, respectively.
Availability of data and materials
The transcriptome assembly, normalized expression matrix, functional annotation results, and the custom Shell, Python, and R scripts of all the bioinformatics analyses described in this article are available at the GitHub repository: https://github.com/zengxiaofei/monolignol-biosynthesis.
Cinnamic acid 4-hydroxylase
Caffeic acid/5-hydroxyconiferaldehyde 3/5-O-methyltransferase
4-Hydroxycinnamate: CoA ligase
Hydroxycinnamoyl CoA: shikimate hydroxycinnamoyl transferase
Coumaroyl shikimate 3′-hydroxylase
Caffeoyl CoA 3-O-methyltransferase
Cinnamoyl CoA reductase
Cinnamyl alcohol dehydrogenase
Ferulic acid/coniferaldehyde 5-hydroxylase
Caffeoyl shikimate esterase
Kyoto Encyclopedia of Genes and Genomes
Transcription factor binding site
Secondary wall NAC
Chae WB, Hong SJ, Gifford JM, Rayburn AL, Sacks EJ, Juvik JA. Plant morphology, genome size, and SSR markers differentiate five distinct taxonomic groups among accessions in the genus Miscanthus. GCB Bioenergy. 2014;6(6):646–60.
Heaton EA, Dohleman FG, Long SP. Meeting US biofuel goals with less land: the potential of Miscanthus. Glob Change Biol. 2008;14(9):2000–14.
Cadoux S, Ferchaud F, Demay C, Boizard H, Machet J-M, Fourdinier E, Preudhomme M, Chabbert B, Gosse G, Mary B. Implications of productivity and nutrient requirements on greenhouse gas balance of annual and perennial bioenergy crops. GCB Bioenergy. 2014;6(4):425–38.
Zhang B, Hastings A, Clifton-Brown JC, Jiang D, Faaij APC. Modeled spatial assessment of biomass productivity and technical potential of Miscanthus × giganteus, Panicumvirgatum L., and Jatropha on marginal land in China. GCB Bioenergy. 2020;12(5):328–45.
Dey P, Pal P, Kevin JD, Das DB. Lignocellulosic bioethanol production: prospects of emerging membrane technologies to improve the process—a critical review. Rev Chem Eng. 2020;36(3):333.
Banerjee S, Mudliar S, Sen R, Giri B, Satpute D, Chakrabarti T, Pandey RA. Commercializing lignocellulosic bioethanol: technology bottlenecks and possible remedies. Biofuels Bioprod Biorefin. 2010;4(1):77–93.
Li M, Pu Y, Ragauskas AJ. Current understanding of the correlation of lignin structure with biomass recalcitrance. Front Chem. 2016;4:45.
Kim D. Physico-chemical conversion of lignocellulose: inhibitor effects and detoxification strategies: a mini review. Molecules. 2018;23(2):309.
Jönsson LJ, Martín C. Pretreatment of lignocellulose: Formation of inhibitory by-products and strategies for minimizing their effects. Bioresour Technol. 2016;199:103–12.
Park S-H, Mei C, Pauly M, Ong RG, Dale BE, Sabzikar R, Fotoh H, Nguyen T, Sticklen M. Downregulation of maize cinnamoyl-coenzyme a reductase via RNA interference technology causes brown midrib and improves ammonia fiber expansion-pretreated conversion into fermentable sugars for biofuels. Crop Sci. 2012;52(6):2687–701.
Fornalé S, Capellades M, Encina A, Wang K, Irar S, Lapierre C, Ruel K, Joseleau JP, Berenguer J, Puigdomènech P, et al. Altered lignin biosynthesis improves cellulosic bioethanol production in transgenic maize plants down-regulated for cinnamyl alcohol dehydrogenase. Mol plant. 2012;5(4):817–30.
Fu C, Mielenz JR, Xiao X, Ge Y, Hamilton CY, Rodriguez M, Chen F, Foston M, Ragauskas A, Bouton J, et al. Genetic manipulation of lignin reduces recalcitrance and improves ethanol production from switchgrass. Proc Natl Acad Sci. 2011;108(9):3803.
Fu C, Xiao X, Xi Y, Ge Y, Chen F, Bouton J, Dixon RA, Wang Z-Y. Downregulation of cinnamyl alcohol dehydrogenase (CAD) leads to improved saccharification efficiency in switchgrass. BioEnergy Res. 2011;4(3):153–64.
Xu B, Escamilla-Treviño LL, Sathitsuksanoh N, Shen Z, Shen H, Percival Zhang Y-H, Dixon RA, Zhao B. Silencing of 4-coumarate:coenzyme A ligase in switchgrass leads to reduced lignin content and improved fermentable sugar yields for biofuel production. New Phytol. 2011;192(3):611–25.
Li M, Si S, Hao B, Zha Y, Wan C, Hong S, Kang Y, Jia J, Zhang J, Li M, et al. Mild alkali-pretreatment effectively extracts guaiacyl-rich lignin for high lignocellulose digestibility coupled with largely diminishing yeast fermentation inhibitors in Miscanthus. Bioresour Technol. 2014;169:447–54.
Xu N, Zhang W, Ren S, Liu F, Zhao C, Liao H, Xu Z, Huang J, Li Q, Tu Y, et al. Hemicelluloses negatively affect lignocellulose crystallinity for high biomass digestibility under NaOH and H2SO4 pretreatments in Miscanthus. Biotechnol Biofuels. 2012;5(1):58.
Dixon RA, Barros J. Lignin biosynthesis: old roads revisited and new roads explored. Open Biol. 2019;9(12):190215.
Barros J, Escamilla-Trevino L, Song L, Rao X, Serrani-Yarce JC, Palacios MD, Engle N, Choudhury FK, Tschaplinski TJ, Venables BJ, et al. 4-Coumarate 3-hydroxylase in the lignin biosynthesis pathway is a cytosolic ascorbate peroxidase. Nat Commun. 2019;10(1):1994.
Rosler J, Krekel F, Amrhein N, Schmid J. Maize phenylalanine ammonia-lyase has tyrosine ammonia-lyase activity. Plant Physiol. 1997;113(1):175–9.
Cass CL, Peraldi A, Dowd PF, Mottiar Y, Santoro N, Karlen SD, Bukhman YV, Foster CE, Thrower N, Bruno LC, et al. Effects of PHENYLALANINE AMMONIA LYASE (PAL) knockdown on cell wall composition, biomass digestibility, and biotic and abiotic stress responses in Brachypodium. J Exp Bot. 2015;66(14):4317–35.
Barros J, Serrani-Yarce JC, Chen F, Baxter D, Venables BJ, Dixon RA. Role of bifunctional ammonia-lyase in grass cell wall biosynthesis. Nat Plants. 2016;2(6):16050.
Jun SY, Sattler SA, Cortez GS, Vermerris W, Sattler SE, Kang C. Biochemical and structural analysis of substrate specificity of a phenylalanine ammonia-lyase. Plant Physiol. 2018;176(2):1452–68.
Ha CM, Escamilla-Trevino L, Yarce JCS, Kim H, Ralph J, Chen F, Dixon RA. An essential role of caffeoyl shikimate esterase in monolignol biosynthesis in Medicago truncatula. Plant J. 2016;86(5):363–75.
Rao X, Chen X, Shen H, Ma Q, Li G, Tang Y, Pena M, York W, Frazier TP, Lenaghan S, et al. Gene regulatory networks for lignin biosynthesis in switchgrass (Panicum virgatum). Plant Biotechnol J. 2019;17(3):580–93.
Shen H, Mazarei M, Hisano H, Escamilla-Trevino L, Fu C, Pu Y, Rudis MR, Tang Y, Xiao X, Jackson L, et al. A genomics approach to deciphering lignin biosynthesis in switchgrass. Plant Cell. 2013;25(11):4342–61.
Guillaumie S, San-Clemente H, Deswarte C, Martinez Y, Lapierre C, Murigneux A, Barrière Y, Pichon M, Goffner D. MAIZEWALL. Database and developmental gene expression profiling of cell wall biosynthesis and assembly in maize. Plant Physiol. 2007;143(1):339–63.
Yang F, Li W, Jiang N, Yu H, Morohashi K, Ouma WZ, Morales-Mantilla DE, Gomez-Cano FA, Mukundi E, Prada-Salcedo LD, et al. A maize gene regulatory network for phenolic metabolism. Mol plant. 2017;10(3):498–515.
Jardim-Messeder D, da Franca ST, Fonseca JP, Junior JN, Barzilai L, Felix-Cordeiro T, Pereira JC, Rodrigues-Ferreira C, Bastos I, da Silva TC, et al. Identification of genes from the general phenylpropanoid and monolignol-specific metabolism in two sugarcane lignin-contrasting genotypes. Mol Genet Genomics. 2020;295(3):717–39.
Glover N, Dessimoz C, Ebersberger I, Forslund SK, Gabaldón T, Huerta-Cepas J, Martin M-J, Muffato M, Patricio M, Pereira C, et al. Advances and applications in the quest for orthologs. Mol Biol Evol. 2019;36(10):2157–64.
Zeng X, Sheng J, Zhu F, Zhao L, Hu X, Zheng X, Zhou F, Hu Z, Diao Y, Jin S. Differential expression patterns reveal the roles of cellulose synthase genes (CesAs) in primary and secondary cell wall biosynthesis in Miscanthus × giganteus. Ind Crops Prod. 2020;145:112129.
Raes J, Rohde A, Christensen JH, Van de Peer Y, Boerjan W. Genome-wide characterization of the lignification toolbox in Arabidopsis. Plant Physiol. 2003;133(3):1051–71.
Pegueroles C, Laurie S, Albà MM. Accelerated evolution after gene duplication: a time-dependent process affecting just one copy. Mol Biol Evol. 2013;30(8):1830–42.
Scannell DR, Wolfe KH. A burst of protein sequence evolution and a prolonged period of asymmetric evolution follow gene duplication in yeast. Genome Res. 2008;18(1):137–47.
Golfier P, Volkert C, He F, Rausch T, Wolf S. Regulation of secondary cell wall biosynthesis by a NAC transcription factor from Miscanthus. Plant Direct. 2017;1(5):e00024.
Ohtani M, Demura T. The quest for transcriptional hubs of lignin biosynthesis: beyond the NAC-MYB-gene regulatory network model. Curr Opin Biotechnol. 2019;56:82–7.
Yang L, Zhao X, Yang F, Fan D, Jiang Y, Luo K. PtrWRKY19, a novel WRKY transcription factor, contributes to the regulation of pith secondary wall formation in Populus trichocarpa. Sci Rep. 2016;6(1):18643.
Gallego-Giraldo L, Shadle G, Shen H, Barros-Rios J, Fresquet Corrales S, Wang H, Dixon RA. Combining enhanced biomass density with reduced lignin level for improved forage quality. Plant Biotechnol J. 2016;14(3):895–904.
Wang H, Avci U, Nakashima J, Hahn MG, Chen F, Dixon RA. Mutation of WRKY transcription factors initiates pith secondary wall formation and increases stem biomass in dicotyledonous plants. Proc Natl Acad Sci USA. 2010;107(51):22338–43.
Guo W, Jin L, Miao Y, He X, Hu Q, Guo K, Zhu L, Zhang X. An ethylene response-related factor, GbERF1-like, from Gossypium barbadense improves resistance to Verticillium dahliae via activating lignin synthesis. Plant Mol Biol. 2016;91(3):305–18.
Liu Y, Wei M, Hou C, Lu T, Liu L, Wei H, Cheng Y, Wei Z. Functional characterization of Populus PsnSHN2 in coordinated regulation of secondary wall components in tobacco. Sci Rep. 2017;7(1):42.
Ma R, Xiao Y, Lv Z, Tan H, Chen R, Li Q, Chen J, Wang Y, Yin J, Zhang L, et al. AP2/ERF transcription factor, Ii049, positively regulates lignan biosynthesis in Isatis indigotica through activating salicylic acid signaling and lignan/lignin pathway genes. Front Plant Sci. 2017;8:1361.
Zeng J-K, Li X, Xu Q, Chen J-Y, Yin X-R, Ferguson IB, Chen K-S. EjAP2-1, an AP2/ERF gene, is a novel regulator of fruit lignification induced by chilling injury, via interaction with EjMYB transcription factors. Plant Biotechnol J. 2015;13(9):1325–34.
Wuddineh WA, Mazarei M, Turner GB, Sykes RW, Decker SR, Davis MF, Stewart CN Jr. Identification and molecular characterization of the switchgrass AP2/ERF transcription factor superfamily, and overexpression of PvERF001 for improvement of biomass characteristics for biofuel. Front Bioeng Biotechnol. 2015;3:101–101.
Nakano Y, Yamaguchi M, Endo H, Rejab NA, Ohtani M. NAC-MYB-based transcriptional regulation of secondary cell wall biosynthesis in land plants. Front Plant Sci. 2015;6:288.
Rao X, Dixon RA. Current models for transcriptional regulation of secondary cell wall biosynthesis in grasses. Front Plant Sci. 2018;9:399.
Kumar S, Stecher G, Suleski M, Hedges SB. TimeTree: a resource for timelines, timetrees, and divergence times. Mol Biol Evol. 2017;34(7):1812–9.
Zeng X, Cheng N, Zheng X, Diao Y, Fang G, Jin S, Zhou F, Hu Z. Molecular cloning and characterization of two manganese superoxide dismutases from Miscanthus × giganteus. Plant Cell Rep. 2015;34(12):2137–49.
Sheng J, Zheng X, Wang J, Zeng X, Zhou F, Jin S, Hu Z, Diao Y. Transcriptomics and proteomics reveal genetic and biological basis of superior biomass crop Miscanthus. Sci Rep. 2017;7(1):13777.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. BLAST+: architecture and applications. BMC Bioinf. 2009;10(1):421.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.
Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, Lanfear R. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020;37(5):1530–4.
Tang H, Bowers JE, Wang X, Ming R, Alam M, Paterson AH. Synteny and collinearity in plant genomes. Science. 2008;320(5875):486.
Rychlik W. OLIGO 7 primer analysis software. In: Yuryev A, editor. PCR primer design. Totowa, NJ: Humana Press; 2007. p. 35–59.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29(9):e45–e45.
Wickham H. ggplot2: elegant graphics for data analysis. Springer; 2016.
Kolde R, Kolde MR. Package ‘pheatmap’. R Package 2015.
Barling A, Swaminathan K, Mitros T, James BT, Morris J, Ngamboma O, Hall MC, Kirkpatrick J, Alabady M, Spence AK, et al. A detailed gene expression study of the Miscanthus genus reveals changes in the transcriptome associated with the rejuvenation of spring rhizomes. BMC Genomics. 2013;14(1):864.
Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):i884–90.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.
DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, del Angel G, Rivas MA, Hanna M, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43(5):491–8.
Huerta-Cepas J, Forslund K, Coelho LP, Szklarczyk D, Jensen LJ, von Mering C, Bork P. Fast genome-wide functional annotation through orthology assignment by eggNOG-mapper. Mol Biol Evol. 2017;34(8):2115–22.
Huerta-Cepas J, Szklarczyk D, Heller D, Hernández-Plaza A, Forslund SK, Cook H, Mende DR, Letunic I, Rattei T, Jensen Lars J, et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 2018;47(D1):D309–14.
Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Suyama M, Torrents D, Bork P. PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006;34(suppl_2):W609–12.
Wang D, Zhang Y, Zhang Z, Zhu J, Yu J. KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genomics proteome Bioinform. 2010;8(1):77–80.
Tian F, Yang D-C, Meng Y-Q, Jin J, Gao G. PlantRegMap: charting functional regulatory maps in plants. Nucleic Acids Res. 2019;48(D1):D1104–13.
We thank Danni Liu for the technical assistance of transcriptome analysis. The transcriptome assembly and quantification were supported by the Center for Computational Science and Engineering of Southern University of Science and Technology. The genome data of M. sinensis were produced by the US Department of Energy Joint Genome Institute (DOE-JGI).
This study was financially supported by the National Natural Science Foundation of China [Grant No. 31571740], the National High-tech R&D Program [Grant No. 2012AA101801], and the Natural Science Foundation of Hubei Province [Grant No. 2013CFA103].
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zeng, X., Sheng, J., Zhu, F. et al. Genetic, transcriptional, and regulatory landscape of monolignol biosynthesis pathway in Miscanthus × giganteus. Biotechnol Biofuels 13, 179 (2020). https://doi.org/10.1186/s13068-020-01819-4
- Miscanthus × giganteus
- Monolignol biosynthesis pathway
- Monolignol biosynthetic genes
- Transcription factors
- Regulatory mechanism
- Transcriptome analysis
- Genetic engineering