An Integrated Genotyping-by-Sequencing Polymorphism Map for Over 10,000 Sorghum Genotypes

Milletomics: a metabolomics centered integrated omics approach toward genetic progression

Producing alternative staple foods like millet will be essential to feeding ten billion people by 2050. The increased demand for millet is driving researchers to improve its genetic variation. Millets include protein, dietary fiber, phenolic substances, and flavonoid components. Its climate resilience makes millet an appealing crop for agronomic sustainability. Integrative omics technologies could potentially identify and develop millets with desirable phenotypes that may have high agronomic value. Millets' salinity and drought tolerance have been enhanced using transcriptomics. In foxtail, finger, and pearl millet, proteomics has discovered salt-tolerant protein, phytohormone-focused protein, and drought tolerance. Metabolomics studies have revealed that certain metabolic pathways including those involving lignin, flavonoids, phenylpropanoid, and lysophospholipids are critical for many processes, including seed germination, photosynthesis, energy metabolism, and the synthesis of bioactive chemicals necessary for drought tolerance. Metabolomics integration with other omics revealed metabolome engineering and trait-specific metabolite creation. Integrated metabolomics and ionomics are still in the development stage, but they could potentially assist in comprehending the pathway of ionomers to control nutrient levels and biofortify millet. Epigenomic analysis has shown alterations in DNA methylation patterns and chromatin structure in foxtail and pearl millets in response to abiotic stress. Whole-genome sequencing utilizing next-generation sequencing is the most proficient method for finding stress-induced phytoconstituent genes. New genome sequencing enables novel biotechnological interventions including genome-wide association, mutation-based research, and other omics approaches. Millets can breed more effectively by employing next-generation sequencing and genotyping by sequencing, which may mitigate climate change. Millet marker-assisted breeding has advanced with high-throughput markers and combined genotyping technologies.

Sorghum SbGhd7 is a major regulator of floral transition and directly represses genes crucial for flowering activation

Molecular genetic understanding of flowering time regulation is crucial for sorghum development. GRAIN NUMBER, PLANT HEIGHT AND HEADING DATE 7 (SbGhd7) is one of the six classical loci conferring photoperiod sensitivity of sorghum flowering. However, its functions remain poorly studied. The molecular functions of SbGhd7 were characterized. The gene regulatory network controlled by SbGhd7 was constructed and validated. The biological roles of SbGhd7 and its major targets were studied. SbGhd7 overexpression (OE) completely prevented sorghum flowering. Additionally, we show that SbGhd7 is a major negative regulator of flowering, binding to the promoter motif TGAATG(A/T)(A/T/C) and repressing transcription of the major florigen FLOWERING LOCUS T 10 (SbFT10) and floral activators EARLY HEADING DATE (SbEhd1), FLAVIN-BINDING, KELCH REPEAT, F-BOX1 (SbFKF1) and EARLY FLOWERING 3 (SbELF3). Reinforcing the direct effect of SbGhd7, SbEhd1 OE activated the promoters of three functional florigens (SbFT1, SbFT8 and SbFT10), dramatically accelerating flowering. Our studies demonstrate that SbGhd7 is a major repressor of sorghum flowering by directly and indirectly targeting genes for flowering activation. The mechanism appears ancient. Our study extends the current model of floral transition regulation in sorghum and provides a framework for a comprehensive understanding of sorghum photoperiod response.

Genome-Wide Association Analysis Uncovers Genes Associated with Resistance to Head Smut Pathotype 5 in Senegalese Sorghum Accessions

A newly documented pathotype 5 of the soil-borne fungus Sporisorium reilianum, causing head smut in sorghum, was tested against 153 unexplored Senegalese sorghum accessions. Among the 153 sorghum accessions tested, 63 (41%) exhibited complete resistance, showing no signs of infection by the fungus. The remaining 90 accessions (59%) displayed varying degrees of susceptibility. Sorghum responses against S. reilianum were explored to analyze the potential link with previously known seed morphology-related traits and new phenotype data from 59 lines for seed weight. A genome-wide association study (GWAS) screened 297,876 SNPs and identified highly significant associations (p < 1 × 10-5) with head smut resistance in sorghum. By mapping these significant SNPs to the reference genome, this study revealed 35 novel candidate defense genes potentially involved in disease resistance.

Identification of genetic and environmental factors influencing aerial root traits that support biological nitrogen fixation in sorghum

Plant breeding and genetics play a major role in the adaptation of plants to meet human needs. The current requirement to make agriculture more sustainable can be partly met by a greater reliance on biological nitrogen fixation by symbiotic diazotrophic microorganisms that provide crop plants with ammonium. Select accessions of the cereal crop sorghum (Sorghum bicolor (L.) Moench) form mucilage-producing aerial roots that harbor nitrogen-fixing bacteria. Breeding programs aimed at developing sorghum varieties that support diazotrophs will benefit from a detailed understanding of the genetic and environmental factors contributing to aerial root formation. A genome-wide association study of the sorghum minicore, a collection of 242 landraces, and 30 accessions from the sorghum association panel was conducted in Florida and Wisconsin and under 2 fertilizer treatments to identify loci associated with the number of nodes with aerial roots and aerial root diameter. Sequence variation in genes encoding transcription factors that control phytohormone signaling and root system architecture showed significant associations with these traits. In addition, the location had a significant effect on the phenotypes. Concurrently, we developed F2 populations from crosses between bioenergy sorghums and a landrace that produced extensive aerial roots to evaluate the mode of inheritance of the loci identified by the genome-wide association study. Furthermore, the mucilage collected from aerial roots contained polysaccharides rich in galactose, arabinose, and fucose, whose composition displayed minimal variation among 10 genotypes and 2 fertilizer treatments. These combined results support the development of sorghums with the ability to acquire nitrogen via biological nitrogen fixation.

Parallel tuning of semi-dwarfism via differential splicing of Brachytic1 in commercial maize and smallholder sorghum

In the current genomic era, the search and deployment of new semi-dwarf alleles have continued to develop better plant types in all cereals. We characterized an agronomically optimal semi-dwarf mutation in Zea mays L. and a parallel polymorphism in Sorghum bicolor L. We cloned the maize brachytic1 (br1-Mu) allele by a modified PCR-based Sequence Amplified Insertion Flanking Fragment (SAIFF) approach. Histology and RNA-Seq elucidated the mechanism of semi-dwarfism. GWAS linked a sorghum plant height QTL with the Br1 homolog by resequencing a West African sorghum landraces panel. The semi-dwarf br1-Mu allele encodes an MYB transcription factor78 that positively regulates stalk cell elongation by interacting with the polar auxin pathway. Semi-dwarfism is due to differential splicing and low functional Br1 wild-type transcript expression. The sorghum ortholog, SbBr1, co-segregates with the major plant height QTL qHT7.1 and is alternatively spliced. The high frequency of the Sbbr1 allele in African landraces suggests that African smallholder farmers used the semi-dwarf allele to improve plant height in sorghum long before efforts to introduce Green Revolution-style varieties in the 1960s. Surprisingly, variants for differential splicing of Brachytic1 were found in both commercial maize and smallholder sorghum, suggesting parallel tuning of plant architecture across these systems.

Genome-wide association analysis identified molecular markers and candidate genes for flower traits in Chinese orchid (Cymbidium sinense)

The orchid, the champagne of flowers, brings luxury, elegance, and novelty to nature. Cymbidium sinense is a symbol of gigantic floral variability on account of wavering shapes and sizes of floral organs, although marker-trait association (MTA) has not been studied for its floral traits. We evaluated markers associated with 14 floral traits of C. sinense through a genome-wide association study (GWAS) of 195 accessions. A total of 65 318 522 single-nucleotide polymorphisms (SNPs) and 3 906 176 insertion/deletion (InDel) events were identified through genotyping-by-sequencing. Among these, 4694 potential SNPs and 477 InDels were identified as MTAs at -log10 P > 5. The genes related to these SNPs and InDels were largely associated with floral regulators, hormonal pathways, cell division, and metabolism, playing essential roles in tailoring floral morphology. Moreover, 20 candidate SNPs/InDels linked to 11 genes were verified, 8 of which were situated on exons, one was located in the 5'-UTR and two were positioned in introns. Here, the multitepal trait-related gene RABBIT EARS (RBE) was found to be the most crucial gene. We analyzed the role of CsRBE in the regulation of flower-related genes via efficient transient overexpression in C. sinense protoplasts, and found that the floral homeotic genes CsAP3 and CsPI, as well as organ boundary regulators, including CsCUC and CsTCP genes, were regulated by CsRBE. Thus, we obtained key gene loci for important ornamental traits of orchids using genome-wide association analysis of populations with natural variation. The findings of this study can do a great deal to expedite orchid breeding programs for shape variability.

Genetic diversity, population structure and anthracnose resistance response in a novel sweet sorghum diversity panel

Sweet sorghum is an attractive feedstock for the production of renewable chemicals and fuels due to the readily available fermentable sugars that can be extracted from the juice, and the additional stream of fermentable sugars that can be obtained from the cell wall polysaccharides in the bagasse. An important selection criterion for new sweet sorghum germplasm is resistance to anthracnose, a disease caused by the fungal pathogen Colletotrichum sublineolum. The identification of novel anthracnose-resistance sources present in sweet sorghum germplasm offers a fast track towards the development of new resistant sweet sorghum germplasm. We established a sweet sorghum diversity panel (SWDP) of 272 accessions from the USDA-ARS National Plant Germplasm (NPGS) collection that includes landraces from 22 countries and advanced breeding material, and that represents ~15% of the NPGS sweet sorghum collection. Genomic characterization of the SWDP identified 171,954 single nucleotide polymorphisms (SNPs) with an average of one SNP per 4,071 kb. Population structure analysis revealed that the SWDP could be stratified into four populations and one admixed group, and that this population structure could be aligned to sorghum's racial classification. Results from a two-year replicated trial of the SWDP for anthracnose resistance response in Texas, Georgia, Florida, and Puerto Rico showed 27 accessions to be resistant across locations, while 145 accessions showed variable resistance response against local pathotypes. A genome-wide association study identified 16 novel genomic regions associated with anthracnose resistance. Four resistance loci on chromosomes 3, 6, 8 and 9 were identified against pathotypes from Puerto Rico, and two resistance loci on chromosomes 3 and 8 against pathotypes from Texas. In Georgia and Florida, three resistance loci were detected on chromosomes 4, 5, 6 and four on chromosomes 4, 5 (two loci) and 7, respectively. One resistance locus on chromosome 2 was effective against pathotypes from Texas and Puerto Rico and a genomic region of 41.6 kb at the tip of chromosome 8 was associated with resistance response observed in Georgia, Texas, and Puerto Rico. This publicly available SWDP and the extensive evaluation of anthracnose resistance represent a valuable genomic resource for the improvement of sorghum.

Harnessing the Genetic Basis of Sorghum Biomass-Related Traits to Facilitate Bioenergy Applications

The extensive use of fossil fuels and global climate change have raised ever-increasing attention to sustainable development, global food security and the replacement of fossil fuels by renewable energy. Several C4 monocot grasses have excellent photosynthetic ability, stress tolerance and may rapidly produce biomass in marginal lands with low agronomic inputs, thus representing an important source of bioenergy. Among these grasses, Sorghum bicolor has been recognized as not only a promising bioenergy crop but also a research model due to its diploidy, simple genome, genetic diversity and clear orthologous relationship with other grass genomes, allowing sorghum research to be easily translated to other grasses. Although sorghum molecular genetic studies have lagged far behind those of major crops (e.g., rice and maize), recent advances have been made in a number of biomass-related traits to dissect the genetic loci and candidate genes, and to discover the functions of key genes. However, molecular and/or targeted breeding toward biomass-related traits in sorghum have not fully benefited from these pieces of genetic knowledge. Thus, to facilitate the breeding and bioenergy applications of sorghum, this perspective summarizes the bioenergy applications of different types of sorghum and outlines the genetic control of the biomass-related traits, ranging from flowering/maturity, plant height, internode morphological traits and metabolic compositions. In particular, we describe the dynamic changes of carbohydrate metabolism in sorghum internodes and highlight the molecular regulators involved in the different stages of internode carbohydrate metabolism, which affects the bioenergy utilization of sorghum biomass. We argue the way forward is to further enhance our understanding of the genetic mechanisms of these biomass-related traits with new technologies, which will lead to future directions toward tailored designing sorghum biomass traits suitable for different bioenergy applications.

Genomics and phenomics enabled prebreeding improved early-season chilling tolerance in Sorghum

In temperate climates, earlier planting of tropical-origin crops can provide longer growing seasons, reduce water loss, suppress weeds, and escape post-flowering drought stress. However, chilling sensitivity of sorghum, a tropical-origin cereal crop, limits early planting, and over 50 years of conventional breeding has been stymied by coinheritance of chilling tolerance (CT) loci with undesirable tannin and dwarfing alleles. In this study, phenomics and genomics-enabled approaches were used for prebreeding of sorghum early-season CT. Uncrewed aircraft systems (UAS) high-throughput phenotyping platform tested for improving scalability showed moderate correlation between manual and UAS phenotyping. UAS normalized difference vegetation index values from the chilling nested association mapping population detected CT quantitative trait locus (QTL) that colocalized with manual phenotyping CT QTL. Two of the 4 first-generation Kompetitive Allele Specific PCR (KASP) molecular markers, generated using the peak QTL single nucleotide polymorphisms (SNPs), failed to function in an independent breeding program as the CT allele was common in diverse breeding lines. Population genomic fixation index analysis identified SNP CT alleles that were globally rare but common to the CT donors. Second-generation markers, generated using population genomics, were successful in tracking the donor CT allele in diverse breeding lines from 2 independent sorghum breeding programs. Marker-assisted breeding, effective in introgressing CT allele from Chinese sorghums into chilling-sensitive US elite sorghums, improved early-planted seedling performance ratings in lines with CT alleles by up to 13-24% compared to the negative control under natural chilling stress. These findings directly demonstrate the effectiveness of high-throughput phenotyping and population genomics in molecular breeding of complex adaptive traits.

The potentialities of omics resources for millet improvement

Millets are nutrient-rich (nutri-rich) cereals with climate resilience attributes. However, its full productive potential is not realized due to the lack of a focused yield improvement approach, as evidenced by the available literature. Also, the lack of well-characterized genomic resources significantly limits millet improvement. But the recent availability of genomic data and advancement in omics tools has shown its enormous potential to enhance the efficiency and precision faced by conventional breeding in millet improvement. The development of high throughput genotyping platforms based on next-generation sequencing (NGS) has provided a low-cost method for genomic information, specifically for neglected nutri-rich cereals with the availability of a limited number of reference genome sequences. NGS has created new avenues for millet biotechnological interventions such as mutation-based study, GWAS, GS, and other omics technologies. The simultaneous discovery of high-throughput markers and multiplexed genotyping platform has aggressively aided marker-assisted breeding for millet improvement. Therefore, omics technology offers excellent opportunities to explore and combine useful variations for targeted traits that could impart high nutritional value to high-yielding cultivars under changing climatic conditions. In millet improvement, an in-depth account of NGS, integrating genomics data with different biotechnology tools, is reviewed in this context.

Genome-Wide Association Study of Seed Morphology Traits in Senegalese Sorghum Cultivars

Sorghum is considered the fifth most important crop in the world. Despite the potential value of Senegalese germplasm for various traits, such as resistance to fungal diseases, there is limited information on the study of sorghum seed morphology. In this study, 162 Senegalese germplasms were evaluated for seed area size, length, width, length-to-width ratio, perimeter, circularity, the distance between the intersection of length & width (IS) and center of gravity (CG), and seed darkness and brightness by scanning and analyzing morphology-related traits with SmartGrain software at the USDA-ARS Plant Science Research Unit. Correlations between seed morphology-related traits and traits associated with anthracnose and head smut resistance were analyzed. Lastly, genome-wide association studies were performed on phenotypic data collected from over 16,000 seeds and 193,727 publicly available single nucleotide polymorphisms (SNPs). Several significant SNPs were found and mapped to the reference sorghum genome to uncover multiple candidate genes potentially associated with seed morphology. The results indicate clear correlations among seed morphology-related traits and potential associations between seed morphology and the defense response of sorghum. GWAS analysis listed candidate genes associated with seed morphologies that can be used for sorghum breeding in the future.

Multi-Trait Genome-Wide Association Studies of Sorghum bicolor Regarding Resistance to Anthracnose, Downy Mildew, Grain Mold and Head Smut

Multivariate linear mixed models (mvLMMs) are widely applied for genome-wide association studies (GWAS) to detect genetic variants affecting multiple traits with correlations and/or different plant growth stages. Subsets of multiple sorghum populations, including the Sorghum Association Panel (SAP), the Sorghum Mini Core Collection and the Senegalese sorghum population, have been screened against various sorghum diseases such as anthracnose, downy mildew, grain mold and head smut. Still, these studies were generally performed in a univariate framework. In this study, we performed GWAS based on the principal components of defense-related multi-traits against the fungal diseases, identifying new potential SNPs (S04_51771351, S02_66200847, S09_47938177, S08_7370058, S03_72625166, S07_17951013, S04_66666642 and S08_51886715) associated with sorghum's defense against these diseases.

Characterization of grain carotenoids in global sorghum germplasm to guide genomics-assisted breeding strategies

BACKGROUND

Crop biofortification is a successful strategy to ameliorate Vitamin A deficiency. Sorghum is a good candidate for vitamin A biofortification, as it is a staple food in regions with high prevalence of vitamin A deficiency. β-carotene-the main provitamin A carotenoid-is below the target concentration in sorghum grain, therefore biofortification breeding is required. Previous studies found evidence that sorghum carotenoid variation is oligogenic, suggesting that marker-assisted selection can be an appropriate biofortification method. However, we hypothesize that sorghum carotenoids have both oligogenic and polygenic components of variation. Genomics-assisted breeding could accelerate breeding efforts, but there exists knowledge gaps in the genetics underlying carotenoid variation, as well as appropriate germplasm to serve as donors.

RESULTS

In this study, we characterized carotenoids in 446 accessions from the sorghum association panel and carotenoid panel using high-performance liquid chromatography, finding high carotenoid accessions not previously identified. Genome-wide association studies conducted with 345 accessions, confirmed that zeaxanthin epoxidase is a major gene underlying variation for not only zeaxanthin, but also lutein and β-carotene. High carotenoid lines were found to have limited genetic diversity, and originated predominantly from only one country. Potential novel genetic diversity for carotenoid content was identified through genomic predictions in 2,495 accessions of unexplored germplasm. Oligogenic variation of carotenoids was confirmed, as well as evidence for polygenic variation, suggesting both marker-assisted selection and genomic selection can facilitate breeding efforts.

CONCLUSIONS

Sorghum vitamin A biofortification could be beneficial for millions of people who rely on it as a dietary staple. Carotenoid content in sorghum is low, but high heritability suggests that increasing concentrations through breeding is possible. Low genetic diversity among high carotenoid lines might be the main limitation for breeding efforts, therefore further germplasm characterization is needed to assess the feasibility of biofortification breeding. Based on germplasm here evaluated, most countries' germplasm lacks high carotenoid alleles, thus pre-breeding will be needed. A SNP marker within the zeaxanthin epoxidase gene was identified as a good candidate for use in marker-assisted selection. Due to the oligogenic and polygenic variation of sorghum grain carotenoids, both marker-assisted selection and genomic selection can be employed to accelerate breeding efforts.

GWAS of grain color and tannin content in Chinese sorghum based on whole-genome sequencing

Seventy-three QTL related to grain color and tannin content were identified in Chinese sorghum accessions, and a new recessive allelic variant of TAN2 gene was discovered. Sorghum is mainly used for brewing distilled liquors in China. Since grain tannins play an important role in liquor brewing, accurately understanding the relationship between grain color and tannin content can provide basis for selection standards of tannin sorghum. We resequenced a panel of 242 Chinese sorghum accessions and performed population structure and genome-wide association study (GWAS) to identify quantitative trait locus (QTL) affecting pericarp color, testa pigment, and tannin content. Phylogenetic analysis, principal component analysis (PCA), and admixture model were used to infer population structure. Two distinct genetic sub-populations were identified according to their corresponding northern and southern geographic origin. To investigate the genetic basis of natural variation in sorghum grain color, GWAS with 2,760,264 SNPs was conducted in four environments using multiple models (Blink, FarmCPU, GLM, and MLM). Seventy-three QTL were identified to be associated for the color of exocarp, mesocarp, testa, and tannin content on all chromosomes except chromosome 5, of which 47 might be novel QTL. Some important QTL were found to colocalize with orthologous genes in the flavonoid biosynthetic pathway from other plants, including orthologous of Arabidopsis (Arabidopsis thaliana) TT2, TT7, TT12, TT16 and AT5G41220 (GST), as well as orthologous of rice (Oryza sativa) MYB61 and OsbHLH025. Our investigation of the variation in grain color and tannin content in Chinese sorghum germplasm may help guide future sorghum breeding for liquor brewing.

Novel QTL for chilling tolerance at germination and early seedling stages in sorghum

Sorghum (Sorghum bicolor L.) a drought tolerant staple crop for half a billion people in Africa and Asia, an important source of animal feed throughout the world and a biofuel feedstock of growing importanceorghum’s originated from tropical regions rendering the crop to be cold sensitive. Low temperature stresses such as chilling and frost greatly affect the agronomic performance of sorghum and limit its geographical distribution, posing a major problem in temperate environments when sorghum is planted early. Understanding the genetic basis of wide adaptability and of sorghum would facilitate molecular breeding programs and studies of other C4 crops. The objective of this study is to conduct quantitative trait loci analysis using genotying by sequencing for early seed germination and seedling cold tolerance in two sorghum recombinant inbred lines populations. To accomplish that, we used two populations of recombinant inbred lines (RIL) developed from crosses between cold-tolerant (CT19, ICSV700) and cold-sensitive (TX430, M81E) parents. The derived RIL populations were evaluated for single nucleotide polymorphism (SNP) using genotype-by-sequencing (GBS) in the field and under controlled environments for their response to chilling stress. Linkage maps were constructed with 464 and 875 SNPs for the CT19 X TX430 (C1) and ICSV700 X M81 E (C2) populations respectively. Using quantitative trait loci (QTL) mapping, we identified QTL conferring tolerance to chilling tolerance at the seedling stage. A total of 16 and 39 total QTL were identified in the C1 and C2 populations, respectively. Two major QTL were identified in the C1 population, and three major QTL were mapped in the C2 population. Comparisons between the two populations and with previously identified QTL show a high degree of similarity in QTL locations. Given the amount of co-localization of QTL across traits and the direction of allelic effect supports that these regions have a pleiotropic effect. These QTL regions were also identified to be highly enriched for genes encoding chilling stress and hormonal response genes. This identified QTL can be useful in developing tools for molecular breeding of sorghums with improved low-temperature germinability.

Host Determinants of Fungal Species Composition and Symptom Manifestation in the Sorghum Grain Mold Disease Complex

Sorghum grain mold (SGM) is an important multifungal disease complex affecting sorghum (Sorghum bicolor) production systems worldwide. SGM-affected sorghum grain can be contaminated with potent fumonisin mycotoxins produced by Fusarium verticillioides, a prevalent SGM-associated taxon. Historically, efforts to improve resistance to SGM have achieved only limited success. Classical approaches to evaluating SGM resistance are based solely on disease severity, which offers little insight regarding the distinct symptom manifestations within the disease complex. In this study, three novel phenotypes were developed to facilitate assessment of SGM symptom manifestation. A sorghum diversity panel composed of 390 accessions was inoculated with endogenous strains of F. verticillioides and evaluated for these phenotypes, as well as for the conventional panicle grain mold severity rating phenotype, in South Carolina, U.S.A., in 2017 and 2019. Distributions of phenotype values were examined, broad-sense heritability was estimated, and relationships to botanical race were explored. A typology of SGM symptom manifestations was developed to classify accessions using principal component analysis and k-means clustering, constituting a novel option for basing breeding decisions on SGM outcomes more nuanced than disease severity. Genome-wide association studies were performed using SGM trait data, resulting in the identification of 19 significant single nucleotide polymorphisms in linkage disequilibrium with a total of 86 gene models. Our findings provide a basis of exploratory evidence regarding the genetic architecture of SGM symptom manifestation and indicate that traits other than disease severity could be tractable targets for SGM resistance breeding.

A Genome-Wide Association Study of Senegalese Sorghum Seedlings Responding to Pathotype 5 of Sporisorium reilianum

Sporisorium reilianum is a fungal pathogen that causes head smut in sorghum. In addition to pathotypes (P) 1-4, P5 and P6 were identified recently. In this study, seedlings of Senegalese sorghum, comprising 163 accessions, were evaluated for response to Sporisorium reilianum. Teliospores of pathotype P5 of the pathogen in dilute agar were pipetted onto seedling shoots while still in soil, and inoculated seedlings were submerged under water at 4 days post-inoculation. Signs of infection (noticeable spots) on the first leaf were checked daily up to 6 days post submergence. A genome-wide association study (GWAS) was conducted using 193,727 single-nucleotide polymorphisms (SNPs) throughout the genome based on two types of phenotypic data: whether noticeable spots were shown or not and the average time for an observation of the spots across 163 accessions. When mapped back to the reference sorghum genome, most of the top candidate SNP loci were associated with plant defense or plant stress response-related genes. The identified SNP loci were associated with spot appearance in sorghum seedlings under flooding following inoculation with P5 of Sporisorium reilianum.

Genome-wide analyses using multi-locus models revealed marker-trait associations for major agronomic traits in Sorghum bicolor

Globally, sorghum is the fifth most important cereal crop, and it is a major crop in Ethiopia, where it has a high genetic diversity. The country's sorghum gene pool contributes significantly to sorghum improvement worldwide. This study aimed to identify genomic regions and candidate genes associated with major agronomic traits in sorghum by using its genetic resources in Ethiopia for a genome-wide association study (GWAS). Phenotypic data of days to flowering (DTF), plant height (PH), panicle length (PALH), panicle width (PAWD), panicle weight (PAWT), and grain yield (GY) were collected from a GWAS panel comprising 324 sorghum accessions grown in three environments. SeqSNP, a targeted genotyping method, was used to genotype the panel using 5,000 gene-based single nucleotide polymorphism (SNP) markers. For marker-trait association (MTA) analyses, fixed and random model circulating probability unification (FarmCPU), and Bayesian-information and linkage-disequilibrium iteratively nested keyway (BLINK) models were used. In all traits, high phenotypic variation was observed, with broad-sense heritability ranging from 0.32 (for GY) to 0.90 (for PALH). A population structure, principal component analysis, and kinship analysis revealed that the accessions could be divided into two groups. In total, 54 MTAs were identified, 11 of which were detected by both BLINK and farmCPU. MTAs identified for each trait ranged from five (PAWT and GY) to fourteen (PH) representing both novel and previously identified quantitative trait loci (QTLs). Three SNPs were associated with more than one trait, including a SNP within the Sobic.004G189200 gene that was associated with PH and PAWT. Major effect SNP loci, Sbi2393610 (PVE = 23.3%), Sbi10438246 (PVE = 35.2%), Sbi17789352 (PVE = 11.9%) and Sbi30169733 (PVE = 18.9%) on chromosomes 1, 3, 5 and 9 that showed strong association signals for PAWD, DTF, GY and PALH, respectively, were major findings of this study. The SNP markers and candidate genes identified in this study provide insights into the genetic control of grain yield and related agronomic traits, and once validated, the markers could be used in genomics-led breeding.

G-DIRT: a web server for identification and removal of duplicate germplasms based on identity-by-state analysis using single nucleotide polymorphism genotyping data

Maintaining duplicate germplasms in genebanks hampers effective conservation and utilization of genebank resources. The redundant germplasm adds to the cost of germplasm conservation by requiring a large proportion of the genebank financial resources towards conservation rather than enriching the diversity. Besides, genome-wide-association analysis using an association panel with over-represented germplasms can be biased resulting in spurious marker-trait associations. The conventional methods of germplasm duplicate removal using passport information suffer from incomplete or missing passport information and data handling errors at various stages of germplasm enrichment. This limitation is less likely in the case of genotypic data. Therefore, we developed a web-based tool, Germplasm Duplicate Identification and Removal Tool (G-DIRT), which allows germplasm duplicate identification based on identity-by-state analysis using single-nucleotide polymorphism genotyping information along with pre-processing of genotypic data. A homozygous genotypic difference threshold of 0.1% for germplasm duplicates has been determined using tetraploid wheat genotypic data with 94.97% of accuracy. Based on the genotypic difference, the tool also builds a dendrogram that can visually depict the relationship between genotypes. To overcome the constraint of high-dimensional genotypic data, an offline version of G-DIRT in the interface of R has also been developed. The G-DIRT is expected to help genebank curators, breeders and other researchers across the world in identifying germplasm duplicates from the global genebank collections by only using the easily sharable genotypic data instead of physically exchanging the seeds or propagating materials. The web server will complement the existing methods of germplasm duplicate identification based on passport or phenotypic information being freely accessible at http://webtools.nbpgr.ernet.in/gdirt/.

Sorghum Association Panel whole-genome sequencing establishes cornerstone resource for dissecting genomic diversity

Association mapping panels represent foundational resources for understanding the genetic basis of phenotypic diversity and serve to advance plant breeding by exploring genetic variation across diverse accessions. We report the whole-genome sequencing (WGS) of 400 sorghum (Sorghum bicolor (L.) Moench) accessions from the Sorghum Association Panel (SAP) at an average coverage of 38× (25-72×), enabling the development of a high-density genomic marker set of 43 983 694 variants including single-nucleotide polymorphisms (approximately 38 million), insertions/deletions (indels) (approximately 5 million), and copy number variants (CNVs) (approximately 170 000). We observe slightly more deletions among indels and a much higher prevalence of deletions among CNVs compared to insertions. This new marker set enabled the identification of several novel putative genomic associations for plant height and tannin content, which were not identified when using previous lower-density marker sets. WGS identified and scored variants in 5-kb bins where available genotyping-by-sequencing (GBS) data captured no variants, with half of all bins in the genome falling into this category. The predictive ability of genomic best unbiased linear predictor (GBLUP) models was increased by an average of 30% by using WGS markers rather than GBS markers. We identified 18 selection peaks across subpopulations that formed due to evolutionary divergence during domestication, and we found six F_st peaks resulting from comparisons between converted lines and breeding lines within the SAP that were distinct from the peaks associated with historic selection. This population has served and continues to serve as a significant public resource for sorghum research and demonstrates the value of improving upon existing genomic resources.

Genome-wide association study of Senegalese sorghum seedlings responding to a Texas isolate of Colletotrichum sublineola

Colletotrichum sublineola is a destructive fungal pathogen that causes anthracnose in sorghum. Senegalese sorghum germplasm is currently being considered as an option of sources for genetic resistance. In a recent study, Senegalese sorghum accessions were evaluated for response to a mixture of Texas isolates of C. sublineola at the 8-leaf stage in the greenhouse. As a comparison, 159 Senegalese sorghum accessions at the 1-leaf developmental stage were evaluated against a single Texas isolate of C. sublineola (FSP53) using an excised-leaf assay. A genome-wide association study (GWAS) was conducted based on the phenotypic data acquired to discover genetic variation associated with response to C. sublineola using 193,727 single nucleotide polymorphisms (SNPs) throughout the genome. Sorghum seedlings tended to be more resistant when compared with sorghum plants inoculated at the 8-leaf stage in the greenhouse in previous experiments. Based on the highest score evaluated in the 1-leaf developmental stage excised leaf assay for each accession, 16 accessions were labeled as susceptible. GWAS identified the SNP locus S01_72868925 that is associated with protein kinase domain // Leucine rich repeat N-terminal domain at a level of confidence that surpassed Bonferroni correction. Along with the SNP locus S01_72868925, other top SNP loci were also associated with genes that are known to play critical roles in plant defense or plant stress responses.

A pan-genome data structure induced by pooled sequencing facilitates variant mining in heterogeneous germplasm

Valuable genetic variation lies unused in gene banks due to the difficulty of exploiting heterogeneous germplasm accessions. Advances in molecular breeding, including transgenics and genome editing, present the opportunity to exploit hidden sequence variation directly. Here we describe the pan-genome data structure induced by whole-genome sequencing of pooled individuals from wild populations of Patellifolia spp., a source of disease resistance genes for the related crop species sugar beet (Beta vulgaris). We represent the pan-genome as a map of reads from pooled sequencing of a heterogeneous population sample to a reference genome, plus a BLAST data base of the mapped reads. We show that this basic data structure can be queried by reference genome position or homology to identify sequence variants present in the wild relative, at genes of agronomic interest in the crop, a process known as allele or variant mining. Further we demonstrate the possibility of cataloging variants in all Patellifolia genomic regions that have corresponding single copy orthologous regions in sugar beet. The data structure, termed a "pooled read archive," can be produced, altered, and queried using standard tools to facilitate discovery of agronomically-important sequence variation.

Supplementary information

The online version contains supplementary material available at 10.1007/s11032-022-01308-6.

Genome-wide association mapping of resistance to the sorghum aphid in Sorghum bicolor

Since 2013, the sorghum aphid (SA), Melanaphis sorghi (Theobald), has been a serious pest that hampers all types of sorghum production in the U.S. Known sorghum aphid resistance in sorghum is limited to a few genetic regions on SBI-06. In this study, a subset of the Sorghum Association Panel (SAP) was used along with some additional lines to identify genomic regions that confer sorghum aphid resistance. SAP lines were grown in the field and visually evaluated for SA resistance during the growing seasons of 2019 and 2020 in Tifton, GA. In 2020, the SAP accessions were also evaluated for SA resistance in the field using drone-based high throughput phenotyping (HTP). Flowering time was recorded in the field to confirm that our methods were sufficient for identifying known quantitative trait loci (QTL). This study combined phenotypic data from field-based visual ratings and reflectance data to identify genome-wide associated (GWAS) marker-trait associations (MTA) using genotyping-by-sequencing (GBS) data. Several MTAs were identified for SA-related traits across the genome, with a few common markers that were consistently identified on SBI-08 and SBI-10 for aphid count and plant damage, as well as loci for reflectance-based traits on SBI-02, SBI-03, and SBI-05. Candidate genes encoding leucine-rich repeats (LRR), Avr proteins, lipoxygenases (LOXs), calmodulins (CAM) dependent protein kinase, WRKY transcription factors, flavonoid biosynthesis genes, and 12-oxo-phytodienoic acid reductase were identified near SNPs that had significant associations with different SA traits. In this study, flowering time-related genes were also identified as a positive control for the methods. The total phenotypic variation explained by significant SNPs across SA-scored traits, reflectance data, and flowering time ranged from 6 to 61%, while the heritability value ranged from 4 to 69%. This study identified three new sources of resistant lines to sorghum aphid. These results supported the existing literature, and also revealed several new loci. Markers identified in this study will support marker-assisted breeding for sorghum aphid resistance.

Genome-wide association studies identify putative pleiotropic locus mediating drought tolerance in sorghum

Drought is a key constraint on plant productivity and threat to food security. Sorghum (Sorghum bicolor L. Moench), a global staple food and forage crop, is among the most drought-adapted cereal crops, but its adaptation is not yet well understood. This study aims to better understand the genetic basis of preflowering drought in sorghum and identify loci underlying variation in water use and yield components under drought. A panel of 219 diverse sorghum from West Africa was phenotyped for yield components and water use in an outdoor large-tube lysimeter system under well-watered (WW) versus a preflowering drought water-stressed (WS) treatment. The experimental system was validated based on characteristic drought response in international drought tolerant check genotypes and genome-wide association studies (GWAS) that mapped the major height locus at QHT7.1 and Dw3. GWAS further identified marker trait associations (MTAs) for drought-related traits (plant height, flowering time, forage biomass, grain weight, water use) that each explained 7-70% of phenotypic variance. Most MTAs for drought-related traits correspond to loci not previously reported, but some MTA for forage biomass and grain weight under WS co-localized with staygreen post-flowering drought tolerance loci (Stg3a and Stg4). A globally common allele at S7_50055849 is associated with several yield components under drought, suggesting that it tags a major pleiotropic variant controlling assimilate partitioning to grain versus vegetative biomass. The GWAS revealed oligogenic variants for drought tolerance in sorghum landraces, which could be used as trait predictive markers for improved drought adaptation.

Quantitative and population genomics suggest a broad role of stay-green loci in the drought adaptation of sorghum

Drought is a major constraint on plant productivity globally. Sorghum [Sorghum bicolor (L.) Moench] landraces have evolved in drought-prone regions, but the genetics of their adaptation is poorly understood. Here we sought to identify novel drought-tolerance loci and test hypotheses on the role of known loci including those underlying stay-green (Stg) postflowering drought tolerance. We phenotyped 590 diverse sorghum accessions from West Africa in 10 environments, under field-based managed drought stress [preflowering water stress (WS1), postflowering water stress (WS2), and well-watered (WW)] and rainfed (RF) conditions over 4 yr. Days to 50% flowering (DFLo), aboveground dry biomass (DBM), plant height (PH), and plant grain yield components (including grain weight [GrW], panicle weight [PW] and grain number [GrN] per plant, and 1000-grain weight [TGrW]) were measured, and genome-wide association studies (GWAS) was conducted. Broad-sense heritability for biomass and plant grain yield was high (33-92%) across environments. There was a significant correlation between stress tolerance index (STI) for GrW per plant across WS1 and WS2. Genome-wide association studies revealed that SbZfl1 and SbCN12, orthologs of maize (Zea mays L.) flowering genes, likely underlie flowering time variation under these conditions. Genome-wide association studies further identified associations (n = 134; common between two GWAS models) for STI and drought effects on plant yield components including 16 putative pleiotropic associations. Thirty of the associations colocalized with Stg1, Stg2, Stg3, and Stg4 loci and had large effects. Seven lead associations, including some within Stg1, overlapped with positive selection outliers. Our findings reveal previously undescribed natural genetic variation for drought-tolerance-related traits and suggest a broad role of Stg loci in drought adaptation of sorghum.

Decoding the sorghum methylome: understanding epigenetic contributions to agronomic traits

DNA methylation is a chromatin modification that plays an essential role in regulating gene expression and genome stability and it is typically associated with gene silencing and heterochromatin. Owing to its heritability, alterations in the patterns of DNA methylation have the potential to provide for epigenetic inheritance of traits. Contemporary epigenomic technologies provide information beyond sequence variation and could supply alternative sources of trait variation for improvement in crops such as sorghum. Yet, compared with other species such as maize and rice, the sorghum DNA methylome is far less well understood. The distribution of CG, CHG, and CHH methylation in the genome is different compared with other species. CG and CHG methylation levels peak around centromeric segments in the sorghum genome and are far more depleted in the gene dense chromosome arms. The genes regulating DNA methylation in sorghum are also yet to be functionally characterised; better understanding of their identity and functional analysis of DNA methylation machinery mutants in diverse genotypes will be important to better characterise the sorghum methylome. Here, we catalogue homologous genes encoding methylation regulatory enzymes in sorghum based on genes in Arabidopsis, maize, and rice. Discovering variation in the methylome may uncover epialleles that provide extra information to explain trait variation and has the potential to be applied in epigenome-wide association studies or genomic prediction. DNA methylation can also improve genome annotations and discover regulatory elements underlying traits. Thus, improving our knowledge of the sorghum methylome can enhance our understanding of the molecular basis of traits and may be useful to improve sorghum performance.

Current status and prospects of herbicide-resistant grain sorghum (Sorghum bicolor)

Grain sorghum is a versatile crop, which can thrive under limited water and other inputs. However, crop loss from weed infestation continues to be a major constraint in grain sorghum production. Particularly, post-emergence grass weed control is a great challenge in grain sorghum due to the lack of herbicide options. Unlike in other major crops, such as maize or soybean, herbicide-resistant sorghum technology that can facilitate weed control throughout crop growing season is not available to growers yet. The development of herbicide-resistant sorghum can have potential to improve weed management, including post-emergence grass weed control. One of the major concerns in the development of such technology in sorghum is escape of resistance traits into weedy relatives of sorghum (e.g. shattercane and johnsongrass). This review focuses on sources of herbicide resistance in sorghum, the status of the development of herbicide-resistant sorghum technologies, overview of breeding methods, and limitations in the development of such sorghum technology as well as economic benefits for sorghum growers. © 2021 Society of Chemical Industry.

A sorghum genome-wide association study (GWAS) identifies a WRKY transcription factor as a candidate gene underlying sugarcane aphid (Melanaphis sacchari) resistance

A WRKY transcription factor identified through forward genetics is associated with sorghum resistance to the sugarcane aphid and through heterologous expression reduces aphid populations in multiple plant species. Crop plant resistance to insect pests is based on genetically encoded traits which often display variability across diverse germplasm. In a comparatively recent event, a predominant sugarcane aphid (SCA: Melanaphis sacchari) biotype has become a significant agronomic pest of grain sorghum (Sorghum bicolor). To uncover candidate genes underlying SCA resistance, we used a forward genetics approach combining the genetic diversity present in the Sorghum Association Panel (SAP) and the Bioenergy Association Panel (BAP) for a genome-wide association study, employing an established SCA damage rating. One major association was found on Chromosome 9 within the WRKY transcription factor 86 (SbWRKY86). Transcripts encoding SbWRKY86 were previously identified as upregulated in SCA-resistant germplasm and the syntenic ortholog in maize accumulates following Rhopalosiphum maidis infestation. Analyses of SbWRKY86 transcripts displayed patterns of increased SCA-elicited accumulation in additional SCA-resistant sorghum lines. Heterologous expression of SbWRKY86 in both tobacco (Nicotiana benthamiana) and Arabidopsis resulted in reduced population growth of green peach aphid (Myzus persicae). Comparative RNA-Seq analyses of Arabidopsis lines expressing 35S:SbWRKY86-YFP identified changes in expression for a small network of genes associated with carbon-nitrogen metabolism and callose deposition, both contributing factors to defense against aphids. As a test of altered plant responses, 35S:SbWRKY86-YFP Arabidopsis lines were activated using the flagellin epitope elicitor, flg22, and displayed significant increases in callose deposition. Our findings indicate that both heterologous and increased native expression of the transcription factor SbWRKY86 contributes to reduced aphid levels in diverse plant models.

Sorghum genetic, genomic, and breeding resources

Sorghum research has entered an exciting and fruitful era due to the genetic, genomic, and breeding resources that are now available to researchers and plant breeders. As the world faces the challenges of a rising population and a changing global climate, new agricultural solutions will need to be developed to address the food and fiber needs of the future. To that end, sorghum will be an invaluable crop species as it is a stress-resistant C₄ plant that is well adapted for semi-arid and arid regions. Sorghum has already remained as a staple food crop in many parts of Africa and Asia and is critically important for animal feed and niche culinary applications in other regions, such as the United States. In addition, sorghum has begun to be developed into a promising feedstock for forage and bioenergy production. Due to this increasing demand for sorghum and its potential to address these needs, the continuous development of powerful community resources is required. These resources include vast collections of sorghum germplasm, high-quality reference genome sequences, sorghum association panels for genome-wide association studies of traits involved in food and bioenergy production, mutant populations for rapid discovery of causative genes for phenotypes relevant to sorghum improvement, gene expression atlas, and online databases that integrate all resources and provide the sorghum community with tools that can be used in breeding and genomic studies. Used in tandem, these valuable resources will ensure that the rate, quality, and collaborative potential of ongoing sorghum improvement efforts is able to rival that of other major crops.

Cannabis labelling is associated with genetic variation in terpene synthase genes

Analysis of over 100 Cannabis samples quantified for terpene and cannabinoid content and genotyped for over 100,000 single nucleotide polymorphisms indicated that Sativa- and Indica-labelled samples were genetically indistinct on a genome-wide scale. Instead, we found that Cannabis labelling was associated with variation in a small number of terpenes whose concentrations are controlled by genetic variation at tandem arrays of terpene synthase genes.

Genetic characterization of a Sorghum bicolor multiparent mapping population emphasizing carbon-partitioning dynamics

Sorghum bicolor, a photosynthetically efficient C₄ grass, represents an important source of grain, forage, fermentable sugars, and cellulosic fibers that can be utilized in myriad applications ranging from bioenergy to bioindustrial feedstocks. Sorghum’s efficient fixation of carbon per unit time per unit area per unit input has led to its classification as a preferred biomass crop highlighted by its designation as an advanced biofuel by the U.S. Department of Energy. Due to its extensive genetic diversity and worldwide colonization, sorghum has considerable diversity for a range of phenotypes influencing productivity, composition, and sink/source dynamics. To dissect the genetic basis of these key traits, we present a sorghum carbon-partitioning nested association mapping (NAM) population generated by crossing 11 diverse founder lines with Grassl as the single recurrent female. By exploiting existing variation among cellulosic, forage, sweet, and grain sorghum carbon partitioning regimes, the sorghum carbon-partitioning NAM population will allow the identification of important biomass-associated traits, elucidate the genetic architecture underlying carbon partitioning and improve our understanding of the genetic determinants affecting unique phenotypes within Poaceae. We contrast this NAM population with an existing grain population generated using Tx430 as the recurrent female. Genotypic data are assessed for quality by examining variant density, nucleotide diversity, linkage decay, and are validated using pericarp and testa phenotypes to map known genes affecting these phenotypes. We release the 11-family NAM population along with corresponding genomic data for use in genetic, genomic, and agronomic studies with a focus on carbon-partitioning regimes.

Genome-wide association analysis for response of Senegalese sorghum accessions to Texas isolates of anthracnose

Anthracnose disease of sorghum is caused by Colletotrichum sublineola, a filamentous fungus. The genetic basis of resistance to anthracnose in sorghum is largely unclear, especially in Senegalese sorghum germplasm. In this study, 163 Senegalese sorghum accessions were evaluated for response to C. sublineola, and a genome-wide association study (GWAS) was performed to identify genetic variation associated with response to C. sublineola using 193,727 single nucleotide polymorphisms (SNPs) throughout the genome. Germplasm diversity analysis showed low genetic diversity and slow linkage disequilibrium (LD) decay among the Senegalese accessions. Phenotypic analysis resulted in relatively low differences to C. sublineola among the tested population. Genome-wide association study did not identify any significant association based on a strict threshold for the number of SNPs available. However, individual analysis of the top eight SNPs associated with relative susceptibility and resistance identified candidate genes that have been shown to play important roles in plant stress tolerance in previous studies. This study identifies sorghum genes whose annotated properties have known roles in host defense and thus identify them as candidates for use in breeding for resistance to anthracnose.

Sorghum breeding in the genomic era: opportunities and challenges

The importance and potential of the multi-purpose crop sorghum in global food security have not yet been fully exploited, and the integration of the state-of-art genomics and high-throughput technologies into breeding practice is required. Sorghum, a historically vital staple food source and currently the fifth most important major cereal, is emerging as a crop with diverse end-uses as food, feed, fuel and forage and a model for functional genetics and genomics of tropical grasses. Rapid development in high-throughput experimental and data processing technologies has significantly speeded up sorghum genomic researches in the past few years. The genomes of three sorghum lines are available, thousands of genetic stocks accessible and various genetic populations, including NAM, MAGIC, and mutagenised populations released. Functional and comparative genomics have elucidated key genetic loci and genes controlling agronomical and adaptive traits. However, the knowledge gained has far away from being translated into real breeding practices. We argue that the way forward is to take a genome-based approach for tailored designing of sorghum as a multi-functional crop combining excellent agricultural traits for various end uses. In this review, we update the new concepts and innovation systems in crop breeding and summarise recent advances in sorghum genomic researches, especially the genome-wide dissection of variations in genes and alleles for agronomically important traits. Future directions and opportunities for sorghum breeding are highlighted to stimulate discussion amongst sorghum academic and industrial communities.

Polyphenol oxidase genes as integral part of the evolutionary history of domesticated tetraploid wheat

Studying and understanding the genetic basis of polyphenol oxidases (PPO)-related traits plays a crucial role in genetic improvement of crops. A tetraploid wheat collection (T. turgidum ssp., TWC) was analyzed using the 90K wheat SNP iSelect assay and phenotyped for PPO activity. A total of 21,347 polymorphic SNPs were used to perform genome-wide association analysis (GWA) in TWC and durum wheat sub-groups, detecting 23 and 85 marker-trait associations (MTA). In addition, candidate genes responsible for PPO activity were predicted. Based on the 23 MTAs detected in TWC, two haplotypes associated with low and high PPO activity were identified. Four SNPs were developed and validated providing one reliable marker (IWB75732) for marker assisted selection. The 23 MTAs were used to evaluate the genetic divergence (F_ST > 0.25) between the T. turgidum subspecies, providing new information important for understanding the domestication process of Triticum turgidum ssp. and in particular of ssp. carthlicum.

Population structure, genetic diversity and genomic selection signatures among a Brazilian common bean germplasm

Brazil is the world's largest producer of common bean. Knowledge of the genetic diversity and relatedness of accessions adapted to Brazilian conditions is of great importance for the conservation of germplasm and for directing breeding programs aimed at the development of new cultivars. In this context, the objective of this study was to analyze the genetic diversity, population structure, and linkage disequilibrium (LD) of a diversity panel consisting of 219 common bean accessions, most of which belonging to the Mesoamerican gene pool. Genotyping by sequencing (GBS) of these accessions allowed the identification of 49,817 SNPs with minor allele frequency > 0.05. Of these, 17,149 and 12,876 were exclusive to the Mesoamerican and Andean pools, respectively, and 11,805 SNPs could differentiate the two gene pools. Further the separation according to the gene pool, bayesian analysis of the population structure showed a subdivision of the Mesoamerican accessions based on the origin and color of the seed tegument. LD analysis revealed the occurrence of long linkage blocks and low LD decay with physical distance between SNPs (LD half decay in 249 kb, corrected for population structure and relatedness). The GBS technique could effectively characterize the Brazilian common bean germplasms, and the diversity panel used in this study may be of great use in future genome-wide association studies.

Escalation in the host-pathogen arms race: A host resistance response corresponds to a heightened bacterial virulence response

The zig-zag model of host-pathogen interaction describes the relative strength of defense response across a spectrum of pathogen-induced plant phenotypes. A stronger defense response results in increased resistance. Here, we investigate the strength of pathogen virulence during disease and place these findings in the context of the zig-zag model. Xanthomonas vasicola pv. holcicola (Xvh) causes sorghum bacterial leaf streak. Despite being widespread, this disease has not been described in detail at the molecular level. We divided diverse sorghum genotypes into three groups based on disease symptoms: water-soaked lesions, red lesions, and resistance. Bacterial growth assays confirmed that these three phenotypes represent a range of resistance and susceptibility. To simultaneously reveal defense and virulence responses across the spectrum of disease phenotypes, we performed dual RNA-seq on Xvh-infected sorghum. Consistent with the zig-zag model, the expression of plant defense-related genes was strongest in the resistance interaction. Surprisingly, bacterial virulence genes related to the type III secretion system (T3SS) and type III effectors (T3Es) were also most highly expressed in the resistance interaction. This expression pattern was observed at multiple time points within the sorghum-Xvh pathosystem. Further, a similar expression pattern was observed in Arabidopsis infected with Pseudomonas syringae for effector-triggered immunity via AvrRps4 but not AvrRpt2. Specific metabolites were able to repress the Xvh virulence response in vitro and in planta suggesting a possible signaling mechanism. Taken together, these findings reveal multiple permutations of the continually evolving host-pathogen arms race from the perspective of host defense and pathogen virulence responses.

The arms race between plants and pathogens is a complex process. To dissect the plant defense and pathogen virulence responses simultaneously, we used sorghum and Xanthomonas vasicola pv. holcicola, as a model pathosystem. We performed dual RNA-seq on infected sorghum with a range of disease phenotypes. Our characterization of this pathosystem demonstrates that genes related to the plant defense and pathogen virulence responses are most highly induced during a resistance interaction. We observed a similar pattern of escalation in Arabidopsis infected with Pseudomonas syringae. These observations support a conceptual model of fluidity between the different stages of plant immunity and pathogen virulence.

QTL mapping for soybean (Glycine max L.) leaf chlorophyll-content traits in a genotyped RIL population by using RAD-seq based high-density linkage map

BACKGROUND

Different soybean (Glycine max L.) leaf chlorophyll-content traits are considered to be significantly linked to soybean yield. To map the quantitative trait loci (QTLs) of soybean leaf chlorophyll-content traits, an advanced recombinant inbred line (RIL, ZH, Zhonghuang 24 × Huaxia 3) population was adopted to phenotypic data acquisitions for the target traits across six distinct environments (seasons and soybean growth stages). Moreover, the restriction site-associated DNA sequencing (RAD-seq) based high-density genetic linkage map of the RIL population was utilized for QTL mapping by carrying out the composite interval mapping (CIM) approach.

RESULTS

Correlation analyses showed that most traits were correlated with each other under specific chlorophyll assessing method and were regulated both by hereditary and environmental factors. In this study, 78 QTLs for soybean leaf chlorophyll-content traits were identified. Furthermore, 13 major QTLs and five important QTL hotspots were classified and highlighted from the detected QTLs. Finally, Glyma01g15506, Glyma02g08910, Glyma02g11110, Glyma07g15960, Glyma15g19670 and Glyma15g19810 were predicted from the genetic intervals of the major QTLs and important QTL hotspots.

CONCLUSIONS

The detected QTLs and candidate genes may facilitate to gain a better understanding of the hereditary basis of soybean leaf chlorophyll-content traits and may be valuable to pave the way for the marker-assisted selection (MAS) breeding of the target traits.

Comparisons of sampling methods for assessing intra- and inter-accession genetic diversity in three rice species using genotyping by sequencing

To minimize the cost of sample preparation and genotyping, most genebank genomics studies in self-pollinating species are conducted on a single individual to represent an accession, which may be heterogeneous with larger than expected intra-accession genetic variation. Here, we compared various population genetics parameters among six DNA (leaf) sampling methods on 90 accessions representing a wild species (O. barthii), cultivated and landraces (O. glaberrima, O. sativa), and improved varieties derived through interspecific hybridizations. A total of 1,527 DNA samples were genotyped with 46,818 polymorphic single nucleotide polymorphisms (SNPs) using DArTseq. Various statistical analyses were performed on eleven datasets corresponding to 5 plants per accession individually and in a bulk (two sets), 10 plants individually and in a bulk (two sets), all 15 plants individually (one set), and a randomly sampled individual repeated six times (six sets). Overall, we arrived at broadly similar conclusions across 11 datasets in terms of SNP polymorphism, heterozygosity/heterogeneity, diversity indices, concordance among genetic dissimilarity matrices, population structure, and genetic differentiation; there were, however, a few discrepancies between some pairs of datasets. Detailed results of each sampling method, the concordance in their outputs, and the technical and cost implications of each method were discussed.

Days to heading, controlled by the heading date genes, Hd1 and DTH8, limits rice yield-related traits in Hokkaido, Japan

A key aspect of rice breeding programs is the optimization of days to heading (DTH) for maximizing grain productivity in cultivation areas. Here, the effects of genotypes for heading date on yield-related traits in rice (culm and panicle length (CL and PL), panicle number (PN), and total number of seeds) were investigated. Heading date 1 (Hd1) and Days to heading 8 (DTH8) are the main controllers of the variation in heading date in the rice population of Hokkaido, Japan. Thus, an F₂ population (n = 192) derived from a cross between Kitaibuki (Hd1dth8) and Akage (hd1DTH8) was developed. Significant differences in DTH were found among all combinations. Each genotype for heading date showed variations in the yield-related traits without a significant difference. However, DTH exhibited high positive coefficient values (more than 0.709) with the yield-related traits except for PN, which had a negative coefficient value of -0.431. A later heading date resulted in a longer growth duration and a higher yield with a combination of longer PL and CL and lower PN. These results suggest that DTH limits the yield-related traits rather than the genotype for heading date.

Dissecting Adaptive Traits with Nested Association Mapping: Genetic Architecture of Inflorescence Morphology in Sorghum

In the cereal crop sorghum (Sorghum bicolor) inflorescence morphology variation underlies yield variation and confers adaptation across precipitation gradients, but its genetic basis is poorly understood. We characterized the genetic architecture of sorghum inflorescence morphology using a global nested association mapping (NAM) population (2200 recombinant inbred lines) and 198,000 phenotypic observations from multi-environment trials for four inflorescence morphology traits (upper branch length, lower branch length, rachis length, and rachis diameter). Trait correlations suggest that lower and upper branch length are under somewhat independent control, while lower branch length and rachis diameter are highly pleiotropic. Joint linkage and genome-wide association mapping revealed an oligogenic architecture with 1–22 QTL per trait, each explaining 0.1–5.0% of variation across the entire NAM population. There is a significant enrichment (2.twofold) of QTL colocalizing with grass inflorescence gene homologs, notably with orthologs of maize Ramosa2 and rice Aberrant Panicle Organization1 and TAWAWA1. Still, many QTL do not colocalize with inflorescence gene homologs. In global georeferenced germplasm, allelic variation at the major inflorescence QTL is geographically patterned but only weakly associated with the gradient of annual precipitation. Comparison of NAM with diversity panel association suggests that naive association models may capture some true associations not identified by mixed linear models. Overall, the findings suggest that global inflorescence diversity in sorghum is largely controlled by oligogenic, epistatic, and pleiotropic variation in ancestral regulatory networks. The findings also provide a basis for genomics-enabled breeding of locally-adapted inflorescence morphology.

Genomics of sorghum local adaptation to a parasitic plant

Host-parasite coevolution can maintain high levels of genetic diversity in traits involved in species interactions. In many systems, host traits exploited by parasites are constrained by use in other functions, leading to complex selective pressures across space and time. Here, we study genome-wide variation in the staple crop Sorghum bicolor (L.) Moench and its association with the parasitic weed Striga hermonthica (Delile) Benth., a major constraint to food security in Africa. We hypothesize that geographic selection mosaics across gradients of parasite occurrence maintain genetic diversity in sorghum landrace resistance. Suggesting a role in local adaptation to parasite pressure, multiple independent loss-of-function alleles at sorghum LOW GERMINATION STIMULANT 1 (LGS1) are broadly distributed among African landraces and geographically associated with S. hermonthica occurrence. However, low frequency of these alleles within S. hermonthica-prone regions and their absence elsewhere implicate potential trade-offs restricting their fixation. LGS1 is thought to cause resistance by changing stereochemistry of strigolactones, hormones that control plant architecture and below-ground signaling to mycorrhizae and are required to stimulate parasite germination. Consistent with trade-offs, we find signatures of balancing selection surrounding LGS1 and other candidates from analysis of genome-wide associations with parasite distribution. Experiments with CRISPR-Cas9-edited sorghum further indicate that the benefit of LGS1-mediated resistance strongly depends on parasite genotype and abiotic environment and comes at the cost of reduced photosystem gene expression. Our study demonstrates long-term maintenance of diversity in host resistance genes across smallholder agroecosystems, providing a valuable comparison to both industrial farming systems and natural communities.

Genetic Architecture of Chilling Tolerance in Sorghum Dissected with a Nested Association Mapping Population

Dissecting the genetic architecture of stress tolerance in crops is critical to understand and improve adaptation. In temperate climates, early planting of chilling-tolerant varieties could provide longer growing seasons and drought escape, but chilling tolerance (<15°) is generally lacking in tropical-origin crops. Here we developed a nested association mapping (NAM) population to dissect the genetic architecture of early-season chilling tolerance in the tropical-origin cereal sorghum (Sorghum bicolor [L.] Moench). The NAM resource, developed from reference line BTx623 and three chilling-tolerant Chinese lines, is comprised of 771 recombinant inbred lines genotyped by sequencing at 43,320 single nucleotide polymorphisms. We phenotyped the NAM population for emergence, seedling vigor, and agronomic traits (>75,000 data points from ∼16,000 plots) in multi-environment field trials in Kansas under natural chilling stress (sown 30-45 days early) and normal growing conditions. Joint linkage mapping with early-planted field phenotypes revealed an oligogenic architecture, with 5-10 chilling tolerance loci explaining 20-41% of variation. Surprisingly, several of the major chilling tolerance loci co-localize precisely with the classical grain tannin (Tan1 and Tan2) and dwarfing genes (Dw1 and Dw3) that were under strong directional selection in the US during the 20^th century. These findings suggest that chilling sensitivity was inadvertently selected due to coinheritance with desired nontannin and dwarfing alleles. The characterization of genetic architecture with NAM reveals why past chilling tolerance breeding was stymied and provides a path for genomics-enabled breeding of chilling tolerance.

Genome wide association analysis of sorghum mini core lines regarding anthracnose, downy mildew, and head smut

In previous studies, a sorghum mini core collection was scored over several years for response to Colletotrichum sublineola, Peronosclerospora sorghi, and Sporisorium reilianum, the causal agents of the disease anthracnose, downy mildew, and head smut, respectively. The screening results were combined with over 290,000 Single nucleotide polymorphic (SNP) loci from an updated version of a publicly available genotype by sequencing (GBS) dataset available for the mini core collection. GAPIT (Genome Association and Prediction Integrated Tool) R package was used to identify chromosomal locations that differ in disease response. When the top scoring SNPs were mapped to the most recent version of the published sorghum genome, in each case, a nearby and most often the closest annotated gene has precedence for a role in host defense.