Genome-wide association study identifies favorable SNP alleles and candidate genes for waterlogging tolerance in chrysanthemums

Identification of SNP Markers and Candidate Genes Associated with Major Agronomic Traits in Coffea arabica

Genome-wide association studies (GWASs) allow for inferences about the relationships between genomic variants and phenotypic traits in natural or breeding populations. However, few have used this methodology in Coffea arabica. We aimed to identify chromosomal regions with significant associations between SNP markers and agronomic traits in C. arabica. We used a coffee panel consisting of 195 plants derived from 13 families in F2 generations and backcrosses of crosses between leaf rust-susceptible and -resistant genotypes. The plants were phenotyped for 18 agronomic markers and genotyped for 21,211 SNP markers. A GWAS enabled the identification of 110 SNPs with significant associations (p < 0.05) for several agronomic traits in C. arabica: plant height, plagiotropic branch length, number of vegetative nodes, canopy diameter, fruit size, cercosporiosis incidence, and rust incidence. The effects of each SNP marker associated with the traits were analyzed, such that they can be used for molecular marker-assisted selection. For the first time, a GWAS was used for these important agronomic traits in C. arabica, enabling applications in accelerated coffee breeding through marker-assisted selection and ensuring greater efficiency and time reduction. Furthermore, our findings provide preliminary knowledge to further confirm the genomic loci and potential candidate genes contributing to various structural and disease-related traits of C. arabica.

Multi-locus genome-wide association studies reveal the dynamic genetic architecture of flowering time in chrysanthemum

KEY MESSAGE

The dynamic genetic architecture of flowering time in chrysanthemum was elucidated by GWAS. Thirty-six known genes and 14 candidate genes were identified around the stable QTNs and QEIs, among which ERF-1 was highlighted. Flowering time (FT) adaptation is one of the major breeding goals in chrysanthemum, a multipurpose ornamental plant. In order to reveal the dynamic genetic architecture of FT in chrysanthemum, phenotype investigation of ten FT-related traits was conducted on 169 entries in 2 environments. The broad-sense heritability of five non-conditional FT traits, i.e., budding (FBD), visible coloring (VC), early opening (EO), full-bloom (OF) and decay period (DP), ranged from 56.93 to 84.26%, which were higher than that of the five derived conditional FT traits (38.51-75.13%). The phenotypic variation coefficients of OF_EO and DP_OF were relatively large ranging from 30.59 to 36.17%. Based on 375,865 SNPs, the compressed variance component mixed linear model 3VmrMLM was applied for a multi-locus genome-wide association study (GWAS). As a result, 313 quantitative trait nucleotides (QTNs) were identified for the non-conditional FT traits in single-environment analysis, while 119 QTNs and 67 QTN-by-environment interactions (QEIs) were identified in multi-environment analysis. As for the conditional traits, 343 QTNs were detected in single-environment analysis, and 119 QTNs and 83 QEIs were identified in multi- environment analysis. Among the genes around stable QTNs and QEIs, 36 were orthologs of known FT genes in Arabidopsis and other plants; 14 candidates were mined by combining the transcriptomics data and functional annotation, including ERF-1, ACA10, and FOP1. Furthermore, the haplotype analysis of ERF-1 revealed six elite accessions with extreme FBD. Our findings contribute to the understanding of dynamic genetic architecture of FT and provide valuable resources for future chrysanthemum molecular breeding programs.

Genetic architecture and genomic prediction of plant height-related traits in chrysanthemum

Plant height (PH) is a crucial trait determining plant architecture in chrysanthemum. To better understand the genetic basis of PH, we investigated the variations of PH, internode number (IN), internode length (IL), and stem diameter (SD) in a panel of 200 cut chrysanthemum accessions. Based on 330 710 high-quality SNPs generated by genotyping by sequencing, a total of 42 associations were identified via a genome-wide association study (GWAS), and 16 genomic regions covering 2.57 Mb of the whole genome were detected through selective sweep analysis. In addition, two SNPs, Chr1_339370594 and Chr18_230810045, respectively associated with PH and SD, overlapped with the selective sweep regions from FST and π ratios. Moreover, candidate genes involved in hormones, growth, transcriptional regulation, and metabolic processes were highlighted based on the annotation of homologous genes in Arabidopsis and transcriptomes in chrysanthemum. Finally, genomic selection for four PH-related traits was performed using a ridge regression best linear unbiased predictor model (rrBLUP) and six marker sets. The marker set constituting the top 1000 most significant SNPs identified via GWAS showed higher predictabilities for the four PH-related traits, ranging from 0.94 to 0.97. These findings improve our knowledge of the genetic basis of PH and provide valuable markers that could be applied in chrysanthemum genomic selection breeding programs.

Identification of novel putative alleles related to important agronomic traits of wheat using robust strategies in GWAS

Principal component analysis (PCA) is widely used in various genetics studies. In this study, the role of classical PCA (cPCA) and robust PCA (rPCA) was evaluated explicitly in genome-wide association studies (GWAS). We evaluated 294 wheat genotypes under well-watered and rain-fed, focusing on spike traits. First, we showed that some phenotypic and genotypic observations could be outliers based on cPCA and different rPCA algorithms (Proj, Grid, Hubert, and Locantore). Hubert's method provided a better approach to identifying outliers, which helped to understand the nature of these samples. These outliers led to the deviation of the heritability of traits from the actual value. Then, we performed GWAS with 36,000 single nucleotide polymorphisms (SNPs) based on the traditional approach and two robust strategies. In the conventional approach and using the first three components of cPCA as population structure, 184 and 139 marker-trait associations (MTAs) were identified for five traits in well-watered and rain-fed environments, respectively. In the first robust strategy and when rPCA was used as population structure in GWAS, we observed that the Hubert and Grid methods identified new MTAs, especially for yield and spike weight on chromosomes 7A and 6B. In the second strategy, we followed the classical and robust principal component-based GWAS, where the first two PCs obtained from phenotypic variables were used instead of traits. In the recent strategy, despite the similarity between the methods, some new MTAs were identified that can be considered pleiotropic. Hubert's method provided a better linear combination of traits because it had the most MTAs in common with the traditional approach. Newly identified SNPs, including rs19833 (5B) and rs48316 (2B), were annotated with important genes with vital biological processes and molecular functions. The approaches presented in this study can reduce the misleading GWAS results caused by the adverse effect of outlier observations.

Genome-Wide Association Mapping for Yield and Yield-Related Traits in Rice (Oryza Sativa L.) Using SNPs Markers

Rice (Oryza sativa L.) is a staple food for more than 50% of the world's population. Rice cultivar improvement is critical in order to feed the world's growing population. Improving yield is one of the main aims of rice breeders. However, yield is a complex quantitative trait controlled by many genes. The presence of genetic diversity is the key factor to improve the yield hence, the presence of diversity in any germplasm is important for yield improvement. In the current study, the rice germplasm was collected from Pakistan and the United States of America and a panel of 100 diverse genotypes was utilized to identify important yield and yield-related traits. For this, a genome-wide association study (GWAS) was performed to identify the genetic loci related to yield. The GWAS on the diverse germplasm will lead to the identification of new genes which can be utilized in the breeding program for improvement of yield. For this reason, firstly, the germplasm was phenotypically evaluated in two growing seasons for yield and yield-related traits. The analysis of variance results showed significant differences among traits which showed the presence of diversity in the current germplasm. Secondly, the germplasm was also genotypically evaluated using 10K SNP. Genetic structure analysis showed the presence of four groups which showed that enough genetic diversity was present in the rice germplasm to be used for association mapping analysis. The results of GWAS identified 201 significant marker trait associations (MTAs. 16 MTAs were identified for plant height, 49 for days to flowering, three for days to maturity, four for tillers per plant, four for panicle length, eight for grains per panicle, 20 unfilled grains per panicle, 81 for seed setting %, four for thousand-grain weight, five for yield per plot and seven for yield per hectare. Apart from this, some pleiotropic loci were also identified. The results showed that panicle length (PL) and thousand-grain weight (TGW) were controlled by a pleiotropic locus OsGRb23906 on chromosome 1 at 10,116,371 cM. The loci OsGRb25803 and OsGRb15974 on chromosomes 4 and 8 at the position of 14,321,111 cM and 6,205,816 cM respectively, showed pleiotropic effects for seed setting % (SS) and unfilled grain per panicle (UG/P). A locus OsGRb09180 on chromosome 4 at 19,850,601 cM was significantly linked with SS and yield/ha. Furthermore, gene annotation was performed, and results indicated that the 190 candidate genes or QTLs that closely linked with studied traits. These candidate genes and novel significant markers could be useful in marker-assisted gene selection and QTL pyramiding to improve rice yield and the selection of potential parents, recombinants and MTAs which could be used in rice breeding programs to develop high-yielding rice varieties for sustainable food security.

Whole-genome re-sequencing, diversity analysis, and stress-resistance analysis of 77 grape rootstock genotypes

INTRODUCTION

Grape rootstocks play critical role in the development of the grape industry over the globe for their higher adaptability to various environments, and the evaluation of their genetic diversity among grape genotypes is necessary to the conservation and utility of genotypes.

METHODS

To analyze the genetic diversity of grape rootstocks for a better understanding multiple resistance traits, whole-genome re-sequencing of 77 common grape rootstock germplasms was conducted in the present study.

RESULTS

About 645 billion genome sequencing data were generated from the 77 grape rootstocks at an average depth of ~15.5×, based on which the phylogenic clusters were generated and the domestication of grapevine rootstocks was explored. The results indicated that the 77 rootstocks originated from five ancestral components. Through phylogenetic, principal components, and identity-by-descent (IBD) analyses, these 77 grape rootstocks were assembled into ten groups. It is noticed that the wild resources of V. amurensis and V. davidii, originating from China and being generally considered to have stronger resistance against biotic and abiotic stresses, were sub-divided from the other populations. Further analysis indicated that a high level of linkage disequilibrium was found among the 77 rootstock genotypes, and a total of 2,805,889 single nucleotide polymorphisms (SNPs) were excavated, GWAS analysis among the grape rootstocks located 631, 13, 9, 2, 810, and 44 SNP loci that were responsible to resistances to phylloxera, root-knot nematodes, salt, drought, cold and waterlogging traits.

DISCUSSION

This study generated a significant amount of genomic data from grape rootstocks, thus providing a theoretical basis for further research on the resistance mechanism of grape rootstocks and the breeding of resistant varieties. These findings also reveal that China originated V. amurensis and V. davidii could broaden the genetic background of grapevine rootstocks and be important germplasm used in breeding high stress-resistant grapevine rootstocks.

Genome-wide association mapping of bread wheat genotypes using yield and grain morphology-related traits under different environments

Depleting water resources and increasing global temperature due to climate change are major challenges to agriculture and food security worldwide. Deciphering the underlying mechanisms of traits contributing to grain development and yield is essential for the development of climate-resilient cultivars. Therefore, this study assessed 105 bread wheat genotypes grown under control, drought, and heat-stress conditions for two crop seasons and performed a genome-wide association study (GWAS) using a 90k SNP array. The genotypes showed significant trait differences under all environmental conditions. Highly significant variation was observed, with moderate (50.09%) to high (76.19%) heritability in the studied germplasms. The studied traits were all also significantly positively correlated. A total of 541 significant associations (p ≤ 10^-3) between marker and trait (MTAs) were observed after crossing the FDR <0.05 threshold for all traits. Among these, 195, 179, and 167 significant MTAs were detected under control, drought, and heat-stress conditions, respectively. Under the control and drought conditions, pleiotropic loci BS00010616_51 and BS00010868_51 were observed on chromosomes 7B and 1B situated at 186.24 cM and 35.47 cM, respectively. Pleiotropic loci BS00010868_51, Kukri_c11154_1723, and Ex_c10068_1509 were identified on chromosomes 1B, 5B, and 2A, respectively, under control and heat stress conditions. A stable and consistent locus (Excalibur_c20796_395) on chromosome 7A, located at 372.34 cM, was also linked to grain morphology and yield-related attributes in control, drought, and heat-stress conditions. The results of the current study confirmed several previously reported MTAs for the traits under consideration and identified new MTAs under harsh climatic conditions. These SNPs will aid in the discovery of novel genes in wheat. SNPs showing significant associations may be used in marker-assisted selection and allow the development of drought- and heat-tolerant genotypes with high yields to address global food security concerns.

Bioinformatics in Plant Breeding and Research on Disease Resistance

In the context of plant breeding, bioinformatics can empower genetic and genomic selection to determine the optimal combination of genotypes that will produce a desired phenotype and help expedite the isolation of these new varieties. Bioinformatics is also instrumental in collecting and processing plant phenotypes, which facilitates plant breeding. Robots that use automated and digital technologies to collect and analyze different types of information to monitor the environment in which plants grow, analyze the environmental stresses they face, and promptly optimize suboptimal and adverse growth conditions accordingly, have helped plant research and saved human resources. In this paper, we describe the use of various bioinformatics databases and algorithms and explore their potential applications in plant breeding and for research on plant disease resistance.

The transcriptional coactivator CmMBF1c is required for waterlogging tolerance in Chrysanthemum morifolium

Waterlogging is one of the most serious abiotic stressors affecting Chrysanthemum morifolium during its lifespan. However, the molecular mechanisms underlying the waterlogging tolerance of chrysanthemum remain unclear. In this study, we discovered that the transcriptional coactivator MULTIPROTEIN BRIDGING FACTOR 1c (CmMBF1c) was significantly induced by waterlogging stress in chrysanthemums. Promoter sequence analysis and transient dual-luciferase assay using chrysanthemum protoplasts showed that the waterlogging-tolerant cultivar 'Nannongxuefeng' carried more response elements involved in waterlogging and hypoxia stress compared with the waterlogging-sensitive cultivar 'Qinglu', conferring on 'Nannongxuefeng' a stronger hypoxia responsive activity and higher CmMBF1c expression under waterlogging conditions. Subcellular localization and transcriptional activity assays showed that CmMBF1c protein was localized to the nucleus and had no transcriptional activation activity. Overexpression of CmMBF1c in 'Qinglu' enhanced its waterlogging tolerance by promoting its reactive oxygen species (ROS) scavenging ability and maintaining low ROS levels. However, RNAi-mediated knockdown of CmMBF1c in cultivar 'Nannongxuefeng' resulted in the opposite tendency. Yeast two-hybrid screening and tobacco bimolecular fluorescence complementation assays revealed that CmHRE2, a pivotal regulator of hypoxia response, could interact with CmMBF1c. In summary, this study demonstrates that CmMBF1c improves chrysanthemum waterlogging tolerance by regulating its ROS signaling pathway and interacting with CmHRE2. These findings together offer, to our knowledge, new mechanistic insights into chrysanthemum waterlogging tolerance and provide a rational foundation for future research on the genetic improvement of horticultural crops for waterlogging stress tolerance.

Genome-wide association study, haplotype analysis, and genomic prediction reveal the genetic basis of yield-related traits in soybean (Glycine max L.)

Identifying the genetic components underlying yield-related traits in soybean is crucial for improving its production and productivity. Here, 211 soybean genotypes were evaluated across six environments for four yield-related traits, including seed yield per plant (SYP), number of pods per plant number of seeds per plant and 100-seed weight (HSW). Genome-wide association study (GWAS) and genomic prediction (GP) analyses were performed using 12,617 single nucleotide polymorphism markers from NJAU 355K SoySNP Array. A total of 57 SNPs were significantly associated with four traits across six environments and a combined environment using five Genome-wide association study models. Out of these, six significant SNPs were consistently identified in more than three environments using multiple GWAS models. The genomic regions (±670 kb) flanking these six consistent SNPs were considered stable QTL regions. Gene annotation and in silico expression analysis revealed 15 putative genes underlying the stable QTLs that might regulate soybean yield. Haplotype analysis using six significant SNPs revealed various allelic combinations regulating diverse phenotypes for the studied traits. Furthermore, the GP analysis revealed that accurate breeding values for the studied soybean traits is attainable at an earlier generation. Our study paved the way for increasing soybean yield performance within a short breeding cycle.

Genome-Wide Association Studies of Root-Related Traits in Brassica napus L. under Low-Potassium Conditions

Roots are essential organs for a plant’s ability to absorb water and obtain mineral nutrients, hence they are critical to its development. Plants use root architectural alterations to improve their chances of absorbing nutrients when their supply is low. Nine root traits of a Brassica napus association panel were explored in hydroponic-system studies under low potassium (K) stress to unravel the genetic basis of root growth in rapeseed. The quantitative trait loci (QTL) and candidate genes for root development were discovered using a multilocus genome-wide association study (ML-GWAS). For the nine traits, a total of 453 significant associated single-nucleotide polymorphism (SNP) loci were discovered, which were then integrated into 206 QTL clusters. There were 45 pleiotropic clusters, and qRTA04-4 and qRTC04-7 were linked to TRL, TSA, and TRV at the same time, contributing 5.25–11.48% of the phenotypic variance explained (PVE) to the root traits. Additionally, 1360 annotated genes were discovered by examining genomic regions within 100 kb upstream and downstream of lead SNPs within the 45 loci. Thirty-five genes were identified as possibly regulating root-system development. As per protein–protein interaction analyses, homologs of three genes (BnaC08g29120D, BnaA07g10150D, and BnaC04g45700D) have been shown to influence root growth in earlier investigations. The QTL clusters and candidate genes identified in this work will help us better understand the genetics of root growth traits and could be employed in marker-assisted breeding for rapeseed adaptable to various conditions with low K levels.

Genetic Diversity and Genome-Wide Association Study of Architectural Traits of Spray Cut Chrysanthemum Varieties

The architecture of spray cut chrysanthemum is crucial for the quality and quantity of cut flower production. However, the mechanism underlying plant architecture still needs to be clarified. In this study, we measured nine architecture-related traits of 195 spray cut chrysanthemum varieties during a two-year period. The results showed that the number of upper primary branches, number of lateral flower buds and primary branch length widely varied. Additionally, plant height had a significant positive correlation with number of leaf nodes and total number of lateral buds. Number of upper primary branches had a significant negative correlation with primary branch diameter, primary branch angle and primary branch length. Plant height, total number of lateral buds, number of upper primary branches, stem diameter, primary branch diameter and primary branch length were vulnerable to environmental impacts. All varieties could be divided into five categories according to cluster analysis, and the typical plant architecture of the varieties was summarized. Finally, a genome-wide association study (GWAS) was performed to find potential functional genes.

Genome-wide association mapping of bread wheat genotypes for sustainable food security and yield potential under limited water conditions

Determining the genetic basis of yield and water deficient tolerance in wheat is vital for wheat breeding programs. Herein, a genome-wide association study (GWAS) was performed for water deficient and yield-related attributes on wheat genotypes with high-density Illumina 90K Infinium SNP array. Major yield and drought-related attributes were phenotyped on a panel of Pakistani and foreign accessions grown in non-stressed and water deficient stressed environments during two crop cycles. Among all accessions, highly significant variations were shown in studied environments for examined characters. Water deficient conditions, reduced the wheat yield and had strong and positive correlation among relative water content and grain yield per plant. Population structure analyses based on 90,000 SNP data, classify the accessions into 4 sub-populations. Marker-trait association analyses (MTA) revealed that 134 significant SNPs were linked with yield and drought tolerance attributes. Pleotropic loci RAC875_s117925_244 and RAC875_c16333_340 located on chromosome 5A and 2A respectively, were significantly linked with relative water contents (RWC), cell membrane thermo-stability (CMT), grain per spike (GPS), spikelet per spike (SPS) and grain yield per plant (GYP). The markers Ra_c58279_684, BobWhite_c23828_341 and IAAV3414 located on chromosomes 2A, 6B and 7B respectively, showed pleotropic effects for RWC, GPS and GYP under both environments. The current experiment not only validated several MTAs reported in other studies but also discovered novel MTAs which significant under drought-stressed conditions. A total of 171 candidate genes were recognized that could be cloned and functionally characterized for the respective associated traits. For RWC and CMT, total 11 and 3 associated SNPs were mapped on coding DNA sequence (CDS) of the identified candidate genes. Isolation and characterization of the candidate genes herein mapped SNPs will be useful in discovering novel genes underpinning drought tolerance in bread wheat to fulfill the wheat demand and sustainable food security under limited water conditions.

Genome-Wide Association Study Reveals Marker Trait Associations (MTA) for Waterlogging-Triggered Adventitious Roots and Aerenchyma Formation in Barley

Waterlogging is an environmental stress, which severely affects barley growth and development. Limited availability of oxygen in the root zone negatively affects the metabolism of the whole plant. Adventitious roots (AR) and root cortical aerenchyma (RCA) formation are the most important adaptive traits that contribute to a plant's ability to survive in waterlogged soil conditions. This study used a genome-wide association (GWAS) approach using 18,132 single nucleotide polymorphisms (SNPs) in a panel of 697 barley genotypes to reveal marker trait associations (MTA) conferring the above adaptive traits. Experiments were conducted over two consecutive years in tanks filled with soil and then validated in field experiments. GWAS analysis was conducted using general linear models (GLM), mixed linear models (MLM), and fixed and random model circulating probability unification models (FarmCPU model), with the FarmCPU showing to be the best suited model. Six and five significant (approximately -log₁₀ (p) ≥ 5.5) MTA were identified for AR and RCA formation under waterlogged conditions, respectively. The highest -log₁₀ (p) MTA for adventitious root and aerenchyma formation were approximately 9 and 8 on chromosome 2H and 4H, respectively. The combination of different MTA showed to be more effective in forming RCA and producing more AR under waterlogging stress. Genes from major facilitator superfamily (MFS) transporter and leucine-rich repeat (LRR) families for AR formation, and ethylene responsive factor (ERF) family genes and potassium transporter family genes for RCA formation were the potential candidate genes involved under waterlogging conditions. Several genotypes, which performed consistently well under different conditions, can be used in breeding programs to develop waterlogging-tolerant varieties.

Identification of Chilling Accumulation-Associated Genes for Litchi Flowering by Transcriptome-Based Genome-Wide Association Studies

Litchi is an important Sapindaceae fruit tree. Flowering in litchi is triggered by low temperatures in autumn and winter. It can be divided into early-, medium-, and late-flowering phenotypes according to the time for floral induction. Early-flowering varieties need low chilling accumulation level for floral induction, whereas the late-flowering varieties require high chilling accumulation level. In the present study, RNA-Seq of 87 accessions was performed and transcriptome-based genome-wide association studies (GWAS) was used to identify candidate genes involved in chilling accumulation underlying the time for floral induction. A total of 98,155 high-quality single-nucleotide polymorphism (SNP) sites were obtained. A total of 1,411 significantly associated SNPs and 1,115 associated genes (AGs) were identified, of which 31 were flowering-related, 23 were hormone synthesis-related, and 27 were hormone signal transduction-related. Association analysis between the gene expression of the AGs and the flowering phenotypic data was carried out, and differentially expressed genes (DEGs) in a temperature-controlled experiment were obtained. As a result, 15 flowering-related candidate AGs (CAGs), 13 hormone synthesis-related CAGs, and 11 hormone signal transduction-related CAGs were further screened. The expression levels of the CAGs in the early-flowering accessions were different from those in the late-flowering ones, and also between the flowering trees and non-flowering trees. In a gradient chilling treatment, flowering rates of the trees and the CAGs expression were affected by the treatment. Our present work for the first time provided candidate genes for genetic regulation of flowering in litchi using transcriptome-based GWAS.

Transcriptome Analysis Reveals Genes Respond to Chlorophyll Deficiency in Green and Yellow Leaves of Chrysanthemum morifolium Ramat

Chlorophyll is vital for photosynthesis to produce sugars and other useful biochemical products in green plants. However, the molecular effects of chlorophyll deficiency in Chrysanthemum are largely unknown. In this study, we identified a bud sport mutant chrysanthemum belonging to the variety ‘Nannong Binyun’, which has yellow branches. Plant physiological studies have shown that the yellow color is revealed due to chlorophyll loss. RNA extracts of yellow and green tissues were analyzed using high-throughput RNA-sequencing, and a total of 11,649 tissue enriched unigenes that respond to chlorophyll deficiency were identified, including 4803 unigenes upregulated in yellow tissues and 6846 unigenes in green tissues. GO analysis revealed that these tissue-enriched genes may involve in the physiological processes of chlorophyll accumulation and photosynthesis. In addition, many DEGs from the families of AP2-EREBP, bHLH, MYB, and FAR1 that are associated with plant development and stress response were detected. Our study found that most of the genes from the GRAS family were downregulated in yellow leaves, indicating their putative roles in stem cell maintenance and possible contribution to leaf size determination.

Detection of Reproducible Major Effect QTL for Petal Traits in Garden Roses

The detection of QTL by association genetics depends on the genetic architecture of the trait under study, the size and structure of the investigated population and the availability of phenotypic and marker data of sufficient quality and quantity. In roses, we previously demonstrated that major QTL could already be detected in small association panels. In this study, we analyzed petal number, petal size and fragrance in a small panel of 95 mostly tetraploid garden rose genotypes. After genotyping the panel with the 68 K Axiom WagRhSNP chip we detected major QTL for all three traits. Each trait was significantly influenced by several genomic regions. Some of the QTL span genomic regions that comprise several candidate genes. Selected markers from some of these regions were converted into KASP markers and were validated in independent populations of up to 282 garden rose genotypes. These markers demonstrate the robustness of the detected effects independent of the set of genotypes analyzed. Furthermore, the markers can serve as tools for marker-assisted breeding in garden roses. Over an extended timeframe, they may be used as a starting point for the isolation of the genes underlying the QTL.

The soybean (Glycine max L.) cytokinin oxidase/dehydrogenase multigene family; Identification of natural variations for altered cytokinin content and seed yield

Cytokinins (CKs) play a fundamental role in regulating dynamics of organ source/sink relationships during plant development, including flowering and seed formation stages. As a result, CKs are key drivers of seed yield. The cytokinin oxidase/dehydrogenase (CKX) is one of the critical enzymes responsible for regulating plant CK levels by causing their irreversible degradation. Variation of CKX activity is significantly correlated with seed yield in many crop species while in soybean (Glycine max L.), the possible associations between CKX gene family members (GFMs) and yield parameters have not yet been assessed. In this study, 17 GmCKX GFMs were identified, and natural variations among GmCKX genes were probed among soybean cultivars with varying yield characteristics. The key CKX genes responsible for regulating CK content during seed filling stages of reproductive development were highlighted using comparative phylogenetics, gene expression analysis and CK metabolite profiling. Five of the seventeen identified GmCKX GFMs, showed natural variations in the form of single nucleotide polymorphisms (SNPs). The gene GmCKX7-1, with high expression during critical seed filling stages, was found to have a non-synonymous mutation (H105Q), on one of the active site residues, Histidine 105, previously reported to be essential for co-factor binding to maintain structural integrity of the enzyme. Soybean lines with this mutation had higher CK content and desired yield characteristics. The potential for marker-assisted selection based on the identified natural variation within GmCKX7-1, is discussed in the context of hormonal control that can result in higher soybean yield.

A genome-wide association study uncovers a critical role of the RsPAP2 gene in red-skinned Raphanus sativus L

Radish (Raphanus sativus L.) taproot contains high concentrations of flavonoids, including anthocyanins (ATCs), in red-skinned genotypes. However, little information on the genetic regulation of ATC biosynthesis in radish is available. A genome-wide association study of radish red skin color was conducted using whole-genome sequencing data derived from 179 radish genotypes. The R2R3-MYB transcription factor production of anthocyanin pigment 2 (PAP2) gene was found in the region associated with a leading SNP located on chromosome 2. The amino acid sequence encoded by the RsPAP2 gene was different from those of the other published RsMYB genes responsible for the red skin color of radish. The overexpression of the RsPAP2 gene resulted in ATC accumulation in Arabidopsis and radish, which was accompanied by the upregulation of several ATC-related structural genes. RsPAP2 was found to bind the RsUFGT and RsTT8 promoters, as shown by a dual-luciferase reporter system and a yeast one-hybrid assay. The promoter activities of the RsANS, RsCHI, RsPAL, and RsUFGT genes could be strongly activated by coinfiltration with RsPAP2 and RsTT8. These findings showed the effectiveness of GWAS in identifying candidate genes in radish and demonstrated that RsPAP2 could (either directly or together with its cofactor RsTT8) regulate the transcript levels of ATC-related genes to promote ATC biosynthesis, facilitating the genetic enhancement of ATC contents and other related traits in radish.

Genetic diversity, SNP-trait associations and genomic selection accuracy in a west African collection of Kersting's groundnut [Macrotyloma geocarpum(Harms) Maréchal & Baudet]

Understanding the mechanisms governing complex traits variation is a requirement for efficient crop improvement. In this study, the molecular characterization, marker-trait associations and the possibility for genomic selection in a collection of 281 Kersting's groundnut accessions were carried out. The diversity panel was phenotyped using an Alpha lattice design with two replicates in two contrasting environments. Accessions were genotyped using genotyping by sequencing technology. Genome-wide association analyses were performed between single nucleotide polymorphism markers and yield-related traits across tested environments. SNP markers were used to calculate the observed (Ho) and expected heterozygosity (He), and the total gene diversity (Ht). Genetic differentiation among accessions across ecological regions of origin was analysed. Our results revealed 493 quality SNPs of which 113 had a minor allele frequency>0.05, a total gene diversity of 0.43 and average Ho and He values of 0.04 and 0.22, respectively. Four clusters, highly differentiated by seed coat colour (Fst = 0.79), were identified. The population structure analysis showed two subpopulations with high differentiation across ecological regions (Fst = 0.37). The GWAS revealed 10 significant marker-trait associations, of which six SNPs were consistent across environments. The genomic selection through cross-validation showed moderate to high prediction accuracies for leaflet length, seed dimension traits, 100 seed weight, days to 50% flowering and days to maturity. This demonstrates the existence of genetic variability within Kersting's groundnut and shows the potential for the improvement of the species. The findings also provide a first insight into the phenotype-to-genotype relationships in Kersting's groundnut, using SNP markers.

Genome-wide association studies for inflorescence type and remontancy in Hydrangea macrophylla

Inflorescence type and remontancy are two valuable traits in bigleaf hydrangea (Hydrangea macrophylla L.) and both are recessively inherited. Molecular marker-assisted selection (MAS) can greatly reduce the time necessary to breed cultivars with desired traits. In this study, a genome-wide association study (GWAS) using 5803 single-nucleotide polymorphisms (SNPs) was performed using a panel of 82 bigleaf hydrangea cultivars. One SNP locus (Hy_CAPS_Inflo) associated with inflorescence type was identified with general linear model (GLM) and mixed linear model (MLM) methods that explained 65.5% and 36.1% of the phenotypic variations, respectively. Twenty-three SNPs associated with remontancy were detected in GLM whereas no SNP was detected in MLM. The SNP locus (Hy_CAPS_Inflo) was converted to a cleaved amplified polymorphic sequence (CAPS) marker that showed absolute identification accuracy (100%) of inflorescence type in a validation panel consisting of eighteen H. macrophylla cultivars. The SNP was investigated in 341 F1 progenies using genotyping by sequencing (GBS) and co-segregated with inflorescence type (χ 2 = 0.12; P = 0.73). The SNP was subsequently used for breeding selection using kompetitive allele specific PCR (KASP) technology. Future directions for the use of genomics and MAS in hydrangea breeding improvement are discussed. The results presented in this study provide insights for further research on understanding genetic mechanisms behind inflorescence type and remontancy in H. macrophylla. The CAPS and KASP markers developed here will be immediately useful for applying MAS to accelerate breeding improvement in hydrangea.

Current achievements and future prospects in the genetic breeding of chrysanthemum: a review

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a leading flower with applied value worldwide. Developing new chrysanthemum cultivars with novel characteristics such as new flower colors and shapes, plant architectures, flowering times, postharvest quality, and biotic and abiotic stress tolerance in a time- and cost-efficient manner is the ultimate goal for breeders. Various breeding strategies have been employed to improve the aforementioned traits, ranging from conventional techniques, including crossbreeding and mutation breeding, to a series of molecular breeding methods, including transgenic technology, genome editing, and marker-assisted selection (MAS). In addition, the recent extensive advances in high-throughput technologies, especially genomics, transcriptomics, proteomics, metabolomics, and microbiomics, which are collectively referred to as omics platforms, have led to the collection of substantial amounts of data. Integration of these omics data with phenotypic information will enable the identification of genes/pathways responsible for important traits. Several attempts have been made to use emerging molecular and omics methods with the aim of accelerating the breeding of chrysanthemum. However, applying the findings of such studies to practical chrysanthemum breeding remains a considerable challenge, primarily due to the high heterozygosity and polyploidy of the species. This review summarizes the recent achievements in conventional and modern molecular breeding methods and emerging omics technologies and discusses their future applications for improving the agronomic and horticultural characteristics of chrysanthemum.