SMRT-AgRenSeq-d in potato (Solanum tuberosum) as a method to identify candidates for the nematode resistance Gpa5.

Potato is the third most important food crop in the world. Diverse pathogens threaten sustainable crop production but can be controlled, in many cases, through the deployment of disease resistance genes belonging to the family of nucleotide-binding, leucine-rich-repeat (NLR) genes. To identify effective disease resistance genes in established varieties, we have successfully established SMRT-AgRenSeq in tetraploid potatoes and have further enhanced the methodology by including dRenSeq in an approach that we term SMR-AgRenSeq-d. The inclusion of dRenSeq enables the filtering of candidates after the association analysis by establishing a presence/absence matrix across resistant and susceptible varieties that is translated into an F1 score. Using a SMRT-RenSeq-based sequence representation of the NLRome from the cultivar Innovator, SMRT-AgRenSeq-d analyses reliably identified the late blight resistance benchmark genes Rpi-R1, Rpi-R2-like, Rpi-R3a, and Rpi-R3b in a panel of 117 varieties with variable phenotype penetrations. All benchmark genes were identified with an F1 score of 1, which indicates absolute linkage in the panel. This method also identified nine strong candidates for Gpa5 that controls the potato cyst nematode (PCN) species Globodera pallida (pathotypes Pa2/3). Assuming that NLRs are involved in controlling many types of resistances, SMRT-AgRenSeq-d can readily be applied to diverse crops and pathogen systems.

HISS: Snakemake-based workflows for performing SMRT-RenSeq assembly, AgRenSeq and dRenSeq for the discovery of novel plant disease resistance genes.

BACKGROUND: In the ten years since the initial publication of the RenSeq protocol, the method has proved to be a powerful tool for studying disease resistance in plants and providing target genes for breeding programmes. Since the initial publication of the methodology, it has continued to be developed as new technologies have become available and the increased availability of computing power has made new bioinformatic approaches possible. Most recently, this has included the development of a k-mer based association genetics approach, the use of PacBio HiFi data, and graphical genotyping with diagnostic RenSeq. However, there is not yet a unified workflow available and researchers must instead configure approaches from various sources themselves. This makes reproducibility and version control a challenge and limits the ability to perform these analyses to those with bioinformatics expertise. RESULTS: Here we present HISS, consisting of three workflows which take a user from raw RenSeq reads to the identification of candidates for disease resistance genes. These workflows conduct the assembly of enriched HiFi reads from an accession with the resistance phenotype of interest. A panel of accessions both possessing and lacking the resistance are then used in an association genetics approach (AgRenSeq) to identify contigs positively associated with the resistance phenotype. Candidate genes are then identified on these contigs and assessed for their presence or absence in the panel with a graphical genotyping approach that uses dRenSeq. These workflows are implemented via Snakemake, a python-based workflow manager. Software dependencies are either shipped with the release or handled with conda. All code is freely available and is distributed under the GNU GPL-3.0 license. CONCLUSIONS: HISS provides a user-friendly, portable, and easily customised approach for identifying novel disease resistance genes in plants. It is easily installed with all dependencies handled internally or shipped with the release and represents a significant improvement in the ease of use of these bioinformatics analyses.

Conserved signalling components coordinate epidermal patterning and cuticle deposition in barley.

Faced with terrestrial threats, land plants seal their aerial surfaces with a lipid-rich cuticle. To breathe, plants interrupt their cuticles with adjustable epidermal pores, called stomata, that regulate gas exchange, and develop other specialised epidermal cells such as defensive hairs. Mechanisms coordinating epidermal features remain poorly understood. Addressing this, we studied two loci whose allelic variation causes both cuticular wax-deficiency and misarranged stomata in barley, identifying the underlying genes, Cer-g/ HvYDA1, encoding a YODA-like (YDA) MAPKKK, and Cer-s/ HvBRX-Solo, encoding a single BREVIS-RADIX (BRX) domain protein. Both genes control cuticular integrity, the spacing and identity of epidermal cells, and barley's distinctive epicuticular wax blooms, as well as stomatal patterning in elevated CO2 conditions. Genetic analyses revealed epistatic and modifying relationships between HvYDA1 and HvBRX-Solo, intimating that their products participate in interacting pathway(s) linking epidermal patterning with cuticular properties in barley. This may represent a mechanism for coordinating multiple adaptive features of the land plant epidermis in a cultivated cereal.

The evolutionary patterns of barley pericentromeric chromosome regions, as shaped by linkage disequilibrium and domestication.

The distribution of recombination events along large cereal chromosomes is uneven and is generally restricted to gene-rich telomeric ends. To understand how the lack of recombination affects diversity in the large pericentromeric regions, we analysed deep exome capture data from a final panel of 815 Hordeum vulgare (barley) cultivars, landraces and wild barleys, sampled from across their eco-geographical ranges. We defined and compared variant data across the pericentromeric and non-pericentromeric regions, observing a clear partitioning of diversity both within and between chromosomes and germplasm groups. Dramatically reduced diversity was found in the pericentromeres of both cultivars and landraces when compared with wild barley. We observed a mixture of completely and partially differentiated single-nucleotide polymorphisms (SNPs) between domesticated and wild gene pools, suggesting that domesticated gene pools were derived from multiple wild ancestors. Patterns of genome-wide linkage disequilibrium, haplotype block size and number, and variant frequency within blocks showed clear contrasts among individual chromosomes and between cultivars and wild barleys. Although most cultivar chromosomes shared a single major pericentromeric haplotype, chromosome 7H clearly differentiated the two-row and six-row types associated with different geographical origins. Within the pericentromeric regions we identified 22 387 non-synonymous SNPs, 92 of which were fixed for alternative alleles in cultivar versus wild accessions. Surprisingly, only 29 SNPs found exclusively in the cultivars were predicted to be 'highly deleterious'. Overall, our data reveal an unconventional pericentromeric genetic landscape among distinct barley gene pools, with different evolutionary processes driving domestication and diversification.

BaRTv2: a highly resolved barley reference transcriptome for accurate transcript-specific RNA-seq quantification.

Accurate characterisation of splice junctions (SJs) as well as transcription start and end sites in reference transcriptomes allows precise quantification of transcripts from RNA-seq data, and enables detailed investigations of transcriptional and post-transcriptional regulation. Using novel computational methods and a combination of PacBio Iso-seq and Illumina short-read sequences from 20 diverse tissues and conditions, we generated a comprehensive and highly resolved barley reference transcript dataset from the European 2-row spring barley cultivar Barke (BaRTv2.18). Stringent and thorough filtering was carried out to maintain the quality and accuracy of the SJs and transcript start and end sites. BaRTv2.18 shows increased transcript diversity and completeness compared with an earlier version, BaRTv1.0. The accuracy of transcript level quantification, SJs and transcript start and end sites have been validated extensively using parallel technologies and analysis, including high-resolution reverse transcriptase-polymerase chain reaction and 5'-RACE. BaRTv2.18 contains 39 434 genes and 148 260 transcripts, representing the most comprehensive and resolved reference transcriptome in barley to date. It provides an important and high-quality resource for advanced transcriptomic analyses, including both transcriptional and post-transcriptional regulation, with exceptional resolution and precision.

EORNA, a barley gene and transcript abundance database.

A high-quality, barley gene reference transcript dataset (BaRTv1.0), was used to quantify gene and transcript abundances from 22 RNA-seq experiments, covering 843 separate samples. Using the abundance data we developed a Barley Expression Database (EORNA*) to underpin a visualisation tool that displays comparative gene and transcript abundance data on demand as transcripts per million (TPM) across all samples and all the genes. EORNA provides gene and transcript models for all of the transcripts contained in BaRTV1.0, and these can be conveniently identified through either BaRT or HORVU gene names, or by direct BLAST of query sequences. Browsing the quantification data reveals cultivar, tissue and condition specific gene expression and shows changes in the proportions of individual transcripts that have arisen via alternative splicing. TPM values can be easily extracted to allow users to determine the statistical significance of observed transcript abundance variation among samples or perform meta analyses on multiple RNA-seq experiments. * Eòrna is the Scottish Gaelic word for Barley.

Barley sodium content is regulated by natural variants of the Na+ transporter HvHKT1;5.

During plant growth, sodium (Na+) in the soil is transported via the xylem from the root to the shoot. While excess Na+ is toxic to most plants, non-toxic concentrations have been shown to improve crop yields under certain conditions, such as when soil K+ is low. We quantified grain Na+ across a barley genome-wide association study panel grown under non-saline conditions and identified variants of a Class 1 HIGH-AFFINITY-POTASSIUM-TRANSPORTER (HvHKT1;5)-encoding gene responsible for Na+ content variation under these conditions. A leucine to proline substitution at position 189 (L189P) in HvHKT1;5 disturbs its characteristic plasma membrane localisation and disrupts Na+ transport. Under low and moderate soil Na+, genotypes containing HvHKT1:5P189 accumulate high concentrations of Na+ but exhibit no evidence of toxicity. As the frequency of HvHKT1:5P189 increases significantly in cultivated European germplasm, we cautiously speculate that this non-functional variant may enhance yield potential in non-saline environments, possibly by offsetting limitations of low available K+.

Characterisation of barley landraces from Syria and Jordan for resistance to rhynchosporium and identification of diagnostic markers for Rrs1Rh4.

KEY MESSAGE: Diagnostic markers for Rrs1Rh4 have been identified by testing for associations between SNPs within the Rrs1 interval in 150 barley genotypes and their resistance to Rhynchosporium commune isolates recognised by lines containing Rrs1. Rhynchosporium or barley scald, caused by the destructive fungal pathogen Rhynchosporium commune, is one of the most economically important diseases of barley in the world. Barley landraces from Syria and Jordan demonstrated high resistance to rhynchosporium in the field. Genotyping of a wide range of barley cultivars and landraces, including known sources of different Rrs1 genes/alleles, across the Rrs1 interval, followed by association analysis of this genotypic data with resistance phenotypes to R. commune isolates recognised by Rrs1, allowed the identification of diagnostic markers for Rrs1Rh4. These markers are specific to Rrs1Rh4 and do not detect other Rrs1 genes/alleles. The Rrs1Rh4 diagnostic markers represent a resource that can be exploited by breeders for the sustainable deployment of varietal resistance in new cultivars. Thirteen out of the 55 most resistant Syrian and Jordanian landraces were shown to contain markers specific to Rrs1Rh4. One of these lines came from Jordan, with the remaining 12 lines from different locations in Syria. One of the Syrian landraces containing Rrs1Rh4 was also shown to have Rrs2. The remaining landraces that performed well against rhynchosporium in the field are likely to contain other resistance genes and represent an important novel resource yet to be exploited by European breeders.

BaRTv1.0: an improved barley reference transcript dataset to determine accurate changes in the barley transcriptome using RNA-seq.

BACKGROUND: The time required to analyse RNA-seq data varies considerably, due to discrete steps for computational assembly, quantification of gene expression and splicing analysis. Recent fast non-alignment tools such as Kallisto and Salmon overcome these problems, but these tools require a high quality, comprehensive reference transcripts dataset (RTD), which are rarely available in plants. RESULTS: A high-quality, non-redundant barley gene RTD and database (Barley Reference Transcripts - BaRTv1.0) has been generated. BaRTv1.0, was constructed from a range of tissues, cultivars and abiotic treatments and transcripts assembled and aligned to the barley cv. Morex reference genome (Mascher et al. Nature; 544: 427-433, 2017). Full-length cDNAs from the barley variety Haruna nijo (Matsumoto et al. Plant Physiol; 156: 20-28, 2011) determined transcript coverage, and high-resolution RT-PCR validated alternatively spliced (AS) transcripts of 86 genes in five different organs and tissue. These methods were used as benchmarks to select an optimal barley RTD. BaRTv1.0-Quantification of Alternatively Spliced Isoforms (QUASI) was also made to overcome inaccurate quantification due to variation in 5' and 3' UTR ends of transcripts. BaRTv1.0-QUASI was used for accurate transcript quantification of RNA-seq data of five barley organs/tissues. This analysis identified 20,972 significant differentially expressed genes, 2791 differentially alternatively spliced genes and 2768 transcripts with differential transcript usage. CONCLUSION: A high confidence barley reference transcript dataset consisting of 60,444 genes with 177,240 transcripts has been generated. Compared to current barley transcripts, BaRTv1.0 transcripts are generally longer, have less fragmentation and improved gene models that are well supported by splice junction reads. Precise transcript quantification using BaRTv1.0 allows routine analysis of gene expression and AS.

Exome Capture for Variant Discovery and Analysis in Barley.

Exome capture is a reduced representation approach that selectively captures sequence from only the gene-bearing regions of a genome. It is based on probes targeted at these regions and, compared with whole genome shotgun sequencing, leads to a significant reduction in cost and data processing effort while still providing insights into the most relevant part of a genome. An exome capture array for barley was released in 2013 and this has opened the door to numerous studies that have put this technology to good use. In this chapter we detail the laboratory protocols required for enrichment and sequencing, and provide detailed step-by-step instructions for the bioinformatics analysis of the resulting data.

Characterisation of barley resistance to rhynchosporium on chromosome 6HS.

KEY MESSAGE: Major resistance gene to rhynchosporium, Rrs18, maps close to the telomere on the short arm of chromosome 6H in barley. Rhynchosporium or barley scald caused by a fungal pathogen Rhynchosporium commune is one of the most destructive and economically important diseases of barley in the world. Testing of Steptoe × Morex and CIho 3515 × Alexis doubled haploid populations has revealed a large effect QTL for resistance to R. commune close to the telomere on the short arm of chromosome 6H, present in both populations. Mapping markers flanking the QTL from both populations onto the 2017 Morex genome assembly revealed a rhynchosporium resistance locus independent of Rrs13 that we named Rrs18. The causal gene was fine mapped to an interval of 660 Kb using Steptoe × Morex backcross 1 S2 and S3 lines with molecular markers developed from Steptoe exome capture variant calling. Sequencing RNA from CIho 3515 and Alexis revealed that only 4 genes within the Rrs18 interval were transcribed in leaf tissue with a serine/threonine protein kinase being the most likely candidate for Rrs18.

Identification and rapid mapping of a gene conferring broad-spectrum late blight resistance in the diploid potato species Solanum verrucosum through DNA capture technologies.

KEY MESSAGE: A broad-spectrum late blight disease-resistance gene from Solanum verrucosum has been mapped to potato chromosome 9. The gene is distinct from previously identified-resistance genes. We have identified and characterised a broad-spectrum resistance to Phytophthora infestans from the wild Mexican species Solanum verrucosum. Diagnostic resistance gene enrichment (dRenSeq) revealed that the resistance is not conferred by previously identified nucleotide-binding, leucine-rich repeat genes. Utilising the sequenced potato genome as a reference, two complementary enrichment strategies that target resistance genes (RenSeq) and single/low-copy number genes (Generic-mapping enrichment Sequencing; GenSeq), respectively, were deployed for the rapid, SNP-based mapping of the resistance through bulked-segregant analysis. Both approaches independently positioned the resistance, referred to as Rpi-ver1, to the distal end of potato chromosome 9. Stringent post-enrichment read filtering identified a total of 64 informative SNPs that corresponded to the expected ratio for significant polymorphisms in the parents as well as the bulks. Of these, 61 SNPs are located on potato chromosome 9 and reside within 27 individual genes, which in the sequenced potato clone DM locate to positions 45.9 to 60.9 Mb. RenSeq- and GenSeq-derived SNPs within the target region were converted into allele-specific PCR-based KASP markers and further defined the position of the resistance to a 4.3 Mb interval at the bottom end of chromosome 9 between positions 52.62-56.98 Mb.

ELIGULUM-A Regulates Lateral Branch and Leaf Development in Barley.

The shoot apical and axillary meristems control shoot development, effectively influencing lateral branch and leaf formation. The barley (Hordeum vulgare) uniculm2 (cul2) mutation blocks axillary meristem development, and mutant plants lack lateral branches (tillers) that normally develop from the crown. A genetic screen for cul2 suppressors recovered two recessive alleles of ELIGULUM-A (ELI-A) that partially rescued the cul2 tillering phenotype. Mutations in ELI-A produce shorter plants with fewer tillers and disrupt the leaf blade-sheath boundary, producing liguleless leaves and reduced secondary cell wall development in stems and leaves. ELI-A is predicted to encode an unannotated protein containing an RNaseH-like domain that is conserved in land plants. ELI-A transcripts accumulate at the preligule boundary, the developing ligule, leaf margins, cells destined to develop secondary cell walls, and cells surrounding leaf vascular bundles. Recent studies have identified regulatory similarities between boundary development in leaves and lateral organs. Interestingly, we observed ELI-A transcripts at the preligule boundary, suggesting that ELI-A contributes to boundary formation between the blade and sheath. However, we did not observe ELI-A transcripts at the axillary meristem boundary in leaf axils, suggesting that ELI-A is not involved in boundary development for axillary meristem development. Our results show that ELI-A contributes to leaf and lateral branch development by acting as a boundary gene during ligule development but not during lateral branch development.

Development and Evaluation of a Barley 50k iSelect SNP Array.

High-throughput genotyping arrays continue to be an attractive, cost-effective alternative to sequencing based approaches. We have developed a new 50k Illumina Infinium iSelect genotyping array for barley, a cereal crop species of major international importance. The majority of SNPs on the array have been extracted from variants called in exome capture data of a wide range of European barley germplasm. We used the recently published barley pseudomolecule assembly to map the exome capture data, which allowed us to generate markers with accurate physical positions and detailed gene annotation. Markers from an existing and widely used barley 9k Infinium iSelect array were carried over onto the 50k chip for backward compatibility. The array design featured 49,267 SNP markers that converted into 44,040 working assays, of which 43,461 were scorable in GenomeStudio. Of the working assays, 6,251 are from the 9k iSelect platform. We validated the SNPs by comparing the genotype calls from the new array to legacy datasets. Rates of agreement averaged 98.1 and 93.9% respectively for the legacy 9k iSelect SNP set (Comadran et al., 2012) and the exome capture SNPs. To test the utility of the 50k chip for genetic mapping, we genotyped a segregating population derived from a Golden Promise × Morex cross (Liu et al., 2014) and mapped over 14,000 SNPs to genetic positions which showed a near exact correspondence to their known physical positions. Manual adjustment of the cluster files used by the interpreting software for genotype scoring improved results substantially, but migration of cluster files between sites led to a deterioration of results, suggesting that local adjustment of cluster files is required on a site-per-site basis. Information relating to the markers on the chip is available online at https://ics.hutton.ac.uk/50k.

Association Mapping of Diastatic Power in UK Winter and Spring Barley by Exome Sequencing of Phenotypically Contrasting Variety Sets.

Diastatic Power (DP) is an important quality trait for malt used in adjunct brewing and distilling. Substantial genetic variation for DP exists within UK elite barley cultivars, but breeding progress has been slow due to the limited demand, compared to the overall barley market, and difficulties in assessing DP. Estimates of DP (taken from recommended and national list trials between 1994 and 2012) from a collection of UK elite winter and spring varieties were used to identify contrasting sets of high and low DP varieties. DNA samples were pooled within sets and exome capture sequencing performed. Allele frequency estimates of Single Nucleotide Polymorphisms (SNPs) identified from the sequencing were used to identify genomic locations associated with differences in DP. Individual genotypes were generated from a set of custom KASP assays, both within sets and in a wider germplasm collection, to validate allele frequency estimates and marker associations with DP. QTL identified regions previously linked to variation in DP as well as novel associations. QTL colocalised with a number of genes annotated as having a diastase related function. Results indicate that winter barley is more genetically diverse for genes influencing DP. The marker assays produced by this work represent a resource that is available for immediate use by barley breeders in the production of new high DP varieties.

Construction of a map-based reference genome sequence for barley, Hordeum vulgare L.

Barley (Hordeum vulgare L.) is a cereal grass mainly used as animal fodder and raw material for the malting industry. The map-based reference genome sequence of barley cv. 'Morex' was constructed by the International Barley Genome Sequencing Consortium (IBSC) using hierarchical shotgun sequencing. Here, we report the experimental and computational procedures to (i) sequence and assemble more than 80,000 bacterial artificial chromosome (BAC) clones along the minimum tiling path of a genome-wide physical map, (ii) find and validate overlaps between adjacent BACs, (iii) construct 4,265 non-redundant sequence scaffolds representing clusters of overlapping BACs, and (iv) order and orient these BAC clusters along the seven barley chromosomes using positional information provided by dense genetic maps, an optical map and chromosome conformation capture sequencing (Hi-C). Integrative access to these sequence and mapping resources is provided by the barley genome explorer (BARLEX).

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.

Exome sequencing of geographically diverse barley landraces and wild relatives gives insights into environmental adaptation.

After domestication, during a process of widespread range extension, barley adapted to a broad spectrum of agricultural environments. To explore how the barley genome responded to the environmental challenges it encountered, we sequenced the exomes of a collection of 267 georeferenced landraces and wild accessions. A combination of genome-wide analyses showed that patterns of variation have been strongly shaped by geography and that variant-by-environment associations for individual genes are prominent in our data set. We observed significant correlations of days to heading (flowering) and height with seasonal temperature and dryness variables in common garden experiments, suggesting that these traits were major drivers of environmental adaptation in the sampled germplasm. A detailed analysis of known flowering-associated genes showed that many contain extensive sequence variation and that patterns of single- and multiple-gene haplotypes exhibit strong geographical structuring. This variation appears to have substantially contributed to range-wide ecogeographical adaptation, but many factors key to regional success remain unidentified.

Tablet: Visualizing Next-Generation Sequence Assemblies and Mappings.

This chapter is designed to be a practical guide to using Tablet for the visualization of next/second-generation (NGS) sequencing data. NGS data is being produced more frequently and in greater data volumes every year. As such, it is increasingly important to have tools which enable biologists and bioinformaticians to understand and gain key insights into their data. Visualization can play a key role in the exploration of such data as well as aid in the visual validation of sequence assemblies and features such as single nucleotide polymorphisms (SNPs). We aim to show several use cases which demonstrate Tablet's ability to visually highlight various situations of interest which can arise in NGS data.

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.