An integrated, systems-wide approach is needed for public-private partnerships to drive genetic innovation in crops

Public-private partnerships are key to successfully translate knowledge to products, but current frameworks do not foster the systems-wide approach required to improve crops to meet the agricultural production challenges of the 21st century.

Linking the International Wheat Genome Sequencing Consortium bread wheat reference genome sequence to wheat genetic and phenomic data

The Wheat@URGI portal has been developed to provide the international community of researchers and breeders with access to the bread wheat reference genome sequence produced by the International Wheat Genome Sequencing Consortium. Genome browsers, BLAST, and InterMine tools have been established for in-depth exploration of the genome sequence together with additional linked datasets including physical maps, sequence variations, gene expression, and genetic and phenomic data from other international collaborative projects already stored in the GnpIS information system. The portal provides enhanced search and browser features that will facilitate the deployment of the latest genomics resources in wheat improvement.

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.

Molecular and FISH analyses of a 53-kbp intact DNA fragment inserted by biolistics in wheat (Triticum aestivum L.) genome

A large, 53-kbp, intact DNA fragment was inserted into the wheat ( Triticum aestivum L.) genome. FISH analyses of individual transgenic events revealed multiple insertions of intact fragments. Transferring large intact DNA fragments containing clusters of resistance genes or complete metabolic pathways into the wheat genome remains a challenge. In a previous work, we showed that the use of dephosphorylated cassettes for wheat transformation enabled the production of simple integration patterns. Here, we used the same technology to produce a cassette containing a 44-kb Arabidopsis thaliana BAC, flanked by one selection gene and one reporter gene. This 53-kb linear cassette was integrated in the bread wheat (Triticum aestivum L.) genome by biolistic transformation. Our results showed that transgenic plants harboring the entire cassette were generated. The inheritability of the cassette was demonstrated in the T1 and T2 generation. Surprisingly, FISH analysis performed on T1 progeny of independent events identified double genomic insertions of intact fragments in non-homoeologous positions. Inheritability of these double insertions was demonstrated by FISH analysis of the T1 generation. Relative conclusions that can be drawn from molecular or FISH analysis are discussed along with future prospects of the engineering of large fragments for wheat transformation or genome editing.

Improving global integration of crop research

Field laboratories in realistic crop environments are needed

Exploiting the Repetitive Fraction of the Wheat Genome for High-Throughput Single-Nucleotide Polymorphism Discovery and Genotyping

Transposable elements (TEs) account for more than 80% of the wheat genome. Although they represent a major obstacle for genomic studies, TEs are also a source of polymorphism and consequently of molecular markers such as insertion site-based polymorphism (ISBP) markers. Insertion site-based polymorphisms have been found to be a great source of genome-specific single-nucleotide polymorphism (SNPs) in the hexaploid wheat ( L.) genome. Here, we report on the development of a high-throughput SNP discovery approach based on sequence capture of ISBP markers. By applying this approach to the reference sequence of chromosome 3B from hexaploid wheat, we designed 39,077 SNPs that are evenly distributed along the chromosome. We demonstrate that these SNPs can be efficiently scored with the KASPar (Kompetitive allele-specific polymerase chain reaction) genotyping technology. Finally, through genetic diversity and genome-wide association studies, we also demonstrate that ISBP-derived SNPs can be used in marker-assisted breeding programs.

Fine mapping of a large-effect QTL conferring Fusarium crown rot resistance on the long arm of chromosome 3B in hexaploid wheat

BACKGROUND

Fusarium crown rot (FCR) is a major cereal disease in semi-arid areas worldwide. Of the various QTL reported, the one on chromosome arm 3BL (Qcrs.cpi-3B) has the largest effect that can be consistently detected in different genetic backgrounds. Nine sets of near isogenic lines (NILs) for this locus were made available in a previous study. To identify markers that could be reliably used in tagging the Qcrs.cpi-3B locus, a NIL-derived population consisting of 774 F10 lines were generated and exploited to assess markers selected from the existing linkage map and generated from sequences of the 3B pseudomolecule.

RESULTS

This is the first report on fine mapping a QTL conferring FCR resistance in wheat. By three rounds of linkage mapping using the NILs and the NIL-derived population, the Qcrs.cpi-3B locus was mapped to an interval of 0.7 cM covering a physical distance of about 1.5 Mb. Seven markers co-segregating with the locus were developed. This interval contains a total of 63 gene-coding sequences based on the 3B pseudomolecule, and six of them were known to encode disease resistance proteins. Several of the genes in this interval were among those responsive to FCR infection detected in an earlier study.

CONCLUSIONS

The accurate localization of the Qcrs.cpi-3B locus and the development of the markers co-segregating with it should facilitate the incorporation of this large-effect QTL conferring FCR resistance into breeding programs as well as the cloning of the gene(s) underlying the QTL.

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.

Deep transcriptome sequencing provides new insights into the structural and functional organization of the wheat genome

Background

Because of its size, allohexaploid nature, and high repeat content, the bread wheat genome is a good model to study the impact of the genome structure on gene organization, function, and regulation. However, because of the lack of a reference genome sequence, such studies have long been hampered and our knowledge of the wheat gene space is still limited. The access to the reference sequence of the wheat chromosome 3B provided us with an opportunity to study the wheat transcriptome and its relationships to genome and gene structure at a level that has never been reached before.

Results

By combining this sequence with RNA-seq data, we construct a fine transcriptome map of the chromosome 3B. More than 8,800 transcription sites are identified, that are distributed throughout the entire chromosome. Expression level, expression breadth, alternative splicing as well as several structural features of genes, including transcript length, number of exons, and cumulative intron length are investigated. Our analysis reveals a non-monotonic relationship between gene expression and structure and leads to the hypothesis that gene structure is determined by its function, whereas gene expression is subject to energetic cost. Moreover, we observe a recombination-based partitioning at the gene structure and function level.

Conclusions

Our analysis provides new insights into the relationships between gene and genome structure and function. It reveals mechanisms conserved with other plant species as well as superimposed evolutionary forces that shaped the wheat gene space, likely participating in wheat adaptation.

Electronic supplementary material

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

QTug.sau-3B is a major quantitative trait locus for wheat hexaploidization

Meiotic nonreduction resulting in unreduced gametes is thought to be the predominant mechanism underlying allopolyploid formation in plants. Until now, however, its genetic base was largely unknown. The allohexaploid crop common wheat (Triticum aestivum L.), which originated from hybrids of T. turgidum L. with Aegilops tauschii Cosson, provides a model to address this issue. Our observations of meiosis in pollen mother cells from T. turgidum×Ae. tauschii hybrids indicated that first division restitution, which exhibited prolonged cell division during meiosis I, was responsible for unreduced gamete formation. A major quantitative trait locus (QTL) for this trait, named QTug.sau-3B, was detected on chromosome 3B in two T. turgidum×Ae. tauschii haploid populations. This QTL is situated between markers Xgwm285 and Xcfp1012 and covered a genetic distance of 1 cM in one population. QTug.sau-3B is a haploid-dependent QTL because it was not detected in doubled haploid populations. Comparative genome analysis indicated that this QTL was close to Ttam-3B, a collinear homolog of tam in wheat. Although the relationship between QTug.sau-3B and Ttam requires further study, high frequencies of unreduced gametes may be related to reduced expression of Ttam in wheat.

Structural and functional partitioning of bread wheat chromosome 3B

We produced a reference sequence of the 1-gigabase chromosome 3B of hexaploid bread wheat. By sequencing 8452 bacterial artificial chromosomes in pools, we assembled a sequence of 774 megabases carrying 5326 protein-coding genes, 1938 pseudogenes, and 85% of transposable elements. The distribution of structural and functional features along the chromosome revealed partitioning correlated with meiotic recombination. Comparative analyses indicated high wheat-specific inter- and intrachromosomal gene duplication activities that are potential sources of variability for adaption. In addition to providing a better understanding of the organization, function, and evolution of a large and polyploid genome, the availability of a high-quality sequence anchored to genetic maps will accelerate the identification of genes underlying important agronomic traits.

Characterization of polyploid wheat genomic diversity using a high-density 90,000 single nucleotide polymorphism array

High-density single nucleotide polymorphism (SNP) genotyping arrays are a powerful tool for studying genomic patterns of diversity, inferring ancestral relationships between individuals in populations and studying marker-trait associations in mapping experiments. We developed a genotyping array including about 90,000 gene-associated SNPs and used it to characterize genetic variation in allohexaploid and allotetraploid wheat populations. The array includes a significant fraction of common genome-wide distributed SNPs that are represented in populations of diverse geographical origin. We used density-based spatial clustering algorithms to enable high-throughput genotype calling in complex data sets obtained for polyploid wheat. We show that these model-free clustering algorithms provide accurate genotype calling in the presence of multiple clusters including clusters with low signal intensity resulting from significant sequence divergence at the target SNP site or gene deletions. Assays that detect low-intensity clusters can provide insight into the distribution of presence-absence variation (PAV) in wheat populations. A total of 46 977 SNPs from the wheat 90K array were genetically mapped using a combination of eight mapping populations. The developed array and cluster identification algorithms provide an opportunity to infer detailed haplotype structure in polyploid wheat and will serve as an invaluable resource for diversity studies and investigating the genetic basis of trait variation in wheat.

Fine mapping of LrSV2, a race-specific adult plant leaf rust resistance gene on wheat chromosome 3BS

Fine mapping permits the precise positioning of genes within chromosomes, prerequisite for positional cloning that will allow its rational use and the study of the underlying molecular action mechanism. Three leaf rust resistance genes were identified in the durable leaf rust resistant Argentinean wheat variety Sinvalocho MA: the seedling resistance gene Lr3 on distal 6BL and two adult plant resistance genes, LrSV1 and LrSV2, on chromosomes 2DS and 3BS, respectively. To develop a high-resolution genetic map for LrSV2, 10 markers were genotyped on 343 F2 individuals from a cross between Sinvalocho MA and Gama6. The closest co-dominant markers on both sides of the gene (3 microsatellites and 2 STMs) were analyzed on 965 additional F2s from the same cross. Microsatellite marker cfb5010 cosegregated with LrSV2 whereas flanking markers were found at 1 cM distal and at 0.3 cM proximal to the gene. SSR markers designed from the sequences of cv Chinese Spring BAC clones spanning the LrSV2 genetic interval were tested on the recombinants, allowing the identification of microsatellite swm13 at 0.15 cM distal to LrSV2. This delimited an interval of 0.45 cM around the gene flanked by the SSR markers swm13 and gwm533 at the subtelomeric end of chromosome 3BS.

High-resolution analysis of a QTL for resistance to Stagonospora nodorum glume blotch in wheat reveals presence of two distinct resistance loci in the target interval

Stagonospora nodorum glume blotch (SNG), caused by the necrotrophic fungus Stagonospora nodorum, is one of the economically important diseases of bread wheat (Triticum aestivum L.). Resistance to SNG is known to be quantitative and previous studies of a recombinant inbred line (RIL) population identified a major quantitative trait locus (QTL) for resistance to SNG on the short arm of chromosome 3B. To localize this QTL (QSng.sfr-3BS) with high resolution, we constructed a genetic map for the QTL target region using information from sequenced flow-sorted chromosomes 3B of the two parental cultivars 'Arina' and 'Forno', the physical map of chromosome 3B of cultivar 'Chinese Spring' and BAC-clone sequences. The mapping population of near-isogenic lines (NIL) was evaluated for SNG resistance in field infection tests. NILs segregated for disease resistance as well as for plant height; additionally, we observed a high environmental influence on the trait. Our analysis detected a strong negative correlation of SNG resistance and plant height. Further analysis of the target region identified two linked loci associated with SNG resistance. One of them was also associated with plant height, revealing an effect of QSng.sfr-3BS on plant height that was hidden in the RIL population. This result demonstrates an unexpectedly high genetic complexity of resistance controlled by QSng.sfr-3BS and shows the importance of the study of QTL in mendelized form in NILs.

Organization and evolution of transposable elements along the bread wheat chromosome 3B

BACKGROUND

The 17 Gb bread wheat genome has massively expanded through the proliferation of transposable elements (TEs) and two recent rounds of polyploidization. The assembly of a 774 Mb reference sequence of wheat chromosome 3B provided us with the opportunity to explore the impact of TEs on the complex wheat genome structure and evolution at a resolution and scale not reached so far.

RESULTS

We develop an automated workflow, CLARI-TE, for TE modeling in complex genomes. We delineate precisely 56,488 intact and 196,391 fragmented TEs along the 3B pseudomolecule, accounting for 85% of the sequence, and reconstruct 30,199 nested insertions. TEs have been mostly silent for the last one million years, and the 3B chromosome has been shaped by a succession of bursts that occurred between 1 to 3 million years ago. Accelerated TE elimination in the high-recombination distal regions is a driving force towards chromosome partitioning. CACTAs overrepresented in the high-recombination distal regions are significantly associated with recently duplicated genes. In addition, we identify 140 CACTA-mediated gene capture events with 17 genes potentially created by exon shuffling and show that 19 captured genes are transcribed and under selection pressure, suggesting the important role of CACTAs in the recent wheat adaptation.

CONCLUSION

Accurate TE modeling uncovers the dynamics of TEs in a highly complex and polyploid genome. It provides novel insights into chromosome partitioning and highlights the role of CACTA transposons in the high level of gene duplication in wheat.

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.

Wheat syntenome unveils new evidences of contrasted evolutionary plasticity between paleo- and neoduplicated subgenomes

Bread wheat derives from a grass ancestor structured in seven protochromosomes followed by a paleotetraploidization to reach a 12 chromosomes intermediate and a neohexaploidization (involving subgenomes A, B and D) event that finally shaped the 21 modern chromosomes. Insights into wheat syntenome in sequencing conserved orthologous set (COS) genes unravelled differences in genomic structure (such as gene conservation and diversity) and genetical landscape (such as recombination pattern) between ancestral as well as recent duplicated blocks. Contrasted evolutionary plasticity is observed where the B subgenome appears more sensitive (i.e. plastic) in contrast to A as dominant (i.e. stable) in response to the neotetraploidization and D subgenome as supra-dominant (i.e. pivotal) in response to the neohexaploidization event. Finally, the wheat syntenome, delivered through a public web interface PlantSyntenyViewer at http://urgi.versailles.inra.fr/synteny-wheat, can be considered as a guide for accelerated dissection of major agronomical traits in wheat.

[A comparative genetic and cytogenetic mapping of wheat chromosome 5B using introgression lines]

The genetic map of chromosome 5B has been constructed by using microsatellite (SSR) analysis of 381 plants from the F2 population produced by cross of the Chinese Spring (CS) and Renan cultivars. Initially, 180 SSR markers for the common wheat 5B chromosome have been used for analysis of these cultivars. The 32 markers able to detect polymorphism between these cultivars have been located on the genetic map of chromosome 5B. Cytogenetic mapping has involved a set of CS 5B chromosome deletion lines. Totally, 51 SSR markers have been located in ten regions (deletion bins) of this chromosome by SSR analysis of these deletion lines. Five genes--TaCBFIIIc-B10, Vrn--B1, Chi--B1, Skr, and Ph1--have been integrated into the cytogenetic map of chromosome 5B using the markers either specific of or tightly linked to the genes in question. Com- parison of the genetic and cytogenetic maps suggests that recombination is suppressed in the pericentromeric region of chromosome 5B, especially in the short arm segment. The 18 markers localized to deletion bins 5BL16-0.79-1.00 and 5BL8-0.66-0.79 have been used to analyze common wheat introgression lines L842, L5366-180, L73/00i, and L21-4, carrying fragments of alien genomes in the terminal region of 5B long arm. L5366-180 and L842 lines carry a fragment of the Triticum timopheevii 5GL chromosome, while L73/00i ? L21-4 lines, a fragment of the Aegilopsspeltoides 5SL chromosome. As has been shown, the translocated fragments in these four lines are of different lengths, allowing bin 5BL18-0.66-0.79 to be divided into three shorter regions. The utility of wheat introgression lines carrying alien translocations for increasing the resolution of cytogenetic mapping is discussed.

A physical map of the short arm of wheat chromosome 1A

Bread wheat (Triticum aestivum) has a large and highly repetitive genome which poses major technical challenges for its study. To aid map-based cloning and future genome sequencing projects, we constructed a BAC-based physical map of the short arm of wheat chromosome 1A (1AS). From the assembly of 25,918 high information content (HICF) fingerprints from a 1AS-specific BAC library, 715 physical contigs were produced that cover almost 99% of the estimated size of the chromosome arm. The 3,414 BAC clones constituting the minimum tiling path were end-sequenced. Using a gene microarray containing ∼40 K NCBI UniGene EST clusters, PCR marker screening and BAC end sequences, we arranged 160 physical contigs (97 Mb or 35.3% of the chromosome arm) in a virtual order based on synteny with Brachypodium, rice and sorghum. BAC end sequences and information from microarray hybridisation was used to anchor 3.8 Mbp of Illumina sequences from flow-sorted chromosome 1AS to BAC contigs. Comparison of genetic and synteny-based physical maps indicated that ∼50% of all genetic recombination is confined to 14% of the physical length of the chromosome arm in the distal region. The 1AS physical map provides a framework for future genetic mapping projects as well as the basis for complete sequencing of chromosome arm 1AS.

Sequence-based analysis of translocations and inversions in bread wheat (Triticum aestivum L.)

Structural changes of chromosomes are a primary mechanism of genome rearrangement over the course of evolution and detailed knowledge of such changes in a given species and its close relatives should increase the efficiency and precision of chromosome engineering in crop improvement. We have identified sequences bordering each of the main translocation and inversion breakpoints on chromosomes 4A, 5A and 7B of the modern bread wheat genome. The locations of these breakpoints allow, for the first time, a detailed description of the evolutionary origins of these chromosomes at the gene level. Results from this study also demonstrate that, although the strategy of exploiting sorted chromosome arms has dramatically simplified the efforts of wheat genome sequencing, simultaneous analysis of sequences from homoeologous and non-homoeologous chromosomes is essential in understanding the origins of DNA sequences in polyploid species.

Radiation hybrid QTL mapping of Tdes2 involved in the first meiotic division of wheat

Since the dawn of wheat cytogenetics, chromosome 3B has been known to harbor a gene(s) that, when removed, causes chromosome desynapsis and gametic sterility. The lack of natural genetic diversity for this gene(s) has prevented any attempt to fine map and further characterize it. Here, gamma radiation treatment was used to create artificial diversity for this locus. A total of 696 radiation hybrid lines were genotyped with a custom mini array of 140 DArT markers, selected to evenly span the whole 3B chromosome. The resulting map spanned 2,852 centi Ray with a calculated resolution of 0.384 Mb. Phenotyping for the occurrence of meiotic desynapsis was conducted by measuring the level of gametic sterility as seeds produced per spikelet and pollen viability at booting. Composite interval mapping revealed a single QTL with LOD of 16.2 and r (2) of 25.6 % between markers wmc326 and wPt-8983 on the long arm of chromosome 3B. By independent analysis, the location of the QTL was confirmed to be within the deletion bin 3BL7-0.63-1.00 and to correspond to a single gene located ~1.4 Mb away from wPt-8983. The meiotic behavior of lines lacking this gene was characterized cytogenetically to reveal striking similarities with mutants for the dy locus, located on the syntenic chromosome 3 of maize. This represents the first example to date of employing radiation hybrids for QTL analysis. The success achieved by this approach provides an ideal starting point for the final cloning of this interesting gene involved in meiosis of cereals.

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

BACKGROUND

As for other major crops, achieving a complete wheat genome sequence is essential for the application of genomics to breeding new and improved varieties. To overcome the complexities of the large, highly repetitive and hexaploid wheat genome, the International Wheat Genome Sequencing Consortium established a chromosome-based strategy that was validated by the construction of the physical map of chromosome 3B. Here, we present improved strategies for the construction of highly integrated and ordered wheat physical maps, using chromosome 1BL as a template, and illustrate their potential for evolutionary studies and map-based cloning.

RESULTS

Using a combination of novel high throughput marker assays and an assembly program, we developed a high quality physical map representing 93% of wheat chromosome 1BL, anchored and ordered with 5,489 markers including 1,161 genes. Analysis of the gene space organization and evolution revealed that gene distribution and conservation along the chromosome results from the superimposition of the ancestral grass and recent wheat evolutionary patterns, leading to a peak of synteny in the central part of the chromosome arm and an increased density of non-collinear genes towards the telomere. With a density of about 11 markers per Mb, the 1BL physical map provides 916 markers, including 193 genes, for fine mapping the 40 QTLs mapped on this chromosome.

CONCLUSIONS

Here, we demonstrate that high marker density physical maps can be developed in complex genomes such as wheat to accelerate map-based cloning, gain new insights into genome evolution, and provide a foundation for reference sequencing.

Multi-environment analysis and improved mapping of a yield-related QTL on chromosome 3B of wheat

Improved mapping, multi-environment quantitative trait loci (QTL) analysis and dissection of allelic effects were used to define a QTL associated with grain yield, thousand grain weight and early vigour on chromosome 3BL of bread wheat (Triticum aestivum L.) under abiotic stresses. The QTL had pleiotropic effects and showed QTL x environment interactions across 21 diverse environments in Australia and Mexico. The occurrence and the severity of water deficit combined with high temperatures during the growing season affected the responsiveness of this QTL, resulting in a reversal in the direction of allelic effects. The influence of this QTL can be substantial, with the allele from one parent (RAC875) increasing grain yield by up to 12.5 % (particularly in environments where both heat and drought stress occurred) and the allele from the other parent (Kukri) increasing grain yield by up to 9 % in favourable environments. With the application of additional markers and the genotyping of additional recombinant inbred lines, the genetic map in the QTL region was refined to provide a basis for future positional cloning.

Wheat centromeric retrotransposons: the new ones take a major role in centromeric structure

The physical map of the hexaploid wheat chromosome 3B was screened using centromeric DNA probes. A 1.1-Mb region showing the highest number of positive bacterial artificial chromosome (BAC) clones was fully sequenced and annotated, revealing that 96% of the DNA consisted of transposable elements, mainly long terminal repeat (LTR) retrotransposons (88%). Estimation of the insertion time of the transposable elements revealed that CRW (also called Cereba) and Quinta are the youngest elements at the centromeres of common wheat (Triticum spp.) and its diploid ancestors, with Quinta being younger than CRW in both diploid and hexaploid wheats. Chromatin immunoprecipitation experiments revealed that both CRW and Quinta families are targeted by the centromere-specific histone H3 variant CENH3. Immuno colocalization of retroelements and CENH3 antibody indicated that a higher proportion of Quinta than CRWs was associated with CENH3, although CRWs were more abundant. Long arrays of satellite repeats were also identified in the wheat centromere regions, but they lost the ability to bind with CENH3. In addition to transposons, two functional genes and one pseudogene were identified. The gene density in the centromere appeared to be between three and four times lower than the average gene density of chromosome 3B. Comparisons with related grasses also indicated a loss of microcollinearity in this region. Finally, comparison of centromeric sequences of Aegilops tauschii (DD), Triticum boeoticum (AA) and hexaploid wheat revealed that the centromeres in both the polyploids and diploids are still undergoing dynamic changes, and that the new CRWs and Quintas may have undertaken the core role in kinetochore formation.

Integrating cereal genomics to support innovation in the Triticeae

The genomic resources of small grain cereals that include some of the most important crop species such as wheat, barley, and rye are attaining a level of completion that now is contributing to new structural and functional studies as well as refining molecular marker development and mapping strategies for increasing the efficiency of breeding processes. The integration of new efforts to obtain reference sequences in bread wheat and barley, in particular, is accelerating the acquisition and interpretation of genome-level analyses in both of these major crops.

Integrating cereal genomics to support innovation in the Triticeae

The genomic resources of small grain cereals that include some of the most important crop species such as wheat, barley, and rye are attaining a level of completion that now is contributing to new structural and functional studies as well as refining molecular marker development and mapping strategies for increasing the efficiency of breeding processes. The integration of new efforts to obtain reference sequences in bread wheat and barley, in particular, is accelerating the acquisition and interpretation of genome-level analyses in both of these major crops.

Association mapping and haplotype analysis of a 3.1-Mb genomic region involved in Fusarium head blight resistance on wheat chromosome 3BS

A previous study provided an in-depth understanding of molecular population genetics of European and Asian wheat gene pools using a sequenced 3.1-Mb contig (ctg954) on chromosome 3BS. This region is believed to carry the Fhb1 gene for response to Fusarium head blight. In this study, 266 wheat accessions were evaluated in three environments for Type II FHB response based on the single floret inoculation method. Hierarchical clustering (UPGMA) based on a Manhattan dissimilarity matrix divided the accessions into eight groups according to five FHB-related traits which have a high correlation between them; Group VIII comprised six accessions with FHB response levels similar to variety Sumai 3. Based on the compressed mixed linear model (MLM), association analysis between five FHB-related traits and 42 molecular markers along the 3.1-Mb region revealed 12 significant association signals at a threshold of P<0.05. The highest proportion of phenotypic variation (6.2%) in number of diseased spikelets (NDS) occurred at locus cfb6059, and the physical distance was about 2.9 Kb between umn10 and this marker. Haplotype block (HapB) analysis using a sliding window LD of 5 markers, detected six HapBs in the 3.1-Mb region at r(2)>0.1 and P<0.001 between random closely linked markers. F-tests among Haps with frequencies >0.05 within each HapB at r(2)>0.1 and P<0.001 showed significant differences between the Hap carried by FHB resistant resources, such as Sumai 3 and Wangshuibai, and susceptible genotypes in HapB3 and HapB6. These results suggest that Fhb1 is located within HapB6, with the possibility that another gene is located at or near HapB3. SSR markers and Haps detected in this study will be helpful in further understanding the genetic basis of FHB resistance, and provide useful information for marker-assisted selection of Fhb1 in wheat breeding.

Intraspecific sequence comparisons reveal similar rates of non-collinear gene insertion in the B and D genomes of bread wheat

BACKGROUND

Polyploidization is considered one of the main mechanisms of plant genome evolution. The presence of multiple copies of the same gene reduces selection pressure and permits sub-functionalization and neo-functionalization leading to plant diversification, adaptation and speciation. In bread wheat, polyploidization and the prevalence of transposable elements resulted in massive gene duplication and movement. As a result, the number of genes which are non-collinear to genomes of related species seems markedly increased in wheat.

RESULTS

We used new-generation sequencing (NGS) to generate sequence of a Mb-sized region from wheat chromosome arm 3DS. Sequence assembly of 24 BAC clones resulted in two scaffolds of 1,264,820 and 333,768 bases. The sequence was annotated and compared to the homoeologous region on wheat chromosome 3B and orthologous loci of Brachypodium distachyon and rice. Among 39 coding sequences in the 3DS scaffolds, 32 have a homoeolog on chromosome 3B. In contrast, only fifteen and fourteen orthologs were identified in the corresponding regions in rice and Brachypodium, respectively. Interestingly, five pseudogenes were identified among the non-collinear coding sequences at the 3B locus, while none was found at the 3DS locus.

CONCLUSION

Direct comparison of two Mb-sized regions of the B and D genomes of bread wheat revealed similar rates of non-collinear gene insertion in both genomes with a majority of gene duplications occurring before their divergence. Relatively low proportion of pseudogenes was identified among non-collinear coding sequences. Our data suggest that the pseudogenes did not originate from insertion of non-functional copies, but were formed later during the evolution of hexaploid wheat. Some evidence was found for gene erosion along the B genome locus.

Genome-level identification of cell wall invertase genes in wheat for the study of drought tolerance

In wheat (Triticum aestivum L.) drought-induced pollen sterility is a major contributor to grain yield loss and is caused by the downregulation of the cell wall invertase gene IVR1. The IVR1 gene catalyses the irreversible hydrolysis of sucrose to glucose and fructose, the essential energy substrates which support pollen development. Downregulation of IVR1 in response to drought is isoform specific and shows variation in temporal and tissue-specific expression. IVR1 is now prompting interest as a candidate gene for molecular marker development to screen wheat germplasm for improved drought tolerance. The aim of this study was to define the family of IVR1 genes to enable: (1) individual isoforms to be assayed in gene expression studies; and (2) greater accuracy in IVR1 mapping to the wheat genetic map and drought tolerance QTL analysis. Using a cell wall invertase-specific motif as a probe, wheat genomics platforms were screened for the presence of unidentified IVR1 isoforms. Wheat genomics platforms screened included the IWGSC wheat survey sequence, the wheat D genome donor sequence from Aegilops tauschii Coss, and the CCG wheat chromosome 3B assembly: contig506. Chromosome-specific sequences homologous to the query motif were isolated and characterised. Sequence annotation results showed five previously unidentified IVR1 isoforms exist on multiple chromosome arms and on all three genomes (A, B and D): IVR1-3A, IVR1-4A, IVR1-5B, IVR1.2-3B and IVR1-5D. Including three previously characterised IVR1 isoforms (IVR1.1-1A, IVR1.2-1A and IVR1.1-3B), the total number of isoform gene family members is eight. The IVR1 isoforms contain two motifs common to cell wall invertase (NDPN and WECPDF) and a high degree of conservation in exon 4, suggesting conservation of functionality. Sequence divergence at a primary structure level in other regions of the gene was evident amongst the isoforms, which likely contributes to variation in gene regulation and expression in response to water deficit within this subfamily of IVR1 isoforms in wheat.

DNA repair and crossing over favor similar chromosome regions as discovered in radiation hybrid of Triticum

BACKGROUND

The uneven distribution of recombination across the length of chromosomes results in inaccurate estimates of genetic to physical distances. In wheat (Triticum aestivum L.) chromosome 3B, it has been estimated that 90% of the cross over events occur in distal sub-telomeric regions representing 40% of the chromosome. Radiation hybrid (RH) mapping which does not rely on recombination is a strategy to map genomes and has been widely employed in animal species and more recently in some plants. RH maps have been proposed to provide i) higher and ii) more uniform resolution than genetic maps, and iii) to be independent of the distribution patterns observed for meiotic recombination. An in vivo RH panel was generated for mapping chromosome 3B of wheat in an attempt to provide a complete scaffold for this ~1 Gb segment of the genome and compare the resolution to previous genetic maps.

RESULTS

A high density RH map with 541 marker loci anchored to chromosome 3B spanning a total distance of 1871.9 cR was generated. Detailed comparisons with a genetic map of similar quality confirmed that i) the overall resolution of the RH map was 10.5 fold higher and ii) six fold more uniform. A significant interaction (r = 0.879 at p = 0.01) was observed between the DNA repair mechanism and the distribution of crossing-over events. This observation could be explained by accepting the possibility that the DNA repair mechanism in somatic cells is affected by the chromatin state in a way similar to the effect that chromatin state has on recombination frequencies in gametic cells.

CONCLUSIONS

The RH data presented here support for the first time in vivo the hypothesis of non-casual interaction between recombination hot-spots and DNA repair. Further, two major hypotheses are presented on how chromatin compactness could affect the DNA repair mechanism. Since the initial RH application 37 years ago, we were able to show for the first time that the iii) third hypothesis of RH mapping might not be entirely correct.

DNA repair and crossing over favor similar chromosome regions as discovered in radiation hybrid of Triticum

BACKGROUND

The uneven distribution of recombination across the length of chromosomes results in inaccurate estimates of genetic to physical distances. In wheat (Triticum aestivum L.) chromosome 3B, it has been estimated that 90% of the cross over events occur in distal sub-telomeric regions representing 40% of the chromosome. Radiation hybrid (RH) mapping which does not rely on recombination is a strategy to map genomes and has been widely employed in animal species and more recently in some plants. RH maps have been proposed to provide i) higher and ii) more uniform resolution than genetic maps, and iii) to be independent of the distribution patterns observed for meiotic recombination. An in vivo RH panel was generated for mapping chromosome 3B of wheat in an attempt to provide a complete scaffold for this ~1 Gb segment of the genome and compare the resolution to previous genetic maps.

RESULTS

A high density RH map with 541 marker loci anchored to chromosome 3B spanning a total distance of 1871.9 cR was generated. Detailed comparisons with a genetic map of similar quality confirmed that i) the overall resolution of the RH map was 10.5 fold higher and ii) six fold more uniform. A significant interaction (r = 0.879 at p = 0.01) was observed between the DNA repair mechanism and the distribution of crossing-over events. This observation could be explained by accepting the possibility that the DNA repair mechanism in somatic cells is affected by the chromatin state in a way similar to the effect that chromatin state has on recombination frequencies in gametic cells.

CONCLUSIONS

The RH data presented here support for the first time in vivo the hypothesis of non-casual interaction between recombination hot-spots and DNA repair. Further, two major hypotheses are presented on how chromatin compactness could affect the DNA repair mechanism. Since the initial RH application 37 years ago, we were able to show for the first time that the iii) third hypothesis of RH mapping might not be entirely correct.

Whole Genome Profiling provides a robust framework for physical mapping and sequencing in the highly complex and repetitive wheat genome

Background

Sequencing projects using a clone-by-clone approach require the availability of a robust physical map. The SNaPshot technology, based on pair-wise comparisons of restriction fragments sizes, has been used recently to build the first physical map of a wheat chromosome and to complete the maize physical map. However, restriction fragments sizes shared randomly between two non-overlapping BACs often lead to chimerical contigs and mis-assembled BACs in such large and repetitive genomes. Whole Genome Profiling (WGP™) was developed recently as a new sequence-based physical mapping technology and has the potential to limit this problem.

Results

A subset of the wheat 3B chromosome BAC library covering 230 Mb was used to establish a WGP physical map and to compare it to a map obtained with the SNaPshot technology. We first adapted the WGP-based assembly methodology to cope with the complexity of the wheat genome. Then, the results showed that the WGP map covers the same length than the SNaPshot map but with 30% less contigs and, more importantly with 3.5 times less mis-assembled BACs. Finally, we evaluated the benefit of integrating WGP tags in different sequence assemblies obtained after Roche/454 sequencing of BAC pools. We showed that while WGP tag integration improves assemblies performed with unpaired reads and with paired-end reads at low coverage, it does not significantly improve sequence assemblies performed at high coverage (25x) with paired-end reads.

Conclusions

Our results demonstrate that, with a suitable assembly methodology, WGP builds more robust physical maps than the SNaPshot technology in wheat and that WGP can be adapted to any genome. Moreover, WGP tag integration in sequence assemblies improves low quality assembly. However, to achieve a high quality draft sequence assembly, a sequencing depth of 25x paired-end reads is required, at which point WGP tag integration does not provide additional scaffolding value. Finally, we suggest that WGP tags can support the efficient sequencing of BAC pools by enabling reliable assignment of sequence scaffolds to their BAC of origin, a feature that is of great interest when using BAC pooling strategies to reduce the cost of sequencing large genomes.

TriAnnot: A Versatile and High Performance Pipeline for the Automated Annotation of Plant Genomes

In support of the international effort to obtain a reference sequence of the bread wheat genome and to provide plant communities dealing with large and complex genomes with a versatile, easy-to-use online automated tool for annotation, we have developed the TriAnnot pipeline. Its modular architecture allows for the annotation and masking of transposable elements, the structural, and functional annotation of protein-coding genes with an evidence-based quality indexing, and the identification of conserved non-coding sequences and molecular markers. The TriAnnot pipeline is parallelized on a 712 CPU computing cluster that can run a 1-Gb sequence annotation in less than 5 days. It is accessible through a web interface for small scale analyses or through a server for large scale annotations. The performance of TriAnnot was evaluated in terms of sensitivity, specificity, and general fitness using curated reference sequence sets from rice and wheat. In less than 8 h, TriAnnot was able to predict more than 83% of the 3,748 CDS from rice chromosome 1 with a fitness of 67.4%. On a set of 12 reference Mb-sized contigs from wheat chromosome 3B, TriAnnot predicted and annotated 93.3% of the genes among which 54% were perfectly identified in accordance with the reference annotation. It also allowed the curation of 12 genes based on new biological evidences, increasing the percentage of perfect gene prediction to 63%. TriAnnot systematically showed a higher fitness than other annotation pipelines that are not improved for wheat. As it is easily adaptable to the annotation of other plant genomes, TriAnnot should become a useful resource for the annotation of large and complex genomes in the future.

A 3,000-loci transcription map of chromosome 3B unravels the structural and functional features of gene islands in hexaploid wheat

To improve our understanding of the organization and regulation of the wheat (Triticum aestivum) gene space, we established a transcription map of a wheat chromosome (3B) by hybridizing a newly developed wheat expression microarray with bacterial artificial chromosome pools from a new version of the 3B physical map as well as with cDNA probes derived from 15 RNA samples. Mapping data for almost 3,000 genes showed that the gene space spans the whole chromosome 3B with a 2-fold increase of gene density toward the telomeres due to an increase in the number of genes in islands. Comparative analyses with rice (Oryza sativa) and Brachypodium distachyon revealed that these gene islands are composed mainly of genes likely originating from interchromosomal gene duplications. Gene Ontology and expression profile analyses for the 3,000 genes located along the chromosome revealed that the gene islands are enriched significantly in genes sharing the same function or expression profile, thereby suggesting that genes in islands acquired shared regulation during evolution. Only a small fraction of these clusters of cofunctional and coexpressed genes was conserved with rice and B. distachyon, indicating a recent origin. Finally, genes with the same expression profiles in remote islands (coregulation islands) were identified suggesting long-distance regulation of gene expression along the chromosomes in wheat.

Sequence-based marker development in wheat: advances and applications to breeding

In the past two decades, the wheat community has made remarkable progress in developing molecular resources for breeding. A wide variety of molecular tools has been established to accelerate genetic and physical mapping for facilitating the efficient identification of molecular markers linked to genes and QTL of agronomic interest. Already, wheat breeders are benefiting from a wide range of techniques to follow the introgression of the most favorable alleles in elite material and develop improved varieties. Breeders soon will be able to take advantage of new technological developments based on Next Generation Sequencing. In this paper, we review the molecular toolbox available to wheat scientists and breeders for performing fundamental genomic studies and breeding. Special emphasis is given on the production and detection of single nucleotide polymorphisms (SNPs) that should enable a step change in saturating the wheat genome for more efficient genetic studies and for the development of new selection methods. The perspectives offered by the access to an ordered full genome sequence for further marker development and enhanced precision breeding is also discussed. Finally, we discuss the advantages and limitations of marker-assisted selection for supporting wheat improvement.

Frequent gene movement and pseudogene evolution is common to the large and complex genomes of wheat, barley, and their relatives

All six arms of the group 1 chromosomes of hexaploid wheat (Triticum aestivum) were sequenced with Roche/454 to 1.3- to 2.2-fold coverage and compared with similar data sets from the homoeologous chromosome 1H of barley (Hordeum vulgare). Six to ten thousand gene sequences were sampled per chromosome. These were classified into genes that have their closest homologs in the Triticeae group 1 syntenic region in Brachypodium, rice (Oryza sativa), and/or sorghum (Sorghum bicolor) and genes that have their homologs elsewhere in these model grass genomes. Although the number of syntenic genes was similar between the homologous groups, the amount of nonsyntenic genes was found to be extremely diverse between wheat and barley and even between wheat subgenomes. Besides a small core group of genes that are nonsyntenic in other grasses but conserved among Triticeae, we found thousands of genic sequences that are specific to chromosomes of one single species or subgenome. By examining in detail 50 genes from chromosome 1H for which BAC sequences were available, we found that many represent pseudogenes that resulted from transposable element activity and double-strand break repair. Thus, Triticeae seem to accumulate nonsyntenic genes frequently. Since many of them are likely to be pseudogenes, total gene numbers in Triticeae are prone to pronounced overestimates.

Cross-genome map based dissection of a nitrogen use efficiency ortho-metaQTL in bread wheat unravels concerted cereal genome evolution

Monitoring nitrogen use efficiency (NUE) in plants is becoming essential to maintain yield while reducing fertilizer usage. Optimized NUE application in major crops is essential for long-term sustainability of agriculture production. Here, we report the precise identification of 11 major chromosomal regions controlling NUE in wheat that co-localise with key developmental genes such as Ppd (photoperiod sensitivity), Vrn (vernalization requirement), Rht (reduced height) and can be considered as robust markers from a molecular breeding perspective. Physical mapping, sequencing, annotation and candidate gene validation of an NUE metaQTL on wheat chromosome 3B allowed us to propose that a glutamate synthase (GoGAT) gene that is conserved structurally and functionally at orthologous positions in rice, sorghum and maize genomes may contribute to NUE in wheat and other cereals. We propose an evolutionary model for the NUE locus in cereals from a common ancestral region, involving species specific shuffling events such as gene deletion, inversion, transposition and the invasion of repetitive elements.

Crop genome sequencing: lessons and rationales

2010 marks the 10th anniversary of the completion of the first plant genome sequence (Arabidopsis thaliana). Triggered by advancements in sequencing technologies, many crop genome sequences have been produced, with eight published since 2008. To date, however, only the rice (Oryza sativa) genome sequence has been finished to a quality level similar to that of the Arabidopsis sequence. This trend to produce draft genomes could affect the ability of researchers to address biological questions of speciation and recent evolution or to link sequence variation accurately to phenotypes. Here, we review the current crop genome sequencing activities, discuss how variability in sequence quality impacts utility for different studies and provide a perspective for a paradigm shift in selecting crops for sequencing in the future.

Palaeogenomics in cereals: modeling of ancestors for modern species improvement

During the last decade, technological improvements led to the development of large sets of plant genomic resources permitting the emergence of high-resolution comparative genomic studies. Synteny-based identification of seven shared duplications in cereals led to the modeling of a common ancestral genome structure of 33.6 Mb structured in five protochromosomes containing 9138 protogenes and provided new insights into the evolution of cereal genomes from their extinct ancestors. Recent palaeogenomic data indicate that whole genome duplications were a driving force in the evolutionary success of cereals over the last 50 to 70 millions years. Finally, detailed synteny and duplication relationships led to an improved representation of cereal genomes in concentric circles, thus providing a new reference tool for improved gene annotation and cross-genome markers development.

Considering transposable element diversification in de novo annotation approaches

Transposable elements (TEs) are mobile, repetitive DNA sequences that are almost ubiquitous in prokaryotic and eukaryotic genomes. They have a large impact on genome structure, function and evolution. With the recent development of high-throughput sequencing methods, many genome sequences have become available, making possible comparative studies of TE dynamics at an unprecedented scale. Several methods have been proposed for the de novo identification of TEs in sequenced genomes. Most begin with the detection of genomic repeats, but the subsequent steps for defining TE families differ. High-quality TE annotations are available for the Drosophila melanogaster and Arabidopsis thaliana genome sequences, providing a solid basis for the benchmarking of such methods. We compared the performance of specific algorithms for the clustering of interspersed repeats and found that only a particular combination of algorithms detected TE families with good recovery of the reference sequences. We then applied a new procedure for reconciling the different clustering results and classifying TE sequences. The whole approach was implemented in a pipeline using the REPET package. Finally, we show that our combined approach highlights the dynamics of well defined TE families by making it possible to identify structural variations among their copies. This approach makes it possible to annotate TE families and to study their diversification in a single analysis, improving our understanding of TE dynamics at the whole-genome scale and for diverse species.

Specific patterns of gene space organisation revealed in wheat by using the combination of barley and wheat genomic resources

Background

Because of its size, allohexaploid nature and high repeat content, the wheat genome has always been perceived as too complex for efficient molecular studies. We recently constructed the first physical map of a wheat chromosome (3B). However gene mapping is still laborious in wheat because of high redundancy between the three homoeologous genomes. In contrast, in the closely related diploid species, barley, numerous gene-based markers have been developed. This study aims at combining the unique genomic resources developed in wheat and barley to decipher the organisation of gene space on wheat chromosome 3B.

Results

Three dimensional pools of the minimal tiling path of wheat chromosome 3B physical map were hybridised to a barley Agilent 15K expression microarray. This led to the fine mapping of 738 barley orthologous genes on wheat chromosome 3B. In addition, comparative analyses revealed that 68% of the genes identified were syntenic between the wheat chromosome 3B and barley chromosome 3 H and 59% between wheat chromosome 3B and rice chromosome 1, together with some wheat-specific rearrangements. Finally, it indicated an increasing gradient of gene density from the centromere to the telomeres positively correlated with the number of genes clustered in islands on wheat chromosome 3B.

Conclusion

Our study shows that novel structural genomics resources now available in wheat and barley can be combined efficiently to overcome specific problems of genetic anchoring of physical contigs in wheat and to perform high-resolution comparative analyses with rice for deciphering the organisation of the wheat gene space.

Variation in crossover rates across a 3-Mb contig of bread wheat (Triticum aestivum) reveals the presence of a meiotic recombination hotspot

In bread wheat (Triticum aestivum L.), initial studies using deletion lines indicated that crossover (CO) events occur mainly in the telomeric regions of the chromosomes with a possible correlation with the presence of genes. However, little is known about the distribution of COs at the sequence level. To investigate this, we studied in detail the pattern of COs along a contig of 3.110 Mb using two F2 segregating populations (Chinese Spring × Renan (F2-CsRe) and Chinese Spring × Courtot (F2-CsCt)) each containing ~2,000 individuals. The availability of the sequence of the contig from Cs enabled the development of 318 markers among which 23 co-dominant polymorphic markers (11 SSRs and 12 SNPs) were selected for CO distribution analyses. The distribution of CO events was not homogeneous throughout the contig, ranging from 0.05 to 2.77 cM/Mb, but was conserved between the two populations despite very different contig recombination rate averages (0.82 cM/Mb in F2-CsRe vs 0.35 cM/Mb in F2-CsCt). The CO frequency was correlated with the percentage of coding sequence in Cs and with the polymorphism rate between Cs and Re or Ct in both populations, indicating an impact of these two factors on CO distribution. At a finer scale, COs were found in a region covering 2.38 kb, spanning a gene coding for a glycosyl transferase (Hga3), suggesting the presence of a CO hotspot. A non-crossover event covering at least 453 bp was also identified in the same interval. From these results, we can conclude that gene content could be one of the factors driving recombination in bread wheat.

LTC: a novel algorithm to improve the efficiency of contig assembly for physical mapping in complex genomes

Background

Physical maps are the substrate of genome sequencing and map-based cloning and their construction relies on the accurate assembly of BAC clones into large contigs that are then anchored to genetic maps with molecular markers. High Information Content Fingerprinting has become the method of choice for large and repetitive genomes such as those of maize, barley, and wheat. However, the high level of repeated DNA present in these genomes requires the application of very stringent criteria to ensure a reliable assembly with the FingerPrinted Contig (FPC) software, which often results in short contig lengths (of 3-5 clones before merging) as well as an unreliable assembly in some difficult regions. Difficulties can originate from a non-linear topological structure of clone overlaps, low power of clone ordering algorithms, and the absence of tools to identify sources of gaps in Minimal Tiling Paths (MTPs).

Results

To address these problems, we propose a novel approach that: (i) reduces the rate of false connections and Q-clones by using a new cutoff calculation method; (ii) obtains reliable clusters robust to the exclusion of single clone or clone overlap; (iii) explores the topological contig structure by considering contigs as networks of clones connected by significant overlaps; (iv) performs iterative clone clustering combined with ordering and order verification using re-sampling methods; and (v) uses global optimization methods for clone ordering and Band Map construction. The elements of this new analytical framework called Linear Topological Contig (LTC) were applied on datasets used previously for the construction of the physical map of wheat chromosome 3B with FPC. The performance of LTC vs. FPC was compared also on the simulated BAC libraries based on the known genome sequences for chromosome 1 of rice and chromosome 1 of maize.

Conclusions

The results show that compared to other methods, LTC enables the construction of highly reliable and longer contigs (5-12 clones before merging), the detection of "weak" connections in contigs and their "repair", and the elongation of contigs obtained by other assembly methods.

Genetic structure and ecogeographical adaptation in wild barley (Hordeum chilense Roemer et Schultes) as revealed by microsatellite markers

BACKGROUND

Multi-allelic microsatellite markers have become the markers of choice for the determination of genetic structure in plants. Synteny across cereals has allowed the cross-species and cross-genera transferability of SSR markers, which constitute a valuable and cost-effective tool for the genetic analysis and marker-assisted introgression of wild related species. Hordeum chilense is one of the wild relatives with a high potential for cereal breeding, due to its high crossability (both interspecies and intergenera) and polymorphism for adaptation traits. In order to analyze the genetic structure and ecogeographical adaptation of this wild species, it is necessary to increase the number of polymorphic markers currently available for the species. In this work, the possibility of using syntenic wheat SSRs as a new source of markers for this purpose has been explored.

RESULTS

From the 98 wheat EST-SSR markers tested for transferability and polymorphism in the wild barley genome, 53 primer pairs (54.0%) gave cross-species transferability and 20 primer pairs (20.4%) showed polymorphism. The latter were used for further analysis in the H. chilense germplasm. The H. chilense-Triticum aestivum addition lines were used to test the chromosomal location of the new polymorphic microsatellite markers. The genetic structure and diversity was investigated in a collection of 94 H. chilense accessions, using a set of 49 SSR markers distributed across the seven chromosomes. Microsatellite markers showed a total of 351 alleles over all loci. The number of alleles per locus ranged from two to 27, with a mean of 7.2 alleles per locus and a mean Polymorphic Information Content (PIC) of 0.5.

CONCLUSIONS

According to the results, the germplasm can be divided into two groups, with morphological and ecophysiological characteristics being key determinants of the population structure. Geographic and ecological structuring was also revealed in the analyzed germplasm. A significant correlation between geographical and genetic distance was detected in the Central Chilean region for the first time in the species. In addition, significant ecological influence in genetic distance has been detected for one of the population structure groups (group II) in the Central Chilean region. Finally, the association of the SSR markers with ecogeographical variables was investigated and one marker was found significantly associated with precipitation. These findings have a potential application in cereal breeding.