Skip to main content

Transcriptome analysis and characterization of genes associated to leaf tannin content in foxtail millet [Setaria italica (L.) P. Beauv.]

Abstract

Background

Chinese chestnut is an economically important tree species whose yield and quality are seriously affected by red spider attack. Tannins is one of the most important class secondary metabolites in plants, and is closely associated with plant defense mechanisms against insect and herbivory. In our previous studies, it was revealed that several low-tannin foxtail millet varieties growing under the Chinese chestnut trees could attract red spiders to feed on their leaves and protect the chestnut trees from the infestation of red spiders, meanwhile, the growth and yield of foxtail millet plants themselves were not greatly affected.

Results

To identify genes related to leaf tannin content and selection of foxtail millet germplasm resources with low tannin content for interplanting with Chinese chestnut and preventing the red spider attack, the leaves of 4 varieties with different levels of tannin content were harvested for comparative transcriptome analysis. In total, 335 differentially expressed genes (DEGs) were identified. For acquisition of gene functions and biological pathways they involved in, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses were performed, and several DEGs were found to possibly participate in the tannins biosynthesis pathway and transport processes of precursors. In addition, according to the PlantTFDB database, some transcription factors were predicted among the DEGs, suggesting their role in regulation of tannins biosynthesis pathway.

Conclusion

Our results provide valuable gene resources for understanding the biosynthesis and regulation mechanisms of tannins in foxtail millet, and pave the way for speeding up the breeding of low-tannin varieties through marker-assisted selection, which could be utilized for interplanting with Chinese chestnut trees to confer protection against red spider attack.

Peer Review reports

Introduction

Chinese chestnut (Castanea mollssina Blume) is an economically important tree species which originated from China. In Yanshan Mountain and Taihang Mountain area of Beijing–Tianjin–Hebei region, the major producing areas of Chinese chestnut in China, the planting area covered by Chinese chestnut has reached about 313.33 thousand hm2, and Chinese chestnut plays an irreplaceable role in economic development of this area. However, it was found that the yield and the quality of Chinese chestnut kernel decrease with the incidence of red spider (Tetranychus cinnbarinus), an agricultural insect pest which can damage more than 110 kinds of plants [1], thereby seriously affecting the income of local farmers. Thus, there is urgent need to find an effective approach to control the occurrence of red spiders on Chinese chestnut trees.

Insecticides have been applied to fight against red spiders, however, chemical pesticide residues pose hazards to both humans and the environment. Besides, trapping and killing red spiders, and release of natural enemies of red spiders were also effective measures, which significantly reduce the harm of red spiders to Chinese chestnut [2, 3]. Nevertheless, these biological control methods are time consuming and lead high economic cost. Therefore, it is important to explore a new way to prevent the red spider attack on Chinese chestnut, which takes into consideration the economic, social and environmental aspects.

Plant phenols have been regarded as one of the most common and widespread group of defensive compounds against insects in various studies [4,5,6]. Plant phenols range from simple phenol (MW 94) to polyphenols, such as anthocyanin pigments (MW 2000) and tannins (MW up to 20,000) [7]. Tannins have been reported to have strong deleterious effects on phytophagous insects affecting both growth and development by protein binding, reducing nutrient absorption efficiency and causing midgut lesions [8,9,10]. Mohamed and Abd-El Hameed found that the tannin content in leaves and seeds of faba bean (Vicia faba L.) plants is significantly increased in the insect-resistant genotypes as than in the susceptible genotypes [11], and Rizwan et al. reported that whitefly populations exhibit a significant and negative response to the leaf tannin content of cotton varieties [12], revealing the role of tannins in plant defense against insect damage.

Foxtail millet [Setaria italica (L.) P. Beauv.] is one of the most important commercial crops in Beijing–Tianjin–Hebei region of China, and has a strong resistance to drought and barren soil [13]. In our previous study on stereo interplanting system of foxtail millet with Chinese chestnut tree, we breed several early maturing and shade-tolerant foxtail millet materials, which could grow well under the spreading Chinese chestnut tree, realizing the high-efficient utilization of land resources under the trees. Meanwhile, we found that these materials could attract red spiders to feed on the leaves and protect the trees from the infestation of red spiders, while the growth and yield of foxtail millet plants themselves were not greatly affected. Moreover, the leaf tannin content of foxtail millet is about 10 times lower than the content of Chinese chestnut, which is probably one of the major reasons leading to attraction of foxtail millet plants for red spiders, according to the negative correlation of tannin content with insect infestation mentioned above [11, 12]. For this reason, identification of foxtail millet germplasm resources with low tannin content and using them for interplanting with Chinese chestnut may be an effective method to prevent the red spider attack on Chinese chestnut tree.

Molecular marker that tightly linked to desired genes is a powerful tool to screen genotypes of interest, which saves time and resources, and has been used in several plant species for selection of genotypes with low tannin content or zero tannin content. For example, two complementary recessive genes, zt-1 ant zt-2, control the absence of tannins in faba bean, and development of tannin-free cultivars will enhance the nutritional quality of this legume [14]. Gutierrez et al. developed Random Amplified Polymorphic DNA (RAPD) based Sequence-Characterised Amplified Region (SCAR) markers closely linked to zt-1 and zt-2 genes and validated their effectiveness and applicability in discriminating the genotypes differing in tannin content [15, 16]. Webb et al. mined a Single Nucleotide Polymorphism (SNP) marker for ZT1 [17]. Hou et al. reported Simple Sequence Repeat (SSR) markers for zt1 gene, among which a marker predicted the zt-1 genotypes with absolute accuracy [18]. Recently, Zanotto et al. reported a robust Kompetitive Allele-Specific PCR (KASP) marker associated with zt2 gene [19]. However, no potential candidate genes related to tannin content and molecular markers were reported in foxtail millet so far. Moreover, elucidation of the molecular mechanisms underlying tannins biosynthesis also needs identification and functional analysis of key genes.

In this study, to explore genes related to tannin content in foxtail millet, we compared the transcriptomes of foxtail millet varieties with different levels of tannin content. Through examination of differential gene expression, several potential candidate genes contributed to tannins biosynthesis in foxtail millet were identified, and their GO and KEGG enrichment analyses were performed. Moreover, the transcription factors involved in the regulation of tannins biosynthesis pathway were predicted by using the PlantTFDB database. Our results provide valuable gene resources for development of molecular markers and speeding up the breeding of low-tannin varieties, and will provide a basis for understanding the biological pathways and molecular regulation mechanisms behind tannins biosynthesis in foxtail millet.

Methods

Plant materials

Five foxtail millet hybrid varieties, 56229, 57295, 12950, 1121, HK950, and a conventional variety, JG32, which were preserved by Tangshan Normal University, were used in this study. All these materials were planted in Beiguan Village of Qianxi County, Tangshan, China. Three standard plants were selected for each variety, and middle part of the first, second and third leaves from the top of each plant at booting stage were sampled. Then the collected leaves of 3 standard plants for each variety were mixed and ground for tannin content determination.

Determination of tannin content

The chemical analysis to determine tannin content was carried out at ZKGX Research Institute of Chemical Technology (Plant Lab) of Chinese Institute of Chemical (Beijing, China). Tannin was determined by spectrophotometric method according to the standard NY/T 1600–2008 [20]. Briefly, the reaction mixture contained tannin extract solution, water, Folin-Denis reagent (containing Na2WO4·2H2O, H3Mo12O40P·XH2O and phosphomolybdic acid solution) and Na2CO3 solution was incubated at room temperature for 2 h. The absorbance was measured immediately at 765 nm using a spectrophotometer. Gallic acid was used for standard curve.

The production practice showed that the ability of foxtail millet hybrid HK950 to attract red spiders to feed was strong, while the ability of JG32 was weak. Therefore, the average (0.25) of leaf tannin content of these two varieties was used as a control to screen varieties with high or low leaf tannin content. According to the tannin content listed in Table S1, varieties 57295, 56229 and 12950 were considered as varieties with high tannin content, and 1121 was identified as a variety with low tannin content. These 4 varieties were then used for high-throughput sequencing.

RNA extraction, cDNA library construction and Illumina sequencing

To mine genes related to leaf tannin content, the middle part of the first, second and third leaves from the top of plants at booting stage were sampled from the 4 foxtail millet hybrid varieties, 57295, 56229, 12950 and 1121. Leaves collected from 3 standard plants for each variety were mixed and frozen in liquid nitrogen for total RNA isolation. Three biological replicates were performed for each sample.

Total RNA isolation, cDNA library construction and RNA-seq were performed by Berry Genomics Co., Ltd (Beijing, China). A total of 2 µg RNA with high quality per sample was adopted to construct cDNA libraries using VAHTS mRNA-seq v2 Library Prep Kit for Illumina following manufacturer’s recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using fragmentation buffer. First strand cDNA was synthesized and second strand cDNA synthesis was subsequently performed. Remaining overhangs were converted into blunt ends. After adenylation of 3’ ends of DNA fragments, adaptor with hairpin loop structure were ligated. Then the PCR was performed. At last, Qubit HS quantification, Agilent 2100 Bioanalyzer/Fragment Analyzer 5300 quality control, the final library size of about 350 bp. The libraries were sequenced on an Illumina NovaSeq platform to generate 150 bp paired-end reads (raw data) according to the manufacturer’s instructions.

Assembly of RNA-seq data

The raw data (raw reads) of fastq format were firstly processed through primary quality control. In this step, clean data (clean reads) were obtained by removing read pairs that contain N more than 3 or the proportion of base with quality value below 5 is more than 20%, in any end, or adapter sequence was founded. All the downstream analyses were based on the clean data with high quality.

Paired-end clean reads were aligned to the Setaria italica reference genome (Setaria_italica_v2.0, http://plants.ensembl.org/Setaria_italica/Info/Index) using HISAT2 version 2.0.6 [21] and featureCounts version 2.0.0 [22] was used to count the reads numbers mapped to each gene.

Differential expression genes analysis

The gene expression levels were determined as fragments per kilobase of transcript per million mapped reads (FPKM) value. The software edgeR version 3.3.3 [23] was used to detect the differential expression of genes between foxtail millet varieties with high leaf tannin content and low tannin content, and |log2(fold change)|> 1 and adjusted p-value (padj) < 0.05 were set as the threshold in the differentially expressed genes (DEGs) screening process.

Chromosome location

Chromosome positional data including the chromosome number and position of DEGs were retrieved from Ensembl Plants (http://plants.ensembl.org/index.html). The genes were then plotted individually onto the respective foxtail millet chromosomes and the resultant physical maps were displayed using MapChart software [24].

Quantitative real-time PCR (qRT-PCR) verification

Several DEGs were selected to verify the expression levels of RNA-Seq by qRT-PCR. The total RNA extraction, cDNA synthesis and qRT-PCR analysis were performed as described previously [25]. The sequences of primers applied in qRT-PCR were listed in Table S2. The PCR thermal profile was under the following cycle conditions: an initial 95 °C for 15 s, followed by 40 cycles at 95 °C for 10 s, and 60 °C for 31 s. All experiments were used in triplicates for each sample and the relative expression levels of each gene were calculated using 2−ΔΔCT method [26].

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses

The GO (http://www.geneontology.org) and KEGG (https://www.kegg.jp/) pathway enrichment analyses of DEGs were carried out using topGO [27] and KOBAS 3.0 [28] with default parameters, respectively.

Transcription factor identification

The transcription factor families were identified using the Plant Transcription Factor Database PlantTFDB 5.0 (http://planttfdb.gao-lab.org/prediction.php) [29].

Results

Identification of DEGs

To identify possible candidate genes involved in tannins biosynthesis, estimation of tannin content was carried out among 23 foxtail millet varieties, and varieties 57295, 56229 and 12950 were selected as varieties with high tannin content, while 1121 was considered as a variety with low tannin content (Table S1). Illumina sequencing was performed on leaves of these 4 samples, with 3 replicates on each sample (Table S3 and Figure S1), and a total of 30,141 genes were acquired. According to their FPKM value in different samples, the DEGs were identified.

In total, when compared with variety 1121, 717, 1124 and 2126 up-regulated DEGs were identified in 57295, 56229 and 12950, respectively (Fig. 1A). Among them, 167 DEGs were shared in 3 foxtail millet varieties with high tannin content. Meanwhile, there were 960, 1819 and 1103 down-regulated DEGs in 57295, 56229 and 12950, respectively, among which 168 DEGs were in common (Fig. 1B). The differential expression analysis provided amounts of potential candidate genes with respected to tannins biosynthesis in foxtail millet.

Fig. 1
figure 1

Venn diagram depicting the number of differentially expressed genes (DEGs) between foxtail millet varieties with high and low leaf tannin content. (A) and (B) shows the distribution of up-regulated and down-regulated DEGs in each comparison, respectively

Chromosome location of DEGs

Based on their chromosomal location, the 335 common DEGs shared in all comparisons were plotted onto 9 chromosomes of foxtail millet and the physical map was constructed, through which it was found that these genes were distributed unevenly on the 9 chromosomes (Fig. 2). For example, 71 genes were localized at chromosomes 1, while chromosomes 4 and 6 had 12 and 16 genes, respectively. The chromosome location analysis showed important positional information of candidate genes related to tannin content for further development of molecular markers.

Fig. 2
figure 2

Physical map of differentially expressed genes (DEGs) commonly existed in each comparison showing their chromosomal localization. The scale bar at left indicates chromosome length

qRT-PCR verification

To verify the reliability of gene changes in transcriptome analysis, 6 DEGs were randomly selected to analyze their transcriptional levels by qRT-PCR. As shown in Fig. 3, the alteration pattern of these genes was consistent with that of transcriptome analysis, indicating that the DEGs identified by comparative transcriptome analysis were reliable. Among these genes, SETIT_040859mg, SETIT_030369mg, SETIT_000657mg and SETIT_016714mg, were significantly up-regulated, whereas SETIT_013588mg and SETIT_026138mg were remarkably down-regulated in foxtail millet varieties with high leaf tannin content.

Fig. 3
figure 3

The expression pattern of 6 differentially expressed genes (DEGs) in the leaves of foxtail millet varieties with different levels of leaf tannin content. Transcription levels were verified by qRT-PCR with SiACTIN gene as an internal control. The X axis represents different foxtail millet varieties, and the Y axis (left side) indicates the relative expression level of selected genes determined by qRT-PCR (blue columns). The transcript level of genes in variety 56229 was set as 1.0, and error bars represent standard error (n = 3). The Y axis (right side) depicts the expression level of genes in RNA-seq data (evaluated by FPKM, red lines)

GO functional enrichment of DEGs

As shown in Fig. 4, the GO annotation results of 335 DEGs among different samples were obtained. These DEGs were classified into 3 categories, including biological process (BP), molecular function (MF), cellular component (CC) and 38 subcategories, in which 39.5%, 23.7% and 36.8% subcategories were classified as BP, MF and CC, respectively (Table S4). In BP terms, metabolic process (GO:0008152), cellular process (GO:0009987) and single-organism process (GO:0044699) showed significant enrichment. In MF terms, binding (GO:0005488) and catalytic activity (GO:0003824) indicated significant enrichment. In CC terms, cell (GO:0005623), cell part (GO:0044464), membrane (GO:0016020), membrane part (GO:0044425) and organelle (GO:0043226) had more numbers of genes. These results indicated that these biological processes, molecular functions and cellular components had undergone remarkable changes in foxtail millet varieties with different levels of tannin content.

Fig. 4
figure 4

GO enrichment analysis of differentially expressed genes (DEGs). The GO categories include biological processes, molecular functions, cellular components and 38 subcategories. Blue represents up-regulated DEGs in varieties with high leaf tannin content, and red represents down-regulated DEGs

KEGG pathway enrichment analysis of DEGs

To further identify the biological pathways and their function in the leaf tissues of foxtail millet, the 335 DEGs were subjected to KEGG pathway analysis. A total of 115 DEGs were assigned to 50 KEGG pathways, among which 26 DEGs (22.61%) were involved in metabolic pathways, similarly, 17 DEGs (14.78%) for biosynthesis of secondary metabolites, and 4 DEGs (3.48%) of phenylpropanoid biosynthesis (Table S5). The 50 pathways belong to 4 main categories that could be further divided into 17 subcategories. As shown in Fig. 5, the metabolism category had the highest proportion, followed by genetic information processing, environmental information processing and organismal systems.

Fig. 5
figure 5

KEGG classification of differentially expressed genes (DEGs). The name of KEGG pathway is shown on the Y axis, and the number of DEGs annotated to the pathway is shown on the X axis

Precursor molecules for tannins biosynthesis were reported to be derived from the phenylpropanoid pathway [30]. As per the KEGG analysis, the genes SETIT_009509mg and SETIT_029093mg present in phenylalanine/tyrosine ammonia-lyase (PTAL) [EC:4.3.1.25], SETIT_014043mg present in peroxidase [EC1.11.1.7] and SETIT_006847mg present in cinnamoyl-CoA reductase (CCR) [EC:1.2.1.44] were enriched into the phenylpropanoid biosynthesis pathway (Figure S2). Except for the gene SETIT_029093mg, they all had higher expression in leaf tissues of foxtail millet varieties with high tannin content as compared to gene expression in varieties with low tannin content.

Transcription factor prediction

The GO annotation indicated that regulation of transcription (GO:0006355) showed significant enrichment in biological process, implying that several DEGs may function as transcription factors (TFs) to involve in gene control and regulation of tannins biosynthesis in foxtail millet. Therefore, the potential TF genes were identified among the 335 common DEGs according to the family assignment rules illustrated in PlantTFDB [29]. A total of 20 DEGs were predicted to be involved in the regulation of transcription and were classified into 10 TF families (Fig. 6, Table S6). The most abundant TF families enriched were bHLH followed by FAR1 and B3. Among these TF genes, 13 were up-regulated and 7 were down-regulated.

Fig. 6
figure 6

Transcription factor families enriched in differentially expressed genes (DEGs). The number of members in each transcription factor family is presented within the brackets

Discussion

During growth and maturation period in plants, the secondary metabolites have essential roles. As one of the most important secondary metabolites, tannins is a diverse group of compounds that are related primarily in their ability to complex with proteins [31, 32], and is closely associated with plant defense mechanisms against insect. In this study, 335 DEGs were identified through Illumina sequencing among varieties with different levels of tannin content, which may be in relation to tannins biosynthesis in foxtail millet. Our findings provide valuable resources for investigation on the biosynthesis and regulation mechanisms of tannins in foxtail millet, and for development of molecular markers aimed at marker-assisted selection (MAS) for breeding of varieties with low tannin content, which could be used for interplanting with Chinese chestnut trees to confer protection against red spider attack.

According to their chemical structure and properties, tannins are divided into 2 main groups, condensed tannins (CTs) and hydrolysable tannins (HTs), and CTs are usually found higher concentrations in plants than HTs [33]. CTs, also known as proanthocyanidins (PAs), are synthesized via the phenylpropanoid and flavonoid pathways, and are widespread in the plant kingdom [30, 34]. In the present study, we found several genes were enriched into flavonoid and phenylpropanoid pathways according to the gene analyses. For example, the DEGs, SETIT_009941mg, SETIT_004762mg, SETIT_033256mg and SETIT_006373mg, were annotated to flavonoid biosynthetic process (GO:0009813), and SETIT_017426mg was annotated to flavonoid metabolic process (GO:0009812). As per the KEGG enrichment analysis, the 2 DEGs, SETIT_009509mg and SETIT_029093mg, which present in phenylalanine/tyrosine ammonia-lyase (PTAL) [EC:4.3.1.25], may code for enzymes in the early steps of the phenylpropanoid pathway [30, 35], and significantly affect the CTs biosynthesis and levels in foxtail millet. Our results suggested that these DEGs may be key genes involved in biosynthesis of tannins in foxtail millet, and be utilized for MAS breeding program and generating low-tannin content varieties. However, their functions still need to be further clarified.

The biosynthesis pathway of CTs is under complex control by multiple regulatory genes at the transcriptional level [36]. According to gene structure, these regulatory genes are classified into 6 different families, including bHLH, MYB, WD40, WRKY, MADS and WIP transcription factors, and most of the regulatory elements belong to the top 3 families. For example, the gene TRANSPARENT TESTA8 (VfTT8), encoding a bHLH transcription factor, was identified as the candidate for zt2, a gene controls absence of tannin in faba bean [37]. In tetraploid cotton, 2 R2R3-type MYB transcription factors, GhMYB36 and GhMYB10, were identified to involve in regulating CTs biosynthesis [38]. A WRKY transcription factor, McWRKY71, was identified in Malus crabapple, and overexpression of McWRKY71 promotes the accumulation of CTs in apple calli. Besides, it was found that McWRKY71 directly binds to the W-BOX element in the promoter of McANR, which promotes CTs accumulation, and McWRKY71 interacts with McMYB12, a functional regulator of CTs synthesis [39]. In this study, 20 genes were predicted to encode transcription factors categorized into 10 families, including 4 bHLH members, 1 MYB member and 2 WRKY members, and they were differentially expressed in foxtail millet varieties with different levels of tannin content. Therefore, these genes may be key candidates for involving in regulation of biosynthesis and accumulation of tannins in foxtail millet, and their function and regulation mechanism will be investigated in future work.

It is well known that CTs accumulate in the vacuole, while the known enzymes participate in the biosynthesis pathway mostly locate on the endoplasmic reticulum membranes or in the cytoplasm [40]. Hence, the CTs precursors must be transported from the site of synthesis to the site of storage by 2 models: the vesicle trafficking-mediated model and the membrane transporter-mediated model [41, 42]. The latter is established based on 3 types of transporters, glutathione S-transferase (GST), multidrug and toxic compound extrusion (MATE) and H+-ATPase. For example, a GST gene family member in Arabidopsis thaliana, TT19, was reported to localize both in the cytoplasm and on the tonoplast, and is required for transferring the anthocyanins into the vacuole [43]. In the current study, based on the GO enrichment analysis, we found a DEG, SETIT_009171mg, was annotated to transmembrane transport process (GO:0055085) and ATPase activity, coupled to transmembrane movement of substances (GO:0042626). Furthermore, the DEG SETIT_002832mg was categorized to glutathione metabolic process (GO:0006749) and glutathione transferase activity (GO:0004364). As per the KEGG analysis, this gene was present in glutathione S-transferase [EC:2.5.1.18]. Whether these genes involve in transport processes required for CTs biosynthesis is also an interesting focus of our future work.

Conclusion

In this study, by using RNA-seq data derived from 4 foxtail millet varieties with different levels of tannin content, we identified 335 differentially expressed genes, which were potential candidate genes contributed to tannins biosynthesis in foxtail millet. To understand their functions and biological pathways they involved in, GO and KEGG enrichment analyses were performed. Moreover, the transcription factors involved in the regulation of tannins biosynthesis pathway were predicted by using the PlantTFDB database. These results will provide a basis for investigation on biosynthesis and regulation mechanisms of tannins in foxtail millet, and provide a resource for breeding of low-tannin varieties through marker-assisted selection, which could be utilized for interplanting with Chinese chestnut tree and prevent red spider attack.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available in the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI) under the accession number PRJNA772942.

Abbreviations

RAPD:

Random Amplified Polymorphic DNA

SCAR:

Sequence-Characterised Amplified Region

SNP:

Single Nucleotide Polymorphism

SSR:

Simple Sequence Repeat

KASP:

Kompetitive Allele-Specific PCR

FPKM:

Fragments per kilobase of transcript per million mapped reads

DEG:

Differentially expressed gene

qRT-PCR:

Quantitative real-time PCR

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

BP:

Biological process

MF:

Molecular function

CC:

Cellular component

PTAL:

Phenylalanine/tyrosine ammonia-lyase

CCR:

Cinnamoyl-CoA reductase

TF:

Transcription factor

MAS:

Marker-assisted selection

CTs:

Condensed tannins

HTs:

Hydrolysable tannins

TT8:

TRANSPARENT TESTA8

GST:

Glutathione S-transferase

MATE:

Multidrug and toxic compound extrusion

References

  1. Huang Y, Song H, Ding P, Guo L, Xu Y, Lu J, et al. Efficacy evaluation of 25% abamectin·etoxazole SC against red spider (Tetranychus cinnbarinus) in ornamental rose. Plant Diseases and Pests. 2017;8(6):39–40.

    Google Scholar 

  2. Zhao H, Fan R. Occurrence and control methods of red spiders of Chestnut. Modern Rural Sci Technol. 2019;9:34.

    Google Scholar 

  3. Wang X, Wen X, Lu C. Research advance in controlling insect pests on Castanea mollissima in China. J Hebei Normal Univ Sci Technol. 2010;24(1):39–43.

    Google Scholar 

  4. Barbehenn RV, Martin MM. Tannin sensitivity in larvae of Malacosoma disstria (Lepidoptera): Roles of the peritrophic envelope and midgut oxidation. J Chem Ecol. 1994;20(8):1985–2001.

    Article  CAS  PubMed  Google Scholar 

  5. Barbehenn RV, Martin MM, Hagerman AE. Reassessment of the roles of the peritrophic envelope and hydrolysis in protecting polyphagous grasshoppers from ingested hydrolyzable tannins. J Chem Ecol. 1996;22(10):1901–19.

    Article  CAS  PubMed  Google Scholar 

  6. Henn MW. Adsorption von Tanninen aus Eichenblattern an Kohlen e hydrathaltiger Verbindungen. Mitt Dtsch Entomol Ges. 1997;11:495–9.

    Google Scholar 

  7. Rattan RS. Mechanism of action of insecticidal secondary metabolites of plant origin. Crop Prot. 2010;29(9):913–20.

    Article  CAS  Google Scholar 

  8. Sharma HC, Sujana G, Rao DM. Morphological and chemical components of resistance to pod borer, Helicoverpa armigera in wild relatives of pigeonpea. Arthropod-Plant Inte. 2009;3(3):151–61.

    Article  Google Scholar 

  9. Sharma HC, Agarwal RA. Role of some chemical components and leaf hairs in varietal resistance in cotton to jassid, Amrasca biguttula biguttula Ishida. J Entomol Res. 1983;7(2):145–9.

    CAS  Google Scholar 

  10. Barbehenn RV, Constabel CP. Tannins in plant–herbivore interactions. Phytochemistry. 2011;72(13):1551–65.

    Article  CAS  PubMed  Google Scholar 

  11. Mohamed HI, Abd-El Hameed AG. Molecular and biochemical markers of some Vicia faba L. genotypes in response to storage insect pests infestation. J Plant Interact. 2014;9(1):618–26.

    Article  CAS  Google Scholar 

  12. Rizwan M, Abro S, Asif MU, Hameed A, Mahboob W, Deho ZA, et al. Evaluation of cotton germplasm for morphological and biochemical host plant resistance traits against sucking insect pests complex. J Cotton Res. 2021;4(1):1–8.

    Article  CAS  Google Scholar 

  13. Liu Z, Bai G, Zhang D, Zhu C, Xia X, Cheng R, et al. Genetic diversity and population structure of elite foxtail millet [Setaria italica (L.) P. Beauv.] germplasm in China. Crop Sci. 2011;51(4):1655–63.

    Article  Google Scholar 

  14. Picard J. Aperçu sur l’hérédité du caractère absence de tanins dans les graines de féverole (Vicia faba L.). Annales de l’Ame´lioration des Plantes. 1976;26:101–6.

    Google Scholar 

  15. Gutierrez N, Avila CM, Rodriguez-Suarez C, Moreno MT, Torres AM. Development of SCAR markers linked to a gene controlling absence of tannins in faba bean. Mol Breeding. 2006;19(4):305–14.

    Article  CAS  Google Scholar 

  16. Gutierrez N, Avila CM, Moreno MT, Torres AM. Development of SCAR markers linked to zt-2, one of the genes controlling absence of tannins in faba bean. Aust J Agr Res. 2008;59(1):62–8.

    Article  CAS  Google Scholar 

  17. Webb A, Cottage A, Wood T, Khamassi K, Hobbs D, Gostkiewicz K, et al. A SNP-based consensus genetic map for synteny-based trait targeting in faba bean (Vicia faba L.). Plant Biotechnol J. 2016;14(1):177–85.

    Article  CAS  PubMed  Google Scholar 

  18. Hou W, Zhang X, Yan Q, Li P, Sha W, Tian Y, et al. Linkage map of a gene controlling zero tannins (zt-1) in faba bean (Vicia faba L.) with SSR and ISSR markers. Agronomy. 2018;8(6):80.

    Article  CAS  Google Scholar 

  19. Zanotto S, Vandenberg A, Khazaei H. Development and validation of a robust KASP marker for zt2 locus in faba bean (Vicia faba). Plant Breeding. 2020;139(2):375–80.

    Article  CAS  Google Scholar 

  20. Chinese Ministry of Agriculture. Determination of tannin content in fruit, vegetable and derived product-Spectrophotometry method: NY/T 1600–2008. In: Beijing: China Agriculture Press; 2008.

  21. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.

    Article  CAS  PubMed  Google Scholar 

  23. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009;26(1):139–40.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Voorrips RE. MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002;93(1):77–8.

    Article  CAS  PubMed  Google Scholar 

  25. Liu D, Cui Y, Zhao Z, Li S, Liang D, Wang C, et al. Genome-wide identification and characterization of the BES/BZR gene family in wheat and foxtail millet. BMC Genomics. 2021;22(1):1–14.

    Article  Google Scholar 

  26. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3(6):1101–8.

    Article  CAS  PubMed  Google Scholar 

  27. Yu G, Li F, Qin Y, Bo X, Wu Y, Wang S. GOSemSim: an R package for measuring semantic similarity among GO terms and gene products. Bioinformatics. 2010;26(7):976–8.

    Article  CAS  PubMed  Google Scholar 

  28. Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39(Web Server issue):W316-322.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Jin J, Tian F, Yang DC, Meng YQ, Kong L, Luo J, et al. PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res. 2017;45(D1):D1040–5.

    Article  CAS  PubMed  Google Scholar 

  30. Singh V, Goel R, Pande V, Asif MH, Mohanty CS. De novo sequencing and comparative analysis of leaf transcriptomes of diverse condensed tannin-containing lines of underutilized Psophocarpus tetragonolobus (L.) DC. Sci Rep. 2017;7:1–13.

    CAS  Google Scholar 

  31. Fahey GC Jr, Jung HG. Phenolic compounds in forages and fibrous feedstuffs. Toxicants of plant origin. 1989;4:123–90.

    Google Scholar 

  32. Hassanpour S, Sadaghian M, MaheriSis N, Eshratkhah B, ChaichiSemsari M. Effect of condensed tannin on controlling faecal protein excretion in nematode-infected sheep: in vivo study. J Am Sci. 2011;7(5):896–900.

    Google Scholar 

  33. Hassanpour S, Maheri-Sis N, Eshratkhah B, Mehmandar FB. Plants and secondary metabolites (Tannins): A Review. Int J Forest, Soil and Erosion. 2011;1(1):47–53.

    Google Scholar 

  34. He F, Pan QH, Shi Y, Duan CQ. Biosynthesis and genetic regulation of proanthocyanidins in plants. Molecules. 2008;13(10):2674–703.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Paolocci F, Bovone T, Tosti N, Arcioni S, Damiani F. Light and an exogenous transcription factor qualitatively and quantitatively affect the biosynthetic pathway of condensed tannins in Lotus corniculatus leaves. J Exp Bot. 2005;56(414):1093–103.

    Article  CAS  PubMed  Google Scholar 

  36. Tian L, Pang Y, Dixon RA. Biosynthesis and genetic engineering of proanthocyanidins and (iso) flavonoids. Phytochem Rev. 2008;7(3):445–65.

    Article  CAS  Google Scholar 

  37. Gutierrez N, Avila CM, Torres AM. The bHLH transcription factor VfTT8 underlies zt2, the locus determining zero tannin content in faba bean (Vicia faba L.). Sci Rep. 2020;10(1):1–10.

    Article  CAS  Google Scholar 

  38. Lu N, Roldan M, Dixon RA. Characterization of two TT2-type MYB transcription factors regulating proanthocyanidin biosynthesis in tetraploid cotton, Gossypium hirsutum. Planta. 2017;246(2):323–35.

    CAS  PubMed  Google Scholar 

  39. Zhang J, Wang Y, Mao Z, Liu W, Ding L, Zhang X, et al. Transcription factor McWRKY71 induced by ozone stress regulates anthocyanin and proanthocyanidin biosynthesis in Malus crabapple. Ecotoxicol Environ Saf. 2022;232:113274.

    Article  CAS  PubMed  Google Scholar 

  40. Saslowsky D, Winkel-Shirley B. Localization of flavonoid enzymes in Arabidopsis roots. Plant J. 2001;27(1):37–48.

    Article  CAS  PubMed  Google Scholar 

  41. Grotewold E, Davies K. Trafficking and sequestration of anthocyanins. Nat Prod Commun. 2008;3(8):1251–8.

    CAS  Google Scholar 

  42. Zhao J, Dixon RA. The ‘ins’ and ‘outs’ of flavonoid transport. Trends Plant Sci. 2010;15(2):72–80.

    Article  CAS  PubMed  Google Scholar 

  43. Sun Y, Li H, Huang JR. Arabidopsis TT19 functions as a carrier to transport anthocyanin from the cytosol to tonoplasts. Mol Plant. 2012;5(2):387–400.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by S&T Program of Hebei (21326337D and 21326302D).

Author information

Authors and Affiliations

Authors

Contributions

Z.L. and S.L. contributed to the study conception and design. D.L., Z.Z., S.L. Y.C., and J.Z. performed the experiments. D.L. and Z.Z. analyzed the data. Y.C. drafted the manuscript. Z.L. and S.L. revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Zhengli Liu.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethics approval and consent to participate

The foxtail millet used in this study is not endangered or a protected species in China, and it was preserved by Tangshan Normal University and planted in Beiguan Village of Qianxi County, Tangshan, China. The identification of plant materials was undertaken by Professor Zhengli Liu. The appropriate permissions and/or licenses for sample collection were obtained. We comply with relevant institutional, national, and international guidelines and legislation for plant study.

Consent for publication

Not applicable.

Competing interest

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1:

 Supplementary file Table S1. The tannin content in leaves of different foxtail millet varieties. TableS2. Primer sequences used for quantitative real-time PCR. Table S3. Characteristics of the RNA-sequencing data obtained from analysis of 4 leaf samples of foxtail millet. Table S4. GO annotation of the differentially expressed genes shared in all comparisons. Table S5. KEGG pathway annotation and classification of the differentially expressed genes shared in all comparisons. Table S6. The transcription factor (TF) genes predicted using PlantTFDB. Up and down indicate gene is up-regulated and down-regulated in 3 foxtail millet varieties with high leaf tannin content, respectively. Figure S1. Heatmap of pairwise Pearson correlation coefficients (R2) between samples. Figure S2. KEGG analysis of phenylpropanoid pathway in foxtail millet. The pathway map image is obtained from KEGG [1-3].

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, S., Cui, Y., Liu, D. et al. Transcriptome analysis and characterization of genes associated to leaf tannin content in foxtail millet [Setaria italica (L.) P. Beauv.]. BMC Genomics 23, 512 (2022). https://doi.org/10.1186/s12864-022-08746-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-022-08746-8

Keywords