The coordinated regulation of early meiotic stages is dominated by non-coding RNAs and stage-specific transcription in wheat.

Reproductive success hinges on precisely coordinated meiosis, yet our understanding of how structural rearrangements of chromatin and phase transitions during meiosis are transcriptionally regulated is limited. In crop plants, detailed analysis of the meiotic transcriptome could identify regulatory genes and epigenetic regulators that can be targeted to increase recombination rates and broaden genetic variation, as well as provide a resource for comparison among eukaryotes of different taxa to answer outstanding questions about meiosis. We conducted a meiotic stage-specific analysis of messenger RNA (mRNA), small non-coding RNA (sncRNA), and long intervening/intergenic non-coding RNA (lincRNA) in wheat (Triticum aestivum L.) and revealed novel mechanisms of meiotic transcriptional regulation and meiosis-specific transcripts. Amidst general repression of mRNA expression, significant enrichment of ncRNAs was identified during prophase I relative to vegetative cells. The core meiotic transcriptome was comprised of 9309 meiosis-specific transcripts, 48 134 previously unannotated meiotic transcripts, and many known and novel ncRNAs differentially expressed at specific stages. The abundant meiotic sncRNAs controlled the reprogramming of central metabolic pathways by targeting genes involved in photosynthesis, glycolysis, hormone biosynthesis, and cellular homeostasis, and lincRNAs enhanced the expression of nearby genes. Alternative splicing was not evident in this polyploid species, but isoforms were switched at phase transitions. The novel, stage-specific regulatory controls uncovered here challenge the conventional understanding of this crucial biological process and provide a new resource of requisite knowledge for those aiming to directly modulate meiosis to improve crop plants. The wheat meiosis transcriptome dataset can be queried for genes of interest using an eFP browser located at https://bar.utoronto.ca/efp_wheat/cgi-bin/efpWeb.cgi?dataSource=Wheat_Meiosis.

Long-read sequence assembly: a technical evaluation in barley.

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

Sequencing of Camelina neglecta, a diploid progenitor of the hexaploid oilseed Camelina sativa.

Camelina neglecta is a diploid species from the genus Camelina, which includes the versatile oilseed Camelina sativa. These species are closely related to Arabidopsis thaliana and the economically important Brassica crop species, making this genus a useful platform to dissect traits of agronomic importance while providing a tool to study the evolution of polyploids. A highly contiguous chromosome-level genome sequence of C. neglecta with an N50 size of 29.1 Mb was generated utilizing Pacific Biosciences (PacBio, Menlo Park, CA) long-read sequencing followed by chromosome conformation phasing. Comparison of the genome with that of C. sativa shows remarkable coincidence with subgenome 1 of the hexaploid, with only one major chromosomal rearrangement separating the two. Synonymous substitution rate analysis of the predicted 34 061 genes suggested subgenome 1 of C. sativa directly descended from C. neglecta around 1.2 mya. Higher functional divergence of genes in the hexaploid as evidenced by the greater number of unique orthogroups, and differential composition of resistant gene analogs, might suggest an immediate adaptation strategy after genome merger. The absence of genome bias in gene fractionation among the subgenomes of C. sativa in comparison with C. neglecta, and the complete lack of fractionation of meiosis-specific genes attests to the neopolyploid status of C. sativa. The assembled genome will provide a tool to further study genome evolution processes in the Camelina genus and potentially allow for the identification and exploitation of novel variation for Camelina crop improvement.

Alternative splicing dynamics and evolutionary divergence during embryogenesis in wheat species.

Among polyploid species with complex genomic architecture, variations in the regulation of alternative splicing (AS) provide opportunities for transcriptional and proteomic plasticity and the potential for generating trait diversities. However, the evolution of AS and its influence on grain development in diploid grass and valuable polyploid wheat crops are poorly understood. To address this knowledge gap, we developed a pipeline for the analysis of alternatively spliced transcript isoforms, which takes the high sequence similarity among polyploid wheat subgenomes into account. Through analysis of synteny and detection of collinearity of homoeologous subgenomes, conserved and specific AS events across five wheat and grass species were identified. A global analysis of the regulation of AS in diploid grass and polyploid wheat grains revealed diversity in AS events not only between the endosperm, pericarp and embryo overdevelopment, but also between subgenomes. Analysis of AS in homoeologous triads of polyploid wheats revealed evolutionary divergence between gene-level and transcript-level regulation of embryogenesis. Evolutionary age analysis indicated that the generation of novel transcript isoforms has occurred in young genes at a more rapid rate than in ancient genes. These findings, together with the development of comprehensive AS resources for wheat and grass species, advance understanding of the evolution of regulatory features of AS during embryogenesis and grain development in wheat.

Structural features of two major nucleolar organizer regions (NORs), Nor-B1 and Nor-B2, and chromosome-specific rRNA gene expression in wheat.

The reference genome sequence of wheat 'Chinese Spring' (CS) is now available (IWGSC RefSeq v1.0), but the core sequences defining the nucleolar organizer regions (NORs) have not been characterized. We estimated that the total copy number of the rDNA units in the wheat genome is 11 160, of which 30.5%, 60.9% and 8.6% are located on Nor-B1 (1B), Nor-B2 (6B) and other NORs, respectively. The total length of the NORs is estimated to be 100 Mb, corresponding to approximately 10% of the unassembled portion of the genome not represented in RefSeq v1.0. Four subtypes (S1-S4) of the rDNA units were identified based on differences within the 3' external transcribed spacer regions in Nor-B1 and Nor-B2, and quantitative PCR indicated locus-specific variation in rDNA subtype contents. Expression analyses of rDNA subtypes revealed that S1 was predominantly expressed and S2 weakly expressed, in contrast to the relative abundance of rDNA subtypes in the wheat genome. These results suggest a regulation mechanism of differential rDNA expression based on sequence differences. S3 expression increased in the ditelosomic lines Dt1BL and Dt6BL, suggesting that S3 is subjected to chromosome-mediated silencing. Structural differences were detected in the regions surrounding the NOR among homoeologous chromosomes of groups 1 and 6. The adjacent regions distal to the major NORs were expanded compared with their homoeologous counterparts, and the gene density of these expanded regions was relatively low. We provide evidence that these regions are likely to be important for autoregulation of the associated major NORs as well as silencing of minor NORs.

WheatCRISPR: a web-based guide RNA design tool for CRISPR/Cas9-mediated genome editing in wheat.

BACKGROUND: CRISPR/Cas9 gene editing has become a revolutionary technique for crop improvement as it can facilitate fast and efficient genetic changes without the retention of transgene components in the final plant line. Lack of robust bioinformatics tools to facilitate the design of highly specific functional guide RNAs (gRNAs) and prediction of off-target sites in wheat is currently an obstacle to effective application of CRISPR technology to wheat improvement. DESCRIPTION: We have developed a web-based bioinformatics tool to design specific gRNAs for genome editing and transcriptional regulation of gene expression in wheat. A collaborative study between the Broad Institute and Microsoft Research used large-scale empirical evidence to devise algorithms (Doech et al., 2016, Nature Biotechnology 34, 184-191) for predicting the on-target activity and off-target potential of CRISPR/SpCas9 (Streptococcus pyogenes Cas9). We applied these prediction models to determine on-target specificity and potential off-target activity for individual gRNAs targeting specific loci in the wheat genome. The genome-wide gRNA mappings and the corresponding Doench scores predictive of the on-target and off-target activities were used to create a gRNA database which was used as a data source for the web application termed WheatCRISPR. CONCLUSION: The WheatCRISPR tool allows researchers to browse all possible gRNAs targeting a gene or sequence of interest and select effective gRNAs based on their predicted high on-target and low off-target activity scores, as well as other characteristics such as position within the targeted gene. It is publicly available at https://crispr.bioinfo.nrc.ca/WheatCrispr/ .

Correction to: MeioCapture: an efficient method for staging and isolation of meiocytes in the prophase I sub-stages of meiosis in wheat.

Following publication of the original article [1], a reader spotted an incorrect citation of the reference 14 [2] in the 'Background'. The male meiocyte isolation work described in this article [2] was carried out in rice and not in Brassica as originally stated in the 'Background' [1]. Thus, the following amendment to the Background section should be noted.

Mapping quantitative trait loci associated with common bunt resistance in a spring wheat (Triticum aestivum L.) variety Lillian.

KEY MESSAGE: Based on their consistency over environments, two QTL identified in Lillian on chromosomes 5A and 7A could be useful targets for marker assisted breeding of common bunt resistance. Common bunt of wheat (Triticum aestivum L.) caused by Tilletia tritici and T. laevis is an economically important disease because of losses in grain yield and reduced grain quality. Resistance can be quantitative, under the control of multiple small effect genes. The Canada Western Red Spring wheat variety Lillian is moderately resistant to common bunt races found on the Canadian prairies. This study was conducted to identify and map quantitative trait loci (QTL) conferring resistance against common bunt in Lillian. A doubled haploid population comprising 280 lines was developed from F1 plants of the cross of Lillian by Vesper. The lines were inoculated at seeding with the two races L16 (T. laevis) and T19 (T. tritici), grown in field near Swift Current, SK, in 2014, 2015 and 2016 and assessed for disease incidence. The lines were genotyped with the 90 K iSelect SNP genotyping assay, and a high-density genetic map was constructed. Quantitative trait locus analysis was performed with MapQTL.6® software. Two relatively stable common bunt resistance QTL, detected in two of the 3 years, were identified on chromosomes 5A and 7A from Lillian. In addition, three less stable QTL, appearing in one out of 3 years, were identified: one was contributed by Lillian on chromosome 3D and two were contributed by Vesper on chromosomes 1D and 2A. Epistatic interaction was identified for the bunt incidence between 3D and 7A resulting in greater bunt resistance. Future bunt resistance breeding will benefit from combining these QTL through gene pyramiding.

Genome re-sequencing and simple sequence repeat markers reveal the existence of divergent lineages in the Canadian Puccinia striiformis f. sp. tritici population with extensive DNA methylation.

Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is an important disease in Canada. The worldwide genetic structure of Pst populations have been characterized, excluding Canada. Here, we elucidated the genetic structure of the western Canadian Pst population using molecular markers, revealing the presence of four divergent lineages with predominantly clonal structure. In the worldwide context, two previously reported lineages were identified: PstS0 (22%), representing an old Northwestern-European and PstS1 (35%), an invasive warm-temperature adapted. Additionally, two new, unreported lineages, PstPr (9%) and PstS1-related (35%), were detected, which produced more telia than other lineages and had double the number of unique recombination events. The PstPr was a recent invasion, and likely evolved in a diverse, recombinant population as it was closely related to the PstS5, PstS7/Warrior, PstS8/Kranich, and PstS9 lineages originating from sexually recombining populations in the centre of diversity. The DNA methylation analysis revealed DNA-methyltransferase1-homologs, providing compelling evidence for epigenetic regulation and as a first report, an average of ∼5%, 5hmC in the Puccinia epigenome merits further investigation. The divergent lineages in the Canadian Pst population with the potential for genetic recombination, as well as epigenetic regulation needs consideration in the context of pathogen adaptation and management.

MeioCapture: an efficient method for staging and isolation of meiocytes in the prophase I sub-stages of meiosis in wheat.

BACKGROUND: Molecular analysis of meiosis has been hindered by difficulties in isolating high purity subpopulations of sporogenous cells representing the succeeding stages of meiosis. Isolation of purified male meiocytes from defined meiotic stages is crucial in discovering meiosis specific genes and associated regulatory networks. RESULTS: We describe an optimized method termed MeioCapture for simultaneous isolation of uncontaminated male meiocytes from wheat (Triticum spp.), specifically from the pre-meiotic G2 and the five sub-stages of meiotic prophase I. The MeioCapture protocol builds on the traditional anther squash technique and the capillary collection method, and involves extrusion of intact sporogenous archesporial columns (SACs) containing meiocytes. This improved method exploits the natural meiotic synchrony between anthers of the same floret, the correlation between the length of anthers and meiotic stage, and the occurrence of meiocytes in intact SACs largely free of somatic cells. The main advantage of MeioCapture, compared to previous methods, is that it allows simultaneous collection of meiocytes from different sub-stages of prophase I at a very high level of purity, through correlation of stages with anther sizes. A detailed description is provided for all steps, including the collection of tissue, isolation and size sorting of anthers, extrusion of intact SACs, and staging of meiocytes. Precautions for individual steps throughout the procedure are also provided to facilitate efficient isolation of pure meiocytes. The proof-of-concept was successfully established in wheat, and a light microscopic atlas of meiosis, encompassing all stages from pre-meiosis to telophase II, was developed. CONCLUSION: The MeioCapture method provides an essential technique to study the molecular basis of chromosome pairing and exchange of genetic information in wheat, leading to strategies for manipulating meiotic recombination frequencies. The method also provides a foundation for similar studies in other crop species.

The developmental transcriptome atlas of the biofuel crop Camelina sativa.

Camelina sativa is currently being embraced as a viable industrial bio-platform crop due to a number of desirable agronomic attributes and the unique fatty acid profile of the seed oil that has applications for food, feed and biofuel. The recent completion of the reference genome sequence of C. sativa identified a young hexaploid genome. To complement this work, we have generated a genome-wide developmental transcriptome map by RNA sequencing of 12 different tissues covering major developmental stages during the life cycle of C. sativa. We have generated a digital atlas of this comprehensive transcriptome resource that enables interactive visualization of expression data through a searchable database of electronic fluorescent pictographs (eFP browser). An analysis of this dataset supported expression of 88% of the annotated genes in C. sativa and provided a global overview of the complex architecture of temporal and spatial gene expression patterns active during development. Conventional differential gene expression analysis combined with weighted gene expression network analysis uncovered similarities as well as differences in gene expression patterns between different tissues and identified tissue-specific genes and network modules. A high-quality census of transcription factors, analysis of alternative splicing and tissue-specific genome dominance provided insight into the transcriptional dynamics and sub-genome interplay among the well-preserved triplicated repertoire of homeologous loci. The comprehensive transcriptome atlas in combination with the reference genome sequence provides a powerful resource for genomics research which can be leveraged to identify functional associations between genes and understand the regulatory networks underlying developmental processes.

Highly predictive SNP markers for efficient selection of the wheat leaf rust resistance gene Lr16.

BACKGROUND: Lr16 is a widely deployed leaf rust resistance gene in wheat (Triticum aestivum L.) that is highly effective against the North American Puccinia triticina population when pyramided with the gene Lr34. Lr16 is a seedling leaf rust resistance gene conditioning an incompatible interaction with a distinct necrotic ring surrounding the uredinium. Lr16 was previously mapped to the telomeric region of the short arm of wheat chromosome 2B. The goals of this study were to develop numerous single nucleotide polymorphism (SNP) markers for the Lr16 region and identify diagnostic gene-specific SNP marker assays for marker-assisted selection (MAS). RESULTS: Forty-three SNP markers were developed and mapped on chromosome 2BS tightly linked with the resistance gene Lr16 across four mapping populations representing a total of 1528 gametes. Kompetitive Allele Specific PCR (KASP) assays were designed for all identified SNPs. Resistance gene analogs (RGAs) linked with the Lr16 locus were identified and RGA-based SNP markers were developed. The diagnostic potential of the SNPs co-segregating with Lr16 was evaluated in a diverse set of 133 cultivars and breeding lines. Six SNP markers were consistent with the Lr16 phenotype and are accurately predictive of Lr16 for all wheat lines/cultivars in the panel. CONCLUSIONS: Lr16 was mapped relative to SNP markers in four populations. Six SNP markers exhibited high quality clustering in the KASP assay and are suitable for MAS of Lr16 in wheat breeding programs.

Polyploid evolution of the Brassicaceae during the Cenozoic era.

The Brassicaceae (Cruciferae) family, owing to its remarkable species, genetic, and physiological diversity as well as its significant economic potential, has become a model for polyploidy and evolutionary studies. Utilizing extensive transcriptome pyrosequencing of diverse taxa, we established a resolved phylogeny of a subset of crucifer species. We elucidated the frequency, age, and phylogenetic position of polyploidy and lineage separation events that have marked the evolutionary history of the Brassicaceae. Besides the well-known ancient α (47 million years ago [Mya]) and ß (124 Mya) paleopolyploidy events, several species were shown to have undergone a further more recent (∼7 to 12 Mya) round of genome multiplication. We identified eight whole-genome duplications corresponding to at least five independent neo/mesopolyploidy events. Although the Brassicaceae family evolved from other eudicots at the beginning of the Cenozoic era of the Earth (60 Mya), major diversification occurred only during the Neogene period (0 to 23 Mya). Remarkably, the widespread species divergence, major polyploidy, and lineage separation events during Brassicaceae evolution are clustered in time around epoch transitions characterized by prolonged unstable climatic conditions. The synchronized diversification of Brassicaceae species suggests that polyploid events may have conferred higher adaptability and increased tolerance toward the drastically changing global environment, thus facilitating species radiation.

Quantitative trait loci for resistance to stripe rust of wheat revealed using global field nurseries and opportunities for stacking resistance genes.

KEY MESSAGE: Quantitative trait loci controlling stripe rust resistance were identified in adapted Canadian spring wheat cultivars providing opportunity for breeders to stack loci using marker-assisted breeding. Stripe rust or yellow rust, caused by Puccinia striiformis Westend. f. sp. tritici Erikss., is a devastating disease of common wheat (Triticum aestivum L.) in many regions of the world. The objectives of this research were to identify and map quantitative trait loci (QTL) associated with stripe rust resistance in adapted Canadian spring wheat cultivars that are effective globally, and investigate opportunities for stacking resistance. Doubled haploid (DH) populations from the crosses Vesper/Lillian, Vesper/Stettler, Carberry/Vesper, Stettler/Red Fife and Carberry/AC Cadillac were phenotyped for stripe rust severity and infection response in field nurseries in Canada (Lethbridge and Swift Current), New Zealand (Lincoln), Mexico (Toluca) and Kenya (Njoro), and genotyped with SNP markers. Six QTL for stripe rust resistance in the population of Vesper/Lillian, five in Vesper/Stettler, seven in Stettler/Red Fife, four in Carberry/Vesper and nine in Carberry/AC Cadillac were identified. Lillian contributed stripe rust resistance QTL on chromosomes 4B, 5A, 6B and 7D, AC Cadillac on 2A, 2B, 3B and 5B, Carberry on 1A, 1B, 4A, 4B, 7A and 7D, Stettler on 1A, 2A, 3D, 4A, 5B and 6A, Red Fife on 2D, 3B and 4B, and Vesper on 1B, 2B and 7A. QTL on 1A, 1B, 2A, 2B, 3B, 4A, 4B, 5B, 7A and 7D were observed in multiple parents. The populations are compelling sources of recombination of many stripe rust resistance QTL for stacking disease resistance. Gene pyramiding should be possible with little chance of linkage drag of detrimental genes as the source parents were mostly adapted cultivars widely grown in Canada.

Gene-based SNP discovery in tepary bean (Phaseolus acutifolius) and common bean (P. vulgaris) for diversity analysis and comparative mapping.

BACKGROUND: Common bean (Phaseolus vulgaris) is an important grain legume and there has been a recent resurgence in interest in its relative, tepary bean (P. acutifolius), owing to this species' ability to better withstand abiotic stresses. Genomic resources are scarce for this minor crop species and a better knowledge of the genome-level relationship between these two species would facilitate improvement in both. High-throughput genotyping has facilitated large-scale single nucleotide polymorphism (SNP) identification leading to the development of molecular markers with associated sequence information that can be used to place them in the context of a full genome assembly. RESULTS: Transcript-based SNPs were identified from six common bean and two tepary bean accessions and a subset were used to generate a 768-SNP Illumina GoldenGate assay for each species. The tepary bean assay was used to assess diversity in wild and cultivated tepary bean and to generate the first gene-based map of the tepary bean genome. Genotypic analyses of the diversity panel showed a clear separation between domesticated and cultivated tepary beans, two distinct groups within the domesticated types, and P. parvifolius was confirmed to be distinct. The genetic map of tepary bean was compared to the common bean genome assembly to demonstrate high levels of collinearity between the two species with differences limited to a few intra-chromosomal rearrangements. CONCLUSIONS: The development of the first set of genomic resources specifically for tepary bean has allowed for greater insight into the structure of this species and its relationship to its agriculturally more prominent relative, common bean. These resources will be helpful in the development of efficient breeding strategies for both species and will facilitate the introgression of agriculturally important traits from one crop into the other.

A high-density SNP genotyping array for Brassica napus and its ancestral diploid species based on optimised selection of single-locus markers in the allotetraploid genome.

KEY MESSAGE: The Brassica napus Illumina array provides genome-wide markers linked to the available genome sequence, a significant tool for genetic analyses of the allotetraploid B. napus and its progenitor diploid genomes. A high-density single nucleotide polymorphism (SNP) Illumina Infinium array, containing 52,157 markers, was developed for the allotetraploid Brassica napus. A stringent selection process employing the short probe sequence for each SNP assay was used to limit the majority of the selected markers to those represented a minimum number of times across the highly replicated genome. As a result approximately 60 % of the SNP assays display genome-specificity, resolving as three clearly separated clusters (AA, AB, and BB) when tested with a diverse range of B. napus material. This genome specificity was supported by the analysis of the diploid ancestors of B. napus, whereby 26,504 and 29,720 markers were scorable in B. oleracea and B. rapa, respectively. Forty-four percent of the assayed loci on the array were genetically mapped in a single doubled-haploid B. napus population allowing alignment of their physical and genetic coordinates. Although strong conservation of the two positions was shown, at least 3 % of the loci were genetically mapped to a homoeologous position compared to their presumed physical position in the respective genome, underlying the importance of genetic corroboration of locus identity. In addition, the alignments identified multiple rearrangements between the diploid and tetraploid Brassica genomes. Although mostly attributed to genome assembly errors, some are likely evidence of rearrangements that occurred since the hybridisation of the progenitor genomes in the B. napus nucleus. Based on estimates for linkage disequilibrium decay, the array is a valuable tool for genetic fine mapping and genome-wide association studies in B. napus and its progenitor genomes.

A saturated SNP linkage map for the orange wheat blossom midge resistance gene Sm1.

KEY MESSAGE: SNP markers were developed for the OWBM resistance gene Sm1 that will be useful for MAS. The wheat Sm1 region is collinear with an inverted syntenic interval in B. distachyon. Orange wheat blossom midge (OWBM, Sitodiplosis mosellana Géhin) is an important insect pest of wheat (Triticum aestivum) in many growing regions. Sm1 is the only described OWBM resistance gene and is the foundation of managing OWBM through host genetics. Sm1 was previously mapped to wheat chromosome arm 2BS relative to simple sequence repeat (SSR) markers and the dominant, sequence characterized amplified region (SCAR) marker WM1. The objectives of this research were to saturate the Sm1 region with markers, develop improved markers for marker-assisted selection (MAS), and examine the synteny between wheat, Brachypodium distachyon, and rice (Oryza sativa) in the Sm1 region. The present study mapped Sm1 in four populations relative to single nucleotide polymorphisms (SNPs), SSRs, Diversity Array Technology (DArT) markers, single strand conformation polymorphisms (SSCPs), and the SCAR WM1. Numerous high quality SNP assays were designed that mapped near Sm1. BLAST delineated the syntenic intervals in B. distachyon and rice using gene-based SNPs as query sequences. The Sm1 region in wheat was inverted relative to B. distachyon and rice, which suggests a chromosomal rearrangement within the Triticeae lineage. Seven SNPs were tested on a collection of wheat lines known to carry Sm1 and not to carry Sm1. Sm1-flanking SNPs were identified that were useful for predicting the presence or absence of Sm1 based upon haplotype. These SNPs will be a major improvement for MAS of Sm1 in wheat breeding programs.

Genome wide SNP identification in chickpea for use in development of a high density genetic map and improvement of chickpea reference genome assembly.

BACKGROUND: In the whole genome sequencing, genetic map provides an essential framework for accurate and efficient genome assembly and validation. The main objectives of this study were to develop a high-density genetic map using RAD-Seq (Restriction-site Associated DNA Sequencing) genotyping-by-sequencing (RAD-Seq GBS) and Illumina GoldenGate assays, and to examine the alignment of the current map with the kabuli chickpea genome assembly. RESULTS: Genic single nucleotide polymorphisms (SNPs) totaling 51,632 SNPs were identified by 454 transcriptome sequencing of Cicer arietinum and Cicer reticulatum genotypes. Subsequently, an Illumina GoldenGate assay for 1,536 SNPs was developed. A total of 1,519 SNPs were successfully assayed across 92 recombinant inbred lines (RILs), of which 761 SNPs were polymorphic between the two parents. In addition, the next generation sequencing (NGS)-based GBS was applied to the same population generating 29,464 high quality SNPs. These SNPs were clustered into 626 recombination bins based on common segregation patterns. Data from the two approaches were used for the construction of a genetic map using a population derived from an intraspecific cross. The map consisted of 1,336 SNPs including 604 RAD recombination bins and 732 SNPs from Illumina GoldenGate assay. The map covered 653 cM of the chickpea genome with an average distance between adjacent markers of 0.5 cM. To date, this is the most extensive genetic map of chickpea using an intraspecific population. The alignment of the map with the CDC Frontier genome assembly revealed an overall conserved marker order; however, a few local inconsistencies within the Cicer arietinum pseudochromosome 1 (Ca1), Ca5 and Ca8 were detected. The map enabled the alignment of 215 unplaced scaffolds from the CDC Frontier draft genome assembly. The alignment also revealed varying degrees of recombination rates and hotspots across the chickpea genome. CONCLUSIONS: A high-density genetic map using RAD-Seq GBS and Illumina GoldenGate assay was developed and aligned with the existing kabuli chickpea draft genome sequence. The analysis revealed an overall conserved marker order, although some localized inversions between draft genome assembly and the genetic map were detected. The current analysis provides an insight of the recombination rates and hotspots across the chickpea genome.

High-density single nucleotide polymorphism (SNP) array mapping in Brassica oleracea: identification of QTL associated with carotenoid variation in broccoli florets.

KEY MESSAGE: A high-resolution genetic linkage map of B. oleracea was developed from a B. napus SNP array. The work will facilitate genetic and evolutionary studies in Brassicaceae. A broccoli population, VI-158 × BNC, consisting of 150 F2:3 families was used to create a saturated Brassica oleracea (diploid: CC) linkage map using a recently developed rapeseed (Brassica napus) (tetraploid: AACC) Illumina Infinium single nucleotide polymorphism (SNP) array. The map consisted of 547 non-redundant SNP markers spanning 948.1 cM across nine chromosomes with an average interval size of 1.7 cM. As the SNPs are anchored to the genomic reference sequence of the rapid cycling B. oleracea TO1000, we were able to estimate that the map provides 96 % coverage of the diploid genome. Carotenoid analysis of 2 years data identified 3 QTLs on two chromosomes that are associated with up to half of the phenotypic variation associated with the accumulation of total or individual compounds. By searching the genome sequences of the two related diploid species (B. oleracea and B. rapa), we further identified putative carotenoid candidate genes in the region of these QTLs. This is the first description of the use of a B. napus SNP array to rapidly construct high-density genetic linkage maps of one of the constituent diploid species. The unambiguous nature of these markers with regard to genomic sequences provides evidence to the nature of genes underlying the QTL, and demonstrates the value and impact this resource will have on Brassica research.

Gene-based SNP discovery and genetic mapping in pea.

KEY MESSAGE: Gene-based SNPs were identified and mapped in pea using five recombinant inbred line populations segregating for traits of agronomic importance. Pea (Pisum sativum L.) is one of the world's oldest domesticated crops and has been a model system in plant biology and genetics since the work of Gregor Mendel. Pea is the second most widely grown pulse crop in the world following common bean. The importance of pea as a food crop is growing due to its combination of moderate protein concentration, slowly digestible starch, high dietary fiber concentration, and its richness in micronutrients; however, pea has lagged behind other major crops in harnessing recent advances in molecular biology, genomics and bioinformatics, partly due to its large genome size with a large proportion of repetitive sequence, and to the relatively limited investment in research in this crop globally. The objective of this research was the development of a genome-wide transcriptome-based pea single-nucleotide polymorphism (SNP) marker platform using next-generation sequencing technology. A total of 1,536 polymorphic SNP loci selected from over 20,000 non-redundant SNPs identified using deep transcriptome sequencing of eight diverse Pisum accessions were used for genotyping in five RIL populations using an Illumina GoldenGate assay. The first high-density pea SNP map defining all seven linkage groups was generated by integrating with previously published anchor markers. Syntenic relationships of this map with the model legume Medicago truncatula and lentil (Lens culinaris Medik.) maps were established. The genic SNP map establishes a foundation for future molecular breeding efforts by enabling both the identification and tracking of introgression of genomic regions harbouring QTLs related to agronomic and seed quality traits.