Genome-wide analysis of UDP-glycosyltransferase super family in Brassica rapa and Brassica oleracea reveals its evolutionary history and functional characterization

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.

Cloning, Expression, and Functional Characterization of Two Highly Efficient Flavonoid-di-O-glycosyltransferases ZmUGT84A1 and ZmUGT84A2 from Maize (Zea mays L.)

The maize (Zea mays L.) glycosyltransferase family 1 comprises many uridine diphosphate glycosyltransferase (UGT) members. However, UGT activities and biochemical functions have seldom been revealed. In this study, the genes of two flavonoid di-O-glycosyltransferases ZmUGT84A1 and ZmUGT84A2 were cloned from maize plant and expressed in Escherichia coli. Phylogenetic analysis showed that the two enzymes were homologous to AtUGT84A1 and AtUGT84A3. The two recombinant enzymes showed a high conversion rate of luteolin to its glucosides, mainly 4',7-di-O-glucoside and minorly 3',7-di-O-glucoside in two-step glycosylation reactions in vitro. Moreover, the recombinant ZmUGT84A1 and ZmUGT84A2 had a broad substrate spectrum, converting eriodictyol, naringenin, apigenin, quercetin, and kaempferol to monoglucosides and diglucosides. The highly efficient ZmUGT84A1 and ZmUGT84A2 may be used as a tool for the effective synthesis of various flavonoid O-glycosides and as markers for crop breeding to increase O-glycosyl flavonoid content in food.

Identification of Key Genes Associated with 1,2,6-Tri-O-galloyl-β-D-glucopyranose Accumulation in Camellia sinensis Based on Transcriptome Sequencing

Hydrolyzed tannin 1,2,6-tri-O-galloyl-β-D-glucopyranose (1,2,6-TGGP) possesses significant medicinal properties. However, little is known about its underlying molecular mechanisms. In this study, the levels of 1,2,6-TGGP in tea materials from different cultivars and leaf positions were compared. Additionally, one leaf and one bud sample from six tea cultivars with significant variations in 1,2,6-TGGP levels were analyzed using transcriptome high-throughput sequencing to identify the genes that are responsible for 1,2,6-TGGP accumulation. The sequencing results were mapped to the reference tea genome, revealing a total of 2735 differentially expressed genes (DEGs). This set included four UDP glycosyltransferase (UGTs) and six serine carboxypeptidases-like (SCPLs) genes. Among them, the upregulated SCPLs (CSS0032817) may directly participate in the acylation reaction of 1,2,6-TGGP. In addition, several classes of DEGs, including cytochrome P450, were significantly associated with the 1,2,6-TGGP content, which is potentially involved in their regulation. Overall, these results provide new insights into the molecular mechanism of 1,2,6-TGGP accumulation.

Systematic analysis of the UDP-glucosyltransferase family: discovery of a member involved in rutin biosynthesis in Solanum melongena

Eggplant (Solanum melongena) is an economically important crop and rich in various nutrients, among which rutin that has positive effects on human health is found in eggplant. Glycosylation mediated by UDP-glycosyltransferases (UGTs) is a key step in rutin biosynthesis. However, the UGT gene has not been reported in eggplant to date. Herein, 195 putative UGT genes were identified in eggplant by genome-wide analysis, and they were divided into 17 subgroups (Group A-P and Group R) according to the phylogenetic evolutionary tree. The members of Groups A, B, D, E and L were related to flavonol biosynthesis, and rutin was the typical flavonol. The expression profile showed that the transcriptional levels of SmUGT genes in Clusters 7-10 were closely related to those of rutin biosynthetic pathway genes. Notably, SmUGT89B2 was classified into Cluster 7 and Group B; its expression was consistent with rutin accumulation in different tissues and different leaf stages of eggplant. SmUGT89B2 was located in the nucleus and cell membrane. Virus-induced gene silencing (VIGS) and transient overexpression assays showed that SmUGT89B2 can promote rutin accumulation in eggplant. These findings provide new insights into the UGT genes in eggplant, indicating that SmUGT89B2 is likely to encode the final enzyme in rutin biosynthesis.

Gap-free genome assembly and CYP450 gene family analysis reveal the biosynthesis of anthocyanins in Scutellaria baicalensis

Scutellaria baicalensis Georgi, a member of the Lamiaceae family, is a widely utilized medicinal plant. The flavones extracted from S. baicalensis contribute to numerous health benefits, including anti-inflammatory, antiviral, and anti-tumor activities. However, the incomplete genome assembly hinders biological studies on S. baicalensis. This study presents the first telomere-to-telomere (T2T) gap-free genome assembly of S. baicalensis through the integration of Pacbio HiFi, Nanopore ultra-long and Hi-C technologies. A total of 384.59 Mb of genome size with a contig N50 of 42.44 Mb was obtained, and all sequences were anchored into nine pseudochromosomes without any gap or mismatch. In addition, we analysed the major cyanidin- and delphinidin-based anthocyanins involved in the determination of blue-purple flower using a widely-targeted metabolome approach. Based on the genome-wide identification of Cytochrome P450 (CYP450) gene family, three genes (SbFBH1, 2, and 5) encoding flavonoid 3'-hydroxylases (F3'Hs) and one gene (SbFBH7) encoding flavonoid 3'5'-hydroxylase (F3'5'H) were found to hydroxylate the B-ring of flavonoids. Our studies enrich the genomic information available for the Lamiaceae family and provide a toolkit for discovering CYP450 genes involved in the flavonoid decoration.

ATP binding cassette transporters and uridine diphosphate glycosyltransferases are ancient protein families that evolved roles in herbicide resistance through exaptation

ATP-binding cassette (ABC) transporters actively transport various substances across membranes, while uridine diphosphate (UDP) glycosyltransferases (UGTs) are proteins that catalyse the chemical modification of various organic compounds. Both of these protein superfamilies have been associated with conferring herbicide resistance in weeds. Little is known about the evolutionary history of these protein families in the Archaeplastida. To infer the evolutionary histories of these protein superfamilies, we compared protein sequences collected from 10 species which represent distinct lineages of the Archaeplastida-the lineage including glaucophyte algae, rhodophyte algae, chlorophyte algae and the streptophytes-and generated phylogenetic trees. We show that ABC transporters were present in the last common ancestor of the Archaeplastida which lived 1.6 billion years ago, and the major clades identified in extant plants were already present then. Conversely, we only identified UGTs in members of the streptophyte lineage, which suggests a loss of these proteins in earlier diverging Archaeplastida lineages or arrival of UGTs into a common ancestor of the streptophyte lineage through horizontal gene transfer from a non-Archaeplastida eukaryote lineage. We found that within the streptophyte lineage, most diversification of the UGT protein family occurred in the vascular lineage, with 17 of the 20 clades identified in extant plants present only in vascular plants. Based on our findings, we conclude that ABC transporters and UGTs are ancient protein families which diversified during Archaeplastida evolution, which may have evolved for developmental functions as plants began to occupy new environmental niches and are now being selected to confer resistance to a diverse range of herbicides in weeds.

GWAS with principal component analysis identify QTLs associated with main peanut flavor-related traits

Peanut flavor is a complex and important trait affected by raw material and processing technology owing to its significant impact on consumer preference. In this research, principal component analysis (PCA) on 33 representative traits associated with flavor revealed that total sugars, sucrose, and total tocopherols provided more information related to peanut flavor. Genome-wide association studies (GWAS) using 102 U.S. peanut mini-core accessions were performed to study associations between 12,526 single nucleotide polymorphic (SNP) markers and the three traits. A total of 7 and 22 significant quantitative trait loci (QTLs) were identified to be significantly associated with total sugars and sucrose, respectively. Among these QTLs, four and eight candidate genes for the two traits were mined. In addition, two and five stable QTLs were identified for total sugars and sucrose in both years separately. No significant QTLs were detected for total tocopherols. The results from this research provide useful knowledge about the genetic control of peanut flavor, which will aid in clarifying the molecular mechanisms of flavor research in peanuts.

Recent advances in auxin biosynthesis and homeostasis

The plant proliferation is linked with auxins which in turn play a pivotal role in the rate of growth. Also, auxin concentrations could provide insights into the age, stress, and events leading to flowering and fruiting in the sessile plant kingdom. The role in rejuvenation and plasticity is now evidenced. Interest in plant auxins spans many decades, information from different plant families for auxin concentrations, transcriptional, and epigenetic evidences for gene regulation is evaluated here, for getting an insight into pattern of auxin biosynthesis. This biosynthesis takes place via an tryptophan-independent and tryptophan-dependent pathway. The independent pathway initiated before the tryptophan (trp) production involves indole as the primary substrate. On the other hand, the trp-dependent IAA pathway passes through the indole pyruvic acid (IPyA), indole-3-acetaldoxime (IAOx), and indole acetamide (IAM) pathways. Investigations on trp-dependent pathways involved mutants, namely yucca (1-11), taa1, nit1, cyp79b and cyp79b2, vt2 and crd, and independent mutants of tryptophan, ins are compiled here. The auxin conjugates of the IAA amide and ester-linked mutant gh3, iar, ilr, ill, iamt1, ugt, and dao are remarkable and could facilitate the assimilation of auxins. Efforts are made herein to provide an up-to-date detailed information about biosynthesis leading to plant sustenance. The vast information about auxin biosynthesis and homeostasis is consolidated in this review with a simplified model of auxin biosynthesis with keys and clues for important missing links since auxins can enable the plants to proliferate and override the environmental influence and needs to be probed for applications in sustainable agriculture.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03709-6.

UGT76F1 glycosylates an isomer of the C7-necic acid component of pyrrolizidine alkaloids in Arabidopsis thaliana

Identification of unknown metabolites and their biosynthetic genes is an active research area in plant specialized metabolism. By following a gene-metabolite association from a genome-wide association study of Arabidopsis stem metabolites, we report a previously unknown metabolite, 2-hydroxy-2-(1-hydroxyethyl)pentanoic acid glucoside, and demonstrated that UGT76F1 is responsible for its production in Arabidopsis. The chemical structure of the glucoside was determined by a series of analyses, including tandem MS, acid and base hydrolysis, and NMR spectrometry. T-DNA knockout mutants of UGT76F1 are devoid of the glucoside but accumulate increased levels of the aglycone. 2-hydroxy-2-(1-hydroxyethyl)pentanoic acid is structurally related to the C7-necic acid component of lycopsamine-type pyrrolizidine alkaloids such as trachelantic acid and viridifloric acid. Feeding norvaline greatly enhances the accumulation of 2-hydroxy-2-(1-hydroxyethyl)pentanoic acid glucoside in wild-type but not the UGT76F1 knockout mutant plants, providing evidence for an orthologous C7-necic acid biosynthetic pathway in Arabidopsis despite the apparent lack of pyrrolizidine alkaloids.

Transcriptome Profiling during Sequential Stages of Cryopreservation in Banana (Musa AAA cv Borjahaji) Shoot Meristem

Cryopreservation approaches have been implemented in gene banks as a strategy to back up plant genetic resource collections that are vegetatively propagated. Different strategies have been employed to effectively cryopreserve plant tissue. There is little information on the cellular processes and molecular adjustments that confer resilience to the multiple stresses imposed during a cryoprotocol. In the present work, the cryobionomics of banana (Musa sp.), a non-model species, was investigated through the transcriptomic approach using RNA-Seq. Proliferating meristems of in vitro explants (Musa AAA cv 'Borjahaji') were cryopreserved using the droplet-vitrification technique. Transcriptome profiling analysis of eight cDNA libraries including the bio-replicates for T0 (stock cultures (control tissue), T1 (high sucrose pre-cultured), T2 (vitrification solution-treated) and T3 (liquid nitrogen-treated) meristem tissues was carried out. The raw reads obtained were mapped with a Musa acuminata reference genome sequence. A total of 70 differentially expressed genes (DEGs) comprising 34 upregulated and 36 downregulated were identified in all three phases as compared to control (T0). Among the significant DEGs (>log FC 2.0), during sequential steps, 79 in T1, 3 in T2 and the 4 in T3 were upregulated and 122 in T1, 5 in T2 and 9 in T3 were downregulated. Gene ontology (GO) enrichment analysis showed that these significant DEGs were involved in the upregulation of biological process (BP-170), cellular component (CC-10) and molecular function (MF-94) and downregulation of biological process (BP-61), cellular component (CC-3) and molecular function (MF-56). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that DEGs were involved in the biosynthesis of secondary metabolites, glycolysis/gluconeogenesis, MAPK signaling, EIN 3-lke 1 protein, 3-ketoacy-CoA synthase 6-like, and fatty acid elongation during cryopreservation. For the first time, a comprehensive transcript profiling during four stages of cryopreservation in banana were carried out, which will pave the way for devising an effective cryopreservation protocol.

UDP-Glycosyltransferases in Edible Fungi: Function, Structure, and Catalytic Mechanism

UDP-glycosyltransferases (UGTs) are the most studied glycosyltransferases, and belong to large GT1 family performing the key roles in antibiotic synthesis, the development of bacterial glycosyltransferase inhibitors, and in animal inflammation. They transfer the glycosyl groups from nucleotide UDP-sugars (UDP-glucose, UDP-galactose, UDP-xylose, and UDP-rhamnose) to the acceptors including saccharides, proteins, lipids, and secondary metabolites. The present review summarized the recent of UDP-glycosyltransferases, including their structures, functions, and catalytic mechanism, especially in edible fungi. The future perspectives and new challenges were also summarized to understand of their structure–function relationships in the future. The outputs in this field could provide a reference to recognize function, structure, and catalytic mechanism of UDP-glycosyltransferases for understanding the biosynthesis pathways of secondary metabolites, such as hydrocarbons, monoterpenes, sesquiterpene, and polysaccharides in edible fungi.

Advances in the Biosynthesis and Molecular Evolution of Steroidal Saponins in Plants

Steroidal saponins are an important type of plant-specific metabolite that are essential for plants' responses to biotic and abiotic stresses. Because of their extensive pharmacological activities, steroidal saponins are also important industrial raw materials for the production of steroidal drugs. In recent years, more and more studies have explored the biosynthesis of steroidal saponins in plants, but most of them only focused on the biosynthesis of their molecular skeleton, diosgenin, and their subsequent glycosylation modification mechanism needs to be further studied. In addition, the biosynthetic regulation mechanism of steroidal saponins, their distribution pattern, and their molecular evolution in plants remain unclear. In this review, we summarized and discussed recent studies on the biosynthesis, molecular regulation, and function of steroidal saponins. Finally, we also reviewed the distribution and molecular evolution of steroidal saponins in plants. The elucidation of the biosynthesis, regulation, and molecular evolutionary mechanisms of steroidal saponins is crucial to provide new insights and references for studying their distribution, diversity, and evolutionary history in plants. Furthermore, a deeper understanding of steroidal saponin biosynthesis will contribute to their industrial production and pharmacological applications.

Genome assembly of Musa beccarii shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes

BACKGROUND

Musa beccarii (Musaceae) is a banana species native to Borneo, sometimes grown as an ornamental plant. The basic chromosome number of Musa species is x = 7, 10, or 11; however, M. beccarii has a basic chromosome number of x = 9 (2n = 2x = 18), which is the same basic chromosome number of species in the sister genera Ensete and Musella. Musa beccarii is in the section Callimusa, which is sister to the section Musa. We generated a high-quality chromosome-scale genome assembly of M. beccarii to better understand the evolution and diversity of genomes within the family Musaceae.

FINDINGS

The M. beccarii genome was assembled by long-read and Hi-C sequencing, and genes were annotated using both long Iso-seq and short RNA-seq reads. The size of M. beccarii was the largest among all known Musaceae assemblies (∼570 Mbp) due to the expansion of transposable elements and increased 45S ribosomal DNA sites. By synteny analysis, we detected extensive genome-wide chromosome fusions and fissions between M. beccarii and the other Musa and Ensete species, far beyond those expected from differences in chromosome number. Within Musaceae, M. beccarii showed a reduced number of terpenoid synthase genes, which are related to chemical defense, and enrichment in lipid metabolism genes linked to the physical defense of the cell wall. Furthermore, type III polyketide synthase was the most abundant biosynthetic gene cluster (BGC) in M. beccarii. BGCs were not conserved in Musaceae genomes.

CONCLUSIONS

The genome assembly of M. beccarii is the first chromosome-scale genome assembly in the Callimusa section in Musa, which provides an important genetic resource that aids our understanding of the evolution of Musaceae genomes and enhances our knowledge of the pangenome.

Genome-wide identification, evolution and function analysis of UGTs superfamily in cotton

Glycosyltransferases mainly catalyse the glycosylation reaction in living organisms and widely exists in plants. UGTs have been identified from G. raimondii, G. arboreum and G. hirsutum. However, Genome-wide systematic analysis of UGTs superfamily have not been studied in G. barbadense. 752 UGTs were identified from four cotton species and grouped into 18 clades, of which R was newly discovered clades. Most UGTs were clustered at both ends of the chromosome and showed a heterogeneous distribution. UGT proteins were widely distributed in cells, with the highest distribution in chloroplasts. UGTs of the same clade shared similar intron/exon structural features. During evolution, the gene family has undergone strong selection for purification. UGTs were significantly enriched in “transcriptional activity (GO:0016758)” and “metabolic processes (GO:0008152)”. Genes from the same clade differed in function under various abiotic stresses. The analysis of cis-acting element and qRT–PCR may indicate that GHUGTs play important roles in plant growth, development and abiotic stress. We further found that GHUGT74-2 plays an important role under submergence. The study broadens the understanding of UGTs in terms of gene characteristics, evolutionary processes, and gene function in cotton and provides a new way to systematically and globally understand the structure–function relationship of multigene families in the evolutionary process.

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.

Characterization of Three Paris polyphylla Glycosyltransferases from Different UGT Families for Steroid Functionalization

Plant steroid glycosides, such as phytosterol glycosides, steroidal saponins, and steroidal glycoalkaloids, are natural products with great pharmaceutical values. In this study, we characterized three UDP-glycosyltransferases (UGTs) involved in the glycosylation of steroidal sapogenin from Paris polyphylla. Substrate specificity analysis revealed that UGT73CR1 could glycosylate steroidal sapogenins and steroidal alkaloids, with the highest affinity for diosgenin. The residues His27 and Asp129 of UGT73CR1 are conserved in corresponding positions of plant glycosyltransferases, which are crucial for activating the C-3 OH of the receptor substrates. In comparison, UGT80A33 and UGT80A34 exhibited a higher affinity for cholesterol than other steroids. UGT80s have a larger active pocket, which allows them to accommodate the side chain of sterols. In summary, we assessed three P. polyphylla glycosyltransferases from two UGT families for the functionalization of steroidal molecules, which will provide a basis for the future biomanufacturing of diverse bioactive steroid glycosides.

High Resistance to Quinclorac in Multiple-Resistant Echinochloa colona Associated with Elevated Stress Tolerance Gene Expression and Enriched Xenobiotic Detoxification Pathway

Echinochloa colona and other species in this genus are a threat to global rice production and food security. Quinclorac, an auxin mimic, is a common herbicide for grass weed control in rice, and Echinochloa spp. have evolved resistance to it. The complete mode of quinclorac action and subsequent evolution of resistance is not fully understood. We analyzed the de novo transcriptome of multiple-herbicide-resistant (ECO-R) and herbicide-susceptible genotypes in response to quinclorac. Several biological processes were constitutively upregulated in ECO-R, including carbon metabolism, photosynthesis, and ureide metabolism, indicating improved metabolic efficiency. The transcriptional change in ECO-R following quinclorac treatment indicates an efficient response, with upregulation of trehalose biosynthesis, which is also known for abiotic stress mitigation. Detoxification-related genes were induced in ECO-R, mainly the UDP-glycosyltransferase (UGT) family, most likely enhancing quinclorac metabolism. The transcriptome data also revealed that many antioxidant defense elements were uniquely elevated in ECO-R to protect against the auxin-mediated oxidative stress. We propose that upon quinclorac treatment, ECO-R detoxifies quinclorac utilizing UGT genes, which modify quinclorac using the sufficient supply of UDP-glucose from the elevated trehalose pathway. Thus, we present the first report of upregulation of trehalose synthesis and its association with the herbicide detoxification pathway as an adaptive mechanism to herbicide stress in Echinochloa, resulting in high resistance.

Characterization of UDP-glycosyltransferase family members reveals how major flavonoid glycoside accumulates in the roots of Scutellaria baicalensis

Background

Flavonoid glycosides extracted from roots of Scutellaria baicalensis exhibit strong pharmaceutical antitumor, antioxidative, anti-inflammatory, and antiviral activities. UDP glycosyltransferase (UGT) family members are responsible for the transfer of a glycosyl moiety from UDP sugars to a wide range of acceptor flavonoids. Baicalin is the major flavonoid glycoside found in S. baicalensis roots, and its aglycone baicalein is synthesized from a specially evolved pathway that has been elucidated. However, it is necessary to carry out a genome-wide study of genes involved in 7-O-glucuronidation, the final biosynthesis step of baicalin, which might elucidate the relationship between the enzymes and the metabolic accumulation patterns in this medicinal plant.

Results

We reported the phylogenetic analysis, tissue-specific expression, biochemical characterization and evolutionary analysis of glucosyltransferases (SbUGTs) and glucuronosyltransferases (SbUGATs) genes based on the recently released genome of S. baicalensis. A total of 124 UGTs were identified, and over one third of them were highly expressed in roots. In vitro enzyme assays showed that 6 SbUGTs could use UDP-glucose as a sugar donor and convert baicalein to oroxin A (baicalein 7-O-glucoside), while 4 SbUGATs used only UDP-glucuronic acid as the sugar donor and catalyzed baicalein to baicalin. SbUGAT4 and SbUGT2 are the most highly expressed SbUGAT and SbUGT genes in root tissues, respectively. Kinetic measurements revealed that SbUGAT4 had a lower Km value and higher Vmax/Km ratio to baicalein than those of SbUGT2. Furthermore, tandem duplication events were detected in SbUGTs and SbUGATs.

Conclusions

This study demonstrated that glucosylation and glucuronidation are two major glycosylated decorations in the roots of S. baicalensis. Higher expression level and affinity to substrate of SbUGAT4, and expansion of this gene family contribute high accumulation of baicalin in the root of S. baicalensis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08391-1.

Genome-Wide Identification of the TIFY Gene Family in Brassiceae and Its Potential Association with Heavy Metal Stress in Rapeseed

The TIFY gene family plays important roles in various plant biological processes and responses to stress and hormones. The chromosome-level genome of the Brassiceae species has been released, but knowledge concerning the TIFY family is lacking in the Brassiceae species. The current study performed a bioinformatics analysis on the TIFY family comparing three diploid (B. rapa, B. nigra, and B. oleracea) and two derived allotetraploid species (B. juncea, and B. napus). A total of 237 putative TIFY proteins were identified from five Brassiceae species, and classified into ten subfamilies (six JAZ types, one PPD type, two TIFY types, and one ZML type) based on their phylogenetic relationships with TIFY proteins in A. thaliana and Brassiceae species. Duplication and synteny analysis revealed that segmental and tandem duplications led to the expansion of the TIFY family genes during the process of polyploidization, and most of these TIFY family genes (TIFYs) were subjected to purifying selection after duplication based on Ka/Ks values. The spatial and temporal expression patterns indicated that different groups of BnaTIFYs have distinct spatiotemporal expression patterns under normal conditions and heavy metal stresses. Most of the JAZIII subfamily members were highest in all tissues, but JAZ subfamily members were strongly induced by heavy metal stresses. BnaTIFY34, BnaTIFY59, BnaTIFY21 and BnaTIFY68 were significantly upregulated mostly under As3+ and Cd2+ treatment, indicating that they could be actively induced by heavy metal stress. Our results may contribute to further exploration of TIFYs, and provided valuable information for further studies of TIFYs in plant tolerance to heavy metal stress.

Integrated Metabolomic and Transcriptomic Analysis Reveals Differential Mechanism of Flavonoid Biosynthesis in Two Cultivars of Angelica sinensis

Angelica sinensis is a traditional Chinese medicinal plant that has been primarily used as a blood tonic. It largely relies on its bioactive metabolites, which include ferulic acid, volatile oils, polysaccharides and flavonoids. In order to improve the yield and quality of A. sinensis, the two cultivars Mingui 1 (M1), with a purple stem, and Mingui 2 (M2), with a green stem, have been selected in the field. Although a higher root yield and ferulic acid content in M1 than M2 has been observed, the differences of flavonoid biosynthesis and stem-color formation are still limited. In this study, the contents of flavonoids and anthocyanins were determined by spectrophotometer, the differences of flavonoids and transcripts in M1 and M2 were conducted by metabolomic and transcriptomic analysis, and the expression level of candidate genes was validated by qRT-PCR. The results showed that the contents of flavonoids and anthocyanins were 1.5- and 2.6-fold greater in M1 than M2, respectively. A total of 26 differentially accumulated flavonoids (DAFs) with 19 up-regulated (UR) and seven down-regulated (DR) were obtained from the 131 identified flavonoids (e.g., flavonols, flavonoid, isoflavones, and anthocyanins) in M1 vs. M2. A total 2210 differentially expressed genes (DEGs) were obtained from the 34,528 full-length isoforms in M1 vs. M2, and 29 DEGs with 24 UR and 5 DR were identified to be involved in flavonoid biosynthesis, with 25 genes (e.g., CHS1, CHI3, F3H, DFR, ANS, CYPs and UGTs) mapped on the flavonoid biosynthetic pathway and four genes (e.g., RL1, RL6, MYB90 and MYB114) belonging to transcription factors. The differential accumulation level of flavonoids is coherent with the expression level of candidate genes. Finally, the network of DAFs regulated by DEGs was proposed. These findings will provide references for flavonoid production and cultivars selection of A. sinensis.

Broadening the roles of UDP-glycosyltransferases in auxin homeostasis and plant development

The levels of the important plant growth regulator indole-3-acetic acid (IAA) are tightly controlled within plant tissues to spatiotemporally orchestrate concentration gradients that drive plant growth and development. Metabolic inactivation of bioactive IAA is known to participate in the modulation of IAA maxima and minima. IAA can be irreversibly inactivated by oxidation and conjugation to aspartate and glutamate. Usually overlooked because of its reversible nature, the most abundant inactive IAA form is the IAA-glucose (IAA-glc) conjugate. Glycosylation of IAA in Arabidopsis thaliana is reported to be carried out by UDP-glycosyltransferase 84B1 (UGT84B1), while UGT74D1 has been implicated in the glycosylation of the irreversibly formed IAA catabolite oxIAA. Here we demonstrated that both UGT84B1 and UGT74D1 modulate IAA levels throughout plant development by dual IAA and oxIAA glycosylation. Moreover, we identified a novel UGT subfamily whose members redundantly mediate the glycosylation of oxIAA and modulate skotomorphogenic growth.

Structure-function relationship of terpenoid glycosyltransferases from plants

Covering: up to 2021Terpenoids are physiologically active substances that are of great importance to humans. Their physicochemical properties are modified by glycosylation, in terms of polarity, volatility, solubility and reactivity, and their bioactivities are altered accordingly. Significant scientific progress has been made in the functional study of glycosylated terpenes and numerous plant enzymes involved in regio- and enantioselective glycosylation have been characterized, a reaction that remains chemically challenging. Crucial clues to the mechanism of terpenoid glycosylation were recently provided by the first crystal structures of a diterpene glycosyltransferase UGT76G1. Here, we review biochemically characterized terpenoid glycosyltransferases, compare their functions and primary structures, discuss their acceptor and donor substrate tolerance and product specificity, and elaborate features of the 3D structures of the first terpenoid glycosyltransferases from plants.

Comprehensive analysis and expression profiles of cassava UDP-glycosyltransferases (UGT) family reveal their involvement in development and stress responses in cassava

UDP-glycosyltransferases (UGTs) are widely involved in plant growth and stress responses. However, UGT family are not well understood in cassava. Here, we identified 121 MeUGT genes and classified them into 14 subfamilies by phylogenetic analysis. All MeUGT proteins have typical feature of the UGTs family. Tandem duplications are the crucial driving force for the expansion of MeUGT family. Cis-Acting elements analysis uncovered those 14 kinds of cis-elements associated with biotic and abiotic stress responses. Transcriptomic and qRT-PCR analyses indicated that MeUGT genes participate in postharvest physiological deterioration of storage root and the responses of biotic and abiotic stresses. Of which, MeUGT-14/41 were significantly induced after Xam treatment. Silencing of MeUGT-14 or MeUGT-41 reduced cassava resistance to Xam, verifying the accuracy of transcriptomic data for function prediction. Together, this study characterized the MeUGTs family and revealed their potential functions, which build a solid foundation for MeUGTs associated genetic improvement of cassava.

De novo genome assembly of the potent medicinal plant Rehmannia glutinosa using nanopore technology

Rehmannia glutinosa is a potent medicinal plant with a significant importance in traditional Chinese medicine. Its root is enriched with various bioactive molecules mainly iridoids, possessing important pharmaceutical properties. However, the molecular biology and evolution of R. glutinosa have been largely unexplored. Here, we report a reference genome of R. glutinosa using Nanopore technology, Illumina and Hi-C sequencing. The assembly genome is 2.49 Gb long with a scaffold N50 length of 70 Mb and high heterozygosity (2%). Since R. glutinosa is an autotetraploid (4n = 56), the difference between each set of chromosomes is very small, and it is difficult to distinguish the two sets of chromosomes using Hi-C. Hence, only one set of the genome size was mounted to the chromosome level. Scaffolds covering 52.61% of the assembled genome were anchored on 14 pseudochromosomes. Over 67% of the genome consists of repetitive sequences dominated by Copia long terminal repeats and 48,475 protein-coding genes were predicted. Phylogenetic analysis corroborates the placement of R. glutinosa in the Orobanchaceae family. Our results indicated an independent and very recent whole genome duplication event that occurred 3.64 million year ago in the R. glutinosa lineage. Comparative genomics analysis demonstrated expansion of the UDP-dependent glycosyltransferases and terpene synthase gene families, known to be involved in terpenoid biosynthesis and diversification. Furthermore, the molecular biosynthetic pathway of iridoids has been clarified in this work. Collectively, the generated reference genome of R. glutinosa will facilitate discovery and development of important pharmacological compounds.

Characterization of bacterial community and flavor differences of different types of Douchi

According to the appearance and technology, traditional fermented Douchi can be divided into dried Douchi and wet Douchi. However, there are few reports on the difference of bacterial community structure between them or the influence of bacterial community on product flavor. In this study, high-throughput sequencing technology and electronic nose were used to measure the bacterial diversity and flavor of 40 Douchi samples, and the correlation between them was explored by multivariate statistical means combined with COG database. Results showed that the cumulative average relative abundance of Firmicutes and Proteobacteria in the samples was as high as 95.93%, and the former was the core bacteria phylum. On the whole, the dominant bacteria in Douchi were Bacillus (50.67%), Staphylococcus (14.07%), Enterococcus (2.54%), Proteus (1.61%), Brevibacillus (1.46%), Providencia (1.26%), Weissella (1.24%), and Ureibacillus (1.19%). LEfSe analysis indicated that Bacillus can be used as a biomarker in dried fermented soybeans. Meanwhile, dried samples contained more intensive aromatic substances, but were significantly lower in W6S (selectivity to hydrogen) and W3S (methane-aliph) compared with the wet samples. Aneurinibacillus and Brevibacillus were helpful to the formation of aromatic flavor in Douchi, but Vagococcus and Corynebacterium were the opposite. Gene and microbial phenotypic prediction showed that microorganisms in dried Douchi use protein more efficiently, while in wet Douchi, microbial energy metabolism was more vigorous. The pathogenic potential of microorganisms in dried samples was higher than that in wet. This study can sound the alarm for improving the safety of home-brewed Douchi and provide guidance for the subsequent screening of strains that enhance the flavor of fermented soybeans.

Genome-Wide Analysis of the YABBY Transcription Factor Family in Rapeseed (Brassica napus L.)

The YABBY family of plant-specific transcription factors play important regulatory roles during the development of leaves and floral organs, but their functions in Brassica species are incompletely understood. Here, we identified 79 YABBY genes from Arabidopsis thaliana and five Brassica species (B. rapa, B. nigra, B. oleracea, B. juncea, and B. napus). A phylogenetic analysis of YABBY proteins separated them into five clusters (YAB1–YAB5) with representatives from all five Brassica species, suggesting a high degree of conservation and similar functions within each subfamily. We determined the gene structure, chromosomal location, and expression patterns of the 21 BnaYAB genes identified, revealing extensive duplication events and gene loss following polyploidization. Changes in exon–intron structure during evolution may have driven differentiation in expression patterns and functions, combined with purifying selection, as evidenced by K_a/K_s values below 1. Based on transcriptome sequencing data, we selected nine genes with high expression at the flowering stage. qRT-PCR analysis further indicated that most BnaYAB family members are tissue-specific and exhibit different expression patterns in various tissues and organs of B. napus. This preliminary study of the characteristics of the YABBY gene family in the Brassica napus genome provides theoretical support and reference for the later functional identification of the family genes.

Roles of Alkaline Phosphatase PhoA in Algal Metabolic Regulation under Phosphorus-replete Conditions

Alkaline phosphatase (AP) in plants and algae is known to hydrolyze dissolved organophosphate (DOP) in order to obtain phosphorus when the preferred dissolved inorganic phosphorus (DIP) is present in limited supply. By conducting comparative analyses of physiologies and transcriptomes on a mutant of PhoA type AP (mPhoA) and wild type (WT) of the marine diatom Phaeodactylum tricornutum CCAP 1055/1 under P-replete and P-depleted conditions, we document other roles of this gene than DOP scavenging. PhoA mutation created by CRISPR/Cas9 diminished its DOP hydrolase activity but led to significant increases in cellular contents of pigment, carbon, and lipids, photosynthetic rate, growth rate, and the transcriptional levels of their corresponding metabolic pathways. All the results in concert indicate that besides P-nutrient scavenging under DIP deficiency, AP also functions, under the P-replete condition, to constrain pigment biosynthesis, photosynthesis, fatty acid biosynthesis, and cell division. These functions have important implications in maintaining metabolic homeostasis and preventing premature cell division.

The genome of the euryhaline rotifer Brachionus paranguensis: Potential use in molecular ecotoxicology

Brachionus spp. rotifers have been proposed as model organisms for ecotoxicological studies. We analyzed the whole-genome sequence of B. paranguensis through NextDenovo, resulting in a total length of 106.2 Mb and 71 contigs. The N50 and the GC content were 4.13 Mb and 28%, respectively. A total of 18,501 genes were predicted within the genome of B. paranguensis. Prominent detoxification-related gene families of phase I and II detoxifications have been investigated. In parallel with other Brachionus rotifers, high gene expansion was observed in CYP clan 3 and GST sigma class in B. paranguensis. Moreover, species-specific expansion of sulfotransferase (SULTs) and gain of UDP-glucuronosyltransferases (UGTs) through horizontal gene transfer has been specifically found within B. plicatilis complex. This whole-genome analysis of B. paranguensis provides a basis for molecular ecotoxicological studies and provides useful information for comparative studies of the evolution of detoxification mechanisms in Brachionus spp.

Analysis of microbial diversity and functional differences in different types of high-temperature Daqu

Bacterial communities that enrich in high-temperature Daqu are important for the Chinese maotai-flavor liquor brewing process. However, the bacterial communities in three different types of high-temperature Daqu (white Daqu, black Daqu, and yellow Daqu) are still undercharacterized. In this study, the bacterial diversity of three different types of high-temperature Daqu was investigated using Illumina MiSeq high-throughput sequencing. The bacterial community of high-temperature Daqu is mainly composed of thermophilic bacteria, and seven bacterial phyla along with 262 bacterial genera were identified in all 30 high-temperature Daqu samples. Firmicutes, Actinobacteria, Proteobacteria, and Acidobacteria were the dominant bacterial phyla in high-temperature Daqu samples, while Thermoactinomyces, Staphylococcus, Lentibacillus, Bacillus, Kroppenstedtia, Saccharopolyspora, Streptomyces, and Brevibacterium were the dominant bacterial genera. The bacterial community structure of three different types of high-temperature Daqu was significantly different (p < .05). In addition, the results of microbiome phenotype prediction by BugBase and bacterial functional potential prediction using PICRUSt show that bacteria from different types of high-temperature Daqu have similar functions as well as phenotypes, and bacteria in high-temperature Daqu have vigorous metabolism in the transport and decomposition of amino acids and carbohydrates. These results offer a reference for the comprehensive understanding of bacterial diversity of high-temperature Daqu.

Bacterial diversity and flavor profile of Zha-Chili, a traditional fermented food in China

Zha-chili is a traditional Chinese fermented food with special flavor, which is often used as an appetizer in condiments and an important energy source. The final quality of zha-chili is closely related to its microbial community structure. However, the differences of bacterial diversity in zha-chili from different regions and how bacterial species affect zha-chili fermentation process and flavor quality have not been reported. In this study, the bacterial diversity and flavor quality of zha-chili samples from different regions were analyzed using Illumina Miseq high-throughput sequencing, electronic nose and electronic tongue technology. Twenty-three bacterial phyla and 665 bacterial genera were identified in all zha-chili samples. Firmicutes, Proteobacteria and Actinobacteria were the dominant bacterial phyla in zha-chili samples, while Lactobacillus, Pseudomonas, Pediococcus, Weissella and Staphylococcus were the dominant bacterial genera. The bacterial community structure of zha-chili samples from different regions was significantly diverse (p < 0.05). The flavor of zha-chili samples also varied in different regions, and the discrepancy of taste was much greater than that of aroma. Moreover, there were significant correlations (p < 0.05) between 6 dominant bacterial genera and 8 flavor indicators (3 aroma indicators, 5 taste indicators). In addition, the results of microbiome phenotypes prediction by BugBase and bacterial functional potential prediction using PICRUSt showed that eight out of nine predicted phenotypic functions of zha-chili samples from different regions were significantly different (p < 0.05), bacterial metabolism was vigorous in the zha-chili samples, and Lactobacillus was the dominant bacterial genus involved in metabolism during fermentation.

Genome-wide characterization, evolution and expression profiling of UDP-glycosyltransferase family in pomelo (Citrus grandis) fruit

BACKGROUND

Pomelo is one of the three major species of citrus. The fruit accumulates a variety of abundant secondary metabolites that affect the flavor. UDP-glycosyltransferases (UGTs) are involved in the glycosylation of secondary metabolites.

RESULTS

In the present study, we performed a genome-wide analysis of pomelo UGT family, a total of 145 UGTs was identified based on the conserved plant secondary product glycosyltransferase (PSPG) motif. These UGT genes were clustered into 16 major groups through phylogenetic analysis of these genes with other plant UGTs (A-P). Pomelo UGTs were distributed unevenly among the chromosomes. At least 10 intron insertion events were observed in these UGT genome sequences, and I-5 was identified to be the highest conserved one. The expression profile analysis of pomelo UGT genes in different fruit tissues during development and ripening was carried out by RNA-seq.

CONCLUSIONS

We identified 145 UGTs in pomelo fruit through transcriptome data and citrus genome database. Our research provides available information on UGTs studies in pomelo, and provides an important research foundation for screening and identification of functional UGT genes.

A Genome-Wide Survey of MATE Transporters in Brassicaceae and Unveiling Their Expression Profiles under Abiotic Stress in Rapeseed

The multidrug and toxic compound extrusion (MATE) protein family is important in the export of toxins and other substrates, but detailed information on this family in the Brassicaceae has not yet been reported compared to Arabidopsis thaliana. In this study, we identified 57, 124, 81, 85, 130, and 79 MATE genes in A. thaliana, Brassica napus, Brassica oleracea, Brassica rapa, Brassica juncea, and Brassica nigra, respectively, which were unevenly distributed on chromosomes owing to both tandem and segmental duplication events. Phylogenetic analysis showed that these genes could be classified into four subgroups, shared high similarity and conservation within each group, and have evolved mainly through purifying selection. Furthermore, numerous B. napusMATE genes showed differential expression between tissues and developmental stages and between plants treated with heavy metals or hormones and untreated control plants. This differential expression was especially pronounced for the Group 2 and 3 BnaMATE genes, indicating that they may play important roles in stress tolerance and hormone induction. Our results provide a valuable foundation for the functional dissection of the different BnaMATE homologs in B. napus and its parental lines, as well as for the breeding of more stress-tolerant B. napus genotypes.

Genome-wide identification AINTEGUMENTA-like (AIL) genes in Brassica species and expression patterns during reproductive development in Brassica napus L

The AINTEGUMENTA-like (AIL) proteins, which belong to the AP2 family, play important roles in regulating the growth and development of plant organs. The AIL family has not yet been comprehensively studied in rapeseed (Brassica napus), an allotetraploid and model organism for the study of polyploid evolution. In the present study, 99 AIL family genes were identified and characterized from B. rapa, B. oleracea, B. napus, B. juncea, and B. nigra using a comprehensive genome-wide study, including analyses of phylogeny, gene structure, chromosomal localization, and expression pattern. Using a phylogenetic analysis, the AIL genes were divided into eight groups, which were closely related to the eight AtAIL genes, and which shared highly conserved structural features within the same subfamily. The non-synonymous/synonymous substitution ratios of the paralogs and orthologs were less than 1, suggesting that the AIL genes mainly experienced purifying selection during evolution. In addition, the RNA sequencing data and qRT-PCR analysis revealed that the B. napus AIL genes exhibited organ- and developmental stage-specific expression patterns. Certain genes were highly expressed in the developing seeds (BnaAIL1, BnaAIL2, BnaAIL5, and BnaAIL6), the roots (BnaANT, BnaAIL5, and BnaAIL6), and the stem (BnaAIL7B). Our results provide valuable information for further functional analysis of the AIL family in B. napus and related Brassica species.

Identification and analysis of CYP450 and UGT supergene family members from the transcriptome of Aralia elata (Miq.) seem reveal candidate genes for triterpenoid saponin biosynthesis

BACKGROUND

Members of the cytochrome P450 (CYP450) and UDP-glycosyltransferase (UGT) gene superfamily have been shown to play essential roles in regulating secondary metabolite biosynthesis. However, the systematic identification of CYP450s and UGTs has not been reported in Aralia elata (Miq.) Seem, a highly valued medicinal plant.

RESULTS

In the present study, we conducted the RNA-sequencing (RNA-seq) analysis of the leaves, stems, and roots of A. elata, yielding 66,713 total unigenes. Following annotation and KEGG pathway analysis, we were able to identify 64 unigenes related to triterpenoid skeleton biosynthesis, 254 CYP450s and 122 UGTs, respectively. A total of 150 CYP450s and 92 UGTs encoding > 300 amino acid proteins were utilized for phylogenetic and tissue-specific expression analyses. This allowed us to cluster 150 CYP450s into 9 clans and 40 families, and then these CYP450 proteins were further grouped into two primary branches: A-type (53%) and non-A-type (47%). A phylogenetic analysis of 92 UGTs and other plant UGTs led to clustering into 16 groups (A-P). We further assessed the expression patterns of these CYP450 and UGT genes across A. elata tissues, with 23 CYP450 and 16 UGT members being selected for qRT-PCR validation, respectively. From these data, we identified CYP716A295 and CYP716A296 as the candidate genes most likely to be associated with oleanolic acid synthesis, while CYP72A763 and CYP72A776 were identified as being the most likely to play roles in hederagenin biosynthesis. We also selected five unigenes as the best candidates for oleanolic acid 3-O-glucosyltransferase. Finally, we assessed the subcellular localization of three CYP450 proteins within Arabidopsis protoplasts, highlighting the fact that they localize to the endoplasmic reticulum.

CONCLUSIONS

This study presents a systematic analysis of the CYP450 and UGT gene family in A. elata and provides a foundation for further functional characterization of these two multigene families.

Genome-Wide Analysis of Phosphorus Transporter Genes in Brassica and Their Roles in Heavy Metal Stress Tolerance

Phosphorus transporter (PHT) genes encode H2PO4-/H+ co-transporters that absorb and transport inorganic nutrient elements required for plant development and growth and protect plants from heavy metal stress. However, little is known about the roles of PHTs in Brassica compared to Arabidopsis thaliana. In this study, we identified and extensively analyzed 336 PHTs from three diploid (B. rapa, B. oleracea, and B. nigra) and two allotetraploid (B. juncea and B. napus) Brassica species. We categorized the PHTs into five phylogenetic clusters (PHT1-PHT5), including 201 PHT1 homologs, 15 PHT2 homologs, 40 PHT3 homologs, 54 PHT4 homologs, and 26 PHT5 homologs, which are unevenly distributed on the corresponding chromosomes of the five Brassica species. All PHT family genes from Brassica are more closely related to Arabidopsis PHTs in the same vs. other clusters, suggesting they are highly conserved and have similar functions. Duplication and synteny analysis revealed that segmental and tandem duplications led to the expansion of the PHT gene family during the process of polyploidization and that members of this family have undergone purifying selection during evolution based on Ka/Ks values. Finally, we explored the expression profiles of BnaPHT family genes in specific tissues, at various developmental stages, and under heavy metal stress via RNA-seq analysis and qRT-PCR. BnaPHTs that were induced by heavy metal treatment might mediate the response of rapeseed to this important stress. This study represents the first genome-wide analysis of PHT family genes in Brassica species. Our findings improve our understanding of PHT family genes and provide a basis for further studies of BnaPHTs in plant tolerance to heavy metal stress.

Identification of Differentially Expressed Proteins in Sugarcane in Response to Infection by Xanthomonas albilineans Using iTRAQ Quantitative Proteomics

Sugarcane can suffer severe yield losses when affected by leaf scald, a disease caused by Xanthomonas albilineans. This bacterial pathogen colonizes the vascular system of sugarcane, which can result in reduced plant growth and plant death. In order to better understand the molecular mechanisms involved in the resistance of sugarcane to leaf scald, a comparative proteomic study was performed with two sugarcane cultivars inoculated with X. albilineans: one resistant (LCP 85-384) and one susceptible (ROC20) to leaf scald. The iTRAQ (isobaric tags for relative and absolute quantification) approach at 0 and 48 h post-inoculation (hpi) was used to identify and annotate differentially expressed proteins (DEPs). A total of 4295 proteins were associated with 1099 gene ontology (GO) terms by GO analysis. Among those, 285 were DEPs during X. albilineans infection in cultivars LCP 85-384 and ROC20. One hundred seventy-two DEPs were identified in resistant cultivar LCP 85-384, and 113 of these proteins were upregulated and 59 were downregulated. One hundred ninety-two DEPs were found in susceptible cultivar ROC20 and half of these (92) were upregulated, whereas the other half corresponded to downregulated proteins. The significantly upregulated DEPs in LCP 85-384 were involved in metabolic pathways, the biosynthesis of secondary metabolites, and the phenylpropanoid biosynthesis pathway. Additionally, the expression of seven candidate genes related to photosynthesis and glycolytic pathways, plant innate immune system, glycosylation process, plant cytochrome P450, and non-specific lipid transfer protein was verified based on transcription levels in sugarcane during infection by X. albilineans. Our findings shed new light on the differential expression of proteins in sugarcane cultivars in response to infection by X. albilineans. The identification of these genes provides important information for sugarcane variety improvement programs using molecular breeding strategies.

Phylogenomic analysis of UDP-dependent glycosyltransferases provides insights into the evolutionary landscape of glycosylation in plant metabolism

Glycosylated metabolites generated by UDP-dependent glycosyltransferases (UGTs) play critical roles in plant interactions with the environment as well as human and animal nutrition. The evolution of plant UGTs has previously been explored, but with a limited taxon sampling. In this study, 65 fully sequenced plant genomes were analyzed, and stringent criteria for selection of candidate UGTs were applied to ensure a more comprehensive taxon sampling and reliable sequence inclusion. In addition to revealing the overall evolutionary landscape of plant UGTs, the phylogenomic analysis also resolved the phylogenetic association of UGTs from free-sporing plants and gymnosperms, and identified an additional UGT group (group R) in seed plants. Furthermore, lineage-specific expansions and contractions of UGT groups were detected in angiosperms, with the total number of UGTs per genome remaining constant generally. The loss of group Q UGTs in Poales and Brassicales, rather than functional convergence in the group Q containing species, was supported by a gene tree of group Q UGTs sampled from many species, and further corroborated by the absence of group Q homologs on the syntenic chromosomal regions in Arabidopsis thaliana (Brassicales). Branch-site analyses of the group Q UGT gene tree allowed for identification of branches and amino acid sites that experienced episodic positive selection. The positively selected sites are located on the surface of a representative group Q UGT (PgUGT95B2), away from the active site, suggesting their role in protein folding/stability or protein-protein interactions.

Molecular cloning and biochemical characterization of a new flavonoid glycosyltransferase from the aquatic plant lotus

Lotus (Nelumbo nucifera Gaertn.), one of the earliest plants in angiosperms, is a perennial aquatic herb widely distributed throughout Eastern Asia. Quercetin and its glycosides are the most abundant phenolic compounds in lotus with multiple pharmacological activities. Although many flavonoid O-glycosyltransferases involved in the biosynthesis of quercetin glycosides have been identified from terrestrial plants, no glycosyltransferase has been identified in aquatic plants. In this study, a new glycosyltransferase (NpUGT6) was identified from the embryo of Nelumbo nucifera (Nelumbinis Plumula). Function characterization demonstrated that NpUGT6 exhibited a robust capability to regio- and stereo-specific O-glycosylation at the 3-hydroxy group of quercetin scaffolds with UDP-glucose. Moreover, the O-glycosylation catalyzed by NpUGT6 was reversible. NpUGT6 is the first identified flavonoid O-glycosyltransferase from aquatic plants. Its sequence will provide useful guidance for the discovery of additional flavonoid glycosyltransferses in Nymphaeaceae and other aquatic plants.

Molecular Regulation of Catalpol and Acteoside Accumulation in Radial Striation and non-Radial Striation of Rehmannia glutinosa Tuberous Root

Rehmannia glutinosa L., a perennial plant of Scrophulariaceae, is one of the most commonly used herbs in traditional Chinese medicine (TCM) that have been widely cultivated in China. However, to date, the biosynthetic pathway of its two quality-control components, catalpol and acteoside, are only partially elucidated and the mechanism for their tissue-specific accumulation remains unknown. To facilitate the basic understanding of the key genes and transcriptional regulators involved in the biosynthesis of catalpol and acteoside, transcriptome sequencing of radial striation (RS) and non-radial striation (nRS) from four R. glutinosa cultivars was performed. A total of 715,158,202 (~107.27 Gb) high quality reads obtained using paired-end Illumina sequencing were de novo assembled into 150,405 transcripts. Functional annotation with multiple public databases identified 155 and 223 unigenes involved in catalpol and acteoside biosynthesis, together with 325 UGTs, and important transcription factor (TF) families. Comparative analysis of the transcriptomes identified 362 unigenes, found to be differentially expressed in all RS vs. nRS comparisons, with 143 upregulated unigenes, including those encoding enzymes of the catalpol and acteoside biosynthetic pathway, such as geranyl diphosphate synthase (RgGPPS), geraniol 8-hydroxylase (RgG10H), and phenylalanine ammonia-lyase (RgPAL). Other differentially expressed unigenes predicted to be related to catalpol and acteoside biosynthesis fall into UDP-dependent glycosyltransferases (UGTs), as well as transcription factors. In addition, 16 differentially expressed genes were selectively confirmed by real-time PCR. In conclusion, a large unigene dataset of R. glutinosa generated in the current study will serve as a resource for the identification of potential candidate genes for investigation of the tuberous root development and biosynthesis of active components.

Elucidation of the Amygdalin Pathway Reveals the Metabolic Basis of Bitter and Sweet Almonds (Prunus dulcis)

Almond (Prunus dulcis) is the principal Prunus species in which the consumed and thus commercially important part of the fruit is the kernel. As a result of continued selection, the vast majority of almonds have a nonbitter kernel. However, in the field, there are trees carrying bitter kernels, which are toxic to humans and, consequently, need to be removed. The toxicity of bitter almonds is caused by the accumulation of the cyanogenic diglucoside amygdalin, which releases toxic hydrogen cyanide upon hydrolysis. In this study, we identified and characterized the enzymes involved in the amygdalin biosynthetic pathway: PdCYP79D16 and PdCYP71AN24 as the cytochrome P450 (CYP) enzymes catalyzing phenylalanine-to-mandelonitrile conversion, PdUGT94AF3 as an additional monoglucosyl transferase (UGT) catalyzing prunasin formation, and PdUGT94AF1 and PdUGT94AF2 as the two enzymes catalyzing amygdalin formation from prunasin. This was accomplished by constructing a sequence database containing UGTs known, or predicted, to catalyze a β(1→6)-O-glycosylation reaction and a Basic Local Alignment Search Tool search of the draft version of the almond genome versus these sequences. Functional characterization of candidate genes was achieved by transient expression in Nicotiana benthamiana Reverse transcription quantitative polymerase chain reaction demonstrated that the expression of PdCYP79D16 and PdCYP71AN24 was not detectable or only reached minute levels in the sweet almond genotype during fruit development, while it was high and consistent in the bitter genotype. Therefore, the basis for the sweet kernel phenotype is a lack of expression of the genes encoding the two CYPs catalyzing the first steps in amygdalin biosynthesis.

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.

Genome-wide analysis of family-1 UDP glycosyltransferases (UGT) and identification of UGT genes for FHB resistance in wheat (Triticum aestivum L.)

BACKGROUND

Fusarium head blight (FHB), a devastating disease in wheat worldwide, results in yield loses and mycotoxin, such as deoxynivalenol (DON), accumulation in infected grains. DON also facilitates the pathogen colonization and spread of FHB symptoms during disease development. UDP-glycosyltransferase enzymes (UGTs) are known to contribute to detoxification and enhance FHB resistance by glycosylating DON into DON-3-glucoside (D3G) in wheat. However, a comprehensive investigation of wheat (Triticum aestivum) UGT genes is still lacking.

RESULTS

In this study, we carried out a genome-wide analysis of family-1 UDP glycosyltransferases in wheat based on the PSPG conserved box that resulted in the identification of 179 putative UGT genes. The identified genes were clustered into 16 major phylogenetic groups with a lack of phylogenetic group K. The UGT genes were invariably distributed among all the chromosomes of the 3 genomes. At least 10 intron insertion events were found in the UGT sequences, where intron 4 was observed as the most conserved intron. The expression analysis of the wheat UGT genes using both online microarray data and quantitative real-time PCR verification suggested the distinct role of UGT genes in different tissues and developmental stages. The expression of many UGT genes was up-regulated after Fusarium graminearum inoculation, and six of the genes were further verified by RT-qPCR.

CONCLUSION

We identified 179 UGT genes from wheat using the available sequenced wheat genome. This study provides useful insight into the phylogenetic structure, distribution, and expression patterns of family-1 UDP glycosyltransferases in wheat. The results also offer a foundation for future work aimed at elucidating the molecular mechanisms underlying the resistance to FHB and DON accumulation.

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.

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.

Association mapping for total polyphenol content, total flavonoid content and antioxidant activity in barley

BACKGROUND

The interest has been increasing on the phenolic compounds in plants because of their nutritive function as food and the roles regulating plant growth. However, their underlying genetic mechanism in barley is still not clear.

RESULTS

A genome-wide association study (GWAS) was conducted for total phenolic content (TPC), total flavonoid content (FLC) and antioxidant activity (AOA) in 67 cultivated and 156 Tibetan wild barley genotypes. Most markers associated with phenolic content were different in cultivated and wild barleys. The markers bPb-0572 and bPb-4531 were identified as the major QTLs controlling phenolic compounds in Tibetan wild barley. Moreover, the marker bPb-4531 was co-located with the UDP- glycosyltransferase gene (HvUGT), which is a homolog to Arabidopsis UGTs and involved in biosynthesis of flavonoid glycosides .

CONCLUSIONS

GWAS is an efficient tool for exploring the genetic architecture of phenolic compounds in the cultivated and Tibetan wild barleys. The DArT markers applied in this study can be used in barley breeding for developing new barley cultivars with higher phenolics content. The candidate gene (HvUGT) provides a potential route for deep understanding of the molecular mechanism of flavonoid synthesis.