Methods for linkage disequilibrium mapping in crops

Omics-driven utilization of wild relatives for empowering pre-breeding in pearl millet

MAIN CONCLUSION

Pearl millet wild relatives harbour novel alleles which could be utilized to broaden genetic base of cultivated species. Genomics-informed pre-breeding is needed to speed up introgression from wild to cultivated gene pool in pearl millet. Rising episodes of intense biotic and abiotic stresses challenge pearl millet production globally. Wild relatives provide a wide spectrum of novel alleles which could address challenges posed by climate change. Pre-breeding holds potential to introgress novel diversity in genetically narrow cultivated Pennisetum glaucum from diverse gene pool. Practical utilization of gene pool diversity remained elusive due to genetic intricacies. Harnessing promising traits from wild pennisetum is limited by lack of information on underlying candidate genes/QTLs. Next-Generation Omics provide vast scope to speed up pre-breeding in pearl millet. Genomic resources generated out of draft genome sequence and improved genome assemblies can be employed to utilize gene bank accessions effectively. The article highlights genetic richness in pearl millet and its utilization with a focus on harnessing next-generation Omics to empower pre-breeding.

Reviewing the essential roles of remote phenotyping, GWAS and explainable AI in practical marker-assisted selection for drought-tolerant winter wheat breeding

Marker-assisted selection (MAS) plays a crucial role in crop breeding improving the speed and precision of conventional breeding programmes by quickly and reliably identifying and selecting plants with desired traits. However, the efficacy of MAS depends on several prerequisites, with precise phenotyping being a key aspect of any plant breeding programme. Recent advancements in high-throughput remote phenotyping, facilitated by unmanned aerial vehicles coupled to machine learning, offer a non-destructive and efficient alternative to traditional, time-consuming, and labour-intensive methods. Furthermore, MAS relies on knowledge of marker-trait associations, commonly obtained through genome-wide association studies (GWAS), to understand complex traits such as drought tolerance, including yield components and phenology. However, GWAS has limitations that artificial intelligence (AI) has been shown to partially overcome. Additionally, AI and its explainable variants, which ensure transparency and interpretability, are increasingly being used as recognised problem-solving tools throughout the breeding process. Given these rapid technological advancements, this review provides an overview of state-of-the-art methods and processes underlying each MAS, from phenotyping, genotyping and association analyses to the integration of explainable AI along the entire workflow. In this context, we specifically address the challenges and importance of breeding winter wheat for greater drought tolerance with stable yields, as regional droughts during critical developmental stages pose a threat to winter wheat production. Finally, we explore the transition from scientific progress to practical implementation and discuss ways to bridge the gap between cutting-edge developments and breeders, expediting MAS-based winter wheat breeding for drought tolerance.

Genome-wide association study and development of molecular markers for yield and quality traits in peanut (Arachis hypogaea L.)

BACKGROUND

This study aims to decipher the genetic basis governing yield components and quality attributes of peanuts, a critical aspect for advancing molecular breeding techniques. Integrating genotype re-sequencing and phenotypic evaluations of seven yield components and two grain quality traits across four distinct environments allowed for the execution of a genome-wide association study (GWAS).

RESULTS

The nine phenotypic traits were all continuous and followed a normal distribution. The broad heritability ranged from 88.09 to 98.08%, and the genotype-environment interaction effects were all significant. There was a highly significant negative correlation between protein content (PC) and oil content (OC). The 10× genome re-sequencing of 199 peanut accessions yielded a total of 631,988 high-quality single nucleotide polymorphisms (SNPs), with 374 significant SNP loci identified in association with the nine traits of interest. Notably, 66 of these pertinent SNPs were detected in multiple environments, and 48 of them were linked to multiple traits of interest. Five loci situated on chromosome 16 (Chr16) exhibited pleiotropic effects on yield traits, accounting for 17.64-32.61% of the observed phenotypic variation. Two loci on Chr08 were found to be strongly associated with protein and oil contents, accounting for 12.86% and 14.06% of their respective phenotypic variations, respectively. Linkage disequilibrium (LD) block analysis of these seven loci unraveled five nonsynonymous variants, leading to the identification of one yield-related candidate gene and two quality-related candidate genes. The correlation between phenotypic variation and SNP loci in these candidate genes was validated by Kompetitive allele-specific PCR (KASP) marker analysis.

CONCLUSIONS

Overall, molecular markers were developed for genetic loci associated with yield and quality traits through a GWAS investigation of 199 peanut accessions across four distinct environments. These molecular tools can aid in the development of desirable peanut germplasm with an equilibrium of yield and quality through marker-assisted breeding.

Genome-wide association analysis of time to heading and maturity in bread wheat using 55K microarrays

To investigate the genetic mechanisms underlying the reproductive traits (time to flowering and maturity) in wheat and identify candidate genes associated, a phenotypic analysis was conducted on 239 wheat accessions (lines) from around the world. A genome-wide association study (GWAS) of wheat heading and maturity phases was performed using the MLM (Q+K) model in the TASSLE software, combined with the Wheat 55K SNP array. The results revealed significant phenotypic variation in heading and maturity among the wheat accessions across different years, with coefficients of variation ranging from 0.96% to 1.97%. The phenotypic data from different years exhibited excellent correlation, with a genome-wide linkage disequilibrium (LD) attenuation distance of 3 Mb. Population structure analysis, evolutionary tree analysis, and principal component analysis indicated that the 239 wheat accessions formed a relatively homogeneous natural population, which could be divided into three subgroups. The GWAS results identified a total of 293 SNP marker loci that were significantly associated with wheat heading and maturity stages (P ≤ 0.001) in different environments. Among them, nine stable SNP marker loci were consistently detected in multiple environments. These marker loci were distributed on wheat chromosomes 1A、1B、2D、3A、5B、6D and 7A. Each individual locus explained 4.03%-16.06% of the phenotypic variation. Furthermore, through careful analysis of the associated loci with large phenotypic effect values and stable inheritance, a total of nine candidate genes related to wheat heading and maturity stages were identified. These findings have implications for molecular marker-assisted selection breeding programs targeting specific wheat traits at the heading and maturity stages. In summary, this study conducted a comprehensive GWAS of wheat heading and maturity phases, revealing significant associations between genetic markers and key developmental stages in wheat. The identification of candidate genes and marker loci provides valuable information for further studies on wheat breeding and genetic improvement targeted at enhancing heading and maturity traits.

Leveraging Single-Cell Populations to Uncover the Genetic Basis of Complex Traits

The ease and throughput of single-cell genomics have steadily improved, and its current trajectory suggests that surveying single-cell populations will become routine. We discuss the merger of quantitative genetics with single-cell genomics and emphasize how this synergizes with advantages intrinsic to plants. Single-cell population genomics provides increased detection resolution when mapping variants that control molecular traits, including gene expression or chromatin accessibility. Additionally, single-cell population genomics reveals the cell types in which variants act and, when combined with organism-level phenotype measurements, unveils which cellular contexts impact higher-order traits. Emerging technologies, notably multiomics, can facilitate the measurement of both genetic changes and genomic traits in single cells, enabling single-cell genetic experiments. The implementation of single-cell genetics will advance the investigation of the genetic architecture of complex molecular traits and provide new experimental paradigms to study eukaryotic genetics.

RecView: an interactive R application for locating recombination positions using pedigree data

BACKGROUND

Recombination reshuffles alleles at linked loci, allowing genes to evolve independently and consequently enhancing the efficiency of selection. This makes quantifying recombination along chromosomes an important goal for understanding how selection and drift are acting on genes and chromosomes.

RESULTS

We present RecView, an interactive R application and its homonymous R package, to facilitate locating recombination positions along chromosomes or scaffolds using whole-genome genotype data of a three-generation pedigree. RecView analyses and plots the grandparent-of-origin of all informative alleles along each chromosome of the offspring in the pedigree, and infers recombination positions with either of two built-in algorithms: one based on change in the proportion of the alleles with specific grandparent-of-origin, and one on the degree of continuity of alleles with the same grandparent-of-origin. RecView handles multiple offspring and chromosomes simultaneously, and all putative recombination positions are reported in base pairs together with an estimated precision based on the local density of informative alleles. We demonstrate RecView using genotype data of a passerine bird with an available reference genome, the great reed warbler (Acrocephalus arundinaceus), and show that recombination events can be located to specific positions.

CONCLUSIONS

RecView is an easy-to-use and highly effective application for locating recombination positions with high precision. RecView is available on GitHub ( https://github.com/HKyleZhang/RecView.git ).

Recombinant inbred lines derived from wide crosses in Pisum

Genomic resources are becoming available for Pisum but to link these to phenotypic diversity requires well marked populations segregating for relevant traits. Here we describe two such resources. Two recombinant inbred populations, derived from wide crosses in Pisum are described. One high resolution mapping population involves cv Caméor, for which the first pea whole genome assembly was obtained, crossed to JI0281, a basally divergent P. sativum sativum landrace from Ethiopia. The other is an inter sub-specific cross between P. s. sativum and the independently domesticated P. s. abyssinicum. The corresponding genetic maps provide information on chromosome level sequence assemblies and identify structural differences between the genomes of these two Pisum subspecies. In order to visualise chromosomal translocations that distinguish the mapping parents, we created a simplified version of Threadmapper to optimise it for interactive 3-dimensional display of multiple linkage groups. The genetic mapping of traits affecting seed coat roughness and colour, plant height, axil ring pigmentation, leaflet number and leaflet indentation enabled the definition of their corresponding genomic regions. The consequence of structural rearrangement for trait analysis is illustrated by leaf serration. These analyses pave the way for identification of the underlying genes and illustrate the utility of these publicly available resources. Segregating inbred populations derived from wide crosses in Pisum, together with the associated marker data, are made publicly available for trait dissection. Genetic analysis of these populations is informative about chromosome scale assemblies, structural diversity in the pea genome and has been useful for the fine mapping of several discrete and quantitative traits.

Multi-year field trials provide a massive repository of trait data on a highly diverse population of tomato and uncover novel determinants of tomato productivity

Tomato (Solanum lycopersicum) is a prominent fruit with rich genetic resources for crop improvement. By using a phenotype-guided screen of over 7900 tomato accessions from around the world, we identified new associations for complex traits such as fruit weight and total soluble solids (Brix). Here, we present the phenotypic data from several years of trials. To illustrate the power of this dataset we use two case studies. First, evaluation of color revealed allelic variation in phytoene synthase 1 that resulted in differently colored or even bicolored fruit. Secondly, in view of the negative relationship between fruit weight and Brix, we pre-selected a subset of the collection that includes high and low Brix values in each category of fruit size. Genome-wide association analysis allowed us to detect novel loci associated with total soluble solid content and fruit weight. In addition, we developed eight F2 biparental intraspecific populations. Furthermore, by taking a phenotype-guided approach we were able to isolate individuals with high Brix values that were not compromised in terms of yield. In addition, the demonstration of novel results despite the high number of previous genome-wide association studies of these traits in tomato suggests that adoption of a phenotype-guided pre-selection of germplasm may represent a useful strategy for finding target genes for breeding.

Uncovering the genetic architecture of parallel evolution

Identifying the genetic architecture underlying adaptive traits is exceptionally challenging in natural populations. This is because associations between traits not only mask the targets of selection but also create correlated patterns of genomic divergence that hinder our ability to isolate causal genetic effects. Here, we examine the repeated evolution of components of the auxin pathway that have contributed to the replicated loss of gravitropism (i.e. the ability of a plant to bend in response to gravity) in multiple populations of the Senecio lautus species complex in Australia. We use a powerful approach which combines parallel population genomics with association mapping in a Multiparent Advanced Generation Inter-Cross (MAGIC) population to break down genetic and trait correlations to reveal how adaptive traits evolve during replicated evolution. We sequenced auxin and shoot gravitropism-related gene regions in 80 individuals from six natural populations (three parallel divergence events) and 133 individuals from a MAGIC population derived from two of the recently diverged natural populations. We show that artificial tail selection on gravitropism in the MAGIC population recreates patterns of parallel divergence in the auxin pathway in the natural populations. We reveal a set of 55 auxin gene regions that have evolved repeatedly during the evolution of the species, of which 50 are directly associated with gravitropism divergence in the MAGIC population. Our work creates a strong link between patterns of genomic divergence and trait variation contributing to replicated evolution by natural selection, paving the way to understand the origin and maintenance of adaptations in natural populations.

Development and evaluation of the utility of GenoBaits Peanut 40K for a peanut MAGIC population

Population and genotype data are essential for genetic mapping. The multi-parent advanced generation intercross (MAGIC) population is a permanent mapping population used for precisely mapping quantitative trait loci. Moreover, genotyping-by-target sequencing (GBTS) is a robust high-throughput genotyping technology characterized by its low cost, flexibility, and limited requirements for information management and support. In this study, an 8-way MAGIC population was constructed using eight elite founder lines. In addition, GenoBaits Peanut 40K was developed and utilized for the constructed MAGIC population. A subset (297 lines) of the MAGIC population at the S2 stage was genotyped using GenoBaits Peanut 40K. Furthermore, these lines and the eight parents were analyzed in terms of pod length, width, area, and perimeter. A total of 27 single nucleotide polymorphisms (SNPs) were revealed to be significantly associated with peanut pod size-related traits according to a genome-wide association study. The GenoBaits Peanut 40K provided herein and the constructed MAGIC population will be applicable for future research to identify the key genes responsible for important peanut traits.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01417-w.

Discovering useful genetic variation in the seed parent gene pool for sorghum improvement

Multi-parent populations contain valuable genetic material for dissecting complex, quantitative traits and provide a unique opportunity to capture multi-allelic variation compared to the biparental populations. A multi-parent advanced generation inter-cross (MAGIC) B-line (MBL) population composed of 708 F6 recombinant inbred lines (RILs), was recently developed from four diverse founders. These selected founders strategically represented the four most prevalent botanical races (kafir, guinea, durra, and caudatum) to capture a significant source of genetic variation to study the quantitative traits in grain sorghum [Sorghum bicolor (L.) Moench]. MBL was phenotyped at two field locations for seven yield-influencing traits: panicle type (PT), days to anthesis (DTA), plant height (PH), grain yield (GY), 1000-grain weight (TGW), tiller number per meter (TN) and yield per panicle (YPP). High phenotypic variation was observed for all the quantitative traits, with broad-sense heritabilities ranging from 0.34 (TN) to 0.84 (PH). The entire population was genotyped using Diversity Arrays Technology (DArTseq), and 8,800 single nucleotide polymorphisms (SNPs) were generated. A set of polymorphic, quality-filtered markers (3,751 SNPs) and phenotypic data were used for genome-wide association studies (GWAS). We identified 52 marker-trait associations (MTAs) for the seven traits using BLUPs generated from replicated plots in two locations. We also identified desirable allelic combinations based on the plant height loci (Dw1, Dw2, and Dw3), which influences yield related traits. Additionally, two novel MTAs were identified each on Chr1 and Chr7 for yield traits independent of dwarfing genes. We further performed a multi-variate adaptive shrinkage analysis and 15 MTAs with pleiotropic effect were identified. The five best performing MBL progenies were selected carrying desirable allelic combinations. Since the MBL population was designed to capture significant diversity for maintainer line (B-line) accessions, these progenies can serve as valuable resources to develop superior sorghum hybrids after validation of their general combining abilities via crossing with elite pollinators. Further, newly identified desirable allelic combinations can be used to enrich the maintainer germplasm lines through marker-assisted backcross breeding.

Transcriptional repression of MdMa1 by MdMYB21 in Ma locus decreases malic acid content in apple fruit

Malic acid is a major organic acid component of apples and a crucial determinant of fruit organoleptic quality. A candidate gene for malic acid content, designated MdMa1, was previously identified in the Ma locus, which is a major quantitative trait locus (QTL) for apple fruit acidity located on the linkage group 16. Region-based association mapping to detect candidate genes in the Ma locus identified MdMa1 and an additional MdMYB21 gene putatively associated with malic acid. MdMYB21 was significantly associated with fruit malic acid content, accounting for ~7.48% of the observed phenotypic variation in the apple germplasm collection. Analyses of transgenic apple calli, fruits and tomatoes demonstrated that MdMYB21 negatively regulated malic acid accumulation. The apple fruit acidity-related MdMa1 and its tomato ortholog, SlALMT9, exhibited lower expression profiles in apple calli, mature fruits and tomatoes in which MdMYB21 was overexpressed, compared with their corresponding wild-type variety. MdMYB21 directly binds to the MdMa1 promoter and represses its expression. Interestingly, a 2-bp variation in the MdMYB21 promoter region altered its expression and regulation of its target gene, MdMa1, expression. Our findings not only demonstrate the efficiency of integrating QTL and association mapping in the identification of candidate genes controlling complex traits in apples, but also provide insights into the complex regulatory mechanism of fruit malic acid accumulation.

Identification of QTLs and Putative Candidate Genes for Plant Architecture of Lotus Revealed by Regional Association Mapping

The lotus (Nelumbo Adans.) is one of the most economically relevant ornamental aquatic plants. Plant architecture (PA) is an important trait for lotus classification, cultivation, breeding, and applications. However, the underlying genetic and molecular basis controlling PA remains poorly understood. In this study, an association study for PA-related traits was performed with 93 genome-wide microsatellite markers (simple sequence repeat, SSR) and 51 insertion-deletion (InDel) markers derived from the candidate regions using a panel of 293 lotus accessions. Phenotypic data analysis of the five PA-related traits revealed a wide normal distribution and high heritability from 2013 to 2016, which indicated that lotus PA-related traits are highly polygenic traits. The population structure (Q-matrix) and the relative kinships (K-matrix) of the association panels were analyzed using 93 SSR markers. The mixed linear model (MLM) taking Q-matrix and K-matrix into account was used to estimate the association between markers and the traits. A total of 26 markers and 65 marker-trait associations were identified by considering associations with p < 0.001 and Q < 0.05. Based on the significant markers, two QTLs on Chromosome 1 were identified, and two candidate genes were preliminarily determined. The results of our study provided useful information for the lotus breeding aiming at different PA phenotypes using a molecular-assisted selection (MAS) method and also laid the foundation for the illustration of the molecular mechanism underlying the major QTL and key markers associated with lotus PA.

A Genome-Wide Association Study and Genomic Prediction for Fiber and Sucrose Contents in a Mapping Population of LCP 85-384 Sugarcane

Sugarcane (Saccharum spp. hybrids) is an economically important crop for both sugar and biofuel industries. Fiber and sucrose contents are the two most critical quantitative traits in sugarcane breeding that require multiple-year and multiple-location evaluations. Marker-assisted selection (MAS) could significantly reduce the time and cost of developing new sugarcane varieties. The objectives of this study were to conduct a genome-wide association study (GWAS) to identify DNA markers associated with fiber and sucrose contents and to perform genomic prediction (GP) for the two traits. Fiber and sucrose data were collected from 237 self-pollinated progenies of LCP 85-384, the most popular Louisiana sugarcane cultivar from 1999 to 2007. The GWAS was performed using 1310 polymorphic DNA marker alleles with three models of TASSEL 5, single marker regression (SMR), general linear model (GLM) and mixed linear model (MLM), and the fixed and random model circulating probability unification (FarmCPU) of R package. The results showed that 13 and 9 markers were associated with fiber and sucrose contents, respectively. The GP was performed by cross-prediction with five models, ridge regression best linear unbiased prediction (rrBLUP), Bayesian ridge regression (BRR), Bayesian A (BA), Bayesian B (BB) and Bayesian least absolute shrinkage and selection operator (BL). The accuracy of GP varied from 55.8% to 58.9% for fiber content and 54.6% to 57.2% for sucrose content. Upon validation, these markers can be applied in MAS and genomic selection (GS) to select superior sugarcane with good fiber and high sucrose contents.

Genome-Wide association analysis of phenotypic traits in Bambara groundnut under drought-stressed and non-stressed conditions based on DArTseq SNP

INTRODUCTION

Bambara groundnut (BG) (Vigna subterranea [L.] Verdc) is an indigenous, resilient, but underutilized leguminous crop that occurs mostly as genetically heterogeneous landraces with limited information on the drought tolerant attributes. This study elucidates the associations between sequencing-based diversity array technology (DArTseq) and phenotypic character as well as differing indices related to drought tolerance in one hundred accessions of Bambara groundnut.

METHODS

The field experiments were conducted at IITA research stations in Kano and Ibadan between 2016 and 2018 planting seasons. The experiments were arranged in randomised complete block design with three replications, under the different water regimes. The phenotypic traits evaluated was further to construct the dendrogram. Genome-wide association mapping was conducted based on 5927 DArTs loci with < 20% missing data.

RESULTS AND DISCUSSIONS

The genome wide association study predicted drought tolerance in Bambara accessions for geometric mean productivity (GMP) and stress tolerance index (STI). TVSu-423 had the highest GMP and STI values (28.50, 2.40), while TVSu-2017 had the lowest at GMP (1.74) and STI (0.01) respectively. The relative water content (%) was significantly higher for accessions; TVSu-266 (60.35, 61.49), TVSu-2 (58.29, 53.94), and TVSu-411 (55.17, 58.92) in 2016/2017 and 2017/2018, respectively. The phenotypic characters studied delineated the accessions into two major clusters and five distinct sub-clusters, indicating variations across all the geographical locations. The 5,927 DArTseq genomic markers in association with STI further grouped the 100 accessions into two main clusters. TVSu-1897 from Botswana (Southern Africa) was in the first cluster, while the remaining 99 accessions from Western, Central, and Eastern Africa made up the second cluster. The eight significant Quantitative Trait Loci (QTLs) (24346377|F|0-22:A>G-22:A>G, 24384105|F|0-56:A>G33 :A> G, 24385643|F|0-53:G>C-53:G>C, 24385696|F|0-43:A>G-43:A>G, 4177257|F|0-44:A>T-44:A>T, 4182070|F|0-66:G>A-66:G>A, 4183483|F|0-24:G>A-24:G>A, 4183904|F|0-11:C>T-11:C>T) identified with Bonferroni threshold was in association with STI, indicative of variations under the drought-stressed condition. The observation of consistent SNPs in the 2016 and 2017 planting seasons, as well as in combination with the 2016 and 2017 planting seasons, led to the designation of these QTLs as significant. The drought selected accessions could form basis for hybridization breeding. The identified quantitative trait loci could be useful in marker-assisted selection in drought molecular breeding programs.

Conventional and new-breeding technologies for improving disease resistance in lentil (Lens culinaris Medik)

Lentil, an important cool season food legume, is a rich source of easily digestible protein, folic acid, bio-available iron, and zinc nutrients. Lentil grows mainly as a sole crop in the winter after harvesting rice in South Asia. However, the annual productivity is low due to its slow growth during the early phase, competitive weed infestation, and disease outbreaks during the crop growth period. Disease resistance breeding has been practiced for a long time to enhance resistance to various diseases. Often the sources of resistance are available in wild crop relatives. Thus, wide hybridization and the ovule rescue technique have helped to introgress the resistance trait into cultivated lentils. Besides hybridization, induced mutagenesis contributed immensely in creating variability for disease tolerance, and several disease-resistant mutant lines have been developed. However, to overcome the limitations of traditional breeding approaches, advancement in molecular marker technologies, and genomics has helped to develop disease-resistant and climate-resilient lentil varieties with more precision and efficiency. This review describes types of diseases, disease screening methods, the role of conventional and new breeding technologies in alleviating disease-incurred damage and progress toward making lentil varieties more resilient to disease outbreaks under the shadow of climate change.

Development of a MAGIC population and high-resolution quantitative trait mapping for nicotine content in tobacco

Multiparent Advanced Generation Inter-Cross (MAGIC) population is an ideal genetic and breeding material for quantitative trait locus (QTL) mapping and molecular breeding. In this study, a MAGIC population derived from eight tobacco parents was developed. Eight parents and 560 homozygous lines were genotyped by a 430K single-nucleotide polymorphism (SNP) chip assay and phenotyped for nicotine content under different conditions. Four QTLs associated with nicotine content were detected by genome-wide association mapping (GWAS), and one major QTL, named qNIC7-1, was mapped repeatedly under different conditions. Furthermore, by combining forward mapping, bioinformatics analysis and gene editing, we identified an ethylene response factor (ERF) transcription factor as a candidate gene underlying the major QTL qNIC7-1 for nicotine content in tobacco. A presence/absence variation (PAV) at qNIC7-1 confers changes in nicotine content. Overall, the large size of this MAGIC population, diverse genetic composition, balanced parental contributions and high levels of recombination all contribute to its value as a genetic and breeding resource. The application of the tobacco MAGIC population for QTL mapping and detecting rare allelic variation was demonstrated using nicotine content as a proof of principle.

Genetic diversity and population structure assessment of Western Canadian barley cooperative trials

Studying the population structure and genetic diversity of historical datasets is a proposed use for association analysis. This is particularly important when the dataset contains traits that are time-consuming or costly to measure. A set of 96 elite barley genotypes, developed from eight breeding programs of the Western Canadian Cooperative Trials were used in the current study. Genetic diversity, allelic variation, and linkage disequilibrium (LD) were investigated using 5063 high-quality SNP markers via the Illumina 9K Barley Infinium iSelect SNP assay. The distribution of SNPs markers across the barley genome ranged from 449 markers on chromosome 1H to 1111 markers on chromosome 5H. The average polymorphism information content (PIC) per locus was 0.275 and ranged from 0.094 to 0.375. Bayesian clustering in STRUCTURE and principal coordinate analysis revealed that the populations are differentiated primarily due to the different breeding program origins and ear-row type into five subpopulations. Analysis of molecular variance based on PhiPT values suggested that high values of genetic diversity were observed within populations and accounted for 90% of the total variance. Subpopulation 5 exhibited the most diversity with the highest values of the diversity indices, which represent the breeding program gene pool of AFC, AAFRD, AU, and BARI. With increasing genetic distance, the LD values, expressed as r², declined to below the critical r^{2 =} 0.18 after 3.91 cM, and the same pattern was observed on each chromosome. Our results identified an important pattern of genetic diversity among the Canadian barley panel that was proposed to be representative of target breeding programs and may have important implications for association mapping in the future. This highlight, that efforts to identify novel variability underlying this diversity may present practical breeding opportunities to develop new barley genotypes.

Analysis of genotype × environment interactions for agronomic traits of soybean (Glycine max [L.] Merr.) using association mapping

The soybean yield is a complex quantitative trait that is significantly influenced by environmental factors. G × E interaction (GEI), which derives the performance of soybean genotypes differentially in various environmental conditions, is one of the main obstacles to increasing the net production. The primary goal of this study is to identify the outperforming genotypes in different latitudes, which can then be used in future breeding programs. A total of 96 soybean genotypes were examined in two different ecological regions: Faisalabad and Tando Jam in Pakistan. The evaluation of genotypes in different environmental conditions showed a substantial amount of genetic diversity for grain yield. We identified 13 environment-specific genotypes showing their maximum grain yield in each environment. Genotype G69 was found to be an ideal genotype with higher grain yield than other genotypes tested in this study and is broadly adapted for environments E1 and E2 and also included in top-yielding genotypes in E3, E4, and E5. G92 is another genotype that is broadly adapted in E1, E3, and E4. In the case of environments, E3 is suggested to be a more ideal environment as it is plotted near the concentric circle and is very informative for the selection of genotypes with high yield. Despite the presence of GEI, advances in DNA technology provided very useful tools to investigate the insight of advanced genotypes. Association mapping is a useful method for swiftly and efficiently investigating the genetic basis of significant plant traits. A total of 26 marker-trait associations were found for six agronomic traits in five environments, with the highest significance (p-value = 2.48 × 10^-08) for plant height and the lowest significance (1.03 × 10^-03) for hundred-grain weight. Soybean genotypes identified in the present study could be a valuable source for future breeding programs as they are adaptable to a wide range of environments. Genetic selection of genotypes with the best yields can be used for gross grain production in a wide range of climatic conditions, and it would give an essential reference in terms of soybean variety selection.

Genomic Approaches for the Study of Flower Development in Floriculture Crops

The advances in genomics and bioinformatics have made possible the study in non-model plants of phenotypes associated to flower development. Floriculture crops are an interesting source of traits associated to flower development such as the transition between zygomorphic and actinomorphic flowers or the production of flowers with double and triple corollas. In this chapter, we summarize the material and methods for the use of floriculture crops to study flower development using genomic tools, from the sequencing and assembly of a reference genome to QTL and RNA-Seq analysis to search candidate genes associated to specific traits.

Association mapping for important agronomic traits in wild and cultivated Vigna species using cross-species and cross-genera simple sequence repeat markers

The genus Vigna is an agronomically important taxon, with many of its species inhabiting a wide range of environments and offering numerous useful genes for the improvement of the cultivated types. The present study aimed to detect the genomic regions associated with yield-attributing traits by genome-wide association mapping. A diverse panel of 98 wild and cultivated Vigna accessions (acc.) belonging to 13 different species was evaluated for yield and related traits during the kharif season of 2017 and 2018. The panel was also genotyped using 92 cross-genera and cross-species simple sequence repeat markers to study the population genetic structure and useful market-trait associations. The PCA and trait correlation established relationships amongst the traits during both seasons while 100-seed weight (HSW) had a positive correlation with pod length (PL), and days to first flowering (DFF) with days to maturity (DM). The population genetic structure analysis grouped different acc. into three genetically distinct sub-populations with SP-1 comprising 34 acc., SP-2 (24 acc.), and SP-3 (33 acc.) and one admixture group (7 acc.). Mixed linear model analysis revealed an association of 13 markers, namely, VR018, VR039, VR022, CEDG033, GMES0337, MBSSR008, CEDG220, VM27, CP1225, CP08695, CEDG100, CEDG008, and CEDG096A with nine traits. Seven of the aforementioned markers, namely, VR018 for plant height (PH) and terminal leaflet length (TLL), VR022 for HSW and pod length (PL), CEDG033 for DFF and DM, MBSSR008 for DFF and DM, CP1225 for CC at 30 days (CC30), DFF and DM, CEDG100 for PH and terminal leaflet length (TLL), and CEDG096A for CC30 and chlorophyll content at 45 days were associated with multiple traits. The marker CEDG100, associated with HSW, PH, and TLL, is co-localized in gene-encoding histone-lysine N-methyltransferase ATX5. Similarly, VR22, associated with PL and HSW, is co-located in gene-encoding SHOOT GRAVITROPISM 5 in mungbean. These associations may be highly useful for marker-assisted genetic improvement of mungbean and other related Vigna species.

Genomic selection in tropical perennial crops and plantation trees: a review

To overcome the multiple challenges currently faced by agriculture, such as climate change and soil deterioration, more efficient plant breeding strategies are required. Genomic selection (GS) is crucial for the genetic improvement of quantitative traits, as it can increase selection intensity, shorten the generation interval, and improve selection accuracy for traits that are difficult to phenotype. Tropical perennial crops and plantation trees are of major economic importance and have consequently been the subject of many GS articles. In this review, we discuss the factors that affect GS accuracy (statistical models, linkage disequilibrium, information concerning markers, relatedness between training and target populations, the size of the training population, and trait heritability) and the genetic gain expected in these species. The impact of GS will be particularly strong in tropical perennial crops and plantation trees as they have long breeding cycles and constrained selection intensity. Future GS prospects are also discussed. High-throughput phenotyping will allow constructing of large training populations and implementing of phenomic selection. Optimized modeling is needed for longitudinal traits and multi-environment trials. The use of multi-omics, haploblocks, and structural variants will enable going beyond single-locus genotype data. Innovative statistical approaches, like artificial neural networks, are expected to efficiently handle the increasing amounts of heterogeneous multi-scale data. Targeted recombinations on sites identified from profiles of marker effects have the potential to further increase genetic gain. GS can also aid re-domestication and introgression breeding. Finally, GS consortia will play an important role in making the best of these opportunities.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01326-4.

Genome-wide association analysis of 101 accessions dissects the genetic basis of shell thickness for genetic improvement in Persian walnut (Juglans regia L.)

BACKGROUND

Understanding the underlying genetic mechanisms that drive phenotypic variations is essential for enhancing the efficacy of crop improvement. Persian walnut (Juglans regia L.), which is grown extensively worldwide, is an important economic tree fruit due to its horticultural, medicinal, and material value. The quality of the walnut fruit is related to the selection of traits such as thinner shells, larger filling rates, and better taste, which is very important for breeding in China. The complex quantitative fruit-related traits are influenced by a variety of physiological and environmental factors, which can vary widely between walnut genotypes.

RESULTS

For this study, a set of 101 Persian walnut accessions were re-sequenced, which generated a total of 906.2 Gb of Illumina sequence data with an average read depth of 13.8× for each accession. We performed the genome-wide association study (GWAS) using 10.9 Mb of high-quality single-nucleotide polymorphisms (SNPs) and 10 agronomic traits to explore the underlying genetic basis of the walnut fruit. Several candidate genes are proposed to be involved in walnut characteristics, including JrPXC1, JrWAKL8, JrGAMYB, and JrFRK1. Specifically, the JrPXC1 gene was confirmed to participate in the regulation of secondary wall cellulose thickening in the walnut shell.

CONCLUSION

In addition to providing considerable available genetic resources for walnut trees, this study revealed the underlying genetic basis involved in important walnut agronomic traits, particularly shell thickness, as well as providing clues for the improvement of genetic breeding and domestication in other perennial economic crops.

Analysis of historical selection in winter wheat

KEY MESSAGE

Modeling of the distribution of allele frequency over year of variety release identifies major loci involved in historical breeding of winter wheat. Winter wheat is a major crop with a rich selection history in the modern era of crop breeding. Genetic gains across economically important traits like yield have been well characterized and are the major force driving its production. Winter wheat is also an excellent model for analyzing historical genetic selection. As a proof of concept, we analyze two major collections of winter wheat varieties that were bred in Western Europe from 1916 to 2010, namely the Triticeae Genome (TG) and WAGTAIL panels, which include 333 and 403 varieties, respectively. We develop and apply a selection mapping approach, Regression of Alleles on Years (RALLY), in these panels, as well as in simulated populations. RALLY maps loci under sustained historical selection by using a simple logistic model to regress allele counts on years of variety release. To control for drift-induced allele frequency change, we develop a hybrid approach of genomic control and delta control. Within the TG panel, we identify 22 significant RALLY quantitative selection loci (QSLs) and estimate the local heritabilities for 12 traits across these QSLs. By correlating predicted marker effects with RALLY regression estimates, we show that alleles whose frequencies have increased over time are heavily biased toward conferring positive yield effect, but negative effects in flowering time, lodging, plant height and grain protein content. Altogether, our results (1) demonstrate the use of RALLY to identify selected genomic regions while controlling for drift, and (2) reveal key patterns in the historical selection in winter wheat and guide its future breeding.

Exploring Chickpea Germplasm Diversity for Broadening the Genetic Base Utilizing Genomic Resourses

Legume crops provide significant nutrition to humans as a source of protein, omega-3 fatty acids as well as specific macro and micronutrients. Additionally, legumes improve the cropping environment by replenishing the soil nitrogen content. Chickpeas are the second most significant staple legume food crop worldwide behind dry bean which contains 17%–24% protein, 41%–51% carbohydrate, and other important essential minerals, vitamins, dietary fiber, folate, β-carotene, anti-oxidants, micronutrients (phosphorus, calcium, magnesium, iron, and zinc) as well as linoleic and oleic unsaturated fatty acids. Despite these advantages, legumes are far behind cereals in terms of genetic improvement mainly due to far less effort, the bottlenecks of the narrow genetic base, and several biotic and abiotic factors in the scenario of changing climatic conditions. Measures are now called for beyond conventional breeding practices to strategically broadening of narrow genetic base utilizing chickpea wild relatives and improvement of cultivars through advanced breeding approaches with a focus on high yield productivity, biotic and abiotic stresses including climate resilience, and enhanced nutritional values. Desirable donors having such multiple traits have been identified using core and mini core collections from the cultivated gene pool and wild relatives of Chickpea. Several methods have been developed to address cross-species fertilization obstacles and to aid in inter-specific hybridization and introgression of the target gene sequences from wild Cicer species. Additionally, recent advances in “Omics” sciences along with high-throughput and precise phenotyping tools have made it easier to identify genes that regulate traits of interest. Next-generation sequencing technologies, whole-genome sequencing, transcriptomics, and differential genes expression profiling along with a plethora of novel techniques like single nucleotide polymorphism exploiting high-density genotyping by sequencing assays, simple sequence repeat markers, diversity array technology platform, and whole-genome re-sequencing technique led to the identification and development of QTLs and high-density trait mapping of the global chickpea germplasm. These altogether have helped in broadening the narrow genetic base of chickpeas.

Genome properties of key oil palm (Elaeis guineensis Jacq.) breeding populations

A good knowledge of the genome properties of the populations makes it possible to optimize breeding methods, in particular genomic selection (GS). In oil palm (Elaeis guineensis Jacq), the world's main source of vegetable oil, this would provide insight into the promising GS results obtained so far. The present study considered two complex breeding populations, Deli and La Mé, with 943 individuals and 7324 single-nucleotide polymorphisms (SNPs) from genotyping-by-sequencing. Linkage disequilibrium (LD), haplotype sharing, effective size (N_e), and fixation index (F_st) were investigated. A genetic linkage map spanning 1778.52 cM and with a recombination rate of 2.85 cM/Mbp was constructed. The LD at r²=0.3, considered the minimum to get reliable GS results, spanned over 1.05 cM/0.22 Mbp in Deli and 0.9 cM/0.21 Mbp in La Mé. The significant degree of differentiation existing between Deli and La Mé was confirmed by the high F_st value (0.53), the pattern of correlation of SNP heterozygosity and allele frequency among populations, and the decrease of persistence of LD and of haplotype sharing among populations with increasing SNP distance. However, the level of resemblance between the two populations over short genomic distances (correlation of r values between populations >0.6 for SNPs separated by <0.5 cM/1 kbp and percentage of common haplotypes >40% for haplotypes <3600 bp/0.20 cM) likely explains the superiority of GS models ignoring the parental origin of marker alleles over models taking this information into account. The two populations had low N_e (<5). Population-specific genetic maps and reference genomes are recommended for future studies.

Using breeding and quantitative genetics to understand the C4 pathway

Reducing photorespiration in C3 crops could significantly increase rates of photosynthesis and yield. One method to achieve this would be to integrate C4 photosynthesis into C3 species. This objective is challenging as it involves engineering incompletely understood traits into C3 leaves, including complex changes to their biochemistry, cell biology, and anatomy. Quantitative genetics and selective breeding offer underexplored routes to identify regulators of these processes. We first review examples of natural intraspecific variation in C4 photosynthesis as well as the potential for hybridization between C3 and C4 species. We then discuss how quantitative genetic approaches including artificial selection and genome-wide association could be used to better understand the C4 syndrome and in so doing guide the engineering of the C4 pathway into C3 crops.

Screening of stable resistant accessions and identification of resistance loci to Barley yellow mosaic virus disease

Background

The disease caused by Barley yellow mosaic virus (BaYMV) infection is a serious threat to autumn-sown barley (Hordeum vulgare L.) production in Europe, East Asia and Iran. Due to the rapid diversification of BaYMV strains, it is urgent to discover novel germplasm and genes to assist breeding new varieties with resistance to different BaYMV strains, thus minimizing the effect of BaYMV disease on barley cropping.

Methods

A natural population consisting of 181 barley accessions with different levels of resistance to BaYMV disease was selected for field resistance identification in two separate locations (Yangzhou and Yancheng, Jiangsu Province, China). Additive main effects and multiplicative interaction (AMMI) analysis was used to identify accessions with stable resistance. Genome-wide association study (GWAS) of BaYMV disease resistance was broadly performed by combining both single nucleotide polymorphisms (SNPs) and specific molecular markers associated with the reported BaYMV disease resistance genes. Furthermore, the viral protein genome linked (VPg) sequences of the virus were amplified and analyzed to assess the differences between the BaYMV strains sourced from the different experimental sites.

Results

Seven barley accessions with lower standardized Area Under the Disease Progress Steps (sAUDPS) index in every environment were identified and shown to have stable resistance to BaYMV disease in each assessed location. Apart from the reported BaYMV disease resistance genes rym4 and rym5, one novel resistance locus explaining 24.21% of the phenotypic variation was identified at the Yangzhou testing site, while two other novel resistance loci that contributed 19.23% and 19.79% of the phenotypic variation were identified at the Yancheng testing site, respectively. Further analysis regarding the difference in the VPg sequence of the predominant strain of BaYMV collected from these two testing sites may explain the difference of resistance loci differentially identified under geographically distinct regions. Our research provides novel genetic resources and resistance loci for breeding barley varieties for BaMYV disease resistance.

Germplasm Resources and Strategy for Genetic Breeding of Lycium Species: A Review

Lycium species (goji), belonging to Solanaceae, are widely spread in the arid to semiarid environments of Eurasia, Africa, North and South America, among which most species have affinal drug and diet functions, resulting in their potential to be a superior healthy food. However, compared with other crop species, scientific research on breeding Lycium species lags behind. This review systematically introduces the present germplasm resources, cytological examination and molecular-assisted breeding progress in Lycium species. Introduction of the distribution of Lycium species around the world could facilitate germplasm collection for breeding. Karyotypes of different species could provide a feasibility analysis of fertility between species. The introduction of mapping technology has discussed strategies for quantitative trait locus (QTL) mapping in Lycium species according to different kinds of traits. Moreover, to extend the number of traits and standardize the protocols of trait detection, we also provide 1,145 potential traits (275 agronomic and 870 metabolic) in different organs based on different reference studies on Lycium, tomato and other Solanaceae species. Finally, perspectives on goji breeding research are discussed and concluded. This review will provide breeders with new insights into breeding Lycium species.

Association Mapping for Yield Attributing Traits and Yellow Mosaic Disease Resistance in Mung Bean [Vigna radiata (L.) Wilczek]

Mung bean [Vigna radiata (L.) Wilczek] is an important short-duration grain legume widely known for its nutritional, soil ameliorative, and cropping system intensification properties. This study aims at evaluating genetic diversity among mung bean genotypes and detecting genomic regions associated with various yield attributing traits and yellow mosaic disease (YMD) resistance by association mapping. A panel of 80 cultivars and advanced breeding lines was evaluated for 10 yield-related and YMD resistance traits during kharif (monsoon) and summer seasons of 2018-2019 and 2019-2020. A total of 164 genome-wide simple sequence repeat (SSR) markers were initially screened, out of which 89 were found polymorphic which generated 317 polymorphic alleles with an average of 3.56 alleles per SSR locus. The number of alleles at each locus varied from 2 to 7. The population genetic structure analysis grouped different genotypes in three major clusters and three genetically distinct subpopulations (SPs) (i.e., SP-1, SP-2, and SP-3) with one admixture subpopulation (SP-4). Both cluster and population genetic structure analysis categorized the advanced mung bean genotypes in a single group/SP and the released varieties in other groups/SPs, suggesting that the studied genotypes may have common ancestral history at some level. The population genetic structure was also in agreement with the genetic diversity analysis. The estimate of the average degree of linkage disequilibrium (LD) present at the genome level in 80 mung bean genotypes unveiled significant LD blocks. Over the four seasons, 10 marker-trait associations were observed significant for YMD and four seed yield (SY)-related traits viz., days to flowering, days to maturity, plant height, and number of pods per plant using the mixed linear model (MLM) method. These associations may be useful for marker-assisted mung bean yield improvement programs and YMD resistance.

Genomics of maize resistance to kernel contamination with fumonisins using a multiparental advanced generation InterCross maize population (MAGIC)

Maize kernel is exposed to several fungal species, most notably Fusarium verticillioides, which can contaminate maize kernels with fumonisins. In an effort to increase genetic gains and avoid the laborious tasks of conventional breeding, the use of marker-assisted selection or genomic selection programs was proposed. To this end, in the present study a Genome Wide Association Study (GWAS) was performed on 339 RILs of a Multiparental Advanced Generation InterCross (MAGIC) population that had previously been used to locate Quantitative Trait Locus (QTL) for resistance to Fusarium Ear Rot (FER). Six QTLs for fumonisin content were detected in the bins 3.08, 4.07, 4.10, 7.03-7.04, 9.04-9.05 and 10.04-10.5. Five of the six QTLs collocate in regions where QTLs for FER were also found. However, the genetic variation for fumonisin content in kernel is conditioned by many other QTLs of small effect that could show QTL x environment interaction effects. Although a genomic selection approach to directly reduce fumonisin content in the kernel could be suitable, improving resistance to fumonisin content by genomic selection for FER would be more advisable.

Nucleotide Diversity of the Maize ZmCNR13 Gene and Association With Ear Traits

The maize (Zea mays L.) ZmCNR13 gene, encoding a protein of fw2.2-like (FWL) family, has been demonstrated to be involved in cell division, expansion, and differentiation. In the present study, the genomic sequences of the ZmCNR13 locus were re-sequenced in 224 inbred lines, 56 landraces and 30 teosintes, and the nucleotide polymorphism and selection signature were estimated. A total of 501 variants, including 415 SNPs and 86 Indels, were detected. Among them, 51 SNPs and 4 Indels were located in the coding regions. Although neutrality tests revealed that this locus had escaped from artificial selection during the process of maize domestication, the population of inbred lines possesses lower nucleotide diversity and decay of linkage disequilibrium. To estimate the association between sequence variants of ZmCNR13 and maize ear characteristics, a total of ten ear-related traits were obtained from the selected inbred lines. Four variants were found to be significantly associated with six ear-related traits. Among them, SNP2305, a non-synonymous mutation in exon 2, was found to be associated with ear weight, ear grain weight, ear diameter and ear row number, and explained 4.59, 4.61, 4.31, and 8.42% of the phenotypic variations, respectively. These results revealed that natural variations of ZmCNR13 might be involved in ear development and can be used in genetic improvement of maize ear-related traits.

Accelerating wheat breeding for end-use quality through association mapping and multivariate genomic prediction

In hard-winter wheat (Triticum aestivum L.) breeding, the evaluation of end-use quality is expensive and time-consuming, being relegated to the final stages of the breeding program after selection for many traits including disease resistance, agronomic performance, and grain yield. In this study, our objectives were to identify genetic variants underlying baking quality traits through genome-wide association study (GWAS) and develop improved genomic selection (GS) models for the quality traits in hard-winter wheat. Advanced breeding lines (n = 462) from 2015-2017 were genotyped using genotyping-by-sequencing (GBS) and evaluated for baking quality. Significant associations were detected for mixograph mixing time and bake mixing time, most of which were within or in tight linkage to glutenin and gliadin loci and could be suitable for marker-assisted breeding. Candidate genes for newly associated loci are phosphate-dependent decarboxylase and lipid transfer protein genes, which are believed to affect nitrogen metabolism and dough development, respectively. The use of GS can both shorten the breeding cycle time and significantly increase the number of lines that could be selected for quality traits, thus we evaluated various GS models for end-use quality traits. As a baseline, univariate GS models had 0.25-0.55 prediction accuracy in cross-validation and from 0 to 0.41 in forward prediction. By including secondary traits as additional predictor variables (univariate GS with covariates) or correlated response variables (multivariate GS), the prediction accuracies were increased relative to the univariate model using only genomic information. The improved genomic prediction models have great potential to further accelerate wheat breeding for end-use quality.

Identification of Novel Quantitative Trait Nucleotides and Candidate Genes for Bacterial Wilt Resistance in Tobacco (Nicotiana tabacum L.) Using Genotyping-by-Sequencing and Multi-Locus Genome-Wide Association Studies

Tobacco bacterial wilt (TBW) is a devastating soil-borne disease threatening the yield and quality of tobacco. However, its genetic foundations are not fully understood. In this study, we identified 126,602 high-quality single-nucleotide polymorphisms (SNPs) in 94 tobacco accessions using genotyping-by-sequencing (GBS) and a 94.56 KB linkage disequilibrium (LD) decay rate for candidate gene selection. The population structure analysis revealed two subpopulations with 37 and 57 tobacco accessions. Four multi-locus genome-wide association study (ML-GWAS) approaches identified 142 quantitative trait nucleotides (QTNs) in E1-E4 and the best linear unbiased prediction (BLUP), explaining 0.49-22.52% phenotypic variance. Of these, 38 novel stable QTNs were identified across at least two environments/methods, and their alleles showed significant TBW-DI differences. The number of superior alleles associated with TBW resistance for each accession ranged from 4 to 24; eight accessions had more than 18 superior alleles. Based on TBW-resistant alleles, the five best cross combinations were predicted, including MC133 × Ruyuan No. 1 and CO258 × ROX28. We identified 52 candidate genes around 38 QTNs related to TBW resistance based on homologous functional annotation and KEGG enrichment analysis, e.g., CYCD3;2, BSK1, Nitab4.5_0000641g0050, Nitab4.5_0000929g0030. To the best of our knowledge, this is the first comprehensive study to identify QTNs, superior alleles, and their candidate genes for breeding TBW-resistant tobacco varieties. The results provide further insight into the genetic architecture, marker-assisted selection, and functional genomics of TBW resistance, improving future breeding efforts to increase crop productivity.

Crafting for a better MAGIC: systematic design and test for Multiparental Advanced Generation Inter-Cross population

Multiparental Advanced Generation Inter-Cross (MAGIC) populations are valuable crop resources with a wide array of research uses including genetic mapping of complex traits, management of genetic resources and breeding of new varieties. Multiple founders are crossed to create a rich mosaic of highly recombined founder genomes in the MAGIC recombinant inbred lines (RILs). Many variations of MAGIC population designs exist; however, a large proportion of the currently available populations have been created empirically and based on similar designs. In our evaluations of five MAGIC populations, we found that the choice of designs has a large impact on the recombination landscape in the RILs. The most popular design used in many MAGIC populations has been shown to have a bias in recombinant haplotypes and low level of unique recombinant haplotypes, and therefore is not recommended. To address this problem and provide a remedy for the future, we have developed the “magicdesign” R package for creating and testing any MAGIC population design via simulation. A Shiny app version of the package is available as well. Our “magicdesign” package provides a unifying tool and a framework for creativity and innovation in MAGIC population designs. For example, using this package, we demonstrate that MAGIC population designs can be found which are very effective in creating haplotype diversity without the requirement for very large crossing programs. Furthermore, we show that interspersing cycles of crossing with cycles of selfing is effective in increasing haplotype diversity. These approaches are applicable in species that are hard to cross or in which resources are limited.

GWAS for genetics of complex quantitative traits: Genome to pangenome and SNPs to SVs and k-mers

The development of improved methods for genome-wide association studies (GWAS) for genetics of quantitative traits has been an active area of research during the last 25 years. This activity initially started with the use of mixed linear model (MLM), which was variously modified. During the last decade, however, with the availability of high throughput next generation sequencing (NGS) technology, development and use of pangenomes and novel markers including structural variations (SVs) and k-mers for GWAS has taken over as a new thrust area of research. Pangenomes and SVs are now available in humans, livestock, and a number of plant species, so that these resources along with k-mers are being used in GWAS for exploring additional genetic variation that was hitherto not available for analysis. These developments have resulted in significant improvement in GWAS methodology for detection of marker-trait associations (MTAs) that are relevant to human healthcare and crop improvement.

Genome-wide association mapping for agronomic traits in an 8-way Upland cotton MAGIC population by SLAF-seq

One sub-MAGIC population was genotyped using SLAF-seq, and QTLs and candidate genes for agronomic traits were identified in Upland cotton. The agronomic traits of Upland cotton have serious impacts on cotton production, as well as economic benefits. To discover the genetic basis of important agronomic traits in Upland cotton, a subset MAGIC (multi-parent advanced generation inter-cross) population containing 372 lines (SMLs) was selected from an 8-way MAGIC population with 960 lines. The 372 lines and 8 parents were phenotyped in six environments and deeply genotyped by SLAF-seq with 60,495 polymorphic SNPs. The genetic diversity indexes of all SNPs were 0.324 and 0.362 for the parents and MAGIC lines, respectively. The LD decay distance of the SMLs was 600 kb (r² = 0.1). Genome-wide association mapping was performed using 60,495 SNPs and the phenotypic data of the SMLs, and 177 SNPs were identified to be significantly associated with 9 stable agronomic traits in multiple environments. The identified SNPs were divided into 117 QTLs (quantitative trait loci) by LD decay distance, explaining 5.44% to 31.64% of the phenotypic variation. Among the 117 QTLs, 3 QTLs were stable in multiple environments, and 11 QTL regions were proven to have pleiotropism associated with multiple traits. Within QTL regions, 154 genes were preferentially expressed in correlated tissues, and 8 genes with known functions were identified as priori candidate genes. Two genes, GhACT1 and GhGASL3, reported to have clear functions, were, respectively, located in qFE-A05-4 and qFE-D04-3, two stable QTLs for FE. This study revealed the genetic basis of important agronomic traits of Upland cotton, and the results will facilitate molecular breeding in cotton.

Toward the identification of molecular markers associated with phytochemical traits in the Iranian sumac (Rhus coriaria L.) population

Sumac (Rhus coriaria L.) is one of the important forest species dispersed in the northwest of Iran. It is one of the important spice in Iran and the Middle East because of active components containing organic acids, phenolic acids, flavonoids, anthocyanins, tannins and terpenoids. This study aimed to investigate population structure and linkage disequilibrium (LD) extent within Rhus coriaria L. genotypes using ISSR markers and identify molecular markers associated with phytochemical traits using association analysis. In the molecular part of the experiment, the genetic diversity of 75 sumac genotypes from five different areas of northwest Iran was assessed by 18 ISSR primers. In the phenotypic assessment, the fruits of the sumac genotypes were analyzed using HPLC-LC/MSMS for determining phytochemical components including maleic acid, ellagic acid, maleic acid hexoside, gallic acid, coumaric acid, quercetin, caftaric acid, and linoleic acid. The phenotypic data analysis revealed the great phenotypic diversity among and within Iranian sumac populations for the studied phytochemical traits. The studied sumac genotypes were divided into two subpopulations based on molecular marker-based structure analysis. A significant level of LD was observed in 11.64% of the ISSR marker pairs (p < .05). The mixed linear model procedure showed that 12 loci had a significant association with investigated traits. The ISSR loci identified in this study can be potentially used in marker-assisted selection in sumac breeding programs.

Genome-Wide Association Study Uncover the Genetic Architecture of Salt Tolerance-Related Traits in Common Wheat (Triticum aestivum L.)

Soil salinity is a serious threat to wheat yield affecting sustainable agriculture. Although salt tolerance is important for plant establishment at seedling stage, its genetic architecture remains unclear. In the present study, we have evaluated eight salt tolerance-related traits at seedling stage and identified the loci for salt tolerance by genome-wide association study (GWAS). This GWAS panel comprised 317 accessions and was genotyped with the wheat 90 K single-nucleotide polymorphism (SNP) chip. In total, 37 SNPs located at 16 unique loci were identified, and each explained 6.3 to 18.6% of the phenotypic variations. Among these, six loci were overlapped with previously reported genes or quantitative trait loci, whereas the other 10 were novel. Besides, nine loci were detected for two or more traits, indicating that the salt-tolerance genetic architecture is complex. Furthermore, five candidate genes were identified for salt tolerance-related traits, including kinase family protein, E3 ubiquitin-protein ligase-like protein, and transmembrane protein. SNPs identified in this study and the accessions with more favorable alleles could further enhance salt tolerance in wheat breeding. Our results are useful for uncovering the genetic mechanism of salt tolerance in wheat at seeding stage.

Genome-wide association study of multiple traits linked to heat tolerance in emmer-derived hexaploid wheat genotypes

Heat stress tolerance in plants is a complex trait controlled by multiple genes of minor effect which are influenced by the environment and this makes breeding and selection complicated. Emmer wheat (Triticum dicoccon Schrank) carries valuable diversity that can be used to improve the heat tolerance of modern bread wheat. A diverse set of emmer-based genotypes was developed by crossing emmer wheat with hexaploid wheat. These materials, along with their hexaploid recurrent parents and commercial cultivars, were evaluated at optimum (E1) and heat stressed (E2) sowing times in the field for three consecutive years (2014-2016). The material was genotyped using the Infinium iSelect SNP 90K SNP Assay. The phenotypic data were combined across years within each sowing time and best linear unbiased estimators calculated for each genotype in each environment. These estimates were used for GWAS analysis. Significant phenotypic and genotypic variation was observed for all traits. A total of 125 and 142 marker-trait associations (MTAs) were identified in E1 and E2, respectively. The highest number of MTAs were observed on the A genome (106), followed by the B (105) and D (56) genomes. MTAs with pleiotropic effects within and across the environments were observed. Many of the MTAs found were reported previously for various traits, and a few significant MTAs under heat stress were new and linked to emmer genome. Genomic regions identified on chromosomes 2B and 3A had a significant positive impact on grain yield under stress with a 7% allelic effect. Genomic regions on chromosomes 1A and 4B contributed 11% and 9% of the variation for thousand kernel weight (TKW) under heat stress respectively. Following fine mapping, these regions could be used for marker-assisted selection to improve heat tolerance in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01222-3.

Drivers of linkage disequilibrium across a species' geographic range

While linkage disequilibrium (LD) is an important parameter in genetics and evolutionary biology, the drivers of LD remain elusive. Using whole-genome sequences from across a species’ range, we assessed the impact of demographic history and mating system on LD. Both range expansion and a shift from outcrossing to selfing in North American Arabidopsis lyrata were associated with increased average genome-wide LD. Our results indicate that range expansion increases short-distance LD at the farthest range edges by about the same amount as a shift to selfing. However, the extent over which LD in genic regions unfolds was shorter for range expansion compared to selfing. Linkage among putatively neutral variants and between neutral and deleterious variants increased to a similar degree with range expansion, providing support that genome-wide LD was positively associated with mutational load. As a consequence, LD combined with mutational load may decelerate range expansions and set range limits. Finally, a small number of genes were identified as LD outliers, suggesting that they experience selection by either of the two demographic processes. These included genes involved in flowering and photoperiod for range expansion, and the self-incompatibility locus for mating system.

Nearby genomic variants are often co-inherited because of limited recombination. The extent of non-random association of alleles at different loci is called linkage disequilibrium (LD) and is commonly used in genomic analyses, for example to detect regions under selection or to determine effective population size. Here we reversed testing and addressed how demographic history may affect LD within a species. Using genomic data from more than a thousand individuals of North American Arabidopsis lyrata from across the entire species’ range, we quantified the effect of postglacial range expansion and a shift in mating system from outcrossing to selfing on LD. We show that both factors lead to increased LD, and that the maximal effect of range expansion is comparable with a shift in mating system to selfing. Heightened LD involves deleterious mutations, and therefore, LD can also serve as an indicator of mutation accumulation. Furthermore, we provide evidence that some genes experienced stronger increases in LD possibly due to selection associated with the two demographic changes. Our results provide a novel and broad view on the evolutionary factors shaping LD that may also apply to the very many species that underwent postglacial range expansion.

Using natural variation to understand plant responses to iron availability

Iron bioavailability varies dramatically between soil types across the globe. This has given rise to high levels of natural variation in plant iron responses, allowing members of even a single species to thrive across a wide range of soil types. In recent years we have seen the use of genome-wide association analysis to identify natural variants underlying plant responses to changes in iron availability in both Arabidopsis and important crop species. These studies have provided insights into which genes have been important in shaping local adaptation to iron availability in different plant species and have allowed the discovery of novel regulators and mechanisms, not previously identified using mutagenesis approaches. Furthermore, these studies have allowed the identification of markers that can be used to accelerate breeding of future elite varieties with increased resilience to iron stress and improved nutritional quality. The studies highlighted here show that, in addition to studying plant responses to iron alone, it is important to consider these responses within the context of plant nutrition more broadly and to also consider iron regulation in relation to additional traits of agronomic importance such as yield and disease resistance.

Genome-Wide Association Study Reveals Marker-Trait Associations for Early Vegetative Stage Salinity Tolerance in Rice

Rice germplasm is a rich resource for discovering genes associated with salt tolerance. In the current study, a set of 96 accessions were evaluated for seedling stage salinity tolerance and its component traits. Significant phenotypic variation was observed among the genotypes for all the measured traits and eleven accessions with high level of salt tolerance at seedling stage were identified. The germplasm set comprised of three sub-populations and genome-wide association study (GWAS) identified a total of 23 marker-trait associations (MTAs) for traits studied. These MTAs were located on rice chromosomes 1, 2, 5, 6, 7, 9, and 12 and explained the trait phenotypic variances ranging from 13.98 to 29.88 %. Twenty-one MTAs identified in this study were located either in or near the previously reported quantitative trait loci (QTLs), while two MTAs namely, qSDW2.1 and qSNC5 were novel. A total of 18 and 13 putative annotated candidate genes were identified in a genomic region spanning ~200 kb around the MTAs qSDW2.1 and qSNC5, respectively. Some of the important genes underlying the novel MTAs were OsFBA1,OsFBL7, and mTERF which are known to be associated with salinity tolerance in crops. These MTAs pave way for combining salinity tolerance with high yield in rice genotypes through molecular breeding.

Association analysis and molecular tagging of phytochemicals in the endangered medicinal plant licorice (Glycyrrhiza glabra L.)

Licorice (Glycyrrhiza glabra L.) is a medicinal plant species valued in many countries in Asia and Europe for its phytochemical characteristics. Licorice biodiversity is becoming threatened nowadays in Iran due to increasing demand and a drastic decline of its natural habitats. Therefore, licorice domestication would be necessary in the near future, and molecular breeding would help to introduce genotypes suitable for cultivation. The present study was carried out with 170 individual licorice plants sampled in the wild in 59 localizations in 21 provinces of Iran. The association of 436 polymorphic AFLP markers, produced by 15 primer combinations (EcoRI/MseI), with six phenotypic phytochemical traits was studied. The AMOVA analysis show gene diversity among and within localizations. The population structure analysis identified two main sub-populations with significant genetic variation. Significant associations were identified between three markers (E3/M40-4, E34/M4-12 and E12/M31-15) and glycyrrhizin concentration, and between four markers (E11/M34-12, E11/M34-15, E9/M7-29, and E9/M7-30) and phenolic compounds contents. Markers detected can be useful in the domestication of licorice as well as in breeding programs. Licorice sampled in four localizations (KBA1, KBA2, SKh2 and Fa1) were found to be superior in terms of glycyrrhizin and antioxidants content, and therefore they can be considered as elite genotypes which could be included in the domestication process.

Pairwise linkage disequilibrium estimation for polyploids

Many tasks in statistical genetics involve pairwise estimation of linkage disequilibrium (LD). The study of LD in diploids is mature. However, in polyploids, the field lacks a comprehensive characterization of LD. Polyploids also exhibit greater levels of genotype uncertainty than diploids, yet no methods currently exist to estimate LD in polyploids in the presence of such genotype uncertainty. Furthermore, most LD estimation methods do not quantify the level of uncertainty in their LD estimates. Our study contains three major contributions. (i) We characterize haplotypic and composite measures of LD in polyploids. These composite measures of LD turn out to be functions of common statistical measures of association. (ii) We derive procedures to estimate haplotypic and composite LD in polyploids in the presence of genotype uncertainty. We do this by estimating LD directly from genotype likelihoods, which may be obtained from many genotyping platforms. (iii) We derive standard errors of all LD estimators that we discuss. We validate our methods on both real and simulated data. Our methods are implemented in the R package ldsep, available on the Comprehensive R Archive Network https://cran.r-project.org/package=ldsep.

Association analysis of physiological traits in spring barley (Hordeum vulgare L.) under water-deficit conditions

In the present study, 148 commercial barley cultivars were assessed by 14 AFLP primer combinations and 32 SSRs primer pairs. Population structure, linkage disequilibrium, and genomic regions associated with physiological traits under drought stress were investigated. The phenotypic results showed a high level of diversity between studied cultivars. The studied barley cultivars were divided into two subgroups. Linkage disequilibrium analysis revealed that r 2 values among all possible marker pairs have an average value of 0.0178. The mixed linear model procedure showed that totally, 207 loci had a significant association with investigated traits. 120 QTLs out of 207 were detected for traits under normal conditions, and 90 QTLs were detected for traits under drought stress conditions. Identified QTLs after validation and transferring to SCAR markers in the case of AFLPs can be used to develop MAS strategies for barley breeding programs. Some common markers were identified for a particular trait or some traits across normal and drought stress conditions. These markers show low interaction with environmental conditions (stable markers); therefore, selection by them for a trait under normal conditions will improve the trait value under stress conditions, too.

Innovative approach to identify multigenomic and environmental interactions associated with birth defects in family-based hybrid designs

Genes, including those with transgenerational effects, work in concert with behavioral, environmental, and social factors via complex biological networks to determine human health. Understanding complex relationships between causal factors underlying human health is an essential step towards deciphering biological mechanisms. We propose a new analytical framework to investigate the interactions between maternal and offspring genetic variants or their surrogate single nucleotide polymorphisms (SNPs) and environmental factors using family-based hybrid study design. The proposed approach can analyze diverse genetic and environmental factors and accommodate samples from a variety of family units, including case/control-parental triads, and case/control-parental dyads, while minimizing potential bias introduced by population admixture. Comprehensive simulations demonstrated that our innovative approach outperformed the log-linear approach, the best available method for case-control family data. The proposed approach had greater statistical power and was capable to unbiasedly estimate the maternal and child genetic effects and the effects of environmental factors, while controlling the Type I error rate against population stratification. Using our newly developed approach, we analyzed the associations between maternal and fetal SNPs and obstructive and conotruncal heart defects, with adjustment for demographic and lifestyle factors and dietary supplements. Fourteen and 11 fetal SNPs were associated with obstructive and conotruncal heart defects, respectively. Twenty-seven and 17 maternal SNPs were associated with obstructive and conotruncal heart defects, respectively. In addition, maternal body mass index was a significant risk factor for obstructive defects. The proposed approach is a powerful tool for interrogating the etiological mechanism underlying complex traits.

Genetic characterization of an almond germplasm collection and volatilome profiling of raw and roasted kernels

Almond is appreciated for its nutraceutical value and for the aromatic profile of the kernels. In this work, an almond collection composed of 96 Sicilian accessions complemented with 10 widely cultivated cultivars was phenotyped for the production of volatile organic compounds using a proton-transfer time-of-flight mass spectrometer and genotyped using the Illumina Infinium^®18 K Peach SNP array. The profiling of the aroma was carried out on fresh and roasted kernels enabling the detection of 150 mass peaks. Sixty eight, for the most related with sulfur compounds, furan containing compounds, and aldehydes formed by Strecker degradation, significantly increased during roasting, while the concentration of fifty-four mass peaks, for the most belonging to alcohols and terpenes, significantly decreased. Four hundred and seventy-one robust SNPs were selected and employed for population genetic studies. Structure analysis detected three subpopulations with the Sicilian accessions characterized by a different genetic stratification compared to those collected in Apulia (South Italy) and the International cultivars. The linkage-disequilibrium (LD) decay across the genome was equal to r² = 0.083. Furthermore, a high level of collinearity (r² = 0.96) between almond and peach was registered confirming the high synteny between the two genomes. A preliminary application of a genome-wide association analysis allowed the detection of significant marker-trait associations for 31 fresh and 33 roasted almond mass peaks respectively. An accurate genetic and phenotypic characterization of novel germplasm can represent a valuable tool for the set-up of marker-assisted selection of novel cultivars with an enhanced aromatic profile.

Multiparental Population in Crops: Methods of Development and Dissection of Genetic Traits

Multiparental populations are located midway between association mapping that relies on germplasm collections and classic linkage analysis, based upon biparental populations. They provide several key advantages such as the possibility to include a higher number of alleles and increased level of recombination with respect to biparental populations, and more equilibrated allelic frequencies than association mapping panels. Moreover, in these populations new allele's combinations arise from recombination that may reveal transgressive phenotypes and make them a useful pre-breeding material. Here we describe the strategies for working with multiparental populations, focusing on nested association mapping populations (NAM) and multiparent advanced generation intercross populations (MAGIC). We provide details from the selection of founders, population development, and characterization to the statistical methods for genetic mapping and quantitative trait detection.

Understanding the classics: the unifying concepts of transgressive segregation, inbreeding depression and heterosis and their central relevance for crop breeding

Transgressive segregation and heterosis are the reasons that plant breeding works. Molecular explanations for both phenomena have been suggested and play a contributing role. However, it is often overlooked by molecular genetic researchers that transgressive segregation and heterosis are most simply explained by dispersion of favorable alleles. Therefore, advances in molecular biology will deliver the most impact on plant breeding when integrated with sources of heritable trait variation - and this will be best achieved within a quantitative genetics framework. An example of the power of quantitative approaches is the implementation of genomic selection, which has recently revolutionized animal breeding. Genomic selection is now being applied to both hybrid and inbred crops and is likely to be the major source of improvement in plant breeding practice over the next decade. Breeders' ability to efficiently apply genomic selection methodologies is due to recent technology advances in genotyping and sequencing. Furthermore, targeted integration of additional molecular data (such as gene expression, gene copy number and methylation status) into genomic prediction models may increase their performance. In this review, we discuss and contextualize a suite of established quantitative genetics themes relating to hybrid vigour, transgressive segregation and their central relevance to plant breeding, with the aim of informing crop researchers outside of the quantitative genetics discipline of their relevance and importance to crop improvement. Better understanding between molecular and quantitative disciplines will increase the potential for further improvements in plant breeding methodologies and so help underpin future food security.