Kmasker--a tool for in silico prediction of single-copy FISH probes for the large-genome species Hordeum vulgare

Kmasker plants - a tool for assessing complex sequence space in plant species

Many plant genomes display high levels of repetitive sequences. The assembly of these complex genomes using short high-throughput sequence reads is still a challenging task. Underestimation or disregard of repeat complexity in these datasets can easily misguide downstream analysis. Detection of repetitive regions by k-mer counting methods has proved to be reliable. Easy-to-use applications utilizing k-mer counting are in high demand, especially in the domain of plants. We present Kmasker plants, a tool that uses k-mer count information as an assistant throughout the analytical workflow of genome data that is provided as a command-line and web-based solution. Beside its core competence to screen and mask repetitive sequences, we have integrated features that enable comparative studies between different cultivars or closely related species and methods that estimate target specificity of guide RNAs for application of site-directed mutagenesis using Cas9 endonuclease. In addition, we have set up a web service for Kmasker plants that maintains pre-computed indices for 10 of the economically most important cultivated plants. Source code for Kmasker plants has been made publically available at https://github.com/tschmutzer/kmasker. The web service is accessible at https://kmasker.ipk-gatersleben.de.

Discovery of multi-megabase polymorphic inversions by chromosome conformation capture sequencing in large-genome plant species

Chromosomal inversions occur in natural populations of many species, and may underlie reproductive isolation and local adaptation. Traditional methods of inversion discovery are labor-intensive and lack sensitivity. Here, we report the use of three-dimensional contact probabilities between genomic loci as assayed by chromosome-conformation capture sequencing (Hi-C) to detect multi-megabase polymorphic inversions in four barley genotypes. Inversions are validated by fluorescence in situ hybridization and Bionano optical mapping. We propose Hi-C as a generally applicable method for inversion discovery in natural populations.

Evolution of Plant B Chromosome Enriched Sequences

B chromosomes are supernumerary chromosomes found in addition to the normal standard chromosomes (A chromosomes). B chromosomes are well known to accumulate several distinct types of repeated DNA elements. Although the evolution of B chromosomes has been the subject of numerous studies, the mechanisms of accumulation and evolution of repetitive sequences are not fully understood. Recently, new genomic approaches have shed light on the origin and accumulation of different classes of repetitive sequences in the process of B chromosome formation and evolution. Here we discuss the impact of repetitive sequences accumulation on the evolution of plant B chromosomes.

How Next-Generation Sequencing Has Aided Our Understanding of the Sequence Composition and Origin of B Chromosomes

Accessory, supernumerary, or-most simply-B chromosomes, are found in many eukaryotic karyotypes. These small chromosomes do not follow the usual pattern of segregation, but rather are transmitted in a higher than expected frequency. As increasingly being demonstrated by next-generation sequencing (NGS), their structure comprises fragments of standard (A) chromosomes, although in some plant species, their sequence also includes contributions from organellar genomes. Transcriptomic analyses of various animal and plant species have revealed that, contrary to what used to be the common belief, some of the B chromosome DNA is protein-encoding. This review summarizes the progress in understanding B chromosome biology enabled by the application of next-generation sequencing technology and state-of-the-art bioinformatics. In particular, a contrast is drawn between a direct sequencing approach and a strategy based on a comparative genomics as alternative routes that can be taken towards the identification of B chromosome sequences.

From plant genomes to phenotypes

Recent advances in sequencing technologies have greatly accelerated the rate of plant genome and applied breeding research. Despite this advancing trend, plant genomes continue to present numerous difficulties to the standard tools and pipelines not only for genome assembly but also gene annotation and downstream analysis. Here we give a perspective on tools, resources and services necessary to assemble and analyze plant genomes and link them to plant phenotypes.

Towards a whole-genome sequence for rye (Secale cereale L.)

We report on a whole-genome draft sequence of rye (Secale cereale L.). Rye is a diploid Triticeae species closely related to wheat and barley, and an important crop for food and feed in Central and Eastern Europe. Through whole-genome shotgun sequencing of the 7.9-Gbp genome of the winter rye inbred line Lo7 we obtained a de novo assembly represented by 1.29 million scaffolds covering a total length of 2.8 Gbp. Our reference sequence represents nearly the entire low-copy portion of the rye genome. This genome assembly was used to predict 27 784 rye gene models based on homology to sequenced grass genomes. Through resequencing of 10 rye inbred lines and one accession of the wild relative S. vavilovii, we discovered more than 90 million single nucleotide variants and short insertions/deletions in the rye genome. From these variants, we developed the high-density Rye600k genotyping array with 600 843 markers, which enabled anchoring the sequence contigs along a high-density genetic map and establishing a synteny-based virtual gene order. Genotyping data were used to characterize the diversity of rye breeding pools and genetic resources, and to obtain a genome-wide map of selection signals differentiating the divergent gene pools. This rye whole-genome sequence closes a gap in Triticeae genome research, and will be highly valuable for comparative genomics, functional studies and genome-based breeding in rye.

Multiplex sequencing of bacterial artificial chromosomes for assembling complex plant genomes

Hierarchical shotgun sequencing remains the method of choice for assembling high-quality reference sequences of complex plant genomes. The efficient exploitation of current high-throughput technologies and powerful computational facilities for large-insert clone sequencing necessitates the sequencing and assembly of a large number of clones in parallel. We developed a multiplexed pipeline for shotgun sequencing and assembling individual bacterial artificial chromosomes (BACs) using the Illumina sequencing platform. We illustrate our approach by sequencing 668 barley BACs (Hordeum vulgare L.) in a single Illumina HiSeq 2000 lane. Using a newly designed parallelized computational pipeline, we obtained sequence assemblies of individual BACs that consist, on average, of eight sequence scaffolds and represent >98% of the genomic inserts. Our BAC assemblies are clearly superior to a whole-genome shotgun assembly regarding contiguity, completeness and the representation of the gene space. Our methods may be employed to rapidly obtain high-quality assemblies of a large number of clones to assemble map-based reference sequences of plant and animal species with complex genomes by sequencing along a minimum tiling path.

transPLANT Resources for Triticeae Genomic Data

The genome sequences of many important Triticeae species, including bread wheat ( L.) and barley ( L.), remained uncharacterized for a long time because their high repeat content, large sizes, and polyploidy. As a result of improvements in sequencing technologies and novel analyses strategies, several of these have recently been deciphered. These efforts have generated new insights into Triticeae biology and genome organization and have important implications for downstream usage by breeders, experimental biologists, and comparative genomicists. transPLANT () is an EU-funded project aimed at constructing hardware, software, and data infrastructure for genome-scale research in the life sciences. Since the Triticeae data are intrinsically complex, heterogenous, and distributed, the transPLANT consortium has undertaken efforts to develop common data formats and tools that enable the exchange and integration of data from distributed resources. Here we present an overview of the individual Triticeae genome resources hosted by transPLANT partners, introduce the objectives of transPLANT, and outline common developments and interfaces supporting integrated data access.

Species-wide genome sequence and nucleotide polymorphisms from the model allopolyploid plant Brassica napus

Brassica napus (oilseed rape, canola) is one of the world's most important sources of vegetable oil for human nutrition and biofuel, and also a model species for studies investigating the evolutionary consequences of polyploidisation. Strong bottlenecks during its recent origin from interspecific hybridisation, and subsequently through intensive artificial selection, have severely depleted the genetic diversity available for breeding. On the other hand, high-throughput genome profiling technologies today provide unprecedented scope to identify, characterise and utilise genetic diversity in primary and secondary crop gene pools. Such methods also enable implementation of genomic selection strategies to accelerate breeding progress. The key prerequisite is availability of high-quality sequence data and identification of high-quality, genome-wide sequence polymorphisms representing relevant gene pools. We present comprehensive genome resequencing data from a panel of 52 highly diverse natural and synthetic B. napus accessions, along with a stringently selected panel of 4.3 million high-confidence, genome-wide SNPs. The data is of great interest for genomics-assisted breeding and for evolutionary studies on the origins and consequences in allopolyploidisation in plants.

Evolution of the Grain Dispersal System in Barley

About 12,000 years ago in the Near East, humans began the transition from hunter-gathering to agriculture-based societies. Barley was a founder crop in this process, and the most important steps in its domestication were mutations in two adjacent, dominant, and complementary genes, through which grains were retained on the inflorescence at maturity, enabling effective harvesting. Independent recessive mutations in each of these genes caused cell wall thickening in a highly specific grain "disarticulation zone," converting the brittle floral axis (the rachis) of the wild-type into a tough, non-brittle form that promoted grain retention. By tracing the evolutionary history of allelic variation in both genes, we conclude that spatially and temporally independent selections of germplasm with a non-brittle rachis were made during the domestication of barley by farmers in the southern and northern regions of the Levant, actions that made a major contribution to the emergence of early agrarian societies.

Cytogenetic mapping with centromeric bacterial artificial chromosomes contigs shows that this recombination-poor region comprises more than half of barley chromosome 3H

Genetic maps are based on the frequency of recombination and often show different positions of molecular markers in comparison to physical maps, particularly in the centromere that is generally poor in meiotic recombinations. To decipher the position and order of DNA sequences genetically mapped to the centromere of barley (Hordeum vulgare) chromosome 3H, fluorescence in situ hybridization with mitotic metaphase and meiotic pachytene chromosomes was performed with 70 genomic single-copy probes derived from 65 fingerprinted bacterial artificial chromosomes (BAC) contigs genetically assigned to this recombination cold spot. The total physical distribution of the centromeric 5.5 cM bin of 3H comprises 58% of the mitotic metaphase chromosome length. Mitotic and meiotic chromatin of this recombination-poor region is preferentially marked by a heterochromatin-typical histone mark (H3K9me2), while recombination enriched subterminal chromosome regions are enriched in euchromatin-typical histone marks (H3K4me2, H3K4me3, H3K27me3) suggesting that the meiotic recombination rate could be influenced by the chromatin landscape.

Chromosome-Specific Painting in Cucumis Species Using Bulked Oligonucleotides

Chromosome-specific painting is a powerful technique in molecular cytogenetic and genome research. We developed an oligonucleotide (oligo)-based chromosome painting technique in cucumber (Cucumis sativus) that will be applicable in any plant species with a sequenced genome. Oligos specific to a single chromosome of cucumber were identified using a newly developed bioinformatic pipeline and then massively synthesized de novo in parallel. The synthesized oligos were amplified and labeled with biotin or digoxigenin for use in fluorescence in situ hybridization (FISH). We developed three different probes with each containing 23,000-27,000 oligos. These probes spanned 8.3-17 Mb of DNA on targeted cucumber chromosomes and had the densities of 1.5-3.2 oligos per kilobases. These probes produced FISH signals on a single cucumber chromosome and were used to paint homeologous chromosomes in other Cucumis species diverged from cucumber for up to 12 million years. The bulked oligo probes allowed us to track a single chromosome in early stages during meiosis. We were able to precisely map the pairing between cucumber chromosome 7 and chromosome 1 of Cucumis hystrix in a F1 hybrid. These two homeologous chromosomes paired in 71% of prophase I cells but only 25% of metaphase I cells, which may provide an explanation of the higher recombination rates compared to the chiasma frequencies between homeologous chromosomes reported in plant hybrids.

Mapping-by-sequencing accelerates forward genetics in barley

Mapping-by-sequencing has emerged as a powerful technique for genetic mapping in several plant and animal species. As this resequencing-based method requires a reference genome, its application to complex plant genomes with incomplete and fragmented sequence resources remains challenging. We perform exome sequencing of phenotypic bulks of a mapping population of barley segregating for a mutant phenotype that increases the rate of leaf initiation. Read depth analysis identifies a candidate gene, which is confirmed by the analysis of independent mutant alleles. Our method illustrates how the genomic resources of barley together with exome resequencing can underpin mapping-by-sequencing.

Mapping-by-sequencing accelerates forward genetics in barley

Mapping-by-sequencing has emerged as a powerful technique for genetic mapping in several plant and animal species. As this resequencing-based method requires a reference genome, its application to complex plant genomes with incomplete and fragmented sequence resources remains challenging. We perform exome sequencing of phenotypic bulks of a mapping population of barley segregating for a mutant phenotype that increases the rate of leaf initiation. Read depth analysis identifies a candidate gene, which is confirmed by the analysis of independent mutant alleles. Our method illustrates how the genomic resources of barley together with exome resequencing can underpin mapping-by-sequencing.

PROTEIN DISULFIDE ISOMERASE LIKE 5-1 is a susceptibility factor to plant viruses

Protein disulfide isomerases (PDIs) catalyze the correct folding of proteins and prevent the aggregation of unfolded or partially folded precursors. Whereas suppression of members of the PDI gene family can delay replication of several human and animal viruses (e.g., HIV), their role in interactions with plant viruses is largely unknown. Here, using a positional cloning strategy we identified variants of PROTEIN DISULFIDE ISOMERASE LIKE 5-1 (HvPDIL5-1) as the cause of naturally occurring resistance to multiple strains of Bymoviruses. The role of wild-type HvPDIL5-1 in conferring susceptibility was confirmed by targeting induced local lesions in genomes for induced mutant alleles, transgene-induced complementation, and allelism tests using different natural resistance alleles. The geographical distribution of natural genetic variants of HvPDIL5-1 revealed the origin of resistance conferring alleles in domesticated barley in Eastern Asia. Higher sequence diversity was correlated with areas with increased pathogen diversity suggesting adaptive selection for bymovirus resistance. HvPDIL5-1 homologs are highly conserved across species of the plant and animal kingdoms implying that orthologs of HvPDIL5-1 or other closely related members of the PDI gene family may be potential susceptibility factors to viruses in other eukaryotic species.