Natural variation in meiotic recombination rate shapes introgression patterns in intraspecific hybrids between wild and domesticated barley.

Meiotic recombination rates vary considerably between species, populations and individuals. The genetic exchange between homologous chromosomes plays a major role in evolution by breaking linkage between advantageous and deleterious alleles in the case of introgressions. Identifying recombination rate modifiers is thus of both fundamental and practical interest to understand and utilize variation in meiotic recombination rates. We investigated recombination rate variation in a large intraspecific hybrid population (named HEB-25) derived from a cross between domesticated barley and 25 wild barley accessions. We observed quantitative variation in total crossover number with a maximum of a 1.4-fold difference between subpopulations and increased recombination rates across pericentromeric regions. The meiosis-specific α-kleisin cohesin subunit REC8 was identified as a candidate gene influencing crossover number and patterning. Furthermore, we quantified wild barley introgression patterns and revealed how local and genome-wide recombination rate variation shapes patterns of introgression. The identification of allelic variation in REC8 in combination with the observed changes in crossover patterning suggest a difference in how chromatin loops are tethered to the chromosome axis, resulting in reduced crossover suppression across pericentromeric regions. Local and genome-wide recombination rate variation is shaping patterns of introgressions and thereby directly influences the consequences of linkage drag.

Genome-wide association of yield traits in a nested association mapping population of barley reveals new gene diversity for future breeding.

To explore wild barley as a source of useful alleles for yield improvement in breeding, we have carried out a genome-wide association scan using the nested association mapping population HEB-25, which contains 25 diverse exotic barley genomes superimposed on an ~70% genetic background of cultivated barley. A total of 1420 HEB-25 lines were trialled for nine yield-related grain traits for 2 years in Germany and Scotland, with varying N fertilizer application. The phenotypic data were related to genotype scores for 5398 gene-based single nucleotide polymorphism (SNP) markers. A total of 96 quantitative trait locus (QTL) regions were identified across all measured traits, the majority of which co-localize with known major genes controlling flowering time (Ppd-H2, HvCEN, HvGI, VRN-H1, and VRN-H3) and spike morphology (VRS3, VRS1, VRS4, and INT-C) in barley. Fourteen QTL hotspots, with at least three traits coinciding, were also identified, several of which co-localize with barley orthologues of genes controlling grain dimensions in rice. Most of the allele effects are specific to geographical location and/or exotic parental genotype. This study shows the existence of beneficial alleles for yield-related traits in exotic barley germplasm and provides candidate alleles for future improvement of these traits by the breeder.

Contrasting genetic regulation of plant development in wild barley grown in two European environments revealed by nested association mapping.

Barley is cultivated more widely than the other major world crops because it adapts well to environmental constraints, such as drought, heat, and day length. To better understand the genetic control of local adaptation in barley, we studied development in the nested association mapping population HEB-25, derived from crossing 25 wild barley accessions with the cultivar 'Barke'. HEB-25 was cultivated in replicated field trials in Dundee (Scotland) and Halle (Germany), differing in regard to day length, precipitation, and temperature. Applying a genome-wide association study, we located 60 and 66 quantitative trait locus (QTL) regions regulating eight plant development traits in Dundee and Halle, respectively. A number of QTLs could be explained by known major genes such as PHOTOPERIOD 1 (Ppd-H1) and FLOWERING LOCUS T (HvFT-1) that regulate plant development. In addition, we observed that developmental traits in HEB-25 were partly controlled via genotype × environment and genotype × donor interactions, defined as location-specific and family-specific QTL effects. Our findings indicate that QTL alleles are available in the wild barley gene pool that show contrasting effects on plant development, which may be deployed to improve adaptation of cultivated barley to future environmental changes.

Chromatin state analysis of the barley epigenome reveals a higher-order structure defined by H3K27me1 and H3K27me3 abundance.

Combinations of histones carrying different covalent modifications are a major component of epigenetic variation. We have mapped nine modified histones in the barley seedling epigenome by chromatin immunoprecipitation next-generation sequencing (ChIP-seq). The chromosomal distributions of the modifications group them into four different classes, and members of a given class also tend to coincide at the local DNA level, suggesting that global distribution patterns reflect local epigenetic environments. We used this peak sharing to define 10 chromatin states representing local epigenetic environments in the barley genome. Five states map mainly to genes and five to intergenic regions. Two genic states involving H3K36me3 are preferentially associated with constitutive gene expression, while an H3K27me3-containing genic state is associated with differentially expressed genes. The 10 states display striking distribution patterns that divide barley chromosomes into three distinct global environments. First, telomere-proximal regions contain high densities of H3K27me3 covering both genes and intergenic DNA, together with very low levels of the repressive H3K27me1 modification. Flanking these are gene-rich interior regions that are rich in active chromatin states and have greatly decreased levels of H3K27me3 and increasing amounts of H3K27me1 and H3K9me2. Lastly, H3K27me3-depleted pericentromeric regions contain gene islands with active chromatin states separated by extensive retrotransposon-rich regions that are associated with abundant H3K27me1 and H3K9me2 modifications. We propose an epigenomic framework for barley whereby intergenic H3K27me3 specifies facultative heterochromatin in the telomere-proximal regions and H3K27me1 is diagnostic for constitutive heterochromatin elsewhere in the barley genome.

An investigation of causes of false positive single nucleotide polymorphisms using simulated reads from a small eukaryote genome.

BACKGROUND: Single Nucleotide Polymorphisms (SNPs) are widely used molecular markers, and their use has increased massively since the inception of Next Generation Sequencing (NGS) technologies, which allow detection of large numbers of SNPs at low cost. However, both NGS data and their analysis are error-prone, which can lead to the generation of false positive (FP) SNPs. We explored the relationship between FP SNPs and seven factors involved in mapping-based variant calling - quality of the reference sequence, read length, choice of mapper and variant caller, mapping stringency and filtering of SNPs by read mapping quality and read depth. This resulted in 576 possible factor level combinations. We used error- and variant-free simulated reads to ensure that every SNP found was indeed a false positive. RESULTS: The variation in the number of FP SNPs generated ranged from 0 to 36,621 for the 120 million base pairs (Mbp) genome. All of the experimental factors tested had statistically significant effects on the number of FP SNPs generated and there was a considerable amount of interaction between the different factors. Using a fragmented reference sequence led to a dramatic increase in the number of FP SNPs generated, as did relaxed read mapping and a lack of SNP filtering. The choice of reference assembler, mapper and variant caller also significantly affected the outcome. The effect of read length was more complex and suggests a possible interaction between mapping specificity and the potential for contributing more false positives as read length increases. CONCLUSIONS: The choice of tools and parameters involved in variant calling can have a dramatic effect on the number of FP SNPs produced, with particularly poor combinations of software and/or parameter settings yielding tens of thousands in this experiment. Between-factor interactions make simple recommendations difficult for a SNP discovery pipeline but the quality of the reference sequence is clearly of paramount importance. Our findings are also a stark reminder that it can be unwise to use the relaxed mismatch settings provided as defaults by some read mappers when reads are being mapped to a relatively unfinished reference sequence from e.g. a non-model organism in its early stages of genomic exploration.

The low-recombining pericentromeric region of barley restricts gene diversity and evolution but not gene expression.

The low-recombining pericentromeric region of the barley genome contains roughly a quarter of the genes of the species, embedded in low-recombining DNA that is rich in repeats and repressive chromatin signatures. We have investigated the effects of pericentromeric region residency upon the expression, diversity and evolution of these genes. We observe no significant difference in average transcript level or developmental RNA specificity between the barley pericentromeric region and the rest of the genome. In contrast, all of the evolutionary parameters studied here show evidence of compromised gene evolution in this region. First, genes within the pericentromeric region of wild barley show reduced diversity and significantly weakened purifying selection compared with the rest of the genome. Second, gene duplicates (ohnolog pairs) derived from the cereal whole-genome duplication event ca. 60MYa have been completely eliminated from the barley pericentromeric region. Third, local gene duplication in the pericentromeric region is reduced by 29% relative to the rest of the genome. Thus, the pericentromeric region of barley is a permissive environment for gene expression but has restricted gene evolution in a sizeable fraction of barley's genes.

Barley whole exome capture: a tool for genomic research in the genus Hordeum and beyond.

Advanced resources for genome-assisted research in barley (Hordeum vulgare) including a whole-genome shotgun assembly and an integrated physical map have recently become available. These have made possible studies that aim to assess genetic diversity or to isolate single genes by whole-genome resequencing and in silico variant detection. However such an approach remains expensive given the 5 Gb size of the barley genome. Targeted sequencing of the mRNA-coding exome reduces barley genomic complexity more than 50-fold, thus dramatically reducing this heavy sequencing and analysis load. We have developed and employed an in-solution hybridization-based sequence capture platform to selectively enrich for a 61.6 megabase coding sequence target that includes predicted genes from the genome assembly of the cultivar Morex as well as publicly available full-length cDNAs and de novo assembled RNA-Seq consensus sequence contigs. The platform provides a highly specific capture with substantial and reproducible enrichment of targeted exons, both for cultivated barley and related species. We show that this exome capture platform provides a clear path towards a broader and deeper understanding of the natural variation residing in the mRNA-coding part of the barley genome and will thus constitute a valuable resource for applications such as mapping-by-sequencing and genetic diversity analyzes.

A unified classification system for eukaryotic transposable elements.

Our knowledge of the structure and composition of genomes is rapidly progressing in pace with their sequencing. The emerging data show that a significant portion of eukaryotic genomes is composed of transposable elements (TEs). Given the abundance and diversity of TEs and the speed at which large quantities of sequence data are emerging, identification and annotation of TEs presents a significant challenge. Here we propose the first unified hierarchical classification system, designed on the basis of the transposition mechanism, sequence similarities and structural relationships, that can be easily applied by non-experts. The system and nomenclature is kept up to date at the WikiPoson web site.

Islands and streams: clusters and gene flow in wild barley populations from the Levant.

The domestication of plants frequently results in a high level of genetic differentiation between domesticated plants and their wild progenitors. This process is counteracted by gene flow between wild and domesticated plants because they are usually able to inter-mate and to exchange genes. We investigated the extent of gene flow between wild barley Hordeum spontaneum and cultivated barley Hordeum vulgare, and its effect on population structure in wild barley by analysing a collection of 896 wild barley accessions (Barley1K) from Israel and all available Israeli H. vulgare accessions from the Israeli gene bank. We compared the performance of simple sequence repeats (SSR) and single nucleotide polymorphisms (SNP) marker data genotyped over a core collection in estimating population parameters. Estimates of gene flow rates with SSR markers indicated a high level of introgression from cultivated barley into wild barley. After removing accessions from the wild barley sample that were recently admixed with cultivated barley, the inference of population structure improved significantly. Both SSR and SNP markers showed that the genetic population structure of wild barley in Israel corresponds to the three major ecogeographic regions: the coast, the Mediterranean north and the deserts in the Jordan valley and the South. Gene flow rates were estimated to be higher from north to south than in the opposite direction. As has been observed in other crop species, there is a significant exchange of alleles between the wild species and domesticated varieties that needs to be accounted for in the population genetic analysis of domestication.

Flapjack--graphical genotype visualization.

SUMMARY: New software tools for graphical genotyping are required that can routinely handle the large data volumes generated by the high-throughput single-nucleotide polymorphism (SNP) platforms, genotyping-by-sequencing and other comparable genotyping technologies. Flapjack has been developed to facilitate analysis of these data, providing real time rendering with rapid navigation and comparisons between lines, markers and chromosomes, with visualization, sorting and querying based on associated data, such as phenotypes, quantitative trait loci or other mappable features. AVAILABILITY: Flapjack is freely available for Microsoft Windows, Mac OS X, Linux and Solaris, and can be downloaded from http://bioinf.scri.ac.uk/flapjack .

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.

The genetic diversity and evolution of field pea (Pisum) studied by high throughput retrotransposon based insertion polymorphism (RBIP) marker analysis.

BACKGROUND: The genetic diversity of crop species is the result of natural selection on the wild progenitor and human intervention by ancient and modern farmers and breeders. The genomes of modern cultivars, old cultivated landraces, ecotypes and wild relatives reflect the effects of these forces and provide insights into germplasm structural diversity, the geographical dimension to species diversity and the process of domestication of wild organisms. This issue is also of great practical importance for crop improvement because wild germplasm represents a rich potential source of useful under-exploited alleles or allele combinations. The aim of the present study was to analyse a major Pisum germplasm collection to gain a broad understanding of the diversity and evolution of Pisum and provide a new rational framework for designing germplasm core collections of the genus. RESULTS: 3020 Pisum germplasm samples from the John Innes Pisum germplasm collection were genotyped for 45 retrotransposon based insertion polymorphism (RBIP) markers by the Tagged Array Marker (TAM) method. The data set was stored in a purpose-built Germinate relational database and analysed by both principal coordinate analysis and a nested application of the Structure program which yielded substantially similar but complementary views of the diversity of the genus Pisum. Structure revealed three Groups (1-3) corresponding approximately to landrace, cultivar and wild Pisum respectively, which were resolved by nested Structure analysis into 14 Sub-Groups, many of which correlate with taxonomic sub-divisions of Pisum, domestication related phenotypic traits and/or restricted geographical locations. Genetic distances calculated between these Sub-Groups are broadly supported by principal coordinate analysis and these, together with the trait and geographical data, were used to infer a detailed model for the domestication of Pisum. CONCLUSIONS: These data provide a clear picture of the major distinct gene pools into which the genus Pisum is partitioned and their geographical distribution. The data strongly support the model of independent domestications for P. sativum ssp abyssinicum and P. sativum. The relationships between these two cultivated germplasms and the various sub-divisions of wild Pisum have been clarified and the most likely ancestral wild gene pools for domesticated P. sativum identified. Lastly, this study provides a framework for defining global Pisum germplasm which will be useful for designing core collections.

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.

Optimizing the procedure of grain nutrient predictions in barley via hyperspectral imaging.

Hyperspectral imaging enables researchers and plant breeders to analyze various traits of interest like nutritional value in high throughput. In order to achieve this, the optimal design of a reliable calibration model, linking the measured spectra with the investigated traits, is necessary. In the present study we investigated the impact of different regression models, calibration set sizes and calibration set compositions on prediction performance. For this purpose, we analyzed concentrations of six globally relevant grain nutrients of the wild barley population HEB-YIELD as case study. The data comprised 1,593 plots, grown in 2015 and 2016 at the locations Dundee and Halle, which have been entirely analyzed through traditional laboratory methods and hyperspectral imaging. The results indicated that a linear regression model based on partial least squares outperformed neural networks in this particular data modelling task. There existed a positive relationship between the number of samples in a calibration model and prediction performance, with a local optimum at a calibration set size of ~40% of the total data. The inclusion of samples from several years and locations could clearly improve the predictions of the investigated nutrient traits at small calibration set sizes. It should be stated that the expansion of calibration models with additional samples is only useful as long as they are able to increase trait variability. Models obtained in a certain environment were only to a limited extent transferable to other environments. They should therefore be successively upgraded with new calibration data to enable a reliable prediction of the desired traits. The presented results will assist the design and conceptualization of future hyperspectral imaging projects in order to achieve reliable predictions. It will in general help to establish practical applications of hyperspectral imaging systems, for instance in plant breeding concepts.

"Wild barley serves as a source for biofortification of barley grains".

The continuing growth of the human population creates an inevitable necessity for higher crop yields, which are mandatory for the supply with adequate amounts of food. However, increasing grain yield may lead to a reduction of grain quality, such as a decline in protein and mineral nutrient concentrations causing the so-called hidden hunger. To assess the interdependence between quantity and quality and to evaluate the biofortification potential of wild barley, we conducted field studies, examining the interplay between plant development, yield, and nutrient concentrations, using HEB-YIELD, a subset of the wild barley nested association mapping population HEB-25. A huge variation of nutrient concentration in grains was obtained, since we identified lines with a more than 50% higher grain protein, iron, and zinc concentration in comparison to the recurrent parent 'Barke'. We observed a negative relationship between grain yield and nutritional value in barley, indicated by predominantly negative correlations between yield and nutrient concentrations. Analyzing the genetic control of nutrient concentration in mature grains indicated that numerous genomic regions determine the final nutritional value of grains and wild alleles were frequently associated with higher nutrient concentrations. The targeted introgression of wild barley alleles may enable biofortification in future barley breeding.

Barley yield formation under abiotic stress depends on the interplay between flowering time genes and environmental cues.

Since the dawn of agriculture, crop yield has always been impaired through abiotic stresses. In a field trial across five locations worldwide, we tested three abiotic stresses, nitrogen deficiency, drought and salinity, using HEB-YIELD, a selected subset of the wild barley nested association mapping population HEB-25. We show that barley flowering time genes Ppd-H1, Sdw1, Vrn-H1 and Vrn-H3 exert pleiotropic effects on plant development and grain yield. Under field conditions, these effects are strongly influenced by environmental cues like day length and temperature. For example, in Al-Karak, Jordan, the day length-sensitive wild barley allele of Ppd-H1 was associated with an increase of grain yield by up to 30% compared to the insensitive elite barley allele. The observed yield increase is accompanied by pleiotropic effects of Ppd-H1 resulting in shorter life cycle, extended grain filling period and increased grain size. Our study indicates that the adequate timing of plant development is crucial to maximize yield formation under harsh environmental conditions. We provide evidence that wild barley alleles, introgressed into elite barley cultivars, can be utilized to support grain yield formation. The presented knowledge may be transferred to related crop species like wheat and rice securing the rising global food demand for cereals.

Gene-based sequence diversity analysis of field pea (Pisum).

Sequence diversity of 39 dispersed gene loci was analyzed in 48 diverse individuals representative of the genus Pisum. The different genes show large variation in diversity parameters, suggesting widely differing levels of selection and a high overall diversity level for the species. The data set yields a genetic diversity tree whose deep branches, involving wild samples, are preserved in a tree derived from a polymorphic retrotransposon insertions in an identical sample set. Thus, gene regions and intergenic "junk DNA" share a consistent picture for the genomic diversity of Pisum, despite low linkage disequilibrium in wild and landrace germplasm, which might be expected to allow independent evolution of these very different DNA classes. Additional lines of evidence indicate that recombination has shuffled gene haplotypes efficiently within Pisum, despite its high level of inbreeding and widespread geographic distribution. Trees derived from individual gene loci show marked differences from each other, and genetic distance values between sample pairs show high standard deviations. Sequence mosaic analysis of aligned sequences identifies nine loci showing evidence for intragenic recombination. Lastly, phylogenetic network analysis confirms the non-treelike structure of Pisum diversity and indicates the major germplasm classes involved. Overall, these data emphasize the artificiality of simple tree structures for representing genomic sequence variation within Pisum and emphasize the need for fine structure haplotype analysis to accurately define the genetic structure of the species.

Involvement of Serotonergic and Relaxin-3 Neuropeptide Systems in the Expression of Anxiety-like Behavior.

Anxiety-related defensive behavior is controlled by a distributed network of brain regions and interconnected neural circuits. The dorsal raphe nucleus (DR), which contains the majority of forebrain-projecting serotonergic neurons, is a key brain region involved in fear states and anxiety-related behavior via modulation of this broad neural network. Evidence suggests that relaxin-3 neurons in the nucleus incertus (NI) may also interact with this network, however, the potential role of the NI in the control of anxiety-related defensive behavior requires further investigation. In this study, we examined the response of an anxiety-related neuronal network, including serotonergic neurons in the DR and relaxin-3-containing neurons in the NI, to administration of an anxiogenic drug and exposure to an aversive environment. We administered an anxiogenic dose of the adenosine receptor antagonist, caffeine (50 mg/kg, i.p.), or vehicle, to adult male Wistar rats and 30 min later exposed them to either an elevated plus-maze (EPM) or a home cage environment. Administration of caffeine and exposure to the EPM activated a broad network of brain regions involved in control of anxiety-like behaviors, including serotonergic neurons in the DR, as measured using c-Fos immunohistochemistry. However, only exposure to the EPM activated relaxin-3-containing neurons in the NI, and activation of these neurons was not correlated with changes in anxiety-like behavior. These data suggest activation of the NI relaxin-3 system is associated with expression of behavior in tests of anxiety, but may not be directly involved in the approach-avoidance conflict inherent in anxiety-related defensive behavior in rodents.

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.

Genome-wide association mapping in winter barley for grain yield and culm cell wall polymer content using the high-throughput CoMPP technique.

A collection of 112 winter barley varieties (Hordeum vulgare L.) was grown in the field for two years (2008/09 and 2009/10) in northern Italy and grain and straw yields recorded. In the first year of the trial, a severe attack of barley yellow mosaic virus (BaYMV) strongly influenced final performances with an average reduction of ~ 50% for grain and straw harvested in comparison to the second year. The genetic determination (GD) for grain yield was 0.49 and 0.70, for the two years respectively, and for straw yield GD was low in 2009 (0.09) and higher in 2010 (0.29). Cell wall polymers in culms were quantified by means of the monoclonal antibodies LM6, LM11, JIM13 and BS-400-3 and the carbohydrate-binding module CBM3a using the high-throughput CoMPP technique. Of these, LM6, which detects arabinan components, showed a relatively high GD in both years and a significantly negative correlation with grain yield (GYLD). Overall, heritability (H2) was calculated for GYLD, LM6 and JIM and resulted to be 0.42, 0.32 and 0.20, respectively. A total of 4,976 SNPs from the 9K iSelect array were used in the study for the analysis of population structure, linkage disequilibrium (LD) and genome-wide association study (GWAS). Marker-trait associations (MTA) were analyzed for grain yield and cell wall determination by LM6 and JIM13 as these were the traits showing significant correlations between the years. A single QTL for GYLD containing three MTAs was found on chromosome 3H located close to the Hv-eIF4E gene, which is known to regulate resistance to BaYMV. Subsequently the QTL was shown to be tightly linked to rym4, a locus for resistance to the virus. GWAs on arabinans quantified by LM6 resulted in the identification of major QTLs closely located on 3H and hypotheses regarding putative candidate genes were formulated through the study of gene expression levels based on bioinformatics tools.