Genome-wide association study identifies candidate genes that affect plant height in Chinese elite maize (Zea mays L.) inbred lines

Identification of Multiple Genetic Loci Related to Low-Temperature Tolerance during Germination in Maize (Zea maize L.) through a Genome-Wide Association Study

Low-temperature stress during the germination stage is an important abiotic stress that affects the growth and development of northern spring maize and seriously restricts maize yield and quality. Although some quantitative trait locis (QTLs) related to low-temperature tolerance in maize have been detected, only a few can be commonly detected, and the QTL intervals are large, indicating that low-temperature tolerance is a complex trait that requires more in-depth research. In this study, 296 excellent inbred lines from domestic and foreign origins (America and Europe) were used as the study materials, and a low-coverage resequencing method was employed for genome sequencing. Five phenotypic traits related to low-temperature tolerance were used to assess the genetic diversity of maize through a genome-wide association study (GWAS). A total of 14 SNPs significantly associated with low-temperature tolerance were detected (-log10(P) > 4), and an SNP consistently linked to low-temperature tolerance in the field and indoors during germination was utilized as a marker. This SNP, 14,070, was located on chromosome 5 at position 2,205,723, which explained 4.84-9.68% of the phenotypic variation. The aim of this study was to enrich the genetic theory of low-temperature tolerance in maize and provide support for the innovation of low-temperature tolerance resources and the breeding of new varieties.

Identification of quantitative trait loci for related traits of stalk lodging resistance using genome-wide association studies in maize (Zea mays L.)

BACKGROUND

Stalk lodging is one of the main factors affecting maize (Zea mays L.) yield and limiting mechanized harvesting. Developing maize varieties with high stalk lodging resistance requires exploring the genetic basis of lodging resistance-associated agronomic traits. Stalk strength is an important indicator to evaluate maize lodging and can be evaluated by measuring stalk rind penetrometer resistance (RPR) and stalk buckling strength (SBS). Along with morphological traits of the stalk for the third internodes length (TIL), fourth internode length (FIL), third internode diameter (TID), and the fourth internode diameter (FID) traits are associated with stalk lodging resistance.

RESULTS

In this study, a natural population containing 248 diverse maize inbred lines genotyped with 83,057 single nucleotide polymorphism (SNP) markers was used for genome-wide association study (GWAS) for six stalk lodging resistance-related traits. The heritability of all traits ranged from 0.59 to 0.72 in the association mapping panel. A total of 85 significant SNPs were identified for the association mapping panel using best linear unbiased prediction (BLUP) values of all traits. Additionally, five candidate genes were associated with stalk strength traits, which were either directly or indirectly associated with cell wall components.

CONCLUSIONS

These findings contribute to our understanding of the genetic basis of maize stalk lodging and provide valuable theoretical guidance for lodging resistance in maize breeding in the future.

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.

Breeding effects on canopy light attenuation in maize: a retrospective and prospective analysis

The light attenuation process within a plant canopy defines energy capture and vertical distribution of light and nitrogen (N). The vertical light distribution can be quantitatively described with the extinction coefficient (k), which associates the fraction of intercepted photosynthetically active radiation (fPARi) with the leaf area index (LAI). Lower values of k correspond to upright leaves and homogeneous vertical light distribution, increasing radiation use efficiency (RUE). Yield gains in maize (Zea mays L.) were accompanied by increases in optimum plant density and leaf erectness. Thus, the yield-driven breeding programs and management changes, such as reduced row spacing, selected a more erect leaf habit under different maize production systems (e.g., China and the USA). In this study, data from Argentina revealed that k decreased at a rate of 1.1% year-1 since 1989, regardless of plant density and in agreement with Chinese reports (1.0% year-1 since 1981). A reliable assessment of changes in k over time is critical for predicting (i) modifications in resource use efficiency (e.g. radiation, water, and N), improving estimations derived from crop simulation models; (ii) differences in productivity caused by management practices; and (iii) limitations to further exploit this trait with breeding.

Genetics of Germination and Seedling Traits under Drought Stress in a MAGIC Population of Maize

Drought is one of the most detrimental abiotic stresses hampering seed germination, development, and productivity. Maize is more sensitive to drought than other cereals, especially at seedling stage. Our objective was to study genetic regulation of drought tolerance at germination and during seedling growth in maize. We evaluated 420 RIL with their parents from a multi-parent advanced generation inter-cross (MAGIC) population with PEG-induced drought at germination and seedling establishment. A genome-wide association study (GWAS) was carried out to identify genomic regions associated with drought tolerance. GWAS identified 28 and 16 SNPs significantly associated with germination and seedling traits under stress and well-watered conditions, respectively. Among the SNPs detected, two SNPs had significant associations with several traits with high positive correlations, suggesting a pleiotropic genetic control. Other SNPs were located in regions that harbored major QTLs in previous studies, and co-located with QTLs for cold tolerance previously published for this MAGIC population. The genomic regions comprised several candidate genes related to stresses and plant development. These included numerous drought-responsive genes and transcription factors implicated in germination, seedling traits, and drought tolerance. The current analyses provide information and tools for subsequent studies and breeding programs for improving drought tolerance.

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.

Validation of Functional Polymorphisms Affecting Maize Plant Height by Unoccupied Aerial Systems (UAS) Discovers Novel Temporal Phenotypes

Plant height (PHT) in maize (Zea mays L.) has been scrutinized genetically and phenotypically due to relationship with other agronomically valuable traits (e.g. yield). Heritable variation of PHT is determined by many discovered quantitative trait loci (QTLs); however, phenotypic effects of such loci often lack validation across environments and genetic backgrounds, especially in the hybrid state grown by farmers rather than the inbred state more often used by geneticists. A previous genome wide association study using a topcrossed hybrid diversity panel identified two novel quantitative trait variants (QTVs) controlling both PHT and grain yield. Here, heterogeneous inbred families demonstrated that these two loci, characterized by two single nucleotide polymorphisms (SNPs), cause phenotypic variation in inbred lines, but that size of these effects were variable across four different genetic backgrounds, ranging from 1 to 10 cm. Weekly unoccupied aerial system flights demonstrated the two SNPs had larger effects, varying from 10 to 25 cm, in early growth while effects decreased towards the end of the season. These results show that allelic effect sizes of economically valuable loci are both dynamic in temporal growth and dynamic across genetic backgrounds, resulting in informative phenotypic variability overlooked following traditional phenotyping methods. Public genotyping data shows recent favorable allele selection in elite temperate germplasm with little change across tropical backgrounds. As these loci remain rarer in tropical germplasm, with effects most visible early in growth, they are useful for breeding and selection to expand the genetic basis of maize.

Mapping QTL/QTN and mining candidate genes for plant height and its response to planting densities in soybean [Glycine max (L.) Merr.] through a FW-RIL population

Plant height (PH) determines the morphology and seed yield of soybean, so it is an important breeding target, which is controlled by multiple genes and affected by plant density. In this research, it was used about a four-way recombinant inbred lines (FW-RIL) with 144 families constructed by double cross (Kenfeng 14 × Kenfeng 15) × (Heinong 48 × Kenfeng 19) as experimental materials, with the purpose to map QTL/QTN associated with PH under densities of 2.2×105 plant/ha (D1) and 3×105 plant/ha (D2) in five environments. The results showed that response of PH to densities varied in accordance to genotypes among environments. A total of 26 QTLs and 13 QTNs were identified specifically in D1; 20 QTLs and 21 QTNs were identified specifically in D2. Nine QTLs and one QTN were discovered commonly in two densities. Fifteen QTLs and 9 QTNs were repeatedly detected by multiple statistical methods, densities, or environments, which could be considered stable. Eighteen QTLs were detected, as well as 7 QTNs underlying responses of PH to density increment. Six QTNs, co-located in the interval of QTL, were detected in more than two environments or methods with a longer genome length over 3000 kb. Based on gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, five genes were predicted as candidates, which were likely to be involved in growth and development of PH. The results will help elucidate the genetic basis and improve molecular assistant selection of PH.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01209-0.

Genome-Wide Association Study of Maize Aboveground Dry Matter Accumulation at Seedling Stage

Dry matter accumulation and partitioning during the early phases of development could significantly affect crop growth and productivity. In this study, the aboveground dry matter (DM), the DM of different organs, and partition coefficients of a maize association mapping panel of 412 inbred lines were evaluated at the third and sixth leaf stages (V3 and V6). Further, the properties of these phenotypic traits were analyzed. Genome-wide association studies (GWAS) were conducted on the total aboveground biomass and the DM of different organs. Analysis of GWAS results identified a total of 1,103 unique candidate genes annotated by 678 significant SNPs (P value < 1.28e-6). A total of 224 genes annotated by SNPs at the top five of each GWAS method and detected by multiple GWAS methods were regarded as having high reliability. Pathway enrichment analysis was also performed to explore the biological significance and functions of these candidate genes. Several biological pathways related to the regulation of seed growth, gibberellin-mediated signaling pathway, and long-day photoperiodism were enriched. The results of our study could provide new perspectives on breeding high-yielding maize varieties.

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.

QTL mapping and GWAS for field kernel water content and kernel dehydration rate before physiological maturity in maize

Kernel water content (KWC) and kernel dehydration rate (KDR) are two main factors affecting maize seed quality and have a decisive influence on the mechanical harvest. It is of great importance to map and mine candidate genes related to KWCs and KDRs before physiological maturity in maize. 120 double-haploid (DH) lines constructed from Si287 with low KWC and JiA512 with high KWC were used as the mapping population. KWCs were measured every 5 days from 10 to 40 days after pollination, and KDRs were calculated. A total of 1702 SNP markers were used to construct a linkage map, with a total length of 1,309.02 cM and an average map distance of 0.77 cM. 10 quantitative trait loci (QTLs) and 27 quantitative trait nucleotides (QTNs) were detected by genome-wide composite interval mapping (GCIM) and multi-locus random-SNP-effect mixed linear model (mrMLM), respectively. One and two QTL hotspot regions were found on Chromosome 3 and 7, respectively. Analysis of the Gene Ontology showed that 2 GO terms of biological processes (BP) were significantly enriched (P ≤ 0.05) and 6 candidate genes were obtained. This study provides theoretical support for marker-assisted breeding of mechanical harvest variety in maize.

Identification of candidate tolerance genes to low-temperature during maize germination by GWAS and RNA-seqapproaches

BACKGROUND

Maize (Zea mays L.) is one of the main agricultural crops with the largest yield and acreage in the world. However, maize germplasm is very sensitive to low temperatures, mainly during germination, and low temperatures significantly affect plant growth and crop yield. Therefore, the identification of genes capable of increasing tolerance to low temperature has become necessary.

RESULTS

In this study, fourteen phenotypic traits related to seed germination were used to assess the genetic diversity of maize through genome-wide association study (GWAS). A total of 30 single-nucleotide polymorphisms (SNPs) linked to low-temperature tolerance were detected (-log10(P) > 4), fourteen candidate genes were found to be directly related to the SNPs, further additional 68 genes were identified when the screen was extended to include a linkage disequilibrium (LD) decay distance of r2 ≥ 0.2 from the SNPs. RNA-sequencing (RNA-seq) analysis was then used to confirm the linkage between the candidate gene and low-temperature tolerance. A total of ten differentially expressed genes (DEGs) (|log2 fold change (FC)| ≥ 0.585, P < 0.05) were found within the set distance of LD decay (r2 ≥ 0.2). Among these genes, the expression of six DEGs was verified using qRT-PCR. Zm00001d039219 and Zm00001d034319 were putatively involved in 'mitogen activated protein kinase (MAPK) signal transduction' and 'fatty acid metabolic process', respectively, based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Thus, these genes appeared to be related to low-temperature signal transduction and cell membrane fluidity.

CONCLUSION

Overall, by integrating the results of our GWAS and DEG analysis of low-temperature tolerance during germination in maize, we were able to identify a total of 30 SNPs and 82 related candidate genes, including 10 DEGs, two of which were involved in the response to tolerance to low temperature. Functional analysis will provide valuable information for understanding the genetic mechanism of low-temperature tolerance during germination in maize.

Genetic Architecture of Early Vigor Traits in Wild Soybean

A worldwide food shortage has been projected as a result of the current increase in global population and climate change. In order to provide sufficient food to feed more people, we must develop crops that can produce higher yields. Plant early vigor traits, early growth rate (EGR), early plant height (EPH), inter-node length, and node count are important traits that are related to crop yield. Glycine soja, the wild counterpart to cultivated soybean, Glycine max, harbors much higher genetic diversity and can grow in diverse environments. It can also cross easily with cultivated soybean. Thus, it holds a great potential in developing soybean cultivars with beneficial agronomic traits. In this study, we used 225 wild soybean accessions originally from diverse environments across its geographic distribution in East Asia. We quantified the natural variation of several early vigor traits, investigated the relationships among them, and dissected the genetic basis of these traits by applying a Genome-Wide Association Study (GWAS) with genome-wide single nucleotide polymorphism (SNP) data. Our results showed positive correlation between all early vigor traits studied. A total of 12 SNPs significantly associated with EPH were identified with 4 shared with EGR. We also identified two candidate genes, Glyma.07G055800.1 and Glyma.07G055900.1, playing important roles in influencing trait variation in both EGR and EPH in G. soja.

Identification of genetic loci and candidate genes related to soybean flowering through genome wide association study

BACKGROUND

As a photoperiod-sensitive and self-pollinated species, the growth periods traits play important roles in the adaptability and yield of soybean. To examine the genetic architecture of soybean growth periods, we performed a genome-wide association study (GWAS) using a panel of 278 soybean accessions and 34,710 single nucleotide polymorphisms (SNPs) with minor allele frequencies (MAF) higher than 0.04 detected by the specific-locus amplified fragment sequencing (SLAF-seq) with a 6.14-fold average sequencing depth. GWAS was conducted by a compressed mixed linear model (CMLM) involving in both relative kinship and population structure.

RESULTS

GWAS revealed that 37 significant SNP peaks associated with soybean flowering time or other growth periods related traits including full bloom, beginning pod, full pod, beginning seed, and full seed in two or more environments at -log10(P) > 3.75 or -log10(P) > 4.44 were distributed on 14 chromosomes, including chromosome 1, 2, 3, 5, 6, 9, 11, 12, 13, 14, 15, 17, 18, 19. Fourteen SNPs were novel loci and 23 SNPs were located within known QTLs or 75 kb near the known SNPs. Five candidate genes (Glyma.05G101800, Glyma.11G140100, Glyma.11G142900, Glyma.19G099700, Glyma.19G100900) in a 90 kb genomic region of each side of four significant SNPs (Gm5_27111367, Gm11_10629613, Gm11_10950924, Gm19_34768458) based on the average LD decay were homologs of Arabidopsis flowering time genes of AT5G48385.1, AT3G46510.1, AT5G59780.3, AT1G28050.1, and AT3G26790.1. These genes encoding FRI (FRIGIDA), PUB13 (plant U-box 13), MYB59, CONSTANS, and FUS3 proteins respectively might play important roles in controlling soybean growth periods.

CONCLUSIONS

This study identified putative SNP markers associated with soybean growth period traits, which could be used for the marker-assisted selection of soybean growth period traits. Furthermore, the possible candidate genes involved in the control of soybean flowering time were predicted.

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.

SNP-based mixed model association of growth- and yield-related traits in popcorn

The identification of the genes responsible for complex traits is highly promising to accelerate crop breeding, but such information is still limited for popcorn. Thus, in the present study, a mixed linear model-based association analysis (MLMA) was applied for six important popcorn traits: plant and ear height, 100-grain weight, popping expansion, grain yield and expanded popcorn volume per hectare. To this end, 196 plants of the open-pollinated popcorn population UENF-14 were sampled, selfed (S1), and then genotyped with a panel of 10,507 single nucleotide polymorphisms (SNPs) markers distributed throughout the genome. The six traits were studied under two environments [Campos dos Goytacazes-RJ (ENV1) and Itaocara-RJ (ENV2)] in an incomplete block design. Based on the phenotypic data of the S1 progenies and on the genetic characteristics of the parents, the MLMA was performed. Thereafter, genes annotated in the MaizeGDB platform were screened for potential linkage disequilibrium with the SNPs associated to the six evaluated traits. Overall, seven and eight genes were identified as associated with the traits in ENV1 and ENV2, respectively, and proteins encoded by these genes were evaluated for their function. The results obtained here contribute to increase knowledge on the genetic architecture of the six evaluated traits and might be used for marker-assisted selection in breeding programs.

Candidate genes and genome-wide association study of grain protein content and protein deviation in durum wheat

Stable QTL for grain protein content co-migrating with nitrogen-related genes have been identified by the candidate genes and genome-wide association mapping approaches useful for marker-assisted selection. Grain protein content (GPC) is one of the most important quality traits in wheat, defining the nutritional and end-use properties and rheological characteristics. Over the years, a number of breeding programs have been developed aimed to improving GPC, most of them having been prevented by the negative correlation with grain yield. To overcome this issue, a collection of durum wheat germplasm was evaluated for both GPC and grain protein deviation (GPD) in seven field trials. Fourteen candidate genes involved in several processes related to nitrogen metabolism were precisely located on two high-density consensus maps of common and durum wheat, and six of them were found to be highly associated with both traits. The wheat collection was genotyped using the 90 K iSelect array, and 11 stable quantitative trait loci (QTL) for GPC were detected in at least three environments and the mean across environments by the genome-wide association mapping. Interestingly, seven QTL were co-migrating with N-related candidate genes. Four QTL were found to be significantly associated to increases of GPD, indicating that selecting for GPC could not affect final grain yield per spike. The combined approaches of candidate genes and genome-wide association mapping led to a better understanding of the genetic relationships between grain storage proteins and grain yield per spike, and provided useful information for marker-assisted selection programs.

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.

Genome-wide association analysis of stalk biomass and anatomical traits in maize

BACKGROUND

Maize stover is an important source of crop residues and a promising sustainable energy source in the United States. Stalk is the main component of stover, representing about half of stover dry weight. Characterization of genetic determinants of stalk traits provide a foundation to optimize maize stover as a biofuel feedstock. We investigated maize natural genetic variation in genome-wide association studies (GWAS) to detect candidate genes associated with traits related to stalk biomass (stalk diameter and plant height) and stalk anatomy (rind thickness, vascular bundle density and area).

RESULTS

Using a panel of 942 diverse inbred lines, 899,784 RNA-Seq derived single nucleotide polymorphism (SNP) markers were identified. Stalk traits were measured on 800 members of the panel in replicated field trials across years. GWAS revealed 16 candidate genes associated with four stalk traits. Most of the detected candidate genes were involved in fundamental cellular functions, such as regulation of gene expression and cell cycle progression. Two of the regulatory genes (Zmm22 and an ortholog of Fpa) that were associated with plant height were previously shown to be involved in regulating the vegetative to floral transition. The association of Zmm22 with plant height was confirmed using a transgenic approach. Transgenic lines with increased expression of Zmm22 showed a significant decrease in plant height as well as tassel branch number, indicating a pleiotropic effect of Zmm22.

CONCLUSION

Substantial heritable variation was observed in the association panel for stalk traits, indicating a large potential for improving useful stalk traits in breeding programs. Genome-wide association analyses detected several candidate genes associated with multiple traits, suggesting common regulatory elements underlie various stalk traits. Results of this study provide insights into the genetic control of maize stalk anatomy and biomass.

Genetic characterization of maize germplasm derived from Suwan population and temperate resources

Background

The Suwan population is a well-known maize germplasm that has greatly contributed to the development of maize breeding in tropical and subtropical regions, especially in southern China. Inbred lines derived from the Suwan population always contain stronger resistance and extensive adaptability in different environments. To evaluate the genetic character of inbred lines derived from the Suwan population, a panel including 226 inbred line derived from the Suwan population and temperate resources was assembled and genotyped by using MaizeSNP50 BeadChip, which contained 56,110 genome-wide single nucleotide polymorphism (SNP) markers. This panel contained 98 temperate inbred line and 128 lines derived from the Suwan population.

Results

The results showed that high genetic diversity was found, with PIC and GD to be 0.67 and 0.60, respectively. In addition, two novel subgroups were identified, with representative inbred lines as HCL645 and Ki32, respectively. One acknowledged heterotic group of Iowa Stiff Stalk Synthetic (SS) was also identified in this study. This study can provide some additional scientific evidence for heterotic group division and use in maize. Additionally, lower linkage disequilibrium (LD) levels and weaker genetic relationships were found, with an average LD level of 41.15 kb that varied from 3.5 to 96 kb. A total of 82.8% of paired relative kinships ranged from 0.05 to 0.28.

Conclusions

These results would not only facilitate maize breeding practices in tropical and subtropical regions, but also revealed that this panel can be used in dissecting the genetic basis of complex quantitative traits’ variations by using genome-wide association studies (GWAS).

Electronic supplementary material

The online version of this article (10.1186/s41065-018-0077-1) contains supplementary material, which is available to authorized users.

Genome-wide association study and candidate gene analysis of rice cadmium accumulation in grain in a diverse rice collection

BACKGROUND

Cadmium (Cd) accumulation in rice followed by transfer to the food chain causes severe health problems in humans. Breeding of low Cd accumulation varieties is one of the most economical ways to solve the problem. However, information on the identity of rice germplasm with low Cd accumulation is limited, particularly in indica, and the genetic basis of Cd accumulation in rice is not well understood.

RESULTS

Screening of 312 diverse rice accessions revealed that the grain Cd concentrations of these rice accessions ranged from 0.12 to 1.23 mg/kg, with 24 accessions less than 0.20 mg/kg. Three of the 24 accessions belong to indica. Japonica accumulated significantly less Cd than indica (p < 0.001), while tropical japonica accumulated significantly less Cd than temperate japonica (p < 0.01). GWAS in all accessions identified 14 QTLs for Cd accumulation, with 7 identified in indica and 7 identified in japonica subpopulations. No common QTL was identified between indica and japonica. The previously identified genes (OsHMA3, OsNRAMP1, and OsNRAMP5) from japonica were colocalized with QTLs identified in japonica instead of indica. Expression analysis of OsNRAMP2, the candidate gene of the novel QTL (qCd3-2) identified in the present study, demonstrated that OsNRAMP2 was mainly induced in the shoots of high Cd accumulation accessions after Cd treatment. Four amino acid differences were found in the open reading frame of OsNRAMP2 between high and low Cd accumulation accessions. The allele from low Cd accumulation accessions significantly increased the Cd sensitivity and accumulation in yeast. Subcellular localization analysis demonstrated OsNRAMP2 expressed in the tonoplast of rice protoplast.

CONCLUSION

The results suggest that grain Cd concentrations are significantly different among subgroups, with Cd concentrations decreasing from indica to temperate japonica to tropical japonica. However, considerable variations exist within subgroups. The fact that no common QTL was identified between indica and japonica implies that there is a different genetic basis for determining Cd accumulation between indica and japonica, or that some QTLs for Cd accumulation in rice are subspecies-specific. Through further integrated analysis, it is speculated that OsNRAMP2 could be a novel functional gene associated with Cd accumulation in rice.

Identification of the Genomic Region Underlying Seed Weight per Plant in Soybean (Glycine max L. Merr.) via High-Throughput Single-Nucleotide Polymorphisms and a Genome-Wide Association Study

Seed weight per plant (SWPP) of soybean (Glycine max (L.) Merr.), a complicated quantitative trait controlled by multiple genes, was positively associated with soybean seed yields. In the present study, a natural soybean population containing 185 diverse accessions primarily from China was used to analyze the genetic basis of SWPP via genome-wide association analysis (GWAS) based on high-throughput single-nucleotide polymorphisms (SNPs) generated by the Specific Locus Amplified Fragment Sequencing (SLAF-seq) method. A total of 33,149 SNPs were finally identified with minor allele frequencies (MAF) > 5% which were present in 97% of all the genotypes. Twenty association signals associated with SWPP were detected via GWAS. Among these signals, eight SNPs were novel loci, and the other twelve SNPs were overlapped or located in the linked genomic regions of the reported QTL from SoyBase database. Several genes belonging to the categories of hormone pathways, RNA regulation of transcription in plant development, ubiquitin, transporting systems, and other metabolisms were considered as candidate genes associated with SWPP. Furthermore, nine genes from the flanking region of Gm07:19488264, Gm08:15768591, Gm08:15768603, or Gm18:23052511 were significantly associated with SWPP and were stable among multiple environments. Nine out of 18 haplotypes from nine genes showed the effect of increasing SWPP. The identified loci along with the beneficial alleles and candidate genes could be of great value for studying the molecular mechanisms underlying SWPP and for improving the potential seed yield of soybean in the future.

Dissecting the Genetic Basis Underlying Combining Ability of Plant Height Related Traits in Maize

Maize plant height related traits including plant height, ear height, and internode number are tightly linked with biomass, planting density, and grain yield in the field. Previous studies have focused on understanding the genetic basis of plant architecture traits per se, but the genetic basis of combining ability remains poorly understood. In this study, 328 recombinant inbred lines were inter-group crossed with two testers to produce 656 hybrids using the North Carolina II mating design. Both of the parental lines and hybrids were evaluated in two summer maize-growing regions of China in 2015 and 2016. QTL mapping highlighted that 7 out of 16 QTL detected for RILs per se could be simultaneously detected for general combining ability (GCA) effects, suggesting that GCA effects and the traits were genetically controlled by different sets of loci. Among the 35 QTL identified for hybrid performance, 57.1% and 28.5% QTL overlapped with additive/GCA and non-additive/SCA effects, suggesting that the small percentage of hybrid variance due to SCA effects in our design. Two QTL hotspots, located on chromosomes 5 and 10 and including the qPH5-1 and qPH10 loci, were validated for plant height related traits by Ye478 derivatives. Notably, the qPH5-1 locus could simultaneously affect the RILs per se and GCA effects while the qPH10, a major QTL (PVE > 10%) with pleiotropic effects, only affected the GCA effects. These results provide evidence that more attention should be focused on loci that influence combining ability directly in maize hybrid breeding.

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.

Multi-Locus Genome-Wide Association Study Reveals the Genetic Architecture of Stalk Lodging Resistance-Related Traits in Maize

Stalk lodging resistance, which is mainly measured by stem diameter (SD), stalk bending strength (SBS), and rind penetrometer resistance (RPR) in maize, seriously affects the yield and quality of maize (Zea mays L.). To dissect its genetic architecture, in this study multi-locus genome-wide association studies for stalk lodging resistance-related traits were conducted in a population of 257 inbred lines, with tropical, subtropical, and temperate backgrounds, genotyped with 48,193 high-quality single nucleotide polymorphisms. The analyses of phenotypic variations for the above traits in three environments showed high broad-sense heritability (0.679, 0.720, and 0.854, respectively). In total, 423 significant Quantitative Trait Nucleotides (QTNs) were identified by mrMLM, FASTmrEMMA, ISIS EM-BLASSO, and pLARmEB methods to be associated with the above traits. Among these QTNs, 29, 34, and 48 were commonly detected by multiple methods or across multiple environments to be related to SD, SBS, and RPR, respectively. The superior allele analyses in 30 elite lines showed that only eight lines contained more than 50% of the superior alleles, indicating that stalk lodging resistance can be improved by the integration of more superior alleles. Among sixty-three candidate genes of the consistently expressed QTNs, GRMZM5G856734 and GRMZM2G116885, encoding membrane steroid-binding protein 1 and cyclin-dependent kinase inhibitor 1, respectively, possibly inhibit cell elongation and division, which regulates lodging resistance. Our results provide the further understanding of the genetic foundation of maize lodging resistance.

Identification of minor effect QTLs for plant architecture related traits using super high density genotyping and large recombinant inbred population in maize (Zea mays)

BACKGROUND

Plant Architecture Related Traits (PATs) are of great importance for maize breeding, and mainly controlled by minor effect quantitative trait loci (QTLs). However, cloning or even fine-mapping of minor effect QTLs is very difficult in maize. Theoretically, large population and high density genetic map can be helpful for increasing QTL mapping resolution and accuracy, but such a possibility have not been actually tested.

RESULTS

Here, we employed a genotyping-by-sequencing (GBS) strategy to construct a linkage map with 16,769 marker bins for 1021 recombinant inbred lines (RILs). Accurately mapping of well studied genes P1, pl1 and r1 underlying silk color demonstrated the map quality. After QTL analysis, a total of 51 loci were mapped for six PATs. Although all of them belong to minor effect alleles, the lengths of the QTL intervals, with a minimum and median of 1.03 and 3.40 Mb respectively, were remarkably reduced as compared with previous reports using smaller size of population or small number of markers. Several genes with known function in maize were shown to be overlapping with or close neighboring to these QTL peaks, including na1, td1, d3 for plant height, ra1 for tassel branch number, and zfl2 for tassel length. To further confirm our mapping results, a plant height QTL, qPH1a, was verified by an introgression lines (ILs).

CONCLUSIONS

We demonstrated a method for high resolution mapping of minor effect QTLs in maize, and the resulted comprehensive QTLs for PATs are valuable for maize molecular breeding in the future.

Genome-wide association study discovered candidate genes of Verticillium wilt resistance in upland cotton (Gossypium hirsutum L.)

Verticillium wilt (VW), caused by infection by Verticillium dahliae, is considered one of the most yield-limiting diseases in cotton. To examine the genetic architecture of cotton VW resistance, we performed a genome-wide association study (GWAS) using a panel of 299 accessions and 85 630 single nucleotide polymorphisms (SNPs) detected using the specific-locus amplified fragment sequencing (SLAF-seq) approach. Trait-SNP association analysis detected a total of 17 significant SNPs at P < 1.17 × 10-5 (P = 1/85 630, -log10 P = 4.93); the peaks of SNPs associated with VW resistance on A10 were continuous and common in three environments (RDIG2015, RDIF2015 and RDIF2016). Haplotype block structure analysis predicted 22 candidate genes for VW resistance based on A10_99672586 with a minimum P-value (-log10 P = 6.21). One of these genes (CG02) was near the significant SNP A10_99672586 (0.26 Mb), located in a 372-kb haplotype block, and its Arabidopsis AT3G25510 homologues contain TIR-NBS-LRR domains that may be involved in disease resistance response. Real-time quantitative PCR and virus-induced gene silencing (VIGS) analysis showed that CG02 was specific to up-regulation in the resistant (R) genotype Zhongzhimian2 (ZZM2) and that silenced plants were more susceptible to V. dahliae. These results indicate that CG02 is likely the candidate gene for resistance against V. dahliae in cotton. The identified locus or gene may serve as a promising target for genetic engineering and selection for improving resistance to VW in cotton.

Brassinosteroid and gibberellin control of seedling traits in maize (Zea mays L.)

In this study, we established two doubled haploid (DH) libraries with a total of 207 DH lines. We applied BR and GA inhibitors to all DH lines at seedling stage and measured seedling BR and GA inhibitor responses. Moreover, we evaluated field traits for each DH line (untreated). We conducted genome-wide association studies (GWAS) with 62,049 genome wide SNPs to explore the genetic control of seedling traits by BR and GA. In addition, we correlate seedling stage hormone inhibitor response with field traits. Large variation for BR and GA inhibitor response and field traits was observed across these DH lines. Seedling stage BR and GA inhibitor response was significantly correlate with yield and flowering time. Using three different GWAS approaches to balance false positive/negatives, multiple SNPs were discovered to be significantly associated with BR/GA inhibitor responses with some localized within gene models. SNPs from gene model GRMZM2G013391 were associated with GA inhibitor response across all three GWAS models. This gene is expressed in roots and shoots and was shown to regulate GA signaling. These results show that BRs and GAs have a great impact for controlling seedling growth. Gene models from GWAS results could be targets for seeding traits improvement.

QTL mapping for haploid male fertility by a segregation distortion method and fine mapping of a key QTL qhmf4 in maize

Four QTL related to haploid male fertility were detected by a segregation distortion method and the key QTL qhmf4 was fine mapped to an interval of ~800 kb. Doubled haploid (DH) technology enables rapid development of homozygous lines in maize breeding programs. However, haploid genome doubling is a bottleneck for the commercialization of DH technology and is limited by haploid male fertility (HMF). This is the first study reporting the quantitative trait locus (QTL) analysis of HMF in maize. Four QTL, qhmf1, qhmf2, qhmf3, and qhmf4, controlling HMF have been identified by segregation distortion (SD) loci detection in the selected haploid population derived from 'Yu87-1/Zheng58'. Three loci, qhmf1, qhmf2, and qhmf4, were also detected in the selected haploid population derived from '4F1/Zheng58'. The QTL qhmf4 showed the strongest SD in both haploid populations. Based on the sequence information of 'Yu87-1' and 'Zheng58', thirteen markers being polymorphic between the two lines were developed to saturate the qhmf4 region. A total of 8168 H₁BC₂ (haploid backcross generation) plants produced from 'Yu87-1' and 'Zheng58' were screened for recombinants. All the 48 recombinants were backcrossed to 'Zheng58' to develop H₁BC₃ progeny. The heterozygous H₁BC₃ individuals were crossed with CAU5 to induce haploids. In each H₁BC₃ progeny, haploids were genotyped and evaluated for anther emergence score (AES). Significant (or no significant) difference (P < 0.05) between haploids with or without 'Yu87-1' donor segment indicated presence or absence of qhmf4 in the donor segment. The analysis of the 48 recombinants narrowed the qhmf4 locus down to an ~800 kb interval flanked by markers IND166 and IND1668.

Analysis of the genetic architecture of maize ear and grain morphological traits by combined linkage and association mapping

Using combined linkage and association mapping, 26 stable QTL and six stable SNPs were detected across multiple environments for eight ear and grain morphological traits in maize. One QTL, PKS2, might play an important role in maize yield improvement. In the present study, one bi-parental population and an association panel were used to identify quantitative trait loci (QTL) for eight ear and grain morphological traits. A total of 108 QTL related to these traits were detected across four environments using an ultra-high density bin map constructed using recombinant inbred lines (RILs) derived from a cross between Ye478 and Qi319, and 26 QTL were identified in more than two environments. Furthermore, 64 single nucleotide polymorphisms (SNPs) were found to be significantly associated with the eight ear and grain morphological traits (-log₁₀(P) > 4) in an association panel of 240 maize inbred lines. Combining the two mapping populations, a total of 17 pleiotropic QTL/SNPs (pQTL/SNPs) were associated with various traits across multiple environments. PKS2, a stable locus influencing kernel shape identified on chromosome 2 in a genome-wide association study (GWAS), was within the QTL confidence interval defined by the RILs. The candidate region harbored a short 13-Kb LD block encompassing four SNPs (SYN11386, PHM14783.16, SYN11392, and SYN11378). In the association panel, 13 lines derived from the hybrid PI78599 possessed the same allele as Qi319 at the PHM14783.16 (GG) locus, with an average value of 0.21 for KS, significantly lower than that of the 34 lines derived from Ye478 that carried a different allele (0.25, P < 0.05). Therefore, further fine mapping of PKS2 will provide valuable information for understanding the genetic components of grain yield and improving molecular marker-assisted selection (MAS) in 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.

Genome-wide Association Studies in Maize: Praise and Stargaze

Genome-wide association study (GWAS) has become a widely accepted strategy for decoding genotype-phenotype associations in many species thanks to advances in next-generation sequencing (NGS) technologies. Maize is an ideal crop for GWAS and significant progress has been made in the last decade. This review summarizes current GWAS efforts in maize functional genomics research and discusses future prospects in the omics era. The general goal of GWAS is to link genotypic variations to corresponding differences in phenotype using the most appropriate statistical model in a given population. The current review also presents perspectives for optimizing GWAS design and analysis. GWAS analysis of data from RNA, protein, and metabolite-based omics studies is discussed, along with new models and new population designs that will identify causes of phenotypic variation that have been hidden to date. The joint and continuous efforts of the whole community will enhance our understanding of maize quantitative traits and boost crop molecular breeding designs.

Gibberellin in plant height control: old player, new story

Height relates to plant architecture, lodging resistance, and yield performance. Growth-promoting phytohormones gibberellins (GAs) play a pivotal role in plant height control. Mutations in GA biosynthesis, metabolism, and signaling cascades influence plant height. Moreover, GA interacts with other phytohormones in the modulation of plant height. Here, we first briefly describe the regulation of plant height by altered GA pathway. Then, we depict effects of the crosstalk between GA and other phytohormones on plant height. We also dissect the co-localization of GA pathway genes and established quantitative genetic loci for plant height. Finally, we suggest ways forward for the application of hormone GA knowledge in breeding of crops with plant height ideotypes.

Gibberellins Promote Brassinosteroids Action and Both Increase Heterosis for Plant Height in Maize (Zea mays L.)

Brassinosteroids (BRs) and Gibberellins (GAs) are two classes of plant hormones affecting plant height (PHT). Thus, manipulation of BR and GA levels or signaling enables optimization of crop grain and biomass yields. We established backcross (BC) families, selected for increased PHT, in two elite maize inbred backgrounds. Various exotic accessions used in the germplasm enhancement in maize project served as donors. BC1-derived doubled haploid lines in the same two elite maize inbred backgrounds established without selection for plant height were included for comparison. We conducted genome-wide association studies to explore the genetic control of PHT by BR and GA. In addition, we used BR and GA inhibitors to compare the relationship between PHT, BR, and GA in inbred lines and heterozygotes from a physiological and biological perspective. A total of 73 genomic loci were discovered to be associated with PHT, with seven co-localized with GA, and two co-localized with BR candidate genes. PHT determined in field trials was significantly correlated with seedling stage BR and GA inhibitor responses. However, this observation was only true for maize heterozygotes, not for inbred lines. Path analysis results suggest that heterozygosity increases GA levels, which in turn promote BR levels. Thus, at least part of heterosis for PHT in maize can be explained by increased GA and BR levels, and seedling stage hormone inhibitor response is promising to predict heterosis for PHT.

Dissection of Recombination Attributes for Multiple Maize Populations Using a Common SNP Assay

Recombination is a vital characteristic for quantitative trait loci mapping and breeding to enhance the yield potential of maize. However, recombination characteristics in globally used segregating populations have never been evaluated at similar genetic marker densities. This study aimed to divulge the characteristics of recombination events, recombinant chromosomal segments, and recombination frequency for four dissimilar populations. These populations were doubled haploid (DH), recombination inbred line (RIL), intermated B73xMo17 (IBM), and multi-parent advanced generation inter-cross (MAGIC), using the Illumina MaizeSNP50 BeadChip to provide markers. Our results revealed that the average number of recombination events was 16, 41, 72, and 86 per line in DH, RIL, IBM, and MAGIC populations, respectively. Accordingly, the average length of recombinant chromosomal segments was 84.8, 47.3, 29.2, and 20.4 Mb in DH, RIL, IBM, and MAGIC populations, respectively. Furtherly, the recombination frequency varied in different genomic regions and population types [DH (0-12.7 cM/Mb), RIL (0-15.5 cM/Mb), IBM (0-24.1 cM/Mb), MAGIC (0-42.3 cM/Mb)]. Utilizing different sub-sets of lines, the recombination bin number and size were analyzed in each population. Additionally, different sub-sets of markers and lines were employed to estimate the recombination bin number and size via formulas for relationship in these populations. The relationship between recombination events and recombination bin length was also examined. Our results contribute to determining the most suitable number of genetic markers, lines in each population, and population type for successful mapping and breeding.

Genome-Wide Association Study and QTL Mapping Reveal Genomic Loci Associated with Fusarium Ear Rot Resistance in Tropical Maize Germplasm

Fusarium ear rot (FER) incited by Fusarium verticillioides is a major disease of maize that reduces grain quality globally. Host resistance is the most suitable strategy for managing the disease. We report the results of genome-wide association study (GWAS) to detect alleles associated with increased resistance to FER in a set of 818 tropical maize inbred lines evaluated in three environments. Association tests performed using 43,424 single-nucleotide polymorphic (SNPs) markers identified 45 SNPs and 15 haplotypes that were significantly associated with FER resistance. Each associated SNP locus had relatively small additive effects on disease resistance and accounted for 1-4% of trait variation. These SNPs and haplotypes were located within or adjacent to 38 candidate genes, 21 of which were candidate genes associated with plant tolerance to stresses, including disease resistance. Linkage mapping in four biparental populations to validate GWAS results identified 15 quantitative trait loci (QTL) associated with F. verticillioides resistance. Integration of GWAS and QTL to the maize physical map showed eight colocated loci on chromosomes 2, 3, 4, 5, 9, and 10. QTL on chromosomes 2 and 9 are new. These results reveal that FER resistance is a complex trait that is conditioned by multiple genes with minor effects. The value of selection on identified markers for improving FER resistance is limited; rather, selection to combine small effect resistance alleles combined with genomic selection for polygenic background for both the target and general adaptation traits might be fruitful for increasing FER resistance in maize.

Bandwagons I, too, have known

Bandwagons come in waves. A plant breeder, just like a surfer, needs to carefully choose which waves to be on. A bandwagon is an idea, activity, or cause that becomes increasingly fashionable as more and more people adopt it. In a 1991 article entitled Bandwagons I Have Known, Professor N. W. Simmonds described several bandwagons that he encountered in his career, beginning with induced polyploidy and mutation breeding and ending with the then-new field of biotechnology. This article reviews and speculates about post-1990 bandwagons in plant improvement, including transgenic cultivars, quantitative trait locus (QTL) mapping, association mapping, genomewide (or genomic) selection, phenomics, envirotyping, and genome editing. The life cycle of a bandwagon includes an excitement phase of hype and funding; a realization phase when the initial hype is either tempered or the initial expectations are found to have been too low; and a reality phase when the useful aspects of a bandwagon become part of mainstream thinking and practice, or when an unsuccessful bandwagon is largely abandoned. During the realization phase, a new bandwagon that draws our attention and gives us renewed optimism typically arises. The most popular bandwagons, such as QTL mapping, are those for which the needed experimental resources are accessible, the required technical knowledge and skills can be easily learned, and the outputs can almost always be reported. The favorite bandwagon of any plant breeder has, in one way or another, resulted from Mendel's seminal discoveries 150 years ago. Our community of plant breeders needs to be continually diligent in welcoming new bandwagons, but also in hopping off from those that do not prove useful.

Combined linkage and association mapping identifies a major QTL (qRtsc8-1), conferring tar spot complex resistance in maize

A major QTL ( qRtsc8 - 1 ) conditioning resistance to tar spot complex of maize and occurring at a frequency of 3.5 % across 890 maize inbred lines. Tar spot complex (TSC) is a highly destructive disease of maize found in some countries in America. Identification of TSC resistant germplasm and elucidating the genetic mechanism of resistance is crucial for the use of host resistance to manage this disease. We evaluated 890 elite maize inbred lines in multiple environments and used genome wide association analysis (GWAS) with genotypic data from Illumina MaizeSNP50 BeadChip containing 56 K SNPs to dissect the genetics of TSC resistance. GWAS results were validated through linkage analysis in three bi-parental populations derived from different resistant and susceptible parents. Through GWAS, three TSC resistance loci were identified on chromosome 2, 7 and 8 (-log10 (p) > 5.99). A major quantitative resistance locus (QTL) designated qRtsc8-1, was detected on maize chromosome bin 8.03. qRtsc8-1, was confirmed in three independent bi-parental populations and it accounted for 18-43 % of the observed phenotypic variation for TSC. A rare haplotype within the qRtsc8-1 region, occurring at a frequency of 3.5 % increased TSC resistance by 14 %. Candidate gene analysis revealed that a leucine-rich repeat receptor-like protein (LRR-RLKs) gene family maybe the candidate gene for qRtsc8-1. Identification and localization of a major locus conditioning TSC resistance provides the foundation for fine mapping qRtsc8-1 and developing functional markers for improving TSC resistance in maize breeding programs. To the best of our knowledge, this is the first report of a major QTL for TSC resistance.

Genetic dissection of maize plant architecture with an ultra-high density bin map based on recombinant inbred lines

Background

Plant architecture attributes, such as plant height, ear height, and internode number, have played an important role in the historical increases in grain yield, lodging resistance, and biomass in maize (Zea mays L.). Analyzing the genetic basis of variation in plant architecture using high density QTL mapping will be of benefit for the breeding of maize for many traits. However, the low density of molecular markers in existing genetic maps has limited the efficiency and accuracy of QTL mapping. Genotyping by sequencing (GBS) is an improved strategy for addressing a complex genome via next-generation sequencing technology. GBS has been a powerful tool for SNP discovery and high-density genetic map construction. The creation of ultra-high density genetic maps using large populations of advanced recombinant inbred lines (RILs) is an efficient way to identify QTL for complex agronomic traits.

Results

A set of 314 RILs derived from inbreds Ye478 and Qi319 were generated and subjected to GBS. A total of 137,699,000 reads with an average of 357,376 reads per individual RIL were generated, which is equivalent to approximately 0.07-fold coverage of the maize B73 RefGen_V3 genome for each individual RIL. A high-density genetic map was constructed using 4183 bin markers (100-Kb intervals with no recombination events). The total genetic distance covered by the linkage map was 1545.65 cM and the average distance between adjacent markers was 0.37 cM with a physical distance of about 0.51 Mb. Our results demonstrated a relatively high degree of collinearity between the genetic map and the B73 reference genome. The quality and accuracy of the bin map for QTL detection was verified by the mapping of a known gene, pericarp color 1 (P1), which controls the color of the cob, with a high LOD value of 80.78 on chromosome 1. Using this high-density bin map, 35 QTL affecting plant architecture, including 14 for plant height, 14 for ear height, and seven for internode number were detected across three environments. Interestingly, pQTL10, which influences all three of these traits, was stably detected in three environments on chromosome 10 within an interval of 14.6 Mb. Two MYB transcription factor genes, GRMZM2G325907 and GRMZM2G108892, which might regulate plant cell wall metabolism are the candidate genes for qPH10.

Conclusions

Here, an ultra-high density accurate linkage map for a set of maize RILs was constructed using a GBS strategy. This map will facilitate identification of genes and exploration of QTL for plant architecture in maize. It will also be helpful for further research into the mechanisms that control plant architecture while also providing a basis for marker-assisted selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2555-z) contains supplementary material, which is available to authorized users.

A comprehensive meta-analysis of plant morphology, yield, stay-green, and virus disease resistance QTL in maize (Zea mays L.)

The meta-QTL and candidate genes will facilitate the elucidation of molecular bases underlying agriculturally important traits and open new avenues for functional markers development and elite alleles introgression in maize breeding program. A large number of QTLs attributed to grain productivity and other agriculturally important traits have been identified and deposited in public repositories. The integration of fruitful QTL becomes a major issue in current plant genomics. To this end, we first collected QTL for six agriculturally important traits in maize, including yield, plant height, ear height, leaf angle, stay-green, and maize rough dwarf disease resistance. The meta-analysis method was then employed to retrieve 113 meta-QTL. Additionally, we also isolated candidate genes for target traits by the bioinformatic technique. Several candidates, including some well-characterized genes, GA3ox2 for plant height, lg1 and lg4 for leaf angle, zfl1 and zfl2 for flowering time, were co-localized with established meta-QTL intervals. Intriguingly, in a relatively narrow meta-QTL region, the maize ortholog of rice yield-related gene GW8/OsSPL16 was believed to be a candidate for yield. Leveraging results presented in this study will provide further insights into the genetic architecture of maize agronomic traits. Moreover, the meta-QTL and candidate genes reported here could be harnessed for the enhancement of stress tolerance and yield performance in maize and translation to other crops.

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

Plant height (PH) and ear height (EH) are two very important agronomic traits related to the population density and lodging in maize. In order to better understand of the genetic basis of nature variation in PH and EH, two bi-parental populations and one genome-wide association study (GWAS) population were used to map quantitative trait loci (QTL) for both traits. Phenotypic data analysis revealed a wide normal distribution and high heritability for PH and EH in the three populations, which indicated that maize height is a highly polygenic trait. A total of 21 QTL for PH and EH in three common genomic regions (bin 1.05, 5.04/05, and 6.04/05) were identified by QTL mapping in the two bi-parental populations under multiple environments. Additionally, 41 single nucleotide polymorphisms (SNPs) were identified for PH and EH by GWAS, of which 29 SNPs were located in 19 unique candidate gene regions. Most of the candidate genes were related to plant growth and development. One QTL on Chromosome 1 was further verified in a near-isogenic line (NIL) population, and GWAS identified a C2H2 zinc finger family protein that maybe the candidate gene for this QTL. These results revealed that nature variation of PH and EH are strongly controlled by multiple genes with low effect and facilitated a better understanding of the underlying mechanism of height in maize.

A genome-wide association study of plant height and primary branch number in rapeseed (Brassica napus)

Crop plant architecture plays a highly important role in its agronomic performance. Plant height (PH) and primary branch number (PB) are two major factors that affect the plant architecture of rapeseed (Brassica napus). Previous studies have shown that these two traits are controlled by multiple quantitative trait loci (QTL); however, QTLs have not been delimited to regions less than 10cM. Genome-wide association study (GWAS) is a highly efficient approach for identifying genetic loci controlling traits at relatively high resolution. In this study, variations in PH and PB of a panel of 472 rapeseed accessions that had previously been analyzed by a 60k SNP array were investigated for three consecutive years and studied by GWAS. Eight QTLs on chromosome A03, A05, A07 and C07 were identified for PH, and five QTLs on A01, A03, A07 and C07 were identified for PB. Although most QTLs have been detected in previous studies based on linkage analyses, the two QTLs of PH on A05 and the QTL of PB on C07 were novel. In the genomic regions close to the GWAS peaks, orthologs of the genes involved in flower development, phytohormone biosynthesis, metabolism and signaling in Arabidopsis were identified.

Genetic Analysis in Maize Foundation Parents with Mapping Population and Testcross Population: Ye478 Carried More Favorable Alleles and Using QTL Information Could Improve Foundation Parents

The development of maize foundation parents is an important part of genetics and breeding research, and applying new genetic information to produce foundation parents has been challenging. In this study, we focused on quantitative trait loci (QTLs) and general combining ability (GCA) of Ye478, a widely used foundation parent in China. We developed three sets of populations for QTL mapping and to analyze the GCA for some agronomic traits. The assessment of 15 traits resulted in the detection of 251 QTLs in six tested environments, with 119 QTLs identified through a joint analysis across all environments. Further, analyses revealed that most favorable alleles for plant type-related traits were from Ye478, and more than half of the favorable alleles for yield-related traits were from R08, another foundation parent used in southwestern China, suggesting that different types of foundation parents carried different favorable alleles. We observed that the GCA for most traits (e.g., plant height and 100-kernel weight) was maintained in the inbred lines descended from the foundation parents. Additionally, the continuous improvement in the GCA of the descendants of the foundation parents was consistent with the main trend in maize breeding programs. We identified three significant genomic regions that were highly conserved in three Ye478 descendants, including the stable QTL for plant height. The GCA for the traits in the F₇ generation revealed that the QTLs for the given traits per se were affected by additive effects in the same way in different populations.

Genome-wide association of drought-related and biomass traits with HapMap SNPs in Medicago truncatula

Improving drought tolerance of crop plants is a major goal of plant breeders. In this study, we characterized biomass and drought-related traits of 220 Medicago truncatula HapMap accessions. Characterized traits included shoot biomass, maximum leaf size, specific leaf weight, stomatal density, trichome density and shoot carbon-13 isotope discrimination (δ(13) C) of well-watered M. truncatula plants, and leaf performance in vitro under dehydration stress. Genome-wide association analyses were carried out using the general linear model (GLM), the standard mixed linear model (MLM) and compressed MLM (CMLM) in TASSEL, which revealed significant overestimation of P-values by CMLM. For each trait, candidate genes and chromosome regions containing SNP markers were found that are in significant association with the trait. For plant biomass, a 0.5 Mbp region on chromosome 2 harbouring a plasma membrane intrinsic protein, PIP2, was discovered that could potentially be targeted to increase dry matter yield. A protein disulfide isomerase-like protein was found to be tightly associated with both shoot biomass and leaf size. A glutamate-cysteine ligase and an aldehyde dehydrogenase family protein with Arabidopsis homologs strongly expressed in the guard cells were two of the top genes identified by stomata density genome-wide association studies analysis.

Association mapping in sunflower (Helianthus annuus L.) reveals independent control of apical vs. basal branching

BACKGROUND

Shoot branching is an important determinant of plant architecture and influences various aspects of growth and development. Selection on branching has also played an important role in the domestication of crop plants, including sunflower (Helianthus annuus L.). Here, we describe an investigation of the genetic basis of variation in branching in sunflower via association mapping in a diverse collection of cultivated sunflower lines.

RESULTS

Detailed phenotypic analyses revealed extensive variation in the extent and type of branching within the focal population. After correcting for population structure and kinship, association analyses were performed using a genome-wide collection of SNPs to identify genomic regions that influence a variety of branching-related traits. This work resulted in the identification of multiple previously unidentified genomic regions that contribute to variation in branching. Genomic regions that were associated with apical and mid-apical branching were generally distinct from those associated with basal and mid-basal branching. Homologs of known branching genes from other study systems (i.e., Arabidopsis, rice, pea, and petunia) were also identified from the draft assembly of the sunflower genome and their map positions were compared to those of associations identified herein. Numerous candidate branching genes were found to map in close proximity to significant branching associations.

CONCLUSIONS

In sunflower, variation in branching is genetically complex and overall branching patterns (i.e., apical vs. basal) were found to be influenced by distinct genomic regions. Moreover, numerous candidate branching genes mapped in close proximity to significant branching associations. Although the sunflower genome exhibits localized islands of elevated linkage disequilibrium (LD), these non-random associations are known to decay rapidly elsewhere. The subset of candidate genes that co-localized with significant associations in regions of low LD represents the most promising target for future functional analyses.

Genome-wide association analysis of seedling root development in maize (Zea mays L.)

Background

Plants rely on the root system for anchorage to the ground and the acquisition and absorption of nutrients critical to sustaining productivity. A genome wide association analysis enables one to analyze allelic diversity of complex traits and identify superior alleles. 384 inbred lines from the Ames panel were genotyped with 681,257 single nucleotide polymorphism markers using Genotyping-by-Sequencing technology and 22 seedling root architecture traits were phenotyped.

Results

Utilizing both a general linear model and mixed linear model, a GWAS study was conducted identifying 268 marker trait associations (p ≤ 5.3×10^-7). Analysis of significant SNP markers for multiple traits showed that several were located within gene models with some SNP markers localized within regions of previously identified root quantitative trait loci. Gene model GRMZM2G153722 located on chromosome 4 contained nine significant markers. This predicted gene is expressed in roots and shoots.

Conclusion

This study identifies putatively associated SNP markers associated with root traits at the seedling stage. Some SNPs were located within or near (<1 kb) gene models. These gene models identify possible candidate genes involved in root development at the seedling stage. These and respective linked or functional markers could be targets for breeders for marker assisted selection of seedling root traits.

Electronic supplementary material

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

Genome-wide association analysis of seedling root development in maize (Zea mays L.)

BACKGROUND

Plants rely on the root system for anchorage to the ground and the acquisition and absorption of nutrients critical to sustaining productivity. A genome wide association analysis enables one to analyze allelic diversity of complex traits and identify superior alleles. 384 inbred lines from the Ames panel were genotyped with 681,257 single nucleotide polymorphism markers using Genotyping-by-Sequencing technology and 22 seedling root architecture traits were phenotyped.

RESULTS

Utilizing both a general linear model and mixed linear model, a GWAS study was conducted identifying 268 marker trait associations (p ≤ 5.3×10(-7)). Analysis of significant SNP markers for multiple traits showed that several were located within gene models with some SNP markers localized within regions of previously identified root quantitative trait loci. Gene model GRMZM2G153722 located on chromosome 4 contained nine significant markers. This predicted gene is expressed in roots and shoots.

CONCLUSION

This study identifies putatively associated SNP markers associated with root traits at the seedling stage. Some SNPs were located within or near (<1 kb) gene models. These gene models identify possible candidate genes involved in root development at the seedling stage. These and respective linked or functional markers could be targets for breeders for marker assisted selection of seedling root traits.