Whole-genome diversity, population structure and linkage disequilibrium analysis of globally diverse wheat genotypes using genotyping-by-sequencing DArTseq platform

Extent and patterns of morphological and molecular genetic diversity and population structure of Nigerian Taro cultivars

BACKGROUND

Genetic diversity is crucial for conservation efforts as well as breeding programs targeted at the development of improved varieties. Taro, a climate-resilient crop, plays a vital role in the nutritional and economic livelihoods of many households in Nigeria, but its yield is very low due to inadequate genetic improvement efforts. A diversity assessment of Nigerian taro is therefore required to create a premise for its improvement in yield, quality and disease tolerance. In this study, the genetic diversity and population structure of 490 taro cultivars comprising two main gene pools: Dasheen (215) and Eddoe (275), collected from farmers and marketers across seven states in Nigeria was assessed using 3047 Diversity Array Technology single nucleotide polymorphism (DArT-SNP) markers. A subset of 114 taro cultivars, comprising 30 Dasheens and 84 Eddoes were further phenotyped using 24 agro-morphological descriptors.

RESULTS

Both phenotypic and molecular characterization revealed higher genetic diversity among the Eddoes than Dasheens. Estimates of gene flow (Nm = 0.353) revealed intermixing of cultivars among the States of collection, with the highest gene flow occurring between cultivars from Anambra and Ondo states and the lowest between Anambra and Kwara states. Population structure and Ward's minimum variance hierarchical cluster based on DArT-SNPs identified four groups, one comprising Dasheen and three comprising Eddoe cultivars. Hierarchical clustering based on phenotypic traits delineated three clusters. Variation between gene pools (49%) was higher than within gene pools (32%). Variation among States of collection was high (41%), while variation among individuals within gene pools (18%) and States of collection (19%) was relatively low. Correlation between phenotypic and genotypic diversity assessments was low (r = 0.01), indicating that both approaches were necessary for assessing genetic diversity in taro. However, genotypic assessment provided better information about genetic diversity of the taro cultivars.

CONCLUSION

This is the first study that represented germplasm collection across the major taro growing regions of Nigeria. The findings from this study based on agro-morphological characterization and DArT-SNP genotyping are critical for genetic characterization, conservation and breeding of taro in Nigeria, mainly initiating hybridization between the two genepools after careful assessment of ploidy levels of the accessions collected in this study. This will facilitate in developing improved taro varieties with desirable traits, such as higher yield, better disease resistance, and improved nutritional quality.

Harnessing genome-wide genetic diversity, population structure and linkage disequilibrium in Ethiopian durum wheat gene pool

Yanyang Liu, Henan Academy of Agricultural Sciences (HNAAS), China; Landraces are an important genetic source for transferring valuable novel genes and alleles required to enhance genetic variation. Therefore, information on the gene pool's genetic diversity and population structure is essential for the conservation and sustainable use of durum wheat genetic resources. Hence, the aim of this study was to assess genetic diversity, population structure, and linkage disequilibrium, as well as to identify regions with selection signature. Five hundred (500) individuals representing 46 landraces, along with 28 cultivars were evaluated using the Illumina Infinium 25K wheat SNP array, resulting in 8,178 SNPs for further analysis. Gene diversity (GD) and the polymorphic information content (PIC) ranged from 0.13-0.50 and 0.12-0.38, with mean GD and PIC values of 0.34 and 0.27, respectively. Linkage disequilibrium (LD) revealed 353,600 pairs of significant SNPs at a cut-off (r2 > 0.20, P < 0.01), with an average r2 of 0.21 for marker pairs. The nucleotide diversity (π) and Tajima's D (TD) per chromosome for the populations ranged from 0.29-0.36 and 3.46-5.06, respectively, with genome level, mean π values of 0.33 and TD values of 4.43. Genomic scan using the Fst outlier test revealed 85 loci under selection signatures, with 65 loci under balancing selection and 17 under directional selection. Putative candidate genes co-localized with regions exhibiting strong selection signatures were associated with grain yield, plant height, host plant resistance to pathogens, heading date, grain quality, and phenolic content. The Bayesian Model (STRUCTURE) and distance-based (principal coordinate analysis, PCoA, and unweighted pair group method with arithmetic mean, UPGMA) methods grouped the genotypes into five subpopulations, where landraces from geographically non-adjoining environments were clustered in the same cluster. This research provides further insights into population structure and genetic relationships in a diverse set of durum wheat germplasm, which could be further used in wheat breeding programs to address production challenges sustainably.

SNP markers identification by genome wide association study for chemical quality traits of coffee (Coffea spp.) Germplasm

BACKGROUND

Coffee quality is an important selection criterion for coffee breeding. Metabolite profiling and Genome-Wide Association Studies (GWAS) effectively dissect the genetic background of complex traits such as metabolites content (caffeine, trigonelline, and 5-caffeoylquinic acid (5-CQA)) in coffee that affect quality. Therefore, it is important to determine the metabolic profiles of Coffea spp. genotypes. This study aimed to identify Single Nucleotide Polymorphisms (SNPs) within Coffea spp. genotypes through GWAS and associate these significant SNPs to the metabolic profiles of the different genotypes.

METHODS AND RESULTS

A total of 1,739 SNP markers were obtained from 80 genotypes using the DArTseq™ method. Caffeine, trigonelline, and 5-CQA content were determined in coffee leaves using Ultra-Performance Liquid Chromatography/tandem mass spectrometry (UPLC-MS/MS) analyses. The GWAS was carried out using the Genome Association and Prediction Integrated Tool (GAPIT) software and a compressed mixed linear model. Finally, a total of three significant SNP markers out of ten were identified. One SNP, located in the coffee chromosome (Chr) 8, was significantly associated with caffeine. The two remaining SNPs, located in Chr 4 and 5, were significantly associated with trigonelline and six SNPs markers were associated with 5-CQA in Chr 1, 5 and 10, but these six markers were not significant.

CONCLUSIONS

These significant SNP sequences were associated with protein ubiquitination, assimilation, and wall receptor kinases. Therefore, these SNPs might be useful hits in subsequent quality coffee breeding programs.

Exploration of the Genetic Diversity of Solina Wheat and Its Implication for Grain Quality

Different Solina wheat accessions (n = 24) collected in the Abruzzo region (Italy) were studied using 45,000 SNP markers generated from the DarTseq platform. The structure of genetic data was analyzed by Principal Component Analysis and Hierarchical Cluster analysis that revealed the existence of two main clusters (Clu1 and Clu2) characterized by samples with different geographical origin. The Solina genetic dataset was further merged and analyzed with a public genetic one provided by CIMMYT containing 25,963 genotypes from all over the world. The Solina accessions occupied a vast space, thus confirming a high heterogeneity of this landrace that, nevertheless, is considerably unique and placed quite far from other clusters. Clu1 and Clu2 divergence were clearly visible. Solina clusters were genetically closer to landraces from Turkey and the central fertile crescent than to the Italian genotypes present in the dataset. Selected commercial quality traits of accessions of the two Solina clusters were analyzed (yield, thousand kernel weight, test weight, and protein content), and significant differences were found between clusters. The results of this investigation did not highlight any relationships of Solina with Italian genotypes, and confirmed its wide genetic diversity by permitting to identify two genetic groups with distinct origin and quality traits.

Assessment of genetic relationships among native and introduced Himalayan balsam (Impatiens glandulifera) plants based on genome profiling

We conducted genomic characterization based on SNP and SilicoDArT markers on the invasive Himalayan balsam (Impatiens glandulifera) plants originating from native and non-native regions of their distribution. When genetic relationships were explored by PCoA using SNP and SilicoDArT marker data, the first, second, and third principal coordinates explained altogether 37.4% and 31.0% of the variability, respectively. Samples from the UK, Canada, and Pakistan were grouped together, while Indian plants were clearly distinct based on SNP markers but relatively close to the UK-Canada-Pakistan group based on SilicoDArT markers. Constructed trees differentiated individuals into clusters resembling the PCoA patterns. The Bayesian BAPS analysis performed for the SNP data revealed that the individuals were distributed in seven clusters, representing samples from each of the four Finnish populations, India, Pakistan, and the combination of the UK and Canada. Similar clustering was visible in the UPGMA tree. The Indian cluster did not display any ancestral gene flow with the others, while the Pakistani cluster showed ancestral gene flow only with the combined UK and Canada cluster. Furthermore, the latter cluster displayed ancestral gene flow with the Finnish populations varying from 0% to 3.1%. The BAPS analyses conducted for the SilicoDArT data differ slightly: The individuals were distributed in nine clusters, and the Indian cluster exhibited ancestral gene flow with the mixed cluster including Canadian, Pakistani, and UK samples, and one Finnish sample. The AMOVA showed that 45% and 26% of variation was present among the I. glandulifera groups/populations and the rest within them based on SNP and SilicoDArT markers, respectively. The Bayesian BAPS analyses and the gene flow networks were the most informative tools for resolving relationships among native and introduced plants. It is notable that the small sample sizes for non-Finnish plant materials may affect the accuracy of the gene flow and other estimates.