Ten Years of the Maize Nested Association Mapping Population: Impact, Limitations, and Future Directions

Early developmental shifts in root exudation profiles of five Zea mays L. genotypes

Root exudates impact soil-plant-microbe interactions and play important roles in ecosystem functioning and plant growth. During early plant development the root rhizosphere may change drastically. For maize (Zea mays L.), one of the world's most important crop species, little is known about root exudation patterns during early plant development. We determined abundance and composition of root exudation among maize genotypes from five inbred lines across three early plant development stages (Emergence, V1-2, and V3-4). We characterized the exudates for non-purgeable organic carbon and performed non-targeted metabolomics with high-performance liquid chromatography tandem mass spectrometry (HPLC-MS/MS). Across all genotypes, plant development stage had a significant influence on both abundance and composition of exudates. Exudation rates (mg C per cm2 root area d-1) were highest in the emergence stage and logarithmically decreased with plant development. In the emergence stage, the roots released predominantly sugars (most indicative: glucose and fructose) and the metabolite richness was generally higher than in later stages. Secondary compounds (e.g. phenolics, benzoxazinoids, or mucilage) increased significantly in later development stages. Differences in the composition of exudates between genotypes may be related to their respective development strategies, with genotypes accumulating more biomass releasing relatively more compounds related to root establishment (growth and rhizosphere development, e.g. mucilage, fatty and organic acids) and slower developing genotypes relatively more metabolites related to maintenance and defense (e.g. phenolics). Our results shed light onto the early dynamics of maize root exudation and rhizosphere establishment, over a phenotypical spectrum of genotypes.

A deficient CP24 allele defines variation for dynamic nonphotochemical quenching and photosystem II efficiency in maize

Maize (Zea mays L.) is a global crop species in which CO2 assimilation occurs via the C4 pathway. C4 photosynthesis is typically more efficient than C3 photosynthesis under warm and dry conditions; however, despite this inherent advantage, considerable variation remains in photosynthetic efficiency for C4 species that could be leveraged to benefit crop performance. Here, we investigate the genetic architecture of nonphotochemical quenching (NPQ) and photosystem II (PSII) efficiency using a combination of high-throughput phenotyping and quantitative trait loci (QTL) mapping in a field-grown Multi-parent Advanced Generation Inter-Cross (MAGIC) mapping population. QTL mapping was followed by the identification of putative candidate genes using a combination of genomics, transcriptomics, protein biochemistry, and targeted physiological phenotyping. We identified four genes with a putative causal role in the observed QTL effects. The highest confidence causal gene was found for a large effect QTL for photosynthetic efficiency on chromosome 10, which was underpinned by allelic variation in the expression of the minor PSII antenna protein light harvesting complex photosystem II subunit (LHCB6 or CP24), mainly driven by poor expression associated with the haplotype of the F7 founder line. The historical role of this line in breeding for early flowering time may suggest that the presence of this deficient allele could be enriched in temperate maize germplasm. These findings advance our understanding of the genetic basis of NPQ and PSII efficiency in C4 plants and highlight the potential for breeding strategies aimed at optimizing photosynthetic efficiency in maize.

Characterization of a G. max × G. soja nested association mapping population and identification of loci controlling seed composition traits from wild soybean

Wild soybean (Glycine soja Siebold & Zucc.) has valuable genetic diversity for improved disease resistance, stress tolerance, seed protein content and seed sulfur-containing amino acid concentrations. Many studies have reported loci controlling seed composition traits based on cultivated soybean populations, but wild soybean has been largely overlooked. In this study, a nested association mapping (NAM) population consisting of 10 families and 1107 recombinant inbred lines was developed by crossing 10 wild accessions with the common cultivar NC-Raleigh. Seed composition of the F6 generation grown at two locations was phenotyped, and genetic markers were identified for each line. The average number of recombination events in the wild soybean-derived population was significantly higher than that in the cultivated soybean-derived population, which resulted in a higher resolution for QTL mapping. Segregation bias in almost all NAM families was significantly biased toward the alleles of the wild soybean parent. Through single-family linkage mapping and association analysis of the entire NAM population, new QTLs with positive allele effects were identified from wild parents, including 5, 6, 18, 9, 16, 17 and 20 for protein content, oil content, total protein and oil content, methionine content, cysteine content, lysine content and threonine content, respectively. Candidate genes associated with these traits were identified based on gene annotations and gene expression levels in different tissues. This is the first study to reveal the genetic characteristics of wild soybean-derived populations, landscapes and the extent of effects of QTLs and candidate genes controlling traits from different wild soybean parents.

Regulatory variation controlling architectural pleiotropy in maize

An early event in plant organogenesis is establishment of a boundary between the stem cell containing meristem and differentiating lateral organ. In maize (Zea mays), evidence suggests a common gene network functions at boundaries of distinct organs and contributes to pleiotropy between leaf angle and tassel branch number, two agronomic traits. To uncover regulatory variation at the nexus of these two traits, we use regulatory network topologies derived from specific developmental contexts to guide multivariate genome-wide association analyses. In addition to defining network plasticity around core pleiotropic loci, we identify new transcription factors that contribute to phenotypic variation in canopy architecture, and structural variation that contributes to cis-regulatory control of pleiotropy between tassel branching and leaf angle across maize diversity. Results demonstrate the power of informing statistical genetics with context-specific developmental networks to pinpoint pleiotropic loci and their cis-regulatory components, which can be used to fine-tune plant architecture for crop improvement.

Using genome-wide associations and host-by-pathogen predictions to identify allelic interactions that control disease resistance

Characterizing the molecular mechanisms underlying disease symptom expression has been used to improve human health and disease resistance in crops and animal breeds. Quantitative trait loci and genome-wide association studies (GWAS) are widely used to identify genomic regions that are involved in disease progression. This study extends traditional GWAS significance tests of host and pathogen marker main effects by utilizing dual-genome reaction norm models to evaluate the importance of host-single nucleotide polymorphism (SNP) by pathogen-SNP interactions. Disease symptom severity data from Fusarium ear rot (FER) on maize (Zea mays L.) is used to demonstrate the use of both genomes in genomic selection models for breeding and the identification of loci that interact across organisms to impact FER disease development. Dual genome prediction models improved heritability estimates, error variances, and model accuracy while providing predictions for host-by-pathogen interactions that may be used to test the significance of SNP-SNP interactions. Independent GWAS for maize and Fusarium populations identified significantly associated loci and predictions that were used to evaluate the importance of interactions using two different association tests. Predictions from dual genome models were used to evaluate the significance of the SNP-SNP interactions that may be associated with population structure or polygenic effects. As well, association tests incorporating host and pathogen markers in models that also included genomic relationship matrices were used to account for population structure. Subsequent evaluation of protein-protein interactions from candidate genes near the interacting SNPs provides a further in silico evaluation method to expedite the identification of interacting genes.

Maize inbreds show allelic variation for diel transcription patterns

Circadian entrainment and external cues can cause gene transcript abundance to oscillate throughout the day, and these patterns of diel transcript oscillation vary across genes and plant species. Less is known about within-species allelic variation for diel patterns of transcript oscillation, or about how regulatory sequence variation influences diel transcription patterns. In this study, we evaluated diel transcript abundance for 24 diverse maize inbred lines. We observed extensive natural variation in diel transcription patterns, with two-fold variation in the number of genes that oscillate over the course of the day. A convolutional neural network trained to predict oscillation from promoter sequence identified sequences previously reported as binding motifs for known circadian clock genes in other plant systems. Genes showing diel transcription patterns that cosegregate with promoter sequence haplotypes are enriched for associations with photoperiod sensitivity and may have been indirect targets of selection as maize was adapted to longer day lengths at higher latitudes. These findings support the idea that cis-regulatory sequence variation influences patterns of gene expression, which in turn can have effects on phenotypic plasticity and local adaptation.

Expanding the BonnMu sequence-indexed repository of transposon induced maize (Zea mays L.) mutations in dent and flint germplasm

The BonnMu resource is a transposon tagged mutant collection designed for functional genomics studies in maize. To expand this resource, we crossed an active Mutator (Mu) stock with dent (B73, Co125) and flint (DK105, EP1, and F7) germplasm, resulting in the generation of 8064 mutagenized BonnMu F2-families. Sequencing of these Mu-tagged families revealed 425 924 presumptive heritable Mu insertions affecting 36 612 (83%) of the 44 303 high-confidence gene models of maize (B73v5). On average, we observed 12 Mu insertions per gene (425 924 total insertions/36 612 affected genes) and 53 insertions per BonnMu F2-family (425 924 total insertions/8064 families). Mu insertions and photos of seedling phenotypes from segregating BonnMu F2-families can be accessed through the Maize Genetics and Genomics Database (MaizeGDB). Downstream examination via the automated Mutant-seq Workflow Utility (MuWU) identified 94% of the presumptive germinal insertion sites in genic regions and only a small fraction of 6% inserting in non-coding intergenic sequences of the genome. Consistently, Mu insertions aligned with gene-dense chromosomal arms. In total, 42% of all BonnMu insertions were located in the 5' untranslated region of genes, corresponding to accessible chromatin. Furthermore, for 38% of the insertions (163 843 of 425 924 total insertions) Mu1, Mu8 and MuDR were confirmed to be the causal Mu elements. Our publicly accessible European BonnMu resource has archived insertions covering two major germplasm groups, thus facilitating both forward and reverse genetics studies.

Differences in stomatal sensitivity to CO2 and light influence variation in water use efficiency and leaf carbon isotope composition in two genotypes of the C4 plant Zea mays

The ratio of net CO2 uptake (Anet) and stomatal conductance (gs) is an intrinsic measurement of leaf water use efficiency (WUEi); however, its measurement can be challenging for large phenotypic screens. Measurements of the leaf carbon isotope composition (δ13Cleaf) may be a scalable tool to approximate WUEi for screening because it in part reflects the competing influences of Anet and gs on the CO2 partial pressure (pCO2) inside the leaf over time. However, in C4 photosynthesis, the CO2-concentrating mechanism complicates the relationship between δ13Cleaf and WUEi. Despite this complicated relationship, several studies have shown genetic variation in δ13Cleaf across C4 plants. Yet there has not been a clear demonstration of if Anet or gs are the causal mechanisms controlling WUEi and δ13Cleaf. Our approach was to characterize leaf photosynthetic traits of two Zea mays recombinant inbred lines (Z007E0067 and Z007E0150) which consistently differ for δ13Cleaf even though they have minimal confounding genetic differences. We demonstrate that these two genotypes contrasted in WUEi driven by differences in the speed of stomatal responses to changes in pCO2 and light that lead to unproductive leaf water loss. These findings provide support that differences in δ13Cleaf in closely related genotypes do reflect greater WUEi and further suggest that differences in stomatal kinetic response to changing environmental conditions is a key target to improve WUEi.

Genome-wide association mapping in exotic × Canadian elite crosses: mining beneficial alleles for agronomic and seed composition traits in soybean

Given the narrow genetic base of North American soybean germplasm, which originates from approximately 35 ancestral lines, discovering and introducing useful diversity for key traits in exotic germplasm could potentially enhance diversity in the current elite gene pool. This study explores the potential of exotic germplasm to enhance yield and agronomic traits in the University of Guelph soybean germplasm. We utilized a nested association mapping (NAM) design to develop a population (n = 294) composed of crosses of high-yielding Canadian elite cultivar, OAC Bruton, with four high-yielding exotic lines developed at USDA (Urbana, IL), and we mapped the genetic architecture of agronomic and seed composition traits using association mapping methods. The analysis across three Southwestern Ontario environments revealed seven unique genomic regions underlying agronomic traits and four for seed composition traits, with both desirable and undesirable alleles from the exotic parents. Notably, a region on chromosome 10, co-locating to the E2 maturity locus, was found to be associated with seed yield and maturity. The allele that increased yield by 166 kg/ha was contributed by all exotic parents and was absent in the Canadian-adapted parent. The study underscores the potential of using exotic germplasm to introduce novel genetic diversity into the Canadian elite soybean breeding pool. By identifying exotic-derived beneficial alleles, our findings offer a pathway for enhancing agronomic traits in Canadian soybeans with novel exotic diversity.

QTL mapping and genome-wide association analysis reveal genetic loci and candidate gene for resistance to gray leaf spot in tropical and subtropical maize germplasm

KEY MESSAGE

Using QTL mapping and GWAS, two candidate genes (Zm00001d051039 and Zm00001d051147) were consistently identified across the three different environments and BLUP values. GWAS analysis identified the candidate gene, Zm00001d044845. These genes were subsequently validated to exhibit a significant association with maize gray leaf spot (GLS) resistance. Gray leaf spot (GLS) is a major foliar disease of maize (Zea mays L.) that causes significant yield losses worldwide. Understanding the genetic mechanisms underlying gray leaf spot resistance is crucial for breeding high-yielding and disease-resistant varieties. In this study, eight tropical and subtropical germplasms were crossed with the temperate germplasm Ye107 to develop a nested association mapping (NAM) population comprising 1,653 F2:8 RILs, consisting of eight recombinant inbred line (RIL) subpopulations, using the single-seed descent method. The NAM population was evaluated for GLS resistance in three different environments, and genotyping by sequencing of the NAM population generated 593,719 high-quality single-nucleotide polymorphisms (SNPs). Linkage analysis and genome-wide association studies (GWASs) were conducted to identify candidate genes regulating GLS resistance in maize. Both analyses identified 25 QTLs and 149 SNPs that were significantly associated with GLS resistance. Candidate genes were screened 20 Kb upstream and downstream of the significant SNPs, and three novel candidate genes (Zm00001d051039, Zm00001d051147, and Zm00001d044845) were identified. Zm00001d051039 and Zm00001d051147 were located on chromosome 4 and co-localized in both linkage (qGLS4-1 and qGLS4-2) and GWAS analyses. SNP-138,153,206 was located 0.499 kb downstream of the candidate gene Zm00001d051039, which encodes the protein IN2-1 homolog B, a homolog of glutathione S-transferase (GST). GSTs and protein IN2-1 homolog B scavenge reactive oxygen species under various stress conditions, and GSTs are believed to protect plants from a wide range of biotic and abiotic stresses by detoxifying reactive electrophilic compounds. Zm00001d051147 encodes a probable beta-1,4-xylosyltransferase involved in the biosynthesis of xylan in the cell wall, enhancing resistance. SNP-145,813,215 was located 2.69 kb downstream of the candidate gene. SNP-5,043,412 was consistently identified in three different environments and BLUP values and was located 8.788 kb downstream of the candidate gene Zm00001d044845 on chromosome 9. Zm00001d044845 encodes the U-box domain-containing protein 4 (PUB4), which is involved in regulating plant immunity. qRT-PCR analysis showed that the relative expression levels of the three candidate genes were significantly upregulated in the leaves of the TML139 (resistant) parent, indicating that these three candidate genes could be associated with resistance to GLS. The findings of this study are significant for marker-assisted breeding aimed at enhancing resistance to GLS in maize and lay the foundation for further elucidation of the genetic mechanisms underlying resistance to gray leaf spot in maize and breeding of new disease-resistant varieties.

Bacterial bioluminescence is an important regulator of multitrophic interactions in the soil

Enormous efforts have been made to understand the functions of bioluminescence; however, its relevance in soil ecosystems has barely been investigated. In addition, our understanding of the biological relevance of bioluminescence is hampered by the scarcity of tools to genetically manipulate this trait. Using the symbionts of entomopathogenic nematodes, Photorhabdus bacteria, we show that bioluminescence plays important regulatory roles in multitrophic interactions in the soil. Through genetic modifications and exploiting natural variability, we provide direct evidence for the multifunctional nature of bioluminescence. It regulates abiotic and biotic stress resistance, impacts other trophic levels, including nematodes, insects, and plants, and contributes to symbiosis. Our study contributes to understanding the factors that have driven the evolution and maintenance of this trait in belowground ecosystems.

Comprehensive identification of genomic and environmental determinants of phenotypic plasticity in maize

Maize phenotypes are plastic, determined by the complex interplay of genetics and environmental variables. Uncovering the genes responsible and understanding how their effects change across a large geographic region are challenging. In this study, we conducted systematic analysis to identify environmental indices that strongly influence 19 traits (including flowering time, plant architecture, and yield component traits) measured in the maize nested association mapping (NAM) population grown in 11 environments. Identified environmental indices based on day length, temperature, moisture, and combinations of these are biologically meaningful. Next, we leveraged a total of more than 20 million SNP and SV markers derived from recent de novo sequencing of the NAM founders for trait prediction and dissection. When combined with identified environmental indices, genomic prediction enables accurate performance predictions. Genome-wide association studies (GWASs) detected genetic loci associated with the plastic response to the identified environmental indices for all examined traits. By systematically uncovering the major environmental and genomic factors underlying phenotypic plasticity in a wide variety of traits and depositing our results as a track on the MaizeGDB genome browser, we provide a community resource as well as a comprehensive analytical framework to facilitate continuing complex trait dissection and prediction in maize and other crops. Our findings also provide a conceptual framework for the genetic architecture of phenotypic plasticity by accommodating two alternative models, regulatory gene model and allelic sensitivity model, as special cases of a continuum.

Differences in activity and stability drive transposable element variation in tropical and temperate maize

Much of the profound interspecific variation in genome content has been attributed to transposable elements (TEs). To explore the extent of TE variation within species, we developed an optimized open-source algorithm, panEDTA, to de novo annotate TEs in a pangenome context. We then generated a unified TE annotation for a maize pangenome derived from 26 reference-quality genomes, which reveals an excess of 35.1 Mb of TE sequences per genome in tropical maize relative to temperate maize. A small number (n = 216) of TE families, mainly LTR retrotransposons, drive these differences. Evidence from the methylome, transcriptome, LTR age distribution, and LTR insertional polymorphisms reveals that 64.7% of the variability is contributed by LTR families that are young, less methylated, and more expressed in tropical maize, whereas 18.5% is driven by LTR families with removal or loss in temperate maize. Additionally, we find enrichment for Young LTR families adjacent to nucleotide-binding and leucine-rich repeat (NLR) clusters of varying copy number across lines, suggesting TE activity may be associated with disease resistance in maize.

Genomic investigation of 18,421 lines reveals the genetic architecture of rice

Understanding how numerous quantitative trait loci (QTL) shape phenotypic variation is an important question in genetics. To address this, we established a permanent population of 18,421 (18K) rice lines with reduced population structure. We generated reference-level genome assemblies of the founders and genotyped all 18K-rice lines through whole-genome sequencing. Through high-resolution mapping, 96 high-quality candidate genes contributing to variation in 16 traits were identified, including OsMADS22 and OsFTL1 verified as causal genes for panicle number and heading date, respectively. We identified epistatic QTL pairs and constructed a genetic interaction network with 19 genes serving as hubs. Overall, 170 masking epistasis pairs were characterized, serving as an important factor contributing to genetic background effects across diverse varieties. The work provides a basis to guide grain yield and quality improvements in rice.

Characterization of adaptation mechanisms in sorghum using a multireference back-cross nested association mapping design and envirotyping

Identifying the genetic factors impacting the adaptation of crops to environmental conditions is of key interest for conservation and selection purposes. It can be achieved using population genomics, and evolutionary or quantitative genetics. Here we present a sorghum multireference back-cross nested association mapping population composed of 3,901 lines produced by crossing 24 diverse parents to 3 elite parents from West and Central Africa-back-cross nested association mapping. The population was phenotyped in environments characterized by differences in photoperiod, rainfall pattern, temperature levels, and soil fertility. To integrate the multiparental and multi-environmental dimension of our data we proposed a new approach for quantitative trait loci (QTL) detection and parental effect estimation. We extended our model to estimate QTL effect sensitivity to environmental covariates, which facilitated the integration of envirotyping data. Our models allowed spatial projections of the QTL effects in agro-ecologies of interest. We utilized this strategy to analyze the genetic architecture of flowering time and plant height, which represents key adaptation mechanisms in environments like West Africa. Our results allowed a better characterization of well-known genomic regions influencing flowering time concerning their response to photoperiod with Ma6 and Ma1 being photoperiod-sensitive and the region of possible candidate gene Elf3 being photoperiod-insensitive. We also accessed a better understanding of plant height genetic determinism with the combined effects of phenology-dependent (Ma6) and independent (qHT7.1 and Dw3) genomic regions. Therefore, we argue that the West and Central Africa-back-cross nested association mapping and the presented analytical approach constitute unique resources to better understand adaptation in sorghum with direct application to develop climate-smart varieties.

Evaluating and Predicting the Performance of Sorghum Lines in an Elite by Exotic Backcross-Nested Association Mapping Population

Maintaining or introducing genetic diversity into plant breeding programs is necessary for continual genetic gain; however, diversity at the cost of reduced performance is not something sought by breeders. To this end, backcross-nested association mapping (BC-NAM) populations, in which the recurrent parent is an elite line, can be employed as a strategy to introgress diversity from unadapted accessions while maintaining agronomic performance. This study evaluates (i) the hybrid performance of sorghum lines from 18 BC1-NAM families and (ii) the potential of genomic prediction to screen lines from BC1-NAM families for hybrid performance prior to phenotypic evaluation. Despite the diverse geographical origins and agronomic performance of the unadapted parents for BC1-NAM families, many BC1-derived lines performed significantly better in the hybrid trials than the elite recurrent parent, R.Tx436. The genomic prediction accuracies for grain yield, plant height, and days to mid-anthesis were acceptable, but the prediction accuracies for plant height were lower than expected. While the prediction accuracies increased when including more individuals in the training set, improvements tended to plateau between two and five lines per family, with larger training sets being required for more complex traits such as grain yield. Therefore, genomic prediction models can be optimized in a large BC1-NAM population with a relatively low fraction of individuals needing to be evaluated. These results suggest that genomic prediction is an effective method of pre-screening lines within BC1-NAM families prior to evaluation in extensive hybrid field trials.

Leveraging prior biological knowledge improves prediction of tocochromanols in maize grain

With an essential role in human health, tocochromanols are mostly obtained by consuming seed oils; however, the vitamin E content of the most abundant tocochromanols in maize (Zea mays L.) grain is low. Several large-effect genes with cis-acting variants affecting messenger RNA (mRNA) expression are mostly responsible for tocochromanol variation in maize grain, with other relevant associated quantitative trait loci (QTL) yet to be fully resolved. Leveraging existing genomic and transcriptomic information for maize inbreds could improve prediction when selecting for higher vitamin E content. Here, we first evaluated a multikernel genomic best linear unbiased prediction (MK-GBLUP) approach for modeling known QTL in the prediction of nine tocochromanol grain phenotypes (12-21 QTL per trait) within and between two panels of 1,462 and 242 maize inbred lines. On average, MK-GBLUP models improved predictive abilities by 7.0-13.6% when compared with GBLUP. In a second approach with a subset of 545 lines from the larger panel, the highest average improvement in predictive ability relative to GBLUP was achieved with a multi-trait GBLUP model (15.4%) that had a tocochromanol phenotype and transcript abundances in developing grain for a few large-effect candidate causal genes (1-3 genes per trait) as multiple response variables. Taken together, our study illustrates the enhancement of prediction models when informed by existing biological knowledge pertaining to QTL and candidate causal genes.

Leveraging Single-Cell Populations to Uncover the Genetic Basis of Complex Traits

The ease and throughput of single-cell genomics have steadily improved, and its current trajectory suggests that surveying single-cell populations will become routine. We discuss the merger of quantitative genetics with single-cell genomics and emphasize how this synergizes with advantages intrinsic to plants. Single-cell population genomics provides increased detection resolution when mapping variants that control molecular traits, including gene expression or chromatin accessibility. Additionally, single-cell population genomics reveals the cell types in which variants act and, when combined with organism-level phenotype measurements, unveils which cellular contexts impact higher-order traits. Emerging technologies, notably multiomics, can facilitate the measurement of both genetic changes and genomic traits in single cells, enabling single-cell genetic experiments. The implementation of single-cell genetics will advance the investigation of the genetic architecture of complex molecular traits and provide new experimental paradigms to study eukaryotic genetics.

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

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

Plant biochemical genetics in the multiomics era

Our understanding of plant biology has been revolutionized by modern genetics and biochemistry. However, biochemical genetics can be traced back to the foundation of Mendelian genetics; indeed, one of Mendel's milestone discoveries of seven characteristics of pea plants later came to be ascribed to a mutation in a starch branching enzyme. Here, we review both current and historical strategies for the elucidation of plant metabolic pathways and the genes that encode their component enzymes and regulators. We use this historical review to discuss a range of classical genetic phenomena including epistasis, canalization, and heterosis as viewed through the lens of contemporary high-throughput data obtained via the array of approaches currently adopted in multiomics studies.

Genomics and phenomics enabled prebreeding improved early-season chilling tolerance in Sorghum

In temperate climates, earlier planting of tropical-origin crops can provide longer growing seasons, reduce water loss, suppress weeds, and escape post-flowering drought stress. However, chilling sensitivity of sorghum, a tropical-origin cereal crop, limits early planting, and over 50 years of conventional breeding has been stymied by coinheritance of chilling tolerance (CT) loci with undesirable tannin and dwarfing alleles. In this study, phenomics and genomics-enabled approaches were used for prebreeding of sorghum early-season CT. Uncrewed aircraft systems (UAS) high-throughput phenotyping platform tested for improving scalability showed moderate correlation between manual and UAS phenotyping. UAS normalized difference vegetation index values from the chilling nested association mapping population detected CT quantitative trait locus (QTL) that colocalized with manual phenotyping CT QTL. Two of the 4 first-generation Kompetitive Allele Specific PCR (KASP) molecular markers, generated using the peak QTL single nucleotide polymorphisms (SNPs), failed to function in an independent breeding program as the CT allele was common in diverse breeding lines. Population genomic fixation index analysis identified SNP CT alleles that were globally rare but common to the CT donors. Second-generation markers, generated using population genomics, were successful in tracking the donor CT allele in diverse breeding lines from 2 independent sorghum breeding programs. Marker-assisted breeding, effective in introgressing CT allele from Chinese sorghums into chilling-sensitive US elite sorghums, improved early-planted seedling performance ratings in lines with CT alleles by up to 13-24% compared to the negative control under natural chilling stress. These findings directly demonstrate the effectiveness of high-throughput phenotyping and population genomics in molecular breeding of complex adaptive traits.

Detecting directional and non-directional epistasis in bi-parental populations using genomic data

Epistasis, commonly defined as interaction effects between alleles of different loci, is an important genetic component of the variation of phenotypic traits in natural and breeding populations. In addition to its impact on variance, epistasis can also affect the expected performance of a population and is then referred to as directional epistasis. Before the advent of genomic data, the existence of epistasis (both directional and non-directional) was investigated based on complex and expensive mating schemes involving several generations evaluated for a trait of interest. In this study, we propose a methodology to detect the presence of epistasis based on simple inbred biparental populations, both genotyped and phenotyped, ideally along with their parents. Thanks to genomic data, parental proportions as well as shared parental proportions between inbred individuals can be estimated. They allow the evaluation of epistasis through a test of the expected performance for directional epistasis or the variance of genetic values. This methodology was applied to two large multiparental populations, i.e. the American maize and soybean nested association mapping populations, evaluated for different traits. Results showed significant epistasis, especially for the test of directional epistasis, e.g. the increase in anthesis to silking interval observed in most maize inbred progenies or the decrease in grain yield observed in several soybean inbred progenies. In general, the effects detected suggested that shuffling allelic associations of both elite parents had a detrimental effect on the performance of their progeny. This methodology is implemented in the EpiTest R-package and can be applied to any bi/multiparental inbred population evaluated for a trait of interest.

Multi-parental fungal mapping population study to detect genomic regions associated with Pyrenophora teres f. teres virulence

In recent years multi-parental mapping populations (MPPs) have been widely adopted in many crops to detect quantitative trait loci (QTLs) as this method can compensate for the limitations of QTL analyses using bi-parental mapping populations. Here we report the first multi-parental nested association mapping (MP-NAM) population study used to detect genomic regions associated with host-pathogenic interactions. MP-NAM QTL analyses were conducted on 399 Pyrenophora teres f. teres individuals using biallelic, cross-specific and parental QTL effect models. A bi-parental QTL mapping study was also conducted to compare the power of QTL detection between bi-parental and MP-NAM populations. Using MP-NAM with 399 individuals detected a maximum of eight QTLs with a single QTL effect model whilst only a maximum of five QTLs were detected with an individual bi-parental mapping population of 100 individuals. When reducing the number of isolates in the MP-NAM to 200 individuals the number of QTLs detected remained the same for the MP-NAM population. This study confirms that MPPs such as MP-NAM populations can be successfully used in detecting QTLs in haploid fungal pathogens and that the power of QTL detection with MPPs is greater than with bi-parental mapping populations.

Natural variation of respiration-related traits in plants

Plant respiration is one of the greatest global metabolic fluxes, but rates of respiration vary massively both within different cell types as well as between different individuals and different species. Whilst this is well known, few studies have detailed population-level variation of respiration until recently. The last 20 years have seen a renaissance in studies of natural variance. In this review, we describe how experimental breeding populations and collections of large populations of accessions can be used to determine the genetic architecture of plant traits. We further detail how these approaches have been used to study the rate of respiration per se as well as traits that are intimately associated with respiration. The review highlights specific breakthroughs in these areas but also concludes that the approach should be more widely adopted in the study of respiration per se as opposed to the more frequently studied respiration-related traits.

Pan-genome inversion index reveals evolutionary insights into the subpopulation structure of Asian rice

Understanding and exploiting genetic diversity is a key factor for the productive and stable production of rice. Here, we utilize 73 high-quality genomes that encompass the subpopulation structure of Asian rice (Oryza sativa), plus the genomes of two wild relatives (O. rufipogon and O. punctata), to build a pan-genome inversion index of 1769 non-redundant inversions that span an average of ~29% of the O. sativa cv. Nipponbare reference genome sequence. Using this index, we estimate an inversion rate of ~700 inversions per million years in Asian rice, which is 16 to 50 times higher than previously estimated for plants. Detailed analyses of these inversions show evidence of their effects on gene expression, recombination rate, and linkage disequilibrium. Our study uncovers the prevalence and scale of large inversions (≥100 bp) across the pan-genome of Asian rice and hints at their largely unexplored role in functional biology and crop performance.

Parental Inbred Lines of the Nested Association Mapping (NAM) Population of Corn Show Sources of Resistance to Tar Spot in Northern Indiana

Tar spot is a major foliar disease of corn caused by the obligate fungal pathogen Phyllachora maydis, first identified in Indiana in 2015. Under conducive weather conditions, P. maydis causes significant yield losses in the United States and other countries, constituting a major threat to corn production. Relatively little is known about resistance to tar spot other than a major quantitative gene that was identified in tropical maize lines. To test for additional sources of resistance against populations of P. maydis in North America, 26 parental inbred lines of the nested associated mapping (NAM) population were evaluated for tar spot resistance in Indiana in replicated field trials under natural infection for 3 years. Tar spot disease severity was scored visually using a 0-to-100% scale. Maximum disease severity (MDS) for tar spot scoring at reproductive growth stage ranged from 0 to 48.3%, with 0% being most resistant and 48.3% being most susceptible. Nine inbred lines were resistant to P. maydis with MDS ranging from 0 to 5.0%, six were moderately resistant (5.2 to 10.6% MDS), two were moderately susceptible (11.7 to 26.0% MDS), and the remaining eight inbred lines were rated as susceptible (30.0 to 48.3% MDS). There was some variability between years, due to higher disease pressure after 2019. Inbred B73, the common parent of the NAM populations, was rated as susceptible, with MDS of 30.0%. The nine highly resistant lines provide a potential source of new genes for genetic analysis and mapping of tar spot resistance in corn.

Co-expression pan-network reveals genes involved in complex traits within maize pan-genome

BACKGROUND

With the advances in the high throughput next generation sequencing technologies, genome-wide association studies (GWAS) have identified a large set of variants associated with complex phenotypic traits at a very fine scale. Despite the progress in GWAS, identification of genotype-phenotype relationship remains challenging in maize due to its nature with dozens of variants controlling the same trait. As the causal variations results in the change in expression, gene expression analyses carry a pivotal role in unraveling the transcriptional regulatory mechanisms behind the phenotypes.

RESULTS

To address these challenges, we incorporated the gene expression and GWAS-driven traits to extend the knowledge of genotype-phenotype relationships and transcriptional regulatory mechanisms behind the phenotypes. We constructed a large collection of gene co-expression networks and identified more than 2 million co-expressing gene pairs in the GWAS-driven pan-network which contains all the gene-pairs in individual genomes of the nested association mapping (NAM) population. We defined four sub-categories for the pan-network: (1) core-network contains the highest represented ~ 1% of the gene-pairs, (2) near-core network contains the next highest represented 1-5% of the gene-pairs, (3) private-network contains ~ 50% of the gene pairs that are unique to individual genomes, and (4) the dispensable-network contains the remaining 50-95% of the gene-pairs in the maize pan-genome. Strikingly, the private-network contained almost all the genes in the pan-network but lacked half of the interactions. We performed gene ontology (GO) enrichment analysis for the pan-, core-, and private- networks and compared the contributions of variants overlapping with genes and promoters to the GWAS-driven pan-network.

CONCLUSIONS

Gene co-expression networks revealed meaningful information about groups of co-regulated genes that play a central role in regulatory processes. Pan-network approach enabled us to visualize the global view of the gene regulatory network for the studied system that could not be well inferred by the core-network alone.

Genome-wide association studies revealed complex genetic architecture and breeding perspective of maize ear traits

BACKGROUND

Maize (Zea Mays) is one of the world's most important crops. Hybrid maize lines resulted a major improvement in corn production in the previous and current centuries. Understanding the genetic mechanisms of the corn production associated traits greatly facilitate the development of superior hybrid varieties.

RESULT

In this study, four ear traits associated with corn production of Nested Association Mapping (NAM) population were analyzed using a full genetic model, and further, optimal genotype combinations and total genetic effects of current best lines, superior lines, and superior hybrids were predicted for each of the traits at four different locations. The analysis identified 21-34 highly significant SNPs (-log₁₀P > 5), with an estimated total heritability of 37.31-62.34%, while large contributions to variations was due to dominance, dominance-related epistasis, and environmental interaction effects ([Formula: see text] 14.06% ~ 49.28%), indicating these factors contributed significantly to phenotypic variations of the ear traits. Environment-specific genetic effects were also discovered to be crucial for maize ear traits. There were four SNPs found for three ear traits: two for ear length and weight, and two for ear row number and length. Using the Enumeration method and the stepwise tuning technique, optimum multi-locus genotype combinations for superior lines were identified based on the information obtained from GWAS.

CONCLUSIONS

Predictions of genetic breeding values showed that different genotype combinations in different geographical regions may be better, and hybrid-line variety breeding with homozygote and heterozygote genotype combinations may have a greater potential to improve ear traits.

Functional Study of Lipoxygenase-Mediated Resistance against Fusarium verticillioides and Aspergillus flavus Infection in Maize

Mycotoxin contamination of maize kernels by fungal pathogens like Fusarium verticillioides and Aspergillus flavus is a chronic global challenge impacting food and feed security, health, and trade. Maize lipoxygenase genes (ZmLOXs) synthetize oxylipins that play defense roles and govern host-fungal interactions. The current study investigated the involvement of ZmLOXs in maize resistance against these two fungi. A considerable intraspecific genetic and transcript variability of the ZmLOX family was highlighted by in silico analysis comparing publicly available maize pan-genomes and pan-transcriptomes, respectively. Then, phenotyping and expression analysis of ZmLOX genes along with key genes involved in oxylipin biosynthesis were carried out in a maize mutant carrying a Mu transposon insertion in the ZmLOX4 gene (named UFMulox4) together with Tzi18, Mo17, and W22 inbred lines at 3- and 7-days post-inoculation with F. verticillioides and A. flavus. Tzi18 showed the highest resistance to the pathogens coupled with the lowest mycotoxin accumulation, while UFMulox4 was highly susceptible to both pathogens with the most elevated mycotoxin content. F. verticillioides inoculation determined a stronger induction of ZmLOXs and maize allene oxide synthase genes as compared to A. flavus. Additionally, oxylipin analysis revealed prevalent linoleic (18:2) peroxidation by 9-LOXs, the accumulation of 10-oxo-11-phytoenoic acid (10-OPEA), and triglyceride peroxidation only in F. verticillioides inoculated kernels of resistant genotypes.

Identification of two new QTLs of maize (Zea mays L.) underlying kernel row number using the HNAU-NAM1 population

BACKGROUND

Maize kernel row number (KRN) is one of the most important yield traits and has changed greatly during maize domestication and selection. Elucidating the genetic basis of KRN will be helpful to improve grain yield in maize.

RESULTS

Here, we measured KRN in four environments using a nested association mapping (NAM) population named HNAU-NAM1 with 1,617 recombinant inbred lines (RILs) that were derived from 12 maize inbred lines with a common parent, GEMS41. Then, five consensus quantitative trait loci (QTLs) distributing on four chromosomes were identified in at least three environments along with the best linear unbiased prediction (BLUP) values by the joint linkage mapping (JLM) method. These QTLs were further validated by the separate linkage mapping (SLM) and genome-wide association study (GWAS) methods. Three KRN genes cloned through the QTL assay were found in three of the five consensus QTLs, including qKRN1.1, qKRN2.1 and qKRN4.1. Two new QTLs of KRN, qKRN4.2 and qKRN9.1, were also identified. On the basis of public RNA-seq and genome annotation data, five genes highly expressed in ear tissue were considered candidate genes contributing to KRN.

CONCLUSIONS

This study carried out a comprehensive analysis of the genetic architecture of KRN by using a new NAM population under multiple environments. The present results provide solid information for understanding the genetic components underlying KRN and candidate genes in qKRN4.2 and qKRN9.1. Single-nucleotide polymorphisms (SNPs) closely linked to qKRN4.2 and qKRN9.1 could be used to improve inbred yield during molecular breeding in maize.

Fine mapping of candidate quantitative trait loci for plant and ear height in a maize nested-association mapping population

Plant height (PH) and ear height (EH) are two important traits in maize (Zea mays L.), as they are closely related to lodging resistance and planting density. Our objectives were to (1) investigate single-nucleotide polymorphisms (SNPs) that are associated with PH and EH for detecting quantitative trait loci (QTL) and new gene that determines PH and EH, (2) explore the value of the QTL in maize breeding, and (3) investigate whether the "triangle heterotic group" theory is applicable for lowering PH and EH to increase yield. Seven inbred female parents were crossed with a common founder male parent Ye 107 to create a nested association mapping (NAM) population. The analysis of phenotypic data on PH and EH revealed wide variation among the parents of the NAM population. Genome-wide association study (GWAS) and high-resolution linkage mapping were conducted using the NAM population, which generated 264,694 SNPs by genotyping-by-sequencing. A total of 105 SNPs and 22 QTL were identified by GWAS and found to be significantly associated with PH and EH. A high-confidence QTL for PH, Qtl-chr1-EP, was identified on chromosome 1 via GWAS and confirmed by linkage analysis in two recombinant inbred line (RIL) populations. Results revealed that the SNP variation in the promoter region of the candidate gene Zm00001d031938, located at Qtl-chr1-EP, which encoded UDP-N-acetylglucosamine-peptide N-acetyl-glucosaminyl-transferase, might decrease PH and EH. Furthermore, the triangle heterotic pattern adopted in maize breeding programs by our team is practicable in selecting high-yield crosses based on the low ratio of EH/PH (EP).

Effects of drought and heat stresses during reproductive stage on pollen germination, yield, and leaf reflectance properties in maize (Zea mays L.)

Drought and heat stresses are the major abiotic stress factors detrimental to maize (Zea mays L.) production. Much attention has been directed toward plant responses to heat or drought stress. However, maize reproductive stage responses to combined heat and drought remain less explored. Therefore, this study aimed to quantify the impact of optimum daytime (30°C, control) and warmer daytime temperatures (35°C, heat stress) on pollen germination, morpho-physiology, and yield potential using two maize genotypes ("Mo17" and "B73") under contrasting soil moisture content, that is, 100% and 40% irrigation during flowering. Pollen germination of both genotypes decreased under combined stresses (42%), followed by heat stress (30%) and drought stress (19%). Stomatal conductance and transpiration were comparable between control and heat stress but significantly decreased under combined stresses (83% and 72%) and drought stress (52% and 47%) compared with the control. Genotype "Mo17" reduced its green leaf area to minimize the water loss, which appears to be one of the adaptive strategies of "Mo17" under stress conditions. The leaf reflectance of both genotypes varied across treatments. Vegetation indices associated with pigments (chlorophyll index of green, chlorophyll index of red edge, and carotenoid index) and plant health (normalized difference red-edge index) were found to be highly sensitive to drought and combined stressors than heat stress. Combined drought and heat stresses caused a significant reduction in yield and yield components in both Mo17 (49%) and B73 (86%) genotypes. The harvest index of genotype "B73" was extremely low, indicating poor partitioning efficiency. At least when it comes to "B73," the cause of yield reduction appears to be the result of reduced sink number rather than the pollen and source size. To the best of our awareness, this is the first study that showed how the leaf-level spectra, yield, and quality parameters respond to the short duration of independent and combined stresses during flowering in inbred maize. Further studies are required to validate the responses of potential traits involving diverse maize genotypes under field conditions. This study suggests the need to develop maize with improved tolerance to combined stresses to sustain production under increasing temperatures and low rainfall conditions.

How Could the Use of Crop Wild Relatives in Breeding Increase the Adaptation of Crops to Marginal Environments?

Alongside the use of fertilizer and chemical control of weeds, pests, and diseases modern breeding has been very successful in generating cultivars that have increased agricultural production several fold in favorable environments. These typically homogeneous cultivars (either homozygous inbreds or hybrids derived from inbred parents) are bred under optimal field conditions and perform well when there is sufficient water and nutrients. However, such optimal conditions are rare globally; indeed, a large proportion of arable land could be considered marginal for agricultural production. Marginal agricultural land typically has poor fertility and/or shallow soil depth, is subject to soil erosion, and often occurs in semi-arid or saline environments. Moreover, these marginal environments are expected to expand with ongoing climate change and progressive degradation of soil and water resources globally. Crop wild relatives (CWRs), most often used in breeding as sources of biotic resistance, often also possess traits adapting them to marginal environments. Wild progenitors have been selected over the course of their evolutionary history to maintain their fitness under a diverse range of stresses. Conversely, modern breeding for broad adaptation has reduced genetic diversity and increased genetic vulnerability to biotic and abiotic challenges. There is potential to exploit genetic heterogeneity, as opposed to genetic uniformity, in breeding for the utilization of marginal lands. This review discusses the adaptive traits that could improve the performance of cultivars in marginal environments and breeding strategies to deploy them.

Functional Allele Validation by Gene Editing to Leverage the Wealth of Genetic Resources for Crop Improvement

Advances in molecular technologies over the past few decades, such as high-throughput DNA marker genotyping, have provided more powerful plant breeding approaches, including marker-assisted selection and genomic selection. At the same time, massive investments in plant genetics and genomics, led by whole genome sequencing, have led to greater knowledge of genes and genetic pathways across plant genomes. However, there remains a gap between approaches focused on forward genetics, which start with a phenotype to map a mutant locus or QTL with the goal of cloning the causal gene, and approaches using reverse genetics, which start with large-scale sequence data and work back to the gene function. The recent establishment of efficient CRISPR-Cas-based gene editing promises to bridge this gap and provide a rapid method to functionally validate genes and alleles identified through studies of natural variation. CRISPR-Cas techniques can be used to knock out single or multiple genes, precisely modify genes through base and prime editing, and replace alleles. Moreover, technologies such as protoplast isolation, in planta transformation, and the use of developmental regulatory genes promise to enable high-throughput gene editing to accelerate crop improvement.

Variation in upstream open reading frames contributes to allelic diversity in maize protein abundance

SignificanceProteins are the machinery which execute essential cellular functions. However, measuring their abundance within an organism can be difficult and resource-intensive. Cells use a variety of mechanisms to control protein synthesis from mRNA, including short open reading frames (uORFs) that lie upstream of the main coding sequence. Ribosomes can preferentially translate uORFs instead of the main coding sequence, leading to reduced translation of the main protein. In this study, we show that uORF sequence variation between individuals can lead to different rates of protein translation and thus variable protein abundances. We also demonstrate that natural variation in uORFs occurs frequently and can be linked to whole-plant phenotypes, indicating that uORF sequence variation likely contributes to plant adaptation.

Genetic dissection of maize plant architecture using a novel nested association mapping population

The leaf angle (LA), plant height (PH), and ear height (EH) are key plant architectural traits influencing maize (Zea mays L.) yield. However, their genetic determinants have not yet been well-characterized. Here, we developed a maize advanced backcross-nested association mapping population in Henan Agricultural University (HNAU-NAM1) comprised of 1,625 BC₁ F₄ /BC₂ F₄ lines. These were obtained by crossing a diverse set of 12 representative inbred lines with the common GEMS41 line, which were then genotyped using the MaizeSNP9.4K array. Genetic diversity and phenotypic distribution analyses showed considerable levels of genetic variation. We obtained 18-88 quantitative trait loci (QTLs) associated with LA, PH, and EH by using three complementary mapping methods, named as separate linkage mapping, joint linkage mapping, and genome-wide association studies. Our analyses enabled the identification of ten QTL hot-spot regions associated with the three traits, which were distributed on nine different chromosomes. We further selected 13 major QTLs that were simultaneously detected by three methods and deduced the candidate genes, of which eight were not reported before. The newly constructed HNAU-NAM1 population in this study will further broaden our insights into understanding of genetic regulation of plant architecture, thus will help to improve maize yield and provide an invaluable resource for maize functional genomics and breeding research.

Back to the wild: mining maize (Zea mays L.) disease resistance using advanced breeding tools

Cultivated modern maize (Zea mays L.) originated through the continuous process of domestication from its wild progenitors. Today, maize is considered as the most important cereal crop which is extensively cultivated in all parts of the world. Maize shows remarkable genotypic and phenotypic diversity which makes it an ideal model species for crop genetic research. However, intensive breeding and artificial selection of desired agronomic traits greatly narrow down the genetic bases of maize. This reduction in genetic diversity among cultivated maize led to increase the chance of more attack of biotic stress as climate changes hampering the maize grain production globally. Maize germplasm requires to integrate both durable multiple-diseases and multiple insect-pathogen resistance through tapping the unexplored resources of maize landraces. Revisiting the landraces seed banks will provide effective opportunities to transfer the resistant genes into the modern cultivars. Here, we describe the maize domestication process and discuss the unique genes from wild progenitors which potentially can be utilized for disease resistant in maize. We also focus on the genetics and disease resistance mechanism of various genes against maize biotic stresses and then considered the different molecular breeding tools for gene transfer and advanced high resolution mapping for gene pyramiding in maize lines. At last, we provide an insight for targeting identified key genes through CRISPR/Cas9 genome editing system to enhance the maize resilience towards biotic stress.

Exploiting High-Throughput Indoor Phenotyping to Characterize the Founders of a Structured B. napus Breeding Population

Phenotyping is considered a significant bottleneck impeding fast and efficient crop improvement. Similar to many crops, Brassica napus, an internationally important oilseed crop, suffers from low genetic diversity, and will require exploitation of diverse genetic resources to develop locally adapted, high yielding and stress resistant cultivars. A pilot study was completed to assess the feasibility of using indoor high-throughput phenotyping (HTP), semi-automated image processing, and machine learning to capture the phenotypic diversity of agronomically important traits in a diverse B. napus breeding population, SKBnNAM, introduced here for the first time. The experiment comprised 50 spring-type B. napus lines, grown and phenotyped in six replicates under two treatment conditions (control and drought) over 38 days in a LemnaTec Scanalyzer 3D facility. Growth traits including plant height, width, projected leaf area, and estimated biovolume were extracted and derived through processing of RGB and NIR images. Anthesis was automatically and accurately scored (97% accuracy) and the number of flowers per plant and day was approximated alongside relevant canopy traits (width, angle). Further, supervised machine learning was used to predict the total number of raceme branches from flower attributes with 91% accuracy (linear regression and Huber regression algorithms) and to identify mild drought stress, a complex trait which typically has to be empirically scored (0.85 area under the receiver operating characteristic curve, random forest classifier algorithm). The study demonstrates the potential of HTP, image processing and computer vision for effective characterization of agronomic trait diversity in B. napus, although limitations of the platform did create significant variation that limited the utility of the data. However, the results underscore the value of machine learning for phenotyping studies, particularly for complex traits such as drought stress resistance.

Natural genetic variation as a tool for discovery in Caenorhabditis nematodes

Over the last 20 years, studies of Caenorhabditis elegans natural diversity have demonstrated the power of quantitative genetic approaches to reveal the evolutionary, ecological, and genetic factors that shape traits. These studies complement the use of the laboratory-adapted strain N2 and enable additional discoveries not possible using only one genetic background. In this chapter, we describe how to perform quantitative genetic studies in Caenorhabditis, with an emphasis on C. elegans. These approaches use correlations between genotype and phenotype across populations of genetically diverse individuals to discover the genetic causes of phenotypic variation. We present methods that use linkage, near-isogenic lines, association, and bulk-segregant mapping, and we describe the advantages and disadvantages of each approach. The power of C. elegans quantitative genetic mapping is best shown in the ability to connect phenotypic differences to specific genes and variants. We will present methods to narrow genomic regions to candidate genes and then tests to identify the gene or variant involved in a quantitative trait. The same features that make C. elegans a preeminent experimental model animal contribute to its exceptional value as a tool to understand natural phenotypic variation.

Sorghum genetic, genomic, and breeding resources

Sorghum research has entered an exciting and fruitful era due to the genetic, genomic, and breeding resources that are now available to researchers and plant breeders. As the world faces the challenges of a rising population and a changing global climate, new agricultural solutions will need to be developed to address the food and fiber needs of the future. To that end, sorghum will be an invaluable crop species as it is a stress-resistant C₄ plant that is well adapted for semi-arid and arid regions. Sorghum has already remained as a staple food crop in many parts of Africa and Asia and is critically important for animal feed and niche culinary applications in other regions, such as the United States. In addition, sorghum has begun to be developed into a promising feedstock for forage and bioenergy production. Due to this increasing demand for sorghum and its potential to address these needs, the continuous development of powerful community resources is required. These resources include vast collections of sorghum germplasm, high-quality reference genome sequences, sorghum association panels for genome-wide association studies of traits involved in food and bioenergy production, mutant populations for rapid discovery of causative genes for phenotypes relevant to sorghum improvement, gene expression atlas, and online databases that integrate all resources and provide the sorghum community with tools that can be used in breeding and genomic studies. Used in tandem, these valuable resources will ensure that the rate, quality, and collaborative potential of ongoing sorghum improvement efforts is able to rival that of other major crops.

Intergenic spaces: a new frontier to improving plant health

To more sustainably mitigate the impact of crop diseases on plant health and productivity, there is a need for broader spectrum, long-lasting resistance traits. Defense response (DR) genes, located throughout the genome, participate in cellular and system-wide defense mechanisms to stave off infection by diverse pathogens. This multigenic resistance avoids rapid evolution of a pathogen to overcome host resistance. DR genes reside within resistance-associated quantitative trait loci (QTL), and alleles of DR genes in resistant varieties are more active during pathogen attack relative to susceptible haplotypes. Differential expression of DR genes results from polymorphisms in their regulatory regions, that includes cis-regulatory elements such as transcription factor binding sites as well as features that influence epigenetic structural changes to modulate chromatin accessibility during infection. Many of these elements are found in clusters, known as cis-regulatory modules (CRMs), which are distributed throughout the host genome. Regulatory regions involved in plant-pathogen interactions may also contain pathogen effector binding elements that regulate DR gene expression, and that, when mutated, result in a change in the plants' response. We posit that CRMs and the multiple regulatory elements that comprise them are potential targets for marker-assisted breeding for broad-spectrum, durable disease resistance.

Haplotype associated RNA expression (HARE) improves prediction of complex traits in maize

Genomic prediction typically relies on associations between single-site polymorphisms and traits of interest. This representation of genomic variability has been successful for predicting many complex traits. However, it usually cannot capture the combination of alleles in haplotypes and it has generated little insight about the biological function of polymorphisms. Here we present a novel and cost-effective method for imputing cis haplotype associated RNA expression (HARE), studied their transferability across tissues, and evaluated genomic prediction models within and across populations. HARE focuses on tightly linked cis acting causal variants in the immediate vicinity of the gene, while excluding trans effects from diffusion and metabolism. Therefore, HARE estimates were more transferrable across different tissues and populations compared to measured transcript expression. We also showed that HARE estimates captured one-third of the variation in gene expression. HARE estimates were used in genomic prediction models evaluated within and across two diverse maize panels–a diverse association panel (Goodman Association panel) and a large half-sib panel (Nested Association Mapping panel)–for predicting 26 complex traits. HARE resulted in up to 15% higher prediction accuracy than control approaches that preserved haplotype structure, suggesting that HARE carried functional information in addition to information about haplotype structure. The largest increase was observed when the model was trained in the Nested Association Mapping panel and tested in the Goodman Association panel. Additionally, HARE yielded higher within-population prediction accuracy as compared to measured expression values. The accuracy achieved by measured expression was variable across tissues, whereas accuracy by HARE was more stable across tissues. Therefore, imputing RNA expression of genes by haplotype is stable, cost-effective, and transferable across populations.

Genomic marker data is widely used in the prediction of many traits. However, prediction has been primarily carried out within populations and without explicit modeling of RNA or protein expression. In this study, we explored the prediction of field traits within and across populations using estimated RNA expression attributable to only the DNA sequence around a gene. We showed that the estimated RNA expression was more transferable across populations and tissues than measured RNA expression. We improved prediction of field traits up to 15% using estimated gene expression as compared to observed expression or gene sequence alone. Overall, these findings indicate that structural and functional information in the gene sequence is highly transferable.

The identification of grain size genes by RapMap reveals directional selection during rice domestication

Cloning quantitative trait locus (QTL) is time consuming and laborious, which hinders the understanding of natural variation and genetic diversity. Here, we introduce RapMap, a method for rapid multi-QTL mapping by employing F2 gradient populations (F2GPs) constructed by minor-phenotypic-difference accessions. The co-segregation standard of the single-locus genetic models ensures simultaneous integration of a three-in-one framework in RapMap i.e. detecting a real QTL, confirming its effect, and obtaining its near-isogenic line-like line (NIL-LL). We demonstrate the feasibility of RapMap by cloning eight rice grain-size genes using 15 F2GPs in three years. These genes explain a total of 75% of grain shape variation. Allele frequency analysis of these genes using a large germplasm collection reveals directional selection of the slender and long grains in indica rice domestication. In addition, major grain-size genes have been strongly selected during rice domestication. We think application of RapMap in crops will accelerate gene discovery and genomic breeding.

Development of an Aus-Derived Nested Association Mapping (Aus-NAM) Population in Rice

A genetic resource for studying genetic architecture of agronomic traits and environmental adaptation is essential for crop improvements. Here, we report the development of a rice nested association mapping population (aus-NAM) using 7 aus varieties as diversity donors and T65 as the common parent. Aus-NAM showed broad phenotypic variations. To test whether aus-NAM was useful for quantitative trait loci (QTL) mapping, known flowering genes (Ehd1, Hd1, and Ghd7) in rice were characterized using single-family QTL mapping, joint QTL mapping, and the methods based on genome-wide association study (GWAS). Ehd1 was detected in all the seven families and all the methods. On the other hand, Hd1 and Ghd7 were detected in some families, and joint QTL mapping and GWAS-based methods resulted in weaker and uncertain peaks. Overall, the high allelic variations in aus-NAM provide a valuable genetic resource for the rice community.

Meta-Analysis Identifies Pleiotropic Loci Controlling Phenotypic Trade-offs in Sorghum

Community association populations are composed of phenotypically and genetically diverse accessions. Once these populations are genotyped, the resulting marker data can be reused by different groups investigating the genetic basis of different traits. Because the same genotypes are observed and scored for a wide range of traits in different environments, these populations represent a unique resource to investigate pleiotropy. Here we assembled a set of 234 separate trait datasets for the Sorghum Association Panel, a group of 406 sorghum genotypes widely employed by the sorghum genetics community. Comparison of genome wide association studies conducted with two independently generated marker sets for this population demonstrate that existing genetic marker sets do not saturate the genome and likely capture only 35-43% of potentially detectable loci controlling variation for traits scored in this population. While limited evidence for pleiotropy was apparent in cross-GWAS comparisons, a multivariate adaptive shrinkage approach recovered both known pleiotropic effects of existing loci and new pleiotropic effects, particularly significant impacts of known dwarfing genes on root architecture. In addition, we identified new loci with pleiotropic effects consistent with known trade-offs in sorghum development. These results demonstrate the potential for mining existing trait datasets from widely used community association populations to enable new discoveries from existing trait datasets as new, denser genetic marker datasets are generated for existing community association populations.

The carbon isotopic signature of C4 crops and its applicability in breeding for climate resilience

KEY MESSAGE

Carbon isotope discrimination is a promising trait for indirect screening for improved water use efficiency of C4 crops. In the context of a changing climate, drought is one of the major factors limiting plant growth and yield. Hence, breeding efforts are directed toward improving water use efficiency (WUE) as a key factor in climate resilience and sustainability of crop production. As WUE is a complex trait and its evaluation is rather resource consuming, proxy traits, which are easier to screen and reliably reflect variation in WUE, are needed. In C3 crops, a trait established to be indicative for WUE is the carbon isotopic composition (δ13C) of plant material, which reflects the preferential assimilation of the lighter carbon isotope 12C over 13C during photosynthesis. In C4 crops, carbon fixation is more complex and δ13C thus depends on many more factors than in C3 crops. Recent physiological and genetic studies indicate a correlation between δ13C and WUE also in C4 crops, as well as a colocalization of quantitative trait loci for the two traits. Moreover, significant intraspecific variation as well as a medium to high heritability of δ13C has been shown in some of the main C4 crops, such as maize, sorghum and sugarcane, indicating its potential for indirect selection and breeding. Further research on physiological, genetic and environmental components influencing δ13C is needed to support its application in improving WUE and making C4 crops resilient to climate change.

DArTseq-Based High-Throughput SilicoDArT and SNP Markers Applied for Association Mapping of Genes Related to Maize Morphology

Today, agricultural productivity is essential to meet the needs of a growing population, and is also a key tool in coping with climate change. Innovative plant breeding technologies such as molecular markers, phenotyping, genotyping, the CRISPR/Cas method and next-generation sequencing can help agriculture meet the challenges of the 21st century more effectively. Therefore, the aim of the research was to identify single-nucleotide polymorphisms (SNPs) and SilicoDArT markers related to select morphological features determining the yield in maize. The plant material consisted of ninety-four inbred lines of maize of various origins. These lines were phenotyped under field conditions. A total of 14 morphological features was analyzed. The DArTseq method was chosen for genotyping because this technique reduces the complexity of the genome by restriction enzyme digestion. Subsequently, short fragment sequencing was used. The choice of a combination of restrictases allowed the isolation of highly informative low copy fragments of the genome. Thanks to this method, 90% of the obtained DArTseq markers are complementary to the unique sequences of the genome. All the observed features were normally distributed. Analysis of variance indicated that the main effect of lines was statistically significant (p < 0.001) for all 14 traits of study. Thanks to the DArTseq analysis with the use of next-generation sequencing (NGS) in the studied plant material, it was possible to identify 49,911 polymorphisms, of which 33,452 are SilicoDArT markers and the remaining 16,459 are SNP markers. Among those mentioned, two markers associated with four analyzed traits deserved special attention: SNP (4578734) and SilicoDArT (4778900). SNP marker 4578734 was associated with the following features: anthocyanin coloration of cob glumes, number of days from sowing to anthesis, number of days from sowing to silk emergence and anthocyanin coloration of internodes. SilicoDArT marker 4778900 was associated with the following features: number of days from sowing to anthesis, number of days from sowing to silk emergence, tassel: angle between the axis and lateral branches and plant height. Sequences with a length of 71 bp were used for physical mapping. The BLAST and EnsemblPlants databases were searched against the maize genome to identify the positions of both markers. Marker 4578734 was localized on chromosome 7, the closest gene was Zm00001d022467, approximately 55 Kb apart, encoding anthocyanidin 3-O-glucosyltransferase. Marker 4778900 was located on chromosome 7, at a distance of 45 Kb from the gene Zm00001d045261 encoding starch synthase I. The latter observation indicated that these flanking SilicoDArT and SNP markers were not in a state of linkage disequilibrium.

Diversifying maize genomic selection models

Genomic selection (GS) is one of the most powerful tools available for maize breeding. Its use of genome-wide marker data to estimate breeding values translates to increased genetic gains with fewer breeding cycles. In this review, we cover the history of GS and highlight particular milestones during its adaptation to maize breeding. We discuss how GS can be applied to developing superior maize inbreds and hybrids. Additionally, we characterize refinements in GS models that could enable the encapsulation of non-additive genetic effects, genotype by environment interactions, and multiple levels of the biological hierarchy, all of which could ultimately result in more accurate predictions of breeding values. Finally, we suggest the stages in a maize breeding program where it would be beneficial to apply GS. Given the current sophistication of high-throughput phenotypic, genotypic, and other -omic level data currently available to the maize community, now is the time to explore the implications of their incorporation into GS models and thus ensure that genetic gains are being achieved as quickly and efficiently as possible.

Development of an Australian Bread Wheat Nested Association Mapping Population, a New Genetic Diversity Resource for Breeding under Dry and Hot Climates

Genetic diversity, knowledge of the genetic architecture of the traits of interest and efficient means of transferring the desired genetic diversity into the relevant genetic background are prerequisites for plant breeding. Exotic germplasm is a rich source of genetic diversity; however, they harbor undesirable traits that limit their suitability for modern agriculture. Nested association mapping (NAM) populations are valuable genetic resources that enable incorporation of genetic diversity, dissection of complex traits and providing germplasm to breeding programs. We developed the OzNAM by crossing and backcrossing 73 diverse exotic parents to two Australian elite varieties Gladius and Scout. The NAM parents were genotyped using the iSelect wheat 90K Infinium SNP array, and the progeny were genotyped using a custom targeted genotyping-by-sequencing assay based on molecular inversion probes designed to target 12,179 SNPs chosen from the iSelect wheat 90K Infinium SNP array of the parents. In total, 3535 BC1F4:6 RILs from 125 families with 21 to 76 lines per family were genotyped and we found 4964 polymorphic and multi-allelic haplotype markers that spanned the whole genome. A subset of 530 lines from 28 families were evaluated in multi-environment trials over three years. To demonstrate the utility of the population in QTL mapping, we chose to map QTL for maturity and plant height using the RTM-GWAS approach and we identified novel and known QTL for maturity and plant height.