Uncovering the Genetic Architecture of Seed Weight and Size in Intermediate Wheatgrass through Linkage and Association Mapping

Origin of current intermediate wheatgrass germplasm being developed for Kernza grain production

Intermediate wheatgrass (IWG, Thinopyrum intermedium [Host] Barkworth & D. R. Dewey) has been developed as a perennial grain crop to provide ecosystem services, environmental benefits, and human food. Grain and products derived from IWG varieties improved for food production have been marketed under the registered trademark, Kernza. In the 1980s, a joint breeding effort between the Rodale Institute (RI) and the Big Flats Plant Material Center used IWG plant introductions (PI) from the National Plant Germplasm System (NPGS) and recurrent phenotypic selection to improve populations of IWG with the goal of developing a perennial grain. Initial selections were provided to The Land Institute where they were subsequently improved for grain production, yet the identity of the founder material of improved, food-grade IWG has not been publicly documented. Recently recovered original documents have been used to reconstruct the early breeding program to identify the most likely 20 PIs that form the founders of modern food-grade IWG. Molecular data using genotyping-by-sequencing in current elite breeding material, remnant seed from the initial RI selections, and preserved sample material have provided supporting evidence for the historical records. The genetic origin for food-grade IWG is focused between the Black Sea and Caspian Sea in the Stavropol region of Russia, with smaller contributions likely from collections as distant as Kazakhstan in the east to Turkey in the west. This work connects the flow of germplasm and utility of NPGS PIs to present day IWG grain cultivars being developed in multiple breeding programs around the world.

Identification of quantitative trait loci (QTL) and meta-QTL analysis for kernel size-related traits in wheat (Triticum aestivum L.)

BACKGROUND

Kernel size-related traits, including kernel length (KL), kernel width (KW), kernel diameter ratio (KDR) and kernel thickness (KT), are critical determinants for wheat kernel weight and yield and highly governed by a type of quantitative genetic basis. Genome-wide identification of major and stable quantitative trait loci (QTLs) and functional genes are urgently required for genetic improvement in wheat kernel yield. A hexaploid wheat population consisting of 120 recombinant inbred lines was developed to identify QTLs for kernel size-related traits under different water environments. The meta-analysis and transcriptome evaluation were further integrated to identify major genomic regions and putative candidate genes.

RESULTS

The analysis of variance (ANOVA) revealed more significant genotypic effects for kernel size-related traits, indicating the moderate to high heritability of 0.61-0.89. Thirty-two QTLs for kernel size-related traits were identified, explaining 3.06%-14.2% of the phenotypic variation. Eleven stable QTLs were detected in more than three water environments. The 1103 original QTLs from the 34 previous studies and the present study were employed for the MQTL analysis and refined into 58 MQTLs. The average confidence interval of the MQTLs was 3.26-fold less than that of the original QTLs. The 1864 putative candidate genes were mined within the regions of 12 core MQTLs, where 70 candidate genes were highly expressed in spikes and kernels by comprehensive analysis of wheat transcriptome data. They were involved in various metabolic pathways, such as carbon fixation in photosynthetic organisms, carbon metabolism, mRNA surveillance pathway, RNA transport and biosynthesis of secondary metabolites.

CONCLUSIONS

Major genomic regions and putative candidate genes for kernel size-related traits in wheat have been revealed by an integrative strategy with QTL linkage mapping, meta-analysis and transcriptomic assessment. The findings provide a novel insight into understanding the genetic determinants of kernel size-related traits and will be useful for the marker-assisted selection of high yield in wheat breeding.

The use of wheatgrass (<i>Thinopyrum intermedium</i>) in breeding

Wheatgrass (Th. intermedium) has been traditionally used in wheat breeding for obtaining wheat-wheatgrass hybrids and varieties with introgressions of new genes for economically valuable traits. However, in the 1980s in the United States wheatgrass was selected from among perennial plant species as having promise for domestication and the development of dual-purpose varieties for grain (as an alternative to perennial wheat) and hay. The result of this work was the creation of the wheatgrass varieties Kernza (The Land Institute, Kansas) and MN-Clearwater (University of Minnesota, Minnesota). In Omsk State Agrarian University, the variety Sova was developed by mass selection of the most winter-hardy biotypes with their subsequent combination from the population of wheatgrass obtained from The Land Institute. The average grain yield of the variety Sova is 9.2 dt/ha, green mass is 210.0 dt/ ha, and hay is 71.0 dt/ha. Wheatgrass is a crop with a large production potential, benef icial environmental properties, and valuable grain for functional food. Many publications show the advantages of growing the Kernza variety compared to annual crops in reducing groundwater nitrate contamination, increasing soil carbon sequestration, and reducing energy and economic costs. However, breeding programs for domestication of perennial crops are very limited in Russia. This paper presents an overview of main tasks faced by breeders, aimed at enhancing the yield and cultivating wheatgrass eff iciency as a perennial grain and fodder crop. To address them, both traditional and modern biotechnological and molecular cytogenetic approaches are used. The most important task is to transfer target genes of Th. intermedium to modern wheat varieties and decrease the level of chromatin carrying undesirable genes of the wild relative. The f irst consensus map of wheatgrass containing 10,029 markers was obtained, which is important for searching for genes and their introgressions to the wheat genome. The results of research on the nutritional and technological properties of wheatgrass grain for the development of food products as well as the differences in the quality of wheatgrass grain and wheat grain are presented.

Genetic architecture and QTL selection response for Kernza perennial grain domestication traits

Analysis of multi-year breeding program data revealed that the genetic architecture of an intermediate wheatgrass population was highly polygenic for both domestication and agronomic traits, supporting the use of genomic selection for new crop domestication. Perennial grains have the potential to provide food for humans and decrease the negative impacts of annual agriculture. Intermediate wheatgrass (IWG, Thinopyrum intermedium, Kernza®) is a promising perennial grain candidate that The Land Institute has been breeding since 2003. We evaluated four consecutive breeding cycles of IWG from 2016 to 2020 with each cycle containing approximately 1100 unique genets. Using genotyping-by-sequencing markers, quantitative trait loci (QTL) were mapped for 34 different traits using genome-wide association analysis. Combining data across cycles and years, we found 93 marker-trait associations for 16 different traits, with each association explaining 0.8-5.2% of the observed phenotypic variance. Across the four cycles, only three QTL showed an F_ST differentiation > 0.15 with two corresponding to a decrease in floret shattering. Additionally, one marker associated with brittle rachis was 216 bp from an ortholog of the btr2 gene. Power analysis and quantitative genetic theory were used to estimate the effective number of QTL, which ranged from a minimum of 33 up to 558 QTL for individual traits. This study suggests that key agronomic and domestication traits are under polygenic control and that molecular methods like genomic selection are needed to accelerate domestication and improvement of this new crop.

Perennials as Future Grain Crops: Opportunities and Challenges

Perennial grain crops could make a valuable addition to sustainable agriculture, potentially even as an alternative to their annual counterparts. The ability of perennials to grow year after year significantly reduces the number of agricultural inputs required, in terms of both planting and weed control, while reduced tillage improves soil health and on-farm biodiversity. Presently, perennial grain crops are not grown at large scale, mainly due to their early stages of domestication and current low yields. Narrowing the yield gap between perennial and annual grain crops will depend on characterizing differences in their life cycles, resource allocation, and reproductive strategies and understanding the trade-offs between annualism, perennialism, and yield. The genetic and biochemical pathways controlling plant growth, physiology, and senescence should be analyzed in perennial crop plants. This information could then be used to facilitate tailored genetic improvement of selected perennial grain crops to improve agronomic traits and enhance yield, while maintaining the benefits associated with perennialism.

QTL for seed shattering and threshability in intermediate wheatgrass align closely with well-studied orthologs from wheat, barley, and rice

Perennial grain crops have the potential to improve agricultural sustainability but few existing species produce sufficient grain yield to be economically viable. The outcrossing, allohexaploid, and perennial forage species intermediate wheatgrass (IWG) [Thinopyrum intermedium (Host) Barkworth & D. R. Dewey] has shown promise in undergoing direct domestication as a perennial grain crop using phenotypic and genomic selection. However, decades of selection will be required to achieve yields on par with annual small-grain crops. Marker-aided selection could accelerate progress if important genomic regions associated with domestication were identified. Here we use the IWG nested association mapping (NAM) population, with 1,168 F₁ progeny across 10 families to dissect the genetic control of brittle rachis, floret shattering, and threshability. We used a genome-wide association study (GWAS) with 8,003 single nucleotide polymorphism (SNP) markers and linkage mapping-both within-family and combined across families-with a robust phenotypic dataset collected from four unique year-by-location combinations. A total of 29 quantitative trait loci (QTL) using GWAS and 20 using the combined linkage analysis were detected, and most large-effect QTL were in common across the two analysis methods. We reveal that the genetic control of these traits in IWG is complex, with significant QTL across multiple chromosomes, sometimes within and across homoeologous groups and effects that vary depending on the family. In some cases, these QTL align within 216 bp to 31 Mbp of BLAST hits for known domestication genes in related species and may serve as precise targets of selection and directions for further study to advance the domestication of IWG.

Nested association mapping reveals the genetic architecture of spike emergence and anthesis timing in intermediate wheatgrass

Intermediate wheatgrass (Thinopyrum intermedium) is an outcrossing, cool season grass species currently undergoing direct domestication as a perennial grain crop. Though many traits are selection targets, understanding the genetic architecture of those important for local adaptation may accelerate the domestication process. Nested association mapping (NAM) has proven useful in dissecting the genetic control of agronomic traits many crop species, but its utility in primarily outcrossing, perennial species has yet to be demonstrated. Here, we introduce an intermediate wheatgrass NAM population developed by crossing ten phenotypically divergent donor parents to an adapted common parent in a reciprocal manner, yielding 1,168 F₁ progeny from 10 families. Using genotyping by sequencing, we identified 8,003 SNP markers and developed a population-specific consensus genetic map with 3,144 markers across 21 linkage groups. Using both genomewide association mapping and linkage mapping combined across and within families, we characterized the genetic control of flowering time. In the analysis of two measures of maturity across four separate environments, we detected as many as 75 significant QTL, many of which correspond to the same regions in both analysis methods across 11 chromosomes. The results demonstrate a complex genetic control that is variable across years, locations, traits, and within families. The methods were effective at detecting previously identified QTL, as well as new QTL that align closely to the well-characterized flowering time orthologs from barley, including Ppd-H1 and Constans. Our results demonstrate the utility of the NAM population for understanding the genetic control of flowering time and its potential for application to other traits of interest.

Sequenced-based paternity analysis to improve breeding and identify self-incompatibility loci in intermediate wheatgrass (Thinopyrum intermedium)

KEY MESSAGE

Paternity assignment and genome-wide association analyses for fertility were applied to a Thinopyrum intermedium breeding program. A lack of progeny between combinations of parents was associated with loci near self-incompatibility genes. In outcrossing species such as intermediate wheatgrass (IWG, Thinopyrum intermedium), polycrossing is often used to generate novel recombinants through each cycle of selection, but it cannot track pollen-parent pedigrees and it is unknown how self-incompatibility (SI) genes may limit the number of unique crosses obtained. This study investigated the potential of using next-generation sequencing to assign paternity and identify putative SI loci in IWG. Using a reference population of 380 individuals made from controlled crosses of 64 parents, paternity was assigned with 92% agreement using Cervus software. Using this approach, 80% of 4158 progeny (n = 3342) from a polycross of 89 parents were assigned paternity. Of the 89 pollen parents, 82 (92%) were represented with 1633 unique full-sib families representing 42% of all potential crosses. The number of progeny per successful pollen parent ranged from 1 to 123, with number of inflorescences per pollen parent significantly correlated to the number of progeny (r = 0.54, p < 0.001). Shannon's diversity index, assessing the total number and representation of families, was 7.33 compared to a theoretical maximum of 8.98. To test our hypothesis on the impact of SI genes, a genome-wide association study of the number of progeny observed from the 89 parents identified genetic effects related to non-random mating, including marker loci located near putative SI genes. Paternity testing of polycross progeny can impact future breeding gains by being incorporated in breeding programs to optimize polycross methodology, maintain genetic diversity, and reveal genetic architecture of mating patterns.

Genome-wide association analysis of natural variation in seed tocochromanols of barley

Tocochromanols (tocols for short), commonly called Vitamin E, are lipid-soluble plant antioxidants vital for regulating lipid peroxidation in chloroplasts and seeds. Barley (Hordeum vulgare L.) seeds contain all eight different isoforms of tocols; however, the extent of natural variation in their composition and their underlying genetic basis is not known. Tocol levels in barley seeds were quantified in diverse H. vulgare panels comprising 297 wild lines from a diversity panel and 160 cultivated spring-type accessions from the mini-core panel representing the genetic diversity of the USDA barley germplasm collection. Significant differences were observed in the concentration of tocols between the two panels. To identify the genes associated with tocols, genome-wide association analysis was conducted with single nucleotide polymorphisms (SNPs) from Illumina arrays for the mini-core panel and genotyping-by-sequencing for the wild barley panel. Forty unique SNPs in the wild barley and 27 SNPs in the mini-core panel were significantly associated with various tocols. Marker-trait associations (MTAs) were identified on chromosomes 1, 6, and 7 for key genes in the tocol biosynthesis pathway, which have also been reported in other studies. Several novel MTAs were identified on chromosomes 2, 3, 4 and 5 and were found to be in proximity to genes involved in the generation of precursor metabolites required for tocol biosynthesis. This study provides a valuable resource for barley breeding programs targeting specific isoforms of seed tocols and for investigating the physiological roles of these metabolites in seed longevity, dormancy, and germination.

Genome-Wide Association Study of Yield Component Traits in Intermediate Wheatgrass and Implications in Genomic Selection and Breeding

Intermediate wheatgrass (Thinopyrum intermedium, IWG) is a perennial grain crop with high biomass and grain yield, long seeds, and resistance to pests and diseases. It also reduces soil erosion, nitrate and mineral leaching into underground water tables, and sequesters carbon in its roots. The domestication timeline of IWG as a grain crop spans only 3 decades, hence it lags annual grain crops in yield and seed characteristics. One approach to improve its agronomic traits is by using molecular markers to uncover marker-trait associations. In this study, we performed association mapping on IWG breeding germplasm from the third recurrent selection cycle at the University of Minnesota. The IWG population was phenotyped in St Paul, MN in 2017 and 2018, and in Crookston, MN in 2018 for grain yield, seed length, width and weight, spike length and weight, and number of spikelets per spike. Strong positive correlations were observed among most trait pairs, with correlations as high as 0.76. Genotyping using high throughput sequencing identified 8,899 high-quality genome-wide SNPs which were combined with phenotypic data in association mapping to discover regions associated with the yield component traits. We detected 154 genetic loci associated with these traits of which 19 were shared between at least two traits. Prediction of breeding values using significant loci as fixed effects in genomic selection model improved predictive abilities by up to 14%. Genetic mapping of agronomic traits followed by using genomic selection to predict breeding values can assist breeders in selecting superior genotypes to accelerate IWG domestication.

Genome mapping of quantitative trait loci (QTL) controlling domestication traits of intermediate wheatgrass (Thinopyrum intermedium)

Allohexaploid (2n = 6x = 42) intermediate wheatgrass (Thinopyrum intermedium), abbreviated IWG, is an outcrossing perennial grass belonging to the tertiary gene pool of wheat. Perenniality would be valuable option for grain production, but attempts to introgress this complex trait from wheat-Thinopyrum hybrids have not been commercially successful. Efforts to breed IWG itself as a dual-purpose forage and grain crop have demonstrated useful progress and applications, but grain yields are significantly less than wheat. Therefore, genetic and physical maps have been developed to accelerate domestication of IWG. Herein, these maps were used to identify quantitative trait loci (QTLs) and candidate genes associated with IWG grain production traits in a family of 266 full-sib progenies derived from two heterozygous parents, M26 and M35. Transgressive segregation was observed for 17 traits related to seed size, shattering, threshing, inflorescence capacity, fertility, stem size, and flowering time. A total of 111 QTLs were detected in 36 different regions using 3826 genotype-by-sequence markers in 21 linkage groups. The most prominent QTL had a LOD score of 15 with synergistic effects of 29% and 22% over the family means for seed retention and percentage of naked seeds, respectively. Many QTLs aligned with one or more IWG gene models corresponding to 42 possible domestication orthogenes including the wheat Q and RHT genes. A cluster of seed-size and fertility QTLs showed possible alignment to a putative Z self-incompatibility gene, which could have detrimental grain-yield effects when genetic variability is low. These findings elucidate pathways and possible hurdles in the domestication of IWG.

Has Selection for Grain Yield Altered Intermediate Wheatgrass?

Perennial grains are demonstrating a greater probability of occupying land currently dedicated to other agricultural production. Arable land that is currently in use for forage or annual crop production becomes utilized. Breeding materials for the introduction of perennial grains directly into the human food chain has required utilizing existing plant materials in the domestication of species or manufacturing diverse crosses to introduce perenniality into existing crops. In the domestication of intermediate wheatgrass (Thinopyrum intermedium (Host), Barkworth and Dewey), existing forage cultivars or plant accessions were used to develop populations selected for grain production. A comparison of Cycle 3 materials from The Land Institute (TLI), Salina, KS, USA to USDA-Germplasm Resources Information Network (GRIN) accessions took place under space-planted field conditions at Carman, MB, Canada from 2011 to 2014. One hundred plants (75 TLI and 25 GRIN identified in May 2012) were followed through three seed harvests cycles with phenological, morphological and agronomic traits measured throughout. Selection for seed productivity (TLI materials) reduced the importance of biomass plant−1 on seed yield plant−1, leading to an increase in harvest index. Principal component analysis demonstrated the separation of the germplasm sources and the differential impact of years on the performance of all accessions. Path coefficient analysis also indicated that plant biomass production was of less importance on seed yield plant−1 in the TLI materials. Analysis removing area plant−1 as a factor increased both the importance of biomass and heads on seed yield cm−2 in the TLI materials, especially in the first two seed production years. Plant differences due to selection appear to have reduced overall plant area and increased harvest index in the TLI materials, indicating progress for grain yield under selection. However, a greater understanding of the dynamics within a seed production field is needed to provide insight into the development of more effective selection criteria for long-term field level production.