Whole-genome sequencing of diverse wheat accessions uncovers genetic changes during modern breeding in China and the United States.

Breeding has dramatically changed the plant architecture of wheat (Triticum aestivum), resulting in the development of high-yielding varieties adapted to modern farming systems. However, how wheat breeding shaped the genomic architecture of this crop remains poorly understood. Here, we performed a comprehensive comparative analysis of a whole-genome resequencing panel of 355 common wheat accessions (representing diverse landraces and modern cultivars from China and the United States) at the phenotypic and genomic levels. The genetic diversity of modern wheat cultivars was clearly reduced compared to landraces. Consistent with these genetic changes, most phenotypes of cultivars from China and the United States were significantly altered. Of the 21 agronomic traits investigated, 8 showed convergent changes between the 2 countries. Moreover, of the 207 loci associated with these 21 traits, more than half overlapped with genomic regions that showed evidence of selection. The distribution of selected loci between the Chinese and American cultivars suggests that breeding for increased productivity in these 2 regions was accomplished by pyramiding both shared and region-specific variants. This work provides a framework to understand the genetic architecture of the adaptation of wheat to diverse agricultural production environments, as well as guidelines for optimizing breeding strategies to design better wheat varieties.

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.

Introgressions lead to reference bias in wheat RNA-seq analysis.

BACKGROUND: RNA-seq is a fundamental technique in genomics, yet reference bias, where transcripts derived from non-reference alleles are quantified less accurately, can undermine the accuracy of RNA-seq quantification and thus the conclusions made downstream. Reference bias in RNA-seq analysis has yet to be explored in complex polyploid genomes despite evidence that they are often a complex mosaic of wild relative introgressions, which introduce blocks of highly divergent genes. RESULTS: Here we use hexaploid wheat as a model complex polyploid, using both simulated and experimental data to show that RNA-seq alignment in wheat suffers from widespread reference bias which is largely driven by divergent introgressed genes. This leads to underestimation of gene expression and incorrect assessment of homoeologue expression balance. By incorporating gene models from ten wheat genome assemblies into a pantranscriptome reference, we present a novel method to reduce reference bias, which can be readily scaled to capture more variation as new genome and transcriptome data becomes available. CONCLUSIONS: This study shows that the presence of introgressions can lead to reference bias in wheat RNA-seq analysis. Caution should be exercised by researchers using non-sample reference genomes for RNA-seq alignment and novel methods, such as the one presented here, should be considered.

CRISPR-based editing of the ω- and γ-gliadin gene clusters reduces wheat immunoreactivity without affecting grain protein quality.

Wheat immunotoxicity is associated with abnormal reaction to gluten-derived peptides. Attempts to reduce immunotoxicity using breeding and biotechnology often affect dough quality. Here, the multiplexed CRISPR-Cas9 editing of cultivar Fielder was used to modify gluten-encoding genes, specifically focusing on ω- and γ-gliadin gene copies, which were identified to be abundant in immunoreactive peptides based on the analysis of wheat genomes assembled using the long-read sequencing technologies. The whole-genome sequencing of an edited line showed mutation or deletion of nearly all ω-gliadin and half of the γ-gliadin gene copies and confirmed the lack of editing in the α/ß-gliadin genes. The estimated 75% and 64% reduction in ω- and γ-gliadin content, respectively, had no negative impact on the end-use quality characteristics of grain protein and dough. A 47-fold immunoreactivity reduction compared to a non-edited line was demonstrated using antibodies against immunotoxic peptides. Our results indicate that the targeted CRISPR-based modification of the ω- and γ-gliadin gene copies determined to be abundant in immunoreactive peptides by analysing high-quality genome assemblies is an effective mean for reducing immunotoxicity of wheat cultivars while minimizing the impact of editing on protein quality.

Quantitative trait loci for rolled leaf in a wheat EMS mutant from Jagger.

KEY MESSAGE: Two QTLs with major effects on rolled leaf trait were consistently detected on chromosomes 1A (QRl.hwwg-1AS) and 5A (QRl.hwwg-5AL) in the field experiments. Rolled leaf (RL) is a morphological strategy to protect plants from dehydration under stressed field conditions. Identification of quantitative trait loci (QTLs) underlining RL is essential to breed drought-tolerant wheat cultivars. A mapping population of 154 recombinant inbred lines was developed from the cross between JagMut1095, a mutant of Jagger, and Jagger to identify quantitative trait loci (QTLs) for the RL trait. A linkage map of 3106 cM was constructed with 1003 unique SNPs from 21 wheat chromosomes. Two consistent QTLs were identified for RL on chromosomes 1A (QRl.hwwg-1AS) and 5A (QRl.hwwg-5AL) in all field experiments. QRl.hwwg-1AS explained 24-56% of the phenotypic variation and QRl.hwwg-5AL explained up to 20% of the phenotypic variation. The combined percent phenotypic variation associated with the two QTLs was up to 61%. Analyses of phenotypic and genotypic data of recombinants generated from heterogeneous inbred families of JagMut1095 × Jagger delimited QRl.hwwg-1AS to a 6.04 Mb physical interval. This work lays solid foundation for further fine mapping and map-based cloning of QRl.hwwg-1AS.

Multiplexed promoter and gene editing in wheat using a virus-based guide RNA delivery system.

The low efficiency of genetic transformation and gene editing across diverse cultivars hinder the broad application of CRISPR technology for crop improvement. The development of virus-based methods of CRISPR-Cas system delivery into the plant cells holds great promise to overcome these limitations. Here, we perform direct inoculation of wheat leaves with the barley stripe mosaic virus (BSMV) transcripts to deliver guide RNAs (sgRNA) into the Cas9-expressing wheat. We demonstrate that wheat inoculation with the pool of BSMV-sgRNAs could be used to generate heritable precise deletions in the promoter region of a transcription factor and to perform multiplexed editing of agronomic genes. We transfer the high-expressing locus of Cas9 into adapted spring and winter cultivars by marker-assisted introgression and use of the BSMV-sgRNAs to edit two agronomic genes. A strategy presented in our study could be applied to any adapted cultivar for creating new cis-regulatory diversity or large-scale editing of multiple genes in biological pathways or QTL regions, opening possibilities for the effective engineering of crop genomes, and accelerating gene discovery and trait improvement efforts.

Identification and characterization of a natural polymorphism in FT-A2 associated with increased number of grains per spike in wheat.

KEY MESSAGE: We discovered a natural FT-A2 allele that increases grain number per spike in both pasta and bread wheat with limited effect on heading time. Increases in wheat grain yield are necessary to meet future global food demands. A previous study showed that loss-of-function mutations in FLOWERING LOCUS T2 (FT2) increase spikelet number per spike (SNS), an important grain yield component. However, these mutations were also associated with reduced fertility, offsetting the beneficial effect of the increases in SNS on grain number. Here, we report a natural mutation resulting in an aspartic acid to alanine change at position 10 (D10A) associated with significant increases in SNS and no negative effects on fertility. Using a high-density genetic map, we delimited the SNS candidate region to a 5.2-Mb region on chromosome 3AS including 28 genes. Among them, only FT-A2 showed a non-synonymous polymorphism (D10A) present in two different populations segregating for the SNS QTL on chromosome arm 3AS. These results, together with the known effect of the ft-A2 mutations on SNS, suggest that variation in FT-A2 is the most likely cause of the observed differences in SNS. We validated the positive effects of the A10 allele on SNS, grain number, and grain yield per spike in near-isogenic tetraploid wheat lines and in an hexaploid winter wheat population. The A10 allele is present at very low frequency in durum wheat and at much higher frequency in hexaploid wheat, particularly in winter and fall-planted spring varieties. These results suggest that the FT-A2 A10 allele may be particularly useful for improving grain yield in durum wheat and fall-planted common wheat varieties.

Expanding the range of editable targets in the wheat genome using the variants of the Cas12a and Cas9 nucleases.

The development of CRISPR-based editors recognizing distinct protospacer-adjacent motifs (PAMs), or having different spacer length/structure requirements broadens the range of possible genomic applications. We evaluated the natural and engineered variants of Cas12a (FnCas12a and LbCas12a) and Cas9 for their ability to induce mutations in endogenous genes controlling important agronomic traits in wheat. Unlike FnCas12a, LbCas12a-induced mutations in the wheat genome, even though with a lower rate than that reported for SpCas9. The eight-fold improvement in the gene editing efficiency was achieved for LbCas12a by using the guides flanked by ribozymes and driven by the RNA polymerase II promoter from switchgrass. The efficiency of multiplexed genome editing (MGE) using LbCas12a was mostly similar to that obtained using the simplex RNA guides and showed substantial increase after subjecting transgenic plants to high-temperature treatment. We successfully applied LbCas12a-MGE for generating heritable mutations in a gene controlling grain size and weight in wheat. We showed that the range of editable loci in the wheat genome could be further expanded by using the engineered variants of Cas12a (LbCas12a-RVR) and Cas9 (Cas9-NG and xCas9) that recognize the TATV and NG PAMs, respectively, with the Cas9-NG showing higher editing efficiency on the targets with atypical PAMs compared to xCas9. In conclusion, our study reports a set of validated natural and engineered variants of Cas12a and Cas9 editors for targeting loci in the wheat genome not amenable to modification using the original SpCas9 nuclease.

QTL mapping for grain yield and three yield components in a population derived from two high-yielding spring wheat cultivars.

KEY MESSAGE: Four genomic regions on chromosomes 4A, 6A, 7B, and 7D were discovered, each with multiple tightly linked QTL (QTL clusters) associated with two to three yield components. The 7D QTL cluster was associated with grain yield, fertile spikelet number per spike, thousand kernel weight, and heading date. It was located in the flanking region of FT-D1, a homolog gene of Arabidopsis FLOWERING LOCUS T, a major gene that regulates wheat flowering. Genetic manipulation of yield components is an important approach to increase grain yield in wheat (Triticum aestivum). The present study used a mapping population comprised of 181 doubled haploid lines derived from two high-yielding spring wheat cultivars, UI Platinum and LCS Star. The two cultivars and the derived population were assessed for six traits in eight field trials primarily in Idaho in the USA. The six traits were grain yield, fertile spikelet number per spike, productive tiller number per unit area, thousand kernel weight, heading date, and plant height. Quantitative Trait Locus (QTL) analysis of the six traits was conducted using 14,236 single-nucleotide polymorphism (SNP) markers generated from the wheat 90 K SNP and the exome and promoter capture arrays. Of the 19 QTL detected, 14 were clustered in four chromosomal regions on 4A, 6A, 7B and 7D. Each of the four QTL clusters was associated with multiple yield component traits, and these traits were often negatively correlated with one another. As a result, additional QTL dissection studies are needed to optimize trade-offs among yield component traits for specific production environments. Kompetitive allele-specific PCR markers for the four QTL clusters were developed and assessed in an elite spring wheat panel of 170 lines, and eight of the 14 QTL were validated. The two parents contain complementary alleles for the four QTL clusters, suggesting the possibility of improving grain yield via genetic recombination of yield component loci.

Gene editing of the wheat homologs of TONNEAU1-recruiting motif encoding gene affects grain shape and weight in wheat.

Grain size and weight are important components of a suite of yield-related traits in crops. Here, we showed that the CRISPR-Cas9 gene editing of TaGW7, a homolog of rice OsGW7 encoding a TONNEAU1-recruiting motif (TRM) protein, affects grain shape and weight in allohexaploid wheat. By editing the TaGW7 homoeologs in the B and D genomes, we showed that mutations in either of the two or both genomes increased the grain width and weight but reduced the grain length. The effect sizes of mutations in the TaGW7 gene homoeologs coincided with the relative levels of their expression in the B and D genomes. The effects of gene editing on grain morphology and weight traits were dosage dependent with the double-copy mutant showing larger effect than the respective single copy mutants. The TaGW7-centered gene co-expression network indicated that this gene is involved in the pathways regulating cell division and organ growth, also confirmed by the cellular co-localization of TaGW7 with α- and ß-tubulin proteins, the building blocks of microtubule arrays. The analyses of exome capture data in tetraploid domesticated and wild emmer, and hexaploid wheat revealed the loss of diversity around TaGW7-associated with domestication selection, suggesting that TaGW7 is likely to play an important role in the evolution of yield component traits in wheat. Our study showed how integrating CRISPR-Cas9 system with cross-species comparison can help to uncover the function of a gene fixed in wheat for allelic variants targeted by domestication selection and select targets for engineering new gene variants for crop improvement.

Dissection of Cell Death Induction by Wheat Stem Rust Resistance Protein Sr35 and Its Matching Effector AvrSr35.

Nucleotide-binding leucine-rich repeat receptors (NLRs) are the most abundant type of immune receptors in plants and can trigger a rapid cell-death (hypersensitive) response upon sensing pathogens. We previously cloned the wheat NLR Sr35, which encodes a coiled-coil (CC) NLR that confers resistance to the virulent wheat stem rust race Ug99. Here, we investigated Sr35 signaling after Agrobacterium-mediated transient expression in Nicotiana benthamiana. Expression of Sr35 in N. benthamiana leaves triggered a mild cell-death response, which is enhanced in the autoactive mutant Sr35 D503V. The N-terminal tagging of Sr35 with green fluorescent protein (GFP) blocked the induction of cell death, whereas a C-terminal GFP tag did not. No domain truncations of Sr35 generated cell-death responses as strong as the wild type, but a truncation including the NB-ARC (nucleotide binding adaptor) shared by APAF-1, R proteins, and CED-4 domains in combination with the D503V autoactive mutation triggered cell death. In addition, coexpression of Sr35 with the matching pathogen effector protein AvrSr35 resulted in robust cell death and electrolyte leakage levels that were similar to autoactive Sr35 and significantly higher than Sr35 alone. Coexpression of Sr35-CC-NB-ARC and AvrSr35 did not induce cell death, confirming the importance of the leucine-rich repeat (LRR) domain for AvrSr35 recognition. These findings were confirmed through Agrobacterium-mediated transient expression in barley. Taken together, these results implicate the CC-NB-ARC domains of Sr35 in inducing cell death and the LRR domain in AvrSr35 recognition.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Sequence-based mapping identifies a candidate transcription repressor underlying awn suppression at the B1 locus in wheat.

Awns are stiff, hair-like structures which grow from the lemmas of wheat (Triticum aestivum) and other grasses that contribute to photosynthesis and play a role in seed dispersal. Variation in awn length in domesticated wheat is controlled primarily by three major genes, most commonly the dominant awn suppressor Tipped1 (B1). This study identifies a transcription repressor responsible for awn inhibition at the B1 locus. Association mapping was combined with analysis in biparental populations to delimit B1 to a distal region of 5AL colocalized with QTL for number of spikelets per spike, kernel weight, kernel length, and test weight. Fine-mapping located B1 to a region containing only two predicted genes, including C2H2 zinc finger transcriptional repressor TraesCS5A02G542800 upregulated in developing spikes of awnless individuals. Deletions encompassing this candidate gene were present in awned mutants of an awnless wheat. Sequence polymorphisms in the B1 coding region were not observed in diverse wheat germplasm whereas a nearby polymorphism was highly predictive of awn suppression. Transcriptional repression by B1 is the major determinant of awn suppression in global wheat germplasm. It is associated with increased number of spikelets per spike and decreased kernel size.

The genetic architecture of genome-wide recombination rate variation in allopolyploid wheat revealed by nested association mapping.

Recombination affects the fate of alleles in populations by imposing constraints on the reshuffling of genetic information. Understanding the genetic basis of these constraints is critical for manipulating the recombination process to improve the resolution of genetic mapping, and reducing the negative effects of linkage drag and deleterious genetic load in breeding. Using sequence-based genotyping of a wheat nested association mapping (NAM) population of 2,100 recombinant inbred lines created by crossing 29 diverse lines, we mapped QTL affecting the distribution and frequency of 102 000 crossovers (CO). Genome-wide recombination rate variation was mostly defined by rare alleles with small effects together explaining up to 48.6% of variation. Most QTL were additive and showed predominantly trans-acting effects. The QTL affecting the proximal COs also acted additively without increasing the frequency of distal COs. We showed that the regions with decreased recombination carry more single nucleotide polymorphisms (SNPs) with possible deleterious effects than the regions with a high recombination rate. Therefore, our study offers insights into the genetic basis of recombination rate variation in wheat and its effect on the distribution of deleterious SNPs across the genome. The identified trans-acting additive QTL can be utilized to manipulate CO frequency and distribution in the large polyploid wheat genome opening the possibility to improve the efficiency of gene pyramiding and reducing the deleterious genetic load in the low-recombining pericentromeric regions of chromosomes.

Analysis of recombinant inbred line populations derived from wheat landraces to identify new genes for wheat stem sawfly resistance.

Wheat landrace accessions were chosen from areas of the world with historical European wheat stem sawfly (Cephus pygmaeus L.) selection pressure to develop six recombinant inbred line (RIL) populations. Molecular maps were constructed, and resistance due to antibiosis and antixenosis was assessed at sites in Montana naturally infested by Cephus cinctus Norton, the wheat stem sawfly (WSS). Novel QTLs were identified along with QTL previously identified in elite germplasm. A newly identified QTL on chromosome 1B provided a new source for pith-filled solid stems. An allele for resistance on chromosome 4A unrelated to solid stems was identified in four of the six RIL populations. A landrace from Turkey, PI 166471, contained alleles at three QTLs causing high levels of larval mortality. None of the QTLs were related to stem solidness, but their combined effect provided resistance similar to that observed in a solid-stemmed check cultivar. These results show the utility of genetic populations derived from geographically targeted landrace accessions to identify new alleles for insect resistance. New PCR-based molecular markers were developed for introgression of novel alleles for WSS resistance into elite lines. Comparison of results with previous analysis of elite cultivars addresses changes in allele frequencies during the wheat breeding process.

Identification of a candidate gene for a QTL for spikelet number per spike on wheat chromosome arm 7AL by high-resolution genetic mapping.

KEY MESSAGE: A high-resolution genetic map combined with haplotype analyses identified a wheat ortholog of rice gene APO1 as the best candidate gene for a 7AL locus affecting spikelet number per spike. A better understanding of the genes controlling differences in wheat grain yield components can accelerate the improvements required to satisfy future food demands. In this study, we identified a promising candidate gene underlying a quantitative trait locus (QTL) on wheat chromosome arm 7AL regulating spikelet number per spike (SNS). We used large heterogeneous inbred families ( > 10,000 plants) from two crosses to map the 7AL QTL to an 87-kb region (674,019,191-674,106,327 bp, RefSeq v1.0) containing two complete and two partial genes. In this region, we found three major haplotypes that were designated as H1, H2 and H3. The H2 haplotype contributed the high-SNS allele in both H1 × H2 and H2 × H3 segregating populations. The ancestral H3 haplotype is frequent in wild emmer (48%) but rare (~ 1%) in cultivated wheats. By contrast, the H1 and H2 haplotypes became predominant in modern cultivated durum and common wheat, respectively. Among the four candidate genes, only TraesCS7A02G481600 showed a non-synonymous polymorphism that differentiated H2 from the other two haplotypes. This gene, designated here as WHEAT ORTHOLOG OF APO1 (WAPO1), is an ortholog of the rice gene ABERRANT PANICLE ORGANIZATION 1 (APO1), which affects spikelet number. Taken together, the high-resolution genetic map, the association between polymorphisms in the different mapping populations with differences in SNS, and the known role of orthologous genes in other grass species suggest that WAPO-A1 is the most likely candidate gene for the 7AL SNS QTL among the four genes identified in the candidate gene region.

Major structural genomic alterations can be associated with hybrid speciation in Aegilops markgrafii (Triticeae).

During evolutionary history many grasses from the tribe Triticeae have undergone interspecific hybridization, resulting in allopolyploidy; whereas homoploid hybrid speciation was found only in rye. Homoeologous chromosomes within the Triticeae preserved cross-species macrocolinearity, except for a few species with rearranged genomes. Aegilops markgrafii, a diploid wild relative of wheat (2n = 2x = 14), has a highly asymmetrical karyotype that is indicative of chromosome rearrangements. Molecular cytogenetics and next-generation sequencing were used to explore the genome organization. Fluorescence in situ hybridization with a set of wheat cDNAs allowed the macrostructure and cross-genome homoeology of the Ae. markgrafii chromosomes to be established. Two chromosomes maintained colinearity, whereas the remaining were highly rearranged as a result of inversions and inter- and intrachromosomal translocations. We used sets of barley and wheat orthologous gene sequences to compare discrete parts of the Ae. markgrafii genome involved in the rearrangements. Analysis of sequence identity profiles and phylogenic relationships grouped chromosome blocks into two distinct clusters. Chromosome painting revealed the distribution of transposable elements and differentiated chromosome blocks into two groups consistent with the sequence analyses. These data suggest that introgressive hybridization accompanied by gross chromosome rearrangements might have had an impact on karyotype evolution and homoploid speciation in Ae. markgrafii.

Identification and validation of QTL for grain yield and plant water status under contrasting water treatments in fall-sown spring wheats.

KEY MESSAGE: Chromosome regions affecting grain yield, grain yield components and plant water status were identified and validated in fall-sown spring wheats grown under full and limited irrigation. Increases in wheat production are required to feed a growing human population. To understand the genetic basis of grain yield in fall-sown spring wheats, we performed a genome-wide association study (GWAS) including 262 photoperiod-insensitive spring wheat accessions grown under full and limited irrigation treatments. Analysis of molecular variance showed that 4.1% of the total variation in the panel was partitioned among accessions originally developed under fall-sowing or spring-sowing conditions, 11.7% among breeding programs within sowing times and 84.2% among accessions within breeding programs. We first identified QTL for grain yield, yield components and plant water status that were significant in at least three environments in the GWAS, and then selected those that were also significant in at least two environments in a panel of eight biparental mapping populations. We identified and validated 14 QTL for grain yield, 15 for number of spikelets per spike, one for kernel number per spike, 11 for kernel weight and 9 for water status, which were not associated with differences in plant height or heading date. We detected significant correlations among traits and colocated QTL that were consistent with those correlations. Among those, grain yield and plant water status were negatively correlated in all environments, and six QTL for these traits were colocated or tightly linked (< 1 cM). QTL identified and validated in this study provide useful information for the improvement of fall-sown spring wheats under full and limited irrigation.

Gene editing and mutagenesis reveal inter-cultivar differences and additivity in the contribution of TaGW2 homoeologues to grain size and weight in wheat.

KEY MESSAGE: CRISPR-Cas9-based genome editing and EMS mutagenesis revealed inter-cultivar differences and additivity in the contribution of TaGW2 homoeologues to grain size and weight in wheat. The TaGW2 gene homoeologues have been reported to be negative regulators of grain size (GS) and thousand grain weight (TGW) in wheat. However, the contribution of each homoeologue to trait variation among different wheat cultivars is not well documented. We used the CRISPR-Cas9 system and TILLING to mutagenize each homoeologous gene copy in cultivars Bobwhite and Paragon, respectively. Plants carrying single-copy nonsense mutations in different genomes showed different levels of GS/TGW increase, with TGW increasing by an average of 5.5% (edited lines) and 5.3% (TILLING mutants). In any combination, the double homoeologue mutants showed higher phenotypic effects than the respective single-genome mutants. The double mutants had on average 12.1% (edited) and 10.5% (TILLING) higher TGW with respect to wild-type lines. The highest increase in GS and TGW was shown for triple mutants of both cultivars, with increases in 16.3% (edited) and 20.7% (TILLING) in TGW. The additive effects of the TaGW2 homoeologues were also demonstrated by the negative correlation between the functional gene copy number and GS/TGW in Bobwhite mutants and an F2 population. The highest single-genome increases in GS and TGW in Paragon and Bobwhite were obtained by mutations in the B and D genomes, respectively. These inter-cultivar differences in the phenotypic effects between the TaGW2 gene homoeologues coincide with inter-cultivar differences in the homoeologue expression levels. These results indicate that GS/TGW variation in wheat can be modulated by the dosage of homoeologous genes with inter-cultivar differences in the magnitude of the individual homoeologue effects.

Identification of the VERNALIZATION 4 gene reveals the origin of spring growth habit in ancient wheats from South Asia.

Wheat varieties with a winter growth habit require long exposures to low temperatures (vernalization) to accelerate flowering. Natural variation in four vernalization genes regulating this requirement has favored wheat adaptation to different environments. The first three genes (VRN1-VRN3) have been cloned and characterized before. Here we show that the fourth gene, VRN-D4, originated by the insertion of a ∼290-kb region from chromosome arm 5AL into the proximal region of chromosome arm 5DS. The inserted 5AL region includes a copy of VRN-A1 that carries distinctive mutations in its coding and regulatory regions. Three lines of evidence confirmed that this gene is VRN-D4: it cosegregated with VRN-D4 in a high-density mapping population; it was expressed earlier than other VRN1 genes in the absence of vernalization; and induced mutations in this gene resulted in delayed flowering. VRN-D4 was found in most accessions of the ancient subspecies Triticum aestivum ssp. sphaerococcum from South Asia. This subspecies showed a significant reduction of genetic diversity and increased genetic differentiation in the centromeric region of chromosome 5D, suggesting that VRN-D4 likely contributed to local adaptation and was favored by positive selection. Three adjacent SNPs in a regulatory region of the VRN-D4 first intron disrupt the binding of GLYCINE-RICH RNA-BINDING PROTEIN 2 (TaGRP2), a known repressor of VRN1 expression. The same SNPs were identified in VRN-A1 alleles previously associated with reduced vernalization requirement. These alleles can be used to modulate vernalization requirements and to develop wheat varieties better adapted to different or changing environments.

A whole-genome, radiation hybrid mapping resource of hexaploid wheat.

Generating a contiguous, ordered reference sequence of a complex genome such as hexaploid wheat (2n = 6x = 42; approximately 17 GB) is a challenging task due to its large, highly repetitive, and allopolyploid genome. In wheat, ordering of whole-genome or hierarchical shotgun sequencing contigs is primarily based on recombination and comparative genomics-based approaches. However, comparative genomics approaches are limited to syntenic inference and recombination is suppressed within the pericentromeric regions of wheat chromosomes, thus, precise ordering of physical maps and sequenced contigs across the whole-genome using these approaches is nearly impossible. We developed a whole-genome radiation hybrid (WGRH) resource and tested it by genotyping a set of 115 randomly selected lines on a high-density single nucleotide polymorphism (SNP) array. At the whole-genome level, 26 299 SNP markers were mapped on the RH panel and provided an average mapping resolution of approximately 248 Kb/cR1500 with a total map length of 6866 cR1500 . The 7296 unique mapping bins provided a five- to eight-fold higher resolution than genetic maps used in similar studies. Most strikingly, the RH map had uniform bin resolution across the entire chromosome(s), including pericentromeric regions. Our research provides a valuable and low-cost resource for anchoring and ordering sequenced BAC and next generation sequencing (NGS) contigs. The WGRH developed for reference wheat line Chinese Spring (CS-WGRH), will be useful for anchoring and ordering sequenced BAC and NGS based contigs for assembling a high-quality, reference sequence of hexaploid wheat. Additionally, this study provides an excellent model for developing similar resources for other polyploid species.