Genetic structure, linkage disequilibrium and signature of selection in Sorghum: lessons from physically anchored DArT markers

Genetic loci associated with sorghum drought tolerance in multiple environments and their sensitivity to environmental covariables

KEY MESSAGE

Climate change can limit yields of naturally resilient crops, like sorghum, challenging global food security. Agriculture under an erratic climate requires tapping into a reservoir of flexible adaptive loci that can lead to lasting yield stability under multiple abiotic stress conditions. Domesticated in the hot and dry regions of Africa, sorghum is considered a harsh crop, which is adapted to important stress factors closely related to climate change. To investigate the genetic basis of drought stress adaptation in sorghum, we used a multi-environment multi-locus genome-wide association study (MEML-GWAS) in a subset of a diverse sorghum association panel (SAP) phenotyped for performance both under well-watered and water stress conditions. We selected environments in Brazil that foreshadow agriculture where both drought and temperature stresses coincide as in many tropical agricultural frontiers. Drought reduced average grain yield (Gy) by up to 50% and also affected flowering time (Ft) and plant height (Ph). We found 15 markers associated with Gy on all sorghum chromosomes except for chromosomes 7 and 9, in addition to loci associated with phenology traits. Loci associated with Gy strongly interacted with the environment in a complex way, while loci associated with phenology traits were less affected by G × E. Studying environmental covariables potentially underpinning G × E, increases in relative humidity and evapotranspiration favored and disfavored grain yield, respectively. High temperatures influenced G × E and reduced sorghum yields, with a ~ 100 kg ha-1 average decrease in grain yield for each unit increase in maximum temperature between 29 and 38 °C. Extreme G × E for sorghum stress resilience poses an additional challenge to breed crops for moving, erratic weather conditions.

Association mapping in bambara groundnut [Vigna subterranea (L.) Verdc.] reveals loci associated with agro-morphological traits

BACKGROUND

Genome-wide association studies (GWAS) are important for the acceleration of crop improvement through knowledge of marker-trait association (MTA). This report used DArT SNP markers to successfully perform GWAS on agro-morphological traits using 270 bambara groundnut [Vigna subterranea (L.) Verdc.] landraces sourced from diverse origins. The study aimed to identify marker traits association for nine agronomic traits using GWAS and their candidate genes. The experiment was conducted at two different locations laid out in alpha lattice design. The cowpea [Vigna unguiculata (L.) Walp.] reference genome (i.e. legume genome most closely related to bambara groundnut) assisted in the identification of candidate genes.

RESULTS

The analyses showed that linkage disequilibrium was found to decay rapidly with an average genetic distance of 148 kb. The broadsense heritability was relatively high and ranged from 48.39% (terminal leaf length) to 79.39% (number of pods per plant). The GWAS identified a total of 27 significant marker-trait associations (MTAs) for the nine studied traits explaining 5.27% to 24.86% of phenotypic variations. Among studied traits, the highest number of MTAs was obtained from seed coat colour (6) followed by days to flowering (5), while the least is days to maturity (1), explaining 5.76% to 11.03%, 14.5% to 19.49%, and 11.66% phenotypic variations, respectively. Also, a total of 17 candidate genes were identified, varying in number for different traits; seed coat colour (6), days to flowering (3), terminal leaf length (2), terminal leaf width (2), number of seed per pod (2), pod width (1) and days to maturity (1).

CONCLUSION

These results revealed the prospect of GWAS in identification of SNP variations associated with agronomic traits in bambara groundnut. Also, its present new opportunity to explore GWAS and marker assisted strategies in breeding of bambara groundnut for acceleration of the crop improvement.

Exploring the genetic diversity and population structure of upland cotton germplasm by iPBS-retrotransposons markers

BACKGROUND

Upland cotton is one of the utmost significant strategic fiber crops, and play a vital role in the global textile industry.

METHODS AND RESULTS

A total of 128 genotypes comprised Gossypium hirsutum L, Gossypium barbadense L., and pure lines were used to examine genetic diversity using iPBS-retrotransposon markers system. Eleven highly polymorphic primers yielded 287 bands and 99.65% polymorphism was recorded. The mean polymorphism information content was estimated at 0.297 and the average diversity indices for the effective number of alleles, Shannon's information index, and overall gene diversity were 1.481, 0.443, and 0.265, respectively. The analysis of molecular variance (AMOVA) revealed that 69% of the genetic variation was within the population. A model-based STRUCTURE algorithm divided the entire germplasm into four populations and one un-classified population, the genotypes G42 (originating in Egypt) and G128 (originating in the United States), showed the highest genetic distance (0.996) so these genotypes could be suggested for breeding programs as parental lines.

CONCLUSIONS

This is the first investigation using an iPBS-retrotransposon marker system to examine the genetic diversity and population structure of upland cotton germplasm. The rich diversity found in upland cotton germplasm could be exploited as a genetic resource when developing breeding programs and could also help with efforts to breed cotton around the world. These findings also show the applicability and effectiveness of iPBS-retrotransposons for the molecular characterization of cotton germplasm.

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.

Deciphering the genetic diversity and population structure of Turkish bread wheat germplasm using iPBS-retrotransposons markers

BACKGROUND

Research activities aiming to investigate the genetic diversity are very crucial because they provide information for the breeding and germplasm conservation activities. Wheat is one of the most important cereal crops globally by feeding more than a third of the human population around the world.

METHODS AND RESULTS

During present investigation, a total of 74 Turkish bread wheat accessions (54 landraces and 20 cultivars) were used as plant material and iPBS-retrotransposons marker system was used for the molecular characterization. 13 polymorphic primers used for molecular characterization resulted a total of 152 bands. Range of calculated diversity indices like polymorphism information content (0.11-0.702), effective numbers of alleles (1.026-1.526), Shannon's information index (0.101-0.247) and gene diversity (0.098-0.443) confirmed higher genetic variations in studied germplasm. Bread wheat landraces reflected higher genetic variations compared to commercial cultivars. Analysis of molecular variance resulted that higher (98%) genetic variations are present within populations. The model-based structure algorithm separated 74 bread wheat accessions in to two populations. Diversity indices based on structure evaluated population's revealed population B as a more diverse population. The principal coordinate analysis and neighbor-joining analysis separated 74 bread wheat accessions according to their collection points. Genetic distance for 74 Turkish bread wheat accessions explored Bingol and Asure accessions as genetically diverse that can be used as parents for breeding activities.

CONCLUSIONS

The extensive diversity of bread wheat in Turkish germplasm might be used as genetic resource for the exhaustive wheat breeding program. For instance, accessions Bingol and Asure were found genetically diverse and can be used as parents for future breeding activities.

A genomics resource for genetics, physiology, and breeding of West African sorghum

Local landrace and breeding germplasm is a useful source of genetic diversity for regional and global crop improvement initiatives. Sorghum (Sorghum bicolor L. Moench) in western Africa (WA) has diversified across a mosaic of cultures and end uses and along steep precipitation and photoperiod gradients. To facilitate germplasm utilization, a West African sorghum association panel (WASAP) of 756 accessions from national breeding programs of Niger, Mali, Senegal, and Togo was assembled and characterized. Genotyping-by-sequencing (GBS) was used to generate 159,101 high-quality biallelic single nucleotide polymorphisms (SNPs), with 43% in intergenic regions and 13% in genic regions. High genetic diversity was observed within the WASAP (π = .00045), only slightly less than in a global diversity panel (GDP) (π = .00055). Linkage disequilibrium (LD) decayed to background level (r² < .1) by ∼50 kb in the WASAP. Genome-wide diversity was structured both by botanical type and by populations within botanical type with eight ancestral populations identified. Most populations were distributed across multiple countries, suggesting several potential common gene pools across the national programs. Genome-wide association studies (GWAS) of days to flowering (DFLo) and plant height (PH) revealed eight and three significant quantitative trait loci (QTL), respectively, with major height QTL at canonical height loci Dw3 and SbHT7.1. Colocalization of two of eight major flowering time QTL with flowering genes previously described in U.S. germplasm (Ma6 and SbCN8) suggests that photoperiodic flowering in West African sorghum is conditioned by both known and novel genes. This genomic resource provides a foundation for genomics-enabled breeding of climate-resilient varieties in WA.

Meta-Analysis Identifies Pleiotropic Loci Controlling Phenotypic Trade-offs in Sorghum

Community association populations are composed of phenotypically and genetically diverse accessions. Once these populations are genotyped, the resulting marker data can be reused by different groups investigating the genetic basis of different traits. Because the same genotypes are observed and scored for a wide range of traits in different environments, these populations represent a unique resource to investigate pleiotropy. Here we assembled a set of 234 separate trait datasets for the Sorghum Association Panel, a group of 406 sorghum genotypes widely employed by the sorghum genetics community. Comparison of genome wide association studies conducted with two independently generated marker sets for this population demonstrate that existing genetic marker sets do not saturate the genome and likely capture only 35-43% of potentially detectable loci controlling variation for traits scored in this population. While limited evidence for pleiotropy was apparent in cross-GWAS comparisons, a multivariate adaptive shrinkage approach recovered both known pleiotropic effects of existing loci and new pleiotropic effects, particularly significant impacts of known dwarfing genes on root architecture. In addition, we identified new loci with pleiotropic effects consistent with known trade-offs in sorghum development. These results demonstrate the potential for mining existing trait datasets from widely used community association populations to enable new discoveries from existing trait datasets as new, denser genetic marker datasets are generated for existing community association populations.

Genomic patterns of structural variation among diverse genotypes of Sorghum bicolor and a potential role for deletions in local adaptation

Genomic structural mutations, especially deletions, are an important source of variation in many species and can play key roles in phenotypic diversification and evolution. Previous work in many plant species has identified multiple instances of structural variations (SVs) occurring in or near genes related to stress response and disease resistance, suggesting a possible role for SVs in local adaptation. Sorghum [Sorghum bicolor (L.) Moench] is one of the most widely grown cereal crops in the world. It has been adapted to an array of different climates as well as bred for multiple purposes, resulting in a striking phenotypic diversity. In this study, we identified genome-wide SVs in the Biomass Association Panel, a collection of 347 diverse sorghum genotypes collected from multiple countries and continents. Using Illumina-based, short-read whole-genome resequencing data from every genotype, we found a total of 24,648 SVs, including 22,359 deletions. The global site frequency spectrum of deletions and other types of SVs fit a model of neutral evolution, suggesting that the majority of these mutations were not under any types of selection. Clustering results based on single nucleotide polymorphisms separated the genotypes into eight clusters which largely corresponded with geographic origins, with many of the large deletions we uncovered being unique to a single cluster. Even though most deletions appeared to be neutral, a handful of cluster-specific deletions were found in genes related to biotic and abiotic stress responses, supporting the possibility that at least some of these deletions contribute to local adaptation in sorghum.

Genetic diversity, population structure and selection signature in Ethiopian Sorghum (Sorghum bicolor L. [Moench]) germplasm

Ethiopia, the probable center of origin and diversity for sorghum [Sorghum bicolor L. (Moench)] and with unique ecogeographic features, possesses a large number of sorghum landraces that have not been well studied. Increased knowledge of this diverse germplasm through large-scale genomic characterization may contribute for understanding of evolutionary biology, and adequate use of these valuable resources from the center of origin. In this study, we characterized genetic diversity, population structure and selection signature in 304 sorghum accessions collected from diverse sorghum growing regions of Ethiopia using genotyping-by-sequencing. We identified a total of 108,107 high-quality single-nucleotide polymorphism (SNPs) markers that were evenly distributed across the sorghum genome. The average gene diversity among accessions was high (H_e = 0.29). We detected a relatively low frequency of rare alleles (26%), highlighting the potential of this germplasm for subsequent allele mining studies through genome-wide association studies. Although we found no evidence of genetic differentiation among administrative regions (F_ST = 0.02, P = 0.12), population structure and cluster analyses showed clear differentiation among six Ethiopian sorghum populations (F_ST = 0.28, P = 0.01) adapting to different environments. Analysis of SNP differentiation between the identified genetic groups revealed a total of 40 genomic regions carrying signatures of selection. These regions harbored candidate genes potentially involved in a variety of biological processes, including abiotic stress tolerance, pathogen defense and reproduction. Overall, a high level of untapped diversity for sorghum improvement remains available in Ethiopia, with patterns of diversity consistent with divergent selection on a range of adaptive characteristics.

In-Depth Genetic Diversity and Population Structure of Endangered Peruvian Amazon Rosewood Germplasm Using Genotyping by Sequencing (GBS) Technology

Research studies on conservative genetics of endangered plants are very important to establish the management plans for the conservation of biodiversity. Rosewood is an evergreen tree of the Amazon region and its essential oil has great acceptance in the medical and cosmetic industry. The present study aimed to explore the genetic diversity and population structure of 90 rosewood accessions collected from eight localities of Peruvian Amazon territory through DArTseq markers. A total of 7485 informative markers resulted from genotyping by sequencing (GBS) analysis were used for the molecular characterization of rosewood germplasm. Mean values of various calculated diversity parameters like observed number of alleles (1.962), the effective number of alleles (1.669), unbiased expected heterozygosity (0.411), and percent polymorphism (93.51%) over the entire germplasm showed the existence of a good level of genetic variations. Our results showed that the Mairiricay population was more diverse compared to the rest of the populations. Tamshiyacu-2 and Mairiricay-15 accessions were found genetically distinct accessions. The analysis of molecular variance (AMOVA) reflected maximum variations (75%) are due to differences within populations. The implemented clustering algorithms, i.e., STRUCTURE, neighbor-joining analysis and principal coordinate analysis (PCoA) separated the studied germplasm on the basis of their geographical locations. Diversity indices for STRUCTURE-based populations showed that subpopulation A is more diverse population than the rest of the populations, for such reason, individuals belonging to this subpopulation should be used for reintroduction or reinforcement plans of rosewood conservation. We envisage that molecular characterization of Peruvian rosewood germplasm with DArTseq markers will provide a platform for the conservation, management and restoration of endangered rosewood in upcoming years.

Association mapping and genomic selection for sorghum adaptation to tropical soils of Brazil in a sorghum multiparental random mating population

A multiparental random mating population used in sorghum breeding is amenable for the detection of QTLs related to tropical soil adaptation, fine mapping of underlying genes and genomic selection approaches. Tropical soils where low phosphorus (P) and aluminum (Al) toxicity limit sorghum [Sorghum bicolor (L.) Moench] production are widespread in the developing world. We report on BRP13R, a multiparental random mating population (MP-RMP), which is commonly used in sorghum recurrent selection targeting tropical soil adaptation. Recombination dissipated much of BRP13R's likely original population structure and average linkage disequilibrium (LD) persisted up to 2.5 Mb, establishing BRP13R as a middle ground between biparental populations and sorghum association panels. Genome-wide association mapping (GWAS) identified conserved QTL from previous studies, such as for root morphology and grain yield under low-P, and indicated the importance of dominance in the genetic architecture of grain yield. By overlapping consensus QTL regions, we mapped two candidate P efficiency genes to a ~ 5 Mb region on chromosomes 6 (ALMT) and 9 (PHO2). Remarkably, we find that only 200 progeny genotyped with ~ 45,000 markers in BRP13R can lead to GWAS-based positional cloning of naturally rare, subpopulation-specific alleles, such as for SbMATE-conditioned Al tolerance. Genomic selection was found to be useful in such MP-RMP, particularly if markers in LD with major genes are fitted as fixed effects into GBLUP models accommodating dominance. Shifts in allele frequencies in progeny contrasting for grain yield indicated that intermediate to minor-effect genes on P efficiency, such as SbPSTOL1 genes, can be employed in pre-breeding via allele mining in the base population. Therefore, MP-RMPs such as BRP13R emerge as multipurpose resources for efficient gene discovery and deployment for breeding sorghum cultivars adapted to tropical soils.

Genetic Diversity, Population Structure and Marker-Trait Association for 100-Seed Weight in International Safflower Panel Using SilicoDArT Marker Information

Safflower is an important oilseed crop mainly grown in the arid and semi-arid regions of the world. The aim of this study was to explore phenotypic and genetic diversity, population structure, and marker-trait association for 100-seed weight in 94 safflower accessions originating from 26 countries using silicoDArT markers. Analysis of variance revealed statistically significant genotypic effects (p < 0.01), while Turkey samples resulted in higher 100-seed weight compared to Pakistan samples. A Constellation plot divided the studied germplasm into two populations on the basis of their 100-seed weight. Various mean genetic diversity parameters including observed number of alleles (1.99), effective number of alleles (1.54), Shannon's information index (0.48), expected heterozygosity (0.32), and unbiased expected heterozygosity (0.32) for the entire population exhibited sufficient genetic diversity using 12232 silicoDArT markers. Analysis of molecular variance (AMOVA) revealed that most of the variations (91%) in world safflower panel are due to differences within country groups. A model-based structure grouped the 94 safflower accessions into populations A, B, C and an admixture population upon membership coefficient. Neighbor joining analysis grouped the safflower accessions into two populations (A and B). Principal coordinate analysis (PCoA) also clustered the safflower accessions on the basis of geographical origin. Three accessions; Egypt-5, Egypt-2, and India-2 revealed the highest genetic distance and hence might be recommended as candidate parental lines for safflower breeding programs. The mixed linear model i.e., the Q + K model, demonstrated that two DArTseq markers (DArT-45483051 and DArT-15672391) had significant association (p < 0.01) for 100-seed weight. We envisage that identified DArTseq markers associated with 100-seed weight will be helpful to develop high-yielding cultivars of safflower through marker-assisted breeding in the near future.

Whole Genome Scan Reveals Molecular Signatures of Divergence and Selection Related to Important Traits in Durum Wheat Germplasm

The first breeding program in the world for durum wheat was conceived in Italy in the early 1900s. Over the decades, pressure exerted by natural and artificial selection could have progressively reduced the genetic diversity of the durum wheat germplasm. In the present study, a large panel of Italian durum wheat accessions that includes landraces, old and modern cultivars was subjected to genotyping using the Illumina iSelect 15K wheat SNP array. The aim was to assess the impact that selection has in shaping Italian durum wheat genetic diversity and to exploit the patterns of genetic diversity between populations to identify molecular signatures of divergence and selection. Relatively small differences in genetic diversity have been observed among accessions, which have been selected and cultivated in Italy over the past 150 years. Indeed, directional selection combined with that operated by farmers/breeders resulted in the increase of linkage disequilibrium (LD) and in changes of the allelic frequencies in DNA regions that control important agronomic traits. Results from this study also show that major well-known genes and/or QTLs affecting plant height (RHT), earliness (VRN, PPD) and grain quality (GLU, PSY, PSD, LYC, PPO, LOX3) co-localized with outlier SNP loci. Interestingly, many of these SNPs fall in genomic regions where genes involved in nitrogen metabolism are. This finding highlights the key role these genes have played in the transition from landraces to modern cultivars. Finally, our study remarks on the need to fully exploit the genetic diversity of Italian landraces by intense pre-breeding activities aimed at introducing a new source of adaptability and resistance in the genetic background of modern cultivars, to contrast the effect of climate change. The list of divergent loci and loci under selection associated with useful agronomic traits represents an invaluable resource to detect new allelic variants for target genes and for guiding new genomic selection programs in durum wheat.

Genetic Diversity, Population Structure, and Linkage Disequilibrium of Pearl Millet

Mapping of GBS reads of 398 accessions to the draft genome sequence identified 82,112 SNPs Model-based clustering analysis revealed a hierarchical genetic structure of six subgroups Greater LD decay in the west-African subpopulation is likely due to long history of recombination Genetic differentiation analysis among subpopulations revealed variation in selection signatures Pearl millet [Cenchrus americanus (L.) Morrone syn. Pennisetum glaucum (L.) R. Br.] is one of the most extensively cultivated cereals in the world, after wheat (Triticum aestivum L.), maize (Zea mays L.), rice (Oryza sativa L.), barley (Hordeum vulgare L.), and sorghum [Sorghum bicolor (L.) Moench]. It is the main component of traditional farming systems and a staple food in the arid and semiarid regions of Africa and southern Asia. However, its genetic improvement is lagging behind other major cereals and the yield is still low. Genotyping-by-sequencing (GBS)-based single-nucleotide polymorphism (SNP) markers were screened on a total of 398 accessions from different geographic regions to assess genetic diversity, population structure, and linkage disequilibrium (LD). By mapping the GBS reads to the reference genome sequence, 82,112 genome-wide SNPs were discovered. The telomeric regions of the chromosomes have the higher SNP density than in pericentromeric regions. Model-based clustering analysis of the population revealed a hierarchical genetic structure of six subgroups that mostly overlap with the geographic origins or sources of the genotypes but with differing levels of admixtures. A neighbor-joining phylogeny analysis revealed that germplasm from western Africa rooted the dendrogram with much diversity within each subgroup. Greater LD decay was observed in the west-African subpopulation than in the other subpopulations, indicating a long history of recombination among landraces. Also, genome scan of genetic differentiatation detected different selection histories among subpopulations. These results have potential application in the development of genomic-assisted breeding in pearl millet and heterotic grouping of the lines for improved hybrid performance.

New candidate loci and marker genes on chromosome 7 for improved chilling tolerance in sorghum

Sorghum is often exposed to suboptimal low temperature stress under field conditions, particularly at the seedling establishment stage. Enhancing chilling tolerance will facilitate earlier planting and so minimize the negative impacts of other stresses experienced at later growth stages. Genome-wide association mapping was performed on a sorghum association panel grown under control (30/20 °C; day/night) and chilling (20/10 °C) conditions. Genomic regions on chromosome 7, controlling the emergence index and seedling (root and shoot) vigor, were associated with increased chilling tolerance but they did not co-localize with undesirable tannin content quantitative trait loci (QTLs). Shoot and root samples from highly contrasting haplotype pairs expressing differential responses to chilling stress were used to identify candidate genes. Three candidate genes (an alpha/beta hydrolase domain protein, a DnaJ/Hsp40 motif-containing protein, and a YTH domain-containing RNA-binding protein) were expressed at significantly higher levels under chilling stress in the tolerant haplotype compared with the sensitive haplotype and BTx623. Moreover, two CBF/DREB1A transcription factors on chromosome 2 showed a divergent response to chilling in the contrasting haplotypes. These studies identify haplotype differences on chromosome 7 that modulate chilling tolerance by either regulating CBF or feeding back into this signaling pathway. We have identified new candidate genes that will be useful markers in ongoing efforts to develop tannin-free chilling-tolerant sorghum hybrids.

Genotyping by sequencing provides new insights into the diversity of Napier grass (Cenchrus purpureus) and reveals variation in genome-wide LD patterns between collections

Napier grass is an important tropical forage-grass and of growing potential as an energy crop. One-hundred-five Napier grass accessions, encompassing two independent collections, were subjected to genotyping by sequencing which generated a set of high-density genome-wide markers together with short sequence reads. The reads, averaging 54 nucleotides, were mapped to the pearl millet genome and the closest genes and annotation information were used to select candidate genes linked to key forage traits. 980 highly polymorphic SNP markers, distributed across the genome, were used to assess population structure and diversity with seven-subgroups identified. A few representative accessions were selected with the objective of distributing subsets of a manageable size for further evaluation. Genome-wide linkage disequilibrium (LD) analyses revealed a fast LD-decay, on average 2.54 kbp, in the combined population with a slower LD-decay in the ILRI collection compared with the EMBRAPA collection, the significance of which is discussed. This initiative generated high-density markers with a good distribution across the genome. The diversity analysis revealed the existence of a substantial amount of variation in the ILRI collection and identified some unique materials from the EMBRAPA collection, demonstrating the potential of the overall population for further genetic and marker-trait-association studies.

Genomic signatures of adaptation to Sahelian and Soudanian climates in sorghum landraces of Senegal

Uncovering the genomic basis of climate adaptation in traditional crop varieties can provide insight into plant evolution and facilitate breeding for climate resilience. In the African cereal sorghum (Sorghum bicolor L. [Moench]), the genomic basis of adaptation to the semiarid Sahelian zone versus the subhumid Soudanian zone is largely unknown. To address this issue, we characterized a large panel of 421 georeferenced sorghum landrace accessions from Senegal and adjacent locations at 213,916 single-nucleotide polymorphisms (SNPs) using genotyping-by-sequencing. Seven subpopulations distributed along the north-south precipitation gradient were identified. Redundancy analysis found that climate variables explained up to 8% of SNP variation, with climate collinear with space explaining most of this variation (6%). Genome scans of nucleotide diversity suggest positive selection on chromosome 2, 4, 5, 7, and 10 in durra sorghums, with successive adaptation during diffusion along the Sahel. Putative selective sweeps were identified, several of which colocalize with stay-green drought tolerance (Stg) loci, and a priori candidate genes for photoperiodic flowering and inflorescence morphology. Genome-wide association studies of photoperiod sensitivity and panicle compactness identified 35 and 13 associations that colocalize with a priori candidate genes, respectively. Climate-associated SNPs colocalize with Stg3a, Stg1, Stg2, and Ma6 and have allelic distribution consistent with adaptation across Sahelian and Soudanian zones. Taken together, the findings suggest an oligogenic basis of adaptation to Sahelian versus Soudanian climates, underpinned by variation in conserved floral regulatory pathways and other systems that are less understood in cereals.

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

Mining crop genomic variation can facilitate the genetic research of complex traits and molecular breeding. In sorghum [ L. (Moench)], several large-scale single nucleotide polymorphism (SNP) datasets have been generated using genotyping-by-sequencing of KI reduced representation libraries. However, data reuse has been impeded by differences in reference genome coordinates among datasets. To facilitate reuse of these data, we constructed and characterized an integrated 459,304-SNP dataset for 10,323 sorghum genotypes on the version 3.1 reference genome. The SNP distribution showed high enrichment in subtelomeric chromosome arms and in genic regions (48% of SNPs) and was highly correlated ( = 0.82) to the distribution of KI restriction sites. The genetic structure reflected population differences by botanical race, as well as familial structure among recombinant inbred lines (RILs). Faster linkage disequilibrium decay was observed in the diversity panel than in the RILs, as expected, given the greater opportunity for recombination in diverse populations. To validate the quality and utility of the integrated SNP dataset, we used genome-wide association studies (GWAS) of genebank phenotype data, precisely mapping several known genes (e.g and ) and identifying novel associations for other traits. We further validated the dataset with GWAS of new and published plant height and flowering time data in a nested association mapping population, precisely mapping known genes and identifying epistatic interactions underlying both traits. These findings validate this integrated SNP dataset as a useful genomics resource for sorghum genetics and breeding.

Mobile genomic element diversity in world collection of safflower (Carthamus tinctorius L.) panel using iPBS-retrotransposon markers

Safflower (Carthamus tinctorius L.) is a multipurpose crop of dry land yielding very high quality of edible oil. Present study was aimed to investigate the genetic diversity and population structure of 131 safflower accessions originating from 28 different countries using 13 iPBS-retrotransposon markers. A total of 295 iPBS bands were observed among which 275 (93.22%) were found polymorphic. Mean Polymorphism information content (0.48) and diversity parameters including mean effective number of alleles (1.33), mean Shannon's information index (0.33), overall gene diversity (0.19), Fstatistic (0.21), and inbreeding coefficient (1.00) reflected the presence of sufficient amount of genetic diversity in the studied plant materials. Analysis of molecular variance (AMOVA) showed that more than 40% of genetic variation was derived from populations. Model-based structure, principal coordinate analysis (PCoA) and unweighted pair-group method with arithmetic means (UPGMA) algorithms clustered the 131 safflower accessions into four main populations A, B, C, D and an unclassified population, with no meaningful geographical origin. Most diverse accessions originated from Asian countries including Afghanistan, Pakistan, China, Turkey, and India. Four accessions, Turkey3, Afghanistan4, Afghanistan2, and Pakistan24 were found most genetically distant and might be recommended as a candidate parents for breeding purposes. The findings of this study are most probably supported by the seven similarity centers hypothesis of safflower. This is a first study to explore the genetic diversity and population structure in safflower accessions using the iPBS-retrotransposon markers. The information provided in this work will therefore be helpful for scientists interested in safflower breeding.

Semiautomated Feature Extraction from RGB Images for Sorghum Panicle Architecture GWAS

Because structural variation in the inflorescence architecture of cereal crops can influence yield, it is of interest to identify the genes responsible for this variation. However, the manual collection of inflorescence phenotypes can be time consuming for the large populations needed to conduct genome-wide association studies (GWAS) and is difficult for multidimensional traits such as volume. A semiautomated phenotyping pipeline, TIM (Toolkit for Inflorescence Measurement), was developed and used to extract unidimensional and multidimensional features from images of 1,064 sorghum (Sorghum bicolor) panicles from 272 genotypes comprising a subset of the Sorghum Association Panel. GWAS detected 35 unique single-nucleotide polymorphisms associated with variation in inflorescence architecture. The accuracy of the TIM pipeline is supported by the fact that several of these trait-associated single-nucleotide polymorphisms (TASs) are located within chromosomal regions associated with similar traits in previously published quantitative trait locus and GWAS analyses of sorghum. Additionally, sorghum homologs of maize (Zea mays) and rice (Oryza sativa) genes known to affect inflorescence architecture are enriched in the vicinities of TASs. Finally, our TASs are enriched within genomic regions that exhibit high levels of divergence between converted tropical lines and cultivars, consistent with the hypothesis that these chromosomal intervals were targets of selection during modern breeding.

Evaluation of Genetic Variation among Sorghum Varieties from Southwest China via Genome Resequencing

Little is known regarding genomic variation among glutinous sorghum [ (L.) Moench] varieties grown in southwest China, which are primarily used to brew the popular Jiang-flavor liquor. This study evaluated genomic variation among six representative sorghum accessions via whole-genome resequencing. The evaluation revealed 2365,363 single-nucleotide polymorphisms (SNPs), 394,365 insertions and deletions, and 47,567 copy number variations among the six genomes. Chromosomes 5 and 10 showed relatively high SNP densities, whereas whole-genome diversity in this population was low. In addition, some chromosomal loci exhibited obvious selection during the breeding process. Sorghum accessions from southwest China formed an elite germplasm population compared with the findings of other geographic populations, and the elite variety 'Hongyingzi' contained 79 unique genes primarily involved in basic metabolism. The six sorghum lines contained a large number of high-confidence genes, with Hongyingzi in particular possessing 104 unique genes. These findings advance our understanding of domestication of the sorghum genome, and Chinese sorghum accessions will be valuable resources for further research and breeding improvements.

Characterization of genetic diversity in Turkish common bean gene pool using phenotypic and whole-genome DArTseq-generated silicoDArT marker information

Turkey presents a great diversity of common bean landraces in farmers' fields. We collected 183 common bean accessions from 19 different Turkish geographic regions and 5 scarlet runner bean accessions to investigate their genetic diversity and population structure using phenotypic information (growth habit, and seed weight, flower color, bracteole shape and size, pod shape and leaf shape and color), geographic provenance and 12,557 silicoDArT markers. A total of 24.14% markers were found novel. For the entire population (188 accessions), the expected heterozygosity was 0.078 and overall gene diversity, Fst and Fis were 0.14, 0.55 and 1, respectively. Using marker information, model-based structure, principal coordinate analysis (PCoA) and unweighted pair-group method with arithmetic means (UPGMA) algorithms clustered the 188 accessions into two main populations A (predominant) and B, and 5 unclassified genotypes, representing 3 meaningful heterotic groups for breeding purposes. Phenotypic information clearly distinguished these populations; population A and B, respectively, were bigger (>40g/100 seeds) and smaller (<40g/100 seeds) seed-sized. The unclassified population was pure and only contained climbing genotypes with 100 seed weight 2-3 times greater than populations A and B. Clustering was mainly based on A: seed weight, B: growth habit, C: geographical provinces and D: flower color. Mean kinship was generally low, but population B was more diverse than population A. Overall, a useful level of gene and genotypic diversity was observed in this work and can be used by the scientific community in breeding efforts to develop superior common bean strains.

Genetic architecture of photosynthesis in Sorghum bicolor under non-stress and cold stress conditions

Sorghum (Sorghum bicolor L. Moench) is a C4 species sensitive to the cold spring conditions that occur at northern latitudes, especially when coupled with excessive light, and that greatly affect the photosynthetic rate. The objective of this study was to discover genes/genomic regions that control the capacity to cope with excessive energy under low temperature conditions during the vegetative growth period. A genome-wide association study (GWAS) was conducted for seven photosynthetic gas exchange and chlorophyll fluorescence traits under three consecutive temperature treatments: control (28 °C/24 °C), cold (15 °C/15 °C), and recovery (28 °C/24 °C). Cold stress significantly reduced the rate of photosynthetic CO2 uptake of sorghum plants, and a total of 143 unique genomic regions were discovered associated with at least one trait in a particular treatment or with derived variables. Ten regions on chromosomes 3, 4, 6, 7, and 8 that harbor multiple significant markers in linkage disequilibrium (LD) were consistently identified in gas exchange and chlorophyll fluorescence traits. Several candidate genes within those intervals have predicted functions related to carotenoids, phytohormones, thioredoxin, components of PSI, and antioxidants. These regions represent the most promising results for future validation and with potential application for the improvement of crop productivity under cold stress.

Genome-Wide Association Study of Developing Leaves' Heat Tolerance during Vegetative Growth Stages in a Sorghum Association Panel

Heat stress reduces grain yield and quality worldwide. Enhancing heat tolerance of crops at all developmental stages is one of the essential strategies required for sustaining agricultural production especially as frequency of temperature extremes escalates in response to climate change. Although heat tolerance mechanisms have been studied extensively in model plant species, little is known about the genetic control underlying heat stress responses of crop plants at the vegetative stage under field conditions. To dissect the genetic basis of heat tolerance in sorghum [ (L.) Moench], we performed a genome-wide association study (GWAS) for traits responsive to heat stress at the vegetative stage in an association panel. Natural variation in leaf firing (LF) and leaf blotching (LB) were evaluated separately for 3 yr in experimental fields at three locations where sporadic heat waves occurred throughout the sorghum growing season. We identified nine single-nucleotide polymorphisms (SNPs) that were significantly associated with LF and five SNPs that were associated with LB. Candidate genes near the SNPs were investigated and 14 were directly linked to biological pathways involved in plant stress responses including heat stress response. The findings of this study provide new knowledge on the genetic control of leaf traits responsive to heat stress in sorghum, which could aid in elucidating the genetic and molecular mechanisms of vegetative stage heat tolerance in crops. The results also provide candidate markers for molecular breeding of enhanced heat tolerance in cereal and bioenergy crops.

DArT whole genome profiling provides insights on the evolution and taxonomy of edible Banana (Musa spp.)

BACKGROUND AND AIMS

Dessert and cooking bananas are vegetatively propagated crops of great importance for both the subsistence and the livelihood of people in developing countries. A wide diversity of diploid and triploid cultivars including AA, AB, AS, AT, AAA, AAB, ABB, AAS and AAT genomic constitutions exists. Within each of this genome groups, cultivars are classified into subgroups that are reported to correspond to varieties clonally derived from each other after a single sexual event. The number of those founding events at the basis of the diversity of bananas is a matter of debate.

METHODS

We analysed a large panel of 575 accessions, 94 wild relatives and 481 cultivated accessions belonging to the section Musa with a set of 498 DArT markers previously developed.

KEY RESULTS

DArT appeared successful and accurate to describe Musa diversity and help in the resolution of cultivated banana genome constitution and taxonomy, and highlighted discrepancies in the acknowledged classification of some accessions. This study also argues for at least two centres of domestication corresponding to South-East Asia and New Guinea, respectively. Banana domestication in New Guinea probably followed different schemes that those previously reported where hybridization underpins the emergence of edible banana. In addition, our results suggest that not all wild ancestors of bananas are known, especially in M. acuminata subspecies. We also estimate the extent of the two consecutive bottlenecks in edible bananas by evaluating the number of sexual founding events underlying our sets of edible diploids and triploids, respectively.

CONCLUSIONS

The attribution of clone identity to each sample of the sets allowed the detection of subgroups represented by several sets of clones. Although morphological characterization of some of the accessions is needed to correct potentially erroneous classifications, some of the subgroups seem polyclonal.

Genome-Wide Association Studies of Grain Yield Components in Diverse Sorghum Germplasm

Grain yield and its primary determinants, grain number and weight, are important traits in cereal crops that have been well studied; however, the genetic basis of and interactions between these traits remain poorly understood. Characterization of grain yield per primary panicle (YPP), grain number per primary panicle (GNP), and 1000-grain weight (TGW) in sorghum [ (L.) Moench], a hardy C cereal with a genome size of ∼730 Mb, was implemented in a diversity panel containing 390 accessions. These accessions were genotyped to obtain 268,830 single-nucleotide polymorphisms (SNPs). Genome-wide association studies (GWAS) were performed to identify loci associated with each grain yield component and understand the genetic interactions between these traits. Genome-wide association studies identified associations across the genome with YPP, GNP, and TGW that were located within previously mapped sorghum QTL for panicle weight, grain yield, and seed size, respectively. There were no significant associations between GNP and TGW that were within 100 kb, much greater than the average linkage disequilibrium (LD) in sorghum. The identification of nonoverlapping loci for grain number and weight suggests these traits may be manipulated independently to increase the grain yield of sorghum. Following GWAS, genomic regions surrounding each associated SNP were mined for candidate genes. Previously published expression data indicated several TGW candidate genes, including an ethylene receptor homolog, were primarily expressed within developing seed tissues to support GWAS. Furthermore, maize ( L.) homologs of identified TGW candidates were differentially expressed within the seed between small- and large-kernel lines from a segregating maize population.

Molecular Approaches to Understand Nutritional Potential of Coarse Cereals

Coarse grains are important group of crops that constitutes staple food for large population residing primarily in the arid and semi-arid regions of the world. Coarse grains are designated as nutri-cereals as they are rich in essential amino acids, minerals and vitamins. In spite of having several nutritional virtues in coarse grain as mentioned above, there is still scope for improvement in quality parameters such as cooking qualities, modulation of nutritional constituents and reduction or elimination of anti-nutritional factors. Besides its use in traditional cooking, coarse grains have been used mainly in the weaning food preparation and other malted food production. Improvement in quality parameters will certainly increase consumer's preference for coarse grains and increase their demand. The overall genetic gain in quality traits of economic importance in the cultivated varieties will enhance their industrial value and simultaneously increase income of farmers growing these varieties. The urgent step for improvement of quality traits in coarse grains requires a detailed understanding of molecular mechanisms responsible for varied level of different nutritional contents in different genotypes of these crops. In this review we have discussed the progresses made in understanding of coarse grain biology with various omics tool coupled with modern breeding approaches and the current status with regard to our effort towards dissecting traits related to improvement of quality and nutritional constituents of grains.

Association mapping of germinability and seedling vigor in sorghum under controlled low-temperature conditions

Sorghum is one of the world's most important food, feed, and fiber crops as well as a potential feedstock for lignocellulosic bioenergy. Early-season planting extends sorghum's growing season and increases yield in temperate regions. However, sorghum's sensitivity to low soil temperatures adversely impacts seed germination. In this study, we evaluated the 242 accessions of the ICRISAT sorghum mini core collection for seed germination and seedling vigor at 12 °C as a measure of cold tolerance. Genome-wide association analysis was performed with approximately 162,177 single nucleotide polymorphism markers. Only one marker locus (Locus 7-2) was significantly associated with low-temperature germination and none with vigor. The linkage of Locus 7-2 to low-temperature germination was supported by four lines of evidence: strong association in three independent experiments, co-localization with previously mapped cold tolerance quantitative trait loci (QTL) in sorghum, a candidate gene that increases cold tolerance and germination rate when its wheat homolog is overexpressed in tobacco, and its syntenic region in rice co-localized with two cold tolerance QTL in rice. This locus may be useful in developing tools for molecular breeding of sorghums with improved low-temperature germinability.

Identification of Genetic Differentiation between Waxy and Common Maize by SNP Genotyping

Waxy maize (Zea mays L. var. ceratina) is an important vegetable and economic crop that is thought to have originated from cultivated flint maize and most recently underwent divergence from common maize. In this study, a total of 110 waxy and 110 common maize inbred lines were genotyped with 3072 SNPs to evaluate the genetic diversity, population structure, and linkage disequilibrium decay as well as identify putative loci that are under positive selection. The results revealed abundant genetic diversity in the studied panel and that genetic diversity was much higher in common than in waxy maize germplasms. Principal coordinate analysis and neighbor-joining cluster analysis consistently classified the 220 accessions into two major groups and a mixed group with mixed ancestry. Subpopulation structure in both waxy and common maize sets were associated with the germplasm origin and corresponding heterotic groups. The LD decay distance (1500–2000 kb) in waxy maize was lower than that in common maize. Fourteen candidate loci were identified as under positive selection between waxy and common maize at the 99% confidence level. The information from this study can assist waxy maize breeders by enhancing parental line selection and breeding program design.

Evaluating the Influence of the Microsatellite Marker Set on the Genetic Structure Inferred in Pyrus communis L

Fingerprinting information can be used to elucidate in a robust manner the genetic structure of germplasm collections, allowing a more rational and fine assessment of genetic resources. Bayesian model-based approaches are nowadays majorly preferred to infer genetic structure, but it is still largely unresolved how marker sets should be built in order to obtain a robust inference. The objective was to evaluate, in Pyrus germplasm collections, the influence of the SSR marker set size on the genetic structure inferred, also evaluating the influence of the criterion used to select those markers. Inferences were performed considering an increasing number of SSR markers that ranged from just two up to 25, incorporated one at a time into the analysis. The influence of the number of SSR markers used was evaluated comparing the number of populations and the strength of the signal detected, and also the similarity of the genotype assignments to populations between analyses. In order to test if those results were influenced by the criterion used to select the SSRs, several choosing scenarios based on the discrimination power or the fixation index values of the SSRs were tested. Our results indicate that population structure could be inferred accurately once a certain SSR number threshold was reached, which depended on the underlying structure within the genotypes, but the method used to select the markers included on each set appeared not to be very relevant. The minimum number of SSRs required to provide robust structure inferences and adequate measurements of the differentiation, even when low differentiation levels exist within populations, was proved similar to that of the complete list of recommended markers for fingerprinting. When a SSR set size similar to the minimum marker sets recommended for fingerprinting it is used, only major divisions or moderate (FST>0.05) differentiation of the germplasm are detected.

Subgenomic Diversity Patterns Caused by Directional Selection in Bread Wheat Gene Pools

Genetic diversity represents the fundamental key to breeding success, providing the basis for breeders to select varieties with constantly improving yield performance. On the other hand, strong selection during domestication and breeding have eliminated considerable genetic diversity in the breeding pools of major crops, causing erosion of genetic potential for adaptation to emerging challenges like climate change. High-throughput genomic technologies can address this dilemma by providing detailed knowledge to characterize and replenish genetic diversity in breeding programs. In hexaploid bread wheat (Triticum aestivum L.), the staple food for 35% of the world's population, bottlenecks during allopolyploidisation followed by strong artificial selection have considerably narrowed diversity to the extent that yields in many regions appear to be unexpectedly stagnating. In this study, we used a 90,000 single nucleotide polymorphism (SNP) wheat genotyping array to assay high-frequency, polymorphic SNP markers in 460 accessions representing different phenological diversity groups from Asian, Australian, European, and North American bread wheat breeding materials. Detailed analysis of subgroup diversity at the chromosome and subgenome scale revealed highly distinct patterns of conserved linkage disequilibrium between different gene pools. The data enable identification of genome regions in most need of rejuvenation with novel diversity and provide a high-resolution molecular basis for genomic-assisted introgression of new variation into chromosome segments surrounding directionally selected metaloci conferring important adaptation and quality traits.

Genetic analysis of inflorescence and plant height components in sorghum (Panicoidae) and comparative genetics with rice (Oryzoidae)

BACKGROUND

Domestication has played an important role in shaping characteristics of the inflorescence and plant height in cultivated cereals. Taking advantage of meta-analysis of QTLs, phylogenetic analyses in 502 diverse sorghum accessions, GWAS in a sorghum association panel (n = 354) and comparative data, we provide insight into the genetic basis of the domestication traits in sorghum and rice.

RESULTS

We performed genome-wide association studies (GWAS) on 6 traits related to inflorescence morphology and 6 traits related to plant height in sorghum, comparing the genomic regions implicated in these traits by GWAS and QTL mapping, respectively. In a search for signatures of selection, we identify genomic regions that may contribute to sorghum domestication regarding plant height, flowering time and pericarp color. Comparative studies across taxa show functionally conserved 'hotspots' in sorghum and rice for awn presence and pericarp color that do not appear to reflect corresponding single genes but may indicate co-regulated clusters of genes. We also reveal homoeologous regions retaining similar functions for plant height and flowering time since genome duplication an estimated 70 million years ago or more in a common ancestor of cereals. In most such homoeologous QTL pairs, only one QTL interval exhibits strong selection signals in modern sorghum.

CONCLUSIONS

Intersections among QTL, GWAS and comparative data advance knowledge of genetic determinants of inflorescence and plant height components in sorghum, and add new dimensions to comparisons between sorghum and rice.

PAV markers in Sorghum bicolour: genome pattern, affected genes and pathways, and genetic linkage map construction

KEY MESSAGE

5,511 genic small-size PAVs in sorghum were identified and examined, including the pattern and the function enrichment of PAV genes. 325 PAV markers were developed to construct a genetic map. Presence/absence variants (PAVs) correlate closely to the phenotypic variation, by impacting plant genome sizes and the adaption to the environment. To shed more light on their genome-wide patterns, functions and the possibility of using them as molecular markers, we generated next generation genome sequencing data for four sorghum inbred lines and used associated bioinformatic pipelines to identify small-size PAVs (40-10 kb). Five thousand five hundreds and eleven genic PAVs (40-10 kb) were identified and found to affect 3,238 genes. These PAVs were mainly distributed on the sub-telomeric regions, but the highest proportions occurred in the vicinity of the centromeric regions. One of the prominent features of the PAVs is the high occurrence of long terminal repeats retrotransposons and DNA transposons. PAVs caused various alterations to gene structure, primarily including the coding sequence variants, intron variants, transcript ablation, and initiator codon changes. The genes affected by PAVs were significantly enriched in those involved in stress responses and protein modification. We used 325 PAVs polymorphic between two sorghum inbred lines Ji2731 and E-Tian, together with 49 SSR markers, and constructed a genetic map, which consisted of 10 linkage groups corresponding to the 10 chromosomes of sorghum and spanned 1,430.3 cM in length covering 97% of the physical genome. The resources reported here should be useful for genetic study and breeding of sorghum and related species.

Genetic variation of temperature-regulated curd induction in cauliflower: elucidation of floral transition by genome-wide association mapping and gene expression analysis

Cauliflower (Brassica oleracea var. botrytis) is a vernalization-responsive crop. High ambient temperatures delay harvest time. The elucidation of the genetic regulation of floral transition is highly interesting for a precise harvest scheduling and to ensure stable market supply. This study aims at genetic dissection of temperature-dependent curd induction in cauliflower by genome-wide association studies and gene expression analysis. To assess temperature-dependent curd induction, two greenhouse trials under distinct temperature regimes were conducted on a diversity panel consisting of 111 cauliflower commercial parent lines, genotyped with 14,385 SNPs. Broad phenotypic variation and high heritability (0.93) were observed for temperature-related curd induction within the cauliflower population. GWA mapping identified a total of 18 QTL localized on chromosomes O1, O2, O3, O4, O6, O8, and O9 for curding time under two distinct temperature regimes. Among those, several QTL are localized within regions of promising candidate flowering genes. Inferring population structure and genetic relatedness among the diversity set assigned three main genetic clusters. Linkage disequilibrium (LD) patterns estimated global LD extent of r(2) = 0.06 and a maximum physical distance of 400 kb for genetic linkage. Transcriptional profiling of flowering genes FLOWERING LOCUS C (BoFLC) and VERNALIZATION 2 (BoVRN2) was performed, showing increased expression levels of BoVRN2 in genotypes with faster curding. However, functional relevance of BoVRN2 and BoFLC2 could not consistently be supported, which probably suggests to act facultative and/or might evidence for BoVRN2/BoFLC-independent mechanisms in temperature-regulated floral transition in cauliflower. Genetic insights in temperature-regulated curd induction can underpin genetically informed phenology models and benefit molecular breeding strategies toward the development of thermo-tolerant cultivars.

Sub-genomic selection patterns as a signature of breeding in the allopolyploid Brassica napus genome

Background

High-density single-nucleotide polymorphism (SNP) genotyping arrays are a powerful tool for genome-wide association studies and can give valuable insight into patterns of population structure and linkage disequilibrium (LD). In this study we used the Brassica 60kSNP Illumina consortium genotyping array to assess the influence of selection and breeding for important agronomic traits on LD and haplotype structure in a diverse panel of 203 Chinese semi-winter rapeseed (Brassica napus) breeding lines.

Results

Population structure and principal coordinate analysis, using a subset of the SNPs, revealed diversification into three subpopulations and one mixed population, reflecting targeted introgressions from external gene pools during breeding. Pairwise LD analysis within the A- and C-subgenomes of allopolyploid B. napus revealed that mean LD, at a threshold of r² = 0.1, decayed on average around ten times more rapidly in the A-subgenome (0.25-0.30 Mb) than in the C-subgenome (2.00-2.50 Mb). A total of 3,097 conserved haplotype blocks were detected over a total length of 182.49 Mb (15.17% of the genome). The mean size of haplotype blocks was considerably longer in the C-subgenome (102.85 Kb) than in the A-subgenome (33.51 Kb), and extremely large conserved haplotype blocks were found on a number of C-genome chromosomes. Comparative sequence analysis revealed conserved blocks containing homoloeogous quantitative trait loci (QTL) for seed erucic acid and glucosinolate content, two key seed quality traits under strong agronomic selection. Interestingly, C-subgenome QTL were associated with considerably greater conservation of LD than their corresponding A-subgenome homoeologues.

Conclusions

The data we present in this paper provide evidence for strong selection of large chromosome regions associated with important rapeseed seed quality traits conferred by C-subgenome QTL. This implies that an increase in genetic diversity and recombination within the C-genome is particularly important for breeding. The resolution of genome-wide association studies is also expected to vary greatly across different genome regions.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1170) contains supplementary material, which is available to authorized users.

Population genetics and origin of the native North American agricultural weed waterhemp (Amaranthus tuberculatus; Amaranthaceae)

•

PREMISE OF THE STUDY

The evolution of invasiveness has been extensively studied in natural ecosystems; however, far less is known about the evolution of agricultural invasiveness, despite the major economic impact of weeds on crop productivity. Examining the population structure of recently arisen weeds can provide insights into evolutionary avenues to invasion of agroecosystems. Weeds that originate from wild plants are the most common yet least frequently studied type of agricultural invasive. Here we address several questions about the origin of the native North American agricultural weed waterhemp (Amaranthus tuberculatus), which invaded corn and soy fields in the midwestern United States in the 20th century.•

METHODS

We genotyped 38 populations from across the species range with 10 microsatellite markers and used these data to assess genetic diversity and population structure within and outside the geographical region where waterhemp is agriculturally problematic.•

KEY RESULTS

We found evidence for two ancestral genetic lineages in our data, supporting the hypothesis that A. tuberculatus was diverging into two evolutionary lineages prior to the 20th century. However, we found no support for the hypothesis that agricultural weed populations arose from admixture of these two lineages after secondary contact. Our data suggest that eastward movement of the western genetic lineage, facilitated by changing agricultural practices, is the source of the agricultural invasion of waterhemp.•

CONCLUSIONS

This research demonstrates that agricultural invasion by native, wild plant species can proceed via different evolutionary trajectories from weeds related to domesticated plants, which has implications for evolutionary biology and weed control.

Duplicate and conquer: multiple homologs of PHOSPHORUS-STARVATION TOLERANCE1 enhance phosphorus acquisition and sorghum performance on low-phosphorus soils

Low soil phosphorus (P) availability is a major constraint for crop production in tropical regions. The rice (Oryza sativa) protein kinase, PHOSPHORUS-STARVATION TOLERANCE1 (OsPSTOL1), was previously shown to enhance P acquisition and grain yield in rice under P deficiency. We investigated the role of homologs of OsPSTOL1 in sorghum (Sorghum bicolor) performance under low P. Association mapping was undertaken in two sorghum association panels phenotyped for P uptake, root system morphology and architecture in hydroponics and grain yield and biomass accumulation under low-P conditions, in Brazil and/or in Mali. Root length and root surface area were positively correlated with grain yield under low P in the soil, emphasizing the importance of P acquisition efficiency in sorghum adaptation to low-P availability. SbPSTOL1 alleles reducing root diameter were associated with enhanced P uptake under low P in hydroponics, whereas Sb03g006765 and Sb03g0031680 alleles increasing root surface area also increased grain yield in a low-P soil. SbPSTOL1 genes colocalized with quantitative trait loci for traits underlying root morphology and dry weight accumulation under low P via linkage mapping. Consistent allelic effects for enhanced sorghum performance under low P between association panels, including enhanced grain yield under low P in the soil in Brazil, point toward a relatively stable role for Sb03g006765 across genetic backgrounds and environmental conditions. This study indicates that multiple SbPSTOL1 genes have a more general role in the root system, not only enhancing root morphology traits but also changing root system architecture, which leads to grain yield gain under low-P availability in the soil.

Genetic dissection of temperature-dependent sorghum growth during juvenile development

Promising genome regions for improving cold tolerance of sorghum were identified on chromosomes SBI-01, SBI-03, SBI-07, and SBI-10. Chlorophyll fluorescence had no major effect on growth rates at low temperatures. Developing fast growing sorghum seedlings is an important breeding goal for temperate climates since low springtime temperatures are resulting in a prolonged juvenile development. The adaptation of sorghum to tropical and subtropical highlands gives hint for certain genetic variation. The goals of the present study were to detect marker-trait associations for leaf and dry matter growth rate and for chlorophyll fluorescence and content (SPAD) in relation to temperature. A diversity set comprising 194 genotypes was tested in eight controlled environments with temperatures ranging from 9.4 to 20.8 °C. Significant marker-trait associations (p < 0.05) were identified for each individual temperature regime and on the parameters of regression analyses describing the responses of growth or chlorophyll related traits to temperatures. The diversity set was fingerprinted with 171 diversity array technology (DArT) and 31 simple-sequence repeat (SSR) markers. SSRs were used to analyze the population structure while association studies were performed on DArT markers. Promising marker-trait associations for growth rates in relation to temperature were detected on chromosomes SBI-01, SBI-03, SBI-07, and SBI-10. Many promising loci were also significantly associated to the results obtained in individual low-temperature environments. Marker-trait associations for chlorophyll content and fluorescence did occasionally co-locate to those for growth during juvenile development but there was no evidence supporting our hypothesis that seedling growth at low temperatures is largely influenced by SPAD or fluorescence.

Two in one sweep: aluminum tolerance and grain yield in P-limited soils are associated to the same genomic region in West African sorghum

BACKGROUND

Sorghum (Sorghum bicolor L. Moench) productivity is severely impeded by low phosphorus (P) and aluminum (Al) toxic soils in sub-Saharan Africa and especially West Africa (WA). Improving productivity of this staple crop under these harsh conditions is crucial to improve food security and farmer's incomes in WA.

RESULTS

This is the first study to examine the genetics underlying sorghum adaptation to phosphorus limitation in a wide range of WA growing conditions. A set of 187 diverse sorghum genotypes were grown in 29 -P and + P field experiments from 2006-2012 in three WA countries. Sorghum grain yield performance under -P and + P conditions was highly correlated (r = 0.85***). Significant genotype-by-phosphorus interaction was detected but with small magnitude compared to the genotype variance component. We observed high genetic diversity within our panel, with rapid linkage disequilibrium decay, confirming recent sequence based studies in sorghum. Using genome wide association mapping based on 220 934 SNPs we identified one genomic region on chromosome 3 that was highly associated to grain yield production. A major Al-tolerance gene in sorghum, SbMATE, was collocated in this region and SbMATE specific SNPs showed very high associations to grain yield production, especially under -P conditions, explaining up to 16% of the genotypic variance.

CONCLUSION

The results suggest that SbMATE has a possible pleiotropic role in providing tolerance to two of the most serious abiotic stresses for sorghum in WA, Al toxicity and P deficiency. The identified SNPs can help accelerate breeding for increased sorghum productivity under unfavorable soil conditions and contribute to assuring food security in WA.

Whole-genome sequencing reveals untapped genetic potential in Africa's indigenous cereal crop sorghum

Sorghum is a food and feed cereal crop adapted to heat and drought and a staple for 500 million of the world’s poorest people. Its small diploid genome and phenotypic diversity make it an ideal C₄ grass model as a complement to C₃ rice. Here we present high coverage (16–45 × ) resequenced genomes of 44 sorghum lines representing the primary gene pool and spanning dimensions of geographic origin, end-use and taxonomic group. We also report the first resequenced genome of S. propinquum, identifying 8 M high-quality SNPs, 1.9 M indels and specific gene loss and gain events in S. bicolor. We observe strong racial structure and a complex domestication history involving at least two distinct domestication events. These assembled genomes enable the leveraging of existing cereal functional genomics data against the novel diversity available in sorghum, providing an unmatched resource for the genetic improvement of sorghum and other grass species.

Sorghum is a drought-resistant food and feed cereal crop used by over half a billion of the world’s poorest people. Here the authors present high-coverage resequencing genome data of 44 sorghum lines of varying geographic and taxonomic origin, which include a number of sorghum wild relatives.

Association mapping provides insights into the origin and the fine structure of the sorghum aluminum tolerance locus, AltSB

Root damage caused by aluminum (Al) toxicity is a major cause of grain yield reduction on acid soils, which are prevalent in tropical and subtropical regions of the world where food security is most tenuous. In sorghum, Al tolerance is conferred by SbMATE, an Al-activated root citrate efflux transporter that underlies the major Al tolerance locus, AltSB, on sorghum chromosome 3. We used association mapping to gain insights into the origin and evolution of Al tolerance in sorghum and to detect functional variants amenable to allele mining applications. Linkage disequilibrium across the AltSB locus decreased much faster than in previous reports in sorghum, and reached basal levels at approximately 1000 bp. Accordingly, intra-locus recombination events were found to be extensive. SNPs and indels highly associated with Al tolerance showed a narrow frequency range, between 0.06 and 0.1, suggesting a rather recent origin of Al tolerance mutations within AltSB. A haplotype network analysis suggested a single geographic and racial origin of causative mutations in primordial guinea domesticates in West Africa. Al tolerance assessment in accessions harboring recombinant haplotypes suggests that causative polymorphisms are localized to a ∼6 kb region including intronic polymorphisms and a transposon (MITE) insertion, whose size variation has been shown to be positively correlated with Al tolerance. The SNP with the strongest association signal, located in the second SbMATE intron, recovers 9 of the 14 highly Al tolerant accessions and 80% of all the Al tolerant and intermediately tolerant accessions in the association panel. Our results also demonstrate the pivotal importance of knowledge on the origin and evolution of Al tolerance mutations in molecular breeding applications. Allele mining strategies based on associated loci are expected to lead to the efficient identification, in diverse sorghum germplasm, of Al tolerant accessions able maintain grain yields under Al toxicity.

High-throughput genomics in sorghum: from whole-genome resequencing to a SNP screening array

With its small, diploid and completely sequenced genome, sorghum (Sorghum bicolor L. Moench) is highly amenable to genomics-based breeding approaches. Here, we describe the development and testing of a robust single-nucleotide polymorphism (SNP) array platform that enables polymorphism screening for genome-wide and trait-linked polymorphisms in genetically diverse S. bicolor populations. Whole-genome sequences with 6× to 12× coverage from five genetically diverse S. bicolor genotypes, including three sweet sorghums and two grain sorghums, were aligned to the sorghum reference genome. From over 1 million high-quality SNPs, we selected 2124 Infinium Type II SNPs that were informative in all six source genomes, gave an optimal Assay Design Tool (ADT) score, had allele frequencies of 50% in the six genotypes and were evenly spaced throughout the S. bicolor genome. Furthermore, by phenotype-based pool sequencing, we selected an additional 876 SNPs with a phenotypic association to early-stage chilling tolerance, a key trait for European sorghum breeding. The 3000 attempted bead types were used to populate half of a dual-species Illumina iSelect SNP array. The array was tested using 564 Sorghum spp. genotypes, including offspring from four unrelated recombinant inbred line (RIL) and F2 populations and a genetic diversity collection. A high call rate of over 80% enabled validation of 2620 robust and polymorphic sorghum SNPs, underlining the efficiency of the array development scheme for whole-genome SNP selection and screening, with diverse applications including genetic mapping, genome-wide association studies and genomic selection.

SNP genotyping reveals genetic diversity between cultivated landraces and contemporary varieties of tomato

Background

The tomato (Solanum lycopersium L.) is the most widely grown vegetable in the world. It was domesticated in Latin America and Italy and Spain are considered secondary centers of diversification. This food crop has experienced severe genetic bottlenecks and modern breeding activities have been characterized by trait introgression from wild species and divergence in different market classes.

Results

With the aim to examine patterns of polymorphism, characterize population structure and identify putative loci under positive selection, we genotyped 214 tomato accessions (which include cultivated landraces, commercial varieties and wild relatives) using a custom-made Illumina SNP-panel. Most of the 175 successfully scored SNP loci were found to be polymorphic. Population structure analysis and estimates of genetic differentiation indicated that landraces constitute distinct sub-populations. Furthermore, contemporary varieties could be separated in groups (processing, fresh and cherry) that are consistent with the recent breeding aimed at market-class specialization. In addition, at the 95% confidence level, we identified 30, 34 and 37 loci under positive selection between landraces and each of the groups of commercial variety (cherry, processing and fresh market, respectively). Their number and genomic locations imply the presence of some extended regions with high genetic variation between landraces and contemporary varieties.

Conclusions

Our work provides knowledge concerning the level and distribution of genetic variation within cultivated tomato landraces and increases our understanding of the genetic subdivision of contemporary varieties. The data indicate that adaptation and selection have led to a genomic signature in cultivated landraces and that the subpopulation structure of contemporary varieties is shaped by directed breeding and largely of recent origin. The genomic characterization presented here is an essential step towards a future exploitation of the available tomato genetic resources in research and breeding programs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-835) contains supplementary material, which is available to authorized users.

Assessment of genetic diversity in the sorghum reference set using EST-SSR markers

Selection and use of genetically diverse genotypes are key factors in any crop breeding program to develop cultivars with a broad genetic base. Molecular markers play a major role in selecting diverse genotypes. In the present study, a reference set representing a wide range of sorghum genetic diversity was screened with 40 EST-SSR markers to validate both the use of these markers for genetic structure analyses and the population structure of this set. Grouping of accessions is identical in distance-based and model-based clustering methods. Genotypes were grouped primarily based on race within the geographic origins. Accessions derived from the African continent contributed 88.6 % of alleles confirming the African origin of sorghum. In total, 360 alleles were detected in the reference set with an average of 9 alleles per marker. The average PIC value was 0.5230 with a range of 0.1379-0.9483. Sub-race, guinea margaritiferum (Gma) from West Africa formed a separate cluster in close proximity to wild accessions suggesting that the Gma group represents an independent domestication event. Guineas from India and Western Africa formed two distinct clusters. Accessions belongs to the kafir race formed the most homogeneous group as observed in earlier studies. This analysis suggests that the EST-SSR markers used in the present study have greater discriminating power than the genomic SSRs. Genetic variance within the subpopulations was very high (71.7 %) suggesting that the germplasm lines included in the set are more diverse. Thus, this reference set representing the global germplasm is an ideal material for the breeding community, serving as a community resource for trait-specific allele mining as well as genome-wide association mapping.

Identification of genetic markers linked to anthracnose resistance in sorghum using association analysis

Anthracnose in sorghum caused by Colletotrichum sublineolum is one of the most destructive diseases affecting sorghum production under warm and humid conditions. Markers and genes linked to resistance to the disease are important for plant breeding. Using 14,739 SNP markers, we have mapped eight loci linked to resistance in sorghum through association analysis of a sorghum mini-core collection consisting of 242 diverse accessions evaluated for anthracnose resistance for 2 years in the field. The mini-core was representative of the International Crops Research Institute for the Semi-Arid Tropics' world-wide sorghum landrace collection. Eight marker loci were associated with anthracnose resistance in both years. Except locus 8, disease resistance-related genes were found in all loci based on their physical distance from linked SNP markers. These include two NB-ARC class of R genes on chromosome 10 that were partially homologous to the rice blast resistance gene Pib, two hypersensitive response-related genes: autophagy-related protein 3 on chromosome 1 and 4 harpin-induced 1 (Hin1) homologs on chromosome 8, a RAV transcription factor that is also part of R gene pathway, an oxysterol-binding protein that functions in the non-specific host resistance, and homologs of menthone:neomenthol reductase (MNR) that catalyzes a menthone reduction to produce the antimicrobial neomenthol. These genes and markers may be developed into molecular tools for genetic improvement of anthracnose resistance in sorghum.

Polymorphism of TaSAP1-A1 and its association with agronomic traits in wheat

TaSAP1, a member of the stress association protein (SAP) gene family from wheat, is involved in response to several abiotic stresses, including drought, salt and cold. In this study, TaSAP1-A1, a TaSAP1 member on chromosome 7A, and its flanking sequences were isolated. Polymorphism analysis indicated that the average nucleotide diversity (π) of the whole region was 0.00296. The highest nucleotide diversity occurred in the promoter region (π = 0.00631) and no polymorphism was identified in the coding region. Three markers T7AM5, T7AM2606 and T7AM39 located in the promoter region, were developed from sequence variations (InDel5-1810, SNP-2606 and InDel39-1637). Six haplotypes were identified among 300 accessions based on the three markers. TaSAP1-A1 was located on chromosome 7A using marker T7AM39 and was flanked by markers Xwmc530 and Xbarc174. QTL for yield-related traits, including 1,000-grain weight, number of grains per spike and grain yield, were located in the same region. In marker- and haplotype-trait association analyses, TaSAP1-A1 was significantly associated with 1,000-grain weight, number of grains per spike, spike length, peduncle length and total number of spikelets per spike in multiple environments. These results provide useful information for marker-assisted selection for yield-related traits under well-watered and drought-stressed conditions.

Genetic structure and linkage disequilibrium in a diverse, representative collection of the C4 model plant, Sorghum bicolor

To facilitate the mapping of genes in sorghum [Sorghum bicolor (L.) Moench] underlying economically important traits, we analyzed the genetic structure and linkage disequilibrium in a sorghum mini core collection of 242 landraces with 13,390 single-nucleotide polymorphims. The single-nucleotide polymorphisms were produced using a highly multiplexed genotyping-by-sequencing methodology. Genetic structure was established using principal component, Neighbor-Joining phylogenetic, and Bayesian cluster analyses. These analyses indicated that the mini-core collection was structured along both geographic origin and sorghum race classification. Examples of the former were accessions from Southern Africa, East Asia, and Yemen. Examples of the latter were caudatums with widespread geographical distribution, durras from India, and guineas from West Africa. Race bicolor, the most primitive and the least clearly defined sorghum race, clustered among other races and formed only one clear bicolor-centric cluster. Genome-wide linkage disequilibrium analyses showed linkage disequilibrium decayed, on average, within 10-30 kb, whereas the short arm of SBI-06 contained a linkage disequilibrium block of 20.33 Mb, confirming a previous report of low recombination on this chromosome arm. Four smaller but equally significant linkage disequilibrium blocks of 3.5-35.5 kb were detected on chromosomes 1, 2, 9, and 10. We examined the genes encoded within each block to provide a first look at candidates such as homologs of GS3 and FT that may indicate a selective sweep during sorghum domestication.

Massive sorghum collection genotyped with SSR markers to enhance use of global genetic resources

Large ex situ collections require approaches for sampling manageable amounts of germplasm for in-depth characterization and use. We present here a large diversity survey in sorghum with 3367 accessions and 41 reference nuclear SSR markers. Of 19 alleles on average per locus, the largest numbers of alleles were concentrated in central and eastern Africa. Cultivated sorghum appeared structured according to geographic regions and race within region. A total of 13 groups of variable size were distinguished. The peripheral groups in western Africa, southern Africa and eastern Asia were the most homogeneous and clearly differentiated. Except for Kafir, there was little correspondence between races and marker-based groups. Bicolor, Caudatum, Durra and Guinea types were each dispersed in three groups or more. Races should therefore better be referred to as morphotypes. Wild and weedy accessions were very diverse and scattered among cultivated samples, reinforcing the idea that large gene-flow exists between the different compartments. Our study provides an entry to global sorghum germplasm collections. Our reference marker kit can serve to aggregate additional studies and enhance international collaboration. We propose a core reference set in order to facilitate integrated phenotyping experiments towards refined functional understanding of sorghum diversity.

Population genomic and genome-wide association studies of agroclimatic traits in sorghum

Accelerating crop improvement in sorghum, a staple food for people in semiarid regions across the developing world, is key to ensuring global food security in the context of climate change. To facilitate gene discovery and molecular breeding in sorghum, we have characterized ~265,000 single nucleotide polymorphisms (SNPs) in 971 worldwide accessions that have adapted to diverse agroclimatic conditions. Using this genome-wide SNP map, we have characterized population structure with respect to geographic origin and morphological type and identified patterns of ancient crop diffusion to diverse agroclimatic regions across Africa and Asia. To better understand the genomic patterns of diversification in sorghum, we quantified variation in nucleotide diversity, linkage disequilibrium, and recombination rates across the genome. Analyzing nucleotide diversity in landraces, we find evidence of selective sweeps around starch metabolism genes, whereas in landrace-derived introgression lines, we find introgressions around known height and maturity loci. To identify additional loci underlying variation in major agroclimatic traits, we performed genome-wide association studies (GWAS) on plant height components and inflorescence architecture. GWAS maps several classical loci for plant height, candidate genes for inflorescence architecture. Finally, we trace the independent spread of multiple haplotypes carrying alleles for short stature or long inflorescence branches. This genome-wide map of SNP variation in sorghum provides a basis for crop improvement through marker-assisted breeding and genomic selection.