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 ( Setaria italica 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.     

Genome-wide development of transposable elements-based markers in foxtail millet and construction of an integrated database

Transposable elements (TEs) are major components of plant genome and are reported to play significant roles in functional genome diversity and phenotypic variations. Several TEs are highly polymorphic for insert location in the genome and this facilitates development of TE-based markers for various genotyping purposes. Considering this, a genome-wide analysis was performed in the model plant foxtail millet. A total of 30,706 TEs were identified and classified as DNA transposons (24,386), full-length Copia type (1,038), partial or solo Copia type (10,118), full-length Gypsy type (1,570), partial or solo Gypsy type (23,293) and Long- and Short-Interspersed Nuclear Elements (3,659 and 53, respectively). Further, 20,278 TE-based markers were developed, namely Retrotransposon-Based Insertion Polymorphisms (4,801, ∼24%), Inter-Retrotransposon Amplified Polymorphisms (3,239, ∼16%), Repeat Junction Markers (4,451, ∼22%), Repeat Junction-Junction Markers (329, ∼2%), Insertion-Site-Based Polymorphisms (7,401, ∼36%) and Retrotransposon-Microsatellite Amplified Polymorphisms (57, 0.2%). A total of 134 Repeat Junction Markers were screened in 96 accessions of Setaria italica and 3 wild Setaria accessions of which 30 showed polymorphism. Moreover, an open access database for these developed resources was constructed (Foxtail millet Transposable Elements-based Marker Database; http://59.163.192.83/ltrdb/index.html). Taken together, this study would serve as a valuable resource for large-scale genotyping applications in foxtail millet and related 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.     

Comparative genetic maps of foxtail millet (Setaria italica) and rice (Oryza sativa)

 A foxtail millet-rice comparative genetic map was constructed using mapped rice RFLP markers and wheat genomic and cDNA clones with known map position in rice. About 74% and 37% of the cDNA and genomic clones, respectively, were transferable to foxtail millet, confirming that conservation at the DNA level is greatest in genic regions. A high degree of conserved colinearity was observed between the two genomes. Five entire foxtail millet chromosomes appear to be colinear with five entire rice chromosomes. The remaining four foxtail millet linkage groups each show colinearity with segments of two rice chromosomes. The rearrangements of rice chromosomes 3 and 10 to form foxtail millet chromosome IX, and 7 and 9 to form chromosome II are very similar to those required to form maize chromosomes 1 and 7 and sorghum linkage groups C and B, indicating Setaria’s clear taxonomic position within the subfamily of the Panicoideae. Received: 18 December 1996 / Accepted: 4 August 1997     

Phytolith analysis for differentiating between foxtail millet (Setaria italica) and green foxtail (Setaria viridis)

Foxtail millet (Setaria italica) is one of the oldest domesticated cereal crops in Eurasia, but identifying foxtail millets, especially in charred grains, and differentiating it from its wild ancestor, green foxtail (Setaria viridis), in the archaeobotanical remains, is still problematic. Phytolithic analysis provides a meaningful method for identifying this important crop. In this paper, the silicon structure patterns in the glumes, lemmas, and paleas from inflorescence bracts in 16 modern plants of foxtail millet and green foxtail from China and Europe are examined using light microscopy with phase-contrast and a microscopic interferometer. Our research shows that the silicon structure of ΩIII from upper lemmas and paleas in foxtail millet and green foxtail can be correspondingly divided into two groups. The size of ΩIII type phytolith of foxtail millet is bigger than that from green foxtail. Discriminant function analysis reveals that 78.4% of data on foxtail millet and 76.9% of data on green foxtail are correctly classified. This means certain morphotypes of phytoliths are relatively reliable tools for distinguishing foxtail millet from green foxtail. Our results also revealed that the husk phytolith morphologies of foxtail millets from China and Eastern Europe are markedly different from those from Western Europe. Our research gives a meaningful method of separating foxtail millet and green foxtail. The implications of these findings for understanding the history of foxtail millet domestication and cultivation in ancient civilizations are significant.     

Construction of a foxtail millet linkage map and mapping of spikelet-tipped bristles 1(stb1) by using transposon display markers and simple sequence repeat markers with genome sequence information

We attempted genetic analysis and mapping of a gene responsible for the trait “spikelet-tipped bristles” (stb) in foxtail millet, Setaria italica (L.) P.Beauv., as the first step in positional cloning of the gene. This trait is important not only in grain yield such as grain number per panicle of this millet but also in the evolutionary development of the “bristle grass” clade including genera Setaria, Pennisetum and Cenchrus in subfamily Panicoideae. First of all, we confirmed that this trait is controlled by a single recessive gene, using two populations of F2 plants; one was a cross combination between two Taiwanese landraces and the other was a combination between a Taiwanese landrace and a Japanese landrace. Using the latter of the two F2 populations, with transposon display (TD) markers and simple sequence repeat (SSR) markers developed previously, we constructed a genetic map with 13 linkage groups and mapped the responsible gene (stb1) on chromosome 2. We also developed novel SSR markers by using foxtail millet genome sequence information, and we finally constructed nine linkage groups corresponding to nine chromosomes with a total length of 1287.5 cM, and mapped stb1 more precisely on chromosome 2. This work suggests that the foxtail millet genome sequences recently published are useful for developing genome-wide SSR markers for constructing linkage maps and mapping genes in this millet.     

Static and Dynamic Factors Limit Chromosomal Replication Complexity in Escherichia coli,Avoiding Dangers of Runaway Overreplication

We define chromosomal replication complexity (CRC) as the ratio of the copy number of the most replicated regions to that of unreplicated regions on the same chromosome. Although a typical CRC of eukaryotic or bacterial chromosomes is 2, rapidly growing Escherichia coli cells induce an extra round of replication in their chromosomes (CRC = 4). There are also E. coli mutants with stable CRC∼6. We have investigated the limits and consequences of elevated CRC in E. coli and found three limits: the “natural” CRC limit of ∼8 (cells divide more slowly); the “functional” CRC limit of ∼22 (cells divide extremely slowly); and the “tolerance” CRC limit of ∼64 (cells stop dividing). While the natural limit is likely maintained by the eclipse system spacing replication initiations, the functional limit might reflect the capacity of the chromosome segregation system, rather than dedicated mechanisms, and the tolerance limit may result from titration of limiting replication factors. Whereas recombinational repair is beneficial for cells at the natural and functional CRC limits, we show that it becomes detrimental at the tolerance CRC limit, suggesting recombinational misrepair during the runaway overreplication and giving a rationale for avoidance of the latter.     

FmMDb: A Versatile Database of Foxtail Millet Markers for Millets and Bioenergy Grasses Research

The prominent attributes of foxtail millet (Setaria italica L.) including its small genome size, short life cycle, inbreeding nature, and phylogenetic proximity to various biofuel crops have made this crop an excellent model system to investigate various aspects of architectural, evolutionary and physiological significances in Panicoid bioenergy grasses. After release of its whole genome sequence, large-scale genomic resources in terms of molecular markers were generated for the improvement of both foxtail millet and its related species. Hence it is now essential to congregate, curate and make available these genomic resources for the benefit of researchers and breeders working towards crop improvement. In view of this, we have constructed the Foxtail millet Marker Database (FmMDb; http://www.nipgr.res.in/foxtail.html), a comprehensive online database for information retrieval, visualization and management of large-scale marker datasets with unrestricted public access. FmMDb is the first database which provides complete marker information to the plant science community attempting to produce elite cultivars of millet and bioenergy grass species, thus addressing global food insecurity.     

Genetic and Genomic Resources of Small Millets

Small millets are very promising agricultural entity to ensure global food security. They gained remarkable importance in agriculture due to their resilience to climatic changes and increasing demand for nutritious food and feed. The genetic variability in the core and mini-core germplasm of small millets was characterized for nutritional composition and capacity to tolerate abiotic stresses that can be infused in breeding programs. Other than the foxtail millet, availability of genomic information in small millets is far below the mark for use in marker-assisted breeding and other genetic improvement programs. The genome sequence of foxtail millet has recently triggered a plethora of post-genomic analysis and envisaged foxtail millet as a model organism for the C₄ grasses and bioenergy research. Recent developments in the next-generation sequencing technologies enabled us, with the simultaneous discovery of high-throughput markers and multiplexed genotyping of germplasm, to speedup marker-assisted breeding. In this context, an in-depth analysis of the wealth of diverse germplasm resources and future perspectives of integrating genomics in genome-wide marker-trait association and breeding in small millets is worthy.     

Genome-wide analysis of tandem repeats in Tribolium castaneum genome reveals abundant and highly dynamic tandem repeat families with satellite DNA features in euchromatic chromosomal arms

Although satellite DNAs are well-explored components of heterochromatin and centromeres, little is known about emergence, dispersal and possible impact of comparably structured tandem repeats (TRs) on the genome-wide scale. Our bioinformatics analysis of assembled Tribolium castaneum genome disclosed significant contribution of TRs in euchromatic chromosomal arms and clear predominance of satellite DNA-typical 170 bp monomers in arrays of ≥5 repeats. By applying different experimental approaches, we revealed that the nine most prominent TR families Cast1–Cast9 extracted from the assembly comprise ∼4.3% of the entire genome and reside almost exclusively in euchromatic regions. Among them, seven families that build ∼3.9% of the genome are based on ∼170 and ∼340 bp long monomers. Results of phylogenetic analyses of 2500 monomers originating from these families show high-sequence dynamics, evident by extensive exchanges between arrays on non-homologous chromosomes. In addition, our analysis shows that concerted evolution acts more efficiently on longer than on shorter arrays. Efficient genome-wide distribution of nine TR families implies the role of transposition only in expansion of the most dispersed family, and involvement of other mechanisms is anticipated. Despite similarities in sequence features, FISH experiments indicate high-level compartmentalization of centromeric and euchromatic tandem repeats.     

Development and genetic mapping of SSR markers in foxtail millet [Setaria italica (L.) P. Beauv.]

SSR markers are desirable markers in analysis of genetic diversity, quantitative trait loci mapping and gene locating. In this study, SSR markers were developed from two genomic libraries enriched for (GA)n and (CA)n of foxtail millet [Setaria italica (L.) P. Beauv.], a crop of historical importance in China. A total of 100 SSR markers among the 193 primer pairs detected polymorphism between two mapping parents of an F₂ population, i.e. “B100” of cultivated S. italica and “A10” of wild S. viridis. Excluding 14 markers with unclear amplifications, and five markers unlinked with any linkage group, a foxtail millet SSR linkage map was constructed by integrating 81 new developed SSR markers with 20 RFLP anchored markers. The 81 SSRs covered nine chromosomes of foxtail millet. The length of the map was 1,654 cM, with an average interval distance between markers of 16.4 cM. The 81 SSR markers were not evenly distributed throughout the nine chromosomes, with Ch.8 harbouring the least (3 markers) and Ch.9 harbouring the most (18 markers). To verify the usefulness of the SSR markers developed, 37 SSR markers were randomly chosen to analyze genetic diversity of 40 foxtail millet accessions. Totally 228 alleles were detected, with an average 6.16 alleles per locus. Polymorphism information content (PIC) value for each locus ranged from 0.413 to 0.847, with an average of 0.697. A positive correlation between PIC and number of alleles and between PIC and number of repeat unit were found [0.802 and 0.429, respectively (P < 0.01)]. UPGMA analysis revealed that the 40 foxtail millet cultivars could be grouped into five clusters in which the landraces’ grouping was largely consistent with ecotypes while the breeding varieties from different provinces in China tended to be grouped together. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genome-Wide Investigation and Expression Analyses of WD40 Protein Family in the Model Plant Foxtail Millet (Setaria italica L.)

WD40 proteins play a crucial role in diverse protein-protein interactions by acting as scaffolding molecules and thus assisting in the proper activity of proteins. Hence, systematic characterization and expression profiling of these WD40 genes in foxtail millet would enable us to understand the networks of WD40 proteins and their biological processes and gene functions. In the present study, a genome-wide survey was conducted and 225 potential WD40 genes were identified. Phylogenetic analysis categorized the WD40 proteins into 5 distinct sub-families (I–V). Gene Ontology annotation revealed the biological roles of the WD40 proteins along with its cellular components and molecular functions. In silico comparative mapping with sorghum, maize and rice demonstrated the orthologous relationships and chromosomal rearrangements including duplication, inversion and deletion of WD40 genes. Estimation of synonymous and non-synonymous substitution rates revealed its evolutionary significance in terms of gene-duplication and divergence. Expression profiling against abiotic stresses provided novel insights into specific and/or overlapping expression patterns of SiWD40 genes. Homology modeling enabled three-dimensional structure prediction was performed to understand the molecular functions of WD40 proteins. Although, recent findings had shown the importance of WD40 domains in acting as hubs for cellular networks during many biological processes, it has invited a lesser research attention unlike other common domains. Being a most promiscuous interactors, WD40 domains are versatile in mediating critical cellular functions and hence this genome-wide study especially in the model crop foxtail millet would serve as a blue-print for functional characterization of WD40s in millets and bioenergy grass species. In addition, the present analyses would also assist the research community in choosing the candidate WD40s for comprehensive studies towards crop improvement of millets and biofuel grasses.     

High-Throughput Discovery of Mutations in Tef Semi-Dwarfing Genes by Next-Generation Sequencing Analysis

Tef (Eragrostis tef) is a major cereal crop in Ethiopia. Lodging is the primary constraint to increasing productivity in this allotetraploid species, accounting for losses of ∼15–45% in yield each year. As a first step toward identifying semi-dwarf varieties that might have improved lodging resistance, an ∼6× fosmid library was constructed and used to identify both homeologues of the dw3 semi-dwarfing gene of Sorghum bicolor. An EMS mutagenized population, consisting of ∼21,210 tef plants, was planted and leaf materials were collected into 23 superpools. Two dwarfing candidate genes, homeologues of dw3 of sorghum and rht1 of wheat, were sequenced directly from each superpool with 454 technology, and 120 candidate mutations were identified. Out of 10 candidates tested, six independent mutations were validated by Sanger sequencing, including two predicted detrimental mutations in both dw3 homeologues with a potential to improve lodging resistance in tef through further breeding. This study demonstrates that high-throughput sequencing can identify potentially valuable mutations in under-studied plant species like tef and has provided mutant lines that can now be combined and tested in breeding programs for improved lodging resistance.     

Population structure and association mapping of yield contributing agronomic traits in foxtail millet

Key message

Association analyses accounting for population structure and relative kinship identified eight SSR markers ( p < 0.01) showing significant association ( R ²  = 18 %) with nine agronomic traits in foxtail millet.

Abstract

Association mapping is an efficient tool for identifying genes regulating complex traits. Although association mapping using genomic simple sequence repeat (SSR) markers has been successfully demonstrated in many agronomically important crops, very few reports are available on marker-trait association analysis in foxtail millet. In the present study, 184 foxtail millet accessions from diverse geographical locations were genotyped using 50 SSR markers representing the nine chromosomes of foxtail millet. The genetic diversity within these accessions was examined using a genetic distance-based and a general model-based clustering method. The model-based analysis using 50 SSR markers identified an underlying population structure comprising five sub-populations which corresponded well with distance-based groupings. The phenotyping of plants was carried out in the field for three consecutive years for 20 yield contributing agronomic traits. The linkage disequilibrium analysis considering population structure and relative kinship identified eight SSR markers (p < 0.01) on different chromosomes showing significant association (R ² = 18 %) with nine agronomic traits. Four of these markers were associated with multiple traits. The integration of genetic and physical map information of eight SSR markers with their functional annotation revealed strong association of two markers encoding for phospholipid acyltransferase and ubiquitin carboxyl-terminal hydrolase located on the same chromosome (5) with flag leaf width and grain yield, respectively. Our findings on association mapping is the first report on Indian foxtail millet germplasm and this could be effectively applied in foxtail millet breeding to further uncover marker-trait associations with a large number of markers.     

Feeding behavior and development response of the perilla seed bugs (Nysius sp.) (Heteroptera: Lygaeidae) on multiple crop seeds combinations

The perilla seed bugs, Nysius plebeius Distant and Nysius hidakai Nakatani, are emerging pests of several field crop species including perilla crop in Korea. Here, we investigated the oviposition, nymphal mortality, development, adult emergence, adult longevity, adult weight and female ratio of N. plebeius and N. hidakai on multiple (MC) and double (DC) combinations of five crop seed species: perilla [Perilla frutescens (L.) Britton.], sesame (Sesamum indicum L.), sorghum [Sorghum bicolor (L.)], proso millet (Panicum miliaceum L.), and foxtail millet (Setaria italic L.) in no-choice test. The life-history variables of N. plebeius and N. hidakai were significantly influenced by the crop seed combinations. The most preferred crop seed combinations for oviposition were MC (perilla + sorghum + sesame + proso millet + foxtail millet), followed by DC (perilla + foxtail millet and perilla + sorghum), whereas the least preferred seed combination was DC (proso millet + foxtail millet). Both Nysius sp. were developed faster in MC (perilla + sorghum + sesame + proso millet + foxtail millet), with minimal nymphal mortality, higher adult emergence, longer adult longevity and heavier adult weight. The positive adult bug’s performance on the MC (perilla + sorghum + sesame + proso millet + foxtail millet) may be due to the nutrient contents and physicochemical components associated with the crop seeds. Further, the findings of this study are discussed in terms of bug's behavioral response to crop seed species combinations and management strategies of Nysius sp. as an alternative to chemical management.