Shifting the limits in wheat research and breeding using a fully annotated reference genome

TaCER1-1A is involved in cuticular wax alkane biosynthesis in hexaploid wheat and responds to plant abiotic stresses

To protect above-ground plant organs from excessive water loss, their surfaces are coated by waxes. The genes involved in wax formation have been investigated in detail in Arabidopsis but scarcely in crop species. Here, we aimed to isolate and characterize a CER1 enzyme responsible for formation of the very long-chain alkanes present in high concentrations especially during late stages of wheat development. On the basis of comparative wax and transcriptome analyses of various wheat organs, we selected TaCER1-1A as a primary candidate and demonstrated that it was located to the endoplasmic reticulum, the subcellular compartment for wax biosynthesis. A wheat nullisomic-tetrasomic substitution line lacking TaCER1-1A had significantly reduced amounts of C₃₃ alkane, whereas rice plants overexpressing TaCER1-1A showed substantial increases of C₂₅ -C₃₃ alkanes relative to wild type control. Similarly, heterologous expression of TaCER1-1A in Arabidopsis wild type and the cer1 mutant resulted in increased levels of unbranched alkanes, iso-branched alkanes and alkenes. Finally, the expression of TaCER1-1A was found activated by abiotic stresses and abscisic acid treatment, resulting in increased production of alkanes in wheat. Taken together, our results demonstrate that TaCER1-1A plays an important role in wheat wax alkane biosynthesis and involved in responding to drought and other environmental stresses.

In-silico prediction of novel genes responsive to drought and salinity stress tolerance in bread wheat (Triticum aestivum)

Common wheat (Triticum aestivum) is the most widely grown cereal crop and is cultivated extensively in dry regions. Water shortage, resulting from either drought or salinity, leads to slow growth and loss of wheat yield. In order to predict new genes responsive to the drought and salt stresses in wheat, 6,717 expressed sequence tags (ESTs), expressed in drought and salinity stress conditions were collected from the National Center for Biotechnology Information (NCBI). The downloaded ESTs were clustered and assembled into 354 contigs; 14 transcription factor families in 29 contigs were identified. In addition, 119 contigs were organized in five enzyme classes. Biological functions were obtained for only 324 of the 354 contigs using gene ontology. In addition, using Kyoto Encyclopedia of Genes and Genomes database, 191 metabolic pathways were identified. The remaining contigs were used for further analysis and the search for new genes responsive to drought and salt stresses. These contigs were mapped on the International Wheat Genome Sequencing Consortium RefSeq v1.0 assembly, the most complete version of the reference sequence of the bread wheat variety Chinese Spring. They were found to have from one to three locations on the subgenomes A, B, and D. Full-length gene sequences were designed for these contigs, which were further validated using promoter analysis. These predicted genes may have applications in molecular breeding programs and wheat drought and salinity research.

Early Root Transcriptomic Changes in Wheat Seedlings Colonized by Trichoderma harzianum Under Different Inorganic Nitrogen Supplies

Wheat is one of the most important crops worldwide. The use of plant growth promoting microorganisms, such as those of the genus Trichoderma, constitutes an alternative to chemical fertilizers, since they are cheaper and are not detrimental to the environment. However, the interaction between Trichoderma and wheat plants has been scarcely studied, at least at a molecular level. In the present work, a microarray approach was used to study the early transcriptomic changes induced in wheat roots by Trichoderma harzianum, applied alone or in combination with different concentrations of calcium nitrate [Ca(NO₃)₂], which was last used as nitrogen (N) source. Our results show that T. harzianum causes larger transcriptomic changes than Ca(NO₃)₂ in wheat roots, and such changes are different when plants are challenged with Trichoderma alone or treated with a combination of T. harzianum and Ca(NO₃)₂. Overall, T. harzianum activates the expression of defense-related genes at early stages of the interaction with the roots, while this fungus reduces the expression of genes related to plant growth and development. Moreover, the current study in wheat roots, subjected to the different T. harzianum and Ca(NO₃)₂ combinations, reveals that the number of transcriptomic changes was higher when compared against those caused by the different Ca(NO₃)₂ concentrations than when it was compared against those caused by T. harzianum. N metabolism gene expression changes were in agreement with the levels of nitrate reductase activity measured in plants from Trichoderma plus Ca(NO₃)₂ conditions. Results were also concordant with plant phenotypes, which showed reduced growth at early interaction stages when inoculated with T. harzianum or with its combination with Ca(NO₃)₂ at the lowest dosage. These results were in a good agreement with the recognized role of Trichoderma as an inducer of plant defense.

Phenotypic and molecular characterization of Hessian fly resistance in diploid wheat, Aegilops tauschii

BACKGROUND

The Hessian fly (Mayetiola destructor), belonging to the gall midge family (Cecidomyiidae), is a devastating pest of wheat (Triticum aestivum) causing significant yield losses. Despite identification and characterization of numerous Hessian fly-responsive genes and associated biological pathways involved in wheat defense against this dipteran pest, their functional validation has been challenging. This is largely attributed to the large genome, polyploidy, repetitive DNA, and limited genetic resources in hexaploid wheat. The diploid progenitor Aegilops tauschii, D-genome donor of modern-day hexaploid wheat, offers an ideal surrogate eliminating the need to target all three homeologous chromosomes (A, B and D) individually, and thereby making the functional validation of candidate Hessian fly-responsive genes plausible. Furthermore, the well-annotated sequence of Ae. tauschii genome and availability of genetic resources amenable to manipulations makes the functional assays less tedious and time-consuming. However, prior to utilization of this diploid genome for downstream studies, it is imperative to characterize its physical and molecular responses to Hessian fly.

RESULTS

In this study we screened five Ae. tauschii accessions for their response to the Hessian fly biotypes L and vH13. Two lines were identified that exhibited a homozygous resistance response to feeding by both Hessian fly biotypes. Studies using physical measurements and neutral red staining showed that the resistant Ae. tauschii accessions resembled hexaploid wheat in their phenotypic responses to Hessian fly, that included similarities in larval developmental stages, leaf and plant growth, and cell wall permeability. Furthermore, molecular responses, characterized by gene expression profiling using quantitative real-time PCR, in select resistant Ae. tauschii lines also revealed similarities with resistant hexaploid wheat.

CONCLUSIONS

Phenotypic and molecular characterization of Ae. tauschii to Hessian fly infestation revealed resistant accessions that shared similarities to hexaploid wheat. Resembling the resistant hexaploid wheat, the Ae. tauschii accessions mount an early defense strategy involving defense proteins including lectins, secondary metabolites and reactive oxygen species (ROS) radicals. Our results reveal the suitability of the diploid progenitor for use as an ideal tool for functional genomics research in deciphering the wheat-Hessian fly molecular interactions.

Characterization of genetic diversity and population structure in wheat using array based SNP markers

Genetic diversity is crucial for successful adaptation and sustained improvement in crops. India is bestowed with diverse agro-climatic conditions which makes it rich in wheat germplasm adapted to various niches. Germplasm repository consists of local landraces, trait specific genetic stocks including introgressions from wild relatives, exotic collections, released varieties, and improved germplasm. Characterization of genetic diversity is done using morpho-physiological characters as well as by analyzing variations at DNA level. However, there are not many reports on array based high throughput SNP markers having characteristics of genome wide coverage employed in Indian spring wheat germplasm. Amongst wheat SNP arrays, 35K Axiom Wheat Breeder's Array has the highest SNP polymorphism efficiency suitable for genetic mapping and genetic diversity characterization. Therefore, genotyping was done using 35K in 483 wheat genotypes resulting in 14,650 quality filtered SNPs, that were distributed across the B (~ 50%), A (~ 39%), and D (~ 10%) genomes. The total genetic distance coverage was 4477.85 cM with 3.27 SNP/cM and 0.49 cM/SNP as average marker density and average inter-marker distance, respectively. The PIC ranged from 0.09 to 0.38 with an average of 0.29 across genomes. Population structure and Principal Coordinate Analysis resulted in two subpopulations (SP1 and SP2). The analysis of molecular variance revealed the genetic variation of 2% among and 98% within subpopulations indicating high gene flow between SP1 and SP2. The subpopulation SP2 showed high level of genetic diversity based on genetic diversity indices viz. Shannon's information index (I) = 0.648, expected heterozygosity (He) = 0.456 and unbiased expected heterozygosity (uHe) = 0.456. To the best of our knowledge, this study is the first to include the largest set of Indian wheat genotypes studied exclusively for genetic diversity. These findings may serve as a potential source for the identification of uncharacterized QTL/gene using genome wide association studies and marker assisted selection in wheat breeding programs.

A Co-Expression Network in Hexaploid Wheat Reveals Mostly Balanced Expression and Lack of Significant Gene Loss of Homeologous Meiotic Genes Upon Polyploidization

Polyploidization has played an important role in plant evolution. However, upon polyploidization, the process of meiosis must adapt to ensure the proper segregation of increased numbers of chromosomes to produce balanced gametes. It has been suggested that meiotic gene (MG) duplicates return to a single copy following whole genome duplication to stabilize the polyploid genome. Therefore, upon the polyploidization of wheat, a hexaploid species with three related (homeologous) genomes, the stabilization process may have involved rapid changes in content and expression of MGs on homeologous chromosomes (homeologs). To examine this hypothesis, sets of candidate MGs were identified in wheat using co-expression network analysis and orthology informed approaches. In total, 130 RNA-Seq samples from a range of tissues including wheat meiotic anthers were used to define co-expressed modules of genes. Three modules were significantly correlated with meiotic tissue samples but not with other tissue types. These modules were enriched for GO terms related to cell cycle, DNA replication, and chromatin modification and contained orthologs of known MGs. Overall, 74.4% of genes within these meiosis-related modules had three homeologous copies which was similar to other tissue-related modules. Amongst wheat MGs identified by orthology, rather than co-expression, the majority (93.7%) were either retained in hexaploid wheat at the same number of copies (78.4%) or increased in copy number (15.3%) compared to ancestral wheat species. Furthermore, genes within meiosis-related modules showed more balanced expression levels between homeologs than genes in non-meiosis-related modules. Taken together, our results do not support extensive gene loss nor changes in homeolog expression of MGs upon wheat polyploidization. The construction of the MG co-expression network allowed identification of hub genes and provided key targets for future studies.

Identification of quantitative trait loci for kernel traits in a wheat cultivar Chuannong16

Background

Kernel length (KL), kernel width (KW) and thousand-kernel weight (TKW) are key agronomic traits in wheat breeding. Chuannong16 (‘CN16’) is a commercial cultivar with significantly longer kernels than the line ‘20828’. To identify and characterize potential alleles from CN16 controlling KL, the previously developed recombinant inbred line (RIL) population derived from the cross ‘20828’ × ‘CN16’ and the genetic map constructed by the Wheat55K SNP array and SSR markers were used to perform quantitative trait locus/loci (QTL) analyses for kernel traits.

Results

A total of 11 putative QTL associated with kernel traits were identified and they were located on chromosomes 1A (2 QTL), 2B (2 QTL), 2D (3 QTL), 3D, 4A, 6A, and 7A, respectively. Among them, three major QTL, QKL.sicau-2D, QKW.sicau-2D and QTKW.sicau-2D, controlling KL, KW and TKW, respectively, were detected in three different environments. Respectively, they explained 10.88–18.85%, 17.21–21.49% and 10.01–23.20% of the phenotypic variance. Further, they were genetically mapped in the same interval on chromosome 2DS. A previously developed kompetitive allele-specific PCR (KASP) marker KASP-AX-94721936 was integrated in the genetic map and QTL re-mapping finally located the three major QTL in a 1- cM region flanked by AX-111096297 and KASP-AX-94721936. Another two co-located QTL intervals for KL and TKW were also identified. A few predicted genes involved in regulation of kernel growth and development were identified in the intervals of these identified QTL. Significant relationships between kernel traits and spikelet number per spike and anthesis date were detected and discussed.

Conclusions

Three major and stably expressed QTL associated with KL, KW, and TKW were identified. A KASP marker tightly linked to these three major QTL was integrated. These findings provide information for subsequent fine mapping and cloning the three co-localized major QTL for kernel traits.

Preparation and Characterization of Avenin-Enriched Oat Protein by Chill Precipitation for Feeding Trials in Celiac Disease

The safety of oats for people with celiac disease remains unresolved. While oats have attractive nutritional properties that can improve the quality and palatability of the restrictive, low fiber gluten-free diet, rigorous feeding studies to address their safety in celiac disease are needed. Assessing the oat prolamin proteins (avenins) in isolation and controlling for gluten contamination and other oat components such as fiber that can cause non-specific effects and symptoms is crucial. Further, the avenin should contain all reported immunogenic T cell epitopes, and be deliverable at a dose that enables biological responses to be correlated with clinical effects. To date, isolation of a purified food-grade avenin in sufficient quantities for feeding studies has not been feasible. Here, we report a new gluten isolation technique that enabled 2 kg of avenin to be extracted from 400 kg of wheat-free oats under rigorous gluten-free and food grade conditions. The extract consisted of 85% protein of which 96% of the protein was avenin. The concentration of starch (1.8% dry weight), β-glucan (0.2% dry weight), and free sugars (1.8% dry weight) were all low in the final avenin preparation. Other sugars including oligosaccharides, small fructans, and other complex sugars were also low at 2.8% dry weight. Liquid chromatography tandem mass spectrometry (LC-MS/MS) analysis of the proteins in these preparations showed they consisted only of oat proteins and were uncontaminated by gluten containing cereals including wheat, barley or rye. Proteomic analysis of the avenin enriched samples detected more avenin subtypes and fewer other proteins compared to samples obtained using other extraction procedures. The identified proteins represented five main groups, four containing known immune-stimulatory avenin peptides. All five groups were identified in the 50% (v/v) ethanol extract however the group harboring the epitope DQ2.5-ave-1b was less represented. The avenin-enriched protein fractions were quantitatively collected by reversed phase HPLC and analyzed by MALDI-TOF mass spectrometry. Three reverse phase HPLC peaks, representing ~40% of the protein content, were enriched in proteins containing DQ2.5-ave-1a epitope. The resultant high quality avenin will facilitate controlled and definitive feeding studies to establish the safety of oat consumption by people with celiac disease.

Genome-Wide Identification and Analysis of the AP2 Transcription Factor Gene Family in Wheat (Triticum aestivum L.)

The AP2 transcription factors play important roles in regulating plant growth and development. However, limited data are available on the contributions of AP2 transcription factors in wheat (Triticum aestivum L.). In the present study, a total of 62 AP2 genes were identified in wheat from a genome-wide search against the latest wheat genome data. Phylogenetic and sequence alignment analyses divided the wheat AP2 genes into 3 clusters, euAP2, euANT, and basalANT. Chromosomal distribution, gene structure and duplication, and motif composition were subsequently investigated. The 62 TaAP2 genes were unevenly distributed on 21 chromosomes. Twenty-four homologous gene sets among A, B, and D sub-genomes were detected, which contributed to the expansion of the wheat AP2 gene family. The expression levels of TaAP2 genes were examined using the WheatExp database; most detected genes exhibited tissue-specific expression patterns. The transcript levels of 9 randomly selected TaAP2 genes were validated through qPCR analyses. Overexpression of TaAP2-10-5D, the most likely homolog of Arabidopsis ANT gene, increased organ sizes in Arabidopsis. Our results extend our knowledge of the AP2 gene family in wheat, and contribute to further functional characterization of AP2s during wheat development with the ultimate goal of improving crop production.

Novel Informatic Tools to Support Functional Annotation of the Durum Wheat Genome

Seed mutagenesis is one strategy to create a population with thousands of useful mutations for the direct selection of desirable traits, to introduce diversity into varietal improvement programs, or to generate a mutant collection to support gene functional analysis. However, phenotyping such large collections, where each individual may carry many mutations, is a bottleneck for downstream analysis. Targeting Induced Local Lesions in Genomes (TILLinG), when coupled with next-generation sequencing allows high-throughput mutation discovery and selection by genotyping. We mutagenized an advanced durum breeding line, UAD0951096_F2:5 and performed short-read (2x125 bp) Illumina sequencing of the exome of 100 lines using an available exome capture platform. To improve variant calling, we generated a consolidated exome reference using the recently available genome sequences of the cultivars Svevo and Kronos to facilitate the alignment of reads from the UAD0951096_F2:5 derived mutants. The resulting exome reference was 484.4 Mbp. We also developed a user-friendly, searchable database and bioinformatic analysis pipeline that allowed us to predict zygosity of the mutations discovered and extracts flanking sequences for rapid marker development. Here, we present these tools with the aim of allowing researchers fast and accurate downstream selection of mutations discovered by TILLinG by sequencing to support functional annotation of the durum wheat genome.

Wheat AGAMOUS LIKE 6 transcription factors function in stamen development by regulating the expression of Ta APETALA3

Previous studies have revealed the functions of rice and maize AGAMOUS LIKE 6 (AGL6) genes OsMADS6 and ZAG3, respectively, in floral development; however, the functions of three wheat (Triticum aestivum) AGL6 genes are still unclear. Here, we report the main functions of wheat AGL6 homoeologous genes in stamen development. In RNAi plants, stamens showed abnormality in number and morphology, and a tendency to transform into carpels. Consistently, the expression of the B-class gene TaAPETALA3 (AP3) and the auxin-responsive gene TaMGH3 was downregulated, whereas the wheat ortholog of the rice carpel identity gene DROOPING LEAF was ectopically expressed in RNAi stamens. TaAGL6 proteins bind to the promoter of TaAP3 directly. Yeast one-hybrid and transient expression assays further showed that TaAGL6 positively regulates the expression of TaAP3 in vivo. Wheat AGL6 transcription factors interact with TaAP3, TaAGAMOUS and TaMADS13. Our findings indicate that TaAGL6 transcription factors play an essential role in stamen development through transcriptional regulation of TaAP3 and other related genes. We propose a model to illustrate the function and probable mechanism of this regulation. This study extends our understanding of AGL6 genes.

Quantitative trait loci for seeding root traits and the relationships between root and agronomic traits in common wheat

A completely developed and vigorous root system can provide a stable platform for aboveground plant organs. To identify loci controlling root traits that could be used in wheat (Triticum aestivum L.) breeding, 199 recombinant inbred lines were used to measure and analyze eight root traits. A total of 18 quantitative trait loci (QTL) located on chromosomes 1A, 2A, 2B, 2D, 4B, 4D, 6A, 7A, and 7B were identified. The phenotypic variation explained by these 18 QTL ranged from 3.27% to 11.75%, and the logarithm of odds scores ranged from 2.50 to 6.58. A comparison of physical intervals indicated several new QTL for root traits were identified. In addition, significant correlations between root and agronomic traits were detected and discussed. The results presented in this study, along with those of previous reports, suggest that chromosomes 2 and 7 likely play important roles in the growth and development of wheat roots.

Using environmental clustering to identify specific drought tolerance QTLs in bread wheat (T. aestivum L.)

Environmental clustering helps to identify QTLs associated with grain yield in different water stress scenarios. These QTLs could be useful for breeders to improve grain yields and increase genetic resilience in marginal environments. Drought is one of the main abiotic stresses limiting winter bread wheat growth and productivity around the world. The acquisition of new high-yielding and stress-tolerant varieties is therefore necessary and requires improved understanding of the physiological and genetic bases of drought resistance. A panel of 210 elite European varieties was evaluated in 35 field trials. Grain yield and its components were scored in each trial. A crop model was then run with detailed climatic data and soil water status to assess the dynamics of water stress in each environment. Varieties were registered from 1992 to 2011, allowing us to test timewise genetic progress. Finally, a genome-wide association study (GWAS) was carried out using genotyping data from a 280 K SNP chip. The crop model simulation allowed us to group the environments into four water stress scenarios: an optimal condition with no water stress, a post-anthesis water stress, a moderate-anthesis water stress and a high pre-anthesis water stress. Compared to the optimal water condition, grain yield losses in the stressed conditions were 3.3%, 12.4% and 31.2%, respectively. This environmental clustering improved understanding of the effect of drought on grain yields and explained 20% of the G × E interaction. The greatest genetic progress was obtained in the optimal condition, mostly represented in France. The GWAS identified several QTLs, some of which were specific of the different water stress patterns. Our results make breeding for improved drought resistance to specific environmental scenarios easier and will facilitate genetic progress in future environments, i.e., water stress environments.

A Roadmap for Lowering Crop Nitrogen Requirement

Increasing nitrogen fertilizer applications have sustained a growing world population in the 20th century. However, to avoid any further associated environmental damage, new sustainable agronomic practices together with new cultivars must be developed. To date the concept of nitrogen use efficiency (NUE) has been useful in quantifying the processes of nitrogen uptake and utilization, but we propose a shift in focus to consider nitrogen responsiveness as a more appropriate trait to select varieties with lower nitrogen requirements. We provide a roadmap to integrate the regulation of nitrogen uptake and assimilation into varietal selection and crop breeding programs. The overall goal is to reduce nitrogen inputs by farmers growing crops in contrasting cropping systems around the world, while sustaining yields and reducing greenhouse gas (GHG) emissions.

Plasticity in Triticeae centromere DNA sequences: a wheat × tall wheatgrass (decaploid) model

Centromeres mediate chromosome attachment to microtubules and maintain the integrity of chromosomes for proper segregation of the sister chromatids during cell division. Advances in the assembly of Triticeae genome sequences combined with the capacity to recover hybrid species derived from very distantly related species provides potential experimental systems for linking retrotransposon amplification and repositioning of centromeres via non-mendelian inheritance in partial amphiploid breeds. The decaploid tall wheatgrass (Thinopyrum ponticum) is one of the most successfully used perennial species in wheat breeding for generating translocation lines with valuable agronomic traits. We found that wheat centromere retrotransposons CRW and Quinta widely occur within the tall wheatgrass genome. In addition, one of the genome donors to Th. ponticum, Pseudoroegneria stipifolia (StSt), has been shown to have Abigail and a satellite repeat, CentSt. We also found two other centromeric retrotransposons, Abia and CL135 in Th. ponticum by ChIP-seq. Examination of partial amphiploid lines that were generated in the 1970s demonstrated extensive modification in centromere sequences using CentSt, Abigail and Abia as probes. We also detected that St-genome chromosomes were more enriched with Abigail and CentSt, whereas E-genome chromosomes were enriched with CRW and Quinta in tall wheatgrass and its closer relatives. It can be concluded that bursts of transposition of retrotransposons and repositioning of centromeres via non-mendelian segregation are common in partial amphiploids derived from interspecific hybrids. Practically speaking, our study reveals that the existence of homologous centromere functional sequences in both a donor and its receptor can substantially contribute to the successful transfer of alien genes into crop species. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/sra/SRR9089557; https://www.ncbi.nlm.nih.gov/sra/SRR9089558; https://www.ncbi.nlm.nih.gov/sra/SRR9089559; https://www.ncbi.nlm.nih.gov/sra/SRR9089560; https://www.ncbi.nlm.nih.gov/sra/SRR9089561; https://www.ncbi.nlm.nih.gov/sra/SRR9089562; https://www.ncbi.nlm.nih.gov/sra/SRR9089563; https://www.ncbi.nlm.nih.gov/sra/SRR9089564; https://www.ncbi.nlm.nih.gov/nuccore/MK999394; https://www.ncbi.nlm.nih.gov/nuccore/MK999395; https://www.ncbi.nlm.nih.gov/nuccore/MK999396.

Uncovering homeologous relationships between tetraploid Agropyron cristatum and bread wheat genomes using COS markers

Using COS markers, the study reveals homeologous relationships between tetraploid Agropyron cristatum and bread wheat to support alien introgression breeding of wheat. Crested wheatgrass (Agropyron cristatum L. Gaertn.) is a wild relative of wheat that possesses many genes that are potentially useful in wheat improvement. The species comprises a complex of diploid, tetraploid and hexaploid forms. In this study, wheat-A. cristatum chromosome, telosome and translocation lines were used to characterize syntenic relationships between tetraploid A. cristatum and bread wheat. Prior to mapping COS markers, the cytogenetic stock lines were characterized for fertility and by FISH and GISH for karyotype stability. Out of 328 COS markers selected for the study, 279 consistently amplified products in tetraploid A. cristatum, and, out of these, 139 were polymorphic between tetraploid crested wheatgrass and wheat. Sixty-nine markers were found to be suitable for the detection of tetraploid A. cristatum chromosomes 1P-6P in wheat, ranging from 6 to 17 markers per chromosome. BLASTn of the source ESTs resulted in significant hits for 67 markers on the wheat pseudomolecules. Generally, COS markers of the same homeologous group were detected on similar arms in both Agropyron and wheat. However, some intragenomic duplications and chromosome rearrangements were detected in tetraploid A. cristatum. These results provide new insights into the structure and evolution of the tetraploid A. cristatum genome and will facilitate the exploitation of the wild species for introgression breeding of bread wheat.

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.

Re-acquisition of the brittle rachis trait via a transposon insertion in domestication gene Q during wheat de-domestication

De-domestication is a unique evolutionary process during which crops re-acquire wild-like traits to survive and persist in agricultural fields without the need for human cultivation. The re-acquisition of seed dispersal mechanisms is crucial for crop de-domestication. Common wheat is an important cereal crop worldwide. Tibetan semi-wild wheat is a potential de-domesticated common wheat subspecies. However, the crucial genes responsible for its brittle rachis trait have not been identified. Genetic mapping, functional analyses and phylogenetic analyses were completed to identify the gene associated with Qbr.sau-5A, which is a major locus for the brittle rachis trait of Tibetan semi-wild wheat. The cloned Qbr.sau-5A gene is a new Q allele (Q^t ) with a 161-bp transposon insertion in exon 5. Although Q^t is expressed normally, its encoded peptide lacks some key features of the APETALA2 family. The abnormal functions of Q^t in developing wheat spikes result in brittle rachises. Phylogenetic and genotyping analyses confirmed that Q^t originated from Q in common wheat and is naturally distributed only in Tibetan semi-wild wheat populations. The identification of Q^t provides new evidence regarding the origin of Tibetan semi-wild wheat, and new insights into the re-acquisition of wild traits during crop de-domestication.

TaAPO-A1, an ortholog of rice ABERRANT PANICLE ORGANIZATION 1, is associated with total spikelet number per spike in elite European hexaploid winter wheat (Triticum aestivum L.) varieties

We dissected the genetic basis of total spikelet number (TSN) along with other traits, viz. spike length (SL) and flowering time (FT) in a panel of 518 elite European winter wheat varieties. Genome-wide association studies (GWAS) based on 39,908 SNP markers revealed highly significant quantitative trait loci (QTL) for TSN on chromosomes 2D, 7A, and 7B, for SL on 5A, and FT on 2D, with 2D-QTL being the functional marker for the gene Ppd-D1. The physical region of the 7A-QTL for TSN revealed the presence of a wheat ortholog (TaAPO-A1) to APO1–a rice gene that positively controls the spikelet number on the panicles. Interspecific analyses of the TaAPO-A1 orthologs showed that it is a highly conserved gene important for floral development and present in a wide range of terrestrial plants. Intraspecific studies of the TaAPO-A1 across wheat genotypes revealed a polymorphism in the conserved F-box domain, defining two haplotypes. A KASP marker developed on the polymorphic site showed a highly significant association of TaAPO-A1 with TSN, explaining 23.2% of the total genotypic variance. Also, the TaAPO-A1 alleles showed weak but significant differences for SL and grain yield. Our results demonstrate the importance of wheat sequence resources to identify candidate genes for important traits based on genetic analyses.

Population Genomic Approaches for Weed Science

Genomic approaches are opening avenues for understanding all aspects of biological life, especially as they begin to be applied to multiple individuals and populations. However, these approaches typically depend on the availability of a sequenced genome for the species of interest. While the number of genomes being sequenced is exploding, one group that has lagged behind are weeds. Although the power of genomic approaches for weed science has been recognized, what is needed to implement these approaches is unfamiliar to many weed scientists. In this review we attempt to address this problem by providing a primer on genome sequencing and provide examples of how genomics can help answer key questions in weed science such as: (1) Where do agricultural weeds come from; (2) what genes underlie herbicide resistance; and, more speculatively, (3) can we alter weed populations to make them easier to control? This review is intended as an introduction to orient weed scientists who are thinking about initiating genome sequencing projects to better understand weed populations, to highlight recent publications that illustrate the potential for these methods, and to provide direction to key tools and literature that will facilitate the development and execution of weed genomic projects.

Genome wide genetic dissection of wheat quality and yield related traits and their relationship with grain shape and size traits in an elite × non-adapted bread wheat cross

The genetic gain in yield and quality are two major targets of wheat breeding programs around the world. In this study, a high density genetic map consisting of 10,172 SNP markers identified a total of 43 genomic regions associated with three quality traits, three yield traits and two agronomic traits in hard red spring wheat (HRSW). When compared with six grain shape and size traits, the quality traits showed mostly independent genetic control (~18% common loci), while the yield traits showed moderate association (~53% common loci). Association of genomic regions for grain area (GA) and thousand-grain weight (TGW), with yield suggests that targeting an increase in GA may help enhancing wheat yield through an increase in TGW. Flour extraction (FE), although has a weak positive phenotypic association with grain shape and size, they do not share any common genetic loci. A major contributor to plant height was the Rht8 locus and the reduced height allele was associated with significant increase in grains per spike (GPS) and FE, and decrease in number of spikes per square meter and test weight. Stable loci were identified for almost all the traits. However, we could not find any QTL in the region of major known genes like GPC-B1, Ha, Rht-1, and Ppd-1. Epistasis also played an important role in the genetics of majority of the traits. In addition to enhancing our knowledge about the association of wheat quality and yield with grain shape and size, this study provides novel loci, genetic information and pre-breeding material (combining positive alleles from both parents) to enhance the cultivated gene pool in wheat germplasm. These resources are valuable in facilitating molecular breeding for wheat quality and yield improvement.

The Use of Pentaploid Crosses for the Introgression of Amblyopyrum muticum and D-Genome Chromosome Segments Into Durum Wheat

The wild relatives of wheat provide an important source of genetic variation for wheat improvement. Much of the work in the past aimed at transferring genetic variation from wild relatives into wheat has relied on the exploitation of the ph1b mutant, located on the long arm of chromosome 5B. This mutation allows homologous recombination to occur between chromosomes from related but different genomes, e.g. between the chromosomes of wheat and related chromosomes from a wild relative resulting in the generation of interspecific recombinant chromosomes. However, the ph1b mutant also enables recombination to occur between the homologous genomes of wheat, e.g. A/B, A/D, B/D, resulting in the generation of wheat intergenomic recombinant chromosomes. In this work we report on the presence of wheat intergenomic recombinants in the genomic background of hexaploid wheat/Amblyopyrum muticum introgression lines. The transfer of genomic rearrangements involving the D-genome through pentaploid crosses provides a strategy by which the D-genome of wheat can be introgressed into durum wheat. Hence, a pentaploid crossing strategy was used to transfer D-genome segments, introgressed with either the A- and/or the B-genome, into the tetraploid background of two durum wheat genotypes Karim and Om Rabi 5 in either the presence or absence of different Am. muticum (2n = 2x = 14, TT) introgressions. Introgressions were monitored in backcross generations to the durum wheat parents via multi-color genomic in situ hybridization (mc-GISH). Tetraploid lines carrying homozygous D-genome introgressions, as well as simultaneous homozygous D- and T-genome introgressions, were developed. Introgression lines were characterized via Kompetitive Allele-Specific PCR (KASP) markers and multi-color fluorescence in situ hybridization (FISH). Results showed that new wheat sub-genomic translocations were generated at each generation in progeny that carried any Am. muticum chromosome introgression irrespective of the linkage group that the segment was derived from. The highest frequencies of homologous recombination were observed between the A- and the D-genomes. Results indicated that the genotype Karim had a higher tolerance to genomic rearrangements and T-genome introgressions compared to Om Rabi 5. This indicates the importance of the selection of the parental genotype when attempting to transfer/develop introgressions into durum wheat from pentaploid crosses.

Conserved residues in the wheat (Triticum aestivum) NAM-A1 NAC domain are required for protein binding and when mutated lead to delayed peduncle and flag leaf senescence

BACKGROUND

NAC transcription factors contain five highly conserved subdomains which are required for protein dimerisation and DNA binding. Few residues within these subdomains have been identified as essential for protein function, and fewer still have been shown to be of biological relevance in planta. Here we use a positive regulator of senescence in wheat, NAM-A1, to test the impact of missense mutations at specific, highly conserved residues of the NAC domain on protein function.

RESULTS

We identified missense mutations in five highly conserved residues of the NAC domain of NAM-A1 in a tetraploid TILLING population. TILLING lines containing these mutations, alongside synonymous and non-conserved mutation controls, were grown under glasshouse conditions and scored for senescence. Four of the five mutations showed a significant and consistent delay in peduncle senescence but had no consistent effects on flag leaf senescence. All four mutant alleles with the delayed senescence phenotype also lost the ability to interact with the homoeolog NAM-B1 in a yeast two-hybrid assay. Two of these residues were previously shown to be involved in NAC domain function in Arabidopsis, suggesting conservation of residue function between species. Three of these four alleles led to an attenuated cell death response compared to wild-type NAM-A1 when transiently over-expressed in Nicotiana benthamiana. One of these mutations was further tested under field conditions, in which there was a significant and consistent delay in both peduncle and leaf senescence.

CONCLUSIONS

We combined field and glasshouse studies of a series of mutant alleles with biochemical analyses to identify four residues of the NAC domain which are required for NAM-A1 function and protein interaction. We show that mutations in these residues lead to a gradient of phenotypes, raising the possibility of developing allelic series of mutations for traits of agronomic importance. We also show that mutations in NAM-A1 more severely impact peduncle senescence, compared to the more commonly studied flag leaf senescence, highlighting this as an area deserving of further study. The results from this integrated approach provide strong evidence that conserved residues within the functional domains of NAC transcription factors have biological significance in planta.

Conserved residues in the wheat (Triticum aestivum) NAM-A1 NAC domain are required for protein binding and when mutated lead to delayed peduncle and flag leaf senescence

BACKGROUND

NAC transcription factors contain five highly conserved subdomains which are required for protein dimerisation and DNA binding. Few residues within these subdomains have been identified as essential for protein function, and fewer still have been shown to be of biological relevance in planta. Here we use a positive regulator of senescence in wheat, NAM-A1, to test the impact of missense mutations at specific, highly conserved residues of the NAC domain on protein function.

RESULTS

We identified missense mutations in five highly conserved residues of the NAC domain of NAM-A1 in a tetraploid TILLING population. TILLING lines containing these mutations, alongside synonymous and non-conserved mutation controls, were grown under glasshouse conditions and scored for senescence. Four of the five mutations showed a significant and consistent delay in peduncle senescence but had no consistent effects on flag leaf senescence. All four mutant alleles with the delayed senescence phenotype also lost the ability to interact with the homoeolog NAM-B1 in a yeast two-hybrid assay. Two of these residues were previously shown to be involved in NAC domain function in Arabidopsis, suggesting conservation of residue function between species. Three of these four alleles led to an attenuated cell death response compared to wild-type NAM-A1 when transiently over-expressed in Nicotiana benthamiana. One of these mutations was further tested under field conditions, in which there was a significant and consistent delay in both peduncle and leaf senescence.

CONCLUSIONS

We combined field and glasshouse studies of a series of mutant alleles with biochemical analyses to identify four residues of the NAC domain which are required for NAM-A1 function and protein interaction. We show that mutations in these residues lead to a gradient of phenotypes, raising the possibility of developing allelic series of mutations for traits of agronomic importance. We also show that mutations in NAM-A1 more severely impact peduncle senescence, compared to the more commonly studied flag leaf senescence, highlighting this as an area deserving of further study. The results from this integrated approach provide strong evidence that conserved residues within the functional domains of NAC transcription factors have biological significance in planta.

Dynamic QTL for adult plant resistance to powdery mildew in common wheat (Triticum aestivum L.)

Agriculture will benefit from a rigorous characterization of genes for adult plant resistance (APR) since this gene class was recognized to provide more durable protection from plant diseases. The present study reports the identification of APR loci to powdery mildew in German winter wheat cultivars Cortez and Atlantis. Cortez was previously shown to carry all-stage resistance gene Pm3e. To avoid interference of Pm3e in APR studies, line 6037 that lacked Pm3e but showed field resistance from doubled-haploid (DH) population Atlantis/Cortez was used in two backcrosses to Atlantis for the establishment of DH population 6037/Atlantis//Atlantis. APR was assessed in the greenhouse 10, 15, and 20 days after inoculation (dai) from the 4-leaf stage onwards and combined with single-nucleotide polymorphism data in a genome-wide association study (GWAS) and a linkage map-based quantitative trait loci (QTL) analysis. In GWAS, two QTL were detected: one on chromosome 1BL 10 dai, the other on chromosome 2BL 20 dai. In conventional QTL analysis, both QTL were detected with all three disease ratings: the QTL on chromosome 1BL explained a maximum of 35.2% of the phenotypic variation 10 dai, whereas the QTL on chromosome 2BL explained a maximum of 43.5% of the phenotypic variation 20 dai. Compared with GWAS, linkage map-based QTL analysis allowed following the dynamics of QTL action. The two large-effect QTL for APR to powdery mildew with dynamic gene action can be useful for the enhancement of wheat germplasm.

Genome-wide identification, characterization and expression pattern analysis of APYRASE family members in response to abiotic and biotic stresses in wheat

APYRASEs, which directly regulate intra- and extra-cellular ATP homeostasis, play a pivotal role in the regulation of various stress adaptations in mammals, bacteria and plants. In the present study, we identified and characterized wheat APYRASE family members at the genomic level in wheat. The results identified a total of nine APY homologs with conserved ACR domains. The sequence alignments, phylogenetic relations and conserved motifs of wheat APYs were bioinformatically analyzed. Although they share highly conserved secondary and tertiary structures, the wheat APYs could be mainly categorized into three groups, according to phylogenetic and structural analysis. Additionally, these APYs exhibited similar expression patterns in the root and shoot, among which TaAPY3-1, TaAPY3-3 and TaAPY3-4 had the highest expression levels. The time-course expression patterns of the eight APYs in response to biotic and abiotic stress in the wheat seedlings were also investigated. TaAPY3-2, TaAPY3-3, TaAPY3-4 and TaAPY6 exhibited strong sensitivity to all kinds of stresses in the leaves. Some APYs showed specific expression responses, such as TaAPY6 to heavy metal stress, and TaAPY7 to heat and salt stress. These results suggest that the stress-inducible APYs could have potential roles in the regulation of environmental stress adaptations. Moreover, the catalytic activity of TaAPY3-1 was further analyzed in the in vitro system. The results showed that TaAPY3-1 protein exhibited high catalytic activity in the degradation of ATP and ADP, but with low activity in degradation of TTP and GTP. It also has an extensive range of temperature adaptability, but preferred relatively acidic pH conditions. In this study, the genome-wide identification and characterization of APYs in wheat were suggested to be useful for further genetic modifications in the generation of high-stress-tolerant wheat cultivars.

Genome-Wide Identification and Expression Profiles of Late Embryogenesis-Abundant (LEA) Genes during Grain Maturation in Wheat (Triticum aestivum L.)

Late embryogenesis-abundant (LEA) genes play important roles in plant growth and development, especially the cellular dehydration tolerance during seed maturation. In order to comprehensively understand the roles of LEA family members in wheat, we carried out a series of analyses based on the latest genome sequence of the bread wheat Chinese Spring. 121 Triticum aestivum L. LEA (TaLEA) genes, classified as 8 groups, were identified and characterized. TaLEA genes are distributed in all chromosomes, most of them with a low number of introns (≤3). Expression profiles showed that most TaLEA genes expressed specifically in grains. By qRT-PCR analysis, we confirmed that 12 genes among them showed high expression levels during late stage grain maturation in two spring wheat cultivars, Yangmai16 and Yangmai15. For most genes, the peak of expression appeared earlier in Yangmai16. Statistical analysis indicated that expression level of 8 genes in Yangmai 16 were significantly higher than Yangmai 15 at 25 days after anthesis. Taken together, our results provide more knowledge for future functional analysis and potential utilization of TaLEA genes in wheat breeding.

Wheat Seed Proteins: Factors Influencing Their Content, Composition, and Technological Properties, and Strategies to Reduce Adverse Reactions

Wheat is the primary source of nutrition for many, especially those living in developing countries, and wheat proteins are among the most widely consumed dietary proteins in the world. However, concerns about disorders related to the consumption of wheat and/or wheat gluten proteins have increased sharply in the last 20 years. This review focuses on wheat gluten proteins and amylase trypsin inhibitors, which are considered to be responsible for eliciting most of the intestinal and extraintestinal symptoms experienced by susceptible individuals. Although several approaches have been proposed to reduce the exposure to gluten or immunogenic peptides resulting from its digestion, none have proven sufficiently effective for general use in coeliac-safe diets. Potential approaches to manipulate the content, composition, and technological properties of wheat proteins are therefore discussed, as well as the effects of using gluten isolates in various food systems. Finally, some aspects of the use of gluten-free commodities are discussed.

The Highly Conserved Barley Powdery Mildew Effector BEC1019 Confers Susceptibility to Biotrophic and Necrotrophic Pathogens in Wheat

Effector proteins secreted by plant pathogens play important roles in promoting colonization. Blumeria effector candidate (BEC) 1019, a highly conserved metalloprotease of Blumeria graminis f. sp. hordei (Bgh), is essential for fungal haustorium formation, and silencing BEC1019 significantly reduces Bgh virulence. In this study, we found that BEC1019 homologs in B. graminis f. sp. tritici (Bgt) and Gaeumannomyces graminis var. tritici (Ggt) have complete sequence identity with those in Bgh, prompting us to investigate their functions. Transcript levels of BEC1019 were abundantly induced concomitant with haustorium formation in Bgt and necrosis development in Ggt-infected plants. BEC1019 overexpression considerably increased wheat susceptibility to Bgt and Ggt, whereas silencing this gene using host-induced gene silencing significantly enhanced wheat resistance to Bgt and Ggt, which was associated with hydrogen peroxide accumulation, cell death, and pathogenesis-related gene expression. Additionally, we found that the full and partial sequences of BEC1019 can trigger cell death in Nicotiana benthamiana leaves. These results indicate that Bgt and Ggt can utilize BEC1019 as a virulence effector to promote plant colonization, and thus these genes represent promising new targets in breeding wheat cultivars with broad-spectrum resistance.

An ancestral NB-LRR with duplicated 3'UTRs confers stripe rust resistance in wheat and barley

Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is a global threat to wheat production. Aegilops tauschii, one of the wheat progenitors, carries the YrAS2388 locus for resistance to Pst on chromosome 4DS. We reveal that YrAS2388 encodes a typical nucleotide oligomerization domain-like receptor (NLR). The Pst-resistant allele YrAS2388R has duplicated 3’ untranslated regions and is characterized by alternative splicing in the nucleotide-binding domain. Mutation of the YrAS2388R allele disrupts its resistance to Pst in synthetic hexaploid wheat; transgenic plants with YrAS2388R show resistance to eleven Pst races in common wheat and one race of P. striiformis f. sp. hordei in barley. The YrAS2388R allele occurs only in Ae. tauschii and the Ae. tauschii-derived synthetic wheat; it is absent in 100% (n = 461) of common wheat lines tested. The cloning of YrAS2388R will facilitate breeding for stripe rust resistance in wheat and other Triticeae species.

Stripe rust is a serious threat to wheat production. Here, the authors reveal that the resistance gene, only present in the wheat progenitor Aegilops tauschii and its derived synthetic wheat, encodes a nucleotide oligomerization domain-like receptor and confers resistance in common wheat and barley.

Genetic Screening for Mutants with Altered Seminal Root Numbers in Hexaploid Wheat Using a High-Throughput Root Phenotyping Platform

Roots are the main channel for water and nutrient uptake in plants. Optimization of root architecture provides a viable strategy to improve nutrient and water uptake efficiency and maintain crop productivity under water-limiting and nutrient-poor conditions. We know little, however, about the genetic control of root development in wheat, a crop supplying 20% of global calorie and protein intake. To improve our understanding of the genetic control of seminal root development in wheat, we conducted a high-throughput screen for variation in seminal root number using an exome-sequenced mutant population derived from the hexaploid wheat cultivar Cadenza. The screen identified seven independent mutants with homozygous and stably altered seminal root number phenotypes. One mutant, Cadenza0900, displays a recessive extra seminal root number phenotype, while six mutants (Cadenza0062, Cadenza0369, Cadenza0393, Cadenza0465, Cadenza0818 and Cadenza1273) show lower seminal root number phenotypes most likely originating from defects in the formation and activation of seminal root primordia. Segregation analysis in F2 populations suggest that the phenotype of Cadenza0900 is controlled by multiple loci whereas the Cadenza0062 phenotype fits a 3:1 mutant:wild-type segregation ratio characteristic of dominant single gene action. This work highlights the potential to use the sequenced wheat mutant population as a forward genetic resource to uncover novel variation in agronomic traits, such as seminal root architecture.

Application of image-based phenotyping tools to identify QTL for in-field winter survival of winter wheat (Triticum aestivum L.)

Genome-wide association on winter survival was conducted using data from image-based phenotyping method. Nine QTL were observed and three of them with candidate gene identified. Winter survival is an essential trait of winter wheat (Triticum aestivum L.) grown in regions with high risk of winterkill. We characterized a diversity panel of 450 Canadian wheat varieties that included mostly winter-growth habit wheats to identify key genetic factors that contribute to higher winter survival under field conditions. To more accurately quantify winter survival differences among varieties, image-based phenotyping methods, captured by unmanned aerial vehicle (UAV) and on ground level, were used to estimate the winter survival of each varieties. Winter survival index was developed to correct for emergence when evaluating winter survival. Winter survival measurement estimated by visual estimation, UAV imagery and ground imagery showed strong correlation with each other and had comparable broad-sense heritability. Genome-wide association studies resulted in the identification of seven quantitative trait loci (QTL) for winter survival including Vrn-A1. By using the recently released annotated sequence of the wheat genome and the available RNA-Seq data, two putative candidate genes underlying the QTL for winter survival were identified. However, our study showed that certain QTL was unique to specific winter survival measurement. Collectively, our study demonstrated the feasibility of using UAV-based imagery for the identification of loci associated with winter survival in wheat. The complexity of in-field condition make our result a valuable complement to indoor frost-tolerance studies in the identification of genetic factors not directly linked to freezing tolerance.

Chromosome engineering-mediated introgression and molecular mapping of novel Aegilops speltoides-derived resistance genes for tan spot and Septoria nodorum blotch diseases in wheat

We identified, mapped and introduced novel Aegilops speltoides-derived resistance genes for tan spot and SNB diseases into wheat, enhancing understanding and utilization of host resistance to both diseases in wheat. Tan spot and Septoria nodorum blotch (SNB) are two important fungal diseases of wheat. Resistance to these diseases is often observed as the lack of sensitivity to the necrotrophic effectors (NE) produced by the fungal pathogens and thus exhibits a recessive inheritance pattern. In this study, we identified novel genes for resistance to tan spot and SNB on Aegilops speltoides (2n = 2x = 14, genome SS) chromosome 2S. These genes confer dominant resistance in the wheat background, indicating a distinct NE-independent mechanism of resistance. Ae. speltoides chromosome 2S was engineered for resistance gene introgression and molecular mapping by inducing meiotic homoeologous recombination with wheat chromosome 2B. Twenty representative 2B-2S recombinants were evaluated for reaction to tan spot and SNB and were delineated by genomic in situ hybridization and high-throughput wheat 90 K SNP assay. The resistance genes physically mapped to the sub-telomeric region (~ 8 Mb) on the short arm of chromosome 2S and designated TsrAes1 for tan spot resistance and SnbAes1 for SNB resistance. In addition, we developed SNP-derived PCR markers closely linked to TsrAes1/SnbAes1 for marker-assisted selection in wheat breeding. TsrAes1 and SnbAes1 are the first set of NE-independent tan spot, and SNB resistance genes are identified from Ae. speltoides. The 2SS-2BS·2BL recombinants with minimal amounts of Ae. speltoides chromatin containing TsrAes1/SnbAes1 were produced for germplasm development, making the wild species-derived resistance genes usable in wheat breeding. This will strengthen and diversify resistance of wheat to tan spot and SNB and facilitate understanding of resistance to these two diseases.

Amelioration in gliadin antigenicity and maintenance of viscoelastic properties of wheat (Triticum aestivum L.) cultivars with mixed probiotic fermentation

Sourdough fermentation of twenty wheat cultivars was carried out using mixed probiotic culture (Lactobacillus acidophilus UNI, Lactobacillus brevis LR/5 and Lactobacillus plantarum ATCC 8014). The gliadin antigenicity was expressed in terms of its content in twenty different wheat cultivars. The gliadin proteins were characterized by SDS-PAGE and structural changes analyzed on FTIR spectrophotometer. Moreover, changes in the viscoelastic character of fermented and non-fermented dough were studied by rheometry. The results showed a remarkable reduction in antigenicity by 60% (average) in all wheat cultivars on sourdough fermentation. This reduction may be due to the synergistic effect of protease secretion by mixed lactobacilli, responsible for gliadin degradation. These changes in gliadins by mixed culture proteolysis were confirmed on SDS-PAGE on observing new gliadin-derived low molecular weight peptides. The results were further validated by FTIR spectroscopy where structural changes of gliadins were analyzed in the fermented dough. The rheological data indicated a higher storage modulus (G') compared to loss modulus (G″) in both control and fermented flour of all wheat cultivars, however, with a lower efficacy in sourdoughs. The present study thus establishes that mixed culture sourdough fermentation decreases the antigenic potential of gliadins without any change in the rheology and thereby maintaining the baking or viscoelastic properties of the wheat flour.

Genetic dissection of drought and heat-responsive agronomic traits in wheat

High yield and wide adaptation are principal targets of wheat breeding but are hindered by limited knowledge on genetic basis of agronomic traits and abiotic stress tolerances. In this study, 277 wheat accessions were phenotyped across 30 environments with non-stress, drought-stressed, heat-stressed, and drought-heat-stressed treatments and were subjected to genome-wide association study using 395 681 single nucleotide polymorphisms. We detected 295 associated loci including consistent loci for agronomic traits across different treatments and eurytopic loci for multiple abiotic stress tolerances. A total of 22 loci overlapped with quantitative trait loci identified by biparental quantitative trait loci mapping. Six loci were simultaneously associated with agronomic traits and abiotic stress tolerance, four of which fell within selective sweep regions. Selection in Chinese wheat has increased the frequency of superior marker alleles controlling yield-related traits in the four loci during past decades, which conversely diminished favourable genetic variation controlling abiotic stress tolerance in the same loci; two promising candidate paralogous genes colocalized with such loci, thereby providing potential targets for studying the molecular mechanism of stress tolerance-productivity trade-off. These results uncovering promising alleles controlling agronomic traits and/or multiple abiotic stress tolerances, providing insights into heritable covariation between yield and abiotic stress tolerance, will accelerate future efforts for wheat improvement.

Characterization of Pm65, a new powdery mildew resistance gene on chromosome 2AL of a facultative wheat cultivar

A new powdery mildew resistance gene that can be readily used in wheat breeding, Pm65, was identified in the facultative wheat cultivar Xinmai 208 and mapped to the terminal region of chromosome 2AL. Wheat powdery mildew, a widely occurring disease caused by the biotrophic fungus Blumeriagraminis f. sp. tritici (Bgt), poses a serious threat to wheat production. A high breeding priority is to identify powdery mildew resistance genes that can be readily used either alone or in gene complexes involving other disease resistance genes. An F₂ population and 227 F_2:3 families derived from the cross Xinmai 208 × Stardust were generated to map a powdery mildew resistance gene in Xinmai 208, a high-yielding Chinese wheat cultivar. Genetic analysis indicated that Xinmai 208 carries a single dominant powdery mildew resistance gene, designated herein Pm65, and linkage analysis delimited Pm65 to a 0.5 cM interval covering 531.8 Kb (763,289,667-763,821,463 bp) on chromosome 2AL in the Chinese Spring reference sequence. An allelism test indicated that Pm65 is a new gene about 10.3 cM distal to the Pm4 locus. Pm65 was 0.3 cM proximal to Xstars355 and 0.2 cM distal to Xstars356. It conferred near-immunity to 19 of 20 Bgt isolates collected from different wheat-growing regions of the USA. Coming from a high-yield potential cultivar, Pm65 can be directly used to enhance powdery mildew resistance in wheat. The newly developed SSR markers Xstars355 and Xstars356 have the potential to tag Pm65 for wheat improvement.

Genomic Prediction and Genome-Wide Association Studies of Flour Yield and Alveograph Quality Traits Using Advanced Winter Wheat Breeding Material

Use of genetic markers and genomic prediction might improve genetic gain for quality traits in wheat breeding programs. Here, flour yield and Alveograph quality traits were inspected in 635 F6 winter wheat breeding lines from two breeding cycles. Genome-wide association studies revealed single nucleotide polymorphisms (SNPs) on chromosome 5D significantly associated with flour yield, Alveograph P (dough tenacity), and Alveograph W (dough strength). Additionally, SNPs on chromosome 1D were associated with Alveograph P and W, SNPs on chromosome 1B were associated with Alveograph P, and SNPs on chromosome 4A were associated with Alveograph L (dough extensibility). Predictive abilities based on genomic best linear unbiased prediction (GBLUP) models ranged from 0.50 for flour yield to 0.79 for Alveograph W based on a leave-one-out cross-validation strategy. Predictive abilities were negatively affected by smaller training set sizes, lower genetic relationship between lines in training and validation sets, and by genotype-environment (G×E) interactions. Bayesian Power Lasso models and genomic feature models resulted in similar or slightly improved predictions compared to GBLUP models. SNPs with the largest effects can be used for screening large numbers of lines in early generations in breeding programs to select lines that potentially have good quality traits. In later generations, genomic predictions might be used for a more accurate selection of high quality wheat lines.

Genetic Contribution of Synthetic Hexaploid Wheat to CIMMYT's Spring Bread Wheat Breeding Germplasm

Synthetic hexaploid (SH) wheat (AABBD'D') is developed by artificially generating a fertile hybrid between tetraploid durum wheat (Triticum turgidum, AABB) and diploid wild goat grass (Aegilops tauschii, D'D'). Over three decades, the International Maize and Wheat Improvement Center (CIMMYT) has developed and utilized SH wheat to bridge gene transfer from Ae. tauschii and durum wheat to hexaploid bread wheat. This is a unique example of success utilizing wild relatives in mainstream breeding at large scale worldwide. Our study aimed to determine the genetic contribution of SH wheat to CIMMYT's global spring bread wheat breeding program. We estimated the theoretical and empirical contribution of D' to synthetic derivative lines using the ancestral pedigree and marker information using over 1,600 advanced lines and their parents. The average marker-estimated D' contribution was 17.5% with difference in genome segments suggesting application of differential selection pressure. The pedigree-based contribution was correlated with marker-based estimates without providing chromosome segment specific variation. Results from international yield trials showed that 20% of the lines were synthetic derived with an average D' contribution of 15.6%. Our results underline the importance of SH wheat in maintaining and enhancing genetic diversity and genetic gain over years and is important for development of a more targeted introgression strategy. The study provides retrospective view into development and utilization of SH in the CIMMYT Global Wheat Program.

Genetic diversity and population structure analysis of synthetic and bread wheat accessions in Western Siberia

Recurrent selection and intercrossing between best of the best parents in each generation of breeding cycle resulted in a narrower genetic diversity in elite wheat (Triticum aestivum L.) germplasm. Therefore, we investigated diverse source of 143 synthetic and bread wheat accessions for identifying potentially rich genetic resources for improving the genetic diversity in wheat. This study identified 47,526 genotyping-by-sequencing-derived SNP markers that were nearly evenly distributed across three genomes of wheat. The population structure analysis identified three distinct clusters (Japan synthetics, CIMMYT synthetics, and bread wheat) of wheat genotypes on the basis of type and geographical origin of wheat accessions. Population differentiation using analysis of molecular variance indicated 21% of the total genetic variance among subgroups and the remainder within subgroups. This study also identified that the Japan synthetic group was the most divergent group compared with other subgroups. The genetic diversity comparisons between synthetic and bread wheat lines showed that the gene diversity of synthetic wheat was 33% higher than bread wheat accessions, indicating the potential use of these lines for broadening the genetic diversity of modern wheat cultivars. The results from this study will be helpful in further understanding genomic features of wheat and facilitate their use in wheat breeding programs.

Genome-wide association study reveals new loci for yield-related traits in Sichuan wheat germplasm under stripe rust stress

BACKGROUND

As one of the most important food crops in the world, increasing wheat (Triticum aestivum L.) yield is an urgent task for global food security under the continuous threat of stripe rust (caused by Puccinia striiformis f. sp. tritici) in many regions of the world. Molecular marker-assisted breeding is one of the most efficient ways to increase yield. Here, we identified loci associated to multi-environmental yield-related traits under stripe rust stress in 244 wheat accessions from Sichuan Province through genome-wide association study (GWAS) using 44,059 polymorphic markers from the 55 K single nucleotide polymorphism (SNP) chip.

RESULTS

A total of 13 stable quantitative trait loci (QTLs) were found to be highly associating to yield-related traits, including 6 for spike length (SL), 3 for thousand-kernel weight (TKW), 2 for kernel weight per spike (KWPS), and 2 for both TKW and KWPS, in at least two test environments under stripe rust stress conditions. Of them, ten QTLs were overlapped or very close to the reported QTLs, three QTLs, QSL.sicau-1AL, QTKW.sicau-4AL, and QKWPS.sicau-4AL.1, were potentially novel through the physical location comparison with previous QTLs. Further, 21 candidate genes within three potentially novel QTLs were identified, they were mainly involved in the regulation of phytohormone, cell division and proliferation, meristem development, plant or organ development, and carbohydrate transport.

CONCLUSIONS

QTLs and candidate genes detected in our study for yield-related traits under stripe rust stress will facilitate elucidating genetic basis of yield-related trait and could be used in marker-assisted selection in wheat yield breeding.

Serpins: Genome-Wide Characterisation and Expression Analysis of the Serine Protease Inhibitor Family in Triticum aestivum

The serine protease inhibitor (serpin) gene family is the largest family of protease inhibitors. Serine protease inhibitors have an active, but under-characterized, role in grain development and defense against pathogen attack in cereal crops. By exploiting publicly available genomic, transcriptomic and proteomic data for wheat (Triticum aestivum), we have identified and annotated the entire ’serpinome’ of wheat and constructed a high-quality and robust phylogenetic tree of the gene family, identifying paralogous and homeologous clades from the hexaploid wheat genome, including the Serpin-Z group that have been well characterized in barley. Using publicly available RNAseq data (http://www.wheat-expression.com/), expression profiles of the wheat serpins were explored across a variety of tissues from the developing grain, spikelet and spike. We show that the SERPIN-Z clade, among others, are highly expressed during grain development, and that there is homeologous and paralogous functional redundancy in this gene family. Further to their role in grain development, serpins play an important but under-explored role in response to fungal pathogens. Using 13 RNAseq datasets of wheat tissues infected by fungal pathogens, we identified 37 serpins with a significant disease response. The majority of the disease-responsive serpins were upregulated by Fusarium graminearum, a destructive fungal pathogen that attacks the spike and developing grain of wheat. As serpins are ubiquitous in wheat grain, the genes encoding serpins may be linked to grain development, with their disease response a result of pleiotropy.

isqg: A Binary Framework for in Silico Quantitative Genetics

The dna is the fundamental basis of genetic information, just as bits are for computers. Whenever computers are used to represent genetic data, the computational encoding must be efficient to allow the representation of processes driving the inheritance and variability. This is especially important across simulations in view of the increasing complexity and dimensions brought by genomics. This paper introduces a new binary representation of genetic information. Algorithms as bitwise operations that mimic the inheritance of a wide range of polymorphisms are also presented. Different kinds and mixtures of polymorphisms are discussed and exemplified. Proposed algorithms and data structures were implemented in C++ programming language and is available to end users in the R package “isqg” which is available at the R repository (cran). Supplementary data are available online.

Genetic Mapping and Prediction Analysis of FHB Resistance in a Hard Red Spring Wheat Breeding Population

Fusarium head blight (FHB) is one of the most destructive diseases in wheat worldwide. Breeding for FHB resistance is hampered by its complex genetic architecture, large genotype by environment interaction, and high cost of phenotype screening. Genomic selection (GS) is a powerful tool to enhance improvement of complex traits such as FHB resistance. The objectives of this study were to (1) investigate the genetic architecture of FHB resistance in a North Dakota State University (NDSU) hard red spring wheat breeding population, (2) test if the major QTL Fhb1 and Fhb5 play an important role in this breeding population; and (3) assess the potential of GS to enhance breeding efficiency of FHB resistance. A total of 439 elite spring wheat breeding lines from six breeding cycles were genotyped using genotyping-by-sequencing (GBS) and 102,147 SNP markers were obtained. Evaluation of FHB severity was conducted in 10 unbalanced field trials across multiple years and locations. One QTL for FHB resistance was identified and located on chromosome arm 1AL, explaining 5.3% of total phenotypic variation. The major type II resistance QTL Fhb1 only explained 3.1% of total phenotypic variation and the QTL Fhb5 was not significantly associated with FHB resistance in this breeding population. Our results suggest that integration of many genes with medium/minor effects in this breeding population should provide stable FHB resistance. Genomic prediction accuracies of 0.22-0.44 were obtained when predicting over breeding cycles in this study, indicating the potential of GS to enhance the improvement of FHB resistance.

Chromosome Location Contributing to Ozone Tolerance in Wheat

Breeding wheat for higher grain yield can contribute to global food security and sustainable production on less land. Tropospheric ozone can injure wheat plants and subsequently reduce grain yield. Identification of ozone tolerance in the wheat genome can assist plant breeders in developing new sources of tolerant germplasm. Our objective was to use the ‘Chinese Spring’ monosomic lines to screen for ozone response and identify the chromosomic locations contributing to ozone tolerance based on foliar injury. Two methodologies, Continuous Stirred Tank Reactors and Outdoor Plant Environment Chambers, were used to expose wheat monosomic lines to varying concentrations and durations of ozone. Each wheat monosomic line in ‘Chinese Spring’ has a missing chromosome in each of the wheat subgenomes (A, B, and D). In both methodologies, we found significant and repeatable data to identify chromosome 7A as a major contributor to tolerance to ozone injury in ‘Chinese Spring’. In every experiment, the absence of chromosome 7A resulted in significant injury to wheat due to ozone. This was not the case when any other chromosome was missing.

A breeding strategy targeting the secondary gene pool of bread wheat: introgression from a synthetic hexaploid wheat

Introgressing one-eighth of synthetic hexaploid wheat genome through a double top-cross plus a two-phase selection is an effective strategy to develop high-yielding wheat varieties. The continued expansion of the world population and the likely onset of climate change combine to form a major crop breeding challenge. Genetic advances in most crop species to date have largely relied on recombination and reassortment within a relatively narrow gene pool. Here, we demonstrate an efficient wheat breeding strategy for improving yield potentials by introgression of multiple genomic regions of de novo synthesized wheat. The method relies on an initial double top-cross (DTC), in which one parent is synthetic hexaploid wheat (SHW), followed by a two-phase selection procedure. A genotypic analysis of three varieties (Shumai 580, Shumai 969 and Shumai 830) released from this program showed that each harbors a unique set of genomic regions inherited from the SHW parent. The first two varieties were generated from very small populations, whereas the third used a more conventional scale of selection since one of bread wheat parents was a pre-breeding material. The three varieties had remarkably enhanced yield potential compared to those developed by conventional breeding. A widely accepted consensus among crop breeders holds that introducing unadapted germplasm, such as landraces, as parents into a breeding program is a risky proposition, since the size of the breeding population required to overcome linkage drag becomes too daunting. However, the success of the proposed DTC strategy has demonstrated that novel variation harbored by SHWs can be accessed in a straightforward, effective manner. The strategy is in principle generalizable to any allopolyploid crop species where the identity of the progenitor species is known.

Wheat ms5 male-sterility is induced by recessive homoeologous A and D genome non-specific lipid transfer proteins

Nuclear male-sterile mutants with non-conditional, recessive and strictly monogenic inheritance are useful for both hybrid and conventional breeding systems, and have long been a research focus for many crops. In allohexaploid wheat, however, genic redundancy results in rarity of such mutants, with the ethyl methanesulfonate-induced mutant ms5 among the few reported to date. Here, we identify TaMs5 as a glycosylphosphatidylinositol-anchored lipid transfer protein required for normal pollen exine development, and by transgenic complementation demonstrate that TaMs5-A restores fertility to ms5. We show ms5 locates to a centromere-proximal interval and has a sterility inheritance pattern modulated by TaMs5-D but not TaMs5-B. We describe two allelic forms of TaMs5-D, one of which is non-functional and confers mono-factorial inheritance of sterility. The second form is functional but shows incomplete dominance. Consistent with reduced functionality, transcript abundance in developing anthers was found to be lower for TaMs5-D than TaMs5-A. At the 3B homoeolocus, we found only non-functional alleles among 178 diverse hexaploid and tetraploid wheats that include landraces and Triticum dicoccoides. Apparent ubiquity of non-functional TaMs5-B alleles suggests loss-of-function arose early in wheat evolution and, therefore, at most knockout of two homoeoloci is required for sterility. This work provides genetic information, resources and tools required for successful implementation of ms5 sterility in breeding systems for bread and durum wheats.

The anther-specific CYP704B is potentially responsible for MSG26 male sterility in barley

Plant male sterility is a valuable trait in breeding and hybrid seed production. The barley male-sterility gene msg26 was mapped to a 0.02-cM region that anchors to a 506-kb low-quality assembly between two cleaved amplified polymorphic sequence (CAPS) markers, SP1M14 and SP1M49. The barley gene HORVU4Hr1G074840, which encodes a putative cytochrome P450 CYP704B protein, appears to be a strong candidate for the MSG26 trait. Barley (Hordeum vulgare L.) is an important cereal crop worldwide. Traditional breeding in barley is time-consuming and labor-intensive. The use of male-sterile genotypes may significantly improve the efficacy of hybrid breeding and seed production. The barley accession 'GSHO745' is a spontaneous male-sterile mutant from the barley variety, 'Unitan'. The male sterility in 'GSHO745' is controlled by the recessive gene, msg26 (originally named as ms-u). We revealed that the barley plants homozygous for msg26 proceeded normally through Meiosis II until the tetrad stage, but became fully defective in the late uninucleate microspores and developed pollen-less anthers. Using seven barley F₂ populations, we mapped MSG26 to a 0.02-cM region that anchored to a 506-kb low-quality assembly between two cleaved amplified polymorphic sequence markers, SP1M14 and SP1M49. The HORVU4Hr1G074840 gene that encodes a putative cytochrome P450 protein (CYP704B) was identified as the most plausible candidate for MSG26. First, HORVU4Hr1G074840 is located in a collinear region of the rice CYP704B2 and the maize CYP704B1. Both of these genes are essential for male gamete production. Second, the male-sterile allele of HORVU4Hr1G074840 in GSHO745 contained a 4-bp deletion in the last exon. The resulting frame shift causes a Gly436Gln substitution, scrambles the sequence of the remainder of the protein, and forms a new termination site at the 70th triplet of the shifted reading frame. We thus called the variant protein CYP704B:p.G436Qfs*70. Third, the barley HORVU4Hr1G074840 gene was specifically expressed in anthers. Altogether, HORVU4Hr1G074840 represents a strong candidate for MSG26 in barley.

Cell surface immune receptors: the guardians of the plant's extracellular spaces

Since the original 'Zigzag model', several iterations have been proposed to reconcile both the Pattern Triggered Immunity (PTI) and the Effector Triggered Immunity (ETI) branches of the plant immune system. The recent cloning of new disease resistance genes, functioning in gene-for-gene interactions, which structurally resemble cell surface broad spectrum Pattern Recognition Receptors, have further blurred the distinctions between PTI and ETI in plant immunity. In an attempt to simplify further the existing conceptual models, we, herein, propose a scheme based on the spatial localization of the key proteins (receptors) which function to induce plant immune responses. We believe this 'Spatial Invasion model' will prove useful for understanding how immune receptors interact with different pathogen types which peripherally or totally invade plant cells, colonize solely extracellularly or switch locations during a successful infection.