Application of discrete Fourier inter-coefficient difference for assessing genetic sequence similarity

Digital signal processing (DSP) techniques for biological sequence analysis continue to grow in popularity due to the inherent digital nature of these sequences. DSP methods have demonstrated early success for detection of coding regions in a gene. Recently, these methods are being used to establish DNA gene similarity. We present the inter-coefficient difference (ICD) transformation, a novel extension of the discrete Fourier transformation, which can be applied to any DNA sequence. The ICD method is a mathematical, alignment-free DNA comparison method that generates a genetic signature for any DNA sequence that is used to generate relative measures of similarity among DNA sequences. We demonstrate our method on a set of insulin genes obtained from an evolutionarily wide range of species, and on a set of avian influenza viral sequences, which represents a set of highly similar sequences. We compare phylogenetic trees generated using our technique against trees generated using traditional alignment techniques for similarity and demonstrate that the ICD method produces a highly accurate tree without requiring an alignment prior to establishing sequence similarity.     

The Molecular Clock Runs at Different Rates Among Closely Related Members of a Gene Family

The serum albumin gene family is composed of four members that have arisen by a series of duplications from a common ancestor. From sequence differences between members of the gene family, we infer that a gene duplication some 580 Myr ago gave rise to the vitamin D–binding protein (DBP) gene and a second lineage, which reduplicated about 295 Myr ago to give the albumin (ALB) gene and a common precursor to α-fetoprotein (AFP) and α-albumin (ALF). This precursor itself duplicated about 250 Myr ago, giving rise to the youngest family members, AFP and ALF. It should be possible to correlate these dates with the phylogenetic distribution of members of the gene family among different species. All four genes are found in mammals, but AFP and ALF are not found in amphibia, which diverged from reptiles about 360 Myr ago, before the divergence of the AFP-ALF progenitor from albumin. Although individual family members display an approximate clock-like evolution, there are significant deviations—the rates of divergence for AFP differ by a factor of 7, the rates for ALB differ by a factor of 2.1. Since the progenitor of this gene family itself arose by triplication of a smaller gene, the rates of evolution of individual domains were also calculated and were shown to vary within and between family members. The great variation in the rates of the molecular clock raises questions concerning whether it can be used to infer evolutionary time from contemporary sequence differences. Received: 28 February 1995 / Accepted: 6 October 1997     

Isolation and analysis of a cDNA clone that encodes an alfalfa (Medicago sativa) aspartate aminotransferase

Summary Nitrate reductase (NR) assays revealed a bi-specific NAD(P)H-NR (EC 1.6.6.2.) to be the only nitrate-reducing enzyme in leaves of hydroponically grown birches. To obtain the primary structure of the NAD(P)H-NR, leaf poly(A)⁺ mRNA was used to construct a cDNA library in the lambda gt11 phage. Recombinant clones were screened with heterologous gene probes encoding NADH-NR from tobacco and squash. A 3.0 kb cDNA was isolated which hybridized to a 3.2 kb mRNA whose level was significantly higher in plants grown on nitrate than in those grown on ammonia. The nucleotide sequence of the cDNA comprises a reading frame encoding a protein of 898 amino acids which reveals 67%–77% identity with NADH-nitrate reductase sequences from higher plants. To identify conserved and variable regions of the multicentre electron-transfer protein a graphical evaluation of identities found in NR sequence alignments was carried out. Thirteen well-conserved sections exceeding a size of 10 amino acids were found in higher plant nitrate reductases. Sequence comparisons with related redox proteins indicate that about half of the conserved NR regions are involved in cofactor binding. The most striking difference in the birch NAD(P)H-NR sequence in comparison to NADH-NR sequences was found at the putative pyridine nucleotide binding site. Southern analysis indicates that the bi-specific NR is encoded by a single copy gene in birch. These sequence data appeared in the EMBL/GenBank/DDBJ nucleotide sequence data bases under the accession number X54097     

Sequence instability in the long terminal repeats of avian spleen necrosis virus and reticuloendotheliosis virus

Summary Sequence divergence between the 3 long terminal repeats (LTR) of avian reticuloendotheliosis virus (REV), deletion variant proviral clone 2-20-4, and spleen necrosis virus (SNV)—proviral clones 14-44, 60, and 70—was found to involve two classes of base substitutions: low-frequency interspersed and high-frequency clustered substitutions. Clones 2-20-4 and 14-44 have diverged 4.4% owing to low-frequency substitutions. In contrast, two high-frequency substitution segments have diverged by 30% and 29%, respectively. Clustered substitutions appear to be located either within or next to tandem repeats, suggesting their introduction concomitant with sequence deletions and duplications commonly associated with such repeats. A new 19-bp tandem repeat is found in clone 2-20-4. Its sequence could have evolved from the 26-bp repeats found in the SNV clones.     

HDAC inhibition upregulates the expression of angiostatic ADAMTS1

HDAC inhibitors are promising anticancer agents that induce cell cycle arrest and apoptosis. However, the role of HDACs in cancer progression, such as angiogenesis and metastasis, remains largely unexplored. Among various HDAC inhibitors, we demonstrate that TSA and SAHA upregulated the expression of angiostatic ADAMTS1 in A549 cells. HDAC6 inhibitor tubacin, and knockdown of HDAC6, also lead to ADAMTS1 upregulation. By reporter, DAPA, and ChIP assays, the proximal GC boxes were demonstrated to be essential for ADAMTS1 induction. Decreased binding of SP1 and HDAC6 to the ADAMTS1 promoter after TSA treatment was also seen. These data suggest the involvement of HDAC6 and SP1 in the HDACi-induced expression of angiostatic ADAMTS1.     

Expression and Regulation of the Arabidopsis thaliana Cel1 Endo 1,4 �� Glucanase Gene During Compatible Plant-Nematode Interactions

The root-knot nematode Meloidogyne incognita is an obligate endoparasite of plant roots and stimulates elaborate modifications of selected root vascular cells to form giant cells for feeding. An Arabidopsis thaliana endoglucanase (Atcel1) promoter is activated in giant cells that were formed in Atcel1::UidA transgenic tobacco and Arabidopsis plants. Activity of the full-length Atcel1 promoter was detected in root and shoot elongation zones and in the lateral root primordia. Different 5’ and internal deletions of regions of the 1,673 bp Atcel1 promoter were each fused to the UidA reporter gene and transformed in tobacco, and roots of the transformants were inoculated with M. incognita to assay for GUS expression in giant cells and noninfected plant tissues. Comparison of the Atcel1 promoter deletion constructs showed that the region between −1,673 and −1,171 (fragment 1) was essential for Atcel1 promoter activity in giant cells and roots. Fragment 1 alone, however, was not sufficient for Atcel1 expression in giant cells or roots, suggesting that cis-acting elements in fragment 1 may function in consort with other elements within the Atcel1 promoter. Root-knot nematodes and giant cells developed normally within roots of Arabidopsis that expressed a functional antisense construct to Atcel1, suggesting that a functional redundancy in endoglucanase activity may represent another level of regulatory control of cell wall-modifying activity within nematode feeding cells.     

Molecular cloning of Plasmodium vivax calcium-dependent protein kinase 4

A family of calcium-dependent protein kinases (CDPKs) is a unique enzyme which plays crucial roles in intracellular calcium signaling in plants, algae, and protozoa. CDPKs of malaria parasites are known to be key regulators for stage-specific cellular responses to calcium, a widespread secondary messenger that controls the progression of the parasite. In our study, we identified a gene encoding Plasmodium vivax CDPK4 (PvCDPK4) and characterized its molecular property and cellular localization. PvCDPK4 was a typical CDPK which had well-conserved N-terminal kinase domain and C-terminal calmodulin-like structure with 4 EF hand motifs for calcium-binding. The recombinant protein of EF hand domain of PvCDPK4 was expressed in E. coli and a 34 kDa product was obtained. Immunofluorescence assay by confocal laser microscopy revealed that the protein was expressed at the mature schizont of P. vivax. The expression of PvCDPK4-EF in schizont suggests that it may participate in the proliferation or egress process in the life cycle of this parasite.     

OPUS-Dom: Applying the Folding-Based Method VECFOLD to Determine Protein Domain Boundaries

In this article, we present a de novo method for predicting protein domain boundaries, called OPUS-Dom. The core of the method is a novel coarse-grained folding method, VECFOLD, which constructs low-resolution structural models from a target sequence by folding a chain of vectors representing the predicted secondary-structure elements. OPUS-Dom generates a large ensemble of folded structure decoys by VECFOLD and labels the domain boundaries of each decoy by a domain parsing algorithm. Consensus domain boundaries are then derived from the statistical distribution of the putative boundaries and three empirical sequence-based domain profiles. OPUS-Dom generally outperformed several state-of-the-art domain prediction algorithms over various benchmark protein sets. Even though each VECFOLD-generated structure contains large errors, collectively these structures provide a more robust delineation of domain boundaries. The success of OPUS-Dom suggests that the arrangement of protein domains is more a consequence of limited coordination patterns per domain arising from tertiary packing of secondary-structure segments, rather than sequence-specific constraints.