Genetic architecture of phenotypic means and plasticities of kernel size and weight in maize

Genome-Wide Association Study for Maize Hybrid Performance in a Typical Breeder Population

Maize is one of the major crops that has demonstrated success in the utilization of heterosis. Developing high-yield hybrids is a crucial part of plant breeding to secure global food demand. In this study, we conducted a genome-wide association study (GWAS) for 10 agronomic traits using a typical breeder population comprised 442 single-cross hybrids by evaluating additive, dominance, and epistatic effects. A total of 49 significant single nucleotide polymorphisms (SNPs) and 69 significant pairs of epistasis were identified, explaining 26.2% to 64.3% of the phenotypic variation across the 10 traits. The enrichment of favorable genotypes is significantly correlated to the corresponding phenotype. In the confident region of the associated site, 532 protein-coding genes were discovered. Among these genes, the Zm00001d044211 candidate gene was found to negatively regulate starch synthesis and potentially impact yield. This typical breeding population provided a valuable resource for dissecting the genetic architecture of yield-related traits. We proposed a novel mating strategy to increase the GWAS efficiency without utilizing more resources. Finally, we analyzed the enrichment of favorable alleles in the Shaan A and Shaan B groups, as well as in each inbred line. Our breeding practice led to consistent results. Not only does this study demonstrate the feasibility of GWAS in F1 hybrid populations, it also provides a valuable basis for further molecular biology and breeding research.

Genetic interrogation of phenotypic plasticity informs genome-enabled breeding in cotton

Phenotypic plasticity, or the ability to adapt to and thrive in changing climates and variable environments, is essential for developmental programs in plants. Despite its importance, the genetic underpinnings of phenotypic plasticity for key agronomic traits remain poorly understood in many crops. In this study, we aim to fill this gap by using genome-wide association studies to identify genetic variations associated with phenotypic plasticity in upland cotton (Gossypium hirsutum L.). We identified 73 additive quantitative trait loci (QTLs), 32 dominant QTLs, and 6799 epistatic QTLs associated with 20 traits. We also identified 117 additive QTLs, 28 dominant QTLs, and 4691 epistatic QTLs associated with phenotypic plasticity in 19 traits. Our findings reveal new genetic factors, including additive, dominant, and epistatic QTLs, that are linked to phenotypic plasticity and agronomic traits. Meanwhile, we find that the genetic factors controlling the mean phenotype and phenotypic plasticity are largely independent in upland cotton, indicating the potential for simultaneous improvement. Additionally, we envision a genomic design strategy by utilizing the identified QTLs to facilitate cotton breeding. Taken together, our study provides new insights into the genetic basis of phenotypic plasticity in cotton, which should be valuable for future breeding.

A loss-of-function mutant allele of a glycosyl hydrolase gene has been co-opted for seed weight control during soybean domestication

The resultant DNA from loss-of-function mutation can be recruited in biological evolution and development. Here, we present such a rare and potential case of "to gain by loss" as a neomorphic mutation during soybean domestication for increasing seed weight. Using a population derived from a chromosome segment substitution line of Glycine max (SN14) and Glycine soja (ZYD06), a quantitative trait locus (QTL) of 100-seed weight (qHSW) was mapped on chromosome 11, corresponding to a truncated β-1, 3-glucosidase (βGlu) gene. The novel gene hsw results from a 14-bp deletion, causing a frameshift mutation and a premature stop codon in the βGlu. In contrast to HSW, the hsw completely lost βGlu activity and function but acquired a novel function to promote cell expansion, thus increasing seed weight. Overexpressing hsw instead of HSW produced large soybean seeds, and surprisingly, truncating hsw via gene editing further increased the seed size. We further found that the core 21-aa peptide of hsw and its variants acted as a promoter of seed size. Transcriptomic variation in these transgenic soybean lines substantiated the integration hsw into cell and seed size control. Moreover, the hsw allele underwent selection and expansion during soybean domestication and improvement. Our work cloned a likely domesticated QTL controlling soybean seed weight, revealed a novel genetic variation and mechanism in soybean domestication, and provided new insight into crop domestication and breeding, and plant evolution.

Genetic analysis of phenotypic plasticity identifies BBX6 as the candidate gene for maize adaptation to temperate regions

Introduction

Climate changes pose a significant threat to crop adaptation and production. Dissecting the genetic basis of phenotypic plasticity and uncovering the responsiveness of regulatory genes to environmental factors can significantly contribute to the improvement of climate- resilience in crops.

Methods

We established a BC1F3:4 population using the elite inbred lines Zheng58 and PH4CV and evaluated plant height (PH) across four environments characterized by substantial variations in environmental factors. Then, we quantified the correlation between the environmental mean of PH (the mean performance in each environment) and the environmental parameters within a specific growth window. Furthermore, we performed GWAS analysis of phenotypic plasticity, and identified QTLs and candidate gene that respond to key environment index. After that, we constructed the coexpression network involving the candidate gene, and performed selective sweep analysis of the candidate gene.

Results

We found that the environmental parameters demonstrated substantial variation across the environments, and genotype by environment interaction contributed to the variations of PH. Then, we identified PTT(35-48) (PTT is the abbreviation for photothermal units), the mean PTT from 35 to 48 days after planting, as the pivotal environmental index that closely correlated with environmental mean of PH. Leveraging the slopes of the response of PH to both the environmental mean and PTT(35-48), we successfully pinpointed QTLs for phenotypic plasticity on chromosomes 1 and 2. Notably, the PH4CV genotypes at these two QTLs exhibited positive contributions to phenotypic plasticity. Furthermore, our analysis demonstrated a direct correlation between the additive effects of each QTL and PTT(35-48). By analyzing transcriptome data of the parental lines in two environments, we found that the 1009 genes responding to PTT(35-48) were enriched in the biological processes related to environmental sensitivity. BBX6 was the prime candidate gene among the 13 genes in the two QTL regions. The coexpression network of BBX6 contained other genes related to flowering time and photoperiod sensitivity. Our investigation, including selective sweep analysis and genetic differentiation analysis, suggested that BBX6 underwent selection during maize domestication.

Discussion

Th is research substantially advances our understanding of critical environmental factors influencing maize adaptation while simultaneously provides an invaluable gene resource for the development of climate-resilient maize hybrid varieties.

Complex genetic architecture underlying the plasticity of maize agronomic traits

Phenotypic plasticity is the ability of a given genotype to produce multiple phenotypes in response to changing environmental conditions. Understanding the genetic basis of phenotypic plasticity and establishing a predictive model is highly relevant to future agriculture under a changing climate. Here we report findings on the genetic basis of phenotypic plasticity for 23 complex traits using a diverse maize population planted at five sites with distinct environmental conditions. We found that latitude-related environmental factors were the main drivers of across-site variation in flowering time traits but not in plant architecture or yield traits. For the 23 traits, we detected 109 quantitative trait loci (QTLs), 29 for mean values, 66 for plasticity, and 14 for both parameters, and 80% of the QTLs interacted with latitude. The effects of several QTLs changed in magnitude or sign, driving variation in phenotypic plasticity. We experimentally validated one plastic gene, ZmTPS14.1, whose effect was likely mediated by the compensation effect of ZmSPL6 from a downstream pathway. By integrating genetic diversity, environmental variation, and their interaction into a joint model, we could provide site-specific predictions with increased accuracy by as much as 9.9%, 2.2%, and 2.6% for days to tassel, plant height, and ear weight, respectively. This study revealed a complex genetic architecture involving multiple alleles, pleiotropy, and genotype-by-environment interaction that underlies variation in the mean and plasticity of maize complex traits. It provides novel insights into the dynamic genetic architecture of agronomic traits in response to changing environments, paving a practical way toward precision agriculture.

Identifying yield-related genes in maize based on ear trait plasticity

BACKGROUND

Phenotypic plasticity is defined as the phenotypic variation of a trait when an organism is exposed to different environments, and it is closely related to genotype. Exploring the genetic basis behind the phenotypic plasticity of ear traits in maize is critical to achieve climate-stable yields, particularly given the unpredictable effects of climate change. Performing genetic field studies in maize requires development of a fast, reliable, and automated system for phenotyping large numbers of samples.

RESULTS

Here, we develop MAIZTRO as an automated maize ear phenotyping platform for high-throughput measurements in the field. Using this platform, we analyze 15 common ear phenotypes and their phenotypic plasticity variation in 3819 transgenic maize inbred lines targeting 717 genes, along with the wild type lines of the same genetic background, in multiple field environments in two consecutive years. Kernel number is chosen as the primary target phenotype because it is a key trait for improving the grain yield and ensuring yield stability. We analyze the phenotypic plasticity of the transgenic lines in different environments and identify 34 candidate genes that may regulate the phenotypic plasticity of kernel number.

CONCLUSIONS

Our results suggest that as an integrated and efficient phenotyping platform for measuring maize ear traits, MAIZTRO can help to explore new traits that are important for improving and stabilizing the yield. This study indicates that genes and alleles related with ear trait plasticity can be identified using transgenic maize inbred populations.

Partitioning the forms of genotype-by-environment interaction in the reaction norm analysis of stability

KEY MESSAGE

The reaction norm analysis of stability can be enhanced by partitioning the contribution of different types of G × E to the variation in slope. The slope of regression in a reaction norm model, where the performance of a genotype is regressed over an environmental covariable, is often used as a measure of stability of genotype performance. This method could be developed further by partitioning variation in the slope of regression into the two sources of genotype-by-environment interaction (G × E) which cause it: scale-type G × E (heterogeneity of variance) and rank-type G × E (heterogeneity of correlation). Because the two types of G × E have very different properties, separating their effect would enable a clearer understanding of stability. The aim of this paper was to demonstrate two methods which seek to achieve this in reaction norm models. Reaction norm models were fit to yield data from a multi-environment trial in Barley (Hordeum vulgare), with the adjusted mean yield from each environment used as the environmental covariable. Stability estimated from factor-analytic models, which can disentangle the two types of G × E and estimate stability based on rank-type G × E, was used for comparison. Adjusting the reaction norm slope to account for scale-type G × E using a genetic regression more than tripled the correlation with factor-analytic estimates of stability (0.24-0.26 to 0.80-0.85), indicating that it removed variation in the reaction norm slope that originated from scale-type G × E. A standardisation procedure had a more modest increase (055-0.59) but could be useful when curvilinear reaction norms are required. Analyses which use reaction norms to explore the stability of genotypes could gain additional insight into the mechanisms of stability by applying the methods outlined in this study.

Genome-Wide Association Analysis for Candidate Genes Contributing to Kernel-Related Traits in Maize

Maize grain size is the main factor determining grain yield. Dissecting the genetic basis of maize grain size may help reveal the regulatory mechanism of maize seed development and yield formation. In this study, two associated populations were used for genome-wide association analysis of kernel length, kernel width, kernel thickness, and hundred-kernel weight from multiple locations in AM122 and AM180, respectively. Then, genome-wide association mapping was performed based on the maize 6H90K SNP chip. A total of 139 loci were identified as associated with the four traits with p < 1 × 10^-4 using two models (FarmCPU and MLM). The transcriptome data showed that 15 of them were expressed differentially in two maize-inbred lines KB182 (small kernel) and KB020 (big kernel) during kernel development. These candidate genes were enriched in regulating peroxidase activity, oxidoreductase, and leaf senescence. The molecular function was major in binding and catalytic activity. This study provided important reference information for exploring maize kernel development mechanisms and applying molecular markers in high-yield breeding.

Phenotypic Plasticity Contributes to Maize Adaptation and Heterosis

Plant phenotypic plasticity describes altered phenotypic performance of an individual when grown in different environments. Exploring genetic architecture underlying plant plasticity variation may help mitigate the detrimental effects of a rapidly changing climate on agriculture, but little research has been done in this area to date. In the present study, we established a population of 976 maize F₁ hybrids by crossing 488 diverse inbred lines with two elite testers. Genome-wide association study identified hundreds of quantitative trait loci associated with phenotypic plasticity variation across diverse F₁ hybrids, the majority of which contributed very little variance, in accordance with the polygenic nature of these traits. We identified several quantitative trait locus regions that may have been selected during the tropical-temperate adaptation process. We also observed heterosis in terms of phenotypic plasticity, in addition to the traditional genetic value differences measured between hybrid and inbred lines, and the pattern of which was affected by genetic background. Our results demonstrate a landscape of phenotypic plasticity in maize, which will aid in the understanding of its genetic architecture, its contribution to adaptation and heterosis, and how it may be exploited for future maize breeding in a rapidly changing environment.

Genetic determinants of seed protein plasticity in response to the environment in Medicago truncatula

As the frequency of extreme environmental events is expected to increase with climate change, identifying candidate genes for stabilizing the protein composition of legume seeds or optimizing this in a given environment is increasingly important. To elucidate the genetic determinants of seed protein plasticity, major seed proteins from 200 ecotypes of Medicago truncatula grown in four contrasting environments were quantified after one-dimensional electrophoresis. The plasticity index of these proteins was recorded for each genotype as the slope of Finlay and Wilkinson's regression and then used for genome-wide association studies (GWASs), enabling the identification of candidate genes for determining this plasticity. This list was enriched in genes related to transcription, DNA repair and signal transduction, with many of them being stress responsive. Other over-represented genes were related to sulfur and aspartate family pathways leading to the synthesis of the nutritionally essential amino acids methionine and lysine. By placing these genes in metabolic pathways, and using a M. truncatula mutant impaired in regenerating methionine from S-methylmethionine, we discovered that methionine recycling pathways are major contributors to globulin composition establishment and plasticity. These data provide a unique resource of genes that can be targeted to mitigate negative impacts of environmental stresses on seed protein composition.

Association Mapping and Transcriptome Analysis Reveal the Genetic Architecture of Maize Kernel Size

Kernel length, kernel width, and kernel thickness are important traits affecting grain yield and product quality. Here, the genetic architecture of the three kernel size traits was dissected in an association panel of 309 maize inbred lines using four statistical methods. Forty-two significant single nucleotide polymorphisms (SNPs; p < 1.72E-05) and 70 genes for the three traits were identified under five environments. One and eight SNPs were co-detected in two environments and by at least two methods, respectively, and they explained 5.87-9.59% of the phenotypic variation. Comparing the transcriptomes of two inbred lines with contrasting seed size, three and eight genes identified in the association panel showed significantly differential expression between the two genotypes at 15 and 39 days after pollination, respectively. Ten and 17 genes identified by a genome-wide association study were significantly differentially expressed between the two development stages in the two genotypes. Combining environment-/method-stable SNPs and differential expression analysis, ribosomal protein L7, jasmonate-regulated gene 21, serine/threonine-protein kinase RUNKEL, AP2-EREBP-transcription factor 16, and Zm00001d035222 (cell wall protein IFF6-like) were important candidate genes for maize kernel size and development.

Wheat individual grain-size variance originates from crop development and from specific genetic determinism

Wheat grain yield is usually decomposed in the yield components: number of spikes / m2, number of grains / spike, number of grains / m2 and thousand kernel weight (TKW). These are correlated one with another due to yield component compensation. Under optimal conditions, the number of grains per m2 has been identified as the main determinant of yield. However, with increasing occurrences of post-flowering abiotic stress associated with climate change, TKW may become severely limiting and hence a target for breeding. TKW is usually studied at the plot scale as it represents the average mass of a grain. However, this view disregards the large intra-genotypic variance of individual grain mass and its effect on TKW. The aim of this study is to investigate the determinism of the variance of individual grain size. We measured yield components and individual grain size variances of two large genetic wheat panels grown in two environments. We also carried out a genome-wide association study using a dense SNPs array. We show that the variance of individual grain size partly originates from the pre-flowering components of grain yield; in particular it is driven by canopy structure via its negative correlation with the number of spikes per m2. But the variance of final grain size also has a specific genetic basis. The genome-wide analysis revealed the existence of QTL with strong effects on the variance of individual grain size, independently from the other yield components. Finally, our results reveal some interesting drivers for manipulating individual grain size variance either through canopy structure or through specific chromosomal regions.