Skip to main content


Metabolic pathways and genes identified by RNA-seq analysis of barley near-isogenic lines differing by allelic state of the Black lemma and pericarp (Blp) gene



Some plant species have ‘melanin-like’ black seed pigmentation. However, the chemical and genetic nature of this ‘melanin-like’ black pigment have not yet been fully explored due to its complex structure and ability to withstand almost all solvents. Nevertheless, identification of genetic networks participating in trait formation is key to understanding metabolic processes involved in the expression of ‘melanin-like’ black seed pigmentation. The aim of the current study was to identify differentially expressed genes (DEGs) in barley near-isogenic lines (NILs) differing by allelic state of the Blp (black lemma and pericarp) locus.


RNA-seq analysis of six libraries (three replicates for each line) was performed. A total of 957 genome fragments had statistically significant changes in expression levels between lines BLP and BW, with 632 fragments having increased expression levels in line BLP and 325 genome fragments having decreased expression. Among identified DEGs, 191 genes were recognized as participating in known pathways. Among these were metabolic pathways including ‘suberin monomer biosynthesis’, ‘diterpene phytoalexins precursors biosynthesis’, ‘cutin biosynthesis’, ‘cuticular wax biosynthesis’, and ‘phenylpropanoid biosynthesis, initial reactions’. Differential expression was confirmed by real-time PCR analysis of selected genes.


Metabolic pathways and genes presumably associated with black lemma and pericarp colour as well as Blp-associated resistance to oxidative stress and pathogens, were revealed. We suggest that the black pigmentation of lemmas and pericarps is related to increased level of phenolic compounds and their oxidation. The effect of functional Blp on the synthesis of ferulic acid and other phenolic compounds can explain the increased antioxidant capacity and biotic and abiotic stress tolerance of black-grained cereals. Their drought tolerance and resistance to diseases affecting the spike may also be related to cuticular wax biosynthesis. In addition, upregulated synthesis of phytoalexins, suberin and universal stress protein (USP) in lemmas and pericarps of the Blp carriers may contribute to their increased disease resistance. Further description of the DEGs haplotypes and study of their association with physiological characteristics may be useful for future application in barley pre-breeding.


Grain colour is one of the stable taxonomic characters used for the description and classification of barley (Hordeum vulgare L., 2n =  = 14, HH) accessions. Barley grain at maturity may have white, yellow, purple, blue, or black pigmentation. In white grains, the pigments are absent, whereas two different classes of compounds determine the other types of colours. Yellow, purple, and blue are caused by a diverse subgroup of flavonoid compounds, synthesized in different grain tissues [1, 2]. Unlike the pigmentations caused by flavonoids, studies of black coloration are very scant from biochemical and molecular genetic points of view.

Black coloration is caused by ‘melanin-like’ pigments, which are groups of high molecular weight irregular polymers, arising in the course of oxidation and polymerization of phenolic compounds [3]. Because of their complex structure and resistance to almost all solvents, their chemical nature and the metabolic pathways leading to these black pigments have not yet been determined [4].

In barley, the ‘melanin-like’ pigments accumulate in the lemma and pericarp of grains during the later stages of growth, slightly before maturation of the spike [1, 5]. This trait is controlled by the monogenically inherited Blp1 (black lemma and pericarp 1) locus that has been mapped to chromosome 1HL [6, 7]. Three dominant alleles of the Blp1 locus conferring intensity of grain pigmentation have been reported: Blp1.b (B), Blp1.mb (B mb), and Blp1.g (B g) determine extreme black, medium black, and light black or grey colours, respectively [5, 8, 9].

Barley landraces with black grains were found in southwestern Asia, and sparsely in Tibet and China [10]. In Syria, the black-grained locally adopted landrace ‘Arabi Aswad’ was grown on approximately 75% of the barley growing area, particularly in the driest part of the country, unlike the white-seeded ‘Arabi Abiad’ landrace which was grown in more favourable conditions [11,12,13]. Compared to white-seeded varieties, black-seeded barley types are considered to be more drought tolerant, more vigorous in early stages of growth, more cold-tolerant, more prostrate, taller and faster-maturing [11, 12].

In addition to better viability under abiotic stress conditions, black-grained barley demonstrate higher resistance to fungal infection. FHB screening showed that among black-grained accessions, 20% were resistant, while among non-coloured accessions only 5% of genotypes were found to be resistant [14, 15]. Barley recombinant inbred lines (RILs) with black grains demonstrated lower Fusarium head blight incidence and lower accumulation of mycotoxin deoxynivalenol than yellow RILs, at least at some disease pressure levels [16]. Increased resistance to Fusarium disease has also been reported in oat genotypes having dark-coloured floral glumes in comparison with light-coloured ones [17].

Some studies of the phytochemical composition of barley grains with black colour have revealed that real ‘melanin-like’ black pigments may be present simultaneously with anthocyanins and other related co-pigments that contribute to the total phenol content [18, 19].

Comparative analysis of phenolic profiles and transcription of flavonoid biosynthesis pathway genes in black- and white-grained barley near-isogenic lines (NILs) demonstrated that none of the key flavonoid biosynthesis genes were upregulated in the black NIL, whereas in the presence of a dominant Blp1 allele, the total amount of phenolic compounds was enhanced in the black-grained line in comparison with the white-grained one [20]. These data suggest that neither coloured (proanthocyanidins, anthocyanins, phlobaphenes) nor uncoloured flavonoids participate in ‘melanin-like’ pigment biosynthesis, and instead, other classes of phenolic compounds are involved.

In the current study, to reveal genetic networks (and hence metabolic pathways) participating in formation of the ‘melanin-like’ black pigments, differentially expressed genes (DEGs) in black- vs white-grained NILs were determined by RNA-seq.


Six pooled libraries containing a total of 23 million short reads were produced as raw sequence data. After filtering, 2.49 million (15%) reads were removed from libraries before further analysis. Alignment of the libraries to the H. vulgare reference genome resulted in an average of 40% reads being mapped for each library. Metrics of library preprocessing and mapping are shown in Table 1.

Table 1 Metrics of libraries preprocessing and mapping

Among all genome fragments, 22,976 were identified as having non-zero expression levels. The fragments, their locations, annotations and respective RPKM values are listed in Additional file 1 (Additional file 1, Table ‘Genes’). These fragments include genes listed in genome assembly annotations as well as genome fragments for which expression has been detected by the Cufflinks pipeline and are further referred to as ‘transcripts’. RPKM values for the transcripts were used for library clustering (Fig. 1). Libraries BLP (2) and BLP (3) are highly similar in terms of gene expression, as are the BW (2) and BW (3) libraries. BLP (1) differs somewhat from the other two libraries belonging to the BLP line. Library BW (1) is further separated from the two other libraries belonging to the Bowman line (Fig. 1). Nevertheless, two clusters can be distinguished on the map, each containing all libraries of one of the two biological lines.

Fig. 1

Clustering of mapped libraries based on RPKM values for the transcripts (Additional file 1)

EdgeR detected 648 genes with lower and 1155 genes with higher expression levels in line BLP (Additional file 1, Table ‘EdgeR’), whereas Cufflinks detected 567 and 938 genes with lower and higher levels of expression in line BLP, respectively (Additional file 1, Table ‘Cufflinks’). In the list of genes mutually identified by both tools, 325 genes with decreased and 632 with increased expression levels in BLP were detected (Fig. 2; Additional file 1, Tables ‘UpReg DEGs Cufflinks & EdgeR’, ‘DownReg DEGs Cufflinks & EdgeR’). Thus, 957 transcripts with differential expression were detected in total, comprising 4.1% of all transcripts expressed in either of the biological lines.

Fig. 2

Number of differently expressed genes (DEGs) revealed by Cufflinks and EdgeR

PlantCyc database (BarleyCyc division) annotations include 196 differentially expressed genes, identified by Cufflinks and EdgeR, as participating in known BarleyCyc pathways (Additional file 2). The five pathways with the highest DEG PPS scores are ‘superpathway of pyrimidine deoxyribonucleotides de novo biosynthesis’, ‘UTP and CTP dephosphorylation II’, ‘flavonoid biosynthesis (in equisetum)’, ‘flavonol biosynthesis’ and ‘UTP and CTP dephosphorylation I’. Pathways with high numbers of DEGs include ‘suberin monomer biosynthesis’, ‘cutin biosynthesis’, ‘phenylpropanoid biosynthesis, initial reactions’, ‘superpathway of scopolin and esculin biosynthesis’ and ‘eugenol and isoeugenol biosynthesis’. Among additional pathways potentially related to the differing phenotypes of the NILs were ‘cuticular wax biosynthesis’ and ‘diterpene phytoalexins precursors biosynthesis’.

Using Blastp, homologies among Affymetrix Barley Genome Array sequences were found for 508 upregulated and 270 downregulated transcripts in BLP. Gene Ontology enrichment analysis showed that the GO term ‘fatty acid biosynthesis process’ is significantly (FDR = 7,5·10−4) enriched in transcripts with higher expression in BLP (Additional file 3), while cellular component GO terms ‘nucleoplasm’ (FDR = 0,03) and ‘chloroplast thylakoid’ (FDR = 0,029) and molecular function GO terms ‘transporter activity’ (FDR = 3,4·10−3) and ‘transcription factor activity’ (FDR = 5,5·10−3) were significantly overrepresented among transcripts with lower expression in line BLP (Additional files 4-5).

qPCR differential expression verification

qRT-PCR was used to evaluate the transcription of eight genes in seeds of the BLP and Bowman NILs, and results are summarized in Fig. 3. Six genes with higher expression levels in BLP and two genes with higher expression levels in the Bowman line were tested to verify the RNA-seq results.

Fig. 3

qRT-PCR validation of DEGs (relative mRNA levels of 8 genes obtained using gene-specific primers and cDNA of pericarp and lemma of BLP and BW lines)

Among the genes from the RNA-seq data showing increased expression levels in BLP are MLOC_62322 encoding phenylalanine ammonia-lyase (PAL), MLOC_76018 encoding Caffeic acid O-methyltransferase, MLOC_55029 encoding lipoxygenase B, MLOC_51393 encoding isoflavone reductase, MLOC_3802 encoding universal stress protein and MLOC_80301 encoding polyphenol oxidase. We also experimentally verified two genes from the RNA-seq data showing increased expression in the Bowman line: MLOC_54913 encoding lys-63-specific deubiquitinase and MLOC_16300 encoding apyrase 3. The log2(FC) values, p-values from RNA-seq data corrected with the Benjamini-Hochberg procedure and the qRT-PCR experimental fold changes and associated p-values for verified genes are shown in Table 2. The obtained experimental data are consistent with the RNA-seq predictions.

Table 2 The log2 (FC) values and corrected with Benjamini-Hochberg procedure p-value predicted by RNA-seq data and experimental fold changes and p-value for verified genes


The diversity of metabolic pathways up- and down-regulated in the presence of Blp (Blp1) suggests a pleiotropic nature of this gene. This finding is in agreement with previous observations at the phenotypic and physiological levels, suggesting that in addition to conferring black lemma and pericarp colour, the presence of the dominant Blp allele may contribute to total antioxidant level in barley grains (Additional file 6), abiotic stress tolerance [11,12,13] and disease resistance [14,15,16]. Investigation of DEG functions allowed us to identify genes and metabolic pathways presumably related to known phenotypic and physiological effects of the Blp locus.

Black ‘melanin-like’ pigments in plants have complex structures and can withstand almost all solvents of different chemical natures [4]. It has been proposed that these pigments are high molecular weight products of oxidation and polymerization of phenolic compounds [3, 21]. Indeed, an increase of total phenolic compounds in BLP in comparison to BW has been observed [20]. Here, we report significantly increased levels of PAL mRNA in the BLP line (Table 2, Fig. 3). PAL encodes phenylalanine ammonia-lyase, the key enzyme of the phenylpropanoid pathway, underlying biosynthesis of phenolic compounds in plants. Nevertheless, these quantitative changes in phenolic metabolites and PAL transcripts do not explain the qualitative differences between genotypes carrying dominant and recessive Blp alleles. Investigating additional DEGs, we propose that black pigmentation may be associated with high expression level of polyphenol oxidase (PPO) in BLP (Table 2, Fig. 3). PPO catalyses the oxidation of phenols to quinones. The products of the polymerization of quinones cause the darkening of some vegetables and fruits upon contact with air [22]. In humans, a similar enzyme, tyrosinase, is involved in the synthesis of melanin [23].

In wheat, black-spiked genotypes carry dominant allele Bg (Rg-A1c) on chromosome 1A, however on the short arm [24, 25], so barley Blp and wheat Bg are unlikely to be orthologues. The origin of black glume colour in wheat remains unknown, although it has been suggested to be related to phlobaphenes, a class of flavonoid compounds [26]. In barley, however, we did not observe phlobaphene biosynthesis structural genes among the identified DEGs.

Higher resistance of black grains to oxidative stress may also associated with increased phenolic content, since many phenolic compounds are known to be antioxidants [27]. Among the DEGs was a gene related to phenolic metabolism, which is of particular interest in relation to mechanisms increasing the antioxidant capacity of grains in the BLP line. This gene, with increased mRNA level in BLP, encodes Caffeic acid O-methyltransferase (Table 2, Fig. 3), which converts caffeic acid into ferulic acid [28]. Ferulic acid is a ubiquitous plant constituent that arises from the metabolism of phenylalanine and tyrosine. Ferulic acid, which accumulates primarily in seeds and leaves (both in its free form as well as covalently linked to lignin or other biopolymers), is known for its ability to terminate free radical chain reactions and hence serves an important antioxidant function to preserve the physiological integrity of cells in unfavourable conditions [29].

Higher resistance of black grains to oxidative stress may partially explain earlier findings that black seeded types are more drought tolerant, more vigorous in early stages of growth and more cold-tolerant than white seeded types [11]. Higher drought tolerance may also be related to cuticular wax biosynthesis, as revealed in our study as one of the metabolic pathways induced in the presence of Blp.

Antioxidant capacity can contribute to higher biotic stress resistance.

Black grain colour is associated with higher resistance to Fusarium in barley [14,15,16] and oat [17]. Fusarium head blight is a fungal disease favoured by humid conditions during flowering and early stages of kernel development. Due to better heat accumulation and accelerated drying, black-grained cereals may have an advantage. In addition to this possible effect of black pigmentation, the elevated resistance of dominant Blp allele carriers to Fusarium graminearum may also be explained by intensive production of phenolic compounds. Gunnaiah et al. applied an integrated metabolo-proteomic approach to decipher the mechanisms by which a wheat QTL (Fhb1) contributes to resistance against F. graminearum [30]. Fhb1 is not related to Blp or black grain pigmentation. However, it confers resistance to F. graminearum by inducing the phenylpropanoid pathway, leading to secondary cell wall thickening in rachises of wheat NILs with a resistant Fhb1 allele. Finally, increased transcriptional activity of the PAL gene can affect the biosynthesis of suberin and phytoalexins. Suberin, a waterproof waxy substance found in higher plants, can benefit black-grained cereals in humid conditions by preventing F. graminearum development, while phytoalexins (a wide class of plant molecules known for their antimicrobial properties) can act as toxins for the fungi. In addition to the PAL gene, we identified another DEG important for phenolic phytoalexin synthesis, encoding Isoflavone reductase (IFR) (Table 2, Fig. 3). Furthermore, DEGs necessary for the synthesis of terpenoid phytoalexins were identified (Additional file 2; ‘diterpene phytoalexins precursors biosynthesis’).

Improved resistance to F. graminearum in Blp carriers can also be explained by cuticular wax formation. Cuticular wax biosynthesis was revealed in our study as a metabolic pathway induced in the presence of Blp. Another mechanism of resistance can putatively be related to increased expression of the gene encoding universal stress protein (USP), which is usually expressed under environmental stress and stimulates further stress tolerance [31]. The BLP line has an increased USP level in optimal conditions, so it may potentially have increased tolerance to environmental stresses even in the absence of preliminary stimulation.


Metabolic pathways and genes presumably associated with black lemma and pericarp colour as well as with the pleiotropic effects of the Blp locus (i.e., resistance to oxidative stress and pathogens) were revealed. We suggest that black pigmentation of lemmas and pericarps may be related to increased level of phenolic compounds and their oxidation. The effect of Blp presence on the synthesis of ferulic acid and other phenolic compounds may be one of the explanations for the increased antioxidant capacity and biotic and abiotic stress tolerance of black-grained cereals. Their drought tolerance and resistance to diseases affecting the spike may also be related to cuticular wax biosynthesis. In addition, upregulated synthesis of phytoalexins, suberin and universal stress protein (USP) in lemmas and pericarps of the Blp carriers may also contribute to their increased disease resistance. Further description of the DEGs haplotypes and study of their association with physiological characteristics may be useful for future application in barley pre-breeding.


Plant materials, RNA extraction, genotyping and TAC measurements

Two H. vulgare cultivar Bowman [32] near-isogenic lines (NILs) were used in the study: ‘Black lemma and pericarp 1.b’, BLP (NGB20470) with black lemma and pericarp, and ‘Bowman From Fargo’, BW (NGB22812) with green lemma and uncoloured pericarp. These lines were kindly provided by the Nordic Gene Bank (NGB, Plants were grown at the ICG Greenhouse Core Facilities (Novosibirsk, Russia) under 12 h of light per day at 20–25 °C. Plants were genotyped and used for RNA-seq analysis. DNA was extracted for genotyping with microsatellite markers from leaf material (during leaf development stage, BBCH code 12) using a procedure described by Plaschke et al. [33]. For amplification of the barley 1H microsatellite loci, primers [34, 35] were used in PCRs conducted according to Röder et al. [36]. The resulting amplicons were visualized after separation through a 5% agarose (ACTGene, Inc., Piscataway, NJ, USA) gel. The scheme of donor fragments in chromosome 1H of BLP NIL is presented in Additional file 7.

Total RNA was extracted from lemma and pericarp tissue (early dough stage maturity, BBCH code 83) from BLP and Bowman lines using an RNeasy Plant Mini Kit (QIAGEN, Hilden, Germany). RNA extracted from several plants was pooled to exclude possible errors introduced by deviations of biological material. Three biological replicates were prepared for each genotype.

Extracts for TAC measurements were obtained by incubating 1 g of whole grain powder in either 10 ml of 1% HCl for 1 h at 37 °C (method 1) or 10 ml of 1% HCl in 40% ethanol for 30 min at 37 °C (method 2). Three replicates for each line were analysed in case of each methods). Measurements were performed with the antioxidant activity analyser ‘Blizar’ (Interlab Ltd., Russia) according to the manufacturer’s instructions. Results are presented in Additional file 6.

RNA preparation and sequencing

RNA was pooled into six libraries (three replicates for each line). ERCC spike-in mix was added to each library to serve as an internal control. Libraries were incubated with poly-T-tailed beads for poly-A enrichment. Library RNA was fragmented by incubation with nuclease enzymes. Sequencing was carried out with the IonTorrent platform.

Read preprocessing and mapping

Prinseq [37] tool v 0.20.4 was used to assess sequence quality and filter the libraries. Sequences with length below 50 nucleotides and mean Phred quality less than 20 were excluded from further analysis. Filtered libraries were mapped to H. vulgare genome assembly version 32,608 v 1.33 from the Ensembl Plants database ( Mapping was performed with TopHat2 tool v 2.1.0 [38] after constructing genome indexes using Bowtie2 v 2.2.6 software [39]. TopHat2 parameters ‘allowed mismatches’, ‘read edit distance’ and ‘read gap length’ were set at 4.

After mapping libraries to the reference genome, the resulting alignments were processed with the Cufflinks v 2.2.1 pipeline [40]. The number of reads mapped to each genome fragment, either expressed or annotated in the genome assembly (referred to here as ‘transcripts’), and respective RPKM (Reads Per Kilobase per Million mapped reads) [41] values were then used to detect differential expression of genes between the two studied barley lines. Furthermore, dendrogram of library clustering was plotted using the cummeRbund [42] package for R. Whole list of transcripts identified by cufflinks pipeline was used as a subset.

Gene expression analysis

Differential expression search was carried out using the Cufflinks pipeline cuffdiff utility and EdgeR v 3.8.6 [43] package for R v 3.1.3. Likelihood ratio tests were performed with EdgeR. Transcripts with total RPKM values below 15 in three libraries belonging to either of the biological lines were discarded from analysis. Transcripts with two-fold and higher differences in expression levels (|logFC| > 1) and FDR < 0.05 for EdgeR or q < 0.05 for Cufflinks were considered differentially expressed between the two lines. Transcripts with higher and lower levels of expression in BLP were analysed separately.

Transcript functional annotation

Blastp utility from the BLAST suite v 2.6.0 [44] was used to align sequenced of differentially expressed transcripts to sequences of the barley Affymetrix genome array downloaded from the AgriGO database [45]. For each transcript, the sequence from the barley Affymetrix genome array with the highest homology was identified. Lists of Affymetrix sequences corresponding to the transcripts with heightened and lowered levels of expression in BLP were processed with the AgriGO database Singular Enrichment Analysis tool.

Annotated genes with differential expression were analysed using the Plant Metabolic Network (PMN) database, BarleyCyc division (


For qPCR, isolated RNA was treated with DNase (QIAGEN RNase-Free DNase Set). Then, 0.7 μg of RNA was used to prepare single-stranded cDNA by reverse transcription, based on a RevertAid™ kit (Thermo Fisher Scientific Inc., Waltham, MA, USA) and a (dT)15 primer.

Primers were designed using IDT PrimerQuest software ( for eight DEGs (Table 3).

Table 3 Primers sequences constructed and used in qRT-PCRs

The Ubc (ubiquitin) gene sequence was used as a reference [46]. The subsequent qRT-PCR was based on a SYNTOL SYBR Green I kit (Syntol, Moscow, Russia). Three technical replicates of each reaction were performed.



Black lemma and pericarp




Differently expressed genes


Fusarium head blight


Gene ontology


Isoflavone reductase


Phenylalanine ammonia-lyase


Polyphenol oxidase


Pathway Perturbation Scores


Quantitative RT-PCR


Rids Per Kilobase per Million of reads


Reverse transcription


The total antioxidant content


  1. 1.

    Harlan HV. Some distinctions in our cultivated barleys with reference to their use in plant breeding. US Dept Agric Bul (Bulletin of the US Department of Agriculture). 1914;137:38.

  2. 2.

    Aastrup S, Outtrup H, Erdal K. Location of the proanthocyanidins in the barley grain. Carlsb Res Commun. 1984;49:105–9.

  3. 3.

    Britton G. The biochemistry of natural pigments. Cambridge: Cambridge University Press; 1983.

  4. 4.

    Jana BK, Mukherjee SK. Notes on the distribution of phytomelanin layer in higher plants – a short communication. J Pharm Biol. 2014;4:131–2.

  5. 5.

    Franckowiak JD, Lundqvist U. BGS 203. Barley Genet Newsl. 1997;26:209.

  6. 6.

    Buckley GFH. Inheritance in barley with special reference to the color of the caryopsis and lemma. Sci Agric. 1930;10:460–92.

  7. 7.

    Robertson DW, Wiebe GA, Shands RG, Hagberg A. A summary of linkage studies in cultivated barley, Hordeum species. Crop Sci. 1965;5(Suppl 3):33–43. 1954-1963

  8. 8.

    Woodward RW. Inheritance of melanin-like pigment in the glumes and caryopses of barley. J Agric Res. 1941;63:21–8.

  9. 9.

    Woodward RW. Linkage relationships between the allelomorphic series, B, Bmb, Bg, and Atat factors in barley. J Amer Soc Agron. 1942;34:659–61.

  10. 10.

    Takahashi R. Genetic features of east Asian barleys. In: Yasuda S, Konishi T, editors. Barley genetics V. Proc. fifth Int. barley genetics Symp. (Okayama;1986). Okayama: Okayama. Sanyo Press Co; 1987. p. 7–20.

  11. 11.

    Ceccarelli S, Grando S, van Leur JAG. Genetic diversity in barley landraces from Syria and Jordan. Euphytica. 1987;36:389–405.

  12. 12.

    Weltzien E. Evaluation of barley (Hordeum vulgare L) landrace populations originating from different growing regions in the near east. Plant Breed. 1988;101:95–106.

  13. 13.

    Bishaw Z, Struik PC, van Gastel AJG. Wheat and barley seed system in Syria: farmers, varietal perceptions, seed sources and seed management. Int J Plant Prod. 2011;5:323–47.

  14. 14.

    Zhou X, Chao M, Liang X. Screening and testing of barley varieties for scab resistance (in Chinese with English abstract). Acta Phytophylacica Sinica. 1991;18:261–5.

  15. 15.

    Choo TM. Breeding barley for resistance to fusarium head blight and mycotoxin accumulation. Plant Breed Rev. 2006;26:125–69.

  16. 16.

    Choo TM, Vigier BJ, Savard ME, Blackwell BA, Martinn RA, Wang JM, Yang J, Abdel-Aal ESM. Black barley as a means of mitigating deoxynivalenol contamination. Crop Sci. 2015;55:1096–103.

  17. 17.

    Loskutov IG, Blinova EV, Gavrilova OP, Gagkaeva TY. The valuable characteristics of oats genotypes and resistance to Fusarium disease. Vavilov Zh Genet Selekt. 2016;20:286–94.

  18. 18.

    Kim MJ, Hyun JN, Kim JA, Park JC, Kim MY, Kim JG, Lee SJ, Chun SC, Chung IM. Relationship between phenolic compounds, anthocyanins content and antioxidant activity in colored barley germplasm. J Agric Food Chem. 2007;55:4802–9.

  19. 19.

    Siebenhandl S, Grausgruber H, Pellegrini N, Del Rio D, Fogliano V, Pernice R, Berghofer E. Phytochemical profile of main antioxidants in different fractions of purple and blue wheat, and black barley. J Agric Food Chem. 2007;55:8541–7.

  20. 20.

    Shoeva OY, Mock H-P, Kukoeva TV, Börner A, Khlestkina EK. Regulation of the flavonoid biosynthesis pathway genes in purple and black grains of Hordeum vulgare. PLoS One. 2016;11:e0163782.

  21. 21.

    Varga M, Berkesi O, Darula Z, May NV, Palágyi A. Structural characterization of allomelanin from black oat. Phytochemistry. 2016;130:313–20.

  22. 22.

    Mayer AM. Polyphenol oxidases in plants and fungi: going places? A review. Phytochemistry. 2006;67:2318–31.

  23. 23.

    Sullivan ML. Beyond Brown: Polyphenol oxidases as enzymes of plant specialized metabolism. Frontiers Plant Sci. 2014;5:783.

  24. 24.

    Khlestkina EK, Pshenichnikova TA, Röder MS, Arbuzova VS, Salina EA, Börner A. Comparative mapping of genes for glume colouration and pubescence in hexaploid wheat (Triticum aestivum L.). Theor Appl Genet. 2006;113:801–7.

  25. 25.

    Khlestkina EK, Salina EA, Pshenichnikova TA, Röder MS, Börner A. Glume coloration in wheat: allelism test, consensus mapping and its association with specific microsatellite allele. Cereal Res Commun. 2009;37:37–43.

  26. 26.

    Khlestkina EK. Regulatory-target gene relationships in allopolyploid and hybrid genomes. In: Urbano KV, editor. Advances in genetics research, vol. 3. USA: NOVA Science Publishers, Inc; 2010. p. 311–28.

  27. 27.

    Rice-Evans C, Miller N, Paganga G. Antioxidant properties of phenolic compounds. Trends Plant Sci. 1997;2:152–9.

  28. 28.

    Shahidi F, Naczk M. Phenolics in food and Nutraceuticals. Florida: CRC Press LLC; 2004.

  29. 29.

    Graf E. Antioxidant potential of ferulic acid. Free Radic Biol Med. 1992;13(4):435–48.

  30. 30.

    Gunnaiah R, Kushalappa AC, Duggavathi R, Fox S, Somers DJ. Integrated metabolo-proteomic approach to decipher the mechanisms by which wheat QTL (Fhb1) contributes to resistance against Fusarium graminearum. PLoS One. 2012;7:e40695.

  31. 31.

    Siegele DA. Universal stress proteins in Escherichia coli. J Bacteriol. 2005;187:6253–4.

  32. 32.

    Franckowiak JD, Foster AE, Pederson VD, Pyler RE. Registration of Bowman barley. Crop Sci. 1985;25:883.

  33. 33.

    Plaschke J, Ganal MW, Röder MS. Detection of genetic diversity in closely related bread wheat using microsatellite markers. Theor Appl Genet. 1995;91:1001–7.

  34. 34.

    Ramsay L, Macaulay M, Degli Ivanissivich S, MacLean K, Cardle L, Fuller J, Edwards K, Tuvensson S, Morgante M, Massari A, Maesti E, Marmiroli N, Sjakste T, Ganal M, Powell W, Waugh R. A simple sequence repeat-based linkage map of barley. Genetics. 2000;156:1997–2005.

  35. 35.

    Li JZ, Sjakste TG, Röder MS, Ganal MW. Development and genetic mapping of 127 new microsatellite markers in barley. Theor Appl Genet. 2003;107:1021–7.

  36. 36.

    Röder MS, Korzun V, Wendehake K, Plaschke J, Tixier M-H, Leroy P, Ganal MW. A microsatellite map of wheat. Genetics. 1998;149:2007–23.

  37. 37.

    Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27:863–4.

  38. 38.

    Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.

  39. 39.

    Langmead B, Salzberg S. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–9.

  40. 40.

    Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Tanscript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotech. 2010;28:511–5.

  41. 41.

    Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;7:621–8.

  42. 42.

    Trapnell K, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Tanscript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechnology. 2010;28:511–5.

  43. 43.

    Robinson M, McCarthy D, Smyth G. EdgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.

  44. 44.

    Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421.

  45. 45.

    Zhou D, Xin Z, Yi L, Zhang Z, Zhen S. agriGO: a GO analysis toolkit for the agricultural community. Nucl Acids Res. 2010;68:64–70.

  46. 46.

    Himi E, Nisar A, Noda K. Colour genes (R and Rc) for grain and coleoptile upregulate flavonoid biosynthesis genes in wheat. Genome. 2005;48:747–54.

Download references


We thank Ms. Galina Generalova, Mrs. Tatiana Kukoeva and Mrs. Olga Zakharova (ICG, Novosibirsk, Russia) for technical assistance, and Drs Elena Gordeeva and Rimma Yudina for the help with analysis of TAC.


Publication costs were funded by the Russian Science Foundation (№ 16–14-00086). Growing of barley plants in ICG Plant Growth Core Facility is supported by ICG project 0324–2016-0001.

Availability of data and materials

Raw RNA-seq data is available at NCBI as BioProject PRJNA399215.

About this supplement

This article has been published as part of BMC Plant Biology Volume 17 Supplement 1, 2017: Selected articles from PlantGen 2017. The full contents of the supplement are available online at

Author information

AYG participated in extraction of RNA, preparation and sequencing of libraries, performed qRT-PCRs and statistical analysis, participated in the drafting the manuscript. NAS performed in silico analysis of the sequencing data and participated in the drafting the manuscript. OYS contributed to interpretation of data and participated in the drafting the manuscript. GVV prepared libraries, carried out sequencing of libraries on IonTorrent platform, participated in study coordination and interpretation of data. NVS participated in preparation and sequencing of libraries. AB provided plant material, contributed to conception of the study and interpretation of data. DAA participated in study coordination and interpretation of data. EKK conceived of the study, contributed to its conception, design and coordination, to interpretation of data and to revising the manuscript critically. All authors read and approved the final manuscript.

Correspondence to Anastasiya Y. Glagoleva.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Tables “Genes” (List of genome fragments identified by TopHat2, their genomic location, RPKM values in six samples and annotation), “Cufflinks” (Differentially expressed genes identified by Cufflinks), “EdgeR” (Differentially expressed genes identified by EdgeR), Table “UpRegDEGS_Cufflinks & EdgeR” (Differentially expressed genes identified by Cufflinks and EdgeR, upregulated in BLP line), and “DownRegDEGS_Cufflinks & EdgeR” (Differentially expressed genes identified by Cufflinks and EdgeR, downregulated in BLP line). (XLSX 2557 kb)

Additional file 2:

Links of differentially expressed genes identified by Cufflinks and EdgeR programs to BarleyCyc pathways. (XLSX 25 kb)

Additional file 3:

Ontology terms associated with biological process of protein products for genes with higher expression level in BLP line. (PNG 130 kb)

Additional file 4:

Ontology terms associated with cellular localization of protein products for genes with lower expression level in BLP line. (PNG 189 kb)

Additional file 5:

Ontology terms associated with molecular functions of genes with lower expression level in BLP line. (PNG 37 kb)

Additional file 6:

The total antioxidant content (mg/g) in grains of Bowman (BW) and BLP lines (equivalent to Gallic acid). Method 1: extracts were prepared by incubation of 1 g of whole grain powder in 10 ml of 1% HCl for 1 h at 37 °C. Method 2: 1 g of whole grain powder was incubated in 10 ml of 1% HCl in 40% ethanol for 30 min at 37 °C. Measurements were performed with the antioxidant activity analyser ‘Blizar’ (Interlab Ltd., Russia) according to the manufacturer’s instructions. * - difference between BW and BLP is significant (p ≤ 0.01; U-test). (PDF 183 kb)

Additional file 7:

Chromosome 1H scheme of Bowman (BW) and BLP lines. Blp donor segment on chromosome 1H remaining in the BLP NIL, revealed by microsatellite genotyping, is in gray. (PDF 49 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Glagoleva, A.Y., Shmakov, N.A., Shoeva, O.Y. et al. Metabolic pathways and genes identified by RNA-seq analysis of barley near-isogenic lines differing by allelic state of the Black lemma and pericarp (Blp) gene. BMC Plant Biol 17, 182 (2017) doi:10.1186/s12870-017-1124-1

Download citation


  • Barley
  • Near-isogenic lines
  • Phytomelanin
  • Black lemma
  • Black pericarp
  • RNA-seq
  • Differential expression