Identification of QTN-by-environment interactions for yield related traits in maize under multiple abiotic stresses

Genome-wide association mapping and genomic prediction analyses reveal the genetic architecture of grain yield and agronomic traits under drought and optimum conditions in maize

BACKGROUND

Drought is a major abiotic stress in sub-Saharan Africa, impacting maize growth and development leading to severe yield loss. Drought tolerance is a complex trait regulated by multiple genes, making direct grain yield selection ineffective. To dissect the genetic architecture of grain yield and flowering traits under drought stress, a genome-wide association study (GWAS) was conducted on a panel of 236 maize lines testcrossed and evaluated under managed drought and optimal growing conditions in multiple environments using seven multi-locus GWAS models (mrMLM, FASTmrMLM, FASTmrEMMA, pLARmEB, pKWmEB, ISIS EM-BLASSO, and FARMCPU) from mrMLM and GAPIT R packages. Genomic prediction with RR-BLUP model was applied on BLUEs across locations under optimum and drought conditions.

RESULTS

A total of 172 stable and reliable quantitative trait nucleotides (QTNs) were identified, of which 77 are associated with GY, AD, SD, ASI, PH, EH, EPO and EPP under drought and 95 are linked to GY, AD, SD, ASI, PH, EH, EPO and EPP under optimal conditions. Among these QTNs, 17 QTNs explained over 10% of the phenotypic variation (R2 ≥ 10%). Furthermore, 43 candidate genes were discovered and annotated. Two major candidate genes, Zm00001eb041070 closely associated with grain yield near peak QTN, qGY_DS1.1 (S1_216149215) and Zm00001eb364110 closely related to anthesis-silking interval near peak QTN, qASI_DS8.2 (S8_167256316) were identified, encoding AP2-EREBP transcription factor 60 and TCP-transcription factor 20, respectively under drought stress. Haplo-pheno analysis identified superior haplotypes for qGY_DS1.1 (S1_216149215) associated with the higher grain yield under drought stress. Genomic prediction revealed moderate to high prediction accuracies under optimum and drought conditions.

CONCLUSION

The lines carrying superior haplotypes can be used as potential donors in improving grain yield under drought stress. Integration of genomic selection with GWAS results leads not only to an increase in the prediction accuracy but also to validate the function of the identified candidate genes as well increase in the accumulation of favorable alleles with minor and major effects in elite breeding lines. This study provides valuable insight into the genetic architecture of grain yield and secondary traits under drought stress.

Integrating GWAS with a gene co-expression network better prioritizes candidate genes associated with root metaxylem phenes in maize

Root metaxylems are phenotypically diverse structures whose function is particularly important under drought stress. Significant research has dissected the genetic machinery underlying metaxylem phenotypes in dicots, but that of monocots are relatively underexplored. In maize (Zea mays), a robust pipeline integrated a genome-wide association study (GWAS) of root metaxylem phenes under well-watered and water-stress conditions with a gene co-expression network to prioritize the strongest gene candidates. We identified 244 candidate genes by GWAS, of which 103 reside in gene co-expression modules most relevant to xylem development. Several candidate genes may be involved in biosynthetic processes related to the cell wall, hormone signaling, oxidative stress responses, and drought responses. Of those, six gene candidates were detected in multiple root metaxylem phenes in both well-watered and water-stress conditions. We posit that candidate genes that are more essential to network function based on gene co-expression (i.e., hubs or bottlenecks) should be prioritized and classify 33 essential genes for further investigation. Our study demonstrates a new strategy for identifying promising gene candidates and presents several gene candidates that may enhance our understanding of vascular development and responses to drought in cereals.

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.

Pilot-scale genome-wide association mapping in diverse sorghum germplasms identified novel genetic loci linked to major agronomic, root and stomatal traits

This genome-wide association studies (GWAS) used a subset of 96 diverse sorghum accessions, constructed from a large collection of 219 accessions for mining novel genetic loci linked to major agronomic, root morphological and physiological traits. The subset yielded 43,452 high quality single nucleotide polymorphic (SNP) markers exhibiting high allelic diversity. Population stratification showed distinct separation between caudatum and durra races. Linkage disequilibrium (LD) decay was rapidly declining with increasing physical distance across all chromosomes. The initial 50% LD decay was ~ 5 Kb and background level was within ~ 80 Kb. This study detected 42 significant quantitative trait nucleotide (QTNs) for different traits evaluated using FarmCPU, SUPER and 3VmrMLM which were in proximity with candidate genes related and were co-localized in already reported quantitative trait loci (QTL) and phenotypic variance (R2) of these QTNs ranged from 3 to 20%. Haplotype validation of the candidate genes from this study resulted nine genes showing significant phenotypic difference between different haplotypes. Three novel candidate genes associated with agronomic traits were validated including Sobic.001G499000, a potassium channel tetramerization domain protein for plant height, Sobic.010G186600, a nucleoporin-related gene for dry biomass, and Sobic.002G022600 encoding AP2-like ethylene-responsive transcription factor for plant yield. Several other candidate genes were validated and associated with different root and physiological traits including Sobic.005G104100, peroxidase 13-related gene with root length, Sobic.010G043300, homologous to Traes_5BL_8D494D60C, encoding inhibitor of apoptosis with iWUE, and Sobic.010G125500, encoding zinc finger, C3HC4 type domain with Abaxial stomatal density. In this study, 3VmrMLM was more powerful than FarmCPU and SUPER for detecting QTNs and having more breeding value indicating its reliable output for validation. This study justified that the constructed subset of diverse sorghums can be used as a panel for mapping other key traits to accelerate molecular breeding in sorghum.