A sequence-ready physical map of barley anchored genetically by two million single-nucleotide polymorphisms

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Amand P, Diggle AJ, Rengel Z. Genetic characterization of root architectural traits in barley (Hordeum vulgare L.) using SNP markers

Increasing attention is paid to providing new tools to breeders for targeted breeding for specific root traits that are beneficial in low-fertility, drying soils; however, such information is not available for barley (Hordeum vulgare L.). A panel of 191 barley accessions (originating from Australia, Europe, and Africa) was phenotyped for 26 root and shoot traits using the semi-hydroponic system and genotyped using 21 062 high-quality single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing (GBS). The population structure analysis of the barley panel identified six distinct groups. We detected 1199 significant (P<0.001) marker-trait associations (MTAs) with r2 values up to 0.41. The strongest MTAs were found for root diameter in the top 20 cm and the longest root length. Based on the physical locations of these MTAs in the barley reference genome, we identified 37 putative QTLs for the root traits, and three QTLs for shoot traits, with nine QTLs located in the same physical regions. The genomic region 640-653 Mb on chromosome 7H was significant for five root length-related traits, where 440 annotated genes were located. The putative QTLs for various root traits identified in this study may be useful for genetic improvement regarding the adaptation of new barley cultivars to suboptimal environments and abiotic stresses.

Scripting Analyses of Genomes in Ensembl Plants

Ensembl Plants ( http://plants.ensembl.org ) offers genome-scale information for plants, with four releases per year. As of release 47 (April 2020) it features 79 species and includes genome sequence, gene models, and functional annotation. Comparative analyses help reconstruct the evolutionary history of gene families, genomes, and components of polyploid genomes. Some species have gene expression baseline reports or variation across genotypes. While the data can be accessed through the Ensembl genome browser, here we review specifically how our plant genomes can be interrogated programmatically and the data downloaded in bulk. These access routes are generally consistent across Ensembl for other non-plant species, including plant pathogens, pests, and pollinators.

A high-resolution consensus linkage map for barley based on GBS-derived genotypes

As genotyping-by-sequencing (GBS) is widely used in barley genetic studies, the translation of the physical position of GBS-derived SNPs into accurate genetic positions has become relevant. The main aim of this study was to develop a high-resolution consensus linkage map based on GBS-derived SNPs. The construction of this integrated map involved 11 bi-parental populations composed of 3743 segregating progenies. We adopted a uniform set of SNP-calling and filtering conditions to identify 50 875 distinct SNPs segregating in at least one population. These SNPs were grouped into 18 580 non-redundant SNPs (bins). The resulting consensus linkage map spanned 1050.1 cM, providing an average density of 17.7 bins and 48.4 SNPs per cM. The consensus map is characterized by the absence of large intervals devoid of marker coverage (significant gaps), the largest interval between bins was only 3.7 cM and the mean distance between adjacent bins was 0.06 cM. This high-resolution linkage map will contribute to several applications in genomic research, such as providing useful information on the recombination landscape for QTLs/genes identified via GWAS or ensuring a uniform distribution of SNPs when developing low-cost genotyping tools offering a limited number of markers.

A linkage disequilibrium-based approach to position unmapped SNPs in crop species

BACKGROUND

High-density SNP arrays are now available for a wide range of crop species. Despite the development of many tools for generating genetic maps, the genome position of many SNPs from these arrays is unknown. Here we propose a linkage disequilibrium (LD)-based algorithm to allocate unassigned SNPs to chromosome regions from sparse genetic maps. This algorithm was tested on sugarcane, wheat, and barley data sets. We calculated the algorithm's efficiency by masking SNPs with known locations, then assigning their position to the map with the algorithm, and finally comparing the assigned and true positions.

RESULTS

In the 20-fold cross-validation, the mean proportion of masked mapped SNPs that were placed by the algorithm to a chromosome was 89.53, 94.25, and 97.23% for sugarcane, wheat, and barley, respectively. Of the markers that were placed in the genome, 98.73, 96.45 and 98.53% of the SNPs were positioned on the correct chromosome. The mean correlations between known and new estimated SNP positions were 0.97, 0.98, and 0.97 for sugarcane, wheat, and barley. The LD-based algorithm was used to assign 5920 out of 21,251 unpositioned markers to the current Q208 sugarcane genetic map, representing the highest density genetic map for this species to date.

CONCLUSIONS

Our LD-based approach can be used to accurately assign unpositioned SNPs to existing genetic maps, improving genome-wide association studies and genomic prediction in crop species with fragmented and incomplete genome assemblies. This approach will facilitate genomic-assisted breeding for many orphan crops that lack genetic and genomic resources.

Efficient construction of a linkage map and haplotypes for Mentha suaveolens using sequence capture

The sustainability of many crops is hindered by the lack of genomic resources and a poor understanding of natural genetic diversity. Particularly, application of modern breeding requires high-density linkage maps that are integrated into a highly contiguous reference genome. Here, we present a rapid method for deriving haplotypes and developing linkage maps, and its application to Mentha suaveolens, one of the diploid progenitors of cultivated mints. Using sequence-capture via DNA hybridization to target single nucleotide polymorphisms (SNPs), we successfully genotyped ∼5000 SNPs within the genome of >400 individuals derived from a self cross. After stringent quality control, and identification of nonredundant SNPs, 1919 informative SNPs were retained for linkage map construction. The resulting linkage map defined a total genetic space of 942.17 cM divided among 12 linkage groups, ranging from 56.32 to 122.61 cM in length. The linkage map is in good agreement with pseudomolecules from our preliminary genome assembly, proving this resource effective for the correction and validation of the reference genome. We discuss the advantages of this method for the rapid creation of linkage maps.

Long-read sequence assembly: a technical evaluation in barley

Sequence assembly of large and repeat-rich plant genomes has been challenging, requiring substantial computational resources and often several complementary sequence assembly and genome mapping approaches. The recent development of fast and accurate long-read sequencing by circular consensus sequencing (CCS) on the PacBio platform may greatly increase the scope of plant pan-genome projects. Here, we compare current long-read sequencing platforms regarding their ability to rapidly generate contiguous sequence assemblies in pan-genome studies of barley (Hordeum vulgare). Most long-read assemblies are clearly superior to the current barley reference sequence based on short-reads. Assemblies derived from accurate long reads excel in most metrics, but the CCS approach was the most cost-effective strategy for assembling tens of barley genomes. A downsampling analysis indicated that 20-fold CCS coverage can yield very good sequence assemblies, while even five-fold CCS data may capture the complete sequence of most genes. We present an updated reference genome assembly for barley with near-complete representation of the repeat-rich intergenic space. Long-read assembly can underpin the construction of accurate and complete sequences of multiple genomes of a species to build pan-genome infrastructures in Triticeae crops and their wild relatives.

[NGS sequencing in barley breeding and genetic studies]

Barley (Hordeum vulgare L.) is the one of the most important cereal species used as food and feed crops, as well as for malting and alcohol production. At the end of the last century, traditional breeding techniques were complemented by the use of DNA markers. Molecular markers have also been used extensively for molecular genetic mapping and QTL analysis. In 2012, the barley genome sequencing was completed, which provided a broad range of new opportunities - from a more efficient search for candidate genes controlling economically important traits to genomic selection. The review summarizes the results of the studies performed after barley genome sequencing, which discovered new areas of barley genetics and breeding with high throughput screening and genotyping methods. During this period, intensive studies aimed at identification of barley genomic loci associated with economically important traits have been carried out; online databases and tools for working with barley genomic data and their deposition have appeared and are being replenished. In recent years, GWAS analysis has been used for large-scale phenotypegenotype association studies, which has been widely used in barley since 2010 due to the developed SNP-arrays, as well as genotyping methods based on direct NGS sequencing of selected fractions of the genome. To date, more than 80 papers have been published that describe the results of the GWAS analysis in barley. SNP identification associated with economically important traits and their transformation into CAPS or KASP markers convenient for screening selection material significantly expands the possibilities of marker-assisted selection of barley. In addition, the currently available information on potential target genes and the quality of the whole barley genome sequence provides a good base for applying genome editing technologies to create material for the creation of varieties with desired properties.

Gramene 2021: harnessing the power of comparative genomics and pathways for plant research

Gramene (http://www.gramene.org), a knowledgebase founded on comparative functional analyses of genomic and pathway data for model plants and major crops, supports agricultural researchers worldwide. The resource is committed to open access and reproducible science based on the FAIR data principles. Since the last NAR update, we made nine releases; doubled the genome portal's content; expanded curated genes, pathways and expression sets; and implemented the Domain Informational Vocabulary Extraction (DIVE) algorithm for extracting gene function information from publications. The current release, #63 (October 2020), hosts 93 reference genomes-over 3.9 million genes in 122 947 families with orthologous and paralogous classifications. Plant Reactome portrays pathway networks using a combination of manual biocuration in rice (320 reference pathways) and orthology-based projections to 106 species. The Reactome platform facilitates comparison between reference and projected pathways, gene expression analyses and overlays of gene-gene interactions. Gramene integrates ontology-based protein structure-function annotation; information on genetic, epigenetic, expression, and phenotypic diversity; and gene functional annotations extracted from plant-focused journals using DIVE. We train plant researchers in biocuration of genes and pathways; host curated maize gene structures as tracks in the maize genome browser; and integrate curated rice genes and pathways in the Plant Reactome.

Gramene 2021: harnessing the power of comparative genomics and pathways for plant research

Gramene (http://www.gramene.org), a knowledgebase founded on comparative functional analyses of genomic and pathway data for model plants and major crops, supports agricultural researchers worldwide. The resource is committed to open access and reproducible science based on the FAIR data principles. Since the last NAR update, we made nine releases; doubled the genome portal's content; expanded curated genes, pathways and expression sets; and implemented the Domain Informational Vocabulary Extraction (DIVE) algorithm for extracting gene function information from publications. The current release, #63 (October 2020), hosts 93 reference genomes-over 3.9 million genes in 122 947 families with orthologous and paralogous classifications. Plant Reactome portrays pathway networks using a combination of manual biocuration in rice (320 reference pathways) and orthology-based projections to 106 species. The Reactome platform facilitates comparison between reference and projected pathways, gene expression analyses and overlays of gene-gene interactions. Gramene integrates ontology-based protein structure-function annotation; information on genetic, epigenetic, expression, and phenotypic diversity; and gene functional annotations extracted from plant-focused journals using DIVE. We train plant researchers in biocuration of genes and pathways; host curated maize gene structures as tracks in the maize genome browser; and integrate curated rice genes and pathways in the Plant Reactome.

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

Diagnostic markers for Rrs1_Rh4 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 Rrs1_Rh4. These markers are specific to Rrs1_Rh4 and do not detect other Rrs1 genes/alleles. The Rrs1_Rh4 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 Rrs1_Rh4. One of these lines came from Jordan, with the remaining 12 lines from different locations in Syria. One of the Syrian landraces containing Rrs1_Rh4 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.

Pairing and Exchanging between Daypyrum villosum Chromosomes 6V#2 and 6V#4 in the Hybrids of Two Different Wheat Alien Substitution Lines

Normal pairing and exchanging is an important basis to evaluate the genetic relationship between homologous chromosomes in a wheat background. The pairing behavior between 6V#2 and 6V#4, two chromosomes from different Dasypyrum villosum accessions, is still not clear. In this study, two wheat alien substitution lines, 6V#2 (6A) and 6V#4 (6D), were crossed to obtain the F₁ hybrids and F₂ segregating populations, and the testcross populations were obtained by using the F₁ as a parent crossed with wheat variety Wan7107. The chromosomal behavior at meiosis in pollen mother cells (PMCs) of the F₁ hybrids was observed using a genomic in situ hybridization (GISH) technique. Exchange events of two alien chromosomes were investigated in the F₂ populations using nine polymerase chain reaction (PCR) markers located on the 6V short arm. The results showed that the two alien chromosomes could pair with each other to form ring- or rod-shaped bivalent chromosomes in 79.76% of the total PMCs, and most were pulled to two poles evenly at anaphase I. Investigation of the F₂ populations showed that the segregation ratios of seven markers were consistent with the theoretical values 3:1 or 1:2:1, and recombinants among markers were detected. A genetic linkage map of nine PCR markers for 6VS was accordingly constructed based on the exchange frequencies and compared with the physical maps of wheat and barley based on homologous sequences of the markers, which showed that conservation of sequence order compared to 6V was 6H and 6B > 6A > 6D. In the testcross populations with 482 plants, seven showed susceptibility to powdery mildew (PM) and lacked amplification of alien chromosomal bands. Six other plants had amplification of specific bands of both the alien chromosomes at multiple sites, which suggested that the alien chromosomes had abnormal separation behavior in about 1.5% of the PMCs in F₁, which resulted in some gametes containing two alien chromosomes. In addition, three new types of chromosome substitution were developed. This study lays a foundation for alien allelism tests and further assessment of the genetic relationship among 6V#2, 6V#4, and their wheat homoeologous chromosomes.

Construction of Genetic Linkage Maps in Multiparental Populations

Construction of genetic linkage maps has become a routine step for mapping quantitative trait loci (QTL), particularly in animal and plant breeding populations. Many multiparental populations have recently been produced to increase genetic diversity and QTL mapping resolution. However, few software packages are available for map construction in these populations. In this paper, we build a general framework for the construction of genetic linkage maps from genotypic data in diploid populations, including bi- and multiparental populations, cross-pollinated (CP) populations, and breeding pedigrees. The framework is implemented as an automatic pipeline called magicMap, where the maximum multilocus likelihood approach utilizes genotypic information efficiently. We evaluate magicMap by extensive simulations and eight real datasets: one biparental, one CP, four multiparent advanced generation intercross (MAGIC), and two nested association mapping (NAM) populations, the number of markers ranging from a few hundred to tens of thousands. Not only is magicMap the only software capable of accommodating all of these designs, it is more accurate and robust to missing genotypes and genotyping errors than commonly used packages.

Leaf Variegation and Impaired Chloroplast Development Caused by a Truncated CCT Domain Gene in albostrians Barley

Chloroplasts fuel plant development and growth by converting solar energy into chemical energy. They mature from proplastids through the concerted action of genes in both the organellar and the nuclear genome. Defects in such genes impair chloroplast development and may lead to pigment-deficient seedlings or seedlings with variegated leaves. Such mutants are instrumental as tools for dissecting genetic factors underlying the mechanisms involved in chloroplast biogenesis. Characterization of the green-white variegated albostrians mutant of barley (Hordeum vulgare) has greatly broadened the field of chloroplast biology, including the discovery of retrograde signaling. Here, we report identification of the ALBOSTRIANS gene HvAST (also known as Hordeum vulgare CCT Motif Family gene 7, HvCMF7) by positional cloning as well as its functional validation based on independently induced mutants by Targeting Induced Local Lesions in Genomes (TILLING) and RNA-guided clustered regularly interspaced short palindromic repeats-associated protein 9 endonuclease-mediated gene editing. The phenotypes of the independent HvAST mutants imply residual activity of HvCMF7 in the original albostrians allele conferring an imperfect penetrance of the variegated phenotype even at homozygous state of the mutation. HvCMF7 is a homolog of the Arabidopsis (Arabidopsis thaliana) CONSTANS, CO-like, and TOC1 (CCT) Motif transcription factor gene CHLOROPLAST IMPORT APPARATUS2, which was reported to be involved in the expression of nuclear genes essential for chloroplast biogenesis. Notably, in barley we localized HvCMF7 to the chloroplast, without any clear evidence for nuclear localization.

Integrating a newly developed BAC-based physical mapping resource for Lolium perenne with a genome-wide association study across a L. perenne European ecotype collection identifies genomic contexts associated with agriculturally important traits

BACKGROUND AND AIMS

Lolium perenne (perennial ryegrass) is the most widely cultivated forage and amenity grass species in temperate areas worldwide and there is a need to understand the genetic architectures of key agricultural traits and crop characteristics that deliver wider environmental services. Our aim was to identify genomic regions associated with agriculturally important traits by integrating a bacterial artificial chromosome (BAC)-based physical map with a genome-wide association study (GWAS).

METHODS

BAC-based physical maps for L. perenne were constructed from ~212 000 high-information-content fingerprints using Fingerprint Contig and Linear Topology Contig software. BAC clones were associated with both BAC-end sequences and a partial minimum tiling path sequence. A panel of 716 L. perenne diploid genotypes from 90 European accessions was assessed in the field over 2 years, and genotyped using a Lolium Infinium SNP array. The GWAS was carried out using a linear mixed model implemented in TASSEL, and extended genomic regions associated with significant markers were identified through integration with the physical map.

KEY RESULTS

Between ~3600 and 7500 physical map contigs were derived, depending on the software and probability thresholds used, and integrated with ~35 k sequenced BAC clones to develop a resource predicted to span the majority of the L. perenne genome. From the GWAS, eight different loci were significantly associated with heading date, plant width, plant biomass and water-soluble carbohydrate accumulation, seven of which could be associated with physical map contigs. This allowed the identification of a number of candidate genes.

CONCLUSIONS

Combining the physical mapping resource with the GWAS has allowed us to extend the search for candidate genes across larger regions of the L. perenne genome and identified a number of interesting gene model annotations. These physical maps will aid in validating future sequence-based assemblies of the L. perenne genome.

Resistance to Rhynchosporium commune in a collection of European spring barley germplasm

Association analyses of resistance to Rhynchosporium commune in a collection of European spring barley germplasm detected 17 significant resistance quantitative trait loci. The most significant association was confirmed as Rrs1. Rhynchosporium commune is a fungal pathogen of barley which causes a highly destructive and economically important disease known as rhynchosporium. Genome-wide association mapping was used to investigate the genetic control of host resistance to R. commune in a collection of predominantly European spring barley accessions. Multi-year disease nursery field trials revealed 8 significant resistance quantitative trait loci (QTL), whilst a separate association mapping analysis using historical data from UK national and recommended list trials identified 9 significant associations. The most significant association identified in both current and historical data sources, collocated with the known position of the major resistance gene Rrs1. Seedling assays with R. commune single-spore isolates expressing the corresponding avirulence protein NIP1 confirmed that this locus is Rrs1. These results highlight the significant and continuing contribution of Rrs1 to host resistance in current elite spring barley germplasm. Varietal height was shown to be negatively correlated with disease severity, and a resistance QTL was identified that co-localised with the semi-dwarfing gene sdw1, previously shown to contribute to disease escape. The remaining QTL represent novel resistances that are present within European spring barley accessions. Associated markers to Rrs1 and other resistance loci, identified in this study, represent a set of tools that can be exploited by breeders for the sustainable deployment of varietal resistance in new cultivars.

Construction of A High-Density Genetic Map and Mapping of Fruit Traits in Watermelon (Citrullus Lanatus L.) Based on Whole-Genome Resequencing

Watermelon (Citrullus lanatus L.) is an important horticultural crop that is grown worldwide and has a high economic value. To dissect the loci associated with important horticultural traits and to analyze the genetic and genomic information of this species, a high-density genetic map was constructed based on whole-genome resequencing (WGR), a powerful high-resolution method for single-nucleotide polymorphism (SNP) marker development, genetic map construction, and gene mapping. Resequencing of both parental lines and 126 recombinant inbred lines (RIL) resulted in the detection of 178,762 single-nucleotide polymorphism (SNP) markers in the parental lines at a sequencing depth greater than four-fold. Additionally, 2132 recombination bin markers comprising 103,029 SNP markers were mapped onto 11 linkage groups (LGs). Substantially more SNP markers were mapped to the genetic map compared with other recent studies. The total length of the linkage map was 1508.94 cM, with an average distance of 0.74 cM between adjacent bin markers. Based on this genetic map, one locus for fruit bitterness, one locus for rind color, and one locus for seed coat color with high LOD scores (58.361, 18.353, 26.852) were identified on chromosome 1, chromosome 8, and chromosome 3, respectively. These prominent loci were identified in a region of 6.16 Mb, 2.07 Mb, and 0.37 Mb, respectively. On the basis of current research, the high-density map and mapping results will provide a valuable tool for identifying candidate genes, map-based gene cloning, comparative mapping, and marker-assisted selection (MAS) in watermelon breeding.

Genomic approaches for studying crop evolution

Understanding how crop plants evolved from their wild relatives and spread around the world can inform about the origins of agriculture. Here, we review how the rapid development of genomic resources and tools has made it possible to conduct genetic mapping and population genetic studies to unravel the molecular underpinnings of domestication and crop evolution in diverse crop species. We propose three future avenues for the study of crop evolution: establishment of high-quality reference genomes for crops and their wild relatives; genomic characterization of germplasm collections; and the adoption of novel methodologies such as archaeogenetics, epigenomics, and genome editing.

QTL Mapping for Fiber Quality and Yield Traits Based on Introgression Lines Derived from Gossypium hirsutum × G. tomentosum

The tetraploid species Gossypium hirsutum is cultivated widely throughout the world with high yield and moderate fiber quality, but its genetic basis is narrow. A set of 107 introgression lines (ILs) was developed with an interspecific cross using G. hirsutumacc. 4105 as the recurrent parent and G. tomentosum as the donor parent. A specific locus amplified fragment sequencing (SLAF-seq) strategy was used to obtain high-throughput single nucleotide polymorphism (SNP) markers. In total, 3157 high-quality SNP markers were obtained and further used for identification of quantitative trait loci (QTLs) for fiber quality and yield traits evaluated in multiple environments. In total, 74 QTLs were detected that were associated with five fiber quality traits (30 QTLs) and eight yield traits (44 QTLs), with 2.02-30.15% of the phenotypic variance explained (PVE), and 69 markers were found to be associated with these thirteen traits. Eleven chromosomes in the A sub-genome (At) harbored 47 QTLs, and nine chromosomes in the D sub-genome (Dt) harbored 27 QTLs. More than half (44 QTLs = 59.45%) showed positive additive effects for fiber and yield traits. Five QTL clusters were identified, with three in the At, comprised of thirteen QTLs, and two in the Dt comprised of seven QTLs. The ILs developed in this study and the identified QTLs will facilitate further molecular breeding for improvement of Upland cotton in terms of higher yield with enhanced fiber quality.

Fine mapping of the restorer gene Rfp3 from an Iranian primitive rye (Secale cereale L.)

A comparative genetics approach allowed to precisely determine the map position of the restorer gene Rfp3 in rye and revealed that Rfp3 and the restorer gene Rfm1 in barley reside at different positions in a syntenic 4RL/6HS segment. Cytoplasmic male sterility (CMS) is a reliable and striking genetic mechanism for hybrid seed production. Breeding of CMS-based hybrids in cereals requires the use of effective restorer genes as an indispensable pre-requisite. We report on the fine mapping of a restorer gene for the Pampa cytoplasm in winter rye that has been tapped from the Iranian primitive rye population Altevogt 14160. For this purpose, we have mapped 41 gene-derived markers to a 38.8 cM segment in the distal part of the long arm of chromosome 4R, which carries the restorer gene. Male fertility restoration was comprehensively analyzed in progenies of crosses between a male-sterile tester genotype and 21 recombinant as well as six non-recombinant BC₄S₂ lines. This approach allowed us to validate the position of this restorer gene, which we have designated Rfp3, on chromosome 4RL. Rfp3 was mapped within a 2.5 cM interval and cosegregated with the EST-derived marker c28385. The gene-derived conserved ortholog set (COS) markers enabled us to investigate the orthology of restorer genes originating from different genetic resources of rye as well as barley. The observed localization of Rfp3 and Rfm1 in a syntenic 4RL/6HS segment asks for further efforts towards cloning of both restorer genes as an option to study the mechanisms of male sterility and fertility restoration in cereals.

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).

Fine mapping of a major QTL for awn length in barley using a multiparent mapping population

Awn length was mapped using a multiparent population derived from cv. Morex and four wild accessions. One QTL was fine mapped and candidate genes were identified in NILs by RNA-seq. Barley awns are photosynthetically active and contribute to grain yield. Awn length is variable among both wild and cultivated barley genotypes and many mutants with alterations in awn length have been identified. Here, we used a multiparent mapping population derived from cv. Morex and four genetically diverse wild barley lines to detect quantitative trait loci (QTLs) for awn length. Twelve QTLs, distributed over the barley genome, were identified with the most significant one located on chromosome arm 7HL (QTL AL7.1). The effect of AL7.1 was confirmed using near isogenic lines (NILs) and fine-mapped in two independent heterogeneous inbred families to a < 0.9 cM interval. With exception of a small effect on grain width, no other traits such as plant height or flowering time were affected by AL7.1. Variant calling on transcripts obtained from RNA sequencing reads in NILs was used to narrow down the list of candidate genes located in the interval. This data may be used for further characterization and unravelling of the mechanisms underlying natural variation in awn length.

Identification of nuclear genes controlling chlorophyll synthesis in barley by RNA-seq

BACKGROUND

Albinism in plants is characterized by lack of chlorophyll and results in photosynthesis impairment, abnormal plant development and premature death. These abnormalities are frequently encountered in interspecific crosses and tissue culture experiments. Analysis of albino mutant phenotypes with full or partial chlorophyll deficiency can shed light on genetic determinants and molecular mechanisms of albinism. Here we report analysis of RNA-seq transcription profiling of barley (Hordeum vulgare L.) near-isogenic lines, one of which is a carrier of mutant allele of the Alm gene for albino lemma and pericarp phenotype (line i:BwAlm).

RESULTS

1221 genome fragments have statistically significant changes in expression levels between lines i:BwAlm and Bowman, with 148 fragments having increased expression levels in line i:BwAlm, and 1073 genome fragments, including 42 plastid operons, having decreased levels of expression in line i:BwAlm. We detected functional dissimilarity between genes with higher and lower levels of expression in i:BwAlm line. Genes with lower level of expression in the i:BwAlm line are mostly associated with photosynthesis and chlorophyll synthesis, while genes with higher expression level are functionally associated with vesicle transport. Differentially expressed genes are shown to be involved in several metabolic pathways; the largest fraction of such genes was observed for the Calvin-Benson-Bassham cycle. Finally, de novo assembly of transcriptome contains several transcripts, not annotated in current H. vulgare genome version.

CONCLUSIONS

Our results provide the new information about genes which could be involved in formation of albino lemma and pericarp phenotype. They demonstrate the interplay between nuclear and chloroplast genomes in this physiological process.

Marker development using SLAF-seq and whole-genome shotgun strategy to fine-map the semi-dwarf gene ari-e in barley

Background

Barley semi-dwarf genes have been extensively explored and widely used in barley breeding programs. The semi-dwarf gene ari-e from Golden Promise is an important gene associated with some agronomic traits and salt tolerance. While ari-e has been mapped on barley chromosome 5H using traditional markers and next-generation sequencing technologies, it has not yet been finely located on this chromosome.

Results

We integrated two methods to develop molecular markers for fine-mapping the semi-dwarf gene ari-e: (1) specific-length amplified fragment sequencing (SLAF-seq) with bulked segregant analysis (BSA) to develop SNP markers, and (2) the whole-genome shotgun sequence to develop InDels. Both SNP and InDel markers were developed in the target region and used for fine-mapping the ari-e gene. Linkage analysis showed that ari-e co-segregated with marker InDel-17 and was delimited by two markers (InDel-16 and DGSNP21) spanning 6.8 cM in the doubled haploid (DH) Dash × VB9104 population. The genetic position of ari-e was further confirmed in the Hindmarsh × W1 DH population which was located between InDel-7 and InDel-17. As a result, the overlapping region of the two mapping populations flanked by InDel-16 and InDel-17 was defined as the candidate region spanning 0.58 Mb on the POPSEQ physical map.

Conclusions

The current study demonstrated the SLAF-seq for SNP discovery and whole-genome shotgun sequencing for InDel development as an efficient approach to map complex genomic region for isolation of functional gene. The ari-e gene was fine mapped from 10 Mb to 0.58 Mb interval.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3247-4) contains supplementary material, which is available to authorized users.

Genetic Mapping of Millions of SNPs in Safflower (Carthamus tinctorius L.) via Whole-Genome Resequencing

Accurate assembly of complete genomes is facilitated by very high density genetic maps. We performed low-coverage, whole-genome shotgun sequencing on 96 F₆ recombinant inbred lines (RILs) of a cross between safflower (Carthamus tinctorius L.) and its wild progenitor (C. palaestinus Eig). We also produced a draft genome assembly of C. tinctorius covering 866 million bp (∼two-thirds) of the expected 1.35 Gbp genome after sequencing a single, short insert library to ∼21 × depth. Sequence reads from the RILs were mapped to this genome assembly to facilitate SNP identification, and the resulting polymorphisms were used to construct a genetic map. The resulting map included 2,008,196 genetically located SNPs in 1178 unique positions. A total of 57,270 scaffolds, each containing five or more mapped SNPs, were anchored to the map. This resulted in the assignment of sequence covering 14% of the expected genome length to a genetic position. Comparison of this safflower map to genetic maps of sunflower and lettuce revealed numerous chromosomal rearrangements, and the resulting patterns were consistent with a whole-genome duplication event in the lineage leading to sunflower. This sequence-based genetic map provides a powerful tool for the assembly of a low-cost draft genome of safflower, and the same general approach is expected to work for other species.

Fine Mapping of the Barley Chromosome 6H Net Form Net Blotch Susceptibility Locus

Net form net blotch, caused by the necrotrophic fungal pathogen Pyrenophora teres f. teres, is a destructive foliar disease of barley with the potential to cause significant yield loss in major production regions throughout the world. The complexity of the host-parasite genetic interactions in this pathosystem hinders the deployment of effective resistance in barley cultivars, warranting a deeper understanding of the interactions. Here, we report on the high-resolution mapping of the dominant susceptibility locus near the centromere of chromosome 6H in the barley cultivars Rika and Kombar, which are putatively targeted by necrotrophic effectors from P. teres f. teres isolates 6A and 15A, respectively. Utilization of progeny isolates derived from a cross of P. teres f. teres isolates 6A × 15A harboring single major virulence loci (VK1, VK2, and VR2) allowed for the Mendelization of single inverse gene-for-gene interactions in a high-resolution population consisting of 2976 Rika × Kombar recombinant gametes. Brachypodium distachyon synteny was exploited to develop and saturate the susceptibility region with markers, delimiting it to ∼0.24 cM and a partial physical map was constructed. This genetic and physical characterization further resolved the dominant susceptibility locus, designated Spt1 (susceptibility to P. teres f. teres). The high-resolution mapping and cosegregation of the Spt1.R and Spt1.K gene/s indicates tightly linked genes in repulsion or alleles possibly targeted by different necrotrophic effectors. Newly developed barley genomic resources greatly enhance the efficiency of positional cloning efforts in barley, as demonstrated by the Spt1 fine mapping and physical contig identification reported here.

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.

A Cluster of Nucleotide-Binding Site-Leucine-Rich Repeat Genes Resides in a Barley Powdery Mildew Resistance Quantitative Trait Loci on 7HL

Powdery mildew causes severe yield losses in barley production worldwide. Although many resistance genes have been described, only a few have already been cloned. A strong QTL (quantitative trait locus) conferring resistance to a wide array of powdery mildew isolates was identified in a Spanish barley landrace on the long arm of chromosome 7H. Previous studies narrowed down the QTL position, but were unable to identify candidate genes or physically locate the resistance. In this study, the exome of three recombinant lines from a high-resolution mapping population was sequenced and analyzed, narrowing the position of the resistance down to a single physical contig. Closer inspection of the region revealed a cluster of closely related NBS-LRR (nucleotide-binding site-leucine-rich repeat containing protein) genes. Large differences were found between the resistant lines and the reference genome of cultivar Morex, in the form of PAV (presence-absence variation) in the composition of the NBS-LRR cluster. Finally, a template-guided assembly was performed and subsequent expression analysis revealed that one of the new assembled candidate genes is transcribed. In summary, the results suggest that NBS-LRR genes, absent from the reference and the susceptible genotypes, could be functional and responsible for the powdery mildew resistance. The procedure followed is an example of the use of NGS (next-generation sequencing) tools to tackle the challenges of gene cloning when the target gene is absent from the reference genome.

A Homolog of Blade-On-Petiole 1 and 2 (BOP1/2) Controls Internode Length and Homeotic Changes of the Barley Inflorescence

Inflorescence architecture in small-grain cereals has a direct effect on yield and is an important selection target in breeding for yield improvement. We analyzed the recessive mutation laxatum-a (lax-a) in barley (Hordeum vulgare), which causes pleiotropic changes in spike development, resulting in (1) extended rachis internodes conferring a more relaxed inflorescence, (2) broadened base of the lemma awns, (3) thinner grains that are largely exposed due to reduced marginal growth of the palea and lemma, and (4) and homeotic conversion of lodicules into two stamenoid structures. Map-based cloning enforced by mapping-by-sequencing of the mutant lax-a locus enabled the identification of a homolog of BLADE-ON-PETIOLE1 (BOP1) and BOP2 as the causal gene. Interestingly, the recently identified barley uniculme4 gene also is a BOP1/2 homolog and has been shown to regulate tillering and leaf sheath development. While the Arabidopsis (Arabidopsis thaliana) BOP1 and BOP2 genes act redundantly, the barley genes contribute independent effects in specifying the developmental growth of vegetative and reproductive organs, respectively. Analysis of natural genetic diversity revealed strikingly different haplotype diversity for the two paralogous barley genes, likely affected by the respective genomic environments, since no indication for an active selection process was detected.

The big five of the monocot genomes

Monocots represent a monophyletic clade of the angiosperms that - based on fossil and molecular records - originated at around the Early Cretaceous from aquatic and wetland ancestors. Among their members are important crops including maize, wheat, rice, sorghum and barley, accounting for the major source for the daily calorie uptake by humans. Reflecting this importance, the partly large and complex genomes of these plants were major targets for ambitious and innovative sequencing projects, which will be discussed in this review article.

Isolation and fine mapping of Rps6: an intermediate host resistance gene in barley to wheat stripe rust

We uncouple host and nonhost resistance in barley to Puccinia striiformis ff. spp. hordei and tritici . We isolate, fine map, and physically anchor Rps6 to chromosome 7H in barley. A plant may be considered a nonhost of a pathogen if all known genotypes of a plant species are resistant to all known isolates of a pathogen species. However, if a small number of genotypes are susceptible to some known isolates of a pathogen species this plant may be considered an intermediate host. Barley (Hordeum vulgare) is an intermediate host for Puccinia striiformis f. sp. tritici (Pst), the causal agent of wheat stripe rust. We wanted to understand the genetic architecture underlying resistance to Pst and to determine whether any overlap exists with resistance to the host pathogen, Puccinia striiformis f. sp. hordei (Psh). We mapped Pst resistance to chromosome 7H and show that host and intermediate host resistance is genetically uncoupled. Therefore, we designate this resistance locus Rps6. We used phenotypic and genotypic selection on F2:3 families to isolate Rps6 and fine mapped the locus to a 0.1 cM region. Anchoring of the Rps6 locus to the barley physical map placed the region on a single fingerprinted contig spanning a physical region of 267 kb. Efforts are now underway to sequence the minimal tiling path and to delimit the physical region harboring Rps6. This will facilitate additional marker development and permit identification of candidate genes in the region.

Identification and Fine Mapping of a White Husk Gene in Barley (Hordeum vulgare L.)

Barley is the only crop in the Poaceae family with adhering husks at maturity. The color of husk at barely development stage could influence the agronomic traits and malting qualities of grains. A barley mutant with a white husk was discovered from the malting barley cultivar Supi 3 and designated wh (white husk). Morphological changes and the genetics of white husk barley were investigated. Husks of the mutant were white at the heading and flowering stages but yellowed at maturity. The diastatic power and α-amino nitrogen contents also significantly increased in wh mutant. Transmission electron microscopy examination showed abnormal chloroplast development in the mutant. Genetic analysis of F₂ and BC₁F₁ populations developed from a cross of wh and Yangnongpi 5 (green husk) showed that the white husk was controlled by a single recessive gene (wh). The wh gene was initially mapped between 49.64 and 51.77 cM on chromosome 3H, which is syntenic with rice chromosome 1 where a white husk gene wlp1 has been isolated. The barley orthologous gene of wlp1 was sequenced from both parents and a 688 bp deletion identified in the wh mutant. We further fine-mapped the wh gene between SSR markers Bmac0067 and Bmag0508a with distances of 0.36 cM and 0.27 cM in an F₂ population with 1115 individuals of white husk. However, the wlp1 orthologous gene was mapped outside the interval. New candidate genes were identified based on the barley genome sequence.

Mitogen-Activated Protein Kinase Kinase 3 Regulates Seed Dormancy in Barley

Seed dormancy has fundamental importance in plant survival and crop production; however, the mechanisms regulating dormancy remain unclear [1-3]. Seed dormancy levels generally decrease during domestication to ensure that crops successfully germinate in the field. However, reduction of seed dormancy can cause devastating losses in cereals like wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.) due to pre-harvest sprouting, the germination of mature seed (grain) on the mother plant when rain occurs before harvest. Understanding the mechanisms of dormancy can facilitate breeding of crop varieties with the appropriate levels of seed dormancy [4-8]. Barley is a model crop [9, 10] and has two major seed dormancy quantitative trait loci (QTLs), SD1 and SD2, on chromosome 5H [11-19]. We detected a QTL designated Qsd2-AK at SD2 as the single major determinant explaining the difference in seed dormancy between the dormant cultivar "Azumamugi" (Az) and the non-dormant cultivar "Kanto Nakate Gold" (KNG). Using map-based cloning, we identified the causal gene for Qsd2-AK as Mitogen-activated Protein Kinase Kinase 3 (MKK3). The dormant Az allele of MKK3 is recessive; the N260T substitution in this allele decreases MKK3 kinase activity and appears to be causal for Qsd2-AK. The N260T substitution occurred in the immediate ancestor allele of the dormant allele, and the established dormant allele became prevalent in barley cultivars grown in East Asia, where the rainy season and harvest season often overlap. Our findings show fine-tuning of seed dormancy during domestication and provide key information for improving pre-harvest sprouting tolerance in barley and wheat.

Are PECTIN ESTERASE INHIBITOR Genes Involved in Mediating Resistance to Rhynchosporium commune in Barley?

A family of putative PECTIN ESTERASE INHIBITOR (PEI) genes, which were detected in the genomic region co-segregating with the resistance gene Rrs2 against scald caused by Rhynchosporium commune in barley, were characterized and tested for their possible involvement in mediating resistance to the pathogen by complementation and overexpression analysis. The sequences of the respective genes were derived from two BAC contigs originating from the susceptible cultivar ‘Morex’. For the genes HvPEI2, HvPEI3, HvPEI4 and HvPEI6, specific haplotypes for 18 resistant and 23 susceptible cultivars were detected after PCR-amplification and haplotype-specific CAPS-markers were developed. None of the tested candidate genes HvPEI2, HvPEI3 and HvPEI4 alone conferred a high resistance level in transgenic over-expression plants, though an improvement of the resistance level was observed especially with OE-lines for gene HvPEI4. These results do not confirm but also do not exclude an involvement of the PEI gene family in the response to the pathogen. A candidate for the resistance gene Rrs2 could not be identified yet. It is possible that Rrs2 is a PEI gene or another type of gene which has not been detected in the susceptible cultivar ‘Morex’ or the full resistance reaction requires the presence of several PEI genes.

High-density genetic map construction and gene mapping of pericarp color in wax gourd using specific-locus amplified fragment (SLAF) sequencing

Background

High-density map is a valuable tool for genetic and genomic analysis. Although wax gourd is a widely distributed vegetable of Cucurbitaceae and has important medicinal and health value, no genetic map has been constructed because of the lack of efficient markers. Specific-locus amplified fragment sequencing (SLAF-seq) is a newly developed high-throughput strategy for large-scale single nucleotide polymorphism (SNP) discovery and genotyping.

Results

In our present study, we constructed a high-density genetic map by using SLAF-seq and identified a locus controlling pericarp color in wax gourd. An F2 population of 140 individuals and their two parents were subjected to SLAF-seq. A total of 143.38 M pair-end reads were generated. The average sequencing depth was 26.51 in the maternal line (B214), 27.01 in the parental line (B227), and 5.11 in each F2 individual. When filtering low-depth SLAF tags, a total of 142,653 high-quality SLAFs were detected, and 22,151 of them were polymorphic, with a polymorphism rate of 15.42 %. And finally, 4,607 of the polymorphic markers were selected for genetic map construction, and 12 linkage groups (LGs) were generated. The map spanned 2,172.86 cM with an average distance between adjacent markers for 0.49 cM. The inheritance of pericarp color was also studied, which showed that the pericarp color was controlled by one single gene. And based on the newly constructed high-density map, a single locus locating on chromosome 5 was identified for controlling the pericarp color of wax gourd.

Conclusions

This is the first report of high-density genetic map construction and gene mapping in wax gourd, which will be served as an invaluable tool for gene mapping, marker assisted breeding, map-based gene cloning, comparative mapping and draft genome assembling of wax gourd.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2220-y) contains supplementary material, which is available to authorized users.

The Relationships between Development and Low Temperature Tolerance in Barley Near Isogenic Lines Differing for Flowering Behavior

Flowering time, vernalization requirement, photoperiod sensitivity and low temperature tolerance are key traits in the Triticeae. We characterized a set of isogenic genetic stocks-representing single and pairwise substitutions of spring alleles at the VRN-H1, VRN-H2 and VRN-H3 loci in a winter barley background-at the structural, functional and phenotypic levels. High density mapping with reference to the barley genome sequence confirmed that in all cases target VRN alleles were present in the near isogenic lines (NILs) and allowed estimates of introgression size (at the genetic and physical levels) and gene content. Expression data corroborated the structural and phenotypic results. The latter confirmed that substitution of a spring allele at any of the VRN loci is sufficient to eliminate vernalization requirement. There was no significant change in low temperature tolerance with substitution of a spring allele at VRN-H2, but there were significant losses in cold tolerance with substitutions at VRN-H1 and VRN-H3. Reductions in cold tolerance are ascribed to an accelerated transition from the vegetative to reproductive state. The set of NILs will be a rich resource for understanding the genetics of vernalization, low temperature tolerance and other traits encoded/regulated by genes within the introgressed intervals.

Improvement of barley genome annotations by deciphering the Haruna Nijo genome

Full-length (FL) cDNA sequences provide the most reliable evidence for the presence of genes in genomes. In this report, detailed gene structures of barley, whole genome shotgun (WGS) and additional transcript data of the cultivar Haruna Nijo were quality controlled and compared with the published Morex genome information. Haruna Nijo scaffolds have longer total sequence length with much higher N50 and fewer sequences than those in Morex WGS contigs. The longer Haruna Nijo scaffolds provided efficient FLcDNA mapping, resulting in high coverage and detection of the transcription start sites. In combination with FLcDNAs and RNA-Seq data from four different tissue samples of Haruna Nijo, we identified 51,249 gene models on 30,606 loci. Overall sequence similarity between Haruna Nijo and Morex genome was 95.99%, while that of exon regions was higher (99.71%). These sequence and annotation data of Haruna Nijo are combined with Morex genome data and released from a genome browser. The genome sequence of Haruna Nijo may provide detailed gene structures in addition to the current Morex barley genome information.

Application of Population Sequencing (POPSEQ) for Ordering and Imputing Genotyping-by-Sequencing Markers in Hexaploid Wheat

The advancement of next-generation sequencing technologies in conjunction with new bioinformatics tools enabled fine-tuning of sequence-based, high-resolution mapping strategies for complex genomes. Although genotyping-by-sequencing (GBS) provides a large number of markers, its application for association mapping and genomics-assisted breeding is limited by a large proportion of missing data per marker. For species with a reference genomic sequence, markers can be ordered on the physical map. However, in the absence of reference marker order, the use and imputation of GBS markers is challenging. Here, we demonstrate how the population sequencing (POPSEQ) approach can be used to provide marker context for GBS in wheat. The utility of a POPSEQ-based genetic map as a reference map to create genetically ordered markers on a chromosome for hexaploid wheat was validated by constructing an independent de novo linkage map of GBS markers from a Synthetic W7984 × Opata M85 recombinant inbred line (SynOpRIL) population. The results indicated that there is strong agreement between the independent de novo linkage map and the POPSEQ mapping approach in mapping and ordering GBS markers for hexaploid wheat. After ordering, a large number of GBS markers were imputed, thus providing a high-quality reference map that can be used for QTL mapping for different traits. The POPSEQ-based reference map and whole-genome sequence assemblies are valuable resources that can be used to order GBS markers and enable the application of highly accurate imputation methods to leverage the application GBS markers in wheat.

Physical Mapping of Bread Wheat Chromosome 5A: An Integrated Approach

The huge size, redundancy, and highly repetitive nature of the bread wheat [Triticum aestivum (L.)] genome, makes it among the most difficult species to be sequenced. To overcome these limitations, a strategy based on the separation of individual chromosomes or chromosome arms and the subsequent production of physical maps was established within the frame of the International Wheat Genome Sequence Consortium (IWGSC). A total of 95,812 bacterial artificial chromosome (BAC) clones of short-arm chromosome 5A (5AS) and long-arm chromosome 5A (5AL) arm-specific BAC libraries were fingerprinted and assembled into contigs by complementary analytical approaches based on the FingerPrinted Contig (FPC) and Linear Topological Contig (LTC) tools. Combined anchoring approaches based on polymerase chain reaction (PCR) marker screening, microarray, and sequence homology searches applied to several genomic tools (i.e., genetic maps, deletion bin map, neighbor maps, BAC end sequences (BESs), genome zipper, and chromosome survey sequences) allowed the development of a high-quality physical map with an anchored physical coverage of 75% for 5AS and 53% for 5AL with high portions (64 and 48%, respectively) of contigs ordered along the chromosome. In the genome of grasses, Brachypodium [Brachypodium distachyon (L.) Beauv.], rice (Oryza sativa L.), and sorghum [Sorghum bicolor (L.) Moench] homologs of genes on wheat chromosome 5A were separated into syntenic blocks on different chromosomes as a result of translocations and inversions during evolution. The physical map presented represents an essential resource for fine genetic mapping and map-based cloning of agronomically relevant traits and a reference for the 5A sequencing projects.

Assessing the Barley Genome Zipper and Genomic Resources for Breeding Purposes

The aim of this study was to estimate the accuracy and convergence of newly developed barley (Hordeum vulgare L.) genomic resources, primarily genome zipper (GZ) and population sequencing (POPSEQ), at the genome-wide level and to assess their usefulness in applied barley breeding by analyzing seven known loci. Comparison of barley GZ and POPSEQ maps to a newly developed consensus genetic map constructed with data from 13 individual linkage maps yielded an accuracy of 97.8% (GZ) and 99.3% (POPSEQ), respectively, regarding the chromosome assignment. The percentage of agreement in marker position indicates that on average only 3.7% GZ and 0.7% POPSEQ positions are not in accordance with their centimorgan coordinates in the consensus map. The fine-scale comparison involved seven genetic regions on chromosomes 1H, 2H, 4H, 6H, and 7H, harboring major genes and quantitative trait loci (QTL) for disease resistance. In total, 179 GZ loci were analyzed and 64 polymorphic markers were developed. Entirely, 89.1% of these were allocated within the targeted intervals and 84.2% followed the predicted order. Forty-four markers showed a match to a POPSEQ-anchored contig, the percentage of collinearity being 93.2%, on average. Forty-four markers allowed the identification of twenty-five fingerprinted contigs (FPCs) and a more clear delimitation of the physical regions containing the traits of interest. Our results demonstrate that an increase in marker density of barley maps by using new genomic data significantly improves the accuracy of GZ. In addition, the combination of different barley genomic resources can be considered as a powerful tool to accelerate barley breeding.

Characterization and fine mapping of a novel barley Stage Green-Revertible Albino Gene (HvSGRA) by Bulked Segregant Analysis based on SSR assay and Specific Length Amplified Fragment Sequencing

Background

Leaf color variations are common in plants. Herein we describe a natural mutant of barley cultivar Edamai No.6, whs18, whose leaf color showed stable and inheritable stage-green-revertible-albino under field condition.

Methods

Bulked Segregant Analysis (BSA) based on SSR assay and Specific Length Amplified Fragment Sequencing (SLAF-seq) was used to map the candidate gene for this trait.

Results

We found that leaf color of whs18 was green at seedling stage, while the seventh or eighth leaf began to show etiolation, and albino leaves emerged after a short period. The newly emerged leaves began to show stripe white before jointing stage, and normal green leaves emerged gradually. The duration of whs18 with abnormal leaf color lasted for about 3 months, which had some negative impacts on yield-related-traits. Further investigations showed that the variation was associated with changes in chlorophyII content and chloroplast development. Genetic analysis revealed that the trait was controlled by a single recessive nuclear gene, and was designed as HvSGRA in this study. Based on the F₂ population derived from Edamai No.9706 and whs18, we initially mapped the HvSGRA gene on the short arm of chromosome 2H using SSR and BSA. GBMS247 on 2HS showed co-segregation with HvSGRA. The genetic distance between the other marker GBM1187 and HvSGRA was 1.2 cM. Further analysis using BSA with SLAF-seq also identified this region as candidate region. Finally, HvSGRA interval was narrowed to 0.4 cM between morex_contig_160447 and morex_contig_92239, which were anchored to two adjacent FP contigs, contig_34437 and contig_46434, respectively. Furthermore, six putative genes with high-confidence in this interval were identified by POPSEQ. Further analysis showed that the substitution from C to A in the third exon of fructokinase-1-like gene generated a premature stop codon in whs18, which may lead to loss function of this gene.

Conclusions

Using SSR and SLAF-seq in conjunction with BSA, we mapped HvSGRA within two adjacent FP contigs of barley. The mutation of fructokinase-1-like gene in whs18 may cause the stage green-revertible albino of barley. The current study lays foundation for hierarchical map-based cloning of HvSGRA and utilizing the gene/trait as a visualized maker in molecular breeding in future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2015-1) contains supplementary material, which is available to authorized users.

Development of genome-wide InDel markers and their integration with SSR, DArT and SNP markers in single barley map

Background

Development of molecular markers such as SSR (simple sequence repeat), DArT (diversity arrays technology) and SNP (single nucleotide polymorphism) is fundamental for linkage map construction and QTL mapping. However, DArT and SNP genotyping require special tools, and detection of SSR polymorphisms requires time-consuming polyacrylamide electrophoresis. Furthermore, many markers have been mapped in different populations such that their genetic positions are inconsistent. Recently, InDel (insertion and deletion) markers have become popular in genetic map construction and map-based cloning.

Results

Aligning genomic DNA sequences in two barley cultivars (Morex and Barke) identified 436,640 InDels. We designed 1140 InDel markers across the barley genome with an average genetic distance of 1 cM, each having a unique location in the barley genome. High-resolution melting (HRM) technology was used to genotype 55 InDel markers; those PCR amplicons with melting temperature differences >0.3 °C were ideal for HRM genotyping. The 1140 InDel markers together with 383 SSRs, 3909 gene-based SNPs and 1544 DArT markers were integrated into single barley genetic map according to their physical map positions.

Conclusions

High-density InDel markers with specific genome locations were developed, with 6976 molecular markers (SSRs, DArTs, SNPs and InDels) integrated into single barley genetic map. HRM genotyping of the InDel markers each with single PCR band will facilitate quick map construction and gene fine-mapping.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2027-x) contains supplementary material, which is available to authorized users.

Sequencing of 15 622 gene-bearing BACs clarifies the gene-dense regions of the barley genome

Barley (Hordeum vulgare L.) possesses a large and highly repetitive genome of 5.1 Gb that has hindered the development of a complete sequence. In 2012, the International Barley Sequencing Consortium released a resource integrating whole-genome shotgun sequences with a physical and genetic framework. However, because only 6278 bacterial artificial chromosome (BACs) in the physical map were sequenced, fine structure was limited. To gain access to the gene-containing portion of the barley genome at high resolution, we identified and sequenced 15 622 BACs representing the minimal tiling path of 72 052 physical-mapped gene-bearing BACs. This generated ~1.7 Gb of genomic sequence containing an estimated 2/3 of all Morex barley genes. Exploration of these sequenced BACs revealed that although distal ends of chromosomes contain most of the gene-enriched BACs and are characterized by high recombination rates, there are also gene-dense regions with suppressed recombination. We made use of published map-anchored sequence data from Aegilops tauschii to develop a synteny viewer between barley and the ancestor of the wheat D-genome. Except for some notable inversions, there is a high level of collinearity between the two species. The software HarvEST:Barley provides facile access to BAC sequences and their annotations, along with the barley-Ae. tauschii synteny viewer. These BAC sequences constitute a resource to improve the efficiency of marker development, map-based cloning, and comparative genomics in barley and related crops. Additional knowledge about regions of the barley genome that are gene-dense but low recombination is particularly relevant.

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.

Application of whole genome re-sequencing data in the development of diagnostic DNA markers tightly linked to a disease-resistance locus for marker-assisted selection in lupin (Lupinus angustifolius)

BACKGROUND

Molecular marker-assisted breeding provides an efficient tool to develop improved crop varieties. A major challenge for the broad application of markers in marker-assisted selection is that the marker phenotypes must match plant phenotypes in a wide range of breeding germplasm. In this study, we used the legume crop species Lupinus angustifolius (lupin) to demonstrate the utility of whole genome sequencing and re-sequencing on the development of diagnostic markers for molecular plant breeding.

RESULTS

Nine lupin cultivars released in Australia from 1973 to 2007 were subjected to whole genome re-sequencing. The re-sequencing data together with the reference genome sequence data were used in marker development, which revealed 180,596 to 795,735 SNP markers from pairwise comparisons among the cultivars. A total of 207,887 markers were anchored on the lupin genetic linkage map. Marker mining obtained an average of 387 SNP markers and 87 InDel markers for each of the 24 genome sequence assembly scaffolds bearing markers linked to 11 genes of agronomic interest. Using the R gene PhtjR conferring resistance to phomopsis stem blight disease as a test case, we discovered 17 candidate diagnostic markers by genotyping and selecting markers on a genetic linkage map. A further 243 candidate diagnostic markers were discovered by marker mining on a scaffold bearing non-diagnostic markers linked to the PhtjR gene. Nine out from the ten tested candidate diagnostic markers were confirmed as truly diagnostic on a broad range of commercial cultivars. Markers developed using these strategies meet the requirements for broad application in molecular plant breeding.

CONCLUSIONS

We demonstrated that low-cost genome sequencing and re-sequencing data were sufficient and very effective in the development of diagnostic markers for marker-assisted selection. The strategies used in this study may be applied to any trait or plant species. Whole genome sequencing and re-sequencing provides a powerful tool to overcome current limitations in molecular plant breeding, which will enable plant breeders to precisely pyramid favourable genes to develop super crop varieties to meet future food demands.