SIFGD: Setaria italica Functional Genomics Database

C4 cereal and biofuel crop microbiomes

Microbiomes provide multiple life-support functions for plants, including nutrient acquisition and tolerance to abiotic and biotic stresses. Considering the importance of C4 cereal and biofuel crops for food security under climate change conditions, more attention has been given recently to C4 plant microbiome assembly and functions. Here, we review the current status of C4 cereal and biofuel crop microbiome research with a focus on beneficial microbial traits for crop growth and health. We highlight the importance of environmental factors and plant genetics in C4 crop microbiome assembly and pinpoint current knowledge gaps. Finally, we discuss the potential of foxtail millet as a C4 model species and outline future perspectives of C4 plant microbiome research.

Amaranth Genomic Resource Database: an integrated database resource of Amaranth genes and genomics

Amaranth (Amaranthus L.) is native to Mexico and North America, where it was cultivated thousands of years ago, but now amaranth is grown worldwide. Amaranth is one of the most promising food crops with high nutritional value and belongs to the family Amaranthaceae. The high-quality genome assembly of cultivated amaranth species (A. hypochondriacus, A. cruentus) and wild/weedy species (A. tuberculatus, A. hybridus, and A. palmeri) has already been reported; therefore, we developed an Amaranth Genomic Resource Database (AGRDB) to provide access to all the genomic information such as genes, SSRs, SNPs, TFs, miRNAs, and transporters in one place. The AGRDB database contains functionally annotated gene information with their sequence details, genic as well as genomic SSRs with their three sets of primers, transcription factors classified into different families with their sequence information and annotation details, putative miRNAs with their family, sequences, and targeted gene details, transporter genes with their superfamily, trans-membrane domain details, and details of genic as well as nongenic SNPs with 3' and 5' flanking sequence information of five amaranth species. A database search can be performed using the gene ID, sequence ID, sequence motif, motif repeat, family name, annotation keyword, scaffold or chromosome numbers, etc. This resource also includes some useful tools, including JBrowse for the visualization of genes, SSRs, SNPs, and TFs on the respective amaranth genomes and BLAST search to perform a BLAST search of the user's query sequence against the amaranth genome as well as protein sequences. The AGRDB database will serve as a potential platform for genetic improvement and characterization of this futuristic crop. The AGRDB database will be accessible via the link: http://www.nbpgr.ernet.in:8080/AmaranthGRD/.

MDSi: Multi-omics Database for Setaria italica

BACKGROUND

Foxtail millet (Setaria italica) harbors the small diploid genome (~ 450 Mb) and shows the high inbreeding rate and close relationship to several major foods, feed, fuel and bioenergy grasses. Previously, we created a mini foxtail millet, xiaomi, with an Arabidopsis-like life cycle. The de novo assembled genome data with high-quality and an efficient Agrobacterium-mediated genetic transformation system made xiaomi an ideal C₄ model system. The mini foxtail millet has been widely shared in the research community and as a result there is a growing need for a user-friendly portal and intuitive interface to perform exploratory analysis of the data.

RESULTS

Here, we built a Multi-omics Database for Setaria italica (MDSi, http://sky.sxau.edu.cn/MDSi.htm ), that contains xiaomi genome of 161,844 annotations, 34,436 protein-coding genes and their expression information in 29 different tissues of xiaomi (6) and JG21 (23) samples that can be showed as an Electronic Fluorescent Pictograph (xEFP) in-situ. Moreover, the whole-genome resequencing (WGS) data of 398 germplasms, including 360 foxtail millets and 38 green foxtails and the corresponding metabolic data were available in MDSi. The SNPs and Indels of these germplasms were called in advance and can be searched and compared in an interactive manner. Common tools including BLAST, GBrowse, JBrowse, map viewer, and data downloads were implemented in MDSi.

CONCLUSION

The MDSi constructed in this study integrated and visualized data from three levels of genomics, transcriptomics and metabolomics, and also provides information on the variation of hundreds of germplasm resources that can satisfies the mainstream requirements and supports the corresponding research community.

Comparison of SWEET gene family between maize and foxtail millet through genomic, transcriptomic, and proteomic analyses

Foxtail millet [Setaria italica (L.) P. Beauv.] does not show high yield and biomass compared with maize (Zea mays L.) although it is a C₄ crop with the potential for high productivity. Because SWEET genes, which are important for sugar transport in plants, play critical roles in biomass production and seed filling in crops, genome-wide, transcriptomic, and proteomic comparison on SWEET gene family between these two species would provide some clues for unlocking this issue. In our study, 24 SWEET genes were identified in foxtail millet and maize. Sequence-based bioinformatics combined with gene expression analyses identified several candidate functional orthologs in these two species. A comparative analysis on expression characteristics of SWEET genes and proteins between maize and foxtail millet indicate that not only some critical major SWEET proteins show significant upregulation in maize compared with their orthologs in foxtail millet, but also there are more quantities of maize SWEET genes showing high expressions than that of foxtail millet genes, suggesting that compared with foxtail millet, maize possesses higher capacity of sugar transport, the crucial determinant for crop yield and biomass. These results provide a basis on revealing why foxtail millet exhibits low yield and biomass although it is a C₄ crop with the potential for high productivity.

Multi-omics intervention in Setaria to dissect climate-resilient traits: Progress and prospects

Millets constitute a significant proportion of underutilized grasses and are well known for their climate resilience as well as excellent nutritional profiles. Among millets, foxtail millet (Setaria italica) and its wild relative green foxtail (S. viridis) are collectively regarded as models for studying broad-spectrum traits, including abiotic stress tolerance, C₄ photosynthesis, biofuel, and nutritional traits. Since the genome sequence release, the crop has seen an exponential increase in omics studies to dissect agronomic, nutritional, biofuel, and climate-resilience traits. These studies have provided first-hand information on the structure, organization, evolution, and expression of several genes; however, knowledge of the precise roles of such genes and their products remains elusive. Several open-access databases have also been instituted to enable advanced scientific research on these important crops. In this context, the current review enumerates the contemporary trend of research on understanding the climate resilience and other essential traits in Setaria, the knowledge gap, and how the information could be translated for the crop improvement of related millets, biofuel crops, and cereals. Also, the review provides a roadmap for studying other underutilized crop species using Setaria as a model.

MicroRNAs expression profiles in early responses to different levels of water deficit in Setaria viridis

Water deficit is a major constraint for crops of economic importance in almost all agricultural regions. However, plants have an active defense system to adapt to these adverse conditions, acting in the reprogramming of gene expression responsible for encoding microRNAs (miRNAs). These miRNAs promote the regulation to the target gene expression by the post-transcriptional (PTGS) and transcriptional gene silencing (TGS), modulating several pathways including defense response to water deficit. The broader knowledge of the miRNA expression profile and its regulatory networks in response to water deficit can provide evidence for the development of new biotechnological tools for genetic improvement of several important crops. In this study, we used Setaria viridis accession A10.1 as a C4 model plant to widely investigate the miRNA expression profile in early responses to different levels of water deficit. Ecophysiological studies in Setaria viridis under water deficit and after rewatering demonstrated a drought tolerant accession, capable of a rapid recovery from the stress. Deep small RNA sequencing and degradome studies were performed in plants submitted to drought to identify differentially expressed miRNA genes and their predicted targets, using in silico analysis. Our findings showed that several miRNAs were differentially modulated in response to distinctive levels of water deficit and after rewatering. The predicted mRNA targets mainly corresponded to genes related to cell wall remodeling, antioxidant system and drought-related transcription factors, indicating that these genes are rapidly regulated in early responses to drought stress. The implications of these modulations are extensively discussed, and higher-effect miRNAs are suggested as major players for potential use in genetic engineering to improve drought tolerance in economically important crops, such as sugarcane, maize, and sorghum.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01226-z.

Evolutionary analysis and functional characterization of SiBRI1 as a Brassinosteroid receptor gene in foxtail millet

Brassinosteroids (BRs) play important roles in plant growth and development. Although BR receptors have been intensively studied in Arabidopsis, those in foxtail millet remain largely unknown. Here, we show that the BR signaling function of BRASSINOSTEROID INSENSITIVE 1 (BRI1) is conserved between Arabidopsis and foxtail millet, a new model species for C4 and Panicoideae grasses. We identified four putative BR receptor genes in the foxtail millet genome: SiBRI1, SiBRI1-LIKE RECEPTOR KINASE 1 (SiBRL1), SiBRL2 and SiBRL3. Phylogenetic analysis was used to classify the BR receptors in dicots and monocots into three branches. Analysis of their expression patterns by quantitative real-time PCR (qRT-PCR) showed that these receptors were ubiquitously expressed in leaves, stems, dark-grown seedlings, roots and non-flowering spikelets. GFP fusion experiments verified that SiBRI1 localized to the cell membrane. We also explored the SiBRI1 function in Arabidopsis through complementation experiments. Ectopic overexpression of SiBRI1 in an Arabidopsis BR receptor loss-of-function mutant, bri1-116, mostly reversed the developmental defects of the mutant. When SiBRI1 was overexpressed in foxtail millet, the plants showed a drooping leaf phenotype and root development inhibition, lateral root initiation inhibition, and the expression of BR synthesis genes was inhibited. We further identified BRI1-interacting proteins by immunoprecipitation (IP)-mass spectrometry (MS). Our results not only demonstrate that SiBRI1 plays a conserved role in BR signaling in foxtail millet but also provide insight into the molecular mechanism of SiBRI1.

Identification of QTL for resistance to leaf blast in foxtail millet by genome re-sequencing analysis

Three QTL for resistance to leaf blast were identified on chromosomes 1, 2, and 8 of the foxtail millet cultivar Yugu 5. Leaf blast disease of foxtail millet (Setaria italica) is caused by Pyricularia spp., can infect all the aboveground parts of plants, and is the most frequently observed blast disease in China. Lack of information on genetic control of disease resistance impedes developing leaf blast-resistant cultivars. An F₆ recombinant inbred line (RIL) population from the cross Yugu 5 × Jigu 31 was phenotyped for its reactions to leaf blast in six field trials in the naturally diseased nurseries. An ultra-density genetic linkage map was constructed using 35,065 single nucleotide polymorphism (SNP) markers generated by sequencing of the RIL population. Three QTL, QLB-czas1, QLB-czas2, and QLB-cazas8, were detected in the genomic intervals of 276.6 kb, 1.62 Mb, and 1.75 Mb on chromosomes 1, 2, and 8 of Yugu 5, which explained 14-17% (2 environments), 9% (5 environments), and 12-20% (6 environments) of the phenotypic variations. Bulked segregant analysis (BSA) and RNA sequencing (BSR-Seq) method identified common SNPs that fell in the genomic region of QLB-czas8, providing additional evidence of localization of this QTL. Three and 19 predicted genes were annotated to be associated with disease resistance in the genomic intervals for QLB-czas2 and QLB-czas8. Due to their unique positions, these QTL appear to be new loci conferring resistance to leaf blast. The identification of these new resistance QTL will be useful in cultivar development and the study of the genetic control of blast resistance in foxtail millet.

HpeNet: Co-expression Network Database for de novo Transcriptome Assembly of Paeonia lactiflora Pall

The herbaceous peony (Paeonia lactiflora Pall.) is a well-known ornamental flowering and pharmaceutical plant found in China. Its high medicinal value has long been recognized by traditional Chinese medicine (as Radix paeoniae Alba and Radix paeoniae Rubra), and it has become economically valued for its oilseed in recent years; like other Paeonia species, it has been identified as a novel resource for the α-linolenic acid used in seed oil production. However, its genome has not yet been sequenced, and little transcriptome data on Paeonia lactiflora are available. To obtain a comprehensive transcriptome for Paeonia lactiflora, RNAs from 10 tissues of the Paeonia lactiflora Pall. cv Shaoyou17C were used for de novo assembly, and 416,062 unigenes were obtained. Using a homology search, it was found that 236,222 (approximately 57%) unigenes had at least one BLAST hit in one or more public data resources. The construction of co-expression networks is a feasible means for improving unigene annotation. Using in-house transcriptome data, we obtained a co-expression network covering 95.13% of the unigenes. Then we integrated co-expression network analyses and lipid-related pathway genes to study lipid metabolism in Paeonia lactiflora cultivars. Finally, we constructed the online database HpeNet (http://bioinformatics.cau.edu.cn/HpeNet) to integrate transcriptome data, gene information, the co-expression network, and so forth. The database can also be searched for gene details, gene functions, orthologous matches, and other data. Our online database may help the research community identify functional genes and perform research on Paeonia lactiflora more conveniently. We hope that de novo transcriptome assembly, combined with co-expression networks, can provide a feasible means to predict the gene function of species that do not have a reference genome.

croFGD: Catharanthus roseus Functional Genomics Database

Catharanthus roseus is a medicinal plant, which can produce monoterpene indole alkaloid (MIA) metabolites with biological activity and is rich in vinblastine and vincristine. With release of the scaffolded genome sequence of C. roseus, it is necessary to annotate gene functions on the whole-genome level. Recently, 53 RNA-seq datasets are available in public with different tissues (flower, root, leaf, seedling, and shoot) and different treatments (MeJA, PnWB infection and yeast elicitor). We used in-house data process pipeline with the combination of PCC and MR algorithms to construct a co-expression network exploring multi-dimensional gene expression (global, tissue preferential, and treat response) through multi-layered approaches. In the meanwhile, we added miRNA-target pairs, predicted PPI pairs into the network and provided several tools such as gene set enrichment analysis, functional module enrichment analysis, and motif analysis for functional prediction of the co-expression genes. Finally, we have constructed an online croFGD database (http://bioinformatics.cau.edu.cn/croFGD/). We hope croFGD can help the communities to study the C. roseus functional genomics and make novel discoveries about key genes involved in some important biological processes.

PlantEAR: Functional Analysis Platform for Plant EAR Motif-Containing Proteins

The Ethylene-responsive element binding factor-associated Amphiphilic Repression (EAR) motifs, which were initially identified in members of the Arabidopsis ethylene response factor (ERF) family, are transcriptional repression motifs in plants and are defined by the consensus sequence patterns of either LxLxL or DLNxxP. EAR motif-containing proteins can function as transcription repressors, thus interacting with co-repressors, such as TOPLESS and AtSAP18, affecting the structure of chromatin by histone modifications and thereby repressing gene transcription. EAR motif-containing proteins are highly conserved across diverse plant species and play important roles in hormone signal transduction, stress responses and development, but they have not been identified in most plants. In this study, we identified 20,542 EAR motif-containing proteins from 71 plant species based on a Hidden Markov Model and orthologous gene search, and then we constructed a functional analysis platform for plant EAR motif-containing proteins (PlantEAR, http://structuralbiology.cau.edu.cn/plantEAR) by integrating a variety of functional annotations and processed data. Several tools were provided as functional support for EAR motif-containing proteins, such as browse, search, co-expression and protein-protein interaction (PPI) network analysis as well as cis-element analysis and gene set enrichment analysis (GSEA). In addition, basing on the identified EAR motif-containing proteins, we also explored their distribution in various species and found that the numbers of EAR motif-containing proteins showed an increasing trend in evolution from algae to angiosperms.

Genomic, expressional, protein-protein interactional analysis of Trihelix transcription factor genes in Setaria italia and inference of their evolutionary trajectory

Background

Trihelix transcription factors (TTF) play important roles in plant growth and response to adversity stress. Until now, genome-wide identification and analysis of this gene family in foxtail millet has not been available. Here, we identified TTF genes in the foxtail millet and its grass relatives, and characterized their functional domains.

Results

As to sequence divergence, TTF genes were previously divided into five subfamilies, I-V. We found that Trihelix family members in foxtail millet and other grasses mostly preserved their ancestral chromosomal locations during millions of years’ evolution. Six amino acid sites of the SIP1 subfamily possibly were likely subjected to significant positive selection. Highest expression level was observed in the spica, with the SIP1 subfamily having highest expression level. As to the origination and expansion of the gene family, notably we showed that a subgroup of subfamily IV was the oldest, and therefore was separated to define a new subfamily O. Overtime, starting from the subfamily O, certain genes evolved to form subfamilies III and I, and later from subfamily I to develop subfamilies II and V. The oldest gene, Si1g016284, has the most structural changes, and a high expression in different tissues. What’s more interesting is that it may have bridge the interaction with different proteins.

Conclusions

By performing phylogenetic analysis using non-plant species, notably we showed that a subgroup of subfamily IV was the oldest, and therefore was separated to define a new subfamily O. Starting from the subfamily O, certain genes evolved to form other subfamilies. Our work will contribute to understanding the structural and functional innovation of Trihelix transcription factor, and the evolutionary trajectory.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5051-9) contains supplementary material, which is available to authorized users.

Genome-wide identification, phylogeny and expression analyses of SCARECROW-LIKE(SCL) genes in millet (Setaria italica)

As a member of the GRAS gene family, SCARECROW-LIKE (SCL) genes encode transcriptional regulators that are involved in plant information transmission and signal transduction. In this study, 44 SCL genes including two SCARECROW genes in millet were identified to be distributed on eight chromosomes, except chromosome 6. All the millet genes contain motifs 6-8, indicating that these motifs are conserved during the evolution. SCL genes of millet were divided into eight groups based on the phylogenetic relationship and classification of Arabidopsis SCL genes. Several putative millet orthologous genes in Arabidopsis, maize and rice were identified. High throughput RNA sequencing revealed that the expressions of millet SCL genes in root, stem, leaf, spica, and along leaf gradient varied greatly. Analyses combining the gene expression patterns, gene structures, motif compositions, promoter cis-elements identification, alternative splicing of transcripts and phylogenetic relationship of SCL genes indicate that the these genes may play diverse functions. Functionally characterized SCL genes in maize, rice and Arabidopsis would provide us some clues for future characterization of their homologues in millet. To the best of our knowledge, this is the first study of millet SCL genes at the genome wide level. Our work provides a useful platform for functional analysis of SCL genes in millet, a model crop for C₄ photosynthesis and bioenergy studies.

SesameFG: an integrated database for the functional genomics of sesame

Sesame (Sesamum indicum L.) has high oil content, a small diploid genome and a short growth period, making it an attractive species for genetic studies on oilseed crops. With the advancement of next-generation sequencing technology, genomics and functional genomics research of sesame has developed quickly in the last few years, and large amounts of data have been generated. However, these results are distributed in many different publications, and there is a lack of integration. To promote functional genomics research of sesame, we collected genetic information combined with comprehensive phenotypic information and integrated them in the web-based database named SesameFG. The current version of SesameFG contains phenotypic information on agronomic traits of 705 sesame accessions, de novo assembled genomes of three sesame varieties, massive numbers of identified SNPs, gene expression profiles of five tissues, gene families, candidate genes for the important agronomic traits and genomic-SSR markers. All phenotypic and genotypic information in SesameFG is available for online queries and can be downloaded freely. SesameFG provides useful search functions and data mining tools, including Genome Browser and local BLAST services. SesameFG is freely accessible at http://ncgr.ac.cn/SesameFG/. SesameFG provides valuable resources and tools for functional genomics research and the molecular breeding of sesame.

Co-expression network analyses identify functional modules associated with development and stress response in Gossypium arboreum

Cotton is an economically important crop, essential for the agriculture and textile industries. Through integrating transcriptomic data, we discovered that multi-dimensional co-expression network analysis was powerful for predicting cotton gene functions and functional modules. Here, the recently available transcriptomic data on Gossypium arboreum, including data on multiple growth stages of tissues and stress treatment samples were applied to construct a co-expression network exploring multi-dimensional expression (development and stress) through multi-layered approaches. Based on differential gene expression and network analysis, a fibre development regulatory module of the gene GaKNL1 was found to regulate the second cell wall through repressing the activity of REVOLUTA, and a tissue-selective module of GaJAZ1a was examined in response to water stress. Moreover, comparative genomics analysis of the JAZ1-related regulatory module revealed high conservation across plant species. In addition, 1155 functional modules were identified through integrating the co-expression network, module classification and function enrichment tools, which cover functions such as metabolism, stress responses, and transcriptional regulation. In the end, an online platform was built for network analysis (http://structuralbiology.cau.edu.cn/arboreum), which could help to refine the annotation of cotton gene function and establish a data mining system to identify functional genes or modules with important agronomic traits.

ccNET: Database of co-expression networks with functional modules for diploid and polyploid Gossypium

Plant genera with both diploid and polyploid species are a common evolutionary occurrence. Polyploids, especially allopolyploids such as cotton and wheat, are a great model system for heterosis research. Here, we have integrated genome sequences and transcriptome data of Gossypium species to construct co-expression networks and identified functional modules from different cotton species, including 1155 and 1884 modules in G. arboreum and G. hirsutum, respectively. We overlayed the gene expression results onto the co-expression network. We further provided network comparison analysis for orthologous genes across the diploid and allotetraploid Gossypium. We also constructed miRNA-target networks and predicted PPI networks for both cotton species. Furthermore, we integrated in-house ChIP-seq data of histone modification (H3K4me3) together with cis-element analysis and gene sets enrichment analysis tools for studying possible gene regulatory mechanism in Gossypium species. Finally, we have constructed an online ccNET database (http://structuralbiology.cau.edu.cn/gossypium) for comparative gene functional analyses at a multi-dimensional network and epigenomic level across diploid and polyploid Gossypium species. The ccNET database will be beneficial for community to yield novel insights into gene/module functions during cotton development and stress response, and might be useful for studying conservation and diversity in other polyploid plants, such as T. aestivum and Brassica napus.

Histone modifications facilitate the coexpression of bidirectional promoters in rice

Background

Bidirectional gene pairs are highly abundant and mostly co-regulated in eukaryotic genomes. The structural features of bidirectional promoters (BDPs) have been well studied in yeast, humans and plants. However, the underlying mechanisms responsible for the coexpression of BDPs remain understudied, especially in plants.

Results

Here, we characterized chromatin features associated with rice BDPs. Several unique chromatin features were present in rice BDPs but were missing from unidirectional promoters (UDPs), including overrepresented active histone marks, canonical nucleosomes and underrepresented H3K27me3. In particular, overrepresented active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me2 and H3K36me3) were truly overrepresented in type I BDPs but not in the other two BDPs, based on a Kolmogorov-Smirnov test.

Conclusions

Our analyses indicate that active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me3, H3K9ac and H3K27ac) may coordinate with repressive marks (H3K27me3 and H3K9me1/3) to build a unique chromatin structure that favors the coregulation of bidirectional gene pairs. Thus, our findings help to enhance the understanding of unique epigenetic mechanisms that regulate bidirectional gene pairs and may improve the manipulation of gene pairs for crop bioengineering.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3125-0) contains supplementary material, which is available to authorized users.

SorghumFDB: sorghum functional genomics database with multidimensional network analysis

Sorghum (Sorghum bicolor [L.] Moench) has excellent agronomic traits and biological properties, such as heat and drought-tolerance. It is a C₄ grass and potential bioenergy-producing plant, which makes it an important crop worldwide. With the sorghum genome sequence released, it is essential to establish a sorghum functional genomics data mining platform. We collected genomic data and some functional annotations to construct a sorghum functional genomics database (SorghumFDB). SorghumFDB integrated knowledge of sorghum gene family classifications (transcription regulators/factors, carbohydrate-active enzymes, protein kinases, ubiquitins, cytochrome P450, monolignol biosynthesis related enzymes, R-genes and organelle-genes), detailed gene annotations, miRNA and target gene information, orthologous pairs in the model plants Arabidopsis, rice and maize, gene loci conversions and a genome browser. We further constructed a dynamic network of multidimensional biological relationships, comprised of the co-expression data, protein–protein interactions and miRNA-target pairs. We took effective measures to combine the network, gene set enrichment and motif analyses to determine the key regulators that participate in related metabolic pathways, such as the lignin pathway, which is a major biological process in bioenergy-producing plants.

Database URL:http://structuralbiology.cau.edu.cn/sorghum/index.html.

Gene mapping and functional analysis of the novel leaf color gene SiYGL1 in foxtail millet [Setaria italica (L.) P. Beauv]

Setaria italica and its wild ancestor Setaria viridis are emerging as model systems for genetics and functional genomics research. However, few systematic gene mapping or functional analyses have been reported in these promising C4 models. We herein isolated the yellow-green leaf mutant (siygl1) in S. italica using forward genetics approaches. Map-based cloning revealed that SiYGL1, which is a recessive nuclear gene encoding a magnesium-chelatase D subunit (CHLD), is responsible for the mutant phenotype. A single Phe to Leu amino acid change occurring near the ATPase-conserved domain resulted in decreased chlorophyll (Chl) accumulation and modified chloroplast ultrastructure. However, the mutation enhanced the light-use efficiency of the siygl1 mutant, suggesting that the mutated CHLD protein does not completely lose its original activity, but instead, gains novel features. A transcriptional analysis of Chl a oxygenase revealed that there is a strong negative feedback control of Chl b biosynthesis in S. italica. The SiYGL1 mRNA was expressed in all examined tissues, with higher expression observed in the leaves. Comparison of gene expression profiles in wild-type and siygl1 mutant plants indicated that SiYGL1 regulates a subset of genes involved in photosynthesis (rbcL and LHCB1), thylakoid development (DEG2) and chloroplast signaling (SRP54CP). These results provide information regarding the mutant phenotype at the transcriptional level. This study demonstrated that the genetic material of a Setaria species could be ideal for gene discovery investigations using forward genetics approaches and may help to explain the molecular mechanisms associated with leaf color variation.