GWAS: Fast-forwarding gene identification and characterization in temperate Cereals: lessons from Barley - A review

Drought-Tolerance QTLs Associated with Grain Yield and Related Traits in Spring Bread Wheat

The present research aims to identify the efficient combination of drought-tolerance selection criteria and associated quantitative trait loci. A panel of 197 bread wheat genotypes was evaluated for yield- and drought-tolerance-related traits in two environments (favorable and semiarid) for 2 years (2015-2016). Grain number, biomass, number of fertile spikes per plant and ground cover exhibited a significant correlation with grain yield and constitute potential secondary selection criteria for yield under drought conditions. About 73 significant marker-trait associations were detected along various chromosomal positions. The markers "wsnp_Ex_Rep_c67786_66472676" and "ExcalibuR_c24593_1217" exhibited important genetic gains associated with yield increase under drought (11 and 7%, respectively). The markers "KukRi_c94792_127" and "wsnp_Ex_c298_580660" showed a significant correlation with grain yield, biomass and grain number and were associated with a significant increase in yield performance at the semiarid site (+6 and +7%, respectively). The ground cover was found associated with grain yield and biomass through the markers "wsnp_Ex_Rep_c67786_66472676" (+11%) and "KukRi_c49927_151" (+10%). One marker "TduRuM_contig25432_1377" on chromosome 5B at 20 cM was consistently correlated with the number of fertile spikes across both environments. Further research should be considered to validate the efficiency of these markers to undertake selection for drought tolerance under various environments and genetic backgrounds.

High-LD SNP markers exhibiting pleiotropic effects on salt tolerance at germination and seedlings stages in spring wheat

Salt tolerance at germination and seedling growth stages was investigated. GWAS revealed nine genomic regions with pleiotropic effects on salt tolerance. Salt tolerant genotypes were identified for future breeding program. With 20% of the irrigated land worldwide affected by it, salinity is a serious threat to plant development and crop production. While wheat is the most stable food source worldwide, it has been classified as moderately tolerant to salinity. In several crop plants; such as barley, maize and rice, it has been shown that salinity tolerance at seed germination and seedling establishment is under polygenic control. As yield was the ultimate goal of breeders and geneticists, less attention has been paid to understanding the genetic architecture of salt tolerance at early stages. Thus, the genetic control of salt tolerance at these stages is poorly understood relative to the late stages. In the current study, 176 genotypes of spring wheat were tested for salinity tolerance at seed germination and seedling establishment. Genome-Wide Association Study (GWAS) has been used to identify the genomic regions/genes conferring salt tolerance at seed germination and seedling establishment. Salinity stress negatively impacted all germination and seedling development parameters. A set of 137 SNPs showed significant association with the traits of interest. Across the whole genome, 33 regions showed high linkage disequilibrium (LD). These high LD regions harbored 15 SNPs with pleiotropic effect (i.e. SNPs that control more than one trait). Nine genes belonging to different functional groups were found to be associated with the pleiotropic SNPs. Noteworthy, chromosome 2B harbored the gene TraesCS2B02G135900 that acts as a potassium transporter. Remarkably, one SNP marker, reported in an early study, associated with salt tolerance was validated in this study. Our findings represent potential targets of genetic manipulation to understand and improve salinity tolerance in wheat.

Genetic analysis toward more nutritious barley grains for a food secure world

BACKGROUND

Understanding the relationships between nutrition, human health and plant food source is among the highest priorities for public health. Therefore, enhancing the minerals content such as iron (Fe), zinc (Zn) and selenium (Se) in barley (Hordeum vulgare L.) grains is an urgent need to improve the nutritive value of barley grains in overcoming malnutrition and its potential consequencing. This study aimed to expedite biofortification of barley grains by elucidating the genetic basis of Zn, Fe, and Se accumulation in the grains, which will contribute to improved barley nutritional quality.

RESULTS

A genome-wide association study (GWAS) was conducted to detect the genetic architecture for grain Zn, Fe, and Se accumulations in 216 spring barley accessions across two years. All the accessions were genotyped by single nucleotide polymorphisms (SNPs) molecular markers. Mineral heritability values ranging from moderate to high were revealed in both environments. Remarkably, there was a high natural phenotypic variation for all micronutrient accumulation in the used population. High-LD SNP markers (222 SNPs) were detected to be associated with all micronutrients in barley grains across the two environments plus BLUEs. Three genomic regions were detected based on LD, which were identified for the most effective markers that had associations with more than one trait. The strongest SNP-trait associations were found to be physically located within genes that may be involved in grain Zn and Fe homeostasis. Two putative candidate genes were annotated as Basic helix loop helix (BHLH) family transcription factor and Squamosa promoter binding-like protein, respectively, and have been suggested as candidates for increased grain Zn, Fe, and Se accumulation.

CONCLUSIONS

These findings shed a light on the genetic basis of Zn, Fe, and Se accumulation in barley grains and have the potential to assist plant breeders in selecting accessions with high micronutrient concentrations to enhance grain quality and, ultimately human health.

Genome-wide association study for lignocellulosic compounds and fermentable sugar in rice straw

Cellulose and lignin are the two main components of secondary plant cell walls with substantial impact on stalk in the field and on straw during industrial processing. The amount of fermentable sugar that can be accessed is another important parameter affecting various industrial applications. In the present study, genetic variability of rice (Oryza sativa L.) genotypes for cellulose, lignin, and fermentable sugars contents was analyzed in rice straw. A genome-wide association study of 33,484 single nucleotide polymorphisms (SNPs) with a minor allele frequency (MAF) >0.05 was performed. The genome-wide association study identified seven, three, and three genomic regions to be significantly associated with cellulose, lignin, and fermentable sugar contents, respectively. Candidate genes in the associated genomic regions were enzymes mainly involved in cell wall metabolism. Novel SNP markers associated with cellulose were tagged to GH16, peroxidase, GT6, GT8, and CSLD2. For lignin content, Villin protein, OsWAK1/50/52/53, and GH16 were identified. For fermentable sugar content, UTP-glucose-1-phosphate uridylyltransferase, BRASSINOSTEROID INSENSITIVE 1, and receptor-like protein kinase 5 were found. The results of this study should improve our understanding of the genetic basis of the factors that might be involved in biosynthesis, turnover, and modification of major cell wall components and saccharides in rice straw.

Genome-wide association study reveals white lupin candidate gene involved in anthracnose resistance

GWAS identifies candidate gene controlling resistance to anthracnose disease in white lupin. White lupin (Lupinus albus L.) is a promising grain legume to meet the growing demand for plant-based protein. Its cultivation, however, is severely threatened by anthracnose disease caused by the fungal pathogen Colletotrichum lupini. To dissect the genetic architecture for anthracnose resistance, genotyping by sequencing was performed on white lupin accessions collected from the center of domestication and traditional cultivation regions. GBS resulted in 4611 high-quality single-nucleotide polymorphisms (SNPs) for 181 accessions, which were combined with resistance data observed under controlled conditions to perform a genome-wide association study (GWAS). Obtained disease phenotypes were shown to highly correlate with overall three-year disease assessments under Swiss field conditions (r > 0.8). GWAS results identified two significant SNPs associated with anthracnose resistance on gene Lalb_Chr05_g0216161 encoding a RING zinc-finger E3 ubiquitin ligase which is potentially involved in plant immunity. Population analysis showed a remarkably fast linkage disequilibrium decay, weak population structure and grouping of commercial varieties with landraces, corresponding to the slow domestication history and scarcity of modern breeding efforts in white lupin. Together with 15 highly resistant accessions identified in the resistance assay, our findings show promise for further crop improvement. This study provides the basis for marker-assisted selection, genomic prediction and studies aimed at understanding anthracnose resistance mechanisms in white lupin and contributes to improving breeding programs worldwide.

SoySNP618K array: A high-resolution single nucleotide polymorphism platform as a valuable genomic resource for soybean genetics and breeding

Innovations in genomics have enabled the development of low-cost, high-resolution, single nucleotide polymorphism (SNP) genotyping arrays that accelerate breeding progress and support basic research in crop science. Here, we developed and validated the SoySNP618K array (618,888 SNPs) for the important crop soybean. The SNPs were selected from whole-genome resequencing data containing 2,214 diverse soybean accessions; 29.34% of the SNPs mapped to genic regions representing 86.85% of the 56,044 annotated high-confidence genes. Identity-by-state analyses of 318 soybeans revealed 17 redundant accessions, highlighting the potential of the SoySNP618K array in supporting gene bank management. The patterns of population stratification and genomic regions enriched through domestication were highly consistent with previous findings based on resequencing data, suggesting that the ascertainment bias in the SoySNP618K array was largely compensated for. Genome-wide association mapping in combination with reported quantitative trait loci enabled fine-mapping of genes known to influence flowering time, E2 and GmPRR3b, and of a new candidate gene, GmVIP5. Moreover, genomic prediction of flowering and maturity time in 502 recombinant inbred lines was highly accurate (>0.65). Thus, the SoySNP618K array is a valuable genomic tool that can be used to address many questions in applied breeding, germplasm management, and basic crop research.

Genome-Wide Association Study Reveals Genetic Architecture and Candidate Genes for Yield and Related Traits under Terminal Drought, Combined Heat and Drought in Tropical Maize Germplasm

Maize (Zea mays L.) production is constrained by drought and heat stresses. The combination of these two stresses is likely to be more detrimental. To breed for maize cultivars tolerant of these stresses, 162 tropical maize inbred lines were evaluated under combined heat and drought (CHD) and terminal drought (TD) conditions. The mixed linear model was employed for the genome-wide association study using 7834 SNP markers and several phenotypic data including, days to 50% anthesis (AD) and silking (SD), husk cover (HUSKC), and grain yield (GY). In total, 66, 27, and 24 SNPs were associated with the traits evaluated under CHD, TD, and their combined effects, respectively. Of these, four single nucleotide polymorphism (SNP) markers (SNP_161703060 on Chr01, SNP_196800695 on Chr02, SNP_195454836 on Chr05, and SNP_51772182 on Chr07) had pleiotropic effects on both AD and SD under CHD conditions. Four SNPs (SNP_138825271 (Chr03), SNP_244895453 (Chr04), SNP_168561609 (Chr05), and SNP_62970998 (Chr06)) were associated with AD, SD, and HUSKC under TD. Twelve candidate genes containing phytohormone cis-acting regulating elements were implicated in the regulation of plant responses to multiple stress conditions including heat and drought. The SNPs and candidate genes identified in the study will provide invaluable information for breeding climate smart maize varieties under tropical conditions following validation of the SNP markers.

Genetic characterization and genome-wide association mapping for stem rust resistance in spring bread wheat

Background

Emerging wheat stem rust races have become a major threat to global wheat production. Finding additional loci responsible for resistance to these races and incorporating them into currently cultivated varieties is the most economic and environmentally sound strategy to combat this problem. Thus, this study was aimed at characterizing the genetic diversity and identifying the genetic loci conferring resistance to the stem rust of wheat. To accomplish this, 245 elite lines introduced from the International Center for Agricultural Research in the Dry Areas (ICARDA) were evaluated under natural stem rust pressure in the field at the Debre Zeit Agricultural Research Center, Ethiopia. The single nucleotide polymorphisms (SNP) marker data was retrieved from a 15 K SNP wheat array. A mixed linear model was used to investigate the association between SNP markers and the best linear unbiased prediction (BLUP) values of the stem rust coefficient of infection (CI).

Results

Phenotypic analysis revealed that 46% of the lines had a coefficient of infection (CI) in a range of 0 to 19. Genome-wide average values of 0.38, 0.20, and 0.71 were recorded for Nei’s gene diversity, polymorphism information content, and major allele frequency, respectively. A total of 46 marker-trait associations (MTAs) encompassed within eleven quantitative trait loci (QTL) were detected on chromosomes 1B, 3A, 3B, 4A, 4B, and 5A for CI. Two major QTLs with –log₁₀ (p) ≥ 4 (EWYP1B.1 and EWYP1B.2) were discovered on chromosome 1B.

Conclusions

This study identified several novel markers associated with stem rust resistance in wheat with the potential to facilitate durable rust resistance development through marker-assisted selection. It is recommended that the resistant wheat genotypes identified in this study be used in the national wheat breeding programs to improve stem rust resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01030-4.

Genetic Insight Into the Insect Resistance in Bread Wheat Exploiting the Untapped Natural Diversity

Climate change is an undeniable threat to sustainable wheat production in the future as an increased temperature will significantly increase grain loss due to the increased number of generations per season of multivoltine species that are detrimental to plants. Among insects, orange wheat blossom midge (OWBM), yellow wheat blossom midge (YWBM), saddle gall midge (SGM), thrips, and frit fly (FF) are important wheat pests in the European environments, which can be managed by the development of resistant cultivars. This involves the identification, confirmation, and incorporation of insect resistance sources into new high-yielding cultivars. We used two diverse and unrelated wheat [winter wheat (WW) and spring wheat (SW)] panels to associate single-nucleotide polymorphism (SNP) markers with the mentioned pests using the tools of association mapping. All in all, a total of 645 and 123 significant associations were detected in WW and SW, respectively, which were confined to 246 quantitative trait loci. Many candidate genes were identified using the BLAST analysis of the sequences of associated SNPs. Some of them are involved in controlling the physical structures of plants such as stomatal immunity and closure, cuticular wax in leaf blade, whereas others are involved in the production of certain enzymes in response to biotic and abiotic stresses. To our knowledge, this is the first detailed investigation that deals with YWBM, SGM, thrips, and FF resistance genetics using the natural variation in wheat. The reported germplasm is also readily available to breeders across the world that can make rational decisions to breed for the pest resilience of their interest by including the resistant genotypes being reported.

Genome-wide association study reveals a genomic region on 5AL for salinity tolerance in wheat

Soil salinity is a major threat to crop productivity and quality worldwide. In order to reduce the negative effects of salinity stress, it is important to understand the genetic basis of salinity tolerance. Identifying new salinity tolerance QTL or genes is crucial for breeders to pyramid different tolerance mechanisms to improve crop adaptability to salinity. Being one of the major cereal crops, wheat is known as a salt-sensitive glycophyte and subject to substantial yield losses when grown in the presence of salt. In this study, both pot and tank experiments were conducted to investigate the genotypic variation present in 328 wheat varieties in their salinity tolerance at the vegetative stage. A Genome-Wide Association Studies (GWAS) were carried out to identify QTL conferring salinity tolerance through a mixed linear model. Six, five and eight significant marker-trait associations (MTAs) were identified from pot experiments, tank experiments and average damage scores, respectively. These markers are located on the wheat chromosomes 1B, 2B, 2D, 3A, 4B, and 5A. These tolerance alleles were additive in their effects and, when combined, increased tolerance to salinity. Candidate genes identified in these QTL regions encoded a diverse class of proteins involved in salinity tolerance in plants. A Na⁺/H⁺ exchanger and a potassium transporter on chromosome 5A (IWB30519) will be of a potential value for improvement of salt tolerance of wheat cultivars using marker assisted selection programs. Some useful genotypes, which showed consistent tolerance in different trials, can also be effectively used in breeding programs.

Dissection of the polygenic architecture of neuronal Aβ production using a large sample of individual iPSC lines derived from Alzheimer's disease patients

Genome-wide association studies have demonstrated that polygenic risks shape Alzheimer's disease (AD). To elucidate the polygenic architecture of AD phenotypes at a cellular level, we established induced pluripotent stem cells from 102 patients with AD, differentiated them into cortical neurons and conducted a genome-wide analysis of the neuronal production of amyloid β (Aβ). Using such a cellular dissection of polygenicity (CDiP) approach, we identified 24 significant genome-wide loci associated with alterations in Aβ production, including some loci not previously associated with AD, and confirmed the influence of some of the corresponding genes on Aβ levels by the use of small interfering RNA. CDiP genotype sets improved the predictions of amyloid positivity in the brains and cerebrospinal fluid of patients in the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort. Secondary analyses of exome sequencing data from the Japanese ADNI and the ADNI cohorts focused on the 24 CDiP-derived loci associated with alterations in Aβ led to the identification of rare AD variants in KCNMA1.

Genetic characterization of carrot root shape and size using genome-wide association analysis and genomic-estimated breeding values

The principal phenotypic determinants of market class in carrot-the size and shape of the root-are under primarily additive, but also highly polygenic, genetic control. The size and shape of carrot roots are the primary determinants not only of yield, but also market class. These quantitative phenotypes have historically been challenging to objectively evaluate, and thus subjective visual assessment of market class remains the primary method by which selection for these traits is performed. However, advancements in digital image analysis have recently made possible the high-throughput quantification of size and shape attributes. It is therefore now feasible to utilize modern methods of genetic analysis to investigate the genetic control of root morphology. To this end, this study utilized both genome wide association analysis (GWAS) and genomic-estimated breeding values (GEBVs) and demonstrated that the components of market class are highly polygenic traits, likely under the influence of many small effect QTL. Relatively large proportions of additive genetic variance for many of the component phenotypes support high predictive ability of GEBVs; average prediction ability across underlying market class traits was 0.67. GWAS identified multiple QTL for four of the phenotypes which compose market class: length, aspect ratio, maximum width, and root fill, a previously uncharacterized trait which represents the size-independent portion of carrot root shape. By combining digital image analysis with GWAS and GEBVs, this study represents a novel advance in our understanding of the genetic control of market class in carrot. The immediate practical utility and viability of genomic selection for carrot market class is also described, and concrete guidelines for the design of training populations are provided.

Genomic Analysis of Resistance to Fall Armyworm (Spodoptera frugiperda) in CIMMYT Maize Lines

The recent invasion, rapid spread, and widescale destruction of the maize crop by the fall armyworm (FAW; Spodoptera frugiperda (J.E. Smith)) is likely to worsen the food insecurity situation in Africa. In the present study, a set of 424 maize lines were screened for their responses to FAW under artificial infestation to dissect the genetic basis of resistance. All lines were evaluated for two seasons under screen houses and genotyped with the DArTseq platform. Foliar damage was rated on a scale of 1 (highly resistant) to 9 (highly susceptible) and scored at 7, 14, and 21 days after artificial infestation. Analyses of variance revealed significant genotypic and genotype by environment interaction variances for all traits. Heritability estimates for leaf damage scores were moderately high and ranged from 0.38 to 0.58. Grain yield was negatively correlated with a high magnitude to foliar damage scores, ear rot, and ear damage traits. The genome-wide association study (GWAS) revealed 56 significant marker–trait associations and the predicted functions of the putative candidate genes varied from a defense response to several genes of unknown function. Overall, the study revealed that native genetic resistance to FAW is quantitative in nature and is controlled by many loci with minor effects.

Broadening the horizon of crop research: a decade of advancements in plant molecular genetics to divulge phenotype governing genes

Advancements in sequencing, genotyping, and computational technologies during the last decade (2011-2020) enabled new forward-genetic approaches, which subdue the impediments of precise gene mapping in varied crops. The modern crop improvement programs rely heavily on two major steps-trait-associated QTL/gene/marker's identification and molecular breeding. Thus, it is vital for basic and translational crop research to identify genomic regions that govern the phenotype of interest. Until the advent of next-generation sequencing, the forward-genetic techniques were laborious and time-consuming. Over the last 10 years, advancements in the area of genome assembly, genotyping, large-scale data analysis, and statistical algorithms have led faster identification of genomic variations regulating the complex agronomic traits and pathogen resistance. In this review, we describe the latest developments in genome sequencing and genotyping along with a comprehensive evaluation of the last 10-year headways in forward-genetic techniques that have shifted the focus of plant research from model plants to diverse crops. We have classified the available molecular genetic methods under bulk-segregant analysis-based (QTL-seq, GradedPool-Seq, QTG-Seq, Exome QTL-seq, and RapMap), target sequence enrichment-based (RenSeq, AgRenSeq, and TACCA), and mutation-based groups (MutMap, NIKS algorithm, MutRenSeq, MutChromSeq), alongside improvements in classical mapping and genome-wide association analyses. Newer methods for outcrossing, heterozygous, and polyploid plant genetics have also been discussed. The use of k-mers has enriched the nature of genetic variants which can be utilized to identify the phenotype-causing genes, independent of reference genomes. We envisage that the recent methods discussed herein will expand the repertoire of useful alleles and help in developing high-yielding and climate-resilient crops.

Mini review: Targeting below-ground plant performance to improve nitrogen use efficiency (NUE) in barley

Nitrogen (N) fertilizer is one of the major inputs for grain crops including barley and its usage is increasing globally. However, N use efficiency (NUE) is low in cereal crops, leading to higher production costs, unfulfilled grain yield potential and environmental hazards. N uptake is initiated from plant root tips but a very limited number of studies have been conducted on roots relevant to NUE specifically. In this review, we used barley, the fourth most important cereal crop, as the primary study plant to investigate this topic. We first highlighted the recent progress and study gaps in genetic analysis results, primarily, the genome-wide association study (GWAS) regarding both biological and statistical considerations. In addition, different factors contributing to NUE are discussed in terms of root morphological and anatomical traits, as well as physiological mechanisms such as N transporter activities and hormonal regulation.

Preparation and Curation of Multiyear, Multilocation, Multitrait Datasets

Genome-wide association studies (GWAS) are a powerful approach to dissect genotype-phenotype associations and identify causative regions. However, this power is highly influenced by the accuracy of the phenotypic data. To obtain accurate phenotypic values, the phenotyping should be achieved through multienvironment trials (METs). In order to avoid any technical errors, the required time needs to be spent on exploring, understanding, curating and adjusting the phenotypic data in each trial before combining them using an appropriate linear mixed model (LMM). The LMM is chosen to minimize as much as possible any effect that can lead to misestimation of the phenotypic values. The purpose of this chapter is to explain a series of important steps to explore and analyze data from METs used to characterize an association panel. Two datasets are used to illustrate two different scenarios.

Two decades of association mapping: Insights on disease resistance in major crops

Climate change across the globe has an impact on the occurrence, prevalence, and severity of plant diseases. About 30% of yield losses in major crops are due to plant diseases; emerging diseases are likely to worsen the sustainable production in the coming years. Plant diseases have led to increased hunger and mass migration of human populations in the past, thus a serious threat to global food security. Equipping the modern varieties/hybrids with enhanced genetic resistance is the most economic, sustainable and environmentally friendly solution. Plant geneticists have done tremendous work in identifying stable resistance in primary genepools and many times other than primary genepools to breed resistant varieties in different major crops. Over the last two decades, the availability of crop and pathogen genomes due to advances in next generation sequencing technologies improved our understanding of trait genetics using different approaches. Genome-wide association studies have been effectively used to identify candidate genes and map loci associated with different diseases in crop plants. In this review, we highlight successful examples for the discovery of resistance genes to many important diseases. In addition, major developments in association studies, statistical models and bioinformatic tools that improve the power, resolution and the efficiency of identifying marker-trait associations. Overall this review provides comprehensive insights into the two decades of advances in GWAS studies and discusses the challenges and opportunities this research area provides for breeding resistant varieties.

High Resolution Genome Wide Association Studies Reveal Rich Genetic Architectures of Grain Zinc and Iron in Common Wheat (Triticum aestivum L.)

Biofortification is a sustainable strategy to alleviate micronutrient deficiency in humans. It is necessary to improve grain zinc (GZnC) and iron concentrations (GFeC) in wheat based on genetic knowledge. However, the precise dissection of the genetic architecture underlying GZnC and GFeC remains challenging. In this study, high-resolution genome-wide association studies were conducted for GZnC and GFeC by three different models using 166 wheat cultivars and 373,106 polymorphic markers from the wheat 660K and 90K single nucleotide polymorphism (SNP) arrays. Totally, 25 and 16 stable loci were detected for GZnC and GFeC, respectively. Among them, 17 loci for GZnC and 8 for GFeC are likely to be new quantitative trait locus/loci (QTL). Based on gene annotations and expression profiles, 28 promising candidate genes were identified for Zn/Fe uptake (8), transport (11), storage (3), and regulations (6). Of them, 11 genes were putative wheat orthologs of known Arabidopsis and rice genes related to Zn/Fe homeostasis. A brief model, such as genes related to Zn/Fe homeostasis from root uptake, xylem transport to the final seed storage was proposed in wheat. Kompetitive allele-specific PCR (KASP) markers were successfully developed for two major QTL of GZnC on chromosome arms 3AL and 7AL, respectively, which were independent of thousand kernel weight and plant height. The 3AL QTL was further validated in a bi-parental population under multi-environments. A wheat multidrug and toxic compound extrusion (MATE) transporter TraesCS3A01G499300, the ortholog of rice gene OsPEZ2, was identified as a potential candidate gene. This study has advanced our knowledge of the genetic basis underlying GZnC and GFeC in wheat and provides valuable markers and candidate genes for wheat biofortification.

Advanced high-throughput plant phenotyping techniques for genome-wide association studies: A review

Linking phenotypes and genotypes to identify genetic architectures that regulate important traits is crucial for plant breeding and the development of plant genomics. In recent years, genome-wide association studies (GWASs) have been applied extensively to interpret relationships between genes and traits. Successful GWAS application requires comprehensive genomic and phenotypic data from large populations. Although multiple high-throughput DNA sequencing approaches are available for the generation of genomics data, the capacity to generate high-quality phenotypic data is lagging far behind. Traditional methods for plant phenotyping mostly rely on manual measurements, which are laborious, inaccurate, and time-consuming, greatly impairing the acquisition of phenotypic data from large populations. In contrast, high-throughput phenotyping has unique advantages, facilitating rapid, non-destructive, and high-throughput detection, and, in turn, addressing the shortcomings of traditional methods. Aim of Review: This review summarizes the current status with regard to the integration of high-throughput phenotyping and GWAS in plants, in addition to discussing the inherent challenges and future prospects. Key Scientific Concepts of Review: High-throughput phenotyping, which facilitates non-contact and dynamic measurements, has the potential to offer high-quality trait data for GWAS and, in turn, to enhance the unraveling of genetic structures of complex plant traits. In conclusion, high-throughput phenotyping integration with GWAS could facilitate the revealing of coding information in plant genomes.

Alternative splicing triggered by the insertion of a CACTA transposon attenuates LsGLK and leads to the development of pale-green leaves in lettuce

Lettuce (Lactuca sativa) is one of the most important vegetable crops in the world. As a leafy vegetable, the polymorphism of lettuce leaves from dark to pale green is an important trait. However, the genetic and molecular mechanisms underlying such variations remain poorly understood. In this study, one major locus controlling the polymorphism of dark- and pale-green leaves in lettuce was identified using genome-wide association studies (GWAS). This locus was then fine mapped to an interval of 5375 bp on chromosome 4 using a segregating population containing 2480 progeny. Only one gene, homologous to the GLK genes in Arabidopsis and other plants, is present in the candidate region. A complementation test confirmed that the candidate gene, LsGLK, contributes to the variation of dark- and pale-green leaves. Sequence analysis showed that a CACTA transposon of 7434 bp was inserted 10 bp downstream of the stop codon of LsGLK, followed by a duplication of a 1826-bp fragment covering exons 3-6 of the LsGLK gene. The transposon insertion did not change the expression level of the LsGLK gene. However, because of alternative splicing, only 6% of the transcripts produced from the transposon insertion were wild-type transcripts, which led to the production of pale-green leaves. An evolutionary analysis revealed that the insertion of the CACTA transposon occurred in cultivated lettuce and might have been selected in particular cultivars to satisfy the diverse demands of consumers. In this study, we demonstrated that a transposon insertion near a gene may affect its splicing and consequently generate phenotypic variations.

Discovery of a Novel Leaf Rust (Puccinia recondita) Resistance Gene in Rye (Secale cereale L.) Using Association Genomics

Leaf rust constitutes one of the most important foliar diseases in rye (Secale cereale L.). To discover new sources of resistance, we phenotyped 180 lines belonging to a less well-characterized Gülzow germplasm at three field trial locations in Denmark and Northern Germany in 2018 and 2019. We observed lines with high leaf rust resistance efficacy at all locations in both years. A genome-wide association study using 261,406 informative single-nucleotide polymorphisms revealed two genomic regions associated with resistance on chromosome arms 1RS and 7RS, respectively. The most resistance-associated marker on chromosome arm 1RS physically co-localized with molecular markers delimiting Pr3. In the reference genomes Lo7 and Weining, the genomic region associated with resistance on chromosome arm 7RS contained a large number of nucleotide-binding leucine-rich repeat (NLR) genes. Residing in close proximity to the most resistance-associated marker, we identified a cluster of NLRs exhibiting close protein sequence similarity with the wheat leaf rust Lr1 gene situated on chromosome arm 5DL in wheat, which is syntenic to chromosome arm 7RS in rye. Due to the close proximity to the most resistance-associated marker, our findings suggest that the considered leaf rust R gene, provisionally denoted Pr6, could be a Lr1 ortholog in rye.

Portuguese Common Bean Natural Variation Helps to Clarify the Genetic Architecture of the Legume's Nutritional Composition and Protein Quality

Common bean is a nutritious food legume widely appreciated by consumers worldwide. It is a staple food in Latin America, and a component of the Mediterranean diet, being an affordable source of protein with high potential as a gourmet food. Breeding for nutritional quality, including both macro and micronutrients, and meeting organoleptic consumers' preferences is a difficult task which is facilitated by uncovering the genetic basis of related traits. This study explored the diversity of 106 Portuguese common bean accessions, under two contrasting environments, to gain insight into the genetic basis of nutritional composition (ash, carbohydrates, fat, fiber, moisture, protein, and resistant starch contents) and protein quality (amino acid contents and trypsin inhibitor activity) traits through a genome-wide association study. Single-nucleotide polymorphism-trait associations were tested using linear mixed models accounting for the accessions' genetic relatedness. Mapping resolution to the gene level was achieved in 56% of the cases, with 102 candidate genes proposed for 136 genomic regions associated with trait variation. Only one marker-trait association was stable across environments, highlighting the associations' environment-specific nature and the importance of genotype × environment interaction for crops' local adaptation and quality. This study provides novel information to better understand the molecular mechanisms regulating the nutritional quality in common bean and promising molecular tools to aid future breeding efforts to answer consumers' concerns.

Pyramiding of scald resistance genes in four spring barley MAGIC populations

Genome-Wide Association Studies (GWAS) of four Multi-parent Advanced Generation Inter-Cross (MAGIC) populations identified nine regions on chromosomes 1H, 3H, 4H, 5H, 6H and 7H associated with resistance against barley scald disease. Three of these regions are putatively novel resistance Quantitative Trait Loci (QTL). Barley scald is caused by Rhynchosporium commune, one of the most important barley leaf diseases that are prevalent in most barley-growing regions. Up to 40% yield losses can occur in susceptible barley cultivars. Four MAGIC populations were generated in a Nordic Public-Private Pre-breeding of spring barley project (PPP Barley) to introduce resistance to several important diseases. Here, these MAGIC populations consisting of six to eight founders each were tested for scald resistance in field trials in Finland and Iceland. Eight different model covariate combinations were compared for GWAS studies, and the models that deviated the least from the expected p-values were selected. For all QTL, candidate genes were identified that are predicted to be involved in pathogen defence. The MAGIC progenies contained new haplotypes of significant SNP-markers with high resistance levels. The lines with successfully pyramided resistance against scald and mildew and the significant markers are now distributed among Nordic plant breeders and will benefit development of disease-resistant cultivars.

Omics-Facilitated Crop Improvement for Climate Resilience and Superior Nutritive Value

Novel crop improvement approaches, including those that facilitate for the exploitation of crop wild relatives and underutilized species harboring the much-needed natural allelic variation are indispensable if we are to develop climate-smart crops with enhanced abiotic and biotic stress tolerance, higher nutritive value, and superior traits of agronomic importance. Top among these approaches are the "omics" technologies, including genomics, transcriptomics, proteomics, metabolomics, phenomics, and their integration, whose deployment has been vital in revealing several key genes, proteins and metabolic pathways underlying numerous traits of agronomic importance, and aiding marker-assisted breeding in major crop species. Here, citing several relevant examples, we appraise our understanding on the recent developments in omics technologies and how they are driving our quest to breed climate resilient crops. Large-scale genome resequencing, pan-genomes and genome-wide association studies are aiding the identification and analysis of species-level genome variations, whilst RNA-sequencing driven transcriptomics has provided unprecedented opportunities for conducting crop abiotic and biotic stress response studies. Meanwhile, single cell transcriptomics is slowly becoming an indispensable tool for decoding cell-specific stress responses, although several technical and experimental design challenges still need to be resolved. Additionally, the refinement of the conventional techniques and advent of modern, high-resolution proteomics technologies necessitated a gradual shift from the general descriptive studies of plant protein abundances to large scale analysis of protein-metabolite interactions. Especially, metabolomics is currently receiving special attention, owing to the role metabolites play as metabolic intermediates and close links to the phenotypic expression. Further, high throughput phenomics applications are driving the targeting of new research domains such as root system architecture analysis, and exploration of plant root-associated microbes for improved crop health and climate resilience. Overall, coupling these multi-omics technologies to modern plant breeding and genetic engineering methods ensures an all-encompassing approach to developing nutritionally-rich and climate-smart crops whose productivity can sustainably and sufficiently meet the current and future food, nutrition and energy demands.

Genome-Wide Association Mapping for Heat and Drought Adaptive Traits in Pea

Heat and drought, individually or in combination, limit pea productivity. Fortunately, substantial genetic diversity exists in pea germplasm for traits related to abiotic stress resistance. Understanding the genetic basis of resistance could accelerate the development of stress-adaptive cultivars. We conducted a genome-wide association study (GWAS) in pea on six stress-adaptive traits with the aim to detect the genetic regions controlling these traits. One hundred and thirty-five genetically diverse pea accessions were phenotyped in field studies across three or five environments under stress and control conditions. To determine marker trait associations (MTAs), a total of 16,877 valuable single nucleotide polymorphisms (SNPs) were used in association analysis. Association mapping detected 15 MTAs that were significantly (p ≤ 0.0005) associated with the six stress-adaptive traits averaged across all environments and consistent in multiple individual environments. The identified MTAs were four for lamina wax, three for petiole wax, three for stem thickness, two for the flowering duration, one for the normalized difference vegetation index (NDVI), and two for the normalized pigment and chlorophyll index (NPCI). Sixteen candidate genes were identified within a 15 kb distance from either side of the markers. The detected MTAs and candidate genes have prospective use towards selecting stress-hardy pea cultivars in marker-assisted selection.

Genome-Wide Association Studies of Soybean Yield-Related Hyperspectral Reflectance Bands Using Machine Learning-Mediated Data Integration Methods

In conjunction with big data analysis methods, plant omics technologies have provided scientists with cost-effective and promising tools for discovering genetic architectures of complex agronomic traits using large breeding populations. In recent years, there has been significant progress in plant phenomics and genomics approaches for generating reliable large datasets. However, selecting an appropriate data integration and analysis method to improve the efficiency of phenome-phenome and phenome-genome association studies is still a bottleneck. This study proposes a hyperspectral wide association study (HypWAS) approach as a phenome-phenome association analysis through a hierarchical data integration strategy to estimate the prediction power of hyperspectral reflectance bands in predicting soybean seed yield. Using HypWAS, five important hyperspectral reflectance bands in visible, red-edge, and near-infrared regions were identified significantly associated with seed yield. The phenome-genome association analysis of each tested hyperspectral reflectance band was performed using two conventional genome-wide association studies (GWAS) methods and a machine learning mediated GWAS based on the support vector regression (SVR) method. Using SVR-mediated GWAS, more relevant QTL with the physiological background of the tested hyperspectral reflectance bands were detected, supported by the functional annotation of candidate gene analyses. The results of this study have indicated the advantages of using hierarchical data integration strategy and advanced mathematical methods coupled with phenome-phenome and phenome-genome association analyses for a better understanding of the biology and genetic backgrounds of hyperspectral reflectance bands affecting soybean yield formation. The identified yield-related hyperspectral reflectance bands using HypWAS can be used as indirect selection criteria for selecting superior genotypes with improved yield genetic gains in large breeding populations.

The Relationships between Plant Developmental Traits and Winter Field Survival in Rye (Secale cereale L.)

Overwintering cereals accumulate low temperature tolerance (LTT) during cold acclimation in the autumn. Simultaneously, the plants adjust to the colder season by making developmental changes at the shoot apical meristem. These processes lead to higher winter hardiness in winter rye varieties (Secale cereale L.) adapted to Northern latitudes as compared to other cereal crops. To dissect the winter-hardiness trait in rye, a panel of 96 genotypes of different origins and growth habits was assessed for winter field survival (WFS), LTT, and six developmental traits. Best Linear Unbiased Estimates for WFS determined from five field trials correlated strongly with LTT (r = 0.90, p < 0.001); thus, cold acclimation efficiency was the major contributor to WFS. WFS also correlated strongly (p < 0.001) with final leaf number (r = 0.80), prostrate growth habit (r = 0.61), plant height (r = 0.34), but showed weaker associations with top internode length (r = 0.30, p < 0.01) and days to anthesis (r = 0.25, p < 0.05). The heritability estimates (h²) for WFS-associated traits ranged from 0.45 (prostrate growth habit) to 0.81 (final leaf number) and were overall higher than for WFS (h² = 0.48). All developmental traits associated with WFS and LTT are postulated to be regulated by phytohormone levels at shoot apical meristem.

Identification of genomic regions affecting grain peroxidase activity in bread wheat using genome-wide association study

BACKGROUND

Peroxidase (POD) activity plays an important role in flour-based product quality, which is mainly associated with browning and bleaching effects of flour. Here, we performed a genome-wide association study (GWAS) on POD activity using an association population consisted with 207 wheat world-wide collected varieties. Our study also provide basis for the genetic improvement of flour color-based quality in wheat.

RESULTS

Twenty quantitative trait loci (QTLs) were detected associated with POD activity, explaining 5.59-12.67% of phenotypic variation. Superior alleles were positively correlated with POD activity. In addition, two SNPs were successfully developed to KASP (Kompetitive Allele-Specific PCR) markers. Two POD genes, TraesCS2B02G615700 and TraesCS2D02G583000, were aligned near the QTLs flanking genomic regions, but only TraesCS2D02G583000 displayed significant divergent expression levels (P < 0.001) between high and low POD activity varieties in the investigated association population. Therefore, it was deduced to be a candidate gene. The expression level of TraesCS2D02G583000 was assigned as a phenotype for expression GWAS (eGWAS) to screen regulatory elements. In total, 505 significant SNPs on 20 chromosomes (excluding 4D) were detected, and 9 of them located within 1 Mb interval of TraesCS2D02G583000.

CONCLUSIONS

To identify genetic loci affecting POD activity in wheat grain, we conducted GWAS on POD activity and the candidate gene TraesCS2D02G583000 expression. Finally, 20 QTLs were detected for POD activity, whereas two QTLs associated SNPs were converted to KASP markers that could be used for marker-assisted breeding. Both cis- and trans-acting elements were revealed by eGWAS of TraesCS2D02G583000 expression. The present study provides genetic loci for improving POD activity across wide genetic backgrounds and largely improved the selection efficiency for breeding in wheat.

Genomic regions associated with physiological, biochemical and yield-related responses under water deficit in diploid potato at the tuber initiation stage revealed by GWAS

Water deficit, which is increasing with climate change, is a serious threat to agricultural sustainability worldwide. Dissection of the genetic architecture of water deficit responses is highly desirable for developing water-deficit tolerant potato cultivars and enhancing the resilience of existing cultivars. This study examined genetic variation in response to water deficit in a panel of diploid potato and identified the QTL governing this trait via a genome-wide association study (GWAS). A panel of 104 diploid potato accessions were evaluated under both well-watered and water deficit treatments at tuber initiation stage. Drought stress index (DTI) was calculated to assess tolerance of the diploid potato genotypes to water deficit. The GWAS was conducted using a matrix of 47K single nucleotide polymorphisms (SNP), recently available for this population. We are reporting 38 QTL, seven for well-watered conditions, twenty-two for water deficit conditions and nine for DTI which explain between 12.6% and 44.1% of the phenotypic variance. A set of 6 QTL were found to be associated with more than one variable. Marker WDP-9.21 was found associated with tuber fresh weigh under WD and gene annotation analysis revealed co-localization with the Glucan/water dikinase (GWD) gene. Of the nine QTL detected from DTI on chromosomes 2,3,5,8,10 and 12, three candidate genes with a feasible role in water deficit response were identified. The findings of this study can be used in marker-assisted selection (MAS) for water- deficit tolerance breeding in potato.

Genetic factors controlling nTiO2 nanoparticles stress tolerance in barley (Hordeum vulgare) during seed germination and seedling development

Titanium dioxide nanoparticle (nTiO2) is one of the most produced nanoparticles worldwide. Its mechanism on crop development and performance is unclear as it is hard to predict their toxicity or benefit. Therefore, understanding the genetics of crop development under nTiO2 is a prerequisite for their applications in agriculture and crop improvement. Here, we aimed to examine the influnce of 300ppm nTiO2 on seed germination, seedling morphology, root-related traits in 121 worldwide spring barley (Hordeum vulgare L.) accessions. Results show that nTiO2 significantley affected all traits scored in this study. Response to nTiO2 treatment, clear wide natural variation among accesions was detected. Remarkably, 10 genotypes showed increased root length under nTiO2 at the seedling stage indicating that nTiO2 enhanced the root elongation. Genome-wide association scan (GWAS) was applied using 9K single nucleotide polymorphism (SNPs) in a mixed-linear model that revealed 86 significant marker-trait associations with all traits scored in this study. Many significant SNPs were physically located near candidate genes, of which 191 genes were detected within the linkage disequilibrium and distributed over all barley chromosomes. Mostly, the genes harboured by chromosome 2H, specially calcium-binding genes family, regulate the variation of seedling length-related traits. Candidate genes on 7H encode zinc finger protein that controls the rate of germination. Therefore, these genomic regions at 2H and 7H can be targeted to select for improved seedling development and seed germination under nTiO2 stress in soils. These results improve understanding the genetic control of seed germination and seedling development under high levels of nTiO2 that can support plant breeding and crop improvement programmes.

A Genome-Wide Association Study Pinpoints Quantitative Trait Genes for Plant Height, Heading Date, Grain Quality, and Yield in Rye (Secale cereale L.)

Rye is the only cross-pollinating Triticeae crop species. Knowledge of rye genes controlling complex-inherited traits is scarce, which, currently, largely disables the genomics assisted introgression of untapped genetic variation from self-incompatible germplasm collections in elite inbred lines for hybrid breeding. We report on the first genome-wide association study (GWAS) in rye based on the phenotypic evaluation of 526 experimental hybrids for plant height, heading date, grain quality, and yield in 2 years and up to 19 environments. We established a cross-validated NIRS calibration model as a fast, effective, and robust analytical method to determine grain quality parameters. We observed phenotypic plasticity in plant height and tiller number as a resource use strategy of rye under drought and identified increased grain arabinoxylan content as a striking phenotype in osmotically stressed rye. We used DArTseq™ as a genotyping-by-sequencing technology to reduce the complexity of the rye genome. We established a novel high-density genetic linkage map that describes the position of almost 19k markers and that allowed us to estimate a low genome-wide LD based on the assessed genetic diversity in elite germplasm. We analyzed the relationship between plant height, heading date, agronomic, as well as grain quality traits, and genotype based on 20k novel single-nucleotide polymorphism markers. In addition, we integrated the DArTseq™ markers in the recently established 'Lo7' reference genome assembly. We identified cross-validated SNPs in 'Lo7' protein-coding genes associated with all traits studied. These include associations of the WUSCHEL-related homeobox transcription factor DWT1 and grain yield, the DELLA protein gene SLR1 and heading date, the Ethylene overproducer 1-like protein gene ETOL1 and thousand-grain weight, protein and starch content, as well as the Lectin receptor kinase SIT2 and plant height. A Leucine-rich repeat receptor protein kinase and a Xyloglucan alpha-1,6-xylosyltransferase count among the cross-validated genes associated with water-extractable arabinoxylan content. This study demonstrates the power of GWAS, hybrid breeding, and the reference genome sequence in rye genetics research to dissect and identify the function of genes shaping genetic diversity in agronomic and grain quality traits of rye. The described links between genetic causes and phenotypic variation will accelerate genomics-enabled rye improvement.

Identification and Validation of High LD Hotspot Genomic Regions Harboring Stem Rust Resistant Genes on 1B, 2A (Sr38), and 7B Chromosomes in Wheat

Stem rust caused by Puccinia graminis f. sp. tritici Eriks. is an important disease of common wheat globally. The production and cultivation of genetically resistant cultivars are one of the most successful and environmentally friendly ways to protect wheat against fungal pathogens. Seedling screening and genome-wide association study (GWAS) were used to determine the genetic diversity of wheat genotypes obtained on stem rust resistance loci. At the seedling stage, the reaction of the common stem rust race QFCSC in Nebraska was measured in a set of 212 genotypes from F_3:6 lines. The results indicated that 184 genotypes (86.8%) had different degrees of resistance to this common race. While 28 genotypes (13.2%) were susceptible to stem rust. A set of 11,911 single-nucleotide polymorphism (SNP) markers was used to perform GWAS which detected 84 significant marker-trait associations (MTAs) with SNPs located on chromosomes 1B, 2A, 2B, 7B and an unknown chromosome. Promising high linkage disequilibrium (LD) genomic regions were found in all chromosomes except 2B which suggested they include candidate genes controlling stem rust resistance. Highly significant LD was found among these 59 significant SNPs on chromosome 2A and 12 significant SNPs with an unknown chromosomal position. The LD analysis between SNPs located on 2A and Sr38 gene reveal high significant LD genomic regions which was previously reported. To select the most promising stem rust resistant genotypes, a new approach was suggested based on four criteria including, phenotypic selection, number of resistant allele(s), the genetic distance among the selected parents, and number of the different resistant allele(s) in the candidate crosses. As a result, 23 genotypes were considered as the most suitable parents for crossing to produce highly resistant stem rust genotypes against the QFCSC.

Metabolomics for Crop Breeding: General Considerations

The development of new, more productive varieties of agricultural crops is becoming an increasingly difficult task. Modern approaches for the identification of beneficial alleles and their use in elite cultivars, such as quantitative trait loci (QTL) mapping and marker-assisted selection (MAS), are effective but insufficient for keeping pace with the improvement of wheat or other crops. Metabolomics is a powerful but underutilized approach that can assist crop breeding. In this review, basic methodological information is summarized, and the current strategies of applications of metabolomics related to crop breeding are explored using recent examples. We briefly describe classes of plant metabolites, cellular localization of metabolic pathways, and the strengths and weaknesses of the main metabolomics technique. Among the commercialized genetically modified crops, about 50 with altered metabolic enzyme activities have been identified in the International Service for the Acquisition of Agri-biotech Applications (ISAAA) database. These plants are reviewed as encouraging examples of the application of knowledge of biochemical pathways. Based on the recent examples of metabolomic studies, we discuss the performance of metabolic markers, the integration of metabolic and genomic data in metabolic QTLs (mQTLs) and metabolic genome-wide association studies (mGWAS). The elucidation of metabolic pathways and involved genes will help in crop breeding and the introgression of alleles of wild relatives in a more targeted manner.

Genome-Wide Association Study to Map Genomic Regions Related to the Initiation Time of Four Growth Stage Traits in Soybean

The time to flowering (DF), pod beginning (DPB), seed formation (DSF), and maturity initiation (DMI) in soybean (Glycine max [L.] Merr) are important characteristics of growth stage traits (GSTs) in Chinese summer-sowing soybean, and are influenced by genetic as well as environmental factors. To better understand the molecular mechanism underlying the initiation times of GSTs, we investigated four GSTs of 309 diverse soybean accessions in six different environments and Best Linear Unbiased Prediction values. Furthermore, the genome-wide association study was conducted by a Fixed and random model Circulating Probability Unification method using over 60,000 single nucleotide polymorphism (SNP) markers to identify the significant quantitative trait nucleotide (QTN) regions with phenotypic data. As a result, 212 SNPs within 102 QTN regions were associated with four GSTs. Of which, eight stable regions were repeatedly detected in least three datasets for one GST. Interestingly, half of the QTN regions overlapped with previously reported quantitative trait loci or well-known soybean growth period genes. The hotspots associated with all GSTs were concentrated on chromosome 10. E2 (Glyma10g36600), a gene with a known function in regulating flowering and maturity in soybean, is also found on this chromosome. Thus, this genomic region may account for the strong correlation among the four GSTs. All the significant SNPs in the remaining 7 QTN regions could cause the significant phenotypic variation with both the major and minor alleles. Two hundred and seventy-five genes in soybean and their homologs in Arabidopsis were screened within ± 500 kb of 7 peak SNPs in the corresponding QTN regions. Most of the genes are involved in flowering, response to auxin stimulus, or regulation of seed germination, among others. The findings reported here provide an insight for genetic improvement which will aid in breeding of soybean cultivars that can be adapted to the various summer sowing areas in China and beyond.

Toward identification of a putative candidate gene for nutrient mineral accumulation in wheat grains for human nutrition purposes

A multilocus genome-wide association study of a panel of 369 diverse wheat (Triticum aestivum) genotypes was carried out in order to examine the genetic basis of variations in nutrient mineral concentrations in the grains. The panel was grown under field conditions for three consecutive years and the concentrations of Ca, K, Mg, Mn, P, and S were determined. Wide ranges of natural variation were detected among the genotypes. Strong positive correlations were found among the minerals except for K, which showed negative correlation trends with the other minerals. Genetic association analysis detected 86 significant marker-trait associations (MTAs) underlying the natural variations in mineral concentrations in grains. The major MTA was detected on the long arm of chromosome 5A and showed a pleiotropic effect on Ca, K, Mg, Mn, and S. Further significant MTAs were distributed among the whole genome except for chromosomes 3D and 6D. We identified putative candidate genes that are potentially involved in metal uptake, transport, and assimilation, including TraesCS5A02G542600 on chromosome 5A, which was annotated as a Major Facilitator Superfamily transporter and acted on all the minerals except K. TraesCS5A02G542600 was highly expressed in seed coat, and to a lesser extent in the peduncle, awns, and lemma. Our results provide important insights into the genetic basis of enhancement of nutrient mineral concentrations that can help to inform future breeding studies in order to improve human nutrition.

Analysis of 427 genomes reveals moso bamboo population structure and genetic basis of property traits

Moso bamboo (Phyllostachys edulis) is an economically and ecologically important nontimber forestry species. Further development of this species as a sustainable bamboo resource has been hindered by a lack of population genome information. Here, we report a moso bamboo genomic variation atlas of 5.45 million single-nucleotide polymorphisms (SNPs) from whole-genome resequencing of 427 individuals covering 15 representative geographic areas. We uncover low genetic diversity, high genotype heterozygosity, and genes under balancing selection underlying moso bamboo population adaptation. We infer its demographic history with one bottleneck and its recently small population without a rebound. We define five phylogenetic groups and infer that one group probably originated by a single-origin event from East China. Finally, we conduct genome-wide association analysis of nine important property-related traits to identify candidate genes, many of which are involved in cell wall, carbohydrate metabolism, and environmental adaptation. These results provide a foundation and resources for understanding moso bamboo evolution and the genetic mechanisms of agriculturally important traits.

Identification of genetic loci associated with major agronomic traits of wheat (Triticum aestivum L.) based on genome-wide association analysis

BACKGROUND

Bread wheat (Triticum aestivum L.) is one of the most widely consumed cereal crops, but its complex genome makes it difficult to investigate the genetic effect on important agronomic traits. Genome-wide association (GWA) analysis is a useful method to identify genetic loci controlling complex phenotypic traits. With the RNA-sequencing based gene expression analysis, putative candidate genes governing important agronomic trait can be suggested and also molecular markers can be developed.

RESULTS

We observed major quantitative agronomic traits of wheat; the winter survival rate (WSR), days to heading (DTH), days to maturity (DTM), stem length (SL), spike length (SPL), awn length (AL), liter weight (LW), thousand kernel weight (TKW), and the number of seeds per spike (SPS), of 287 wheat accessions from diverse country origins. A significant correlation was observed between the observed traits, and the wheat genotypes were divided into three subpopulations according to the population structure analysis. The best linear unbiased prediction (BLUP) values of the genotypic effect for each trait under different environments were predicted, and these were used for GWA analysis based on a mixed linear model (MLM). A total of 254 highly significant marker-trait associations (MTAs) were identified, and 28 candidate genes closely located to the significant markers were predicted by searching the wheat reference genome and RNAseq data. Further, it was shown that the phenotypic traits were significantly affected by the accumulation of favorable or unfavorable alleles.

CONCLUSIONS

From this study, newly identified MTA and putative agronomically useful genes will help to study molecular mechanism of each phenotypic trait. Further, the agronomically favorable alleles found in this study can be used to develop wheats with superior agronomic traits.

Genetic dissection of grain architecture-related traits in a winter wheat population

BACKGROUND

The future productivity of wheat (T. aestivum L.) as the most grown crop worldwide is of utmost importance for global food security. Thousand kernel weight (TKW) in wheat is closely associated with grain architecture-related traits, e.g. kernel length (KL), kernel width (KW), kernel area (KA), kernel diameter ratio (KDR), and factor form density (FFD). Discovering the genetic architecture of natural variation in these traits, identifying QTL and candidate genes are the main aims of this study. Therefore, grain architecture-related traits in 261 worldwide winter accessions over three field-year experiments were evaluated.

RESULTS

Genome-wide association analysis using 90K SNP array in FarmCPU model revealed several interesting genomic regions including 17 significant SNPs passing false discovery rate threshold and strongly associated with the studied traits. Four of associated SNPs were physically located inside candidate genes within LD interval e.g. BobWhite_c5872_589 (602,710,399 bp) found to be inside TraesCS6A01G383800 (602,699,767-602,711,726 bp). Further analysis reveals the four novel candidate genes potentially involved in more than one grain architecture-related traits with a pleiotropic effects e.g. TraesCS6A01G383800 gene on 6A encoding oxidoreductase activity was associated with TKW and KA. The allelic variation at the associated SNPs showed significant differences betweeen the accessions carying the wild and mutated alleles e.g. accessions carying C allele of BobWhite_c5872_589, TraesCS6A01G383800 had significantly higher TKW than the accessions carying T allele. Interestingly, these genes were highly expressed in the grain-tissues, demonstrating their pivotal role in controlling the grain architecture.

CONCLUSIONS

These results are valuable for identifying regions associated with kernel weight and dimensions and potentially help breeders in improving kernel weight and architecture-related traits in order to increase wheat yield potential and end-use quality.

Exploring natural diversity reveals alleles to enhance antioxidant system in barley under salt stress

Soil salinity stress causes osmotic/ionic imbalances and induces oxidative stress that causes cellular structure damage, perturbs metabolism, antioxidant system (comprising enzymatic and non-enzymatic components) and hence inhibits plant growth performance. In this study, we used genome-wide association scan (GWAS) in 174 diverse spring barley accessions which were exposed to salt stress under field conditions at the vegetative stage to uncover the genetic basis of antioxidant components and agronomic traits. High activities of enzymatic and content of non-enzymatic antioxidants were observed under salt stress compared to control conditions. Under salt stress, all the agronomic and yield-related traits were significantly reduced. Six genomic regions were associated with antioxidants and agronomic traits under salt stress conditions which were found to be linked with candidate genes. Several significant associations were physically located inside or near genes which are potentially involved in antioxidants. Two candidate genes at 2H (40,659,364 bp) and 7H (416,743,127 bp) were found to be involved in Dihydroflavonol 4-reductase/flavanone protein and Glyceraldehyde-3-phosphate dehydrogenase, respectively. The allelic variation at SNP of BK_07 at 7H inside the GAPDH gene demonstrates a negative selection of accessions carrying A allele. This allele appears in cultivars with lower activity of enzymatic antioxidants e.g. superoxide dismutase and catalases under salt stress conditions. These accessions are predominantly two-rowed, cultivars, originated from Europe, and carrying photoperiod sensitive alleles. The detected associated molecular markers in this work are considered as an important source for selection of increased amount of antioxidant compounds in barley under stress conditions.

Multi-locus genome-wide association mapping for spike-related traits in bread wheat (Triticum aestivum L.)

BACKGROUND

Bread wheat (Triticum aestivum L.) is one of the most important cereal food crops for the global population. Spike-layer uniformity (the consistency of the spike distribution in the vertical space)-related traits (SLURTs) are quantitative and have been shown to directly affect yield potential by modifying the plant architecture. Therefore, these parameters are important breeding targets for wheat improvement. The present study is the first genome-wide association study (GWAS) targeting SLURTs in wheat. In this study, a set of 225 diverse spring wheat accessions were used for multi-locus GWAS to evaluate SLURTs, including the number of spikes per plant (NSPP), spike length (SL), number of spikelets per spike (NSPS), grain weight per spike (GWPS), lowest tiller height (LTH), spike-layer thickness (SLT), spike-layer number (SLN) and spike-layer uniformity (SLU).

RESULTS

In total, 136 significant marker trait associations (MTAs) were identified when the analysis was both performed individually and combined for two environments. Twenty-nine MTAs were detected in environment one, 48 MTAs were discovered in environment two and 59 MTAs were detected using combined data from the two environments. Altogether, 15 significant MTAs were found for five traits in one of the two environments, and four significant MTAs were detected for the two traits, LTH and SLU, in both environments i.e. E1, E2 and also in combined data from the two environments. In total, 279 candidate genes (CGs) were identified, including Chaperone DnaJ, ABC transporter-like, AP2/ERF, SWEET sugar transporter, as well as genes that have previously been associated with wheat spike development, seed development and grain yield.

CONCLUSIONS

The MTAs detected through multi-locus GWAS will be useful for improving SLURTs and thus yield in wheat production through marker-assisted and genomic selection.

An Overview of Abiotic Stress in Cereal Crops: Negative Impacts, Regulation, Biotechnology and Integrated Omics

Abiotic stresses (AbS), such as drought, salinity, and thermal stresses, could highly affect the growth and development of plants. For decades, researchers have attempted to unravel the mechanisms of AbS for enhancing the corresponding tolerance of plants, especially for crop production in agriculture. In the present communication, we summarized the significant factors (atmosphere, soil and water) of AbS, their regulations, and integrated omics in the most important cereal crops in the world, especially rice, wheat, sorghum, and maize. It has been suggested that using systems biology and advanced sequencing approaches in genomics could help solve the AbS response in cereals. An emphasis was given to holistic approaches such as, bioinformatics and functional omics, gene mining and agronomic traits, genome-wide association studies (GWAS), and transcription factors (TFs) family with respect to AbS. In addition, the development of omics studies has improved to address the identification of AbS responsive genes and it enables the interaction between signaling pathways, molecular insights, novel traits and their significance in cereal crops. This review compares AbS mechanisms to omics and bioinformatics resources to provide a comprehensive view of the mechanisms. Moreover, further studies are needed to obtain the information from the integrated omics databases to understand the AbS mechanisms for the development of large spectrum AbS-tolerant crop production.

Genome-wide association mapping reveals key genomic regions for physiological and yield-related traits under salinity stress in wheat (Triticum aestivum L.)

A genome-wide association study (GWAS) was conducted using six different multi-locus GWAS models and 35K SNP array to demarcate genomic regions underlying reproductive stage salinity tolerance. Marker-trait association analysis was performed for salt tolerance indices (STI) of 11 morpho-physiological traits, and the actual concentrations of Na⁺ and K⁺, and the Na⁺/K⁺ ratio in flag leaf. A total of 293 significantly associated quantitative trait nucleotides (QTNs) for 14 morpho-physiological traits were identified. Of these 293 QTNs, 12 major QTNs with R² ≥ 10.0% were detected in three or more GWAS models. Novel major QTNs were identified for plant height, number of effective tillers, biomass, grain yield, thousand grain weight, Na⁺ and K⁺ content, and the Na⁺/K⁺ ratio in flag leaf. Moreover, 48 candidate genes were identified from the associated genomic regions. The QTNs identified in this study could potentially be targeted for improving salinity tolerance in wheat.

Optimizing sampling design and sequencing strategy for the genomic analysis of quantitative traits in natural populations

Mapping the genes underlying ecologically relevant traits in natural populations is fundamental to develop a molecular understanding of species adaptation. Current sequencing technologies enable the characterization of a species' genetic diversity across the landscape or even over its whole range. The relevant capture of the genetic diversity across the landscape is critical for a successful genetic mapping of traits and there are no clear guidelines on how to achieve an optimal sampling and which sequencing strategy to implement. Here we determine, through simulation, the sampling scheme that maximizes the power to map the genetic basis of a complex trait in an outbreeding species across an idealized landscape and draw genomic predictions for the trait, comparing individual and pool sequencing strategies. Our results show that quantitative trait locus detection power and prediction accuracy are higher when more populations over the landscape are sampled and this is more cost-effectively done with pool sequencing than with individual sequencing. Additionally, we recommend sampling populations from areas of high genetic diversity. As progress in sequencing enables the integration of trait-based functional ecology into landscape genomics studies, these findings will guide study designs allowing direct measures of genetic effects in natural populations across the environment.

Crop breeding for a changing climate: integrating phenomics and genomics with bioinformatics

Safeguarding crop yields in a changing climate requires bioinformatics advances in harnessing data from vast phenomics and genomics datasets to translate research findings into climate smart crops in the field. Climate change and an additional 3 billion mouths to feed by 2050 raise serious concerns over global food security. Crop breeding and land management strategies will need to evolve to maximize the utilization of finite resources in coming years. High-throughput phenotyping and genomics technologies are providing researchers with the information required to guide and inform the breeding of climate smart crops adapted to the environment. Bioinformatics has a fundamental role to play in integrating and exploiting this fast accumulating wealth of data, through association studies to detect genomic targets underlying key adaptive climate-resilient traits. These data provide tools for breeders to tailor crops to their environment and can be introduced using advanced selection or genome editing methods. To effectively translate research into the field, genomic and phenomic information will need to be integrated into comprehensive clade-specific databases and platforms alongside accessible tools that can be used by breeders to inform the selection of climate adaptive traits. Here we discuss the role of bioinformatics in extracting, analysing, integrating and managing genomic and phenomic data to improve climate resilience in crops, including current, emerging and potential approaches, applications and bottlenecks in the research and breeding pipeline.

Advances in Cereal Crop Genomics for Resilience under Climate Change

Adapting to climate change, providing sufficient human food and nutritional needs, and securing sufficient energy supplies will call for a radical transformation from the current conventional adaptation approaches to more broad-based and transformative alternatives. This entails diversifying the agricultural system and boosting productivity of major cereal crops through development of climate-resilient cultivars that can sustainably maintain higher yields under climate change conditions, expanding our focus to crop wild relatives, and better exploitation of underutilized crop species. This is facilitated by the recent developments in plant genomics, such as advances in genome sequencing, assembly, and annotation, as well as gene editing technologies, which have increased the availability of high-quality reference genomes for various model and non-model plant species. This has necessitated genomics-assisted breeding of crops, including underutilized species, consequently broadening genetic variation of the available germplasm; improving the discovery of novel alleles controlling important agronomic traits; and enhancing creation of new crop cultivars with improved tolerance to biotic and abiotic stresses and superior nutritive quality. Here, therefore, we summarize these recent developments in plant genomics and their application, with particular reference to cereal crops (including underutilized species). Particularly, we discuss genome sequencing approaches, quantitative trait loci (QTL) mapping and genome-wide association (GWAS) studies, directed mutagenesis, plant non-coding RNAs, precise gene editing technologies such as CRISPR-Cas9, and complementation of crop genotyping by crop phenotyping. We then conclude by providing an outlook that, as we step into the future, high-throughput phenotyping, pan-genomics, transposable elements analysis, and machine learning hold much promise for crop improvements related to climate resilience and nutritional superiority.

Genome-wide association mapping in maize: status and prospects

Genome-wide association study (GWAS) provides a robust and potent tool to retrieve complex phenotypic traits back to their underlying genetics. Maize is an excellent crop for performing GWAS due to diverse genetic variability, rapid decay of linkage disequilibrium, availability of distinct sub-populations and abundant SNP information. The application of GWAS in maize has resulted in successful identification of thousands of genomic regions associated with many abiotic and biotic stresses. Many agronomic and quality traits of maize are severely affected by such stresses and, significantly affecting its growth and productivity. To improve productivity of maize crop in countries like India which contribute only 2% to the world's total production in 2019-2020, it is essential to understand genetic complexity of underlying traits. Various DNA markers and trait associations have been revealed using conventional linkage mapping methods. However, it has achieved limited success in improving polygenic complex traits due to lower resolution of trait mapping. The present review explores the prospects of GWAS in improving yield, quality and stress tolerance in maize besides, strengths and challenges of using GWAS for molecular breeding and genomic selection. The information gathered will facilitate elucidation of genetic mechanisms of complex traits and improve efficiency of marker-assisted selection in maize breeding.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02799-4.

QTL Analysis of Adult Plant Resistance to Stripe Rust in a Winter Wheat Recombinant Inbred Population

Stripe rust, caused by the fungus Puccinia striiformis f. sp. tritici, is a worldwide disease of wheat that causes devastating crop losses. Resistant cultivars have been developed over the last 40 years that have significantly reduced the economic impact of the disease on growers, but in heavy infection years it is mostly controlled through the intensive application of fungicides. The Pacific Northwest of the United States has an ideal climate for stripe rust and has one of the most diverse race compositions in the country. This has resulted in many waves of epidemics that have overcome most of the resistance genes traditionally used in elite germplasm. The best way to prevent high yield losses, reduce production costs to growers, and reduce the heavy application of fungicides is to pyramid multiple stripe rust resistance genes into new cultivars. Using genotyping-by-sequencing, we identified 4662 high quality variant positions in a recombinant inbred line population of 196 individuals derived from a cross between Skiles, a highly resistant winter wheat cultivar, and Goetze, a moderately to highly susceptible winter wheat cultivar, both developed at Oregon State University. A subsequent genome wide association study identified two quantitative trait loci (QTL) on chromosomes 3B and 3D within the predicted locations of stripe rust resistance genes. Resistance QTL, when combined together, conferred high levels of stripe rust resistance above the level of Skiles in some locations, indicating that these QTL would be important additions to future breeding efforts of Pacific Northwest winter wheat cultivars.

Genomic regions associated with stripe rust resistance against the Egyptian race revealed by genome-wide association study

BACKGROUND

Wheat stripe rust (caused by Puccinia striiformis f. sp. Tritici), is a major disease that causes huge yield damage. New pathogen races appeared in the last few years and caused a broke down in the resistant genotypes. In Egypt, some of the resistant genotypes began to be susceptible to stripe rust in recent years. This situation increases the need to produce new genotypes with durable resistance. Besides, looking for a new resistant source from the available wheat genotypes all over the world help in enhancing the breeding programs.

RESULTS

In the recent study, a set of 103-spring wheat genotypes from different fourteen countries were evaluated to their field resistant to stripe rust for two years. These genotypes included 17 Egyptian genotypes from the old and new cultivars. The 103-spring wheat genotypes were reported to be well adapted to the Egyptian environmental conditions. Out of the tested genotypes, eight genotypes from four different countries were found to be resistant in both years. Genotyping was carried out using genotyping-by-sequencing and a set of 26,703 SNPs were used in the genome-wide association study. Five SNP markers, located on chromosomes 2A and 4A, were found to be significantly associated with the resistance in both years. Three gene models associated with disease resistance and underlying these significant SNPs were identified. One immune Iranian genotype, with the highest number of different alleles from the most resistant Egyptian genotypes, was detected.

CONCLUSION

the high variation among the tested genotypes in their resistance to the Egyptian stripe rust race confirming the possible improvement of stripe rust resistance in the Egyptian wheat genotypes. The identified five SNP markers are stable and could be used in marker-assisted selection after validation in different genetic backgrounds. Crossing between the immune Iranian genotype and the Egyptian genotypes will improve stripe rust resistance in Egypt.