ABSTRACT
Durable host plant resistance (HPR) to insect pests is critical for sustainable agriculture. Natural variation exists for aphid HPR in sorghum (Sorghum bicolor), but the genetic architecture and phenotype have not been clarified and characterized for most sources. In order to assess the current threat of a sorghum aphid (Melanaphis sorghi) biotype shift, we characterized the phenotype of Resistance to Melanaphis sorghi 1 (RMES1) and additional HPR architecture in globally admixed populations selected under severe sorghum aphid infestation in Haiti. We found RMES1 reduces sorghum aphid fecundity but not bird cherry-oat aphid (Rhopalosiphum padi) fecundity, suggesting a discriminant HPR response typical of gene-for-gene interaction. A second resistant gene, Resistance to Melanaphis sorghi 2 (RMES2), was more frequent than RMES1 resistant alleles in landraces and historic breeding lines. RMES2 contributes early and mid-season aphid resistance in a segregating F2 population; however, RMES1 was only significant with mid-season fitness. In a fixed population with high sorghum aphid resistance, RMES1 and RMES2 were selected for demonstrating a lack of severe antagonistic pleiotropy. Associations with resistance colocated with cyanogenic glucoside biosynthesis genes support additional HPR sources. Globally, therefore, an HPR source vulnerable to biotype shift via selection pressure (RMES1) is bolstered by a second common source of resistance in breeding programs (RMES2), which may be staving off a biotype shift and is critical for sustainable sorghum production.
Subject(s)
Aphids , Sorghum , Aphids/physiology , Aphids/genetics , Sorghum/genetics , Sorghum/parasitology , Animals , Genes, Plant , Phenotype , Plant Diseases/genetics , Plant Diseases/parasitologyABSTRACT
Rapid environmental change can lead to population extinction or evolutionary rescue. The global staple crop sorghum (Sorghum bicolor) has recently been threatened by a global outbreak of an aggressive new biotype of sugarcane aphid (SCA; Melanaphis sacchari). We characterized genomic signatures of adaptation in a Haitian breeding population that had rapidly adapted to SCA infestation, conducting evolutionary population genomics analyses on 296 Haitian lines versus 767 global accessions. Genome scans and geographic analyses suggest that SCA adaptation has been conferred by a globally rare East African allele of RMES1, which spread to breeding programs in Africa, Asia, and the Americas. De novo genome sequencing revealed potential causative variants at RMES1. Markers developed from the RMES1 sweep predicted resistance in eight independent commercial and public breeding programs. These findings demonstrate the value of evolutionary genomics to develop adaptive trait technology and highlight the benefits of global germplasm exchange to facilitate evolutionary rescue.
ABSTRACT
Successful management and utilization of increasingly large genomic datasets is essential for breeding programs to accelerate cultivar development. To help with this, we developed a Sorghum bicolor Practical Haplotype Graph (PHG) pangenome database that stores haplotypes and variant information. We developed two PHGs in sorghum that were used to identify genome-wide variants for 24 founders of the Chibas sorghum breeding program from 0.01x sequence coverage. The PHG called single nucleotide polymorphisms (SNPs) with 5.9% error at 0.01x coverage-only 3% higher than PHG error when calling SNPs from 8x coverage sequence. Additionally, 207 progenies from the Chibas genomic selection (GS) training population were sequenced and processed through the PHG. Missing genotypes were imputed from PHG parental haplotypes and used for genomic prediction. Mean prediction accuracies with PHG SNP calls range from .57-.73 and are similar to prediction accuracies obtained with genotyping-by-sequencing or targeted amplicon sequencing (rhAmpSeq) markers. This study demonstrates the use of a sorghum PHG to impute SNPs from low-coverage sequence data and shows that the PHG can unify genotype calls across multiple sequencing platforms. By reducing input sequence requirements, the PHG can decrease the cost of genotyping, make GS more feasible, and facilitate larger breeding populations. Our results demonstrate that the PHG is a useful research and breeding tool that maintains variant information from a diverse group of taxa, stores sequence data in a condensed but readily accessible format, unifies genotypes across genotyping platforms, and provides a cost-effective option for genomic selection.
Subject(s)
Sorghum , Cost-Benefit Analysis , Genome , Genomics , Haplotypes , Sorghum/geneticsABSTRACT
Young breeding programs in developing countries, like the Chibas sorghum breeding program in Haiti, face the challenge of increasing genetic gain with limited resources. Implementing genomic selection (GS) could increase genetic gain, but optimization of GS is needed to account for these programs' unique challenges and advantages. Here, we used simulations to identify conditions under which genomic-assisted recurrent selection (GARS) would be more effective than phenotypic recurrent selection (PRS) in small new breeding programs. We compared genetic gain, cost per unit gain, genetic variance, and prediction accuracy of GARS (two or three cycles per year) vs. PRS (one cycle per year) assuming various breeding population sizes and trait genetic architectures. For oligogenic architecture, the maximum relative genetic gain advantage of GARS over PRS was 12-88%, which was observed only during the first few cycles. For the polygenic architecture, GARS provided maximum relative genetic gain advantage of 26-165%, and was always superior to PRS. Average prediction accuracy declines substantially after several cycles of selection, suggesting the prediction models should be updated regularly. Updating prediction models every year increased the genetic gain by up to 33-39% compared to no-update scenarios. For small populations and oligogenic traits, cost per unit gain was lower in PRS than GARS. However, with larger populations and polygenic traits cost per unit gain was up to 67% lower in GARS than PRS. Collectively, the simulations suggest that GARS could increase the genetic gain in small young breeding programs by accelerating the breeding cycles and enabling evaluation of larger populations.
Subject(s)
Genome-Wide Association Study/methods , Models, Genetic , Plant Breeding/methods , Selection, Genetic , Sorghum/genetics , Computer Simulation , Haiti , Multifactorial Inheritance , Selective BreedingABSTRACT
Flowering time is a complex trait that controls adaptation of plants to their local environment in the outcrossing species Zea mays (maize). We dissected variation for flowering time with a set of 5000 recombinant inbred lines (maize Nested Association Mapping population, NAM). Nearly a million plants were assayed in eight environments but showed no evidence for any single large-effect quantitative trait loci (QTLs). Instead, we identified evidence for numerous small-effect QTLs shared among families; however, allelic effects differ across founder lines. We identified no individual QTLs at which allelic effects are determined by geographic origin or large effects for epistasis or environmental interactions. Thus, a simple additive model accurately predicts flowering time for maize, in contrast to the genetic architecture observed in the selfing plant species rice and Arabidopsis.
Subject(s)
Flowers/genetics , Quantitative Trait Loci , Zea mays/genetics , Alleles , Chromosome Mapping , Chromosomes, Plant/genetics , Epistasis, Genetic , Flowers/growth & development , Gene Frequency , Genes, Plant , Genetic Variation , Geography , Inbreeding , Phenotype , Polymorphism, Single Nucleotide , Quantitative Trait, Heritable , Recombination, Genetic , Time Factors , Zea mays/growth & development , Zea mays/physiologyABSTRACT
Maize genetic diversity has been used to understand the molecular basis of phenotypic variation and to improve agricultural efficiency and sustainability. We crossed 25 diverse inbred maize lines to the B73 reference line, capturing a total of 136,000 recombination events. Variation for recombination frequencies was observed among families, influenced by local (cis) genetic variation. We identified evidence for numerous minor single-locus effects but little two-locus linkage disequilibrium or segregation distortion, which indicated a limited role for genes with large effects and epistatic interactions on fitness. We observed excess residual heterozygosity in pericentromeric regions, which suggested that selection in inbred lines has been less efficient in these regions because of reduced recombination frequency. This implies that pericentromeric regions may contribute disproportionally to heterosis.
Subject(s)
Chromosome Mapping , Chromosomes, Plant/genetics , Genetic Variation , Quantitative Trait, Heritable , Zea mays/genetics , Alleles , Centromere/genetics , Crosses, Genetic , Epistasis, Genetic , Flowers/genetics , Flowers/growth & development , Genome, Plant , Heterozygote , Hybrid Vigor , Inbreeding , Linkage Disequilibrium , Phenotype , Polymorphism, Single Nucleotide , Quantitative Trait Loci , Recombination, Genetic , Selection, Genetic , Zea mays/classification , Zea mays/physiologyABSTRACT
We characterized allelic variation at barren inflorescence2 (bif2), a maize co-ortholog of the Arabidopsis PINOID protein kinase (PID), and tested for trait associations with bif2 in both an association mapping population of 277 diverse maize inbreds and in the inter-mated B73 x Mo17 (IBM) linkage population. Results from the quantitative analyses were compared with previous reports of bif2 phenotypes in mutagenesis studies. All three approaches (association, linkage, and mutagenesis) detect a significant effect of bif2 on tassel architecture. Association mapping implicates bif2 in an unexpectedly wide range of traits including plant height, node number, leaf length, and flowering time. Linkage mapping finds a significant interaction effect for node number between bif2 and other loci, in keeping with previous reports that bif2;spi1 and Bif2;Bif1 double mutants produce fewer phytomers. The Mo17 allele is associated with a reduced tassel branch zone and shows lower expression than the B73 allele in hybrid B73-Mo17 F(1) inflorescences, consistent with the complete absence of tassel branches in the bif2 knockout mutant. Overall, these data suggest that allelic variation at bif2 affects maize architecture by modulating auxin transport during vegetative and inflorescence development.
Subject(s)
Flowers/growth & development , Plant Proteins/metabolism , Protein Serine-Threonine Kinases/metabolism , Zea mays/genetics , Alleles , Arabidopsis Proteins , Base Sequence , DNA, Plant/genetics , Flowers/genetics , Flowers/metabolism , Gene Expression Regulation, Plant , Genes, Plant , Genetic Variation , Indoleacetic Acids/metabolism , Linkage Disequilibrium , Molecular Sequence Data , Phenotype , Plant Proteins/genetics , Protein Serine-Threonine Kinases/genetics , Quantitative Trait, Heritable , Sequence Alignment , Sequence Analysis, DNA , Zea mays/enzymology , Zea mays/growth & developmentABSTRACT
As population structure can result in spurious associations, it has constrained the use of association studies in human and plant genetics. Association mapping, however, holds great promise if true signals of functional association can be separated from the vast number of false signals generated by population structure. We have developed a unified mixed-model approach to account for multiple levels of relatedness simultaneously as detected by random genetic markers. We applied this new approach to two samples: a family-based sample of 14 human families, for quantitative gene expression dissection, and a sample of 277 diverse maize inbred lines with complex familial relationships and population structure, for quantitative trait dissection. Our method demonstrates improved control of both type I and type II error rates over other methods. As this new method crosses the boundary between family-based and structured association samples, it provides a powerful complement to currently available methods for association mapping.
Subject(s)
Genetic Techniques , Heredity/genetics , Models, Genetic , Zea mays/genetics , Gene Expression , Genetic Variation , Humans , Phenotype , Quantitative Trait, Heritable , Research DesignABSTRACT
Crop improvement and the dissection of complex genetic traits require germplasm diversity. Although this necessary phenotypic variability exists in diverse maize, most research is conducted using a small subset of inbred lines. An association population of 302 lines is now available--a valuable research tool that captures a large proportion of the alleles in cultivated maize. Provided that appropriate statistical models correcting for population structure are included, this tool can be used in association analyses to provide high-resolution evaluation of multiple alleles. This study describes the population structure of the 302 lines, and investigates the relationship between population structure and various measures of phenotypic and breeding value. On average, our estimates of population structure account for 9.3% of phenotypic variation, roughly equivalent to a major quantitative trait locus (QTL), with a high of 35%. Inclusion of population structure in association models is critical to meaningful analyses. This new association population has the potential to identify QTL with small effects, which will aid in dissecting complex traits and in planning future projects to exploit the rich diversity present in maize.