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

Bread wheat: a role model for plant domestication and breeding

BACKGROUND

Bread wheat is one of the most important crops in the world. Its domestication coincides with the beginning of agriculture and since then, it has been constantly under selection by humans. Its breeding has followed millennia of cultivation, sometimes with unintended selection on adaptive traits, and later by applying intentional but empirical selective pressures. For more than one century, wheat breeding has been based on science, and has been constantly evolving due to on farm agronomy and breeding program improvements. The aim of this work is to briefly review wheat breeding, with emphasis on the current advances.

DISCUSSION

Improving yield potential, resistance/tolerance to biotic and abiotic stresses, and baking quality, have been priorities for breeding this cereal, however, new objectives are arising, such as biofortification enhancement. The narrow genetic diversity and complexity of its genome have hampered the breeding progress and the application of biotechnology. Old approaches, such as the introgression from relative species, mutagenesis, and hybrid breeding are strongly reappearing, motivated by an accumulation of knowledge and new technologies. A revolution has taken place regarding the use of molecular markers whereby thousands of plants can be routinely genotyped for thousands of loci. After 13 years, the wheat reference genome sequence and annotation has finally been completed, and is currently available to the scientific community. Transgenics, an unusual approach for wheat improvement, still represents a potential tool, however it is being replaced by gene editing, whose technology along with genomic selection, speed breeding, and high-throughput phenotyping make up the most recent frontiers for future wheat improvement.

FINAL CONSIDERATION

Agriculture and plant breeding are constantly evolving, wheat has played a major role in these processes and will continue through decades to come.

iTRAQ-Based Quantitative Analysis of Responsive Proteins Under PEG-Induced Drought Stress in Wheat Leaves

Drought is an important abiotic stress that seriously restricts crop productivity. An understanding of drought tolerance mechanisms offers guidance for cultivar improvement. In order to understand how a well-known wheat genotype Jinmai 47 responds to drought, we adopted the iTRAQ and LC/MS approaches and conducted proteomics analysis of leaves after exposure to 20% of polyethylene glycol-6000 (PEG)-induced stress for 4 days. The study identified 176 differentially expressed proteins (DEPs), with 65 (36.5%) of them being up-regulated, and 111 (63.5%) down-regulated. DEPs, located in cellular membranes and cytosol mainly, were involved in stress and redox regulation (51), carbohydrate and energy metabolism (36), amino acid metabolism (24), and biosynthesis of other secondary metabolites (20) primarily. Under drought stress, TCA cycle related proteins were up-regulated. Antioxidant system, signaling system, and nucleic acid metabolism etc. were relatively weakened. In comparison, the metabolism pathways that function in plasma dehydration protection and protein structure protection were strongly enhanced, as indicated by the improved biosynthesis of 2 osmolytes, sucrose and Proline, and strongly up-regulated protective proteins, LEA proteins and chaperones. SUS4, P5CSs, OAT, Rab protein, and Lea14-A were considered to be important candidate proteins, which deserve to be further investigated.

Integrative Analysis of the Wheat PHT1 Gene Family Reveals A Novel Member Involved in Arbuscular Mycorrhizal Phosphate Transport and Immunity

Phosphorus (P) deficiency is one of the main growth-limiting factors for plants. However, arbuscular mycorrhizal (AM) symbiosis can significantly promote P uptake. Generally, PHT1 transporters play key roles in plants' P uptake, and thus, PHT1 genes have been investigated in some plants, but the regulation and functions of these genes in wheat (TaPHT1) during AM symbiosis have not been studied in depth. Therefore, a comprehensive analysis of TaPHT1 genes was performed, including sequence, phylogeny, cis-elements, expression, subcellular localization and functions, to elucidate their roles in AM-associated phosphate transport and immunity. In total, 35 TaPHT1s were identified in the latest high-quality bread wheat genome, 34 of which were unevenly distributed on 13 chromosomes, and divided into five groups. Sequence analysis indicated that there are 11 types of motif architectures and five types of exon-intron structures in the TaPHT1 family. Duplication mode analysis indicated that the TaPHT1 family has expanded mainly through segmental and tandem duplication events, and that all duplicated gene pairs have been under purifying selection. Transcription analysis of the 35 TaPHT1s revealed that not only known the mycorrhizal-specific genes TaPht-myc, TaPT15-4B (TaPT11) and TaPT19-4D (TaPT10), but also four novel mycorrhizal-specific/inducible genes (TaPT3-2D, TaPT11-4A, TaPT29-6A, and TaPT31-7A) are highly up-regulated in AM wheat roots. Furthermore, the mycorrhizal-specific/inducible genes are significantly induced in wheat roots at different stages of infection by colonizing fungi. Transient Agrobacterium tumefaciens-mediated transformation expression in onion epidermal cells showed that TaPT29-6A is a membrane-localized protein. In contrast to other AM-specific/inducible PHT1 genes, TaPT29-6A is apparently required for the symbiotic and direct Pi pathway. TaPT29-6A-silenced lines exhibited reduced levels of AM fungal colonization and arbuscules, but increased susceptibility to biotrophic, hemi-biotrophic and necrotrophic pathogens. In conclusion, TaPT29-6A was not only essential for the AM symbiosis, but also played vital roles in immunity.

Identification and mapping of expressed genes associated with the 2DL QTL for fusarium head blight resistance in the wheat line Wuhan 1

Background

Fusarium head blight (FHB) is a problem of great concern in small grain cereals, especially wheat. A quantitative trait locus (QTL) for FHB resistance (FHB_SFI) located on the long arm of chromosome 2D in the spring wheat genotype Wuhan 1 is a resistance locus which has potential to improve the FHB resistance of bread wheat since it confers effective resistance to wheat breeding lines. Recently, differentially expressed genes (DEG) have been identified by comparing near isogenic lines (NIL) carrying the susceptible and resistant alleles for the 2DL QTL, using RNA-Seq. In the present study, we aimed to identify candidate genes located within the genetic interval for the 2DL QTL for FHB resistance, as assessed by single floret inoculation (FHB_SFI), and possibly contributing to it.

Results

Combining previous and additional bioinformatics analyses, 26 DEG that were located on chromosome arm 2DL were selected for further characterization of their expression profile by RT-qPCR. Seven of those DEG showed a consistent differential expression profile between either three pairs of near isogenic lines or other genotypes carrying the R and S alleles for the 2DL QTL for FHB resistance. UN25696, which was identified in previous expression work using microarray was also confirmed to have a differential expression pattern. Those eight candidate genes were further characterized in 85 lines of a double haploid mapping population derived from the cross Wuhan 1/Nyubai, the population where the 2DL QTL was originally identified. The expression QTL for gene Traes_2DL_179570792 overlapped completely with the mapping interval for the 2DL QTL for FHB_SFI while the expression QTL for UN25696 mapped near the QTL, but did not overlap with it. None of the other genes had a significant eQTL on chromosome 2DL. Higher expression of Traes_2DL_179570792 and UN25696 was associated with the resistant allele at that locus.

Conclusions

Of the 26 DEG from the 2DL chromosome further characterized in this study, only two had an expression QTL located in or near the interval for the 2DL QTL. Traes_2DL_179570792 is the first expression marker identified as associated with the 2DL QTL.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0748-6) contains supplementary material, which is available to authorized users.

Global transcriptome analysis uncovers the gene co-expression regulation network and key genes involved in grain development of wheat (Triticum aestivum L.)

Wheat grain development is a robust biological process that largely determines grain quality and yield. In this study, we investigated the grain transcriptome of winter wheat cv. Xiaoyan-6 at four developmental stages (5, 10, 15, and 20 days post-anthesis), using high-throughput RNA sequencing (RNA-Seq). We identified 427 grain-specific transcription factors (TFs) and 1653 differentially expressed TFs during grain development as well as a grain co-expression regulation network (GrainNet) of the TFs and their predicted co-expressed genes. Our study identified ten putative key TFs and the predicted regulatory genes of these TFs in wheat grain development of Xiaoyan-6. The analysis was given a firm basis through the study of additional wheat tissues, including root, stem, leaf, flag leaf, grain, spikes (from wheat plants at booting or heading stages) to provide a dataset of 92,478 high-confidence protein-coding genes that were mostly evenly distributed among subgenomes, but unevenly distributed across each of the chromosomes or each of the seven homeologous groups. Within this larger framework of the transcriptomes, we identified 4659 grain-specific genes (SEGs) and 26,500 differentially expressed genes (DEGs) throughout grain development stages tested. The SEGs identified mainly associate with regulation and signaling-related biological processes, while the DEGs mainly involve in cellular component organization or biogenesis and nutrient reservoir activity during grain development of Xiaoyan-6. This study establishes new targets for modifying genes related to grain development and yield, to fine-tune expression in different varieties and environments.

The Dark Matter of Large Cereal Genomes: Long Tandem Repeats

Reference genomes of important cereals, including barley, emmer wheat and bread wheat, were released recently. Their comparison with genome size estimates obtained by flow cytometry indicated that the assemblies represent not more than 88-98% of the complete genome. This work is aimed at identifying the missing parts in two cereal genomes and proposing techniques to make the assemblies more complete. We focused on tandemly organised repetitive sequences, known to be underrepresented in genome assemblies generated from short-read sequence data. Our study found arrays of three tandem repeats with unit sizes of 1242 to 2726 bp present in the bread wheat reference genome generated from short reads. However, this and another wheat genome assembly employing long PacBio reads failed in integrating correctly the 2726-bp repeat in the pseudomolecule context. This suggests that tandem repeats of this size, frequently incorporated in unassigned scaffolds, may contribute to shrinking of pseudomolecules without reducing size of the entire assembly. We demonstrate how this missing information may be added to the pseudomolecules with the aid of nanopore sequencing of individual BAC clones and optical mapping. Using the latter technique, we identified and localised a 470-kb long array of 45S ribosomal DNA absent from the reference genome of barley.

Transcriptional reference map of hormone responses in wheat spikes

Background

Phytohormones are key regulators of plant growth, development, and signalling networks involved in responses to diverse biotic and abiotic stresses. Transcriptional reference maps of hormone responses have been reported for several model plant species such as Arabidopsis thaliana, Oryza sativa, and Brachypodium distachyon. However, because of species differences and the complexity of the wheat genome, these transcriptome data are not appropriate reference material for wheat studies.

Results

We comprehensively analysed the transcriptomic responses in wheat spikes to seven phytohormones, including indole acetic acid (IAA), gibberellic acid (GA), abscisic acid (ABA), ethylene (ET), cytokinin (CK), salicylic acid (SA), and methyl jasmonic acid (MeJA). A total of 3386 genes were differentially expressed at 24 h after the hormone treatments. Furthermore, 22.7% of these genes exhibited overlapping transcriptional responses for at least two hormones, implying there is crosstalk among phytohormones. We subsequently identified genes with expression levels that were significantly and differentially induced by a specific phytohormone (i.e., hormone-specific responses). The data for these hormone-responsive genes were then compared with the transcriptome data for wheat spikes exposed to biotic (Fusarium head blight) and abiotic (water deficit) stresses.

Conclusion

Our data were used to develop a transcriptional reference map of hormone responses in wheat spikes.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5726-x) contains supplementary material, which is available to authorized users.

Fine Mapping of the Wheat Leaf Rust Resistance Gene Lr42

Leaf rust caused by Puccinia triticina Eriks is one of the most problematic diseases of wheat throughout the world. The gene Lr42 confers effective resistance against leaf rust at both seedling and adult plant stages. Previous studies had reported Lr42 to be both recessive and dominant in hexaploid wheat; however, in diploid Aegilops tauschii (TA2450), we found Lr42 to be dominant by studying segregation in two independent F₂ and their F_2:3 populations. We further fine-mapped Lr42 in hexaploid wheat using a KS93U50/Morocco F₅ recombinant inbred line (RIL) population to a 3.7 cM genetic interval flanked by markers TC387992 and WMC432. The 3.7 cM Lr42 region physically corresponds to a 3.16 Mb genomic region on chromosome 1DS based on the Chinese Spring reference genome (RefSeq v.1.1) and a 3.5 Mb genomic interval on chromosome 1 in the Ae. tauschii reference genome. This region includes nine nucleotide-binding domain leucine-rich repeat (NLR) genes in wheat and seven in Ae. tauschii, respectively, and these are the likely candidates for Lr42. Furthermore, we developed two kompetitive allele-specific polymorphism (KASP) markers (SNP113325 and TC387992) flanking Lr42 to facilitate marker-assisted selection for rust resistance in wheat breeding programs.

Evolving Therapy for Celiac Disease

Gluten is known to be the main triggering factor for celiac disease (CeD), an immune-mediated disorder. CeD is therefore managed using a strict and lifelong gluten-free diet (GFD), the only effective treatment available currently. However, the GFD is restrictive. Hence, efforts are being made to explore alternative therapies. Based on their mechanisms of action on various molecular targets involved in the pathogenesis of CeD, these therapies may be classified into one of the following five broad approaches. The first approach focuses on decreasing the immunogenic content of gluten, using strategies like genetically modified wheat, intra-intestinal gluten digestion using glutenases, microwave thermal treatment of hydrated wheat kernels, and gluten pretreatment with either bacterial/ fungal derived endopeptidases or microbial transglutaminase. The second approach involves sequestering gluten in the gut lumen before it is digested into immunogenic peptides and absorbed, using binder drugs like polymer p(HEMA-co-SS), single chain fragment variable (scFv), and anti- gluten antibody AGY. The third approach aims to prevent uptake of digested gluten through intestinal epithelial tight junctions, using a zonulin antagonist. The fourth approach involves tissue transglutaminase (tTG) inhibitors to prevent the enhancement of immunogenicity of digested gluten by the intestinal tTG enzyme. The fifth approach seeks to prevent downstream immune activation after uptake of gluten immunogenic peptides through the intestinal mucosal epithelial layer. Examples include HLA-DQ2 blockers that prevent presentation of gluten derived- antigens by dendritic cells to T cells, immune- tolerizing therapies like the vaccine Nexvax2 and TIMP-Glia, cathepsin inhibitors, immunosuppressants like corticosteroids, azathioprine etc., and anti-cytokine agents targeting TNF-α and interleukin-15. Apart from these approaches, research is being done to evaluate the effectiveness of probiotics/prebiotics, helminth therapy using Necator americanus, low FODMAP diet, and pancreatic enzyme supplementation in CeD symptom control; however, the mechanisms by which they play a beneficial role in CeD are yet to be clearly established. Overall, although many therapies being explored are still in the pre-clinical phase, some like the zonulin antagonist, immune tolerizing therapies and glutenases have reached phase II/III clinical trials. While these potential options appear exciting, currently they may at best be used to supplement rather than supplant the GFD.

Genetic Dissection of Resistance to the Three Fungal Plant Pathogens Blumeria graminis, Zymoseptoria tritici, and Pyrenophora tritici-repentis Using a Multiparental Winter Wheat Population

Bread wheat (Triticum aestivum L.) is one of the world’s most important crop species. The development of new varieties resistant to multiple pathogens is an ongoing task in wheat breeding, especially in times of increasing demand for sustainable agricultural practices. Despite this, little is known about the relations between various fungal disease resistances at the genetic level, and the possible consequences for wheat breeding strategies. As a first step to fill this gap, we analyzed the genetic relations of resistance to the three fungal diseases – powdery mildew (PM), septoria tritici blotch (STB), and tan spot (TS) – using a winter wheat multiparent advanced generation intercross population. Six, seven, and nine QTL for resistance to PM, STB, and TS, respectively, were genetically mapped. Additionally, 15 QTL were identified for the three agro-morphological traits plant height, ear emergence time, and leaf angle distribution. Our results suggest that resistance to STB and TS on chromosome 2B is conferred by the same genetic region. Furthermore, we identified two genetic regions on chromosome 1AS and 7AL, which are associated with all three diseases, but not always in a synchronal manner. Based on our results, we conclude that parallel marker-assisted breeding for resistance to the fungal diseases PM, STB, and TS appears feasible. Knowledge of the genetic co-localization of alleles with contrasting effects for different diseases, such as on chromosome 7AL, allows the trade-offs of selection of these regions to be better understood, and ultimately determined at the genic level.

Development and characterisation of interspecific hybrid lines with genome-wide introgressions from Triticum timopheevii in a hexaploid wheat background

BACKGROUND

Triticum timopheevii (2n = 4x = 28; AtAtGG), is an important source for new genetic variation for wheat improvement with genes for potential disease resistance and salt tolerance. By generating a range of interspecific hybrid lines, T. timopheevii can contribute to wheat's narrow gene-pool and be practically utilised in wheat breeding programmes. Previous studies that have generated such introgression lines between wheat and its wild relatives have been unable to use high-throughput methods to detect the presence of wild relative segments in such lines.

RESULTS

A whole genome introgression approach, exploiting homoeologous recombination in the absence of the Ph1 locus, has resulted in the transfer of different chromosome segments from both the At and G genomes of T. timopheevii into wheat. These introgressions have been detected and characterised using single nucleotide polymorphism (SNP) markers present on a high-throughput Axiom® Genotyping Array. The analysis of these interspecific hybrid lines has resulted in the detection of 276 putative unique introgressions from T. timopheevii, thereby allowing the generation of a genetic map of T. timopheevii containing 1582 SNP markers, spread across 14 linkage groups representing each of the seven chromosomes of the At and G genomes of T. timopheevii. The genotyping of the hybrid lines was validated through fluorescence in situ hybridisation (FISH). Comparative analysis of the genetic map of T. timopheevii and the physical map of the hexaploid wheat genome showed that synteny between the two species is highly conserved at the macro-level and confirmed the presence of inter- and intra-genomic translocations within the At and G genomes of T. timopheevii that have been previously only detected through cytological techniques.

CONCLUSIONS

In this work, we report a set of SNP markers present on a high-throughput genotyping array, able to detect the presence of T. timopheevii in a hexaploid wheat background making it a potentially valuable tool for marker assisted selection (MAS) in wheat pre-breeding programs. These valuable resources of high-density molecular markers and wheat-T. timopheevii hybrid lines will greatly enhance the work being undertaken for wheat improvement through wild relative introgressions.

Comparative proteomic analysis of multi-ovary wheat under heterogeneous cytoplasm suppression

BACKGROUND

DUOII is a multi-ovary wheat (Triticum aestivum L.) line with two or three pistils and three stamens in each floret. The multi-ovary trait of DUOII is controlled by a dominant gene, whose expression can be suppressed by the heterogeneous cytoplasm of TeZhiI (TZI), a line with the nucleus of common wheat and the cytoplasm of Aegilops. Crosses between female DUOII plants and male TZI plants resulted in multi-ovary F₁s; whereas, the reciprocal crosses resulted in mono-ovary F₁s. Although the multi-ovary trait is inherited as single trait controlled by a dominant allele in lines with a Triticum cytoplasm, the mechanism by which the special heterogeneous cytoplasm suppresses the expression of multi-ovary is not well understood.

RESULTS

Observing the developmental process, we found that the critical stage of additional pistil primordium development was when the young spikes were 2-6 mm long. Then, we compared the quantitative proteomic profiles of 2-6 mm long young spikes obtained from the reciprocal crosses between DUOII and TZI. A total of 90 differentially expressed proteins were identified and analyzed based on their biological functions. These proteins had obvious functional pathways mainly implicated in chloroplast metabolism, nuclear and cell division, plant respiration, protein metabolism, and flower development. Importantly, we identified two key proteins, Flowering Locus K Homology Domain and PEPPER, which are known to play an essential role in the specification of pistil organ identity. By drawing relationships between the 90 differentially expressed proteins, we found that these proteins revealed a complex network which is associated with multi-ovary gene expression under heterogeneous cytoplasmic suppression.

CONCLUSIONS

Our proteomic analysis has identified certain differentially expressed proteins in 2-6 mm long young spikes, which was the critical stage of additional primordium development. This paper provided a universal proteomic profiling involved in the cytoplasmic suppression of wheat floral meristems; and our findings have laid a solid foundation for further mechanistic studies on the underlying mechanisms that control the heterogeneous cytoplasm-induced suppression of the nuclear multi-ovary gene in wheat.

Broadening the bread wheat D genome

Although Ae. tauschii has been extensively utilised for wheat breeding, the D-genome-containing allopolyploids have largely remained unexploited. In this review, we discuss approaches that can be used to exploit the D genomes of the different Aegilops species for the improvement of bread wheat. The D genome of allohexaploid bread wheat (Triticum aestivum, 2n = AABBDD) is the least diverse of the three wheat genomes and is unarguably less diverse than that of diploid progenitor Aegilops tauschii (2n = DD). Useful genetic variation and phenotypic traits also exist within each of the wheat group species containing a copy of the D genome: allopolyploid Aegilops species Ae. cylindrica (2n = D^cD^cC^cC^c), Ae. crassa 4x (2n = D¹D¹X^crX^cr), Ae. crassa 6x (2n = D¹D¹X^crX^crD^crD^cr), Ae. ventricosa (2n = D^vD^vN^vN^v), Ae. vavilovii (2n = D¹D¹X^crX^crS^vS^v) and Ae. juvenalis (2n = D¹D¹X^crX^crU^jU^j). Although Ae. tauschii has been extensively utilised for wheat breeding, the D-genome-containing allopolyploids have largely remained unexploited. Some of these D genomes appear to be modified relative to the bread wheat and Ae. tauschii D genomes, and others present in the allopolyploids may also contain useful variation as a result of adaptation to an allopolyploid, multi-genome environment. We summarise the genetic relationships, karyotypic variation and phenotypic traits known to be present in each of the D genome species that could be of relevance for bread wheat improvement and discuss approaches that can be used to exploit the D genomes of the different Aegilops species for the improvement of bread wheat. Better understanding of factors controlling chromosome inheritance and recombination in wheat group interspecific hybrids, as well as effective utilisation of new and developing genetics and genomics technologies, have great potential to improve the agronomic potential of the bread wheat D genome.

Development and validation of an exome-based SNP marker set for identification of the St, Jr and Jvs genomes of Thinopyrym intermedium in a wheat background

KEY MESSAGE

Cytogenetic analysis and array-based SNP genotyping of wheat- Th. intermedium introgression lines allowed identification of 634 chromosome-specific SNP markers across all twenty-one chromosomes of Th. intermedium (StJ r J vs , 2 n  = 6 x  = 42). Thinopyrum intermedium (2n = 6x = 42, StJrJvs) is one of the most promising reservoirs of useful genes including tolerance to abiotic stresses, perenniality and disease resistance not available in the cultivated bread wheat. The transfer of genetic diversity from wild species to wheat offers valuable responses to the effects of climate change. The new array-based single-nucleotide polymorphism (SNP) marker technology provides cheap and easy-to-use molecular markers for marker-assisted selection (MAS) in wheat breeding programmes. Here, we focus on the generation of a new chromosome-specific SNP marker set that can be used to characterize and identify the Th. intermedium chromosomes or chromosome segments transferred into wheat. A progressive investigation of marker development was conducted using 187 various newly developed wheat-Th. intermedium introgression lines and the Axiom® Wheat-Relative Genotyping array. We employed molecular cytogenetic techniques to clarify the genome constitution of the Th. intermedium parental lines and validated 634 chromosome-specific SNPs. Our data confirmed the allohexaploid nature of Th. intermedium and demonstrated that the St genome-specific GISH signal and markers are present at the centromeric regions of chromosomes 1Jvs, 2Jvs, 3Jvs and 7Jvs. The SNP markers presented here will be introduced into current wheat improvement programmes, offering a significant speed-up in wheat breeding and making it possible to deal with the transfer of the full genetic potential of Th. intermedium into wheat.

Identification and validation of a novel major QTL for all-stage stripe rust resistance on 1BL in the winter wheat line 20828

A major, likely novel stripe rust resistance QTL for all-stage resistance on chromosome arm 1BL identified in a 1.76-cM interval using a saturated linkage map was validated in four populations with different genetic backgrounds. Stripe rust is a globally important disease of wheat. Identification and utilization of new resistance genes are essential for breeding resistant cultivars. Wheat line 20828 has exhibited high levels of stripe rust resistance for over a decade. However, the genetics of stripe rust resistance in this line has not been studied. A set of 199 recombinant inbred lines (RILs) were developed from a cross between 20828 and a susceptible cultivar Chuannong 16. The RIL population was genotyped with the Wheat55K SNP (single nucleotide polymorphism) array and SSR (simple sequence repeat) markers and evaluated in four environments with current predominant Puccinia striiformis f. sp. tritici t races including CYR32, CYR33 and CYR34. Four stable QTL were located on chromosomes 1B (2 QTL), 4A and 6A. Among them, the major QTL, QYr.sicau-1B.1 (LOD = 23-28, PVE = 16-39%), was localized to a 1.76-cM interval flanked by SSR markers Xwmc216 and Xwmc156 on chromosome 1BL. Eight resistance genes were previously identified in the physical interval of QYr.sicau-1B.1. Compared with previous studies, QYr.sicau-1B.1 is a new gene for resistant to stripe rust. It was further verified by analysis of the closely linked SSR markers Xwmc216 and Xwmc156 in four other populations with different genetic backgrounds. QYr.sicau-1B.1 reduced the stripe rust disease index by up to 82.8%. Three minor stable QTL (located on chromosomes 1B, 4A and 6A, respectively) also added to the resistance level of QYr.sicau-1B.1. Our results provide valuable information for further fine mapping and cloning as well as molecular-assisted breeding with QYr.sicau-1B.1.

Fine mapping of the wheat powdery mildew resistance gene Pm52 using comparative genomics analysis and the Chinese Spring reference genomic sequence

A high-resolution genetic linkage map was constructed using the comparative genomics analysis approach and the wheat reference genome, which placed wheat powdery mildew resistance gene Pm52 in a 0.21-cM genetic interval on chromosome arm 2BL. The gene Pm52 confers resistance to powdery mildew and has been previously mapped on chromosome arm 2BL in winter wheat cultivar Liangxing 99. Because of its effectiveness against the disease, this study was initiated to finely map Pm52 using the comparative genomics analysis approach and the wheat reference genomic sequence. Based on the EST sequences that were located in the chromosome region flanking Pm52, four EST-SSR markers were developed, and another nine SSR markers were developed using the comparative genomics technology. These thirteen markers were integrated into a genetic linkage map using an F_2:3 subpopulation of the Liangxing 99 × Zhongzuo 9504 cross. Pm52 was mapped within a 3.2-cM genetic interval in the subpopulation that corresponded to a ~40-Mb genomic interval on chromosome arm 2BL of the Chinese Spring reference genome. The Pm52-flanking markers Xicsl163 and Xicsl62 identified 344 recombinant individuals from 8820 F₂ plants. Nine SSR markers generated from the Chinese Spring genomic interval were incorporated into a high-resolution genetic linkage map, which placed Pm52 in a 0.21-cM genetic interval corresponding to 5.6-Mb genomic region. The constructed high-resolution genetic linkage map will facilitate the map-based cloning of Pm52 and its marker-assisted selection.

A major QTL co-localized on chromosome 6BL and its epistatic interaction for enhanced wheat stripe rust resistance

Co-localization of a major QTL for wheat stripe rust resistance to a 3.9-cM interval on chromosome 6BL across both populations and another QTL on chromosome 2B with epistatic interaction. Cultivars with diverse resistance are the optimal strategy to minimize yield losses caused by wheat stripe rust (Puccinia striiformis f. sp. tritici). Two wheat populations involving resistant wheat lines P10078 and Snb"S" from CIMMYT were evaluated for stripe rust response in multiple environments. Pool analysis by Wheat660K SNP array showed that the overlapping interval on chromosome 6B likely harbored a major QTL between two populations. Then, linkage maps were constructed using KASP markers, and a co-localized locus with large effect on chromosome 6BL was detected using QTL analysis in both populations. The coincident QTL, named QYr.nwafu-6BL.2, explained 59.7% of the phenotypic maximum variation in the Mingxian 169 × P10078 and 52.5% in the Zhengmai 9023 × Snb"S" populations, respectively. This co-localization interval spanning 3.9 cM corresponds to ~ 30.5-Mb genomic region of the newest common wheat reference genome (IWGSC RefSeq v.1.0). In addition, another QTL was also detected on chromosome 2B in Zhengmai 9023 × Snb"S" population and it can accelerate expression of QYr.nwafu-6BL.2 to enhance resistance with epistatic interaction. Allowing for Pst response, marker genotypes, pedigree analysis and relative genetic distance, QYr.nwafu-6BL.2 is likely to be a distinct adult plant resistance QTL. Haplotype analysis of QYr.nwafu-6BL.2 revealed specific SNPs or alleles in the target region from a diversity panel of 176 unrelated wheat accessions. This QTL region provides opportunity for further map-based cloning, and haplotypes analysis enables pyramiding favorable alleles into commercial cultivars by marker-assisted selection.

Fine mapping of Aegilops peregrina co-segregating leaf and stripe rust resistance genes to distal-most end of 5DS

Novel rust resistance genes LrP and YrP from Ae. peregrina identified on chromosome 5D and the linked markers will aid deployment of these genes in combination with other major/minor genes. Aegilops peregrina, a wild tetraploid relative of wheat with genome constitution UUSS, displays genetic variation for resistance to leaf and stripe (yellow) rust. The wheat Ae. peregrina introgression line, IL pau16058, harbouring leaf and stripe rust resistance, was crossed with wheat cv. WL711 to generate an F_2:3 mapping population. Inheritance studies on this population indicated the transfer of dominant co-segregating resistance to leaf and stripe rust. Ethyl methane sulphonate mutagenesis of IL pau16058 identified independent loss-of-function mutants for leaf and stripe rust resistance, indicating that the leaf and stripe rust resistance is controlled by independent genes, herein designated LrP and YrP, respectively. A high-resolution genetic map of LrP and YrP was constructed using the Illumina Infinium iSelect 90K wheat array and resistance gene enrichment sequencing (RenSeq) markers. The map spans 4.19 cM on the distal-most region of the short arm of chromosome 5D, consisting of eight SNP markers and one microsatellite marker. LrP and YrP co-segregated with markers BS00163889 and 5DS44573_snp and was flanked distally by the SNP marker BS00129707 and proximally by 5DS149010, defining a 15.71 Mb region in the RefSeq v1.0 genome assembly.

Worldwide phylogeography and history of wheat genetic diversity

Since its domestication in the Fertile Crescent ~8000 to 10,000 years ago, wheat has undergone a complex history of spread, adaptation, and selection. To get better insights into the wheat phylogeography and genetic diversity, we describe allele distribution through time using a set of 4506 landraces and cultivars originating from 105 different countries genotyped with a high-density single-nucleotide polymorphism array. Although the genetic structure of landraces is collinear to ancient human migration roads, we observe a reshuffling through time, related to breeding programs, with the appearance of new alleles enriched with structural variations that may be the signature of introgressions from wild relatives after 1960.

Exome sequencing highlights the role of wild-relative introgression in shaping the adaptive landscape of the wheat genome

Introgression is a potential source of beneficial genetic diversity. The contribution of introgression to adaptive evolution and improvement of wheat as it was disseminated worldwide remains unknown. We used targeted re-sequencing of 890 diverse accessions of hexaploid and tetraploid wheat to identify wild-relative introgression. Introgression, and selection for improvement and environmental adaptation, each reduced deleterious allele burden. Introgression increased diversity genome wide and in regions harboring major agronomic genes, and contributed alleles explaining a substantial proportion of phenotypic variation. These results suggest that historic gene flow from wild relatives made a substantial contribution to the adaptive diversity of modern bread wheat.

Insights on the Proteases Involved in Barley and Wheat Grain Germination

Seed storage proteins must be hydrolyzed by proteases to deliver the amino acids essential for embryo growth and development. Several groups of proteases involved in this process have been identified in both the monocot and the dicot species. This review focuses on the implication of proteases during germination in two cereal species, barley and wheat, where proteolytic control during the germination process has considerable economic importance. Formerly, the participation of proteases during grain germination was inferred from reports of proteolytic activities, the expression of individual genes, or the presence of individual proteins and showed a prominent role for papain-like and legumain-like cysteine proteases and for serine carboxypeptidases. Nowadays, the development of new technologies and the release of the genomic sequences of wheat and barley have permitted the application of genome-scale approaches, such as those used in functional genomics and proteomics. Using these approaches, the repertoire of proteases known to be involved in germination has increased and includes members of distinct protease families. The development of novel techniques based on shotgun proteomics, activity-based protein profiling, and comparative and structural genomics will help to achieve a general view of the proteolytic process during germination.

Uncovering Genomic Regions Associated With 36 Agro-Morphological Traits in Indian Spring Wheat Using GWAS

Wheat genetic improvement by integration of advanced genomic technologies is one way of improving productivity. To facilitate the breeding of economically important traits in wheat, SNP loci and underlying candidate genes associated with the 36 agro-morphological traits were studied in a diverse panel of 404 genotypes. By using Breeders' 35K Axiom array in a comprehensive genome-wide association study covering 4364.79 cM of the wheat genome and applying a compressed mixed linear model, a total of 146 SNPs (-log10 P ≥ 4) were found associated with 23 traits out of 36 traits studied explaining 3.7-47.0% of phenotypic variance. To reveal this a subset of 260 genotypes was characterized phenotypically for six quantitative traits [days to heading (DTH), days to maturity (DTM), plant height (PH), spike length (SL), awn length (Awn_L), and leaf length (Leaf_L)] under five environments. Gene annotations mined ∼38 putative candidate genes which were confirmed using tissue and stage specific gene expression data from RNA Seq. We observed strong co-localized loci for four traits (glume pubescence, SL, PH, and awn color) on chromosome 1B (24.64 cM) annotated five putative candidate genes. This study led to the discovery of hitherto unreported loci for some less explored traits (such as leaf sheath wax, awn attitude, and glume pubescence) besides the refined chromosomal regions of known loci associated with the traits. This study provides valuable information of the genetic loci and their potential genes underlying the traits such as awn characters which are being considered as important contributors toward yield enhancement.

Genome-wide association study of resistance to stripe rust (Puccinia striiformis f. sp. tritici) in Sichuan wheat

BACKGROUND

Stripe rust (also called yellow rust) is a common and serious fungal disease of wheat (Triticum aestivum L.) caused by Puccinia striiformis f. sp. tritici. The narrow genetic basis of modern wheat cultivars and rapid evolution of the rust pathogen have been responsible for periodic and devastating epidemics of wheat rust diseases. In this study, we conducted a genome-wide association study with 44,059 single nucleotide polymorphism markers to identify loci associated with resistance to stripe rust in 244 Sichuan wheat accessions, including 79 landraces and 165 cultivars, in six environments.

RESULTS

In all the field assessments, 24 accessions displayed stable high resistance to stripe rust. Significant correlations among environments were observed for both infection (IT) and disease severity (DS), and high heritability levels were found for both IT and DS. Using mixed linear models, 12 quantitative trait loci (QTLs) significantly associated with IT and/or DS were identified. Two QTLs were mapped on chromosomes 5AS and 5AL and were distant from previously identified stripe rust resistance genes or QTL regions, indicating that they may be novel resistance loci.

CONCLUSIONS

Our results revealed that resistance alleles to stripe rust were accumulated in Sichuan wheat germplasm, implying direct or indirect selection for improved stripe rust resistance in elite wheat breeding programs. The identified stable QTLs or favorable alleles could be important chromosome regions in Sichuan wheat that controlled the resistance to stripe rust. These markers can be used molecular marker-assisted breeding of Sichuan wheat cultivars, and will be useful in the ongoing effort to develop new wheat cultivars with strong resistance to stripe rust.

Analysis of the recombination landscape of hexaploid bread wheat reveals genes controlling recombination and gene conversion frequency

BACKGROUND

Sequence exchange between homologous chromosomes through crossing over and gene conversion is highly conserved among eukaryotes, contributing to genome stability and genetic diversity. A lack of recombination limits breeding efforts in crops; therefore, increasing recombination rates can reduce linkage drag and generate new genetic combinations.

RESULTS

We use computational analysis of 13 recombinant inbred mapping populations to assess crossover and gene conversion frequency in the hexaploid genome of wheat (Triticum aestivum). We observe that high-frequency crossover sites are shared between populations and that closely related parents lead to populations with more similar crossover patterns. We demonstrate that gene conversion is more prevalent and covers more of the genome in wheat than in other plants, making it a critical process in the generation of new haplotypes, particularly in centromeric regions where crossovers are rare. We identify quantitative trait loci for altered gene conversion and crossover frequency and confirm functionality for a novel RecQ helicase gene that belongs to an ancient clade that is missing in some plant lineages including Arabidopsis.

CONCLUSIONS

This is the first gene to be demonstrated to be involved in gene conversion in wheat. Harnessing the RecQ helicase has the potential to break linkage drag utilizing widespread gene conversions.

Genome-wide identification and expression profiling of trihelix gene family under abiotic stresses in wheat

BACKGROUND

The trihelix gene family is a plant-specific transcription factor family that plays important roles in plant growth, development, and responses to abiotic stresses. However, to date, no systemic characterization of the trihelix genes has yet been conducted in wheat and its close relatives.

RESULTS

We identified a total of 94 trihelix genes in wheat, as well as 22 trihelix genes in Triticum urartu, 29 in Aegilops tauschii, and 31 in Brachypodium distachyon. We analyzed the chromosomal locations and orthology relations of the identified trihelix genes, and no trihelix gene was found to be located on chromosome 7A, 7B, or 7D of wheat, thereby reflecting the uneven distributions of wheat trihelix genes. Phylogenetic analysis indicated that the 186 identified trihelix proteins in wheat, rice, B. distachyon, and Arabidopsis were clustered into five major clades. The trihelix genes belonging to the same clades usually shared similar motif compositions and exon/intron structural patterns. Five pairs of tandem duplication genes and three pairs of segmental duplication genes were identified in the wheat trihelix gene family, thereby validating the supposition that more intrachromosomal gene duplication events occur in the genome of wheat than in that of other grass species. The tissue-specific expression and differential expression profiling of the identified genes under cold and drought stresses were analyzed by using RNA-seq data. qRT-PCR was also used to confirm the expression profiles of ten selected wheat trihelix genes under multiple abiotic stresses, and we found that these genes mainly responded to salt and cold stresses.

CONCLUSIONS

In this study, we identified trihelix genes in wheat and its close relatives and found that gene duplication events are the main driving force for trihelix gene evolution in wheat. Our expression profiling analysis demonstrated that wheat trihelix genes responded to multiple abiotic stresses, especially salt and cold stresses. The results of our study built a basis for further investigation of the functions of wheat trihelix genes and provided candidate genes for stress-resistant wheat breeding programs.

A Major Root Architecture QTL Responding to Water Limitation in Durum Wheat

The optimal root system architecture (RSA) of a crop is context dependent and critical for efficient resource capture in the soil. Narrow root growth angle promoting deeper root growth is often associated with improved access to water and nutrients in deep soils during terminal drought. RSA, therefore is a drought-adaptive trait that could minimize yield losses in regions with limited rainfall. Here, GWAS for seminal root angle (SRA) identified seven marker-trait associations clustered on chromosome 6A, representing a major quantitative trait locus (qSRA-6A) which also displayed high levels of pairwise LD (r 2 = 0.67). Subsequent haplotype analysis revealed significant differences between major groups. Candidate gene analysis revealed loci related to gravitropism, polar growth and hormonal signaling. No differences were observed for root biomass between lines carrying hap1 and hap2 for qSRA-6A, highlighting the opportunity to perform marker-assisted selection for the qSRA-6A locus and directly select for wide or narrow RSA, without influencing root biomass. Our study revealed that the genetic predisposition for deep rooting was best expressed under water-limitation, yet the root system displayed plasticity producing root growth in response to water availability in upper soil layers. We discuss the potential to deploy root architectural traits in cultivars to enhance yield stability in environments that experience limited rainfall.

Genomic tools for durum wheat breeding: de novo assembly of Svevo transcriptome and SNP discovery in elite germplasm

BACKGROUND

The tetraploid durum wheat (Triticum turgidum L. ssp. durum Desf. Husnot) is an important crop which provides the raw material for pasta production and a valuable source of genetic diversity for breeding hexaploid wheat (Triticum aestivum L.). Future breeding efforts to enhance yield potential and climate resilience will increasingly rely on genomics-based approaches to identify and select beneficial alleles. A deeper characterisation of the molecular and functional diversity of the durum wheat transcriptome will be instrumental to more effectively harness its genetic diversity.

RESULTS

We report on the de novo transcriptome assembly of durum wheat cultivar 'Svevo'. The transcriptome of four tissues/organs (shoots and roots at the seedling stage, reproductive organs and developing grains) was assembled de novo, yielding 180,108 contigs, with a N50 length of 1121 bp and mean contig length of 883 bp. Alignment against the transcriptome of nine plant species identified 43% of transcripts with homology to at least one reference transcriptome. The functional annotation was completed by means of a combination of complementary software. The presence of differential expression between the A- and B-homoeolog copies of the durum wheat tetraploid genome was ascertained by phase reconstruction of polymorphic sites based on the T. urartu transcripts and inferring homoeolog-specific sequences. We observed greater expression divergence between A and B homoeologs in grains rather than in leaves and roots. The transcriptomes of 13 durum wheat cultivars spanning the breeding period from 1969 to 2005 were analysed for SNP diversity, leading to 95,358 non-rare, hemi-SNPs shared among two or more cultivars and 33,747 locus-specific (diploid inheritance) SNPs.

CONCLUSIONS

Our study updates and expands the de novo transcriptome reference assembly available for durum wheat. Out of 180,108 assembled transcripts, 13,636 were specific to the Svevo cultivar as compared to the only other reference transcriptome available for durum, thus contributing to the identification of the tetraploid wheat pan-transcriptome. Additionally, the analysis of 13 historically relevant hallmark varieties produced a SNP dataset that could successfully validate the genotyping in tetraploid wheat and provide a valuable resource for genomics-assisted breeding of both tetraploid and hexaploid wheats.

Durum wheat genome highlights past domestication signatures and future improvement targets

The domestication of wild emmer wheat led to the selection of modern durum wheat, grown mainly for pasta production. We describe the 10.45 gigabase (Gb) assembly of the genome of durum wheat cultivar Svevo. The assembly enabled genome-wide genetic diversity analyses revealing the changes imposed by thousands of years of empirical selection and breeding. Regions exhibiting strong signatures of genetic divergence associated with domestication and breeding were widespread in the genome with several major diversity losses in the pericentromeric regions. A locus on chromosome 5B carries a gene encoding a metal transporter (TdHMA3-B1) with a non-functional variant causing high accumulation of cadmium in grain. The high-cadmium allele, widespread among durum cultivars but undetected in wild emmer accessions, increased in frequency from domesticated emmer to modern durum wheat. The rapid cloning of TdHMA3-B1 rescues a wild beneficial allele and demonstrates the practical use of the Svevo genome for wheat improvement.

Assessing yield gap in high productive countries by designing wheat ideotypes

Designing crop ideotypes in silico is a powerful tool to explore the crop yield potential and yield gap. We defined yield gap as the difference between yield potential of a crop ideotype optimized under local environment and yield of an existing cultivar under optimal management. Wheat ideotypes were designed for the current climate using the Sirius model for both water-limited and irrigated conditions in two high wheat-productive countries viz. the United Kingdom (UK) and New Zealand (NZ) with the objective of estimating yield gap. The mean ideotype yields of 15.0-19.0 t ha-1 were achieved in water-limited conditions in the UK and NZ, whereas 15.6-19.5 t ha-1 under irrigated conditions. Substantial yield gaps were found in both water-limited, 28-31% (4-6 t ha-1), and irrigated conditions, 30-32% (5-6 t ha-1) in the UK and NZ. Both yield potential (25-27%) and yield gap (32-38%) were greater in NZ than the UK. Ideotype design is generic and could apply globally for estimating yield gap. Despite wheat breeding efforts, the considerable yield gap still potentially exists in high productive countries such as the UK and NZ. To accelerate breeding, wheat ideotypes can provide the key traits for wheat improvement and closing the yield gap.

CRISPR-based genome editing in wheat: a comprehensive review and future prospects

CRISPR technology has vividly increased its applications in last five years for genome editing in a wide range of organisms from bacteria to plants. It is mostly applied in the field of mammalian research. This emerging versatile tool can be utilized in crop improvement by targeting various traits to increase economic value and adaptability of the crop species under changing climate. In plants, Arabidopsis and rice are the most studied plant species in genome editing through CRISPR technology. Wheat is lagging behind in the utilization of CRISPR based genome modifications. The hexaploid, large genome size and the recalcitrant nature in terms of tissue culture are the major obstacles for CRISPR utilization in wheat. Recently, the IWGSC released the high quality of reference genome for wheat which will greatly accelerate the application of CRISPR-based genome engineering in wheat and helps to resolve the global issue of food security in coming decades. The exogenous DNA-free improved mutants with CRISPR technology having desired traits will increase the productivity under biotic and abiotic stress conditions. To address complex traits involving multigene, recently developed multiplex genome editing toolkits can be used. This is a first review of its kind in which the practical utilization and updates on CRISPR validation in wheat along with its future prospects for use of this technology in wheat improvement are comprehensively discussed. Thus, the compiled information will immensely benefit the researchers for utilization of CRISPR system in wheat improvement across the globe.

Integrating genomic resources to present full gene and putative promoter capture probe sets for bread wheat

BACKGROUND

Whole-genome shotgun resequencing of wheat is expensive because of its large, repetitive genome. Moreover, sequence data can fail to map uniquely to the reference genome, making it difficult to unambiguously assign variation. Resequencing using target capture enables sequencing of large numbers of individuals at high coverage to reliably identify variants associated with important agronomic traits. Previous studies have implemented complementary DNA/exon or gene-based probe sets in which the promoter and intron sequence is largely missing alongside newly characterized genes from the recent improved reference sequences.

RESULTS

We present and validate 2 gold standard capture probe sets for hexaploid bread wheat, a gene and a putative promoter capture, which are designed using recently developed genome sequence and annotation resources. The captures can be combined or used independently. We demonstrate that the capture probe sets effectively enrich the high-confidence genes and putative promoter regions that were identified in the genome alongside a large proportion of the low-confidence genes and associated promoters. Finally, we demonstrate successful sample multiplexing that allows generation of adequate sequence coverage for single-nucleotide polymorphism calling while significantly reducing cost per sample for gene and putative promoter capture.

CONCLUSIONS

We show that a capture design employing an "island strategy" can enable analysis of the large gene/putative promoter space of wheat with only 2 × 160 Mbp probe sets. Furthermore, these assays extend the regions of the wheat genome that are amenable to analyses beyond its exome, providing tools for detailed characterization of these regulatory regions in large populations.

CropSNPdb: a database of SNP array data for Brassica crops and hexaploid bread wheat

Advances in sequencing technology have led to a rapid rise in the genomic data available for plants, driving new insights into the evolution, domestication and improvement of crops. Single nucleotide polymorphisms (SNPs) are a major component of crop genomic diversity, and are invaluable as genetic markers in research and breeding programs. High-throughput SNP arrays, or 'SNP chips', can generate reproducible sets of informative SNP markers and have been broadly adopted. Although there are many public repositories for sequencing data, which are routinely uploaded, there are no formal repositories for crop SNP array data. To make SNP array data more easily accessible, we have developed CropSNPdb (http://snpdb.appliedbioinformatics.com.au), a database for SNP array data produced by the Illumina Infinium™ hexaploid bread wheat (Triticum aestivum) 90K and Brassica 60K arrays. We currently host SNPs from datasets covering 526 Brassica lines and 309 bread wheat lines, and provide search, download and upload utilities for users. CropSNPdb provides a useful repository for these data, which can be applied for a range of genomics and molecular crop-breeding activities.

Establishment of an In Vitro Fertilization System in Wheat (Triticum aestivum L.)

In vitro fertilization (IVF) systems using isolated gametes have been utilized to dissect post-fertilization events in angiosperms, since the female gametophytes of plants are deeply embedded within ovaries. In addition, IVF systems have been used to produce hybrid and polyploid zygotes. Complete IVF systems have been established in maize and rice, two of three major crop species providing the majority of human energy intake. Among those crop species, gametes of wheat have not been used to establish a complete IVF system successfully. In this study, a wheat IVF system was developed to introduce the advantages of this technology to wheat research. Fusion of gametes was performed via a modified electrofusion method, and the fusion product, a zygote, formed a cell wall and two nucleoli. The first division of zygotes was observed 19-27 h after fusion, and the resulting two-celled embryo developed into 10-20-celled globular-like embryos and multicellular club-shaped embryos by 3 and 7-10 d after fusion, respectively. Such zygotic division profiles were mostly consistent with those in the embryo sac, suggesting that the division profile of IVF-produced early embryos reflects that of early embryos in planta. Although the IVF-produced club-shaped embryos did not develop into differentiated embryos but into compact embryonic calli, fertile plants could be regenerated from the embryonic calli, and the seeds harvested from those plants grew normally into seedlings. The IVF system described here might become an important technique for generating new genotypes of wheat and/or new hybrids as well as for investigating fertilization-induced events in wheat.

The role of transposable elements and DNA damage repair mechanisms in gene duplications and gene fusions in plant genomes

Plant genomes are shaped by structural variation. Gene-size insertions and among most prominent events and can have significant effects on amplification of gene families as well as facilitate new gene fusions. Transposable elements as well as plant DNA repair machinery have overlapping contributions to these events, and often work in synergy. Activity of transposable elements is often lineage specific and can preferentially affect specific gene families, such as disease resistance genes. Once duplicated, genes themselves can serve templates for additional variation that can arise from non-allelic homologous recombination. Non-homologous DNA repair mechanisms contribute to additional variation and diversify the mechanisms of gene movement, such as through ligation of extra-chromosomal DNA fragments. Genomic processes that generate structural variation can be induced by stress and, therefore, can provide adaptive potential. This review describes mechanisms that contribute to gene-size structural variation in plants, result in gene duplication and generation of new plant genes through gene fusion.

Divergence between bread wheat and Triticum militinae in the powdery mildew resistance QPm.tut-4A locus and its implications for cloning of the resistance gene

A segment of Triticum militinae chromosome 7G harbors a gene(s) conferring powdery mildew resistance which is effective at both the seedling and the adult plant stages when transferred into bread wheat (T. aestivum). The introgressed segment replaces a piece of wheat chromosome arm 4AL. An analysis of segregating materials generated to positionally clone the gene highlighted that in a plant heterozygous for the introgression segment, only limited recombination occurs between the introgressed region and bread wheat 4A. Nevertheless, 75 genetic markers were successfully placed within the region, thereby confining the gene to a 0.012 cM window along the 4AL arm. In a background lacking the Ph1 locus, the localized rate of recombination was raised 33-fold, enabling the reduction in the length of the region containing the resistance gene to a 480 kbp stretch harboring 12 predicted genes. The substituted segment in the reference sequence of bread wheat cv. Chinese Spring is longer (640 kbp) and harbors 16 genes. A comparison of the segments' sequences revealed a high degree of divergence with respect to both their gene content and nucleotide sequence. Of the 12 T. militinae genes, only four have a homolog in cv. Chinese Spring. Possible candidate genes for the resistance have been identified based on function predicted from their sequence.

Integrating genomic resources to present full gene and putative promoter capture probe sets for bread wheat

BACKGROUND

Whole-genome shotgun resequencing of wheat is expensive because of its large, repetitive genome. Moreover, sequence data can fail to map uniquely to the reference genome, making it difficult to unambiguously assign variation. Resequencing using target capture enables sequencing of large numbers of individuals at high coverage to reliably identify variants associated with important agronomic traits. Previous studies have implemented complementary DNA/exon or gene-based probe sets in which the promoter and intron sequence is largely missing alongside newly characterized genes from the recent improved reference sequences.

RESULTS

We present and validate 2 gold standard capture probe sets for hexaploid bread wheat, a gene and a putative promoter capture, which are designed using recently developed genome sequence and annotation resources. The captures can be combined or used independently. We demonstrate that the capture probe sets effectively enrich the high-confidence genes and putative promoter regions that were identified in the genome alongside a large proportion of the low-confidence genes and associated promoters. Finally, we demonstrate successful sample multiplexing that allows generation of adequate sequence coverage for single-nucleotide polymorphism calling while significantly reducing cost per sample for gene and putative promoter capture.

CONCLUSIONS

We show that a capture design employing an "island strategy" can enable analysis of the large gene/putative promoter space of wheat with only 2 × 160 Mbp probe sets. Furthermore, these assays extend the regions of the wheat genome that are amenable to analyses beyond its exome, providing tools for detailed characterization of these regulatory regions in large populations.

Plant genome sequences: past, present, future

The green plants (Viridiplantae) are an essential kingdom of life, responsible via photosynthesis for the majority of global primary production, and directly utilized by humankind for nutrition, animal feed, fuel, clothing, medicine and other purposes. There are an estimated 391 000 species of land plants, in addition to 8000 species of green algae. Their genomes are unusually diverse compared to those of other kingdoms, ranging in size from ∼10 Mb to over 100 Gb. Knowledge of plant genomes initially lagged behind those of other kingdoms but has greatly increased with the development of new technologies for DNA sequencing; bioinformatic analysis, rather than data production, is increasingly the bottleneck to further knowledge. Recent proposals are now contemplating the sequencing, assembly and annotation of the genomes of all of the world's plant species; meanwhile, low coverage sequencing to measure diversity across collections and wild populations has already become commonplace for many species, especially those utilized as crops.

Assessing the genetic diversity and characterizing genomic regions conferring Tan Spot resistance in cultivated rye

Rye (Secale cereale L.) is known for its wide adaptation due to its ability to tolerate harsh environments in semiarid areas. To assess the diversity in rye we genotyped a panel of 178 geographically diverse accessions of four Secale sp. from U.S. National Small Grains Collection using 4,037 high-quality SNPs (single nucleotide polymorphisms) developed by genotyping-by-sequencing (GBS). PCA and STRUCTURE analysis revealed three major clusters that separate S. cereale L. from S. strictum and S. sylvestre, however, genetic clusters did not correlate with geographic origins and growth habit (spring/winter). The panel was evaluated for response to Pyrenophora tritici-repentis race 5 (PTR race 5) and nearly 59% accessions showed resistance or moderate resistance. Genome-wide association study (GWAS) was performed on S. cereale subsp. cereale using the 4,037 high-quality SNPs. Two QTLs (QTs.sdsu-5R and QTs.sdsu-2R) on chromosomes 5R and 2R were identified conferring resistance to PTR race 5 (p < 0.001) that explained 13.1% and 11.6% of the phenotypic variation, respectively. Comparative analysis showed a high degree of synteny between rye and wheat with known rearrangements as expected. QTs.sdsu-2R was mapped in the genomic region corresponding to wheat chromosome group 2 and QTs.sdsu-5R was mapped to a small terminal region on chromosome 4BL. Based on the genetic diversity, a set of 32 accessions was identified to represents more than 99% of the allelic diversity with polymorphic information content (PIC) of 0.25. This set can be utilized for genetic characterization of useful traits and genetic improvement of rye, triticale, and wheat.

Mal de Río Cuarto virus infection causes hormone imbalance and sugar accumulation in wheat leaves

BACKGROUND

Mal de Río Cuarto virus (MRCV) infects several monocotyledonous species including maize and wheat. Infected plants show shortened internodes, partial sterility, increased tillering and reduced root length. To better understand the molecular basis of the plant-virus interactions leading to these symptoms, we combined RNA sequencing with metabolite and hormone measurements.

RESULTS

More than 3000 differentially accumulated transcripts (DATs) were detected in MRCV-infected wheat plants at 21 days post inoculation compared to mock-inoculated plants. Infected plants exhibited decreased levels of TaSWEET13 transcripts, which are involved in sucrose phloem loading. Soluble sugars, starch, trehalose 6-phosphate (Tre6P), and organic and amino acids were all higher in MRCV-infected plants. In addition, several transcripts related to plant hormone metabolism, transport and signalling were increased upon MRCV infection. Transcripts coding for GA20ox, D14, MAX2 and SMAX1-like proteins involved in gibberellin biosynthesis and strigolactone signalling, were reduced. Transcripts involved in jasmonic acid, ethylene and brassinosteroid biosynthesis, perception and signalling and in auxin transport were also altered. Hormone measurements showed that jasmonic acid, brassinosteroids, abscisic acid and indole-3-acetic acid were significantly higher in infected leaves.

CONCLUSIONS

Our results indicate that MRCV causes a profound hormonal imbalance that, together with alterations in sugar partitioning, could account for the symptoms observed in MRCV-infected plants.

Viruses and Evolution - Viruses First? A Personal Perspective

The discovery of exoplanets within putative habitable zones revolutionized astrobiology in recent years. It stimulated interest in the question about the origin of life and its evolution. Here, we discuss what the roles of viruses might have been at the beginning of life and during evolution. Viruses are the most abundant biological entities on Earth. They are present everywhere, in our surrounding, the oceans, the soil and in every living being. Retroviruses contributed to about half of our genomic sequences and to the evolution of the mammalian placenta. Contemporary viruses reflect evolution ranging from the RNA world to the DNA-protein world. How far back can we trace their contribution? Earliest replicating and evolving entities are the ribozymes or viroids fulfilling several criteria of life. RNA can perform many aspects of life and influences our gene expression until today. The simplest structures with non-protein-coding information may represent models of life built on structural, not genetic information. Viruses today are obligatory parasites depending on host cells. Examples of how an independent lifestyle might have been lost include mitochondria, chloroplasts, Rickettsia and others, which used to be autonomous bacteria and became intracellular parasites or endosymbionts, thereby losing most of their genes. Even in vitro the loss of genes can be recapitulated all the way from coding to non-coding RNA. Furthermore, the giant viruses may indicate that there is no sharp border between living and non-living entities but an evolutionary continuum. Here, it is discussed how viruses can lose and gain genes, and that they are essential drivers of evolution. This discussion may stimulate the thinking about viruses as early possible forms of life. Apart from our view “viruses first”, there are others such as “proteins first” and “metabolism first.”

Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton

Allotetraploid cotton is an economically important natural-fiber-producing crop worldwide. After polyploidization, Gossypium hirsutum L. evolved to produce a higher fiber yield and to better survive harsh environments than Gossypium barbadense, which produces superior-quality fibers. The global genetic and molecular bases for these interspecies divergences were unknown. Here we report high-quality de novo-assembled genomes for these two cultivated allotetraploid species with pronounced improvement in repetitive-DNA-enriched centromeric regions. Whole-genome comparative analyses revealed that species-specific alterations in gene expression, structural variations and expanded gene families were responsible for speciation and the evolutionary history of these species. These findings help to elucidate the evolution of cotton genomes and their domestication history. The information generated not only should enable breeders to improve fiber quality and resilience to ever-changing environmental conditions but also can be translated to other crops for better understanding of their domestication history and use in improvement.

Novel Insights into Plant Genome Evolution and Adaptation as Revealed through Transposable Elements and Non-Coding RNAs in Conifers

Plant genomes are punctuated by repeated bouts of proliferation of transposable elements (TEs), and these mobile bursts are followed by silencing and decay of most of the newly inserted elements. As such, plant genomes reflect TE-related genome expansion and shrinkage. In general, these genome activities involve two mechanisms: small RNA-mediated epigenetic repression and long-term mutational decay and deletion, that is, genome-purging. Furthermore, the spatial relationships between TE insertions and genes are an important force in shaping gene regulatory networks, their downstream metabolic and physiological outputs, and thus their phenotypes. Such cascading regulations finally set up a fitness differential among individuals. This brief review demonstrates factual evidence that unifies most updated conceptual frameworks covering genome size, architecture, epigenetic reprogramming, and gene expression. It aims to give an overview of the impact that TEs may have on genome and adaptive evolution and to provide novel insights into addressing possible causes and consequences of intimidating genome sizes (20⁻30 Gb) in a taxonomic group, conifers.

Potential of Aegilops sp. for Improvement of Grain Processing and Nutritional Quality in Wheat (Triticum aestivum)

Wheat is one of the most important staple crops in the world and good source of calories and nutrition. Its flour and dough have unique physical properties and can be processed to make unique products like bread, cakes, biscuits, pasta, noodles etc., which is not possible from other staple crops. Due to domestication, the genetic variability of the genes coding for different economically important traits in wheat is narrow. This genetic variability can be increased by utilizing its wild relatives. Its closest relative, genus Aegilops can be an important source of new alleles. Aegilops has played a very important role in evolution of tetraploid and hexaploid wheat. It consists of 22 species with C, D, M, N, S, T and U genomes with high allelic diversity relative to wheat. Its utilization for wheat improvement for various abiotic and biotic stresses has been reported by various scientific publications. Here in, for the first time, we review the potential of Aegilops for improvement of processing and nutritional traits in wheat. Among processing quality related gluten proteins; high molecular weight glutenins (HMW GS), being easiest to study have been explored in highest number of accessions or lines i.e., 681 belonging to 13 species and selected ones like Ae. searsii, Ae. geniculata and Ae. longissima have been linked with improved bread making quality of wheat. Gliadins and low molecular weight glutenins (LMW GS) have also been extensively explored for wheat improvement and Ae. umbellulata specific LMW GS have been linked with wheat bread making quality improvement. Aegilops has been explored for seed texture diversity and proteins like puroindolins (Pin) and grain softness proteins (GSP). For nutrition quality improvement, it has been screened for essential micronutrients like Fe, Zn, phytochemicals like carotenoids and dietary fibers like arabinoxylan and β-glucan. Ae. kotschyi and Ae. biuncialis transfer in wheat have been associated with higher Fe, Zn content. In this article we have tried to compile information available on exploration of nutritional and processing quality related traits in Aegilops section and their utilization for wheat improvement by different approaches.

Genetic Modification for Wheat Improvement: From Transgenesis to Genome Editing

To feed the growing human population, global wheat yields should increase to approximately 5 tonnes per ha from the current 3.3 tonnes by 2050. To reach this goal, existing breeding practices must be complemented with new techniques built upon recent gains from wheat genome sequencing, and the accumulated knowledge of genetic determinants underlying the agricultural traits responsible for crop yield and quality. In this review we primarily focus on the tools and techniques available for accessing gene functions which lead to clear phenotypes in wheat. We provide a view of the development of wheat transformation techniques from a historical perspective, and summarize how techniques have been adapted to obtain gain-of-function phenotypes by gene overexpression, loss-of-function phenotypes by expressing antisense RNAs (RNA interference or RNAi), and most recently the manipulation of gene structure and expression using site-specific nucleases, such as CRISPR/Cas9, for genome editing. The review summarizes recent successes in the application of wheat genetic manipulation to increase yield, improve nutritional and health-promoting qualities in wheat, and enhance the crop's resistance to various biotic and abiotic stresses.

A Low Resolution Epistasis Mapping Approach To Identify Chromosome Arm Interactions in Allohexaploid Wheat

Epistasis is an important contributor to genetic variance. In inbred populations, pairwise epistasis is present as additive by additive interactions. Testing for epistasis presents a multiple testing problem as the pairwise search space for modest numbers of markers is large. Single markers do not necessarily track functional units of interacting chromatin as well as haplotype based methods do. To harness the power of multiple markers while minimizing the number of tests conducted, we present a low resolution test for epistatic interactions across whole chromosome arms. Epistasis covariance matrices were constructed from the additive covariances of individual chromosome arms. These covariances were subsequently used to estimate an epistatic variance parameter while correcting for background additive and epistatic effects. We find significant epistasis for 2% of the interactions tested for four agronomic traits in a winter wheat breeding population. Interactions across homeologous chromosome arms were identified, but were less abundant than other chromosome arm pair interactions. The homeologous chromosome arm pair 4BL/4DL showed a strong negative relationship between additive and interaction effects that may be indicative of functional redundancy. Several chromosome arms appeared to act as hubs in an interaction network, suggesting that they may contain important regulatory factors. The differential patterns of epistasis across different traits demonstrate that detection of epistatic interactions is robust when correcting for background additive and epistatic effects in the population. The low resolution epistasis mapping method presented here identifies important epistatic interactions with a limited number of statistical tests at the cost of low precision.

Prediction of Subgenome Additive and Interaction Effects in Allohexaploid Wheat

Whole genome duplications have played an important role in the evolution of angiosperms. These events often occur through hybridization between closely related species, resulting in an allopolyploid with multiple subgenomes. With the availability of affordable genotyping and a reference genome to locate markers, breeders of allopolyploids now have the opportunity to manipulate subgenomes independently. This also presents a unique opportunity to investigate epistatic interactions between homeologous orthologs across subgenomes. We present a statistical framework for partitioning genetic variance to the subgenomes of an allopolyploid, predicting breeding values for each subgenome, and determining the importance of inter-genomic epistasis. We demonstrate using an allohexaploid wheat breeding population evaluated in Ithaca, NY and an important wheat dataset from CIMMYT previously shown to demonstrate non-additive genetic variance. Subgenome covariance matrices were constructed and used to calculate subgenome interaction covariance matrices for variance component estimation and genomic prediction. We propose a method to extract population structure from all subgenomes at once before covariances are calculated to reduce collinearity between subgenome estimates. Variance parameter estimation was shown to be reliable for additive subgenome effects, but was less reliable for subgenome interaction components. Predictive ability was equivalent to current genomic prediction methods. Including only inter-genomic interactions resulted in the same increase in accuracy as modeling all pairwise marker interactions. Thus, we provide a new tool for breeders of allopolyploid crops to characterize the genetic architecture of existing populations, determine breeding goals, and develop new strategies for selection of additive effects and fixation of inter-genomic epistasis.

Genome-wide analysis, expansion and expression of the NAC family under drought and heat stresses in bread wheat (T. aestivum L.)

The NAC family is one of the largest plant-specific transcription factor families, and some of its members are known to play major roles in plant development and response to biotic and abiotic stresses. Here, we inventoried 488 NAC members in bread wheat (Triticum aestivum). Using the recent release of the wheat genome (IWGS RefSeq v1.0), we studied duplication events focusing on genomic regions from 4B-4D-5A chromosomes as an example of the family expansion and neofunctionalization of TaNAC members. Differentially expressed TaNAC genes in organs and in response to abiotic stresses were identified using publicly available RNAseq data. Expression profiling of 23 selected candidate TaNAC genes was studied in leaf and grain from two bread wheat genotypes at two developmental stages in field drought conditions and revealed insights into their specific and/or overlapping expression patterns. This study showed that, of the 23 TaNAC genes, seven have a leaf-specific expression and five have a grain-specific expression. In addition, the grain-specific genes profiles in response to drought depend on the genotype. These genes may be considered as potential candidates for further functional validation and could present an interest for crop improvement programs in response to climate change. Globally, the present study provides new insights into evolution, divergence and functional analysis of NAC gene family in bread wheat.

Genome-wide Association Mapping and Prediction of Adult Stage Septoria tritici Blotch Infection in European Winter Wheat via High-Density Marker Arrays

blotch (STB) caused by the fungus is a devastating foliar disease of wheat ( L.) that can lead to substantial yield losses. Quantitative genetic resistance has been proposed as a durable strategy for STB control. In this study, we dissected the genetic basis of STB infection in 371 European wheat varieties based on 35k and 90k single nucleotide polymorphism marker arrays. The phenotypic data analyses suggested that large genetic variance exists for STB infection with a broad-sense heritability of 0.78. Genome-wide association studies (GWAS) propose the highly quantitative nature of STB infection with potential associations on chromosomes 1A, 1B, 2D, 4A, 5A, 6A, 6D, 7A, and 7B. Increased marker density in GWAS by combining markers from both arrays helped to detect additional markers explaining increased genotypic variance. Linkage disequilibrium analyses revealed genes with a possible role in disease resistance. The potential of genomic prediction (GP) assessed via two models accounting for additive effects and additive plus epistatic interactions among the loci suggested the possibility of genomic selection for improved STB resistance. Genomic prediction results also indicated that the higher-order epistatic interactions are not abundant and that both marker platforms are equally suitable for GP of STB infection. Our results provide further understanding of the quantitative genetic nature of STB infection, serve as a resource for marker-assisted breeding, and highlight the potential of genomic selection for improved STB resistance.