Genome-wide analysis of Family-1 UDP-glycosyltransferases in soybean confirms their abundance and varied expression during seed development

Identification and analysis of UGT genes associated with triterpenoid saponin in soapberry (Sapindus mukorossi Gaertn.)

BACKGROUND

Soapberry (Sapindus mukorossi) is an economically important multifunctional tree species. Triterpenoid saponins have many functions in soapberry. However, the types of uridine diphosphate (UDP) glucosyltransferases (UGTs) involved in the synthesis of triterpenoid saponins in soapberry have not been clarified.

RESULTS

In this study, 42 SmUGTs were identified in soapberry, which were unevenly distributed on 12 chromosomes and had sequence lengths of 450 bp to 1638 bp, with an average of 1388 bp. The number of amino acids in SmUGTs was 149 to 545, with an average of 462. Most SmUGTs were acidic and hydrophilic unstable proteins, and their secondary structures were mainly α-helices and random coils. All had conserved UDPGT and PSPG-box domains. Phylogenetic analysis divided them into four subclasses, which glycosylated different carbon atoms. Prediction of cis-acting elements suggested roles of SmUGTs in plant development and responses to environmental stresses. The expression patterns of SmUGTs differed according to the developmental stage of fruits, as determined by transcriptomics and RT-qPCR. Co-expression network analysis of SmUGTs and related genes/transcription factors in the triterpenoid saponin synthesis pathway was also performed. The results indicated potential roles for many transcription factors, such as SmERFs, SmGATAs and SmMYBs. A correlation analysis showed that 42 SmUGTs were crucial in saponin synthesis in soapberry.

CONCLUSIONS

Our findings suggest optimal targets for manipulating glycosylation in soapberry triterpenoid saponin biosynthesis; they also provide a theoretical foundation for further evaluation of the functions of SmUGTs and analyses of their biosynthetic mechanisms.

Identification and Functional Verification of the Glycosyltransferase Gene Family Involved in Flavonoid Synthesis in Rubus chingii Hu

Glycosylation is catalyzed by UDP-glycosyltransferase (UGT) and plays an important role in enriching the diversity of flavonoids. Rubus plants contain a lot of natural flavonoid glycosides, which are important plants with a homology of medicine and food. However, information about the Rubus UGT gene family is very limited. In this study, we carried out genome-wide analysis and identified the 172, 121, 130, 121 UGT genes in R. chingii, R. corchorifolius, R. idaeus, and R. occidentalis, respectively, and divided them into 18 groups. The analysis of the protein motif and gene structure showed that there were structural and functional conservations in the same group, but there were differences among different groups. Gene replication analysis showed that raspberry and dicotyledons had a higher homology. The expansion of the UGTs gene family was mainly driven by tandem replication events, and experienced purified selection during the long evolution of the raspberry. Cis-acting element analysis showed that they were related to plant growth and development, hormone regulation, and stress response. In addition, according to a comprehensive analysis of the co-expression network constructed by transcriptome data and phylogenetic homology, RchUGT169 was identified as a flavonoid glucosyltransferase. Through the transient expression in tobacco, it was verified that RchUGT169 could catalyze the conversion of kaempferol and quercetin to the corresponding flavonoid glycosides. In conclusion, this research enriched the understanding of the diversity of UGTs in Rubus and determined that RcUGT169 can catalyze flavonoids.

Genome-wide analysis of UDP-glycosyltransferases family and identification of UGT genes involved in drought stress of Platycodon grandiflorus

Introduction

The uridine diphosphate (UDP)-glycosyltransferase (UGT) family is the largest glycosyltransferase family, which is involved in the biosynthesis of natural plant products and response to abiotic stress. UGT has been studied in many medicinal plants, but there are few reports on Platycodon grandiflorus. This study is devoted to genome-wide analysis of UGT family and identification of UGT genes involved in drought stress of Platycodon grandiflorus (PgUGTs).

Methods

The genome data of Platycodon grandiflorus was used for genome-wide identification of PgUGTs, online website and bioinformatics analysis software was used to conduct bioinformatics analysis of PgUGT genes and the genes highly responsive to drought stress were screened out by qRT-PCR, these genes were cloned and conducted bioinformatics analysis.

Results

A total of 75 PgUGT genes were identified in P.grandiflorus genome and clustered into 14 subgroups. The PgUGTs were distributed on nine chromosomes, containing multiple cis-acting elements and 22 pairs of duplicate genes were identified. Protein-protein interaction analysis was performed to predict the interaction between PgUGT proteins. Additionally, six genes were upregulated after 3d under drought stress and three genes (PGrchr09G0563, PGrchr06G0523, PGrchr06G1266) responded significantly to drought stress, as confirmed by qRT-PCR. This was especially true for PGrchr06G1266, the expression of which increased 16.21-fold after 3d of treatment. We cloned and conducted bioinformatics analysis of three candidate genes, both of which contained conserved motifs and several cis-acting elements related to stress response, PGrchr06G1266 contained the most elements.

Discussion

PgGT1 was confirmed to catalyze the C-3 position of platycodin D and only eight amino acids showed differences between gene PGr008G1527 and PgGT1, which means PGr008G1527 may be able to catalyze the C-3 position of platycodin D in the same manner as PgGT1. Seven genes were highly expressed in the roots, stems, and leaves, these genes may play important roles in the development of the roots, stems, and leaves of P. grandiflorus. Three genes were highly responsive to drought stress, among which the expression of PGrchr06G1266 was increased 16.21-fold after 3d of drought stress treatment, indicating that PGrchr06G1266 plays an important role in drought stress tolerance. To summarize, this study laied the foundation to better understand the molecular bases of responses to drought stress and the biosynthesis of platycodin.

Genome-wide investigation of UDP-Glycosyltransferase family in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum) belongs to Polygonaceae family and has attracted increasing attention owing to its high nutritional value. UDP-glycosyltransferases (UGTs) glycosylate a variety of plant secondary metabolites to control many metabolic processes during plant growth and development. However, there have been no systematic reports of UGT superfamily in F. tataricum.

RESULTS

We identified 173 FtUGTs in F. tataricum based on their conserved UDPGT domain. Phylogenetic analysis of FtUGTs with 73 Arabidopsis UGTs clustered them into 21 families. FtUGTs from the same family usually had similar gene structure and motif compositions. Most of FtUGTs did not contain introns or had only one intron. Tandem repeats contributed more to FtUGTs amplification than segmental duplications. Expression analysis indicates that FtUGTs are widely expressed in various tissues and likely play important roles in plant growth and development. The gene expression analysis response to different abiotic stresses showed that some FtUGTs were involved in response to drought and cadmium stress. Our study provides useful information on the UGTs in F. tataricum, and will facilitate their further study to better understand their function.

CONCLUSIONS

Our results provide a theoretical basis for further exploration of the functional characteristics of FtUGTs and for understanding the growth, development, and metabolic model in F. tataricum.

Genome-wide analysis of UDP-glycosyltransferase gene family and identification of a flavonoid 7-O-UGT (AhUGT75A) enhancing abiotic stress in peanut (Arachis hypogaea L.)

BACKGROUND

Glycosylation, catalyzed by UDP-glycosyltransferase (UGT), was important for enhancing solubility, bioactivity, and diversity of flavonoids. Peanut (Arachis hypogaea L.) is an important oilseed and cash crop worldwide. In addition to provide high quality of edible oils and proteins, peanut seeds contain a rich source of flavonoid glycosides that benefit human health. However, information of UGT gene family was quite limited in peanut.

RESULTS

In present study, a total of 267 AhUGTs clustered into 15 phylogenetic groups were identified in peanut genome. Group I has greatly expanded to contain the largest number of AhUGT genes. Segmental duplication was the major driving force for AhUGT gene family expansion. Transcriptomic analysis of gene expression profiles in various tissues and under different abiotic stress treatments indicated AhUGTs were involved in peanut growth and abiotic stress response. AhUGT75A (UGT73CG33), located in mitochondria, was characterized as a flavonoid 7-O-UGT by in vitro enzyme assays. The transcript level of AhUGT75A was strongly induced by abiotic stress. Overexpression of AhUGT75A resulted in accumulating less amount of malondialdehyde (MDA) and superoxide, and enhancing tolerance against drought and/or salt stress in transgenic Arabidopsis. These results indicated AhUGT75A played important roles in conferring abiotic stress tolerance through reactive oxygen species scavenging.

CONCLUSIONS

Our research only not provides valuable information for functional characterization of UGTs in peanut, but also gives new insights into potential applications in breeding new cultivars with both desirable stress tolerance and health benefits.

Protein Modelling Highlighted Key Catalytic Sites Involved in Position-Specific Glycosylation of Isoflavonoids

Uridine diphosphate glycosyltransferases (UGTs) are known for promiscuity towards sugar acceptors, a valuable characteristic for host plants but not desirable for heterologous biosynthesis. UGTs characterized for the O-glycosylation of isoflavonoids have shown a variable efficiency, substrate preference, and OH site specificity. Thus, 22 UGTs with reported isoflavonoid O-glycosylation activity were analyzed and ranked for OH site specificity and catalysis efficiency. Multiple-sequence alignment (MSA) showed a 33.2% pairwise identity and 4.5% identical sites among selected UGTs. MSA and phylogenetic analysis highlighted a comparatively higher amino acid substitution rate in the N-terminal domain that likely led to a higher specificity for isoflavonoids. Based on the docking score, OH site specificity, and physical and chemical features of active sites, selected UGTs were divided into three groups. A significantly high pairwise identity (67.4%) and identical sites (31.7%) were seen for group 1 UGTs. The structural and chemical composition of active sites highlighted key amino acids that likely define substrate preference, OH site specificity, and glycosylation efficiency towards selected (iso)flavonoids. In conclusion, physical and chemical parameters of active sites likely control the position-specific glycosylation of isoflavonoids. The present study will help the heterologous biosynthesis of glycosylated isoflavonoids and protein engineering efforts to improve the substrate and site specificity of UGTs.

Genome-wide analysis of UDP-glycosyltransferases family and identification of UGT genes involved in abiotic stress and flavonol biosynthesis in Nicotiana tabacum

BACKGROUND

Uridine disphosphate (UDP) glycosyltransferases (UGTs) act upon a huge variety of highly diverse and complex substrates, such as phytohormones and specialized metabolites, to regulate plant growth, development, disease resistance, and environmental interactions. However, a comprehensive investigation of UGT genes in tobacco has not been conducted.

RESULTS

In this study, we carried out a genome-wide analysis of family-1 UDP glycosyltransferases in Nicotiana tabacum. We predicted 276 NtUGT genes, which were classified into 18 major phylogenetic subgroups. The NtUGT genes were invariably distributed among all the 24 chromosomes with structural diversity in exon/intron structure, conserved motifs, and cis-acting elements of promoters. Three groups of proteins which involved in flavonoid biosynthesis, plant growth and development, transportation and modification were identified that interact with NtUGT proteins using the PPI analysis. Expression analysis of NtUGT genes in cold stress, drought stress and different flower color using both online RNA-Seq data and the realtime PCR analysis, suggested the distinct role of NtUGT genes in resistance of cold, drought and in flavonoid biosynthesis. The enzymatic activities of seven NtUGT proteins that potentially involved in flavonoid glycosylation were analyzed, and found that all seven exhibited activity on myricetin; six (NtUGT108, NtUGT123, NtUGT141, NtUGT155, NtUGT179, and NtUGT195) showed activity on cyanidin; and three (NtUGT108, NtUGT195, and NtUGT217) were active on the flavonol aglycones kaempferol and quercetin, which catalyzing the substrates (myricetin, cyanidin or flavonol) to form new products. We further investigated the enzymatic products and enzymatic properties of NtUGT108, NtUGT195, and NtUGT217, suggested their diverse enzymatic activity toward flavonol, and NtUGT217 showed the highest catalyzed efficient toward quercetin. Overexpression of NtUGT217 significantly increase the content levels of the quercetin-3-O-glucoside, quercetin-3-O-rutinoside and kaempferol-3-O-rutinoside in transgenic tobacco leaves.

CONCLUSION

We identified 276 UGT genes in Nicotiana tabacum. Our study uncovered valuable information about the phylogenetic structure, distribution, genomic characters, expression patterns and enzymatic activity of NtUGT genes in tobacco. We further identified three NtUGT genes involved in flavonoid biosynthesis, and overexpressed NtUGT217 to validate its function in catalyze quercetin. The results provide key candidate NtUGT genes for future breeding of cold and drought resistance and for potential metabolic engineering of flavonoid compounds.

Genome-Wide Association Study Revealed SNP Alleles Associated with Seed Size Traits in African Yam Bean (Sphenostylis stenocarpa (Hochst ex. A. Rich.) Harms)

Seed size is an important yield and quality-determining trait in higher plants and is also crucial to their evolutionary fitness. In African yam bean (AYB), seed size varies widely among different accessions. However, the genetic basis of such variation has not been adequately documented. A genome-wide marker-trait association study was conducted to identify genomic regions associated with four seed size traits (seed length, seed width, seed thickness, and 100-seed weight) in a panel of 195 AYB accessions. A total of 5416 SNP markers were generated from the diversity array technology sequence (DArTseq) genotype-by-sequencing (GBS)- approach, in which 2491 SNPs were retained after SNP quality control and used for marker-trait association analysis. Significant phenotypic variation was observed for the traits. Broad-sense heritability ranged from 50.0% (seed width) to 66.4% (seed length). The relationships among the traits were positive and significant. Genome-wide association study (GWAS) using the general linear model (GLM) and the mixed linear model (MLM) approaches identified 12 SNP markers significantly associated with seed size traits across the six test environments. The 12 makers explained 6.5-10.8% of the phenotypic variation. Two markers (29420334|F|0-52:C>G-52:C>G and 29420736|F|0-57:G>T-57:G>T) with pleiotropic effects associated with seed width and seed thickness were found. A candidate gene search identified five significant markers (100026424|F|0-37:C>T-37:C>T, 100041049|F|0-42:G>C-42:G>C, 100034480|F|0-31:C>A-31:C>A, 29420365|F|0-55:C>G-55:C>G, and 29420736|F|0-57:G>T-57:G>T) located close to 43 putative genes whose encoding protein products are known to regulate seed size traits. This study revealed significant makers not previously reported for seed size in AYB and could provide useful information for genomic-assisted breeding in AYB.

Identification and bioinformatic analysis of the CaCesA/Csls family members and the expression of the CaCslD1 in the flower buds of CMS/Rf system in pepper

The cellulose synthase gene superfamily contains cellulose synthase (CesA) and cellulose synthase-like (Csl) gene families, which synthesize cellulose and hemicellulose in plant cell walls and play a crucial role in plant growth and development. However, the CesA/Csl gene family has not been reported in pepper. Therefore, the genome-wide research of the CaCesA/CaCsl gene family was conducted in pepper. In this study, a total of 39 CaCesA/CaCsls genes (10 CesAs genes and 29 Csls genes) were identified in pepper and unevenly distributed on 11 chromosomes. These CaCesA/Csls were divided into seven subfamilies (CesAs, CslAs, CslBs, CslCs, CslDs, CslEs, CslGs), and most of CaCesA/Csls genes are closely related to AtCesA/Csls genes. The cis-acting elements in the promoters of CaCesA/Csls genes are mainly related to hormone response and stress response. There are ten collinear gene pairs between the CesA/Csls gene family of pepper and Arabidopsis, and four fragment duplication gene pairs of the CaCesA/Csls genes were discovered. RNA-seq analysis shows that the majority of CaCesA/Csls are expressed in a variety of plant tissues, indicating that most CaCesA/Csls gene expression patterns are not organ-specific, and CaCslD1/D4 have the highest expression in anthers, followed by petal, ovary, and F9. RNA-seq analysis shows that most CaCesA/Csls are responsive to five hormones (IAA, GA₃, ABA, SA, and MeJA). The tissue-specific expression analysis of the CaCslD1 gene shows that the CaCslD1 gene is expressed specifically in flowers. In the flower buds IV of cytoplasmic male sterility (CMS) and its restoration of fertility (Rf) system, CaCslD1 reach the highest expression respectively. However, the relative expression level of CaCslD1 in the fertile accessions is extremely significantly higher than in the sterile accessions. This study shows an overall understanding of the CaCesA/Csls gene family and provides a new insight for understanding the function of CaCslD1 in pollen development and exploring the fertility restoration of CMS in pepper.

Genome-Wide Analysis of UGT Genes in Petunia and Identification of PhUGT51 Involved in the Regulation of Salt Resistance

UDP-glycosyltransferase (UGT) plays an essential role in regulating the synthesis of hormones and secondary metabolites in plants. In this study, 129 members of the Petunia UGT family were identified and classified into 16 groups (A-P) based on phylogenetic analysis. The same subgroups have conserved motif compositions and intron/exon arrangement. In the promoters of the Petunia UGT genes, several cis-elements associated with plant hormones, growth and development, and abiotic stress have been discovered. Their expression profiles in five tissues were revealed by tissue expression based on RNA-seq data. Subcellular localization analysis showed that PhUGT51 was located in the nucleus and cell membrane. Salt stress caused an increase in the expression level of PhUGT51, but the expression level remained stable with the growth over time. In addition, the overexpression of PhUGT51 caused a significant increase in salt resistance. Our study systematically analyses the UGT gene family in Petunia for the first time and provides some valuable clues for the further functional studies of UGT genes.

Genome Wide Analysis of Family-1 UDP Glycosyltransferases in Populus trichocarpa Specifies Abiotic Stress Responsive Glycosylation Mechanisms

Populus trichocarpa (Black cottonwood) is a dominant timber-yielding tree that has become a notable model plant for genome-level insights in forest trees. The efficient transport and solubility of various glycoside-associated compounds is linked to Family-1 UDP-glycosyltransferase (EC 2.4.1.x; UGTs) enzymes. These glycosyltransferase enzymes play a vital role in diverse plant functions, such as regulation of hormonal homeostasis, growth and development (seed, flower, fiber, root, etc.), xenobiotic detoxification, stress response (salt, drought, and oxidative), and biosynthesis of secondary metabolites. Here, we report a genome-wide analysis of the P. trichocarpa genome that identified 191 putative UGTs distributed across all chromosomes (with the exception of chromosome 20) based on 44 conserved plant secondary product glycosyltransferase (PSPG) motif amino acid sequences. Phylogenetic analysis of the 191 Populus UGTs together with 22 referenced UGTs from Arabidopsis and maize clustered the putative UGTs into 16 major groups (A–P). Whole-genome duplication events were the dominant pattern of duplication among UGTs in Populus. A well-conserved intron insertion was detected in most intron-containing UGTs across eight examined eudicots, including Populus. Most of the UGT genes were found preferentially expressed in leaf and root tissues in general. The regulation of putative UGT expression in response to drought, salt and heat stress was observed based on microarray and available RNA sequencing datasets. Up- and down-regulated UGT expression models were designed, based on transcripts per kilobase million values, confirmed their maximally varied expression under drought, salt and heat stresses. Co-expression networking of putative UGTs indicated their maximum co-expression with cytochrome P450 genes involved in triterpenoid biosynthesis. Our results provide an important resource for the identification of functional UGT genes to manipulate abiotic stress responsive glycosylation in Populus.

Glutathione Promotes Degradation and Metabolism of Residual Fungicides by Inducing UDP-Glycosyltransferase Genes in Tomato

Reduced glutathione (GSH) is a key antioxidant, which plays a crucial role in the detoxification of xenobiotics in plants. In the present study, glutathione could reduce chlorothalonil (CHT) residues in tomatoes by inducing the expression of the UDP-glycosyltransferase (UGT) gene. In plants, UGT is an important glycosylation catalyst, which can respond to stresses in time by activating plant hormones and defense compounds. Given the importance of plant growth and development, the genome-wipe analyses of Arabidopsis and soybean samples have been carried out, though not on the tomato, which is a vital vegetable crop. In this study, we identified 143 UGT genes in the tomato that were unevenly distributed on 12 chromosomes and divided into 16 subgroups and found that a variety of plant hormones and stress response cis-elements were discovered in the promoter region of the SlUGT genes, indicating that the UGT genes were involved in several aspects of the tomato stress response. Transcriptome analysis and results of qRT-PCR showed that most SlUGT genes could be induced by CHT, and the expression of these genes was regulated by glutathione. In addition, we found that SlUGT genes could participate in plant detoxification through interaction with transcription factors. These findings further clarify the potential function of the UGT gene family in the detoxification of exogenous substances in tomatoes and provide valuable information for the future study of functional genomics of tomatoes.

Genome-Wide Analysis and Profile of UDP-Glycosyltransferases Family in Alfalfa (Medicago sativa L.) under Drought Stress

Drought stress is one of the major constraints that decreases global crop productivity. Alfalfa, planted mainly in arid and semi-arid areas, is of crucial importance in sustaining the agricultural system. The family 1 UDP-glycosyltransferases (UGT) is indispensable because it takes part in the regulation of plant growth and stress resistance. However, a comprehensive insight into the participation of the UGT family in adaptation of alfalfa to drought environments is lacking. In the present study, a genome-wide analysis and profiling of the UGT in alfalfa were carried out. A total of 409 UGT genes in alfalfa (MsUGT) were identified and they are clustered into 13 groups. The expression pattern of MsUGT genes were analyzed by RNA-seq data in six tissues and under different stresses. The quantitative real-time PCR verification genes suggested the distinct role of the MsUGT genes under different drought stresses and abscisic acid (ABA) treatment. Furthermore, the function of MsUGT003 and MsUGT024, which were upregulated under drought stress and ABA treatment, were characterized by heterologous expression in yeast. Taken together, this study comprehensively analyzed the UGT gene family in alfalfa for the first time and provided useful information for improving drought tolerance and in molecular breeding of alfalfa.

Glycosyltransferases: Mining, engineering and applications in biosynthesis of glycosylated plant natural products

UDP-Glycosyltransferases (UGTs) catalyze the transfer of nucleotide-activated sugars to specific acceptors, among which the GT1 family enzymes are well-known for their function in biosynthesis of natural product glycosides. Elucidating GT function represents necessary step in metabolic engineering of aglycone glycosylation to produce drug leads, cosmetics, nutrients and sweeteners. In this review, we systematically summarize the phylogenetic distribution and catalytic diversity of plant GTs. We also discuss recent progress in the identification of novel GT candidates for synthesis of plant natural products (PNPs) using multi-omics technology and deep learning predicted models. We also highlight recent advances in rational design and directed evolution engineering strategies for new or improved GT functions. Finally, we cover recent breakthroughs in the application of GTs for microbial biosynthesis of some representative glycosylated PNPs, including flavonoid glycosides (fisetin 3-O-glycosides, astragalin, scutellarein 7-O-glucoside), terpenoid glycosides (rebaudioside A, ginsenosides) and polyketide glycosides (salidroside, polydatin).

Functional analysis of a novel C-glycosyltransferase in the orchid Dendrobium catenatum

Flavonoids, which are a diverse class of phytonutrients, are used by organisms to respond to nearly all abiotic stresses and are beneficial for human health. Glycosyltransferase, used during the last step of flavonoid biosynthesis, is important in flavonoid enrichment. However, little is known about glycosyltransferase in the orchid Dendrobium catenatum (D. officinale). In this study, we isolated a novel C-glycosyltransferase (designated DcaCGT) from the orchid D. catenatum by identifying and analyzing 82 putative genes in the GT1 family. DcaCGT could specifically catalyze not only di-C-glycosylation but also O-glycosylation. Apart from the normal function of catalyzing 2-hydroxynaringenin and phloretin to the respective di-C-glycosides, DcaCGT also catalyzes apigenin to cosmosiin. Targeted metabolic profiling of the substrates (2-hydroxynaringenin, phloretin, and apigenin) and products (vitexin, isovitexin, vicenin-2, nothofagin, 3',5'-di-C-glucosylphloretin, and cosmosiin) in different tissues showed that vicenin-2 was the most abundant product of this novel enzyme. Cosmosiin was detected in flowers and flower buds. We also established that DcaCGT functions expanded throughout the evolution of D. catenatum. Residual OGT activity may help D. catenatum resist drought stress. Our study illustrates the function, origin, and differentiation of DcaCGT and provides insights into glycosylation and molecular propagation processes, which can be used to improve the production of flavonoids by the cultivated medicinal plant D. catenatum.

Reconstruction of the Evolutionary Histories of UGT Gene Superfamily in Legumes Clarifies the Functional Divergence of Duplicates in Specialized Metabolism

Plant uridine 5'-diphosphate glycosyltransferases (UGTs) influence the physiochemical properties of several classes of specialized metabolites including triterpenoids via glycosylation. To uncover the evolutionary past of UGTs of soyasaponins (a group of beneficial triterpene glycosides widespread among Leguminosae), the UGT gene superfamily in Medicago truncatula, Glycine max, Phaseolus vulgaris, Lotus japonicus, and Trifolium pratense genomes were systematically mined. A total of 834 nonredundant UGTs were identified and categorized into 98 putative orthologous loci (POLs) using tree-based and graph-based methods. Major key findings in this study were of, (i) 17 POLs represent potential catalysts for triterpene glycosylation in legumes, (ii) UGTs responsible for the addition of second (UGT73P2: galactosyltransferase and UGT73P10: arabinosyltransferase) and third (UGT91H4: rhamnosyltransferase and UGT91H9: glucosyltransferase) sugars of the C-3 sugar chain of soyasaponins were resulted from duplication events occurred before and after the hologalegina-millettoid split, respectively, and followed neofunctionalization in species-/ lineage-specific manner, and (iii) UGTs responsible for the C-22-O glycosylation of group A (arabinosyltransferase) and DDMP saponins (DDMPtransferase) and the second sugar of C-22 sugar chain of group A saponins (UGT73F2: glucosyltransferase) may all share a common ancestor. Our findings showed a way to trace the evolutionary history of UGTs involved in specialized metabolism.

Genome-Wide Identification and Characterization of Glycosyltransferase Family 47 in Cotton

The glycosyltransferase (GT) 47 family is involved in the biosynthesis of xylose, pectin and xyloglucan and plays a significant role in maintaining the normal morphology of the plant cell wall. However, the functions of GT47s are less well known in cotton. In the present study, a total of 53, 53, 105 and 109 GT47 genes were detected by genome-wide identification in Gossypium arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. All the GT47s were classified into six major groups via phylogenetic analysis. The exon/intron structure and protein motifs indicated that each branch of the GT47 genes was highly conserved. Collinearity analysis showed that GT47 gene family expansion occurred in Gossypium spp. mainly through whole-genome duplication and that segmental duplication mainly promoted GT47 gene expansion within the A and D subgenomes. The Ka/Ks values suggested that the GT47 gene family has undergone purifying selection during the long-term evolutionary process. Transcriptomic data and qRT-PCR showed that GhGT47 genes exhibited different expression patterns in each tissue and during fiber development. Our results suggest that some genes in the GhGT47 family might be associated with fiber development and the abiotic stress response, which could promote further research involving functional analysis of GT47 genes in cotton.

Utilization of Interspecific High-Density Genetic Map of RIL Population for the QTL Detection and Candidate Gene Mining for 100-Seed Weight in Soybean

Seed-weight is one of the most important traits determining soybean yield. Hence, it is prerequisite to have detailed understanding of the genetic basis regulating seed-weight for the development of improved cultivars. In this regard, the present study used high-density interspecific linkage map of NJIR4P recombinant inbred population evaluated in four different environments to detect stable Quantitative trait loci (QTLs) as well as mine candidate genes for 100-seed weight. In total, 19 QTLs distributed on 12 chromosomes were identified in all individual environments plus combined environment, out of which seven were novel and eight are stable identified in more than one environment. However, all the novel QTLs were minor (R 2 < 10%). The remaining 12 QTLs detected in this study were co-localized with the earlier reported QTLs with narrow genomic regions, and out of these only 2 QTLs were major (R 2 > 10%) viz., qSW-17-1 and qSW-17-4. Beneficial alleles of all identified QTLs were derived from cultivated soybean parent (Nannong493-1). Based on Protein ANalysis THrough Evolutionary Relationships, gene annotation information, and literature search, 29 genes within 5 stable QTLs were predicted to be possible candidate genes that might regulate seed-weight/size in soybean. However, it needs further validation to confirm their role in seed development. In conclusion, the present study provides better understanding of trait genetics and candidate gene information through the use high-density inter-specific bin map, and also revealed considerable scope for genetic improvement of 100-seed weight in soybean using marker-assisted breeding.

Family-1 UDP glycosyltransferases in pear (Pyrus bretschneideri): Molecular identification, phylogenomic characterization and expression profiling during stone cell formation

Stone cells are a characteristic trait of pear fruits, and excessive stone cell formation has a significant negative impact on the texture and flavour of the pulp. Lignin is one of the main components of stone cells. Family-1 uridine diphosphate-glycosyltransferases (UGTs) are responsible for the glycosylation modification of monolignols. However, information remains limited regarding the relationship between UGTs and stone cell formation. To address this problem, we identified 139 UGTs from the pear genome, which were distributed in 15 phylogenetic groups (A-M, O, and P). We also performed a collinearity analysis of UGTs among four Rosaceae plants (pear, peach, mei, and strawberry). Phylogenetic analysis suggested that 13 PbUGTs might be related to the glycosylation of monolignols. Analysis of expression patterns demonstrated that most putative monolignol glycosylation-related PbUGTs not only showed high expression levels in flowers and buds but were also induced by exogenous ABA, SA, and MeJA. In addition, the transcript level of Pbr005014.1 (named PbUGT72AJ2) was consistent with the changing trend of lignin content in pear fruit, and the transcript level was also higher in 'Dangshan Su' pear with higher lignin and stone cell contents. Subcellular localization results showed that PbUGT72AJ2 was located mainly in the cytomembrane and cytoplasm. Based on our study, PbUGT72AJ2 is considered to be a monolignol glycosylation-related UGT. Our results provide an important source for the identification of UGTs and a foundation for the future understanding and manipulation of lignin metabolism and stone cell formation in pear fruit.

Comparative Genomic and Transcriptomic Analyses of Family-1 UDP Glycosyltransferase in Prunus Mume

Glycosylation mediated by Family-1 UDP-glycosyltransferases (UGTs) plays crucial roles in plant growth and adaptation to various stress conditions. Prunus mume is an ideal crop for analyzing flowering for its early spring flowering characteristics. Revealing the genomic and transcriptomic portfolio of the UGT family in P. mume, a species in which UGTs have not yet been investigated, is therefore important. In this study, 130 putative UGT genes were identified and phylogenetically clustered into 14 groups. These PmUGTs were distributed unevenly across eight chromosomes and 32 tandem duplication and 8 segmental duplication pairs were revealed. A highly conserved intron insertion event was revealed on the basis of intron/exon patterns within PmUGTs. According to RNA-seq data, these PmUGTs were specifically expressed in different tissues and during the bud dormancy process. In addition, we confirmed the differential expression of some representative genes in response to abscisic acid treatment. Our results will provide important information on the UGT family in P. mume that should aid further characterization of their biological roles in response to environmental stress.

Characterization of Cellulose Synthase A (CESA) Gene Family in Eudicots

Cellulose synthase A (CESA) is a key enzyme involved in the complex process of plant cell wall biosynthesis, and it remains a productive subject for research. We employed systems biology approaches to explore structural diversity of eudicot CESAs by exon-intron organization, mode of duplication, synteny, and splice site analyses. Using a combined phylogenetics and comparative genomics approach coupled with co-expression networks we reconciled the evolution of cellulose synthase gene family in eudicots and found that the basic forms of CESA proteins are retained in angiosperms. Duplications have played an important role in expansion of CESA gene family members in eudicots. Co-expression networks showed that primary and secondary cell wall modules are duplicated in eudicots. We also identified 230 simple sequence repeat markers in 103 eudicot CESAs. The 13 identified conserved motifs in eudicots will provide a basis for gene identification and functional characterization in other plants. Furthermore, we characterized (in silico) eudicot CESAs against senescence and found that expression levels of CESAs decreased during leaf senescence.

Evolution of intron-poor clades and expression patterns of the glycosyltransferase family 47

A large-scale bioinformatics analysis revealed the origin and evolution of GT47 gene family, and identified two clades of intron-poor genes with putative functions in drought stress responses and seed development in maize. Glycosyltransferase family 47 (GT47) genes encode β-galactosyltransferases and β-glucuronyltransferases that synthesize pectin, xyloglucans and xylan, which are important components of the plant cell wall. In this study, we performed a systematic and large-scale bioinformatics analysis of GT47 gene family using 352 GT47 proteins from 15 species ranging from cyanobacteria to seed plants. The analysis results showed that GT47 family may originate in cyanobacteria and expand along the evolutionary trajectory to moss. Further analysis of 47 GT47 genes in maize revealed that they can divide into five clades with diverse exon-intron structures. Among these five clades, two were mainly composed with intron-poor genes, which may originate in the moss. Gene duplication analysis revealed that the expansion of GT47 gene family in maize was significantly driven from tandem duplication events and segmental duplication events. Significantly, almost all duplicated genes are intron-poor genes. Expression analysis indicated that several intron-poor GT47 genes may be involved in the drought stress response and seed development in maize. This work provides insight into the origin and evolutionary process, expansion mechanisms and expression patterns of GT47 genes, thus facilitating their functional investigations in the future.

Functional characterization of naturally occurring wild soybean mutant (sg-5) lacking astringent saponins using whole genome sequencing approach

Triterpenoid saponins are one of the most highly accumulated groups of functional components in soybean (Glycine max) and the oxidative reactions during their biosynthesis are required for their aglycone diversity. Natural mutants of soyasaponins in wild soybean (Glycine soja) are valuable resources for establishing the soyasaponin biosynthesis pathway and breeding new soybean varieties. In this study, we investigated the genetic mechanism behind the absence of group A saponins in a Korean wild soybean mutant, CWS5095. Whole genome sequencing (WGS) of CWS5095 identified four point mutations [Val6 → Asp, Ile231 → Thr, His294 → Gln, and Arg376 → Lys] in CYP72A69 (Glyma15g39090), which oxygenate the C-21 position of soyasapogenol B or other intermediates to produce soyasapogenol A, leading to group A saponin production. An in vitro enzyme activity assay of single-sited mutated clones indicated that the Arg376 > Lys mutation (a highly conserved mutation based on a nucleotide change from G → A at the 1,127th position) may lead to loss of gene function in the sg-5 mutant. A very high normalized expression value of 377 reads per kilo base per million (RPKM) of Glyma15g39090 in the hypocotyl axis at the early maturation seed-development stage confirmed their abundant presence in seed hypocotyls. A molecular dynamics analysis of the Arg376 > Lys mutation based on the CYP3A4 (a human CYP450) protein structure found that it was responsible for the increase in axis length toward the heme (active site), which is critically important for biological activity and ligand binding. Our results provide important information on how to eradicate bitter and astringent saponins in soybean by utilizing the reported mutation in Glyma15g39090, and its importance for seed hypocotyl development based on transcript abundance.

Comparative genomic and transcriptomic analyses of Family-1 UDP glycosyltransferase in three Brassica species and Arabidopsis indicates stress-responsive regulation

In plants, UGTs (UDP-glycosyltransferases) glycosylate various phytohormones and metabolites in response to biotic and abiotic stresses. Little is known about stress-responsive glycosyltransferases in plants. Therefore, it is important to understand the genomic and transcriptomic portfolio of plants with regard to biotic and abiotic stresses. Here, we identified 140, 154, and 251 putative UGTs in Brassica rapa, Brassica oleracea, and Brassica napus, respectively, and clustered them into 14 major phylogenetic groups (A–N). Fourteen major KEGG pathways and 24 biological processes were associated with the UGTs, highlighting them as unique modulators against environmental stimuli. Putative UGTs from B. rapa and B. oleracea showed a negative selection pressure and biased gene fractionation pattern during their evolution. Polyploidization increased the intron proportion and number of UGT-containing introns among Brassica. The putative UGTs were preferentially expressed in developing tissues and at the senescence stage. Differential expression of up- and down-regulated UGTs in response to phytohormone treatments, pathogen responsiveness and abiotic stresses, inferred from microarray and RNA-Seq data in Arabidopsis and Brassica broaden the glycosylation impact at the molecular level. This study identifies unique candidate UGTs for the manipulation of biotic and abiotic stress pathways in Brassica and Arabidopsis.

De Novo Assembly and Analysis of Polygonatum sibiricum Transcriptome and Identification of Genes Involved in Polysaccharide Biosynthesis

Polygonatum sibiricum polysaccharides (PSPs) are used to improve immunity, alleviate dryness, promote the secretion of fluids, and quench thirst. However, the PSP biosynthetic pathway is largely unknown. Understanding the genetic background will help delineate that pathway at the molecular level so that researchers can develop better conservation strategies. After comparing the PSP contents among several different P. sibiricum germplasms, we selected two groups with the largest contrasts in contents and subjected them to HiSeq2500 transcriptome sequencing to identify the candidate genes involved in PSP biosynthesis. In all, 20 kinds of enzyme-encoding genes were related to PSP biosynthesis. The polysaccharide content was positively correlated with the expression patterns of β-fructofuranosidase (sacA), fructokinase (scrK), UDP-glucose 4-epimerase (GALE), Mannose-1-phosphate guanylyltransferase (GMPP), and UDP-glucose 6-dehydrogenase (UGDH), but negatively correlated with the expression of Hexokinase (HK). Through qRT-PCR validation and comprehensive analysis, we determined that sacA, HK, and GMPP are key genes for enzymes within the PSP metabolic pathway in P. sibiricum. Our results provide a public transcriptome dataset for this species and an outline of pathways for the production of polysaccharides in medicinal plants. They also present more information about the PSP biosynthesis pathway at the molecular level in P. sibiricum and lay the foundation for subsequent research of gene functions.

Genome and transcriptome-wide analyses of cellulose synthase gene superfamily in soybean

The plant cellulose synthase gene superfamily belongs to the category of type-2 glycosyltransferases, and is involved in cellulose and hemicellulose biosynthesis. These enzymes are vital for maintaining cell-wall structural integrity throughout plant life. Here, we identified 78 putative cellulose synthases (CS) in the soybean genome. Phylogenetic analysis against 40 reference Arabidopsis CS genes clustered soybean CSs into seven major groups (CESA, CSL A, B, C, D, E and G), located on 19 chromosomes (except chromosome 18). Soybean CS expansion occurred in 66 duplication events. Additionally, we identified 95 simple sequence repeat makers related to 44 CSs. We next performed digital expression analysis using publically available datasets to understand potential CS functions in soybean. We found that CSs were highly expressed during soybean seed development, a pattern confirmed with an Affymatrix soybean IVT array and validated with RNA-seq profiles. Within CS groups, CESAs had higher relative expression than CSLs. Soybean CS models were designed based on maximum average RPKM values. Gene co-expression networks were developed to explore which CSs could work together in soybean. Finally, RT-PCR analysis confirmed the expression of 15 selected CSs during all four seed developmental stages.

In-Depth Genomic and Transcriptomic Analysis of Five K+ Transporter Gene Families in Soybean Confirm Their Differential Expression for Nodulation

Plants have evolved a sophisticated network of K⁺ transport systems to regulate growth and development. Limited K⁺ resources are now forcing us to investigate how plant demand can be satisfied. To answer this complex question, we must understand the genomic and transcriptomic portfolio of K⁺ transporters in plants. Here, we have identified 70 putative K⁺ transporter genes from soybean, including 29 HAK/KT/KUP genes, 16 genes encoding voltage-gated K⁺ channels, 9 TPK/KCO genes, 4 HKT genes, and 12 KEA genes. To clarify the molecular evolution of each family in soybean, we analyzed their phylogeny, mode of duplication, exon structures and splice sites, and paralogs. Additionally, ortholog clustering and syntenic analysis across five other dicots further explored the evolution of these gene families and indicated that the soybean data is suitable as a model for all other legumes. Available microarray data sets from Genevestigator about nodulation was evaluated and further confirmed with the RNA sequencing data available by a web server. For each family, expression models were designed based on Transcripts Per Kilobase Million (TPM) values; the outcomes indicated differential expression linked to nodulation and confirmed the genes' putative roles. In-depth studies such as ours provides the basis for understanding K⁺ inventories in all other plants.