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.

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.

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.