Automatic Identification of Players in the Flavonoid Biosynthesis with Application on the Biomedicinal Plant Croton tiglium

Isoform-resolved genome annotation enables mapping of tissue-specific betalain regulation in amaranth

Betalains are coloring pigments produced in some families of the order Caryophyllales, where they replace anthocyanins as coloring pigments. While the betalain pathway itself is well studied, the tissue-specific regulation of the pathway remains mostly unknown. We enhance the high-quality Amaranthus hypochondriacus reference genome and produce a substantially more complete genome annotation, incorporating isoform details. We annotate betalain and anthocyanin pathway genes along with their regulators in amaranth and map the genetic control and tissue-specific regulation of the betalain pathway. Our improved genome annotation allowed us to identify causal mutations that lead to a knock-out of red betacyanins in natural accessions of amaranth. We reveal the tissue-specific regulation of flower color via a previously uncharacterized MYB transcription factor, AhMYB2. Downregulation of AhMYB2 in the flower leads to reduced expression of key betalain enzyme genes and loss of red flower color. Our improved amaranth reference genome represents the most complete genome of amaranth to date and is a valuable resource for betalain and amaranth research. High similarity of the flower betalain regulator AhMYB2 to anthocyanin regulators and a partially conserved interaction motif support the co-option of anthocyanin regulators for the betalain pathway as a possible reason for the mutual exclusiveness of the two pigments.

Genetic factors explaining anthocyanin pigmentation differences

BACKGROUND

Anthocyanins are important contributors to coloration across a wide phylogenetic range of plants. Biological functions of anthocyanins span from reproduction to protection against biotic and abiotic stressors. Owing to a clearly visible phenotype of mutants, the anthocyanin biosynthesis and its sophisticated regulation have been studied in numerous plant species. Genes encoding the anthocyanin biosynthesis enzymes are regulated by a transcription factor complex comprising MYB, bHLH and WD40 proteins.

RESULTS

A systematic comparison of anthocyanin-pigmented vs. non-pigmented varieties was performed within numerous plant species covering the taxonomic diversity of flowering plants. The literature was screened for cases in which genetic factors causing anthocyanin loss were reported. Additionally, transcriptomic data sets from four previous studies were reanalyzed to determine the genes possibly responsible for color variation based on their expression pattern. The contribution of different structural and regulatory genes to the intraspecific pigmentation differences was quantified. Differences concerning transcription factors are by far the most frequent explanation for pigmentation differences observed between two varieties of the same species. Among the transcription factors in the analyzed cases, MYB genes are significantly more prone to account for pigmentation differences compared to bHLH or WD40 genes. Among the structural genes, DFR genes are most often associated with anthocyanin loss.

CONCLUSIONS

These findings support previous assumptions about the susceptibility of transcriptional regulation to evolutionary changes and its importance for the evolution of novel coloration phenotypes. Our findings underline the particular significance of MYBs and their apparent prevalent role in the specificity of the MBW complex.

Multiple mechanisms explain loss of anthocyanins from betalain-pigmented Caryophyllales, including repeated wholesale loss of a key anthocyanidin synthesis enzyme

In this study, we investigate the genetic mechanisms responsible for the loss of anthocyanins in betalain-pigmented Caryophyllales, considering our hypothesis of multiple transitions to betalain pigmentation. Utilizing transcriptomic and genomic datasets across 357 species and 31 families, we scrutinize 18 flavonoid pathway genes and six regulatory genes spanning four transitions to betalain pigmentation. We examined evidence for hypotheses of wholesale gene loss, modified gene function, altered gene expression, and degeneration of the MBW (MYB-bHLH-WD40) trasnscription factor complex, within betalain-pigmented lineages. Our analyses reveal that most flavonoid synthesis genes remain conserved in betalain-pigmented lineages, with the notable exception of TT19 orthologs, essential for the final step in anthocyanidin synthesis, which appear to have been repeatedly and entirely lost. Additional late-stage flavonoid pathway genes upstream of TT19 also manifest strikingly reduced expression in betalain-pigmented species. Additionally, we find repeated loss and alteration in the MBW transcription complex essential for canonical anthocyanin synthesis. Consequently, the loss and exclusion of anthocyanins in betalain-pigmented species appear to be orchestrated through several mechanisms: loss of a key enzyme, downregulation of synthesis genes, and degeneration of regulatory complexes. These changes have occurred iteratively in Caryophyllales, often coinciding with evolutionary transitions to betalain pigmentation.

Automatic annotation of the bHLH gene family in plants

BACKGROUND

The bHLH transcription factor family is named after the basic helix-loop-helix (bHLH) domain that is a characteristic element of their members. Understanding the function and characteristics of this family is important for the examination of a wide range of functions. As the availability of genome sequences and transcriptome assemblies has increased significantly, the need for automated solutions that provide reliable functional annotations is emphasised.

RESULTS

A phylogenetic approach was adapted for the automatic identification and functional annotation of the bHLH transcription factor family. The bHLH_annotator, designed for the automated functional annotation of bHLHs, was implemented in Python3. Sequences of bHLHs described in literature were collected to represent the full diversity of bHLH sequences. Previously described orthologs form the basis for the functional annotation assignment to candidates which are also screened for bHLH-specific motifs. The pipeline was successfully deployed on the two Arabidopsis thaliana accessions Col-0 and Nd-1, the monocot species Dioscorea dumetorum, and a transcriptome assembly of Croton tiglium. Depending on the applied search parameters for the initial candidates in the pipeline, species-specific candidates or members of the bHLH family which experienced domain loss can be identified.

CONCLUSIONS

The bHLH_annotator allows a detailed and systematic investigation of the bHLH family in land plant species and classifies candidates based on bHLH-specific characteristics, which distinguishes the pipeline from other established functional annotation tools. This provides the basis for the functional annotation of the bHLH family in land plants and the systematic examination of a wide range of functions regulated by this transcription factor family.

KIPEs3: Automatic annotation of biosynthesis pathways

Flavonoids and carotenoids are pigments involved in stress mitigation and numerous other processes. Both pigment classes can contribute to flower and fruit coloration. Flavonoid aglycones and carotenoids are produced by a pathway that is largely conserved across land plants. Glycosylations, acylations, and methylations of the flavonoid aglycones can be species-specific and lead to a plethora of biochemically diverse flavonoids. We previously developed KIPEs for the automatic annotation of biosynthesis pathways and presented an application on the flavonoid aglycone biosynthesis. KIPEs3 is an improved version with additional features and the potential to identify not just the core biosynthesis players, but also candidates involved in the decoration steps and in the transport of flavonoids. Functionality of KIPEs3 is demonstrated through the analysis of the flavonoid biosynthesis in Arabidopsis thaliana Nd-1, Capsella grandiflora, and Dioscorea dumetorum. We demonstrate the applicability of KIPEs to other pathways by adding the carotenoid biosynthesis to the repertoire. As a technical proof of concept, the carotenoid biosynthesis was analyzed in the same species and Daucus carota. KIPEs3 is available as an online service to enable access without prior bioinformatics experience. KIPEs3 facilitates the automatic annotation and analysis of biosynthesis pathways with a consistent and high quality in a large number of plant species. Numerous genome sequencing projects are generating a huge amount of data sets that can be analyzed to identify evolutionary patterns and promising candidate genes for biotechnological and breeding applications.

Identification of key genes responsible for green and white colored spathes in Anthurium andraeanum (Hort.)

Modern anthuriums, Anthurium andraeanum (Hort.) are among the most popular flowering plants and widely used for interior decoration. Their popularity is largely attributed to the exotic spathes with different colors. Previous studies have reported color development in red spathe cultivars, but limited information is available on key genes regulating white and green colored spathes. This study analyzed anthocyanin, chlorophyll, and carotenoid contents as well as transcript differences in spathes of eight cultivars that differed in spathe colors ranging from red to white and green. Results showed that increased expression of a transcription factor AaMYB2 was associated with elevated levels of anthocyanin in spathes, but decreased expression of AaMYB2 and increased expression of AaLAR (leucoanthocyanidin reductase) and AaANR (anthocyanidin reductase) were accompanied with the accumulation of colorless proanthocyanidin, thus the white spathe. As to the green colored spathe, chlorophyll content in the green spathe cultivar was substantially higher than the other cultivars. Correspondingly, transcripts of chlorophyll biosynthesis-related genes AaHemB (porphobilinogen synthase) and AaPor (protochlorophyllide oxidoreductase) were highly upregulated but almost undetectable in white and red spathes. The increased expression of AaHemB and AaPor was correlated with the expression of transcription factor AaMYB124. Subsequently, qRT-PCR analysis confirmed their expression levels in nine additional cultivars with red, white, and green spathes. A working model for the formation of white and green spathes was proposed. White colored spathes are likely due to the decreased expression of AaMYB2 which results in increased expression of AaLAR and AaANR, and the green spathes are attributed to AaMYB124 enhanced expression of AaHemB and AaPor. Further research is warranted to test this working model.

Lanthanum Silicate Nanomaterials Enhance Sheath Blight Resistance in Rice: Mechanisms of Action and Soil Health Evaluation

In the current study, foliar spray with lanthanum (La) based nanomaterials (La10Si6O27 nanorods, La10Si6O27 nanoparticle, La(OH)3 nanorods, and La2O3 nanoparticle) suppressed the occurrence of sheath blight (Rhizoctonia solani) in rice. The beneficial effects were morphology-, composition-, and concentration-dependent. Foliar application of La10Si6O27 nanorods (100 mg/L) yielded the greatest disease suppression, significantly decreasing the disease severity by 62.4% compared with infected controls; this level of control was 2.7-fold greater than the commercially available pesticide (Thifluzamide). The order of efficacy was as follows: La10Si6O27 nanorods > La10Si6O27 nanoparticle > La(OH)3 nanorods > La2O3 nanoparticle. Mechanistically, (1) La10Si6O27 nanorods had greater bioavailability, slower dissolution, and simultaneous Si nutrient benefits; (2) transcriptomic and metabolomic analyses revealed that La10Si6O27 nanorods simultaneously strengthened rice systemic acquired resistance, physical barrier formation, and antioxidative systems. Additionally, La10Si6O27 nanorods improved rice yield by 35.4% and promoted the nutritional quality of the seeds as compared with the Thifluzamide treatment. A two-year La10Si6O27 nanorod exposure had no effect on soil health based on the evaluated chemical, physical, and biological soil properties. These findings demonstrate that La based nanomaterials can serve as an effective and sustainable strategy to safeguard crops and highlight the importance of nanomaterial composition and morphology in terms of optimizing benefit.

Integrated transcriptomics and metabolomics analysis provide insight into anthocyanin biosynthesis for sepal color formation in Heptacodium miconioides

Heptacodium miconioides Rehd., commonly known as "seven-son flower," is an ornamental species with a beautiful flower pattern and persistent sepals. Its sepals are of horticultural value, turning bright red and elongating in the autumn; however, the molecular mechanisms that cause sepal color change remain unclear. We analyzed the dynamic changes in anthocyanin composition in the sepal of H. miconioides at four developmental stages (S1-S4). A total of 41 anthocyanins were detected and classified into 7 major anthocyanin aglycones. High levels of the pigments cyanidin-3,5-O-diglucoside, cyanidin-3-O-galactoside, cyanidin-3-O-glucoside, and pelargonidin-3-O-glucoside were responsible for sepal reddening. Transcriptome analysis revealed 15 differentially expressed genes involved in anthocyanin biosynthesis that were detected between 2 developmental stages. Of these, the high expression of HmANS was considered critical structural gene related to anthocyanin biosynthesis pathway in the sepal through co-expression analysis with anthocyanin content. In addition, a transcription factor (TF)-metabolite correlation analysis revealed that three HmMYB, two HmbHLH, two HmWRKY, and two HmNAC TFs exhibited a strong positive role in the regulation of the anthocyanin structural genes (Pearson's correlation coefficient > 0.90). Luciferase activity assay showed that HmMYB114, HmbHLH130, HmWRKY6, and HmNAC1 could activate the promoters of HmCHS4 and HmDFR1 genes in vitro. These findings increase our understanding of anthocyanin metabolism in the sepal of H. miconioides and provide a guide for studies involving sepal color conversion and regulation.

Apiaceae FNS I originated from F3H through tandem gene duplication

BACKGROUND

Flavonoids are specialized metabolites with numerous biological functions in stress response and reproduction of plants. Flavones are one subgroup that is produced by the flavone synthase (FNS). Two distinct enzyme families evolved that can catalyze the biosynthesis of flavones. While the membrane-bound FNS II is widely distributed in seed plants, one lineage of soluble FNS I appeared to be unique to Apiaceae species.

RESULTS

We show through phylogenetic and comparative genomic analyses that Apiaceae FNS I evolved through tandem gene duplication of flavanone 3-hydroxylase (F3H) followed by neofunctionalization. Currently available datasets suggest that this event happened within the Apiaceae in a common ancestor of Daucus carota and Apium graveolens. The results also support previous findings that FNS I in the Apiaceae evolved independent of FNS I in other plant species.

CONCLUSION

We validated a long standing hypothesis about the evolution of Apiaceae FNS I and predicted the phylogenetic position of this event. Our results explain how an Apiaceae-specific FNS I lineage evolved and confirm independence from other FNS I lineages reported in non-Apiaceae species.

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.

The evidence for anthocyanins in the betalain-pigmented genus Hylocereus is weak

Here we respond to Zhou (BMC Genomics 21:734, 2020) “Combined Transcriptome and Metabolome analysis of Pitaya fruit unveiled the mechanisms underlying peel and pulp color formation” published in BMC Genomics. Given the evolutionary conserved anthocyanin biosynthesis pathway in betalain-pigmented species, we are open to the idea that species with both anthocyanins and betalains might exist. However, in absence of LC-MS/MS spectra, apparent lack of biological replicates, and no comparison to authentic standards, the findings of Zhou (BMC Genomics 21:734, 2020) are not a strong basis to propose the presence of anthocyanins in betalain-pigmented pitaya. In addition, our re-analysis of the datasets indicates the misidentification of important genes and the omission of key flavonoid and anthocyanin synthesis genes ANS and DFR. Finally, our re-analysis of the RNA-Seq dataset reveals no correlation between anthocyanin biosynthesis gene expression and pigment status.

Dissecting the genetic basis of bioactive metabolites and fruit quality traits in blueberries (Vaccinium corymbosum L.)

Blueberry is well-recognized as a healthy fruit with functionality derived largely from anthocyanin and chlorogenic acid. Despite their importance, no study to date has evaluated the genetic basis of these bioactives in blueberries and their relationship with fruit quality traits. Hence, to fill this gap, a mapping population including 196 F₁ individuals was phenotyped for anthocyanin and chlorogenic acid concentration and fruit quality traits (titratable acidity, pH, and total soluble solids) over 3 years and data were used for QTL mapping and correlation analysis. Total soluble solids and chlorogenic acid were positively correlated with glycosylated anthocyanin and total anthocyanin, respectively, indicating that parallel selection for these traits is possible. Across all the traits, a total of 188 QTLs were identified on chromosomes 1, 2, 4, 8, 9, 11 and 12. Notably, four major regions with overlapping major-effect QTLs were identified on chromosomes 1, 2, 4 and 8, and were responsible for acylation and glycosylation of anthocyanins in a substrate and sugar donor specific manner. Through comparative transcriptome analysis, multiple candidate genes were identified for these QTLs, including glucosyltransferases and acyltransferases. Overall, the study provides the first insights into the genetic basis controlling anthocyanins accumulation and composition, chlorogenic acid and fruit quality traits, and establishes a framework to advance genetic studies and molecular breeding for anthocyanins in blueberry.

Analysis of flavonol regulator evolution in the Brassicaceae reveals MYB12, MYB111 and MYB21 duplications and MYB11 and MYB24 gene loss

Background

Flavonols are the largest subgroup of flavonoids, possessing multiple functions in plants including protection against ultraviolet radiation, antimicrobial activities, and flower pigmentation together with anthocyanins. They are of agronomical and economical importance because the major off-taste component in rapeseed protein isolates is a flavonol derivative, which limits rapeseed protein use for human consumption. Flavonol production in Arabidopsis thaliana is mainly regulated by the subgroup 7 (SG7) R2R3-MYB transcription factors MYB11, MYB12, and MYB111. Recently, the SG19 MYBs MYB21, MYB24, and MYB57 were shown to regulate flavonol accumulation in pollen and stamens. The members of each subgroup are closely related, showing gene redundancy and tissue-specific expression in A. thaliana. However, the evolution of these flavonol regulators inside the Brassicaceae, especially inside the Brassiceae, which include the rapeseed crop species, is not fully understood.

Results

We studied the SG7 and SG19 MYBs in 44 species, including 31 species of the Brassicaceae, by phylogenetic analyses followed by synteny and gene expression analyses. Thereby we identified a deep MYB12 and MYB111 duplication inside the Brassicaceae, which likely occurred before the divergence of Brassiceae and Thelypodieae. These duplications of SG7 members were followed by the loss of MYB11 after the divergence of Eruca vesicaria from the remaining Brassiceae species. Similarly, MYB21 experienced duplication before the emergence of the Brassiceae tribe, where the gene loss of MYB24 is also proposed to have happened. The members of each subgroup revealed frequent overlapping spatio-temporal expression patterns in the Brassiceae member B. napus, which are assumed to compensate for the loss of MYB11 and MYB24 in the analysed tissues.

Conclusions

We identified a duplication of MYB12, MYB111, and MYB21 inside the Brassicaceae and MYB11 and MYB24 gene loss inside the tribe Brassiceae. We propose that polyploidization events have shaped the evolution of the flavonol regulators in the Brassicaceae, especially in the Brassiceae.

Transcriptomics and Metabolomics Analyses Reveal Defensive Responses and Flavonoid Biosynthesis of Dracaena cochinchinensis (Lour.) S. C. Chen under Wound Stress in Natural Conditions

Dracaena cochinchinensis has special defensive reactions against wound stress. Under wound stress, D. cochinchinensis generates a resin that is an important medicine known as dragon’s blood. However, the molecular mechanism underlying the defensive reactions is unclear. Metabolomics and transcriptomics analyses were performed on stems of D. cochinchinensis at different timepoints from the short term to the long term after wounding. According to the 378 identified compounds, wound-induced secondary metabolic processes exhibited three-phase characteristics: short term (0–5 days), middle term (10 days–3 months), and long term (6–17 months). The wound-induced transcriptome profile exhibited characteristics of four stages: within 24 h, 1–5 days, 10–30 days, and long term. The metabolic regulation in response to wound stress mainly involved the TCA cycle, glycolysis, starch and sucrose metabolism, phenylalanine biosynthesis, and flavonoid biosynthesis, along with some signal transduction pathways, which were all well connected. Flavonoid biosynthesis and modification were the main reactions against wound stress, mainly comprising 109 flavonoid metabolites and 93 wound-induced genes. A group of 21 genes encoding CHS, CHI, DFR, PPO, OMT, LAR, GST, and MYBs were closely related to loureirin B and loureirin C. Wound-induced responses at the metabolome and transcriptome level exhibited phase characteristics. Complex responses containing primary metabolism and flavonoid biosynthesis are involved in the defense mechanism against wound stress in natural conditions, and flavonoid biosynthesis and modification are the main strategies of D. cochinchinensis in the long-term responses to wound stress.

Mapping‑by‑Sequencing Reveals Genomic Regions Associated with Seed Quality Parameters in Brassica napus

Rapeseed (Brassica napus L.) is an important oil crop and has the potential to serve as a highly productive source of protein. This protein exhibits an excellent amino acid composition and has high nutritional value for humans. Seed protein content (SPC) and seed oil content (SOC) are two complex quantitative and polygenic traits which are negatively correlated and assumed to be controlled by additive and epistatic effects. A reduction in seed glucosinolate (GSL) content is desired as GSLs cause a stringent and bitter taste. The goal here was the identification of genomic intervals relevant for seed GSL content and SPC/SOC. Mapping by sequencing (MBS) revealed 30 and 15 new and known genomic intervals associated with seed GSL content and SPC/SOC, respectively. Within these intervals, we identified known but also so far unknown putatively causal genes and sequence variants. A 4 bp insertion in the MYB28 homolog on C09 shows a significant association with a reduction in seed GSL content. This study provides insights into the genetic architecture and potential mechanisms underlying seed quality traits, which will enhance future breeding approaches in B. napus.

Integrative Analysis of Metabolome and Transcriptome Identifies Potential Genes Involved in the Flavonoid Biosynthesis in Entada phaseoloides Stem

Entada phaseoloides stem is known for its high medicinal benefits and ornamental value. Flavonoids are one of the main active constituents in E. phaseoloides stem. However, the regulatory mechanism of flavonoids accumulation in E. phaseoloides is lacking. Here, phytochemical compounds and transcripts from stems at different developmental stages in E. phaseoloides were investigated by metabolome and transcriptome analysis. The metabolite profiling of the oldest stem was obviously different from young and older stem tissues. A total of 198 flavonoids were detected, and flavones, flavonols, anthocyanins, isoflavones, and flavanones were the main subclasses. The metabolome data showed that the content of acacetin was significantly higher in the young stem and older stem than the oldest stem. Rutin and myricitrin showed significantly higher levels in the oldest stem. A total of 143 MYBs and 143 bHLHs were identified and classified in the RNA-seq data. Meanwhile, 34 flavonoid biosynthesis structural genes were identified. Based on the expression pattern of structural genes involved in flavonoid biosynthesis, it indicated that flavonol, anthocyanin, and proanthocyanin biosynthesis were first active during the development of E. phaseoloides stem, and the anthocyanin or proanthocyanin biosynthesis branch was dominant; the flavone biosynthesis branch was active at the late developmental stage of the stem. Through the correlation analysis of transcriptome and metabolome data, the potential candidate genes related to regulating flavonoid synthesis and transport were identified. Among them, the MYBs, bHLH, and TTG1 are coregulated biosynthesis of flavonols and structural genes, bHLH and transporter genes are coregulated biosynthesis of anthocyanins. In addition, the WDR gene TTG1-like (AN11) may regulate dihydrochalcones and flavonol biosynthesis in specific combinations with IIIb bHLH and R2R3-MYB proteins. Furthermore, the transport gene protein TRANSPARENT TESTA 12-like gene is positively regulated the accumulation of rutin, and the homolog of ABC transporter B family member gene is positively correlated with the content of flavone acacetin. This study offered candidate genes involved in flavonoid biosynthesis, information of flavonoid composition and characteristics of flavonoids accumulation, improved our understanding of the MYBs and bHLHs-related regulation networks of flavonoid biosynthesis in E. phaseoloides stem, and provided references for the metabolic engineering of flavonoid biosynthesis in E. phaseoloides stem.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

Biochemistry and Molecular Basis of Intracellular Flavonoid Transport in Plants

Flavonoids are a biochemically diverse group of specialized metabolites in plants that are derived from phenylalanine. While the biosynthesis of the flavonoid aglycone is highly conserved across species and well characterized, numerous species-specific decoration steps and their relevance remained largely unexplored. The flavonoid biosynthesis takes place at the cytosolic side of the endoplasmatic reticulum (ER), but accumulation of various flavonoids was observed in the central vacuole. A universal explanation for the subcellular transport of flavonoids has eluded researchers for decades. Current knowledge suggests that a glutathione S-transferase-like protein (ligandin) protects anthocyanins and potentially proanthocyanidin precursors during the transport to the central vacuole. ABCC transporters and to a lower extend MATE transporters sequester anthocyanins into the vacuole. Glycosides of specific proanthocyanidin precursors are sequestered through MATE transporters. A P-ATPase in the tonoplast and potentially other proteins generate the proton gradient that is required for the MATE-mediated antiport. Vesicle-mediated transport of flavonoids from the ER to the vacuole is considered as an alternative or additional route.

Automatic identification and annotation of MYB gene family members in plants

BACKGROUND

MYBs are among the largest transcription factor families in plants. Consequently, members of this family are involved in a plethora of processes including development and specialized metabolism. The MYB families of many plant species were investigated in the last two decades since the first investigation looked at Arabidopsis thaliana. This body of knowledge and characterized sequences provide the basis for the identification, classification, and functional annotation of candidate sequences in new genome and transcriptome assemblies.

RESULTS

A pipeline for the automatic identification and functional annotation of MYBs in a given sequence data set was implemented in Python. MYB candidates are identified, screened for the presence of a MYB domain and other motifs, and finally placed in a phylogenetic context with well characterized sequences. In addition to technical benchmarking based on existing annotation, the transcriptome assembly of Croton tiglium and the annotated genome sequence of Castanea crenata were screened for MYBs. Results of both analyses are presented in this study to illustrate the potential of this application. The analysis of one species takes only a few minutes depending on the number of predicted sequences and the size of the MYB gene family. This pipeline, the required bait sequences, and reference sequences for a classification are freely available on github: https://github.com/bpucker/MYB_annotator .

CONCLUSIONS

This automatic annotation of the MYB gene family in novel assemblies makes genome-wide investigations consistent and paves the way for comparative studies in the future. Candidate genes for in-depth analyses are presented based on their orthology to previously characterized sequences which allows the functional annotation of the newly identified MYBs with high confidence. The identification of orthologs can also be harnessed to detect duplication and deletion events.

The UV-B-Induced Transcription Factor HY5 Regulated Anthocyanin Biosynthesis in Zanthoxylum bungeanum

Pericarp color is an important economic characteristic of Zanthoxylum bungeanum. Anthocyanins are the main reason for the pericarp’s red appearance in Z. bungeanum. In this study, through the combined analysis of the metabolome and transcriptome, HY5, whose expression is highly correlated to changes in the anthocyanin content, was screened and identified. Under natural ripening conditions, the Z. bungeanum fruit gradually changed in color from green to red, while bagging resulted in the fruit maintaining its green color. After unbagging, the fruit gradually turned red, and the ZbHY5 expression and anthocyanin content increased. In addition, the leaves changed from green to red after exposure to UV-B radiation, and the ZbHY5 expression and anthocyanin content increased. The transient overexpression of ZbHY5 deepened the redness of the Z. bungeanum leaves and promoted the expression of ZbHY5 and ZbMYB113 as well as anthocyanin accumulation. Bimolecular fluorescence complementation (BIFC) showed that there was an interaction between ZbHY5 and ZbMYB113. These results revealed that under UV-B irradiation, ZbHY5 might regulate the expression levels of the structural genes related to anthocyanin biosynthesis through combination with ZbMYB113, thereby affecting anthocyanin accumulation. This finding provides useful insights for further studies focusing on UV-B-induced anthocyanin accumulation in Z. bungeanum.

PlantTribes2: Tools for comparative gene family analysis in plant genomics

Plant genome-scale resources are being generated at an increasing rate as sequencing technologies continue to improve and raw data costs continue to fall; however, the cost of downstream analyses remains large. This has resulted in a considerable range of genome assembly and annotation qualities across plant genomes due to their varying sizes, complexity, and the technology used for the assembly and annotation. To effectively work across genomes, researchers increasingly rely on comparative genomic approaches that integrate across plant community resources and data types. Such efforts have aided the genome annotation process and yielded novel insights into the evolutionary history of genomes and gene families, including complex non-model organisms. The essential tools to achieve these insights rely on gene family analysis at a genome-scale, but they are not well integrated for rapid analysis of new data, and the learning curve can be steep. Here we present PlantTribes2, a scalable, easily accessible, highly customizable, and broadly applicable gene family analysis framework with multiple entry points including user provided data. It uses objective classifications of annotated protein sequences from existing, high-quality plant genomes for comparative and evolutionary studies. PlantTribes2 can improve transcript models and then sort them, either genome-scale annotations or individual gene coding sequences, into pre-computed orthologous gene family clusters with rich functional annotation information. Then, for gene families of interest, PlantTribes2 performs downstream analyses and customizable visualizations including, (1) multiple sequence alignment, (2) gene family phylogeny, (3) estimation of synonymous and non-synonymous substitution rates among homologous sequences, and (4) inference of large-scale duplication events. We give examples of PlantTribes2 applications in functional genomic studies of economically important plant families, namely transcriptomics in the weedy Orobanchaceae and a core orthogroup analysis (CROG) in Rosaceae. PlantTribes2 is freely available for use within the main public Galaxy instance and can be downloaded from GitHub or Bioconda. Importantly, PlantTribes2 can be readily adapted for use with genomic and transcriptomic data from any kind of organism.

Characterization of the Brassica napus Flavonol Synthase Gene Family Reveals Bifunctional Flavonol Synthases

Flavonol synthase (FLS) is a key enzyme for the formation of flavonols, which are a subclass of the flavonoids. FLS catalyzes the conversion of dihydroflavonols to flavonols. The enzyme belongs to the 2-oxoglutarate-dependent dioxygenases (2-ODD) superfamily. We characterized the FLS gene family of Brassica napus that covers 13 genes, based on the genome sequence of the B. napus cultivar Express 617. The goal was to unravel which BnaFLS genes are relevant for seed flavonol accumulation in the amphidiploid species B. napus. Two BnaFLS1 homeologs were identified and shown to encode bifunctional enzymes. Both exhibit FLS activity as well as flavanone 3-hydroxylase (F3H) activity, which was demonstrated in vivo and in planta. BnaFLS1-1 and -2 are capable of converting flavanones into dihydroflavonols and further into flavonols. Analysis of spatio-temporal transcription patterns revealed similar expression profiles of BnaFLS1 genes. Both are mainly expressed in reproductive organs and co-expressed with the genes encoding early steps of flavonoid biosynthesis. Our results provide novel insights into flavonol biosynthesis in B. napus and contribute information for breeding targets with the aim to modify the flavonol content in rapeseed.

Functional Characterisation of Banana (Musa spp.) 2-Oxoglutarate-Dependent Dioxygenases Involved in Flavonoid Biosynthesis

Bananas (Musa) are non-grass, monocotyledonous, perennial plants that are well known for their edible fruits. Their cultivation provides food security and employment opportunities in many countries. Banana fruits contain high levels of minerals and phytochemicals, including flavonoids, which are beneficial for human nutrition. To broaden the knowledge on flavonoid biosynthesis in this major crop plant, we aimed to identify and functionally characterise selected structural genes encoding 2-oxoglutarate-dependent dioxygenases, involved in the formation of the flavonoid aglycon. Musa candidates genes predicted to encode flavanone 3-hydroxylase (F3H), flavonol synthase (FLS) and anthocyanidin synthase (ANS) were assayed. Enzymatic functionalities of the recombinant proteins were confirmed in vivo using bioconversion assays. Moreover, transgenic analyses in corresponding Arabidopsis thaliana mutants showed that MusaF3H, MusaFLS and MusaANS were able to complement the respective loss-of-function phenotypes, thus verifying functionality of the enzymes in planta. Knowledge gained from this work provides a new aspect for further research towards genetic engineering of flavonoid biosynthesis in banana fruits to increase their antioxidant activity and nutritional value.

Metabolite Profiling and Transcriptome Analysis Provide Insight into Seed Coat Color in Brassica juncea

The allotetraploid species Brassica juncea (mustard) is grown worldwide as oilseed and vegetable crops; the yellow seed-color trait is particularly important for oilseed crops. Here, to examine the factors affecting seed coat color, we performed a metabolic and transcriptomic analysis of yellow- and dark-seeded B. juncea seeds. In this study, we identified 236 compounds, including 31 phenolic acids, 47 flavonoids, 17 glucosinolates, 38 lipids, 69 other hydroxycinnamic acid compounds, and 34 novel unknown compounds. Of these, 36 compounds (especially epicatechin and its derivatives) accumulated significantly different levels during the development of yellow- and dark-seeded B. juncea. In addition, the transcript levels of BjuDFR, BjuANS,BjuBAN, BjuTT8, and BjuTT19 were closely associated with changes to epicatechin and its derivatives during seed development, implicating this pathway in the seed coat color determinant in B. juncea. Furthermore, we found numerous variations of sequences in the TT8A genes that may be associated with the stability of seed coat color in B. rapa, B. napus, and B. juncea, which might have undergone functional differentiation during polyploidization in the Brassica species. The results provide valuable information for understanding the accumulation of metabolites in the seed coat color of B. juncea and lay a foundation for exploring the underlying mechanism.

The report of anthocyanins in the betalain-pigmented genus Hylocereus is not well evidenced and is not a strong basis to refute the mutual exclusion paradigm

Here we respond to the paper entitled "Contribution of anthocyanin pathways to fruit flesh coloration in pitayas" (Fan et al., BMC Plant Biol 20:361, 2020). In this paper Fan et al. 2020 propose that the anthocyanins can be detected in the betalain-pigmented genus Hylocereus, and suggest they are responsible for the colouration of the fruit flesh. We are open to the idea that, given the evolutionary maintenance of fully functional anthocyanin synthesis genes in betalain-pigmented species, anthocyanin pigmentation might co-occur with betalain pigments, as yet undetected, in some species. However, in absence of the LC-MS/MS spectra and co-elution/fragmentation of the authentic standard comparison, the findings of Fan et al. 2020 are not credible. Furthermore, our close examination of the paper, and re-analysis of datasets that have been made available, indicate numerous additional problems. Namely, the failure to detect betalains in an untargeted metabolite analysis, accumulation of reported anthocyanins that does not correlate with the colour of the fruit, absence of key anthocyanin synthesis genes from qPCR data, likely mis-identification of key anthocyanin genes, unreproducible patterns of correlated RNAseq data, lack of gene expression correlation with pigmentation accumulation, and putative transcription factors that are weak candidates for transcriptional up-regulation of the anthocyanin pathway.

Expression Characterization of Flavonoid Biosynthetic Pathway Genes and Transcription Factors in Peanut Under Water Deficit Conditions

Drought is one of the hostile environmental stresses that limit the yield production of crop plants by modulating their growth and development. Peanut (Arachis hypogaea) has a wide range of adaptations to arid and semi-arid climates, but its yield is prone to loss due to drought. Other than beneficial fatty acids and micronutrients, peanut harbors various bioactive compounds including flavonoids that hold a prominent position as antioxidants in plants and protect them from oxidative stress. In this study, understanding of the biosynthesis of flavonoids in peanut under water deficit conditions was developed through expression analysis and correlational analysis and determining the accumulation pattern of phenols, flavonols, and anthocyanins. Six peanut varieties (BARD479, BARI2011, BARI2000, GOLDEN, PG1102, and PG1265) having variable responses against drought stress have been selected. Higher water retention and flavonoid accumulation have been observed in BARI2011 but downregulation has been observed in the expression of genes and transcription factors (TFs) which indicated the maintenance of normal homeostasis. ANOVA revealed that the expression of flavonoid genes and TFs is highly dependent upon the genotype of peanut in a spatiotemporal manner. Correlation analysis between expression of flavonoid biosynthetic genes and TFs indicated the role of AhMYB111 and AhMYB7 as an inhibitor for AhF3H and AhFLS, respectively, and AhMYB7, AhTTG1, and AhCSU2 as a positive regulator for the expression of Ah4CL, AhCHS, and AhF3H, respectively. However, AhbHLH and AhGL3 revealed nil-to-little relation with the expression of flavonoid biosynthetic pathway genes. Correlational analysis between the expression of TFs related to the biosynthesis of flavonoids and the accumulation of phenolics, flavonols, and anthocyanins indicated coregulation of flavonoid synthesis by TFs under water deficit conditions in peanut. This study would provide insight into the role of flavonoid biosynthetic pathway in drought response in peanut and would aid to develop drought-tolerant varieties of peanut.

Transcriptome Sequence Reveals Candidate Genes Involving in the Post-Harvest Hardening of Trifoliate Yam Dioscorea dumetorum

Storage ability of trifoliate yam (Dioscorea dumetorum) is restricted by a severe post-harvest hardening (PHH) phenomenon, which starts within the first 24 h after harvest and renders tubers inedible. Previous work has only focused on the biochemical changes affecting PHH in D. dumetorum. To the best of our knowledge, the candidate genes responsible for the hardening of D. dumetorum have not been identified. Here, transcriptome analyses of D. dumetorum tubers were performed in yam tubers of four developmental stages: 4 months after emergence (4MAE), immediately after harvest (AH), 3 days after harvest (3DAH) and 14 days after harvest (14DAH) of four accessions (Bangou 1, Bayangam 2, Fonkouankem 1, and Ibo sweet 3) using RNA-Seq. In total, between AH and 3DAH, 165, 199, 128 and 61 differentially expressed genes (DEGs) were detected in Bayangam 2, Fonkouankem 1, Bangou 1 and Ibo sweet 3, respectively. Functional analysis of DEGs revealed that genes encoding for CELLULOSE SYNTHASE A (CESA), XYLAN O-ACETYLTRANSFERASE (XOAT), CHLOROPHYLL A/B BINDING PROTEIN1, 2, 3, 4 (LHCB1, LHCB2, LHCB3, and LCH4) and an MYB transcription factor were predominantly and significantly up-regulated 3DAH, implying that these genes were potentially involved in the PHH as confirmed by qRT-PCR. A hypothetical mechanism of this phenomenon and its regulation has been proposed. These findings provide the first comprehensive insights into gene expression in yam tubers after harvest and valuable information for molecular breeding against the PHH.

The Expression of IbMYB1 Is Essential to Maintain the Purple Color of Leaf and Storage Root in Sweet Potato [Ipomoea batatas (L.) Lam]

IbMYB1 was one of the major anthocyanin biosynthesis regulatory genes that has been identified and utilized in purple-fleshed sweet potato breeding. At least three members of this gene, namely, IbMYB1-1, -2a, and -2b, have been reported. We found that IbMYB1-2a and -2b are not necessary for anthocyanin accumulation in a variety of cultivated species (hexaploid) with purple shoots or purplish rings/spots of flesh. Transcriptomic and quantitative reverse transcription PCR (RT-qPCR) analyses revealed that persistent and vigorous expression of IbMYB1 is essential to maintain the purple color of leaves and storage roots in this type of cultivated species, which did not contain IbMYB1-2 gene members. Compared with IbbHLH2, IbMYB1 is an early response gene of anthocyanin biosynthesis in sweet potato. It cannot exclude the possibility that other MYBs participate in this gene regulation networks. Twenty-two MYB-like genes were identified from 156 MYBs to be highly positively or negatively correlated with the anthocyanin content in leaves or flesh. Even so, the IbMYB1 was most coordinately expressed with anthocyanin biosynthesis genes. Differences in flanking and coding sequences confirm that IbMYB2s, the highest similarity genes of IbMYB1, are not the members of IbMYB1. This phenomenon indicates that there may be more members of IbMYB1 in sweet potato, and the genetic complementation of these members is involved in the regulation of anthocyanin biosynthesis. The 3' flanking sequence of IbMYB1-1 is homologous to the retrotransposon sequence of TNT1-94. Transposon movement is involved in the formation of multiple members of IbMYB1. This study provides critical insights into the expression patterns of IbMYB1, which are involved in the regulation of anthocyanin biosynthesis in the leaf and storage root. Notably, our study also emphasized the presence of a multiple member of IbMYB1 for genetic improvement.