SUN gene is a group of key genes regulating plant growth and development. Here, SUN gene families of strawberry were identified from the genome of the diploid Fragaria vesca, and their physicochemical properties, genes structure, evolution and genes expression were also analyzed. Our results showed that there were thirty-one FvSUN genes in F. vesca and the FvSUNs encoded proteins were classified into seven groups, and the members in the same group showed high similarity in gene structures and conservative motifs. The electronic subcellular localization of FvSUNs was mainly in the nucleus. Collinearity analysis showed that the members of FvSUN gene family were mainly expanded by segmental duplication in F. vesca, and Arabidopsis and F. vesca shared twenty-three pairs of orthologous SUN genes. According to the expression pattern in different tissues shown by the transcriptome data of F. vesca, the FvSUNs gene can be divided into three types: (1) expressed in nearly all tissues, (2) hardly expressed in any tissues, and (3) expressed in special tissues. The gene expression pattern of FvSUNs was further verified by quantitative real-time polymerase chain reaction (qRT-PCR). Additionally, the seedlings of F. vesca were treated by different abiotic stresses, and the expression level of 31 FvSUNs genes were assayed by qRT-PCR. The expression of most of the tested genes was induced by cold, high salt or drought stress. Our studies may facilitate revealing the biological function and molecular mechanism of SUN genes in strawberry.
Fragaria/metabolism* ; Genes, Plant ; Stress, Physiological/genetics* ; Arabidopsis/genetics* ; Plant Development ; Gene Expression Regulation, Plant ; Plant Proteins/metabolism*

MADS-box gene family is a significant transcription factor family that plays a crucial role in regulating plant growth, development, signal transduction, and other processes. In order to study the characteristics of MADS-box gene family in Docynia delavayi (Franch.) Schneid. and its expression during different stages of seed germination, this study used seedlings at different stages of germination as materials and screened MADS-box transcription factors from the transcriptome database of D. delavayi using bioinformatics methods based on transcriptome sequencing. The physical and chemical properties, protein conservative motifs, phylogenetic evolution, and expression patterns of the MADS-box transcription factors were analyzed. Quantitative real-time PCR (qRT-PCR) was used to verify the expression of MADS-box gene family members during different stages of seed germination in D. delavayi. The results showed that 81 genes of MADS-box gene family were identified from the transcriptome data of D. delavayi, with the molecular weight distribution ranged of 6 211.34-173 512.77 Da and the theoretical isoelectric point ranged from 5.21 to 10.97. Phylogenetic analysis showed that the 81 genes could be divided into 15 subgroups, among which DdMADS27, DdMADS42, DdMADS45, DdMADS46, DdMADS53, DdMADS61, DdMADS76, DdMADS77 and DdMADS79 might be involved in the regulation of ovule development in D. delavayi. The combination of the transcriptome data and the qRT-PCR analysis results of D. delavayi seeds indicated that DdMADS25 and DdMADS42 might be involved in the regulation of seed development, and that DdMADS37 and DdMADS38 might have negative regulation effects on seed dormancy. Previous studies have reported that the MIKC^* subgroup is mainly involved in regulating flower organ development. For the first time, we found that the transcription factors of the MIKC^* subgroup exhibited a high expression level at the early stage of seed germination, so we speculated that the MIKC^* subgroup played a regulatory role in the process of seed germination. To verify the accuracy of this speculation, we selected DdMADS60 and DdMADS75 from the MIKC^* subgroup for qRT-PCR experiments, and the experimental results were consistent with the expression trend of transcriptome sequencing. This study provides a reference for further research on the biological function of D. delavayi MADS-box gene family from the perspective of molecular evolution.
MADS Domain Proteins/metabolism* ; Phylogeny ; Gene Expression Regulation, Plant ; Genes, Plant ; Transcription Factors/genetics* ; Plant Proteins/metabolism* ; Gene Expression Profiling

Dracaena marginata is a widely cultivated horticultural plant in the world, which has high ornamental and medicinal value. In this study, the whole genome of leaves from D. marginata was sequenced by Illumina HiSeq 4000 platform. The chloroplast genome were assembled for functional annotation, sequence characteristics and phylogenetic analysis. The results showed that the chloroplast genome of D. marginata composed of four regions with a size of 154 926 bp, which was the smallest chloroplast genome reported for Dracaena species to date. A total of 132 genes were identified, including 86 coding genes, 38 tRNA genes and 8 rRNA genes. Codon bias analysis found that the codon usage bias was weak and there was a bias for using A/U base endings. 46 simple sequence repeat and 54 repeats loci were detected in the chloroplast genome, with the maximum detection rate in the large single copy region and inverted repeat region, respectively. The inverted repeats boundaries of D. marginata and Dracaena were highly conserved, whereas gene location differences occurred. Phylogenetic analysis revealed that D. serrulata and D. cinnabari form a monophyletic clade, which was the closest relationship and conformed to the morphological classification characteristics. The analysis of the chloroplast genome of D. marginata provides important data basis for species identification, genetic diversity and chloroplast genome engineering of Dracaena.
Phylogeny ; Dracaena ; Genome, Chloroplast/genetics* ; Base Sequence ; Genes, Plant

This paper aimed to study the role of asparagine endopeptidase(AEP) gene in the biosynthesis mechanism of cyclic peptide compounds in Pseudostellaria heterophylla. The transcriptome database of P. heterophylla was systematically mined and screened, and an AEP gene, tentatively named PhAEP, was successfully cloned. The heterologous function verification by Nicotiana benthamiana showed that the expression of the gene played a role in the biosynthesis of heterophyllin A in P. heterophylla. Bioinformatics analysis showed that the cDNA of PhAEP was 1 488 bp in length, encoding 495 amino acids with a molecular weight of 54.72 kDa. The phylogenetic tree showed that the amino acid sequence encoded by PhAEP was highly similar to that of Butelase-1 in Clitoria ternatea, reaching 80%. The sequence homology and cyclase active site analysis revealed that the PhAEP enzyme may specifically hydrolyse the C-terminal Asn/Asp(Asx) site of the core peptide in the HA linear precursor peptide of P. heterophylla, thereby participating in the ring formation of the linear precursor peptide. The results of real-time quantitative polymerase chain reaction(RT-qPCR) showed that the expression level of PhAEP was the highest in fruits, followed by in roots, and the lowest in leaves. The heterophyllin A of P. heterophylla was detected in N. benthamiana that co-expressed PrePhHA and PhAEP genes instantaneously. In this study, the PhAEP gene, a key enzyme in the biosynthesis of heterophyllin A in P. heterophylla, has been successfully cloned, which lays a foundation for further analysis of the molecular mechanism of PhAEP enzyme in the biosynthesis of heterophyllin A in P. heterophylla and has important significance for the study of synthetic biology of cyclic peptide compounds in P. heterophylla.
Genes, vif ; Phylogeny ; Plant Leaves/genetics* ; Peptides, Cyclic ; Cloning, Molecular ; Caryophyllaceae/genetics*

The plant growth, development, and secondary metabolism are regulated by R2 R3-MYB transcription factors. This study identified the R2 R3-MYB genes in the genome of Andrographis paniculata and analyzed the chromosomal localization, gene structure, and conserved domains, phylogenetic relationship, and promoter cis-acting elements of these R2 R3-MYB genes. Moreover, the gene expression profiles of R2 R3-MYB genes under abiotic stress and hormone treatments were generated by RNA-seq and validated by qRT-PCR. The results showed that A. paniculata contained 73 R2 R3-MYB genes on 21 chromosomes. These members belonged to 34 subfamilies, 19 of which could be classified into the known subfamilies in Arabidopsis thaliana. The 73 R2 R3-MYB members included 36 acidic proteins and 37 basic proteins, with the lengths of 148-887 aa. The domains, motifs, and gene structures of R2 R3-MYBs in A. paniculata were conserved. The promoter regions of these genes contains a variety of cis-acting elements related to the responses to environmental factors and plant hormones including light, ABA, MeJA, and drought. Based on the similarity of functions of R2 R3-MYBs in the same subfamily and the transcription profiles, ApMYB13/21/35/67/73(S22) may regulate drought stress through ABA pathway; ApMYB20(S11) and ApMYB55(S2) may play a role in the response of A. paniculata to high temperature and UV-C stress; ApMYB5(S7) and ApMYB33(S20) may affect the accumulation of andrographolide by regulating the expression of key enzymes in the MEP pathway. This study provides theoretical reference for further research on the functions of R2 R3-MYB genes in A. paniculata and breeding of A. paniculata varieties with high andrographolide content.
Andrographis paniculata ; Gene Expression Regulation, Plant ; Genes, myb ; Multigene Family ; Phylogeny ; Plant Proteins/metabolism*

Artemisia Argyi Folium, a traditional Chinese medicine of important medicinal and economic value, sees increasing demand in medicinal and moxibustion product market. Screening stable and reliable reference genes for quantitative real-time PCR(qRT-PCR) is a prerequisite for the analysis of gene expression in Artemisia argyi. In this study, eight commonly used reference genes, Actin, 18s, EF-1α, GAPDH, SAND, PAL, TUA, and TUB, from the transcriptome of A. argyi, were selected as candidate genes. The expression of each gene in different tissues(roots, stems, and leaves) of A. argyi and in leaves of A. argyi after treatment with methyl jasmonate(MeJA) for different time(0, 4, 8, 12 h) was detected by qRT-PCR. Then, geNorm, NormFinder, BestKeeper, ΔCT, and RefFinder were employed to evaluate their expression stability. The results demonstrated that Actin was the most stable reference gene in different tissues and in leaves treated with MeJA, and coming in the second was SAND. Furthermore, the expression of DXS and MCT which are involved in terpenoid backbone biosynthesis was detected in different tissues and after MeJA treatment. The results showed that the expression patterns of DXS and MCT in different tissues and under MeJA treatment calculated with Actin and SAND as internal reference genes were consistent, which validated the screening results. In conclusion, Actin is the most suitable reference gene for the analysis of gene expression in different tissues of A. argyi and after MeJA treatment. This study provides valuable information for gene expression analysis in A. argyi and lays a foundation for further research on molecular mechanism of quality formation of Artemisia Argyi Folium.
Artemisia/genetics* ; Gene Expression Profiling ; Genes, Plant/genetics* ; Plant Leaves/genetics* ; Real-Time Polymerase Chain Reaction ; Reference Standards ; Transcriptome

To select suitable references gene of Polygonum multiflorum for gene expression analysis in different tissues, five candidate reference genes like Actin,GAPDH,SAND,PP2A,TIP41 were selected from the transcriptome data of P. multiflorum, then the specific primers were designed. The expression stability of the five reference genes in different tissues of P. multiflorum was analyzed by Real-time quantitative PCR through avilable analysis methods such as geNorm, NormFinder, BestKeeper, Delta CT and RefFinder, to ensure the reliability of the analysis results. The results showed that there were significant differences in the expression levels and stability of candidate genes in different tissues of P. multiflorum. Ct distribution analysis of the expression levels of candidate genes showed that the expression levels of Actin and GAPDH genes were relatively high in different tissues, while the expression levels of SAND, PP2A and TIP41 were lower. The stability of each candidate gene was analyzed by different methods. The results of geNorm analysis showed that the expression of PP2A and GAPDH was the most stable, the expression stability of SAND was the worst, the stability of PP2A was the highest in both NormFinder and Delta CT, the stability of SAND was the lowest, and the stability of Actin was the most stable in BestKeeper analysis. Through the comprehensive evaluation and analysis of the stability of candidate genes by RefFinder, it is concluded that the stability of PP2A gene is the highest, followed by GAPDH, Actin, TIP41, SAND, and SAND gene is the worst. Therefore, the PP2A gene is an ideal reference gene for the analysis of gene expression in different tissues of P. multiflorum.
Fallopia multiflora ; Gene Expression Profiling ; Gene Expression Regulation, Plant ; Genes, Plant/genetics* ; Real-Time Polymerase Chain Reaction ; Reference Standards ; Reproducibility of Results

Amana edulis is a traditional Chinese medicinal plant with low propagation coefficient. In recent years, the increasing demands of A. edulis lead to a shortage of its wild resources. In order to analyze the expression of related functional genes in A. edulis, the selection of suitable internal reference genes is crucial to improve the accuracy of experimental results. Eight genes(ACT, TUA, CYP, GAPDH, UBQ, UBI, EF1a, UBC)were chosen as candidate reference genes based on the RNA-Seq. Real-time fluorescence quantitative technique was used to detect the expression level of candidate internal reference genes in different organs(bulb, leaf, flo-wer) and stolons at different development stages of A. edulis. Then GeNorm, NormFinder, BestKeeper softwares and RefFinder website were used for a comprehensive analysis of the expression stability of the candidate genes.The results showed that among the 8 candidate reference genes, the variation range of Ct value of UBC was the smallest, and the expression level was stable, which was suitable for an reference gene. GeNorm and NormFinder software analysis showed that UBC and UBI were the optimal reference genes. BestKeeper analysis showed that CYP and UBC expression were relatively stable. Comprehensive evaluation of RefFinder website showed that UBC and UBI were the most stable genes, and ACT displayed the lowest stability in all software evaluation, indicating UBC and UBI were suitable for reference genes. Additionally, the most stable UBC, UBI and the most unstable ACT were used as internal reference genes to detect the expression of GBSS gene in A. edulis, and expression pattern of GBSS gene was the same under the calibration of UBC and UBI. The expression data of GBSS gene confirmed that UBC and UBI genes were reliable for A. edulis qRT-PCR as internal reference genes. The results would benefit future studies on related gene expression of A. edulis.
Gene Expression Profiling ; Gene Expression Regulation, Plant ; Genes, Plant/genetics* ; Real-Time Polymerase Chain Reaction ; Reference Standards

Podophyllotoxin (PTOX) is an aryl-tetralin lignan of plant origin found in some species of Podophyllum such as Dysosma versipellis, Diphylleia sinensis, and Sinopodophyllum hexandrum. Etoposide and teniposide are produced semisynthetically from PTOX and used clinically to treat several forms of cancer. As a typical representative of new drug discovery from natural products, the production of PTOX solely depends on extraction from plants, resulting in severe contradiction between supply and demand. With the advantages of unconstrained resources and eco-friendly reaction conditions, biosynthesis method has become a trend in the production of PTOX and its derivatives. In this review, we summarize the research progress of PTOX biosynthesis in plants and expound the functions of the key enzymes as well as their subcellular location. The synthetic biology for production of PTOX intermediates in a tobacco chassis is also introduced. Finally, the heterologous expression and biotransformation of PTOX in microorganisms is summarized, which sets the foundation for the efficient microbial production of PTOX using cell factories.
Genes, Plant ; Podophyllotoxin/biosynthesis* ; Podophyllum/genetics*

Screening suitable reference genes is the premise of quantitative Real-time PCR(qRT-PCR)for gene expression analysis. To provide stable reference genes for expression analysis of genes in Aconitum vilmorinianum, this study selected 19 candidate re-ference genes(ACT1, ACT2, ACT3, aTUB1, aTUB2, bTUB, 18S rRNA, UBQ, eIF2, eIF3, eIF4, eIF5, CYP, GAPDH1, GAPDH2, PP2A1, PP2A2, ACP, and EF1α) based on the transcriptome data of A. vilmorinianum. qRT-PCR was conducted to profile the expression of these genes in the root, stem, leaf, and flower of A. vilmorinianum. The Ct values showed that 18S rRNA with high expression level and GAPDH2 with large expression difference among organs were not suitable as the reference genes. NormFinder and geNorm showed similar results of the expression stability of the other candidate reference genes and demonstrated PP2A1, EF1α, and CYP as the highly stable ones. However, BestKeeper suggested EF1α, ACT3, and PP2A1 as the top stable genes. In view of the different results from different softwares, the geometric mean method was employed to analyze the expression stability of the candidate re-ference genes, the results of which indicated that PP2A1, EF1α, and ACT3 were the most stable. Based on the comprehensive analysis results of geNorm, NormFinder, BestKeeper, and geometric mean method, PP2A1 and EF1α presented the most stable expression in different organs of A. vilmorinianum. PP2A1 and EF1α were the superior reference genes for gene expression profiling in different organs of A. vilmorinianum.
Aconitum ; Gene Expression Profiling ; Genes, Plant/genetics* ; Real-Time Polymerase Chain Reaction ; Reference Standards ; Reverse Transcriptase Polymerase Chain Reaction