Wheat functional genomics in the era of next generation sequencing: An update

Genome-wide association mapping in elite winter wheat breeding for yield improvement

Increased grain yield (GY) is the primary breeding target of wheat breeders. We performed the genome-wide association study (GWAS) on 168 elite winter wheat lines from an ongoing breeding program to identify the main determinants of grain yield. Sequencing of Diversity Array Technology fragments (DArTseq) resulted in 19,350 single-nucleotide polymorphism (SNP) and presence-absence variation (PAV) markers. We identified 15 main genomic regions located in ten wheat chromosomes (1B, 2B, 2D, 3A, 3D, 5A, 5B, 6A, 6B, and 7B) that explained from 7.9 to 20.3% of the variation in grain yield and 13.3% of the yield stability. Loci identified in the reduced genepool are important for wheat improvement using marker-assisted selection. We found marker-trait associations between three genes involved in starch biosynthesis and grain yield. Two starch synthase genes (TraesCS2B03G1238800 and TraesCS2D03G1048800) and a sucrose synthase gene (TraesCS3D03G0024300) were found in regions of QGy.rut-2B.2, QGy.rut-2D.1, and QGy.rut-3D, respectively. These loci and other significantly associated SNP markers found in this study can be used for pyramiding favorable alleles in high-yielding varieties or to improve the accuracy of prediction in genomic selection.

Plant Tolerance to Drought Stress with Emphasis on Wheat

Environmental stresses, such as drought, have negative effects on crop yield. Drought is a stress whose impact tends to increase in some critical regions. However, the worldwide population is continuously increasing and climate change may affect its food supply in the upcoming years. Therefore, there is an ongoing effort to understand the molecular processes that may contribute to improving drought tolerance of strategic crops. These investigations should contribute to delivering drought-tolerant cultivars by selective breeding. For this reason, it is worthwhile to review regularly the literature concerning the molecular mechanisms and technologies that could facilitate gene pyramiding for drought tolerance. This review summarizes achievements obtained using QTL mapping, genomics, synteny, epigenetics, and transgenics for the selective breeding of drought-tolerant wheat cultivars. Synthetic apomixis combined with the msh1 mutation opens the way to induce and stabilize epigenomes in crops, which offers the potential of accelerating selective breeding for drought tolerance in arid and semi-arid regions.

Genome resequencing reveals independent domestication and breeding improvement of naked oat

As an important cereal crop, common oat, has attracted more and more attention due to its healthy nutritional components and bioactive compounds. Here, high-depth resequencing of 115 oat accessions and closely related hexaploid species worldwide was performed. Based on genetic diversity and linkage disequilibrium analysis, it was found that hulled oat (Avena sativa) experienced a more severe bottleneck than naked oat (Avena sativa var. nuda). Combined with the divergence time of ∼51,200 years ago, the previous speculation that naked oat was a variant of hulled oat was rejected. It was found that the common segments that hulled oat introgressed to naked oat cultivars contained 444 genes, mainly enriched in photosynthetic efficiency-related pathways. Selective sweeps during environmental adaptation and breeding improvement were identified in the naked oat genome. Candidate genes associated with smut resistance and the days to maturity phenotype were also identified. Our study provides genomic resources and new insights into naked oat domestication and breeding.

Conditional QTL mapping for seed germination and seedling traits under salt stress and candidate gene prediction in wheat

Breeding new wheat varieties with salt resistance is one of the best ways to solve a constraint on the sustainability and expansion of wheat cultivation. Therefore, understanding the molecular components or genes related to salt tolerance must contribute to the cultivation of salt-tolerant varieties. The present study used a recombinant inbred line (RIL) population to genetically dissect the effects of different salt stress concentrations on wheat seed germination and seedling traits using two quantitative trait locus (QTL) mapping methods. A total of 31 unconditional and 11 conditional QTLs for salt tolerance were identified on 11 chromosomes explaining phenotypic variation (PVE) ranging from 2.01 to 65.76%. Of these, 15 major QTLs were found accounting for more than 10% PVE. QTL clusters were detected on chromosomes 2A and 3B in the marker intervals 'wPt-8328 and wPt-2087' and 'wPt-666008 and wPt-3620', respectively, involving more than one salt tolerance trait. QRdw3B and QSfw3B.2 were most consistent in two or more salt stress treatments. 16 candidate genes associated with salt tolerance were predicted in wheat. These results could be useful to improve salt tolerance by marker-assisted selection (MAS) and shed new light on understanding the genetic basis of salt tolerance in wheat.

Insights into breeding history, hotspot regions of selection, and untapped allelic diversity for bread wheat breeding

Breeding has increasingly altered the genetics of crop plants since the domestication of their wild progenitors. It is postulated that the genetic diversity of elite wheat breeding pools is too narrow to cope with future challenges. In contrast, plant genetic resources (PGRs) of wheat stored in genebanks are valuable sources of unexploited genetic diversity. Therefore, to ensure breeding progress in the future, it is of prime importance to identify the useful allelic diversity available in PGRs and to transfer it into elite breeding pools. Here, a diverse collection consisting of modern winter wheat cultivars and genebank accessions was investigated based on reduced-representation genomic sequencing and an iSelect single nucleotide polymorphism (SNP) chip array. Analyses of these datasets provided detailed insights into population structure, levels of genetic diversity, sources of new allelic diversity, and genomic regions affected by breeding activities. We identified 57 regions representing genomic signatures of selection and 827 regions representing private alleles associated exclusively with genebank accessions. The presence of known functional wheat genes, quantitative trait loci, and large chromosomal modifications, i.e., introgressions from wheat wild relatives, provided initial evidence for putative traits associated within these identified regions. These findings were supported by the results of ontology enrichment analyses. The results reported here will stimulate further research and promote breeding in the future by allowing for the targeted introduction of novel allelic diversity into elite wheat breeding pools.

Genomic Regions, Molecular Markers, and Flanking Genes of Metribuzin Tolerance in Wheat (Triticum aestivum L.)

Understanding the genetics of metribuzin (a group C herbicide) tolerance in wheat is vital in developing tolerant cultivars to improve wheat productivity in dryland farming systems. This study investigated metribuzin tolerance in wheat by conducting a Genome-wide Association Studies (GWAS) with a panel of 150 wheat genotypes of diverse genetic backgrounds and genotyped them with the wheat 90 K SNP genotyping assay. The phenotyping was conducted in a temperature-controlled glasshouse at the University of Western Australia (UWA). Genotypes were sprayed with a metribuzin dose of 400 grams of active ingredient (g. a.i.) ha^-1 as pre-emergent in a specialized spraying cabinet and transferred to the glasshouse where the tolerance level of the genotypes was assessed by measuring the relative reduction in chlorophyll content of the leaves. The decrease in chlorophyll content of the treated plants compared to the control was regarded as the phytotoxic effects of metribuzin. GWAS analysis following a mixed linear model revealed 19 genomic regions with significant marker-trait associations (MTAs), including ten on chromosome 6A, three on chromosome 2B, and one on chromosomes 3A, 5B, 6B 6D, 7A, and 7B, respectively. Sequences of the significant markers were blasted against the wheat genome, IWGSC RefSeq V1.0, and candidate genes having annotations related to herbicide tolerance in wheat, especially in pathways reported to be involved in metribuzin tolerance, such as cytochrome P450 pathways and ATP Binding Cassette (ABC) superfamilies, were identified in these genomic regions. These included TraesCS6A01G028800, TraesCS6A02G353700, TraesCS6A01G326200, TraesCS7A02G331000, and TraesCS2B01G465200. These genomic regions were validated on 30 top tolerant and 30 most susceptible genotypes using the five closest SSR makers to the flanked SNPs. Sufficient polymorphism was detected on two markers (wms193 and barc1036) that were found to differentiate between the susceptible and tolerant alleles and a t-test analysis of the phenotypic data shows a significant (value of p < 0.001) difference suggesting that these markers can be used for marker-assisted selection (MAS) in metribuzin studies and wheat breeding programs.

Progenitor species hold untapped diversity for potential climate-responsive traits for use in wheat breeding and crop improvement

Climate change will have numerous impacts on crop production worldwide necessitating a broadening of the germplasm base required to source and incorporate novel traits. Major variation exists in crop progenitor species for seasonal adaptation, photosynthetic characteristics, and root system architecture. Wheat is crucial for securing future food and nutrition security and its evolutionary history and progenitor diversity offer opportunities to mine favourable functional variation in the primary gene pool. Here we provide a review of the status of characterisation of wheat progenitor variation and the potential to use this knowledge to inform the use of variation in other cereal crops. Although significant knowledge of progenitor variation has been generated, we make recommendations for further work required to systematically characterise underlying genetics and physiological mechanisms and propose steps for effective use in breeding. This will enable targeted exploitation of useful variation, supported by the growing portfolio of genomics and accelerated breeding approaches. The knowledge and approaches generated are also likely to be useful across wider crop improvement.

The wheat AGL6-like MADS-box gene is a master regulator for floral organ identity and a target for spikelet meristem development manipulation

The AGAMOUS-LIKE6 (AGL6)-like genes are ancient MADS-box genes and are functionally studied in a few model plants. The knowledge of these genes in wheat remains limited. Here, by studying a 'double homoeolog mutant' of the AGL6 gene in tetraploid wheat, we showed that AGL6 was required for the development of all four whorls of floral organs with dosage-dependent effect on floret fertility. Yeast two-hybrid analyses detected interactions of AGL6 with all classes of MADS-box proteins in the ABCDE model for floral organ development. AGL6 was found to interact with several additional proteins, including the G protein β and γ (DEP1) subunits. Analysis of the DEP1-B mutant showed a significant reduction in spikelet number per spike in tetraploid wheat, while overexpression of AGL6 in common wheat increased the spikelet number per spike and hence the grain number per spike. RNA-seq analysis identified the regulation of several meristem activity genes by AGL6, such as FUL2 and TaMADS55. Our work therefore extensively updated the wheat ABCDE model and proposed an alternative approach to improve wheat grain yield by manipulating the AGL6 gene.

Challenges for a Massive Implementation of Phenomics in Plant Breeding Programs

Due to climate change and expected food shortage in the coming decades, not only will it be necessary to develop cultivars with greater tolerance to environmental stress, but it is also imperative to reduce breeding cycle time. In addition to yield evaluation, plant breeders resort to many sensory assessments and some others of intermediate complexity. However, to develop cultivars better adapted to current/future constraints, it is necessary to incorporate a new set of traits, such as morphophysiological and physicochemical attributes, information relevant to the successful selection of genotypes or parents. Unfortunately, because of the large number of genotypes to be screened, measurements with conventional equipment are unfeasible, especially under field conditions. High-throughput plant phenotyping (HTPP) facilitates collecting a significant amount of data quickly; however, it is necessary to transform all this information (e.g., plant reflectance) into helpful descriptors to the breeder. To the extent that a holistic characterization of the plant (phenomics) is performed in challenging environments, it will be possible to select the best genotypes (forward phenomics) objectively but also understand why the said individual differs from the rest (reverse phenomics). Unfortunately, several elements had prevented phenomics from developing as desired. Consequently, a new set of prediction/validation methodologies, seasonal ambient information, and the fusion of data matrices (e.g., genotypic and phenotypic information) need to be incorporated into the modeling. In this sense, for the massive implementation of phenomics in plant breeding, it will be essential to count an interdisciplinary team that responds to the urgent need to release material with greater capacity to tolerate environmental stress. Therefore, breeding programs should (i) be more efficient (e.g., early discarding of unsuitable material), (ii) have shorter breeding cycles (fewer crosses to achieve the desired cultivar), and (iii) be more productive, increasing the probability of success at the end of the breeding process (percentage of cultivars released to the number of initial crosses).

Systems biology for crop improvement

In recent years, generation of large-scale data from genome, transcriptome, proteome, metabolome, epigenome, and others, has become routine in several plant species. Most of these datasets in different crop species, however, were studied independently and as a result, full insight could not be gained on the molecular basis of complex traits and biological networks. A systems biology approach involving integration of multiple omics data, modeling, and prediction of the cellular functions is required to understand the flow of biological information that underlies complex traits. In this context, systems biology with multiomics data integration is crucial and allows a holistic understanding of the dynamic system with the different levels of biological organization interacting with external environment for a phenotypic expression. Here, we present recent progress made in the area of various omics studies-integrative and systems biology approaches with a special focus on application to crop improvement. We have also discussed the challenges and opportunities in multiomics data integration, modeling, and understanding of the biology of complex traits underpinning yield and stress tolerance in major cereals and legumes.

Myb10-D confers PHS-3D resistance to pre-harvest sprouting by regulating NCED in ABA biosynthesis pathway of wheat

Pre-harvest sprouting (PHS), the germination of grain before harvest, is a serious problem resulting in wheat yield and quality losses. Here, we mapped the PHS resistance gene PHS-3D from synthetic hexaploid wheat to a 2.4 Mb presence-absence variation (PAV) region and found that its resistance effect was attributed to the pleiotropic Myb10-D by integrated omics and functional analyses. Three haplotypes were detected in this PAV region among 262 worldwide wheat lines and 16 Aegilops tauschii, and the germination percentages of wheat lines containing Myb10-D was approximately 40% lower than that of the other lines. Transcriptome and metabolome profiling indicated that Myb10-D affected the transcription of genes in both the flavonoid and abscisic acid (ABA) biosynthesis pathways, which resulted in increases in flavonoids and ABA in transgenic wheat lines. Myb10-D activates 9-cis-epoxycarotenoid dioxygenase (NCED) by biding the secondary wall MYB-responsive element (SMRE) to promote ABA biosynthesis in early wheat seed development stages. We revealed that the newly discovered function of Myb10-D confers PHS resistance by enhancing ABA biosynthesis to delay germination in wheat. The PAV harboring Myb10-D associated with grain color and PHS will be useful for understanding and selecting white grained PHS resistant wheat cultivars.

Evaluation of genetic structure in European wheat cultivars and advanced breeding lines using high-density genotyping-by-sequencing approach

BACKGROUND

The genetic diversity and gene pool characteristics must be clarified for efficient genome-wide association studies, genomic selection, and hybrid breeding. The aim of this study was to evaluate the genetic structure of 509 wheat accessions representing registered varieties and advanced breeding lines via the high-density genotyping-by-sequencing approach.

RESULTS

More than 30% of 13,499 SNP markers representing 2162 clusters were mapped to genes, whereas 22.50% of 26,369 silicoDArT markers overlapped with coding sequences and were linked in 3527 blocks. Regarding hexaploidy, perfect sequence matches following BLAST searches were not sufficient for the unequivocal mapping to unique loci. Moreover, allelic variations in homeologous loci interfered with heterozygosity calculations for some markers. Analyses of the major genetic changes over the last 27 years revealed the selection pressure on orthologs of the gibberellin biosynthesis-related GA2 gene and the senescence-associated SAG12 gene. A core collection representing the wheat population was generated for preserving germplasm and optimizing breeding programs.

CONCLUSIONS

Our results confirmed considerable differences among wheat subgenomes A, B and D, with D characterized by the lowest diversity but the highest LD. They revealed genomic regions that have been targeted by breeding.

Transcriptomic profiling of wheat near-isogenic lines reveals candidate genes on chromosome 3A for pre-harvest sprouting resistance

BACKGROUND

Pre-harvest sprouting (PHS) in wheat can cause severe damage to both grain yield and quality. Resistance to PHS is a quantitative trait controlled by many genes located across all 21 wheat chromosomes. The study targeted a large-effect quantitative trait locus (QTL) QPhs.ccsu-3A.1 for PHS resistance using several sets previously developed near-isogenic lines (NILs). Two pairs of NILs with highly significant phenotypic differences between the isolines were examined by RNA sequencing for their transcriptomic profiles on developing seeds at 15, 25 and 35 days after pollination (DAP) to identify candidate genes underlying the QTL and elucidate gene effects on PHS resistance. At each DAP, differentially expressed genes (DEGs) between the isolines were investigated.

RESULTS

Gene ontology and KEGG pathway enrichment analyses of key DEGs suggested that six candidate genes underlie QPhs.ccsu-3A.1 responsible for PHS resistance in wheat. Candidate gene expression was further validated by quantitative RT-PCR. Within the targeted QTL interval, 16 genetic variants including five single nucleotide polymorphisms (SNPs) and 11 indels showed consistent polymorphism between resistant and susceptible isolines.

CONCLUSIONS

The targeted QTL is confirmed to harbor core genes related to hormone signaling pathways that can be exploited as a key genomic region for marker-assisted selection. The candidate genes and SNP/indel markers detected in this study are valuable resources for understanding the mechanism of PHS resistance and for marker-assisted breeding of the trait in wheat.

Genome-wide association study of yield and related traits in common wheat under salt-stress conditions

BACKGROUND

Soil salinization is a major threat to wheat production. It is essential to understand the genetic basis of salt tolerance for breeding and selecting new salt-tolerant cultivars that have the potential to increase wheat yield.

RESULT

In this study, a panel of 191 wheat accessions was subjected to genome wide association study (GWAS) to identify SNP markers linked with adult-stage characters. The population was genotyped by Wheat660K SNP array and eight phenotype traits were investigated under low and high salinity environments for three consecutive years. A total of 389 SNPs representing 11 QTLs were significantly associated with plant height, spike number, spike length, grain number, thousand kernels weight, yield and biological mass under different salt treatments, with the phenotypic explanation rate (R²) ranging from 9.14 to 50.45%. Of these, repetitive and pleiotropic loci on chromosomes 4A, 5A, 5B and 7A were significantly linked to yield and yield related traits such as thousand kernels weight, spike number, spike length, grain number and so on under low salinity conditions. Spike length-related loci were mainly located on chromosomes 1B, 3B, 5B and 7A under different salt treatments. Two loci on chromosome 4D and 5A were related with plant height in low and high salinity environment, respectively. Three salt-tolerant related loci were confirmed to be important in two bi-parental populations. Distribution of favorable haplotypes indicated that superior haplotypes of pleiotropic loci on group-5 chromosomes were strongly selected and had potential for increasing wheat salt tolerance. A total of 14 KASP markers were developed for nine loci associating with yield and related traits to improve the selection efficiency of wheat salt-tolerance breeding.

CONCLUSION

Utilizing a Wheat660K SNPs chip, QTLs for yield and its related traits were detected under salt treatment in a natural wheat population. Important salt-tolerant related loci were validated in RIL and DH populations. This study provided reliable molecular markers that could be crucial for marker-assisted selection in wheat salt tolerance breeding programs.

An Efficient Brome mosaic virus-Based Gene Silencing Protocol for Hexaploid Wheat (Triticum aestivum L.)

Virus-induced gene silencing (VIGS) is a rapid and powerful method to evaluate gene function, especially for species like hexaploid wheat that have large, redundant genomes and are difficult and time-consuming to transform. The Brome mosaic virus (BMV)-based VIGS vector is widely used in monocotyledonous species but not wheat. Here we report the establishment of a simple and effective VIGS procedure in bread wheat using BMVCP5, the most recently improved BMV silencing vector, and wheat genes PHYTOENE DESATURASE (TaPDS) and PHOSPHATE2 (TaPHO2) as targets. Time-course experiments revealed that smaller inserts (~100 nucleotides, nt) were more stable in BMVCP5 and conferred higher silencing efficiency and longer silencing duration, compared with larger inserts. When using a 100-nt insert and a novel coleoptile inoculation method, BMVCP5 induced extensive silencing of TaPDS transcript and a visible bleaching phenotype in the 2nd to 5th systemically-infected leaves from nine to at least 28 days post inoculation (dpi). For TaPHO2, the ability of BMVCP5 to simultaneously silence all three homoeologs was demonstrated. To investigate the feasibility of BMV VIGS in wheat roots, ectopically expressed enhanced GREEN FLUORESCENT PROTEIN (eGFP) in a transgenic wheat line was targeted for silencing. Silencing of eGFP fluorescence was observed in both the maturation and elongation zones of roots. BMVCP5 mediated significant silencing of eGFP and TaPHO2 mRNA expression in roots at 14 and 21 dpi, and TaPHO2 silencing led to the doubling of inorganic phosphate concentration in the 2nd through 4th systemic leaves. All 54 wheat cultivars screened were susceptible to BMV infection. BMVCP5-mediated TaPDS silencing resulted in the expected bleaching phenotype in all eight cultivars examined, and decreased TaPDS transcript was detected in all three cultivars examined. This BMVCP5 VIGS technology may serve as a rapid and effective functional genomics tool for high-throughput gene function studies in aerial and root tissues and in many wheat cultivars.

An improved and robust method to efficiently deplete repetitive elements from complex plant genomes

Genome size and complexity often present major challenges to genome-based approaches in crop plants and other agricultural species. For instance, repetitive sequences comprise 80% to 90% of the genome of hexaploid wheat, which has a haploid genome size of approximately 17 Gb. In this study, we developed an improved design and procedure for short-read library preparation that uses a modified adaptor and duplex-specific nuclease (DSN) for the efficient elimination of highly repeated sequence elements within genomes. The improved adapter, which has a hairpin-like form for stability, was constructed from truncated sequences adjacent to the original Illumina TruSeq adapter and can be converted to a full-length adapter structure during PCR amplification. Using the hairpin-structured adaptor, we prepared randomly sheared genomic libraries from rice and diploid, tetraploid, and hexaploid wheat cultivars and evaluated the efficiency of DSN for the enzymatic depletion of repetitive elements. According to real-time quantitative PCR analysis, the relative abundances of 18S and 25S ribosomal DNA decreased respectively to 1.15% and 3.54% in rice and 1.70%-1.95% and 14.71%-20.01% in the three wheat cultivars. Whole-genome sequencing analysis of a diploid wheat cultivar, KU104-1, indicated that DSN treatment with the designed hairpin-structured adapter dramatically reduced highly repetitive elements, such as Ty1-Copia and Ty3-Gypsy retrotransposons and DNA transposons, within the genome, while sequencing reads derived from low-copy genes and protein coding sequences increased more than 50%. Our new procedure should be useful not only for wheat genomes but also for other agricultural plant species with relatively large and complex genomes.