Sequencing 4.3 million mutations in wheat promoters to understand and modify gene expression.

Wheat is an important contributor to global food security, and further improvements are required to feed a growing human population. Functional genetics and genomics tools can help us to understand the function of different genes and to engineer beneficial changes. In this study, we used a promoter capture assay to sequence 2-kb regions upstream of all high-confidence annotated genes from 1,513 mutagenized plants from the tetraploid wheat variety Kronos. We identified 4.3 million induced mutations with an accuracy of 99.8%, resulting in a mutation density of 41.9 mutations per kb. We also remapped Kronos exome capture reads to Chinese Spring RefSeq v1.1, identified 4.7 million mutations, and predicted their effects on annotated genes. Using these predictions, we identified 59% more nonsynonymous substitutions and 49% more truncation mutations than in the original study. To show the biological value of the promoter dataset, we selected two mutations within the promoter of the VRN-A1 vernalization gene. Both mutations, located within transcription factor binding sites, significantly altered VRN-A1 expression, and one reduced the number of spikelets per spike. These publicly available sequenced mutant datasets provide rapid and inexpensive access to induced variation in the promoters and coding regions of most wheat genes. These mutations can be used to understand and modulate gene expression and phenotypes for both basic and commercial applications, where limited governmental regulations can facilitate deployment. These mutant collections, together with gene editing, provide valuable tools to accelerate functional genetic studies in this economically important crop.

Distribution of Cereal Cyst Nematodes (Heterodera avenae and H. filipjevi) in Eastern Washington State.

Cereal cyst nematodes (CCN; Heterodera avenae and H. filipjevi), cause substantial worldwide yield loss in small grain cereals such as wheat, barley, and oat. H. avenae was first detected in the United States in western Oregon in 1974 and had spread to northeast Oregon by the mid-1980s. Although H. avenae was detected in eastern Washington in 1984, extensive infestations were not recognized until 2010. H. filipjevi, first detected in Oregon in 2008, was found in eastern Washington in 2014. To gain more information about the distribution of these two species, an extensive survey was undertaken in eastern Washington, and methods were developed to distinguish species using DNA sequencing of single cysts. In this study, we surveyed 356 wheat and barley fields in eastern Washington from 2007 to 2017. CCN from the infested locations were identified to species level by sequencing the ribosomal internal transcribed spacers (ITS) and/or 28S ribosomal RNA (rRNA) genes. The sequences were compared in the GenBank database in the National Center for Biotechnology Information (NCBI) to identify species. The results show that H. filipjevi is primarily confined to southern Whitman County, WA; and H. avenae has a wider distribution across the higher precipitation annual cropping area of eastern Whitman County. Knowledge of species identification is critical for deployment of host resistance as an effective means of management, since resistance genes for one species of CCN may not be effective against the other.

Copy number and haplotype variation at the VRN-A1 and central FR-A2 loci are associated with frost tolerance in hexaploid wheat.

KEY MESSAGE: The interaction between VRN - A1 and FR - A2 largely affect the frost tolerance of hexaploid wheat. Frost tolerance is critical for wheat survival during cold winters. Natural variation for this trait is mainly associated with allelic differences at the VERNALIZATION 1 (VRN1) and FROST RESISTANCE 2 (FR2) loci. VRN1 regulates the transition between vegetative and reproductive stages and FR2, a locus including several tandemly duplicated C-REPEAT BINDING FACTOR (CBF) transcription factors, regulates the expression of Cold-regulated genes. We identified sequence and copy number variation at these two loci among winter and spring wheat varieties and characterized their association with frost tolerance. We identified two FR-A2 haplotypes-'FR-A2-S' and 'FR-A2-T'-distinguished by two insertion/deletions and ten single nucleotide polymorphisms within the CBF-A12 and CBF-A15 genes. Increased copy number of CBF-A14 was frequently associated with the FR-A2-T haplotype and with higher CBF14 transcript levels in response to cold. Factorial ANOVAs revealed significant interactions between VRN1 and FR-A2 for frost tolerance in both winter and spring panels suggesting a crosstalk between vernalization and cold acclimation pathways. The model including these two loci and their interaction explained 32.0 and 20.7 % of the variation in frost tolerance in the winter and spring panels, respectively. The interaction was validated in a winter wheat F 4:5 population segregating for both genes. Increased VRN-A1 copy number was associated with improved frost tolerance among varieties carrying the FR-A2-T allele but not among those carrying the FR-A2-S allele. These results suggest that selection of varieties carrying the FR-A2-T allele and three copies of the recessive vrn-A1 allele would be a good strategy to improve frost tolerance in wheat.

Effect of climate on the distribution of Fusarium spp. causing crown rot of wheat in the Pacific Northwest of the United States.

Fusarium crown rot (FCR) is one of the most widespread root and crown diseases of wheat in the Pacific Northwest (PNW) of the United States. Our objectives were to characterize crown rot severity and distribution throughout the PNW by conducting a survey of 210 fields covering the diverse dryland wheat-producing areas of Washington and Oregon and to utilize a factor analysis statistical approach to determine the effects of climate and geography on species distribution and disease severity. Climatic variables were based on 30-year averages and 2008 and 2009 separately (the 2 years of the survey). Mean annual temperature, mean temperature in the coldest month, mean temperature in the warmest month, mean annual precipitation, snowfall, elevation, soil type, and cropping intensity were highly intercorrelated. The factor analysis of the climate variables resulted in the development of two latent factors that could be used as predictor variables in logistic regression models for the presence or absence of Fusarium spp. and of FCR disease scores. Isolates of Fusarium spp. were recovered from 99% of 105 fields sampled in 2008 and 97% of fields in 2009. There were differences between years for responses of FCR and nodes scores, and isolations of Fusarium pseudograminearum with more significant results in 2008, due to warmer drier weather. Results of the factor analysis showed that the distribution of F. pseudograminearum occurred in a greater frequency in areas of the PNW at lower elevations with lower moisture and higher temperatures in 2008, whereas F. culmorum occurred in greater frequency in areas at higher elevations with moderate to high moisture and cooler temperatures consistently across both years. Disease scores increased with increasing levels of factors 1 (primarily temperature) and 2 (primarily precipitation). Both the frequency of pathogen species and disease scores were influenced by the year, indicating that soilborne pathogens are responsive to short-term changes in environment. This factor analysis approach can be utilized in studies to determine the effects of climate and other environmental (soil, cropping system, and so on) factors on the distribution and severity of root diseases.

Relationship Between Climatic Factors and Distribution of Pratylenchus spp. in the Dryland Wheat-Production Areas of Eastern Washington.

Field surveys were conducted by collecting soil samples to estimate nematode densities in soil from winter wheat, spring wheat, spring barley, and spring legumes (lentil, chickpea, and pea) fields during 2010 and 2011. Pratylenchus spp. were observed in 60% of sampled fields. However, nematodes were detected in nearly all of the survey fields in high numbers where crops were grown every year. To identify climatic variables associated with density of Pratylenchus spp. in soil, correlation and regression analyses were performed using climate data of survey sites from 1979 to 2010. Fifty-seven climate variables were significantly correlated with densities of Pratylenchus spp. All precipitation variables were significantly positively correlated with nematode abundance. Summer maximum air temperature was negatively correlated and winter minimum air temperature was positively correlated with nematode densities. In addition, both years' nematode densities were significantly correlated with cropping intensity. Five multivariate regression models for 2010 and seven models for 2011 nematode abundance levels were developed. The majority of the climate variables selected in the models were related to precipitation. Knowledge of root-lesion nematode distribution in the dryland region of eastern Washington and associated climate variables may be helpful to determine risk and apply management practices to minimize crop damage.

Spring Wheat Tolerance and Resistance to Heterodera avenae in the Pacific Northwest.

The cereal cyst nematode Heterodera avenae reduces wheat yields in the Pacific Northwest. Previous evaluations of cultivar resistance had been in controlled environments. Cultivar tolerance had not been evaluated. Seven spring wheat trials were conducted in naturally infested fields in three states over 2 years. A split-plot design was used for all trials. Five trials evaluated both tolerance and resistance in 1.8-by-9-m plots treated or not treated with nematicides. Two trials evaluated resistance in 1-m head rows where each wheat entry was paired with an adjacent row of a susceptible cultivar. Cultivars with the Cre1 resistance gene ('Ouyen' and 'Chara') reduced the postharvest density of H. avenae under field conditions, confirming Cre1 parents as useful for germplasm development. Ouyen was resistant but it was also intolerant, producing significantly lower grain yield in controls than in plots treated with nematicides. Susceptible cultivars varied in tolerance. Undefined resistance was identified in one commercial cultivar ('WB-Rockland) and four breeding lines (UC1711, SO900163, SY-B041418, and SY-97621-05). This research was the first systematic field demonstration of potential benefits to be derived through development and deployment of cultivars with resistance plus tolerance to cereal cyst nematode in North America.

Identification of genetic factors controlling kernel hardness and related traits in a recombinant inbred population derived from a soft × 'extra-soft' wheat (Triticum aestivum L.) cross.

Kernel hardness or texture, used to classify wheat (Triticum aestivum L.) into soft and hard classes, is a major determinant of milling and baking quality. Wheat genotypes in the soft class that are termed 'extra-soft' (with kernel hardness in the lower end of the spectrum) have been associated with superior end-use quality. In order to better understand the relationship between kernel hardness, milling yield, and various agronomic traits, we performed quantitative trait mapping using a recombinant inbred line population derived from a cross between a common soft wheat line and a genotype classified as an 'extra-soft' line. A total of 47 significant quantitative trait loci (QTL) (LOD ≥ 3.0) were identified for nine traits with the number of QTL affecting each trait ranging from three to nine. The percentage of phenotypic variance explained by these QTL ranged from 3.7 to 50.3%. Six QTL associated with kernel hardness and break flour yield were detected on chromosomes 1BS, 4BS, 5BS, 2DS, 4DS, and 5DL. The two most important QTL were mapped onto orthologous regions on chromosomes 4DS (Xbarc1118-Rht-D1) and 4BS (Xwmc617-Rht-B1). These results indicated that the 'extra-soft' characteristic was not controlled by the Hardness (Ha) locus on chromosome 5DS. QTL for eight agronomic traits occupied two genomic regions near semi-dwarf genes Rht-D1 on chromosome 4DS and Rht-B1 on chromosome 4BS. The clustering of these QTL is either due to the pleiotropic effects of single genes or tight linkage of genes controlling these various traits.

Characterizing and Validating Stripe Rust Resistance Loci in US Pacific Northwest Winter Wheat Accessions (Triticum aestivum L.) by Genome-wide Association and Linkage Mapping.

Stripe rust resistance is a critical need for wheat cultivars in the US Pacific Northwest (PNW). Our previous genome-wide association study (GWAS) for stripe rust resistance in a set of PNW winter wheat accessions (Panel-2) identified multiple marker-trait associations (MTAs) for both all-stage and field resistance. In this study, we conducted additional GWAS using a different set of PNW winter wheat accessions (Panel-1) that contained recently bred soft white winter wheat breeding lines and cultivars. A total of 12 all-stage resistance MTAs and eight field resistance MTAs were identified. Within these MTAs, nine MTAs for all-stage resistance and two MTAs for field resistance were located distinctly from previously characterized genes and likely represent novel loci. Markers IWB60567 (1B), IWB24342 (2A), and IWB46564 (2B) explained the largest phenotypic variances for disease responses. The analysis confirmed that MTAs on chromosome 1B were indeed the same as identified in Panel-2 and that MTAs on chromosome 2A were likely and closely linked to another field resistance QTL, (Panel-2). Haplotypes for MTAs on chromosome 1B, , and linked loci on chromosome 2A provide useful information for marker development and introgression of these QTL into wheat breeding programs.

Genome-wide association mapping for eyespot disease in US Pacific Northwest winter wheat.

Eyespot, caused by the soil-borne necrotrophic fungi Oculimacula yallundae and O. acuformis, is a disease of major economic significance for wheat, barley and rye. Pacific Northwest (PNW) winter wheat (Triticum aestivum L.) grown in areas of high rainfall and moderate winters is most vulnerable to infection. The objective of this research was to identify novel genomic regions associated with eyespot resistance in winter wheat adapted to the PNW. Two winter wheat panels of 469 and 399 lines were compiled for one of the first genome-wide association studies (GWAS) of eyespot resistance in US winter wheat germplasm. These panels were genotyped with the Infinium 9K and 90K iSelect SNP arrays. Both panels were phenotyped for disease resistance in a two-year field study and in replicated growth chamber trials. Growth chamber trials were used to evaluate the genetic resistance of O. acuformis and O. yallundae species separately. Best linear unbiased predictors (BLUPs) were calculated across all field and growth chamber environments. A total of 73 marker-trait associations (MTAs) were detected on nine different chromosomes (1A, 2A, 2B, 4A, 5A, 5B, 7A, 7B and 7D) that were significantly associated (p-value <0.001) with eyespot resistance in Panel A, and 19 MTAs on nine different chromosomes (1A, 1B, 2A, 2D, 3B, 5A, 5B, 7A, and 7B) in Panel B. The most significant SNPs were associated with Pch1 and Pch2 resistance genes on the long arms of chromosome 7D and 7A. Most of the novel MTAs appeared to have a minor effect on reducing eyespot disease. Nevertheless, eyespot disease scores decreased as the number of resistance alleles increased. Seven SNP markers, significantly associated with reducing eyespot disease across environments and in the absence and presence of Pch1 were identified. These markers were located on chromosomes 2A (IWB8331), 5A (IWB73709), 5B (IWB47298), 7AS (IWB47160), 7B (IWB45005) and two SNPs (Ex_c44379_2509 and IAAV4340) had unknown map positions. The additive effect of the MTAs explained most of the remaining phenotypic variation not accounted for by Pch1 or Pch2. This study provides breeders with adapted germplasm and novel sources of eyespot resistance to be used in the development of superior cultivars with increased eyespot resistance.

Genetic Dissection of End-Use Quality Traits in Adapted Soft White Winter Wheat.

Soft white wheat is used in domestic and foreign markets for various end products requiring specific quality profiles. Phenotyping for end-use quality traits can be costly, time-consuming and destructive in nature, so it is advantageous to use molecular markers to select experimental lines with superior traits. An association mapping panel of 469 soft white winter wheat cultivars and advanced generation breeding lines was developed from regional breeding programs in the U.S. Pacific Northwest. This panel was genotyped on a wheat-specific 90 K iSelect single nucleotide polymorphism (SNP) chip. A total of 15,229 high quality SNPs were selected and combined with best linear unbiased predictions (BLUPs) from historical phenotypic data of the genotypes in the panel. Genome-wide association mapping was conducted using the Fixed and random model Circulating Probability Unification (FarmCPU). A total of 105 significant marker-trait associations were detected across 19 chromosomes. Potentially new loci for total flour yield, lactic acid solvent retention capacity, flour sodium dodecyl sulfate sedimentation and flour swelling volume were also detected. Better understanding of the genetic factors impacting end-use quality enable breeders to more effectively discard poor quality germplasm and increase frequencies of favorable end-use quality alleles in their breeding populations.

Genome-Wide Association Mapping for Tolerance to Preharvest Sprouting and Low Falling Numbers in Wheat.

Preharvest sprouting (PHS), the germination of grain on the mother plant under cool and wet conditions, is a recurring problem for wheat farmers worldwide. α-amylase enzyme produced during PHS degrades starch resulting in baked good with poor end-use quality. The Hagberg-Perten Falling Number (FN) test is used to measure this problem in the wheat industry, and determines how much a farmer's wheat is discounted for PHS damage. PHS tolerance is associated with higher grain dormancy. Thus, breeding programs use germination-based assays such as the spike-wetting test to measure PHS susceptibility. Association mapping identified loci associated with PHS tolerance in U.S. Pacific Northwest germplasm based both on FN and on spike-wetting test data. The study was performed using a panel of 469 white winter wheat cultivars and elite breeding lines grown in six Washington state environments, and genotyped for 15,229 polymorphic markers using the 90k SNP Illumina iSelect array. Marker-trait associations were identified using the FarmCPU R package. Principal component analysis was directly and a kinship matrix was indirectly used to account for population structure. Nine loci were associated with FN and 34 loci associated with PHS based on sprouting scores. None of the QFN.wsu loci were detected in multiple environments, whereas six of the 34 QPHS.wsu loci were detected in two of the five environments. There was no overlap between the QTN detected based on FN and PHS, and there was little correlation between the two traits. However, both traits appear to be PHS-related since 19 of the 34 QPHS.wsu loci and four of the nine QFN.wsu loci co-localized with previously published dormancy and PHS QTL. Identification of these loci will lead to a better understanding of the genetic architecture of PHS and will help with the future development of genomic selection models.

Genomic Regions Associated with Tolerance to Freezing Stress and Snow Mold in Winter Wheat.

Plants grown through the winter are subject to selective pressures that vary with each year's unique conditions, necessitating tolerance of numerous abiotic and biotic stress factors. The objective of this study was to identify molecular markers in winter wheat (Triticum aestivum L.) associated with tolerance of two of these stresses, freezing temperatures and snow mold-a fungal disease complex active under snow cover. A population of 155 F2:5 recombinant inbred lines from a cross between soft white wheat cultivars "Finch" and "Eltan" was evaluated for snow mold tolerance in the field, and for freezing tolerance under controlled conditions. A total of 663 molecular markers was used to construct a genetic linkage map and identify marker-trait associations. One quantitative trait locus (QTL) associated with both freezing and snow mold tolerance was identified on chromosome 5A. A second, distinct, QTL associated with freezing tolerance also was found on 5A, and a third on 4B. A second QTL associated with snow mold tolerance was identified on chromosome 6B. The QTL on 5A associated with both traits was closely linked with the Fr-A2 (Frost-Resistance A2) locus; its significant association with both traits may have resulted from pleiotropic effects, or from greater low temperature tolerance enabling the plants to better defend against snow mold pathogens. The QTL on 4B associated with freezing tolerance, and the QTL on 6B associated with snow mold tolerance have not been reported previously, and may be useful in the identification of sources of tolerance for these traits.

Measuring freezing tolerance: survival and regrowth assays.

Screening plants for freezing tolerance under tightly controlled conditions is an invaluable technique for studying freezing tolerance and selecting for improved winterhardiness. Artificial freezing tests of cereal plants historically have used isolated crown and stem tissue prepared by "removing all plant parts 3 cm above and 0.5 cm below the crown tissue" (Fowler et al., Crop Sci 21:896-901, 1981). Here, we describe a method of conducting freezing tolerance tests using intact plants grown in small horticultural containers, including suggested methods for collecting and analyzing the data.

Mapping stripe rust resistance in a BrundageXCoda winter wheat recombinant inbred line population.

A recombinant inbred line (RIL) mapping population developed from a cross between winter wheat (Triticum aestivum L.) cultivars Coda and Brundage was evaluated for reaction to stripe rust (caused by Puccinia striiformis f. sp. tritici). Two hundred and sixty eight RIL from the population were evaluated in replicated field trials in a total of nine site-year locations in the U.S. Pacific Northwest. Seedling reaction to stripe rust races PST-100, PST-114 and PST-127 was also examined. A linkage map consisting of 2,391 polymorphic DNA markers was developed covering all chromosomes of wheat with the exception of 1D. Two QTL on chromosome 1B were associated with adult plant and seedling reaction and were the most significant QTL detected. Together these QTL reduced adult plant infection type from a score of seven to a score of two reduced disease severity by an average of 25% and provided protection against race PST-100, PST-114 and PST-127 in the seedling stage. The location of these QTL and the race specificity provided by them suggest that observed effects at this locus are due to a complementation of the previously known but defeated resistances of the cultivar Tres combining with that of Madsen (the two parent cultivars of Coda). Two additional QTL on chromosome 3B and one on 5B were associated with adult plant reaction only, and a single QTL on chromosome 5D was associated with seedling reaction to PST-114. Coda has been resistant to stripe rust since its release in 2000, indicating that combining multiple resistance genes for stripe rust provides durable resistance, especially when all-stage resistance genes are combined in a fashion to maximize the number of races they protect against. Identified molecular markers will allow for an efficient transfer of these genes into other cultivars, thereby continuing to provide excellent resistance to stripe rust.

Registration of Zak ERA8 Soft White Spring Wheat Germplasm with Enhanced Response to ABA and Increased Seed Dormancy.

Zak ERA8 (ENHANCED RESPONSE to ABA8) (Reg. No. GP-966, PI 669443) is a unique line derived from soft white spring wheat (Triticum aestivum L.) cultivar Zak that has increased seed dormancy but after-ripens within 10 to 16 wk. The goal in developing this germplasm was to use increased seed dormancy to improve tolerance to preharvest sprouting, a problem that can cause severe economic losses. This germplasm was developed by USDA-ARS, Pullman, WA, in collaboration with Washington State University. Zak ERA8was tested under experimental number 60.1.27.10. The ERA8mutation was generated by chemical mutagenesis followed by selection for the inability to germinate on abscisic acid (ABA) concentrations too low to inhibit wild-type Zak seed germination. The semidominant Zak ERA8 line has been backcrossed twice to wild-type Zak. Following the first backcross, Zak ERA8 showed similar morphological and grain quality traits to the original Zak cultivar.