Overexpression of the WAPO-A1 gene increases the number of spikelets per spike in bread wheat.

Two homoeologous QTLs for number of spikelets per spike (SPS) were mapped on chromosomes 7AL and 7BL using two wheat MAGIC populations. Sets of lines contrasting for the QTL on 7AL were developed which allowed for the validation and fine mapping of the 7AL QTL and for the identification of a previously described candidate gene, WHEAT ORTHOLOG OF APO1 (WAPO1). Using transgenic overexpression in both a low and a high SPS line, we provide a functional validation for the role of this gene in determining SPS also in hexaploid wheat. We show that the expression levels of this gene positively correlate with SPS in multiple MAGIC founder lines under field conditions as well as in transgenic lines grown in the greenhouse. This work highlights the potential use of WAPO1 in hexaploid wheat for further yield increases. The impact of WAPO1 and SPS on yield depends on other genetic and environmental factors, hence, will require a finely balanced expression level to avoid the development of detrimental pleiotropic phenotypes.

Optical and physical mapping with local finishing enables megabase-scale resolution of agronomically important regions in the wheat genome.

BACKGROUND: Numerous scaffold-level sequences for wheat are now being released and, in this context, we report on a strategy for improving the overall assembly to a level comparable to that of the human genome. RESULTS: Using chromosome 7A of wheat as a model, sequence-finished megabase-scale sections of this chromosome were established by combining a new independent assembly using a bacterial artificial chromosome (BAC)-based physical map, BAC pool paired-end sequencing, chromosome-arm-specific mate-pair sequencing and Bionano optical mapping with the International Wheat Genome Sequencing Consortium RefSeq v1.0 sequence and its underlying raw data. The combined assembly results in 18 super-scaffolds across the chromosome. The value of finished genome regions is demonstrated for two approximately 2.5 Mb regions associated with yield and the grain quality phenotype of fructan carbohydrate grain levels. In addition, the 50 Mb centromere region analysis incorporates cytological data highlighting the importance of non-sequence data in the assembly of this complex genome region. CONCLUSIONS: Sufficient genome sequence information is shown to now be available for the wheat community to produce sequence-finished releases of each chromosome of the reference genome. The high-level completion identified that an array of seven fructosyl transferase genes underpins grain quality and that yield attributes are affected by five F-box-only-protein-ubiquitin ligase domain and four root-specific lipid transfer domain genes. The completed sequence also includes the centromere.

Zebularine treatment is associated with deletion of FT-B1 leading to an increase in spikelet number in bread wheat.

The number of rachis nodes (spikelets) on a wheat spike is a component of grain yield that correlates with flowering time. The genetic basis regulating flowering in cereals is well understood, but there are reports that flowering time can be modified at a high frequency by selective breeding, suggesting that it may be regulated by both epigenetic and genetic mechanisms. We investigated the role of DNA methylation in regulating spikelet number and flowering time by treating a semi-spring wheat with the demethylating agent, Zebularine. Three lines with a heritable increase in spikelet number were identified. The molecular basis for increased spikelet number was not determined in 2 lines, but the phenotype showed non-Mendelian inheritance, suggesting that it could have an epigenetic basis. In the remaining line, the increased spikelet phenotype behaved as a Mendelian recessive trait and late flowering was associated with a deletion encompassing the floral promoter, FT-B1. Deletion of FT-B1 delayed the transition to reproductive growth, extended the duration of spike development, and increased spikelet number under different temperature regimes and photoperiod. Transiently disrupting DNA methylation can generate novel flowering behaviour in wheat, but these changes may not be sufficiently stable for use in breeding programs.

TEOSINTE BRANCHED1 Regulates Inflorescence Architecture and Development in Bread Wheat (Triticum aestivum).

The flowers of major cereals are arranged on reproductive branches known as spikelets, which group together to form an inflorescence. Diversity for inflorescence architecture has been exploited during domestication to increase crop yields, and genetic variation for this trait has potential to further boost grain production. Multiple genes that regulate inflorescence architecture have been identified by studying alleles that modify gene activity or dosage; however, little is known in wheat. Here, we show TEOSINTE BRANCHED1 (TB1) regulates inflorescence architecture in bread wheat (Triticum aestivum) by investigating lines that display a form of inflorescence branching known as "paired spikelets." We show that TB1 interacts with FLOWERING LOCUS T1 and that increased dosage of TB1 alters inflorescence architecture and growth rate in a process that includes reduced expression of meristem identity genes, with allelic diversity for TB1 found to associate genetically with paired spikelet development in modern cultivars. We propose TB1 coordinates formation of axillary spikelets during the vegetative to floral transition and that alleles known to modify dosage or function of TB1 could help increase wheat yields.

Increased accuracy of starch granule type quantification using mixture distributions.

BACKGROUND: The proportion of granule types in wheat starch is an important characteristic that can affect its functionality. It is widely accepted that granule types are either large, disc-shaped A-type granules or small, spherical B-type granules. Additionally, there are some reports of the tiny C-type granules. The differences between these granule types are due to its carbohydrate composition and crystallinity which is highly, but not perfectly, correlated with the granule size. A majority of the studies that have considered granule types analyse them based on a size threshold rather than chemical composition. This is understandable due to the expense of separating starch into different types. While the use of a size threshold to classify granule type is a low-cost measure, this results in misclassification. We present an alternative, statistical method to quantify the proportion of granule types by a fit of the mixture distribution, along with an R package, a web based app and a video tutorial for how to use the web app to enable its straightforward application. RESULTS: Our results show that the reliability of the genotypic effects increase approximately 60% using the proportions of the A-type and B-type granule estimated by the mixture distribution over the standard size-threshold measure. Although there was a marginal drop in reliability for C-type granules. The latter is likely due to the low observed genetic variance for C-type granules. CONCLUSIONS: The determination of the proportion of granule types from size-distribution is better achieved by using the mixing probabilities from the fit of the mixture distribution rather than using a size-threshold.

Transcriptomic analysis of wheat near-isogenic lines identifies PM19-A1 and A2 as candidates for a major dormancy QTL.

BACKGROUND: Next-generation sequencing technologies provide new opportunities to identify the genetic components responsible for trait variation. However, in species with large polyploid genomes, such as bread wheat, the ability to rapidly identify genes underlying quantitative trait loci (QTL) remains non-trivial. To overcome this, we introduce a novel pipeline that analyses, by RNA-sequencing, multiple near-isogenic lines segregating for a targeted QTL. RESULTS: We use this approach to characterize a major and widely utilized seed dormancy QTL located on chromosome 4AL. It exploits the power and mapping resolution afforded by large multi-parent mapping populations, whilst reducing complexity by using multi-allelic contrasts at the targeted QTL region. Our approach identifies two adjacent candidate genes within the QTL region belonging to the ABA-induced Wheat Plasma Membrane 19 family. One of them, PM19-A1, is highly expressed during grain maturation in dormant genotypes. The second, PM19-A2, shows changes in sequence causing several amino acid alterations between dormant and non-dormant genotypes. We confirm that PM19 genes are positive regulators of seed dormancy. CONCLUSIONS: The efficient identification of these strong candidates demonstrates the utility of our transcriptomic pipeline for rapid QTL to gene mapping. By using this approach we are able to provide a comprehensive genetic analysis of the major source of grain dormancy in wheat. Further analysis across a diverse panel of bread and durum wheats indicates that this important dormancy QTL predates hexaploid wheat. The use of these genes by wheat breeders could assist in the elimination of pre-harvest sprouting in wheat.

The genetics of rhizosheath size in a multiparent mapping population of wheat.

Rhizosheaths comprise soil that adheres to plant roots and, in some species, are indicative of root hair length. In this study, the genetics of rhizosheath size in wheat was investigated by screening the progeny of multiparent advanced generation intercrosses (MAGIC). Two MAGIC populations were screened for rhizosheath size using a high throughput method. One MAGIC population was developed from intercrosses between four parents (4-way) and the other from intercrosses between eight parents (8-way). Transgressive segregation for rhizosheath size was observed in both the 4-way and 8-way MAGIC populations. A quantitative trait loci (QTL) analysis of the 4-way population identified six major loci located on chromosomes 2B, 4D, 5A, 5B, 6A, and 7A together accounting for 42% of the variation in rhizosheath size. Rhizosheath size was strongly correlated with root hair length and was robust across different soil types in the absence of chemical constraints. Rhizosheath size in the MAGIC populations was a reliable surrogate for root hair length and, therefore, the QTL identified probably control root hair elongation. Members of the basic helix-loop-helix family of transcription factors have previously been identified to regulate root hair length in Arabidopsis and rice. Since several wheat members of the basic helix-loop-helix family of genes are located within or near the QTL, these genes are candidates for controlling the long root hair trait. The QTL for rhizosheath size identified in this study provides the opportunity to implement marker-assisted selection to increase root hair length for improved phosphate acquisition in wheat.

MAGIC populations in crops: current status and future prospects.

KEY MESSAGE: MAGIC populations present novel challenges and opportunities in crops due to their complex pedigree structure. They offer great potential both for dissecting genomic structure and for improving breeding populations. The past decade has seen the rise of multiparental populations as a study design offering great advantages for genetic studies in plants. The genetic diversity of multiple parents, recombined over several generations, generates a genetic resource population with large phenotypic diversity suitable for high-resolution trait mapping. While there are many variations on the general design, this review focuses on populations where the parents have all been inter-mated, typically termed Multi-parent Advanced Generation Intercrosses (MAGIC). Such populations have already been created in model animals and plants, and are emerging in many crop species. However, there has been little consideration of the full range of factors which create novel challenges for design and analysis in these populations. We will present brief descriptions of large MAGIC crop studies currently in progress to motivate discussion of population construction, efficient experimental design, and genetic analysis in these populations. In addition, we will highlight some recent achievements and discuss the opportunities and advantages to exploit the unique structure of these resources post-QTL analysis for gene discovery.

Genome-wide association mapping in plants.

KEY MESSAGE: We present new association mapping methods which address the unique challenges of analyzing genome-wide data from multi-environment plant studies. Association studies on a genome-wide scale are being performed in plants. Unlike human studies, plant studies contain replicates whose data may be recorded across different environments. Plant studies also often employ elaborate experimental designs for controlling extraneous phenotypic variation. As a result, the genome-wide analysis of data from plant studies can be challenging. In this paper, we present QK-based association mapping for the analysis of data from plant association studies. In doing so, we have developed: (a) a general multivariate QK framework for association mapping in plant studies of arbitrary complexity; (b) a new weighted two-stage analysis approach for QK-based association mapping; (c) a heuristic procedure for determining when two-stage analysis is appropriate; and (d) a Monte Carlo sampling procedure for controlling the genome-wide type I error rate. We conduct a simulation study to evaluate the performance of our genome-wide mapping technique. We also analyze data from a multi-environment association study in wheat.

A genetic strategy generating wheat with very high amylose content.

Resistant starch (RS), a type of dietary fibre, plays an important role in human health; however, the content of RS in most modern processed starchy foods is low. Cereal starch, when structurally manipulated through a modified starch biosynthetic pathway to greatly increase the amylose content, could be an important food source of RS. Transgenic studies have previously revealed the requirement of simultaneous down-regulation of two starch branching enzyme (SBE) II isoforms both located on the long arm of chromosome 2, namely SBEIIa and SBEIIb, to elevate the amylose content in wheat from ~25% to ~75%. The current study revealed close proximity of genes encoding SBEIIa and SBEIIb isoforms in wheat with a genetic distance of 0.5 cM on chromosome 2B. A series of deletion and single nucleotide polymorphism (SNP) loss of function alleles in SBEIIa, SBEIIb or both was isolated from two different wheat populations. A breeding strategy to combine deletions and SNPs generated wheat genotypes with altered expression levels of SBEIIa and SBEIIb, elevating the amylose content to an unprecedented ~85%, with a marked concomitant increase in RS content. Biochemical assays were used to confirm the complete absence in the grain of expression of SBEIIa from all three genomes in combination with the absence of SBEIIb from one of the genomes.

Ppd-1 is a key regulator of inflorescence architecture and paired spikelet development in wheat.

The domestication of cereal crops such as wheat, maize, rice and barley has included the modification of inflorescence architecture to improve grain yield and ease harvesting(1). Yield increases have often been achieved through modifying the number and arrangement of spikelets, which are specialized reproductive branches that form part of the inflorescence. Multiple genes that control spikelet development have been identified in maize, rice and barley(2-5). However, little is known about the genetic underpinnings of this process in wheat. Here, we describe a modified spikelet arrangement in wheat, termed paired spikelets. Combining comprehensive QTL and mutant analyses, we show that Photoperiod-1 (Ppd-1), a pseudo-response regulator gene that controls photoperiod-dependent floral induction, has a major inhibitory effect on paired spikelet formation by regulating the expression of FLOWERING LOCUS T (FT)(6,7). These findings show that modulated expression of the two important flowering genes, Ppd-1 and FT, can be used to form a wheat inflorescence with a more elaborate arrangement and increased number of grain producing spikelets.

Transcriptome and allele specificity associated with a 3BL locus for Fusarium crown rot resistance in bread wheat.

Fusarium pathogens cause two major diseases in cereals, Fusarium crown rot (FCR) and head blight (FHB). A large-effect locus conferring resistance to FCR disease was previously located to chromosome arm 3BL (designated as Qcrs-3B) and several independent sets of near isogenic lines (NILs) have been developed for this locus. In this study, five sets of the NILs were used to examine transcriptional changes associated with the Qcrs-3B locus and to identify genes linked to the resistance locus as a step towards the isolation of the causative gene(s). Of the differentially expressed genes (DEGs) detected between the NILs, 12.7% was located on the single chromosome 3B. Of the expressed genes containing SNP (SNP-EGs) detected, 23.5% was mapped to this chromosome. Several of the DEGs and SNP-EGs are known to be involved in host-pathogen interactions, and a large number of the DEGs were among those detected for FHB in previous studies. Of the DEGs detected, 22 were mapped in the Qcrs-3B interval and they included eight which were detected in the resistant isolines only. The enrichment of DEG, and not necessarily those containing SNPs between the resistant and susceptible isolines, around the Qcrs-3B locus is suggestive of local regulation of this region by the resistance allele. Functions for 13 of these DEGs are known. Of the SNP-EGs, 28 were mapped in the Qcrs-3B interval and biological functions for 16 of them are known. These results provide insights into responses regulated by the 3BL locus and identify a tractable number of target genes for fine mapping and functional testing to identify the causative gene(s) at this QTL.

Computationally efficient map construction in the presence of segregation distortion.

KEY MESSAGE: We present a novel estimator for map construction in the presence of segregation distortion which is highly computationally efficient. For multi-parental designs this estimator outperforms methods that do not account for segregation distortion, at no extra computational cost. Inclusion of genetic markers exhibiting segregation distortion in a linkage map can result in biased estimates of genetic distance and distortion of map positions. Removal of distorted markers is hence a typical filtering criterion; however, this may result in exclusion of biologically interesting regions of the genome such as introgressions and translocations. Estimation of additional parameters characterizing the distortion is computationally slow, as it relies on estimation via the Expectation Maximization algorithm or a higher dimensional numerical optimisation. We propose a robust M-estimator (RM) capable of handling tens of thousands of distorted markers from a single linkage group. We show via simulation that for multi-parental designs the RM estimator can perform much better than uncorrected estimation, at no extra computational cost. We then apply the RM estimator to chromosome 2B in wheat in a multi-parent population segregating for the Sr36 introgression, a known transmission distorter. The resulting map contains over 700 markers, and is consistent with maps constructed from crosses which do not exhibit segregation distortion.

Characterizing uncertainty in high-density maps from multiparental populations.

Multiparental populations are of considerable interest in high-density genetic mapping due to their increased levels of polymorphism and recombination relative to biparental populations. However, errors in map construction can have significant impact on QTL discovery in later stages of analysis, and few methods have been developed to quantify the uncertainty attached to the reported order of markers or intermarker distances. Current methods are computationally intensive or limited to assessing uncertainty only for order or distance, but not both simultaneously. We derive the asymptotic joint distribution of maximum composite likelihood estimators for intermarker distances. This approach allows us to construct hypothesis tests and confidence intervals for simultaneously assessing marker-order instability and distance uncertainty. We investigate the effects of marker density, population size, and founder distribution patterns on map confidence in multiparental populations through simulations. Using these data, we provide guidelines on sample sizes necessary to map markers at sub-centimorgan densities with high certainty. We apply these approaches to data from a bread wheat Multiparent Advanced Generation Inter-Cross (MAGIC) population genotyped using the Illumina 9K SNP chip to assess regions of uncertainty and validate them against the recently released pseudomolecule for the wheat chromosome 3B.

Whole-genome analysis of multienvironment or multitrait QTL in MAGIC.

Multiparent Advanced Generation Inter-Cross (MAGIC) populations are now being utilized to more accurately identify the underlying genetic basis of quantitative traits through quantitative trait loci (QTL) analyses and subsequent gene discovery. The expanded genetic diversity present in such populations and the amplified number of recombination events mean that QTL can be identified at a higher resolution. Most QTL analyses are conducted separately for each trait within a single environment. Separate analysis does not take advantage of the underlying correlation structure found in multienvironment or multitrait data. By using this information in a joint analysis-be it multienvironment or multitrait - it is possible to gain a greater understanding of genotype- or QTL-by-environment interactions or of pleiotropic effects across traits. Furthermore, this can result in improvements in accuracy for a range of traits or in a specific target environment and can influence selection decisions. Data derived from MAGIC populations allow for founder probabilities of all founder alleles to be calculated for each individual within the population. This presents an additional layer of complexity and information that can be utilized to identify QTL. A whole-genome approach is proposed for multienvironment and multitrait QTL analysis in MAGIC. The whole-genome approach simultaneously incorporates all founder probabilities at each marker for all individuals in the analysis, rather than using a genome scan. A dimension reduction technique is implemented, which allows for high-dimensional genetic data. For each QTL identified, sizes of effects for each founder allele, the percentage of genetic variance explained, and a score to reflect the strength of the QTL are found. The approach was demonstrated to perform well in a small simulation study and for two experiments, using a wheat MAGIC population.

GrainScan: a low cost, fast method for grain size and colour measurements.

BACKGROUND: Measuring grain characteristics is an integral component of cereal breeding and research into genetic control of seed development. Measures such as thousand grain weight are fast, but do not give an indication of variation within a sample. Other methods exist for detailed analysis of grain size, but are generally costly and very low throughput. Grain colour analysis is generally difficult to perform with accuracy, and existing methods are expensive and involved. RESULTS: We have developed a software method to measure grain size and colour from images captured with consumer level flatbed scanners, in a robust, standardised way. The accuracy and precision of the method have been demonstrated through screening wheat and Brachypodium distachyon populations for variation in size and colour. CONCLUSION: By using GrainScan, cheap and fast measurement of grain colour and size will enable plant research programs to gain deeper understanding of material, where limited or no information is currently available.

Whole-genome QTL analysis for MAGIC.

KEY MESSAGE: An efficient whole genome method of QTL analysis is presented for Multi-parent advanced generation integrated crosses. Multi-parent advanced generation inter-cross (MAGIC) populations have been developed for mice and several plant species and are useful for the genetic dissection of complex traits. The analysis of quantitative trait loci (QTL) in these populations presents some additional challenges compared with traditional mapping approaches. In particular, pedigree and marker information need to be integrated and founder genetic data needs to be incorporated into the analysis. Here, we present a method for QTL analysis that utilizes the probability of inheriting founder alleles across the whole genome simultaneously, either for intervals or markers. The probabilities can be found using three-point or Hidden Markov Model (HMM) methods. This whole-genome approach is evaluated in a simulation study and it is shown to be a powerful method of analysis. The HMM probabilities lead to low rates of false positives and low bias of estimated QTL effect sizes. An implementation of the approach is available as an R package. In addition, we illustrate the approach using a bread wheat MAGIC population.

Characterization of polyploid wheat genomic diversity using a high-density 90,000 single nucleotide polymorphism array.

High-density single nucleotide polymorphism (SNP) genotyping arrays are a powerful tool for studying genomic patterns of diversity, inferring ancestral relationships between individuals in populations and studying marker-trait associations in mapping experiments. We developed a genotyping array including about 90,000 gene-associated SNPs and used it to characterize genetic variation in allohexaploid and allotetraploid wheat populations. The array includes a significant fraction of common genome-wide distributed SNPs that are represented in populations of diverse geographical origin. We used density-based spatial clustering algorithms to enable high-throughput genotype calling in complex data sets obtained for polyploid wheat. We show that these model-free clustering algorithms provide accurate genotype calling in the presence of multiple clusters including clusters with low signal intensity resulting from significant sequence divergence at the target SNP site or gene deletions. Assays that detect low-intensity clusters can provide insight into the distribution of presence-absence variation (PAV) in wheat populations. A total of 46 977 SNPs from the wheat 90K array were genetically mapped using a combination of eight mapping populations. The developed array and cluster identification algorithms provide an opportunity to infer detailed haplotype structure in polyploid wheat and will serve as an invaluable resource for diversity studies and investigating the genetic basis of trait variation in wheat.

Use of a large multiparent wheat mapping population in genomic dissection of coleoptile and seedling growth.

Identification of alleles towards the selection for improved seedling vigour is a key objective of many wheat breeding programmes. A multiparent advanced generation intercross (MAGIC) population developed from four commercial spring wheat cultivars (cvv. Baxter, Chara, Westonia and Yitpi) and containing ca. 1000 F(2) -derived, F(6:7) RILs was assessed at two contrasting soil temperatures (12 and 20 °C) for shoot length and coleoptile characteristics length and thickness. Narrow-sense heritabilities were high for coleoptile and shoot length (h(2) = 0.68-0.70), indicating a strong genetic basis for the differences among progeny. Genotypic variation was large, and distributions of genotype means were approximately Gaussian with evidence for transgressive segregation for all traits. A number of significant QTL were identified for all early growth traits, and these were commonly repeatable across the different soil temperatures. The largest negative effects on coleoptile lengths were associated with Rht-B1b (-8.2%) and Rht-D1b (-10.9%) dwarfing genes varying in the population. Reduction in coleoptile length with either gene was particularly large at the warmer soil temperature. Other large QTL for coleoptile length were identified on chromosomes 1A, 2B, 4A, 5A and 6B, but these were relatively smaller than allelic effects at the Rht-B1 and Rht-D1 loci. A large coleoptile length effect allele (a = 5.3 mm at 12 °C) was identified on chromosome 1AS despite the relatively shorter coleoptile length of the donor Yitpi. Strong, positive genetic correlations for coleoptile and shoot lengths (r(g) = 0.85-0.90) support the co-location of QTL for these traits and suggest a common physiological basis for both. The multiparent population has enabled the identification of promising shoot and coleoptile QTL despite the potential for the confounding of large effect dwarfing gene alleles present in the commercial parents. The incidence of these alleles in commercial wheat breeding programmes should facilitate their ready implementation in selection of varieties with improved establishment and early growth.

Genome-wide comparative diversity uncovers multiple targets of selection for improvement in hexaploid wheat landraces and cultivars.

Domesticated crops experience strong human-mediated selection aimed at developing high-yielding varieties adapted to local conditions. To detect regions of the wheat genome subject to selection during improvement, we developed a high-throughput array to interrogate 9,000 gene-associated single-nucleotide polymorphisms (SNP) in a worldwide sample of 2,994 accessions of hexaploid wheat including landraces and modern cultivars. Using a SNP-based diversity map we characterized the impact of crop improvement on genomic and geographic patterns of genetic diversity. We found evidence of a small population bottleneck and extensive use of ancestral variation often traceable to founders of cultivars from diverse geographic regions. Analyzing genetic differentiation among populations and the extent of haplotype sharing, we identified allelic variants subjected to selection during improvement. Selective sweeps were found around genes involved in the regulation of flowering time and phenology. An introgression of a wild relative-derived gene conferring resistance to a fungal pathogen was detected by haplotype-based analysis. Comparing selective sweeps identified in different populations, we show that selection likely acts on distinct targets or multiple functionally equivalent alleles in different portions of the geographic range of wheat. The majority of the selected alleles were present at low frequency in local populations, suggesting either weak selection pressure or temporal variation in the targets of directional selection during breeding probably associated with changing agricultural practices or environmental conditions. The developed SNP chip and map of genetic variation provide a resource for advancing wheat breeding and supporting future population genomic and genome-wide association studies in wheat.