Crop Breeding Chips and Genotyping Platforms: Progress, Challenges, and Perspectives

Semi-Thermal Asymmetric Reverse PCR (STARP) Genotyping

PCR-based individual Single nucleotide polymorphism (SNP) genotyping methods are preferred due to their flexibility, high-throughput, and improved accuracy. Semi-thermal asymmetric reverse PCR (STARP) is one of the SNP genotyping methods developed to reduce operational cost with improved platform compatibility. STARP is a unique method which can be used either as a gel-free SNP genotyping by detection of fluorescent signals or polyacrylamide gel-based size separation. SNP assay designing using sequence information and detection methods of STARP are described in detail.

Transcriptomic Approach for Global Distribution of SNP/Indel and Plant Genotyping

Single Nucleotide Polymorphisms (SNPs) are the most common structural variants found in any genome. They have been used for different genetic studies, from the understanding of genetic structure of populations to the development of breeding selection markers. In this chapter we present the use of transcriptomic data obtained from contrasting phenotypes for a target trait, in searching of SNPs and insertions/deletions (InDels). This approach has the advantage that the identified markers are in or close to differentially expressed genes, and so they have higher chances to tag the genes underlying the phenotypic expression of a particular trait.

3D-GBS: a universal genotyping-by-sequencing approach for genomic selection and other high-throughput low-cost applications in species with small to medium-sized genomes

Despite the increased efficiency of sequencing technologies and the development of reduced-representation sequencing (RRS) approaches allowing high-throughput sequencing (HTS) of multiplexed samples, the per-sample genotyping cost remains the most limiting factor in the context of large-scale studies. For example, in the context of genomic selection (GS), breeders need genome-wide markers to predict the breeding value of large cohorts of progenies, requiring the genotyping of thousands candidates. Here, we introduce 3D-GBS, an optimized GBS procedure, to provide an ultra-high-throughput and ultra-low-cost genotyping solution for species with small to medium-sized genome and illustrate its use in soybean. Using a combination of three restriction enzymes (PstI/NsiI/MspI), the portion of the genome that is captured was reduced fourfold (compared to a "standard" ApeKI-based protocol) while reducing the number of markers by only 40%. By better focusing the sequencing effort on limited set of restriction fragments, fourfold more samples can be genotyped at the same minimal depth of coverage. This GBS protocol also resulted in a lower proportion of missing data and provided a more uniform distribution of SNPs across the genome. Moreover, we investigated the optimal number of reads per sample needed to obtain an adequate number of markers for GS and QTL mapping (500-1000 markers per biparental cross). This optimization allows sequencing costs to be decreased by ~ 92% and ~ 86% for GS and QTL mapping studies, respectively, compared to previously published work. Overall, 3D-GBS represents a unique and affordable solution for applications requiring extremely high-throughput genotyping where cost remains the most limiting factor.

Development of breeder chip for gene detection and molecular-assisted selection by target sequencing in wheat

Wheat is an essential food crop and its high and stable yield is suffering from great challenges due to the limitations of current breeding technology and various stresses. Accelerating molecularly assisted stress-resistance breeding is critical. Through a meta-analysis of published loci in wheat over the last two decades, we selected 60 loci with main breeding objectives, high heritability, and reliable genotyping, such as stress resistance, yield, plant height, and resistance to spike germination. Then, using genotyping by target sequencing (GBTS) technology, we developed a liquid phase chip based on 101 functional or closely linked markers. The genotyping of 42 loci was confirmed in an extensive collection of Chinese wheat cultivars, indicating that the chip can be used in molecular-assisted selection (MAS) for target breeding goals. Besides, we can perform the preliminary parentage analysis with the genotype data. The most significant contribution of this work lies in translating a large number of molecular markers into a viable chip and providing reliable genotypes. Breeders can quickly screen germplasm resources, parental breeding materials, and intermediate materials for the presence of excellent allelic variants using the genotyping data by this chip, which is high throughput, convenient, reliable, and cost-efficient.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01359-3.

Mining legume germplasm for genetic gains: An Indian perspective

Legumes play a significant role in food and nutritional security and contribute to environmental sustainability. Although legumes are highly beneficial crops, it has not yet been possible to enhance their yield and production to a satisfactory level. Amid a rising population and low yield levels, per capita average legume consumption in India has fallen by 71% over the last 50 years, and this has led to protein-related malnutrition in a large segment of the Indian population, especially women and children. Several factors have hindered attempts to achieve yield enhancement in grain legumes, including biotic and abiotic pressures, a lack of good ideotypes, less amenability to mechanization, poorer responsiveness to fertilizer input, and a poor genetic base. Therefore, there is a need to mine the approximately 0.4 million ex situ collections of legumes that are being conserved in gene banks globally for identification of ideal donors for various traits. The Indian National Gene Bank conserves over 63,000 accessions of legumes belonging to 61 species. Recent initiatives have been undertaken in consortia mode with the aim of unlocking the genetic potential of ex situ collections and conducting large-scale germplasm characterization and evaluation analyses. We assume that large-scale phenotyping integrated with omics-based science will aid the identification of target traits and their use to enhance genetic gains. Additionally, in cases where the genetic base of major legumes is narrow, wild relatives have been evaluated, and these are being exploited through pre-breeding. Thus far, >200 accessions of various legumes have been registered as unique donors for various traits of interest.

Wheat Omics: Advancements and Opportunities

Plant omics, which includes genomics, transcriptomics, metabolomics and proteomics, has played a remarkable role in the discovery of new genes and biomolecules that can be deployed for crop improvement. In wheat, great insights have been gleaned from the utilization of diverse omics approaches for both qualitative and quantitative traits. Especially, a combination of omics approaches has led to significant advances in gene discovery and pathway investigations and in deciphering the essential components of stress responses and yields. Recently, a Wheat Omics database has been developed for wheat which could be used by scientists for further accelerating functional genomics studies. In this review, we have discussed various omics technologies and platforms that have been used in wheat to enhance the understanding of the stress biology of the crop and the molecular mechanisms underlying stress tolerance.

The Pyramiding of Elite Allelic Genes Related to Grain Number Increases Grain Number per Panicle Using the Recombinant Lines Derived from Indica-japonica Cross in Rice

Indica(xian)-japonica(geng) hybrid rice has many heterosis traits that can improve rice yield. However, the traditional hybrid technology will struggle to meet future needs for the development of higher-yield rice. Available genomics resources can be used to efficiently understand the gene-trait association trait for rice breeding. Based on the previously constructed high-density genetic map of 272 high-generation recombinant inbred lines (RILs) originating from the cross of Luohui 9 (indica, as female) and RPY geng (japonica, as male) and high-quality genomes of parents, here, we further explore the genetic basis for an important complex trait: possible causes of grain number per panicle (GNPP). A total of 20 genes related to grains number per panicle (GNPP) with the differences of protein amino acid between LH9 and RPY were used to analyze genotype combinations, and PCA results showed a combination of PLY1, LAX1, DTH8 and OSH1 from the RPY geng with PYL4, SP1, DST and GNP1 from Luohui 9 increases GNPP. In addition, we also found that the combination of LAX1-T2 and GNP1-T3 had the most significant increase in GNPP. Notably, Molecular Breeding Knowledgebase (MBK) showed a few aggregated rice cultivars, LAX1-T2 and GNP1-T3, which may be a result of the natural geographic isolation between the two gene haplotypes. Therefore, we speculate that the pyramiding of japonica-type LAX-T2 with indica-type GNP1-T3 via hybridization can significantly improve rice yield by increasing GNPP.

Development and Evaluation of an AxiomTM 60K SNP Array for Almond (Prunus dulcis)

A high-density single nucleotide polymorphism (SNP) array is essential to enable faster progress in plant breeding for new cultivar development. In this regard, we have developed an Axiom 60K almond SNP array by resequencing 81 almond accessions. For the validation of the array, a set of 210 accessions were genotyped and 82.8% of the SNPs were classified in the best recommended SNPs. The rate of missing data was between 0.4% and 2.7% for the almond accessions and less than 15.5% for the few peach and wild accessions, suggesting that this array can be used for peach and interspecific peach × almond genetic studies. The values of the two SNPs linked to the RMja (nematode resistance) and SK (bitterness) genes were consistent. We also genotyped 49 hybrids from an almond F2 progeny and could build a genetic map with a set of 1159 SNPs. Error rates, less than 1%, were evaluated by comparing replicates and by detection of departures from Mendelian inheritance in the F2 progeny. This almond array is commercially available and should be a cost-effective genotyping tool useful in the search for new genes and quantitative traits loci (QTL) involved in the control of agronomic traits.

BREEDIT: a multiplex genome editing strategy to improve complex quantitative traits in maize

Ensuring food security for an ever-growing global population while adapting to climate change is the main challenge for agriculture in the 21st century. Although new technologies are being applied to tackle this problem, we are approaching a plateau in crop improvement using conventional breeding. Recent advances in CRISPR/Cas9-mediated gene engineering have paved the way to accelerate plant breeding to meet this increasing demand. However, many traits are governed by multiple small-effect genes operating in complex interactive networks. Here, we present the gene discovery pipeline BREEDIT, which combines multiplex genome editing of whole gene families with crossing schemes to improve complex traits such as yield and drought tolerance. We induced gene knockouts in 48 growth-related genes into maize (Zea mays) using CRISPR/Cas9 and generated a collection of over 1,000 gene-edited plants. The edited populations displayed (on average) 5%-10% increases in leaf length and up to 20% increases in leaf width compared with the controls. For each gene family, edits in subsets of genes could be associated with enhanced traits, allowing us to reduce the gene space to be considered for trait improvement. BREEDIT could be rapidly applied to generate a diverse collection of mutants to identify promising gene modifications for later use in breeding programs.

Genome-wide association study for grain zinc concentration in bread wheat (Triticum aestivum L.)

Introduction

Zinc (Zn) deficiency causes serious diseases in people who rely on cereals as their main food source. However, the grain zinc concentration (GZnC) in wheat is low. Biofortification is a sustainable strategy for reducing human Zn deficiency.

Methods

In this study, we constructed a population of 382 wheat accessions and determined their GZnC in three field environments. Phenotype data was used for a genome-wide association study (GWAS) using a 660K single nucleotide polymorphism (SNP) array, and haplotype analysis identified an important candidate gene for GZnC.

Results

We found that GZnC of the wheat accessions showed an increasing trend with their released years, indicating that the dominant allele of GZnC was not lost during the breeding process. Nine stable quantitative trait loci (QTLs) for GZnC were identified on chromosomes 3A, 4A, 5B, 6D, and 7A. And an important candidate gene for GZnC, namely, TraesCS6D01G234600, and GZnC between the haplotypes of this gene showed, significant difference (P ≤ 0.05) in three environments.

Discussion

A novel QTL was first identified on chromosome 6D, this finding enriches our understanding of the genetic basis of GZnC in wheat. This study provides new insights into valuable markers and candidate genes for wheat biofortification to improve GZnC.

Polymorphisms in alternative oxidase genes from ecotypes of Arabidopsis and rice revealed an environment-induced linkage to altitude and rainfall

We investigated SNPs in alternative oxidase (AOX) genes and their connection to ecotype origins (climate, altitude, and rainfall) by using genomic data sets of Arabidopsis and rice populations from 1190 and 90 ecotypes, respectively. Parameters were defined to detect non-synonymous SNPs in the AOX ORF, which revealed amino acid (AA) changes in AOX1c, AOX1d, and AOX2 from Arabidopsis and AOX1c from rice in comparison to AOX references from Columbia-0 and Japonica ecotypes, respectively. Among these AA changes, Arabidopsis AOX1c_A161E&G165R and AOX1c_R242S revealed a link to high rainfall and high altitude, respectively, while all other changes in Arabidopsis and rice AOX was connected to high altitude and rainfall. Comparative 3D modeling showed that all mutant AOX presented structural differences in relation to the respective references. Molecular docking analysis uncovered lower binding affinity values between AOX and the substrate ubiquinol for most of the identified structures compared to their reference, indicating better enzyme-substrate binding affinities. Thus, our in silico data suggest that the majority of the AA changes found in the available ecotypes will confer better enzyme-subtract interactions and thus indicate environment-related, more efficient AOX activity.

Molecular mapping of drought-responsive QTLs during the reproductive stage of rice using a GBS (genotyping-by-sequencing) based SNP linkage map

BACKGROUND

In rice, drought stress at reproductive stage drastically reduces yield, which in turn hampers farmer's efforts towards crop production. The majority of the rice varieties have resistance genes against several abiotic and biotic stresses. Therefore, the traditional landraces were studied to identify QTLs/candidate genes associated with drought tolerance.

METHODS AND RESULTS

A high-density SNP-based genetic map was constructed using a Genotyping-by-sequencing (GBS) approach. The recombinant inbred lines (RILs) derived from crossing 'Banglami × Ranjit' were used for QTL analysis. A total map length of 1306.424 cM was constructed, which had an average inter-marker distance of 0.281 cM. The phenotypic evaluation of F₆ and F₇ RILs were performed under drought stress and control conditions. A total of 42 QTLs were identified under drought stress and control conditions for yield component traits explaining 1.95-13.36% of the total phenotypic variance (PVE). Among these, 19 QTLs were identified under drought stress conditions, whereas 23 QTLs were located under control conditions. A total of 4 QTLs explained a PVE ≥ 10% which are considered as the major QTLs. Moreover, bioinformatics analysis revealed the presence of 6 candidate genes, which showed differential expression under drought and control conditions.

CONCLUSION

These QTLs/genes may be deployed for marker-assisted pyramiding to improve drought tolerance in the existing rice varieties.

Development and Complex Application of Methods for the Identification of Mutations in the FAD3A and FAD3B Genes Resulting in the Reduced Content of Linolenic Acid in Flax Oil

Flax is grown worldwide for seed and fiber production. Linseed varieties differ in their oil composition and are used in pharmaceutical, food, feed, and industrial production. The field of application primarily depends on the content of linolenic (LIN) and linoleic (LIO) fatty acids. Inactivating mutations in the FAD3A and FAD3B genes lead to a decrease in the LIN content and an increase in the LIO content. For the identification of the three most common low-LIN mutations in flax varieties (G-to-A in exon 1 of FAD3A substituting tryptophan with a stop codon, C-to-T in exon 5 of FAD3A leading to arginine to a stop codon substitution, and C-to-T in exon 2 of FAD3B resulting in histidine to tyrosine substitution), three approaches were proposed: (1) targeted deep sequencing, (2) high resolution melting (HRM) analysis, (3) cleaved amplified polymorphic sequences (CAPS) markers. They were tested on more than a thousand flax samples of various types and showed promising results. The proposed approaches can be used in marker-assisted selection to choose parent pairs for crosses, separate heterogeneous varieties into biotypes, and select genotypes with desired homozygous alleles of the FAD3A and FAD3B genes at the early stages of breeding for the effective development of varieties with a particular LIN and LIO content, as well as in basic studies of the molecular mechanisms of fatty acid synthesis in flax seeds to select genotypes adequate to the tasks.

Conditional QTL mapping for seed germination and seedling traits under salt stress and candidate gene prediction in wheat

Breeding new wheat varieties with salt resistance is one of the best ways to solve a constraint on the sustainability and expansion of wheat cultivation. Therefore, understanding the molecular components or genes related to salt tolerance must contribute to the cultivation of salt-tolerant varieties. The present study used a recombinant inbred line (RIL) population to genetically dissect the effects of different salt stress concentrations on wheat seed germination and seedling traits using two quantitative trait locus (QTL) mapping methods. A total of 31 unconditional and 11 conditional QTLs for salt tolerance were identified on 11 chromosomes explaining phenotypic variation (PVE) ranging from 2.01 to 65.76%. Of these, 15 major QTLs were found accounting for more than 10% PVE. QTL clusters were detected on chromosomes 2A and 3B in the marker intervals 'wPt-8328 and wPt-2087' and 'wPt-666008 and wPt-3620', respectively, involving more than one salt tolerance trait. QRdw3B and QSfw3B.2 were most consistent in two or more salt stress treatments. 16 candidate genes associated with salt tolerance were predicted in wheat. These results could be useful to improve salt tolerance by marker-assisted selection (MAS) and shed new light on understanding the genetic basis of salt tolerance in wheat.

Identification of arsenic-tolerant varieties and candidate genes of tolerance in spring wheat (Triticum aestivum L.)

Despite the growing concerns about arsenic toxicity, information on tolerance and responsible genetic factors in wheat remains elusive. To address that, the present study aimed to screen the wheat varieties against arsenic based on growth parameters, yield, grain accumulation, and associated genes. A total of 110 wheat varieties were grown in arsenic-contaminated regions to record physio-morphological traits. The wheat 90K Infinium iSelect SNP array was used for the genome-wide association model to identify genomic regions. Wheat varieties such as Punjab-81, AARI-11, and Daman showed arsenic concentrations >45 μg/kg in similar conditions as well as the impact on grain yield, chlorophyll, Thousand Kernel Weight, and plant height. Contrastingly, varieties like Kohistan-97, As-2002, Barani-70, and Pari-73 showed grain concentrations <5 μg/kg grown under highly contaminated conditions. Three significant loci associated with arsenic accumulation in grain were identified on chromosomes 6A (qASG1-6A) and 6B (qASG3-6B and qASG4-6B). Annotation at these loci identified 39 wheat genes among which several were important for growth and tolerance against stress. The candidate gene (TraesCS6B02G429400) responsible for Glutathione-S-transferase was identified in the present study and must be investigated further using a transcriptomic approach. The present study provided background information for breeding prospects to improve wheat yield and tolerance against arsenic.

Development of the maize 5.5K loci panel for genomic prediction through genotyping by target sequencing

Genotyping platforms are important for genetic research and molecular breeding. In this study, a low-density genotyping platform containing 5.5K SNP markers was successfully developed in maize using genotyping by target sequencing (GBTS) technology with capture-in-solution. Two maize populations (Pop1 and Pop2) were used to validate the GBTS panel for genetic and molecular breeding studies. Pop1 comprised 942 hybrids derived from 250 inbred lines and four testers, and Pop2 contained 540 hybrids which were generated from 123 new-developed inbred lines and eight testers. The genetic analyses showed that the average polymorphic information content and genetic diversity values ranged from 0.27 to 0.38 in both populations using all filtered genotyping data. The mean missing rate was 1.23% across populations. The Structure and UPGMA tree analyses revealed similar genetic divergences (76-89%) in both populations. Genomic prediction analyses showed that the prediction accuracy of reproducing kernel Hilbert space (RKHS) was slightly lower than that of genomic best linear unbiased prediction (GBLUP) and three Bayesian methods for general combining ability of grain yield per plant and three yield-related traits in both populations, whereas RKHS with additive effects showed superior advantages over the other four methods in Pop1. In Pop1, the GBLUP and three Bayesian methods with additive-dominance model improved the prediction accuracies by 4.89-134.52% for the four traits in comparison to the additive model. In Pop2, the inclusion of dominance did not improve the accuracy in most cases. In general, low accuracies (0.33-0.43) were achieved for general combing ability of the four traits in Pop1, whereas moderate-to-high accuracies (0.52-0.65) were observed in Pop2. For hybrid performance prediction, the accuracies were moderate to high (0.51-0.75) for the four traits in both populations using the additive-dominance model. This study suggests a reliable genotyping platform that can be implemented in genomic selection-assisted breeding to accelerate maize new cultivar development and improvement.

Unlocking the hidden variation from wild repository for accelerating genetic gain in legumes

The fluctuating climates, rising human population, and deteriorating arable lands necessitate sustainable crops to fulfil global food requirements. In the countryside, legumes with intriguing but enigmatic nitrogen-fixing abilities and thriving in harsh climatic conditions promise future food security. However, breaking the yield plateau and achieving higher genetic gain are the unsolved problems of legume improvement. Present study gives emphasis on 15 important legume crops, i.e., chickpea, pigeonpea, soybean, groundnut, lentil, common bean, faba bean, cowpea, lupin, pea, green gram, back gram, horse gram, moth bean, rice bean, and some forage legumes. We have given an overview of the world and India's area, production, and productivity trends for all legume crops from 1961 to 2020. Our review article investigates the importance of gene pools and wild relatives in broadening the genetic base of legumes through pre-breeding and alien gene introgression. We have also discussed the importance of integrating genomics, phenomics, speed breeding, genetic engineering and genome editing tools in legume improvement programmes. Overall, legume breeding may undergo a paradigm shift once genomics and conventional breeding are integrated in the near future.

Identification of quantitative trait loci of agronomic traits in bread wheat using a Pamyati Azieva × Paragon mapping population harvested in three regions of Kazakhstan

Background

Although genome-wide association studies (GWAS) are an increasingly informative tool in the mining of new quantitative trait loci (QTLs), a classical biparental mapping approach is still a powerful, widely used method to search the unique genetic factors associated with important agronomic traits in bread wheat.

Methods

In this study, a newly constructed mapping population of Pamyati Azieva (Russian Federation) × Paragon (UK), consisting of 94 recombinant inbred lines (RILs), was tested in three different regions of Kazakhstan with the purpose of QTL identification for key agronomic traits. The RILs were tested in 11 environments of two northern breeding stations (Petropavlovsk, North Kazakhstan region, and Shortandy, Aqmola region) and one southeastern station (Almalybak, Almaty region). The following eight agronomic traits were studied: heading days, seed maturation days, plant height, spike length, number of productive spikes, number of kernels per spike, thousand kernel weight, and yield per square meter. The 94 RILs of the PAxP cross were genotyped using Illumina's iSelect 20K single nucleotide polymorphism (SNP) array and resulted in the identification of 4595 polymorphic SNP markers.

Results

The application of the QTL Cartographer statistical package allowed the identification of 53 stable QTLs for the studied traits. A survey of published studies related to common wheat QTL identification suggested that 28 of those 53 QTLs were presumably novel genetic factors. The SNP markers for the identified QTLs of the analyzed agronomic traits of common wheat can be efficiently applied in ongoing breeding activities in the wheat breeding community using a marker-assisted selection approach.

Development of an inclusive 580K SNP array and its application for genomic selection and genome-wide association studies in rice

Rice is a globally cultivated crop and is primarily a staple food source for more than half of the world's population. Various single-nucleotide polymorphism (SNP) arrays have been developed and utilized as standard genotyping methods for rice breeding research. Considering the importance of SNP arrays with more inclusive genetic information for GWAS and genomic selection, we integrated SNPs from eight different data resources: resequencing data from the Korean World Rice Collection (KRICE) of 475 accessions, 3,000 rice genome project (3 K-RGP) data, 700 K high-density rice array, Affymetrix 44 K SNP array, QTARO, Reactome, and plastid and GMO information. The collected SNPs were filtered and selected based on the breeder's interest, covering all key traits or research areas to develop an integrated array system representing inclusive genomic polymorphisms. A total of 581,006 high-quality SNPs were synthesized with an average distance of 200 bp between adjacent SNPs, generating a 580 K Axiom Rice Genotyping Chip (580 K _ KNU chip). Further validation of this array on 4,720 genotypes revealed robust and highly efficient genotyping. This has also been demonstrated in genome-wide association studies (GWAS) and genomic selection (GS) of three traits: clum length, heading date, and panicle length. Several SNPs significantly associated with cut-off, -log10 p-value >7.0, were detected in GWAS, and the GS predictabilities for the three traits were more than 0.5, in both rrBLUP and convolutional neural network (CNN) models. The Axiom 580 K Genotyping array will provide a cost-effective genotyping platform and accelerate rice GWAS and GS studies.

A high-throughput skim-sequencing approach for genotyping, dosage estimation and identifying translocations

The development of next-generation sequencing (NGS) enabled a shift from array-based genotyping to directly sequencing genomic libraries for high-throughput genotyping. Even though whole-genome sequencing was initially too costly for routine analysis in large populations such as breeding or genetic studies, continued advancements in genome sequencing and bioinformatics have provided the opportunity to capitalize on whole-genome information. As new sequencing platforms can routinely provide high-quality sequencing data for sufficient genome coverage to genotype various breeding populations, a limitation comes in the time and cost of library construction when multiplexing a large number of samples. Here we describe a high-throughput whole-genome skim-sequencing (skim-seq) approach that can be utilized for a broad range of genotyping and genomic characterization. Using optimized low-volume Illumina Nextera chemistry, we developed a skim-seq method and combined up to 960 samples in one multiplex library using dual index barcoding. With the dual-index barcoding, the number of samples for multiplexing can be adjusted depending on the amount of data required, and could be extended to 3,072 samples or more. Panels of doubled haploid wheat lines (Triticum aestivum, CDC Stanley x CDC Landmark), wheat-barley (T. aestivum x Hordeum vulgare) and wheat-wheatgrass (Triticum durum x Thinopyrum intermedium) introgression lines as well as known monosomic wheat stocks were genotyped using the skim-seq approach. Bioinformatics pipelines were developed for various applications where sequencing coverage ranged from 1 × down to 0.01 × per sample. Using reference genomes, we detected chromosome dosage, identified aneuploidy, and karyotyped introgression lines from the skim-seq data. Leveraging the recent advancements in genome sequencing, skim-seq provides an effective and low-cost tool for routine genotyping and genetic analysis, which can track and identify introgressions and genomic regions of interest in genetics research and applied breeding programs.

Genetic dissection of fruit maturity date in apricot (P. armeniaca L.) through a Single Primer Enrichment Technology (SPET) approach

Background

Single primer enrichment technology (SPET) is an emerging and increasingly popular solution for high-throughput targeted genotyping in plants. Although SPET requires a priori identification of polymorphisms for probe design, this technology has potentially higher reproducibility and transferability compared to other reduced representation sequencing (RRS) approaches, also enabling the discovery of closely linked polymorphisms surrounding the target one.

Results

The potential for SPET application in fruit trees was evaluated by developing a 25K target SNPs assay to genotype a panel of apricot accessions and progenies. A total of 32,492 polymorphic sites were genotyped in 128 accessions (including 8,188 accessory non-target SNPs) with extremely low levels of missing data and a significant correlation of allelic frequencies compared to whole-genome sequencing data used for array design. Assay performance was further validated by estimating genotyping errors in two biparental progenies, resulting in an overall 1.8% rate. SPET genotyping data were used to infer population structure and to dissect the architecture of fruit maturity date (MD), a quantitative reproductive phenological trait of great agronomical interest in apricot species. Depending on the year, GWAS revealed loci associated to MD on several chromosomes. The QTLs on chromosomes 1 and 4 (the latter explaining most of the phenotypic variability in the panel) were the most consistent over years and were further confirmed by linkage mapping in two segregating progenies.

Conclusions

Besides the utility for marker assisted selection and for paving the way to in-depth studies to clarify the molecular bases of MD trait variation in apricot, the results provide an overview of the performance and reliability of SPET for fruit tree genetics.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08901-1.

Selection and Validation of 48 KASP Markers for Variety Identification and Breeding Guidance in Conventional and Hybrid Rice (Oryza sativa L.)

BACKGROUND

Breeding of conventional and hybrid rice (Oryza sativa L.) have solved hunger problems and increased farmers' income in the world. Molecular markers have been widely used in marker-assisted breeding and identification of larger numbers of different bred varieties in the past decades. The recently developed SNP markers are applied for more stable and detectable compared with other markers. But the cost of genotyping lots SNPs is high. So, it is essential to select less representative SNPs and inexpensive detecting methods to lower the cost and accelerate variety identification and breeding process. KASP (Kompetitive Allele-Specific PCR) is a flexible method to detect the SNPs, and large number of KASP markers have been widely used in variety identification and breeding. However, the ability of less KASP markers on massive variety identification and breeding remains unknown.

RESULTS

Here, 48 KASP markers were selected from 378 markers to classify and analyze 518 varieties including conventional and hybrid rice. Through analyzing the population structure, the 48 markers could almost represent the 378 markers. In terms of variety identification, the 48 KASP markers had a 100% discrimination rate in 53 conventional indica varieties and 193 hybrid varieties, while they could distinguish 89.1% conventional japonica rice from different breeding institutes. Two more markers added would increase the ratio from 68.38 to 77.94%. Additionally, the 48 markers could be used for classification of subpopulations in the bred variety. Also, 8 markers had almost completely different genotypes between japonica and indica, and 3 markers were found to be very important for japonica hybrid rice. In hybrid varieties, the heterozygosity of chromosomes 3, 6 and 11 was relatively higher than others.

CONCLUSIONS

Our results showed that 48 KASP markers could be used to identify rice varieties, and the panel we tested could provide a database for breeders to identify new breeding lines. Also, the specific markers we found were useful for marker-assisted breeding in rice, including conventional and hybrid.

Genome-wide characterization and sequence polymorphism analyses of cysteine-rich poly comb-like protein in Glycine max

Cysteine-rich poly comb-like protein (CPP) is a member of cysteine-rich transcription factors that regulates plant growth and development. In the present work, we characterized twelve CPP transcription factors encoding genes in soybean (Glycine max). Phylogenetic analyses classified CPP genes into six clades. Sequence logos analyses between G. max and G. soja amino acid residues exhibited high conservation. The presence of growth and stress-related cis-acting elements in the upstream regions of GmCPPs highlight their role in plant development and tolerance against abiotic stress. Ka/Ks levels showed that GmCPPs experienced limited selection pressure with limited functional divergence arising from segmental or whole genome duplication events. By using the PAN-genome of soybean, a single nucleotide polymorphism was identified in GmCPP-6. To perform high throughput genotyping, a kompetitive allele-specific PCR (KASP) marker was developed. Association analyses indicated that GmCPP-6-T allele of GmCPP-6 (in exon region) was associated with higher thousand seed weight under both water regimes (well-water and water-limited). Taken together, these results provide vital information to further decipher the biological functions of CPP genes in soybean molecular breeding.

Pleiotropic effect analysis and marker development for grain zinc and iron concentrations in spring wheat

Wheat (Triticum aestivum L.) is one of the main food crops in the world and a primary source of zinc (Zn) and iron (Fe) in the human body. The genetic mechanisms underlying related traits have been clarified, thereby providing a molecular theoretical foundation for the development of germplasm resources. In this study, a total of 23,536 high-quality DArT markers was used to map quantitative trait loci (QTL) of grain Zn (GZn) and grain Fe (GFe) concentrations in recombinant inbred lines crossed by Avocet/Chilero. A total of 17 QTLs was located on chromosomes 1BL, 2BL, 3BL, 4AL, 4BS, 5AL, 5DL, 6AS, 6BS, 6DS, and 7AS accounting for 0.38-16.62% of the phenotypic variance. QGZn.haust-4AL, QGZn.haust-7AS.1, and QGFe.haust-6BS were detected on chromosomes 4AL, 6BS, and 7AS, accounting for 10.63-16.62% of the phenotypic variance. Four stable QTLs, QGZn.haust-4AL, QGFe.haust-1BL, QGFe.haust-4AL, and QGFe.haust-5DL, were located on chromosomes 1BL, 4AL, and 5DL. Three pleiotropic effects loci for GZn and GFe concentrations were located on chromosomes 1BL, 4AL, and 5DL. Two high-throughput Kompetitive Allele Specific PCR markers were developed by closely linking single-nucleotide polymorphisms on chromosomes 4AL and 5DL, which were validated by a germplasm panel. Therefore, it is the most important that quantitative trait loci and KASP marker for grain zinc and iron concentrations were developed for utilizing in marker-assisted breeding and biofortification of wheat grain in breeding programs.

Genetic networks underlying salinity tolerance in wheat uncovered with genome-wide analyses and selective sweeps

A genetic framework underpinning salinity tolerance at reproductive stage was revealed by genome-wide SNP markers and major adaptability genes in synthetic-derived wheats, and trait-associated loci were used to predict phenotypes. Using wild relatives of crops to identify genes related to improved productivity and resilience to climate extremes is a prioritized area of crop genetic improvement. High salinity is a widespread crop production constraint, and development of salt-tolerant cultivars is a sustainable solution. We evaluated a panel of 294 wheat accessions comprising synthetic-derived wheat lines (SYN-DERs) and modern bread wheat advanced lines under control and high salinity conditions at two locations. The GWAS analysis revealed a quantitative genetic framework of more than 200 loci with minor effect underlying salinity tolerance at reproductive stage. The significant trait-associated SNPs were used to predict phenotypes using a GBLUP model, and the prediction accuracy (r²) ranged between 0.57 and 0.74. The r² values for flag leaf weight, days to flowering, biomass, and number of spikes per plant were all above 0.70, validating the phenotypic effects of the loci discovered in this study. Furthermore, the germplasm sets were compared to identify selection sweeps associated with salt tolerance loci in SYN-DERs. Six loci associated with salinity tolerance were found to be differentially selected in the SYN-DERs (12.4 Mb on chromosome (chr)1B, 7.1 Mb on chr2A, 11.2 Mb on chr2D, 200 Mb on chr3D, 600 Mb on chr6B, and 700.9 Mb on chr7B). A total of 228 reported markers and genes, including 17 well-characterized genes, were uncovered using GWAS and EigenGWAS. A linkage disequilibrium (LD) block on chr5A, including the Vrn-A1 gene at 575 Mb and its homeologs on chr5D, were strongly associated with multiple yield-related traits and flowering time under salinity stress conditions. The diversity panel was screened with more than 68 kompetitive allele-specific PCR (KASP) markers of functional genes in wheat, and the pleiotropic effects of superior alleles of Rht-1, TaGASR-A1, and TaCwi-A1 were revealed under salinity stress. To effectively utilize the extensive genetic information obtained from the GWAS analysis, a genetic interaction network was constructed to reveal correlations among the investigated traits. The genetic network data combined with GWAS, selective sweeps, and the functional gene survey provided a quantitative genetic framework for identifying differentially retained loci associated with salinity tolerance in wheat.

Genome wide association study to detect genetic regions related to isoflavone content in a mutant soybean population derived from radiation breeding

Isoflavones are major secondary metabolites that are exclusively produced by legumes, including soybean. Soy isoflavones play important roles in human health as well as in the plant defense system. The isoflavone content is influenced by minor-effect quantitative trait loci, which interact with polygenetic and environmental factors. It has been difficult to clarify the regulation of isoflavone biosynthesis because of its complex heritability and the influence of external factors. Here, using a genotype-by-sequencing-based genome-wide association mapping study, 189 mutant soybean genotypes (the mutant diversity pool, MDP) were genotyped on the basis of 25,646 high-quality single nucleotide polymorphisms (SNPs) with minor allele frequency of >0.01 except for missing data. All the accessions were phenotyped by determining the contents of 12 isoflavones in the soybean seeds in two consecutive years (2020 and 2021). Then, quantitative trait nucleotides (QTNs) related to isoflavone contents were identified and validated using multi-locus GWAS models. A total of 112 and 46 QTNs related to isoflavone contents were detected by multiple MLM-based models in 2020 and 2021, respectively. Of these, 12 and 5 QTNs were related to more than two types of isoflavones in 2020 and 2021, respectively. Forty-four QTNs were detected within the 441-Kb physical interval surrounding Gm05:38940662. Of them, four QTNs (Gm05:38936166, Gm05:38936167, Gm05:38940662, and Gm05:38940717) were located at Glyma.05g206900 and Glyma.05g207000, which encode glutathione S-transferase THETA 1 (GmGSTT1), as determined from previous quantitative trait loci annotations and the literature. We detected substantial differences in the transcript levels of GmGSTT1 and two other core genes (IFS1 and IFS2) in the isoflavone biosynthetic pathway between the original cultivar and its mutant. The results of this study provide new information about the factors affecting isoflavone contents in soybean seeds and will be useful for breeding soybean lines with high and stable concentrations of isoflavones.

Integrated Approach in Genomic Selection to Accelerate Genetic Gain in Sugarcane

Marker-assisted selection (MAS) has been widely used in the last few decades in plant breeding programs for the mapping and introgression of genes for economically important traits, which has enabled the development of a number of superior cultivars in different crops. In sugarcane, which is the most important source for sugar and bioethanol, marker development work was initiated long ago; however, marker-assisted breeding in sugarcane has been lagging, mainly due to its large complex genome, high levels of polyploidy and heterozygosity, varied number of chromosomes, and use of low/medium-density markers. Genomic selection (GS) is a proven technology in animal breeding and has recently been incorporated in plant breeding programs. GS is a potential tool for the rapid selection of superior genotypes and accelerating breeding cycle. However, its full potential could be realized by an integrated approach combining high-throughput phenotyping, genotyping, machine learning, and speed breeding with genomic selection. For better understanding of GS integration, we comprehensively discuss the concept of genetic gain through the breeder’s equation, GS methodology, prediction models, current status of GS in sugarcane, challenges of prediction accuracy, challenges of GS in sugarcane, integrated GS, high-throughput phenotyping (HTP), high-throughput genotyping (HTG), machine learning, and speed breeding followed by its prospective applications in sugarcane improvement.

Genotyping by Sequencing Advancements in Barley

Barley is considered an ideal crop to study cereal genetics due to its close relationship with wheat and diploid ancestral genome. It plays a crucial role in reducing risks to global food security posed by climate change. Genetic variations in the traits of interest in crops are vital for their improvement. DNA markers have been widely used to estimate these variations in populations. With the advancements in next-generation sequencing, breeders could access different types of genetic variations within different lines, with single-nucleotide polymorphisms (SNPs) being the most common type. However, genotyping barley with whole genome sequencing (WGS) is challenged by the higher cost and computational demand caused by the large genome size (5.5GB) and a high proportion of repetitive sequences (80%). Genotyping-by-sequencing (GBS) protocols based on restriction enzymes and target enrichment allow a cost-effective SNP discovery by reducing the genome complexity. In general, GBS has opened up new horizons for plant breeding and genetics. Though considered a reliable alternative to WGS, GBS also presents various computational difficulties, but GBS-specific pipelines are designed to overcome these challenges. Moreover, a robust design for GBS can facilitate the imputation to the WGS level of crops with high linkage disequilibrium. The complete exploitation of GBS advancements will pave the way to a better understanding of crop genetics and offer opportunities for the successful improvement of barley and its close relatives.

Developing and deploying an efficient genotyping workflow for accelerating maize improvement in developing countries

Background: Molecular breeding is an essential tool for accelerating genetic gain in crop improvement towards meeting the need to feed an ever-growing world population. Establishing low-cost, flexible genotyping platforms in small, public and regional laboratories can stimulate the application of molecular breeding in developing countries. These laboratories can serve plant breeding projects requiring low- to medium-density markers for marker-assisted selection (MAS) and quality control (QC) activities. Methods: We performed two QC and MAS experiments consisting of 637 maize lines, using an optimised genotyping workflow involving an in-house competitive allele-specific PCR (KASP) genotyping system with an optimised sample collection, preparation, and DNA extraction and quantitation process. A smaller volume of leaf-disc size plant samples was collected directly in 96-well plates for DNA extraction, using a slightly modified CTAB-based DArT DNA extraction protocol. DNA quality and quantity analyses were performed using a microplate reader, and the KASP genotyping and data analysis was performed in our laboratory. Results: Applying the optimized genotyping workflow expedited the QC and MAS experiments from over five weeks (when outsourcing) to two weeks and eliminated the shipping cost. Using a set of 28 KASP single nucleotide polymorphisms (SNPs) validated for maize, the QC experiment revealed the genetic identity of four maize varieties taken from five seed sources. Another set of 10 KASP SNPs was sufficient in verifying the parentage of 390 F1 lines. The KASP-based MAS was successfully applied to a maize pro-vitamin A (PVA) breeding program and for introgressing the aflatoxin resistance gene into elite tropical maize lines. Conclusion: This improved workflow has helped accelerate maize improvement activities of IITA's Maize Improvement Program and facilitated DNA fingerprinting for tracking improved crop varieties. National Agricultural Research Systems (NARS) in developing countries can adopt this workflow to fast-track molecular marker-based genotyping for crop improvement.

Identifying Parameters for Defining “Essentially Derived Varieties” of Maize Inbred Lines Using High-Throughput Genome-Wide SNP Markers

Well-developed maize reference genomes and genotyping technology along with fast decreasing detection costs have enabled the chance of shifting essentially derived varieties (EDV) identification to high-throughput SNP genotyping technology. However, attempts of using high-throughput technologies such as SNP array on EDV identification and the essential baseline parameters such as genetic homozygosity and/or stability in EDV practices have not been characterized. Here, we selected 28 accessions of 21 classical maize inbreds, which definitely form a pedigree network from initial founders to derivatives that had made huge contribution to corn production, to demonstrate these fundamental analyses. Our data showed that average residual heterozygosity (RH) rate of these 28 accessions across genome was about 1.03%. However, the RH rate of some accessions was higher than 3%. In addition, some inbreds were found to have an overall RH rate lower than 2% but over 8% level at certain chromosomes. Genetic drift (GD) between two accessions from different years or breeding programs varied from 0.13% to 13.16%. Accessions with low GD level showed cluster distribution pattern and compared with RH distributions indicated that RH was not the only resource of GD. Both RH and GD data suggested that genetic purity analysis is an essential procedure before determining EDV. Eleven derivative lines were characterized with regard to their genome compositions and were inferred as their breeding histories. The backcross, bi-parental recycling, and mutation breeding records could be identified. The data provide insights of underlining fundamental parameters for defining EDV threshold and the results demonstrate the EDV identification process.

Epigenomics as Potential Tools for Enhancing Magnitude of Breeding Approaches for Developing Climate Resilient Chickpea

Epigenomics has become a significant research interest at a time when rapid environmental changes are occurring. Epigenetic mechanisms mainly result from systems like DNA methylation, histone modification, and RNA interference. Epigenetic mechanisms are gaining importance in classical genetics, developmental biology, molecular biology, cancer biology, epidemiology, and evolution. Epigenetic mechanisms play important role in the action and interaction of plant genes during development, and also have an impact on classical plant breeding programs, inclusive of novel variation, single plant heritability, hybrid vigor, plant-environment interactions, stress tolerance, and performance stability. The epigenetics and epigenomics may be significant for crop adaptability and pliability to ambient alterations, directing to the creation of stout climate-resilient elegant crop cultivars. In this review, we have summarized recent progress made in understanding the epigenetic mechanisms in plant responses to biotic and abiotic stresses and have also tried to provide the ways for the efficient utilization of epigenomic mechanisms in developing climate-resilient crop cultivars, especially in chickpea, and other legume crops.

Characterization of Histone H3 Gene Family Reveals That GmHH3-3 is Associated With Higher Seed Weight in Glycine max

The main function of histone protein is to provide support to the structure of chromosomes. It helps in binding a long thread of DNA into a more condensed shape to fit into the nucleus. From histone variants, histone H3 (HH3) plays a crucial role in plant growth and development. Characterization of histones has not been reported in Glycine max till now. The objective of this study was to characterize the HH3 gene family for molecular breeding of G. max. In this study, 17 HH3 members in G. max were identified by performing local BLASTp using HH3 members from Arabidopsis as a query. Phylogenetic analysis classified HH3 genes in seven clades. Sequence logo analysis among Arabidopsis thaliana, Oryza sativa, and Glycine max showed a higher level of similarity in amino acids. Furthermore, conserveness of G. max HH3 genes was also confirmed by Gene Structure Display. Ten paralogous gene pairs were identified in GmHH3 genes in the Glycine max genome by conducting collinearity analysis. G. max HH3 genes have experienced strong purifying selection pressure, with limited functional divergence originating from the segmental and whole-genome duplication, as evidenced by the Ka/Ks ratio. The KASP marker was developed for GmHH3-3 gene. Genotyping was performed on 46 G. max genotypes. This differentiation was based upon the presence of either GmHH3-3-C or GmHH3-3-T allele in the CDS region. The results showed that G. max accessions containing the GmHH3-3-T allele at respective locus showed higher thousand seed weight than that of those accessions that contain the GmHH3-3-C allele. This research provides the basic information to further decipher the function of HH3 in soybean.

HTPdb and HTPtools: Exploiting maize haplotype-tag polymorphisms for germplasm resource analyses and genomics-informed breeding

Along with rapid advances in high-throughput-sequencing technology, the development and application of molecular markers has been critical for the progress that has been made in crop breeding and genetic research. Desirable molecular markers should be able to rapidly genotype tens of thousands of breeding accessions with tens to hundreds of markers. In this study, we developed a multiplex molecular marker, the haplotype-tag polymorphism (HTP), that integrates Maize6H-60K array data from 3,587 maize inbred lines with 6,375 blocks from the recombination block map. After applying strict filtering criteria, we obtained 6,163 highly polymorphic HTPs, which were evenly distributed in the genome. Furthermore, we developed a genome-wide HTP analysis toolkit, HTPtools, which we used to establish an HTP database (HTPdb) covering the whole genomes of 3,587 maize inbred lines commonly used in breeding. A total of 172,921 non-redundant HTP allelic variations were obtained. Three major HTPtools modules combine seven algorithms (e.g., chain Bayes probability and the heterotic-pattern prediction algorithm) and a new plotting engine named "BCplot" that enables rapid visualization of the background information of multiple backcross groups. HTPtools was designed for big-data analyses such as complex pedigree reconstruction and maize heterotic-pattern prediction. The HTP-based analytical strategy and the toolkit developed in this study are applicable for high-throughput genotyping and for genetic mapping, germplasm resource analyses, and genomics-informed breeding in maize.

The Role of Hydroxycinnamic Acid Amide Pathway in Plant Immunity

The compounds involved in the hydroxycinnamic acid amide (HCAA) pathway are an important class of metabolites in plants. Extensive studies have reported that a variety of plant hydroxycinnamamides exhibit pivotal roles in plant-pathogen interactions, such as p-coumaroylagmatine and ferulic acid. The aim of this review is to discuss the emerging findings on the functions of hydroxycinnamic acid amides (HCAAs) accumulation associated with plant defenses against plant pathologies, antimicrobial activity of HCAAs, and the mechanism of HCAAs involved in plant immune responses (such as reactive oxygen species (ROS), cell wall response, plant defense hormones, and stomatal immunity). However, these advances have also revealed the complexity of HCAAs participation in plant defense reactions, and many mysteries remain to be revealed. This review provides an overview of the mechanistic and conceptual insights obtained so far and highlights areas for future exploration of phytochemical defense metabolites.

Developing and deploying an efficient genotyping workflow for accelerating maize improvement in developing countries

Background: Molecular breeding is an essential tool for accelerating genetic gain in crop improvement, towards meeting the need to feed an ever-growing world population. Establishing low-cost, flexible genotyping platforms in small, public and regional laboratories can stimulate the application of molecular breeding in developing countries. These laboratories can serve plant breeding projects requiring low- to medium-density markers for marker-assisted selection (MAS) and quality control (QC) activities. Methods: We performed two QC and MAS experiments consisting of 637 maize lines, using an optimised genotyping workflow involving an in-house competitive allele-specific PCR (KASP) genotyping system with an optimised sample collection, preparation, and DNA extraction and quantitation process. A smaller volume of leaf-disc size plant samples was collected directly in 96-well plates for DNA extraction, using a slightly modified CTAB-based DArT DNA extraction protocol. DNA quality and quantity analyses were performed using a microplate reader, and the KASP genotyping and data analysis was performed in our laboratory. Results: Applying the optimized genotyping workflow expedited the QC and MAS experiments from over five weeks (when outsourcing) to two weeks and eliminated the shipping cost. Using a set of 28 KASP single nucleotide polymorphisms (SNPs) validated for maize, the QC experiment revealed the genetic identity of four maize varieties taken from five seed sources. Another set of 10 KASP SNPs was sufficient in verifying the parentage of 390 F1 lines. The KASP-based MAS was successfully applied to a maize pro-vitamin A (PVA) breeding program and for introgressing the aflatoxin resistance gene into elite tropical maize lines. Conclusion: This improved workflow has helped accelerate maize improvement activities of IITA's Maize Improvement Program and facilitated DNA fingerprinting for tracking improved crop varieties. National Agricultural Research Systems (NARS) in developing countries can adopt this workflow to fast-track molecular marker-based genotyping for crop improvement.

Development of Breeder-Friendly KASP Markers from Genome-Wide Association Studies Results

Array-based SNP markers are commonly used in genome-wide association studies (GWAS) to identify genomic regions involved in important agronomical traits. However, conversion of these SNP markers into breeder-friendly kompetitive allele-specific PCR (KASP) markers for use in marker-assisted selection is often challenging. In this chapter we describe general considerations and successfully applied protocols for the conversion of Illumina array SNP markers into locus-specific KASP markers with a special emphasis and examples on how to overcome difficulties in polyploid wheat.

Asexually propagated Agave tequilana var. azul exhibits variation in genetic markers and defence responses to Fusarium solani

Agave (Agave tequilana var. azul) is considered a crop with low genetic diversity because it has been propagated vegetatively for centuries for commercial purposes, and consequently, it could be equally susceptible to pests and diseases. However, the present study employs plant material derived from field-grown plants exhibiting phenotypic variability in susceptibility to agave wilt. The offshoots from rhizomes of these plants were reproduced in vitro and classified as potentially resistant or susceptible. Amplified fragment length polymorphism analysis confirmed wide genetic differences among individuals, but these differences were not correlated with the observed phenotypic variability in resistance. Propagated plantlets were inoculated with Fusarium solani in two time-lapse confrontations for 72 h and 30 days. The early biochemical response showed statistically superior levels in the accumulation of shikimic acid, phenolic compounds, and chitinase activity in potentially resistant plantlets. There was an inverse correlation of these early biochemical responses and salicylic acid and the incidence of diseased root cells in isogenic plantlets in the 30-day confrontation with F. solani, suggesting that these activities and accumulation of molecules were primordial in the defence against this pathogen.

Five plants per RIL for phenotyping traits of high or moderate heritability ensure the power of QTL mapping in a rice MAGIC population

Currently, the power of QTL mapping is mainly dependent on the quality of phenotypic data in a given population, regardless of the statistical method, as the quality of genotypic data is easily guaranteed in the laboratory. Increasing the sample size per line used for phenotyping is a good way to improve the quality of phenotypic data. However, accommodating a large-scale mapping population takes a large area of rice field, which frequently results in high costs and extra environmental noises. To acquire a reasonable small sample size without a penalty in mapping power, we conducted three experiments with a 4-way MAGIC population and measured phenotypes of 5, 10, and 20 plants per RIL. Three traits including heading date, plant height, and tillers per plant were focused. With SNP- and bin-based QTL mapping, 3 major and 3 minor QTLs for heading date with high heritability and 2 major QTLs for plant height with moderate heritability were commonly detected across the three experiments, but no QTL for tillers per plant with low heritability were commonly identified. In addition, bin-based QTL mapping was more powerful than SNP-based mapping and able to rank the genetic effects of parental alleles. Thus, 5 plants per RIL for phenotyping ensure the power of QTL mapping for traits of high or moderate heritability, and bin-based QTL mapping is recommended for multiparent populations.

Drought-Tolerance QTLs Associated with Grain Yield and Related Traits in Spring Bread Wheat

The present research aims to identify the efficient combination of drought-tolerance selection criteria and associated quantitative trait loci. A panel of 197 bread wheat genotypes was evaluated for yield- and drought-tolerance-related traits in two environments (favorable and semiarid) for 2 years (2015-2016). Grain number, biomass, number of fertile spikes per plant and ground cover exhibited a significant correlation with grain yield and constitute potential secondary selection criteria for yield under drought conditions. About 73 significant marker-trait associations were detected along various chromosomal positions. The markers "wsnp_Ex_Rep_c67786_66472676" and "ExcalibuR_c24593_1217" exhibited important genetic gains associated with yield increase under drought (11 and 7%, respectively). The markers "KukRi_c94792_127" and "wsnp_Ex_c298_580660" showed a significant correlation with grain yield, biomass and grain number and were associated with a significant increase in yield performance at the semiarid site (+6 and +7%, respectively). The ground cover was found associated with grain yield and biomass through the markers "wsnp_Ex_Rep_c67786_66472676" (+11%) and "KukRi_c49927_151" (+10%). One marker "TduRuM_contig25432_1377" on chromosome 5B at 20 cM was consistently correlated with the number of fertile spikes across both environments. Further research should be considered to validate the efficiency of these markers to undertake selection for drought tolerance under various environments and genetic backgrounds.

Dissection of the practical soybean breeding pipeline by developing ZDX1, a high-throughput functional array

KEY MESSAGE

We developed the ZDX1 high-throughput functional soybean array for high accuracy evaluation and selection of both parents and progeny, which can greatly accelerate soybean breeding. Microarray technology facilitates rapid, accurate, and economical genotyping. Here, using resequencing data from 2214 representative soybean accessions, we developed the high-throughput functional array ZDX1, containing 158,959 SNPs, covering 90.92% of soybean genes and sites related to important traits. By application of the array, a total of 817 accessions were genotyped, including three subpopulations of candidate parental lines, parental lines and their progeny from practical breeding. The fixed SNPs were identified in progeny, indicating artificial selection during the breeding process. By identifying functional sites of target traits, novel soybean cyst nematode-resistant progeny and maturity-related novel sources were identified by allele combinations, demonstrating that functional sites provide an efficient method for the rapid screening of desirable traits or gene sources. Notably, we found that the breeding index (BI) was a good indicator for progeny selection. Superior progeny were derived from the combination of distantly related parents, with at least one parent having a higher BI. Furthermore, new combinations based on good performance were proposed for further breeding after excluding redundant and closely related parents. Genomic best linear unbiased prediction (GBLUP) analysis was the best analysis method and achieved the highest accuracy in predicting four traits when comparing SNPs in genic regions rather than whole genomic or intergenic SNPs. The prediction accuracy was improved by 32.1% by using progeny to expand the training population. Collectively, a versatile assay demonstrated that the functional ZDX1 array provided efficient information for the design and optimization of a breeding pipeline for accelerated soybean breeding.

Deciphering Haplotypic Variation and Gene Expression Dynamics Associated with Nutritional and Cooking Quality in Rice

Nutritional quality improvement of rice is the key to ensure global food security. Consequently, enormous efforts have been made to develop genomics and transcriptomics resources for rice. The available omics resources along with the molecular understanding of trait development can be utilized for efficient exploration of genetic resources for breeding programs. In the present study, 80 genes known to regulate the nutritional and cooking quality of rice were extensively studied to understand the haplotypic variability and gene expression dynamics. The haplotypic variability of selected genes were defined using whole-genome re-sequencing data of ~4700 diverse genotypes. The analytical workflow identified 133 deleterious single-nucleotide polymorphisms, which are predicted to affect the gene function. Furthermore, 788 haplotype groups were defined for 80 genes, and the distribution and evolution of these haplotype groups in rice were described. The nucleotide diversity for the selected genes was significantly reduced in cultivated rice as compared with that in wild rice. The utility of the approach was successfully demonstrated by revealing the haplotypic association of chalk5 gene with the varying degree of grain chalkiness. The gene expression atlas was developed for these genes by analyzing RNA-Seq transcriptome profiling data from 102 independent sequence libraries. Subsequently, weighted gene co-expression meta-analyses of 11,726 publicly available RNAseq libraries identified 19 genes as the hub of interactions. The comprehensive analyses of genetic polymorphisms, allelic distribution, and gene expression profiling of key quality traits will help in exploring the most desired haplotype for grain quality improvement. Similarly, the information provided here will be helpful to understand the molecular mechanism involved in the development of nutritional and cooking quality traits in rice.

Identification and validation of QTLs for kernel number per spike and spike length in two founder genotypes of wheat

BACKGROUND

Kernel number per spike (KNS) and spike length (SL) are important spike-related traits in wheat variety improvement. Discovering genetic loci controlling these traits is necessary to elucidate the genetic basis of wheat yield traits and is very important for marker-assisted selection breeding.

RESULTS

In the present study, we used a recombinant inbred line population with 248 lines derived from the two founder genotypes of wheat, Bima4 and BainongAK58, to construct a high-density genetic map using wheat 55 K genotyping assay. The final genetic linkage map consists of 2356 bin markers (14,812 SNPs) representing all 21 wheat chromosomes, and the entire map spanned 4141.24 cM. A total of 7 and 18 QTLs were identified for KNS and SL, respectively, and they were distributed on 11 chromosomes. The allele effects of the flanking markers for 12 stable QTLs, including four QTLs for KNS and eight QTLs for SL, were estimated based on phenotyping data collected from 15 environments in a diverse wheat panel including 384 elite cultivars and breeding lines. The positive alleles at seven loci, namely, QKns.his-7D2-1, QKns.his-7D2-2, QSl.his-4A-1, QSl.his-5D1, QSl.his-4D2-2, QSl.his-5B and QSl.his-5A-2, significantly increased KNS or SL in the diverse panel, suggesting they are more universal in their effects and are valuable for gene pyramiding in breeding programs. The transmission of Bima4 allele indicated that the favorite alleles at five loci (QKns.his-7D2-1, QSl.his-5A-2, QSl.his-2D1-1, QSl.his-3A-2 and QSl.his-3B) showed a relatively high frequency or an upward trend following the continuity of generations, suggesting that they underwent rigorous selection during breeding. At two loci (QKns.his-7D2-1 and QSl.his-5A-2) that the positive effects of the Bima4 alleles have been validated in the diverse panel, two and one kompetitive allele-specific PCR (KASP) markers were further developed, respectively, and they are valuable for marker-assisted selection breeding.

CONCLUSION

Important chromosome regions controlling KNS and SL were identified in the founder parents. Our results are useful for knowing the molecular mechanisms of founder parents and future molecular breeding in wheat.

SoySNP618K array: A high-resolution single nucleotide polymorphism platform as a valuable genomic resource for soybean genetics and breeding

Innovations in genomics have enabled the development of low-cost, high-resolution, single nucleotide polymorphism (SNP) genotyping arrays that accelerate breeding progress and support basic research in crop science. Here, we developed and validated the SoySNP618K array (618,888 SNPs) for the important crop soybean. The SNPs were selected from whole-genome resequencing data containing 2,214 diverse soybean accessions; 29.34% of the SNPs mapped to genic regions representing 86.85% of the 56,044 annotated high-confidence genes. Identity-by-state analyses of 318 soybeans revealed 17 redundant accessions, highlighting the potential of the SoySNP618K array in supporting gene bank management. The patterns of population stratification and genomic regions enriched through domestication were highly consistent with previous findings based on resequencing data, suggesting that the ascertainment bias in the SoySNP618K array was largely compensated for. Genome-wide association mapping in combination with reported quantitative trait loci enabled fine-mapping of genes known to influence flowering time, E2 and GmPRR3b, and of a new candidate gene, GmVIP5. Moreover, genomic prediction of flowering and maturity time in 502 recombinant inbred lines was highly accurate (>0.65). Thus, the SoySNP618K array is a valuable genomic tool that can be used to address many questions in applied breeding, germplasm management, and basic crop research.

Chromosome-specific KASP markers for detecting Amblyopyrum muticum segments in wheat introgression lines

Many wild-relative species are being used in prebreeding programs to increase the genetic diversity of wheat (Triticum aestivum L.). Genotyping tools such as single nucleotide polymorphism (SNP)-based arrays and molecular markers have been widely used to characterize wheat-wild relative introgression lines. However, due to the polyploid nature of the recipient wheat genome, it is difficult to develop SNP-based Kompetitive allele-specific polymerase chain reaction (KASP) markers that are codominant to track the introgressions from the wild species. Previous attempts to develop KASP markers have involved both exome- and polymerase chain reaction (PCR)-amplicon-based sequencing of the wild species. But chromosome-specific KASP assays have been hindered by homoeologous SNPs within the wheat genome. This study involved whole genome sequencing of the diploid wheat wild relative Amblyopyrum muticum (Boiss.) Eig and development of a de novo SNP discovery pipeline that generated ∼38,000 SNPs in unique wheat genome sequences. New assays were designed to increase the density of Am. muticum polymorphic KASP markers. With a goal of one marker per 60 Mbp, 335 new KASP assays were validated as diagnostic for Am. muticum in a wheat background. Together with assays validated in previous studies, 498 well distributed chromosome-specific markers were used to recharacterize previously genotyped wheat-Am. muticum doubled haploid (DH) introgression lines. The chromosome-specific nature of the KASP markers allowed clarification of which wheat chromosomes were involved with recombination events or substituted with Am. muticum chromosomes and the higher density of markers allowed detection of new small introgressions in these DH lines.

Identification and DNA Marker Development for a Wheat-Leymus mollis 2Ns (2D) Disomic Chromosome Substitution

Leymus mollis (2n = 4x = 28, NsNsXmXm), a wild relative of common wheat (Triticum aestivum L.), carries numerous loci which could potentially be used in wheat improvement. In this study, line 17DM48 was isolated from the progeny of a wheat and L. mollis hybrid. This line has 42 chromosomes forming 21 bivalents at meiotic metaphase I. Genomic in situ hybridization (GISH) demonstrated the presence of a pair chromosomes from the Ns genome of L. mollis. This pair substituted for wheat chromosome 2D, as shown by fluorescence in situ hybridization (FISH), DNA marker analysis, and hybridization to wheat 55K SNP array. Therefore, 17DM48 is a wheat-L. mollis 2Ns (2D) disomic substitution line. It shows longer spike and a high level of stripe rust resistance. Using specific-locus amplified fragment sequencing (SLAF-seq), 13 DNA markers were developed to identify and trace chromosome 2Ns of L. mollis in wheat background. This line provides a potential bridge germplasm for genetic improvement of wheat stripe rust resistance.

Genome-Wide Association Study on Root System Architecture and Identification of Candidate Genes in Wheat (Triticum aestivum L.)

The root tissues play important roles in water and nutrient acquisition, environmental adaptation, and plant development. In this study, a diversity panel of 388 wheat accessions was collected to investigate nine root system architecture (RSA) traits at the three-leaf stage under two growing environments: outdoor pot culture (OPC) and indoor pot culture (IPC). Phenotypic analysis revealed that root development was faster under OPC than that under IPC and a significant correlation was observed between the nine RSA traits. The 660K single-nucleotide polymorphism (SNP) chip was used for a genome-wide association study (GWAS). Significant SNPs with a threshold of −log10 (p-value) ≥ 4 were considered. Thus, 36 quantitative trait loci (QTLs), including 13 QTL clusters that were associated with more than one trait, were detected, and 31 QTLs were first identified. The QTL clusters on chromosomes 3D and 5B were associated with four and five RSA traits, respectively. Two candidate genes, TraesCS2A01G516200 and TraesCS7B01G036900, were found to be associated with more than one RSA trait using haplotype analysis, and preferentially expressed in the root tissues. These favourable alleles for RSA traits identified in this study may be useful to optimise the root system in wheat.

Broadening the horizon of crop research: a decade of advancements in plant molecular genetics to divulge phenotype governing genes

Advancements in sequencing, genotyping, and computational technologies during the last decade (2011-2020) enabled new forward-genetic approaches, which subdue the impediments of precise gene mapping in varied crops. The modern crop improvement programs rely heavily on two major steps-trait-associated QTL/gene/marker's identification and molecular breeding. Thus, it is vital for basic and translational crop research to identify genomic regions that govern the phenotype of interest. Until the advent of next-generation sequencing, the forward-genetic techniques were laborious and time-consuming. Over the last 10 years, advancements in the area of genome assembly, genotyping, large-scale data analysis, and statistical algorithms have led faster identification of genomic variations regulating the complex agronomic traits and pathogen resistance. In this review, we describe the latest developments in genome sequencing and genotyping along with a comprehensive evaluation of the last 10-year headways in forward-genetic techniques that have shifted the focus of plant research from model plants to diverse crops. We have classified the available molecular genetic methods under bulk-segregant analysis-based (QTL-seq, GradedPool-Seq, QTG-Seq, Exome QTL-seq, and RapMap), target sequence enrichment-based (RenSeq, AgRenSeq, and TACCA), and mutation-based groups (MutMap, NIKS algorithm, MutRenSeq, MutChromSeq), alongside improvements in classical mapping and genome-wide association analyses. Newer methods for outcrossing, heterozygous, and polyploid plant genetics have also been discussed. The use of k-mers has enriched the nature of genetic variants which can be utilized to identify the phenotype-causing genes, independent of reference genomes. We envisage that the recent methods discussed herein will expand the repertoire of useful alleles and help in developing high-yielding and climate-resilient crops.

Breeding for Economically and Environmentally Sustainable Wheat Varieties: An Integrated Approach from Genomics to Selection

There is currently a strong societal demand for sustainability, quality, and safety in bread wheat production. To address these challenges, new and innovative knowledge, resources, tools, and methods to facilitate breeding are needed. This starts with the development of high throughput genomic tools including single nucleotide polymorphism (SNP) arrays, high density molecular marker maps, and full genome sequences. Such powerful tools are essential to perform genome-wide association studies (GWAS), to implement genomic and phenomic selection, and to characterize the worldwide diversity. This is also useful to breeders to broaden the genetic basis of elite varieties through the introduction of novel sources of genetic diversity. Improvement in varieties particularly relies on the detection of genomic regions involved in agronomical traits including tolerance to biotic (diseases and pests) and abiotic (drought, nutrient deficiency, high temperature) stresses. When enough resolution is achieved, this can result in the identification of candidate genes that could further be characterized to identify relevant alleles. Breeding must also now be approached through in silico modeling to simulate plant development, investigate genotype × environment interactions, and introduce marker-trait linkage information in the models to better implement genomic selection. Breeders must be aware of new developments and the information must be made available to the world wheat community to develop new high-yielding varieties that can meet the challenge of higher wheat production in a sustainable and fluctuating agricultural context. In this review, we compiled all knowledge and tools produced during the BREEDWHEAT project to show how they may contribute to face this challenge in the coming years.