Computational and Transcriptomic Analysis Unraveled OsMATE34 as a Putative Anthocyanin Transporter in Black Rice (Oryza sativa L.) Caryopsis

Transcriptional convergence after repeated duplication of an amino acid transporter gene leads to the independent emergence of the black husk/pericarp trait in barley and rice

The repeated emergence of the same trait (convergent evolution) in distinct species is an interesting phenomenon and manifests visibly the power of natural selection. The underlying genetic mechanisms have important implications to understand how the genome evolves under environmental challenges. In cereal crops, both rice and barley can develop black-coloured husk/pericarp due to melanin accumulation. However, it is unclear if this trait shares a common origin. Here, we fine-mapped the barley HvBlp gene controlling the black husk/pericarp trait and confirmed its function by gene silencing. The result was further supported by a yellow husk/pericarp mutant with deletion of the HvBlp gene, derived from gamma ray radiation of the wild-type W1. HvBlp encodes a putative tyrosine transporter homologous to the black husk gene OsBh4 in rice. Surprisingly, synteny and phylogenetic analyses showed that HvBlp and OsBh4 belonged to different lineages resulted from dispersed and tandem duplications, respectively, suggesting that the black husk/pericarp trait has emerged independently. The dispersed duplication (dated at 21.23 MYA) yielding HvBlp occurred exclusively in the common ancestor of Triticeae. HvBlp and OsBh4 displayed converged transcription in husk/pericarp tissues, contributing to the black husk/pericarp trait. Further transcriptome and metabolome data identified critical candidate genes and metabolites related to melanin production in barley. Taken together, our study described a compelling case of convergent evolution resulted from transcriptional convergence after repeated gene duplication, providing valuable genetic insights into phenotypic evolution. The identification of the black husk/pericarp genes in barley also has great potential in breeding for stress-resilient varieties with higher nutritional values.

OsGSTU34, a Bz2-like anthocyanin-related glutathione transferase transporter, is essential for rice (Oryza sativa L.) organs coloration

Anthocyanins are a flavonoid compound known as one of the most important chromogenic substances. They play several functions, including health promotion and sustaining plants during adverse conditions. They are synthesized at the endoplasmic reticulum and sequestered in the vacuole. In this work, we generated knock-out lines of OsGSTU34, a glutathione transporter's tau gene family, with no transgene line and off-target through CRISPR/Cas9 mutagenesis and highlighted the loss of pigmentation in rice flowers, leaves, stems, shoots, and caryopsis. The anthocyanin quantification in the wild-type BLWT and mutant line BLG34-8 caryopsis showed that cyanidin-3-O-glucoside (C3G) and peonidin-3-O-glucoside (P3G) were almost undetectable in the mutant line. A tandem mass tag (TMT) labeling proteomic analysis was conducted to elucidate the proteomic changes in the BLWT and BLG34-8. The result revealed that 1175 proteins were altered, including 408 that were down-regulated and 767 that were upregulated. The accumulation of the OsGSTU34-related protein (Q8L576), along with several anthocyanin-related proteins, was down-regulated. The enrichment analysis showed that the down-regulated proteins were enriched in different pathways, among which the phenylpropanoid biosynthesis pathway, flavonoid biosynthesis metabolites, and anthocyanin biosynthesis pathway. Protein interaction network prediction revealed that glutathione-S-transferase (Q8L576) was connected to the proteins involved in the flavonoid and anthocyanin biosynthesis pathways, such as flavanone 3-dioxygenase 1 (Q7XM21), leucoanthocyanidin dioxygenase 1 (Q93VC3), 4-coumarate-CoA ligase 2 (Q42982), phenylalanine ammonia-lyase (P14717), chalcone synthase 1 (Q2R3A1), and 4-coumarate-CoA ligase 5 (Q6ZAC1). However, the expression of the most important anthocyanin biosynthesis gene was not altered, suggesting that only the transport mechanism was affected. Our findings highlight new insight into the anthocyanin pigmentation in black rice and provide new perspectives for future research.

Integrative HPLC profiling and transcriptome analysis revealed insights into anthocyanin accumulation and key genes at three developmental stages of black rice (Oryza sativa. L) caryopsis

Introduction

Anthocyanins are plants' secondary metabolites belonging to the flavonoid class with potential health-promoting properties. They are greatly employed in the food industry as natural alternative food colorants for dairy and ready-to-eat desserts and pH indicators. These tremendous advantages make them economically important with increasing market trends. Black rice is a rich source of anthocyanin that can be used to ensure food and nutritional security around the world. However, research on anthocyanin accumulation and gene expression during rice caryopsis development is lacking.

Methods

In this study, we combined high-performance liquid chromatography (HPLC) and transcriptome analysis to profile the changes in anthocyanin content and gene expression dynamics at three developmental stages (milky, doughy, and mature).

Results

Our results showed that anthocyanin accumulation started to be visible seven days after flowering (DAF), increased rapidly from milky (11 DAF) to dough stage, then started decreasing after the peak was attained at 18 DAF. RNA-seq showed that 519 out of 14889, 477 out of 17914, and 1614 out of 18810 genes were uniquely expressed in the milky, doughy, and mature stages, respectively. We performed three pairwise comparisons: milky vs. dough, milky vs. mature, and dough vs. mature, and identified 6753, 9540, and 2531 DEGs, respectively. The DEGs' abundance was higher in milky vs. mature, with 5527 up-regulated genes and 4013 down-regulated genes, while it was smaller in dough vs. mature, with 1419 up-regulated genes and 1112 down-regulated DEGs. This result was consistent with the changes in anthocyanin profiling, and the expression of structural, regulatory, and transporter genes involved in anthocyanin biosynthesis showed their highest expression at the dough stage. Through the gene expression profile and protein interaction network, we deciphered six main contributors of the anthocyanin peak observed at dough stage, including OsANS, OsDFR, OsGSTU34, OsMYB3, OsbHLH015, and OsWD40-50.

Discussion

This study is the first to report the investigation of anthocyanin and gene expression at three developmental stages of black rice caryopsis. The findings of this study could aid in predicting the best harvesting time to achieve maximum anthocyanin content and the best time to collect samples for various gene expression analysis, laying the groundwork for future research into the molecular mechanisms underlying rice caryopsis coloration.

Genome-Wide Association Study Reveals the Genetic Basis of Total Flavonoid Content in Brown Rice

Flavonoids have anti-inflammatory, antioxidative, and anticarcinogenic effects. Breeding rice varieties rich in flavonoids can prevent chronic diseases such as cancer and cardio-cerebrovascular diseases. However, most of the genes reported are known to regulate flavonoid content in leaves or seedlings. To further elucidate the genetic basis of flavonoid content in rice grains and identify germplasm rich in flavonoids in grains, a set of rice core collections containing 633 accessions from 32 countries was used to determine total flavonoid content (TFC) in brown rice. We identified ten excellent germplasms with TFC exceeding 300 mg/100 g. Using a compressed mixed linear model, a total of 53 quantitative trait loci (QTLs) were detected through a genome-wide association study (GWAS). By combining linkage disequilibrium (LD) analysis, location of significant single nucleotide polymorphisms (SNPs), gene expression, and haplotype analysis, eight candidate genes were identified from two important QTLs (qTFC1-6 and qTFC9-7), among which LOC_Os01g59440 and LOC_Os09g24260 are the most likely candidate genes. We also analyzed the geographic distribution and breeding utilization of favorable haplotypes of the two genes. Our findings provide insights into the genetic basis of TFC in brown rice and could facilitate the breeding of flavonoid-rich varieties, which may be a prevention and adjuvant treatment for cancer and cardio-cerebrovascular diseases.

HvGST plays a key role in anthocyanin accumulation in colored barley

Blue aleurone of barley is caused by the accumulation of delphinidin-based derivatives. Although these compounds are ideal nutrients for human health, they are undesirable contaminants in malt brewing. Therefore, the ability to add and remove this trait easily would facilitate breeding barley for different purposes. Here we identified a glutathione S-transferase gene (HvGST) that was responsible for the blue aleurone trait in Tibetan qingke barley by performing a genome-wide association study and RNA-sequencing analysis. Gene variation and expression analysis indicated that HvGST also participates in the transport and accumulation of anthocyanin in purple barley. Haplotype and the geographic distribution analyses of HvGST alleles revealed two independent natural variants responsible for the emergence of white aleurone: a 203-bp deletion causing premature termination of translation in qingke barley and two key single nucleotide polymorphisms in the promoter resulting in low transcription in Western barley. This study contributes to a better understanding of mechanisms of colored barley formation, and provides a comprehensive reference for marker-assisted barley breeding.

Rice co-expression network analysis identifies gene modules associated with agronomic traits

Identifying trait-associated genes is critical for rice (Oryza sativa) improvement, which usually relies on map-based cloning, quantitative trait locus analysis, or genome-wide association studies. Here we show that trait-associated genes tend to form modules within rice gene co-expression networks, a feature that can be exploited to discover additional trait-associated genes using reverse genetics. We constructed a rice gene co-expression network based on the graphical Gaussian model using 8,456 RNA-seq transcriptomes, which assembled into 1,286 gene co-expression modules functioning in diverse pathways. A number of the modules were enriched with genes associated with agronomic traits, such as grain size, grain number, tiller number, grain quality, leaf angle, stem strength, and anthocyanin content, and these modules are considered to be trait-associated gene modules. These trait-associated gene modules can be used to dissect the genetic basis of rice agronomic traits and to facilitate the identification of trait genes. As an example, we identified a candidate gene, OCTOPUS-LIKE 1 (OsOPL1), a homolog of the Arabidopsis (Arabidopsis thaliana) OCTOPUS gene, from a grain size module and verified it as a regulator of grain size via functional studies. Thus, our network represents a valuable resource for studying trait-associated genes in rice.

RNA-seq reveals transcriptional differences in anthocyanin and vitamin biosynthetic pathways between black and white rice

Anthocyanins and vitamins in black rice are the micronutrients vital to human health, both of which predominantly accumulate in the bran fraction. Some studies have demonstrated that black rice contains more vitamins compared with common white rice, indicating potential association between anthocyanin and vitamin accumulation. In this study, transcriptomes of pericarps collected from 27 black rice accessions and 49 white rice accessions at 10 days after flowering (DAF) were sequenced and analyzed. We identified 830 differentially expressed genes (DEGs) including 58 transcription factors (TFs) between black and white rice. Among 58 differentially expressed transcription factors, OsTTG1 was confirmed to be the one and only WD40 repeat protein regulating anthocyanin biosynthesis in the pericarp. Moreover, we identified 53 differentially expressed synthetic-related genes among 42 main synthesis enzymes in the biosynthesis pathway of seven vitamins including β-carotene, vitamin B1, vitamin B2, vitamin B5, vitamin B7, vitamin B9 and vitamin E. Collectively, our results provide valuable insights into the molecular mechanism of biosynthesis of anthocyanins and vitamins and the potential effect of anthocyanin biosynthesis on vitamin biosynthesis in black rice.

Combined Analysis of BSA-Seq Based Mapping, RNA-Seq, and Metabolomic Unraveled Candidate Genes Associated with Panicle Grain Number in Rice (Oryza sativa L.)

Rice grain yield is a complex and highly variable quantitative trait consisting of several key components, including the grain weight, the effective panicles per unit area, and the grain number per panicle (GNPP). The GNPP is a significant contributor to grain yield controlled by multiple genes (QTL) and is crucial for improvement. Attempts have been made to find genes for this trait, which has always been a challenging and arduous task through conventional methods. We combined a BSA analysis, RNA profiling, and a metabolome analysis in the present study to identify new candidate genes involved in the GNPP. The F₂ population from crossing R4233 (high GNPP) and Ce679 (low GNPP) revealed a frequency distribution fitting two segregated genes. Three pools, including low, middle, and high GNPP, were constructed and a BSA analysis revealed six candidate regions spanning 5.38 Mb, containing 739 annotated genes. Further, a conjunctive analysis of BSA-Seq and RNA-Seq showed 31 differentially expressed genes (DEGs) in the candidate intervals. Subsequently, a metabolome analysis showed 1024 metabolites, with 71 significantly enriched, including 44 up and 27 downregulated in Ce679 vs. R4233. A KEGG enrichment analysis of these 31 DEGs and 71 differentially enriched metabolites (DEMs) showed two genes, Os12g0102100 and Os01g0580500, significantly enriched in the metabolic pathways’ biosynthesis of secondary metabolites, cysteine and methionine metabolism, and fatty acid biosynthesis. Os12g0102100, which encodes for the alcohol dehydrogenase superfamily and a zinc-containing protein, is a novel gene whose contribution to the GNPP is not yet elucidated. This gene coding for mitochondrial trans-2-enoyl-CoA reductase is involved in the biosynthesis of myristic acid, also known as tetradecanoic acid. The Os01g0580500 coding for the enzyme 1-aminoclopropane-1-carboxylate oxidase (OsACO7) is responsible for the final step of the ethylene biosynthesis pathway through the conversion of 1-aminocyclopropane-1-carboxylic acid (ACC) into ethylene. Unlike Os12g0102100, this gene was significantly upregulated in R4233, downregulated in Ce679, and significantly enriched in two of the three metabolite pathways. This result pointed out that these two genes are responsible for the difference in the GNPP in the two cultivars, which has never been identified. Further validation studies may disclose the physiological mechanisms through which they regulate the GNPP in rice.

Developing Germplasm and Promoting Consumption of Anthocyanin-Rich Grains for Health Benefits

Malnutrition, unhealthy diets, and lifestyle changes are the major risk factors for overweight and obesity-linked chronic diseases in humans adversely impact achieving sustainable development goals. Colored grains are a source of anthocyanins, a group of flavonoids, that contribute positively to human health. This review focuses on genetic variation harnessed through breeding and biotechnology tools for developing anthocyanin-rich grain crops. Agronomic practices, genotype × environment interactions, different stresses, seed development and seed maturity are factors that impact the content and composition of anthocyanins. Significant progress has been made in characterizing genes associated with anthocyanin biosynthesis in cereal and other crops. Breeding has led to the development and release of grain anthocyanin-rich crop cultivars in Europe, America and in some countries in Asia. Notably, genetic engineering utilizing specific transcription factors and gene editing has led to the development of anthocyanin-rich genetic variants without any significant yield penalty. A variety of food products derived from colored grains or flours are now available in grocery stores and supermarkets worldwide. The public perception about anthocyanin-rich food is positive, but availability, affordability, and willingness to pay a higher price than before limit consumption. Together with other seed nutrition traits in breeding programs the inclusion of anthocyanins can ensure the development of cultivars that meet nutrition needs of humans, especially in the developing world.

Integrative Analysis of Metabolome and Transcriptome Reveals the Mechanism of Color Formation in Liriope spicata Fruit

Liriope spicata is an important ornamental ground cover plant, with a fruit color that turns from green to black during the development and ripening stages. However, the material basis and regulatory mechanism of the color variation remains unclear. In this study, a total of 31 anthocyanins and 2 flavonols were identified from the skin of L. spicata fruit via integrative analysis on the metabolome and transcriptome of three developmental stages. The pigments of black/mature fruits are composed of five common anthocyanin compounds, of which Peonidin 3–O–rutinoside and Delphinidin 3–O–glucoside are the most differential metabolites for color conversion. Using dual-omics joint analysis, the mechanism of color formation was obtained as follows. The expression of structural genes including 4CL, F3H, F3′H, F3′5′H and UFGT were activated due to the upregulation of transcription factor genes MYB and bHLH. As a result, a large amount of precursor substances for the synthesis of flavonoids accumulated. After glycosylation, stable pigments were generated which promoted the accumulation of anthocyanins and the formation of black skin.

Plant Secondary Metabolite Transporters: Diversity, Functionality, and Their Modulation

Secondary metabolites (SMs) play crucial roles in the vital functioning of plants such as growth, development, defense, and survival via their transportation and accumulation at the required site. However, unlike primary metabolites, the transport mechanisms of SMs are not yet well explored. There exists a huge gap between the abundant presence of SM transporters, their identification, and functional characterization. A better understanding of plant SM transporters will surely be a step forward to fulfill the steeply increasing demand for bioactive compounds for the formulation of herbal medicines. Thus, the engineering of transporters by modulating their expression is emerging as the most viable option to achieve the long-term goal of systemic metabolic engineering for enhanced metabolite production at minimum cost. In this review article, we are updating the understanding of recent advancements in the field of plant SM transporters, particularly those discovered in the past two decades. Herein, we provide notable insights about various types of fully or partially characterized transporters from the ABC, MATE, PUP, and NPF families including their diverse functionalities, structural information, potential approaches for their identification and characterization, several regulatory parameters, and their modulation. A novel perspective to the concept of "Transporter Engineering" has also been unveiled by highlighting its potential applications particularly in plant stress (biotic and abiotic) tolerance, SM accumulation, and removal of anti-nutritional compounds, which will be of great value for the crop improvement program. The present study creates a roadmap for easy identification and a better understanding of various transporters, which can be utilized as suitable targets for transporter engineering in future research.