Genetic design and statistical power of nested association mapping in maize

Agile Genetics: Single gene resolution without the fuss

Gene discovery reveals new biology, expands the utility of marker-assisted selection, and enables targeted mutagenesis. Still, such discoveries can take over a decade. We present a general strategy, "Agile Genetics," that uses nested, structured populations to overcome common limits on gene resolution. Extensive simulation work on realistic genetic architectures shows that, at population sizes of >5000 samples, single gene-resolution can be achieved using bulk segregant pools. At this scale, read depth and technical replication become major drivers of resolution. Emerging enrichment methods to address coverage are on the horizon; we describe one possibility - iterative depth sequencing (ID-seq). In addition, graph-based pangenomics in experimental populations will continue to maximize accuracy and improve interpretation. Based on this merger of agronomic scale with molecular and bioinformatic innovation, we predict a new age of rapid gene discovery.

Linkage and association mapping in multi-parental populations reveal the genetic basis of carotenoid variation in maize kernels

Carotenoids are indispensable to plants and critical components of the human diet. The carotenoid metabolic pathway is conserved across plant species, but our understanding of the genetic basis of carotenoid variation remains limited for the seeds of most cereal crops. To address this issue, we systematically performed linkage and association mapping for eight carotenoid traits using six recombinant inbred line (RIL) populations. Single linkage mapping (SLM) and joint linkage mapping (JLM) identified 77 unique additive QTLs and 104 pairs of epistatic QTLs. Among these QTLs, we identified 22 overlapping hotspots of additive and epistatic loci, highlighting the important contributions of some QTLs to carotenoid levels through additive or epistatic mechanisms. A genome-wide association study based on all RILs detected 244 candidate genes significantly associated with carotenoid traits, 23 of which were annotated as carotenoid pathway genes. Effect comparisons suggested that a small number of loci linked to pathway genes have substantial effects on carotenoid variation in our tested populations, but many loci not associated with pathway genes also make important contributions to carotenoid variation. We identified ZmPTOX as the causal gene for a QTL hotspot (Q10/JLM10/GWAS019); this gene encodes a putative plastid terminal oxidase that produces plastoquinone-9 used by two enzymes in the carotenoid pathway. Natural variants in the promoter and second exon of ZmPTOX were found to alter carotenoid levels. This comprehensive assessment of the genetic mechanisms underlying carotenoid variation establishes a foundation for rewiring carotenoid metabolism and accumulation for efficient carotenoid biofortification.

Harnessing landrace diversity empowers wheat breeding

Harnessing genetic diversity in major staple crops through the development of new breeding capabilities is essential to ensure food security1. Here we examined the genetic and phenotypic diversity of the A. E. Watkins landrace collection2 of bread wheat (Triticum aestivum), a major global cereal, by whole-genome re-sequencing of 827 Watkins landraces and 208 modern cultivars and in-depth field evaluation spanning a decade. We found that modern cultivars are derived from two of the seven ancestral groups of wheat and maintain very long-range haplotype integrity. The remaining five groups represent untapped genetic sources, providing access to landrace-specific alleles and haplotypes for breeding. Linkage disequilibrium-based haplotypes and association genetics analyses link Watkins genomes to the thousands of identified high-resolution quantitative trait loci and significant marker-trait associations. Using these structured germplasm, genotyping and informatics resources, we revealed many Watkins-unique beneficial haplotypes that can confer superior traits in modern wheat. Furthermore, we assessed the phenotypic effects of 44,338 Watkins-unique haplotypes, introgressed from 143 prioritized quantitative trait loci in the context of modern cultivars, bridging the gap between landrace diversity and current breeding. This study establishes a framework for systematically utilizing genetic diversity in crop improvement to achieve sustainable food security.

Genotype-Specific Activation of Autophagy during Heat Wave in Wheat

Recycling of unnecessary or dysfunctional cellular structures through autophagy plays a critical role in cellular homeostasis and environmental resilience. Therefore, the autophagy trait may have been unintentionally selected in wheat breeding programs for higher yields in arid climates. This hypothesis was tested by measuring the response of three common autophagy markers, ATG7, ATG8, and NBR1, to a heat wave under reduced soil moisture content in 16 genetically diverse spring wheat landraces originating from different geographical locations. We observed in the greenhouse trials that ATG8 and NBR1 exhibited genotype-specific responses to a 1 h, 40 °C heat wave, while ATG7 did not show a consistent response. Three genotypes from Uruguay, Mozambique, and Afghanistan showed a pattern consistent with higher autophagic activity: decreased or stable abundance of both ATG8 and NBR1 proteins, coupled with increased transcription of ATG8 and NBR1. In contrast, three genotypes from Pakistan, Ethiopia, and Egypt exhibited elevated ATG8 protein levels alongside reduced or unaltered ATG8 transcript levels, indicating a potential suppression or no change in autophagic activity. Principal component analysis demonstrated a correlation between lower abundance of ATG8 and NBR1 proteins and higher yield in the field trials. We found that (i) the combination of heat and drought activated autophagy only in several genotypes, suggesting that despite being a resilience mechanism, autophagy is a heat-sensitive process; (ii) higher autophagic activity correlates positively with greater yield; (iii) the lack of autophagic activity in some high-yielding genotypes suggests contribution of alternative stress-resilient mechanisms; and (iv) enhanced autophagic activity in response to heat and drought was independently selected by wheat breeding programs in different geographic locations.

Genetic architecture of heritable leaf microbes

Host-associated microbiomes are shaped by both their environment and host genetics, and often impact host performance. The scale of host genetic variation important to microbes is largely unknown yet fundamental to the community assembly of host-associated microbiomes, with implications for the eco-evolutionary dynamics of microbes and hosts. Using Ipomoea hederacea, ivyleaf morning glory, we generated matrilines differing in quantitative genetic variation and leaf shape, which is controlled by a single Mendelian locus. We then investigated the relative roles of Mendelian and quantitative genetic variation in structuring the leaf microbiome and how these two sources of genetic variation contributed to microbe heritability. We found that despite large effects of the environment, both Mendelian and quantitative genetic host variation contribute to microbe heritability and that the cumulative small effect genomic differences due to matriline explained as much or more microbial variation than a single large effect Mendelian locus. Furthermore, our results are the first to suggest that leaf shape itself contributes to variation in the abundances of some phyllosphere microbes.IMPORTANCEWe investigated how host genetic variation affects the assembly of Ipomoea hederacea's natural microbiome. We found that the genetic architecture of leaf-associated microbiomes involves both quantitative genetic variation and Mendelian traits, with similar contributions to microbe heritability. The existence of Mendelian and quantitative genetic variation for host-associated microbes means that plant evolution at the leaf shape locus or other quantitative genetic loci has the potential to shape microbial abundance and community composition.

Beyond pathways: Accelerated flavonoids candidate identification and novel exploration of enzymatic properties using combined mapping populations of wheat

Although forward-genetics-metabolomics methods such as mGWAS and mQTL have proven effective in providing myriad loci affecting metabolite contents, they are somehow constrained by their respective constitutional flaws such as the hidden population structure for GWAS and insufficient recombinant rate for QTL. Here, the combination of mGWAS and mQTL was performed, conveying an improved statistical power to investigate the flavonoid pathways in common wheat. A total of 941 and 289 loci were, respectively, generated from mGWAS and mQTL, within which 13 of them were co-mapped using both approaches. Subsequently, the mGWAS or mQTL outputs alone and their combination were, respectively, utilized to delineate the metabolic routes. Using this approach, we identified two MYB transcription factor encoding genes and five structural genes, and the flavonoid pathway in wheat was accordingly updated. Moreover, we have discovered some rare-activity-exhibiting flavonoid glycosyl- and methyl-transferases, which may possess unique biological significance, and harnessing these novel catalytic capabilities provides potentially new breeding directions. Collectively, we propose our survey illustrates that the forward-genetics-metabolomics approaches including multiple populations with high density markers could be more frequently applied for delineating metabolic pathways in common wheat, which will ultimately contribute to metabolomics-assisted wheat crop improvement.

A Quantitative Trait Locus on Maize Chromosome 5 Is Associated with Root-Knot Nematode Resistance

This study provides the first report of a quantitative trait locus (QTL) in maize (Zea mays) for resistance to the southern root-knot nematode (SRKN) (Meloidogyne incognita). The SRKN can feed on the roots of maize in the U.S. Southern Coastal Plain region and can cause yield losses of 30% or more in heavily infested fields. Increases in SRKN density in the soil may reduce the yield for subsequently planted susceptible crops. The use of maize hybrids with resistance to SRKN could prevent an increase in SRKN density, yet no genetic regions have been identified that confer host resistance. In this study, a B73 (susceptible) × Ky21 (resistant) S5 recombinant inbred line (RIL) population was phenotyped for total number of eggs (TE) and root weight. This population had been genotyped using single-nucleotide polymorphisms (SNPs). By utilizing the SNP data with the phenotype data, a single QTL was identified on chromosome 5 that explained 15% of the phenotypic variation (PV) for the number of eggs and 11% of the PV for the number of eggs per gram of root (EGR). Plants that were homozygous for the Ky21 allele for the most associated marker PZA03172.3 had fewer eggs and fewer EGR than the plants that were homozygous or heterozygous for the B73 allele. Thus, the first QTL for SRKN resistance in maize has been identified and could be incorporated into maize hybrids.

Comparative genomics analysis and transposon mutagenesis provides new insights into high menaquinone-7 biosynthetic potential of Bacillus subtilis natto

This research combined Whole-Genome sequencing, intraspecific comparative genomics and transposon mutagenesis to investigate the menaquinone-7 (MK-7) synthesis potential in Bacillus subtilis natto. First, Whole-Genome sequencing showed that Bacillus subtilis natto BN-P15-11-1 contains one single circular chromosome in size of 3,982,436 bp with a GC content of 43.85 %, harboring 4,053 predicted coding genes. Next, the comparative genomics analysis among strain BN-P15-11-1 with model Bacillus subtilis 168 and four typical Bacillus subtilis natto strains proves that the closer evolutionary relationship Bacillus subtilis natto BN-P15-11-1 and Bacillus subtilis 168 both exhibit strong biosynthetic potential. To further dig for MK-7 biosynthesis latent capacity of BN-P15-11-1, we constructed a mutant library using transposons and a high throughput screening method using microplates. We obtained a YqgQ deficient high MK-7 yield strain F4 with a yield 3.02 times that of the parent strain. Experiments also showed that the high yield mutants had defects in different transcription and translation regulatory factor genes, indicating that regulatory factor defects may affect the biosynthesis and accumulation of MK-7 by altering the overall metabolic level. The findings of this study will provide more novel insights on the precise identification and rational utilization of the Bacillus subtilis subspecies for biosynthesis latent capacity.

Dynamic metabolite QTL analyses provide novel biochemical insights into kernel development and nutritional quality improvement in common wheat

Despite recent advances in crop metabolomics, the genetic control and molecular basis of the wheat kernel metabolome at different developmental stages remain largely unknown. Here, we performed widely targeted metabolite profiling of kernels from three developmental stages (grain-filling kernels [FKs], mature kernels [MKs], and germinating kernels [GKs]) using a population of 159 recombinant inbred lines. We detected 625 annotated metabolites and mapped 3173, 3143, and 2644 metabolite quantitative trait loci (mQTLs) in FKs, MKs, and GKs, respectively. Only 52 mQTLs were mapped at all three stages, indicating the high stage specificity of the wheat kernel metabolome. Four candidate genes were functionally validated by in vitro enzymatic reactions and/or transgenic approaches in wheat, three of which mediated the tricin metabolic pathway. Metabolite flux efficiencies within the tricin pathway were evaluated, and superior candidate haplotypes were identified, comprehensively delineating the tricin metabolism pathway in wheat. Finally, additional wheat metabolic pathways were re-constructed by updating them to incorporate the 177 candidate genes identified in this study. Our work provides new information on variations in the wheat kernel metabolome and important molecular resources for improvement of wheat nutritional quality.

Enhanced pan-genomic resources at the maize genetics and genomics database

Pan-genomes, encompassing the entirety of genetic sequences found in a collection of genomes within a clade, are more useful than single reference genomes for studying species diversity. This is especially true for a species like Zea mays, which has a particularly diverse and complex genome. Presenting pan-genome data, analyses, and visualization is challenging, especially for a diverse species, but more so when pan-genomic data is linked to extensive gene model and gene data, including classical gene information, markers, insertions, expression and proteomic data, and protein structures as is the case at MaizeGDB. Here, we describe MaizeGDB's expansion to include the genic subset of the Zea pan-genome in a pan-gene data center featuring the maize genomes hosted at MaizeGDB, and the outgroup teosinte Zea genomes from the Pan-Andropoganeae project. The new data center offers a variety of browsing and visualization tools, including sequence alignment visualization, gene trees and other tools, to explore pan-genes in Zea that were calculated by the pipeline Pandagma. Combined, these data will help maize researchers study the complexity and diversity of Zea, and to use the comparative functions to validate pan-gene relationships for a selected gene model.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

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.

A new winter wheat genetic resource harbors untapped diversity from synthetic hexaploid wheat

KEY MESSAGE

The NIAB_WW_SHW_NAM population, a large nested association mapping panel, is a useful resource for mapping QTL from synthetic hexaploid wheat that can improve modern elite wheat cultivars. The allelic richness harbored in progenitors of hexaploid bread wheat (Triticum aestivum L.) is a useful resource for addressing the genetic diversity bottleneck in modern cultivars. Synthetic hexaploid wheat (SHW) is created through resynthesis of the hybridisation events between the tetraploid (Triticum turgidum subsp. durum Desf.) and diploid (Aegilops tauschii Coss.) bread wheat progenitors. We developed a large and diverse winter wheat nested association mapping (NAM) population (termed the NIAB_WW_SHW_NAM) consisting of 3241 genotypes derived from 54 nested back-cross 1 (BC1) populations, each formed via back-crossing a different primary SHW into the UK winter wheat cultivar 'Robigus'. The primary SHW lines were created using 15 T. durum donors and 47 Ae. tauschii accessions that spanned the lineages and geographical range of the species. Primary SHW parents were typically earlier flowering, taller and showed better resistance to yellow rust infection (Yr) than 'Robigus'. The NIAB_WW_SHW_NAM population was genotyped using a single nucleotide polymorphism (SNP) array and 27 quantitative trait loci (QTLs) were detected for flowering time, plant height and Yr resistance. Across multiple field trials, a QTL for Yr resistance was found on chromosome 4D that corresponded to the Yr28 resistance gene previously reported in other SHW lines. These results demonstrate the value of the NIAB_WW_SHW_NAM population for genetic mapping and provide the first evidence of Yr28 working in current UK environments and genetic backgrounds. These examples, coupled with the evidence of commercial wheat breeders selecting promising genotypes, highlight the potential value of the NIAB_WW_SHW_NAM to variety improvement.

Genetic basis of maize stalk strength decoded via linkage and association mapping

Stalk lodging is a severe problem that limits maize production worldwide, although little attention has been given to its genetic basis. Here we measured rind penetrometer resistance (RPR), an effective index for stalk lodging, in a multi-parent population of 1948 recombinant inbred lines (RILs) and an association population of 508 inbred lines (AMP508). Linkage and association mapping identified 53 and 29 single quantitative trait loci (QTLs) and 50 and 19 pairs of epistatic interactions for RPR in the multi-parent population and AMP508 population, respectively. Phenotypic variation explained by all identified epistatic QTLs (up to ~5%) was much less than that explained by all single additive QTLs (up to ~33% in the multi-parent population and ~ 60% in the AMP508 population). Among all detected QTLs, only eight single QTLs explained >10% of phenotypic variation in single RIL populations. Alleles that increased RPR were enriched in tropical/subtropical (TST) groups from the AMP508 population. Based on genome-wide association studies in both populations, we identified 137 candidate genes affecting RPR, which were assigned to multiple biological processes, such as the biosynthesis of cell wall components. Sixty-six candidate genes were cross-validated by multiple methods or populations. Most importantly, 23 candidate genes were upregulated or downregulated in high-RPR lines relative to low-RPR lines, supporting the associations between candidate genes and RPR. These findings reveal the complex nature of the genetic basis underlying RPR and provide loci or candidate genes for developing elite varieties that are resistant to stalk lodging via molecular breeding.

Whole genome resequencing and phenotyping of MAGIC population for high resolution mapping of drought tolerance in chickpea

Terminal drought is one of the major constraints to crop production in chickpea (Cicer arietinum L.). In order to map drought tolerance related traits at high resolution, we sequenced multi-parent advanced generation intercross (MAGIC) population using whole genome resequencing approach and phenotyped it under drought stress environments for two consecutive years (2013-14 and 2014-15). A total of 52.02 billion clean reads containing 4.67 TB clean data were generated on the 1136 MAGIC lines and eight parental lines. Alignment of clean data on to the reference genome enabled identification of a total, 932,172 of SNPs, 35,973 insertions, and 35,726 deletions among the parental lines. A high-density genetic map was constructed using 57,180 SNPs spanning a map distance of 1606.69 cM. Using compressed mixed linear model, genome-wide association study (GWAS) enabled us to identify 737 markers significantly associated with days to 50% flowering, days to maturity, plant height, 100 seed weight, biomass, and harvest index. In addition to the GWAS approach, an identity-by-descent (IBD)-based mixed model approach was used to map quantitative trait loci (QTLs). The IBD-based mixed model approach detected major QTLs that were comparable to those from the GWAS analysis as well as some exclusive QTLs with smaller effects. The candidate genes like FRIGIDA and CaTIFY4b can be used for enhancing drought tolerance in chickpea. The genomic resources, genetic map, marker-trait associations, and QTLs identified in the study are valuable resources for the chickpea community for developing climate resilient chickpeas.

Revitalizing agriculture: next-generation genotyping and -omics technologies enabling molecular prediction of resilient traits in the Solanaceae family

This review highlights -omics research in Solanaceae family, with a particular focus on resilient traits. Extensive research has enriched our understanding of Solanaceae genomics and genetics, with historical varietal development mainly focusing on disease resistance and cultivar improvement but shifting the emphasis towards unveiling resilience mechanisms in genebank-preserved germplasm is nowadays crucial. Collecting such information, might help researchers and breeders developing new experimental design, providing an overview of the state of the art of the most advanced approaches for the identification of the genetic elements laying behind resilience. Building this starting point, we aim at providing a useful tool for tackling the global agricultural resilience goals in these crops.

Postzygotic barriers persist despite ongoing introgression in hybridizing Mimulus species

The evolution of postzygotic isolation is thought to be a key step in maintaining species boundaries upon secondary contact, yet the dynamics and persistence of hybrid incompatibilities in naturally hybridizing species are not well understood. Here, we explore these issues using genetic mapping in three independent populations of recombinant inbred lines between naturally hybridizing monkeyflowers, Mimulus guttatus and Mimulus nasutus, from the sympatric Catherine Creek population. We discover that the three M. guttatus founders differ dramatically in admixture history, with nearly a quarter of one founder's genome introgressed from M. nasutus. Comparative genetic mapping in the three RIL populations reveals three new putative inversions, each one segregating among the M. guttatus founders, two due to admixture. We find strong, genome-wide transmission ratio distortion in all RILs, but patterns are highly variable among the three populations. At least some of this distortion appears to be explained by epistatic selection favouring parental genotypes, but tests of inter-chromosomal linkage disequilibrium also reveal multiple candidate Dobzhansky-Muller incompatibilities. We also map several genetic loci for hybrid pollen viability, including two interacting pairs that coincide with peaks of distortion. Remarkably, even with this limited sample of three M. guttatus lines, we discover abundant segregating variation for hybrid incompatibilities with M. nasutus, suggesting this population harbours diverse contributors to postzygotic isolation. Moreover, even with substantial admixture, hybrid incompatibilities between Mimulus species persist, suggesting postzygotic isolation might be a potent force in maintaining species barriers in this system.

Gene function annotations for the maize NAM founder lines

OBJECTIVES

We annotated the latest published sequences of the 26 Zea mays Nested Association Mapping (NAM) founder lines using GOMAP, the Gene Ontology Meta Annotator for Plants. The maize NAM panel enables researchers to understand and identify the genetic basis of complex traits. Annotations of predicted functions for genes can help researchers investigate gene-phenotype associations, prioritize candidate genes for phenotypes of interest, and formulate testable hypotheses about gene function/phenotype associations. The creation and release of high-confidence, high-coverage gene function annotation sets for the NAM founder lines is critical to accelerate the generation of knowledge in maize genetics research. GOMAP is a high-throughput computational pipeline that annotates gene functions genome-wide in plant genomes using Gene Ontology functional class terms. Here we report and share GOMAP-generated functional annotations for the NAM founder lines.

DATA DESCRIPTION

Datasets include the protein sequences used as input, GOMAP-generated annotation files, scripts used to update obsolete terms, and GAF-formatted tab-delimited text files of gene function annotations along with README files that describe formatting, content, and how files relate to each other.

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.

Recombinant inbred lines derived from wide crosses in Pisum

Genomic resources are becoming available for Pisum but to link these to phenotypic diversity requires well marked populations segregating for relevant traits. Here we describe two such resources. Two recombinant inbred populations, derived from wide crosses in Pisum are described. One high resolution mapping population involves cv Caméor, for which the first pea whole genome assembly was obtained, crossed to JI0281, a basally divergent P. sativum sativum landrace from Ethiopia. The other is an inter sub-specific cross between P. s. sativum and the independently domesticated P. s. abyssinicum. The corresponding genetic maps provide information on chromosome level sequence assemblies and identify structural differences between the genomes of these two Pisum subspecies. In order to visualise chromosomal translocations that distinguish the mapping parents, we created a simplified version of Threadmapper to optimise it for interactive 3-dimensional display of multiple linkage groups. The genetic mapping of traits affecting seed coat roughness and colour, plant height, axil ring pigmentation, leaflet number and leaflet indentation enabled the definition of their corresponding genomic regions. The consequence of structural rearrangement for trait analysis is illustrated by leaf serration. These analyses pave the way for identification of the underlying genes and illustrate the utility of these publicly available resources. Segregating inbred populations derived from wide crosses in Pisum, together with the associated marker data, are made publicly available for trait dissection. Genetic analysis of these populations is informative about chromosome scale assemblies, structural diversity in the pea genome and has been useful for the fine mapping of several discrete and quantitative traits.

A genomic analysis of the University of Nebraska Replicated Recurrent Selection program

The inbred-hybrid system of maize breeding closely resembles a reciprocal full-sib (RFS) selection program. Studying changes in genetic variation as a result of RFS selection can help illuminate long-standing questions regarding the relative roles of selection and genetic drift and help understand the nature of adaptation occurring in selection programs. The University of Nebraska-Lincoln Replicated Recurrent Selection (UNL-RpRS) program underwent eight cycles of replicated RFS and S1-progeny selection, making it a powerful system to study genomic changes accompanying selection for inter-population performance. The objectives of this study were to identify regions of the genome under selection after eight cycles of selection and evaluate the effect eight cycles of selection for inter-population full-sib performance had in expanding genome-wide and localized population structure. We address these questions with a large set of individuals sampled from the UNL-RpRS program with dense genotyping-by-sequence data. We found evidence of parallel selection signatures in the UNL-RpRS program, with a region on chromosome 7 being implicated in three of the four selection systems studied. Regions that displayed selection signatures across independently run selection programs represent regions likely to be capitalizing on standing genetic variation and support a soft sweep model of adaptation. We did not find selection to be a strong force in diverging populations undergoing RFS. This could be due to the nature of adaptation occurring in these populations, underlying gene action, or a result of unstable genetic topographies.

Genetic and transcriptomic dissection of host defense to Goss's bacterial wilt and leaf blight of maize

Goss's wilt, caused by the Gram-positive actinobacterium Clavibacter nebraskensis, is an important bacterial disease of maize. The molecular and genetic mechanisms of resistance to the bacterium, or, in general, Gram-positive bacteria causing plant diseases, remain poorly understood. Here, we examined the genetic basis of Goss's wilt through differential gene expression, standard genome-wide association mapping (GWAS), extreme phenotype (XP) GWAS using highly resistant (R) and highly susceptible (S) lines, and quantitative trait locus (QTL) mapping using 3 bi-parental populations, identifying 11 disease association loci. Three loci were validated using near-isogenic lines or recombinant inbred lines. Our analysis indicates that Goss's wilt resistance is highly complex and major resistance genes are not commonly present. RNA sequencing of samples separately pooled from R and S lines with or without bacterial inoculation was performed, enabling identification of common and differential gene responses in R and S lines. Based on expression, in both R and S lines, the photosynthesis pathway was silenced upon infection, while stress-responsive pathways and phytohormone pathways, namely, abscisic acid, auxin, ethylene, jasmonate, and gibberellin, were markedly activated. In addition, 65 genes showed differential responses (up- or down-regulated) to infection in R and S lines. Combining genetic mapping and transcriptional data, individual candidate genes conferring Goss's wilt resistance were identified. Collectively, aspects of the genetic architecture of Goss's wilt resistance were revealed, providing foundational data for mechanistic studies.

A rapid method to quantify vein density in C4 plants using starch staining

C4 photosynthesis has evolved multiple times in the angiosperms and typically involves alterations to the biochemistry, cell biology and development of leaves. One common modification found in C4 plants compared with the ancestral C3 state is an increase in vein density such that the leaf contains a larger proportion of bundle sheath cells. Recent findings indicate that there may be significant intraspecific variation in traits such as vein density in C4 plants but to use such natural variation for trait-mapping, rapid phenotyping would be required. Here we report a high-throughput method to quantify vein density that leverages the bundle sheath-specific accumulation of starch found in C4 species. Starch staining allowed high-contrast images to be acquired permitting image analysis with MATLAB- and Python-based programmes. The method works for dicotyledons and monocotolydons. We applied this method to Gynandropsis gynandra where significant variation in vein density was detected between natural accessions, and Zea mays where no variation was apparent in the genotypically diverse lines assessed. We anticipate this approach will be useful to map genes controlling vein density in C4 species demonstrating natural variation for this trait.

A key to totipotency: Wuschel-like homeobox 2a unlocks embryogenic culture response in maize (Zea mays L.)

The ability of plant somatic cells to dedifferentiate, form somatic embryos and regenerate whole plants in vitro has been harnessed for both clonal propagation and as a key component of plant genetic engineering systems. Embryogenic culture response is significantly limited, however, by plant genotype in most species. This impedes advancements in both plant transformation-based functional genomics research and crop improvement efforts. We utilized natural variation among maize inbred lines to genetically map somatic embryo generation potential in tissue culture and identify candidate genes underlying totipotency. Using a series of maize lines derived from crosses involving the culturable parent A188 and the non-responsive parent B73, we identified a region on chromosome 3 associated with embryogenic culture response and focused on three candidate genes within the region based on genetic position and expression pattern. Two candidate genes showed no effect when ectopically expressed in B73, but the gene Wox2a was found to induce somatic embryogenesis and embryogenic callus proliferation. Transgenic B73 cells with strong constitutive expression of the B73 and A188 coding sequences of Wox2a were found to produce somatic embryos at similar frequencies, demonstrating that sufficient expression of either allele could rescue the embryogenic culture phenotype. Transgenic B73 plants were regenerated from the somatic embryos without chemical selection and no pleiotropic effects were observed in the Wox2a overexpression lines in the regenerated T0 plants or in the two independent events which produced T1 progeny. In addition to linking natural variation in tissue culture response to Wox2a, our data support the utility of Wox2a in enabling transformation of recalcitrant genotypes.

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.

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.

Combining multiplex gene editing and doubled haploid technology in maize

A major advantage of using CRISPR/Cas9 for gene editing is multiplexing, that is, the simultaneous targeting of many genes. However, primary transformants typically contain hetero-allelic mutations or are genetic mosaic, while genetically stable lines that are homozygous are desired for functional analysis. Currently, a dedicated and labor-intensive effort is required to obtain such higher-order mutants through several generations of genetic crosses and genotyping. We describe the design and validation of a rapid and efficient strategy to produce lines of genetically identical plants carrying various combinations of homozygous edits, suitable for replicated analysis of phenotypical differences. This approach was achieved by combining highly multiplex gene editing in Zea mays (maize) with in vivo haploid induction and efficient in vitro generation of doubled haploid plants using embryo rescue doubling. By combining three CRISPR/Cas9 constructs that target in total 36 genes potentially involved in leaf growth, we generated an array of homozygous lines with various combinations of edits within three generations. Several genotypes show a reproducible 10% increase in leaf size, including a septuple mutant combination. We anticipate that our strategy will facilitate the study of gene families via multiplex CRISPR mutagenesis and the identification of allele combinations to improve quantitative crop 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.

The NLRomes of Zea mays NAM founder lines and Zea luxurians display presence-absence variation, integrated domain diversity, and mobility

Plant pathogens cause significant crop loss worldwide, and new resistance genes deployed to combat diseases can be overcome quickly. Understanding the existing resistance gene diversity within the germplasm of major crops, such as maize, is crucial for the development of new disease-resistant varieties. We analysed the nucleotide-binding leucine-rich repeat receptors (NLRs) of 26 recently sequenced diverse founder lines from the maize nested association mapping (NAM) population and compared them to the R gene complement present in a wild relative of maize, Zea luxurians. We found that NLRs in both species contain a large diversity of atypical integrated domains, including many domains that have not previously been found in the NLRs of other species. Additionally, the single Z. luxurians genome was found to have greater integrated atypical domain diversity than all 26 NAM founder lines combined, indicating that this species may represent a rich source of novel resistance genes. NLRs were also found to have very high sequence diversity and presence-absence variation among the NAM founder lines, with a large NLR cluster on Chr10 representing a diversity hotspot. Additionally, NLRs were shown to be mobile within maize genomes, with several putative interchromosomal translocations identified.

A simulation framework for reciprocal recurrent selection-based hybrid breeding under transparent and opaque simulators

Hybrid breeding is an established and effective process to improve offspring performance, while it is resource-intensive and time-consuming for the recurrent process in reality. To enable breeders and researchers to evaluate the effectiveness of competing decision-making strategies, we present a modular simulation framework for reciprocal recurrent selection-based hybrid breeding. Consisting of multiple modules such as heterotic separation, genomic prediction, and genomic selection, this simulation framework allows breeders to efficiently simulate the hybrid breeding process with multiple options of simulators and decision-making strategies. We also integrate the recently proposed concepts of transparent and opaque simulators into the framework in order to reflect the breeding process more realistically. Simulation results show the performance comparison among different breeding strategies under the two simulators.

The potentialities of omics resources for millet improvement

Millets are nutrient-rich (nutri-rich) cereals with climate resilience attributes. However, its full productive potential is not realized due to the lack of a focused yield improvement approach, as evidenced by the available literature. Also, the lack of well-characterized genomic resources significantly limits millet improvement. But the recent availability of genomic data and advancement in omics tools has shown its enormous potential to enhance the efficiency and precision faced by conventional breeding in millet improvement. The development of high throughput genotyping platforms based on next-generation sequencing (NGS) has provided a low-cost method for genomic information, specifically for neglected nutri-rich cereals with the availability of a limited number of reference genome sequences. NGS has created new avenues for millet biotechnological interventions such as mutation-based study, GWAS, GS, and other omics technologies. The simultaneous discovery of high-throughput markers and multiplexed genotyping platform has aggressively aided marker-assisted breeding for millet improvement. Therefore, omics technology offers excellent opportunities to explore and combine useful variations for targeted traits that could impart high nutritional value to high-yielding cultivars under changing climatic conditions. In millet improvement, an in-depth account of NGS, integrating genomics data with different biotechnology tools, is reviewed in this context.

Identification of Candidate QTLs and Genes for Ear Diameter by Multi-Parent Population in Maize

Ear diameter (ED) is a critical component of grain yield (GY) in maize (Zea mays L.). Studying the genetic basis of ED in maize is of great significance in enhancing maize GY. Against this backdrop, this study was framed to (1) map the ED-related quantitative trait locus (QTL) and SNPs associated with ED; and (2) identify putative functional genes that may affect ED in maize. To accomplish this, an elite maize inbred line, Ye107, which belongs to the Reid heterotic group, was used as a common parent and crossed with seven elite inbred lines from three different heterotic groups (Suwan1, Reid, and nonReid) that exhibited abundant genetic variation in ED. This led to the construction of a multi-parent population consisting of 1215 F₇ recombinant inbred lines (F₇RILs). A genome-wide association study (GWAS) and linkage analysis were then conducted for the multi-parent population using 264,694 high-quality SNPs generated via the genotyping-by-sequencing method. Our study identified a total of 11 SNPs that were significantly associated with ED through the GWAS, and three QTLs were revealed by the linkage analysis for ED. The major QTL on chromosome 1 was co-identified in the region by the GWAS at SNP_143985532. SNP_143985532, located upstream of the Zm00001d030559 gene, encodes a callose synthase that is expressed in various tissues, with the highest expression level in the maize ear primordium. Haplotype analysis indicated that the haplotype B (allele AA) of Zm00001d030559 was positively correlated with ED. The candidate genes and SNPs identified in this study provide crucial insights for future studies on the genetic mechanism of maize ED formation, cloning of ED-related genes, and genetic improvement of ED. These results may help develop important genetic resources for enhancing maize yield through marker-assisted breeding.

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.

Fast integration and accumulation of beneficial breeding alleles through an AB-NAMIC strategy in wheat

Wheat (Triticum aestivum) is among the most important staple crops for safeguarding the food security of the growing world population. To bridge the gap between genebank diversity and breeding programs, we developed an advanced backcross-nested association mapping plus inter-crossed population (AB-NAMIC) by crossing three popular wheat cultivars as recurrent founders to 20 germplasm lines from a mini core collection. Selective backcrossing combined with selection against undesirable traits and extensive crossing within and between sub-populations created new opportunities to detect unknown genes and increase the frequency of beneficial alleles in the AB-NAMIC population. We performed phenotyping of 590 AB-NAMIC lines and a natural panel of 476 cultivars for six consecutive growing seasons and genotyped these 1066 lines with a 660K SNP array. Genome-wide association studies of both panels for plant development and yield traits demonstrated improved power to detect rare alleles and loci with medium genetic effects in AB-NAMIC. Notably, genome-wide association studies in AB-NAMIC detected the candidate gene TaSWEET6-7B (TraesCS7B03G1216700), which has high homology to the rice SWEET6b gene and exerts strong effects on adaptation and yield traits. The commercial release of two derived AB-NAMIC lines attests to its direct applicability in wheat improvement. Valuable information on genome-wide association study mapping, candidate genes, and their haplotypes for breeding traits are available through WheatGAB. Our research provides an excellent framework for fast-tracking exploration and accumulation of beneficial alleles stored in genebanks.

Hybrid allele-specific ChIP-seq analysis identifies variation in brassinosteroid-responsive transcription factor binding linked to traits in maize

BACKGROUND

Genetic variation in regulatory sequences that alter transcription factor (TF) binding is a major cause of phenotypic diversity. Brassinosteroid is a growth hormone that has major effects on plant phenotypes. Genetic variation in brassinosteroid-responsive cis-elements likely contributes to trait variation. Pinpointing such regulatory variations and quantitative genomic analysis of the variation in TF-target binding, however, remains challenging. How variation in transcriptional targets of signaling pathways such as the brassinosteroid pathway contributes to phenotypic variation is an important question to be investigated with innovative approaches.

RESULTS

Here, we use a hybrid allele-specific chromatin binding sequencing (HASCh-seq) approach and identify variations in target binding of the brassinosteroid-responsive TF ZmBZR1 in maize. HASCh-seq in the B73xMo17 F1s identifies thousands of target genes of ZmBZR1. Allele-specific ZmBZR1 binding (ASB) has been observed for 18.3% of target genes and is enriched in promoter and enhancer regions. About a quarter of the ASB sites correlate with sequence variation in BZR1-binding motifs and another quarter correlate with haplotype-specific DNA methylation, suggesting that both genetic and epigenetic variations contribute to the high level of variation in ZmBZR1 occupancy. Comparison with GWAS data shows linkage of hundreds of ASB loci to important yield and disease-related traits.

CONCLUSION

Our study provides a robust method for analyzing genome-wide variations of TF occupancy and identifies genetic and epigenetic variations of the brassinosteroid response transcription network in maize.

A Combination of QTL Mapping and GradedPool-Seq to Dissect Genetic Complexity for Gibberella Ear Rot Resistance in Maize Using an IBM Syn10 DH Population

Gibberella ear rot (GER) caused by Fusarium graminearum (teleomorph Gibberella zeae) is one of the most devastating maize diseases that reduces grain yield and quality worldwide. Utilization of host genetic resistance has become one of the most suitable strategies to control GER. In this study, a set of 246 diverse inbred lines derived from the intermated B 73 × Mo 17 doubled haploid population (IBM Syn10 DH) were used to detect quantitative trait loci (QTL) associated with resistance to GER. Meanwhile, a GradedPool-Seq (GPS) approach was performed to identify genomic variations involved in GER resistance. Using artificial inoculation across multiple environments, GER severity of the population was observed with wide phenotypic variation. Based on the linkage mapping, a total of 14 resistant QTLs were detected, accounting for 5.11 to 14.87% of the phenotypic variation, respectively. In GPS analysis, five significant single nucleotide polymorphisms (SNPs) associated with GER resistance were identified. Combining QTL mapping and GPS analysis, a peak-value SNP on chromosome 4 from GPS was overlapped with the QTL qGER4.2, suggesting that the colocalized region could be the most possible target location conferring resistance to GER. Subsequently, seven candidate genes were identified within the peak SNP, linking them to GER resistance. These findings are useful for exploring the complicated genetic variations in maize GER resistance. The genomic regions and genes identified herein provide a list of candidate targets for further investigation, in addition to the combined strategy that can be used for quantitative traits in plant species.

A multi-locus linear mixed model methodology for detecting small-effect QTLs for quantitative traits in MAGIC, NAM, and ROAM populations

Although multi-parent populations (MPPs) integrate the advantages of linkage and association mapping populations in the genetic dissection of complex traits and especially combine genetic analysis with crop breeding, it is difficult to detect small-effect quantitative trait loci (QTL) for complex traits in multiparent advanced generation intercross (MAGIC), nested association mapping (NAM), and random-open-parent association mapping (ROAM) populations. To address this issue, here we proposed a multi-locus linear mixed model method, namely mppQTL, to detect QTLs, especially small-effect QTLs, in these MPPs. The new method includes two steps. The first is genome-wide scanning based on a single-locus linear mixed model; the P-values are obtained from likelihood-ratio test, the peaks of negative logarithm P-value curve are selected by group-lasso, and all the selected peaks are regarded as potential QTLs. In the second step, all the potential QTLs are placed on a multi-locus linear mixed model, all the effects are estimated using expectation-maximization empirical Bayes algorithm, and all the non-zero effect vectors are further evaluated via likelihood-ratio test for significant QTLs. In Monte Carlo simulation studies, the new method has higher power in QTL detection, lower false positive rate, lower mean absolute deviation for QTL position estimate, and lower mean squared error for the estimate of QTL size (r²) than existing methods because the new method increases the power of detecting small-effect QTLs. In real dataset analysis, the new method (19) identified five more known genes than the existing three methods (14). This study provides an effective method for detecting small-effect QTLs in any MPPs.

Genetic architecture of protein and oil content in soybean seed and meal

Soybean is grown primarily for the protein and oil extracted from its seed and its value is influenced by these components. The objective of this study was to map marker-trait associations (MTAs) for the concentration of seed protein, oil, and meal protein using the soybean nested association mapping (SoyNAM) population. The composition traits were evaluated on seed harvested from over 5000 inbred lines of the SoyNAM population grown in 10 field locations across 3 years. Estimated heritabilities were at least 0.85 for all three traits. The genotyping of lines with single nucleotide polymorphism markers resulted in the identification of 107 MTAs for the three traits. When MTAs for the three traits that mapped within 5 cM intervals were binned together, the MTAs were mapped to 64 intervals on 19 of the 20 soybean chromosomes. The majority of the MTA effects were small and of the 107 MTAs, 37 were for protein content, 39 for meal protein, and 31 for oil content. For cases where a protein and oil MTAs mapped to the same interval, most (94%) significant effects were opposite for the two traits, consistent with the negative correlation between these traits. A coexpression analysis identified candidate genes linked to MTAs and 18 candidate genes were identified. The large number of small effect MTAs for the composition traits suggest that genomic prediction would be more effective in improving these traits than marker-assisted selection.

Genetic control of morphological traits useful for improving sorghum

Global climate change and global warming, coupled with the growing population, have raised concerns about sustainable food supply and bioenergy demand. Sorghum [Sorghum bicolor (L.) Moench] ranks fifth among cereals produced worldwide; it is a C₄ crop with a higher stress tolerance than other major cereals and has a wide range of uses, such as grains, forage, and biomass. Therefore, sorghum has attracted attention as a promising crop for achieving sustainable development goals (SDGs). In addition, sorghum is a suitable genetic model for C₄ grasses because of its high morphological diversity and relatively small genome size compared to other C₄ grasses. Although sorghum breeding and genetic studies have lagged compared to other crops such as rice and maize, recent advances in research have identified several genes and many quantitative trait loci (QTLs) that control important agronomic traits in sorghum. This review outlines traits and genetic information with a focus on morphogenetic aspects that may be useful in sorghum breeding for grain and biomass utilization.

High-Density Linkage Maps from Japanese Rice japonica Recombinant Inbred Lines Using Genotyping by Random Amplicon Sequencing-Direct (GRAS-Di)

The genetic dissection of agronomically important traits in closely related Japanese rice cultivars is still in its infancy mainly because of the narrow genetic diversity within japonica rice cultivars. In an attempt to unveil potential polymorphism between closely related Japanese rice cultivars, we used a next-generation-sequencing-based genotyping method: genotyping by random amplicon sequencing-direct (GRAS-Di) to develop genetic linkage maps. In this study, four recombinant inbred line (RIL) populations and their parents were used. A final RIL number of 190 for RIL71, 96 for RIL98, 95 for RIL16, and 94 for RIL91 derived from crosses between a common leading Japanese rice cultivar Koshihikari and Yamadanishiki, Taichung 65, Fujisaka 5, and Futaba, respectively, and the parent plants were subjected to GRAS-Di library construction and sequencing. Approximately 438.7 Mbp, 440 Mbp, 403.1 Mbp, and 392 Mbp called bases covering 97.5%, 97.3%, 98.3%, and 96.1%, respectively, of the estimated rice genome sequence at average depth of 1× were generated. Analysis of genotypic data identified 1050, 1285, 1708, and 1704 markers for each of the above RIL populations, respectively. Markers generated by GRAS-Di were organized into linkage maps and compared with those generated by GoldenGate SNP assay of the same RIL populations; the average genetic distance between markers showed a clear decrease in the four RIL populations when we integrated markers of both linkage maps. Genetic studies using these markers successfully localized five QTLs associated with heading date on chromosomes 3, 6, and 7 and which previously were identified as Hd1, Hd2, Hd6, Hd16, and Hd17. Therefore, GRAS-Di technology provided a low cost and efficient genotyping to overcome the narrow genetic diversity in closely related Japanese rice cultivars and enabled us to generate a high density linkage map in this germplasm.

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.

Important wheat diseases in the US and their management in the 21st century

Wheat is a crop of historical significance, as it marks the turning point of human civilization 10,000 years ago with its domestication. Due to the rapid increase in population, wheat production needs to be increased by 50% by 2050 and this growth will be mainly based on yield increases, as there is strong competition for scarce productive arable land from other sectors. This increasing demand can be further achieved using sustainable approaches including integrated disease pest management, adaption to warmer climates, less use of water resources and increased frequency of abiotic stress tolerances. Out of 200 diseases of wheat, 50 cause economic losses and are widely distributed. Each year, about 20% of wheat is lost due to diseases. Some major wheat diseases are rusts, smut, tan spot, spot blotch, fusarium head blight, common root rot, septoria blotch, powdery mildew, blast, and several viral, nematode, and bacterial diseases. These diseases badly impact the yield and cause mortality of the plants. This review focuses on important diseases of the wheat present in the United States, with comprehensive information of causal organism, economic damage, symptoms and host range, favorable conditions, and disease management strategies. Furthermore, major genetic and breeding efforts to control and manage these diseases are discussed. A detailed description of all the QTLs, genes reported and cloned for these diseases are provided in this review. This study will be of utmost importance to wheat breeding programs throughout the world to breed for resistance under changing environmental conditions.

Incorporating genome-wide and transcriptome-wide association studies to identify genetic elements of longissimus dorsi muscle in Huaxi cattle

Locating the genetic variation of important livestock and poultry economic traits is essential for genetic improvement in breeding programs. Identifying the candidate genes for the productive ability of Huaxi cattle was one crucial element for practical breeding. Based on the genotype and phenotype data of 1,478 individuals and the RNA-seq data of 120 individuals contained in 1,478 individuals, we implemented genome-wide association studies (GWAS), transcriptome-wide association studies (TWAS), and Fisher's combined test (FCT) to identify the candidate genes for the carcass trait, the weight of longissimus dorsi muscle (LDM). The results indicated that GWAS, TWAS, and FCT identified seven candidate genes for LDM altogether: PENK was located by GWAS and FCT, PPAT was located by TWAS and FCT, and XKR4, MTMR3, FGFRL1, DHRS4, and LAP3 were only located by one of the methods. After functional analysis of these candidate genes and referring to the reported studies, we found that they were mainly functional in the progress of the development of the body and the growth of muscle cells. Combining advanced breeding techniques such as gene editing with our study will significantly accelerate the genetic improvement for the future breeding of Huaxi cattle.

Breeding crops for drought-affected environments and improved climate resilience

Breeding climate-resilient crops with improved levels of abiotic and biotic stress resistance as a response to climate change presents both opportunities and challenges. Applying the framework of the "breeder's equation," which is used to predict the response to selection for a breeding program cycle, we review methodologies and strategies that have been used to successfully breed crops with improved levels of drought resistance, where the target population of environments (TPEs) is a spatially and temporally heterogeneous mixture of drought-affected and favorable (water-sufficient) environments. Long-term improvement of temperate maize for the US corn belt is used as a case study and compared with progress for other crops and geographies. Integration of trait information across scales, from genomes to ecosystems, is needed to accurately predict yield outcomes for genotypes within the current and future TPEs. This will require transdisciplinary teams to explore, identify, and exploit novel opportunities to accelerate breeding program outcomes; both improved germplasm resources and improved products (cultivars, hybrids, clones, and populations) that outperform and replace the products in use by farmers, in combination with modified agronomic management strategies suited to their local environments.

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.

Identification of Drought-Tolerance Genes in the Germination Stage of Soybean

Drought stress influences the vigor of plant seeds and inhibits seed germination, making it one of the primary environmental factors adversely affecting food security. The seed germination stage is critical to ensuring the growth and productivity of soybeans in soils prone to drought conditions. We here examined the genetic diversity and drought-tolerance phenotypes of 410 accessions of a germplasm diversity panel for soybean and conducted quantitative genetics analyses to identify loci associated with drought tolerance of seed germination. We uncovered significant differences among the diverse genotypes for four growth indices and five drought-tolerance indices, which revealed abundant variation among genotypes, upon drought stress, and for genotype × treatment effects. We also used 158,327 SNP markers and performed GWAS for the drought-related traits. Our data met the conditions (PCA + K) for using a mixed linear model in TASSEL, and we thus identified 26 SNPs associated with drought tolerance indices for germination stage distributed across 10 chromosomes. Nine SNP sites, including, for example, Gm20_34956219 and Gm20_36902659, were associated with two or more phenotypic indices, and there were nine SNP markers located in or adjacent to (within 500 kb) previously reported drought tolerance QTLs. These SNPs led to our identification of 41 candidate genes related to drought tolerance in the germination stage. The results of our study contribute to a deeper understanding of the genetic mechanisms underlying drought tolerance in soybeans at the germination stage, thereby providing a molecular basis for identifying useful soybean germplasm for breeding new drought-tolerant varieties.

The double round-robin population unravels the genetic architecture of grain size in barley

Grain number, size and weight primarily determine the yield of barley. Although the genes regulating grain number are well studied in barley, the genetic loci and the causal gene for sink capacity are poorly understood. Therefore, the primary objective of our work was to dissect the genetic architecture of grain size and weight in barley. We used a multi-parent population developed from a genetic cross between 23 diverse barley inbreds in a double round-robin design. Seed size-related parameters such as grain length, grain width, grain area and thousand-grain weight were evaluated in the HvDRR population comprising 45 recombinant inbred line sub-populations. We found significant genotypic variation for all seed size characteristics, and observed 84% or higher heritability across four environments. The quantitative trait locus (QTL) detection results indicate that the genetic architecture of grain size is more complex than previously reported. In addition, both cultivars and landraces contributed positive alleles at grain size QTLs. Candidate genes identified using genome-wide variant calling data for all parental inbred lines indicated overlapping and potential novel regulators of grain size in cereals. Furthermore, our results indicated that sink capacity was the primary determinant of grain weight in barley.

Nested association mapping-based GWAS for grain yield and related traits in wheat grown under diverse Australian environments

Utilising a nested association mapping (NAM) population-based GWAS, 98 stable marker-trait associations with 127 alleles unique to the exotic parents were detected for grain yield and related traits in wheat. Grain yield, thousand-grain weight, screenings and hectolitre weight are important wheat yield traits. An understanding of their genetic basis is crucial for improving grain yield in breeding programmes. Nested association mapping (NAM) populations are useful resources for the dissection of the genetic basis of complex traits such as grain yield and related traits in wheat. Coupled with phenotypic data collected from multiple environments, NAM populations have the power to detect quantitative trait loci and their multiple alleles, providing germplasm that can be incorporated into breeding programmes. In this study, we evaluated a large-scale wheat NAM population with two recurrent parents in unbalanced trials in nine diverse Australian field environments over three years. By applying a single-stage factor analytical linear mixed model (FALMM) to the NAM multi-environment trials (MET) data and conducting a genome-wide association study (GWAS), we detected 98 stable marker-trait associations (MTAs) with their multiple alleles. 74 MTAs had 127 alleles that were derived from the exotic parents and were absent in either of the two recurrent parents. The exotic alleles had favourable effects on 46 MTAs of the 74 MTAs, for grain yield, thousand-grain weight, screenings and hectolitre weight. Two NAM RILs with consistently high yield in multiple environments were also identified, highlighting the potential of the NAM population in supporting plant breeding through provision of germplasm that can be readily incorporated into breeding programmes. The identified beneficial exotic alleles introgressed into the NAM population provide potential target alleles for the genetic improvement of wheat and further studies aimed at pinpointing the underlying genes.

Improving drought tolerance in maize: Tools and techniques

Drought is an important constraint to agricultural productivity worldwide and is expected to worsen with climate change. To assist farmers, especially in sub-Saharan Africa (SSA), to adapt to climate change, continuous generation of stress-tolerant and farmer-preferred crop varieties, and their adoption by farmers, is critical to curb food insecurity. Maize is the most widely grown staple crop in SSA and plays a significant role in food security. The aim of this review is to present an overview of a broad range of tools and techniques used to improve drought tolerance in maize. We also present a summary of progress in breeding for maize drought tolerance, while incorporating research findings from disciplines such as physiology, molecular biology, and systems modeling. The review is expected to complement existing knowledge about breeding maize for climate resilience. Collaborative maize drought tolerance breeding projects in SSA emphasize the value of public-private partnerships in increasing access to genomic techniques and useful transgenes. To sustain the impact of maize drought tolerance projects in SSA, there must be complementary efforts to train the next generation of plant breeders and crop scientists.