Gene-set association and epistatic analyses reveal complex gene interaction networks affecting flowering time in a worldwide barley collection

Investigating the genetic control of plant development in spring barley under speed breeding conditions

KEY MESSAGE

This study found that the genes, PPD-H1 and ELF3, control the acceleration of plant development under speed breeding, with important implications for optimizing the delivery of climate-resilient crops. Speed breeding is a tool to accelerate breeding and research programmes. Despite its success and growing popularity with breeders, the genetic basis of plant development under speed breeding remains unknown. This study explored the developmental advancements of barley genotypes under different photoperiod regimes. A subset of the HEB-25 Nested Association Mapping population was evaluated for days to heading and maturity under two contrasting photoperiod conditions: (1) Speed breeding (SB) consisting of 22 h of light and 2 h of darkness, and (2) normal breeding (NB) consisting of 16 h of light and 8 h of darkness. GWAS revealed that developmental responses under both conditions were largely controlled by two loci: PPDH-1 and ELF3. Allelic variants at these genes determine whether plants display early flowering and maturity under both conditions. At key QTL regions, domesticated alleles were associated with late flowering and maturity in NB and early flowering and maturity in SB, whereas wild alleles were associated with early flowering under both conditions. We hypothesize that this is related to the dark-dependent repression of PPD-H1 by ELF3 which might be more prominent in NB conditions. Furthermore, by comparing development under two photoperiod regimes, we derived an estimate of plasticity for the two traits. Interestingly, plasticity in development was largely attributed to allelic variation at ELF3. Our results have important implications for our understanding and optimization of speed breeding protocols particularly for introgression breeding and the design of breeding programmes to support the delivery of climate-resilient crops.

Unlocking the genetic diversity and population structure of the newly introduced two-row spring European HerItage Barley collecTion (ExHIBiT)

In the last century, breeding programs have traditionally favoured yield-related traits, grown under high-input conditions, resulting in a loss of genetic diversity and an increased susceptibility to stresses in crops. Thus, exploiting understudied genetic resources, that potentially harbour tolerance genes, is vital for sustainable agriculture. Northern European barley germplasm has been relatively understudied despite its key role within the malting industry. The European Heritage Barley collection (ExHIBiT) was assembled to explore the genetic diversity in European barley focusing on Northern European accessions and further address environmental pressures. ExHIBiT consists of 363 spring-barley accessions, focusing on two-row type. The collection consists of landraces (~14%), old cultivars (~18%), elite cultivars (~67%) and accessions with unknown breeding history (~1%), with 70% of the collection from Northern Europe. The population structure of the ExHIBiT collection was subdivided into three main clusters primarily based on the accession's year of release using 26,585 informative SNPs based on 50k iSelect single nucleotide polymorphism (SNP) array data. Power analysis established a representative core collection of 230 genotypically and phenotypically diverse accessions. The effectiveness of this core collection for conducting statistical and association analysis was explored by undertaking genome-wide association studies (GWAS) using 24,876 SNPs for nine phenotypic traits, four of which were associated with SNPs. Genomic regions overlapping with previously characterised flowering genes (HvZTLb) were identified, demonstrating the utility of the ExHIBiT core collection for locating genetic regions that determine important traits. Overall, the ExHIBiT core collection represents the high level of untapped diversity within Northern European barley, providing a powerful resource for researchers and breeders to address future climate scenarios.

Compressed variance component mixed model reveals epistasis associated with flowering in Arabidopsis

Introduction

Epistasis is currently a topic of great interest in molecular and quantitative genetics. Arabidopsis thaliana, as a model organism, plays a crucial role in studying the fundamental biology of diverse plant species. However, there have been limited reports about identification of epistasis related to flowering in genome-wide association studies (GWAS). Therefore, it is of utmost importance to conduct epistasis in Arabidopsis.

Method

In this study, we employed Levene's test and compressed variance component mixed model in GWAS to detect quantitative trait nucleotides (QTNs) and QTN-by-QTN interactions (QQIs) for 11 flowering-related traits of 199 Arabidopsis accessions with 216,130 markers.

Results

Our analysis detected 89 QTNs and 130 pairs of QQIs. Around these loci, 34 known genes previously reported in Arabidopsis were confirmed to be associated with flowering-related traits, such as SPA4, which is involved in regulating photoperiodic flowering, and interacts with PAP1 and PAP2, affecting growth of Arabidopsis under light conditions. Then, we observed significant and differential expression of 35 genes in response to variations in temperature, photoperiod, and vernalization treatments out of unreported genes. Functional enrichment analysis revealed that 26 of these genes were associated with various biological processes. Finally, the haplotype and phenotypic difference analysis revealed 20 candidate genes exhibiting significant phenotypic variations across gene haplotypes, of which the candidate genes AT1G12990 and AT1G09950 around QQIs might have interaction effect to flowering time regulation in Arabidopsis.

Discussion

These findings may offer valuable insights for the identification and exploration of genes and gene-by-gene interactions associated with flowering-related traits in Arabidopsis, that may even provide valuable reference and guidance for the research of epistasis in other species.

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.

Sense in sensitivity: difference in the meaning of photoperiod-insensitivity between wheat and barley

The description of long photoperiod sensitivity in wheat and barley is a cause of confusion for researchers working in these crops, usually accustomed to free exchange of physiological and genetic knowledge of such similar crops. Indeed, wheat and barley scientists customarily quote studies of either crop species when researching one of them. Among their numerous similarities the main gene controlling that response is the same in both crops (PPD1; PPD-H1 in barley and PPD-D1 in hexaploid wheat). However, the photoperiod responses are different: (i) the main dominant allele inducing shorter time to anthesis is the insensitive allele in wheat (Ppd-D1a) but the sensitive allele in barley (Ppd-H1) (i.e. sensitivity to photoperiod produces opposite effects on time to heading in wheat and barley), (ii) the main "insensitive" allele in wheat, Ppd-D1a, does confer insensitivity, whilst that of barley reduces the sensitivity but still responds to photoperiod. The different behaviour of PPD1 genes in wheat and barley is put in a common framework based on the similarities and differences of the molecular bases of their mutations, which include polymorphism at gene expression levels, copy number variation, and sequence of coding regions. This common perspective sheds light on a source on confusion for cereal researchers, and prompts us to recommend accounting for the photoperiod sensitivity status of the plant materials when doing research on genetic control of phenology. Finally, we provide advice to facilitate the management of natural PPD1 diversity in breeding programs and suggest targets for further modification through gene editing, based on mutual knowledge on the two crops.

New semi-dwarfing alleles with increased coleoptile length by gene editing of gibberellin 3-oxidase 1 using CRISPR-Cas9 in barley (Hordeum vulgare L.)

The green revolution was based on genetic modification of the gibberellin (GA) hormone system with "dwarfing" gene mutations that reduces GA signals, conferring shorter stature, thus enabling plant adaptation to modern farming conditions. Strong GA-related mutants with shorter stature often have reduced coleoptile length, discounting yield gain due to their unsatisfactory seedling emergence under drought conditions. Here we present gibberellin (GA) 3-oxidase1 (GA3ox1) as an alternative semi-dwarfing gene in barley that combines an optimal reduction in plant height without restricting coleoptile and seedling growth. Using large-scale field trials with an extensive collection of barley accessions, we showed that a natural GA3ox1 haplotype moderately reduced plant height by 5-10 cm. We used CRISPR/Cas9 technology, generated several novel GA3ox1 mutants and validated the function of GA3ox1. We showed that altered GA3ox1 activities changed the level of active GA isoforms and consequently increased coleoptile length by an average of 8.2 mm, which could provide essential adaptation to maintain yield under climate change. We revealed that CRISPR/Cas9-induced GA3ox1 mutations increased seed dormancy to an ideal level that could benefit the malting industry. We conclude that selecting HvGA3ox1 alleles offers a new opportunity for developing barley varieties with optimal stature, longer coleoptile and additional agronomic traits.

Pleiotropy Structures Plant Height and Seed Weight Scaling in Barley despite Long History of Domestication and Breeding Selection

Size scaling describes the relative growth rates of different body parts of an organism following a positive correlation. Domestication and crop breeding often target the scaling traits in the opposite directions. The genetic mechanism of the size scaling influencing the pattern of size scaling remains unexplored. Here, we revisited a diverse barley (Hordeum vulgare L.) panel with genome-wide single-nucleotide polymorphisms (SNPs) profile and the measurement of their plant height and seed weight to explore the possible genetic mechanisms that may lead to a correlation of the two traits and the influence of domestication and breeding selection on the size scaling. Plant height and seed weight are heritable and remain positively correlated in domesticated barley regardless of growth type and habit. Genomic structural equation modeling systematically evaluated the pleiotropic effect of individual SNP on the plant height and seed weight within a trait correlation network. We discovered seventeen novel SNPs (quantitative trait locus) conferring pleiotropic effect on plant height and seed weight, involving genes with function in diverse traits related to plant growth and development. Linkage disequilibrium decay analysis revealed that a considerable proportion of genetic markers associated with either plant height or seed weight are closely linked in the chromosome. We conclude that pleiotropy and genetic linkage likely form the genetic bases of plant height and seed weight scaling in barley. Our findings contribute to understanding the heritability and genetic basis of size scaling and open a new venue for seeking the underlying mechanism of allometric scaling in plants.

Genome-wide association study of agronomical and root-related traits in spring barley collection grown under field conditions

The root system is a key component for plant survival and productivity. In particular, under stress conditions, developing plants with a better root architecture can ensure productivity. The objectives of this study were to investigate the phenotypic variation of selected root- and yield-related traits in a diverse panel of spring barley genotypes. By performing a genome-wide association study (GWAS), we identified several associations underlying the variations occurring in root- and yield-related traits in response to natural variations in soil moisture. Here, we report the results of the GWAS based on both individual single-nucleotide polymorphism markers and linkage disequilibrium (LD) blocks of markers for 11 phenotypic traits related to plant morphology, grain quality, and root system in a group of spring barley accessions grown under field conditions. We also evaluated the root structure of these accessions by using a nondestructive method based on electrical capacitance. The results showed the importance of two LD-based blocks on chromosomes 2H and 7H in the expression of root architecture and yield-related traits. Our results revealed the importance of the region on the short arm of chromosome 2H in the expression of root- and yield-related traits. This study emphasized the pleiotropic effect of this region with respect to heading time and other important agronomic traits, including root architecture. Furthermore, this investigation provides new insights into the roles played by root traits in the yield performance of barley plants grown under natural conditions with daily variations in soil moisture content.

Genome architecture and diverged selection shaping pattern of genomic differentiation in wild barley

Divergent selection of populations in contrasting environments leads to functional genomic divergence. However, the genomic architecture underlying heterogeneous genomic differentiation remains poorly understood. Here, we de novo assembled two high-quality wild barley (Hordeum spontaneum K. Koch) genomes and examined genomic differentiation and gene expression patterns under abiotic stress in two populations. These two populations had a shared ancestry and originated in close geographic proximity but experienced different selective pressures due to their contrasting micro-environments. We identified structural variants that may have played significant roles in affecting genes potentially associated with well-differentiated phenotypes such as flowering time and drought response between two wild barley genomes. Among them, a 29-bp insertion into the promoter region formed a cis-regulatory element in the HvWRKY45 gene, which may contribute to enhanced tolerance to drought. A single SNP mutation in the promoter region may influence HvCO5 expression and be putatively linked to local flowering time adaptation. We also revealed significant genomic differentiation between the two populations with ongoing gene flow. Our results indicate that SNPs and small SVs link to genetic differentiation at the gene level through local adaptation and are maintained through divergent selection. In contrast, large chromosome inversions may have shaped the heterogeneous pattern of genomic differentiation along the chromosomes by suppressing chromosome recombination and gene flow. Our research offers novel insights into the genomic basis underlying local adaptation and provides valuable resources for the genetic improvement of cultivated barley.

Population origin determines the adaptive potential for the advancement of flowering onset in Lupinus angustifolius L. (Fabaceae)

In the present framework of global warming, it is unclear whether evolutionary adaptation can happen quick enough to preserve the persistence of many species. Specifically, we lack knowledge about the adaptive potential of the different populations in relation to the various constraints that may hamper particular adaptations. There is evidence indicating that early flowering often provides an adaptive advantage to plants in temperate zones in response to global warming. Thus, the objective of this study was to assess the adaptive potential for advancing flowering onset in Lupinus angustifolius L. (Fabaceae). Seeds from four populations from two contrasting latitudes in Spain were collected and sown in a common garden environment. Selecting the 25% of the individuals that flowered earlier in the first generation, over three generations, three different early flowering selection lines were established, involving both self-crosses and outcrosses. All artificial selection lines advanced their flowering significantly with respect to the control line in the northernmost populations, but not in the southern ones. Selection lines obtained from outcrossing had a greater advancement in flowering than those from self-crossing. No differences were found in the number or weight of the seeds produced between control and artificial selection lines, probably because plants in the common garden were drip irrigated. These results suggest that northern populations may have a greater adaptive potential and that southern populations may be more vulnerable in the context of climate warming. However, earlier flowering was also associated with changes in other traits (height, biomass, shoot growth, specific leaflet area, and leaflet dry matter content), and the effects of these changes varied greatly depending on the latitude of the population and selection line. Assessments of the ability of populations to cope with climate change through this and other approaches are essential to manage species and populations in a more efficient way.

Multi-environment genome -wide association mapping of culm morphology traits in barley

In cereals with hollow internodes, lodging resistance is influenced by morphological characteristics such as internode diameter and culm wall thickness. Despite their relevance, knowledge of the genetic control of these traits and their relationship with lodging is lacking in temperate cereals such as barley. To fill this gap, we developed an image analysis-based protocol to accurately phenotype culm diameters and culm wall thickness across 261 barley accessions. Analysis of culm trait data collected from field trials in seven different environments revealed high heritability values (>50%) for most traits except thickness and stiffness, as well as genotype-by-environment interactions. The collection was structured mainly according to row-type, which had a confounding effect on culm traits as evidenced by phenotypic correlations. Within both row-type subsets, outer diameter and section modulus showed significant negative correlations with lodging (<-0.52 and <-0.45, respectively), but no correlation with plant height, indicating the possibility of improving lodging resistance independent of plant height. Using 50k iSelect SNP genotyping data, we conducted multi-environment genome-wide association studies using mixed model approach across the whole panel and row-type subsets: we identified a total of 192 quantitative trait loci (QTLs) for the studied traits, including subpopulation-specific QTLs and 21 main effect loci for culm diameter and/or section modulus showing effects on lodging without impacting plant height. Providing insights into the genetic architecture of culm morphology in barley and the possible role of candidate genes involved in hormone and cell wall-related pathways, this work supports the potential of loci underpinning culm features to improve lodging resistance and increase barley yield stability under changing environments.

Leaf thickness of barley: genetic dissection, candidate genes prediction and its relationship with yield-related traits

In this first genetic study on assessing leaf thickness directly in cereals, major and environmentally stable QTL were detected in barley and candidate genes underlying a major locus were identified. Leaf thickness (LT) is an important characteristic affecting leaf functions which have been intensively studied. However, as LT has a small dimension in many plant species and technically difficult to measure, previous studies on this characteristic are often based on indirect estimations. In the first study of detecting QTL controlling LT by directly measuring the characteristic in barley, large and stable loci were detected from both field and glasshouse trials conducted in different cropping seasons by assessing a population of 201 recombinant inbred lines. Four loci (locating on chromosome arms 2H, 3H, 5H and 6H, respectively) were consistently detected for flag leaf thickness (FLT) in each of these trials. The one on 6H had the largest effect, with a maximum LOD 9.8 explaining up to 20.9% of phenotypic variance. FLT does not only show strong interactions with flag leaf width and flag leaf area but has also strong correlations with fertile tiller number, spike row types, kernel number per spike and heading date. Though with reduced efficiency, these loci were also detectable from assessing second last leaf of fully grown plants or even from assessing the third leaves of seedlings. Taking advantage of the high-quality genome assemblies for both parents of the mapping population used in this study, three candidate genes underlying the 6H QTL were predicted based on orthologous analysis. These results do not only broaden our understanding on genetic basis of LT and its relationship with other traits in cereal crops but also form the bases for cloning and functional analysis of genes regulating LT in barley.

Comparative de novo transcriptome analysis identifies salinity stress responsive genes and metabolic pathways in sugarcane and its wild relative Erianthus arundinaceus [Retzius] Jeswiet

Erianthus arundinaceus [Retzius] Jeswiet, a wild relative of sugarcane has a high biomass production potential and a reservoir of many genes for superior agronomic traits and tolerance to biotic and abiotic stresses. A comparative physiological, anatomical and root transcriptome analysis were carried out to identify the salt-responsive genes and metabolic pathways associated with salt-tolerant E. arundinaceus genotype IND99-907 and salinity-sensitive sugarcane genotype Co 97010. IND99-907 recorded growth of young leaves, higher proline content, higher relative water content, intact root anatomical structures and lower Na+/K+, Ca2+/K+ and Mg2+/K+ ratio as compared to the sugarcane genotype Co 97010. We have generated four de novo transcriptome assemblies between stressed and control root samples of IND99-907 and Co 97010. A total of 649 and 501 differentially expressed genes (FDR<0.01) were identified from the stressed and control libraries of IND99-907 and Co 97010 respectively. Genes and pathways related to early stress-responsive signal transduction, hormone signalling, cytoskeleton organization, cellular membrane stabilization, plasma membrane-bound calcium and proton transport, sodium extrusion, secondary metabolite biosynthesis, cellular transporters related to plasma membrane-bound trafficking, nucleobase transporter, clathrin-mediated endocytosis were highly enriched in IND99-907. Whereas in Co 97010, genes related to late stress-responsive signal transduction, electron transport system, senescence, protein degradation and programmed cell death, transport-related genes associated with cellular respiration and mitochondrial respiratory chain, vesicular trafficking, nitrate transporter and fewer secondary metabolite biosynthetic genes were highly enriched. A total of 27 pathways, 24 biological processes, three molecular functions and one cellular component were significantly enriched (FDR≤ 0.05) in IND99-907 as compared to 20 pathways, two biological processes without any significant molecular function and cellular components in Co 97010, indicates the unique and distinct expression pattern of genes and metabolic pathways in both genotypes. The genomic resources developed from this study is useful for sugarcane crop improvement through development of genic SSR markers and genetic engineering approaches.

Genomic structural equation modelling provides a whole-system approach for the future crop breeding

Using genomic structural equation modelling, this research demonstrates an efficient way to identify genetically correlating traits and provides an effective proxy for multi-trait selection to consider the joint genetic architecture of multiple interacting traits in crop breeding. Breeding crop cultivars with optimal value across multiple traits has been a challenge, as traits may negatively correlate due to pleiotropy or genetic linkage. For example, grain yield and grain protein content correlate negatively with each other in cereal crops. Future crop breeding needs to be based on practical yet accurate evaluation and effective selection of beneficial trait to retain genes with the best agronomic score for multiple traits. Here, we test the framework of whole-system-based approach using structural equation modelling (SEM) to investigate how one trait affects others to guide the optimal selection of a combination of agronomically important traits. Using ten traits and genome-wide SNP profiles from a worldwide barley panel and SEM analysis, we revealed a network of interacting traits, in which tiller number contributes positively to both grain yield and protein content; we further identified common genetic factors affecting multiple traits in the network of interaction. Our method demonstrates an efficient way to identify genetically correlating traits and underlying pleiotropic genetic factors and provides an effective proxy for multi-trait selection within a whole-system framework that considers the joint genetic architecture of multiple interacting traits in crop breeding. Our findings suggest the promise of a whole-system approach to overcome challenges such as the negative correlation of grain yield and protein content to facilitating quantitative and objective breeding decisions in future crop breeding.

Major flowering time genes of barley: allelic diversity, effects, and comparison with wheat

This review summarizes the allelic series, effects, interactions between genes and with the environment, for the major flowering time genes that drive phenological adaptation of barley. The optimization of phenology is a major goal of plant breeding addressing the production of high-yielding varieties adapted to changing climatic conditions. Flowering time in cereals is regulated by genetic networks that respond predominately to day length and temperature. Allelic diversity at these genes is at the basis of barley wide adaptation. Detailed knowledge of their effects, and genetic and environmental interactions will facilitate plant breeders manipulating flowering time in cereal germplasm enhancement, by exploiting appropriate gene combinations. This review describes a catalogue of alleles found in QTL studies by barley geneticists, corresponding to the genetic diversity at major flowering time genes, the main drivers of barley phenological adaptation: VRN-H1 (HvBM5A), VRN-H2 (HvZCCTa-c), VRN-H3 (HvFT1), PPD-H1 (HvPRR37), PPD-H2 (HvFT3), and eam6/eps2 (HvCEN). For each gene, allelic series, size and direction of QTL effects, interactions between genes and with the environment are presented. Pleiotropic effects on agronomically important traits such as grain yield are also discussed. The review includes brief comments on additional genes with large effects on phenology that became relevant in modern barley breeding. The parallelisms between flowering time allelic variation between the two most cultivated Triticeae species (barley and wheat) are also outlined. This work is mostly based on previously published data, although we added some new data and hypothesis supported by a number of studies. This review shows the wide variety of allelic effects that provide enormous plasticity in barley flowering behavior, which opens new avenues to breeders for fine-tuning phenology of the barley crop.

A global barley panel revealing genomic signatures of breeding in modern Australian cultivars

The future of plant cultivar improvement lies in the evaluation of genetic resources from currently available germplasm. Today's gene pool of crop genetic diversity has been shaped during domestication and more recently by breeding. Recent efforts in plant breeding have been aimed at developing new and improved varieties from poorly adapted crops to suit local environments. However, the impact of these breeding efforts is poorly understood. Here, we assess the contributions of both historical and recent breeding efforts to local adaptation and crop improvement in a global barley panel by analysing the distribution of genetic variants with respect to geographic region or historical breeding category. By tracing the impact that breeding had on the genetic diversity of Hordeum vulgare (barley) released in Australia, where the history of barley production is relatively young, we identify 69 candidate regions within 922 genes that were under selection pressure. We also show that modern Australian barley varieties exhibit 12% higher genetic diversity than historical cultivars. Finally, field-trialling and phenotyping for agriculturally relevant traits across a diverse range of Australian environments suggests that genomic regions under strong breeding selection and their candidate genes are closely associated with key agronomic traits. In conclusion, our combined data set and germplasm collection provide a rich source of genetic diversity that can be applied to understanding and improving environmental adaptation and enhanced yields.

Haplotype-Based, Genome-Wide Association Study Reveals Stable Genomic Regions for Grain Yield in CIMMYT Spring Bread Wheat

We untangled key regions of the genetic architecture of grain yield (GY) in CIMMYT spring bread wheat by conducting a haplotype-based, genome-wide association study (GWAS), together with an investigation of epistatic interactions using seven large sets of elite yield trials (EYTs) consisting of a total of 6,461 advanced breeding lines. These lines were phenotyped under irrigated and stress environments in seven growing seasons (2011-2018) and genotyped with genotyping-by-sequencing markers. Genome-wide 519 haplotype blocks were constructed, using a linkage disequilibrium-based approach covering 14,036 Mb in the wheat genome. Haplotype-based GWAS identified 7, 4, 10, and 15 stable (significant in three or more EYTs) associations in irrigated (I), mild drought (MD), severe drought (SD), and heat stress (HS) testing environments, respectively. Considering all EYTs and the four testing environments together, 30 stable associations were deciphered with seven hotspots identified on chromosomes 1A, 1B, 2B, 4A, 5B, 6B, and 7B, where multiple haplotype blocks were associated with GY. Epistatic interactions contributed significantly to the genetic architecture of GY, explaining variation of 3.5-21.1%, 3.7-14.7%, 3.5-20.6%, and 4.4- 23.1% in I, MD, SD, and HS environments, respectively. Our results revealed the intricate genetic architecture of GY, controlled by both main and epistatic effects. The importance of these results for practical applications in the CIMMYT breeding program is discussed.

Cryptochromes in the field: how blue light influences crop development

Light is a pivotal environmental element capable of influencing multiple physiological processes across the entire plant life cycle. Over the course of their evolution, plants have developed several families of photoreceptors such as phytochromes, phototropins, ultraviolet (UV) resistance locus 8 and cryptochromes (crys), in order to sense light stimuli and respond to their changes. Numerous genetic studies have demonstrated that functional alterations to these photoreceptors cause a change in important agronomical traits. In particular, crys, which absorb UVA/blue light, can influence seed germination, flowering induction, plant architecture, fruit metabolic content and resistance to biotic and abiotic stresses. In the years to come, the rising temperatures and alterations to precipitation patterns generated by climate change will present a dramatic challenge for our agricultural system, with its few varieties characterized by a narrow genetic pool derived from artificial selection. Here, we review the main roles of crys in determining important agronomic traits in crops, we discuss the opportunities of using these photoreceptors as genetic targets for tuning plant physiological responses to environmental change, and the molecular strategies used so far to manipulate this family of photoreceptors.

The role of epistasis and its interaction with environment in fine-tuning heading time in barley

This article comments on: Afsharyan NP, Sannemann W, Léon J, Ballvora A. 2020. Effect of epistasis and environment on flowering time of barley reveals novel flowering-delaying QTL allele. Journal of Experimental Botany 71, 893–906.

Genome-wide association and epistatic interactions of flowering time in soybean cultivar

Genome-wide association studies (GWAS) have enabled the discovery of candidate markers that play significant roles in various complex traits in plants. Recently, with increased interest in the search for candidate markers, studies on epistatic interactions between single nucleotide polymorphism (SNP) markers have also increased, thus enabling the identification of more candidate markers along with GWAS on single-variant-additive-effect. Here, we focused on the identification of candidate markers associated with flowering time in soybean (Glycine max). A large population of 2,662 cultivated soybean accessions was genotyped using the 180k Axiom® SoyaSNP array, and the genomic architecture of these accessions was investigated to confirm the population structure. Then, GWAS was conducted to evaluate the association between SNP markers and flowering time. A total of 93 significant SNP markers were detected within 59 significant genes, including E1 and E3, which are the main determinants of flowering time. Based on the GWAS results, multilocus epistatic interactions were examined between the significant and non-significant SNP markers. Two significant and 16 non-significant SNP markers were discovered as candidate markers affecting flowering time via interactions with each other. These 18 candidate SNP markers mapped to 18 candidate genes including E1 and E3, and the 18 candidate genes were involved in six major flowering pathways. Although further biological validation is needed, our results provide additional information on the existing flowering time markers and present another option to marker-assisted breeding programs for regulating flowering time of soybean.