dna - based genetic markers are pivotal tools for applications as diverse as qtl mapping , marker assisted selection , association mapping , and fine mapping for the detection of genes linked to a particular phenotype . among the variety of genetic markers that have been developed ,
those derived from single nucleotide polymorphisms ( snps ) have become the marker of choice for many mapping applications because of their abundance and the availability of high - throughput and cost - effective technologies for detection and diagnostics [ 24 ] .
one popular tool for snp identification and detection has been the construction of reduced - representation libraries ( rrl ) and their sequencing with massively parallel sequencing platforms , in species as varied as cattle , worm , soybean , rice , maize , or common bean [ 510 ] . however , one major limitation of such platforms is the relatively short length of individual sequencing reads .
while the availability of a high quality reference sequence may render short reads sufficient for alignment and subsequent snp detection , this limitation may be further compounded in crop species due to ( 1 ) the inherent complexity of genomes ( and transcriptomes ) in economically important species , such as maize , soybean , or canola , due to an elevation in ploidy and/or the frequent expansion of paralogous sequences and gene families , and the need to generate very long sequencing reads for resolving highly duplicated sequences within a single genome , and ( 2 ) the potentially large number of polymorphisms between lines in those same plant species ( including indels and in regions flanking a polymorphism of interest ) and the need to provide long line - specific sequences for identifying variants , open reading frames , or other biologically active regions for lines whose genome sequence has a significantly altered composition in comparison to the reference assembly . because of those limitations , the roche 454 flx platform   is often used as the instrument of choice for providing long sequencing reads and generating an appropriate sequencing scaffold .
however , its relatively lower sequencing throughput when compared to other second - generation sequencing instruments often means that a large number of runs are needed to generate a high quality sequencing assembly from very large genomes .
more recently , the pacbio rs platform from pacific biosciences   has garnered a lot of attention due to its ability to generate very long sequencing reads , in the kilobase range , but its low raw accuracy rate makes it , for now , unpractical to use as a variant detection platform . in a previous study
, we developed a methodology for rapid snp detection in rice and soybean that can be applied to a wide range of moderately or highly complex plant genomes where sufficient genomic reference sequences are available .
the methodology , based on digesting the genome with a methylation - sensitive restriction endonuclease ( re ) followed by a secondary digestion with the 4-bp restriction endonuclease dpnii , generates short dna fragments that are sequenced at one end using the illumina genome analyzer ii .
thus , they can be directly compared and aligned to a reference assembly for snp detection . however , one problem inherent to this approach is the lack of sequencing information flanking the 32  bp read of interest and the need to rely on reference genomic sequences from a different line or species for primer design and snp assay development .
the relative short length of the sequence also can lead to relatively high false positive discovery rates in species lacking robust reference sequences .
several studies already have shown that dna fragments sequenced at both ends on a massively parallel sequencing platforms can be locally assembled into contigs hundreds or even thousands of bases long .
hiatt et al .   developed a methodology to generate long consensus sequences in pseudomonas aeruginosa from short second - generation sequencing reads .
in there , ~500  bps dna fragments are sheared randomly to produce fragments shorter than the original fragment set .
those shorter fragments then are end - sequenced at both ends with the illumina genome analyzer ii and the resulting read pairs are grouped together and assembled to recreate the sequence of the original ~500  bp dna fragment .
more recently , etter et al .   developed a paired - end restriction - site associated dna ( rad - pe ) method , where ~350850  bps rad fragments from three - spined sticklebacks ( gasterosteus aculeatus ) , with one end anchored by a restriction endonuclease site and the other end randomly sheared , were sequenced on the illumina genome analyzer ii , to locally assemble staggered randomly sheared ends associated to a particular restriction site into restriction site - specific contigs several hundred bases in length .
the rad - pe protocol requires several steps , including the digestion of genomic dna , ligation of two independent adapters , and the selective pcr amplification of rad fragments flanked only by the two distinct adapters . a variation in protocol , namely , partial digestion of the genomic dna rather than full digestion , creates overlapping rad fragments , which in turn can be assembled into overlapping contigs covering sequences several thousand bases in length .
on two guppy ( poecilia reticulata ) populations , generating consensus ~200400  bps sequences that were used for polymorphism discovery , checking for the presence or absence of restriction - site associated sequences in the two populations to identify them as polymorphic ( as the absence of a sequence suggests the presence of a polymorphism in the restriction site ) . in this study , a revised version of the rad - pe method was tested and implemented in maize , in order to assess the feasibility of using such a strategy for snp detection on the large and complex genome of an economically important crop .
the revised rad - pe protocol shown here extends the concept presented in deschamps et al . , where snp detection in complex crop genomes had been successfully attempted .
briefly , the digestion of dna fragments with dpnii is replaced by random shearing , and the sheared dna fragments containing the methylation - sensitive re site at one end are recovered , size - selected , enriched , and sequenced at both ends on the illumina genome analyzer iix .
two steps are used to select for dna fragments prior to paired - end sequencing , namely , a biotin selection step performed after the initial digestion with a methylation - sensitive re , and a second selection step when dna fragments are selectively amplified via pcr prior to sequencing .
the resulting individual read pairs then are assembled de novo to create larger contigs that can be used both for snp detection and snp assay design for marker development .
these contigs can be several hundreds bps long and their expected length can be modulated simply by selecting appropriate dna fragment size ranges on gel .
in addition to snp discovery and assay design , they can be used in multiple species for applications whose purpose is facilitated by long sequencing reads , such as copy number variant detection , open reading frame discovery , or metagenomics .
  maize genomic dna samples were extracted according to the protocol described in deschamps et al . .  
b73 and b14 genomic dna ( 10  g ) was digested for 4  h at 37c with 15  u of psti ( promega ) in a total volume of 20  l of 1x psti buffer ( promega ) containing 1x acetylated bsa ( promega ) .
the enzyme was inactivated at 65c for 20  min then the digested genomic dna was purified with qiaquick pcr purification spin columns ( qiagen ) .
biotinylated psti - specific adapters were created by mixing 1500  pmol each of two synthetic oligonucleotides ( upper strand , /5bio
/ ggttgacatgctggattgagacctgcaggtgc*a , where * is a phosphorothioate bond ; lower strand , /5po4-/cctgcaggtctcaatccagcatgtc ) in 100  l water , heating them up at 95c for 2  min , and then allowing them to cool slowly to room temperature . a ~50-fold excess of psti - specific adapters ( 37.5  pmol ) relative to available psti ends was ligated to the digested dna , in a total volume of 20  l of 1x ligase buffer ( new england biolabs ) containing 10,000  u of t4 dna ligase ( new england biolabs ) .
the reaction was incubated overnight at 16c then 10  min at 70c , and then purified with qiaquick pcr purification spin columns ( qiagen ) .
the digested dna was randomly sheared via nebulization at 32  psi for 6  min then purified with qiaquick pcr purification spin columns ( qiagen ) .   
after random shearing , a biotin selection step allowed the capture of smaller dna fragments containing a psti site at one end . for biotin selection , 100  l of streptavidin - dynabeads m-280 ( invitrogen ) were washed twice with 1  ml te buffer then resuspended in 100  l 2x b&w buffer ( 10  mmol  l tris - hcl , ph 7.5 ; 1  mmol  l edta ; 2  mol  l nacl ) .
the dna fragments then were added to the beads and incubated for 30  min at 30c with gentle horizontal mixing . after withdrawing the supernatant ,
the beads were collected and washed three times with 1  ml 1x b&w buffer then three times with 1  ml te buffer .
the washed beads were resuspended in 200  l of 1x nebuffer 4 ( new england biolabs ) containing 20  u sbfi ( new england biolabs ) and incubated for 90  min at 37c with gentle horizontal mixing . an sbfi recognition site ( 5-cctgcagg-3 ) located
immediately upstream from the psti site on the biotinylated adapter allowed for the cleavage and recovery of the biotin - bound dna fragments .
the supernatant containing the released dna fragments was then transferred to a new tube , extracted with 100  l phenol / chloroform / isoamyl alcohol ( 25  :  24  :  1 ) and precipitated with ethanol and 3  mol  l naoac .
end - repair ,  a  base addition and ligation of adapters were performed according to the protocol developed by illumina for preparing dna samples for paired - end sequencing .
buffer and enzymatic reagents were obtained directly from illumina 's pe sample prep kit ( catalog number pe-102 - 1002 ) .
briefly , dna was end repaired for 30  min at 20c in a total volume of 100  l of 1x t4 dna ligase buffer with 10  mm atp containing 4  l 10  mm dntp mix , 5  l t4 dna polymerase , 1  l klenow enzyme , and t4 pnk .
after incubation , dna was purified with a qiaquick pcr purification kit ( qiagen ) and resuspended in 32  l eb buffer ( qiagen ) . 
a  bases then were added to the 3 ends of blunt - ended dna fragments .
dna was incubated at 37c for 30  min in 1x klenow buffer containing 10  l 1  mm datp and 3  l klenow exo ( 3 to 5 exo minus ) .
after incubation , dna was purified with a qiaquick minelute purification kit ( qiagen ) and resuspended in 10  l eb buffer ( qiagen ) .
illumina pe adapters then were ligated to dna fragments at 20c for 15  min in a total volume of 50  l of 1x dna ligase buffer containing 1  l pe adapter oligo mix and 5  l dna ligase . after incubation
, dna was purified with a qiaquick pcr purification kit ( qiagen ) and resuspended in 30  l eb buffer ( qiagen ) . ligated dna then were loaded on an e - gel 2% with sybr safe agarose gel ( invitrogen ) in the presence of low molecular weight dna ladder ( new england biolabs ) for size selection . after running the gel for 26  min in an e - gel powerbase v4 electrophoresis unit ( invitrogen ) ,
a 200500  bp dna smear was cut out and purified with the qiaquick gel extraction kit ( qiagen ) .
after ligation of illumina adapters and gel - based size selection of the 200500  bps ligated dna fragments , a pcr amplification step was performed using the illumina paired - end ( pe ) pcr primer 1.0 and a modified pe pcr primer 2.0 ending in 5-tgcaggtgca-3 ( matching adjacent sbfi and psti recognition sites ) . the use of illumina pe adapters and of a modified pe pcr primer 2.0 further allowed the selective amplification of dna fragments containing a psti recognition site at one end . for pcr amplification ,
1  l of size - selected dna was incubated in 50  l of 1x phusion hf master mix ( finnzymes ) in the presence of pcr primer pe 1.0 ( illumina ) and the modified pcr primer pe 2.0 ( caagcagaagacggcatacgagatcggtctcggcattcctgctgaaccgctcttccgatcttgcaggtgc*a , where * is a phosphorothioate bond and  tgcaggtgc*a  is a signature 3 sequence containing leftovers from sbfi & psti restriction sites ) . after
pcr amplification ( 30  s at 98c , followed by 20 rounds of 10  s at 98c , 30  s at 65c , 30  s at 72c , and a final extension step of 5  min at 72c ) , dna was purified with a qiaquick pcr purification kit ( qiagen ) then subjected to a second round of gel - based size selection using conditions similar to the ones described above .
the 200500  bps amplified dna fragments were size - selected on gel and the presence of illumina adapters at both ends suggests that the actual sizes of the fragments that were sequenced were between ~70  bps and ~370  bps .   cluster generation and paired - end sequencing were performed on an illumina cluster station and genome analyzer iix , respectively , according to protocols and recipes developed by illumina .
sequencing was performed for 76 cycles at both ends of the clustered dna fragments using pe sequencing primers for read 1 and read 2 ( illumina ) , generating read pairs ,  reads 1  and  reads 2  , mapping the ends of pcr products generated with pe pcr primer 1.0 and the modified pe pcr primer 2.0 , respectively .
the resulting read 1 and read 2 sequences were grouped into  read pairs  according to the x and y coordinates of the corresponding dna cluster on the flow cell .
sequencing reads and quality scores were generated in a real - time fashion with the illumina data collection software v2.6 .
after initial base calling , additional custom filtering was performed using calibrated quality scores generated by the illumina pipeline .
reads generated from both ends of dna fragments ( reads 1 and reads 2 ) were trimmed by removing from the 3 ends bases with a phred - equivalent quality score below 10 .
a length threshold of 24 was applied to filtering , indicating that all bases < 24 bases in length after trimming were removed from further analysis .
the use of a modified pe pcr primer 2.0 led to a vast majority of read 2 starting with the signature sequence 5-tgcaggtgca-3. however , in spite of a significant bias in base composition for the first 10 bases of reads 2 , a majority of all read 2 data were high quality reads that paired well with their read 1 counterparts .
nonetheless , it must be noted that since illumina uses images of the first few cycles to calculate phasing , such low 5 sequence variation could cause base calling errors in subsequent sequencing cycles if a different illumina instrument or upgraded base calling software were to be used for sequencing . in that particular case , increasing the amount of phix control dna spiked into each lane ( from ~1% to ~25% ) would be recommended .
the replacement of sbfi ( and of its corresponding recognition sequences in the biotinylated psti - specific adapter ) by a type iis enzyme ( such as bseri ) , thus removing the 10  bp  tgcaggtgca  signature sequence from the 5 ends of reads 2 , also is being explored on separate sets of genomic dna samples .
the 10  bp  tgcaggtgca  signature sequence at the 5 end of read 2 sequences was identified by calculating the hamming distance with the first 10 bases on read 2 . only read 2 sequences that have a hamming distance less than or equal to 2 were considered to have a pst1 site .
local assemblies for each distinct read 2 sequences and their corresponding read 1 sequences were generated using the velvet software   with k - mer values ranging from 23 to 63 with increments of 4 and coverage cutoff values of 4 and 8 for each k - mer value .
of the 22 assemblies generated with varying k - mer and coverage cutoff values , the assembly with the largest contig size was considered the final one for that distinct read 2 sequence and its paired read 1 sequences .  
sequences from 100 random b14 and b73 contigs with at least 1 polymorphism to the b73 reference genome v1.0 assembly were used to design pcr primers ( length 1824  bps , tm 6065c ) using a local version of the primer3 primer design software tool .
crossmatch analysis was performed on contig sequences using a local maize repeat database to mask repetitive dna sequences prior to primer design and m13 forward and reverse  tails  were added to the 5 ends of the pcr primer sequences before ordering .
b73 and b14 genomic dna was subjected to pcr amplification ( 15  min at 95c , followed by 40 rounds of 30  s at 95c , 30  s at 60c , 1  min at 72c , and a final extension step of 10  min at 72c ) by using 5  pmol each of the  tailed  pcr primers in a total volume of 10  l of 1x hotstartaq master mix ( qiagen ) containing 25  u hotstartaq dna polymerase ( qiagen ) . 1.5  mmol  l mgcl2 and 200  mol  l dntps .
pcr cleanup reactions were performed by mixing 2  l of pcr products with 0.75  l of exosap - it ( usb corporation ) in a total volume of 17  l with sterile distilled water , and incubating at 37c for 25  min then 80c for 25  min .
5  l of the cleaned - up amplified dna then were end - sequenced using m13 forward and reverse oligonucleotides and the abi bigdye version 3.1 prism sequencing kit .
after ethanol - based cleanup , cycle sequencing reaction products were resolved and detected on life technologies ( carlsbad , ca , usa ) abi3730xl automated capillary dna sequences .
individual sequences from each genotype were combined into a single project ( one project per amplified fragment ) and assembled with the phred / phrap / consed package ( see http://www.phrap.org/phredphrapconsed.html ) .
confirmed snps identified on the genome analyzer were validated by comparison to single - base mismatches between the b73 and b14 genotypes located in regions that matched the original contig sequence generated by assembling individual read pair sequences with velvet .
library construction and massively parallel sequencing were performed on two public maize inbred lines , b73 and b14 , according to the method described in figure 1 .
dna samples first were digested with the methylation - sensitive re psti ( 5-c*tgcag-3 ) .
the psti activity is blocked by 5-c methylation of the cytosine in the first position ( c * ) in the recognition sequence .
the choice of psti was guided by ( 1 ) the intention of enriching for genic regions and avoiding the capture of the repeated fraction of the genome [ 18 , 19 ] , and ( 2 ) the potential number of unique sites digested by this enzyme in both samples .   
after trimming and filtering , a total of 63.9 million and 94.9 million high - quality read pairs ( reads 1 and reads 2 ) were obtained for b73 and b14 , respectively , from one run on the illumina genome analyzer ii ( table 1 ) .
the disparity between those numbers is due in part to the fact that 3 lanes and 4 lanes of a flow cell were used for b73 and b14 , respectively .
a total of 61.5 million ( 96.2% of total ) and 92.2 million ( 97.1% of total ) reads 2 for b73 and b14 , respectively , contained the signature 5-tgcaggtgca-3 sequence at their 5 ends .   
since reads 2 are anchored by a methylation - sensitive re cut site , each read pairs can be grouped based on the identity of their respective read 2 sequences . in order
to assign read pairs to specific regions of the genome , and perform local assemblies of such regions , the following analysis was performed : ( 1 ) assessing the number of read 2 sequences aligning uniquely to the b73 reference genome assembly , ( 2 ) determining the degree of stacking ( and the number of distinct reads 2 ) for all uniquely aligned paired reads 2 , and ( 3 ) assembling each stacks of distinct read 2 regions covered by more than 100 reads 2 with their corresponding read 1 sequences to generate contigs in targeted regions of the genome .   only read pairs containing the 10  bp signature 5-tgcaggtgca-3 sequence at the 5 end of reads 2
the 10  bp signature sequence was removed and paired reads 2 were aligned to the b73 reference genome v2.0 assembly using bowtie , allowing up to 2 mismatches per individual reads .
31.1 million b73 paired reads 2 ( 50.5% of all paired reads 2 containing the 10  bp signature sequence ) and 57.0 million b14 paired reads 2 ( 61.8% ) were aligned uniquely to the b73 reference genome ( table 1 ) .
in addition to uniquely aligned reads 2 , 13.3 million paired reads 2 in b73 and 22.8 million paired reads 2 in b14 were aligned to multiple regions of the genome , indicating that 72.2% and 86.5% of all paired reads 2 for b73 and b14 , respectively , containing the 10  bp signature sequence , align at least once to the b73 reference genome .
finally , a relatively significant fraction of b73 ( 17 million ) and b14 ( 12.3 million ) reads 2 did not align to the b73 reference genome ( table 1 ) .
it is possible that those unaligned reads were excluded from the alignment due to the presence of more than 2 mismatches per individual reads that were caused by natural polymorphisms between lines , sequencing errors , or assembly errors in the b73 reference genome v2.0 assembly .
additional blast search of a subset of the remaining unaligned b73 and b14 read pairs indicates that at least a fraction of the unaligned sequences contain adapter sequences , suggesting the possibility of dna fragments ligating to adapters or adapter - adapter ligations during the library construction process .  
sequencing coverage for all uniquely aligned regions of the b73 genome reference v2.0 assembly was assessed by determining distinct read 2 coverage information ( figure 2 ) .
distinct reads 2 ( or  regions  ) were generated from all uniquely aligned reads 2 , gathering all similar and uniquely aligned read 2 sequences into a unique sequence entry . in both lines , 312,924 and 401,379 regions for b73 and b14 , respectively , are covered by at least 11 uniquely aligned reads 2 and a majority of uniquely aligned reads 2 ( 76.7% for b73 and 88.2% for b14 ) are contained within regions covered by 11 or more reads ( and by extension read pairs ) ( figure 2 ) , potentially enabling high quality genome - wide de novo assembly .
the relatively large percentage of distinct reads 2 covered by one to ten reads ( 110 in figure 2 ) suggests the presence of sequencing singlets with unchecked sequencing errors and/or the possibility of extraneous contamination during the library construction procedure .
nonetheless , that number is low enough not to decrease the overall quality and usefulness of the analysis .
the disparity in number of covered regions between the two lines can be explained at least in part by the higher number of read pairs sequenced for b14 in comparison to b73 and suggest that additional sequencing could further increase the coverage at each individual region adjacent to psti sites in their respective genomes .
the number of regions covered by 11 or more reads 2 in b73 in comparison to b14 also suggests that additional sequencing mostly would increase the number of regions covered by 11 reads 2 or more , rather than generating contigs aligning to still unknown regions of interest , although that number also may vary depending on genome content and organization .  
genomic position overlap between the b73 and b14 distinct read 2 datasets was determined by aligning all distinct reads 2 against the b73 reference genome , using bowtie and allowing up to two mismatches .
the resulting genomic coordinates then were compared between b73 and b14 to assess genomic position overlap .
redundant and nonredundant genomic positions are listed on table 2 ( where distinct reads 2 aligning to more than one region of the genome are considered as redundant ) . interestingly ,
68.8% and 79.8% of distinct reads 2 in b73 and b14 , respectively , and sequenced at least once in one or both genotypes ( i.e. , sequencing coverage of 1x , or one read 2 sequence ) also are sequenced in the other genotype ( table 2 ) . after including reads aligning to multiple regions (
up to two mismatches were allowed when aligning distinct reads using bowtie ) , 78.2% and 84.9% of distinct reads 2 in b73 and b14 , respectively , overlap at the same genomic position .
these data suggests that read 2 sequence data , or contigs resulting from assembling paired read 1 and read 2 data , can be used for direct comparison between genotypes for snp detection and genotyping , in a manner similar to the approach described in deschamps et al . .  
individual read pairs were assembled de novo into contigs with the velvet short read de novo assembler . only distinct reads 2 generated from uniquely aligned reads 2 and with a sequencing coverage greater than 100 were considered and grouped with their corresponding read 1 sequences for de novo assembly .
a total of 183,609 and 129,018 contigs were generated for b14 and b73 , respectively .
read usages in both lines as well as the possibility of multiple contigs generated within the same regions likely explain why a higher number of contigs were generated in comparison to the number of regions with a sequencing coverage greater than 100 . as suggested by figure 3 , 77% of the b73 contigs ( 99,431 contigs ) and 76.2% of the b14 contigs ( 139,968 contigs ) are > 200  bps in length , including several contigs ( 176 for b73 and 310 for b14 ) above 500  bps in length , likely artifacts from the gel size - selection process but an indication nonetheless of the potential of the method if larger dna fragments were to be size - selected on gel .  
contigs were aligned with bowtie to the b73 reference genome v2.0 assembly , allowing up to two mismatches . out of those , 30,763 for b73 and 62,524 for b14
the lower number of uniquely aligned contigs compared to the total number of contigs shown above may be explained by the fact that most contigs did not extend to the stacked and unique read 2 sequences . as shown on table 3 ,
bowtie alignments indicate that a significant fraction of these contigs align perfectly to the assembly ( 0 mismatch ) . as expected , a larger fraction of b73 uniquely mapped contigs ( 72.9% ) align perfectly to the assembly , when compared to b14 , where only 67.6% align with 0 mismatches .
this is due in part to the different genome organization between the two inbred lines and the likely presence of various polymorphisms , including indels in the b14 contigs when compared to the b73 reference assembly .
interestingly , 13.4% of the b73 contigs align to the b73 v2.0 assembly with 1 mismatch and 13.6% align with 2 mismatches , which could be an indication of sequencing or assembly errors in our assembled data or assembly errors in the public data .  
due to the palindromic nature of psti recognition sites , one valuable side effect of this method is the possibility of generating contigs immediately adjacent to the 5 end and the 3 end of the same psti site . the expected outcome would be that such adjacent contigs are in most cases separated by a variable number of bases , as determined by the reference sequence they are aligned against , assuming that contigs are extended from ends opposite the psti ends in dna fragments assembled together , creating 
contig pairs  spanning larger regions of the genome useful for marker development and de novo   assembly applications .
alignments of contig sequences to the b73 reference genome sequence suggest that 6,495 of the > 100  bps uniquely mapped b73 contigs and 28,605 of the > 100  bps uniquely mapped b14 contigs generate contig pairs separated by a distance less than 1,000  bps , and centered on a given psti site in their respective genome .
in addition , out of those , 680 b73 contigs and 3356 b14 contigs are immediately adjacent contigs centered on a unique psti site .
it is expected that deeper sequencing would increase the number of contig pairs per experiment by increasing sequencing coverage and contig length for each discrete genomic regions .  
one application of generating longer contigs from individual read pairs , anchored at a given psti site is the possibility of detecting new snp in complex genomes and using long contig sequences for primer design and snp assay development . to determine the feasibility of using contig assemblies for snp detection , a subset of b73 and b14 contigs with one or more mismatch to the b73 reference genome was selected randomly .
mismatches detected by bowtie ( putative  snps  ) then were confirmed by sanger sequencing of pcr products generated using primers derived from their respective contig sequences .
a total of 41 amplicons , containing 50 putative snps , were sequenced in b14 .
a total of 45 snps were confirmed via sanger sequencing , while the remaining 5 showed sequencing peaks typical of a heterozygous call , possibly the results of assembly errors in the b73 reference genome of highly conserved regions differing by only one base , assembly errors in the b73 contigs generated by assembling read 2 data from different regions grouped by similar read 1 sequences , or duplication of the corresponding regions in the b14 genome in relation to b73 .
these data nonetheless confirm that the putative snps are present in the dna fragments of interest , and the quality and value of contig sequences for detecting snps in a complex genome - like maize .   to determine the origin of mismatches in the b73 contig sequences aligned to the b73 reference genome , a similar approach was used , where 40 amplicons , containing 52 putative snps , were sequenced via sanger sequencing of pcr products . out of the 52 putative snps ,
only 20 were confirmed , while 17 were not confirmed , exhibiting a monomorphic pattern at the base position of interest , and 15 exhibit sequencing peaks characteristic of a heterozygous call .
the confirmation of 20 snps and the existence of 15 putative snps exhibiting sequencing peaks characteristic of a heterozygous call suggests assembly errors in the b73 reference genome or natural variations between dna materials .
alternatively , errors in the b73 reference genome assembly could have led to artificial grouping of reads 1 data , led by the fact that uniquely aligned read 2 sequences actually correspond in reality to separate regions of the genomes .
a rapid blast search of the b73 contig sequences carrying unconfirmed snps or snps exhibiting sequencing peaks characteristic of a heterozygous call indicates that all or parts of a majority of the contigs ( 63% , data not shown ) align with more than one b73 bac clone sequence ( mapped to different chromosomes ) , confirming the potential value of filtering such contig sequences prior to snp discovery or detection .
a revised rad - pe protocol for generating long contig sequences mapping to discrete regions of the maize genome has been tested and implemented .
contig length essentially varies with the size of the dna fragments selected during the library construction procedure and is expected to be limited by the maximum length of the dna fragments sequenced on the illumina platform ( separate experiments have shown that dna fragments up to ~1  kbps can be sequenced effectively on the illumina genome analyzer ) .
resulting contigs can be used for several applications , including snp discovery and marker assay development .
they also are expected to facilitate snp discovery efforts in complex genomes lacking robust reference sequence information , not only in plants but also any other species of interest .
 the use of next - generation dna sequencing technologies has greatly facilitated reference - guided variant detection in complex plant genomes . 
 however , complications may arise when regions adjacent to a read of interest are used for marker assay development , or when reference sequences are incomplete , as short reads alone may not be long enough to ascertain their uniqueness . here 
 , the possibility of generating longer sequences in discrete regions of the large and complex genome of maize is demonstrated , using a modified version of a paired - end rad library construction strategy . 
 reads are generated from dna fragments first digested with a methylation - sensitive restriction endonuclease , sheared , enriched with biotin and a selective pcr amplification step , and then sequenced at both ends . 
 sequences are locally assembled into contigs by subgrouping pairs based on the identity of the read anchored by the restriction site . 
 this strategy applied to two maize inbred lines ( b14 and b73 ) generated 183,609 and 129,018 contigs , respectively , out of which at least 76% were > 200  bps in length . 
 a subset of putative single nucleotide polymorphisms from contigs aligning to the b73 reference genome with at least one mismatch was resequenced , and 90% of those in b14 were confirmed , indicating that this method is a potent approach for variant detection and marker development in species with complex genomes or lacking extensive reference sequences . 
1. Introduction
2. Materials and Methods
3. Results and Discussion
4. Conclusions