The genetic basis of cytoplasmic male sterility and fertility restoration in wheat.

Hybrid wheat varieties give higher yields than conventional lines but are difficult to produce due to a lack of effective control of male fertility in breeding lines. One promising system involves the Rf1 and Rf3 genes that restore fertility of wheat plants carrying Triticum timopheevii-type cytoplasmic male sterility (T-CMS). Here, by genetic mapping and comparative sequence analyses, we identify Rf1 and Rf3 candidates that can restore normal pollen production in transgenic wheat plants carrying T-CMS. We show that Rf1 and Rf3 bind to the mitochondrial orf279 transcript and induce cleavage, preventing expression of the CMS trait. The identification of restorer genes in wheat is an important step towards the development of hybrid wheat varieties based on a CMS-Rf system. The characterisation of their mode of action brings insights into the molecular basis of CMS and fertility restoration in plants.

Genome-Wide Association Analysis of Grain Yield-Associated Traits in a Pan-European Barley Cultivar Collection.

A collection of 379 Hordeum vulgare cultivars, comprising all combinations of spring and winter growth habits with two and six row ear type, was screened by genome wide association analysis to discover alleles controlling traits related to grain yield. Genotypes were obtained at 6,810 segregating gene-based single nucleotide polymorphism (SNP) loci and corresponding field trial data were obtained for eight traits related to grain yield at four European sites in three countries over two growth years. The combined data were analyzed and statistically significant associations between the traits and regions of the barley genomes were obtained. Combining this information with the high resolution gene map for barley allowed the identification of candidate genes underlying all scored traits and superposition of this information with the known genomics of grain trait genes in rice resulted in the assignation of 13 putative barley genes controlling grain traits in European cultivated barley. Several of these genes are associated with grain traits in both winter and spring barley.

Loss of pollen-specific phospholipase NOT LIKE DAD triggers gynogenesis in maize.

Gynogenesis is an asexual mode of reproduction common to animals and plants, in which stimuli from the sperm cell trigger the development of the unfertilized egg cell into a haploid embryo. Fine mapping restricted a major maize QTL (quantitative trait locus) responsible for the aptitude of inducer lines to trigger gynogenesis to a zone containing a single gene NOT LIKE DAD (NLD) coding for a patatin-like phospholipase A. In all surveyed inducer lines, NLD carries a 4-bp insertion leading to a predicted truncated protein. This frameshift mutation is responsible for haploid induction because complementation with wild-type NLD abolishes the haploid induction capacity. Activity of the NLD promoter is restricted to mature pollen and pollen tube. The translational NLD::citrine fusion protein likely localizes to the sperm cell plasma membrane. In Arabidopsis roots, the truncated protein is no longer localized to the plasma membrane, contrary to the wild-type NLD protein. In conclusion, an intact pollen-specific phospholipase is required for successful sexual reproduction and its targeted disruption may allow establishing powerful haploid breeding tools in numerous crops.

Genomic Selection in Multi-environment Crop Trials.

Genomic selection in crop breeding introduces modeling challenges not found in animal studies. These include the need to accommodate replicate plants for each line, consider spatial variation in field trials, address line by environment interactions, and capture nonadditive effects. Here, we propose a flexible single-stage genomic selection approach that resolves these issues. Our linear mixed model incorporates spatial variation through environment-specific terms, and also randomization-based design terms. It considers marker, and marker by environment interactions using ridge regression best linear unbiased prediction to extend genomic selection to multiple environments. Since the approach uses the raw data from line replicates, the line genetic variation is partitioned into marker and nonmarker residual genetic variation (i.e., additive and nonadditive effects). This results in a more precise estimate of marker genetic effects. Using barley height data from trials, in 2 different years, of up to 477 cultivars, we demonstrate that our new genomic selection model improves predictions compared to current models. Analyzing single trials revealed improvements in predictive ability of up to 5.7%. For the multiple environment trial (MET) model, combining both year trials improved predictive ability up to 11.4% compared to a single environment analysis. Benefits were significant even when fewer markers were used. Compared to a single-year standard model run with 3490 markers, our partitioned MET model achieved the same predictive ability using between 500 and 1000 markers depending on the trial. Our approach can be used to increase accuracy and confidence in the selection of the best lines for breeding and/or, to reduce costs by using fewer markers.

The USDA barley core collection: genetic diversity, population structure, and potential for genome-wide association studies.

New sources of genetic diversity must be incorporated into plant breeding programs if they are to continue increasing grain yield and quality, and tolerance to abiotic and biotic stresses. Germplasm collections provide a source of genetic and phenotypic diversity, but characterization of these resources is required to increase their utility for breeding programs. We used a barley SNP iSelect platform with 7,842 SNPs to genotype 2,417 barley accessions sampled from the USDA National Small Grains Collection of 33,176 accessions. Most of the accessions in this core collection are categorized as landraces or cultivars/breeding lines and were obtained from more than 100 countries. Both STRUCTURE and principal component analysis identified five major subpopulations within the core collection, mainly differentiated by geographical origin and spike row number (an inflorescence architecture trait). Different patterns of linkage disequilibrium (LD) were found across the barley genome and many regions of high LD contained traits involved in domestication and breeding selection. The genotype data were used to define 'mini-core' sets of accessions capturing the majority of the allelic diversity present in the core collection. These 'mini-core' sets can be used for evaluating traits that are difficult or expensive to score. Genome-wide association studies (GWAS) of 'hull cover', 'spike row number', and 'heading date' demonstrate the utility of the core collection for locating genetic factors determining important phenotypes. The GWAS results were referenced to a new barley consensus map containing 5,665 SNPs. Our results demonstrate that GWAS and high-density SNP genotyping are effective tools for plant breeders interested in accessing genetic diversity in large germplasm collections.

Identification of crop cultivars with consistently high lignocellulosic sugar release requires the use of appropriate statistical design and modelling.

BACKGROUND: In this study, a multi-parent population of barley cultivars was grown in the field for two consecutive years and then straw saccharification (sugar release by enzymes) was subsequently analysed in the laboratory to identify the cultivars with the highest consistent sugar yield. This experiment was used to assess the benefit of accounting for both the multi-phase and multi-environment aspects of large-scale phenotyping experiments with field-grown germplasm through sound statistical design and analysis. RESULTS: Complementary designs at both the field and laboratory phases of the experiment ensured that non-genetic sources of variation could be separated from the genetic variation of cultivars, which was the main target of the study. The field phase included biological replication and plot randomisation. The laboratory phase employed re-randomisation and technical replication of samples within a batch, with a subset of cultivars chosen as duplicates that were randomly allocated across batches. The resulting data was analysed using a linear mixed model that incorporated field and laboratory variation and a cultivar by trial interaction, and ensured that the cultivar means were more accurately represented than if the non-genetic variation was ignored. The heritability detected was more than doubled in each year of the trial by accounting for the non-genetic variation in the analysis, clearly showing the benefit of this design and approach. CONCLUSIONS: The importance of accounting for both field and laboratory variation, as well as the cultivar by trial interaction, by fitting a single statistical model (multi-environment trial, MET, model), was evidenced by the changes in list of the top 40 cultivars showing the highest sugar yields. Failure to account for this interaction resulted in only eight cultivars that were consistently in the top 40 in different years. The correspondence between the rankings of cultivars was much higher at 25 in the MET model. This approach is suited to any multi-phase and multi-environment population-based genetic experiment.

Variation in the interaction between alleles of HvAPETALA2 and microRNA172 determines the density of grains on the barley inflorescence.

Within the cereal grasses, variation in inflorescence architecture results in a conspicuous morphological diversity that in crop species influences the yield of cereal grains. Although significant progress has been made in identifying some of the genes underlying this variation in maize and rice, in the temperate cereals, a group that includes wheat, barley, and rye, only the dosage-dependent and highly pleiotropic Q locus in hexaploid wheat has been molecularly characterized. Here we show that the characteristic variation in the density of grains along the inflorescence, or spike, of modern cultivated barley (Hordeum vulgare) is largely the consequence of a perturbed interaction between microRNA172 and its corresponding binding site in the mRNA of an APELATA2 (AP2)-like transcription factor, HvAP2. We used genome-wide association and biparental mapping to identify HvAP2. By comparing inflorescence development and HvAP2 transcript abundance in an extreme dense-spike mutant and its nearly isogenic WT line, we show that HvAP2 turnover driven by microRNA 172 regulates the length of a critical developmental window that is required for elongation of the inflorescence internodes. Our data indicate that this heterochronic change, an altered timing of developmental events caused by specific temporal variation in the efficiency of HvAP2 turnover, leads to the striking differences in the size and shape of the barley spike.

Genome-wide association mapping of frost tolerance in barley (Hordeum vulgare L.).

BACKGROUND: Frost tolerance is a key trait with economic and agronomic importance in barley because it is a major component of winter hardiness, and therefore limits the geographical distribution of the crop and the effective transfer of quality traits between spring and winter crop types. Three main frost tolerance QTL (Fr-H1, Fr-H2 and Fr-H3) have been identified from bi-parental genetic mapping but it can be argued that those mapping populations only capture a portion of the genetic diversity of the species. A genetically broad dataset consisting of 184 genotypes, representative of the barley gene pool cultivated in the Mediterranean basin over an extended time period, was genotyped with 1536 SNP markers. Frost tolerance phenotype scores were collected from two trial sites, Foradada (Spain) and Fiorenzuola (Italy) and combined with the genotypic data in genome wide association analyses (GWAS) using Eigenstrat and kinship approaches to account for population structure. RESULTS: GWAS analyses identified twelve and seven positive SNP associations at Foradada and Fiorenzuola, respectively, using Eigenstrat and six and four, respectively, using kinship. Linkage disequilibrium analyses of the significant SNP associations showed they are genetically independent. In the kinship analysis, two of the significant SNP associations were tightly linked to the Fr-H2 and HvBmy loci on chromosomes 5H and 4HL, respectively. The other significant kinship associations were located in genomic regions that have not previously been associated with cold stress. CONCLUSIONS: Haplotype analysis revealed that most of the significant SNP loci are fixed in the winter or facultative types, while they are freely segregating within the un-adapted spring barley genepool. Although there is a major interest in detecting new variation to improve frost tolerance of available winter and facultative types, from a GWAS perspective, working within the un-adapted spring germplasm pool is an attractive alternative strategy which would minimize statistical issues, simplify the interpretation of the data and identify phenology independent genetic determinants of frost tolerance.

Natural variation in a homolog of Antirrhinum CENTRORADIALIS contributed to spring growth habit and environmental adaptation in cultivated barley.

As early farming spread from the Fertile Crescent in the Near East around 10,000 years before the present, domesticated crops encountered considerable ecological and environmental change. Spring-sown crops that flowered without the need for an extended period of cold to promote flowering and day length-insensitive crops able to exploit the longer, cooler days of higher latitudes emerged and became established. To investigate the genetic consequences of adaptation to these new environments, we identified signatures of divergent selection in the highly differentiated modern-day spring and winter barleys. In one genetically divergent region, we identify a natural variant of the barley homolog of Antirrhinum CENTRORADIALIS (HvCEN) as a contributor to successful environmental adaptation. The distribution of HvCEN alleles in a large collection of wild and landrace accessions indicates that this involved selection and enrichment of preexisting genetic variants rather than the acquisition of mutations after domestication.

Tibet is one of the centers of domestication of cultivated barley.

The Near East Fertile Crescent is well recognized as a primary center of barley origin, diversity, and domestication. A large number of wild barleys have been collected from the Tibetan Plateau, which is characterized by an extreme environment. We used genome-wide diversity array technology markers to analyze the genotypic division between wild barley from the Near East and Tibet. Our results confirmed the existence of Tibetan wild barley and suggested that the split between the wild barleys in the Near East and those in Tibet occurred around 2.76 million years ago (Mya). To test the concept of polyphyletic domestication of barley, we characterized a set of worldwide cultivated barley. Some Chinese hulless and six-rowed barleys showed a close relationship with Tibetan wild barley but showed no common ancestor with other cultivated barley. Our data support the concept of polyphyletic domestication of cultivated barley and indicate that the Tibetan Plateau and its vicinity is one of the centers of domestication of cultivated barley. The current results may be highly significant in exploring the elite germplasm for barley breeding, especially against cold and drought stresses.

Genome-wide association mapping of agronomic and morphologic traits in highly structured populations of barley cultivars.

Genome-wide association study (GWAS) has become an obvious general approach for studying traits of agricultural importance in higher plants, especially crops. Here, we present a GWAS of 32 morphologic and 10 agronomic traits in a collection of 615 barley cultivars genotyped by genome-wide polymorphisms from a recently developed barley oligonucleotide pool assay. Strong population structure effect related to mixed sampling based on seasonal growth habit and ear row number is present in this barley collection. Comparison of seven statistical approaches in a genome-wide scan for significant associations with or without correction for confounding by population structure, revealed that in reducing false positive rates while maintaining statistical power, a mixed linear model solution outperforms genomic control, structured association, stepwise regression control and principal components adjustment. The present study reports significant associations for sixteen morphologic and nine agronomic traits and demonstrates the power and feasibility of applying GWAS to explore complex traits in highly structured plant samples.

INTERMEDIUM-C, a modifier of lateral spikelet fertility in barley, is an ortholog of the maize domestication gene TEOSINTE BRANCHED 1.

The domestication of cereals has involved common changes in morphological features, such as seed size, seed retention and modification of vegetative and inflorescence architecture that ultimately contributed to an increase in harvested yield. In barley, this process has resulted in two different cultivated types, two-rowed and six-rowed forms, both derived from the wild two-rowed ancestor, with archaeo-botanical evidence indicating the origin of six-rowed barley early in the domestication of the species, some 8,600-8,000 years ago. Variation at SIX-ROWED SPIKE 1 (VRS1) is sufficient to control this phenotype. However, phenotypes imposed by VRS1 alleles are modified by alleles at the INTERMEDIUM-C (INT-C) locus. Here we show that INT-C is an ortholog of the maize domestication gene TEOSINTE BRANCHED 1 (TB1) and identify 17 coding mutations in barley TB1 correlated with lateral spikelet fertility phenotypes.

Patterns of polymorphism and linkage disequilibrium in cultivated barley.

We carried out a genome-wide analysis of polymorphism (4,596 SNP loci across 190 elite cultivated accessions) chosen to represent the available genetic variation in current elite North West European and North American barley germplasm. Population sub-structure, patterns of diversity and linkage disequilibrium varied considerably across the seven barley chromosomes. Gene-rich and rarely recombining haplotype blocks that may represent up to 60% of the physical length of barley chromosomes extended across the 'genetic centromeres'. By positioning 2,132 bi-parentally mapped SNP markers with minimum allele frequencies higher than 0.10 by association mapping, 87.3% were located to within 5 cM of their original genetic map position. We show that at this current marker density genetically diverse populations of relatively small size are sufficient to fine map simple traits, providing they are not strongly stratified within the sample, fall outside the genetic centromeres and population sub-structure is effectively controlled in the analysis. Our results have important implications for association mapping, positional cloning, physical mapping and practical plant breeding in barley and other major world cereals including wheat and rye that exhibit comparable genome and genetic features.

Genome-wide association mapping to candidate polymorphism resolution in the unsequenced barley genome.

Although commonplace in human disease genetics, genome-wide association (GWA) studies have only relatively recently been applied to plants. Using 32 phenotypes in the inbreeding crop barley, we report GWA mapping of 15 morphological traits across ∼500 cultivars genotyped with 1,536 SNPs. In contrast to the majority of human GWA studies, we observe high levels of linkage disequilibrium within and between chromosomes. Despite this, GWA analysis readily detected common alleles of high penetrance. To investigate the potential of combining GWA mapping with comparative analysis to resolve traits to candidate polymorphism level in unsequenced genomes, we fine-mapped a selected phenotype (anthocyanin pigmentation) within a 140-kb interval containing three genes. Of these, resequencing the putative anthocyanin pathway gene HvbHLH1 identified a deletion resulting in a premature stop codon upstream of the basic helix-loop-helix domain, which was diagnostic for lack of anthocyanin in our association and biparental mapping populations. The methodology described here is transferable to species with limited genomic resources, providing a paradigm for reducing the threshold of map-based cloning in unsequenced crops.

Whole-genome association mapping in elite inbred crop varieties.

We have previously shown that linkage disequilibrium (LD) in the elite cultivated barley (Hordeum vulgare) gene pool extends, on average, for <1-5 cM. Based on this information, we have developed a platform for whole genome association studies that comprises a collection of elite lines that we have characterized at 3060 genome-wide single nucleotide polymorphism (SNP) marker loci. Interrogating this data set shows that significant population substructure is present within the elite gene pool and that diversity and LD vary considerably across each of the seven barley chromosomes. However, we also show that a subpopulation comprised of only the two-rowed spring germplasm is less structured and well suited to whole genome association studies without the need for extensive statistical intervention to account for structure. At the current marker density, the two-rowed spring population is suited for fine mapping simple traits that are located outside of the genetic centromeres with a resolution that is sufficient for candidate gene identification by exploiting conservation of synteny with fully sequenced model genomes and the emerging barley physical map.

Patterns of genetic diversity and linkage disequilibrium in a highly structured Hordeum vulgare association-mapping population for the Mediterranean basin.

Population structure and genome-wide linkage disequilibrium (LD) were investigated in 192 Hordeum vulgare accessions providing a comprehensive coverage of past and present barley breeding in the Mediterranean basin, using 50 nuclear microsatellite and 1,130 DArT((R)) markers. Both clustering and principal coordinate analyses clearly sub-divided the sample into five distinct groups centred on key ancestors and regions of origin of the germplasm. For given genetic distances, large variation in LD values was observed, ranging from closely linked markers completely at equilibrium to marker pairs at 50 cM separation still showing significant LD. Mean LD values across the whole population sample decayed below r (2) of 0.15 after 3.2 cM. By assaying 1,130 genome-wide DArT((R)) markers, we demonstrated that, after accounting for population substructure, current genome coverage of 1 marker per 1.5 cM except for chromosome 4H with 1 marker per 3.62 cM is sufficient for whole genome association scans. We show, by identifying associations with powdery mildew that map in genomic regions known to have resistance loci, that associations can be detected in strongly stratified samples provided population structure is effectively controlled in the analysis. The population we describe is, therefore, shown to be a valuable resource, which can be used in basic and applied research in barley.