A high density physical map of chromosome 1BL supports evolutionary studies, map-based cloning and sequencing in wheat

Physical mapping of the wheat genes in low-recombination regions: radiation hybrid mapping of the C-locus

KEY MESSAGE

This work reports the physical mapping of an important gene affecting spike compactness located in a low-recombination region of hexaploid wheat. This work paves the way for the eventual isolation and characterization of the factor involved but also opens up possibilities to use this approach to precisely map other wheat genes located on proximal parts of wheat chromosomes that show highly reduced recombination. Mapping wheat genes, in the centromeric and pericentromeric regions (~ 2/3rd of a given chromosome), poses a formidable challenge due to highly suppressed recombination. Using an example of compact spike locus (C-locus), this study provides an approach to precisely map wheat genes in the pericentromeric and centromeric regions that house ~ 30% of wheat genes. In club-wheat, spike compactness is controlled by the dominant C-locus, but previous efforts have failed to localize it, on a particular arm of chromosome 2D. We integrated radiation hybrid (RH) and high-resolution genetic mapping to locate C-locus on the short arm of chromosome 2D. Flanking markers of the C-locus span a physical distance of 11.0 Mb (231.0-242 Mb interval) and contain only 11 high-confidence annotated genes. This work demonstrates the value of this integrated strategy in mapping dominant genes in the low-recombination regions of the wheat genome. A comparison of the mapping resolutions of the RH and genetic maps using common anchored markers indicated that the RH map provides ~ 9 times better resolution that the genetic map even with much smaller population size. This study provides a broadly applicable approach to fine map wheat genes in regions of suppressed recombination.

Breeding for Economically and Environmentally Sustainable Wheat Varieties: An Integrated Approach from Genomics to Selection

There is currently a strong societal demand for sustainability, quality, and safety in bread wheat production. To address these challenges, new and innovative knowledge, resources, tools, and methods to facilitate breeding are needed. This starts with the development of high throughput genomic tools including single nucleotide polymorphism (SNP) arrays, high density molecular marker maps, and full genome sequences. Such powerful tools are essential to perform genome-wide association studies (GWAS), to implement genomic and phenomic selection, and to characterize the worldwide diversity. This is also useful to breeders to broaden the genetic basis of elite varieties through the introduction of novel sources of genetic diversity. Improvement in varieties particularly relies on the detection of genomic regions involved in agronomical traits including tolerance to biotic (diseases and pests) and abiotic (drought, nutrient deficiency, high temperature) stresses. When enough resolution is achieved, this can result in the identification of candidate genes that could further be characterized to identify relevant alleles. Breeding must also now be approached through in silico modeling to simulate plant development, investigate genotype × environment interactions, and introduce marker-trait linkage information in the models to better implement genomic selection. Breeders must be aware of new developments and the information must be made available to the world wheat community to develop new high-yielding varieties that can meet the challenge of higher wheat production in a sustainable and fluctuating agricultural context. In this review, we compiled all knowledge and tools produced during the BREEDWHEAT project to show how they may contribute to face this challenge in the coming years.

The effects of crop attributes, selection, and recombination on Canadian bread wheat molecular variation

Cultivated germplasm provides an opportunity to investigate how crop agronomic traits, selection for major genes, and differences in crossing-over rates drive patterns of allelic variation. To identify how these factors correlated with allelic variation within a collection of cultivated bread wheat (Triticum aestivum L.), we generated genotypes for 388 accessions grown in Canada over the past 170 yr using filtered single nucleotide polymorphism (SNP) calls from an Illumina Wheat iSelect 90K SNP-array. Entries' breeding program, era of release, grain texture, kernel color, and growth habit contributed to allelic differentiation. Allelic diversity and linkage disequilibrium (LD) of markers flanking some major loci known to affect traits such as gluten strength, growth habit, and grain color were consistent with selective sweeps. Nonetheless, some flanking markers of major loci had low LD and high allelic diversity. Positive selection may have acted upon homoeologous genes that had significant enrichment for the gene ontology terms 'response-to-auxin' and 'response-to-wounding.' Long regions of LD, spanning approximately one-third the length of entire chromosomes, were associated with many pericentromeric regions. These regions were also characterized by low diversity. Enhancing recombination across these regions could generate novel allele combinations to accelerate Canadian wheat improvement.

Chromosome genomics uncovers plant genome organization and function

The identification of causal genomic loci and their interactions underlying various traits in plants has been greatly aided by progress in understanding the organization of the nuclear genome. This provides clues to the responses of plants to environmental stimuli at the molecular level. Apart from other uses, these insights are needed to fully explore the potential of new breeding techniques that rely on genome editing. However, genome analysis and sequencing is not straightforward in the many agricultural crops and their wild relatives that possess large and complex genomes. Chromosome genomics streamlines this task by dissecting the genome to single chromosomes whose DNA is then used instead of nuclear DNA. This results in a massive and lossless reduction in DNA sample complexity, reduces the time and cost of the experiment, and simplifies data interpretation. Flow cytometric sorting of condensed mitotic chromosomes makes it possible to purify single chromosomes in large quantities, and as the DNA remains intact this process can be coupled successfully with many techniques in molecular biology and genomics. Since the first experiments with flow cytometric sorting in the late 1980s, numerous applications have been developed, and chromosome genomics has been having a significant impact in many areas of research, including the sequencing of complex genomes of important crops and gene cloning. This review discusses these applications, describes their contribution to advancements in plant genome analysis and gene cloning, and outlines future directions.

Soybean BARCSoySNP6K: An assay for soybean genetics and breeding research

The limited number of recombinant events in recombinant inbred lines suggests that for a biparental population with a limited number of recombinant inbred lines, it is unnecessary to genotype the lines with many markers. For genomic prediction and selection, previous studies have demonstrated that only 1000-2000 genome-wide common markers across all lines/accessions are needed to reach maximum efficiency of genomic prediction in populations. Evaluation of too many markers will not only increase the cost but also generate redundant information. We developed a soybean (Glycine max) assay, BARCSoySNP6K, containing 6000 markers, which were carefully chosen from the SoySNP50K assay based on their position in the soybean genome and haplotype block, polymorphism among accessions and genotyping quality. The assay includes 5000 single nucleotide polymorphisms (SNPs) from euchromatic and 1000 from heterochromatic regions. The percentage of SNPs with minor allele frequency >0.10 was 95% and 91% in the euchromatic and heterochromatic regions, respectively. Analysis of progeny from two large families genotyped with SoySNP50K versus BARCSoySNP6K showed that the position of the common markers and number of unique bins along linkage maps were consistent based on the SNPs genotyped with the two assays; however, the rate of redundant markers was dramatically reduced with the BARCSoySNP6K. The BARCSoySNP6K assay is proven as an excellent tool for detecting quantitative trait loci, genomic selection and assessing genetic relationships. The assay is commercialized by Illumina Inc. and being used by soybean breeders and geneticists and the list of SNPs in the assay is an ideal resource for SNP genotyping by targeted amplicon sequencing.

The chromosome-scale genome reveals the evolution and diversification after the recent tetraploidization event in tea plant

Tea is one of the most popular nonalcoholic beverages due to its characteristic secondary metabolites with numerous health benefits. Although two draft genomes of tea plant (Camellia sinensis) have been published recently, the lack of chromosome-scale assembly hampers the understanding of the fundamental genomic architecture of tea plant and potential improvement. Here, we performed a genome-wide chromosome conformation capture technique (Hi-C) to obtain a chromosome-scale assembly based on the draft genome of C. sinensis var. sinensis and successfully ordered 2984.7 Mb (94.7%) scaffolds into 15 chromosomes. The scaffold N50 of the improved genome was 218.1 Mb, ~157-fold higher than that of the draft genome. Collinearity comparison of genome sequences and two genetic maps validated the high contiguity and accuracy of the chromosome-scale assembly. We clarified that only one Camellia recent tetraploidization event (CRT, 58.9-61.7 million years ago (Mya)) occurred after the core-eudicot common hexaploidization event (146.6-152.7 Mya). Meanwhile, 9243 genes (28.6%) occurred in tandem duplication, and most of these expanded after the CRT event. These gene duplicates increased functionally divergent genes that play important roles in tea-specific biosynthesis or stress response. Sixty-four catechin- and caffeine-related quantitative trait loci (QTLs) were anchored to chromosome assembly. Of these, two catechin-related QTL hotspots were derived from the CRT event, which illustrated that polyploidy has played a dramatic role in the diversification of tea germplasms. The availability of a chromosome-scale genome of tea plant holds great promise for the understanding of genome evolution and the discovery of novel genes contributing to agronomically beneficial traits in future breeding programs.

Genomic insights from the first chromosome-scale assemblies of oat (Avena spp.) diploid species

BACKGROUND

Cultivated hexaploid oat (Common oat; Avena sativa) has held a significant place within the global crop community for centuries; although its cultivation has decreased over the past century, its nutritional benefits have garnered increased interest for human consumption. We report the development of fully annotated, chromosome-scale assemblies for the extant progenitor species of the As- and Cp-subgenomes, Avena atlantica and Avena eriantha respectively. The diploid Avena species serve as important genetic resources for improving common oat's adaptive and food quality characteristics.

RESULTS

The A. atlantica and A. eriantha genome assemblies span 3.69 and 3.78 Gb with an N50 of 513 and 535 Mb, respectively. Annotation of the genomes, using sequenced transcriptomes, identified ~ 50,000 gene models in each species-including 2965 resistance gene analogs across both species. Analysis of these assemblies classified much of each genome as repetitive sequence (~ 83%), including species-specific, centromeric-specific, and telomeric-specific repeats. LTR retrotransposons make up most of the classified elements. Genome-wide syntenic comparisons with other members of the Pooideae revealed orthologous relationships, while comparisons with genetic maps from common oat clarified subgenome origins for each of the 21 hexaploid linkage groups. The utility of the diploid genomes was demonstrated by identifying putative candidate genes for flowering time (HD3A) and crown rust resistance (Pc91). We also investigate the phylogenetic relationships among other A- and C-genome Avena species.

CONCLUSIONS

The genomes we report here are the first chromosome-scale assemblies for the tribe Poeae, subtribe Aveninae. Our analyses provide important insight into the evolution and complexity of common hexaploid oat, including subgenome origin, homoeologous relationships, and major intra- and intergenomic rearrangements. They also provide the annotation framework needed to accelerate gene discovery and plant breeding.

Physical location of tandem repeats in the wheat genome and application for chromosome identification

A general distribution of tandem repeats (TRs) in the wheat genome was predicted and a new web page combined with fluorescence in situ hybridization experiments, and the newly developed Oligo probes will improve the resolution for wheat chromosome identification. Comprehensive sequence analysis of tandem repeats (TR) in the wheat reference genome permits discovery and application of TRs for chromosome identification. Genome-wide localization of TRs was identified in the reference sequences of Chinese Spring using Tandem Repeat Finder (TRF). A database of repeats unit size, array number, and physical coverage length of TRs in the wheat genome was built. The distribution of TRs occupied 3-5% of the wheat chromosomes, with non-random dispersal across the A, B, and D genomes. Three classes of TRs surrounding the predicted genes were compared. An optimized computer-assisted website page B2DSC was constructed for the general distribution and chromosomally enriched zones of TR sequences to be displayed graphically. The physical distribution of predicted TRs in the wheat genome by B2DSC matched well with the corresponding hybridization signals obtained with fluorescence in situ hybridization (FISH). We developed 20 oligonucleotide probes representing 20-60 bp lengths of high copy number of TRs and verified by FISH. An integrated physical map of TR-Oligo probes for wheat chromosome identification was constructed. Our results suggest that the combination of both molecular cytogenetics and genomic research will significantly benefit wheat breeding through chromosome manipulation and engineering.

Shifting the limits in wheat research and breeding using a fully annotated reference genome

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.

SNP-based pool genotyping and haplotype analysis accelerate fine-mapping of the wheat genomic region containing stripe rust resistance gene Yr26

NGS-assisted super pooling emerging as powerful tool to accelerate gene mapping and haplotype association analysis within target region uncovering specific linkage SNPs or alleles for marker-assisted gene pyramiding. Conventional gene mapping methods to identify genes associated with important agronomic traits require significant amounts of financial support and time. Here, a single nucleotide polymorphism (SNP)-based mapping approach, RNA-Seq and SNP array assisted super pooling analysis, was used for rapid mining of a candidate genomic region for stripe rust resistance gene Yr26 that has been widely used in wheat breeding programs in China. Large DNA and RNA super-pools were genotyped by Wheat SNP Array and sequenced by Illumina HiSeq, respectively. Hundreds of thousands of SNPs were identified and then filtered by multiple filtering criteria. Among selected SNPs, over 900 were found within an overlapping interval of less than 30 Mb as the Yr26 candidate genomic region in the centromeric region of chromosome arm 1BL. The 235 chromosome-specific SNPs were converted into KASP assays to validate the Yr26 interval in different genetic populations. Using a high-resolution mapping population (> 30,000 gametes), we confined Yr26 to a 0.003-cM interval. The Yr26 target region was anchored to the common wheat IWGSC RefSeq v1.0 and wild emmer WEWSeq v.1.0 sequences, from which 488 and 454 kb fragments were obtained. Several candidate genes were identified in the target genomic region, but there was no typical resistance gene in either genome region. Haplotype analysis identified specific SNPs linked to Yr26 and developed robust and breeder-friendly KASP markers. This integration strategy can be applied to accelerate generating many markers closely linked to target genes/QTL for a trait of interest in wheat and other polyploid species.

Chromosome-scale pseudomolecules refined by optical, physical and genetic maps in flax

Genomes of varying sizes have been sequenced with next-generation sequencing platforms. However, most reference sequences include draft unordered scaffolds containing chimeras caused by mis-scaffolding. A BioNano genome (BNG) optical map was constructed to improve the previously sequenced flax genome (Linum usitatissimum L., 2n = 30, about 373 Mb), which consisted of 3852 scaffolds larger than 1 kb and totalling 300.6 Mb. The high-resolution BNG map of cv. CDC Bethune totalled 317 Mb and consisted of 251 BNG contigs with an N50 of 2.15 Mb. A total of 622 scaffolds (286.6 Mb, 94.9%) aligned to 211 BNG contigs (298.6 Mb, 94.2%). Of those, 99 scaffolds, diagnosed to contain assembly errors, were refined into 225 new scaffolds. Using the newly refined scaffold sequences and the validated bacterial artificial chromosome-based physical map of CDC Bethune, the 211 BNG contigs were scaffolded into 94 super-BNG contigs (N50 of 6.64 Mb) that were further assigned to the 15 flax chromosomes using the genetic map. The pseudomolecules total about 316 Mb, with individual chromosomes of 15.6 to 29.4 Mb, and cover 97% of the annotated genes. Evidence from the chromosome-scale pseudomolecules suggests that flax has undergone palaeopolyploidization and mesopolyploidization events, followed by rearrangements and deletions or fusion of chromosome arms from an ancient progenitor with a haploid chromosome number of eight.

Features of the organization of bread wheat chromosome 5BS based on physical mapping

BACKGROUND

The IWGSC strategy for construction of the reference sequence of the bread wheat genome is based on first obtaining physical maps of the individual chromosomes. Our aim is to develop and use the physical map for analysis of the organization of the short arm of wheat chromosome 5B (5BS) which bears a number of agronomically important genes, including genes conferring resistance to fungal diseases.

RESULTS

A physical map of the 5BS arm (290 Mbp) was constructed using restriction fingerprinting and LTC software for contig assembly of 43,776 BAC clones. The resulting physical map covered ~ 99% of the 5BS chromosome arm (111 scaffolds, N50 = 3.078 Mb). SSR, ISBP and zipper markers were employed for anchoring the BAC clones, and from these 722 novel markers were developed based on previously obtained data from partial sequencing of 5BS. The markers were mapped using a set of Chinese Spring (CS) deletion lines, and F2 and RICL populations from a cross of CS and CS-5B dicoccoides. Three approaches have been used for anchoring BAC contigs on the 5BS chromosome, including clone-by-clone screening of BACs, GenomeZipper analysis, and comparison of BAC-fingerprints with in silico fingerprinting of 5B pseudomolecules of T. dicoccoides. These approaches allowed us to reach a high level of BAC contig anchoring: 96% of 5BS BAC contigs were located on 5BS. An interesting pattern was revealed in the distribution of contigs along the chromosome. Short contigs (200-999 kb) containing markers for the regions interrupted by tandem repeats, were mainly localized to the 5BS subtelomeric block; whereas the distribution of larger 1000-3500 kb contigs along the chromosome better correlated with the distribution of the regions syntenic to rice, Brachypodium, and sorghum, as detected by the Zipper approach.

CONCLUSION

The high fingerprinting quality, LTC software and large number of BAC clones selected by the informative markers in screening of the 43,776 clones allowed us to significantly increase the BAC scaffold length when compared with the published physical maps for other wheat chromosomes. The genetic and bioinformatics resources developed in this study provide new possibilities for exploring chromosome organization and for breeding applications.

Sequencing flow-sorted short arm of Haynaldia villosa chromosome 4V provides insights into its molecular structure and virtual gene order

BACKGROUND

Haynaldia villosa (H. villosa) has been recognized as a species potentially useful for wheat improvement. The availability of its genomic sequences will boost its research and application.

RESULTS

In this work, the short arm of H. villosa chromosome 4V (4VS) was sorted by flow cytometry and sequenced using Illumina platform. About 170.6 Mb assembled sequences were obtained. Further analysis showed that repetitive elements accounted for about 64.6% of 4VS, while the coding fraction, which is corresponding to 1977 annotated genes, represented 1.5% of the arm. The syntenic regions of the 4VS were searched and identified on wheat group 4 chromosomes 4AL, 4BS, 4DS, Brachypodium chromosomes 1 and 4, rice chromosomes 3 and 11, and sorghum chromosomes 1, 5 and 8. Based on genome-zipper analysis, a virtual gene order comprising 735 gene loci on 4VS genome was built by referring to the Brachypodium genome, which was relatively consistent with the scaffold order determined for Ae. tauschii chromosome 4D. The homologous alleles of several cloned genes on wheat group 4 chromosomes including Rht-1 gene were identified.

CONCLUSIONS

The sequences provided valuable information for mapping and positional-cloning genes located on 4VS, such as the wheat yellow mosaic virus resistance gene Wss1. The work on 4VS provided detailed insights into the genome of H. villosa, and may also serve as a model for sequencing the remaining parts of H. villosa genome.

Single-molecule sequencing and Hi-C-based proximity-guided assembly of amaranth (Amaranthus hypochondriacus) chromosomes provide insights into genome evolution

Background

Amaranth (Amaranthus hypochondriacus) was a food staple among the ancient civilizations of Central and South America that has recently received increased attention due to the high nutritional value of the seeds, with the potential to help alleviate malnutrition and food security concerns, particularly in arid and semiarid regions of the developing world. Here, we present a reference-quality assembly of the amaranth genome which will assist the agronomic development of the species.

Results

Utilizing single-molecule, real-time sequencing (Pacific Biosciences) and chromatin interaction mapping (Hi-C) to close assembly gaps and scaffold contigs, respectively, we improved our previously reported Illumina-based assembly to produce a chromosome-scale assembly with a scaffold N50 of 24.4 Mb. The 16 largest scaffolds contain 98% of the assembly and likely represent the haploid chromosomes (n = 16). To demonstrate the accuracy and utility of this approach, we produced physical and genetic maps and identified candidate genes for the betalain pigmentation pathway. The chromosome-scale assembly facilitated a genome-wide syntenic comparison of amaranth with other Amaranthaceae species, revealing chromosome loss and fusion events in amaranth that explain the reduction from the ancestral haploid chromosome number (n = 18) for a tetraploid member of the Amaranthaceae.

Conclusions

The assembly method reported here minimizes cost by relying primarily on short-read technology and is one of the first reported uses of in vivo Hi-C for assembly of a plant genome. Our analyses implicate chromosome loss and fusion as major evolutionary events in the 2n = 32 amaranths and clearly establish the homoeologous relationship among most of the subgenome chromosomes, which will facilitate future investigations of intragenomic changes that occurred post polyploidization.

Electronic supplementary material

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

Coordination of Bactericidal and Iron Regulatory Functions of Hepcidin in Innate Antimicrobial Immunity in a Zebrafish Model

Hepcidin acts as both an antimicrobial peptide and a hormonal regulator of iron homeostasis; however, the biological significance of this dual-function in immune reactions remains elusive. In this study, we provide experimental evidence regarding the coordination of this dual-function in the innate antimicrobial immunity using a zebrafish model. The transcription of hepcidin gene was significantly upregulated in liver by Aeromonas hydrophila (A.h) DNA stimulation, which was accompanied by an increase of hepcidin protein and a decrease of iron concentration in serum. Thus, an enhanced bactericidal activity against A.h and Escherichia coli and inhibitory effects on A.h growth and OmpA expression were observed in A.h cells, the latter of which made the bacterium more susceptible to complement attack. The enhanced bacteriostatic activities in serum following the stimulation were dramatically impaired by neutralizing hepcidin or restoring iron to the samples. Immuno-protection assay showed that zebrafish administrated with A.h DNA or designed CpG-ODNs had a significantly enhanced defence against A.h and Vibrio alginolyticus infections, which was also eliminated by the neutralization of hepcidin. Results indicate that the induction of hepcidin leads to the decrease of iron in circulation, which eventually limits iron availability to invading microorganisms, thus contributing to host defence.

Solar Radiation-Associated Adaptive SNP Genetic Differentiation in Wild Emmer Wheat, Triticum dicoccoides

Whole-genome scans with large number of genetic markers provide the opportunity to investigate local adaptation in natural populations and identify candidate genes under positive selection. In the present study, adaptation genetic differentiation associated with solar radiation was investigated using 695 polymorphic SNP markers in wild emmer wheat originated in a micro-site at Yehudiyya, Israel. The test involved two solar radiation niches: (1) sun, in-between trees; and (2) shade, under tree canopy, separated apart by a distance of 2-4 m. Analysis of molecular variance showed a small (0.53%) but significant portion of overall variation between the sun and shade micro-niches, indicating a non-ignorable genetic differentiation between sun and shade habitats. Fifty SNP markers showed a medium (0.05 ≤ FST ≤ 0.15) or high genetic differentiation (FST > 0.15). A total of 21 outlier loci under positive selection were identified by using four different FST -outlier testing algorithms. The markers and genome locations under positive selection are consistent with the known patterns of selection. These results suggested that genetic differentiation between sun and shade habitats is substantial, radiation-associated, and therefore ecologically determined. Hence, the results of this study reflected effects of natural selection through solar radiation on EST-related SNP genetic diversity, resulting presumably in different adaptive complexes at a micro-scale divergence. The present work highlights the evolutionary theory and application significance of solar radiation-driven natural selection in wheat improvement.

A High Resolution Radiation Hybrid Map of Wheat Chromosome 4A

Bread wheat has a large and complex allohexaploid genome with low recombination level at chromosome centromeric and peri-centromeric regions. This significantly hampers ordering of markers, contigs of physical maps and sequence scaffolds and impedes obtaining of high-quality reference genome sequence. Here we report on the construction of high-density and high-resolution radiation hybrid (RH) map of chromosome 4A supported by high-density chromosome deletion map. A total of 119 endosperm-based RH lines of two RH panels and 15 chromosome deletion bin lines were genotyped with 90K iSelect single nucleotide polymorphism (SNP) array. A total of 2316 and 2695 markers were successfully mapped to the 4A RH and deletion maps, respectively. The chromosome deletion map was ordered in 19 bins and allowed precise identification of centromeric region and verification of the RH panel reliability. The 4A-specific RH map comprises 1080 mapping bins and spans 6550.9 cR with a resolution of 0.13 Mb/cR. Significantly higher mapping resolution in the centromeric region was observed as compared to recombination maps. Relatively even distribution of deletion frequency along the chromosome in the RH panel was observed and putative functional centromere was delimited within a region characterized by two SNP markers.

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.

BioNano genome mapping of individual chromosomes supports physical mapping and sequence assembly in complex plant genomes

The assembly of a reference genome sequence of bread wheat is challenging due to its specific features such as the genome size of 17 Gbp, polyploid nature and prevalence of repetitive sequences. BAC-by-BAC sequencing based on chromosomal physical maps, adopted by the International Wheat Genome Sequencing Consortium as the key strategy, reduces problems caused by the genome complexity and polyploidy, but the repeat content still hampers the sequence assembly. Availability of a high-resolution genomic map to guide sequence scaffolding and validate physical map and sequence assemblies would be highly beneficial to obtaining an accurate and complete genome sequence. Here, we chose the short arm of chromosome 7D (7DS) as a model to demonstrate for the first time that it is possible to couple chromosome flow sorting with genome mapping in nanochannel arrays and create a de novo genome map of a wheat chromosome. We constructed a high-resolution chromosome map composed of 371 contigs with an N50 of 1.3 Mb. Long DNA molecules achieved by our approach facilitated chromosome-scale analysis of repetitive sequences and revealed a ~800-kb array of tandem repeats intractable to current DNA sequencing technologies. Anchoring 7DS sequence assemblies obtained by clone-by-clone sequencing to the 7DS genome map provided a valuable tool to improve the BAC-contig physical map and validate sequence assembly on a chromosome-arm scale. Our results indicate that creating genome maps for the whole wheat genome in a chromosome-by-chromosome manner is feasible and that they will be an affordable tool to support the production of improved pseudomolecules.

A whole-genome, radiation hybrid mapping resource of hexaploid wheat

Generating a contiguous, ordered reference sequence of a complex genome such as hexaploid wheat (2n = 6x = 42; approximately 17 GB) is a challenging task due to its large, highly repetitive, and allopolyploid genome. In wheat, ordering of whole-genome or hierarchical shotgun sequencing contigs is primarily based on recombination and comparative genomics-based approaches. However, comparative genomics approaches are limited to syntenic inference and recombination is suppressed within the pericentromeric regions of wheat chromosomes, thus, precise ordering of physical maps and sequenced contigs across the whole-genome using these approaches is nearly impossible. We developed a whole-genome radiation hybrid (WGRH) resource and tested it by genotyping a set of 115 randomly selected lines on a high-density single nucleotide polymorphism (SNP) array. At the whole-genome level, 26 299 SNP markers were mapped on the RH panel and provided an average mapping resolution of approximately 248 Kb/cR1500 with a total map length of 6866 cR1500 . The 7296 unique mapping bins provided a five- to eight-fold higher resolution than genetic maps used in similar studies. Most strikingly, the RH map had uniform bin resolution across the entire chromosome(s), including pericentromeric regions. Our research provides a valuable and low-cost resource for anchoring and ordering sequenced BAC and next generation sequencing (NGS) contigs. The WGRH developed for reference wheat line Chinese Spring (CS-WGRH), will be useful for anchoring and ordering sequenced BAC and NGS based contigs for assembling a high-quality, reference sequence of hexaploid wheat. Additionally, this study provides an excellent model for developing similar resources for other polyploid species.

Genomics as the key to unlocking the polyploid potential of wheat

Polyploidy has played a central role in plant genome evolution and in the formation of new species such as tetraploid pasta wheat and hexaploid bread wheat. Until recently, the high sequence conservation between homoeologous genes, together with the large genome size of polyploid wheat, had hindered genomic analyses in this important crop species. In the past 5 yr, however, the advent of next-generation sequencing has radically changed the wheat genomics landscape. Here, we review a series of advances in genomic resources and tools for functional genomics that are shifting the paradigm of what is possible in wheat molecular genetics and breeding. We discuss how understanding the relationship between homoeologues can inform approaches to modulate the response of quantitative traits in polyploid wheat; we also argue that functional redundancy has 'locked up' a wide range of phenotypic variation in wheat. We explore how genomics provides key tools to inform targeted manipulation of multiple homoeologues, thereby allowing researchers and plant breeders to unlock the full polyploid potential of wheat.

Exploring the tertiary gene pool of bread wheat: sequence assembly and analysis of chromosome 5M(g) of Aegilops geniculata

Next-generation sequencing (NGS) provides a powerful tool for the discovery of important genes and alleles in crop plants and their wild relatives. Despite great advances in NGS technologies, whole-genome shotgun sequencing is cost-prohibitive for species with complex genomes. An attractive option is to reduce genome complexity to a single chromosome prior to sequencing. This work describes a strategy for studying the genomes of distant wild relatives of wheat by isolating single chromosomes from addition or substitution lines, followed by chromosome sorting using flow cytometry and sequencing of chromosomal DNA by NGS technology. We flow-sorted chromosome 5M(g) from a wheat/Aegilops geniculata disomic substitution line [DS5M(g) (5D)] and sequenced it using an Illumina HiSeq 2000 system at approximately 50 × coverage. Paired-end sequences were assembled and used for structural and functional annotation. A total of 4236 genes were annotated on 5M(g) , in close agreement with the predicted number of genes on wheat chromosome 5D (4286). Single-gene FISH indicated no major chromosomal rearrangements between chromosomes 5M(g) and 5D. Comparing chromosome 5M(g) with model grass genomes identified synteny blocks in Brachypodium distachyon, rice (Oryza sativa), sorghum (Sorghum bicolor) and barley (Hordeum vulgare). Chromosome 5M(g) -specific SNPs and cytogenetic probe-based resources were developed and validated. Deletion bin-mapped and ordered 5M(g) SNP markers will be useful to track 5M-specific introgressions and translocations. This study provides a detailed sequence-based analysis of the composition of a chromosome from a distant wild relative of bread wheat, and opens up opportunities to develop genomic resources for wild germplasm to facilitate crop improvement.

Dissection of barley chromosomes 1H and 6H by the gametocidal system

We dissected barley chromosomes 1H and 6H added to common wheat by the gametocidal system and identified structural changes of the chromosomes by fluorescence in situ hybridization and genomic in situ hybridization. We found five aberrations of chromosome 1H, all of which lacked the long arm: one small fragment with the subtelomeric HvT01 sequence, one terminal deletion, and three telocentric chromosomes of the short arm. We established 33 dissection lines carrying single aberrant 6H chromosomes, of which 15 were deletions, 16 were translocations and two were isochromosomes. We conducted PCR analysis of the aberrant barley chromosomes using 75 and 81 EST markers specific to chromosomes 1H and 6H, respectively. This enabled us to construct a cytological map of chromosome 6H and to compare it to the previously reported genetic map and also to the physical map, which were released by the International Barley Genome Sequencing Consortium. The marker orders on the three maps were largely in agreement. The cytological map had better resolution in the proximal region of chromosome 6H than the corresponding genetic map. We discuss some of the discrepancies in marker order between the three maps that might be due to intraspecific polymorphism and gene duplication, as well as to technical problems inherent in the physical mapping process.

A high-resolution physical map integrating an anchored chromosome with the BAC physical maps of wheat chromosome 6B

Background

A complete genome sequence is an essential tool for the genetic improvement of wheat. Because the wheat genome is large, highly repetitive and complex due to its allohexaploid nature, the International Wheat Genome Sequencing Consortium (IWGSC) chose a strategy that involves constructing bacterial artificial chromosome (BAC)-based physical maps of individual chromosomes and performing BAC-by-BAC sequencing. Here, we report the construction of a physical map of chromosome 6B with the goal of revealing the structural features of the third largest chromosome in wheat.

Results

We assembled 689 informative BAC contigs (hereafter reffered to as contigs) representing 91 % of the entire physical length of wheat chromosome 6B. The contigs were integrated into a radiation hybrid (RH) map of chromosome 6B, with one linkage group consisting of 448 loci with 653 markers. The order and direction of 480 contigs, corresponding to 87 % of the total length of 6B, were determined. We also characterized the contigs that contained a part of the nucleolus organizer region or centromere based on their positions on the RH map and the assembled BAC clone sequences. Analysis of the virtual gene order along 6B using the information collected for the integrated map revealed the presence of several chromosomal rearrangements, indicating evolutionary events that occurred on chromosome 6B.

Conclusions

We constructed a reliable physical map of chromosome 6B, enabling us to analyze its genomic structure and evolutionary progression. More importantly, the physical map should provide a high-quality and map-based reference sequence that will serve as a resource for wheat chromosome 6B.

Electronic supplementary material

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

The physical map of wheat chromosome 5DS revealed gene duplications and small rearrangements

BACKGROUND

The substantially large bread wheat genome, organized into highly similar three sub-genomes, renders genomic research challenging. The construction of BAC-based physical maps of individual chromosomes reduces the complexity of this allohexaploid genome, enables elucidation of gene space and evolutionary relationships, provides tools for map-based cloning, and serves as a framework for reference sequencing efforts. In this study, we constructed the first comprehensive physical map of wheat chromosome arm 5DS, thereby exploring its gene space organization and evolution.

RESULTS

The physical map of 5DS was comprised of 164 contigs, of which 45 were organized into 21 supercontigs, covering 176 Mb with an N50 value of 2,173 kb. Fifty-eight of the contigs were larger than 1 Mb, with the largest contig spanning 6,649 kb. A total of 1,864 molecular markers were assigned to the map at a density of 10.5 markers/Mb, anchoring 100 of the 120 contigs (>5 clones) that constitute ~95 % of the cumulative length of the map. Ordering of 80 contigs along the deletion bins of chromosome arm 5DS revealed small-scale breaks in syntenic blocks. Analysis of the gene space of 5DS suggested an increasing gradient of genes organized in islands towards the telomere, with the highest gene density of 5.17 genes/Mb in the 0.67-0.78 deletion bin, 1.4 to 1.6 times that of all other bins.

CONCLUSIONS

Here, we provide a chromosome-specific view into the organization and evolution of the D genome of bread wheat, in comparison to one of its ancestors, revealing recent genome rearrangements. The high-quality physical map constructed in this study paves the way for the assembly of a reference sequence, from which breeding efforts will greatly benefit.

High-throughput physical map anchoring via BAC-pool sequencing

BACKGROUND

Physical maps created from large insert DNA libraries, typically cloned in BAC vector, are valuable resources for map-based cloning and de novo genome sequencing. The maps are most useful if contigs of overlapping DNA clones are anchored to chromosome(s), and ordered along them using molecular markers. Here we present a novel approach for anchoring physical maps, based on sequencing three-dimensional pools of BAC clones from minimum tilling path.

RESULTS

We used physical map of wheat chromosome arm 3DS to validate the method with two different DNA sequence datasets. The first comprised 567 genes ordered along the chromosome arm based on syntenic relationship of wheat with the sequenced genomes of Brachypodium, rice and sorghum. The second dataset consisted of 7,136 SNP-containing sequences, which were mapped genetically in Aegilops tauschii, the donor of the wheat D genome. Mapping of sequence reads from individual BAC pools to the first and the second datasets enabled unambiguous anchoring 447 and 311 3DS-specific sequences, respectively, or 758 in total.

CONCLUSIONS

We demonstrate the utility of the novel approach for BAC contig anchoring based on mass parallel sequencing of three-dimensional pools prepared from minimum tilling path of physical map. The existing genetic markers as well as any other DNA sequence could be mapped to BAC clones in a single in silico experiment. The approach reduces significantly the cost and time needed for anchoring and is applicable to any genomic project involving the construction of anchored physical map.

New insights into the wheat chromosome 4D structure and virtual gene order, revealed by survey pyrosequencing

Survey sequencing of the bread wheat (Triticum aestivum L.) genome (AABBDD) has been approached through different strategies delivering important information. However, the current wheat sequence knowledge is not complete. The aim of our study is to provide different and complementary set of data for chromosome 4D. A survey sequence was obtained by pyrosequencing of flow-sorted 4DS (7.2×) and 4DL (4.1×) arms. Single ends (SE) and long mate pairs (LMP) reads were assembled into contigs (223Mb) and scaffolds (65Mb) that were aligned to Aegilops tauschii draft genome (DD), anchoring 34Mb to chromosome 4. Scaffolds annotation rendered 822 gene models. A virtual gene order comprising 1973 wheat orthologous gene loci and 381 wheat gene models was built. This order was largely consistent with the scaffold order determined based on a published high density map from the Ae. tauschii chromosome 4, using bin-mapped 4D ESTs as a common reference. The virtual order showed a higher collinearity with homeologous 4B compared to 4A. Additionally, a virtual map was constructed and ∼5700 genes (∼2200 on 4DS and ∼3500 on 4DL) predicted. The sequence and virtual order obtained here using the 454 platform were compared with the Illumina one used by the IWGSC, giving complementary information.

Utilization of deletion bins to anchor and order sequences along the wheat 7B chromosome

A total of 3,671 sequence contigs and scaffolds were mapped to deletion bins on wheat chromosome 7B providing a foundation for developing high-resolution integrated physical map for this chromosome. Bread wheat (Triticum aestivum L.) has a large, complex and highly repetitive genome which is challenging to assemble into high quality pseudo-chromosomes. As part of the international effort to sequence the hexaploid bread wheat genome by the international wheat genome sequencing consortium (IWGSC) we are focused on assembling a reference sequence for chromosome 7B. The successful completion of the reference chromosome sequence is highly dependent on the integration of genetic and physical maps. To aid the integration of these two types of maps, we have constructed a high-density deletion bin map of chromosome 7B. Using the 270 K Nimblegen comparative genomic hybridization (CGH) array on a set of cv. Chinese spring deletion lines, a total of 3,671 sequence contigs and scaffolds (~7.8 % of chromosome 7B physical length) were mapped into nine deletion bins. Our method of genotyping deletions on chromosome 7B relied on a model-based clustering algorithm (Mclust) to accurately predict the presence or absence of a given genomic sequence in a deletion line. The bin mapping results were validated using three different approaches, viz. (a) PCR-based amplification of randomly selected bin mapped sequences (b) comparison with previously mapped ESTs and (c) comparison with a 7B genetic map developed in the present study. Validation of the bin mapping results suggested a high accuracy of the assignment of 7B sequence contigs and scaffolds to the 7B deletion bins.

The wheat-Septoria conflict: a new front opening up?

In the utopic absence of abiotic and/or biotic stressors, attaining the predicted increase (up to 70%) in wheat demand by 2050 in response to global population trends is a challenge. This objective becomes daunting, however, when one factors in the continuous constraint on global wheat production posed by Septoria tritici blotch (STB) disease. This is because, despite resistant loci being identified, a deficit of commercially relevant STB-resistant wheat germplasm remains. The issue is further compounded for growers by the emergence and prevalence of fungicide-resistant/insensitive strains of the causative pathogen Zymoseptoria tritici (formerly known as Mycosphaerella graminicola/Septoria tritici). However, biotechnology-based research is providing new opportunities in this struggle. As the exome response of wheat to STB attack begins to be deciphered, genes intrinsic to resistant and susceptible phenotypes will be identified. Combined with the application of genome-editing techniques and a growing appreciation of the complexity of wheat's and the dynamism of Z. tritici's genome, the generation of resulting STB-resistant wheat varieties will counter the prevalent threat of STB disease in wheat-production systems.

Identification and independent validation of a stable yield and thousand grain weight QTL on chromosome 6A of hexaploid wheat (Triticum aestivum L.)

BACKGROUND

Grain yield in wheat is a polygenic trait that is influenced by environmental and genetic interactions at all stages of the plant's growth. Yield is usually broken down into three components; number of spikes per area, grain number per spike, and grain weight (TGW). In polyploid wheat, studies have identified quantitative trait loci (QTL) which affect TGW, yet few have been validated and fine-mapped using independent germplasm, thereby having limited impact in breeding.

RESULTS

In this study we identified a major QTL for TGW, yield and green canopy duration on wheat chromosome 6A of the Spark x Rialto population, across 12 North European environments. Using independent germplasm in the form of BC2 and BC4 near isogenic lines (NILs), we validated the three QTL effects across environments. In four of the five experiments the Rialto 6A introgression gave significant improvements in yield (5.5%) and TGW (5.1%), with morphometric measurements showing that the increased grain weight was a result of wider grains. The extended green canopy duration associated with the high yielding/TGW Rialto allele was comprised of two independent effects; earlier flowering and delayed final maturity, and was expressed stably across the five environments. The wheat homologue (TaGW2) of a rice gene associated with increased TGW and grain width was mapped within the QTL interval. However, no polymorphisms were identified in the coding sequence between the parents.

CONCLUSION

The discovery and validation through near-isogenic lines of robust QTL which affect yield, green canopy duration, thousand grain weight, and grain width on chromosome 6A of hexaploid wheat provide an important first step to advance our understanding of the genetic mechanisms regulating the complex processes governing grain size and yield in polyploid wheat.

Whole-genome profiling and shotgun sequencing delivers an anchored, gene-decorated, physical map assembly of bread wheat chromosome 6A

Bread wheat (Triticum aestivum L.) is the most important staple food crop for 35% of the world's population. International efforts are underway to facilitate an increase in wheat production, of which the International Wheat Genome Sequencing Consortium (IWGSC) plays an important role. As part of this effort, we have developed a sequence-based physical map of wheat chromosome 6A using whole-genome profiling (WGP™). The bacterial artificial chromosome (BAC) contig assembly tools fingerprinted contig (fpc) and linear topological contig (ltc) were used and their contig assemblies were compared. A detailed investigation of the contigs structure revealed that ltc created a highly robust assembly compared with those formed by fpc. The ltc assemblies contained 1217 contigs for the short arm and 1113 contigs for the long arm, with an L50 of 1 Mb. To facilitate in silico anchoring, WGP™ tags underlying BAC contigs were extended by wheat and wheat progenitor genome sequence information. Sequence data were used for in silico anchoring against genetic markers with known sequences, of which almost 79% of the physical map could be anchored. Moreover, the assigned sequence information led to the 'decoration' of the respective physical map with 3359 anchored genes. Thus, this robust and genetically anchored physical map will serve as a framework for the sequencing of wheat chromosome 6A, and is of immediate use for map-based isolation of agronomically important genes/quantitative trait loci located on this chromosome.

The physical map of wheat chromosome 1BS provides insights into its gene space organization and evolution

Background

The wheat genome sequence is an essential tool for advanced genomic research and improvements. The generation of a high-quality wheat genome sequence is challenging due to its complex 17 Gb polyploid genome. To overcome these difficulties, sequencing through the construction of BAC-based physical maps of individual chromosomes is employed by the wheat genomics community. Here, we present the construction of the first comprehensive physical map of chromosome 1BS, and illustrate its unique gene space organization and evolution.

Results

Fingerprinted BAC clones were assembled into 57 long scaffolds, anchored and ordered with 2,438 markers, covering 83% of chromosome 1BS. The BAC-based chromosome 1BS physical map and gene order of the orthologous regions of model grass species were consistent, providing strong support for the reliability of the chromosome 1BS assembly. The gene space for chromosome 1BS spans the entire length of the chromosome arm, with 76% of the genes organized in small gene islands, accompanied by a two-fold increase in gene density from the centromere to the telomere.

Conclusions

This study provides new evidence on common and chromosome-specific features in the organization and evolution of the wheat genome, including a non-uniform distribution of gene density along the centromere-telomere axis, abundance of non-syntenic genes, the degree of colinearity with other grass genomes and a non-uniform size expansion along the centromere-telomere axis compared with other model cereal genomes. The high-quality physical map constructed in this study provides a solid basis for the assembly of a reference sequence of chromosome 1BS and for breeding applications.

Unlocking Triticeae genomics to sustainably feed the future

The tribe Triticeae includes the major crops wheat and barley. Within the last few years, the whole genomes of four Triticeae species-barley, wheat, Tausch's goatgrass (Aegilops tauschii) and wild einkorn wheat (Triticum urartu)-have been sequenced. The availability of these genomic resources for Triticeae plants and innovative analytical applications using next-generation sequencing technologies are helping to revitalize our approaches in genetic work and to accelerate improvement of the Triticeae crops. Comparative genomics and integration of genomic resources from Triticeae plants and the model grass Brachypodium distachyon are aiding the discovery of new genes and functional analyses of genes in Triticeae crops. Innovative approaches and tools such as analysis of next-generation populations, evolutionary genomics and systems approaches with mathematical modeling are new strategies that will help us discover alleles for adaptive traits to future agronomic environments. In this review, we provide an update on genomic tools for use with Triticeae plants and Brachypodium and describe emerging approaches toward crop improvements in Triticeae.