DNA methylation plays an important role in the growth and development of plant, and would change under different environments. In this study, 5-methyl cytosine (5mC) content and methylation level exhibited tissue specificity in genomic DNA of wheat seedling, and increased significantly in leaf along with the increase of water deficit, which was especially significant in leaf of wheat AK58. Full-methylation might dominate in genomic DNA of wheat seedling, the increase of full-methylation level under water deficit was significantly higher than that of hemi-methylation level. Under water deficit, DNA methylation of wheat seedling showed significant polymorphism, this polymorphism was always higher in root, especially was higher in root of wheat AK58. Further analysis appeared that changes of DNA methylation in wheat seedling took methylation as principle and demethylation as supplement under water deficit. Therefore, under water deficit, the degree, level and polymorphism of DNA methylation in wheat seedling showed tissue specificity and species specificity, and were higher in wheat AK58 compared with those of wheat XM13, perhaps wheat AK58 could more rapidly respond to water deficit by changes of DNA methylation, which would contribute to reveal molecular mechanism of wheat adapting to water deficit.

In this study, ITS, ITS2, matK, rbcL and psbA-trnH in Rehmannia were successfully amplified and sequenced, but some ITS sequences need to be proofread according to ITS2 sequences. Compared with rbcL, matK and psbA-trnH, ITS and ITS2 had higher mutation rate and more information sites, and ITS2 had higher interspecific diversity and lower intraspecific variation in Rehmannia, but the interspecific genetic variation of rbcL and matK was lower. Furthermore, the obvious barcoding gap was found in psbA-trnH or ITS2 + psbA-trnH, and the overlap between interspecific and intraspecific variation of ITS, ITS2 or matK was less. In addition, the phylogenetic tree based on ITS or ITS2 indicated that R. glutinosa, R. chingii or R. henryi with obvious monophyly could be successfully identified, but R. piasezkii and R. elata were clustered into one branch, R. solanifolia could not be distinguished from R. glutinosa, and R. chingii was closer to R. henryi. In phylogenetic tree based on psbA-trnH or ITS2 + psbA-trnH, cultivars and wild varieties of R. glutinosa could be distinguished, were clearly separated from other Rehmannia species, and cultivars or wild varieties of R. glutinosa could be also distinguished by matK. Taken together, ITS2 has great potential in systematic study and species identification of Rehmannia, the combination of ITS2 and psbA-trnH might be the most suitable DNA barcode for Rehmannia species.

Differential metabolites between tuberous roots from cultivated variety (ZP) and wild variety (YS) of Rehmannia glutinosa were analyzed by widely targeted metabolomics, and annotated to KEGG pathways. 228 secondary metabolites (SM) in ZP and YS were detected, of which 58 were differential metabolites (DM), including 41 flavonoids, 10 phenolic acids, 3 terpenoids, 2 alkaloids and 2 others, and 170 were unchanged; Among 58 DMs, 44 (75.9%) were up-regulated in YS, of which 30 were unique to YS, while 14 (24.1%) were down-regulated in YS, of which 10 were unique to ZP; Among flavonoids, 33 (80.5%) were more highly expressed in YS than in ZP; Among phenolic acids, 7 (70%) were more highly expressed in YS than in ZP; 12 of 58 DMs were annotated into 17 types of KEGG pathways. Among them, benzoic acid and p-Coumaryl alcohol were up-regulated in YS, and annotated into 10 pathways (58.8%) and 4 pathways (23.5%), respectively. In addition, much of DMs possess various pharmacological effects. These results indicated better quality of YS than ZP and the necessity of YS domestication. Taken together, this study will provide a reference for the scientific introduction, comprehensive development and utilization of wild Rehmannia glutinosa.

Rehmannia glutinosa has important medicinal value; terpenoid is one of the main active components in R. glutinosa. In this study, iTRAQ technique was used to analyze the relative abundance of proteins in roots of R. glutinosa, and 6,752 reliable proteins were quantified. GO enrichment results indicated that most proteins were involved in metabolic process or cellular process, 57.63% proteins had catalytic activity, and 65.80% proteins were enriched in membrane-bounded organelle. In roots of R. glutinosa, there were 38 KEGG enrichments with significance, more DEPs were found in some pathways, especially the proteasome pathway and TCA cycle with 15.0% DEPs between elongation stage and expansion stage of roots. Furthermore, five KEGG pathways of terpenoid synthesis were found. Most prenyltransferases belong to FPP/GGPP synthase family, involved in terpenoid backbone biosynthesis, and all interacted with biotin carboxylase CAC2. Compared with that at the elongation stage, many prenyltransferases exhibited higher expression at the expansion stage or maturation stage of roots. In addition, eight FPP/GGPP synthase encoding genes were cloned from R. glutinosa, namely FPPS, FPPS1, GGPS, GGPS3, GGPS4, GGPS5, GPPS and GPPS2, introns were also found in FPPS, FPPS1, GGPS5 and GGPS2, and FPP/GPP synthases were more conservative in organisms, especially in viridiplantae, in which the co-occurrence of GPPS or GPPS2 was significantly higher in plants. Further analysis found that FPP/GGPP synthases of R. glutinosa were divided into three kinds, GGPS, GPPS and FPPS, and their gene expression was significantly diverse in different varieties, growth periods, or tissues of R. glutinosa. Compared with that of GGPS, the expression of GPPS and FPPS was much higher in R. glutinosa, especially at the expansion stage and maturation stage. Thus, the synthesis of terpenoids in roots of R. glutinosa is intricately regulated and needs to be further studied.

Rehmannia glutinosa root contains many compounds with important medicinal properties and nutritional benefits, but only more than 140 compounds have been reported so far. Many other compounds and their accumulation and metabolic networks during its development remain unclear. In order to clarify them, its metabolic profiles at three different developmental stages were analyzed using untargeted LC-MS analysis. Multivariate analysis revealed that 434 metabolites differently accumulated in its different stages, suggesting different change trends. The metabolites having the same trend share common metabolic pathways, the metabolites showing increasing contents during its development have medical and nutritional values, and some mature root-specific metabolites may be better candidates for its quality control; 434 metabolites were mapped onto 111 KEGG pathways including 62 enzymes, whose increasing and decreasing patterns were shown during its development. Some metabolites complicatedly interacted with some enzymes and the top-10 pathways enriched from 111 KEGG pathways in network analysis. These findings extended the dataset of its identified compounds, and revealed that its development and quality were associated with the accumulation of different metabolites. Our work will lay the foundation for the better understanding of its chemical constituents, quality and developmental mechanism.

The objective of this research was to establish an efficient system of genetic transformation and plant regeneration from hairy roots by infecting the leaf sections and stem segments of in vitro Rehmannia glutinosa Libosch. f. hueichingensis Hsiao plantlets. Hairy roots were induced from them after co-culturing with Agrobacterium rhizogenes strain 15834 at a frequency of 32 and 29.4%, respectively. The calluses were induced from hairy roots on half-strength Murashige and Skoog medium containing 0.2 mg/l kinetin and 3.0 mg/l benzyladenine at a frequency of 100%, from which transgenic shoots and plantlets were developed. Transgenic plantlets did not have differences in morphology except the shortened internodes and an increase in adventitious root formation compared to wild-type plants. PCR and Southern-blot analyses confirmed that rolB gene of TL-DNA was inserted in the genome of transformed hairy roots and plantlets. RT-PCR analysis and opine paper electrophoresis revealed that rolB gene was expressed in the transformed hairy roots and plantlets. Conclusively, transgenic hairy roots and transgenic plants of Rehmannia glutinosa Libosch. f. hueichingensis Hsiao were developed for the first time.