the ability of certain dna sequences to adopt alternative conformations , in addition to the canonical watson 
indeed , a large number of studies have documented the formation of alternative ( non - b ) dna structures by biophysical methods , including x - ray crystallography ( 24 ) , nuclear magnetic resonance ( nmr ) spectroscopy ( 5 ) and circular dichroism ( 6 ) .
other methods , such as the detection of single - stranded bases upon non - b dna structure formation by chemical and enzymatic probes and the relaxation of negative supercoiling by two - dimensional gel electrophoresis have played a major role in revealing the formation of non - b dna conformations in biological systems ( 79 ) .
, inverted repeats can adopt cruciform structures , runs of alternating purine  pyrimidine bases are able to switch from the right - handed b- to the left - handed z - dna helix , homo(purinepyrimidine ) tracts with mirror repeat symmetry may fold into several types of intramolecular triplexes , four sets of three , four or five guanines , each interrupted by 17 bases , can form highly stable , polymorphic , quadruplex structures and direct repeats can give rise to loops or hairpins through the misalignment of complementary strands , also known as slipped structures ( 10 ) .
a number of bioinformatic searches have been conducted with the aim of identifying the biological relevance of putative non - b dna structures in mammalian and other genomes ( 1 ) .
these studies support the notion that the secondary structure conformational domain , rather than the underlying sequence symmetry , often contributes to the control of diverse biological functions , including replication , transcription , immune response ( 11 ) , recombination and antigenic variation in human pathogens ( 1,12 ) .
concomitant to this notion , a number of studies have provided circumstantial evidence for the involvement of dna secondary structures in inducing genetic instability , both in model systems ( 1315 ) and in association with human genetic disease ( 1620 ) , including genomic regions that do not contain known genes , suggesting that deeper functional annotation across these regions is warranted .
therefore , the need has arisen to provide the scientific community with a tool that offers a systematic cataloguing of all predicted sequences currently known to potentially form alternative dna conformations .
the non - b db database bridges this gap by providing a resource for searching , mapping and comparing non - b dna - forming motifs among various mammalian species .
to date , several reports have detailed methods aimed at enumerating and evaluating predicted non - b dna - forming elements from genomic sequences , including quadbase ( 21 ) , tts ( 22 ) , trf ( 23 ) and others ( documented at http://nonb.abcc.ncifcrf.gov/resources/ ) .
these reports use various consensus - based scanning methods for identifying one specific class of predicted non - b dna structure . in some cases ,
the identified motifs are screened for the presence of other overlapping functional motifs , such as sp1 binding sites and cpg islands ( 24 ) .
in other cases , the resulting motifs can be searched by genomic position and scanned for the presence of other nearby non - b dna predicted features [ e.g. triplex sequences near quadruplexes ( 22 ) ] .
more recently , analyses that incorporate thermodynamic values into the overall scoring method ( 2527 ) have been reported .
together , these resources provide an important , yet partial , view into the complexities of locating and characterizing the many different sequence motifs that have the potential of forming non - b dna structures .
our database expands on these functionalities by including all classes of predicted non - b dna - forming sequences and by using the latest genome assemblies of human , mouse and other mammalian species .
the non - b dna data are available with current genomic annotation data and polymorphism information .
importantly , non - b db provides the capacity to visualize the data in a genomic context that is fully integrated with other genomic features , such as genes and single - nucleotide polymorphisms ( snps ) .
the same interface allows for the users to upload their own annotation data , which are displayed alongside the in - house data through the polybrowse and ucsc interfaces .
one of the main difficulties in developing and evaluating algorithms that predict the likely candidates for each class of non - b structures is the lack of large collections of experimental data that have validated their formation in vivo .
although most non - b dna structures can be formed under in vitro conditions , the identification of such conformations in vivo and the elucidation of parameters that govern their b to non - b equilibria have presented formidable challenges .
in addition , these equilibria are influenced by local superhelical density , the presence of nearby dna unwinding element complexes ( dues ) ( 28 ) , the transcriptional status , nucleosome assembly and other tissue / temporally regulated biological processes . in light of these considerations , we have taken the approach of using rather broad and general identification methods based exclusively on sequence features ; thus , although subsequent filtering of the sampled data is straightforward because of the flexibility provided by the database , our current criteria are expected to include a subset of both false positive and negative hits .
we have previously reported the construction of a database containing information on mouse indel polymorphisms ( 30 ) .
herein , we have extended that system to include motifs with the potential to form non - b dna structures .
a number of studies in vitro ( 3134 ) and in vivo ( 29,3538 ) have indicated that the structural transition from b to non - b dna is assisted by unrestrained negative supercoiling . in mammalian cells ,
the global steady - state levels of negative supercoiling vary depending on chromosomal location ( 39 ) , but are expected to increase transiently by processes , such as transcription , replication and repair , that entail separation of the complementary strands and thus affect nucleosome occupancy ( 29,38,4042 ) . however , because the kinetics of these processes may vary among cell types and various developmental stages , an assessment of the probability that a defined chromosomal sequence might exist in the non - b form is currently not available . indeed , only limited overlap has been reported between the predicted z - dna formation based on in silico thermodynamic predictions and genomic loci bound to the z domain of adar1 , which displays high specificity for z - dna ( 43 ) .
thus , a combination of factors , including nucleosome occupancy , negative supercoiling , matrix attachment sites , replication , transcription and repair may underlie b to non - b equilibria in vivo . in the absence of such information
, our search algorithms were based solely on sequence relationships derived from in vitro data .
the general approach involves running a scanning application for each specific predicted non - b dna class against each chromosome ( table 1 ) , including g - quadruplex motifs , alternating purine  pyrimidine sequences , mirror repeats , inverted repeats and direct repeats . although the  mirror repeat class as a whole has not been reported to form specific non - b dna structures , it is included in the database as it is used as a first step in the identification of triplex - forming motifs , i.e. the subset of mirror repeats with purine / pyrimidine content . 

table 1.criteria for predicting non - b dna - forming motifs in non - b dbdna featuresearch criteriasubset of 
dna feature forming non - b dnasearch criteria for  subset of dna featureinverted repeatrepeat : 10100  ntcruciform motifrepeat : 10100  ntspacer : 0100  ntspacer : 03  ntmirror repeatrepeat : 10100  nttriplex motifrepeat : 10100 ror y ntspacer : 0100  ntspacer : 08  ntdirect repeatrepeat : 1050  ntslipped motifrepeat : 1050  ntspacer : 05  ntspacer : 0  ntz - dna repeat5 units of cg / tg or cg / ca repeatswhole setas per the whole setg - quadruplex forming repeatfour identical blocks of ( 37 ) g nt , each block separated by 17  ntwhole setas per the whole seta - phased repeat3 runs of a - tracts with
10-bp phasingwhole setas per the whole setinverted repeat : a pair of dna sequences , each 10100  nt in length and separated by a spacer of 0100  nt , whose sequence composition on the same strand of dna is such that the bases of the first repeat , when read in the 53 orientation , are complementary to those of the second repeat read in the 35 orientation . the term 
crick hydrogen bonding scheme , whereby a only pairs with t and c only pairs with g. only perfect inverted repeats that conform to this watson 
crick pairing scheme are considered.cruciform motif : the subset of inverted repeat sequences in which the  spacer comprises 03 bases ; due to their proximity , this subset of inverted repeat sequences may fold - back and form intramolecular , antiparallel , double helices stabilized by watson  crick hydrogen bonds , i.e. a cruciform structure ( 1,34).mirror repeat : a pair of dna sequences , each 10100  nt in length and separated by a spacer of 0100  nt , whose sequence composition on the same strand of dna is such that the bases of the first repeat , when read in the 53 orientation , are identical to those of the second repeat read in the 35 orientation ( palindrome ) ; only perfectly matching repeats are included.triplex motif : the subset of mirror repeat sequences comprising only purines ( r  =  a and g ) [ or pyrimidines ( y  =  c and t ) ] on the same strand of dna , and which are separated by few ( 08 ) nt (  spacer ) . these motifs are able to form various intramolecular three - stranded ( triplex , h - dna ) isoforms stabilized by hoogsteen hydrogen bonds ( 1,52,53 ) .
only ry - containing mirror repeats that may yield a : at and g : gc base triplets ( colon indicates hoogsteen hydrogen bonded bases ; dot indicates watson  crick hydrogen bonded bases ) for the r : ry type of intramolecular triplexes and t : at and c : gc triplets for the y : ry type of intramolecular triplexes are included since these are considered the most stable triplet combinations.direct repeat : two tracts of dna , each comprising 1050  nt and separated by 05  nt , having the same sequence composition.slipped motif : the subset of direct repeat sequences without a spacer ( tandem repeats ) ; when aligned in an out - of - register fashion , tandem repeats may give rise to single - stranded loops and/or hairpins ( 1).z - dna motif : five or more tandem repeats , each comprising an alternating pyrimidine 
purine dinucleotide motif , in which the pattern yg is maintained on at least one of the dna strands ; examples include ( cgcg)6 , ( catg)5 and [ ( tg)3(cg)4(cg)4(ca)3 ] ; these motifs may adopt the left - handed z - dna conformation ( 3,54).g - quadruplex - forming repeat : four blocks , each containing the same number ( n ) of g bases ( n can vary from 3 to 7 ) , on the plus or minus strand , separated by 17  nt ; this type of dna sequence may adopt quadruplex structures ( 2 ) ; overlapping tracts of four g - blocks are also considered.a - phased repeat : three runs of a bases ( a - tracts ) in phase with the helical pitch of the dna double - helix , i.e. 10  bp ; an a - tract is defined as a set of at base - pairs without a tpa step ( 47,5557 ) ; three or more tracts of a37 , t37 , aaattt , aaatttt and aaaattt ( in any combination ) on the plus or minus strand , whose centers are separated by 10 bases , are considered ; since a - tracts induce static bends in the dna double helix , the overall dna superhelix is expected to display either a left - handed or a right - handed writhe ( 47,5557 ) ; as mentioned , all the search criteria used herein do not allow for interruptions in the repeats and no thermodynamic information was factored - in in the algorithms used .
criteria for predicting non - b dna - forming motifs in non - b db inverted repeat : a pair of dna sequences , each 10100  nt in length and separated by a spacer of 0100  nt , whose sequence composition on the same strand of dna is such that the bases of the first repeat , when read in the 53 orientation , are complementary to those of the second repeat read in the 35 orientation . the term 
crick hydrogen bonding scheme , whereby a only pairs with t and c only pairs with g. only perfect inverted repeats that conform to this watson 
spacer comprises 03 bases ; due to their proximity , this subset of inverted repeat sequences may fold - back and form intramolecular , antiparallel , double helices stabilized by watson  crick hydrogen bonds , i.e. a cruciform structure ( 1,34 ) . mirror repeat : a pair of dna sequences , each 10100  nt in length and separated by a spacer of 0100  nt , whose sequence composition on the same strand of dna is such that the bases of the first repeat , when read in the 53 orientation , are identical to those of the second repeat read in the 35 orientation ( palindrome ) ; only perfectly matching repeats are included . triplex
motif : the subset of mirror repeat sequences comprising only purines ( r  =  a and g ) [ or pyrimidines ( y  =  c and t ) ] on the same strand of dna , and which are separated by few ( 08 ) nt (  spacer ) .
these motifs are able to form various intramolecular three - stranded ( triplex , h - dna ) isoforms stabilized by hoogsteen hydrogen bonds ( 1,52,53 ) .
only ry - containing mirror repeats that may yield a : at and g : gc base triplets ( colon indicates hoogsteen hydrogen bonded bases ; dot indicates watson  crick hydrogen bonded bases ) for the r : ry type of intramolecular triplexes and t : at and c : gc triplets for the y : ry type of intramolecular triplexes are included since these are considered the most stable triplet combinations .
direct repeat : two tracts of dna , each comprising 1050  nt and separated by 05  nt , having the same sequence composition . slipped motif : the subset of direct repeat sequences without a spacer ( tandem repeats ) ; when aligned in an out - of - register fashion , tandem repeats may give rise to single - stranded loops and/or hairpins ( 1 ) .
z - dna motif : five or more tandem repeats , each comprising an alternating pyrimidine 
purine dinucleotide motif , in which the pattern yg is maintained on at least one of the dna strands ; examples include ( cgcg)6 , ( catg)5 and [ ( tg)3(cg)4(cg)4(ca)3 ] ; these motifs may adopt the left - handed z - dna conformation ( 3,54 ) .
g - quadruplex - forming repeat : four blocks , each containing the same number ( n ) of g bases ( n can vary from 3 to 7 ) , on the plus or minus strand , separated by 17  nt ; this type of dna sequence may adopt quadruplex structures ( 2 ) ; overlapping tracts of four g - blocks are also considered .
a - phased repeat : three runs of a bases ( a - tracts ) in phase with the helical pitch of the dna double - helix , i.e. 10  bp ; an a - tract is defined as a set of at base - pairs without a tpa step ( 47,5557 ) ; three or more tracts of a37 , t37 , aaattt , aaatttt and aaaattt ( in any combination ) on the plus or minus strand , whose centers are separated by 10 bases , are considered ; since a - tracts induce static bends in the dna double helix , the overall dna superhelix is expected to display either a left - handed or a right - handed writhe ( 47,5557 ) ; as mentioned , all the search criteria used herein do not allow for interruptions in the repeats and no thermodynamic information was factored - in in the algorithms used . the output file in gff format ( http://nonb.abcc.ncifcrf.gov/faqs/ )
the data from all such scans are merged and can be queried and displayed using our local instance of gbrowse called polybrowse ( 44 ) at http://pbrowse3.abcc.ncifcrf.gov/cgi-bin/gb2/gbrowse/human_37 and several gff - based query tools at http://nonb.abcc.ncifcrf.gov ( figure 1 ) .
importantly , the result pages produced from the queries contain links that allow the user to switch to the genome browser view of that feature , as well as a view that provides the sequence and other annotations for each feature . 

non - b db user interface ( a ) three query modes are available at the genomic database search tools page of the non - b db web interface : search by feature ( shown ) , search by feature attributes , and search by location ( a feature browser ) . in this example , all non - b dna motifs were queried by the gene symbol myc .
( b ) query results and links to polybrowse , ucsc genome browser , biodbnet ( 51 ) , etc . users may download the results in gff format and tab delimited format .
direct access to the sequence and other annotation information for each of the features is available by clicking on the 
non - b db user interface ( a ) three query modes are available at the genomic database search tools page of the non - b db web interface : search by feature ( shown ) , search by feature attributes , and search by location ( a feature browser ) . in this example , all non - b dna motifs were queried by the gene symbol myc .
( b ) query results and links to polybrowse , ucsc genome browser , biodbnet ( 51 ) , etc .
direct access to the sequence and other annotation information for each of the features is available by clicking on the  a link ( red box ) .
these data represent the basis for the non - b dna annotation information for each species .
the scanning criteria do not allow for mismatches within the repeat segments ; however , this feature may be added as information becomes available as to the acceptable structural tolerances for each mismatch case .
also , currently not included are very large palindromes ( > 100  kb ) , such as those that characterize the y chromosome and whose recombination is known to lead to spermatogenic failure ( 45,46 ) .
nevertheless , some aspects related to the presence of mismatches are presented in the polymorphism analysis described below . after scanning across different mammalian genomes ,
the numbers of each of the predicted classes of non - b dna structure - forming motifs appear to be quite variable ( table 2 ) . 

table 2.statistics for dna repeats and non - b dna - forming motifs in non - b dbdna featurehuman 37mouse 37dog 2chimp 2macaque 1g - quadruplex forming repeat374  545559  280492  535314  171298  142inverted repeat1  044  533801  242814 
080998  249843  889cruciform motif197  910188  532172  032190  736128  334direct repeat871  0451  593  107968  955787  335765  798slipped motif347  969695  150404  750314  516305  285mirror repeat16  51  7233  431  4861  829  86714  85  13514  55  025triplex motif1  79  623618  928336  6421 
05  6401  40  580z - dna repeat294  320690  276261  012278  928280  982a - phased repeat1 
130  7319  09  6531  241  0821  085  5911  098  030for the current releases of the five mammalian genomes indicated , the motif searches were performed and the number of features for each class was counted . according to table 1 , the cruciform motifs represent a subset of the inverted repeat class , the slipped motifs represent a subset of the direct repeat class and the triplex motif represents a subset of the mirror repeat class .
statistics for dna repeats and non - b dna - forming motifs in non - b db for the current releases of the five mammalian genomes indicated , the motif searches were performed and the number of features for each class was counted . according to table 1 , the cruciform motifs represent a subset of the inverted repeat class , the slipped motifs represent a subset of the direct repeat class and the triplex motif represents a subset of the mirror repeat class . as the overall base composition between different mammalian genomes is rather similar ( data not shown ) , the observed differences in the numbers of predicted non - b dna motifs could simply result from the altered arrangement of bases from one species to another .
alternatively , variations in the population of classes of repetitive elements ( sine , line , etc . ) among species , or other unknown features , might also contribute to the observed differences .
this interspecies variability appears to be uniformly distributed along the entire chromosomes , rather than concentrated in large repetitive clusters ( data not shown ) .
whether these differences play any role or contribute to conferring species - specific differences remains to be investigated .
a caveat concerning the simple assessment and comparison of the number of non - b dna - forming repeats among species relates to the criteria used and the counting method .
for example , in the g - quadruplex forming sequences , the pattern of a run of 3gs followed by 17 bases repeated four times can be extended , as long as more runs of gs are encountered , resulting in a single cluster that has the potential to form many substructures .
although our approach identifies this finding as a single cluster in the database , separate database tables are provided , in which all possible permutations of the sequence that satisfies the consensus sequence are reported .
in addition to the non - b dna predicted motifs , the database contains other features of the dna , such as phased a - tracts that impart static bends to the double - helix and may be involved in nucleosome assembly ( 47 ) , simple tandem repeats ( str ) including triplet repeats whose expansions cause a number of neuromuscular disorders ( 20 ) and poly(purinepyrimidine ) tracts , which are characterized by high stacking interactions ( 48 ) .
in addition , ncbi - derived features , such as genes , snps and repeatmasker ( http://www.repeatmasker.org/ ) elements are also included .
this integrated information is critical not only for guiding the user visually , but also for enabling queries that combine 
classes , such as  exons containing predicted  z - dna forming sequences , etc .
one of the main features of the non - b db is the ability to compare different mammalian genomes for the presence of non - b dna - forming motifs .
figure 2 illustrates this capability by comparing the presence of g - quadruplex forming motifs in the region upstream of the myc locus across the human , chimp , macaque , dog and mouse reference genomes . in order to view syntenic regions in other genomes , the liftover application from the ucsc website
areas where a syntenic match failed to be identified ( i.e. that region was absent in the other genome , or mapped redundantly ) do not show a link to that species .
other non - b dna tracks available in polybrowse will be described in more detail elsewhere ( cer et .
figure 2.(a ) polybrowse view of g - quadruplex forming repeats in the myc gene with species syntenic information .
visualization of cross - species information in polybrowse : the genomic region shown in figure 1 is displayed herein to illustrate the syntenic capability . in this case
, the g - quadruplex forming repeat track shows the locations of these features in the region near the beginning of the myc gene . the additional track illustrating the liftover1k blocks is also reported . in that track ,
when the user moves the cursor over the objects , a popup window appears that contains links ( a ) to the syntenic locations in the other available species .
( b  e ) results upon clicking on the links for chimp , macaque , mouse and dog , respectively .
currently , the liftover feature only allows moving from the human to the other species but not among the non - human species .
( a ) polybrowse view of g - quadruplex forming repeats in the myc gene with species syntenic information .
visualization of cross - species information in polybrowse : the genomic region shown in figure 1 is displayed herein to illustrate the syntenic capability .
in this case , the g - quadruplex forming repeat track shows the locations of these features in the region near the beginning of the myc gene . the additional track illustrating the liftover1k blocks is also reported .
in that track , when the user moves the cursor over the objects , a popup window appears that contains links ( a ) to the syntenic locations in the other available species .
( b  e ) results upon clicking on the links for chimp , macaque , mouse and dog , respectively .
currently , the liftover feature only allows moving from the human to the other species but not among the non - human species . the computed non - b dna forming elements are likely to be under - represented in our reference genome as their underlying repeats may be polymorphic among individuals .
because this type of information may be critical in the context of gene regulation or predisposition to disease ( 48 ) , we used a specific parser to scan both the reference human genome as well as additional sequence sources , such as trace reads from the trace archive ( http://www.ncbi.nlm.nih.gov/traces/trace.cgi ) and contigs ( http://www.ncbi.nlm.nih.gov/projects/wgs/wgsprojectlist.cgi ) from personal genome projects ( 49 ) , for matches to the non - b dna motifs .
each match found in either the reference or alternate source is then scored for being polymorphic or not .
of the sites identified as polymorphic , a second evaluation is made to determine whether the polymorphism would affect the motif underlying the putative non - b dna structure .
the results of this scan are incorporated into the database as a series of separate tracks ( figure 3b , trace gplex tracks ) .
additional information can be gathered by extending this type of analysis to sequence alignments using closely related species .
currently , only the g - quadruplex forming motif supports this type of query . 

figure 3.(a ) polybrowse view of some of the non - b dna features in the human myc gene .
links to polybrowse and ucsc browsers : ( a ) display of polybrowse rendering the myc gene and its promoter with some of the non - b dna motifs tracks with the pupy [ poly(purinepyrimidine ) ] and polymorphism tracks turned on .
several other tracks include features computed at the advanced biomedical computing center ( abcc ) , such as strs , base composition , physical dna characteristics , mapping , as well as the ncbi derived features , including genes , snps , cytogenic markers assembly information , repeatmasker elements , etc .
( b ) display of the myc gene in polybrowse showing some of the polymorphism information in non - b db .
below the gene and mrna tracks , the searched motifs in the reference sequence are displayed ( teal ) .
an additional track shows specific g - quadruplex motifs in which all of the observed trace files contained mutations disrupting the g - quadruplex motif .
( c ) display of the myc gene at the ucsc human genome browser ( grch37/hg19 assembly ) as linked to non - b db from the search shown in figure 2b . some of the non - b dna motifs from non - b db can be seen in the red rectangular box .
( a ) polybrowse view of some of the non - b dna features in the human myc gene .
links to polybrowse and ucsc browsers : ( a ) display of polybrowse rendering the myc gene and its promoter with some of the non - b dna motifs tracks with the pupy [ poly(purinepyrimidine ) ] and polymorphism tracks turned on .
several other tracks include features computed at the advanced biomedical computing center ( abcc ) , such as strs , base composition , physical dna characteristics , mapping , as well as the ncbi derived features , including genes , snps , cytogenic markers assembly information , repeatmasker elements , etc .
( b ) display of the myc gene in polybrowse showing some of the polymorphism information in non - b db .
below the gene and mrna tracks , the searched motifs in the reference sequence are displayed ( teal ) .
an additional track shows specific g - quadruplex motifs in which all of the observed trace files contained mutations disrupting the g - quadruplex motif .
( c ) display of the myc gene at the ucsc human genome browser ( grch37/hg19 assembly ) as linked to non - b db from the search shown in figure 2b . some of the non - b dna motifs from non - b db can be seen in the red rectangular box . in order to provide access to the back - end database
the first is a bioperl ( 50 ) set of methods , which is used to query genome databases in various ways , such as by position , by class , or by attribute .
this same set of utilities is used within the context of polybrowse ( 44 ) , so that visualization of the genomic features is made available .
in addition to linking the outputs from the query tools to the browser for visualization , we also provide links allowing the returned data to be displayed in the familiar ucsc ( http://genome.ucsc.edu/ ) genome browser ( figure 3c ) , as well as links to our biodbnet database warehouse ( 51 ) , which contains gene - centric information derived from several sources , and additional links .
to date , several reports have detailed methods aimed at enumerating and evaluating predicted non - b dna - forming elements from genomic sequences , including quadbase ( 21 ) , tts ( 22 ) , trf ( 23 ) and others ( documented at http://nonb.abcc.ncifcrf.gov/resources/ ) .
these reports use various consensus - based scanning methods for identifying one specific class of predicted non - b dna structure . in some cases ,
the identified motifs are screened for the presence of other overlapping functional motifs , such as sp1 binding sites and cpg islands ( 24 ) .
in other cases , the resulting motifs can be searched by genomic position and scanned for the presence of other nearby non - b dna predicted features [ e.g. triplex sequences near quadruplexes ( 22 ) ] .
more recently , analyses that incorporate thermodynamic values into the overall scoring method ( 2527 ) have been reported .
together , these resources provide an important , yet partial , view into the complexities of locating and characterizing the many different sequence motifs that have the potential of forming non - b dna structures .
our database expands on these functionalities by including all classes of predicted non - b dna - forming sequences and by using the latest genome assemblies of human , mouse and other mammalian species .
the non - b dna data are available with current genomic annotation data and polymorphism information .
importantly , non - b db provides the capacity to visualize the data in a genomic context that is fully integrated with other genomic features , such as genes and single - nucleotide polymorphisms ( snps ) .
the same interface allows for the users to upload their own annotation data , which are displayed alongside the in - house data through the polybrowse and ucsc interfaces .
one of the main difficulties in developing and evaluating algorithms that predict the likely candidates for each class of non - b structures is the lack of large collections of experimental data that have validated their formation in vivo .
although most non - b dna structures can be formed under in vitro conditions , the identification of such conformations in vivo and the elucidation of parameters that govern their b to non - b equilibria have presented formidable challenges .
in addition , these equilibria are influenced by local superhelical density , the presence of nearby dna unwinding element complexes ( dues ) ( 28 ) , the transcriptional status , nucleosome assembly and other tissue / temporally regulated biological processes . in light of these considerations , we have taken the approach of using rather broad and general identification methods based exclusively on sequence features ; thus , although subsequent filtering of the sampled data is straightforward because of the flexibility provided by the database , our current criteria are expected to include a subset of both false positive and negative hits .
we have previously reported the construction of a database containing information on mouse indel polymorphisms ( 30 ) .
herein , we have extended that system to include motifs with the potential to form non - b dna structures .
a number of studies in vitro ( 3134 ) and in vivo ( 29,3538 ) have indicated that the structural transition from b to non - b dna is assisted by unrestrained negative supercoiling . in mammalian cells ,
the global steady - state levels of negative supercoiling vary depending on chromosomal location ( 39 ) , but are expected to increase transiently by processes , such as transcription , replication and repair , that entail separation of the complementary strands and thus affect nucleosome occupancy ( 29,38,4042 ) . however , because the kinetics of these processes may vary among cell types and various developmental stages , an assessment of the probability that a defined chromosomal sequence might exist in the non - b form is currently not available . indeed , only limited overlap has been reported between the predicted z - dna formation based on in silico thermodynamic predictions and genomic loci bound to the z domain of adar1 , which displays high specificity for z - dna ( 43 ) .
thus , a combination of factors , including nucleosome occupancy , negative supercoiling , matrix attachment sites , replication , transcription and repair may underlie b to non - b equilibria in vivo . in the absence of such information ,
the general approach involves running a scanning application for each specific predicted non - b dna class against each chromosome ( table 1 ) , including g - quadruplex motifs , alternating purine  pyrimidine sequences , mirror repeats , inverted repeats and direct repeats .
although the  mirror repeat class as a whole has not been reported to form specific non - b dna structures , it is included in the database as it is used as a first step in the identification of triplex - forming motifs , i.e. the subset of mirror repeats with purine / pyrimidine content . 

table 1.criteria for predicting non - b dna - forming motifs in non - b dbdna featuresearch criteriasubset of 
dna feature forming non - b dnasearch criteria for  subset of dna featureinverted repeatrepeat : 10100  ntcruciform motifrepeat : 10100  ntspacer : 0100  ntspacer : 03  ntmirror repeatrepeat : 10100  nttriplex motifrepeat : 10100 ror y ntspacer : 0100  ntspacer : 08  ntdirect repeatrepeat : 1050  ntslipped motifrepeat : 1050  ntspacer : 05  ntspacer : 0  ntz - dna repeat5 units of cg / tg or cg / ca repeatswhole setas per the whole setg - quadruplex forming repeatfour identical blocks of ( 37 ) g nt , each block separated by 17  ntwhole setas per the whole seta - phased repeat3 runs of a - tracts with 10-bp phasingwhole setas per the whole setinverted repeat : a pair of dna sequences , each 10100  nt in length and separated by a spacer of 0100  nt , whose sequence composition on the same strand of dna is such that the bases of the first repeat , when read in the 53 orientation , are complementary to those of the second repeat read in the 35 orientation . the term  complementary refers to the watson 
crick hydrogen bonding scheme , whereby a only pairs with t and c only pairs with g. only perfect inverted repeats that conform to this watson 
crick pairing scheme are considered.cruciform motif : the subset of inverted repeat sequences in which the  spacer comprises 03 bases ; due to their proximity , this subset of inverted repeat sequences may fold - back and form intramolecular , antiparallel , double helices stabilized by watson  crick hydrogen bonds , i.e. a cruciform structure ( 1,34).mirror repeat : a pair of dna sequences , each 10100  nt in length and separated by a spacer of 0100  nt , whose sequence composition on the same strand of dna is such that the bases of the first repeat , when read in the 53 orientation , are identical to those of the second repeat read in the 35 orientation ( palindrome ) ; only perfectly matching repeats are included.triplex motif : the subset of mirror repeat sequences comprising only purines ( r  =  a and g ) [ or pyrimidines ( y  =  c and t ) ] on the same strand of dna , and which are separated by few ( 08 ) nt (  spacer ) .
these motifs are able to form various intramolecular three - stranded ( triplex , h - dna ) isoforms stabilized by hoogsteen hydrogen bonds ( 1,52,53 ) .
only ry - containing mirror repeats that may yield a : at and g : gc base triplets ( colon indicates hoogsteen hydrogen bonded bases ; dot indicates watson  crick hydrogen bonded bases ) for the r : ry type of intramolecular triplexes and t : at and c : gc triplets for the y : ry type of intramolecular triplexes are included since these are considered the most stable triplet combinations.direct repeat : two tracts of dna , each comprising 1050  nt and separated by 05  nt , having the same sequence composition.slipped motif : the subset of direct repeat sequences without a spacer ( tandem repeats ) ; when aligned in an out - of - register fashion , tandem repeats may give rise to single - stranded loops and/or hairpins ( 1).z - dna motif : five or more tandem repeats , each comprising an alternating pyrimidine 
purine dinucleotide motif , in which the pattern yg is maintained on at least one of the dna strands ; examples include ( cgcg)6 , ( catg)5 and [ ( tg)3(cg)4(cg)4(ca)3 ] ; these motifs may adopt the left - handed z - dna conformation ( 3,54).g - quadruplex - forming repeat : four blocks , each containing the same number ( n ) of g bases ( n can vary from 3 to 7 ) , on the plus or minus strand , separated by 17  nt ; this type of dna sequence may adopt quadruplex structures ( 2 ) ; overlapping tracts of four g - blocks are also considered.a - phased repeat : three runs of a bases ( a - tracts ) in phase with the helical pitch of the dna double - helix , i.e. 10  bp ; an a - tract is defined as a set of at base - pairs without a tpa step ( 47,5557 ) ; three or more tracts of a37 , t37 , aaattt , aaatttt and aaaattt ( in any combination ) on the plus or minus strand , whose centers are separated by 10 bases , are considered ; since a - tracts induce static bends in the dna double helix , the overall dna superhelix is expected to display either a left - handed or a right - handed writhe ( 47,5557 ) ; as mentioned , all the search criteria used herein do not allow for interruptions in the repeats and no thermodynamic information was factored - in in the algorithms used .
criteria for predicting non - b dna - forming motifs in non - b db inverted repeat : a pair of dna sequences , each 10100  nt in length and separated by a spacer of 0100  nt , whose sequence composition on the same strand of dna is such that the bases of the first repeat , when read in the 53 orientation , are complementary to those of the second repeat read in the 35 orientation .
crick hydrogen bonding scheme , whereby a only pairs with t and c only pairs with g. only perfect inverted repeats that conform to this watson 
cruciform motif : the subset of inverted repeat sequences in which the  spacer comprises 03 bases ; due to their proximity , this subset of inverted repeat sequences may fold - back and form intramolecular , antiparallel , double helices stabilized by watson  crick hydrogen bonds , i.e. a cruciform structure ( 1,34 ) . mirror repeat : a pair of dna sequences , each 10100  nt in length and separated by a spacer of 0100  nt , whose sequence composition on the same strand of dna is such that the bases of the first repeat , when read in the 53 orientation , are identical to those of the second repeat read in the 35 orientation ( palindrome ) ; only perfectly matching repeats are included . triplex motif : the subset of mirror repeat sequences comprising only purines ( r  =  a and g ) [ or pyrimidines ( y  =  c and t ) ] on the same strand of dna , and which are separated by few ( 08 ) nt (  spacer ) .
these motifs are able to form various intramolecular three - stranded ( triplex , h - dna ) isoforms stabilized by hoogsteen hydrogen bonds ( 1,52,53 ) .
only ry - containing mirror repeats that may yield a : at and g : gc base triplets ( colon indicates hoogsteen hydrogen bonded bases ; dot indicates watson  crick hydrogen bonded bases ) for the r : ry type of intramolecular triplexes and t : at and c : gc triplets for the y : ry type of intramolecular triplexes are included since these are considered the most stable triplet combinations .
direct repeat : two tracts of dna , each comprising 1050  nt and separated by 05  nt , having the same sequence composition . slipped motif : the subset of direct repeat sequences without a spacer ( tandem repeats ) ; when aligned in an out - of - register fashion , tandem repeats may give rise to single - stranded loops and/or hairpins ( 1 ) . z - dna motif : five or more tandem repeats , each comprising an alternating pyrimidine 
purine dinucleotide motif , in which the pattern yg is maintained on at least one of the dna strands ; examples include ( cgcg)6 , ( catg)5 and [ ( tg)3(cg)4(cg)4(ca)3 ] ; these motifs may adopt the left - handed z - dna conformation ( 3,54 ) .
g - quadruplex - forming repeat : four blocks , each containing the same number ( n ) of g bases ( n can vary from 3 to 7 ) , on the plus or minus strand , separated by 17  nt ; this type of dna sequence may adopt quadruplex structures ( 2 ) ; overlapping tracts of four g - blocks are also considered .
a - phased repeat : three runs of a bases ( a - tracts ) in phase with the helical pitch of the dna double - helix , i.e. 10  bp ; an a - tract is defined as a set of at base - pairs without a tpa step ( 47,5557 ) ; three or more tracts of a37 , t37 , aaattt , aaatttt and aaaattt ( in any combination ) on the plus or minus strand , whose centers are separated by 10 bases , are considered ; since a - tracts induce static bends in the dna double helix , the overall dna superhelix is expected to display either a left - handed or a right - handed writhe ( 47,5557 ) ; as mentioned , all the search criteria used herein do not allow for interruptions in the repeats and no thermodynamic information was factored - in in the algorithms used . the output file in gff format ( http://nonb.abcc.ncifcrf.gov/faqs/ )
the data from all such scans are merged and can be queried and displayed using our local instance of gbrowse called polybrowse ( 44 ) at http://pbrowse3.abcc.ncifcrf.gov/cgi-bin/gb2/gbrowse/human_37 and several gff - based query tools at http://nonb.abcc.ncifcrf.gov ( figure 1 ) .
importantly , the result pages produced from the queries contain links that allow the user to switch to the genome browser view of that feature , as well as a view that provides the sequence and other annotations for each feature . 

non - b db user interface ( a ) three query modes are available at the genomic database search tools page of the non - b db web interface : search by feature ( shown ) , search by feature attributes , and search by location ( a feature browser ) . in this example , all non - b dna motifs were queried by the gene symbol myc .
( b ) query results and links to polybrowse , ucsc genome browser , biodbnet ( 51 ) , etc . users may download the results in gff format and tab delimited format .
direct access to the sequence and other annotation information for each of the features is available by clicking on the 
non - b db user interface ( a ) three query modes are available at the genomic database search tools page of the non - b db web interface : search by feature ( shown ) , search by feature attributes , and search by location ( a feature browser ) . in this example , all non - b dna motifs were queried by the gene symbol myc .
( b ) query results and links to polybrowse , ucsc genome browser , biodbnet ( 51 ) , etc . users may download the results in gff format and tab delimited format .
direct access to the sequence and other annotation information for each of the features is available by clicking on the 
these data represent the basis for the non - b dna annotation information for each species .
the scanning criteria do not allow for mismatches within the repeat segments ; however , this feature may be added as information becomes available as to the acceptable structural tolerances for each mismatch case .
also , currently not included are very large palindromes ( > 100  kb ) , such as those that characterize the y chromosome and whose recombination is known to lead to spermatogenic failure ( 45,46 ) .
nevertheless , some aspects related to the presence of mismatches are presented in the polymorphism analysis described below . after scanning across different mammalian genomes ,
the numbers of each of the predicted classes of non - b dna structure - forming motifs appear to be quite variable ( table 2 ) . 

table 2.statistics for dna repeats and non - b dna - forming motifs in non - b dbdna featurehuman 37mouse 37dog 2chimp 2macaque 1g - quadruplex forming repeat374  545559  280492  535314  171298  142inverted repeat1  044  533801  242814  080998  249843  889cruciform motif197  910188 
532172  032190  736128  334direct repeat871  0451  593  107968  955787  335765  798slipped motif347  969695  150404  750314  516305  285mirror repeat16  51  7233  431  4861  829  86714  85  13514  55  025triplex motif1  79  623618  928336  6421 
05  6401  40  580z - dna repeat294  320690  276261  012278  928280  982a - phased repeat1 
130  7319  09  6531  241  0821  085  5911  098  030for the current releases of the five mammalian genomes indicated , the motif searches were performed and the number of features for each class was counted . according to table 1 , the cruciform motifs represent a subset of the inverted repeat class , the slipped motifs represent a subset of the direct repeat class and the triplex motif represents a subset of the mirror repeat class .
statistics for dna repeats and non - b dna - forming motifs in non - b db for the current releases of the five mammalian genomes indicated , the motif searches were performed and the number of features for each class was counted . according to table 1 , the cruciform motifs represent a subset of the inverted repeat class , the slipped motifs represent a subset of the direct repeat class and the triplex motif represents a subset of the mirror repeat class . as the overall base composition between different mammalian genomes is rather similar ( data not shown ) , the observed differences in the numbers of predicted non - b dna motifs could simply result from the altered arrangement of bases from one species to another .
alternatively , variations in the population of classes of repetitive elements ( sine , line , etc . ) among species , or other unknown features , might also contribute to the observed differences .
this interspecies variability appears to be uniformly distributed along the entire chromosomes , rather than concentrated in large repetitive clusters ( data not shown ) .
whether these differences play any role or contribute to conferring species - specific differences remains to be investigated .
a caveat concerning the simple assessment and comparison of the number of non - b dna - forming repeats among species relates to the criteria used and the counting method .
for example , in the g - quadruplex forming sequences , the pattern of a run of 3gs followed by 17 bases repeated four times can be extended , as long as more runs of gs are encountered , resulting in a single cluster that has the potential to form many substructures .
although our approach identifies this finding as a single cluster in the database , separate database tables are provided , in which all possible permutations of the sequence that satisfies the consensus sequence are reported .
in addition to the non - b dna predicted motifs , the database contains other features of the dna , such as phased a - tracts that impart static bends to the double - helix and may be involved in nucleosome assembly ( 47 ) , simple tandem repeats ( str ) including triplet repeats whose expansions cause a number of neuromuscular disorders ( 20 ) and poly(purinepyrimidine ) tracts , which are characterized by high stacking interactions ( 48 ) .
in addition , ncbi - derived features , such as genes , snps and repeatmasker ( http://www.repeatmasker.org/ ) elements are also included .
this integrated information is critical not only for guiding the user visually , but also for enabling queries that combine 
classes , such as  exons containing predicted  z - dna forming sequences , etc .
one of the main features of the non - b db is the ability to compare different mammalian genomes for the presence of non - b dna - forming motifs .
figure 2 illustrates this capability by comparing the presence of g - quadruplex forming motifs in the region upstream of the myc locus across the human , chimp , macaque , dog and mouse reference genomes . in order to view syntenic regions in other genomes , the liftover application from the ucsc website
areas where a syntenic match failed to be identified ( i.e. that region was absent in the other genome , or mapped redundantly ) do not show a link to that species .
other non - b dna tracks available in polybrowse will be described in more detail elsewhere ( cer et .
figure 2.(a ) polybrowse view of g - quadruplex forming repeats in the myc gene with species syntenic information .
visualization of cross - species information in polybrowse : the genomic region shown in figure 1 is displayed herein to illustrate the syntenic capability . in this case
, the g - quadruplex forming repeat track shows the locations of these features in the region near the beginning of the myc gene . the additional track illustrating the liftover1k blocks is also reported . in that track , when the user moves the cursor over the objects , a popup window appears that contains links ( a ) to the syntenic locations in the other available species .
( b  e ) results upon clicking on the links for chimp , macaque , mouse and dog , respectively .
currently , the liftover feature only allows moving from the human to the other species but not among the non - human species .
( a ) polybrowse view of g - quadruplex forming repeats in the myc gene with species syntenic information .
visualization of cross - species information in polybrowse : the genomic region shown in figure 1 is displayed herein to illustrate the syntenic capability . in this case
, the g - quadruplex forming repeat track shows the locations of these features in the region near the beginning of the myc gene .
when the user moves the cursor over the objects , a popup window appears that contains links ( a ) to the syntenic locations in the other available species .
( b  e ) results upon clicking on the links for chimp , macaque , mouse and dog , respectively .
currently , the liftover feature only allows moving from the human to the other species but not among the non - human species .
the computed non - b dna forming elements are likely to be under - represented in our reference genome as their underlying repeats may be polymorphic among individuals .
because this type of information may be critical in the context of gene regulation or predisposition to disease ( 48 ) , we used a specific parser to scan both the reference human genome as well as additional sequence sources , such as trace reads from the trace archive ( http://www.ncbi.nlm.nih.gov/traces/trace.cgi ) and contigs ( http://www.ncbi.nlm.nih.gov/projects/wgs/wgsprojectlist.cgi ) from personal genome projects ( 49 ) , for matches to the non - b dna motifs .
each match found in either the reference or alternate source is then scored for being polymorphic or not .
of the sites identified as polymorphic , a second evaluation is made to determine whether the polymorphism would affect the motif underlying the putative non - b dna structure .
the results of this scan are incorporated into the database as a series of separate tracks ( figure 3b , trace gplex tracks ) .
additional information can be gathered by extending this type of analysis to sequence alignments using closely related species .
currently , only the g - quadruplex forming motif supports this type of query . 

figure 3.(a ) polybrowse view of some of the non - b dna features in the human myc gene .
links to polybrowse and ucsc browsers : ( a ) display of polybrowse rendering the myc gene and its promoter with some of the non - b dna motifs tracks with the pupy [ poly(purinepyrimidine ) ] and polymorphism tracks turned on .
several other tracks include features computed at the advanced biomedical computing center ( abcc ) , such as strs , base composition , physical dna characteristics , mapping , as well as the ncbi derived features , including genes , snps , cytogenic markers assembly information , repeatmasker elements , etc .
( b ) display of the myc gene in polybrowse showing some of the polymorphism information in non - b db .
below the gene and mrna tracks , the searched motifs in the reference sequence are displayed ( teal ) .
an additional track shows specific g - quadruplex motifs in which all of the observed trace files contained mutations disrupting the g - quadruplex motif .
( c ) display of the myc gene at the ucsc human genome browser ( grch37/hg19 assembly ) as linked to non - b db from the search shown in figure 2b .
some of the non - b dna motifs from non - b db can be seen in the red rectangular box .
( a ) polybrowse view of some of the non - b dna features in the human myc gene .
links to polybrowse and ucsc browsers : ( a ) display of polybrowse rendering the myc gene and its promoter with some of the non - b dna motifs tracks with the pupy [ poly(purinepyrimidine ) ] and polymorphism tracks turned on .
several other tracks include features computed at the advanced biomedical computing center ( abcc ) , such as strs , base composition , physical dna characteristics , mapping , as well as the ncbi derived features , including genes , snps , cytogenic markers assembly information , repeatmasker elements , etc .
( b ) display of the myc gene in polybrowse showing some of the polymorphism information in non - b db .
below the gene and mrna tracks , the searched motifs in the reference sequence are displayed ( teal ) .
some of the predicted motifs are found only in the trace sequences . an additional track shows specific g - quadruplex motifs in which all of the observed trace files contained mutations disrupting the g - quadruplex motif .
( c ) display of the myc gene at the ucsc human genome browser ( grch37/hg19 assembly ) as linked to non - b db from the search shown in figure 2b .
some of the non - b dna motifs from non - b db can be seen in the red rectangular box . in order to provide access to the back - end database , we have leveraged two existing tools from the bioinformatics community .
the first is a bioperl ( 50 ) set of methods , which is used to query genome databases in various ways , such as by position , by class , or by attribute .
this same set of utilities is used within the context of polybrowse ( 44 ) , so that visualization of the genomic features is made available .
in addition to linking the outputs from the query tools to the browser for visualization , we also provide links allowing the returned data to be displayed in the familiar ucsc ( http://genome.ucsc.edu/ ) genome browser ( figure 3c ) , as well as links to our biodbnet database warehouse ( 51 ) , which contains gene - centric information derived from several sources , and additional links .
herein , we present a database containing the locations of motifs predicted to adopt the most common non - b dna structures . the database can be used to browse specific genomic regions for the possible contribution of non - b dna - forming elements to inherent biological observations derived from the region .
in addition to the locations of predicted motifs , the database also contains polymorphism information about each of the test sequences , as well as additional candidate sequences not present within the reference genomes .
additional genomes are in the process of being added to the system and will continue to be updated and added as they become available .
input from the community regarding the addition of other tracks , enhanced algorithms for the detection or scoring of the identified motifs or additional query tools are welcome and will be incorporated into the system as appropriate .
further additions , such as a community - based curation capability and the addition of other validation information through literature mining approaches are also under consideration .
we anticipate that significant improvements to our methods will be made in the future by incorporating energetic , and other secondary metrics , to the current predictive algorithms . although significant biological knowledge would be required , such as localized superhelical density , nucleosome positioning , etc .
( see above ) , the overall goal is to associate a likelihood index with each of the predicted locations for each of the non - b dna - forming classes .
finally , as reliable methods are expected to be developed that identify genome - wide data on non - b dna structures in vivo and some of the biological parameters involved , the resulting data sets can be used to train the prediction tools , resulting in improved predictive capabilities for each type of non - b - forming classes .
center for biomedical informatics and information technology ( cbiit)/cancer biomedical informatics grid ( cabig ) isrce yellow task # 09 - 260 to nci - frederick and national cancer institute / national institutes of health contract hhsn261200800001e ( to a.b . ) . funding for open access charge : national cancer institute / national institutes of health contract hhsn261200800001e .
 although the capability of dna to form a variety of non - canonical ( non - b ) structures has long been recognized , the overall significance of these alternate conformations in biology has only recently become accepted en masse . in order to provide access to genome - wide locations of these classes of predicted structures , 
 we have developed non - b db , a database integrating annotations and analysis of non - b dna - forming sequence motifs . 
 the database provides the most complete list of alternative dna structure predictions available , including z - dna motifs , quadruplex - forming motifs , inverted repeats , mirror repeats and direct repeats and their associated subsets of cruciforms , triplex and slipped structures , respectively . 
 the database also contains motifs predicted to form static dna bends , short tandem repeats and homo(purinepyrimidine ) tracts that have been associated with disease . 
 the database has been built using the latest releases of the human , chimp , dog , macaque and mouse genomes , so that the results can be compared directly with other data sources . in order to make the data interpretable in a genomic context , 
 features such as genes , single - nucleotide polymorphisms and repetitive elements ( sine , line , etc . ) have also been incorporated . 
 the database is accessed through query pages that produce results with links to the ucsc browser and a gbrowse - based genomic viewer . 
 it is freely accessible at http://nonb.abcc.ncifcrf.gov . 
INTRODUCTION
RESULTS
Non-B DB versus existing databases
Non-B DB: key features
Cross-species comparisons
Polymorphism analysis
CONCLUSIONS
FUNDING