BACKGROUND
brassica napus  is an allopolyploid species
<cit>  that most probably originated from inter-specific hybridisation of two brassica species that diverged approximately four million years ago
<cit> , brassica rapa  and brassica oleracea .
the genomes of the extant diploid brassica species such as brassica rapa, brassica oleracea and brassica nigra are characterised by a fundamentally triplicated structure
, indicative of a hexaploidisation event in the brassica lineage that post-dated its divergence from the arabidopsis lineage.
the triplicated genome regions remained largely collinear although they have been subjected to structural alterations; interspersed gene loss happened particularly often
.
studies in brassica species greatly benefit from genomic resources that have been assembled.
bac contig maps for the brassica rapa
<cit>  and brassica oleracea genomes were produced
<cit> .
in both studies, high-information-content fingerprinting data of bacs were exploited to establish overlaps between bacs.
integration of the resulting contigs with molecular marker maps was achieved by hybridising sequence-tagged probes to the gridded bac libraries.
for brassica oleracea, almost  <dig>  genome survey sequences  with an average length of 677 bp were established
<cit>  and a draft genome sequence is available for blast analyses
<cit> .
in addition, large collections of expressed sequence tags  and bac end sequences were put together for different brassica species
<cit> .
most importantly, a draft genome sequence was released for brassica rapa
<cit> .
a bac contig map for the brassica napus genome is not yet available, but arabidopsis thaliana genes were used as probes to screen brassica napus bac libraries by colony hybridisation in order to identify bac clones and/or contigs in regions of interest
.
due to the complex structure of the brassica napus genome it is essential that the bac clones resulting from such screens are assigned to different loci before they can be characterised in detail.
since a high-quality reference sequence for the brassica napus genome is pending, studies of specific genomic regions at sequence level relied on the analysis of bac clones and/or contigs.
sequence comparisons between corresponding a genome regions of brassica rapa ssp.
trilocularis and brassica napus var.
tapidor revealed snp frequencies varying from  <dig>  to  <dig> %.
similar values were obtained when c genome copies derived from brassica oleracea ssp.
alboglabra and brassica napus var.
tapidor were studied
<cit> .
small differences with respect to gene content and mobile elements were also detected in such comparisons
<cit> ; nevertheless it is clear that studies in brassica napus can draw on sequence resources that have been assembled for the progenitor genomes, brassica rapa and brassica oleracea.
the identification of single nucleotide polymorphisms  in brassica napus is hampered by the complex genome structure and high level of sequence identity between a and c subgenome sequences
<cit> .
due to polyploidy, two classes of polymorphisms need to be considered.
sites that are polymorphic between accessions represent the so-called intragenomic snps
<cit>  that are especially versatile for genetic mapping and discrimination of accessions.
large snp collections have been put together for brassica napus
.
snps that differentiate between the homoeologous a and c genomes were termed interhomoeolog polymorphism
<cit>  or intergenomic snp
<cit> .
this snp class proved to be valuable for the screening of bac libraries, since intergenomic snps identify and at the same time differentiate homoeologous sequences
<cit> .
compared to intragenomic snps, intergenomic snps are much more common
<cit> .
oligonucleotide extension and ligation assays form the basis for illumina’s goldengate® genotyping assay
.
for each snp locus two allele-specific oligonucleotides that carry the discriminating bases at their 3′-ends are labelled with cy <dig> and cy <dig>  respectively.
additionally, a locus-specific oligonucleotide corresponding to sequences located downstream of the snp is developed that also carries an illumicode address sequence.
the oligonucleotides are hybridised to the dna templates to be analysed and a polymerase fills the gap between one of the allele-specific oligonucleotides and the locus-specific oligonucleotide, subsequently a ligase generates a contiguous template.
the amplified dye-labelled products are annealed to beads of an optical array which are coated with sequences complementary to the different illumicode address sequences.
near pure cy <dig> or cy <dig> fluorescence will be observed for homozygous loci in diploid organisms, in contrast heterozygous loci will elicit both types of fluorescence.
assays which are specific for one of the genomes in an amphidiploid organism will produce cy3/cy <dig> fluorescence ratios indistinguishable from those observed for diploid organisms.
in contrast, assays which match target sequences on the homoeologous chromosomes may either segregate on one of the genomes or more rarely on both genomes.
in such cases cy3/cy <dig> fluorescence ratios will result that are different from those found in diploid organisms
<cit> .
it is also possible to use illumina’s goldengate® genotyping assay for screening of bac pools
<cit> .
if bac pools are analysed for the presence of a particular intergenomic snp with illumina’s goldengate® genotyping assay near pure cy <dig> or cy <dig> fluorescence will be observed for all bac pools containing only one of the homoeologous loci.
in contrast, bac pools harbouring both loci will elicit both types of fluorescence.
residual fluorescence will result if bac pools are evaluated that consist solely of bacs that do not contain any of the homoeologous loci
<cit> .
in a previous study a multidimensional screening platform was developed for a brassica napus bac library.
furthermore, the suitability of illumina’s goldengate® genotyping assay for the screening of this library was evaluated.
by applying a snp calling method especially tailored for the analysis of bac pools it was possible to identify approximately 80% of known bac coordinates from the bac library regardless whether intra- or intergenomic snps were used.
however, it was also recognised that only snp assays that discriminated between homoeologous sequences could be effectively used.
otherwise, even the use of seven or eight screening dimensions did not suffice to identify manageable lists of putative bac coordinates for the library with tenfold genome coverage.
thus, adequate sequence information is an indispensable prerequisite in order to identify snps suitable for bac library screens
<cit> .
in this study it was tested whether intergenomic snps suitable for multiplex pcr screening of brassica napus bac pools can be effectively deduced by drawing on the ever-increasing sequence resources for the progenitor genomes of brassica napus, brassica oleracea and brassica rapa
<cit> .
moreover, it was assessed whether intergenomic snps can be exploited to establish bac contig information for homoeologous regions of the brassica napus genome.
this proof-of-concept study reveals factors that need to be considered in order to apply the described methodology.
RESULTS
identification of regions suitable for development of intergenomic snp assays
a region represented by  <dig> whole genome sequence  contigs and spanning approximately one mbp on brassica rapa chromosome a <dig> and its counterpart on brassica oleracea chromosome c <dig> were chosen for this study.
parts of the brassica rapa region were also represented in six completely sequenced brassica rapa bacs , the brassica oleracea segment included a  <dig> kbp brassica oleracea ssp.
alboglabra sequence contig .figure  <dig>
brassica rapa
region on chromosome a <dig> used for design of intergenomic snps.
positions of bac sequences relative to wgs sequence contigs.
green lines mark bac sequences that were used for assay design, whereas red lines represent sequence contigs from working draft sequences of bac clones that linked the wgs sequence contigs shown as black lines.
designations  <dig> –  <dig> refer to accession numbers aeni <dig>  – aeni <dig> .
position of snp assay sequences used in this study relative to the region of interest.
intergenomic snp assays should be placed in regions in which the 20–25 bp of homoeologous sequences flanking the snp site are identical, moreover they should be distinguished by multiple mismatches and/or indels very close to the assay site when compared to homologous sequences elsewhere in the brassica napus genome
<cit> .
the six brassica rapa bacs were used as blastn
<cit>  queries to identify corresponding sequences in brassica oleracea wgs contigs.
high scoring segment pairs  spanning at least 500 bps and with sequence identities ≥95% were deemed suitable, provided that the hit sequences mapped to the target region on brassica oleracea chromosome c <dig> and did not reveal additional matches ≥150 bps and with sequence identities of ≥90% outside the region of interest.
in total,  <dig> different hsps with an average length of 827 bps and an average sequence identity of  <dig> % were identified.
for the brassica oleracea ssp.
alboglabra sequence contig corresponding brassica rapa wgs sequences were identified in the same manner.
four hsps fulfilled all criteria described above; lowering sequence identity to 92% resulted in eight hsps.
for those regions that were not covered by completely sequenced bacs a different strategy had to be adopted since brassica oleracea o <dig> wgs sequences were only temporarily available for blastn during the time of snp assay design.
sequences corresponding to annotated brassica rapa genes were extracted and used for megablast searches
<cit>  with the brassica oleracea gss dataset
<cit> .
hsps showing ≥90% sequence identity were found for approximately two-thirds of the  <dig> brassica rapa genes tested.
overlapping gss sequences were assembled such that only identical segments were retained to minimise the effect of snps between the different brassica oleracea lines used for gss sequencing and of sequencing errors in single pass sequences.
in total,  <dig> different gss contigs corresponding to  <dig> different genes were established that showed the highest blastn scores for brassica rapa chromosome a <dig> sequences.
the resulting hsps spanned on average 368 bp and showed a sequence identity of  <dig> %.
development of intergenomic snp assays
pcr amplicons were designed based on  <dig> hsps such that the homoeologous a and the c genome sequences should be amplified concomitantly in brassica napus.
sequenced amplification products of brassica napus var.
express were inspected for regions suitable for intergenomic snp assays.
for  <dig>  out of  <dig> potential assay sequences tested the intergenomic snp identified in brassica napus was also found in the relevant hsp.
however, in nine cases the hsps showed one or more additional sequence polymorphisms close to the assay site so that this particular assay would not have been selected based on the hsp sequences.
in conclusion, intergenomic snps for the homoeologous regions of the a and c genomes of brassica napus can be inferred from hsps at reasonable frequency.
only a subset of hsps was chosen for assay design.
for the entire region  <dig> assays were developed ,  <dig> were solely based on hsp sequences, for another  <dig> brassica napus amplicon sequences and hsp sequences were analysed.
screening of the brassica napusbac library with intergenomic snp assays
table
<dig> summarises the results of the bac library screening with illumina’s goldengate® genotyping assay.
the  <dig> assays identified between  <dig> and  <dig> putative bacs.
in total,  <dig> scores corresponding to  <dig> different bacs were found.
approximately a quarter of the bacs were identified with two or more assays, more than  <dig> putative bacs were revealed by four or more assays .table  <dig>
summary of bac library screen
*excluding results for assay # <dig>
assays based on intergenomic snps are expected to identify with the alternative assay nucleotides two different sets of putative bacs.
such results were observed for  <dig> assays and classification with respect to the alternative assay nucleotides was unambiguous for  <dig>  out of  <dig> bacs.
for nine assays putative bacs were found with only one of the assay nucleotides, but in case of # <dig> only a single putative bac was detected.
assay # <dig> did not reveal any putative bac .
assignment of assay nucleotides to the a and c genomes
two different brassica oleracea and brassica rapa accessions each as well as brassica napus var.
express were also analysed with illumina’s goldengate® genotyping assay in order to assign the different snp assay nucleotides experimentally to the a and c genomes.
after normalisation, signal intensities for the different samples  were plotted relative to allele frequencies .
in these plots, plants homozygous for the alternative assay nucleotides are revealing normalised theta values close to  <dig> and  <dig>  respectively.
in case of true intergenomic snps, brassica napus accessions are expected to show clusters of normalised theta values clearly distinct from the clusters of brassica rapa and brassica oleracea accessions that are in the majority of cases homozygous for the alternative assay nucleotides .
this was observed for  <dig>  out of the  <dig> assays tested .
for eight assays only one of the assay nucleotides revealed a signal with brassica napus var.
express .
snp assignment was questionable for brassica napus var.
express for five out of the  <dig> assays tested , among these were two assays, # <dig> and # <dig>  for which very low normalised r values prevented reliable snp scoring .figure  <dig>
analysis of
brassica
accessions with oligonucleotide pool assays.
normalised fluorescence intensities and normalised theta values are indicated on the x- and y-axis, respectively.
values for brassica napus var.
express, brassica rapa and brassica oleracea accessions are displayed in grey, red and blue, respectively.
for brassica rapa and brassica oleracea two accessions were analysed with four replicates each.
two independent plants with six replicates each were examined for brassica napus var.
express.
panels a, b and c show the results for assays #637- <dig>  #578- <dig> and #479-12b, the data for #386-14d, # <dig> and # <dig> are presented in panels d, e and f, respectively.
assignments of snp assay nucleotides to a and c genomes
assay
brassica rapa
brassica oleracea
brassica napus
#386-17c
#386-17a
#386-14d
#386-14c
#386-13
#386-8
#637-3
#637-5
#637-8a
#637-8b
#637-11
#11428
#11426
#11424
#11411
#11410
#11406
#11405
#11396
#11395
#11392a
#11392b
#11385-b
#11385-c
#11382a
#11382b
#11376
#11372
#11356
#11355
#11354
#479-18
#479-13d
#479-13b
#479-12b
#479-12a
#479-9
#479-3
#boli-5
#boli-3
#11308
#11307
#217-1
#217-5
#217-6a
#217-6b
#217-7
#217-9
#578-5
#578-6
assay categories are described in detail in the results section.
no <dig> and no <dig> refer to the brassica rapa and brassica oleracea accessions that were analysed by using illumina’s goldengate® genotyping assay.
for brassica napus var.
express, plant dna samples and bacs that had been allocated to the a and c subgenomes were examined.
the columns designated “snp” show the snp assay nucleotides that were found by using illumina’s goldengate® genotyping assay.
which of the different assay nucleotides matched the whole genome shotgun sequences  and/or genome survey sequences  is displayed in the columns “wgs” and “wgs ”.
data points that were not determined are marked as “nd”.
the two brassica rapa accessions were assigned to the same snp with  <dig> assays , for the two brassica oleracea accessions this was observed for  <dig> assays .
cases in which individual brassica rapa and/or brassica oleracea accessions carry different alleles and/or genes were also observed .
in case of assay #386-14d one of the brassica rapa accessions revealed residual fluorescence, thus it is either missing corresponding sequences or carries an allele that cannot be detected with this particular assay .
sequence alignments of the assay sequences to the gss sequences of brassica oleracea to1000dh <dig> or to1434
<cit>  as well as wgs sequences of brassica rapa chiifu-401
<cit>  and brassica oleracea o <dig> were also exploited to determine which of the assay nucleotides corresponded to the a and c genomes, respectively.
the analysis of  <dig> assay sequences for which wgs and gss contig sequences were available revealed in the majority of cases the same progenitor genome assignment, differences were noted in seven cases .
pcr amplicons aimed at the amplification of the homoeologous regions  were used to establish sequences for brassica napus var.
express across all assay sites.
at  <dig> snp assay sites brassica napus var.
express amplicon sequences revealed both assay nucleotides, but for #386-14d c and g were found instead of the assay nucleotides a and g. for seven snp sites only one of the assay nucleotides was observed in the amplicon sequences, in five of these cases the results were consistent with the snp genotyping of brassica napus var.
express.
in case of #578- <dig>  an intergenomic snp was observed by snp genotyping whereas the amplicon sequence was specific for the c genome copy.
indels differentiating the homoeologous sequences rendered the brassica napus var.
express amplicon sequences non-readable across the assay sites in five cases .
assignment of putative bacs to the a and c subgenomes
for each of the markers bioinformatic and experimental assignments of the assay nucleotides to the brassica rapa and brassica oleracea genomes were assessed in order to allocate tentatively the putative bacs to the a and c subgenomes.
at first, the  <dig> assays were evaluated for which both assay nucleotides had identified sets of bacs .
in case of  <dig> assays one of the assay nucleotides was consistently classified as brassica oleracea and the other as brassica rapa .
all bacs detected by these assays were assigned to the corresponding subgenomes.
for a set of thirteen assays inconsistencies were found for one of the assay nucleotides, but concordant assignments for the other assay nucleotide  permitted to ascribe the putative bacs to the appropriate subgenomes.
thus, for  <dig> out of  <dig> assays which revealed the pattern for an intergenomic snp with brassica napus var.
express, it was possible to assign the assay nucleotides to the a and c subgenomes based on the data established for the progenitor genomes.
in the remaining six cases  the classification of bacs that had been identified with markers of categories  <dig> and  <dig> was used to infer the subgenome assignments for the different assay nucleotides.
eight assays had revealed signals with only one of the assay nucleotides for dna samples of brassica napus var.
express.
such a pattern may result if one of the subgenomes carries one or several sequence polymorphisms that interfere with binding of the assay oligonucleotides.
in such a case all bacs identified will belong to the same subgenome, therefore it is possible to infer the subgenome assignment with the help of markers that allowed for unambiguous assignment .
in contrast, subgenome assignments were not made in cases in which the a and c subgenome copies of brassica napus var.
express are identical in sequence in the assayed region, since such assays identified two different sets of putative bacs that were predicted by other assays to represent the a and c subgenomes, respectively .
subgenome assignments of the assay nucleotides corresponding to assays # <dig> and # <dig> were based on the sequence comparison with wgs sequences only, since poor fluorescence intensity values prevented a meaningful evaluation of the genotyping results .
out of the putative  <dig> bacs that were identified by two or more markers belonging to assay categories  <dig> to  <dig> only  <dig>  showed inconsistencies with respect to subgenome assignment.
for markers of categories  <dig> and  <dig> a value of  <dig> % was found, thus, integrating the scores of categories  <dig> and  <dig> did not result in an increased proportion of bacs with conflicting subgenome assignment.figure
<dig> shows subgenome assignments for all  <dig> putative bacs that were identified with a set of six assays .
the four assays for which both assay nucleotides were unambiguously ascribed to the a and c subgenomes  detected in total  <dig> bacs,  <dig> of which were identified by two or more other assays belonging to categories  <dig> or  <dig>  conflicting assignments were observed for only two putative bacs, 14b <dig> and 92m <dig>  both of which were not truly mapping to the regions of interest.figure  <dig>
subgenome assignment of bacs.
all  <dig> putative bacs that were identified with the six snp assays surrounded by the box are shown as rectangles.
snp assays were ordered according to their sequence in the brassica rapa genome.
two slashes indicate that one or more of the snp assays used in this study map in the brassica rapa genome between the ones shown.
pcr experiments revealed which of the putative bacs truly map to the regions of interest.
black font marks confirmed bacs, whereas bac coordinates that were not verified to map to the homoeologous regions of interest are shown in red font.
red, dark blue and black boxes indicate snp scores that were confirmed by pcr analyses, whereas dotted outlines mark false positive snp scores.
pink or pale blue boxes correspond to scores that were only revealed by pcr but not by illumina’s goldengate® genotyping assay.
red or pink colour represent the a subgenome, c subgenome sequences are displayed in dark or pale blue.
scores of assays for which only one of the assay nucleotides revealed putative bacs are shown as black boxes.
for assays # <dig>  #11392a and # <dig> only one of the assay nucleotides each revealed signals with brassica napus var.
express and the bac pools.
bacs identified by assay #11392a had been consistently allocated to the a subgenome based on the flanking markers # <dig> and #11392b .
analysis of the c genome bacs that had been found with # <dig> and #11392b with pcr amplicons identified seven bacs that harboured sequences matching to #11392a as well as #11392b .
assay #11392a was not suitable to detect bacs carrying the c subgenome sequences because an snp and a seven-bp indel were found in the region of the c subgenome where assay oligonucleotide  <dig> should bind .
assays # <dig> and # <dig> identified bacs that had been allocated to the a and c subgenomes by closely-linked markers .
the bacs representing the a and c subgenomes were identical in sequence at the site of the assay nucleotide in both cases, but at other positions sequence differences differentiated the subgenomes .
without exception, the data derived from bac amplicon sequencing and from snp scores of flanking markers were consistent.figure  <dig>
analysis of bacs with amplicons spanning snp assays.
sequences of pcr amplicons established for representative bacs are shown below six different snp assay sequences.
horizontal arrows represent the extent of the assay nucleotides.
vertical arrows mark the positions of the snp assay nucleotides.
differences between a and c subgenome bacs are highlighted in red and blue, respectively.
bac sequences spanning # <dig>  #11392a and # <dig> are monomorphic at the site of the assay nucleotides, nonetheless a and c subgenome bacs can be distinguished based on sequence differences.
for amplicon  <dig> the region corresponding to the assay is shown and a region mapping approximately 300 bp apart.
for the latter the bac sequences are shown underneath the amplicon sequence of brassica napus var.
express.
selection of bacs for contig building
for the bac library an average insert size of  <dig> to  <dig> kbp had been estimated
<cit> .
the distance between the selected snp assays that were directly adjacent to each other in the brassica rapa genome varied from approximately 100 bp to  <dig> kbp, the median distance was  <dig> kbp.
for snp assays that were not directly adjacent to each other but which were separated by one or two other markers the median values increased to  <dig> kbp and  <dig> kbp, respectively.
the maximum distance between assays which were separated by two additional markers amounted to  <dig> kbp in the brassica rapa genome.
in order to select bacs for validation studies it was reasoned that clones that are truly present in the homoeologous regions of interest should most probably be identified by four or more assays.
the results for assay # <dig> were not considered in the selection procedure since it had revealed more than  <dig> putative bacs.
in total,  <dig> bacs were identified by four or more assays .
notably, all of the bacs shown in figure
<dig> and that are truly mapping to the a subgenome copy met this condition.
previously it was noted that not all assays are equally well suited to identify bac clones in the bac library
<cit> .
similarly, the data for assays #11392a  and #386-14d in this study had also revealed that individual loci had escaped detection in the library screen.
in order to allow for such shortcomings of the screening procedure and/or putative small rearrangements in the brassica napus a and c subgenome copies relative to the progenitor genomes putative bacs were also considered for further analysis if they were identified by any two markers in the region of interest that were separated in the brassica rapa genome by at least two other markers.
bacs 46p <dig> and 177f <dig> for example are truly mapping to one of the homoeologous regions of interest, despite the fact that they were only identified by three snp assays .
among all bacs that had not been selected based on the first criterion this condition was met by  <dig> putative bacs.
in addition, the putative bac coordinates that were identified with the two or three markers directly adjacent to each other and mapping to the very ends of the region of interest in brassica rapa were also taken into account for further studies.
in total,  <dig>  out of the  <dig> putative coordinates that had been identified in the screen were selected.
all but one of these bac clones could be recovered from the library and were subjected to validation studies.
bac contigs corresponding to homoeologous regions of interest
analysis with pcr amplicons  revealed that  <dig> and  <dig> bacs corresponding to the a and c subgenomes, respectively, were truly present in the homoeologous regions of interest.
all assay sequences which were verified to be present in a particular bac showed consistent results regarding subgenome assignment.
out of the  <dig> snp scores that had been obtained for the  <dig> confirmed bacs only  <dig> were not verified by analysis with pcr amplicons.
due to its c-genome specificity pcr amplicon a578- <dig> was not suitable to confirm twelve scores assigned to the a genome with snp assay #578- <dig> , thus only  <dig>  scores represented false positives.
the a and c subgenome bacs were assembled into two and four contigs, respectively .
the deduced order of assays in the a subgenome was generally consistent with their arrangement in the brassica rapa genome.
however,  <dig> of the bacs either suggested altered marker arrangement when compared to the brassica rapa progenitor genome and/or their assay scores were not consistent with the order of snp assays as established by other bacs.
for example, bac 27d <dig> was not identified by any of the four assays mapping between # <dig> and #11382a, although the data for several other bacs suggest a similar order as found in the brassica rapa genome .figure  <dig>
bac contigs of homoeologous regions of interest.
the order of  <dig> snp assays as found in the brassica rapa genome is indicated at the top.
assays labelled in grey font indicate that a particular assay and/or assay nucleotide did not reveal any bac that was confirmed to carry sequences corresponding to this assay.
bacs that were confirmed to carry sequences corresponding to the snp assays are displayed as rectangles.
red and dark blue boxes stand for snp scores that were validated by analyses with pcr amplicons.
pink or pale blue boxes represent false negative snp scores.
scores corresponding to the a subgenome are shown in red or pink, dark or pale blue boxes are indicative of the c subgenome.
coordinates that belong to the  <dig> bacs that had been selected initially for the validation studies are labelled in plain black font.
bac contigs are shown as filled rectangles below the names of the assays.
arrows indicate that assay order cannot be determined based on bac contig data.
coordinates labelled in bold and italics highlight additional bacs that had been integrated at the end of the bac contigs for the refinements of the contigs.
clear rectangles indicate that bac contigs were bridged by such additional bacs.
in order to refine the contig data, the same marker order as observed in brassica rapa was used to order the bacs, then the extent of all bac clones was established by analysing amplicons that spanned assays directly adjacent to the assays that had shown a score for a particular bac.
furthermore, amplicons corresponding to assays that had not yielded a signal in the bac library screen but which were predicted to be present in the bac based on the contig data were evaluated.
in total,  <dig> false negative scores were identified, about half of which were attributed to just four markers, # <dig>  # <dig>  #386-14d and #11392a.
assays # <dig> and # <dig> yielded poor fluorescence values; therefore a high rate of false negative scores was expected.
the false negative scores for #386-14d and #11392a were caused because these assays were only suitable to detect bacs representing one of the subgenomes .
in the remaining  <dig> cases of false negative scores bacs were not identified because one or several of the bac pools did not yield a signal above the applied threshold or a bac pool had been ascribed to the wrong assay nucleotide.
the determination of a bac coordinate required a score in each of seven different screening dimensions, therefore the observed rate of false negative scores has to be regarded as low.
these results are in accordance with a previous study
<cit> .elimination of false positive scores and integration of false negative snp scores resolved all conflicts with respect to marker order, only the order of assays mapping in very close proximity to each other remained ambiguous .
for selected bacs the progenitor genome assignment was also checked.
without exception, the subgenome assignments were verified.
linking of bac contigs
the criteria which were used to select putative bacs for contig building may not be optimal throughout the entire homoeologous regions of interest, due to uneven spacing of assays .
therefore, additional putative bacs which had been identified by assays mapping to the very ends of contigs were analysed in order to evaluate whether some of them may link adjacent contigs.
the gap between contigs a-1a and a-1b was bridged by bac 71d <dig> , likewise bac 54o <dig> linked contigs c-1a and c1-b .
in contrast, it was not possible to find bacs that spanned the region between contigs c- <dig> and c- <dig> in the c subgenome .
the analysis of putative bacs that had been found with assays that were flanking the gap between contigs c- <dig> and c- <dig> showed that  <dig> and  <dig> additional bacs, respectively were truly mapping to the ends of these contigs, but none of these bacs bridged the gap .
in the corresponding region of the brassica rapa genome and the brassica napus a subgenome four markers are found, #479- <dig>  #479- <dig>  #boli- <dig> and #boli- <dig> .
with #boli- <dig> three bacs had been ascribed to the c subgenome, two of which were confirmed to contain sequences corresponding to this assay .
however, the amplicon sequences established for these two bacs perfectly matched a region of brassica oleracea chromosome  <dig> which is very similar but not identical to the region of interest .
with #479- <dig> eight bacs were ascribed to the c genome.
all of them were confirmed to carry sequences corresponding to the assay, but none of these bacs had been identified with any of the other  <dig> snp assays used in this study .
none of the bacs revealed by assays #479- <dig> and #boli- <dig> had been ascribed to the c genome, moreover, the amplicon sequences for brassica napus var.
express only revealed the assay nucleotide corresponding to the a subgenome.
all these data together suggest that this particular region is largely missing from the brassica napus var.
express c subgenome and/or the bac library.
alternatively, sequence divergence may explain why neither amplicon sequences nor assays were suitable to detect bacs representing the c genome copy.
considerations for snp assay design aiming at contig building of homoeologous regions
genome assemblies solely relying on short-read sequences are especially challenging in highly heterozygous and polyploid genomes.
strategies aimed at reducing genome complexity are therefore used to alleviate genome assembly problems.
these include inbreeding, reduction of ploidy, physical separation of chromosomes and the sequencing of diploid progenitor genomes
<cit> .
sequencing of diploid progenitor genomes as surrogates for polyploid species has not only been undertaken in the brassica genus, but for example also in potato
<cit>  and gossypium
<cit> .
due to such efforts it would be possible to apply the methodology described here also in other species.
however, it will be important to consider gene density and degree of sequence identity in homoeologous segments to be analysed because these factors will determine whether hsps suitable for snp assay design can be identified in adequate number.
availability of wgs sequence contigs for both progenitor genomes allows the efficient identification of hsps suitable for the design of intergenomic snps.
using the same criteria that had been applied in this study for the regions covered by completely sequenced bacs, hsps were found for example on average every  <dig>  kbp when the brassica rapa wgs sequence contigs aeni <dig>  – aeni <dig>  were taken as queries to detect corresponding sequences in the brassica oleracea wgs data set .
moreover, the availability of wgs sequence contigs for both progenitor genomes allows elimination of all assay sequences showing highly similar matches in regions other than the target areas.
the analysis of the brassica oleracea wgs sequences at the time of assay design would for example have resulted in the omission of assay #boli- <dig> that matches two highly similar regions on brassica oleracea chromosome c <dig>  nonetheless, it is not essential to have full access to both progenitor genome sequences in order to deduce intergenomic snps for the study of an allopolyploid genome, as a brassica oleracea gss data set with  <dig> -fold estimated genome coverage
<cit>  also yielded a sufficient number of sequence alignments suitable for the design of intergenomic snps in the regions of interest.
genomic survey sequences can be generated in a particularly efficient manner by next generation sequencing technologies, as has recently been demonstrated in the brassica genus for brassica juncea
<cit> .
the establishment of a gss data set for brassica nigra  would enable the analysis of homoeologous regions in the brassica juncea  and brassica carinata  genomes with the methodology described here.
comparison of the brassica napus var.
express sequences with the corresponding sequences in the progenitor genomes revealed many instances in which intergenomic snps were conserved.
such sites are ideally suited for a bac library screen.
cases in which an intergenomic snp was detected in the progenitor genomes but not in brassica napus were also found.
although such assays did not discriminate between homoeologous bacs at least they identified the bac coordinates containing the loci of interest.
however, simultaneous detection of homoeologous sequences by a single assay nucleotide usually leads to a very high number of putative bacs
<cit> ; this was also observed for assays # <dig>  # <dig> and #ac217- <dig> in this study .
the two types of sequence polymorphisms that will result in detection failures of certain loci were rare.
in one case the intergenomic snp found in brassica napus var.
express differed from the one detected in the progenitor genomes .
sequences corresponding to two assays carried in one of the subgenomes additional snps and/or indels in close proximity to the snp site that interfered with proper functioning of the snp assay, in case of #11392a none of the bacs containing the locus of the c genome were detected , for # <dig> the identification of a subgenome bacs was not abolished but comparatively poor.
strategy for bac contig building of homoeologous regions
fingerprinting experiments
<cit>  and whole genome profiling
<cit>  are suitable strategies to build genome-wide bac contig maps.
however, these techniques alone are not suitable to build maps for specific regions of interest.
in contrast, map-based physical-mapping approaches allow the generation of clone contig maps for specific regions of interest.
using molecular markers mapping to the region of interest as hybridisation probes will identify and anchor artificial chromosome clones on the genetic map.
however, in case of polyploid genomes substantial additional efforts are needed because hybridisation probes will usually not discriminate between homoeologous sequences
.
thus, all identified clones need to be assigned to individual loci before placing them into contigs.
strategies aimed at multiplex pcr-based screening of multidimensional bac pools have to consider that for libraries encompassing several genome equivalents usually many putative clones will be detected.
consequently, substantial work has to be invested in order to identify those bacs that are truly present in the region of interest.
considerable less validation effort is necessary if sub-pools with one to twofold genome coverage are used
<cit> , however, studies aiming at contig formation usually employ libraries encompassing at least five to ten genome equivalents, otherwise too many loci will not be detected and/or it will not be possible to determine the order of markers in the contigs.
for the bac library studied here which covers approximately ten genome equivalents it had been shown that out of  <dig> loci tested one was not found in the library and for another one only a single bac was detected
<cit> .
as shown in this study the efficiency of map-based physical-mapping approaches can be fully exploited in polyploid genomes if intergenomic snps are used for the assembly of homoeologous regions, provided that the intergenomic assays are spaced such that each bac spans several assays.
then, bacs spanning three or more assays can be effectively selected from the entire set of putative bacs for validation studies.
in the course of this study  <dig> bacs had been placed into contigs in the homoeologous regions of interest,  <dig>  of which had been selected because they spanned four or more assays.
importantly, this selection procedure even worked well when results of several assays were integrated that did not differentiate between homoeologous sequences and therefore identified several hundred bacs each.
only an assay that detected more than  <dig> putative bacs had to be omitted from the analysis .
even the results of assays that match homologous sequences outside the regions of interest in addition to the target sequences can be integrated, since the bacs that correspond to such homologous sequences will usually not be identified by other assays in the regions of interest.
in contrast, only a small fraction of putative bacs that had been identified by two or three markers that were separated in the brassica rapa genome by at least two other markers were truly present in the homoeologous regions of interest.
out of the  <dig> bacs that were selected based on this criterion, only  <dig>  contained sequences mapping to the regions of interest.
instead of choosing bacs based on this criterion it was more fruitful to select bacs that were identified by two or three assays mapping to the ends of contigs.
for example, out of the  <dig> bacs that were mapping to the very ends of the regions of interest  <dig>  were confirmed to be part of the contigs.
moreover, following such a strategy bacs were identified that bridged contigs a-1a and a-1b as well as c-1a and c-1b.
in the study presented here, all snp scores were validated in order to assess the entire procedure.
taking into account that selection criteria efficiently reduced the number of putative bacs to be tested from more than  <dig> to approximately  <dig>  fingerprinting experiments
<cit>  are the method of choice to place the bacs into contigs.
in this way contig formation is validated by an independent method, moreover, the extent of bacs relative to each other is established more precisely.
this is especially important if minimal tiling paths shall be chosen for sequencing.
bac pools can be established based on the minimal tiling paths and subjected to next generation sequencing technologies as has recently been reported for large homoeologous regions of upland cotton
<cit> .
a comparison of two oryza sativa reference sequence assemblies revealed that the sequence that was established by a wgs strategy showed many more assembly errors than the sequence that was generated by a clone-by-clone sequencing approach
<cit> .
this emphasises that in all cases in which the exact order of sequence contigs on chromosomes needs to be known it may be necessary to validate assemblies.
this will be especially important in cases of highly heterozygous or polyploid genomes that show a high degree of sequence identity in their subgenomes.
in such cases validation studies may benefit from the approach described in this study.
CONCLUSIONS
clone contig maps for genomes can be efficiently generated by bac fingerprinting, but the assembly of contigs for specific chromosome regions in medium or high-throughput fashion remains a challenge, particularly in polyploid genomes in which highly similar regions need to be discriminated.
the screening of multidimensional bac pools with highly parallel snp assays has been established, however, for libraries encompassing several genome equivalents considerable additional work has to be invested to identify among the putative bacs revealed in the screen those that are truly corresponding to the assay sequences.
these limitations are overcome if a screening procedure is implemented in which intergenomic snp assays are spaced throughout homoeologous regions of interest such that any bac truly mapping in the homoeologous regions of interest will be identified by several assays.
taking into account that suitable intergenomic snps can be deduced based on wgs sequence resources of diploid progenitor genomes, the proposed approach enables the efficient assembly of bac contigs for homoeologous chromosome regions of a polyploid genome.

 brassica napus is an allopolyploid species that most probably originated from inter-specific hybridisation of two brassica species that diverged approximately four million years ago, brassica rapa and brassica oleracea. the genomes of the extant diploid brassica species such as brassica rapa, brassica oleracea and brassica nigra are characterised by a fundamentally triplicated structure, indicative of a hexaploidisation event in the brassica lineage that post-dated its divergence from the arabidopsis lineage. the triplicated genome regions remained largely collinear although they have been subjected to structural alterations; interspersed gene loss happened particularly often. studies in brassica species greatly benefit from genomic resources that have been assembled. bac contig maps for the brassica rapa and brassica oleracea genomes were produced. in both studies, high-information-content fingerprinting data of bacs were exploited to establish overlaps between bacs. integration of the resulting contigs with molecular marker maps was achieved by hybridising sequence-tagged probes to the gridded bac libraries. for brassica oleracea, almost genome survey sequences with an average length of 677 bp were established and a draft genome sequence is available for blast analyses. in addition, large collections of expressed sequence tags and bac end sequences were put together for different brassica species. most importantly, a draft genome sequence was released for brassica rapa. a bac contig map for the brassica napus genome is not yet available, but arabidopsis thaliana genes were used as probes to screen brassica napus bac libraries by colony hybridisation in order to identify bac clones and/or contigs in regions of interest. due to the complex structure of the brassica napus genome it is essential that the bac clones resulting from such screens are assigned to different loci before they can be characterised in detail. since a high-quality reference sequence for the brassica napus genome is pending, studies of specific genomic regions at sequence level relied on the analysis of bac clones and/or contigs. sequence comparisons between corresponding a genome regions of brassica rapa ssp. trilocularis and brassica napus var. tapidor revealed snp frequencies varying from to %. similar values were obtained when c genome copies derived from brassica oleracea ssp. alboglabra and brassica napus var. tapidor were studied. small differences with respect to gene content and mobile elements were also detected in such comparisons; nevertheless it is clear that studies in brassica napus can draw on sequence resources that have been assembled for the progenitor genomes, brassica rapa and brassica oleracea. the identification of single nucleotide polymorphisms in brassica napus is hampered by the complex genome structure and high level of sequence identity between a and c subgenome sequences. due to polyploidy, two classes of polymorphisms need to be considered. sites that are polymorphic between accessions represent the so-called intragenomic snps that are especially versatile for genetic mapping and discrimination of accessions. large snp collections have been put together for brassica napus. snps that differentiate between the homoeologous a and c genomes were termed interhomoeolog polymorphism or intergenomic snp. this snp class proved to be valuable for the screening of bac libraries, since intergenomic snps identify and at the same time differentiate homoeologous sequences. compared to intragenomic snps, intergenomic snps are much more common. oligonucleotide extension and ligation assays form the basis for illumina’s goldengate® genotyping assay. for each snp locus two allele-specific oligonucleotides that carry the discriminating bases at their 3′-ends are labelled with cy and cy respectively. additionally, a locus-specific oligonucleotide corresponding to sequences located downstream of the snp is developed that also carries an illumicode address sequence. the oligonucleotides are hybridised to the dna templates to be analysed and a polymerase fills the gap between one of the allele-specific oligonucleotides and the locus-specific oligonucleotide, subsequently a ligase generates a contiguous template. the amplified dye-labelled products are annealed to beads of an optical array which are coated with sequences complementary to the different illumicode address sequences. near pure cy or cy fluorescence will be observed for homozygous loci in diploid organisms, in contrast heterozygous loci will elicit both types of fluorescence. assays which are specific for one of the genomes in an amphidiploid organism will produce cy3/cy fluorescence ratios indistinguishable from those observed for diploid organisms. in contrast, assays which match target sequences on the homoeologous chromosomes may either segregate on one of the genomes or more rarely on both genomes. in such cases cy3/cy fluorescence ratios will result that are different from those found in diploid organisms. it is also possible to use illumina’s goldengate® genotyping assay for screening of bac pools. if bac pools are analysed for the presence of a particular intergenomic snp with illumina’s goldengate® genotyping assay near pure cy or cy fluorescence will be observed for all bac pools containing only one of the homoeologous loci. in contrast, bac pools harbouring both loci will elicit both types of fluorescence. residual fluorescence will result if bac pools are evaluated that consist solely of bacs that do not contain any of the homoeologous loci. in a previous study a multidimensional screening platform was developed for a brassica napus bac library. furthermore, the suitability of illumina’s goldengate® genotyping assay for the screening of this library was evaluated. by applying a snp calling method especially tailored for the analysis of bac pools it was possible to identify approximately 80% of known bac coordinates from the bac library regardless whether intra- or intergenomic snps were used. however, it was also recognised that only snp assays that discriminated between homoeologous sequences could be effectively used. otherwise, even the use of seven or eight screening dimensions did not suffice to identify manageable lists of putative bac coordinates for the library with tenfold genome coverage. thus, adequate sequence information is an indispensable prerequisite in order to identify snps suitable for bac library screens. in this study it was tested whether intergenomic snps suitable for multiplex pcr screening of brassica napus bac pools can be effectively deduced by drawing on the ever-increasing sequence resources for the progenitor genomes of brassica napus, brassica oleracea and brassica rapa. moreover, it was assessed whether intergenomic snps can be exploited to establish bac contig information for homoeologous regions of the brassica napus genome. this proof-of-concept study reveals factors that need to be considered in order to apply the described methodology.  identification of regions suitable for development of intergenomic snp assays a region represented by whole genome sequence contigs and spanning approximately one mbp on brassica rapa chromosome a and its counterpart on brassica oleracea chromosome c were chosen for this study. parts of the brassica rapa region were also represented in six completely sequenced brassica rapa bacs, the brassica oleracea segment included a kbp brassica oleracea ssp. alboglabra sequence contig.figure brassica rapa region on chromosome a used for design of intergenomic snps. positions of bac sequences relative to wgs sequence contigs. green lines mark bac sequences that were used for assay design, whereas red lines represent sequence contigs from working draft sequences of bac clones that linked the wgs sequence contigs shown as black lines. designations – refer to accession numbers aeni – aeni. position of snp assay sequences used in this study relative to the region of interest. intergenomic snp assays should be placed in regions in which the 20–25 bp of homoeologous sequences flanking the snp site are identical, moreover they should be distinguished by multiple mismatches and/or indels very close to the assay site when compared to homologous sequences elsewhere in the brassica napus genome. the six brassica rapa bacs were used as blastn queries to identify corresponding sequences in brassica oleracea wgs contigs. high scoring segment pairs spanning at least 500 bps and with sequence identities ≥95% were deemed suitable, provided that the hit sequences mapped to the target region on brassica oleracea chromosome c and did not reveal additional matches ≥150 bps and with sequence identities of ≥90% outside the region of interest. in total, different hsps with an average length of 827 bps and an average sequence identity of % were identified. for the brassica oleracea ssp. alboglabra sequence contig corresponding brassica rapa wgs sequences were identified in the same manner. four hsps fulfilled all criteria described above; lowering sequence identity to 92% resulted in eight hsps. for those regions that were not covered by completely sequenced bacs a different strategy had to be adopted since brassica oleracea o wgs sequences were only temporarily available for blastn during the time of snp assay design. sequences corresponding to annotated brassica rapa genes were extracted and used for megablast searches with the brassica oleracea gss dataset. hsps showing ≥90% sequence identity were found for approximately two-thirds of the brassica rapa genes tested. overlapping gss sequences were assembled such that only identical segments were retained to minimise the effect of snps between the different brassica oleracea lines used for gss sequencing and of sequencing errors in single pass sequences. in total, different gss contigs corresponding to different genes were established that showed the highest blastn scores for brassica rapa chromosome a sequences. the resulting hsps spanned on average 368 bp and showed a sequence identity of %. development of intergenomic snp assays pcr amplicons were designed based on hsps such that the homoeologous a and the c genome sequences should be amplified concomitantly in brassica napus. sequenced amplification products of brassica napus var. express were inspected for regions suitable for intergenomic snp assays. for out of potential assay sequences tested the intergenomic snp identified in brassica napus was also found in the relevant hsp. however, in nine cases the hsps showed one or more additional sequence polymorphisms close to the assay site so that this particular assay would not have been selected based on the hsp sequences. in conclusion, intergenomic snps for the homoeologous regions of the a and c genomes of brassica napus can be inferred from hsps at reasonable frequency. only a subset of hsps was chosen for assay design. for the entire region assays were developed, were solely based on hsp sequences, for another brassica napus amplicon sequences and hsp sequences were analysed. screening of the brassica napusbac library with intergenomic snp assays table  summarises the results of the bac library screening with illumina’s goldengate® genotyping assay. the assays identified between and putative bacs. in total, scores corresponding to different bacs were found. approximately a quarter of the bacs were identified with two or more assays, more than putative bacs were revealed by four or more assays.table summary of bac library screen *excluding results for assay # assays based on intergenomic snps are expected to identify with the alternative assay nucleotides two different sets of putative bacs. such results were observed for assays and classification with respect to the alternative assay nucleotides was unambiguous for out of bacs. for nine assays putative bacs were found with only one of the assay nucleotides, but in case of # only a single putative bac was detected. assay # did not reveal any putative bac. assignment of assay nucleotides to the a and c genomes two different brassica oleracea and brassica rapa accessions each as well as brassica napus var. express were also analysed with illumina’s goldengate® genotyping assay in order to assign the different snp assay nucleotides experimentally to the a and c genomes. after normalisation, signal intensities for the different samples were plotted relative to allele frequencies. in these plots, plants homozygous for the alternative assay nucleotides are revealing normalised theta values close to and respectively. in case of true intergenomic snps, brassica napus accessions are expected to show clusters of normalised theta values clearly distinct from the clusters of brassica rapa and brassica oleracea accessions that are in the majority of cases homozygous for the alternative assay nucleotides. this was observed for out of the assays tested. for eight assays only one of the assay nucleotides revealed a signal with brassica napus var. express. snp assignment was questionable for brassica napus var. express for five out of the assays tested, among these were two assays, # and # for which very low normalised r values prevented reliable snp scoring.figure analysis of brassica accessions with oligonucleotide pool assays. normalised fluorescence intensities and normalised theta values are indicated on the x- and y-axis, respectively. values for brassica napus var. express, brassica rapa and brassica oleracea accessions are displayed in grey, red and blue, respectively. for brassica rapa and brassica oleracea two accessions were analysed with four replicates each. two independent plants with six replicates each were examined for brassica napus var. express. panels a, b and c show the results for assays #637- #578- and #479-12b, the data for #386-14d, # and # are presented in panels d, e and f, respectively. assignments of snp assay nucleotides to a and c genomes assay brassica rapa brassica oleracea brassica napus #386-17c #386-17a #386-14d #386-14c #386-13 #386-8 #637-3 #637-5 #637-8a #637-8b #637-11 #11428 #11426 #11424 #11411 #11410 #11406 #11405 #11396 #11395 #11392a #11392b #11385-b #11385-c #11382a #11382b #11376 #11372 #11356 #11355 #11354 #479-18 #479-13d #479-13b #479-12b #479-12a #479-9 #479-3 #boli-5 #boli-3 #11308 #11307 #217-1 #217-5 #217-6a #217-6b #217-7 #217-9 #578-5 #578-6 assay categories are described in detail in the results section. no and no refer to the brassica rapa and brassica oleracea accessions that were analysed by using illumina’s goldengate® genotyping assay. for brassica napus var. express, plant dna samples and bacs that had been allocated to the a and c subgenomes were examined. the columns designated “snp” show the snp assay nucleotides that were found by using illumina’s goldengate® genotyping assay. which of the different assay nucleotides matched the whole genome shotgun sequences and/or genome survey sequences is displayed in the columns “wgs” and “wgs ”. data points that were not determined are marked as “nd”. the two brassica rapa accessions were assigned to the same snp with assays, for the two brassica oleracea accessions this was observed for assays. cases in which individual brassica rapa and/or brassica oleracea accessions carry different alleles and/or genes were also observed. in case of assay #386-14d one of the brassica rapa accessions revealed residual fluorescence, thus it is either missing corresponding sequences or carries an allele that cannot be detected with this particular assay. sequence alignments of the assay sequences to the gss sequences of brassica oleracea to1000dh or to1434 as well as wgs sequences of brassica rapa chiifu-401 and brassica oleracea o were also exploited to determine which of the assay nucleotides corresponded to the a and c genomes, respectively. the analysis of assay sequences for which wgs and gss contig sequences were available revealed in the majority of cases the same progenitor genome assignment, differences were noted in seven cases. pcr amplicons aimed at the amplification of the homoeologous regions were used to establish sequences for brassica napus var. express across all assay sites. at snp assay sites brassica napus var. express amplicon sequences revealed both assay nucleotides, but for #386-14d c and g were found instead of the assay nucleotides a and g. for seven snp sites only one of the assay nucleotides was observed in the amplicon sequences, in five of these cases the results were consistent with the snp genotyping of brassica napus var. express. in case of #578- an intergenomic snp was observed by snp genotyping whereas the amplicon sequence was specific for the c genome copy. indels differentiating the homoeologous sequences rendered the brassica napus var. express amplicon sequences non-readable across the assay sites in five cases. assignment of putative bacs to the a and c subgenomes for each of the markers bioinformatic and experimental assignments of the assay nucleotides to the brassica rapa and brassica oleracea genomes were assessed in order to allocate tentatively the putative bacs to the a and c subgenomes. at first, the assays were evaluated for which both assay nucleotides had identified sets of bacs. in case of assays one of the assay nucleotides was consistently classified as brassica oleracea and the other as brassica rapa. all bacs detected by these assays were assigned to the corresponding subgenomes. for a set of thirteen assays inconsistencies were found for one of the assay nucleotides, but concordant assignments for the other assay nucleotide permitted to ascribe the putative bacs to the appropriate subgenomes. thus, for out of assays which revealed the pattern for an intergenomic snp with brassica napus var. express, it was possible to assign the assay nucleotides to the a and c subgenomes based on the data established for the progenitor genomes. in the remaining six cases the classification of bacs that had been identified with markers of categories and was used to infer the subgenome assignments for the different assay nucleotides. eight assays had revealed signals with only one of the assay nucleotides for dna samples of brassica napus var. express. such a pattern may result if one of the subgenomes carries one or several sequence polymorphisms that interfere with binding of the assay oligonucleotides. in such a case all bacs identified will belong to the same subgenome, therefore it is possible to infer the subgenome assignment with the help of markers that allowed for unambiguous assignment. in contrast, subgenome assignments were not made in cases in which the a and c subgenome copies of brassica napus var. express are identical in sequence in the assayed region, since such assays identified two different sets of putative bacs that were predicted by other assays to represent the a and c subgenomes, respectively. subgenome assignments of the assay nucleotides corresponding to assays # and # were based on the sequence comparison with wgs sequences only, since poor fluorescence intensity values prevented a meaningful evaluation of the genotyping results. out of the putative bacs that were identified by two or more markers belonging to assay categories to only showed inconsistencies with respect to subgenome assignment. for markers of categories and a value of % was found, thus, integrating the scores of categories and did not result in an increased proportion of bacs with conflicting subgenome assignment.figure  shows subgenome assignments for all putative bacs that were identified with a set of six assays. the four assays for which both assay nucleotides were unambiguously ascribed to the a and c subgenomes detected in total bacs, of which were identified by two or more other assays belonging to categories or conflicting assignments were observed for only two putative bacs, 14b and 92m both of which were not truly mapping to the regions of interest.figure subgenome assignment of bacs. all putative bacs that were identified with the six snp assays surrounded by the box are shown as rectangles. snp assays were ordered according to their sequence in the brassica rapa genome. two slashes indicate that one or more of the snp assays used in this study map in the brassica rapa genome between the ones shown. pcr experiments revealed which of the putative bacs truly map to the regions of interest. black font marks confirmed bacs, whereas bac coordinates that were not verified to map to the homoeologous regions of interest are shown in red font. red, dark blue and black boxes indicate snp scores that were confirmed by pcr analyses, whereas dotted outlines mark false positive snp scores. pink or pale blue boxes correspond to scores that were only revealed by pcr but not by illumina’s goldengate® genotyping assay. red or pink colour represent the a subgenome, c subgenome sequences are displayed in dark or pale blue. scores of assays for which only one of the assay nucleotides revealed putative bacs are shown as black boxes. for assays # #11392a and # only one of the assay nucleotides each revealed signals with brassica napus var. express and the bac pools. bacs identified by assay #11392a had been consistently allocated to the a subgenome based on the flanking markers # and #11392b. analysis of the c genome bacs that had been found with # and #11392b with pcr amplicons identified seven bacs that harboured sequences matching to #11392a as well as #11392b. assay #11392a was not suitable to detect bacs carrying the c subgenome sequences because an snp and a seven-bp indel were found in the region of the c subgenome where assay oligonucleotide should bind. assays # and # identified bacs that had been allocated to the a and c subgenomes by closely-linked markers. the bacs representing the a and c subgenomes were identical in sequence at the site of the assay nucleotide in both cases, but at other positions sequence differences differentiated the subgenomes. without exception, the data derived from bac amplicon sequencing and from snp scores of flanking markers were consistent.figure analysis of bacs with amplicons spanning snp assays. sequences of pcr amplicons established for representative bacs are shown below six different snp assay sequences. horizontal arrows represent the extent of the assay nucleotides. vertical arrows mark the positions of the snp assay nucleotides. differences between a and c subgenome bacs are highlighted in red and blue, respectively. bac sequences spanning # #11392a and # are monomorphic at the site of the assay nucleotides, nonetheless a and c subgenome bacs can be distinguished based on sequence differences. for amplicon the region corresponding to the assay is shown and a region mapping approximately 300 bp apart. for the latter the bac sequences are shown underneath the amplicon sequence of brassica napus var. express. selection of bacs for contig building for the bac library an average insert size of to kbp had been estimated. the distance between the selected snp assays that were directly adjacent to each other in the brassica rapa genome varied from approximately 100 bp to kbp, the median distance was kbp. for snp assays that were not directly adjacent to each other but which were separated by one or two other markers the median values increased to kbp and kbp, respectively. the maximum distance between assays which were separated by two additional markers amounted to kbp in the brassica rapa genome. in order to select bacs for validation studies it was reasoned that clones that are truly present in the homoeologous regions of interest should most probably be identified by four or more assays. the results for assay # were not considered in the selection procedure since it had revealed more than putative bacs. in total, bacs were identified by four or more assays. notably, all of the bacs shown in figure  and that are truly mapping to the a subgenome copy met this condition. previously it was noted that not all assays are equally well suited to identify bac clones in the bac library. similarly, the data for assays #11392a and #386-14d in this study had also revealed that individual loci had escaped detection in the library screen. in order to allow for such shortcomings of the screening procedure and/or putative small rearrangements in the brassica napus a and c subgenome copies relative to the progenitor genomes putative bacs were also considered for further analysis if they were identified by any two markers in the region of interest that were separated in the brassica rapa genome by at least two other markers. bacs 46p and 177f for example are truly mapping to one of the homoeologous regions of interest, despite the fact that they were only identified by three snp assays. among all bacs that had not been selected based on the first criterion this condition was met by putative bacs. in addition, the putative bac coordinates that were identified with the two or three markers directly adjacent to each other and mapping to the very ends of the region of interest in brassica rapa were also taken into account for further studies. in total, out of the putative coordinates that had been identified in the screen were selected. all but one of these bac clones could be recovered from the library and were subjected to validation studies. bac contigs corresponding to homoeologous regions of interest analysis with pcr amplicons revealed that and bacs corresponding to the a and c subgenomes, respectively, were truly present in the homoeologous regions of interest. all assay sequences which were verified to be present in a particular bac showed consistent results regarding subgenome assignment. out of the snp scores that had been obtained for the confirmed bacs only were not verified by analysis with pcr amplicons. due to its c-genome specificity pcr amplicon a578- was not suitable to confirm twelve scores assigned to the a genome with snp assay #578-, thus only scores represented false positives. the a and c subgenome bacs were assembled into two and four contigs, respectively. the deduced order of assays in the a subgenome was generally consistent with their arrangement in the brassica rapa genome. however, of the bacs either suggested altered marker arrangement when compared to the brassica rapa progenitor genome and/or their assay scores were not consistent with the order of snp assays as established by other bacs. for example, bac 27d was not identified by any of the four assays mapping between # and #11382a, although the data for several other bacs suggest a similar order as found in the brassica rapa genome.figure bac contigs of homoeologous regions of interest. the order of snp assays as found in the brassica rapa genome is indicated at the top. assays labelled in grey font indicate that a particular assay and/or assay nucleotide did not reveal any bac that was confirmed to carry sequences corresponding to this assay. bacs that were confirmed to carry sequences corresponding to the snp assays are displayed as rectangles. red and dark blue boxes stand for snp scores that were validated by analyses with pcr amplicons. pink or pale blue boxes represent false negative snp scores. scores corresponding to the a subgenome are shown in red or pink, dark or pale blue boxes are indicative of the c subgenome. coordinates that belong to the bacs that had been selected initially for the validation studies are labelled in plain black font. bac contigs are shown as filled rectangles below the names of the assays. arrows indicate that assay order cannot be determined based on bac contig data. coordinates labelled in bold and italics highlight additional bacs that had been integrated at the end of the bac contigs for the refinements of the contigs. clear rectangles indicate that bac contigs were bridged by such additional bacs. in order to refine the contig data, the same marker order as observed in brassica rapa was used to order the bacs, then the extent of all bac clones was established by analysing amplicons that spanned assays directly adjacent to the assays that had shown a score for a particular bac. furthermore, amplicons corresponding to assays that had not yielded a signal in the bac library screen but which were predicted to be present in the bac based on the contig data were evaluated. in total, false negative scores were identified, about half of which were attributed to just four markers, # # #386-14d and #11392a. assays # and # yielded poor fluorescence values; therefore a high rate of false negative scores was expected. the false negative scores for #386-14d and #11392a were caused because these assays were only suitable to detect bacs representing one of the subgenomes. in the remaining cases of false negative scores bacs were not identified because one or several of the bac pools did not yield a signal above the applied threshold or a bac pool had been ascribed to the wrong assay nucleotide. the determination of a bac coordinate required a score in each of seven different screening dimensions, therefore the observed rate of false negative scores has to be regarded as low. these results are in accordance with a previous study.elimination of false positive scores and integration of false negative snp scores resolved all conflicts with respect to marker order, only the order of assays mapping in very close proximity to each other remained ambiguous. for selected bacs the progenitor genome assignment was also checked. without exception, the subgenome assignments were verified. linking of bac contigs the criteria which were used to select putative bacs for contig building may not be optimal throughout the entire homoeologous regions of interest, due to uneven spacing of assays. therefore, additional putative bacs which had been identified by assays mapping to the very ends of contigs were analysed in order to evaluate whether some of them may link adjacent contigs. the gap between contigs a-1a and a-1b was bridged by bac 71d, likewise bac 54o linked contigs c-1a and c1-b. in contrast, it was not possible to find bacs that spanned the region between contigs c- and c- in the c subgenome. the analysis of putative bacs that had been found with assays that were flanking the gap between contigs c- and c- showed that and additional bacs, respectively were truly mapping to the ends of these contigs, but none of these bacs bridged the gap. in the corresponding region of the brassica rapa genome and the brassica napus a subgenome four markers are found, #479- #479- #boli- and #boli-. with #boli- three bacs had been ascribed to the c subgenome, two of which were confirmed to contain sequences corresponding to this assay. however, the amplicon sequences established for these two bacs perfectly matched a region of brassica oleracea chromosome which is very similar but not identical to the region of interest. with #479- eight bacs were ascribed to the c genome. all of them were confirmed to carry sequences corresponding to the assay, but none of these bacs had been identified with any of the other snp assays used in this study. none of the bacs revealed by assays #479- and #boli- had been ascribed to the c genome, moreover, the amplicon sequences for brassica napus var. express only revealed the assay nucleotide corresponding to the a subgenome. all these data together suggest that this particular region is largely missing from the brassica napus var. express c subgenome and/or the bac library. alternatively, sequence divergence may explain why neither amplicon sequences nor assays were suitable to detect bacs representing the c genome copy. considerations for snp assay design aiming at contig building of homoeologous regions genome assemblies solely relying on short-read sequences are especially challenging in highly heterozygous and polyploid genomes. strategies aimed at reducing genome complexity are therefore used to alleviate genome assembly problems. these include inbreeding, reduction of ploidy, physical separation of chromosomes and the sequencing of diploid progenitor genomes. sequencing of diploid progenitor genomes as surrogates for polyploid species has not only been undertaken in the brassica genus, but for example also in potato and gossypium. due to such efforts it would be possible to apply the methodology described here also in other species. however, it will be important to consider gene density and degree of sequence identity in homoeologous segments to be analysed because these factors will determine whether hsps suitable for snp assay design can be identified in adequate number. availability of wgs sequence contigs for both progenitor genomes allows the efficient identification of hsps suitable for the design of intergenomic snps. using the same criteria that had been applied in this study for the regions covered by completely sequenced bacs, hsps were found for example on average every kbp when the brassica rapa wgs sequence contigs aeni – aeni were taken as queries to detect corresponding sequences in the brassica oleracea wgs data set. moreover, the availability of wgs sequence contigs for both progenitor genomes allows elimination of all assay sequences showing highly similar matches in regions other than the target areas. the analysis of the brassica oleracea wgs sequences at the time of assay design would for example have resulted in the omission of assay #boli- that matches two highly similar regions on brassica oleracea chromosome c nonetheless, it is not essential to have full access to both progenitor genome sequences in order to deduce intergenomic snps for the study of an allopolyploid genome, as a brassica oleracea gss data set with -fold estimated genome coverage also yielded a sufficient number of sequence alignments suitable for the design of intergenomic snps in the regions of interest. genomic survey sequences can be generated in a particularly efficient manner by next generation sequencing technologies, as has recently been demonstrated in the brassica genus for brassica juncea. the establishment of a gss data set for brassica nigra would enable the analysis of homoeologous regions in the brassica juncea and brassica carinata genomes with the methodology described here. comparison of the brassica napus var. express sequences with the corresponding sequences in the progenitor genomes revealed many instances in which intergenomic snps were conserved. such sites are ideally suited for a bac library screen. cases in which an intergenomic snp was detected in the progenitor genomes but not in brassica napus were also found. although such assays did not discriminate between homoeologous bacs at least they identified the bac coordinates containing the loci of interest. however, simultaneous detection of homoeologous sequences by a single assay nucleotide usually leads to a very high number of putative bacs; this was also observed for assays # # and #ac217- in this study. the two types of sequence polymorphisms that will result in detection failures of certain loci were rare. in one case the intergenomic snp found in brassica napus var. express differed from the one detected in the progenitor genomes. sequences corresponding to two assays carried in one of the subgenomes additional snps and/or indels in close proximity to the snp site that interfered with proper functioning of the snp assay, in case of #11392a none of the bacs containing the locus of the c genome were detected, for # the identification of a subgenome bacs was not abolished but comparatively poor. strategy for bac contig building of homoeologous regions fingerprinting experiments and whole genome profiling are suitable strategies to build genome-wide bac contig maps. however, these techniques alone are not suitable to build maps for specific regions of interest. in contrast, map-based physical-mapping approaches allow the generation of clone contig maps for specific regions of interest. using molecular markers mapping to the region of interest as hybridisation probes will identify and anchor artificial chromosome clones on the genetic map. however, in case of polyploid genomes substantial additional efforts are needed because hybridisation probes will usually not discriminate between homoeologous sequences. thus, all identified clones need to be assigned to individual loci before placing them into contigs. strategies aimed at multiplex pcr-based screening of multidimensional bac pools have to consider that for libraries encompassing several genome equivalents usually many putative clones will be detected. consequently, substantial work has to be invested in order to identify those bacs that are truly present in the region of interest. considerable less validation effort is necessary if sub-pools with one to twofold genome coverage are used, however, studies aiming at contig formation usually employ libraries encompassing at least five to ten genome equivalents, otherwise too many loci will not be detected and/or it will not be possible to determine the order of markers in the contigs. for the bac library studied here which covers approximately ten genome equivalents it had been shown that out of loci tested one was not found in the library and for another one only a single bac was detected. as shown in this study the efficiency of map-based physical-mapping approaches can be fully exploited in polyploid genomes if intergenomic snps are used for the assembly of homoeologous regions, provided that the intergenomic assays are spaced such that each bac spans several assays. then, bacs spanning three or more assays can be effectively selected from the entire set of putative bacs for validation studies. in the course of this study bacs had been placed into contigs in the homoeologous regions of interest, of which had been selected because they spanned four or more assays. importantly, this selection procedure even worked well when results of several assays were integrated that did not differentiate between homoeologous sequences and therefore identified several hundred bacs each. only an assay that detected more than putative bacs had to be omitted from the analysis. even the results of assays that match homologous sequences outside the regions of interest in addition to the target sequences can be integrated, since the bacs that correspond to such homologous sequences will usually not be identified by other assays in the regions of interest. in contrast, only a small fraction of putative bacs that had been identified by two or three markers that were separated in the brassica rapa genome by at least two other markers were truly present in the homoeologous regions of interest. out of the bacs that were selected based on this criterion, only contained sequences mapping to the regions of interest. instead of choosing bacs based on this criterion it was more fruitful to select bacs that were identified by two or three assays mapping to the ends of contigs. for example, out of the bacs that were mapping to the very ends of the regions of interest were confirmed to be part of the contigs. moreover, following such a strategy bacs were identified that bridged contigs a-1a and a-1b as well as c-1a and c-1b. in the study presented here, all snp scores were validated in order to assess the entire procedure. taking into account that selection criteria efficiently reduced the number of putative bacs to be tested from more than to approximately fingerprinting experiments are the method of choice to place the bacs into contigs. in this way contig formation is validated by an independent method, moreover, the extent of bacs relative to each other is established more precisely. this is especially important if minimal tiling paths shall be chosen for sequencing. bac pools can be established based on the minimal tiling paths and subjected to next generation sequencing technologies as has recently been reported for large homoeologous regions of upland cotton. a comparison of two oryza sativa reference sequence assemblies revealed that the sequence that was established by a wgs strategy showed many more assembly errors than the sequence that was generated by a clone-by-clone sequencing approach. this emphasises that in all cases in which the exact order of sequence contigs on chromosomes needs to be known it may be necessary to validate assemblies. this will be especially important in cases of highly heterozygous or polyploid genomes that show a high degree of sequence identity in their subgenomes. in such cases validation studies may benefit from the approach described in this study.  clone contig maps for genomes can be efficiently generated by bac fingerprinting, but the assembly of contigs for specific chromosome regions in medium or high-throughput fashion remains a challenge, particularly in polyploid genomes in which highly similar regions need to be discriminated. the screening of multidimensional bac pools with highly parallel snp assays has been established, however, for libraries encompassing several genome equivalents considerable additional work has to be invested to identify among the putative bacs revealed in the screen those that are truly corresponding to the assay sequences. these limitations are overcome if a screening procedure is implemented in which intergenomic snp assays are spaced throughout homoeologous regions of interest such that any bac truly mapping in the homoeologous regions of interest will be identified by several assays. taking into account that suitable intergenomic snps can be deduced based on wgs sequence resources of diploid progenitor genomes, the proposed approach enables the efficient assembly of bac contigs for homoeologous chromosome regions of a polyploid genome. 
this proof-of-concept study shows that sequence resources of diploid progenitor genomes can be used to deduce intergenomic snps suitable for multiplex polymerase chain reaction -based screening of multidimensional bac pools of a polyploid organism. based on the results of snp assays, a bac contig for the brassica napus a subgenome was established that spanned the entire region of interest. sequence alignments between brassica rapa and brassica oleracea sequences mapping to corresponding regions on chromosomes a and c respectively were used to identify single nucleotide polymorphisms between the a and c genomes. intergenomic snps mapping throughout homoeologous chromosome segments spanning approximately one mbp each were included in illumina’s goldengate® genotyping assay and used to screen multidimensional pools of a brassica napus bacterial artificial chromosome library with tenfold genome coverage. owing to their high abundance and ease of identification, intergenomic snps represent a versatile tool to establish bac contigs for homoeologous regions of a polyploid genome. the c subgenome region was represented in three bac contigs. homoeologous sequences pose a particular challenge if bacterial artificial chromosome contigs shall be established for specific regions of an allopolyploid genome. bac contigbac librarybrassica napushomoeologous regionsmultidimensional poolspcr screeningpolyploidysnpissue-copyright-statement© the author 2014 a large fraction of these polymorphisms was also present in brassica napus, an allopolyploid species that originated from hybridisation of a and c genome species. single nucleotide polymorphisms differentiating between homoeologous genomes may represent a suitable screening tool for such purposes, since they do not only identify homoeologous sequences but also differentiate between them.
 homoeologous sequences pose a particular challenge if bacterial artificial chromosome contigs shall be established for specific regions of an allopolyploid genome. single nucleotide polymorphisms differentiating between homoeologous genomes may represent a suitable screening tool for such purposes, since they do not only identify homoeologous sequences but also differentiate between them.  sequence alignments between brassica rapa and brassica oleracea sequences mapping to corresponding regions on chromosomes a and c respectively were used to identify single nucleotide polymorphisms between the a and c genomes. a large fraction of these polymorphisms was also present in brassica napus, an allopolyploid species that originated from hybridisation of a and c genome species. intergenomic snps mapping throughout homoeologous chromosome segments spanning approximately one mbp each were included in illumina’s goldengate® genotyping assay and used to screen multidimensional pools of a brassica napus bacterial artificial chromosome library with tenfold genome coverage. based on the results of snp assays, a bac contig for the brassica napus a subgenome was established that spanned the entire region of interest. the c subgenome region was represented in three bac contigs.  this proof-of-concept study shows that sequence resources of diploid progenitor genomes can be used to deduce intergenomic snps suitable for multiplex polymerase chain reaction -based screening of multidimensional bac pools of a polyploid organism. owing to their high abundance and ease of identification, intergenomic snps represent a versatile tool to establish bac contigs for homoeologous regions of a polyploid genome. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords bac contigbac librarybrassica napushomoeologous regionsmultidimensional poolspcr screeningpolyploidysnpissue-copyright-statement© the author 2014 