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

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.

Fine mapping and identification of a candidate gene for the barley Un8 true loose smut resistance gene

The candidate gene for the barley Un8 true loose smut resistance gene encodes a deduced protein containing two tandem protein kinase domains. In North America, durable resistance against all known isolates of barley true loose smut, caused by the basidiomycete pathogen Ustilago nuda (Jens.) Rostr. (U. nuda), is under the control of the Un8 resistance gene. Previous genetic studies mapped Un8 to the long arm of chromosome 5 (1HL). Here, a population of 4625 lines segregating for Un8 was used to delimit the Un8 gene to a 0.108 cM interval on chromosome arm 1HL, and assign it to fingerprinted contig 546 of the barley physical map. The minimal tilling path was identified for the Un8 locus using two flanking markers and consisted of two overlapping bacterial artificial chromosomes. One gene located close to a marker co-segregating with Un8 showed high sequence identity to a disease resistance gene containing two kinase domains. Sequence of the candidate gene from the parents of the segregating population, and in an additional 19 barley lines representing a broader spectrum of diversity, showed there was no intron in alleles present in either resistant or susceptible lines, and fifteen amino acid variations unique to the deduced protein sequence in resistant lines differentiated it from the deduced protein sequences in susceptible lines. Some of these variations were present within putative functional domains which may cause a loss of function in the deduced protein sequences within susceptible lines.

Barley: a translational model for adaptation to climate change

Barley (Hordeum vulgare ssp. vulgare) is an excellent model for understanding agricultural responses to climate change. Its initial domestication over 10 millennia ago and subsequent wide migration provide striking evidence of adaptation to different environments, agro-ecologies and uses. A bottleneck in the selection of modern varieties has resulted in a reduction in total genetic diversity and a loss of specific alleles relevant to climate-smart agriculture. However, extensive and well-curated collections of landraces, wild barley accessions (H. vulgare ssp. spontaneum) and other Hordeum species exist and are important new allele sources. A wide range of genomic and analytical tools have entered the public domain for exploring and capturing this variation, and specialized populations, mutant stocks and transgenics facilitate the connection between genetic diversity and heritable phenotypes. These lay the biological, technological and informational foundations for developing climate-resilient crops tailored to specific environments that are supported by extensive environmental and geographical databases, new methods for climate modelling and trait/environment association analyses, and decentralized participatory improvement methods. Case studies of important climate-related traits and their constituent genes - including examples that are indicative of the complexities involved in designing appropriate responses - are presented, and key developments for the future highlighted.

High-resolution genetic mapping and physical map construction for the fertility restorer Rfm1 locus in barley

High-resolution genetic linkage mapping and BAC physical mapping narrowed the fertility restorer locus Rfm1 in barley to a sub-centimorgan genetic interval and a 208-kb physical interval. Rfm1 restores the fertility of msm1 and msm2 male-sterile cytoplasms in barley. The fertility restoration gene is located on the short arm of chromosome 6H (6HS), and we pursued a positional cloning of this gene. Starting from a previous result that has delimited Rfm1 within a 10.8 cM region on 6HS, we developed novel CAPS and SSR markers tightly linked to the gene in barley using the sequence information from the syntenic region of rice and barley genome assemblies. Next, we performed fine mapping of the Rfm1 locus. To isolate recombinants, we surveyed 3,638 F2 plants derived from a cross between the CMS strain and the Rf strain with adjacent markers (NAS2090 and NAS1080). This analysis identified 175 recombinant plants from the F2 population to build a high-resolution map with nine markers tightly linked to the Rfm1 locus. Rfm1 was located within the 0.14 cM region delimited by two markers (NAS9113 and NAS9200). Using these flanking markers as well as marker cosegregating with Rfm1 (NAS9133), we screened the BAC libraries of the cultivar Morex, an rfm1 carrier. We isolated 11 BAC clones and constructed a BAC physical map using their fingerprints. Finally, we delimited the Rfm1 locus encompassing the rfm1 allele on a 208-kb contig composed of three minimally overlapping BAC clones. This precise localization of the Rfm1 locus in the barley genome is expected to greatly accelerate the future map-based cloning of the Rfm1 gene by sequence analysis and its genetic transformation for the complementation of cytoplasmic male-sterile plants.

Construction of a high-density integrated genetic linkage map of rubber tree (Hevea brasiliensis) using genotyping-by-sequencing (GBS)

Construction of linkage maps is crucial for genetic studies and marker-assisted breeding programs. Recent advances in next generation sequencing technologies allow for the generation of high-density linkage maps, especially in non-model species lacking extensive genomic resources. Here, we constructed a high-density integrated genetic linkage map of rubber tree (Hevea brasiliensis), the sole commercial producer of high-quality natural rubber. We applied a genotyping-by-sequencing (GBS) technique to simultaneously discover and genotype single nucleotide polymorphism (SNP) markers in two rubber tree populations. A total of 21,353 single nucleotide substitutions were identified, 55% of which represented transition events. GBS-based genetic maps of populations P and C comprised 1704 and 1719 markers and encompassed 2041 cM and 1874 cM, respectively. The average marker densities of these two maps were one SNP in 1.23-1.25 cM. A total of 1114 shared SNP markers were used to merge the two component maps. An integrated linkage map consisted of 2321 markers and spanned the cumulative length of 2052 cM. The composite map showed a substantial improvement in marker density, with one SNP marker in every 0.89 cM. To our knowledge, this is the most saturated genetic map in rubber tree to date. This integrated map allowed us to anchor 28,965 contigs, covering 135 Mb or 12% of the published rubber tree genome. We demonstrated that GBS is a robust and cost-effective approach for generating a common set of genome-wide SNP data suitable for constructing integrated linkage maps from multiple populations in a highly heterozygous agricultural species.

LAILAPS: the plant science search engine

With the number of sequenced plant genomes growing, the number of predicted genes and functional annotations is also increasing. The association between genes and phenotypic traits is currently of great interest. Unfortunately, the information available today is widely scattered over a number of different databases. Information retrieval (IR) has become an all-encompassing bioinformatics methodology for extracting knowledge from complex, heterogeneous and distributed databases, and therefore can be a useful tool for obtaining a comprehensive view of plant genomics, from genes to traits. Here we describe LAILAPS (http://lailaps.ipk-gatersleben.de), an IR system designed to link plant genomic data in the context of phenotypic attributes for a detailed forward genetic research. LAILAPS comprises around 65 million indexed documents, encompassing >13 major life science databases with around 80 million links to plant genomic resources. The LAILAPS search engine allows fuzzy querying for candidate genes linked to specific traits over a loosely integrated system of indexed and interlinked genome databases. Query assistance and an evidence-based annotation system enable time-efficient and comprehensive information retrieval. An artificial neural network incorporating user feedback and behavior tracking allows relevance sorting of results. We fully describe LAILAPS's functionality and capabilities by comparing this system's performance with other widely used systems and by reporting both a validation in maize and a knowledge discovery use-case focusing on candidate genes in barley.

Genome-Wide Linkage Mapping of QTL for Yield Components, Plant Height and Yield-Related Physiological Traits in the Chinese Wheat Cross Zhou 8425B/Chinese Spring

Identification of genes for yield components, plant height (PH), and yield-related physiological traits and tightly linked molecular markers is of great importance in marker-assisted selection (MAS) in wheat breeding. In the present study, 246 F8 RILs derived from the cross of Zhou 8425B/Chinese Spring were genotyped using the high-density Illumina iSelect 90K single nucleotide polymorphism (SNP) assay. Field trials were conducted at Zhengzhou and Zhoukou of Henan Province, during the 2012-2013 and 2013-2014 cropping season under irrigated conditions, providing data for four environments. Analysis of variance (ANOVA) of agronomic and physiological traits revealed significant differences (P < 0.01) among RILs, environments, and RILs × environments interactions. Broad-sense heritabilities of all traits including thousand kernel weight (TKW), PH, spike length (SL), kernel number per spike (KNS), spike number/m(2) (SN), normalized difference in vegetation index at anthesis (NDVI-A) and at 10 days post-anthesis (NDVI-10), SPAD value of chlorophyll content at anthesis (Chl-A) and at 10 days post-anthesis (Chl-10) ranged between 0.65 and 0.94. A linkage map spanning 3609.4 cM was constructed using 5636 polymorphic SNP markers, with an average chromosome length of 171.9 cM and marker density of 0.64 cM/marker. A total of 866 SNP markers were newly mapped to the hexaploid wheat linkage map. Eighty-six QTL for yield components, PH, and yield-related physiological traits were detected on 18 chromosomes except 1D, 5D, and 6D, explaining 2.3-33.2% of the phenotypic variance. Ten stable QTL were identified across four environments, viz. QTKW.caas-6A.1, QTKW.caas-7AL, QKNS.caas-4AL, QSN.caas-1AL.1, QPH.caas-4BS.2, QPH.caas-4DS.1, QSL.caas-4AS, QSL.caas-4AL.1, QChl-A.caas-5AL, and QChl-10.caas-5BL. Meanwhile, 10 QTL-rich regions were found on chromosome 1BS, 2AL (2), 3AL, 4AL (2), 4BS, 4DS, 5BL, and 7AL exhibiting pleiotropic effects. These QTL or QTL clusters are tightly linked to SNP markers, with genetic distances to the closest SNPs ranging from 0 to 1.5 cM, and could serve as target regions for fine mapping, candidate gene discovery, and MAS in wheat breeding.

An SNP-based saturated genetic map and QTL analysis of fruit-related traits in cucumber using specific-length amplified fragment (SLAF) sequencing

Background

Cucumber, Cucumis sativus L., is an economically important vegetable crop which is processed or consumed fresh worldwide. However, the narrow genetic base in cucumber makes it difficult for constructing high-density genetic maps. The development of massively parallel genotyping methods and next-generation sequencing (NGS) technologies provides an excellent opportunity for developing single nucleotide polymorphisms (SNPs) for linkage map construction and QTL analysis of horticultural traits. Specific-length amplified fragment sequencing (SLAF-seq) is a recent marker development technology that allows large-scale SNP discovery and genotyping at a reasonable cost. In this study, we constructed a high-density SNP map for cucumber using SLAF-seq and detected fruit-related QTLs.

Results

An F₂ population of 148 individuals was developed from an intra-varietal cross between CC3 and NC76. Genomic DNAs extracted from two parents and 148 F₂ individuals were subjected to high-throughput sequencing and SLAF library construction. A total of 10.76 Gb raw data and 75,024,043 pair-end reads were generated to develop 52,684 high-quality SLAFs, out of which 5,044 were polymorphic. 4,817 SLAFs were encoded and grouped into different segregation patterns. A high-resolution genetic map containing 1,800 SNPs was constructed for cucumber spanning 890.79 cM. The average distance between adjacent markers was 0.50 cM. 183 scaffolds were anchored to the SNP-based genetic map covering 46% (168.9 Mb) of the cucumber genome (367 Mb). Nine QTLs for fruit length and weight were detected, a QTL designated fl3.2 explained 44.60% of the phenotypic variance. Alignment of the SNP markers to draft genome scaffolds revealed two mis-assembled scaffolds that were validated by fluorescence in situ hybridization (FISH).

Conclusions

We report herein the development of evenly dispersed SNPs across cucumber genome, and for the first time an SNP-based saturated linkage map. This 1,800-locus map would likely facilitate genetic mapping of complex QTL loci controlling fruit yield, and the orientation of draft genome scaffolds.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1158) contains supplementary material, which is available to authorized users.

chromoWIZ: a web tool to query and visualize chromosome-anchored genes from cereal and model genomes

BACKGROUND

Over the last years reference genome sequences of several economically and scientifically important cereals and model plants became available. Despite the agricultural significance of these crops only a small number of tools exist that allow users to inspect and visualize the genomic position of genes of interest in an interactive manner.

DESCRIPTION

We present chromoWIZ, a web tool that allows visualizing the genomic positions of relevant genes and comparing these data between different plant genomes. Genes can be queried using gene identifiers, functional annotations, or sequence homology in four grass species (Triticum aestivum, Hordeum vulgare, Brachypodium distachyon, Oryza sativa). The distribution of the anchored genes is visualized along the chromosomes by using heat maps. Custom gene expression measurements, differential expression information, and gene-to-group mappings can be uploaded and can be used for further filtering.

CONCLUSIONS

This tool is mainly designed for breeders and plant researchers, who are interested in the location and the distribution of candidate genes as well as in the syntenic relationships between different grass species. chromoWIZ is freely available and online accessible at http://mips.helmholtz-muenchen.de/plant/chromoWIZ/index.jsp.

Unlocking the secondary gene-pool of barley with next-generation sequencing

Crop wild relatives (CWR) provide an important source of allelic diversity for any given crop plant species for counteracting the erosion of genetic diversity caused by domestication and elite breeding bottlenecks. Hordeum bulbosum L. is representing the secondary gene pool of the genus Hordeum. It has been used as a source of genetic introgressions for improving elite barley germplasm (Hordeum vulgare L.). However, genetic introgressions from H. bulbosum have yet not been broadly applied, due to a lack of suitable molecular tools for locating, characterizing, and decreasing by recombination and marker-assisted backcrossing the size of introgressed segments. We applied next-generation sequencing (NGS) based strategies for unlocking genetic diversity of three diploid introgression lines of cultivated barley containing chromosomal segments of its close relative H. bulbosum. Firstly, exome capture-based (re)-sequencing revealed large numbers of single nucleotide polymorphisms (SNPs) enabling the precise allocation of H. bulbosum introgressions. This SNP resource was further exploited by designing a custom multiplex SNP genotyping assay. Secondly, two-enzyme-based genotyping-by-sequencing (GBS) was employed to allocate the introgressed H. bulbosum segments and to genotype a mapping population. Both methods provided fast and reliable detection and mapping of the introgressed segments and enabled the identification of recombinant plants. Thus, the utilization of H. bulbosum as a resource of natural genetic diversity in barley crop improvement will be greatly facilitated by these tools in the future.

Genetic anchoring of whole-genome shotgun assemblies

The recent advances in sequencing throughput and genome assembly algorithms have established whole-genome shotgun (WGS) assemblies as the cornerstone of the genomic infrastructure for many species. WGS assemblies can be constructed with comparative ease and give a comprehensive representation of the gene space even of large and complex genomes. One major obstacle in utilizing WGS assemblies for important research applications such as gene isolation or comparative genomics has been the lack of chromosomal positioning and contextualization of short sequence contigs. Assigning chromosomal locations to sequence contigs required the construction and integration of genome-wide physical maps and dense genetic linkage maps as well as synteny to model species. Recently, methods to rapidly construct ultra-dense linkage maps encompassing millions of genetic markers from WGS sequencing data of segregating populations have made possible the direct assignment of genetic positions to short sequence contigs. Here, we review recent developments in the integration of WGS assemblies and sequence-based linkage maps, discuss challenges for further improvement of the methodology and outline possible applications building on genetically anchored WGS assemblies.

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.

Reduced chlorophyll biosynthesis in heterozygous barley magnesium chelatase mutants

Chlorophyll biosynthesis is initiated by magnesium chelatase, an enzyme composed of three proteins, which catalyzes the insertion of Mg2+ into protoporphyrin IX to produce Mg-protoporphyrin IX. In barley (Hordeum vulgare L.) the three proteins are encoded by Xantha-f, Xantha-g and Xantha-h. Two of the gene products, XanH and XanG, belong to the structurally conserved family of AAA+ proteins (ATPases associated with various cellular activities) and form a complex involving six subunits of each protein. The complex functions as an ATP-fueled motor of the magnesium chelatase that uses XanF as substrate, which is the catalytic subunit responsible for the insertion of Mg2+ into protoporphyrin IX. Previous studies have shown that semi-dominant Xantha-h mutations result in non-functional XanH subunits that participate in the formation of inactive AAA complexes. In the present study, we identify severe mutations in the barley mutants xantha-h.38, -h.56 and -h.57. A truncated form of the protein is seen in xantha-h.38, whereas no XanH is detected in xantha-h.56 and -h.57. Heterozygous mutants show a reduction in chlorophyll content by 14-18% suggesting a slight semi-dominance of xantha-h.38, -h.56 and -h.57, which otherwise have been regarded as recessive mutations.