Peter E. Hedley

A physical, genetic and functional sequence assembly of the barley genome

Barley (Hordeum vulgare L.) is among the world’s earliest domesticated and most important crop plants. It is diploid with a large haploid genome of 5.1 gigabases (Gb). Here we present an integrated and ordered physical, genetic and functional sequence resource that describes the barley gene-space in a structured whole-genome context. We developed a physical map of 4.98 Gb, with more than 3.90 Gb anchored to a high-resolution genetic map. Projecting a deep whole-genome shotgun assembly, complementary DNA and deep RNA sequence data onto this framework supports 79,379 transcript clusters, including 26,159 ‘high-confidence’ genes with homology support from other plant genomes. Abundant alternative splicing, premature termination codons and novel transcriptionally active regions suggest that post-transcriptional processing forms an important regulatory layer. Survey sequences from diverse accessions reveal a landscape of extensive single-nucleotide variation. Our data provide a platform for both genome-assisted research and enabling contemporary crop improvement.

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

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

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

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

A chromosome conformation capture ordered sequence of the barley genome

Cereal grasses of the Triticeae tribe have been the major food source in temperate regions since the dawn of agriculture. Their large genomes are characterized by a high content of repetitive elements and large pericentromeric regions that are virtually devoid of meiotic recombination. Here we present a high-quality reference genome assembly for barley (Hordeum vulgare L.). We use chromosome conformation capture mapping to derive the linear order of sequences across the pericentromeric space and to investigate the spatial organization of chromatin in the nucleus at megabase resolution. The composition of genes and repetitive elements differs between distal and proximal regions. Gene family analyses reveal lineage-specific duplications of genes involved in the transport of nutrients to developing seeds and the mobilization of carbohydrates in grains. We demonstrate the importance of the barley reference sequence for breeding by inspecting the genomic partitioning of sequence variation in modern elite germplasm, highlighting regions vulnerable to genetic erosion.

Single-feature polymorphism discovery in the barley transcriptome

A probe-level model for analysis of GeneChip gene-expression data is presented which identified more than 10,000 single-feature polymorphisms (SFP) between two barley genotypes. The method has good sensitivity, as 67% of known single-nucleotide polymorphisms (SNP) were called as SFPs. This method is applicable to all oligonucleotide microarray data, accounts for SNP effects in gene-expression data and represents an efficient and versatile approach for highly parallel marker identification in large genomes.

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.

The Metabolic and Developmental Roles of Carotenoid Cleavage Dioxygenase4 from Potato

The factors that regulate storage organ carotenoid content remain to be fully elucidated, despite the nutritional and economic importance of this class of compound. Recent findings suggest that carotenoid pool size is determined, at least in part, by the activity of carotenoid cleavage dioxygenases. The aim of this study was to investigate whether Carotenoid Cleavage Dioxygenase4 (CCD4) activity affects potato (Solanum tuberosum) tuber carotenoid content. Microarray analysis revealed elevated expression of the potato CCD4 gene in mature tubers from white-fleshed cultivars compared with higher carotenoid yellow-fleshed tubers. The expression level of the potato CCD4 gene was down-regulated using an RNA interference (RNAi) approach in stable transgenic lines. Down-regulation in tubers resulted in an increased carotenoid content, 2- to 5-fold higher than in control plants. The increase in carotenoid content was mainly due to elevated violaxanthin content, implying that this carotenoid may act as the in vivo substrate. Although transcript level was also reduced in plant organs other than tubers, such as leaves, stems, and roots , there was no change in carotenoid content in these organs. However, carotenoid levels were elevated in flower petals from RNAi lines. As well as changes in tuber carotenoid content, tubers from RNAi lines exhibited phenotypes such as heat sprouting, formation of chain tubers, and an elongated shape. These results suggest that the product of the CCD4 reaction may be an important factor in tuber heat responses.

An evaluation of genotyping by sequencing (GBS) to map the Breviaristatum-e ( ari-e ) locus in cultivated barley

We explored the use of genotyping by sequencing (GBS) on a recombinant inbred line population (GPMx) derived from a cross between the two-rowed barley cultivar ‘Golden Promise’ (ari-e.GP/Vrs1) and the six-rowed cultivar ‘Morex’ (Ari-e/vrs1) to map plant height. We identified three Quantitative Trait Loci (QTL), the first in a region encompassing the spike architecture gene Vrs1 on chromosome 2H, the second in an uncharacterised centromeric region on chromosome 3H, and the third in a region of chromosome 5H coinciding with the previously described dwarfing gene Breviaristatum-e (Ari-e).BackgroundBarley cultivars in North-western Europe largely contain either of two dwarfing genes; Denso on chromosome 3H, a presumed ortholog of the rice green revolution gene OsSd1, or Breviaristatum-e (ari-e) on chromosome 5H. A recessive mutant allele of the latter gene, ari-e.GP, was introduced into cultivation via the cv. ‘Golden Promise’ that was a favourite of the Scottish malt whisky industry for many years and is still used in agriculture today.ResultsUsing GBS mapping data and phenotypic measurements we show that ari-e.GP maps to a small genetic interval on chromosome 5H and that alternative alleles at a region encompassing Vrs1 on 2H along with a region on chromosome 3H also influence plant height. The location of Ari-e is supported by analysis of near-isogenic lines containing different ari-e alleles. We explored use of the GBS to populate the region with sequence contigs from the recently released physically and genetically integrated barley genome sequence assembly as a step towards Ari-e gene identification.ConclusionsGBS was an effective and relatively low-cost approach to rapidly construct a genetic map of the GPMx population that was suitable for genetic analysis of row type and height traits, allowing us to precisely position ari-e.GP on chromosome 5H. Mapping resolution was lower than we anticipated. We found the GBS data more complex to analyse than other data types but it did directly provide linked SNP markers for subsequent higher resolution genetic analysis.

Identification and Characterisation CRN Effectors in Phytophthora capsici Shows Modularity and Functional Diversity

Phytophthora species secrete a large array of effectors during infection of their host plants. The Crinkler (CRN) gene family encodes a ubiquitous but understudied class of effectors with possible but as of yet unknown roles in infection. To appreciate CRN effector function in Phytophthora, we devised a simple Crn gene identification and annotation pipeline to improve effector prediction rates. We predicted 84 full-length CRN coding genes and assessed CRN effector domain diversity in sequenced Oomycete genomes. These analyses revealed evidence of CRN domain innovation in Phytophthora and expansion in the Peronosporales. We performed gene expression analyses to validate and define two classes of CRN effectors, each possibly contributing to infection at different stages. CRN localisation studies revealed that P. capsici CRN effector domains target the nucleus and accumulate in specific sub-nuclear compartments. Phenotypic analyses showed that few CRN domains induce necrosis when expressed in planta and that one cell death inducing effector, enhances P. capsici virulence on Nicotiana benthamiana. These results suggest that the CRN protein family form an important class of intracellular effectors that target the host nucleus during infection. These results combined with domain expansion in hemi-biotrophic and necrotrophic pathogens, suggests specific contributions to pathogen lifestyles. This work will bolster CRN identification efforts in other sequenced oomycete species and set the stage for future functional studies towards understanding CRN effector functions.

An eQTL Analysis of Partial Resistance to Puccinia hordei in Barley

Background Genetic resistance to barley leaf rust caused by Puccinia hordei involves both R genes and quantitative trait loci. The R genes provide higher but less durable resistance than the quantitative trait loci. Consequently, exploring quantitative or partial resistance has become a favorable alternative for controlling disease. Four quantitative trait loci for partial resistance to leaf rust have been identified in the doubled haploid Steptoe (St)/Morex (Mx) mapping population. Further investigations are required to study the molecular mechanisms underpinning partial resistance and ultimately identify the causal genes. Methodology/Principal Findings We explored partial resistance to barley leaf rust using a genetical genomics approach. We recorded RNA transcript abundance corresponding to each probe on a 15K Agilent custom barley microarray in seedlings from St and Mx and 144 doubled haploid lines of the St/Mx population. A total of 1154 and 1037 genes were, respectively, identified as being P. hordei-responsive among the St and Mx and differentially expressed between P. hordei-infected St and Mx. Normalized ratios from 72 distant-pair hybridisations were used to map the genetic determinants of variation in transcript abundance by expression quantitative trait locus (eQTL) mapping generating 15685 eQTL from 9557 genes. Correlation analysis identified 128 genes that were correlated with resistance, of which 89 had eQTL co-locating with the phenotypic quantitative trait loci (pQTL). Transcript abundance in the parents and conservation of synteny with rice allowed us to prioritise six genes as candidates for Rphq11, the pQTL of largest effect, and highlight one, a phospholipid hydroperoxide glutathione peroxidase (HvPHGPx) for detailed analysis. Conclusions/Significance The eQTL approach yielded information that led to the identification of strong candidate genes underlying pQTL for resistance to leaf rust in barley and on the general pathogen response pathway. The dataset will facilitate a systems appraisal of this host-pathogen interaction and, potentially, for other traits measured in this population.