146 research outputs found
GATA: a graphic alignment tool for comparative sequence analysis
BACKGROUND: Several problems exist with current methods used to align DNA sequences for comparative sequence analysis. Most dynamic programming algorithms assume that conserved sequence elements are collinear. This assumption appears valid when comparing orthologous protein coding sequences. Functional constraints on proteins provide strong selective pressure against sequence inversions, and minimize sequence duplications and feature shuffling. For non-coding sequences this collinearity assumption is often invalid. For example, enhancers contain clusters of transcription factor binding sites that change in number, orientation, and spacing during evolution yet the enhancer retains its activity. Dot plot analysis is often used to estimate non-coding sequence relatedness. Yet dot plots do not actually align sequences and thus cannot account well for base insertions or deletions. Moreover, they lack an adequate statistical framework for comparing sequence relatedness and are limited to pairwise comparisons. Lastly, dot plots and dynamic programming text outputs fail to provide an intuitive means for visualizing DNA alignments. RESULTS: To address some of these issues, we created a stand alone, platform independent, graphic alignment tool for comparative sequence analysis (GATA ). GATA uses the NCBI-BLASTN program and extensive post-processing to identify all small sub-alignments above a low cut-off score. These are graphed as two shaded boxes, one for each sequence, connected by a line using the coordinate system of their parent sequence. Shading and colour are used to indicate score and orientation. A variety of options exist for querying, modifying and retrieving conserved sequence elements. Extensive gene annotation can be added to both sequences using a standardized General Feature Format (GFF) file. CONCLUSIONS: GATA uses the NCBI-BLASTN program in conjunction with post-processing to exhaustively align two DNA sequences. It provides researchers with a fine-grained alignment and visualization tool aptly suited for non-coding, 0–200 kb, pairwise, sequence analysis. It functions independent of sequence feature ordering or orientation, and readily visualizes both large and small sequence inversions, duplications, and segment shuffling. Since the alignment is visual and does not contain gaps, gene annotation can be added to both sequences to create a thoroughly descriptive picture of DNA conservation that is well suited for comparative sequence analysis
Recommended from our members
Flexible Promoter Architecture Requirements for Coactivator Recruitment
Background: The spatial organization of transcription factor binding sites in regulatory DNA, and the composition of intersite sequences, influences the assembly of the multiprotein complexes that regulate RNA polymerase recruitment and thereby affects transcription. We have developed a genetic approach to investigate how reporter gene transcription is affected by varying the spacing between transcription factor binding sites. We characterized the components of promoter architecture that govern the yeast transcription factors Cbf1 and Met31/32, which bind independently, but collaboratively recruit the coactivator Met4. Results: A Cbf1 binding site was required upstream of a Met31/32 binding site for full reporter gene expression. Distance constraints on coactivator recruitment were more flexible than those for cooperatively binding transcription factors. Distances from 18 to 50 bp between binding sites support efficient recruitment of Met4, with only slight modulation by helical phasing. Intriguingly, we found that certain sequences located between the binding sites abolished gene expression. Conclusion: These results yield insight to the influence of both binding site architecture and local DNA flexibility on gene expression, and can be used to refine computational predictions of gene expression from promoter sequences. In addition, our approach can be applied to survey promoter architecture requirements for arbitrary combinations of transcription factor binding sites
Large-Scale Turnover of Functional Transcription Factor Binding Sites in Drosophila
The gain and loss of functional transcription factor binding sites has been proposed as a major source of evolutionary change in cis-regulatory DNA and gene expression. We have developed an evolutionary model to study binding-site turnover that uses multiple sequence alignments to assess the evolutionary constraint on individual binding sites, and to map gain and loss events along a phylogenetic tree. We apply this model to study the evolutionary dynamics of binding sites of the Drosophila melanogaster transcription factor Zeste, using genome-wide in vivo (ChIP–chip) binding data to identify functional Zeste binding sites, and the genome sequences of D. melanogaster, D. simulans, D. erecta, and D. yakuba to study their evolution. We estimate that more than 5% of functional Zeste binding sites in D. melanogaster were gained along the D. melanogaster lineage or lost along one of the other lineages. We find that Zeste-bound regions have a reduced rate of binding-site loss and an increased rate of binding-site gain relative to flanking sequences. Finally, we show that binding-site gains and losses are asymmetrically distributed with respect to D. melanogaster, consistent with lineage-specific acquisition and loss of Zeste-responsive regulatory elements
A Novel S -Adenosyl-l-methionine:Arsenic(III) Methyltransferase from Rat Liver Cytosol
S-Adenosyl-l-methionine (AdoMet):arsenic(III) methyltransferase, purified from liver cytosol of adult male Fischer 344 rats, catalyzes transfer of a methyl group from AdoMet to trivalent arsenicals producing methylated and dimethylated arsenicals. The kinetics of production of methylated arsenicals in reaction mixtures containing enzyme, AdoMet, dithiothreitol, glutathione (GSH), and arsenite are consistent with a scheme in which monomethylated arsenical produced from arsenite is the substrate for a second methylation reaction that yields dimethylated arsenical. The mRNA for this protein predicts a 369-amino acid residue protein (molecular mass 41056) that contains common methyltransferase sequence motifs. Its sequence is similar to Cyt19, a putative methyltransferase, expressed in human and mouse tissues. Reverse transcription-polymerase chain reaction detects S-adenosyl-l-methionine:arsenic(III) methyltransferase mRNA in rat tissues and in HepG2 cells, a human cell line that methylates arsenite and methylarsonous acid. S-Adenosyl-l-methionine:arsenic(III) methyltransferase mRNA is not detected in UROtsa cells, an immortalized human urothelial cell line that does not methylate arsenite. Because methylation of arsenic is a critical feature of its metabolism, characterization of this enzyme will improve our understanding of this metalloid's metabolism and its actions as a toxin and a carcinogen
Multi-Jet Event Rates in Deep Inelastic Scattering and Determination of the Strong Coupling Constant
Jet event rates in deep inelastic ep scattering at HERA are investigated
applying the modified JADE jet algorithm. The analysis uses data taken with the
H1 detector in 1994 and 1995. The data are corrected for detector and
hadronization effects and then compared with perturbative QCD predictions using
next-to-leading order calculations. The strong coupling constant alpha_S(M_Z^2)
is determined evaluating the jet event rates. Values of alpha_S(Q^2) are
extracted in four different bins of the negative squared momentum
transfer~\qq in the range from 40 GeV2 to 4000 GeV2. A combined fit of the
renormalization group equation to these several alpha_S(Q^2) values results in
alpha_S(M_Z^2) = 0.117+-0.003(stat)+0.009-0.013(syst)+0.006(jet algorithm).Comment: 17 pages, 4 figures, 3 tables, this version to appear in Eur. Phys.
J.; it replaces first posted hep-ex/9807019 which had incorrect figure 4
Multiplicity Structure of the Hadronic Final State in Diffractive Deep-Inelastic Scattering at HERA
The multiplicity structure of the hadronic system X produced in
deep-inelastic processes at HERA of the type ep -> eXY, where Y is a hadronic
system with mass M_Y< 1.6 GeV and where the squared momentum transfer at the pY
vertex, t, is limited to |t|<1 GeV^2, is studied as a function of the invariant
mass M_X of the system X. Results are presented on multiplicity distributions
and multiplicity moments, rapidity spectra and forward-backward correlations in
the centre-of-mass system of X. The data are compared to results in e+e-
annihilation, fixed-target lepton-nucleon collisions, hadro-produced
diffractive final states and to non-diffractive hadron-hadron collisions. The
comparison suggests a production mechanism of virtual photon dissociation which
involves a mixture of partonic states and a significant gluon content. The data
are well described by a model, based on a QCD-Regge analysis of the diffractive
structure function, which assumes a large hard gluonic component of the
colourless exchange at low Q^2. A model with soft colour interactions is also
successful.Comment: 22 pages, 4 figures, submitted to Eur. Phys. J., error in first
submission - omitted bibliograph
Differential (2+1) Jet Event Rates and Determination of alpha_s in Deep Inelastic Scattering at HERA
Events with a (2+1) jet topology in deep-inelastic scattering at HERA are
studied in the kinematic range 200 < Q^2< 10,000 GeV^2. The rate of (2+1) jet
events has been determined with the modified JADE jet algorithm as a function
of the jet resolution parameter and is compared with the predictions of Monte
Carlo models. In addition, the event rate is corrected for both hadronization
and detector effects and is compared with next-to-leading order QCD
calculations. A value of the strong coupling constant of alpha_s(M_Z^2)=
0.118+- 0.002 (stat.)^(+0.007)_(-0.008) (syst.)^(+0.007)_(-0.006) (theory) is
extracted. The systematic error includes uncertainties in the calorimeter
energy calibration, in the description of the data by current Monte Carlo
models, and in the knowledge of the parton densities. The theoretical error is
dominated by the renormalization scale ambiguity.Comment: 25 pages, 6 figures, 3 tables, submitted to Eur. Phys.
Measurements of Transverse Energy Flow in Deep-Inelastic Scattering at HERA
Measurements of transverse energy flow are presented for neutral current
deep-inelastic scattering events produced in positron-proton collisions at
HERA. The kinematic range covers squared momentum transfers Q^2 from 3.2 to
2,200 GeV^2, the Bjorken scaling variable x from 8.10^{-5} to 0.11 and the
hadronic mass W from 66 to 233 GeV. The transverse energy flow is measured in
the hadronic centre of mass frame and is studied as a function of Q^2, x, W and
pseudorapidity. A comparison is made with QCD based models. The behaviour of
the mean transverse energy in the central pseudorapidity region and an interval
corresponding to the photon fragmentation region are analysed as a function of
Q^2 and W.Comment: 26 pages, 8 figures, submitted to Eur. Phys.
Hadron Production in Diffractive Deep-Inelastic Scattering
Characteristics of hadron production in diffractive deep-inelastic
positron-proton scattering are studied using data collected in 1994 by the H1
experiment at HERA. The following distributions are measured in the
centre-of-mass frame of the photon dissociation system: the hadronic energy
flow, the Feynman-x (x_F) variable for charged particles, the squared
transverse momentum of charged particles (p_T^{*2}), and the mean p_T^{*2} as a
function of x_F. These distributions are compared with results in the gamma^* p
centre-of-mass frame from inclusive deep-inelastic scattering in the
fixed-target experiment EMC, and also with the predictions of several Monte
Carlo calculations. The data are consistent with a picture in which the
partonic structure of the diffractive exchange is dominated at low Q^2 by hard
gluons.Comment: 16 pages, 6 figures, submitted to Phys. Lett.
Measurement of D* Meson Cross Sections at HERA and Determination of the Gluon Density in the Proton using NLO QCD
With the H1 detector at the ep collider HERA, D* meson production cross
sections have been measured in deep inelastic scattering with four-momentum
transfers Q^2>2 GeV2 and in photoproduction at energies around W(gamma p)~ 88
GeV and 194 GeV. Next-to-Leading Order QCD calculations are found to describe
the differential cross sections within theoretical and experimental
uncertainties. Using these calculations, the NLO gluon momentum distribution in
the proton, x_g g(x_g), has been extracted in the momentum fraction range
7.5x10^{-4}< x_g <4x10^{-2} at average scales mu^2 =25 to 50 GeV2. The gluon
momentum fraction x_g has been obtained from the measured kinematics of the
scattered electron and the D* meson in the final state. The results compare
well with the gluon distribution obtained from the analysis of scaling
violations of the proton structure function F_2.Comment: 27 pages, 9 figures, 2 tables, submitted to Nucl. Phys.
- …