Status and Prospects of Association Mapping in Plants

Understanding natural genetic variation for nutritional quality in grain and identification of superior haplotypes in deepwater rice genotypes of Assam, India

Rice grown under deepwater ecosystem is considered to be natural farming and hence they are considered to be input efficient. Thus, to identify gene responsible for nutritional content under natural conditions, a genome-wide association study (GWAS)was performed. GWAS identified single nucleotide polymorphisms (SNPs) significantly associated with various nutritional quality traits such as Zn (mg/kg), Fe (mg/kg), Protein (%), Oil (%), Amylose (%), Starch (%), Phytic acid (%), Phenol (%) and TDF (%) in 184 deepwater rice accessions evaluated over 2 consecutive years. A total of 278 SNPs distributed across 12 chromosomes were found to be significantly associated with Zn, Oil and Phenol content. Among them, eight high confidence SNPs were significant and identified on chr1 (AX-95933712), chr7 (AX-95957036), and chr8 (AX-95965181) for Zn content. Similarly, on chr2 (AX-95945186), chr8 (AX-95964718), and chr11 (AX-95961099) have been found to be associated with Oil content and on chr3 (AX-95922121) and chr4 (AX-95963889) for Phenol content. Genomic regions of ± 220 kb flanking the three consistent lowest p value containing SNPs for each trait were considered for finding superior haplotypes. These SNPs showed significant phenotypic variations with different identified haplotype blocks. The allelic variations with phenotypes were considered to be superior haplotypes i.e., Block 1: Hap 1 (ACCC) for high Zn content, Block 2: Hap 1 (CT) for high Oil content, and Block 2: Hap 1(CGGG) for low Phenol content. The discovered superior haplotype with high nutritional content could be important for understanding the mechanisms involving nutrient use efficiency. Thus, the present study demonstrated that developing rice varieties with appropriate nutritional quality traits will be possible through the incorporation of such superior haplotypes in breeding programs.

Genome-wide association study provides new insight into the underlying mechanism of drought tolerance during seed germination stage in soybean

Drought is one of the major environmental issues that reduce crop yield. Seed germination is a crucial stage of plant development in all crop plants, including soybean. In soybean breeding, information about genetic mechanism of drought tolerance has great importance. However, at germination stage, there is relatively little knowledge on the genetic basis of soybean drought resistance. The objective of this work was to find the quantitative trait nucleotides (QTNs) linked to drought tolerance related three traits using a genome-wide association study (GWAS), viz., germination rate (GR), root length (RL), and whole seedling length (WSL), using germplasm population of 240 soybean PIs with 34,817 SNPs genotype data having MAF > 0.05. It was observed that heritability (H2) for GR, WSL, and RL across both environments (2020, and 2019) were high in the range of 0.76-0.99, showing that genetic factors play a vital role in drought tolerance as compared to environmental factors. A number of 23 and 27 QTNs were found to be linked to three traits using MLM and mrMLM, respectively. Three significant QTNs, qGR8-1, qWSL13-1, and qRL-8, were identified using both MLM and mrMLM methods among these QTNs. QTN8, located on chromosome 8 was consistently linked to two traits (GR and RL). The area (± 100 Kb) associated with this QTN was screened for drought tolerance based on gene annotation. Fifteen candidate genes were found by this screening. Based on the expression data, four candidate genes i.e. Glyma08g156800, Glyma08g160000, Glyma08g162700, and Glyma13g249600 were found to be linked to drought tolerance regulation in soybean. Hence, the current study provides evidence to understand the genetic constitution of drought tolerance during the germination stage and identified QTNs or genes could be utilized in molecular breeding to enhance the yield under drought stress.

Genome-wide association study of fiber quality traits in US upland cotton (Gossypium hirsutum L.)

KEY MESSAGE

A GWAS in an elite diversity panel, evaluated across 10 environments, identified genomic regions regulating six fiber quality traits, facilitating genomics-assisted breeding and gene discovery in upland cotton. In this study, an elite diversity panel of 348 upland cotton accessions was evaluated in 10 environments across the US Cotton Belt and genotyped with the cottonSNP63K array, for a genome-wide association study of six fiber quality traits. All fiber quality traits, upper half mean length (UHML: mm), fiber strength (FS: g tex-1), fiber uniformity (FU: %), fiber elongation (FE: %), micronaire (MIC) and short fiber content (SFC: %), showed high broad-sense heritability (> 60%). All traits except FE showed high genomic heritability. UHML, FS and FU were all positively correlated with each other and negatively correlated with FE, MIC and SFC. GWAS of these six traits identified 380 significant marker-trait associations (MTAs) including 143 MTAs on 30 genomic regions. These 30 genomic regions included MTAs identified in at least three environments, and 23 of them were novel associations. Phenotypic variation explained for the MTAs in these 30 genomic regions ranged from 6.68 to 11.42%. Most of the fiber quality-associated genomic regions were mapped in the D-subgenome. Further, this study confirmed the pleiotropic region on chromosome D11 (UHML, FS and FU) and identified novel co-localized regions on D04 (FU, SFC), D05 (UHML, FU, and D06 UHML, FU). Marker haplotype analysis identified superior combinations of fiber quality-associated genomic regions with high trait values (UHML = 32.34 mm; FS = 32.73 g tex-1; FE = 6.75%). Genomic analyses of traits, haplotype combinations and candidate gene information described in the current study could help leverage genetic diversity for targeted genetic improvement and gene discovery for fiber quality traits in cotton.

Haplotype-based breeding: A new insight in crop improvement

Haplotype-based breeding (HBB) is one of the cutting-edge technologies in the realm of crop improvement due to the increasing availability of Single Nucleotide Polymorphisms identified by Next Generation Sequencing technologies. The complexity of the data can be decreased with fewer statistical tests and a lower probability of spurious associations by combining thousands of SNPs into a few hundred haplotype blocks. The presence of strong genomic regions in breeding lines of most crop species facilitates the use of haplotypes to improve the efficiency of genomic and marker-assisted selection. Haplotype-based breeding as a Genomic Assisted Breeding (GAB) approach harnesses the genome sequence data to pinpoint the allelic variation used to hasten the breeding cycle and circumvent the challenges associated with linkage drag. This review article demonstrates ways to identify candidate genes, superior haplotype identification, haplo-pheno analysis, and haplotype-based marker-assisted selection. The crop improvement strategies that utilize superior haplotypes will hasten the breeding progress to safeguard global food security.

Genetic dissection of value-added quality traits and agronomic parameters through genome-wide association mapping in bread wheat (T. aestivum L.)

Bread wheat (T. aestivum) is one of the world's most widely consumed cereals. Since micronutrient deficiencies are becoming more common among people who primarily depend upon cereal-based diets, a need for better-quality wheat varieties has been felt. An association panel of 154 T. aestivum lines was evaluated for the following quality traits: grain appearance (GA) score, grain hardness (GH), phenol reaction (PR) score, protein percent, sodium dodecyl sulfate (SDS) sedimentation value, and test weight (TWt). In addition, the panel was also phenotyped for grain yield and related traits such as days to heading, days to maturity, plant height, and thousand kernel weight for the year 2017-18 at the Borlaug Institute for South Asia (BISA) Ludhiana and Jabalpur sites. We performed a genome-wide association analysis on this panel using 18,351 genotyping-by-sequencing (GBS) markers to find marker-trait associations for quality and grain yield-related traits. We detected 55 single nucleotide polymorphism (SNP) marker trait associations (MTAs) for quality-related traits on chromosomes 7B (10), 1A (9), 2A (8), 3B (6), 2B (5), 7A (4), and 1B (3), with 3A, 4A, and 6D, having two and the rest, 4B, 5A, 5B, and 1D, having one each. Additionally, 20 SNP MTAs were detected for yield-related traits based on a field experiment conducted in Ludhiana on 7D (4) and 4D (3) chromosomes, while 44 SNP MTAs were reported for Jabalpur on chromosomes 2D (6), 7A (5), 2A (4), and 4A (4). Utilizing these loci in marker-assisted selection will benefit from further validation studies for these loci to improve hexaploid wheat for better yield and grain quality.

Characterization of flag leaf morphology identifies a major genomic region controlling flag leaf angle in the US winter wheat (Triticum aestivum L.)

KEY MESSAGE

Multi-environmental characterization of flag leaf morphology traits in the US winter wheat revealed nine stable genomic regions for different flag leaf-related traits including a major region governing flag leaf angle. Flag leaf in wheat is the primary contributor to accumulating photosynthetic assimilates. Flag leaf morphology (FLM) traits determine the overall canopy structure and capacity to intercept the light, thus influencing photosynthetic efficiency. Hence, understanding the genetic control of these traits could be useful for breeding desirable ideotypes in wheat. We used a panel of 272 accessions from the hard winter wheat (HWW) region of the USA to investigate the genetic architecture of five FLM traits including flag leaf length (FLL), width (FLW), angle (FLANG), length-width ratio, and area using multilocation field experiments. Multi-environment GWAS using 14,537 single-nucleotide polymorphisms identified 36 marker-trait associations for different traits, with nine being stable across environments. A novel and major stable region for FLANG (qFLANG.1A) was identified on chromosome 1A accounting for 9-13% variation. Analysis of spatial distribution for qFLANG.1A in a set of 2354 breeding lines from the HWW region showed a higher frequency of allele associated with narrow leaf angle. A KASP assay was developed for allelic discrimination of qFLANG.1A and was used for its independent validation in a diverse set of spring wheat accessions. Furthermore, candidate gene analysis for two regions associated with FLANG identified seven putative genes of interest for each of the two regions. The present study enhances our understanding of the genetic control of FLM in wheat, particularly FLANG, and these results will be useful for dissecting the genes underlying canopy architecture in wheat facilitating the development of climate-resilient wheat varieties.

Leaf Membrane Stability under High Temperatures as an Indicator of Heat Tolerance in Potatoes and Genome-Wide Association Studies to Understand the Underlying Genetics

High temperatures during the crop growing season are becoming more frequent and unpredictable, resulting in reduced crop productivity and quality. Heat stress disrupts plant metabolic processes that affect cell membrane composition and integrity. Cell membrane permeability, ion leakage, and heat shock proteins have been evaluated to screen for heat tolerance in plants. In potatoes, it is unclear whether leaf membrane stability under heat stress is correlated with underground tuber productivity and quality. The main goal of this study was to evaluate if leaf membrane relative electrolyte conductivity (REC) under high temperatures could be used to identify heat-tolerant potato genotypes. Electrolyte leakage assays, correlation estimations, and genome-wide association studies were carried out in 215 genotypes. Expression levels of small heat shock protein 18 (sHSP18) were evaluated in the heat-sensitive potato variety Russet Burbank and compared with those of the heat-tolerant variety Vanguard Russet using Western blotting. Significant differences were observed among genotypes for leaf membrane REC under extreme heat (50°C); REC values ranged from 47.0-99.5%. Leaf membrane REC was positively correlated with tuber external and internal defects and negatively correlated with yield. REC was negatively correlated with the content of several tuber minerals, such as nitrogen, magnesium, and manganese. Eleven quantitative trait loci (QTLs) were identified for leaf membrane REC, explaining up to 13.8% of the phenotypic variance. Gene annotation in QTL areas indicated associations with genes controlling membrane solute transport and plant responses to abiotic stresses. Vanguard Russet had lower leaf REC and higher expression of sHSP18 under high-temperature stress. Our findings indicate that leaf membrane REC under high temperatures can be used as an indicator of potato heat tolerance.

GWAS elucidated grain yield genetics in Indian spring wheat under diverse water conditions

KEY MESSAGE

Underpinned natural variations and key genes associated with yield under different water regimes, and identified genomic signatures of genetic gain in the Indian wheat breeding program. A novel KASP marker for TKW under water stress was developed and validated. A comprehensive genome-wide association study was conducted on 300 spring wheat genotypes to elucidate the natural variations associated with grain yield and its eleven contributing traits under fully irrigated, restricted water, and simulated no water conditions. Utilizing the 35K Wheat Breeders' Array, we identified 1155 quantitative trait nucleotides (QTNs), with 207 QTNs exhibiting stability across diverse conditions. These QTNs were further delimited into 539 genomic regions using a genome-wide LD value of 3.0 Mbp, revealing pleiotropic control across traits and conditions. Sub-genome A was significantly associated with traits under irrigated conditions, while sub-genome B showed more QTNs under water stressed conditions. Favourable alleles with significantly associated QTNs were delineated, with a notable pyramiding effect for enhancing trait performance. Additionally, allele of only 921 QTNs significantly affected the population mean. Allele profiling highlighted C-306 as a most potential source of drought tolerance. Moreover, 762 genes overlapping significant QTNs were identified, narrowing down to 27 putative candidate genes overlapping 29 novel and functional SNPs expressing (≥ 0.5 tpm) relevance across various growth conditions. A new KASP assay was developed, targeting a gene TraesCS2A03G1123700 regulating thousand kernel weight under severe drought condition. Genomic selection models (GBLUP, BayesB, MxE, and R-Norm) demonstrated an average prediction accuracy of 0.06-0.58 across environments, indicating potential for trait selection. Retrospective analysis of the Indian wheat breeding program supported a genetic gain in GY at the rate of ca. 0.56% per breeding cycle, since 1960, supporting the identification of genomic signatures driving trait selection and genetic gain. These findings offer insight into improving the rate of genetic gain in wheat breeding programs globally.

Genome wide association study of seedling and adult plant leaf rust resistance in two subsets of barley genetic resources

Leaf rust (LR) caused by Puccinia hordei is a serious disease of barley worldwide, causing significant yield losses and reduced grain quality. Discovery and incorporation of new sources of resistance from gene bank accessions into barley breeding programs is essential for the development of leaf rust resistant varieties. To identify Quantitative Trait Loci (QTL) conferring LR resistance in the two barley subsets, the Generation Challenge Program (GCP) reference set of 142 accessions and the leaf rust subset constructed using the Focused Identification of Germplasm Strategy (FIGS) of 76 barley accessions, were genotyped to conduct a genome-wide association study (GWAS). The results revealed a total of 59 QTL in the 218 accessions phenotyped against barley leaf rust at the seedling stage using two P. hordei isolates (ISO-SAT and ISO-MRC), and at the adult plant stage in four environments in Morocco. Out of these 59 QTL, 10 QTL were associated with the seedling resistance (SR) and 49 QTL were associated with the adult plant resistance (APR). Four QTL showed stable effects in at least two environments for APR, whereas two common QTL associated with SR and APR were detected on chromosomes 2H and 7H. Furthermore, 39 QTL identified in this study were potentially novel. Interestingly, the sequences of 27 SNP markers encoded the candidate genes (CGs) with predicted protein functions in plant disease resistance. These results will provide new perspectives on the diversity of leaf rust resistance loci for fine mapping, isolation of resistance genes, and for marker-assisted selection for the LR resistance in barley breeding programs worldwide.

Genetic resources and genes/QTLs for gram pod borer (Helicoverpa armigera Hübner) resistance in chickpea from the Western Himalayas

Helicoverpa armigera (also known as gram pod borer) is a serious threat to chickpea production in the world. A set of 173 chickpea genotypes were evaluated for H. armigera resistance, including mean larval population (MLP), percentage pod damage (PPD), and pest resistance (PR) for 2 consecutive years (year 2020 and 2021). The same core set was also genotyped with 50K Axiom CicerSNP Array. The trait data and 50,000 single nucleotide polymorphism genotypic data were used together to work out marker-trait associations (MTAs) using different genome-wide association studies models. For MLP, a total of 53 MTAs were identified, including 25 MTAs in year 2020 and 28 MTAs in year 2021. A set of three MTAs was found common in both environments. For PPD, two MTAs in year 2020 and five MTAs in year 2021 were identified. A set of two MTAs were common in both environments. Similarly, for PR, only two MTAs common in both environments were identified. Interestingly, a common MTA (Affx_123255526) on chromosome 2 (Ca2) was found to be associated with all the three component traits (MLP, PPD, and PR) of pod borer resistance in chickpea. Further, we report key genes that encode SCAMPs (that facilitates the secretion of defense-related molecules), quinone oxidoreductase (enables the production of reactive oxygen species that promotes diapause of gram pod borer), and NB-LRR proteins that have been implicated in plant defense against H. armigera. The resistant chickpea genotypes, MTAs, and key genes reported in the present study may prove useful in the future for developing pod borer-resistant chickpea varieties.

Genome-Wide Association Study for Identification of Marker-Trait Associations Conferring Resistance to Scald from Globally Collected Barley Germplasm

Scald is one of the major economically important foliar diseases in barley, causing up to 40% yield loss in susceptible varieties. The identification of quantitative trait loci and elite alleles that confer resistance to scald is imperative in reducing the threats to barley production. In this study, genome-wide association studies were conducted using a panel of 697 barley genotypes to identify quantitative trait loci for scald resistance. Field experiments were conducted over three consecutive years. Among different models used for genome-wide association studies analysis, FarmCPU was shown to be the best-suited model. Nineteen significant marker-trait associations related to scald resistance were identified across six different chromosomes. Eleven of these marker-trait associations correspond to previously reported scald resistance genes Rrs1, Rrs4, and Rrs2, respectively. Eight novel marker-trait associations were identified in this study, with the candidate genes encoding a diverse class of proteins, including region leucine-rich repeats, AP2/ERF transcription factor, homeodomain-leucine zipper, and protein kinase family proteins. The combination of identified superior alleles significantly reduces disease severity scores. The results will be valuable for marker-assisted breeding for developing scald-resistant varieties.

Mapping the genetic architecture of low-temperature stress tolerance in citron watermelon

Sweet-fleshed watermelon (Citrullus lanatus) is an important vegetable crop of the tropical origin. It is widely grown and consumed around the world for its hydration and nutritional quality values. Low-temperature stress can affect early planting, seedling establishment, and expansion of crop production to new areas. A collection of 122 citron watermelon (Citrullus amarus) accessions were obtained from the USDA's National Plant Germplasm Repository System gene bank in Griffin, GA. The accessions were genotyped using whole genome resequencing to generate single nucleotide polymorphisms (SNPs) molecular markers and screened under cold-stressed and non-stressed control conditions. Four low-temperature stress tolerance related traits including shoot biomass, vine length, maximum quantum efficiency of photosystem II, and chlorophyll content were measured under cold-stressed and non-stressed control treatment conditions. Correlation analysis revealed the presence of positive relationships among traits. Broad-sense heritability for all traits ranged from 0.35 to 0.73, implying the presence of genetic contributions to the observed phenotypic variation. Genomic regions underlying these traits across several citron watermelon chromosomes were identified. Four low-temperature stress tolerance related putative candidate genes co-located with the peak SNPs from genome-wide association study. These genomic regions and marker information could potentially be used in molecular breeding to accelerate genetic improvements for low-temperature stress tolerance in watermelon.

Genome-Wide Association Study Identifies Quantitative Trait Loci and Candidate Genes Involved in Deep-Sowing Tolerance in Maize (Zea mays L.)

Deep sowing is an efficient strategy for maize to ensure the seedling emergence rate under adverse conditions such as drought or low temperatures. However, the genetic basis of deep-sowing tolerance-related traits in maize remains largely unknown. In this study, we performed a genome-wide association study on traits related to deep-sowing tolerance, including mesocotyl length (ML), coleoptile length (CL), plumule length (PL), shoot length (SL), and primary root length (PRL), using 255 maize inbred lines grown in three different environments. We identified 23, 6, 4, and 4 quantitative trait loci (QTLs) associated with ML, CL, PL, and SL, respectively. By analyzing candidate genes within these QTLs, we found a γ-tubulin-containing complex protein, ZmGCP2, which was significantly associated with ML, PL, and SL. Loss of function of ZmGCP2 resulted in decreased PL, possibly by affecting the cell elongation, thus affecting SL. Additionally, we identified superior haplotypes and allelic variations of ZmGCP2 with a longer PL and SL, which may be useful for breeding varieties with deep-sowing tolerance to improve maize cultivation.

Insights into the genetic architecture of Phytophthora capsici root rot resistance in chile pepper (Capsicum spp.) from multi-locus genome-wide association study

BACKGROUND

Phytophthora root rot, a major constraint in chile pepper production worldwide, is caused by the soil-borne oomycete, Phytophthora capsici. This study aimed to detect significant regions in the Capsicum genome linked to Phytophthora root rot resistance using a panel consisting of 157 Capsicum spp. genotypes. Multi-locus genome wide association study (GWAS) was conducted using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS). Individual plants were separately inoculated with P. capsici isolates, 'PWB-185', 'PWB-186', and '6347', at the 4-8 leaf stage and were scored for disease symptoms up to 14-days post-inoculation. Disease scores were used to calculate disease parameters including disease severity index percentage, percent of resistant plants, area under disease progress curve, and estimated marginal means for each genotype.

RESULTS

Most of the genotypes displayed root rot symptoms, whereas five accessions were completely resistant to all the isolates and displayed no symptoms of infection. A total of 55,117 SNP markers derived from GBS were used to perform multi-locus GWAS which identified 330 significant SNP markers associated with disease resistance. Of these, 56 SNP markers distributed across all the 12 chromosomes were common across the isolates, indicating association with more durable resistance. Candidate genes including nucleotide-binding site leucine-rich repeat (NBS-LRR), systemic acquired resistance (SAR8.2), and receptor-like kinase (RLKs), were identified within 0.5 Mb of the associated markers.

CONCLUSIONS

Results will be used to improve resistance to Phytophthora root rot in chile pepper by the development of Kompetitive allele-specific markers (KASP®) for marker validation, genomewide selection, and marker-assisted 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.

Polymorphisms in two key anthocyanic genes of clivia (Clivia miniata L.) reveal evidence of selection and possible association with flower pigmentation

Members of the genus Clivia show considerable variation in flower pigmentation and morphology. Such variation is affected by mutations that emerge in candidate flower development genes over time. Besides population history, mutations can further illuminate the effects of demographic events in populations in addition to population genetic parameters including selection, recombination, and linkage disequilibrium (LD). The current study aimed to find sequence variants in 2 anthocyanin biosynthetic genes (DFR and bHLH) of Clivia miniata and use the data to assess population genetic factors from a random collection of orange/red- and yellow-flowered specimens. Overall, average nucleotide diversity in the 2 anthocyanin genes was moderate (π = 0.00646), whereas haplotypes differed significantly (Hd ≥ 0.9). Gene evolution was seemingly driven by mutations (CmiDFR) or recombinations (CmibHLH001). LD decayed swiftly within the analyzed gene regions and supported the feasibility of assessing trait-variant associations via the association/linkage mapping approach. In the end, most associations were found to be spurious, but 1 haplotype in CmibHLH001 showed a promising correlation to the orange/red flower phenotype in Clivia specimens. In all, the present study is the first to measure gene-level diversity in C. miniata-data that had never been reported so far. Furthermore, the study also identified allelic and haplotypic variants that may be beneficial in future association genetic studies of Clivia. Such studies, however, consider large diverse populations to control for statistical bias intrinsic to the analysis of small datasets.

Genomic selection in plant breeding: Key factors shaping two decades of progress

Genomic selection, the application of genomic prediction (GP) models to select candidate individuals, has significantly advanced in the past two decades, effectively accelerating genetic gains in plant breeding. This article provides a holistic overview of key factors that have influenced GP in plant breeding during this period. We delved into the pivotal roles of training population size and genetic diversity, and their relationship with the breeding population, in determining GP accuracy. Special emphasis was placed on optimizing training population size. We explored its benefits and the associated diminishing returns beyond an optimum size. This was done while considering the balance between resource allocation and maximizing prediction accuracy through current optimization algorithms. The density and distribution of single-nucleotide polymorphisms, level of linkage disequilibrium, genetic complexity, trait heritability, statistical machine-learning methods, and non-additive effects are the other vital factors. Using wheat, maize, and potato as examples, we summarize the effect of these factors on the accuracy of GP for various traits. The search for high accuracy in GP-theoretically reaching one when using the Pearson's correlation as a metric-is an active research area as yet far from optimal for various traits. We hypothesize that with ultra-high sizes of genotypic and phenotypic datasets, effective training population optimization methods and support from other omics approaches (transcriptomics, metabolomics and proteomics) coupled with deep-learning algorithms could overcome the boundaries of current limitations to achieve the highest possible prediction accuracy, making genomic selection an effective tool in plant breeding.

Genome-wide association study and genomic prediction for resistance to brown planthopper in rice

The brown planthopper (BPH) is the most destructive insect pest that threatens rice production globally. Developing rice varieties incorporating BPH-resistant genes has proven to be an effective control measure against BPH. In this study, we assessed the resistance of a core collection consisting of 502 rice germplasms by evaluating resistance scores, weight gain rates and honeydew excretions. A total of 117 rice varieties (23.31%) exhibited resistance to BPH. Genome-wide association studies (GWAS) were performed on both the entire panel of 502 rice varieties and its subspecies, and 6 loci were significantly associated with resistance scores (P value < 1.0e-8). Within these loci, we identified eight candidate genes encoding receptor-like protein kinase (RLK), nucleotide-binding and leucine-rich repeat (NB-LRR), or LRR proteins. Two loci had not been detected in previous study and were entirely novel. Furthermore, we evaluated the predictive ability of genomic selection for resistance to BPH. The results revealed that the highest prediction accuracy for BPH resistance reached 0.633. As expected, the prediction accuracy increased progressively with an increasing number of SNPs, and a total of 6.7K SNPs displayed comparable accuracy to 268K SNPs. Among various statistical models tested, the random forest model exhibited superior predictive accuracy. Moreover, increasing the size of training population improved prediction accuracy; however, there was no significant difference in prediction accuracy between a training population size of 737 and 1179. Additionally, when there existed close genetic relatedness between the training and validation populations, higher prediction accuracies were observed compared to scenarios when they were genetically distant. These findings provide valuable resistance candidate genes and germplasm resources and are crucial for the application of genomic selection for breeding durable BPH-resistant rice varieties.

Multi-omics analysis reveals novel loci and a candidate regulatory gene of unsaturated fatty acids in soybean (Glycine max (L.) Merr)

BACKGROUND

Soybean is a major oil crop; the nutritional components of soybean oil are mainly controlled by unsaturated fatty acids (FA). Unsaturated FAs mainly include oleic acid (OA, 18:1), linoleic acid (LLA, 18:2), and linolenic acid (LNA, 18:3). The genetic architecture of unsaturated FAs in soybean seeds has not been fully elucidated, although many independent studies have been conducted. A 3 V multi-locus random single nucleotide polymorphism (SNP)-effect mixed linear model (3VmrMLM) was established to identify quantitative trait loci (QTLs) and QTL-by-environment interactions (QEIs) for complex traits.

RESULTS

In this study, 194 soybean accessions with 36,981 SNPs were calculated using the 3VmrMLM model. As a result, 94 quantitative trait nucleotides (QTNs) and 19 QEIs were detected using single-environment (QTN) and multi-environment (QEI) methods. Three significant QEIs, namely rs4633292, rs39216169, and rs14264702, overlapped with a significant single-environment QTN.

CONCLUSIONS

For QTNs and QEIs, further haplotype analysis of candidate genes revealed that the Glyma.03G040400 and Glyma.17G236700 genes were beneficial haplotypes that may be associated with unsaturated FAs. This result provides ideas for the identification of soybean lipid-related genes and provides insights for breeding high oil soybean varieties in the future.

A genome-wide association study for resistance to Fusarium wilt (Fusarium oxysporum f. sp. vasinfectum) race 4 in diploid cotton (Gossypium arboreum) and resistance transfer to tetraploid Gossypium hirsutum

Fusarium wilt, caused by the soilborne fungus Fusarium oxysporum f. sp. vasinfectum (FOV), is a devastating disease affecting cotton (Gossypium spp.) worldwide. Understanding the genetic basis of resistance in diploid cotton and successfully transferring the resistance to tetraploid Upland cotton (G. hirsutum) are crucial for developing resistant cotton cultivars. Although numerous studies have been conducted to investigate the genetic basis of Fusarium wilt in tetraploid cotton, little research has been conducted on diploid species. In this study, an association mapping panel consisting of 246 accessions of G. arboreum, was used to identify chromosomal regions for FOV race 4 (FOV4) resistance based on foliar disease severity ratings in four greenhouse tests. Through a genome-wide association study (GWAS) based on 7,009 single nucleotide polymorphic (SNP) markers, 24 FOV4 resistance QTLs, including three major QTLs on chromosomes A04, A06, and A11, were detected. A validation panel consisting of 97 diploid cotton accessions was employed, confirming the presence of several QTLs. Evaluation of an introgressed BC2F7 population derived from G. hirsutum/G. aridum/G. arboreum showed significant differences in disease incidence and mortality rate, as compared to susceptible and resistant controls, suggesting that the resistance in G. arboreum and/or G. aridum was transferred into Upland cotton for the first time. The identification of novel major resistance QTLs, along with the transfer of resistance from the diploid species, expands our understanding of the genomic regions involved in conferring resistance to FOV4 and contributes to the development of resilient Upland cotton cultivars.

Association study of morpho-phenological traits in quinoa (Chenopodium quinoa Willd.) using SSR markers

In this study, the genetic and molecular diversity of 60 quinoa accessions was assessed using agronomically important traits related to grain yield as well as microsatellite (SSR) markers, and informative markers linked to the studied traits were identified using association study. The results showed that most of the studied traits had a relatively high diversity, but grain saponin and protein content showed the highest diversity. High diversity was also observed in all SSR markers, but KAAT023, KAAT027, KAAT036, and KCAA014 showed the highest values for most of the diversity indices and can be introduced as the informative markers to assess genetic diversity in quinoa. Population structure analysis showed that the studied population probably includes two subclusters, so that out of 60 quinoa accessions, 29 (48%) and 23 (38%) accessions were assigned to the first and second subclusters, respectively, and eight (13%) accessions were considered as the mixed genotypes. The study of the population structure using Structure software showed two possible subgroups (K = 2) in the studied population and the results of the bar plot confirmed it. Association study using the general linear model (GLM) and mixed linear model (MLM) identified the number of 35 and 32 significant marker-trait associations (MTAs) for the first year (2019) and 37 and 35 significant MTAs for the second year (2020), respectively. Among the significant MTAs identified for different traits, the highest number of significant MTAs were obtained for grain yield and 1000-grain weight with six and five MTAs, respectively.

A genomic variation map provides insights into peanut diversity in China and associations with 28 agronomic traits

Peanut (Arachis hypogaea L.) is an important allotetraploid oil and food legume crop. China is one of the world's largest peanut producers and consumers. However, genomic variations underlying the migration and divergence of peanuts in China remain unclear. Here we reported a genome-wide variation map based on the resequencing of 390 peanut accessions, suggesting that peanuts might have been introduced into southern and northern China separately, forming two cultivation centers. Selective sweep analysis highlights asymmetric selection between the two subgenomes during peanut improvement. A classical pedigree from South China offers a context for the examination of the impact of artificial selection on peanut genome. Genome-wide association studies identified 22,309 significant associations with 28 agronomic traits, including candidate genes for plant architecture and oil biosynthesis. Our findings shed light on peanut migration and diversity in China and provide valuable genomic resources for peanut improvement.

Genetic Analysis of Adaptive Traits in Spring Wheat in Northeast China

The dissection of the genetic architecture and the detection of the loci for adaptive traits are important for marker-assisted selection (MAS) for breeding. A spring wheat diversity panel with 251 cultivars, mainly from China, was obtained to conduct a genome-wide association study (GWAS) to detect the new loci, including the heading date (HD), maturating date (MD), plant height (PH), and lodging resistance (LR). In total, 41 loci existing in all 21 chromosomes, except for 4A and 6B, were identified, and each explained 4.3-18.9% of the phenotypic variations existing in two or more environments. Of these, 13 loci are overlapped with the known genes or quantitative trait loci (QTLs), whereas the other 28 are likely to be novel. The 1A locus (296.9-297.7 Mb) is a multi-effect locus for LR and PH, whereas the locus on chromosome 6D (464.5-471.0 Mb) affects both the HD and MD. Furthermore, four candidate genes for adaptive traits were identified, involved in cell division, signal transduction, and plant development. Additionally, two competitive, allele-specific PCR (KASP) markers, Kasp_2D_PH for PH and Kasp_6D_HD for HD, were developed and validated in another 162 spring wheat accessions. Our study uncovered the genetic basis of adaptive traits and provided the associated SNPs and varieties with more favorable alleles for wheat MAS breeding.

Genotyping SNPs in lignin biosynthesis gene (CAD1) and transcription factors (MYB1 and MYB2) exhibits association with wood density in teak (Tectona grandis L.f.)

BACKGROUND

Teak (Tectona grandis L.f.), an important source of tropical timber with immense economic value, is a highly outcrossing forest tree species. 150 unrelated accessions of teak (Tectona grandis L.f.) plus trees assembled as clones at National Teak Germplasm Bank, Chandrapur, Maharashtra, India was investigated for association mapping of candidate lignin biosynthesis gene (CAD1) and transcription factors (MYB1 and MYB2).

METHODS AND RESULTS

The CAD1, MYB1 and MYB2 were amplified using specifically designed primers. The amplified sequences were then sequenced and genotyped for 112 SNPs/11 indels. We evaluated the association between SNPs and wood density in teak accessions using GLM and MLM statistical models, with Bonferroni correction applied. The teak accessions recorded an average wood density of 416.69 kg.m-3 (CV 4.97%) and comprised of three loosely structured admixed sub-populations (K = 3), containing 72.05% genetic variation within sub-populations with low intragenic LD (0-21% SNP pairs) at P < 0.05 and high LD decay (33-934 bp) at R2 = 0.1. GLM and MLM models discounting systematic biases (Q and K matrices) to avoid false discovery revealed five loci at rare variants (MAF 0.003) and three loci at common variants (MAF 0.05) to be significantly (P < 0.05) associated with the wood density. However, the stringent Bonferroni correction (4.06-7.04 × 10-4) yielded only a single associated locus (B1485C/A) from exon of MYB1 transcription factor, contributing to about 10.35% phenotypic variation in wood density trait.

CONCLUSION

Scored SNP locus (B1485C/A) can be developed as a molecular probe for selection of improved planting stock with proven wood density trait for a large-scale teak plantation.

Genome-wide association study of common resistance to rust species in tetraploid wheat

Rusts of the genus Puccinia are wheat pathogens. Stem (black; Sr), leaf (brown; Lr), and stripe (yellow; Yr) rust, caused by Puccinia graminis f. sp. tritici (Pgt), Puccinia triticina (Pt), and Puccinia striiformis f. sp. tritici (Pst), can occur singularly or in mixed infections and pose a threat to wheat production globally in terms of the wide dispersal of their urediniospores. The development of durable resistant cultivars is the most sustainable method for controlling them. Many resistance genes have been identified, characterized, genetically mapped, and cloned; several quantitative trait loci (QTLs) for resistance have also been described. However, few studies have considered resistance to all three rust pathogens in a given germplasm. A genome-wide association study (GWAS) was carried out to identify loci associated with resistance to the three rusts in a collection of 230 inbred lines of tetraploid wheat (128 of which were Triticum turgidum ssp. durum) genotyped with SNPs. The wheat panel was phenotyped in the field and subjected to growth chamber experiments across different countries (USA, Mexico, Morocco, Italy, and Spain); then, a mixed linear model (MLM) GWAS was performed. In total, 9, 34, and 5 QTLs were identified in the A and B genomes for resistance to Pgt, Pt, and Pst, respectively, at both the seedling and adult plant stages. Only one QTL on chromosome 4A was found to be effective against all three rusts at the seedling stage. Six QTLs conferring resistance to two rust species at the adult plant stage were mapped: three on chromosome 1B and one each on 5B, 7A, and 7B. Fifteen QTLs conferring seedling resistance to two rusts were mapped: five on chromosome 2B, three on 7B, two each on 5B and 6A, and one each on 1B, 2A, and 7A. Most of the QTLs identified were specific for a single rust species or race of a species. Candidate genes were identified within the confidence intervals of a QTL conferring resistance against at least two rust species by using the annotations of the durum (cv. 'Svevo') and wild emmer wheat ('Zavitan') reference genomes. The 22 identified loci conferring resistance to two or three rust species may be useful for breeding new and potentially durable resistant wheat cultivars.

Ridge regression and deep learning models for genome-wide selection of complex traits in New Mexican Chile peppers

BACKGROUND

Genomewide prediction estimates the genomic breeding values of selection candidates which can be utilized for population improvement and cultivar development. Ridge regression and deep learning-based selection models were implemented for yield and agronomic traits of 204 chile pepper genotypes evaluated in multi-environment trials in New Mexico, USA.

RESULTS

Accuracy of prediction differed across different models under ten-fold cross-validations, where high prediction accuracy was observed for highly heritable traits such as plant height and plant width. No model was superior across traits using 14,922 SNP markers for genomewide selection. Bayesian ridge regression had the highest average accuracy for first pod date (0.77) and total yield per plant (0.33). Multilayer perceptron (MLP) was the most superior for flowering time (0.76) and plant height (0.73), whereas the genomic BLUP model had the highest accuracy for plant width (0.62). Using a subset of 7,690 SNP loci resulting from grouping markers based on linkage disequilibrium coefficients resulted in improved accuracy for first pod date, ten pod weight, and total yield per plant, even under a relatively small training population size for MLP and random forest models. Genomic and ridge regression BLUP models were sufficient for optimal prediction accuracies for small training population size. Combining phenotypic selection and genomewide selection resulted in improved selection response for yield-related traits, indicating that integrated approaches can result in improved gains achieved through selection.

CONCLUSIONS

Accuracy values for ridge regression and deep learning prediction models demonstrate the potential of implementing genomewide selection for genetic improvement in chile pepper breeding programs. Ultimately, a large training data is relevant for improved genomic selection accuracy for the deep learning models.

Association analysis of production traits of Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing

This study was designed to perform an association analysis and identify SNP markers associated with production traits of Japanese quail using restriction-site-associated DNA sequencing. Weekly body weight data from 805 quail were collected from hatching to 16 weeks of age. A total number of 3990 eggs obtained from 399 female quail were used to assess egg quality traits. Egg-related traits were measured at the beginning of egg production (first stage) and at 12 weeks of age (second stage). Five eggs were analyzed at each stage. Traits, such as egg weight, egg length and short axes, eggshell strength and weight, egg equator thickness, yolk weight, diameter, and colour, albumen weight, age of first egg, total number of laid eggs, and egg production rate, were assessed. A total of 383 SNPs and 1151 associations as well as 734 SNPs and 1442 associations were identified in relation to quail production traits using general linear model (GLM) and mixed linear model (MLM) approaches, respectively. The GLM-identified SNPs were located on chromosomes 1-13, 15, 17-20, 24, 26-28, and Z, underlying phenotypic traits, except for egg and albumen weight at the first stage and yolk yellowness at the second stage. The MLM-identified SNPs were positioned on defined chromosomes associated with phenotypic traits except for the egg long axis at the second stage of egg production. Finally, 35 speculated genes were identified as candidate genes for the targeted traits based on their nearest positions. Our findings provide a deeper understanding and allow a more precise genetic improvement of production traits of Galliformes, particularly in Japanese quail.

Genome-wide association study of soluble solids content, flesh color, and fruit shape in citron watermelon

Fruit quality traits are crucial determinants of consumers' willingness to purchase watermelon produce, making them major goals for breeding programs. There is limited information on the genetic underpinnings of fruit quality traits in watermelon. A total of 125 citron watermelon (Citrullus amarus) accessions were genotyped using single nucleotide polymorphisms (SNPs) molecular markers generated via whole-genome resequencing. A total of 2,126,759 genome-wide SNP markers were used to uncover marker-trait associations using single and multi-locus GWAS models. High broad-sense heritability for fruit quality traits was detected. Correlation analysis among traits revealed positive relationships, with the exception of fruit diameter and fruit shape index (ratio of fruit length to fruit diameter), which was negative. A total of 37 significant SNP markers associated with soluble solids content, flesh color, fruit length, fruit diameter, and fruit shape index traits were uncovered. These peak SNPs accounted for 2.1%-23.4% of the phenotypic variation explained showing the quantitative inheritance nature of the evaluated traits. Candidate genes relevant to fruit quality traits were uncovered on chromosomes Ca01, Ca03, Ca06, and Ca07. These significant molecular markers and candidate genes will be useful in marker-assisted breeding of fruit quality traits in watermelon.

Multi-environment association study highlights candidate genes for robust agronomic quantitative trait loci in a novel worldwide Capsicum core collection

Investigating crop diversity through genome-wide association studies (GWAS) on core collections helps in deciphering the genetic determinants of complex quantitative traits. Using the G2P-SOL project world collection of 10 038 wild and cultivated Capsicum accessions from 10 major genebanks, we assembled a core collection of 423 accessions representing the known genetic diversity. Since complex traits are often highly dependent upon environmental variables and genotype-by-environment (G × E) interactions, multi-environment GWAS with a 10 195-marker genotypic matrix were conducted on a highly diverse subset of 350 Capsicum annuum accessions, extensively phenotyped in up to six independent trials from five climatically differing countries. Environment-specific and multi-environment quantitative trait loci (QTLs) were detected for 23 diverse agronomic traits. We identified 97 candidate genes potentially implicated in 53 of the most robust and high-confidence QTLs for fruit flavor, color, size, and shape traits, and for plant productivity, vigor, and earliness traits. Investigating the genetic architecture of agronomic traits in this way will assist the development of genetic markers and pave the way for marker-assisted selection. The G2P-SOL pepper core collection will be available upon request as a unique and universal resource for further exploitation in future gene discovery and marker-assisted breeding efforts by the pepper community.

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.

Cataloging SCN resistance loci in North American public soybean breeding programs

Soybean cyst nematode (SCN) is a destructive pathogen of soybeans responsible for annual yield loss exceeding $1.5 billion in the United States. Here, we conducted a series of genome-wide association studies (GWASs) to understand the genetic landscape of SCN resistance in the University of Missouri soybean breeding programs (Missouri panel), as well as germplasm and cultivars within the United States Department of Agriculture (USDA) Uniform Soybean Tests-Northern Region (NUST). For the Missouri panel, we evaluated the resistance of breeding lines to SCN populations HG 2.5.7 (Race 1), HG 1.2.5.7 (Race 2), HG 0 (Race 3), HG 2.5.7 (Race 5), and HG 1.3.6.7 (Race 14) and identified seven quantitative trait nucleotides (QTNs) associated with SCN resistance on chromosomes 2, 8, 11, 14, 17, and 18. Additionally, we evaluated breeding lines in the NUST panel for resistance to SCN populations HG 2.5.7 (Race 1) and HG 0 (Race 3), and we found three SCN resistance-associated QTNs on chromosomes 7 and 18. Through these analyses, we were able to decipher the impact of seven major genetic loci, including three novel loci, on resistance to several SCN populations and identified candidate genes within each locus. Further, we identified favorable allelic combinations for resistance to individual SCN HG types and provided a list of available germplasm for integration of these unique alleles into soybean breeding programs. Overall, this study offers valuable insight into the landscape of SCN resistance loci in U.S. public soybean breeding programs and provides a framework to develop new and improved soybean cultivars with diverse plant genetic modes of SCN resistance.

Genome-Wide Association Study Reveals Novel Powdery Mildew Resistance Loci in Bread Wheat

Powdery mildew (PM), caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt), significantly threatens global bread wheat production. Although the use of resistant cultivars is an effective strategy for managing PM, currently available wheat cultivars lack sufficient levels of resistance. To tackle this challenge, we conducted a comprehensive genome-wide association study (GWAS) using a diverse panel of 286 bread wheat genotypes. Over three consecutive years (2020-2021, 2021-2022, and 2022-2023), these genotypes were extensively evaluated for PM severity under field conditions following inoculation with virulent Bgt isolates. The panel was previously genotyped using the Illumina 90K Infinium iSelect assay to obtain genome-wide single-nucleotide polymorphism (SNP) marker coverage. By applying FarmCPU, a multilocus mixed model, we identified a total of 113 marker-trait associations (MTAs) located on chromosomes 1A, 1B, 2B, 3A, 3B, 4A, 4B, 5A, 5B, 6B, 7A, and 7B at a significance level of p ≤ 0.001. Notably, four novel MTAs on chromosome 6B were consistently detected in 2020-2021 and 2021-2022. Furthermore, within the confidence intervals of the identified SNPs, we identified 96 candidate genes belonging to different proteins including 12 disease resistance/host-pathogen interaction-related protein families. Among these, protein kinases, leucine-rich repeats, and zinc finger proteins were of particular interest due to their potential roles in PM resistance. These identified loci can serve as targets for breeding programs aimed at developing disease-resistant wheat cultivars.

Genome wide association in Spanish bread wheat landraces identifies six key genomic regions that constitute potential targets for improving grain yield related traits

KEY MESSAGE

Association mapping conducted in 189 Spanish bread wheat landraces revealed six key genomic regions that constitute stable QTLs for yield and include 15 candidate genes. Genetically diverse landraces provide an ideal population to conduct association analysis. In this study, association mapping was conducted in a collection of 189 Spanish bread wheat landraces whose genomic diversity had been previously assessed. These genomic data were combined with characterization for yield-related traits, including grain size and shape, and phenological traits screened across five seasons. The association analysis revealed a total of 881 significant marker trait associations, involving 434 markers across the genome, that could be grouped in 366 QTLs based on linkage disequilibrium. After accounting for days to heading, we defined 33 high density QTL genomic regions associated to at least four traits. Considering the importance of detecting stable QTLs, 6 regions associated to several grain traits and thousand kernel weight in at least three environments were selected as the most promising ones to harbour targets for breeding. To dissect the genetic cause of the observed associations, we studied the function and in silico expression of the 413 genes located inside these six regions. This identified 15 candidate genes that provide a starting point for future analysis aimed at the identification and validation of wheat yield related genes.

Genome-wide association study of powdery mildew resistance in cultivated soybean from Northeast China

Powdery mildew (PMD), caused by the pathogen Microsphaera diffusa, leads to substantial yield decreases in susceptible soybean under favorable environmental conditions. Effective prevention of soybean PMD damage can be achieved by identifying resistance genes and developing resistant cultivars. In this study, we genotyped 331 soybean germplasm accessions, primarily from Northeast China, using the SoySNP50K BeadChip, and evaluated their resistance to PMD in a greenhouse setting. To identify marker-trait associations while effectively controlling for population structure, we conducted genome-wide association studies utilizing factored spectrally transformed linear mixed models, mixed linear models, efficient mixed-model association eXpedited, and compressed mixed linear models. The results revealed seven single nucleotide polymorphism (SNP) loci strongly associated with PMD resistance in soybean. Among these, one SNP was localized on chromosome (Chr) 14, and six SNPs with low linkage disequilibrium were localized near or in the region of previously mapped genes on Chr 16. In the reference genome of Williams82, we discovered 96 genes within the candidate region, including 17 resistance (R)-like genes, which were identified as potential candidate genes for PMD resistance. In addition, we performed quantitative real-time reverse transcription polymerase chain reaction analysis to evaluate the gene expression levels in highly resistant and susceptible genotypes, focusing on leaf tissues collected at different times after M. diffusa inoculation. Among the examined genes, three R-like genes, including Glyma.16G210800, Glyma.16G212300, and Glyma.16G213900, were identified as strong candidates associated with PMD resistance. This discovery can significantly enhance our understanding of soybean resistance to PMD. Furthermore, the significant SNPs strongly associated with resistance can serve as valuable markers for genetic improvement in breeding M. diffusa-resistant soybean cultivars.

Integrating genome-wide association study into genomic selection for the prediction of agronomic traits in rice (Oryza sativa L.)

Accurately identifying varieties with targeted agronomic traits was thought to contribute to genetic selection and accelerate rice breeding progress. Genomic selection (GS) is a promising technique that uses markers covering the whole genome to predict the genomic-estimated breeding values (GEBV), with the ability to select before phenotypes are measured. To choose the appropriate GS models for breeding work, we analyzed the predictability of nine agronomic traits measured from a population of 459 diverse rice varieties. By the comparison of eight representative GS models, we found that the prediction accuracies ranged from 0.407 to 0.896, with reproducing kernel Hilbert space (RKHS) having the highest predictive ability in most traits. Further results demonstrated the predictivity of GS is altered by several factors. Moreover, we assessed the method of integrating genome-wide association study (GWAS) into various GS models. The predictabilities of GS combined peak-associated markers generated from six different GWAS models were significantly different; a recommendation of Mixed Linear Model (MLM)-RKHS was given for the GWAS-GS-integrated prediction. Finally, based on the above result, we experimented with applying the P-values obtained from optimal GWAS models into ridge regression best linear unbiased prediction (rrBLUP), which benefited the low predictive traits in rice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01423-y.

Exploring genetic diversity and ascertaining genetic loci associated with important fruit quality traits in apple (Malus × domestica Borkh.)

Genetic diversity is the primary source of variability in any crop improvement program, and the diverse germplasm of any crop species represents an important genetic resource for gene or allele mining to meet future needs. Huge genetic and phenotypic diversity is present in the apple gene pool, even though, breeding programs have been mainly focused on a few traits of interests, which have resulted in the reduction of the diversity in the cultivated lines of apple. Therefore, the present study was carried out on 70 diverse apple genotypes with the objective of analyzing the genetic diversity and to identify the genetic loci associated with important fruit quality traits. A total of 140 SSR primers were used to characterize the 70 genotypes of apples, out of which only 88 SSRs were found to be polymorphic. The PIC values varied from 0.03 to 0.75. The value of MI, EMR, and RP varied from 0.03 to 3.5, 0.5 to 5.0, and 1.89 to 6.74, respectively. The dendrogram and structure analysis divided all the genotypes into two main groups. In addition to this, large phenotypic variability was observed for the fruit quality traits under study indicated the suitability of the genotypes for association studies. Altogether 71 novel MTAs were identified for 10 fruit quality traits, of which 15 for fruit length, 15 for fruit diameter, 12 for fruit weight, 2 for total sugar, 2 for TSS, 4 for reducing sugar, 5 for non-reducing sugar, 5 for fruit firmness, 5 for fruit acidity and 6 for anthocyanin, respectively. Consistent with the physicochemical evaluation of traits, there was a significant correlation coefficient among different fruit quality characters, and many common markers were found to be associated with these traits (fruit diameter, length, TSS, total sugar, acidity and anthocyanin, respectively) by using the different modeling techniques (GLM, MLM). The inferred genetic structure, diversity pattern and the identified MTAs will be serving as resourceful grounds for better predictions and understanding of apple genome towards efficient conservation and utilization of apple germplasm for facilitating genetic improvement of fruit quality traits. Furthermore, these findings also suggested that association mapping could be a viable alternative to the conventional QTL mapping approach in apple.

Graphic abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01382-w.

Genetic control of flowering time and fruit yield in citron watermelon

Flowering time and fruit yield are important traits in watermelon crop improvement. There is limited information on the inheritance and genomic loci underlying flowering time and yield performance, especially in citron watermelon. A total of 125 citron watermelon accessions were evaluated in field trials over two growing seasons for days to male and female flowers, fruit count, fruit weight, and fruit yield. The germplasm was genotyped with more than two million single-nucleotide polymorphism (SNP) markers generated via whole-genome resequencing. Trait mapping was conducted using a genome-wide association study (GWAS). Broad-sense heritability for all traits ranged from moderate to high, indicating that genetic improvement through breeding and selection is feasible. Significant marker-trait associations were uncovered for days to female flower (chromosomes Ca04, Ca05, Ca08, and Ca09), fruit count (on Ca02, Ca03, and Ca05), fruit weight (on Ca02, Ca06, Ca08, Ca10, and Ca11), and fruit yield on chromosomes Ca05, Ca07, and Ca09. The phenotypic variation explained by the significant SNPs ranged from 1.6 to 25.4, highlighting the complex genetic architecture of the evaluated traits. Candidate genes relevant to flowering time and fruit yield component traits were uncovered on chromosomes Ca02, Ca04, Ca05, Ca06, Ca09, and Ca11. These results lay a foundation for marker-assisted trait introgression of flowering time and fruit yield component traits in watermelons.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Genome-wide association mapping in a sweet cherry germplasm collection (Prunus avium L.) reveals candidate genes for fruit quality traits

In sweet cherry (Prunus avium L.), large variability exists for various traits related to fruit quality. There is a need to discover the genetic architecture of these traits in order to enhance the efficiency of breeding strategies for consumer and producer demands. With this objective, a germplasm collection consisting of 116 sweet cherry accessions was evaluated for 23 agronomic fruit quality traits over 2-6 years, and characterized using a genotyping-by-sequencing approach. The SNP coverage collected was used to conduct a genome-wide association study using two multilocus models and three reference genomes. We identified numerous SNP-trait associations for global fruit size (weight, width, and thickness), fruit cracking, fruit firmness, and stone size, and we pinpointed several candidate genes involved in phytohormone, calcium, and cell wall metabolisms. Finally, we conducted a precise literature review focusing on the genetic architecture of fruit quality traits in sweet cherry to compare our results with potential colocalizations of marker-trait associations. This study brings new knowledge of the genetic control of important agronomic traits related to fruit quality, and to the development of marker-assisted selection strategies targeted towards the facilitation of breeding efforts.

Rethinking underutilized cereal crops: pan-omics integration and green system biology

MAIN CONCLUSION

Due to harsh lifestyle changes, in the present era, nutritional security is needed along with food security so it is necessary to include underutilized cereal crops (UCCs) in our daily diet to counteract the rising danger of human metabolic illness. We can attain both the goal of zero hunger and nutritional security by developing improved UCCs using advanced pan-omics (genomics, transcriptomics, proteomics, metabolomics, nutrigenomics, phenomics and ionomics) practices. Plant sciences research progressed profoundly since the last few decades with the introduction of advanced technologies and approaches, addressing issues of food demand of the growing population, nutritional security challenges and climate change. However, throughout the expansion and popularization of commonly consumed major cereal crops such as wheat and rice, other cereal crops such as millet, rye, sorghum, and others were impeded, despite their potential medicinal and nutraceutical qualities. Undoubtedly neglected underutilized cereal crops (UCCs) also have the capability to withstand diverse climate change. To relieve the burden of major crops, it is necessary to introduce the new crops in our diet in the way of UCCs. Introgression of agronomically and nutritionally important traits by pan-omics approaches in UCCs could be a defining moment for the population's well-being on the globe. This review discusses the importance of underutilized cereal crops, as well as the application of contemporary omics techniques and advanced bioinformatics tools that could open up new avenues for future study and be valuable assets in the development and usage of UCCs in the perspective of green system biology. The increased and improved use of UCCs is dependent on number of factors that necessitate a concerted research effort in agricultural sciences. The emergence of functional genomics with molecular genetics might gear toward the reawakening of interest in underutilized cereals crops. The need of this era is to focus on potential UCCs in advanced agriculture and breeding programmes. Hence, targeting the UCCs, might provide a bright future for better health and scientific rationale for its use.

Genome-Wide Association Mapping of Processing Quality Traits in Common Wheat (Triticum aestivum L.)

Processing quality is an important economic wheat trait. The marker-assisted selection (MAS) method plays a vital role in accelerating genetic improvement of processing quality. In the present study, processing quality in a panel of 165 cultivars grown in four environments was evaluated by mixograph. An association mapping analysis using 90 K and 660 K single nucleotide polymorphism (SNP) arrays identified 24 loci in chromosomes 1A, 1B (4), 1D, 2A, 2B (2), 3A, 3B, 3D (2), 4A (3), 4B, 5D (2), 6A, 7B (2) and 7D (2), explaining 10.2-42.5% of the phenotypic variances. Totally, 15 loci were stably detected in two or more environments. Nine loci coincided with known genes or QTL, whereas the other fifteen were novel loci. Seven candidate genes encoded 3-ketoacyl-CoA synthase, lipoxygenase, pyridoxal phosphate-dependent decarboxylase, sucrose synthase 3 and a plant lipid transfer protein/Par allergen. SNPs significantly associated with processing quality and accessions with more favorable alleles can be used for marker-assisted selection.

Genome-wide association studies: an intuitive solution for SNP identification and gene mapping in trees

Analysis of natural diversity in wild/cultivated plants can be used to understand the genetic basis for plant breeding programs. Recent advancements in DNA sequencing have expanded the possibilities for genetically altering essential features. There have been several recently disclosed statistical genetic methods for discovering the genes impacting target qualities. One of these useful methods is the genome-wide association study (GWAS), which effectively identifies candidate genes for a variety of plant properties by examining the relationship between a molecular marker (such as SNP) and a target trait. Conventional QTL mapping with highly structured populations has major limitations. The limited number of recombination events results in poor resolution for quantitative traits. Only two alleles at any given locus can be studied simultaneously. Conventional mapping approach fails to work in perennial plants and vegetatively propagated crops. These limitations are sidestepped by association mapping or GWAS. The flexibility of GWAS comes from the fact that the individuals being examined need not be linked to one another, allowing for the use of all meiotic and recombination events to increase resolution. Phenotyping, genotyping, population structure analysis, kinship analysis, and marker-trait association analysis are the fundamental phases of GWAS. With the rapid development of sequencing technologies and computational methods, GWAS is becoming a potent tool for identifying the natural variations that underlie complex characteristics in crops. The use of high-throughput sequencing technologies along with genotyping approaches like genotyping-by-sequencing (GBS) and restriction site associated DNA (RAD) sequencing may be highly useful in fast-forward mapping approach like GWAS. Breeders may use GWAS to quickly unravel the genomes through QTL and association mapping by taking advantage of natural variances. The drawbacks of conventional linkage mapping can be successfully overcome with the use of high-resolution mapping and the inclusion of multiple alleles in GWAS.

Identification of novel candidate loci and genes for seed vigor-related traits in upland cotton (Gossypium hirsutum L.) via GWAS

Seed vigor (SV) is a crucial trait determining the quality of crop seeds. Currently, over 80% of China's cotton-planting area is in Xinjiang Province, where a fully mechanized planting model is adopted, accounting for more than 90% of the total fiber production. Therefore, identifying SV-related loci and genes is crucial for improving cotton yield in Xinjiang. In this study, three seed vigor-related traits, including germination potential, germination rate, and germination index, were investigated across three environments in a panel of 355 diverse accessions based on 2,261,854 high-quality single-nucleotide polymorphisms (SNPs). A total of 26 significant SNPs were detected and divided into six quantitative trait locus regions, including 121 predicted candidate genes. By combining gene expression, gene annotation, and haplotype analysis, two novel candidate genes (Ghir_A09G002730 and Ghir_D03G009280) within qGR-A09-1 and qGI/GP/GR-D03-3 were associated with vigor-related traits, and Ghir_A09G002730 was found to be involved in artificial selection during cotton breeding by population genetic analysis. Thus, understanding the genetic mechanisms underlying seed vigor-related traits in cotton could help increase the efficiency of direct seeding by molecular marker-assisted selection breeding.

Green revolution to genome revolution: driving better resilient crops against environmental instability

Crop improvement programmes began with traditional breeding practices since the inception of agriculture. Farmers and plant breeders continue to use these strategies for crop improvement due to their broad application in modifying crop genetic compositions. Nonetheless, conventional breeding has significant downsides in regard to effort and time. Crop productivity seems to be hitting a plateau as a consequence of environmental issues and the scarcity of agricultural land. Therefore, continuous pursuit of advancement in crop improvement is essential. Recent technical innovations have resulted in a revolutionary shift in the pattern of breeding methods, leaning further towards molecular approaches. Among the promising approaches, marker-assisted selection, QTL mapping, omics-assisted breeding, genome-wide association studies and genome editing have lately gained prominence. Several governments have progressively relaxed their restrictions relating to genome editing. The present review highlights the evolutionary and revolutionary approaches that have been utilized for crop improvement in a bid to produce climate-resilient crops observing the consequence of climate change. Additionally, it will contribute to the comprehension of plant breeding succession so far. Investing in advanced sequencing technologies and bioinformatics will deepen our understanding of genetic variations and their functional implications, contributing to breakthroughs in crop improvement and biodiversity conservation.

Genome-wide association and targeted transcriptomic analyses reveal loci and candidate genes regulating preharvest sprouting in barley

KEY MESSAGE

Genome-wide association study of diverse barley genotypes identified loci, single nucleotide polymorphisms and candidate genes that control seed dormancy and therefore enhance resistance to preharvest sprouting. Preharvest sprouting (PHS) causes significant yield and quality loss in barley and it is strongly associated with the level of seed dormancy. This study performed genome-wide association study using a collection of 255 diverse barley genotypes grown over four environments to identify loci controlling dormancy/PHS. Our phenotypic analysis revealed substantial variation in germination index/dormancy levels among the barley genotypes. Marker-trait association and linkage disequilibrium (LD) decay analyses identified 16 single nucleotide polymorphisms (SNPs) and two QTLs associated with dormancy/PHS, respectively, on chromosome 3H and 5H explaining 6.9% to 11.1% of the phenotypic variation. QTL.5H consist of 14 SNPs of which 12 SNPs satisfy the FDR threshold of α = 0.05, and it may represent the SD2 locus. The QTL on 3H consists of one SNP that doesn't satisfy FDR (α = 0.05). Genes harbouring the significant SNPs were analyzed for their expression pattern in the seeds of selected dormant and non-dormant genotypes. Of these genes, HvRCD1, HvPSRP1 and HvF3H exhibited differential expression between the dormant and non-dormant seed samples, suggesting their role in controlling seed dormancy/PHS. Three SNPs located within the differentially expressed genes residing in QTL.5H explained considerable phenotypic variation (≥ 8.6%), suggesting their importance in regulating PHS resistance. Analysis of the SNP marker data in QTL.5H identified a haplotype for PHS resistance. Overall, the study identified loci, SNPs and candidate genes that control dormancy and therefore play important roles in enhancing PHS resistance in barley through marker-assisted breeding.

A physical map of traits of agronomic importance based on potato and tomato genome sequences

Potato, tomato, pepper, and eggplant are worldwide important crop and vegetable species of the Solanaceae family. Molecular linkage maps of these plants have been constructed and used to map qualitative and quantitative traits of agronomic importance. This research has been undertaken with the vision to identify the molecular basis of agronomic characters on the one hand, and on the other hand, to assist the selection of improved varieties in breeding programs by providing DNA-based markers that are diagnostic for specific agronomic characters. Since 2011, whole genome sequences of tomato and potato became available in public databases. They were used to combine the results of several hundred mapping and map-based cloning studies of phenotypic characters between 1988 and 2022 in physical maps of the twelve tomato and potato chromosomes. The traits evaluated were qualitative and quantitative resistance to pathogenic oomycetes, fungi, bacteria, viruses, nematodes, and insects. Furthermore, quantitative trait loci for yield and sugar content of tomato fruits and potato tubers and maturity or earliness were physically mapped. Cloned genes for pathogen resistance, a few genes underlying quantitative trait loci for yield, sugar content, and maturity, and several hundred candidate genes for these traits were included in the physical maps. The comparison between the physical chromosome maps revealed, in addition to known intrachromosomal inversions, several additional inversions and translocations between the otherwise highly collinear tomato and potato genomes. The integration of the positional information from independent mapping studies revealed the colocalization of qualitative and quantitative loci for resistance to different types of pathogens, called resistance hotspots, suggesting a similar molecular basis. Synteny between potato and tomato with respect to genomic positions of quantitative trait loci was frequently observed, indicating eventual similarity between the underlying genes.

Genetic diversity and genome-wide association study of 13 agronomic traits in 977 Beta vulgaris L. germplasms

BACKGROUND

Sugar beet (Beta vulgaris L.) is an economically essential sugar crop worldwide. Its agronomic traits are highly diverse and phenotypically plastic, influencing taproot yield and quality. The National Beet Medium-term Gene Bank in China maintains more than 1700 beet germplasms with diverse countries of origin. However, it lacks detailed genetic background associated with morphological variability and diversity.

RESULTS

Here, a comprehensive genome-wide association study (GWAS) of 13 agronomic traits was conducted in a panel of 977 sugar beet accessions. Almost all phenotypic traits exhibited wide genetic diversity and high coefficient of variation (CV). A total of 170,750 high-quality single-nucleotide polymorphisms (SNPs) were obtained using the genotyping-by-sequencing (GBS). Neighbour-joining phylogenetic analysis, principal component analysis, population structure and kinship showed no obvious relationships among these genotypes based on subgroups or regional sources. GWAS was carried out using a mixed linear model, and 159 significant associations were detected for these traits. Within the 25 kb linkage disequilibrium decay of the associated markers, NRT1/PTR FAMILY 6.3 (BVRB_5g097760); nudix hydrolase 15 (BVRB_8g182070) and TRANSPORT INHIBITOR RESPONSE 1 (BVRB_8g181550); transcription factor MYB77 (BVRB_2g023500); and ethylene-responsive transcription factor ERF014 (BVRB_1g000090) were predicted to be strongly associated with the taproot traits of root groove depth (RGD); root shape (RS); crown size (CS); and flesh colour (FC), respectively. For the aboveground traits, UDP-glycosyltransferase 79B6 (BVRB_9g223780) and NAC domain-containing protein 7 (BVRB_5g097990); F-box protein At1g10780 (BVRB_6g140760); phosphate transporter PHO1 (BVRB_3g048660); F-box protein CPR1 (BVRB_8g181140); and transcription factor MYB77 (BVRB_2g023500) and alcohol acyltransferase 9 (BVRB_2g023460) might be associated with the hypocotyl colour (HC); plant type (PT); petiole length (PL); cotyledon size (C); and fascicled leaf type (FLT) of sugar beet, respectively. AP-2 complex subunit mu (BVRB_5g106130), trihelix transcription factor ASIL2 (BVRB_2g041790) and late embryogenesis abundant protein 18 (BVRB_5g106150) might be involved in pollen quantity (PQ) variation. The candidate genes extensively participated in hormone response, nitrogen and phosphorus transportation, secondary metabolism, fertilization and embryo maturation.

CONCLUSIONS

The genetic basis of agronomical traits is complicated in heterozygous diploid sugar beet. The putative valuable genes found in this study will help further elucidate the molecular mechanism of each phenotypic trait for beet breeding.

Genome-wide association studies in rice germplasm reveal significant genomic regions for root and yield-related traits under aerobic and irrigated conditions

The development of nutrient-use efficient rice lines is a priority amidst the changing climate and depleting resources viz., water, land, and labor for achieving sustainability in rice cultivation. Along with the traditional transplanted irrigated system of cultivation, the dry direct-seeded aerobic system is gaining ground nationwide. The root-related traits play a crucial role in nutrient acquisition, adaptation and need to be concentrated along with the yield-attributing traits. We phenotyped an association panel of 118 rice lines for seedling vigour index (SVI) traits at 14 and 21 days after sowing (DAS), root-related traits at panicle initiation (PI) stage in polythene bags under controlled aerobic condition, yield and yield-related traits under the irrigated condition at ICAR-IIRR, Hyderabad, Telangana; irrigated and aerobic conditions at ARS, Dhadesugur, Raichur, Karnataka. The panel was genotyped using simple sequence repeats (SSR) markers and genome-wide association studies were conducted for identifying marker-trait associations (MTAs). Significant correlations were recorded for root length, root dry weight with SVI, root volume at the PI stage, number of productive tillers per plant, spikelet fertility, the total number of grains per panicle with grain yield per plant under irrigated conditions, and the total number of grains per panicle with grain yield per plant under aerobic condition. The panel was divided into three sub-groups (K = 3) and correlated with the principal component analysis. The maximum number of MTAs were found on chromosomes 2, 3, and 12 with considerable phenotypic variability. Consistent MTAs were recorded for SVI traits at 14 and 21 DAS (RM25310, RM80, RM22961, RM1385), yield traits under irrigated conditions (RM2584, RM5179, RM410, RM20698, RM14753) across years at ICAR-IIRR, grain yield per plant (RM22961, RM1146) under the aerobic condition, grain yield per plant at irrigated ICAR-IIRR and SVI (RM5501), root traits at PI stage (RM2584, RM80, RM410, RM1146, RM18472). Functionally relevant genes near the MTAs through in-silico expression analysis in root and panicle tissues viz., HBF2 bZIP transcription factor, WD40 repeat-like domain, OsPILS6a auxin efflux carrier, WRKY108, OsSCP42, OsMADS80, nodulin-like domain-containing protein, amino acid transporter using various rice expression databases were identified. The identified MTAs and rice lines having high SVI traits (Langphou, TI-128, Mouli, TI-124, JBB-631-1), high yield under aerobic (Phouren, NPK-43, JBB-684, Ratnamudi, TI-112), irrigated conditions (KR-209, KR-262, Phouren, Keibi-Phou, TI-17), robust root traits like root length (MoirangPhou-Angouba, Wangoo-Phou, JBB-661, Dissi, NPK-45), root volume (Ratnachudi, KJ-221, Mow, Heimang-Phou, PUP-229) can be further employed in breeding programs for the targeted environments aimed at improving seedling vigour, yield-related traits under irrigated condition, aerobic condition as adaptability to water-saving technology.

Genetic diversity and genome-wide association analysis of pine wood nematode populations in different regions of China

Introduction

Pine wilt disease (Bursaphelenchus xylophilus) was recently detected in Liaoning Province, which was previously considered an unfavourable area for B. xylophilus due to its low temperatures. This study aims to compare the reproductivity and genetic variations of B. xylophilus isolates from Liaoning Province and other parts of China to explore their phenotypic and genomic differences.

Methods

The samples from Liaoning, Anhui, Hubei, Henan, Zhejiang and Jiangsu were isolated and purified to obtain the strains. The reproductivity of the strains was determined at 15 °C. The genetic structure was analyzed by using SNP molecular markers, and the whole genome association analysis was carried out by integrating SNP information and feculence traits.

Results

A reproductivity experiment showed that Liaoning isolates have higher reproductive ability at 15 °C. Subsequent SNP profiling and population differentiation analysis revealed obvious genetic differentiation of Liaoning isolates from other isolates. A genome-wide association study showed that SNPs closely related to low-temperature tolerance were mainly located in GPCR, Acyl-CoA, and Cpn10, which are responsible for adaptation to environmental factors, such as temperature change.

Discussion

Pine wood nematodes likely adapted to the climate in Liaoning and maintained a certain reproductive capacity at low temperature via variants of adaptation-related genes. This study provides a theoretical basis for elucidating the prevalence and diffusion status of B. xylophilus in China.

Identification of Quantitative Trait Loci Associated with Powdery Mildew Resistance in Spring Barley under Conditions of Southeastern Kazakhstan

Barley (Hordeum vulgare L.) is one of the most produced cereal crops in the world. It has traditionally been used for the production of animal feed and for malting, as well as for human consumption. However, its production is highly affected by biotic stress factors, particularly the fungal pathogen Blumeria graminis (DC.) f. sp. hordei (Bgh), which causes powdery mildew (PM). In this study, a collection of 406 barley accessions from the USA, Kazakhstan, Europe, and Africa were assessed for resistance to PM over a 3-year period in southeastern Kazakhstan. The collection was grown in the field in 2020, 2021, and 2022 and was genotyped using the 9K SNP Illumina chip. A genome-wide association study (GWAS) was conducted to identify the quantitative trait loci (QTLs) associated with PM resistance. As a result, seven QTLs for PM resistance were detected on chromosomes 4H, 5H, and 7H (FDR p-values < 0.05). Genetic positions of two QTLs were similar to those of PM resistance QTLs previously reported in the scientific literature, suggesting that the five remaining QTLs are novel putative genetic factors for the studied trait. Haplotype analysis for seven QTLs revealed three haplotypes which were associated with total PM resistance and one haplotype associated with the high PM severity in the barley collection. Identified QTLs and haplotypes associated with the PM resistance of barley may be used for further analysis, trait pyramiding, and marker-assisted selection.