Extreme population-dependent linkage disequilibrium detected in an inbreeding plant species, Hordeum vulgare

Farmers without borders-genetic structuring in century old barley (Hordeum vulgare)

The geographic distribution of genetic diversity can reveal the evolutionary history of a species. For crop plants, phylogeographic patterns also indicate how seed has been exchanged and spread in agrarian communities. Such patterns are, however, easily blurred by the intense seed trade, plant improvement and even genebank conservation during the twentieth century, and discerning fine-scale phylogeographic patterns is thus particularly challenging. Using historical crop specimens, these problems are circumvented and we show here how high-throughput genotyping of historical nineteenth century crop specimens can reveal detailed geographic population structure. Thirty-one historical and nine extant accessions of North European landrace barley (Hordeum vulgare L.), in total 231 individuals, were genotyped on a 384 single nucleotide polymorphism assay. The historical material shows constant high levels of within-accession diversity, whereas the extant accessions show more varying levels of diversity and a higher degree of total genotype sharing. Structure, discriminant analysis of principal components and principal component analysis cluster the accessions in latitudinal groups across country borders in Finland, Norway and Sweden. FST statistics indicate strong differentiation between accessions from southern Fennoscandia and accessions from central or northern Fennoscandia, and less differentiation between central and northern accessions. These findings are discussed in the context of contrasting historical records on intense within-country south to north seed movement. Our results suggest that although seeds were traded long distances, long-term cultivation has instead been of locally available, possibly better adapted, genotypes.

Plant systems biology: insights, advances and challenges

Plants dwelling at the base of biological food chain are of fundamental significance in providing solutions to some of the most daunting ecological and environmental problems faced by our planet. The reductionist views of molecular biology provide only a partial understanding to the phenotypic knowledge of plants. Systems biology offers a comprehensive view of plant systems, by employing a holistic approach integrating the molecular data at various hierarchical levels. In this review, we discuss the basics of systems biology including the various 'omics' approaches and their integration, the modeling aspects and the tools needed for the plant systems research. A particular emphasis is given to the recent analytical advances, updated published examples of plant systems biology studies and the future trends.

Linkage disequilibrium and genome-wide association mapping in tetraploid wheat (Triticum turgidum L.)

Association mapping is a powerful tool for the identification of quantitative trait loci through the exploitation of the differential decay of linkage disequilibrium (LD) between marker loci and genes of interest in natural and domesticated populations. Using a sample of 230 tetraploid wheat lines (Triticum turgidum ssp), which included naked and hulled accessions, we analysed the pattern of LD considering 26 simple sequence repeats and 970 mostly mapped diversity array technology loci. In addition, to validate the potential for association mapping in durum wheat, we evaluated the same genotypes for plant height, heading date, protein content, and thousand-kernel weight. Molecular and phenotypic data were used to: (i) investigate the genetic and phenotypic diversity; (ii) study the dynamics of LD across the durum wheat genome, by investigating the patterns of LD decay; and (iii) test the potential of our panel to identify marker–trait associations through the analysis of four quantitative traits of major agronomic importance. Moreover, we compared and validated the association mapping results with outlier detection analysis based on population divergence. Overall, in tetraploid wheat, the pattern of LD is extremely population dependent and is related to the domestication and breeding history of durum wheat. Comparing our data with several other studies in wheat, we confirm the position of many major genes and quantitative trait loci for the traits considered. Finally, the analysis of the selection signature represents a very useful complement to validate marker–trait associations.

Two genomic regions contribute disproportionately to geographic differentiation in wild barley

Genetic differentiation in natural populations is driven by geographic distance and by ecological or physical features within and between natural habitats that reduce migration. The primary population structure in wild barley differentiates populations east and west of the Zagros Mountains. Genetic differentiation between eastern and western populations is uneven across the genome and is greatest on linkage groups 2H and 5H. Genetic markers in these two regions demonstrate the largest difference in frequency between the primary populations and have the highest informativeness for assignment to each population. Previous cytological and genetic studies suggest there are chromosomal structural rearrangements (inversions or translocations) in these genomic regions. Environmental association analyses identified an association with both temperature and precipitation variables on 2H and with precipitation variables on 5H.

Two genomic regions contribute disproportionately to geographic differentiation in wild barley

Genetic differentiation in natural populations is driven by geographic distance and by ecological or physical features within and between natural habitats that reduce migration. The primary population structure in wild barley differentiates populations east and west of the Zagros Mountains. Genetic differentiation between eastern and western populations is uneven across the genome and is greatest on linkage groups 2H and 5H. Genetic markers in these two regions demonstrate the largest difference in frequency between the primary populations and have the highest informativeness for assignment to each population. Previous cytological and genetic studies suggest there are chromosomal structural rearrangements (inversions or translocations) in these genomic regions. Environmental association analyses identified an association with both temperature and precipitation variables on 2H and with precipitation variables on 5H.

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 agronomically important QTL in tetraploid potato cultivars using a marker-trait association analysis

Nineteen tuber quality traits in potato were phenotyped in 205 cultivars and 299 breeder clones. Association analysis using 3364 AFLP loci and 653 SSR-alleles identified QTL for these traits. Two association mapping panels were analysed for marker-trait associations to identify quantitative trait loci (QTL). The first panel comprised 205 historical and contemporary tetraploid potato cultivars that were phenotyped in field trials at two locations with two replicates (the academic panel). The second panel consisted of 299 potato cultivars and included recent breeds obtained from five Dutch potato breeding companies and reference cultivars (the industrial panel). Phenotypic data for the second panel were collected during subsequent clonal selection generations at the individual breeding companies. QTL were identified for 19 agro-morphological and quality traits. Two association mapping models were used: a baseline model without, and a more advanced model with correction for population structure and genetic relatedness. Correction for population structure and genetic relatedness was performed with a kinship matrix estimated from marker information. The detected QTL partly not only confirmed previous studies, e.g. for tuber shape and frying colour, but also new QTL were found like for after baking darkening and enzymatic browning. Pleiotropic effects could be discerned for several QTL.

Genetic diversity and ecological niche modelling of wild barley: refugia, large-scale post-LGM range expansion and limited mid-future climate threats?

Describing genetic diversity in wild barley (Hordeum vulgare ssp. spontaneum) in geographic and environmental space in the context of current, past and potential future climates is important for conservation and for breeding the domesticated crop (Hordeum vulgare ssp. vulgare). Spatial genetic diversity in wild barley was revealed by both nuclear- (2,505 SNP, 24 nSSR) and chloroplast-derived (5 cpSSR) markers in 256 widely-sampled geo-referenced accessions. Results were compared with MaxEnt-modelled geographic distributions under current, past (Last Glacial Maximum, LGM) and mid-term future (anthropogenic scenario A2, the 2080s) climates. Comparisons suggest large-scale post-LGM range expansion in Central Asia and relatively small, but statistically significant, reductions in range-wide genetic diversity under future climate. Our analyses support the utility of ecological niche modelling for locating genetic diversity hotspots and determine priority geographic areas for wild barley conservation under anthropogenic climate change. Similar research on other cereal crop progenitors could play an important role in tailoring conservation and crop improvement strategies to support future human food security.

Genetic structure, linkage disequilibrium and association mapping of Verticillium wilt resistance in elite cotton (Gossypium hirsutum L.) germplasm population

Understanding the population structure and linkage disequilibrium in an association panel can effectively avoid spurious associations and improve the accuracy in association mapping. In this study, one hundred and fifty eight elite cotton (Gossypium hirsutum L.) germplasm from all over the world, which were genotyped with 212 whole genome-wide marker loci and phenotyped with an disease nursery and greenhouse screening method, were assayed for population structure, linkage disequilibrium, and association mapping of Verticillium wilt resistance. A total of 480 alleles ranging from 2 to 4 per locus were identified from all collections. Model-based analysis identified two groups (G1 and G2) and seven subgroups (G1a–c, G2a–d), and differentiation analysis showed that subgroup having a single origin or pedigree was apt to differentiate with those having a mixed origin. Only 8.12% linked marker pairs showed significant LD (P<0.001) in this association panel. The LD level for linked markers is significantly higher than that for unlinked markers, suggesting that physical linkage strongly influences LD in this panel, and LD level was elevated when the panel was classified into groups and subgroups. The LD decay analysis for several chromosomes showed that different chromosomes showed a notable change in LD decay distances for the same gene pool. Based on the disease nursery and greenhouse environment, 42 marker loci associated with Verticillium wilt resistance were identified through association mapping, which widely were distributed among 15 chromosomes. Among which 10 marker loci were found to be consistent with previously identified QTLs and 32 were new unreported marker loci, and QTL clusters for Verticillium wilt resistanc on Chr.16 were also proved in our study, which was consistent with the strong linkage in this chromosome. Our results would contribute to association mapping and supply the marker candidates for marker-assisted selection of Verticillium wilt resistance in cotton.

Candidate gene association mapping of Sclerotinia stalk rot resistance in sunflower (Helianthus annuus L.) uncovers the importance of COI1 homologs

Functional markers for Sclerotinia basal stalk rot resistance in sunflower were obtained using gene-level information from the model species Arabidopsis thaliana. Sclerotinia stalk rot, caused by Sclerotinia sclerotiorum, is one of the most destructive diseases of sunflower (Helianthus annuus L.) worldwide. Markers for genes controlling resistance to S. sclerotiorum will enable efficient marker-assisted selection (MAS). We sequenced eight candidate genes homologous to Arabidopsis thaliana defense genes known to be associated with Sclerotinia disease resistance in a sunflower association mapping population evaluated for Sclerotinia stalk rot resistance. The total candidate gene sequence regions covered a concatenated length of 3,791 bp per individual. A total of 187 polymorphic sites were detected for all candidate gene sequences, 149 of which were single nucleotide polymorphisms (SNPs) and 38 were insertions/deletions. Eight SNPs in the coding regions led to changes in amino acid codons. Linkage disequilibrium decay throughout the candidate gene regions declined on average to an r (2) = 0.2 for genetic intervals of 120 bp, but extended up to 350 bp with r (2) = 0.1. A general linear model with modification to account for population structure was found the best fitting model for this population and was used for association mapping. Both HaCOI1-1 and HaCOI1-2 were found to be strongly associated with Sclerotinia stalk rot resistance and explained 7.4 % of phenotypic variation in this population. These SNP markers associated with Sclerotinia stalk rot resistance can potentially be applied to the selection of favorable genotypes, which will significantly improve the efficiency of MAS during the development of stalk rot resistant cultivars.

Candidate gene association mapping of Sclerotinia stalk rot resistance in sunflower (Helianthus annuus L.) uncovers the importance of COI1 homologs

Functional markers for Sclerotinia basal stalk rot resistance in sunflower were obtained using gene-level information from the model species Arabidopsis thaliana. Sclerotinia stalk rot, caused by Sclerotinia sclerotiorum, is one of the most destructive diseases of sunflower (Helianthus annuus L.) worldwide. Markers for genes controlling resistance to S. sclerotiorum will enable efficient marker-assisted selection (MAS). We sequenced eight candidate genes homologous to Arabidopsis thaliana defense genes known to be associated with Sclerotinia disease resistance in a sunflower association mapping population evaluated for Sclerotinia stalk rot resistance. The total candidate gene sequence regions covered a concatenated length of 3,791 bp per individual. A total of 187 polymorphic sites were detected for all candidate gene sequences, 149 of which were single nucleotide polymorphisms (SNPs) and 38 were insertions/deletions. Eight SNPs in the coding regions led to changes in amino acid codons. Linkage disequilibrium decay throughout the candidate gene regions declined on average to an r (2) = 0.2 for genetic intervals of 120 bp, but extended up to 350 bp with r (2) = 0.1. A general linear model with modification to account for population structure was found the best fitting model for this population and was used for association mapping. Both HaCOI1-1 and HaCOI1-2 were found to be strongly associated with Sclerotinia stalk rot resistance and explained 7.4 % of phenotypic variation in this population. These SNP markers associated with Sclerotinia stalk rot resistance can potentially be applied to the selection of favorable genotypes, which will significantly improve the efficiency of MAS during the development of stalk rot resistant cultivars.

Population structure of barley landrace populations and gene-flow with modern varieties

Landraces are heterogeneous plant varieties that are reproduced by farmers as populations that are subject to both artificial and natural selection. Landraces are distinguished by farmers due to their specific traits, and different farmers often grow different populations of the same landrace. We used simple sequence repeats (SSRs) to analyse 12 barley landrace populations from Sardinia from two collections spanning 10 years. We analysed the population structure, and compared the population diversity of the landraces that were collected at field level (population). We used a representative pool of barley varieties for diversity comparisons and to analyse the effects of gene flow from modern varieties. We found that the Sardinian landraces are a distinct gene pool from those of both two-row and six-row barley varieties. There is also a low, but significant, mean level and population-dependent level of introgression from the modern varieties into the Sardinian landraces. Moreover, we show that the Sardinian landraces have the same level of gene diversity as the representative sample of modern commercial varieties grown in Italy in the last decades, even within population level. Thus, these populations represent crucial sources of germplasm that will be useful for crop improvement and for population genomics studies and association mapping, to identify genes, loci and genome regions responsible for adaptive variations. Our data also suggest that landraces are a source of valuable germplasm for sustainable agriculture in the context of future climate change, and that in-situ conservation strategies based on farmer use can preserve the genetic identity of landraces while allowing adaptation to local environments.

Resequencing data indicate a modest effect of domestication on diversity in barley: a cultigen with multiple origins

The levels of diversity and extent of linkage disequilibrium in cultivated species are largely determined by diversity in their wild progenitors. We report a comparison of nucleotide sequence diversity in wild and cultivated barley (Hordeum vulgare ssp. spontaneum and ssp. vulgare) at 7 nuclear loci totaling 9296bp, using sequence from Hordeum bulbosum to infer the ancestral state of mutations. The sample includes 36 accessions of cultivated barley, including 23 landraces (cultivated forms not subject to modern breeding) and 13 cultivated lines and genetic stocks compared to either 25 or 45 accessions of wild barley for the same loci. Estimates of nucleotide sequence diversity indicate that landraces retain >80% of the diversity in wild barley. The primary population structure in wild barley, which divides the species into eastern and western populations, is reflected in significant differentiation at all loci in wild accessions and at 3 of 7 loci in landraces. "Oriental" landraces have slightly higher diversity than "Occidental" landraces. Genetic assignment suggests more admixture from Occidental landraces into Oriental landraces than the converse, which may explain this difference. Based on θπ for silent sites, modern western cultivars have ~73% of the diversity found in landraces and ~71% of the diversity in wild barley.

Identification and distribution of Puroindoline b-2 variant gene homologs in Hordeum

The barley hordoindoline genes (Hina and Hinb) are homologous to the wheat puroindoline genes (Pina and Pinb). These genes are involved in grain hardness, which is an important quality for barley processing. We identified novel variants of Hina and Hinb in 10 wild Hordeum species (H. bogdanii, H. brachyantherum, H. bulbosum, H. chilense, H. comosum, H. marinum, H. murinum, H. patagonicum, H. pusillum, and H. roshevitzii) covering all Hordeum genomes and preliminarily named them Hinc. These nucleotide sequences were highly similar to those of Puroindoline b-2 variant genes (Pinb-2v) and were located on chromosome 7I in H. chilense. The Hinc genes in H. bogdanii, H. bulbosum, H. patagonicum, and H. roshevitzii were pseudogenes possessing in-frame stop codons. We also found a partial Hinc sequence in H. murinum. This gene was not found in cultivated barley and H. vulgare subsp. spontaneum. The phylogenetic tree of Gsp-1, Hin, and Pin genes demonstrates that Hinc and Pinb-2v genes formed one cluster. Therefore, we considered that Hinc and Pinb-2v genes shared a common ancestral gene and were homologous to each other. We also studied the evolutional process of Gsp-1, Hin, and Pin genes. Our results suggested that Gsp-1 might be the most closely related to a putative ancestral gene on Ha locus.

Association mapping of maturity and plant height using SNP markers with the sorghum mini core collection

Plant height and maturity are two critical traits in sorghum breeding. To develop molecular tools and to identify genes underlying the traits for molecular breeding, we developed 14,739 SNP markers used to genotype the complete sorghum [Sorghum bicolor (L.) Moench] mini core collection. The collection was evaluated in four rainy and three post-rainy season environments for plant height and maturity. Association analysis identified six marker loci linked to height and ten to maturity in at least two environments with at least two SNPs in each locus. Of these, 14 were in close proximity to previously mapped height/maturity QTL in sorghum. Candidate genes for maturity or plant height close to the marker loci include a sugar transporter (SbSUC9), an auxin response factor (SbARF3), an FLC and FT regulator (SbMED12), and a photoperiod response gene (SbPPR1) for maturity and peroxidase 53, and an auxin transporter (SbLAX4) for plant height. Linkage disequilibrium analysis showed that SbPPR1 and SbARF3 were in regions with reduced sequence variation among early-maturing accessions, suggestive of past purifying selection. We also found a linkage disequilibrium block that existed only among the accessions with short plant height in rainy season environments. The block contains a gene homologous to the Arabidopsis flowering time gene, LUMINIDEPENDENS (LD). Functional LD promotes early maturity while mutation delays maturity, affecting plant height. Previous studies also found reduced sequence variations within this gene. These newly-mapped SNP markers will facilitate further efforts to identify plant height or maturity genes in sorghum.

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.

High-density SNP-based genetic map development and linkage disequilibrium assessment in Brassica napus L

Background

High density genetic maps built with SNP markers that are polymorphic in various genetic backgrounds are very useful for studying the genetics of agronomical traits as well as genome organization and evolution. Simultaneous dense SNP genotyping of segregating populations and variety collections was applied to oilseed rape (Brassica napus L.) to obtain a high density genetic map for this species and to study the linkage disequilibrium pattern.

Results

We developed an integrated genetic map for oilseed rape by high throughput SNP genotyping of four segregating doubled haploid populations. A very high level of collinearity was observed between the four individual maps and a large number of markers (>59%) was common to more than two maps. The precise integrated map comprises 5764 SNP and 1603 PCR markers. With a total genetic length of 2250 cM, the integrated map contains a density of 3.27 markers (2.56 SNP) per cM. Genotyping of these mapped SNP markers in oilseed rape collections allowed polymorphism level and linkage disequilibrium (LD) to be studied across the different collections (winter vs spring, different seed quality types) and along the linkage groups. Overall, polymorphism level was higher and LD decayed faster in spring than in “00” winter oilseed rape types but this was shown to vary greatly along the linkage groups.

Conclusions

Our study provides a valuable resource for further genetic studies using linkage or association mapping, for marker assisted breeding and for Brassica napus sequence assembly and genome organization analyses.

Comparative analysis of genome composition in Triticeae reveals strong variation in transposable element dynamics and nucleotide diversity

A 454 sequencing snapshot was utilised to investigate the genome composition and nucleotide diversity of transposable elements (TEs) for several Triticeae taxa, including Triticum aestivum, Hordeum vulgare, Hordeum spontaneum and Secale cereale together with relatives of the A, B and D genome donors of wheat, Triticum urartu (A), Aegilops speltoides (S) and Aegilops tauschii (D). Additional taxa containing the A genome, Triticum monococcum and its wild relative Triticum boeoticum, were also included. The main focus of the analysis was on the genomic composition of TEs as these make up at least 80% of the overall genome content. Although more than 200 TE families were identified in each species, approximately 50% of the overall genome comprised 12-15 TE families. The BARE1 element was the largest contributor to all genomes, contributing more than 10% to the overall genome. We also found that several TE families differ strongly in their abundance between species, indicating that TE families can thrive extremely successfully in one species while going virtually extinct in another. Additionally, the nucleotide diversity of BARE1 populations within individual genomes was measured. Interestingly, the nucleotide diversity in the domesticated barley H. vulgare cv. Barke was found to be twice as high as in its wild progenitor H. spontaneum, suggesting that the domesticated barley gained nucleotide diversity from the addition of different genotypes during the domestication and breeding process. In the rye/wheat lineage, sequence diversity of BARE1 elements was generally higher, suggesting that factors such as geographical distribution and mating systems might play a role in intragenomic TE diversity.

Molecular analysis of the parallel domestication of the common bean (Phaseolus vulgaris) in Mesoamerica and the Andes

We have studied the nucleotide diversity of common bean, Phaseolus vulgaris, which is characterized by two independent domestications in two geographically distinct areas: Mesoamerica and the Andes. This provides an important model, as domestication can be studied as a replicate experiment. We used nucleotide data from five gene fragments characterized by large introns to analyse 214 accessions (102 wild and 112 domesticated). The wild accessions represent a cross-section of the entire geographical distribution of P. vulgaris. A reduction in genetic diversity in both of these gene pools was found, which was three-fold greater in Mesoamerica compared with the Andes. This appears to be a result of a bottleneck that occurred before domestication in the Andes, which strongly impoverished this wild germplasm, leading to the minor effect of the subsequent domestication bottleneck (i.e. sequential bottleneck). These findings show the importance of considering the evolutionary history of crop species as a major factor that influences their current level and structure of genetic diversity. Furthermore, these data highlight a single domestication event within each gene pool. Although the findings should be interpreted with caution, this evidence indicates the Oaxaca valley in Mesoamerica, and southern Bolivia and northern Argentina in South America, as the origins of common bean domestication.

Nucleotide diversity and linkage disequilibrium in antioxidant genes of Brachypodium distachyon

Brachypodium distachyon (Brachypodium) is a powerful model system for studying cereal, bioenergy, forage, and turf grasses. Nucleotide diversity (π) and linkage disequilibrium (LD) in candidate genes involved in the antioxidative pathways in this species are not known. The average π for CAT encoding catalase, GPX encoding glutathione peroxidase, DHAR encoding dehydroascorbate reductase, MDHAR encoding monodehydroascorbate reductase, and APX ecoding ascorbate peroxidase was 0.0027 among 19 accessions contrasting for drought tolerance. The highest value of π was found in APX (0.0046) and the lowest π was in MDHAR (0.0006). The average single nucleotide polymorphism (SNP) frequency across these five genes was one SNP per 131 bp between two randomly sampled sequences for the five genes in the sequence length ranging from 1,447 bp to 1,701 bp. The LD decay was slow and extended to a distance of more than 1.2kb for all genes. The neighbor-joining tree analyses of DHAR, MDHAR, and CAT generally separated accessions differing in drought tolerance. The results indicate a putative role of these candidate genes in increasing general fitness of Brachypodium.

Phenotypic landscapes: phenological patterns in wild and cultivated barley

Phenotypic variation in natural populations is the outcome of the joint effects of environmentally induced adaptations and neutral processes on the genetic architecture of quantitative traits. In this study, we examined the role of adaptation in shaping wild barley phenotypic variation along different environmental gradients. Detailed phenotyping of 164 wild barley (Hordeum spontaneum) accessions from Israel (of the Barley1K collection) and 18 cultivated barley (H. vulgare) varieties was conducted in common garden field trials. Cluster analysis based on phenotypic data indicated that wild barley in this region can be differentiated into three ecotypes in accordance with their ecogeographical distribution: north, coast and desert. Population differentiation (Q(ST) ) for each trait was estimated using a hierarchical Bayesian model and compared to neutral differentiation (F(ST) ) based on 42 microsatellite markers. This analysis indicated that the three clusters diverged in morphological but not in reproductive characteristics. To address the issue of phenotypic variation along environmental gradients, climatic and soil gradients were compared with each of the measured traits given the geographical distance between sampling sites using a partial Mantel test. Flowering time and plant growth were found to be differentially correlated with climatic and soil characteristic gradients, respectively. The H. vulgare varieties were superior to the H. spontaneum accessions in yield components, yet resembled the Mediterranean types in vegetative characteristics and flowering time, which may indicate the geographical origin of domesticated barley.

Patterns of linkage disequilibrium and association mapping in diploid alfalfa (M. sativa L.)

Association mapping enables the detection of marker-trait associations in unstructured populations by taking advantage of historical linkage disequilibrium (LD) that exists between a marker and the true causative polymorphism of the trait phenotype. Our first objective was to understand the pattern of LD decay in the diploid alfalfa genome. We used 89 highly polymorphic SSR loci in 374 unimproved diploid alfalfa (Medicago sativa L.) genotypes from 120 accessions to infer chromosome-wide patterns of LD. We also sequenced four lignin biosynthesis candidate genes (caffeoyl-CoA 3-O-methyltransferase (CCoAoMT), ferulate-5-hydroxylase (F5H), caffeic acid-O-methyltransferase (COMT), and phenylalanine amonialyase (PAL 1)) to identify single nucleotide polymorphisms (SNPs) and infer within gene estimates of LD. As the second objective of this study, we conducted association mapping for cell wall components and agronomic traits using the SSR markers and SNPs from the four candidate genes. We found very little LD among SSR markers implying limited value for genomewide association studies. In contrast, within gene LD decayed within 300 bp below an r (2) of 0.2 in three of four candidate genes. We identified one SSR and two highly significant SNPs associated with biomass yield. Based on our results, focusing association mapping on candidate gene sequences will be necessary until a dense set of genome-wide markers is available for alfalfa.

Genome-wide association mapping in tomato (Solanum lycopersicum) is possible using genome admixture of Solanum lycopersicum var. cerasiforme

Genome-wide association mapping is an efficient way to identify quantitative trait loci controlling the variation of phenotypes, but the approach suffers severe limitations when one is studying inbred crops like cultivated tomato (Solanum lycopersicum). Such crops exhibit low rates of molecular polymorphism and high linkage disequilibrium, which reduces mapping resolution. The cherry type tomato (S. lycopersicum var. cerasiforme) genome has been described as an admixture between the cultivated tomato and its wild ancestor, S. pimpinellifolium. We have thus taken advantage of the properties of this admixture to improve the resolution of association mapping in tomato. As a proof of concept, we sequenced 81 DNA fragments distributed on chromosome 2 at different distances in a core collection of 90 tomato accessions, including mostly cherry type tomato accessions. The 81 Sequence Tag Sites revealed 352 SNPs and indels. Molecular diversity was greatest for S. pimpinellifolium accessions, intermediate for S. l. cerasiforme accessions, and lowest for the cultivated group. We assessed the structure of molecular polymorphism and the extent of linkage disequilibrium over genetic and physical distances. Linkage disequilibrium decreased under r² = 0.3 within 1 cM, and minimal estimated value (r² = 0.13) was reached within 20 kb over the physical regions studied. Associations between polymorphisms and fruit weight, locule number, and soluble solid content were detected. Several candidate genes and quantitative trait loci previously identified were validated and new associations detected. This study shows the advantages of using a collection of S. l. cerasiforme accessions to overcome the low resolution of association mapping in tomato.

Genetic structure and linkage disequilibrium in landrace populations of barley in Sardinia

Multilocus digenic linkage disequilibria (LD) and their population structure were investigated in eleven landrace populations of barley (Hordeum vulgare ssp. vulgare L.) in Sardinia, using 134 dominant simple-sequence amplified polymorphism markers. The analysis of molecular variance for these markers indicated that the populations were partially differentiated (F(ST) = 0.18), and clustered into three geographic areas. Consistent with this population pattern, STRUCTURE analysis allocated individuals from a bulk of all populations into four genetic groups, and these groups also showed geographic patterns. In agreement with other molecular studies in barley, the general level of LD was low (13% of locus pairs, with P < 0.01) in the bulk of 337 lines, and decayed steeply with map distance between markers. The partitioning of multilocus associations into various components indicated that genetic drift and founder effects played a major role in determining the overall genetic makeup of the diversity in these landrace populations, but that epistatic homogenising or diversifying selection was also present. Notably, the variance of the disequilibrium component was relatively high, which implies caution in the pooling of barley lines for association studies. Finally, we compared the analyses of multilocus structure in barley landrace populations with parallel analyses in both composite crosses of barley on the one hand and in natural populations of wild barley on the other. Neither of these serves as suitable mimics of landraces in barley, which require their own study. Overall, the results suggest that these populations can be exploited for LD mapping if population structure is controlled.

Evolution of neutral and flowering genes along pearl millet (Pennisetum glaucum) domestication

Background

Pearl millet landraces display an important variation in their cycle duration. This diversity contributes to the stability of crop production in the Sahel despite inter-annual rainfall fluctuation. Conservation of phenological diversity is important for the future of pearl millet improvement and sustainable use. Identification of genes contributing to flowering time variation is therefore relevant. In this study we focused on three flowering candidate genes, PgHd3a, PgDwarf8 and PgPHYC. We tested for signatures of past selective events within polymorphism patterns of these three genes that could have been associated with pearl millet domestication and/or landraces differentiation. In order to implement ad hoc neutrality tests, a plausible demographic history of pearl millet domestication was inferred through Approximate Bayesian Computation by using eight neutral STS loci.

Results

Domesticated pearl millet exhibited 84% of the nucleotide diversity level found in the wild population. No specific polymorphisms were found either in the wild or in the domestic populations. The Bayesian approach and previous studies suggest that gene flow between wild relatives and domesticated pearl millets is a main factor explaining these results. Early and late landraces did not show significant genetic differentiation at both the neutral and the candidate loci. A positive selection was evidenced in PgHd3a and PgDwarf8 genes of domestic forms but not in the wild population.

Conclusion

Our results strongly suggest that PgHd3a and PgDwarf8 were likely targeted by selection during domestication. However, a potential role of any of the three candidate genes in the phenological differentiation between early and late landraces was not supported by our data. Reasons why these results contrast with previous results that have shown a slight but significant association between PgPHYC polymorphisms and variation in flowering time in pearl millet are discussed.

Genome-wide association studies for agronomical traits in a world wide spring barley collection

BACKGROUND

Genome-wide association studies (GWAS) based on linkage disequilibrium (LD) provide a promising tool for the detection and fine mapping of quantitative trait loci (QTL) underlying complex agronomic traits. In this study we explored the genetic basis of variation for the traits heading date, plant height, thousand grain weight, starch content and crude protein content in a diverse collection of 224 spring barleys of worldwide origin. The whole panel was genotyped with a customized oligonucleotide pool assay containing 1536 SNPs using Illumina's GoldenGate technology resulting in 957 successful SNPs covering all chromosomes. The morphological trait "row type" (two-rowed spike vs. six-rowed spike) was used to confirm the high level of selectivity and sensitivity of the approach. This study describes the detection of QTL for the above mentioned agronomic traits by GWAS.

RESULTS

Population structure in the panel was investigated by various methods and six subgroups that are mainly based on their spike morphology and region of origin. We explored the patterns of linkage disequilibrium (LD) among the whole panel for all seven barley chromosomes. Average LD was observed to decay below a critical level (r2-value 0.2) within a map distance of 5-10 cM. Phenotypic variation within the panel was reasonably large for all the traits. The heritabilities calculated for each trait over multi-environment experiments ranged between 0.90-0.95. Different statistical models were tested to control spurious LD caused by population structure and to calculate the P-value of marker-trait associations. Using a mixed linear model with kinship for controlling spurious LD effects, we found a total of 171 significant marker trait associations, which delineate into 107 QTL regions. Across all traits these can be grouped into 57 novel QTL and 50 QTL that are congruent with previously mapped QTL positions.

CONCLUSIONS

Our results demonstrate that the described diverse barley panel can be efficiently used for GWAS of various quantitative traits, provided that population structure is appropriately taken into account. The observed significant marker trait associations provide a refined insight into the genetic architecture of important agronomic traits in barley. However, individual QTL account only for a small portion of phenotypic variation, which may be due to insufficient marker coverage and/or the elimination of rare alleles prior to analysis. The fact that the combined SNP effects fall short of explaining the complete phenotypic variance may support the hypothesis that the expression of a quantitative trait is caused by a large number of very small effects that escape detection. Notwithstanding these limitations, the integration of GWAS with biparental linkage mapping and an ever increasing body of genomic sequence information will facilitate the systematic isolation of agronomically important genes and subsequent analysis of their allelic diversity.

Distribution of Hordoindoline genes in the genus Hordeum

Hordoindoline (Hin) genes, which are known to comprise Hina, Hinb-1, and Hinb-2, are associated with grain hardness in barley. However, the interspecific variation in the Hin genes in the genus Hordeum has not been studied in detail. We examined the variation in Hin genes and used it to infer the phylogenetic relationships between the genes found in two H. vulgare subspecies (cultivated barley and H. vulgare subsp. spontaneum) and 10 wild relatives (H. bogdanii, H. brachyantherum, H. bulbosum, H. chilense, H. comosum, H. marinum, H. murinum, H. patagonicum, H. pusillum, and H. roshevitzii). The Hina and Hinb genes of these species were amplified by PCR. We found two Hinb genes in three wild species (H. bogdanii, H. brachyantherum, and H. roshevitzii) and preliminarily named them Hinb-A and Hinb-B. Cluster analysis showed that the 17 Hinb genes present in Hordeum formed two distinct clusters (named A and B). Seven Hinb genes were included in Cluster-A, and 10 Hinb genes were included in Cluster-B. All Hinb-A genes were included in Cluster-A, while all of the Hinb-B genes were included in Cluster-B. In contrast, the Hinb-1 and Hinb-2 genes in H. vulgare were included in Cluster-B. These results suggest that the Hinb genes duplicated during the early stages of diversification in the genus Hordeum. On the other hand, the Hinb-1 and Hinb-2 genes in H. vulgare seem to have been generated by a duplication of the Hinb gene after the split of the lineages leading to H. vulgare and H. bulbosum.

Investigation of the domestication of common bean (Phaseolus vulgaris) using multilocus sequence data

Multilocus sequence data collected from domesticated and related wild relatives provides a rich source of information on the effect of human selection on the diversity and adaptability of a species to complex environments. To evaluate the domestication history of common bean (Phaseolus vulgaris L.), multilocus sequence data from landraces representing the various races within the Middle American (MA) and Andean gene pools was evaluated. Across 13 loci, nucleotide diversity was similar between landraces and wild germplasm in both gene pools. The diversity data were evaluated using the approximate Bayesian computation approach to test multiple domestication models and estimate population demographic parameters. A model with a single domestication event coupled with bidirectional migration between wild and domesticated genotypes fitted the data better than models consisting of two or three domestication events in each genepool. The effective bottleneck population size was ~50% of the base population in each genepool. The bottleneck began ~8200 and ~8500 years before present and ended at ~6300 and ~7000 years before present in MA and Andean gene pools respectively. Linkage disequilibrium decayed to a greater extent in the MA genepool. Given the (1) geographical adaptation bottleneck in each wild gene pool, (2) a subsequent domestication bottleneck within each gene pool, (3) differentiation into gene-pool specific races and (4) variable extents of linkage disequilibrium, association mapping experiments for common bean would more appropriately be performed within each genepool.

Nucleotide diversity of a genomic sequence similar to SHATTERPROOF (PvSHP1) in domesticated and wild common bean (Phaseolus vulgaris L.)

Evolutionary studies in plant and animal breeding are aimed at understanding the structure and organization of genetic variations of species. We have identified and characterized a genomic sequence in Phaseolus vulgaris of 1,200 bp (PvSHP1) that is homologous to SHATTERPROOF-1 (SHP1), a gene involved in control of fruit shattering in Arabidopsis thaliana. The PvSHP1 fragment was mapped to chromosome Pv06 in P. vulgaris and is linked to the flower and seed color gene V. Amplification of the PvSHP1 sequence from the most agronomically important legume species showed a high degree of interspecies diversity in the introns within the Phaseoleae, while the coding region was conserved across distant taxa. Sequencing of the PvSHP1 sequence in a sample of 91 wild and domesticated genotypes that span the geographic distribution of this species in the centers of origin showed that PvSHP1 is highly polymorphic and, therefore, particularly useful to further investigate the origin and domestication history of P. vulgaris. Our data confirm the gene pool structure seen in P. vulgaris along with independent domestication processes in the Andes and Mesoamerica; they provide additional evidence for a single domestication event in Mesoamerica. Moreover, our results support the Mesoamerican origin of this species. Finally, we have developed three indel-spanning markers that will be very useful for bean germplasm characterization, and particularly to trace the distribution of the domesticated Andean and Mesoamerican gene pools.

Chromosome 5H of Hordeum species involved in reduction in grain hardness in wheat genetic background

Grain hardness is an important factor affecting end-use quality in wheat. Mutations of the puroindoline genes, which are located on chromosome 5DS, control a majority of grain texture variations. Hordoindoline genes, which are the puroindoline gene homologs in barley, are located on chromosome 5HS and are also responsible for grain texture variation. In this study, we used three types of wheat-barley species (Hordeum vulgare, H. vulgare ssp. spontaneum, and H. chilense) chromosome addition lines and studied the effect of chromosome 5H of these species on wheat grain characteristics. The 5H chromosome addition lines showed significantly lower grain hardness and higher grain weight than the corresponding wheat parents. The effect of enhancing grain softness was largest in the wheat-H. chilense line regardless of having an increase in grain weight similar to those in the wheat-H. vulgare and wheat-H. spontaneum lines. Our results indicated that chromosome 5H of the Hordeum species plays a role in enhancing grain softness and increasing grain weight in the wheat genetic background, and the extent of effect on grain hardness depends on the type of Hordeum species. Protein analysis of hordoindolines indicated that profiles of 2D-electrophoresis of hordoindolines were different among Hordeum species and hordoindolines in the addition lines appeared to be most abundant in wheat-H. chilense line. The differences in enhancing grain softness among the Hordeum species might be attributed to the quantity of hordoindolines expressed in the 5H chromosome addition lines. These results suggested that the barley hordoindolines located on chromosome 5HS play a role in reducing grain hardness in the wheat genetic background.

The heterogeneous levels of linkage disequilibrium in white spruce genes and comparative analysis with other conifers

In plants, knowledge about linkage disequilibrium (LD) is relevant for the design of efficient single-nucleotide polymorphism arrays in relation to their use in population and association genomics studies. Previous studies of conifer genes have shown LD to decay rapidly within gene limits, but exceptions have been reported. To evaluate the extent of heterogeneity of LD among conifer genes and its potential causes, we examined LD in 105 genes of white spruce (Picea glauca) by sequencing a panel of 48 haploid megagametophytes from natural populations and further compared it with LD in other conifer species. The average pairwise r(2) value was 0.19 (s.d.=0.19), and LD dropped quickly with a half-decay being reached at a distance of 65 nucleotides between sites. However, LD was significantly heterogeneous among genes. A first group of 29 genes had stronger LD (mean r(2)=0.28), and a second group of 38 genes had weaker LD (mean r(2)=0.12). While a strong relationship was found with the recombination rate, there was no obvious relationship between LD and functional classification. The level of nucleotide diversity, which was highly heterogeneous across genes, was also not significantly correlated with LD. A search for selection signatures highlighted significant deviations from the standard neutral model, which could be mostly attributed to recent demographic changes. Little evidence was seen for hitchhiking and clear relationships with LD. When compared among conifer species, on average, levels of LD were similar in genes from white spruce, Norway spruce and Scots pine, whereas loblolly pine and Douglas fir genes exhibited a significantly higher LD.

From forest to field: perennial fruit crop domestication

PREMISE OF THE STUDY

Archaeological and genetic analyses of seed-propagated annual crops have greatly advanced our understanding of plant domestication and evolution. Comparatively little is known about perennial plant domestication, a relevant topic for understanding how genes and genomes evolve in long-lived species, and how perennials respond to selection pressures operating on a relatively short time scale. Here, we focus on long-lived perennial crops (mainly trees and other woody plants) grown for their fruits.

KEY RESULTS

We reviewed (1) the basic biology of long-lived perennials, setting the stage for perennial domestication by considering how these species evolve in nature; (2) the suite of morphological features associated with perennial fruit crops undergoing domestication; (3) the origins and evolution of domesticated perennials grown for their fruits; and (4) the genetic basis of domestication in perennial fruit crops.

CONCLUSIONS

Long-lived perennials have lengthy juvenile phases, extensive outcrossing, widespread hybridization, and limited population structure. Under domestication, these features, combined with clonal propagation, multiple origins, and ongoing crop-wild gene flow, contribute to mild domestication bottlenecks in perennial fruit crops. Morphological changes under domestication have many parallels to annual crops, but with key differences for mating system evolution and mode of reproduction. Quantitative trait loci associated with domestication traits in perennials are mainly of minor effect and may not be stable across years. Future studies that take advantage of genomic approaches and consider demographic history will elucidate the genetics of agriculturally and ecologically important traits in perennial fruit crops and their wild relatives.

Genetic variation of HvCBF genes and their association with salinity tolerance in Tibetan annual wild barley

The evaluation of both the genetic variation and the identification of salinity tolerant accessions of Tibetan annual wild barley (hereafter referred to as Tibetan barley) (Hordeum vulgare L. ssp. Spontaneum and H. vulgare L. ssp. agriocrithum) are essential for discovering and exploiting novel alleles involved in salinity tolerance. In this study, we examined tissue dry biomass and the Na⁺ and K⁺ contents of 188 Tibetan barley accessions in response to salt stress. We investigated the genetic variation of transcription factors HvCBF1, HvCBF3 and HvCBF4 within these accessions, conducting association analysis between these three genes and the respective genotypic salt tolerance. Salt stress significantly reduced shoot and root dry weight by 27.6% to 73.1% in the Tibetan barley lines. HvCBF1, HvCBF3 and HvCBF4 showed diverse sequence variation in amplicon as evident by the identification of single nucleotide polymorphisms (SNPs) and 3, 8 and 13 haplotypes, respectively. Furthermore, the decay of Linkage disequilibrium (LD) of chromosome 5 was 8.9 cM (r²<0.1). Marker bpb-4891 and haplotype 13 (Ps 610) of the HvCBF4 gene were significantly (P<0.05) and highly significantly (P<0.001) associated with salt tolerance. However, HvCBF1 and HvCBF3 genes were not associated with salinity tolerance. The accessions from haplotype 13 of the HvCBF4 gene showed high salinity tolerance, maintaining significantly lower Na⁺/K⁺ ratios and higher dry weight. It is thus proposed that these Tibetan barley accessions could be of value for enhancing salinity tolerance in cultivated barley.

Analysis of >1000 single nucleotide polymorphisms in geographically matched samples of landrace and wild barley indicates secondary contact and chromosome-level differences in diversity around domestication genes

Barley is a model species for the investigation of the evolution, adaptation and spread of the world's important crops. In this article, we describe the first application of an oligonucleotide pool assay single nucleotide polymorphism (SNP) platform to assess the evolution of barley in a portion of the Fertile Crescent, a key region in the development of farming. A large collection of >1000 genetically mapped, genome-wide SNPs was assayed in geographically matched landrace and wild barley accessions (N=448) from Jordan and Syria. Landrace and wild barley categories were clearly genetically differentiated, but a limited degree of secondary contact was evident. Significant chromosome-level differences in diversity between barley types were observed around genes known to be involved in the evolution of cultivars. The region of Jordan and southern Syria, compared with the north of Syria, was supported by SNP data as a more likely domestication origin. Our data provide evidence for hybridization as a possible mechanism for the continued adaptation of landrace barley under cultivation, indicate regions of the genome that may be subject to selection processes and suggest limited origins for the development of the cultivated crop.

Genetic analysis of abiotic stress tolerance in crops

Abiotic stress tolerance is complex, but as phenotyping technologies improve, components that contribute to abiotic stress tolerance can be quantified with increasing ease. In parallel with these phenomics advances, genetic approaches with more complex genomes are becoming increasingly tractable as genomic information in non-model crops increases and even whole crop genomes can be re-sequenced. Thus, genetic approaches to elucidating the molecular basis to abiotic stress tolerance in crops are becoming more easily achievable.

Mixed model association scans of multi-environmental trial data reveal major loci controlling yield and yield related traits in Hordeum vulgare in Mediterranean environments

An association panel consisting of 185 accessions representative of the barley germplasm cultivated in the Mediterranean basin was used to localise quantitative trait loci (QTL) controlling grain yield and yield related traits. The germplasm set was genotyped with 1,536 SNP markers and tested for associations with phenotypic data gathered over 2 years for a total of 24 year × location combinations under a broad range of environmental conditions. Analysis of multi-environmental trial (MET) data by fitting a mixed model with kinship estimates detected from two to seven QTL for the major components of yield including 1000 kernel weight, grains per spike and spikes per m(2), as well as heading date, harvest index and plant height. Several of the associations involved SNPs tightly linked to known major genes determining spike morphology in barley (vrs1 and int-c). Similarly, the largest QTL for heading date co-locates with SNPs linked with eam6, a major locus for heading date in barley for autumn sown conditions. Co-localization of several QTL related to yield components traits suggest that major developmental loci may be linked to most of the associations. This study highlights the potential of association genetics to identify genetic variants controlling complex traits.

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.

High levels of nucleotide diversity and fast decline of linkage disequilibrium in rye (Secale cereale L.) genes involved in frost response

BACKGROUND

Rye (Secale cereale L.) is the most frost tolerant cereal species. As an outcrossing species, rye exhibits high levels of intraspecific diversity, which makes it well-suited for allele mining in genes involved in the frost responsive network. For investigating genetic diversity and the extent of linkage disequilibrium (LD) we analyzed eleven candidate genes and 37 microsatellite markers in 201 lines from five Eastern and Middle European rye populations.

RESULTS

A total of 147 single nucleotide polymorphisms (SNPs) and nine insertion-deletion polymorphisms were found within 7,639 bp of DNA sequence from eleven candidate genes, resulting in an average SNP frequency of 1 SNP/52 bp. Nucleotide and haplotype diversity of candidate genes were high with average values π = 5.6 × 10(-3) and Hd = 0.59, respectively. According to an analysis of molecular variance (AMOVA), most of the genetic variation was found between individuals within populations. Haplotype frequencies varied markedly between the candidate genes. ScCbf14, ScVrn1, and ScDhn1 were dominated by a single haplotype, while the other 8 genes (ScCbf2, ScCbf6, ScCbf9b, ScCbf11, ScCbf12, ScCbf15, ScIce2, and ScDhn3) had a more balanced haplotype frequency distribution. Intra-genic LD decayed rapidly, within approximately 520 bp on average. Genome-wide LD based on microsatellites was low.

CONCLUSIONS

The Middle European population did not differ substantially from the four Eastern European populations in terms of haplotype frequencies or in the level of nucleotide diversity. The low LD in rye compared to self-pollinating species promises a high resolution in genome-wide association mapping. SNPs discovered in the promoters or coding regions, which attribute to non-synonymous substitutions, are suitable candidates for association mapping.

Population genetics of genomics-based crop improvement methods

Many genome-wide association studies (GWAS) in humans are concluding that, even with very large sample sizes and high marker densities, most of the genetic basis of complex traits may remain unexplained. At the same time, recent research in plant GWAS is showing much greater success with fewer resources. Both GWAS and genomic selection (GS), a method for predicting phenotypes by the use of genome-wide marker data, are receiving considerable attention among plant breeders. In this review we explore how differences in population genetic histories, as well as past selection for traits of interest, have produced trait architectures and patterns of linkage disequilibrium (LD) that frequently differ dramatically between domesticated plants and humans, making detection of quantitative trait loci (QTL) effects in crops more rewarding and less costly than in humans.

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.

Rapid linkage disequilibrium decay in the Lr10 gene in wild emmer wheat (Triticum dicoccoides) populations

INTRODUCTION

Recombination is a key evolutionary factor enhancing diversity. However, the effect of recombination on diversity in inbreeding species is expected to be low. To estimate this effect, recombination and diversity patterns of Lr10 gene were studied in natural populations of the inbreeder species, wild emmer wheat (Triticum dicoccoides). Wild emmer wheat is the progenitor of most cultivated wheats and it harbors rich genetic resources for disease resistance. Lr10 is a leaf rust resistance gene encoding three domains: a coiled-coil, nucleotide-binding site, and leucine-rich repeat (CC-NBS-LRR).

RESULTS

Lr10 was sequenced from 58 accessions representing 12 diverse habitats in Israel. Diversity analysis revealed a high rate of synonymous and non-synonymous substitutions (d (S) = 0.029, d (N) = 0.018, respectively) in the NBS-LRR domains. Moreover, in contrast to other resistance genes, in Lr10 the CC domain was more diverse than the NBS-LRR domains (d (S) = 0.069 vs. 0.029, d (N) = 0.094 vs. 0.018) and was subjected to positive selection in some of the populations. Seventeen recombination events were detected between haplotypes, especially in the CC domain. Linkage disequilibrium (LD) analysis has shown a rapid decay from r (2) = 0.5 to r (2) = 0.1 within a 2-kb span.

CONCLUSION

These results suggest that recombination is a diversifying force for the R-gene, Lr10, in the selfing species T. dicoccoides. This is the first report of a short-range LD decay in wild emmer wheat.

Chromosomal loci associated with endosperm hardness in a malting barley cross

A breeding objective for the malting barley industry is to produce lines with softer, plumper grain containing moderate protein content (9-12%) as they are more likely to imbibe water readily and contain more starch per grain, which in turn produces higher levels of malt extract. In a malting barley mapping population, 'Arapiles' × 'Franklin', the most significant and robust quantitative trait locus (QTL) for endosperm hardness was observed on the short arm of chromosome 1H, across three environments over two growing seasons. This accounted for 22.6% (Horsham 2000), 26.8% (Esperance 2001), and 12.0% (Tarranyurk 2001) of the genetic variance and significantly increased endosperm hardness by 2.06-3.03 SKCS hardness units. Interestingly, Arapiles and Franklin do not vary in Ha locus alleles. Therefore, this region, near the centromere on chromosome 1H, may be of great importance when aiming to manipulate endosperm hardness and malting quality. Interestingly, this region, close to the centromere on chromosome 1H, in our study, aligns with the region of the genome that includes the HvCslF9 and the HvGlb1 genes. Potentially, one or both of these genes could be considered to be candidate genes that influence endosperm hardness in the barley grain. Additional QTLs for endosperm hardness were detected on chromosomes 2H, 3H, 6H and 7H, confirming that the hardness trait in barley is complex and multigenic, similar to many malting quality traits of interest.

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.

Whole-genome association mapping in elite inbred crop varietiesThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”

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.

Genetic analysis and ecological association of Hina genes based on single nucleotide polymorphisms (SNPs) in wild barley, Hordeum spontaneum

Specific primers were designed to amplify the sequences of Hina genes from 121 wild barley (Hordeum spontaneum) accessions belonging to 18 populations from Iran, Israel and Turkey. Forty-nine single nucleotide polymorphisms (SNPs), nine indels, and 26 haplotypes were determined by sequence analysis. The genetic polymorphism (P), genetic diversity (He), and Shannon's information index (I) in the 18 populations were 0.486, 0.181 and 0.269, respectively. Approximately 2/3 genetic variations of Hina genes were presented within populations, while approximately 1/3 genetic variations were observed between populations. Broad gene flow (Nm= 3.31) and low genetic variation (Gst= 0.0702) were detected. However, the genetic differentiation between populations was independent of geographical distances according to the Mantel test (p = 0.478). The result of Spearman rank correlations (r(s)) showed that the genetic indices (P, He and I) of Hina were not significantly correlated with ecological factors. Only eight SNP positions correlated significantly with ecological factors. Of the eight SNP positions that positively correlated with ecological factors, only one SNP (769, T-C) was located in the coding region; however, it was not responsible for the amino acid change.

Association mapping of spot blotch resistance in wild barley

Spot blotch, caused by Cochliobolus sativus, is an important foliar disease of barley. The disease has been controlled for over 40 years through the deployment of cultivars with durable resistance derived from the line NDB112. Pathotypes of C. sativus with virulence for the NDB112 resistance have been detected in Canada; thus, many commercial cultivars are vulnerable to spot blotch epidemics. To increase the diversity of spot blotch resistance in cultivated barley, we evaluated 318 diverse wild barley accessions comprising the Wild Barley Diversity Collection (WBDC) for reaction to C. sativus at the seedling stage and utilized an association mapping (AM) approach to identify and map resistance loci. A high frequency of resistance was found in the WBDC as 95% (302/318) of the accessions exhibited low infection responses. The WBDC was genotyped with 558 Diversity Array Technology (DArT((R))) and 2,878 single nucleotide polymorphism (SNP) markers and subjected to structure analysis before running the AM procedure. Thirteen QTL for spot blotch resistance were identified with DArT and SNP markers. These QTL were found on chromosomes 1H, 2H, 3H, 5H, and 7H and explained from 2.3 to 3.9% of the phenotypic variance. Nearly half of the identified QTL mapped to chromosome bins where spot blotch resistance loci were previously reported, offering some validation for the AM approach. The other QTL mapped to unique genomic regions and may represent new spot blotch resistance loci. This study demonstrates that AM is an effective technique for identifying and mapping QTL for disease resistance in a wild crop progenitor. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-010-9402-8) contains supplementary material, which is available to authorized users.

Population genomic analysis of Tunisian Medicago truncatula reveals candidates for local adaptation

Genome-wide association studies rely upon segregating natural genetic variation, particularly the patterns of polymorphism and correlation between adjacent markers. To facilitate association studies in the model legume Medicago truncatula, we present a genome-scale polymorphism scan using existing Affymetrix microarrays. We develop and validate a method that uses a simple information-criteria algorithm to call polymorphism from microarray data without reliance on a reference genotype. We genotype 12 inbred M. truncatula lines sampled from four wild Tunisian populations and find polymorphisms at approximately 7% of features, comprising 31 419 probes. Only approximately 3% of these markers assort by population, and of these only 10% differentiate between populations from saline and non-saline sites. Fifty-two differentiated probes with unique genome locations correspond to 18 distinct genome regions. Sanger resequencing was used to characterize a subset of maker loci and develop a single nucleotide polymorphism (SNP)-typing assay that confirmed marker assortment by habitat in an independent sample of 33 individuals from the four populations. Genome-wide linkage disequilibrium (LD) extends on average for approximately 10 kb, falling to background levels by approximately 500 kb. A similar range of LD decay was observed in the 18 genome regions that assort by habitat; these LD blocks delimit candidate genes for local adaptation, many of which encode proteins with predicted functions in abiotic stress tolerance and are targets for functional genomic studies. Tunisian M. truncatula populations contain substantial amounts of genetic variation that is structured in relatively small LD blocks, suggesting a history of migration and recombination. These populations provide a strong resource for genome-wide association studies.