Differential expression of three flavanone 3-hydroxylase genes in grains and coleoptiles of wheat

Genome–Transcriptome Transition Approaches to Characterize Anthocyanin Biosynthesis Pathway Genes in Blue, Black and Purple Wheat

Colored wheat has gained enormous attention from the scientific community, but the information available on the anthocyanin biosynthetic genes is very minimal. The study involved their genome-wide identification, in silico characterization and differential expression analysis among purple, blue, black and white wheat lines. The recently released wheat genome mining putatively identified eight structural genes in the anthocyanin biosynthesis pathway with a total of 1194 isoforms. Genes showed distinct exon architecture, domain profile, regulatory elements, chromosome emplacement, tissue localization, phylogeny and synteny, indicative of their unique function. RNA sequencing of developing seeds from colored (black, blue and purple) and white wheats identified differential expressions in 97 isoforms. The F3H on group two chromosomes and F3′5′H on 1D chromosomes could be significant influencers in purple and blue color development, respectively. Apart from a role in anthocyanin biosynthesis, these putative structural genes also played an important role in light, drought, low temperature and other defense responses. The information can assist in targeted anthocyanin production in the wheat seed endosperm.

Proanthocyanidin biosynthesis in the developing wheat seed coat investigated by chemical and RNA-Seq analysis

The composition of proanthocyanidins in the testa (seed coat) of bread wheat was analyzed by thiolysis of PA oligomers from developing grain and found to consist of (+)-catechin monomers, with a small amount of (+)-gallocatechin. The average chain length of soluble PA stayed relatively constant between 10 and 20 days post-anthesis, whereas that of unextractable PA increased over the same period, suggesting that increases in chain length might account for the insolubility of PAs from mature wheat grain. We carried out RNA-Seq followed by differential expression analysis from dissected tissues of developing grain from red- and white-grained near-isogenic lines differing in the presence of an active R gene that encodes a MYB transcription factor involved in control of PA biosynthesis. In addition to genes already identified encoding chalcone synthase, chalcone isomerase, flavanone 3-hydroxylase, and dihydroxyflavonoid 4-reductase, we showed that wheat genes encoding phenylalanine ammonia lyase, flavonoid 3',5'-hydroxylase, leucoanthocyanidin reductase, and a glutathione S-transferase (the orthologue of maize Bronze-2) were more highly expressed in the red NIL. We also identified candidate orthologues of other catalytic and regulatory components of flavonoid biosynthesis in wheat.

Flavonoid Biosynthesis Genes in Triticum aestivum L.: Methylation Patterns in Cis-Regulatory Regions of the Duplicated CHI and F3H Genes

Flavonoids are a diverse group of secondary plant metabolites that play an important role in the regulation of plant development and protection against stressors. The biosynthesis of flavonoids occurs through the activity of several enzymes, including chalcone isomerase (CHI) and flavanone 3-hydroxylase (F3H). A functional divergence between some copies of the structural TaCHI and TaF3H genes was previously shown in the allohexaploid bread wheat Triticum aestivum L. (BBAADD genome). We hypothesized that the specific nature of TaCHI and TaF3H expression may be induced by the methylation of the promoter. It was found that the predicted position of CpG islands in the promoter regions of the analyzed genes and the actual location of methylation sites did not match. We found for the first time that differences in the methylation status could affect the expression of TaCHI copies, but not the expression of TaF3Hs. At the same time, we revealed significant differences in the structure of the promoters of only the TaF3H genes, while the TaCHI promoters were highly homologous. We assume that the promoter structure in TaF3Hs primarily affects the change in the nature of gene expression. The data obtained are important for understanding the mechanisms that regulate the synthesis of flavonoids in allopolyploid wheat and show that differences in the structure of promoters have a key effect on gene expression.

Genome-Wide Association Study of Pericarp Color in Rice Using Different Germplasm and Phenotyping Methods Reveals Different Genetic Architectures

Pericarp colors (PC) in rice are determined by the types and content of flavonoids in the pericarp. The flavonoid compounds have strong antioxidant activities and are beneficial to human health. However, the genetic basis of PC in rice is still not well-understood. In this study, a genome-wide association study (GWAS) of PC was performed in a diverse rice collection consisting of 442 accessions using different phenotyping methods in two locations over 2 years. In the whole population consisting of white and colored pericarp rice, a total of 11 quantitative trait loci (QTLs) were identified using two phenotyping methods. Among these QTLs, nine were identified using the phenotypes represented by the presence and absence of pigmentation in pericarp, while 10 were identified using phenotypes of the degree of PC (DPC), in which eight are common QTLs identified using the two phenotyping methods. Using colored rice accessions and phenotypes based on DPC, four QTLs were identified, and they were totally different from the QTLs identified using the whole population, suggesting the masking effects of major genes on minor genes. Compared with the previous studies, 10 out of the 15 QTLs are first reported in this study. Based on the differential expression analysis of the predicted genes within the QTL region by both RNA-seq and real-time PCR (RT-PCR) and the gene functions in previous studies, LOC_Os01g49830, encoding a RAV transcription factor was considered as the candidate gene underlying qPC-1, a novel QTL with a large effect in this study. Our results provide a new insight into the genetic basis of PC in rice and contribute to developing the value-added rice with optimized flavonoid content through molecular breeding.

Phenotypic Evaluation and Genetic Analysis of Seedling Emergence in a Global Collection of Wheat Genotypes (Triticum aestivum L.) Under Limited Water Availability

The challenge in establishing an early-sown wheat crop in southern Australia is the need for consistently high seedling emergence when sowing deep in subsoil moisture (>10 cm) or into dry top-soil (4 cm). However, the latter is strongly reliant on a minimum soil water availability to ensure successful seedling emergence. This study aimed to: (1) evaluate 233 Australian and selected international wheat genotypes for consistently high seedling emergence under limited soil water availability when sown in 4 cm of top-soil in field and glasshouse (GH) studies; (2) ascertain genetic loci associated with phenotypic variation using a genome-wide association study (GWAS); and (3) compare across loci for traits controlling coleoptile characteristics, germination, dormancy, and pre-harvest sprouting. Despite significant (P < 0.001) environment and genotype-by-environment interactions within and between field and GH experiments, eight genotypes that included five cultivars, two landraces, and one inbred line had consistently high seedling emergence (mean value > 85%) across nine environments. Moreover, 21 environment-specific quantitative trait loci (QTL) were detected in GWAS analysis on chromosomes 1B, 1D, 2B, 3A, 3B, 4A, 4B, 5B, 5D, and 7D, indicating complex genetic inheritance controlling seedling emergence. We aligned QTL for known traits and individual genes onto the reference genome of wheat and identified 16 QTL for seedling emergence in linkage disequilibrium with coleoptile length, width, and cross-sectional area, pre-harvest sprouting and dormancy, germination, seed longevity, and anthocyanin development. Therefore, it appears that seedling emergence is controlled by multifaceted networks of interrelated genes and traits regulated by different environmental cues.

Population genomics and haplotype analysis in spelt and bread wheat identifies a gene regulating glume color

The cloning of agriculturally important genes is often complicated by haplotype variation across crop cultivars. Access to pan-genome information greatly facilitates the assessment of structural variations and rapid candidate gene identification. Here, we identified the red glume 1 (Rg-B1) gene using association genetics and haplotype analyses in ten reference grade wheat genomes. Glume color is an important trait to characterize wheat cultivars. Red glumes are frequent among Central European spelt, a dominant wheat subspecies in Europe before the 20th century. We used genotyping-by-sequencing to characterize a global diversity panel of 267 spelt accessions, which provided evidence for two independent introductions of spelt into Europe. A single region at the Rg-B1 locus on chromosome 1BS was associated with glume color in the diversity panel. Haplotype comparisons across ten high-quality wheat genomes revealed a MYB transcription factor as candidate gene. We found extensive haplotype variation across the ten cultivars, with a particular group of MYB alleles that was conserved in red glume wheat cultivars. Genetic mapping and transient infiltration experiments allowed us to validate this particular MYB transcription factor variants. Our study demonstrates the value of multiple high-quality genomes to rapidly resolve copy number and haplotype variations in regions controlling agriculturally important traits.

Validation of barley 2OGO gene as a functional orthologue of Arabidopsis DMR6 gene in Fusarium head blight susceptibility

Fusarium head blight (FHB) caused by Fusarium graminearum (Fg) is a devastating disease of crops, especially wheat and barley, resulting in significant yield loss and reduced grain quality. Fg infection leads to the production of mycotoxins, whose consumption is toxic to humans and livestock. The Arabidopsis DMR6 gene encodes a putative 2-oxoglutarate Fe(II)-dependent oxygenase (2OGO) and has been identified as a susceptibility factor to downy mildew. We generated site-specific mutations in Arabidopsis At2OGO by CRISPR/Cas9 gene editing. The resulting At2OGO knock-out (KO) mutants display enhanced resistance to Fg in a detached inflorescence infection assay. Expression profiling of defense genes revealed that impairment of At2OGO function resulted in the upregulation of defense genes that are regulated by salicylic acid (SA), jasmonic acid (JA) and ethylene (ET) pathways. Complementation of the At2OGO-KO lines with a barley (cv. Conlon) orthologue, Hv2OGO, restored susceptibility to Fg. This result indicates that the Hv2OGO gene is functionally equivalent to its Arabidopsis counterpart and, hence, may have a similar role in conditioning susceptibility to FHB in barley. These results provide a molecular basis for proposing 2OGO as a plant immunity suppressor in Arabidopsis and potentially in barley plants and establish a rationale and strategy for enhancing FHB resistance in barley.

Evolutionary and functional analyses of the 2-oxoglutarate-dependent dioxygenase genes involved in the flavonoid biosynthesis pathway in tobacco

MAIN CONCLUSION

This study illustrates the differences in the gene structure of 2-oxoglutarate-dependent oxygenase involved in flavonoid biosynthesis (2ODD-IFB), and their potential roles in regulating tobacco flavonoid biosynthesis and plant growth. Flavonol synthase (FLS), anthocyanidin synthase (ANS), and flavanone 3β-hydroxylase belong to the 2-oxoglutarate-dependent (2ODD) oxygenase family, and each performs crucial functions in the biosynthesis of flavonoids. We identified two NtFLS genes, two NtANS genes, and four NtF3H genes from Nicotiana tabacum genome, as well as their homologous genes in the N. sylvestris and N. tomentosiformis genomes. Our phylogenetic analysis indicated that these three types of genes split from each other before the divergence of gymnosperms and angiosperms. FLS evolved faster in the eudicot plants, whereas ANS evolved faster in the monocot plants. Gene structure analysis revealed two fragment insertions occurred at different times in the intron one position of tobacco FLS genes. Homologous protein modeling revealed distinct structures in the N terminus of the tobacco 2ODD oxygenases. We found that the expression patterns of genes encoding tobacco 2ODD oxygenases in flavonoids biosynthesis (2ODD-IFB) did not determine the accumulation patterns of flavonoids among various tobacco tissues, but strongly affected the concentration of flavonoids in the tissues, where they were biosynthesized. More carbon resource flowed to the flavonol biosynthesis when NtANS gene was silenced, otherwise more anthocyanidin accumulated when NtFLS gene was repressed. This study illustrates the 2ODD-IFB gene structure evolution, differences among their protein structures, and provides a foundation for regulating plant development and altering flavonoid content and/or composition through the manipulation of plant 2ODD-IFB genes.

Functional Analysis of Two Flavanone-3-Hydroxylase Genes from Camellia sinensis: A Critical Role in Flavonoid Accumulation

Flavonoids are major secondary metabolites in Camellia sinensis. Flavanone-3-hydroxylase (F3H) is a key enzyme in flavonoid biosynthesis in plants. However, its role in the flavonoid metabolism in C. sinensis has not been well studied. In this study, we cloned two F3Hs from C. sinensis, named CsF3Ha and CsF3Hb, where CsF3Ha containing 1107 bases encoded 368 amino acids, and CsF3Hb containing 1071 bases encoded 357 amino acids. Enzymatic activity analysis showed both recombinant CsF3H enzymes in Escherichia coli could convert naringenin and eriodictyol into dihydrokaempferol (DHK) and dihydroquercetin (DHQ), respectively. The expression profiles showed that CsF3Ha and CsF3Hb were highly expressed in the tender leaves of tea plants. Under different abiotic stresses, the two CsF3Hs were induced remarkably by ultraviolet (UV) radiation, sucrose, and abscisic acid (ABA). In the seeds of CsF3Hs transgenic Arabidopsis thaliana, the concentration of most flavonol glycosides and oligomeric proanthocyanidins increased significantly, while the content of monocatechin derivatives decreased. The present study revealed that CsF3Hs played critical roles in flavonoid biosynthesis in tea plants.

Exploiting Phenylpropanoid Derivatives to Enhance the Nutraceutical Values of Cereals and Legumes

Phenylpropanoids are a diverse chemical class with immense health benefits that are biosynthesized from the aromatic amino acid L-phenylalanine. This article reviews the progress for accessing variation in phenylpropanoids in germplasm collections, the genetic and molecular basis of phenylpropanoid biosynthesis, and the development of cultivars dense in seed-phenylpropanoids. Progress is also reviewed on high-throughput assays, factors that influence phenylpropanoids, the site of phenylpropanoids accumulation in seed, Genotype × Environment interactions, and on consumer attitudes for the acceptance of staple foods rich in phenylpropanoids. A paradigm shift was noted in barley, maize, rice, sorghum, soybean, and wheat, wherein cultivars rich in phenylpropanoids are grown in Europe and North and Central America. Studies have highlighted some biological constraints that need to be addressed for development of high-yielding cultivars that are rich in phenylpropanoids. Genomics-assisted breeding is expected to facilitate rapid introgression into improved genetic backgrounds by minimizing linkage drag. More research is needed to systematically characterize germplasm pools for assessing variation to support crop genetic enhancement, and assess consumer attitudes to foods rich in phenylpropanoids.

Transcriptome Analysis of Differentially Expressed Genes Involved in Proanthocyanidin Accumulation in the Rhizomes of Fagopyrum dibotrys and an Irradiation-Induced Mutant

The rhizome of Fagopyrum dibotrys is a traditional Chinese medicine that has recently gained attention due to substantial findings regarding its bioactive proanthocyanidin (PA) compounds. However, the molecular mechanism underlying PA accumulation in F. dibotrys remains elusive. We previously obtained an irradiation-induced mutant (RM_R) of F. dibotrys that had a higher PA content compared to that of the wild-type (CK_R). The present study aimed to elucidate the molecular mechanism underlying PA accumulation in F. dibotrys by comparing the rhizome transcriptomes of the irradiation-induced mutant and wild-type using RNA-seq analysis. A total of 53,540 unigenes were obtained, of which 29,901 (55.84%) were annotated based on BLAST searches against public databases, and 501 unique sequences were differentially expressed between the two samples, which consisted of 204 up-regulated and 297 down-regulated unigenes. Further analysis showed that the expression patterns of some unigenes encoding enzymes involved in PAs biosynthesis in F. dibotrys rhizomes differed between RM_R and CK_R. In addition, we identified transcription factor families and several cytochrome P450s that may be involved in PA regulation in F. dibotrys. Finally, 12 unigenes that encode PA biosynthetic enzymes were confirmed by qRT-PCR analysis. This study sheds light on the molecular mechanism underlying radiation-mediated flavonoid accumulation and regulation in F. dibotrys rhizomes. These results will also provide a platform for further functional genomic research on this particular species.

Population Genetics and Structure of a Global Foxtail Millet Germplasm Collection

Foxtail millet [Setaria italica (L.) P. Beauv.] is one among the most ancient crops of dryland agriculture. It is the second most important crop among millets grown for grains or forage. Foxtail millet germplasm resources provide reservoirs of novel alleles and genes for crop improvement that have remained mostly unexplored. We genotyped a set of 190 foxtail millet germplasm accessions (including 155 accessions of the foxtail millet core collection) using genotyping-by-sequencing (GBS) for rapid single nucleotide polymorphisms (SNP) characterization to study population genetics and structure, which enable allele mining through association mapping approaches. After filtering a total 350,000 raw SNPs identified across 190 germplasm accessions for minor allele frequency (MAF), coverage for samples and coverage for sites, we retained 181 accessions with 17,714 high-quality SNPs with ≥5% MAF. Genetic structure analyses revealed that foxtail millet germplasm accessions are structured along both on the basis of races and geographic origin, and the maximum proportion of variation was due to among individuals within populations. Accessions of race indica were less diverse and are highly differentiated from those of maxima and moharia. Genome-wide linkage disequilibrium (LD) analysis showed on an average LD extends up to ∼150 kbp and varied with individual chromosomes. The utility of the data for performing genome-wide association studies (GWASs) was tested with plant pigmentation and days to flowering and identified significant marker-trait associations. This SNP data provides a foundation for exploration of foxtail millet diversity and for mining novel alleles and mapping genes for economically important traits.

Barley Ant17, encoding flavanone 3-hydroxylase (F3H), is a promising target locus for attaining anthocyanin/proanthocyanidin-free plants without pleiotropic reduction of grain dormancy

Preharvest sprouting is a serious problem in grain crop production because it causes quality deterioration and economic losses. It is well known that grain colour is closely associated with grain dormancy in wheat; white-grained lines without accumulating proanthocyanidins in testa tend to be more susceptible to preharvest sprouting than red ones. All available white-grained wheat lines are restricted to triple recessive mutations at the R loci (R-A1, R-B1, and R-D1), but barley is known to have 11 independent loci conferring the proanthocyanidin-free grain phenotype. In this study, we evaluated the dormancy levels of anthocyanin/proanthocyanidin-free ant17 mutants. Three ant17 mutants showed the same levels of dormancy as their respective wild types. Sequencing of three independent ant17 alleles detected a point mutation within the coding regions of flavanone-3-hydroxylase (F3H), which are predicted to cause a premature stop codon at different sites. The F3H locus completely cosegregated with the Ant17 position on the chromosome arm 2HL. Expression of the barley F3H gene was observed in pigmented tissues, but not in nonpigmented roots and stems. This result indicates that wheat F3H may be a promising new target locus for breeding white-grained lines with a practical level of preharvest sprouting resistance.

Gene-to-metabolite network for biosynthesis of lignans in MeJA-elicited Isatis indigotica hairy root cultures

Root and leaf tissue of Isatis indigotica shows notable anti-viral efficacy, and are widely used as "Banlangen" and "Daqingye" in traditional Chinese medicine. The plants' pharmacological activity is attributed to phenylpropanoids, especially a group of lignan metabolites. However, the biosynthesis of lignans in I. indigotica remains opaque. This study describes the discovery and analysis of biosynthetic genes and AP2/ERF-type transcription factors involved in lignan biosynthesis in I. indigotica. MeJA treatment revealed differential expression of three genes involved in phenylpropanoid backbone biosynthesis (IiPAL, IiC4H, Ii4CL), five genes involved in lignan biosynthesis (IiCAD, IiC3H, IiCCR, IiDIR, and IiPLR), and 112 putative AP2/ERF transcription factors. In addition, four intermediates of lariciresinol biosynthesis were found to be induced. Based on these results, a canonical correlation analysis using Pearson's correlation coefficient was performed to construct gene-to-metabolite networks and identify putative key genes and rate-limiting reactions in lignan biosynthesis. Over-expression of IiC3H, identified as a key pathway gene, was used for metabolic engineering of I. indigotica hairy roots, and resulted in an increase in lariciresinol production. These findings illustrate the utility of canonical correlation analysis for the discovery and metabolic engineering of key metabolic genes in plants.

Comparative transcriptional profiling of two wheat genotypes, with contrasting levels of minerals in grains, shows expression differences during grain filling

Wheat is one of the most important cereal crops in the world. To identify the candidate genes for mineral accumulation, it is important to examine differential transcriptome between wheat genotypes, with contrasting levels of minerals in grains. A transcriptional comparison of developing grains was carried out between two wheat genotypes- Triticum aestivum Cv. WL711 (low grain mineral), and T. aestivum L. IITR26 (high grain mineral), using Affymetrix GeneChip Wheat Genome Array. The study identified a total of 580 probe sets as differentially expressed (with log2 fold change of ≥2 at p≤0.01) between the two genotypes, during grain filling. Transcripts with significant differences in induction or repression between the two genotypes included genes related to metal homeostasis, metal tolerance, lignin and flavonoid biosynthesis, amino acid and protein transport, vacuolar-sorting receptor, aquaporins, and stress responses. Meta-analysis revealed spatial and temporal signatures of a majority of the differentially regulated transcripts.

Transcriptome sequencing of a chimaera reveals coordinated expression of anthocyanin biosynthetic genes mediating yellow formation in herbaceous peony (Paeonia lactiflora Pall.)

BACKGROUND

Herbaceous peony (Paeonia lactiflora Pall.) is a traditional flower in China and a wedding attractive flower in worldwide. In its flower colour, yellow is the rarest which is ten times the price of the other colours. However, the breeding of new yellow P. lactiflora varieties using genetic engineering is severely limited due to the little-known biochemical and molecular mechanisms underlying its characteristic formation.

RESULTS

In this study, two cDNA libraries generated from P. lactiflora chimaera with red outer-petal and yellow inner-petal were sequenced using an Illumina HiSeq™ 2000 platform. 66,179,398 and 65,481,444 total raw reads from red outer-petal and yellow inner-petal cDNA libraries were generated, which were assembled into 61,431 and 70,359 Unigenes with an average length of 628 and 617 nt, respectively. Moreover, 61,408 non-redundant All-unigenes were obtained, with 37,511 All-unigenes (61.08%) annotated in public databases. In addition, 6,345 All-unigenes were differentially expressed between the red outer-petal and yellow inner-petal, with 3,899 up-regulated and 2,446 down-regulated All-unigenes, and the flavonoid metabolic pathway related to colour development was identified using the Kyoto encyclopedia of genes and genomes database (KEGG). Subsequently, the expression patterns of 10 candidate differentially expressed genes (DEGs) involved in the flavonoid metabolic pathway were examined, and flavonoids were qualitatively and quantitatively analysed. Numerous anthoxanthins (flavone and flavonol) and a few anthocyanins were detected in the yellow inner-petal, which were all lower than those in the red outer-petal due to the low expression levels of the phenylalanine ammonialyase gene (PlPAL), flavonol synthase gene (PlFLS), dihydroflavonol 4-reductase gene (PlDFR), anthocyanidin synthase gene (PlANS), anthocyanidin 3-O-glucosyltransferase gene (Pl3GT) and anthocyanidin 5-O-glucosyltransferase gene (Pl5GT).

CONCLUSION

Transcriptome sequencing (RNA-Seq) analysis based on the high throughput sequencing technology was an efficient approach to identify critical genes in P. lactiflora and other non-model plants. The flavonoid metabolic pathway and glucide metabolic pathway were identified as relatived yellow formation in P. lactiflora, PlPAL, PlFLS, PlDFR, PlANS, Pl3GT and Pl5GT were selected as potential candidates involved in flavonoid metabolic pathway, which inducing inhibition of anthocyanin biosynthesis mediated yellow formation in P. lactiflora. This study could lay a theoretical foundation for breeding new yellow P. lactiflora varieties.