Novel Expressed Sequence Tag-Derived and Other Genomic Simple Sequence Repeat Markers Revealed Genetic Diversity in Ethiopian Finger Millet Landrace Populations and Cultivars

Finger millet RNA-seq reveals differential gene expression associated with tolerance to aluminum toxicity and provides novel genomic resources

Eleusine coracana, finger millet, is a multipurpose crop cultivated in arid and semi-arid regions of Africa and Asia. RNA sequencing (RNA-seq) was used in this study to obtain valuable genomic resources and identify genes differentially expressed between Al-tolerant and Al-susceptible genotypes. Two groups of finger millet genotypes were used: Al-tolerant (215836, 215845, and 229722) and Al-susceptible (212462, 215804 and 238323). The analysis of the RNA-seq data resulted in 198,546 unigenes, 56.5% of which were annotated with significant hits in one or more of the following six databases: NR (48.8%), GO (29.7%), KEGG (45%), PlantTFDB (19.0%), Uniprot (49.2%), and NT (46.2%). It is noteworthy that only 220 unigenes in the NR database had significant hits against finger millet sequences suggesting that finger millet's genomic resources are scarce. The gene expression analysis revealed that 322 genes were significantly differentially expressed between the Al-tolerant and Al-susceptible genotypes, of which 40.7% were upregulated while 59.3% were downregulated in Al-tolerant genotypes. Among the significant DEGs, 54.7% were annotated in the GO database with the top hits being ATP binding (GO:0005524) and DNA binding (GO:0003677) in the molecular function, DNA integration (GO:0015074) and cell redox homeostasis in the biological process, as well as cellular anatomical entity and intracellular component in the cellular component GO classes. Several of the annotated DEGs were significantly enriched for their corresponding GO terms. The KEGG pathway analysis resulted in 60 DEGs that were annotated with different pathway classes, of which carbohydrate metabolism and signal transduction were the most prominent. The homologs of a number of significant DEGs have been previously reported as being associated with Al or other abiotic stress responses in various crops, including carboxypeptidase SOL1, HMA3, AP2, bZIP, C3H, and WRKY TF genes. A more detailed investigation of these and other DEGs will enable genomic-led breeding for Al tolerance in finger millet. RNA-seq data analysis also yielded 119,073 SNP markers, the majority of which had PIC values above 0.3, indicating that they are highly informative. Additionally, 3,553 single-copy SSR markers were identified, of which trinucleotide SSRs were the most prevalent. These genomic resources contribute substantially to the enrichment of genomic databases for finger millet, and facilitate future research on this crop.

Novel GBS-Based SNP Markers for Finger Millet and Their Use in Genetic Diversity Analyses

Eleusine coracana (L.) Gaertn., commonly known as finger millet, is a multipurpose crop used for food and feed. Genomic tools are required for the characterization of crop gene pools and their genomics-led breeding. High-throughput sequencing-based characterization of finger millet germplasm representing diverse agro-ecologies was considered an effective method for determining its genetic diversity, thereby suggesting potential candidates for breeding. In this study, the genotyping-by-sequencing (GBS) method was used to simultaneously identify novel single nucleotide polymorphism (SNP) markers and genotype 288 finger millet accessions collected from Ethiopia and Zimbabwe. The accessions were characterized at individual and group levels using 5,226 bi-allelic SNPs, with a minimum allele frequency (MAF) of above 0.05, distributed across 2,500 scaffolds of the finger millet reference genome. The polymorphism information content (PIC) of the SNPs was 0.23 on average, and a quarter of them have PIC values over 0.32, making them highly informative. The grouping of the 288 accessions into seven populations based on geographic proximity and the potential for germplasm exchange revealed a narrow range of observed heterozygosity (Ho; 0.09–0.11) and expected heterozygosity (He) that ranged over twofold, from 0.11 to 0.26. Alleles unique to the different groups were also identified, which merit further investigation for their potential association with desirable traits. The analysis of molecular variance (AMOVA) revealed a highly significant genetic differentiation among groups of accessions classified based on the geographic region, country of origin, days to flowering, panicle type, and Al tolerance (p < 0.01). The high genetic differentiation between Ethiopian and Zimbabwean accessions was evident in the AMOVA, cluster, principal coordinate, and population structure analyses. The level of genetic diversity of finger millet accessions varies moderately among locations within Ethiopia, with accessions from the northern region having the lowest level. In the neighbor-joining cluster analysis, most of the improved cultivars included in this study were closely clustered, probably because they were developed using genetically less diverse germplasm and/or selected for similar traits, such as grain yield. The recombination of alleles via crossbreeding genetically distinct accessions from different regions of the two countries can potentially lead to the development of superior cultivars.