Deciphering the microbial and molecular responses of geographically diverse Setaria accessions grown in a nutrient-poor soil.

The microbial and molecular characterization of the ectorhizosphere is an important step towards developing a more complete understanding of how the cultivation of biofuel crops can be undertaken in nutrient poor environments. The ectorhizosphere of Setaria is of particular interest because the plant component of this plant-microbe system is an important agricultural grain crop and a model for biofuel grasses. Importantly, Setaria lends itself to high throughput molecular studies. As such, we have identified important intra- and interspecific microbial and molecular differences in the ectorhizospheres of three geographically distant Setaria italica accessions and their wild ancestor S. viridis. All were grown in a nutrient-poor soil with and without nutrient addition. To assess the contrasting impact of nutrient deficiency observed for two S. italica accessions, we quantitatively evaluated differences in soil organic matter, microbial community, and metabolite profiles. Together, these measurements suggest that rhizosphere priming differs with Setaria accession, which comes from alterations in microbial community abundances, specifically Actinobacteria and Proteobacteria populations. When globally comparing the metabolomic response of Setaria to nutrient addition, plants produced distinctly different metabolic profiles in the leaves and roots. With nutrient addition, increases of nitrogen containing metabolites were significantly higher in plant leaves and roots along with significant increases in tyrosine derived alkaloids, serotonin, and synephrine. Glycerol was also found to be significantly increased in the leaves as well as the ectorhizosphere. These differences provide insight into how C4 grasses adapt to changing nutrient availability in soils or with contrasting fertilization schemas. Gained knowledge could then be utilized in plant enhancement and bioengineering efforts to produce plants with superior traits when grown in nutrient poor soils.

Phylogenetic analysis of the Si7PPO gene in foxtail millet, Setaria italica, provides further evidence for multiple origins of the negative phenol color reaction phenotype.

To elucidate the diversity and evolution of the Si7PPO gene that controls phenol color reaction (Phr) in foxtail millet, Setaria italica, we analyzed sequence polymorphisms of the Si7PPO gene in 39 accessions consisting of foxtail millet landraces (32 accessions) and their wild ancestor ssp. viridis (seven accessions) collected from various regions in Europe and Asia. The accessions included wild type (positive Phr) and three different types of loss-of-function phenotype (negative Phr), "stop codon type", "TE1-insertion type" and "6-bp duplication type", found in our previous study. We constructed a phylogenetic tree of the gene and found that accessions with positive Phr showed higher genetic diversity at the nucleotide sequence level. We also found that the three different loss-of-function types formed different clusters, suggesting that landraces with negative Phr have multiple origins from three different lineages including both landrace and ssp. viridis accessions with positive Phr.

Genomic dissection and expression analysis of stress-responsive genes in C4 panicoid models, Setaria italica and Setaria viridis.

The study reports the identification and expression profiling of five major classes of C4 pathway-specific genes, namely, carbonic anhydrase (CaH), phosphoenolpyruvate carboxylase (PEPC), pyruvate orthophosphate dikinase (PPDK), NADP-dependent malate dehydrogenase (MDH) and NADP-dependent malic enzyme (NADP-ME), in the model species, Setaria italica and Setaria viridis. A total of 42 and 41 genes were identified in S. italica and S. viridis, respectively. Further analysis revealed that segmental and tandem duplications have contributed to the expansion of these gene families. RNA-Seq derived expression profiles of the gene family members showed their differential expression pattern in tissues and dehydration stress. Comparative genome mapping and Ks dating provided insights into their duplication and divergence in the course of evolution. Expression profiling of candidate genes in contrasting S. italica cultivars subjected to abiotic stresses and hormone treatments showed distinct stress-specific upregulation of SiαCaH1, SißCaH5, SiPEPC2, SiPPDK2, SiMDH8, and SiNADP-ME5 in the tolerant cultivar. Overexpression of SiNADP-ME5 in heterologous yeast system enabled the transgenic cells to survive and grow in dehydration stress conditions, which highlights the putative role of SiNADP-ME5 in conferring tolerance to dehydration stress. Altogether, the study highlights key genes that could be potential candidates for elucidating their functional roles in abiotic stress response.

The population genetic structure approach adds new insights into the evolution of plant LTR retrotransposon lineages.

Long terminal repeat retrotransposons (LTR-RTs) in plant genomes differ in abundance, structure and genomic distribution, reflecting the large number of evolutionary lineages. Elements within lineages can be considered populations, in which each element is an individual in its genomic environment. In this way, it would be reasonable to apply microevolutionary analyses to understand transposable element (TE) evolution, such as those used to study the genetic structure of natural populations. Here, we applied a Bayesian method to infer genetic structure of populations together with classical phylogenetic and dating tools to analyze LTR-RT evolution using the monocot Setaria italica as a model species. In contrast to a phylogeny, the Bayesian clusterization method identifies populations by assigning individuals to one or more clusters according to the most probabilistic scenario of admixture, based on genetic diversity patterns. In this work, each LTR-RT insertion was considered to be one individual and each LTR-RT lineage was considered to be a single species. Nine evolutionary lineages of LTR-RTs were identified in the S. italica genome that had different genetic structures with variable numbers of clusters and levels of admixture. Comprehensive analysis of the phylogenetic, clusterization and time of insertion data allowed us to hypothesize that admixed elements represent sequences that harbor ancestral polymorphic sequence signatures. In conclusion, application of microevolutionary concepts in genome evolution studies is suitable as a complementary approach to phylogenetic analyses to address the evolutionary history and functional features of TEs.

Genome-Wide Association Study of Major Agronomic Traits in Foxtail Millet (Setaria italica L.) Using ddRAD Sequencing.

Foxtail millet (Setaria italica), the second largest cultivated millet crop after pearl millet, is utilized for food and forage globally. Further, it is also considered as a model crop for studying agronomic, nutritional and biofuel traits. In the present study, a genome-wide association study (GWAS) was performed for ten important agronomic traits in 142 foxtail millet core eco-geographically diverse genotypes using 10 K SNPs developed through GBS-ddRAD approach. Number of SNPs on individual chromosome ranged from 844 (chromosome 5) to 2153 (chromosome 8) with an average SNP frequency of 25.9 per Mb. The pairwise linkage disequilibrium (LD) estimated using the squared-allele frequency correlations was found to decay rapidly with the genetic distance of 177 Kb. However, for individual chromosome, LD decay distance ranged from 76 Kb (chromosome 6) to 357 Kb (chromosome 4). GWAS identified 81 MTAs (marker-trait associations) for ten traits across the genome. High confidence MTAs for three important agronomic traits including FLW (flag leaf width), GY (grain yield) and TGW (thousand-grain weight) were identified. Significant pyramiding effect of identified MTAs further supplemented its importance in breeding programs. Desirable alleles and superior genotypes identified in the present study may prove valuable for foxtail millet improvement through marker-assisted selection.

[MYB-like transcription factor SiMYB42 from foxtail millet (Setaria italica L.) enhances Arabidopsis tolerance to low-nitrogen stress].

Myeloblastosis (MYB) transcription factors are one of the largest families of transcription factors in higher plants. They play an important role in plant development, defense response processes, and non-biological stresses, i.e., drought stress. Foxtail millet (Setaria italica L.), originated in China, is resistant to drought and low nutrition stresses and has been regarded as an ideal material for studying abiotic stress resistance in monocotyledon. In this study, we ran a transcription profile analysis of zheng 204 under low-nitrogen conditions and identified a MYB-like transcription factor SiMYB42, which was up-regulated under low-nitrogen stress. Phylogenetic tree analysis showed that SiMYB42 belongs to R2R3-MYB subfamily and has two MYB conserved domains. Expression pattern analysis showed that SiMYB42 was significantly up-regulated under various stress conditions, including low-nitrogen stress, high salt, drought and ABA conditions. The results of subcellular localization, quantitative real-time PCR and transcriptional activation analysis indicated that SiMYB42 protein localizes to the nucleus and cell membrane of plant cells, mainly expressed in the leaf or root of foxtail millet, and has transcription activation activity. Functional analysis showed that there was no significant difference between transgenic SiMYB42 Arabidopsis and wild-type (WT) Arabidopsis under normal conditions; however, under low-nitrogen condition, the root length, surface area and seedling fresh weight in transgenic SiMYB42 Arabidopsis, were significantly higher than their counterparts in WT. These results suggest that SiMYB42 transgenic plants exhibit higher tolerance to low-nitrogen stress. Expression levels of nitrate transporters genes NRT2.1, NRT2.4 and NRT2.5, which are the transcriptional targets of SiMYB42, were higher in transgenic SiMYB42 Arabidopsis plants than those in WT; the promoter regions of NRT2.1, NRT2.4 and NRT2.5 all have MYB binding sites. These results indicate that SiMYB42 might enhance foxtail millet tolerance to low-nitrogen condition through regulating the expression of nitrate transporter genes. This study reveals the possible functions of SiMYB42 in a low-nitrogen stress response pathway, and provides a foundation for further understanding the entire regulation network of foxtail millet in response to low-nitrogen stress.

Genome-wide analysis of autophagy-associated genes in foxtail millet (Setaria italica L.) and characterization of the function of SiATG8a in conferring tolerance to nitrogen starvation in rice.

BACKGROUND: Autophagy is a cellular degradation process that is highly evolutionarily-conserved in yeast, plants, and animals. In plants, autophagy plays important roles in regulating intracellular degradation and recycling of amino acids in response to nutrient starvation, senescence, and other environmental stresses. Foxtail millet (Setaria italica) has strong resistance to stresses and has been proposed as an ideal material for use in the study of the physiological mechanisms of abiotic stress tolerance in plants. Although the genome sequence of foxtail millet (Setaria italica) is available, the characteristics and functions of abiotic stress-related genes remain largely unknown for this species. RESULTS: A total of 37 putative ATG (autophagy-associated genes) genes in the foxtail millet genome were identified. Gene duplication analysis revealed that both segmental and tandem duplication events have played significant roles in the expansion of the ATG gene family in foxtail millet. Comparative synteny mapping between the genomes of foxtail millet and rice suggested that the ATG genes in both species have common ancestors, as their ATG genes were primarily located in similar syntenic regions. Gene expression analysis revealed the induced expression of 31 SiATG genes by one or more phytohormone treatments, 26 SiATG genes by drought, salt and cold, 24 SiATG genes by darkness and 25 SiATG genes by nitrogen starvation. Results of qRT-PCR showing that among 37 SiATG genes, the expression level of SiATG8a was the highest after nitrogen starvation treatment 24 h, suggesting its potential role in tolerance to nutrient starvation. Moreover, the heterologous expression of SiATG8a in rice improved nitrogen starvation tolerance. Compared to wild type rice, the transgenic rice performed better and had higher aboveground total nitrogen content when the plants were grown under nitrogen starvation conditions. CONCLUSIONS: Our results deepen understanding about the characteristics and functions of ATG genes in foxtail millet and also identify promising new genetic resources that should be of use in future efforts to develop varieties of foxtail millet and other crop species that have resistance to nitrogen deficiency stress.

Structure, organization and evolution of ADP-ribosylation factors in rice and foxtail millet, and their expression in rice.

ADP-ribosylation factors (ARFs) have been reported to function in diverse physiological and molecular activities. Recent evidences also demonstrate the involvement of ARFs in conferring tolerance to biotic and abiotic stresses in plant species. In the present study, 23 and 25 ARF proteins were identified in C3 model- rice and C4 model- foxtail millet, respectively. These proteins are classified into four classes (I-IV) based on phylogenetic analysis, with ARFs in classes I-III and ARF-like proteins (ARLs) in class IV. Sequence alignment and domain analysis revealed the presence of conserved and additional motifs, which may contribute to neo- and sub-functionalization of these proteins. Promoter analysis showed the presence of several cis-regulatory elements related to stress and hormone response, indicating their role in stress regulatory network. Expression analysis of rice ARFs and ARLs in different tissues, stresses and abscisic acid treatment highlighted temporal and spatial diversification of gene expression. Five rice cultivars screened for allelic variations in OsARF genes showed the presence of allelic polymorphisms in few gene loci. Altogether, the study provides insights on characteristics of ARF/ARL genes in rice and foxtail millet, which could be deployed for further functional analysis to extrapolate their precise roles in abiotic stress responses.

Complete chloroplast genome sequence of green foxtail (Setaria viridis), a promising model system for C4 photosynthesis.

The complete chloroplast genome of green foxtail (Setaria viridis), a promising model system for C4 photosynthesis, is first reported in this study. The genome harbors a large single copy (LSC) region of 81 016 bp and a small single copy (SSC) region of 12 456 bp separated by a pair of inverted repeat (IRa and IRb) regions of 22 315 bp. GC content is 38.92%. The proportion of coding sequence is 57.97%, comprising of 111 (19 duplicated in IR regions) unique genes, 71 of which are protein-coding genes, four are rRNA genes, and 36 are tRNA genes. Phylogenetic analysis indicated that S. viridis was clustered with its cultivated species S. italica in the tribe Paniceae of the family Poaceae. This newly determined chloroplast genome will provide valuable genetic resources to assist future studies on C4 photosynthesis in grasses.

Multiple origins of the phenol reaction negative phenotype in foxtail millet, Setaria italica (L.) P. Beauv., were caused by independent loss-of-function mutations of the polyphenol oxidase (Si7PPO) gene during domestication.

Foxtail millet shows variation in positive phenol color reaction (Phr) and negative Phr in grains, but predominant accessions of this crop are negative reaction type, and the molecular genetic basis of the Phr reaction remains unresolved. In this article, we isolated polyphenol oxidase (PPO) gene responsible for Phr using genome sequence information and investigated molecular genetic basis of negative Phr and crop evolution of foxtail millet. First of all, we searched for PPO gene homologs in a foxtail millet genome database using a rice PPO gene as a query and successfully found three copies of the PPO gene. One of the PPO gene homologs on chromosome 7 showed the highest similarity with PPO genes expressed in hulls (grains) of other cereal species including rice, wheat, and barley and was designated as Si7PPO. Phr phenotypes and Si7PPO genotypes completely co-segregated in a segregating population. We also analyzed the genetic variation conferring negative Phr reaction. Of 480 accessions of the landraces investigated, 87 (18.1 %) showed positive Phr and 393 (81.9 %) showed negative Phr. In the 393 Phr negative accessions, three types of loss-of-function Si7PPO gene were predominant and independently found in various locations. One of them has an SNP in exon 1 resulting in a premature stop codon and was designated as stop codon type, another has an insertion of a transposon (Si7PPO-TE1) in intron 2 and was designated as TE1-insertion type, and the other has a 6-bp duplication in exon 3 resulting in the duplication of 2 amino acids and was designated as 6-bp duplication type. As a rare variant of the stop codon type, one accession additionally has an insertion of a transposon, Si7PPO-TE2, in intron 2 and was designated as "stop codon +TE2 insertion type". The geographical distribution of accessions with positive Phr and those with three major types of negative Phr was also investigated. Accessions with positive Phr were found in subtropical and tropical regions at frequencies of ca. 25-67 % and those with negative Phr were broadly found in Europe and Asia. The stop codon type was found in 285 accessions and was broadly distributed in Europe and Asia, whereas the TE-1 insertion type was found in 99 accessions from Europe and Asia but was not found in India. The 6-bp duplication type was found in only 8 accessions from Nansei Islands (Okinawa Prefecture) of Japan. We also analyzed Phr in the wild ancestor and concluded that the negative Phr type was likely to have originated after domestication of foxtail millet. It was also implied that negative Phr of foxtail millet arose by multiple independent loss of function of PPO gene through dispersal because of some advantages under some environmental conditions and human selection as in rice and barley.

Polyamine metabolism influences antioxidant defense mechanism in foxtail millet (Setaria italica L.) cultivars with different salinity tolerance.

KEY MESSAGE: Polyamines can regulate the expression of antioxidant enzymes and impart plants tolerance to abiotic stresses. A comparative analysis of polyamines, their biosynthetic enzymes at kinetic and at transcriptional level, and their role in regulating the induction of antioxidant defense enzymes under salt stress condition in two foxtail millet (Setaria italica L.) cultivars, namely Prasad, a salt-tolerant, and Lepakshi, a salt-sensitive cultivar was conducted. Salt stress resulted in elevation of free polyamines due to increase in the activity of spermidine synthase and S-adenosyl methionine decarboxylase enzymes in cultivar Prasad compared to cultivar Lepakshi under different levels of NaCl stress. These enzyme activities were further confirmed at the transcript level via qRT-PCR analysis. The cultivar Prasad showed a greater decrease in diamine oxidase and polyamine oxidase activity, which results in the accumulation of polyamine pools over cultivar Lepakshi. Generation of free radicals, such as O 2 (·-) and H2O2, was also analyzed quantitatively. A significant increase in O 2 (·-) and H2O2 in the cultivar Lepakshi compared with cultivar Prasad was recorded in overall pool sizes. Further, histochemical staining showed lesser accumulation of O 2 (·-) and of H2O2 in the leaves of cultivar Prasad than cultivar Lepakshi. Our results also suggest the ability of polyamine oxidation in regulating the induction of antioxidative defense enzymes, which involve in the elimination of toxic levels of O 2 (·-) and H2O2, such as Mn-superoxide dismutase, catalase and ascorbate peroxidase. The contribution of polyamines in modulating antioxidative defense mechanism in NaCl stress tolerance is discussed.

Bioinformatic analysis of the nucleotide binding site-encoding disease-resistance genes in foxtail millet (Setaria italica (L.) Beauv.).

The nucleotide-binding site (NBS) disease-resistance genes are the largest category of plant disease-resistance gene analogs. The complete set of disease-resistant candidate genes, which encode the NBS sequence, was filtered in the genomes of two varieties of foxtail millet (Yugu1 and 'Zhang gu'). This study investigated a number of characteristics of the putative NBS genes, such as structural diversity and phylogenetic relationships. A total of 269 and 281 NBS-coding sequences were identified in Yugu1 and 'Zhang gu', respectively. When the two databases were compared, 72 genes were found to be identical and 164 genes showed more than 90% similarity. Physical positioning and gene family analysis of the NBS disease-resistance genes in the genome revealed that the number of genes on each chromosome was similar in both varieties. The eighth chromosome contained the largest number of genes and the ninth chromosome contained the lowest number of genes. Exactly 34 gene clusters containing the 161 genes were found in the Yugu1 genome, with each cluster containing 4.7 genes on average. In comparison, the 'Zhang gu' genome possessed 28 gene clusters, which had 151 genes, with an average of 5.4 genes in each cluster. The largest gene cluster, located on the eighth chromosome, contained 12 genes in the Yugu1 database, whereas it contained 16 genes in the 'Zhang gu' database. The classification results showed that the CC-NBS-LRR gene made up the largest part of each chromosome in the two databases. Two TIR-NBS genes were also found in the Yugu1 genome.

Population genetics of Setaria viridis, a new model system.

An extensive survey of the standing genetic variation in natural populations is among the priority steps in developing a species into a model system. In recent years, green foxtail (Setaria viridis), along with its domesticated form foxtail millet (S. italica), has rapidly become a promising new model system for C4 grasses and bioenergy crops, due to its rapid life cycle, large amount of seed production and small diploid genome, among other characters. However, remarkably little is known about the genetic diversity in natural populations of this species. In this study, we survey the genetic diversity of a worldwide sample of more than 200 S. viridis accessions, using the genotyping-by-sequencing technique. Two distinct genetic groups in S. viridis and a third group resembling S. italica were identified, with considerable admixture among the three groups. We find the genetic variation of North American S. viridis correlates with both geography and climate and is representative of the total genetic diversity in this species. This pattern may reflect several introduction/dispersal events of S. viridis into North America. We also modelled demographic history and show signal of recent population decline in one subgroup. Finally, we show linkage disequilibrium decay is rapid (<45 kb) in our total sample and slow in genetic subgroups. These results together provide an in-depth understanding of the pattern of genetic diversity of this new model species on a broad geographic scale. They also provide key guidelines for on-going and future work including germplasm preservation, local adaptation, crossing designs and genomewide association studies.

Development of eSSR-Markers in Setaria italica and Their Applicability in Studying Genetic Diversity, Cross-Transferability and Comparative Mapping in Millet and Non-Millet Species.

Foxtail millet (Setariaitalica L.) is a tractable experimental model crop for studying functional genomics of millets and bioenergy grasses. But the limited availability of genomic resources, particularly expressed sequence-based genic markers is significantly impeding its genetic improvement. Considering this, we attempted to develop EST-derived-SSR (eSSR) markers and utilize them in germplasm characterization, cross-genera transferability and in silico comparative mapping. From 66,027 foxtail millet EST sequences 24,828 non-redundant ESTs were deduced, representing ~16 Mb, which revealed 534 (~2%) eSSRs in 495 SSR containing ESTs at a frequency of 1/30 kb. A total of 447 pp were successfully designed, of which 327 were mapped physically onto nine chromosomes. About 106 selected primer pairs representing the foxtail millet genome showed high-level of cross-genera amplification at an average of ~88% in eight millets and four non-millet species. Broad range of genetic diversity (0.02-0.65) obtained in constructed phylogenetic tree using 40 eSSR markers demonstrated its utility in germplasm characterizations and phylogenetics. Comparative mapping of physically mapped eSSR markers showed considerable proportion of sequence-based orthology and syntenic relationship between foxtail millet chromosomes and sorghum (~68%), maize (~61%) and rice (~42%) chromosomes. Synteny analysis of eSSRs of foxtail millet, rice, maize and sorghum suggested the nested chromosome fusion frequently observed in grass genomes. Thus, for the first time we had generated large-scale eSSR markers in foxtail millet and demonstrated their utility in germplasm characterization, transferability, phylogenetics and comparative mapping studies in millets and bioenergy grass species.

Solvent optimization extraction of antioxidants from foxtail millet species' insoluble fibers and their free radical scavenging properties.

In this study, water and 80% of four organic solvents were employed to optimize the extraction of antioxidants from two species of foxtail millet's insoluble fibers under the same temperature, time, and solid/solvent ratio. The results showed that the acetone was able to extract the maximum amount of antioxidants (2.32 mg/g fiber for white specie and 3.86 mg/g fiber for yellow specie) followed by methanol and propanol from both samples. The neutral and the ethanol on the other hand extracted small amount of the antioxidants from the two fiber materials. While considerable level of Total Polyphenols Content (TPC) was recorded in both the water and the organic solvents' extracts, only traces of Total Flavonoid content (TFC) were observed in water, methanol and ethanol extracts. Propanol and acetone extracts was negative to the TFC test. The potency of both white and yellow foxtail millets' insoluble fibers antioxidant extracts was investigated using five different in vitro tests. It was realized that there was a variation in their capacities to quench DPPH and ABTS(+) radicals for the time running of 0-60 min. The samples from the yellow cereal exhibited high inhibition capacity against ABTS(+). No correlation was observed between TPC and radical scavenging capacities for DPPH and ABTS(+). In general, the yellow species contained more antioxidants in comparison with the white one and this accounted for its high antioxidant activity.

Novel genomes and genome constitutions identified by GISH and 5S rDNA and knotted1 genomic sequences in the genus Setaria.

BACKGROUND: The Setaria genus is increasingly of interest to researchers, as its two species, S. viridis and S. italica, are being developed as models for understanding C4 photosynthesis and plant functional genomics. The genome constitution of Setaria species has been studied in the diploid species S. viridis, S. adhaerans and S. grisebachii, where three genomes A, B and C were identified respectively. Two allotetraploid species, S. verticillata and S. faberi, were found to have AABB genomes, and one autotetraploid species, S. queenslandica, with an AAAA genome, has also been identified. The genomes and genome constitutions of most other species remain unknown, even though it was thought there are approximately 125 species in the genus distributed world-wide. RESULTS: GISH was performed to detect the genome constitutions of Eurasia species of S. glauca, S. plicata, and S. arenaria, with the known A, B and C genomes as probes. No or very poor hybridization signal was detected indicating that their genomes are different from those already described. GISH was also performed reciprocally between S. glauca, S. plicata, and S. arenaria genomes, but no hybridization signals between each other were found. The two sets of chromosomes of S. lachnea both hybridized strong signals with only the known C genome of S. grisebachii. Chromosomes of Qing 9, an accession formerly considered as S. viridis, hybridized strong signal only to B genome of S. adherans. Phylogenetic trees constructed with 5S rDNA and knotted1 markers, clearly classify the samples in this study into six clusters, matching the GISH results, and suggesting that the F genome of S. arenaria is basal in the genus. CONCLUSIONS: Three novel genomes in the Setaria genus were identified and designated as genome D (S. glauca), E (S. plicata) and F (S. arenaria) respectively. The genome constitution of tetraploid S. lachnea is putatively CCC'C'. Qing 9 is a B genome species indigenous to China and is hypothesized to be a newly identified species. The difference in genome constitution and origin of S. verticillata and S. faberi is also discussed. The new genomes and the genome constitutions of Setaria species identified in this report provide useful information for Setaria germplasm management, foxtail millet breeding, grass evolution and the development of S. viridis and S. italica as a new model for functional genomics.

Genome-wide development and use of microsatellite markers for large-scale genotyping applications in foxtail millet [Setaria italica (L.)].

The availability of well-validated informative co-dominant microsatellite markers and saturated genetic linkage map has been limited in foxtail millet (Setaria italica L.). In view of this, we conducted a genome-wide analysis and identified 28 342 microsatellite repeat-motifs spanning 405.3 Mb of foxtail millet genome. The trinucleotide repeats (∼48%) was prevalent when compared with dinucleotide repeats (∼46%). Of the 28 342 microsatellites, 21 294 (∼75%) primer pairs were successfully designed, and a total of 15 573 markers were physically mapped on 9 chromosomes of foxtail millet. About 159 markers were validated successfully in 8 accessions of Setaria sp. with ∼67% polymorphic potential. The high percentage (89.3%) of cross-genera transferability across millet and non-millet species with higher transferability percentage in bioenergy grasses (∼79%, Switchgrass and ∼93%, Pearl millet) signifies their importance in studying the bioenergy grasses. In silico comparative mapping of 15 573 foxtail millet microsatellite markers against the mapping data of sorghum (16.9%), maize (14.5%) and rice (6.4%) indicated syntenic relationships among the chromosomes of foxtail millet and target species. The results, thus, demonstrate the immense applicability of developed microsatellite markers in germplasm characterization, phylogenetics, construction of genetic linkage map for gene/quantitative trait loci discovery, comparative mapping in foxtail millet, including other millets and bioenergy grass species.

Development and utilization of novel intron length polymorphic markers in foxtail millet (Setaria italica (L.) P. Beauv.).

Introns are noncoding sequences in a gene that are transcribed to precursor mRNA but spliced out during mRNA maturation and are abundant in eukaryotic genomes. The availability of codominant molecular markers and saturated genetic linkage maps have been limited in foxtail millet (Setaria italica (L.) P. Beauv.). Here, we describe the development of 98 novel intron length polymorphic (ILP) markers in foxtail millet using sequence information of the model plant rice. A total of 575 nonredundant expressed sequence tag (EST) sequences were obtained, of which 327 and 248 unique sequences were from dehydration- and salinity-stressed suppression subtractive hybridization libraries, respectively. The BLAST analysis of 98 EST sequences suggests a nearly defined function for about 64% of them, and they were grouped into 11 different functional categories. All 98 ILP primer pairs showed a high level of cross-species amplification in two millets and two nonmillets species ranging from 90% to 100%, with a mean of ∼97%. The mean observed heterozygosity and Nei's average gene diversity 0.016 and 0.171, respectively, established the efficiency of the ILP markers for distinguishing the foxtail millet accessions. Based on 26 ILP markers, a reasonable dendrogram of 45 foxtail millet accessions was constructed, demonstrating the utility of ILP markers in germplasm characterizations and genomic relationships in millets and nonmillets species.

Genetic structure of landraces in foxtail millet (Setaria italica (L.) P. Beauv.) revealed with transposon display and interpretation to crop evolution of foxtail millet.

Although the origin and domestication process of foxtail millet (Setaria italica subsp. italica (L.) P. Beauv.) has been studied by several groups, the issue is still ambiguous. It is essential to resolve this issue by studying a large number of accessions with sufficient markers covering the entire genome. Genetic structures were analyzed by transposon display (TD) using 425 accessions of foxtail millet and 12 of the wild ancestor green foxtail (Setaria italica subsp. viridis (L.) P. Beauv.). We used three recently active transposons (TSI-1, TSI-7, and TSI-10) as genome-wide markers and succeeded in demonstrating geographical structures of the foxtail millet. A neighbor-joining dendrogram based on TD grouped the foxtail millet accessions into eight major clusters, each of which consisted of accessions collected from adjacent geographical areas. Eleven out of 12 green foxtail accessions were grouped separately from the clusters of foxtail millet. These results indicated strong regional differentiations and a long history of cultivation in each region. Furthermore, we discuss the relationship between foxtail millet and green foxtail and suggest a monophyletic origin of foxtail millet domestication.

Population genetics of foxtail millet and its wild ancestor.

BACKGROUND: Foxtail millet (Setaria italica (L.) P. Beauv.), one of the most ancient domesticated crops, is becoming a model system for studying biofuel crops and comparative genomics in the grasses. However, knowledge on the level of genetic diversity and linkage disequilibrium (LD) is very limited in this crop and its wild ancestor, green foxtail (Setaria viridis (L.) P. Beauv.). Such information would help us to understand the domestication process of cultivated species and will allow further research in these species, including association mapping and identification of agricultural significant genes involved in domestication. RESULTS: In this study, we surveyed DNA sequence for nine loci across 50 accessions of cultivated foxtail millet and 34 of its wild progenitor. We found a low level of genetic diversity in wild green foxtail (θ = 0.0059), θ means Watterson's estimator of θ. Despite of a 55% loss of its wild diversity, foxtail millet still harbored a considerable level of diversity (θ = 0.0027) when compared to rice and sorghum (θ = 0.0024 and 0.0034, respectively). The level of LD in the domesticated foxtail millet extends to 1 kb, while it decayed rapidly to a negligible level within 150 bp in wild green foxtail. Using coalescent simulation, we estimated the bottleneck severity at k = 0.6095 when ρ/θ = 1. These results indicated that the domestication bottleneck of foxtail millet was more severe than that of maize but slightly less pronounced than that of rice. CONCLUSIONS: The results in this study establish a general framework for the domestication history of foxtail millet. The low level of genetic diversity and the increased level of LD in foxtail millet are mainly caused by a population bottleneck, although gene flow from foxtail millet to green foxtail is another factor that may have shaped the pattern of genetic diversity of these two related gene pools. The knowledge provided in this study will benefit future population based studies in foxtail millet.