Combined Linkage and Association Mapping Reveals QTL and Candidate Genes for Plant and Ear Height in Maize

A combination of QTL mapping and genome-wide association study revealed the key gene for husk number in maize

KEY MESSAGE

Two key genes Zm00001d021232 and Zm00001d048138 were identified by QTL mapping and GWAS. Additionally, they were verified to be significantly associated with maize husk number (HN) using gene-based association study. As a by-product of maize production, maize husk is an important industrial raw material. Husk layer number (HN) is an important trait that affects the yield of maize husk. However, the genetic mechanism underlying HN remains unclear. Herein, a total of 13 quantitative trait loci (QTL) controlling HN were identified in an IBM Syn 10 DH population across different locations. Among these, three QTL were individually repeatedly detected in at least two environments. Meanwhile, 26 unique single nucleotide polymorphisms (SNPs) were detected to be significantly (p < 2.15 × 10-6) associated with HN in an association pool. Of these SNPs, three were simultaneously detected across multiple environments or environments and best linear unbiased prediction (BLUP). We focused on these environment-stable and population-common genetic loci for excavating the candidate genes responsible for maize HN. Finally, 173 initial candidate genes were identified, of which 22 were involved in both multicellular organism development and single-multicellular organism process and thus confirmed as the candidate genes for HN. Gene-based association analyses revealed that the variants in four genes were significantly (p < 0.01/N) correlated with HN, of which Zm00001d021232 and Zm00001d048138 were highly expressed in husks and early developing ears among different maize tissues. Our study contributes to the understanding of genetic and molecular mechanisms of maize husk yield and industrial development in the future.

QTL mapping for plant height and ear height using bi-parental immortalized heterozygous populations in maize

Introduction

Plant height (PH) and ear height (EH) are key plant architectural traits in maize, which will affect the photosynthetic efficiency, high plant density tolerance, suitability for mechanical harvesting.

Methods

QTL mapping were conducted for PH and EH using a recombinant inbred line (RIL) population and two corresponding immortalized backcross (IB) populations obtained from crosses between the RIL population and the two parental lines.

Results

A total of 17 and 15 QTL were detected in the RIL and IB populations, respectively. Two QTL, qPH1-1 (qEH1-1) and qPH1-2 (qEH1-4) in the RIL, were simultaneously identified for PH and EH. Combing reported genome-wide association and cloned PH-related genes, co-expression network analyses were constructed, then five candidate genes with high confidence in major QTL were identified including Zm00001d011117 and Zm00001d011108, whose homologs have been confirmed to play a role in determining PH in maize and soybean.

Discussion

QTL mapping used a immortalized backcross population is a new strategy. These identified genes in this study can provide new insights for improving the plant architecture in maize.

Ameliorating the effects of multiple stresses on agronomic traits in crops: modern biotechnological and omics approaches

While global climate change poses a significant environmental threat to agriculture, the increasing population is another big challenge to food security. To address this, developing crop varieties with increased productivity and tolerance to biotic and abiotic stresses is crucial. Breeders must identify traits to ensure higher and consistent yields under inconsistent environmental challenges, possess resilience against emerging biotic and abiotic stresses and satisfy customer demands for safer and more nutritious meals. With the advent of omics-based technologies, molecular tools are now integrated with breeding to understand the molecular genetics of genotype-based traits and develop better climate-smart crops. The rapid development of omics technologies offers an opportunity to generate novel datasets for crop species. Identifying genes and pathways responsible for significant agronomic traits has been made possible by integrating omics data with genetic and phenotypic information. This paper discusses the importance and use of omics-based strategies, including genomics, transcriptomics, proteomics and phenomics, for agricultural and horticultural crop improvement, which aligns with developing better adaptability in these crop species to the changing climate conditions.

Investigating the genetic basis of maize ear characteristics: a comprehensive genome-wide study utilizing high-throughput phenotypic measurement method and system

The morphology of maize ears plays a critical role in the breeding of new varieties and increasing yield. However, the study of traditional ear-related traits alone can no longer meet the requirements of breeding. In this study, 20 ear-related traits, including size, shape, number, and color, were obtained in 407 maize inbred lines at two sites using a high-throughput phenotypic measurement method and system. Significant correlations were found among these traits, particularly the novel trait ear shape (ES), which was correlated with traditional traits: kernel number per row and kernel number per ear. Pairwise comparison tests revealed that the inbred lines of tropical-subtropical were significantly different from other subpopulations in row numbers per ear, kernel numbers per ear, and ear color. A genome-wide association study identified 275, 434, and 362 Single nucleotide polymorphisms (SNPs) for Beijing, Sanya, and best linear unbiased prediction scenarios, respectively, explaining 3.78% to 24.17% of the phenotypic variance. Furthermore, 58 candidate genes with detailed functional descriptions common to more than two scenarios were discovered, with 40 genes being associated with color traits on chromosome 1. After analysis of haplotypes, gene expression, and annotated information, several candidate genes with high reliability were identified, including Zm00001d051328 for ear perimeter and width, zma-MIR159f for ear shape, Zm00001d053080 for kernel width and row number per ear, and Zm00001d048373 for the blue color channel of maize kernels in the red-green-blue color model. This study emphasizes the importance of researching novel phenotypic traits in maize by utilizing high-throughput phenotypic measurements. The identified genetic loci enrich the existing genetic studies related to maize ears.

Silage maize as a potent candidate for sustainable animal husbandry development-perspectives and strategies for genetic enhancement

Maize is recognized as the queen of cereals, with an ability to adapt to diverse agroecologies (from 58^oN to 55^oS latitude) and the highest genetic yield potential among cereals. Under contemporary conditions of global climate change, C₄ maize crops offer resilience and sustainability to ensure food, nutritional security, and farmer livelihood. In the northwestern plains of India, maize is an important alternative to paddy for crop diversification in the wake of depleting water resources, reduced farm diversity, nutrient mining, and environmental pollution due to paddy straw burning. Owing to its quick growth, high biomass, good palatability, and absence of anti-nutritional components, maize is also one of the most nutritious non-legume green fodders. It is a high-energy, low-protein forage commonly used for dairy animals like cows and buffalos, often in combination with a complementary high-protein forage such as alfalfa. Maize is also preferred for silage over other fodders due to its softness, high starch content, and sufficient soluble sugars required for proper ensiling. With a rapid population increase in developing countries like China and India, there is an upsurge in meat consumption and, hence, the requirement for animal feed, which entails high usage of maize. The global maize silage market is projected to grow at a compound annual growth rate of 7.84% from 2021 to 2030. Factors such as increasing demand for sustainable and environment-friendly food sources coupled with rising health awareness are fueling this growth. With the dairy sector growing at about 4%-5% and the increasing shortage faced for fodder, demand for silage maize is expected to increase worldwide. The progress in improved mechanization for the provision of silage maize, reduced labor demand, lack of moisture-related marketing issues as associated with grain maize, early vacancy of farms for next crops, and easy and economical form of feed to sustain household dairy sector make maize silage a profitable venture. However, sustaining the profitability of this enterprise requires the development of hybrids specific for silage production. Little attention has yet been paid to breeding for a plant ideotype for silage with specific consideration of traits such as dry matter yield, nutrient yield, energy in organic matter, genetic architecture of cell wall components determining their digestibility, stalk standability, maturity span, and losses during ensiling. This review explores the available information on the underlying genetic mechanisms and gene/gene families impacting silage yield and quality. The trade-offs between yield and nutritive value in relation to crop duration are also discussed. Based on available genetic information on inheritance and molecular aspects, breeding strategies are proposed to develop maize ideotypes for silage for the development of sustainable animal husbandry.

Mapping of QTLs for morphophysiological and yield traits under water-deficit stress and well-watered conditions in maize

Maize productivity is significantly impacted by drought; therefore, improvement of drought tolerance is a critical goal in maize breeding. To achieve this, a better understanding of the genetic basis of drought tolerance is necessary. Our study aimed to identify genomic regions associated with drought tolerance-related traits by phenotyping a mapping population of recombinant inbred lines (RILs) for two seasons under well-watered (WW) and water-deficit (WD) conditions. We also used single nucleotide polymorphism (SNP) genotyping through genotyping-by-sequencing to map these regions and attempted to identify candidate genes responsible for the observed phenotypic variation. Phenotyping of the RILs population revealed significant variability in most of the traits, with normal frequency distributions, indicating their polygenic nature. We generated a linkage map using 1,241 polymorphic SNPs distributed over 10 chromosomes (chrs), covering a total genetic distance of 5,471.55 cM. We identified 27 quantitative trait loci (QTLs) associated with various morphophysiological and yield-related traits, with 13 QTLs identified under WW conditions and 12 under WD conditions. We found one common major QTL (qCW2-1) for cob weight and a minor QTL (qCH1-1) for cob height that were consistently identified under both water regimes. We also detected one major and one minor QTL for the Normalized Difference Vegetation Index (NDVI) trait under WD conditions on chr 2, bin 2.10. Furthermore, we identified one major QTL (qCH1-2) and one minor QTL (qCH1-1) on chr 1 that were located at different genomic positions to those identified in earlier studies. We found co-localized QTLs for stomatal conductance and grain yield on chr 6 (qg_s6-2 and qGY6-1), while co-localized QTLs for stomatal conductance and transpiration rate were identified on chr 7 (qg_s7-1 and qTR7-1). We also attempted to identify the candidate genes responsible for the observed phenotypic variation; our analysis revealed that the major candidate genes associated with QTLs detected under water deficit conditions were related to growth and development, senescence, abscisic acid (ABA) signaling, signal transduction, and transporter activity in stress tolerance. The QTL regions identified in this study may be useful in designing markers that can be utilized in marker-assisted selection breeding. In addition, the putative candidate genes can be isolated and functionally characterized so that their role in imparting drought tolerance can be more fully understood.

Identification of QTLs and Putative Candidate Genes for Plant Architecture of Lotus Revealed by Regional Association Mapping

The lotus (Nelumbo Adans.) is one of the most economically relevant ornamental aquatic plants. Plant architecture (PA) is an important trait for lotus classification, cultivation, breeding, and applications. However, the underlying genetic and molecular basis controlling PA remains poorly understood. In this study, an association study for PA-related traits was performed with 93 genome-wide microsatellite markers (simple sequence repeat, SSR) and 51 insertion-deletion (InDel) markers derived from the candidate regions using a panel of 293 lotus accessions. Phenotypic data analysis of the five PA-related traits revealed a wide normal distribution and high heritability from 2013 to 2016, which indicated that lotus PA-related traits are highly polygenic traits. The population structure (Q-matrix) and the relative kinships (K-matrix) of the association panels were analyzed using 93 SSR markers. The mixed linear model (MLM) taking Q-matrix and K-matrix into account was used to estimate the association between markers and the traits. A total of 26 markers and 65 marker-trait associations were identified by considering associations with p < 0.001 and Q < 0.05. Based on the significant markers, two QTLs on Chromosome 1 were identified, and two candidate genes were preliminarily determined. The results of our study provided useful information for the lotus breeding aiming at different PA phenotypes using a molecular-assisted selection (MAS) method and also laid the foundation for the illustration of the molecular mechanism underlying the major QTL and key markers associated with lotus PA.

Identification of SNP loci and candidate genes genetically controlling norisoprenoids in grape berry based on genome-wide association study

Obtaining new grapevine varieties with unique aromas has been a long-standing goal of breeders. Norisoprenoids are of particular interest to wine producers and researchers, as these compounds are responsible for the important varietal aromas in wine, characterized by a complex floral and fruity smell, and are likely present in all grape varieties. However, the single-nucleotide polymorphism (SNP) loci and candidate genes genetically controlling the norisoprenoid content in grape berry remain unknown. To this end, in this study, we investigated 13 norisoprenoid traits across two years in an F₁ population consisting of 149 individuals from a hybrid of Vitis vinifera L. cv. Muscat Alexandria and V. vinifera L. cv. Christmas Rose. Based on 568,953 SNP markers, genome-wide association analysis revealed that 27 candidate SNP loci belonging to 18 genes were significantly associated with the concentrations of norisoprenoid components in grape berry. Among them, 13 SNPs were confirmed in a grapevine germplasm population comprising 97 varieties, including two non-synonymous mutations SNPs within the VvDXS1 and VvGGPPS genes, respectively in the isoprenoid metabolic pathway. Genotype analysis showed that the grapevine individuals with the heterozygous genotype C/T at chr5:2987350 of VvGGPPS accumulated higher average levels of 6-methyl-5-hepten-2-one and β-cyclocitral than those with the homozygous genotype C/C. Furthermore, VvGGPPS was highly expressed in individuals with high norisoprenoids concentrations. Transient overexpression of VvGGPPS in the leaves of Vitis quinquangularis and tobacco resulted in an increase in norisoprenoid concentrations. These findings indicate the importance of VvGGPPS in the genetic control of norisoprenoids in grape berries, serving as a potential molecular breeding target for aroma.

Genome-Wide Association Study and Genomic Prediction on Plant Architecture Traits in Sweet Corn and Waxy Corn

Sweet corn and waxy corn has a better taste and higher accumulated nutritional value than regular maize, and is widely planted and popularly consumed throughout the world. Plant height (PH), ear height (EH), and tassel branch number (TBN) are key plant architecture traits, which play an important role in improving grain yield in maize. In this study, a genome-wide association study (GWAS) and genomic prediction analysis were conducted on plant architecture traits of PH, EH, and TBN in a fresh edible maize population consisting of 190 sweet corn inbred lines and 287 waxy corn inbred lines. Phenotypic data from two locations showed high heritability for all three traits, with significant differences observed between sweet corn and waxy corn for both PH and EH. The differences between the three subgroups of sweet corn were not obvious for all three traits. Population structure and PCA analysis results divided the whole population into three subgroups, i.e., sweet corn, waxy corn, and the subgroup mixed with sweet and waxy corn. Analysis of GWAS was conducted with 278,592 SNPs obtained from resequencing data; 184, 45, and 68 significantly associated SNPs were detected for PH, EH, and TBN, respectively. The phenotypic variance explained (PVE) values of these significant SNPs ranged from 3.50% to 7.0%. The results of this study lay the foundation for further understanding the genetic basis of plant architecture traits in sweet corn and waxy corn. Genomic selection (GS) is a new approach for improving quantitative traits in large plant breeding populations that uses whole-genome molecular markers. The marker number and marker quality are essential for the application of GS in maize breeding. GWAS can choose the most related markers with the traits, so it can be used to improve the predictive accuracy of GS.

Characterization and genetic dissection of maize ear leaf midrib acquired by 3D digital technology

The spatial morphological structure of plant leaves is an important index to evaluate crop ideotype. In this study, we characterized the three-dimensional (3D) data of the ear leaf midrib of maize at the grain-filling stage using the 3D digitization technology and obtained the phenotypic values of 15 traits covering four different dimensions of the ear leaf midrib, of which 13 phenotypic traits were firstly proposed for featuring plant leaf spatial structure. Cluster analysis results showed that the 13 traits could be divided into four groups, Group I, -II, -III and -IV. Group I contains HorizontalLength, OutwardGrowthMeasure, LeafAngle and DeviationTip; Group II contains DeviationAngle, MaxCurvature and CurvaturePos; Group III contains LeafLength and ProjectionArea; Group IV contains TipTop, VerticalHeight, UpwardGrowthMeasure, and CurvatureRatio. To investigate the genetic basis of the ear leaf midrib curve, 13 traits with high repeatability were subjected to genome-wide association study (GWAS) analysis. A total of 828 significantly related SNPs were identified and 1365 candidate genes were annotated. Among these, 29 candidate genes with the highest significant and multi-method validation were regarded as the key findings. In addition, pathway enrichment analysis was performed on the candidate genes of traits to explore the potential genetic mechanism of leaf midrib curve phenotype formation. These results not only contribute to further understanding of maize leaf spatial structure traits but also provide new genetic loci for maize leaf spatial structure to improve the plant type of maize varieties.

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

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

GWAS and RNA-seq analysis uncover candidate genes associated with alkaline stress tolerance in maize (Zea mays L.) seedlings

Soil salt-alkalization is a common yet critical environmental stress factor for plant growth and development. Discovering and exploiting genes associated with alkaline tolerance in maize (Zea mays L.) is helpful for improving alkaline resistance. Here, an association panel consisting of 200 maize lines was used to identify the genetic loci responsible for alkaline tolerance-related traits in maize seedlings. A total of nine single-nucleotide polymorphisms (SNPs) and their associated candidate genes were found to be significantly associated with alkaline tolerance using a genome-wide association study (GWAS). An additional 200 genes were identified when the screen was extended to include a linkage disequilibrium (LD) decay distance of r² ≥ 0.2 from the SNPs. RNA-sequencing (RNA-seq) analysis was then conducted to confirm the linkage between the candidate genes and alkali tolerance. From these data, a total of five differentially expressed genes (DEGs; |log2FC| ≥ 0.585, p < 0.05) were verified as the hub genes involved in alkaline tolerance. Subsequently, two candidate genes, Zm00001d038250 and Zm00001d001960, were verified to affect the alkaline tolerance of maize seedlings by qRT-PCR analysis. These genes were putatively involved protein binding and "flavonoid biosynthesis process," respectively, based on Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses. Gene promoter region contains elements related to stress and metabolism. The results of this study will help further elucidate the mechanisms of alkaline tolerance in maize, which will provide the groundwork for future breeding projects.

The Use of DArTseq Technology to Identify New SNP and SilicoDArT Markers Related to the Yield-Related Traits Components in Maize

In the last decade, many scientists have used molecular biology methods in their research to locate the grain-yield-determining loci and yield structure characteristics in maize. Large-scale molecular analyses in maize do not only focus on the identification of new markers and quantitative trait locus (QTL) regions. DNA analysis in the selection of parental components for heterotic crosses is a very important tool for breeders. The aim of this research was to identify and select new markers for maize (SNP and SilicoDArT) linked to genes influencing the size of the yield components in maize. The plant material used for the research was 186 inbred maize lines. The field experiment was established in twolocations. The yield and six yield components were analyzed. For identification of SNP and SilicoDArT markers related to the yield and yield components, next-generation sequencing was used. As a result of the biometric measurements analysis, differentiation in the average elevation of the analyzed traits for the lines in both locations was found. The above-mentioned results indicate the existence of genotype–environment interactions. The analysis of variance for the observed quality between genotypes indicated a statistically significant differentiation between genotypes and a statistically significant differentiation for all the observed properties betweenlocations. A canonical variable analysis was applied to present a multi-trait assessment of the similarity of the tested maize genotypes in a lower number of dimensions with the lowest possible loss of information. No grouping of lines due to the analyzed was observed. As a result of next-generation sequencing, the molecular markers SilicoDArT (53,031) and SNP (28,571) were obtained. The genetic distance between the analyzed lines was estimated on the basis of these markers. Out of 81,602 identified SilicoDArT and SNP markers, 15,409 (1559 SilicoDArT and 13,850 SNPs) significantly related to the analyzed yield components were selected as a result of association mapping. The greatest numbers of molecular markers were associated with cob length (1203), cob diameter (1759), core length (1201) and core diameter (2326). From 15,409 markers significantly related to the analyzed traits of the yield components, 18 DArT markers were selected, which were significant for the same four traits (cob length, cob diameter, core length, core diameter) in both Kobierzyce and Smolice. These markers were used for physical mapping. As a result of the analyses, it was found that 6 out of 18 (1818; 14,506; 2317; 3233; 11,657; 12,812) identified markers are located inside genes. These markers are located on chromosomes 8, 9, 7, 3, 5, and 1, respectively.

Association mapping for general combining ability with yield, plant height and ear height using F1 population in maize

General combining ability (GCA) is an important index for inbred lines breeding of maize. To identify the genetic loci of GCA and associated agronomic traits, an association analysis with 195 SSRs was made in phenotypic traits of 240 F1 derived from 120 elite inbred lines containing current breeding resources of maize crossed with 2 testers (Zheng58 and Chang7-2) in two places in 2018. All of the 20 association loci detected for grain yield (GY), plant height (PH), ear height (EH) and GCA for the three traits in two places could explain a phenotypic variation range of 7.31%-9.29%. Among the 20 association loci, 9 (7.31%-9.04%) were associated with GY, 4 (7.22%-8.91%) were related to GCA of GY, 1 (7.56%) was associated with PH, and 3 (7.53%-8.96%) were related to EH. In addition, 3 loci (9.14%-9.29%) were associated with GCA of PH whereas no locus was identified for GCA of EH. In the comparison of the association loci detected in Baoding and Handan, interestingly, one locus (7.69% and 8.11%) was identified in both environments and one locus (7.52% and 7.82%) was identified for yield and GCA of yield. Therefore, the identification of GY-, PH-, EH- and GCA-related association loci could not only provide references for high yield breeding of maize, but also help us comprehend the relationships among GY, agricultural traits and GCA.

Omics for the Improvement of Abiotic, Biotic, and Agronomic Traits in Major Cereal Crops: Applications, Challenges, and Prospects

Omics technologies, namely genomics, transcriptomics, proteomics, metabolomics, and phenomics, are becoming an integral part of virtually every commercial cereal crop breeding program, as they provide substantial dividends per unit time in both pre-breeding and breeding phases. Continuous advances in omics assure time efficiency and cost benefits to improve cereal crops. This review provides a comprehensive overview of the established omics methods in five major cereals, namely rice, sorghum, maize, barley, and bread wheat. We cover the evolution of technologies in each omics section independently and concentrate on their use to improve economically important agronomic as well as biotic and abiotic stress-related traits. Advancements in the (1) identification, mapping, and sequencing of molecular/structural variants; (2) high-density transcriptomics data to study gene expression patterns; (3) global and targeted proteome profiling to study protein structure and interaction; (4) metabolomic profiling to quantify organ-level, small-density metabolites, and their composition; and (5) high-resolution, high-throughput, image-based phenomics approaches are surveyed in this review.

Genome-wide association study reveals the genetic architecture of root hair length in maize

Background

Root hair, a special type of tubular-shaped cell, outgrows from root epidermal cell and plays important roles in the acquisition of nutrients and water, as well as interactions with biotic and abiotic stress. Although many genes involved in root hair development have been identified, genetic basis of natural variation in root hair growth has never been explored.

Results

Here, we utilized a maize association panel including 281 inbred lines with tropical, subtropical, and temperate origins to decipher the phenotypic diversity and genetic basis of root hair length. We demonstrated significant associations of root hair length with many metabolic pathways and other agronomic traits. Combining root hair phenotypes with 1.25 million single nucleotide polymorphisms (SNPs) via genome-wide association study (GWAS) revealed several candidate genes implicated in cellular signaling, polar growth, disease resistance and various metabolic pathways.

Conclusions

These results illustrate the genetic basis of root hair length in maize, offering a list of candidate genes predictably contributing to root hair growth, which are invaluable resource for the future functional investigation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07961-z.

Dissecting the Genetic Basis of Flowering Time and Height Related-Traits Using Two Doubled Haploid Populations in Maize

In the field, maize flowering time and height traits are closely linked with yield, planting density, lodging resistance, and grain fill. To explore the genetic basis of flowering time and height traits in maize, we investigated six related traits, namely, days to anthesis (AD), days to silking (SD), the anthesis-silking interval (ASI), plant height (PH), ear height (EH), and the EH/PH ratio (ER) in two locations for two years based on two doubled haploid (DH) populations. Based on the two high-density genetic linkage maps, 12 and 22 quantitative trait loci (QTL) were identified, respectively, for flowering time and height-related traits. Of these, ten QTLs had overlapping confidence intervals between the two populations and were integrated into three consensus QTLs (qFT_YZ1a, qHT_YZ5a, and qHT_YZ7a). Of these, qFT_YZ1a, conferring flowering time, is located at 221.1-277.0 Mb on chromosome 1 and explained 10.0-12.5% of the AD and SD variation, and qHT_YZ5a, conferring height traits, is located at 147.4-217.3 Mb on chromosome 5 and explained 11.6-15.3% of the PH and EH variation. These consensus QTLs, in addition to the other repeatedly detected QTLs, provide useful information for further genetic studies and variety improvements in flowering time and height-related traits.

Using chlorate as an analogue to nitrate to identify candidate genes for nitrogen use efficiency in barley

Nitrogen (N) is one of the most important macronutrients for crop growth and development. Large amounts of N fertilizers are applied exogenously to improve grain yield and quality, which has led to environmental pollution and high cost of production. Therefore, improvement of N use efficiency (NUE) is a very important aspect for sustainable agriculture. Here, a pilot experiment was firstly conducted with a set of barley genotypes with confirmed NUE to validate the fast NUE screening, using chlorate as an analogue to nitrate. High NUE genotypes were susceptible to chlorate-induced toxicity whereas the low NUE genotypes were tolerant. A total of 180 barley RILs derived from four parents (Compass, GrangeR, Lockyer and La Trobe) were further screened for NUE. Leaf chlorosis induced by chlorate toxicity was the key parameter observed which was later related to low-N tolerance of the RILs. There was a distinct variation in chlorate susceptibility of the RILs with leaf chlorosis in the oldest leaf ranging from 10 to 80%. A genome-wide association study (GWAS) identified 9 significant marker-trait associations (MTAs) conferring high chlorate sensitivity on chromosomes 2H (2), 3H (1), 4H (4), 5H (1) and Un (1). Genes flanking with these markers were retrieved as potential targets for genetic improvement of NUE. Genes encoding Ferredoxin 3, leucine-rich receptor-like protein kinase family protein and receptor kinase are responsive to N stress. MTA4H5468 which exhibits concordance with high NUE phenotype can further be explored under different genetic backgrounds and successfully applied in marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01239-8.

Different loci control resistance to different isolates of the same race of Colletotrichum lindemuthianum in common bean

Linkage and genome-wide association analyses using high-throughput SNP genotyping revealed different loci controlling resistance to different isolates of race 65 of Colletotrichum lindemuthianum in common bean. Development of varieties with durable resistance to anthracnose is a major challenge in common bean breeding programs because of the extensive virulence diversity of Colletotrichum lindemuthianum fungus. We used linkage and genome-wide association analyses to tap the genomic regions associated with resistance to different isolates of race 65. Linkage mapping was done using an F₂ population derived from the cross between the Mesoamerican common beans BRS Estilo x Ouro Vermelho, inoculated with two different isolates of race 65. Association genetics relied on a diversity common bean panel containing 189 common bean accessions inoculated with five different isolates of race 65 as an attempt to validate the linkage analysis findings and, eventually, identify other genomic regions associated with resistance to race 65. The F₂ population and diversity panel were genotyped with the BARCBean6K_3 Illumina BeadChip containing 5398 SNP markers. Both linkage and genome-wide association analyses identified different loci controlling resistance to different isolates of race 65 on linkage group Pv04. Genome-wide association analysis also detected loci on Pv05, Pv10 and Pv11 associated with resistance to race 65. These findings indicate that resistance to race 65 can be overcome by the virulence diversity among different isolates of the same race and could lead to the loss of resistance after cultivar release. We identified 25 resistant common bean cultivars to all five isolates of race 65 in the diversity panel. The accessions should be useful to develop cultivars combining different resistance genes that favor durable resistance to anthracnose in common bean.

Genetic Dissection of Quantitative Resistance to Common Rust (Puccinia sorghi) in Tropical Maize (Zea mays L.) by Combined Genome-Wide Association Study, Linkage Mapping, and Genomic Prediction

Common rust is one of the major foliar diseases in maize, leading to significant grain yield losses and poor grain quality. To dissect the genetic architecture of common rust resistance, a genome-wide association study (GWAS) panel and a bi-parental doubled haploid (DH) population, DH1, were used to perform GWAS and linkage mapping analyses. The GWAS results revealed six single-nucleotide polymorphisms (SNPs) significantly associated with quantitative resistance of common rust at a very stringent threshold of P-value 3.70 × 10^-6 at bins 1.05, 1.10, 3.04, 3.05, 4.08, and 10.04. Linkage mapping identified five quantitative trait loci (QTL) at bins 1.03, 2.06, 4.08, 7.03, and 9.00. The phenotypic variation explained (PVE) value of each QTL ranged from 5.40 to 12.45%, accounting for the total PVE value of 40.67%. Joint GWAS and linkage mapping analyses identified a stable genomic region located at bin 4.08. Five significant SNPs were only identified by GWAS, and four QTL were only detected by linkage mapping. The significantly associated SNP of S10_95231291 detected in the GWAS analysis was first reported. The linkage mapping analysis detected two new QTL on chromosomes 7 and 10. The major QTL on chromosome 7 in the region between 144,567,253 and 149,717,562 bp had the largest PVE value of 12.45%. Four candidate genes of GRMZM2G328500, GRMZM2G162250, GRMZM2G114893, and GRMZM2G138949 were identified, which played important roles in the response of stress resilience and the regulation of plant growth and development. Genomic prediction (GP) accuracies observed in the GWAS panel and DH1 population were 0.61 and 0.51, respectively. This study provided new insight into the genetic architecture of quantitative resistance of common rust. In tropical maize, common rust could be improved by pyramiding the new sources of quantitative resistance through marker-assisted selection (MAS) or genomic selection (GS), rather than the implementation of MAS for the single dominant race-specific resistance gene.

Dissecting the phenotypic components and genetic architecture of maize stem vascular bundles using high-throughput phenotypic analysis

High-throughput phenotyping is increasingly becoming an important tool for rapid advancement of genetic gain in breeding programmes. Manual phenotyping of vascular bundles is tedious and time-consuming, which lags behind the rapid development of functional genomics in maize. More robust and automated techniques of phenotyping vascular bundles traits at high-throughput are urgently needed for large crop populations. In this study, we developed a standard process for stem micro-CT data acquisition and an automatic CT image process pipeline to obtain vascular bundle traits of stems including geometry-related, morphology-related and distribution-related traits. Next, we analysed the phenotypic variation of stem vascular bundles between natural population subgroup (480 inbred lines) based on 48 comprehensively phenotypic information. Also, the first database for stem micro-phenotypes, MaizeSPD, was established, storing 554 pieces of basic information of maize inbred lines, 523 pieces of experimental information, 1008 pieces of CT scanning images and processed images, and 24 192 pieces of phenotypic data. Combined with genome-wide association studies (GWASs), a total of 1562 significant single nucleotide polymorphism (SNPs) were identified for 30 stem micro-phenotypic traits, and 84 unique genes of 20 traits such as VBNum, VBAvArea and PZVBDensity were detected. Candidate genes identified by GWAS mainly encode enzymes involved in cell wall metabolism, transcription factors, protein kinase and protein related to plant signal transduction and stress response. The results presented here will advance our knowledge about phenotypic trait components of stem vascular bundles and provide useful information for understanding the genetic controls of vascular bundle formation and development.

Genome-Wide Association Study and Identification of Candidate Genes for Nitrogen Use Efficiency in Barley (Hordeum vulgare L.)

Nitrogen (N) fertilizer is largely responsible for barley grain yield potential and quality, yet excessive application leads to environmental pollution and high production costs. Therefore, efficient use of N is fundamental for sustainable agriculture. In the present study, we investigated the performance of 282 barley accessions through hydroponic screening using optimal and low NH4NO3 treatments. Low-N treatment led to an average shoot dry weight reduction of 50%, but there were significant genotypic differences among the accessions. Approximately 20% of the genotypes showed high (>75%) relative shoot dry weight under low-N treatment and were classified as low-N tolerant, whereas 20% were low-N sensitive (≤55%). Low-N tolerant accessions exhibited well-developed root systems with an average increase of 60% in relative root dry weight to facilitate more N absorption. A genome-wide association study (GWAS) identified 66 significant marker trait associations (MTAs) conferring high nitrogen use efficiency, four of which were stable across experiments. These four MTAs were located on chromosomes 1H(1), 3H(1), and 7H(2) and were associated with relative shoot length, relative shoot and root dry weight. Genes corresponding to the significant MTAs were retrieved as candidate genes, including members of the asparagine synthetase gene family, several transcription factor families, protein kinases, and nitrate transporters. Most importantly, the high-affinity nitrate transporter 2.7 (HvNRT2.7) was identified as a promising candidate on 7H for root and shoot dry weight. The identified candidate genes provide new insights into our understanding of the molecular mechanisms driving nitrogen use efficiency in barley and represent potential targets for genetic improvement.

Linkage mapping and genome-wide association study reveals conservative QTL and candidate genes for Fusarium rot resistance in maize

BACKGROUND

Fusarium ear rot (FER) caused by Fusarium verticillioides is a major disease of maize that reduces grain yield and quality globally. However, there have been few reports of major loci for FER were verified and cloned.

RESULT

To gain a comprehensive understanding of the genetic basis of natural variation in FER resistance, a recombinant inbred lines (RIL) population and one panel of inbred lines were used to map quantitative trait loci (QTL) for resistance. As a result, a total of 10 QTL were identified by linkage mapping under four environments, which were located on six chromosomes and explained 1.0-7.1% of the phenotypic variation. Epistatic mapping detected four pairs of QTL that showed significant epistasis effects, explaining 2.1-3.0% of the phenotypic variation. Additionally, 18 single nucleotide polymorphisms (SNPs) were identified across the whole genome by genome-wide association study (GWAS) under five environments. Compared linkage and association mapping revealed five common intervals located on chromosomes 3, 4, and 5 associated with FER resistance, four of which were verified in different near-isogenic lines (NILs) populations. GWAS identified three candidate genes in these consistent intervals, which belonged to the Glutaredoxin protein family, actin-depolymerizing factors (ADFs), and AMP-binding proteins. In addition, two verified FER QTL regions were found consistent with Fusarium cob rot (FCR) and Fusarium seed rot (FSR).

CONCLUSIONS

These results revealed that multi pathways were involved in FER resistance, which was a complex trait that was controlled by multiple genes with minor effects, and provided important QTL and genes, which could be used in molecular breeding for resistance.

Genome-Wide Association Study of Salinity Tolerance During Germination in Barley (Hordeum vulgare L.)

Barley seeds need to be able to germinate and establish seedlings in saline soils in Mediterranean-type climates. Despite being a major cereal crop, barley has few reported quantitative trait loci (QTL) and candidate genes underlying salt tolerance at the germination stage. Breeding programs targeting salinity tolerance at germination require an understanding of genetic loci and alleles in the current germplasm. In this study, we investigated seed-germination-related traits under control and salt stress conditions in 350 diverse barley accessions. A genome-wide association study, using ~24,000 genetic markers, was undertaken to detect marker-trait associations (MTA) and the underlying candidate genes for salinity tolerance during germination. We detected 19 loci containing 52 significant salt-tolerance-associated markers across all chromosomes, and 4 genes belonging to 4 family functions underlying the predicted MTAs. Our results provide new genetic resources and information to improve salt tolerance at germination in future barley varieties via genomic and marker-assisted selection and to open up avenues for further functional characterization of the identified candidate genes.

Combined GWAS and QTL analysis for dissecting the genetic architecture of kernel test weight in maize

Kernel weight in a unit volume is referred to as kernel test weight (KTW) that directly reflects maize (Zea mays L.) grain quality. In this study, an inter-mated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population and an association panel were used to identify loci responsible for KTW of maize across multiple environments. A total of 18 significant KTW-related single-nucleotide polymorphisms (SNPs) were identified using genome-wide association study (GWAS); they were closely linked to 12 candidate genes. In the IBM Syn10 DH population, linkage analysis detected 19 common quantitative trait loci (QTL), five of which were repeatedly detected among multiple environments. Several verified genes that regulate maize seed development were found in the confidence intervals of the mapped QTL and the LD regions of GWAS, such as ZmYUC1, BAP2, ZmTCRR-1, dek36 and ZmSWEET4c. Combined QTL mapping and GWAS identified one significant SNP that was co-identified in the both populations. Based on the co-localized SNP across the both populations, 17 candidate genes were identified. Of them, Zm00001d044075, Zm00001d044086, and Zm00001d044081 were further identified by candidate gene association study for KTW. Zm00001d044081 encodes homeobox-leucine zipper protein ATHB-4, which has been demonstrated to control apical embryo development in Arabidopsis. Our findings provided insights into the mechanism underlying maize KTW and contributed to the application of molecular-assisted selection of high KTW breeding in maize.

Genome-Wide Association Analysis of Plant Height Using the Maize F1 Population

Drastic changes in plant height (PH) are observed when maize adapt to a higher plant density. Most importantly, PH is an important factor affecting maize yield. Although the genetic basis of PH has been extensively studied using different populations during the past decades, genetic basis remains unclear in the F1 population, which was a widely used population in production. In this study, a genome-wide association study (GWAS) was conducted using an F1 population consisting of 300 maize hybrids with 17,652 single nucleotide polymorphisms (SNPs) makers to identify candidate genes for controlling PH. A total of nine significant SNPs makers and two candidate genes were identified for PH. The candidate genes, Zm00001d018617 and Zm00001d023659, were the genes most probable to be involved in the development of PH. Our results provide new insights into the genetic basis of PH in maize.

Single-plant GWAS coupled with bulk segregant analysis allows rapid identification and corroboration of plant-height candidate SNPs

BACKGROUND

Genome wide association studies (GWAS) are a powerful tool for identifying quantitative trait loci (QTL) and causal single nucleotide polymorphisms (SNPs)/genes associated with various important traits in crop species. Typically, GWAS in crops are performed using a panel of inbred lines, where multiple replicates of the same inbred are measured and the average phenotype is taken as the response variable. Here we describe and evaluate single plant GWAS (sp-GWAS) for performing a GWAS on individual plants, which does not require an association panel of inbreds. Instead sp-GWAS relies on the phenotypes and genotypes from individual plants sampled from a randomly mating population. Importantly, we demonstrate how sp-GWAS can be efficiently combined with a bulk segregant analysis (BSA) experiment to rapidly corroborate evidence for significant SNPs.

RESULTS

In this study we used the Shoepeg maize landrace, collected as an open pollinating variety from a farm in Southern Missouri in the 1960's, to evaluate whether sp-GWAS coupled with BSA can efficiently and powerfully used to detect significant association of SNPs for plant height (PH). Plant were grown in 8 locations across two years and in total 768 individuals were genotyped and phenotyped for sp-GWAS. A total of 306 k polymorphic markers in 768 individuals evaluated via association analysis detected 25 significant SNPs (P ≤ 0.00001) for PH. The results from our single-plant GWAS were further validated by bulk segregant analysis (BSA) for PH. BSA sequencing was performed on the same population by selecting tall and short plants as separate bulks. This approach identified 37 genomic regions for plant height. Of the 25 significant SNPs from GWAS, the three most significant SNPs co-localize with regions identified by BSA.

CONCLUSION

Overall, this study demonstrates that sp-GWAS coupled with BSA can be a useful tool for detecting significant SNPs and identifying candidate genes. This result is particularly useful for species/populations where association panels are not readily available.

Genome Wide Association Study of Karnal Bunt Resistance in a Wheat Germplasm Collection from Afghanistan

Karnal bunt disease of wheat, caused by the fungus Neovossia indica, is one of the most important challenges to the grain industry as it affects the grain quality and also restricts the international movement of infected grain. It is a seed-, soil- and airborne disease with limited effect of chemical control. Currently, this disease is contained through the deployment of host resistance but further improvement is limited as only a few genotypes have been found to carry partial resistance. To identify genomic regions responsible for resistance in a set of 339 wheat accessions, genome-wide association study (GWAS) was undertaken using the DArTSeq^® technology, in which 18 genomic regions for Karnal bunt resistance were identified, explaining 5–20% of the phenotypic variation. The identified quantitative trait loci (QTL) on chromosome 2BL showed consistently significant effects across all four experiments, whereas another QTL on 5BL was significant in three experiments. Additional QTLs were mapped on chromosomes 1DL, 2DL, 4AL, 5AS, 6BL, 6BS, 7BS and 7DL that have not been mapped previously, and on chromosomes 4B, 5AL, 5BL and 6BS, which have been reported in previous studies. Germplasm with less than 1% Karnal bunt infection have been identified and can be used for resistance breeding. The SNP markers linked to the genomic regions conferring resistance to Karnal bunt could be used to improve Karnal bunt resistance through marker-assisted selection.

Unravelling the Complex Genetics of Karnal Bunt (Tilletia indica) Resistance in Common Wheat (Triticum aestivum) by Genetic Linkage and Genome-Wide Association Analyses

Karnal bunt caused by Tilletia indica Mitra [syn. Neovossia indica (Mitra) Mundkur] is a significant biosecurity concern for wheat-exporting countries that are free of the disease. It is a seed-, soil-and air-borne disease with no effective chemical control measures. The current study used data from multi-year field experiments of two bi-parental populations and a genome-wide association (GWA) mapping panel to unravel the genetic basis for resistance in common wheat. Broad-sense heritability for Karnal bunt resistance in the populations varied from 0.52 in the WH542×HD29 population, to 0.61 in the WH542×W485 cross and 0.71 in a GWAS panel. Quantitative trait locus (QTL) analysis with seven years of phenotypic data identified a major locus on chromosome 3B (R² = 27.8%) and a minor locus on chromosome 1A (R² = 12.2%), in the WH542×HD29 population, with both parents contributing the high-value alleles. A major locus (R² = 27.8%) and seven minor loci (R² = 4.4–15.8%) were detected in the WH542×W485 population. GWA mapping validated QTL regions in the bi-parent populations, but also identified novel loci not previously associated with Karnal bunt resistance. Meta-QTL analysis aligned the results from this study with those reported in wheat over the last two decades. Two major clusters were detected, the first on chromosome 4B, which clustered with Qkb.ksu-4B, QKb.cimmyt-4BL, Qkb.cim-4BL, and the second on chromosome 3B, which clustered with Qkb.cnl-3B, QKb.cimmyt-3BS and Qkb.cim-3BS1. The results provide definitive chromosomal assignments for QTL/genes controlling Karnal bunt resistance in common wheat, and will be useful in pre-emptive breeding against the pathogen in wheat-producing areas that are free of the disease.

Identification and Characterization of EI (Elongated Internode) Gene in Tomato (Solanum lycopersicum)

Internode length is an important agronomic trait affecting plant architecture and crop yield. However, few genes for internode elongation have been identified in tomato. In this study, we characterized an elongated internode inbred line P502, which is a natural mutant of the tomato cultivar 05T606. The mutant P502 exhibits longer internode and higher bioactive GA concentration compared with wild-type 05T606. Genetic analysis suggested that the elongated internode trait is controlled by quantitative trait loci (QTL). Then, we identified a major QTL on chromosome 2 based on molecular markers and bulked segregant analysis (BSA). The locus was designated as EI (Elongated Internode), which explained 73.6% genetic variance. The EI was further mapped to a 75.8-kb region containing 10 genes in the reference Heinz 1706 genome. One single nucleotide polymorphism (SNP) in the coding region of solyc02g080120.1 was identified, which encodes gibberellin 2-beta-dioxygenase 7 (SlGA2ox7). SlGA2ox7, orthologous to AtGA2ox7 and AtGA2ox8, is involved in the regulation of GA degradation. Overexpression of the wild EI gene in mutant P502 caused a dwarf phenotype with a shortened internode. The difference of EI expression levels was not significant in the P502 and wild-type, but the expression levels of GA biosynthetic genes including CPS, KO, KAO, GA20ox1, GA20ox2, GA20ox4, GA3ox1, GA2ox1, GA2ox2, GA2ox4, and GA2ox5, were upregulated in mutant P502. Our results may provide a better understanding of the genetics underlying the internode elongation and valuable information to improve plant architecture of the tomato.

Genetic analysis of cob resistance to F. verticillioides: another step towards the protection of maize from ear rot

We lay the foundation for further research on maize resistance to Fusarium verticillioides cob rot by identifying a candidate resistance gene. Fusarium verticillioides ear rot is the most common type of maize ear rot in the Huanghuaihai Plain of China. Ear rot resistance includes cob and kernel resistance. Most of the current literature concentrates on kernel resistance, and genetic studies on cob resistance are scarce. We aimed on identifying the QTLs responsible for F. verticillioides cob rot (FCR) resistance. Twenty-eight genes associated with 48 single nucleotide polymorphisms (SNPs) were identified (P < 10^-4) to correlate with FCR resistance using a whole-genome association study. The major quantitative trait locus, qRcfv2, for FCR resistance was identified on chromosome 2 through linkage mapping and was validated in near-isogenic line populations. Two candidate genes associated with two SNPs were detected in the qRcfv2 region with a lower threshold (P < 10^-3). Through real-time fluorescence quantitative PCR, one candidate gene was found to have no expression in the cob but the other was expressed in response to F. verticillioides. These results lay a foundation for research on the resistance mechanisms of cob and provide resources for marker-assisted selection.

Genome-wide haplotype-based association analysis of key traits of plant lodging and architecture of maize identifies major determinants for leaf angle: hapLA4

Traits related to plant lodging and architecture are important determinants of plant productivity in intensive maize cultivation systems. Motivated by the identification of genomic associations with the leaf angle, plant height (PH), ear height (EH) and the EH/PH ratio, we characterized approximately 7,800 haplotypes from a set of high-quality single nucleotide polymorphisms (SNPs), in an association panel consisting of tropical maize inbred lines. The proportion of the phenotypic variations explained by the individual SNPs varied between 7%, for the SNP S1_285330124 (located on chromosome 9 and associated with the EH/PH ratio), and 22%, for the SNP S1_317085830 (located on chromosome 6 and associated with the leaf angle). A total of 40 haplotype blocks were significantly associated with the traits of interest, explaining up to 29% of the phenotypic variation for the leaf angle, corresponding to the haplotype hapLA4.04, which was stable over two growing seasons. Overall, the associations for PH, EH and the EH/PH ratio were environment-specific, which was confirmed by performing a model comparison analysis using the information criteria of Akaike and Schwarz. In addition, five stable haplotypes (83%) and 15 SNPs (75%) were identified for the leaf angle. Finally, approximately 62% of the associated haplotypes (25/40) did not contain SNPs detected in the association study using individual SNP markers. This result confirms the advantage of haplotype-based genome-wide association studies for examining genomic regions that control the determining traits for architecture and lodging in maize plants.

Identification of Loci and Candidate Genes Responsible for Fiber Length in Upland Cotton (Gossypium hirsutum L.) via Association Mapping and Linkage Analyses

Fiber length (FL) is an important fiber quality trait in cotton. Although many fiber quality quantitative trait loci (QTL) responsible for FL have been identified, most cannot be applied to breeding programs, mainly due to unstable environments or large confidence intervals. In this study, we combined a genome-wide association study (GWAS) and linkage mapping to identify and validate high-quality QTLs responsible for FL. For the GWAS, we developed 93,250 high-quality single-nucleotide polymorphism (SNP) markers based on 355 accessions, and the FL was measured in eight different environments. For the linkage mapping, we constructed an F ₂ population from two extreme accessions. The high-density linkage maps spanned 3,848.29 cM, with an average marker interval of 1.41 cM. In total, 14 and 13 QTLs were identified in the association and linkage mapping analyses, respectively. Most importantly, a major QTL on chromosome D03 identified in both populations explained more than 10% of the phenotypic variation (PV). Furthermore, we found that a sucrose synthesis-related gene (Gh_D03G1338) was associated with FL in this QTL region. The RNA-seq data showed that Gh_D03G1338 was highly expressed during the fiber development stage, and the qRT-PCR analysis showed significant expression differences between the long fiber and short fiber varieties. These results suggest that Gh_D03G1338 may determine cotton fiber elongation by regulating the synthesis of sucrose. Favorable QTLs and candidate genes should be useful for increasing fiber quality in cotton breeding.

Combined linkage and association mapping reveal QTL for host plant resistance to common rust (Puccinia sorghi) in tropical maize

BACKGROUND

Common rust, caused by Puccinia sorghi, is an important foliar disease of maize that has been associated with up to 50% grain yield loss. Development of resistant maize germplasm is the ideal strategy to combat P. sorghi.

RESULTS

Association mapping performed using a mixed linear model (MLM), integrating population structure and family relatedness identified 25 QTL (P < 3.12 × 10- 5) that were associated with resistance to common rust and distributed on chromosomes 1, 3, 5, 6, 8, and 10. We identified three QTLs associated with all three disease parameters (final disease rating, mean disease rating, and area under disease progress curve) located on chromosomes 1, 3, and 8. A total of 5 QTLs for resistance to common rust were identified in the RIL population. Nine candidate genes located on chromosomes 1, 5, 6, 8, and 10 for resistance to common rust associated loci were identified through detailed annotation.

CONCLUSIONS

Using a diverse set of inbred lines genotyped with high density markers and evaluated for common rust resistance in multiple environments, it was possible to identify QTL significantly associated with resistance to common rust and several candidate genes. The results point to the need for fine mapping common rust resistance by targeting regions identified in common between this study and others using diverse germplasm.

Discovery of Anthocyanin Acyltransferase1 (AAT1) in Maize Using Genotyping-by-Sequencing (GBS)

The reduced acylation phenotype describes the inability of certain accessions of maize (Zea mays [L.]) to produce significant amounts of acylated anthocyanins, which are typically the most abundant pigments. Acylated anthocyanins are important for their association with stability and are therefore important for the various industries using anthocyanins as natural colorants to replace synthetic dyes. Many anthocyanin acyltransferases have been characterized in other species; however, no anthocyanin acyltransferases have been characterized in maize. Therefore, a mapping population was developed from a cross between mutant stock 707G and wild-type acylation line B73 to identify the locus associated with the reduced acylation trait. High-performance liquid chromatography was used to assay the pigment content and composition of 129 F₂ lines generated in the mapping population. Recessive alleles of Colorless1, Colored1, and the reduced acylation mutant all decreased anthocyanin content while Intensifier1 increased anthocyanin content in aleurone tissue. The association of increased proportions of acylation with increased anthocyanin content indicates acylation may be important for increasing the stability of anthocyanins in vivo. Genotyping-by-sequencing was used to create SNP markers to map the reduced acylation locus. In the QTL analysis, a segment of Chromosome 1 containing transferase family protein GRMZM2G387394 was found to be significant. A UniformMu Mu transposon knockout of GRMZM2G387394 demonstrated this gene has anthocyanidin malonyltransferase activity and will therefore be named Anthocyanin Acyltransferase1 (AAT1). AAT1 is the first anthocyanin acyltransferase characterized in a monocot species and will increase our knowledge of all acyltransferase family members.

Augmentation of crop productivity through interventions of omics technologies in India: challenges and opportunities

With the continuous increase in the population of developing countries and decline of natural resources, there is an urgent need to qualitatively and quantitatively augment crop productivity by using new tools and technologies for improvement of agriculturally important traits. The new scientific and technological omics-based approaches have enabled us to deal with several issues and challenges faced by modern agricultural system and provided us novel opportunities for ensuring food and nutritional security. Recent developments in sequencing techniques have made available huge amount of genomic and transcriptomic data on model and cultivated crop plants including Arabidopsis thaliana, Oryza sativa, Triticum aestivum etc. The sequencing data along with other data generated through several omics platforms have significantly influenced the disciplines of crop sciences. Gene discovery and expression profiling-based technologies are offering enormous opportunities to the scientific community which can now apply marker-assisted selection technology to assess and enhance diversity in their collected germplasm, introgress essential traits from new sources and investigate genes that control key traits of crop plants. Utilization of omics science and technologies for crop productivity, protection and management has recently been receiving a lot of attention; the majority of the efforts have been put into signifying the possible applications of various omics technologies in crop plant sciences. This article highlights the background of challenges and opportunities for augmentation of crop productivity through interventions of omics technologies in India.

Linkage mapping combined with association analysis reveals QTL and candidate genes for three husk traits in maize

Key message Combined linkage and association mapping analyses facilitate the emphasis on the candidate genes putatively involved in maize husk growth. The maize (Zea mays L.) husk consists of multiple leafy layers and plays important roles in protecting the ear from pathogen infection and in preventing grain dehydration. Although husk morphology varies widely among different maize inbred lines, the genetic basis of such variation is poorly understood. In this study, we used three maize recombinant inbred line (RIL) populations to dissect the genetic basis of three husk traits: i.e., husk length (HL), husk width (HW), and the number of husk layers (HN). Three husk traits in all three RIL populations showed wide phenotypic variation and high heritability. The HL showed stronger correlations with ear traits than did HW and HN. A total of 21 quantitative trait loci (QTL) were identified for the three traits in three RIL populations, and some of them were commonly observed for the same trait in different populations. The proportions of total phenotypic variation explained by QTL in three RIL populations were 31.8, 35.3, and 44.5% for HL, HW, and HN, respectively. The highest proportions of phenotypic variation explained by a single QTL were 14.7% for HL in the By815/K22 RIL population (BYK), 13.5% for HW in the By815/DE3 RIL population (BYD), and 19.4% for HN in the BYD population. A combined analysis of linkage mapping with a previous genome-wide association study revealed five candidate genes related to husk morphology situated within three QTL loci. These five genes were related to metabolism, gene expression regulation, and signal transduction.

Quantitative epigenetics and evolution

Epigenetics refers to chemical modifications of chromatin or transcribed DNA that can influence gene activity and expression without changes in DNA sequence. The last 20 years have yielded breakthroughs in our understanding of epigenetic processes that impact many fields of biology. In this review, we discuss how epigenetics relates to quantitative genetics and evolution. We argue that epigenetics is important for quantitative genetics because: (1) quantitative genetics is increasingly being combined with genomics, and therefore we should expand our thinking to include cellular-level mechanisms that can account for phenotypic variance and heritability besides just those that are hard-coded in the DNA sequence; and (2) epigenetic mechanisms change how phenotypic variance is partitioned, and can thereby change the heritability of traits and how those traits are inherited. To explicate these points, we show that epigenetics can influence all aspects of the phenotypic variance formula: VP (total phenotypic variance) = VG (genetic variance) + VE (environmental variance) + VGxE (genotype-by-environment interaction) + 2COVGE (the genotype-environment covariance) + Vɛ (residual variance), requiring new strategies to account for different potential sources of epigenetic effects on phenotypic variance. We also demonstrate how each of the components of phenotypic variance not only can be influenced by epigenetics, but can also have evolutionary consequences. We argue that no sources of epigenetic effects on phenotypic variance can be easily cast aside in a quantitative genetic research program that seeks to understand evolutionary processes.

Identification and genetic mapping for rht-DM, a dominant dwarfing gene in mutant semi-dwarf maize using QTL-seq approach

Semi-dwarfism is an agronomically important trait in breeding for stable high yields and for resistance to damage by wind and rain (lodging resistance). Many QTLs and genes causing dwarf phenotype have been found in maize. However, because of the yield loss associated with these QTLs and genes, they have been difficult to use in breeding for dwarf stature in maize. Therefore, it is important to find the new dwarfing genes or materials without undesirable characters. The objectives of this study were: (1) to figure out the inheritance of semi-dwarfism in mutants; (2) mapping dwarfing gene or QTL. Maize inbred lines '18599' and 'DM173', which is the dwarf mutant derived from the maize inbred line '173' through ⁶⁰Co-γ ray irradiation. F₂ and BC₁F₁ population were used for genetic analysis. Whole genome resequencing-based technology (QTL-seq) were performed to map dwarfing gene and figured out the SNP markers in predicted region using dwarf bulk and tall bulk from F₂ population. Based on the polymorphic SNP markers from QTL-seq, we were fine-mapping the dwarfing gene using F₂ population. In F₂ population, 398 were dwarf plants and 135 were tall plants. Results of χ² tests indicated that the ratio of dwarf plants to tall plants was fitted to 3:1 ratio. Furthermore, the χ² tests of BC₁F₁ population showed that the ratio was fitted to 1:1 ratio. Based on QTL-seq, the dwarfing gene was located at the region from 111.07 to 124.56 Mb of chromosome 9, and we named it rht-DM. Using traditional QTL mapping with SNP markers, the rht-DM was narrowed down to 400 kb region between SNP-21 and SNP-24. The two SNPs were located at 0.43 and 0.11 cM. Segregation analysis of F₂ and BC₁F₁ indicated that the dwarfing gene was likely a dominant gene. This dwarfing gene was located in the region between 115.02 and 115.42 Mb on chromosome 9.

Genetic-based dissection of arsenic accumulation in maize using a genome-wide association analysis method

Understanding the mechanism of arsenic (As) accumulation in plants is important in reducing As's toxicity to plants and its potential risks to human health. Here, we performed a genome-wide association study to dissect the genetic basis of the As contents of different maize tissues in Xixian, which was irrigated with As-rich surface water, and Changge using an association population consisting of 230 representative maize inbred lines. Phenotypic data revealed a wide normal distribution and high repeatability for the As contents in maize tissues. The As concentrations in maize tissues followed the same trend in the two locations: kernels < axes < stems < bracts < leaves. In total, 15, 16 and 15 non-redundant quantitative trait loci (QTLs) associated with As concentrations were identified (P ≤ 2.04 × 10-6 ) in five tissues from Xixian, Changge, and the combination of the locations, respectively, explaining 9.70%-24.65% of the phenotypic variation for each QTL, on average. Additionally, four QTLs [involving 15 single nucleotide polymorphisms (SNPs)] were detected in the single and the combined locations, indicating that these loci/SNPs might be stable across different environments. The candidate genes associated with these four loci were predicted. In addition, four non-redundant QTLs (6 SNPs), including a QTL that was detected in multiple locations according to the genome-wide association study, were found to co-localize with four previously reported QTL intervals. These results are valuable to understand the genetic architecture of As mechanism in maize and facilitate the genetic improvement of varieties without As toxicity.

Genomics-Enabled Next-Generation Breeding Approaches for Developing System-Specific Drought Tolerant Hybrids in Maize

Breeding science has immensely contributed to the global food security. Several varieties and hybrids in different food crops including maize have been released through conventional breeding. The ever growing population, decreasing agricultural land, lowering water table, changing climate, and other variables pose tremendous challenge to the researchers to improve the production and productivity of food crops. Drought is one of the major problems to sustain and improve the productivity of food crops including maize in tropical and subtropical production systems. With advent of novel genomics and breeding tools, the way of doing breeding has been tremendously changed in the last two decades. Drought tolerance is a combination of several component traits with a quantitative mode of inheritance. Rapid DNA and RNA sequencing tools and high-throughput SNP genotyping techniques, trait mapping, functional characterization, genomic selection, rapid generation advancement, and other tools are now available to understand the genetics of drought tolerance and to accelerate the breeding cycle. Informatics play complementary role by managing the big-data generated from the large-scale genomics and breeding experiments. Genome editing is the latest technique to alter specific genes to improve the trait expression. Integration of novel genomics, next-generation breeding, and informatics tools will accelerate the stress breeding process and increase the genetic gain under different production systems.

Genome-wide associations for multiple pest resistances in a Northwestern United States elite spring wheat panel

Northern areas of the western United States are one of the most productive wheat growing regions in the United States. Increasing productivity through breeding is hindered by several biotic stresses which slow and constrain targeted yield improvement. In order to understand genetic variation for stripe rust (Puccinia striiformis f. sp. tritici), Septoria tritici blotch (Mycosphaerella graminicola), and Hessian fly (Mayetiola destructor) in regional germplasm, a panel of 408 elite spring wheat lines was characterized and genotyped with an Illumina 9K wheat single nucleotide polymorphism (SNP) chip to enable genome-wide association study (GWAS) analyses. Significant marker-trait associations were identified for stripe rust (38 loci), Septoria tritici blotch (8) and Hessian fly (9) resistance. Many of the QTL corresponded with previously reported gene locations or QTL, but we also discovered new resistance loci for each trait. We validated one of the stripe rust resistance loci detected by GWAS in a bi-parental mapping population, which confirmed the detection of Yr15 in the panel. This study elucidated well-defined chromosome regions for multiple pest resistances in elite Northwest germplasm. Newly identified resistance loci, along with SNPs more tightly linked to previously reported genes or QTL will help future breeding and marker assisted selection efforts.

Genetic Dissection of Maize Embryonic Callus Regenerative Capacity Using Multi-Locus Genome-Wide Association Studies

The regenerative capacity of the embryonic callus, a complex quantitative trait, is one of the main limiting factors for maize transformation. This trait was decomposed into five traits, namely, green callus rate (GCR), callus differentiating rate (CDR), callus plantlet number (CPN), callus rooting rate (CRR), and callus browning rate (CBR). To dissect the genetic foundation of maize transformation, in this study multi-locus genome-wide association studies (GWAS) for the five traits were performed in a population of 144 inbred lines genotyped with 43,427 SNPs. Using the phenotypic values in three environments and best linear unbiased prediction (BLUP) values, as a result, a total of 127, 56, 160, and 130 significant quantitative trait nucleotides (QTNs) were identified by mrMLM, FASTmrEMMA, ISIS EM-BLASSO, and pLARmEB, respectively. Of these QTNs, 63 QTNs were commonly detected, including 15 across multiple environments and 58 across multiple methods. Allele distribution analysis showed that the proportion of superior alleles for 36 QTNs was <50% in 31 elite inbred lines. Meanwhile, these superior alleles had obviously additive effect on the regenerative capacity. This indicates that the regenerative capacity-related traits can be improved by proper integration of the superior alleles using marker-assisted selection. Moreover, a total of 40 candidate genes were found based on these common QTNs. Some annotated genes were previously reported to relate with auxin transport, cell fate, seed germination, or embryo development, especially, GRMZM2G108933 (WOX2) was found to promote maize transgenic embryonic callus regeneration. These identified candidate genes will contribute to a further understanding of the genetic foundation of maize embryonic callus regeneration.

Genetic Diversity and Genome-Wide Association Study of Major Ear Quantitative Traits Using High-Density SNPs in Maize

Kernel and ear traits are key components of grain yield in maize (Zea mays L.). Investigation of these traits would help to develop high-yield varieties in maize. Genome-wide association study (GWAS) uses the linkage disequilibrium (LD) in the whole genome to determine the genes affecting certain phenotype. In this study, five ear traits (kernel length and width, ear length and diameter, cob diameter) were investigated across multi-environments for 2 years. Combining with the genotype obtained from Maize SNP50 chip, genetic diversity and association mapping in a set of 292 inbred lines were performed. Results showed that maize lines were clustered into seven subgroups and a total of 20 SNPs were found to be associated with ear traits significantly (P < 3.95E-05). The candidate genes identified by GWAS mainly encoded ubiquitin-activation enzymes (GRMZM2G015287), carotenoid cleavage dioxygenase (GRMZM2G446858), MYB-CC type transfactor, and phosphate starvation response protein 3, and they were associated with kernel length (KL) and ear diameter (ED), respectively. Moreover, two novel genes corresponding to RNA processing and fructose metabolism were found. Further, the SNPs detected by GWAS were confirmed by meta-QTL analysis. These genes and SNPs identified in the study would offer essential information for yield-related genes clone and breeding program in maize.

Association mapping unveils favorable alleles for grain iron and zinc concentrations in lentil (Lens culinaris subsp. culinaris)

Lentil is a major cool-season grain legume grown in South Asia, West Asia, and North Africa. Populations in developing countries of these regions have micronutrient deficiencies; therefore, breeding programs should focus more on improving the micronutrient content of food. In the present study, a set of 96 diverse germplasm lines were evaluated at three different locations in India to examine the variation in iron (Fe) and zinc (Zn) concentration and identify simple sequence repeat (SSR) markers that associate with the genetic variation. The genetic variation among genotypes of the association mapping (AM) panel was characterized using a genetic distance-based and a general model-based clustering method. The model-based analysis identified six subpopulations, which satisfactorily explained the genetic structure of the AM panel. AM analysis identified three SSRs (PBALC 13, PBALC 206, and GLLC 563) associated with grain Fe concentration explaining 9% to 11% of phenotypic variation and four SSRs (PBALC 353, SSR 317-1, PLC 62, and PBALC 217) were associated with grain Zn concentration explaining 14%, to 21% of phenotypic variation. These identified SSRs exhibited consistent performance across locations. These candidate SSRs can be used in marker-assisted genetic improvement for developing Fe and Zn fortified lentil varieties. Favorable alleles and promising genotypes identified in this study can be utilized for lentil biofortification.

Assessing and Exploiting Functional Diversity in Germplasm Pools to Enhance Abiotic Stress Adaptation and Yield in Cereals and Food Legumes

There is a need to accelerate crop improvement by introducing alleles conferring host plant resistance, abiotic stress adaptation, and high yield potential. Elite cultivars, landraces and wild relatives harbor useful genetic variation that needs to be more easily utilized in plant breeding. We review genome-wide approaches for assessing and identifying alleles associated with desirable agronomic traits in diverse germplasm pools of cereals and legumes. Major quantitative trait loci and single nucleotide polymorphisms (SNPs) associated with desirable agronomic traits have been deployed to enhance crop productivity and resilience. These include alleles associated with variation conferring enhanced photoperiod and flowering traits. Genetic variants in the florigen pathway can provide both environmental flexibility and improved yields. SNPs associated with length of growing season and tolerance to abiotic stresses (precipitation, high temperature) are valuable resources for accelerating breeding for drought-prone environments. Both genomic selection and genome editing can also harness allelic diversity and increase productivity by improving multiple traits, including phenology, plant architecture, yield potential and adaptation to abiotic stresses. Discovering rare alleles and useful haplotypes also provides opportunities to enhance abiotic stress adaptation, while epigenetic variation has potential to enhance abiotic stress adaptation and productivity in crops. By reviewing current knowledge on specific traits and their genetic basis, we highlight recent developments in the understanding of crop functional diversity and identify potential candidate genes for future use. The storage and integration of genetic, genomic and phenotypic information will play an important role in ensuring broad and rapid application of novel genetic discoveries by the plant breeding community. Exploiting alleles for yield-related traits would allow improvement of selection efficiency and overall genetic gain of multigenic traits. An integrated approach involving multiple stakeholders specializing in management and utilization of genetic resources, crop breeding, molecular biology and genomics, agronomy, stress tolerance, and reproductive/seed biology will help to address the global challenge of ensuring food security in the face of growing resource demands and climate change induced stresses.

Genetic Mapping of the Leaf Number above the Primary Ear and Its Relationship with Plant Height and Flowering Time in Maize

The leaf number above the primary ear (LA) is a major contributing factor to plant architecture in maize. The yield of leafy maize, which has extra LA compared to normal maize, is higher than normal maize in some regions. One major concern is that increasing LA may be accompanied by increased plant height and/or flowering time. Using an F_2:3 population comprising 192 families derived from a leafy maize line and a normal maize line, an association population comprising 437 inbred maize lines, and a pair of near-isogenic maize lines, we mapped the quantitative trait loci (QTL) associated with LA and assessed its genetic relationship with flowering time and plant height. Ten QTL with an additive and dominant effect, 18 pairs of interacting QTL in the F_2:3 population and seventeen significant SNPs in the association population were detected for LA. Two major QTL, qLA3-4 and qLA7-1, were repeatedly detected and explained a large proportion of the phenotypic variation. The qLA3-4 was centered on lfy1, which is a dominant gene underlying extra leaves above the ear in leafy maize. Four LA QTL were found to overlap with flowering time and/or plant height, which suggested that these QTL might have a pleiotropic effect. The pleiotropy of the lfy1 locus on LA, flowering time and plant height were validated by near-isogenic line analysis. These results enhance our understanding of the genetic architecture affecting maize LA and the development of maize hybrids with increased LA.

Genome-Wide Analysis of Tar Spot Complex Resistance in Maize Using Genotyping-by-Sequencing SNPs and Whole-Genome Prediction

Tar spot complex (TSC) is one of the most destructive foliar diseases of maize ( L.) in tropical and subtropical areas of Central and South America, causing significant grain yield losses when weather conditions are conducive. To dissect the genetic architecture of TSC resistance in maize, association mapping, in conjunction with linkage mapping, was conducted on an association-mapping panel and three biparental doubled-haploid (DH) populations using genotyping-by-sequencing (GBS) single-nucleotide polymorphisms (SNPs). Association mapping revealed four quantitative trait loci (QTL) on chromosome 2, 3, 7, and 8. All the QTL, except for the one on chromosome 3, were further validated by linkage mapping in different genetic backgrounds. Additional QTL were identified by linkage mapping alone. A major QTL located on bin 8.03 was consistently detected with the largest phenotypic explained variation: 13% in association-mapping analysis and 13.18 to 43.31% in linkage-mapping analysis. These results indicated that TSC resistance in maize was controlled by a major QTL located on bin 8.03 and several minor QTL with smaller effects on other chromosomes. Genomic prediction results showed moderate-to-high prediction accuracies in different populations using various training population sizes and marker densities. Prediction accuracy of TSC resistance was >0.50 when half of the population was included into the training set and 500 to 1,000 SNPs were used for prediction. Information obtained from this study can be used for developing functional molecular markers for marker-assisted selection (MAS) and for implementing genomic selection (GS) to improve TSC resistance in tropical maize.

Quantitative trait locus analysis of heterosis for plant height and ear height in an elite maize hybrid zhengdan 958 by design III

Background

Plant height (PH) and ear height (EH) are two important agronomic traits in maize selection breeding. F₁ hybrid exhibit significant heterosis for PH and EH as compared to their parental inbred lines. To understand the genetic basis of heterosis controlling PH and EH, we conducted quantitative trait locus (QTL) analysis using a recombinant inbreed line (RIL) based design III population derived from the elite maize hybrid Zhengdan 958 in five environments.

Results

A total of 14 environmentally stable QTLs were identified, and the number of QTLs for Z₁ and Z₂ populations was six and eight, respectively. Notably, all the eight environmentally stable QTLs for Z₂ were characterized by overdominance effect (OD), suggesting that overdominant QTLs were the most important contributors to heterosis for PH and EH. Furthermore, 14 environmentally stable QTLs were anchored on six genomic regions, among which four are trait-specific QTLs, suggesting that the genetic basis for PH and EH is partially different. Additionally, qPH.A-1.3, modifying about 10 centimeters of PH, was further validated in backcross populations.

Conclusions

The genetic basis for PH and EH is partially different, and overdominant QTLs are important factors for heterosis of PH and EH. A major QTL qPH.A-1.3 may be a desired target for genetic improvement of maize plant height.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0503-9) contains supplementary material, which is available to authorized users.

Regional Association Analysis of MetaQTLs Delineates Candidate Grain Size Genes in Rice

Molecular mapping studies which aim to identify genetic basis of diverse agronomic traits are vital for marker-assisted crop improvement. Numerous Quantitative Trait Loci (QTLs) mapped in rice span long genomic intervals with hundreds to thousands of genes, which limits their utilization for marker-assisted genetic enhancement of rice. Although potent, fine mapping of QTLs is challenging task as it requires screening of large number of segregants to identify suitable recombination events. Association mapping offers much higher resolution as compared to QTL mapping, but detects considerable number of spurious QTLs. Therefore, combined use of QTL and association mapping strategies can provide advantages associated with both these methods. In the current study, we utilized meta-analysis approach to identify metaQTLs associated with grain size/weight in diverse Indian indica and aromatic rice accessions. Subsequently, attempt has been made to narrow-down identified grain size/weight metaQTLs through individual SNP- as well as haplotype-based regional association analysis. The study identified six different metaQTL regions, three of which were successfully revalidated, and substantially scaled-down along with GS3 QTL interval (positive control) by regional association analysis. Consequently, two potential candidate genes within two reduced metaQTLs were identified based on their differential expression profiles in different tissues/stages of rice accessions during seed development. The developed strategy has broader practical utility for rapid delineation of candidate genes and natural alleles underlying QTLs associated with complex agronomic traits in rice as well as major crop plants enriched with useful genetic and genomic information.