Exploring the genetic basis of anthracnose resistance in Ethiopian sorghum through a genome-wide association study.

BACKGROUND: Sorghum anthracnose is a major disease that hampers the productivity of the crop globally. The disease is caused by the hemibiotrophic fungal pathogen Colletotrichum sublineola. The identification of anthracnose-resistant sorghum genotypes, defining resistance loci and the underlying genes, and their introgression into adapted cultivars are crucial for enhancing productivity. In this study, we conducted field experiments on 358 diverse accessions of Ethiopian sorghum. Quantitative resistance to anthracnose was evaluated at locations characterized by a heavy natural infestation that is suitable for disease resistance screening. RESULTS: The field-based screening identified 53 accessions that were resistant across locations, while 213 accessions exhibited variable resistance against local pathotypes. Genome-wide association analysis (GWAS) was performed using disease response scores on 329 accessions and 83,861 single nucleotide polymorphisms (SNPs) generated through genotyping-by-sequencing (GBS). We identified 38 loci significantly associated with anthracnose resistance. Interestingly, a subset of these loci harbor genes encoding receptor-like kinases (RLK), nucleotide-binding leucine-rich repeats (NLRs), stress-induced antifungal tyrosine kinase that have been previously implicated in disease resistance. A SNP on chromosome 4 (S04_66140995) and two SNPs on chromosome 2 (S02_75784037, S02_2031925), localized with-in the coding region of genes that encode a putative stress-induced antifungal kinase, an F-Box protein, and Xa21-binding RLK that were strongly associated with anthracnose resistance. We also identified highly significant associations between anthracnose resistance and three SNPs linked to genes (Sobic.002G058400, Sobic.008G156600, Sobic.005G033400) encoding an orthologue of the widely known NLR protein (RPM1), Leucine Rich Repeat family protein, and Heavy Metal Associated domain-containing protein, respectively. Other SNPs linked to predicted immune response genes were also significantly associated with anthracnose resistance. CONCLUSIONS: The sorghum germplasm collections used in the present study are genetically diverse. They harbor potentially useful, yet undiscovered, alleles for anthracnose resistance. This is supported by the identification of novel loci that are enriched for disease resistance regulators such as NLRs, LRKs, Xa21-binding LRK, and antifungal proteins. The genotypic data available for these accessions offer a valuable resource for sorghum breeders to effectively improve the crop. The genomic regions and candidate genes identified can be used to design markers for molecular breeding of sorghum diseases resistance.

Assessing genetic, racial, and geographic diversity among Ethiopian sorghum landraces and implications for heterotic potential for hybrid sorghum breeding.

The wealth of sorghum genetic resources in Africa has not been fully exploited for cultivar development in the continent. Hybrid cultivars developed from locally evolved germplasm are more likely to possess a well-integrated assembly of genes for local adaptation, productivity, quality, as well as for defensive traits and broader stability. A subset of 560 sorghum accessions of known fertility reaction representing the major botanical races and agro-ecologies of Ethiopia were characterized for genetic, agronomic and utilization parameters to lay a foundation for cultivar improvement and parental selection for hybrid breeding. Accessions were genotyped using a genotyping by sequencing (GBS) generating 73,643 SNPs for genetic analysis. Significant genetic variability was observed among accessions with Admixture and Discriminant Analysis of Principal Components where 67% of the accessions fell into K=10 clusters with membership coefficient set to > 0.6. The pattern of aggregation of the accessions partially overlapped with racial category and agro-ecological adaptation. Majority of the non-restorer (B-line) accessions primarily of the bicolor race from the wet highland ecology clustered together away from two clusters of fertility restorer (R-line) accessions. Small members of the B accessions were grouped with the R clusters and in vice-versa while significant numbers of both B and R accessions were spread between the major clusters. Such pattern of diversity along with the complementary agronomic data based information indicate the potential for heterosis providing the foundation for initiating hybrid breeding program based on locally adapted germplasm. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-024-01483-8.

Haplotypes at the sorghum ARG4 and ARG5 NLR loci confer resistance to anthracnose.

Sorghum anthracnose caused by the fungus Colletotrichum sublineola (Cs) is a damaging disease of the crop. Here, we describe the identification of ANTHRACNOSE RESISTANCE GENES (ARG4 and ARG5) encoding canonical nucleotide-binding leucine-rich repeat (NLR) receptors. ARG4 and ARG5 are dominant resistance genes identified in the sorghum lines SAP135 and P9830, respectively, that show broad-spectrum resistance to Cs. Independent genetic studies using populations generated by crossing SAP135 and P9830 with TAM428, fine mapping using molecular markers, comparative genomics and gene expression studies determined that ARG4 and ARG5 are resistance genes against Cs strains. Interestingly, ARG4 and ARG5 are both located within clusters of duplicate NLR genes at linked loci separated by ~1 Mb genomic region. SAP135 and P9830 each carry only one of the ARG genes while having the recessive allele at the second locus. Only two copies of the ARG5 candidate genes were present in the resistant P9830 line while five non-functional copies were identified in the susceptible line. The resistant parents and their recombinant inbred lines carrying either ARG4 or ARG5 are resistant to strains Csgl1 and Csgrg suggesting that these genes have overlapping specificities. The role of ARG4 and ARG5 in resistance was validated through sorghum lines carrying independent recessive alleles that show increased susceptibility. ARG4 and ARG5 are located within complex loci displaying interesting haplotype structures and copy number variation that may have resulted from duplication. Overall, the identification of anthracnose resistance genes with unique haplotype stucture provides a foundation for genetic studies and resistance breeding.

Genome-wide association analysis reveals seed protein loci as determinants of variations in grain mold resistance in sorghum.

KEY MESSAGE: GWAS analysis revealed variations at loci harboring seed storage, late embryogenesis abundant protein, and a tannin biosynthesis gene associated with sorghum grain mold resistance. Grain mold is the most important disease of sorghum [Sorghum bicolor (L.) Moench]. It starts at the early stages of grain development due to concurrent infection by multiple fungal species. The genetic architecture of resistance to grain mold is poorly understood. Using a diverse set of 635 Ethiopian sorghum accessions, we conducted a multi-stage disease rating for resistance to grain mold under natural infestation in the field. Through genome-wide association analyses with 173,666 SNPs and multiple models, two novel loci were identified that were consistently associated with grain mold resistance across environments. Sequence variation at new loci containing sorghum KAFIRIN gene encoding a seed storage protein affecting seed texture and LATE EMBRYOGENESIS ABUNDANT 3 (LEA3) gene encoding a protein that accumulates in seeds, previously implicated in stress tolerance, were significantly associated with grain mold resistance. The KAFIRIN and LEA3 loci were also significant factors in grain mold resistance in accessions with non-pigmented grains. Moreover, we consistently detected the known SNP (S4_62316425) in TAN1 gene, a regulator of tannin accumulation in sorghum grain to be significantly associated with grain mold resistance. Identification of loci associated with new mechanisms of resistance provides fresh insight into genetic control of the trait, while the highly resistant accessions can serve as sources of resistance genes for breeding. Overall, our association data suggest the critical role of loci harboring seed protein genes and implicate grain chemical and physical properties in sorghum grain mold resistance.

A comprehensive phenotypic and genomic characterization of Ethiopian sorghum germplasm defines core collection and reveals rich genetic potential in adaptive traits.

Understanding population genetic structure and diversity of a crop is essential in designing selection strategies in plant breeding. About 2010 Ethiopian sorghum accessions were phenotyped for different traits at multiple locations. A subset of the collection, 1628 accessions, predominantly landraces, some improved varieties, and inbred lines were genotyped by sequencing. Phenotypic data revealed association of important traits with different sorghum growing agro-climatic regions, high genetic diversity and the presence of rare natural variation in the Ethiopian sorghum germplasm. Subsequent genotypic analysis determined optimum number of sub-populations, distinct cluster groups and ancestries of each sorghum accessions. To improve utilization of germplasm, a core subset of 387 lines were selected following posteriori grouping of genotypes based on cluster groups obtained through GBS analysis followed by stratified random sampling using quantitative traits. In order to evaluate how well this new sorghum and millet innovation lab (SMIL) collection from Ethiopia is represented within the largest world sorghum collection at United States Department of Agriculture - National Plant Germplasm System (USDA-NPGS) and the sorghum association panel (SAP), comparisons were conducted based on SNP data. The SMIL collection displayed high genetic diversity with some redundancy with the USDA-NPGS germplasm but SAP showed clear distinction. Furthermore, genome-environment association analysis identified candidate genes associated with adaptation to abiotic factors, that will be important for exploitation of adaptive potential to different environments. In summary, our results described the diversity and relationship of sorghum collections, representativeness of developed core and provide novel insights into candidate genes associated to abiotic stress tolerance.

Global mRNA and microRNA expression dynamics in response to anthracnose infection in sorghum.

BACKGROUND: Anthracnose is a damaging disease of sorghum caused by the fungal pathogen Colletotrichum sublineolum. Genome-wide mRNA and microRNA (miRNA) profiles of resistant and susceptible sorghum genotypes were studied to understand components of immune responses, and fungal induced miRNA and target gene networks. RESULTS: A total of 18 mRNA and 12 miRNA libraries from resistant and susceptible sorghum lines were sequenced prior to and after inoculation with C. sublineolum. Significant differences in transcriptomes of the susceptible and resistant genotypes were observed with dispersion distance and hierarchical cluster tree analyses. Of the total 33,032 genes predicted in the sorghum genome, 19,593 were induced by C. sublineolum, and 15,512 were differentially expressed (DEGs) between the two genotypes. The resistant line was marked by significant reprogramming of the transcriptome at 24 h post inoculation (hpi), and a decrease at 48 hpi, whereas the susceptible line displayed continued changes in gene expression concordant with elevated fungal growth in the susceptible genotype. DEGs encode proteins implicated in diverse functions including photosynthesis, synthesis of tetrapyrrole, carbohydrate and secondary metabolism, immune signaling, and chitin binding. Genes encoding immune receptors, MAPKs, pentatricopeptide repeat proteins, and WRKY transcription factors were induced in the resistant genotype. In a parallel miRNA profiling, the susceptible line displayed greater number of differentially expressed miRNAs than the resistant line indicative of a widespread suppression of gene expression. Interestingly, we found 75 miRNAs, including 36 novel miRNAs, which were differentially expressed in response to fungal inoculation. The expression of 50 miRNAs was significantly different between resistant and susceptible lines. Subsequently, for 35 differentially expressed miRNAs, the corresponding 149 target genes were identified. Expression of 56 target genes were significantly altered after inoculation, showing inverse expression with the corresponding miRNAs. CONCLUSIONS: We provide insights into genome wide dynamics of mRNA and miRNA profiles, biological and cellular processes underlying host responses to fungal infection in sorghum. Resistance is correlated with early transcriptional reprogramming of genes in various pathways. Fungal induced genes, miRNAs and their targets with a potential function in host responses to anthracnose were identified, opening avenues for genetic dissection of resistance mechanisms.

Identification of candidate flowering and sex genes in white Guinea yam (D. rotundata Poir.) by SuperSAGE transcriptome profiling.

Dioecy (distinct male and female individuals) and scarce to non-flowering are common features of cultivated yam (Dioscorea spp.). However, the molecular mechanisms underlying flowering and sex determination in Dioscorea are largely unknown. We conducted SuperSAGE transcriptome profiling of male, female and monoecious individuals to identify flowering and sex-related genes in white Guinea yam (D. rotundata), generating 20,236 unique tags. Of these, 13,901 were represented by a minimum of 10 tags. A total 88 tags were significantly differentially expressed in male, female and monoecious plants, of which 18 corresponded to genes previously implicated in flower development and sex determination in multiple plant species. We validated the SuperSAGE data with quantitative real-time PCR (qRT-PCR)-based analysis of the expression of three candidate genes. We further investigated the flowering patterns of 1938 D. rotundata accessions representing diverse geographical origins over two consecutive years. Over 85% of accessions were either male or non-flowering, less than 15% were female, while monoecious plants were rare. Intensity of flowering varied between male and female plants, with the former flowering more abundantly than the latter. Candidate genes identified in this study can be targeted for further validation and to induce regular flowering in poor to non-flowering cultivars. Findings of the study provide important inputs for further studies aiming to overcome the challenge of flowering in yams and to improve efficiency of yam breeding.

A Large-Scale Genome-Wide Association Analyses of Ethiopian Sorghum Landrace Collection Reveal Loci Associated With Important Traits.

The eastern Africa region, Ethiopia and its surroundings, is considered as the center of origin and diversity for sorghum, and has contributed to global sorghum genetic improvement. The germplasm from this region harbors enormous genetic variation for various traits but little is known regarding the genetic architecture of most traits. Here, 1425 Ethiopian landrace accessions were phenotyped under field conditions for presence or absence of awns, panicle compactness and shape, panicle exsertion, pericarp color, glume cover, plant height and smut resistance under diverse environmental conditions in Ethiopia. In addition, F1 hybrids obtained from a subset of 1341 accessions crossed to an A1 cytoplasmic male sterile line, ATx623, were scored for fertility/sterility reactions. Subsequently, genotyping-by-sequencing generated a total of 879,407 SNPs from which 72,190 robust SNP markers were selected after stringent quality control (QC). Pairwise distance-based hierarchical clustering identified 11 distinct groups. Of the genotypes assigned to either one of the 11 sub-populations, 65% had high ancestry membership coefficient with the likelihood of more than 0.60 and the remaining 35% represented highly admixed accessions. A genome-wide association study (GWAS) identified loci and SNPs associated with aforementioned traits. GWAS based on compressed mixed linear model (CMLM) identified SNPs with significant association (FDR ≤ 0.05) to the different traits studied. The percentage of total phenotypic variation explained with significant SNPs across traits ranged from 2 to 43%. Candidate genes showing significant association with different traits were identified. The sorghum bHLH transcription factor, ABORTED MICROSPORES was identified as a strong candidate gene conditioning male fertility. Notably, sorghum CLAVATA1 receptor like kinase, known for regulation of plant growth, and the ETHYLENE RESPONSIVE TRANSCRIPTION FACTOR gene RAP2-7, known to suppress transition to flowering, were significantly associated with plant height. In addition, the YELLOW SEED1 like MYB transcription factor and TANNIN1 showed strong association with pericarp color validating previous observations. Overall, the genetic architecture of natural variation representing the complex Ethiopian sorghum germplasm was established. The study contributes to the characterization of genes and alleles controlling agronomic traits, and will serve as a source of markers for molecular breeding.

Genetic diversity and population structure of a mini-core subset from the world cowpea (Vigna unguiculata (L.) Walp.) germplasm collection.

The International Institute of Tropical Agriculture maintains the world's largest collection of cowpea germplasm of over 15,000 accessions. A sub-set of 298 lines from the loosely composed mini core collection of 370 landraces were genotyped based on genotyping by sequencing (GBS). Ward's minimum variance hierarchical cluster analysis, model-based ancestry analysis and discriminant analysis of principal component (DAPC) were carried out on this sub-set. Three clusters were identified by the different clustering methods. Principal component analysis further supported the three clusters especially when accessions are scattered along the axes of the first two principal components. The first two principal components explained a total of 22.30% of the variation. Cluster one comprises 115 accessions from the largest number of countries and has the highest gene diversity, heterozygosity and polymorphic information content (PIC) values. Cluster two is made up of 102 accessions, 90 percent of which are from West and Central Africa. Analysis of molecular variance shows that the most variation is among accessions and lowest among clusters. No cluster is made exclusively of accessions from a single country. Based on SNP markers, the sub set of cowpea mini core germplasm collection used in this study encompasses the diversity in the crop.

Genome sequencing of the staple food crop white Guinea yam enables the development of a molecular marker for sex determination.

BACKGROUND: Root and tuber crops are a major food source in tropical Africa. Among these crops are several species in the monocotyledonous genus Dioscorea collectively known as yam, a staple tuber crop that contributes enormously to the subsistence and socio-cultural lives of millions of people, principally in West and Central Africa. Yam cultivation is constrained by several factors, and yam can be considered a neglected "orphan" crop that would benefit from crop improvement efforts. However, the lack of genetic and genomic tools has impeded the improvement of this staple crop. RESULTS: To accelerate marker-assisted breeding of yam, we performed genome analysis of white Guinea yam (Dioscorea rotundata) and assembled a 594-Mb genome, 76.4% of which was distributed among 21 linkage groups. In total, we predicted 26,198 genes. Phylogenetic analyses with 2381 conserved genes revealed that Dioscorea is a unique lineage of monocotyledons distinct from the Poales (rice), Arecales (palm), and Zingiberales (banana). The entire Dioscorea genus is characterized by the occurrence of separate male and female plants (dioecy), a feature that has limited efficient yam breeding. To infer the genetics of sex determination, we performed whole-genome resequencing of bulked segregants (quantitative trait locus sequencing [QTL-seq]) in F1 progeny segregating for male and female plants and identified a genomic region associated with female heterogametic (male = ZZ, female = ZW) sex determination. We further delineated the W locus and used it to develop a molecular marker for sex identification of Guinea yam plants at the seedling stage. CONCLUSIONS: Guinea yam belongs to a unique and highly differentiated clade of monocotyledons. The genome analyses and sex-linked marker development performed in this study should greatly accelerate marker-assisted breeding of Guinea yam. In addition, our QTL-seq approach can be utilized in genetic studies of other outcrossing crops and organisms with highly heterozygous genomes. Genomic analysis of orphan crops such as yam promotes efforts to improve food security and the sustainability of tropical agriculture.

Tracking crop varieties using genotyping-by-sequencing markers: a case study using cassava (Manihot esculenta Crantz).

BACKGROUND: Accurate identification of crop cultivars is crucial in assessing the impact of crop improvement research outputs. Two commonly used identification approaches, elicitation of variety names from farmer interviews and morphological plant descriptors, have inherent uncertainty levels. Genotyping-by-sequencing (GBS) was used in a case study as an alternative method to track released varieties in farmers' fields, using cassava, a clonally propagated root crop widely grown in the tropics, and often disseminated through extension services and informal seed systems. A total of 917 accessions collected from 495 farming households across Ghana were genotyped at 56,489 SNP loci along with a "reference library" of 64 accessions of released varieties and popular landraces. RESULTS: Accurate cultivar identification and ancestry estimation was accomplished through two complementary clustering methods: (i) distance-based hierarchical clustering; and (ii) model-based maximum likelihood admixture analysis. Subsequently, 30% of the identified accessions from farmers' fields were matched to specific released varieties represented in the reference library. ADMIXTURE analysis revealed that the optimum number of major varieties was 11 and matched the hierarchical clustering results. The majority of the accessions (69%) belonged purely to one of the 11 groups, while the remaining accessions showed two or more ancestries. Further analysis using subsets of SNP markers reproduced results obtained from the full-set of markers, suggesting that GBS can be done at higher DNA multiplexing, thereby reducing the costs of variety fingerprinting. A large proportion of discrepancy between genetically unique cultivars as identified by markers and variety names as elicited from farmers were observed. Clustering results from ADMIXTURE analysis was validated using the assumption-free Discriminant Analysis of Principal Components (DAPC) method. CONCLUSION: We show that genome-wide SNP markers from increasingly affordable GBS methods coupled with complementary cluster analysis is a powerful tool for fine-scale population structure analysis and variety identification. Moreover, the ancestry estimation provides a framework for quantifying the contribution of exotic germplasm or older improved varieties to the genetic background of contemporary improved cultivars.

Next-generation sequencing based genotyping, cytometry and phenotyping for understanding diversity and evolution of Guinea yams.

KEY MESSAGE: Genotyping by sequencing (GBS) is used to understand the origin and domestication of guinea yams, including the contribution of wild relatives and polyploidy events to the cultivated guinea yams. Patterns of genetic diversity within and between two cultivated guinea yams (Dioscorea rotundata and D. cayenensis) and five wild relatives (D. praehensilis, D. mangenotiana, D. abyssinica, D. togoensis and D. burkilliana) were investigated using next-generation sequencing (genotyping by sequencing, GBS). Additionally, the two cultivated species were assessed for intra-specific morphological and ploidy variation. In guinea yams, ploidy level is correlated with species identity. Using flow cytometry a single ploidy level was inferred across D. cayenensis (3x, N = 21), D. praehensilis (2x, N = 7), and D. mangenotiana (3x, N = 5) accessions, whereas both diploid and triploid (or aneuploid) accessions were present in D. rotundata (N = 11 and N = 32, respectively). Multi-dimensional scaling and maximum parsimony analyses of 2,215 SNPs revealed that wild guinea yam populations form discrete genetic groupings according to species. D. togoensis and D. burkilliana were most distant from the two cultivated yam species, whereas D. abyssinica, D. mangenotiana, and D. praehensilis were closest to cultivated yams. In contrast, cultivated species were genetically less clearly defined at the intra-specific level. While D. cayenensis formed a single genetic group, D. rotundata comprised three separate groups consisting of; (1) a set of diploid individuals genetically similar to D. praehensilis, (2) a set of diploid individuals genetically similar to D. cayenensis, and (3) a set of triploid individuals. The current study demonstrates the utility of GBS for assessing yam genomic diversity. Combined with morphological and biological data, GBS provides a powerful tool for testing hypotheses regarding the evolution, domestication and breeding of guinea yams.