Variable-length Sequences Research Articles

Variational Upper and Lower Bounds for Probabilistic Graphical Models

Probabilistic phylogenetic models which relax the site independence evolution assumption often face the problem of infeasible likelihood computations, for example, for the task of selecting suitable parameters for the model. We present a new approximation method, applicable for a wide range of probabilistic models, which guarantees to upper and lower bound the true likelihood of data, and apply it to the problem of probabilistic phylogenetic models. The new method is complementary to known variational methods that lower bound the likelihood, and it uses similar methods to optimize the bounds from above and below. We applied our method to aligned DNA sequences of various lengths from human in the region of the CFTR gene and homologous from eight mammals, and found the bounds to be appreciably close to the true likelihood whenever it could be computed. When computing the exact likelihood was not feasible, we demonstrated the proximity of the upper and lower variational bounds, implying a tight approximation of the likelihood.

Elongin B/C Recruitment Regulates Substrate Binding by CIS

SOCS proteins play a major role in the regulation of cytokine signaling. They are recruited to activated receptors and can suppress signaling by different mechanisms including targeting of the receptor complex for proteasomal degradation. The activity of SOCS proteins is regulated at different levels including transcriptional control and posttranslational modification. We describe here a novel regulatory mechanism for CIS, one of the members of this protein family. A CIS mutant deficient in recruitment of the Elongin B/C complex completely failed to suppress STAT5 activation. This deficiency was not caused by altered turnover of CIS but by loss of cytokine receptor interaction. Intriguingly, no such effect was seen for binding to MyD88. The interaction between CIS and the Elongin B/C complex, which depends on the levels of uncomplexed Elongin B/C, was easily disrupted. This regulatory mechanism may be unique for CIS, as similar mutations in SOCS1, -2, -3, -6, and -7 had no functional impact. Our findings indicate that the SOCS box not only plays a role in the formation of E3 ligase complexes but, at least for CIS, can also regulate the binding modus of SOCS box-containing proteins.

Variable length tandem repeat polyglutamine sequences in the flexible tether region of the Tsr chemotaxis receptor of Escherichia coli

Methyl-accepting chemotaxis proteins (MCPs) are receptors that play an important role in bacterial chemotaxis. Methylation of Tsr, the MCP that mediates chemotaxis towards serine in Escherichia coli, is thought to be facilitated by binding of the methyltransferase to a flexible tether region at the C-terminal end of Tsr. This study analysed natural length variants of the tether that occur in E. coli due to genetic instability in tandem repeat DNA sequences that code for glutaminyl (Q) residues, creating polyQ sequences of variable lengths in the tether region. The tsr gene of E. coli K-12 (strain MG1655) codes for 4Q at the beginning of its 35 aa tether region. The tether varies in length from 35 to 47 residues among pathogenic and non-pathogenic strains of Escherichia, Shigella spp., Salmonella, Yersinia and Photorhabdus. Among previous sequences, Escherichia and Shigella mostly have 4Q and 7Q variants, and one strain (E. coli HS) has 10Q. In E. coli isolated from 50 humans and 75 animals (dogs, cats, horses, birds, etc.), polyQ up to 13Q (44 aa tether) were identified (6 strains); relative frequencies were 7Q ( approximately 77 % of the total) >4Q (14 %) >13Q (5 %) >10Q (4 %). Phylogenetic analysis revealed that E. coli strains with 10Q or 13Q largely fell within two clusters. Serine chemotaxis was not significantly different among 7Q, 10Q and 13Q strains, and was comparable to chemotaxis in the frequently studied K-12 strain. These results are consistent with models indicating that polyQ sequences from 7Q to 13Q are flexible, and that longer tethers, within this range, would not change the precision of adaptation mediated by methylation. Studies of this naturally variable polyQ region in E. coli may also have relevance to mechanisms that mediate polyQ instability in human genetic diseases.

Effective DNA inhibitors of cathepsin g by in vitro selection.

Cathepsin G (CatG) is a chymotrypsin-like protease released upon degranulation of neutrophils. In several inflammatory and ischaemic diseases the impaired balance between CatG and its physiological inhibitors leads to tissue destruction and platelet aggregation. Inhibitors of CatG are suitable for the treatment of inflammatory diseases and procoagulant conditions. DNA released upon the death of neutrophils at injury sites binds CatG. Moreover, short DNA fragments are more inhibitory than genomic DNA. Defibrotide, a single stranded polydeoxyribonucleotide with antithrombotic effect is also a potent CatG inhibitor. Given the above experimental evidences we employed a selection protocol to assess whether DNA inhibition of CatG may be ascribed to specific sequences present in defibrotide DNA. A Selex protocol was applied to identify the single-stranded DNA sequences exhibiting the highest affinity for CatG, the diversity of a combinatorial pool of oligodeoxyribonucleotides being a good representation of the complexity found in defibrotide. Biophysical and biochemical studies confirmed that the selected sequences bind tightly to the target enzyme and also efficiently inhibit its catalytic activity. Sequence analysis carried out to unveil a motif responsible for CatG recognition showed a recurrence of alternating TG repeats in the selected CatG binders, adopting an extended conformation that grants maximal interaction with the highly charged protein surface. This unprecedented finding is validated by our results showing high affinity and inhibition of CatG by specific DNA sequences of variable length designed to maximally reduce pairing/folding interactions.

A regulatory SNP of the BICD1 gene contributes to telomere length variation in humans

Telomeres are repetitive sequences of variable length at the ends of chromosomes involved in maintaining their integrity. Telomere dysfunction is associated with increased risk of cancer and other age-related diseases. Telomere length is an important determinant of telomere function and has a strong genetic basis. We previously carried out a genome-wide linkage analysis of mean leukocyte telomere length, and identified a 12 cm quantitative-trait locus affecting telomere length on human chromosome 12. In the present study we confirmed linkage to this locus in an extended sample (380 families, 520 sib-pairs, maximum LOD score 4.3). Fine-mapping identified a 51 kb region of association within intron 1 of the Bicaudal-D homolog 1 (BICD1, MIM 602204) gene. The strongest association (P = 1.9 x 10(-5)) was with SNP rs2630578 where the minor allele C (frequency 0.21) was associated with telomeres that were shorter by 604 (+/-204) base pairs, equivalent to approximately 15-20 years of age-related attrition in telomere length. Subjects carrying the C allele for rs2630778 had 44% lower BICD1 mRNA levels in their leukocytes compared with GG homozygotes (P = 0.004). BICD1 is involved in Golgi-to-endoplasmic reticulum vacuolar transport. Previous studies have implicated vacuolar genes in telomere length homeostasis in yeast. Our study indicates that BICD1 plays a similar role in humans.

Outbreaks and risks of infectious spleen and kidney necrosis virus disease in freshwater ornamental fishes

We examined the distribution of iridoviruses in 10 freshwater ornamental fish species hatched in Korea and imported from other Asian countries using both 1-step and 2-step polymerase chain reation (PCR). None of the 10 fish species analyzed were free of iridovirus as shown by 2-step PCR positive results, and 3 species yielded 1-step PCR positive results with associated mortality. Cloned PCR amplicons of the adenosine triphosphatase (ATPase) and major capsid protein (MCP) genes in genomic DNA of iridovirus showed the same nucleotide sequences as that of infectious spleen and kidney necrosis virus (ISKNV) isolated from the mandarinfish Siniperca chuatsi. These results indicate the presence of ISKNV disease in various ornamental fish as new host species and that the disease is widespread throughout different Asian countries including Korea, Singapore and China. Such infections were either clinical with associated mortality (and 1-step PCR positive) or asymptomatic in fish that were externally healthy (and only positive in 2-step PCR). Molecular analyses of the K2 region performed on iridovirus samples isolated from freshwater ornamental fishes revealed deletion/insertion of repetitive sequences of various lengths (42 to 339 bp), depending on the ISKNV isolates, without substitutions. Experimental infection of pearl gourami Trichogaster leeri and silver gourami T. microlepis with a tissue homogenate of pearl gourami infected by ISKNV induced 70 and 20% cumulative mortalities in the pearl and silver gourami, respectively.

Rate-Distortion Optimized Video Transmission Over DS-CDMA Channels with Auxiliary Vector Filter Single-User Multirate Detection

In this paper, we consider the rate-distortion optimized resource allocation for video transmission over multi-rate wireless direct-sequence code-division-multiple-access (DS-CDMA) channels. We consider the performance of transmitting scalable video over a multipath Rayleigh fading channel via a combination of multi-code multirate CDMA and variable sequence length multirate CDMA channel system. At the receiver, despreading is done using adaptive space-time auxiliary-vector (AV) filters. We propose a new interference cancelling design that uses just a single AV filter for single-user mutirate despreading. Our experimental results show that the proposed interference cancelling design has excellent performance in scalable video transmission over DS-CDMA systems that use a combination of multicode multirate and variable processing gain multirate CDMA. The proposed design takes advantage of the fact that single user's video data is transmitted using two spreading codes, one for the base layer and one for the enhancement layers, and of the fact that these spreading codes can have different processing gains. The proposed interference cancelling design is compared with two conventional single-user multirate CDMA receiver configurations, however now we use an AV filter rather than a simple matched filter. We also propose a resource allocation algorithm for the optimal determination of source coding rate, channel coding rate and processing gain for each scalable layer, in order to minimize the expected distortion at the receiver.

Acoustic Modelling Using Continuous Rational Kernels

Many discriminative classification algorithms are designed for tasks where samples can be represented by fixed-length vectors. However, many examples in the fields of text processing, computational biology and speech recognition are best represented as variable-length sequences of vectors. Although several dynamic kernels have been proposed for mapping sequences of discrete observations into fixed-dimensional feature-spaces, few kernels exist for sequences of continuous observations. This paper introduces continuous rational kernels, an extension of standard rational kernels, as a general framework for classifying sequences of continuous observations. In addition to allowing new task-dependent kernels to be defined, continuous rational kernels allow existing continuous dynamic kernels, such as Fisher and generative kernels, to be calculated using standard weighted finite-state transducer algorithms. Preliminary results on both a large vocabulary continuous speech recognition (LVCSR) task and the TIMIT database are presented.

Autonomous folding of interdomain regions of a modular polyketide synthase

Domains within the multienzyme polyketide synthases are linked by noncatalytic sequences of variable length and unknown function. Recently, the crystal structure was reported of a portion of the linker between the acyltransferase (AT) and ketoreductase (KR) domains from module 1 of the erythromycin synthase (6-deoxyerythronolide B synthase), as a pseudodimer with the adjacent ketoreductase (KR). On the basis of this structure, the homologous linker region between the dehydratase (DH) and enoyl reductase (ER) domains in fully reducing modules has been proposed to occupy a position on the periphery of the polyketide synthases complex, as in porcine fatty acid synthase. We report here the expression and characterization of the same region of the 6-deoxyerythronolide B synthase module 1 AT-KR linker, without the adjacent KR domain (termed DeltaN AT1-KR1), as well as the corresponding section of the DH-ER linker. The linkers fold autonomously and are well structured. However, analytical gel filtration and ultracentrifugation analysis independently show that DeltaN AT1-KR1 is homodimeric in solution; site-directed mutagenesis further demonstrates that linker self-association is compatible with the formation of a linker-KR pseudodimer. Our data also strongly indicate that the DH-ER linker associates with the upstream DH domain. Both of these findings are incompatible with the proposed model for polyketide synthase architecture, suggesting that it is premature to allocate the linker regions to a position in the multienzymes based on the solved structure of animal fatty acid synthase.

Alignment and Topological Accuracy of the Direct Optimization approach via POY and Traditional Phylogenetics via ClustalW + PAUP*

Direct optimization frameworks for simultaneously estimating alignments and phylogenies have recently been developed. One such method, implemented in the program POY, is becoming more common for analyses of variable length sequences (e.g., analyses using ribosomal genes) and for combined evidence analyses (morphology + multiple genes). Simulation of sequences containing insertion and deletion events was performed in order to directly compare a widely used method of multiple sequence alignment (ClustalW) and subsequent parsimony analysis in PAUP* with direct optimization via POY. Data sets were simulated for pectinate, balanced, and random tree shapes under different conditions (clocklike, non-clocklike, and ultrametric). Alignment accuracy scores for the implied alignments from POY and the multiple sequence alignments from ClustalW were calculated and compared. In almost all cases (99.95%), ClustalW produced more accurate alignments than POY-implied alignments, judged by the proportion of correctly identified homologous sites. Topological accuracy (distance to the true tree) for POY topologies and topologies generated under parsimony in PAUP* from the ClustalW alignments were also compared. In 44.94% of the cases, Clustal alignment tree reconstructions via PAUP* were more accurate than POY, whereas in 16.71% of the cases POY reconstructions were more topologically accurate (38.38% of the time they were equally accurate). Comparisons between POY hypothesized alignments and the true alignments indicated that, on average, as alignment error increased, topological accuracy decreased.

Low Energy Peptide Fragmentations in an ESI-Q-Tof Type Mass Spectrometer

Efficient peptide sequencing relies on both high quality MS/MS data acquisition and exhaustive knowledge of gas-phase dissociation mechanisms. We report our contribution to the elaboration of more comprehensive fragmentation models required for efficient automated MS/MS spectra interpretation. Following a statistical approach, various peptides (296 sequences of variable compositions and lengths) were prepared and subjected to low-energy collision-induced dissociations (CID) in an electrospray hybrid instrument (ESI-Q-q-Tof type mass spectrometer) that has retained relatively limited attention so far. Besides, our studies were focused on low molecular weight singly charged peptides that often failed to be identified by sequencing algorithms. Only half of the studied compounds showed charge directed dissociations in accordance with the mobile proton model producing fragment ions directly related to the primary sequence. For the peptides that did not exhibit the expected fragment ion series, alternative dissociation behaviors issued from complex rearrangements were evidenced.

Performance of impulse radio direct sequence ultra-wideband system with variable-length spreading sequences

An adaptive transmission scheme using variable-length spreading sequences (VLSS) is proposed. The system is based on an impulse direct sequence ultra-wideband (DS-UWB) approach for wireless communication where the number of active users is variable. In contrast to the conventional DS-UWB system, which employs a fixed length spreading code, the proposed system changes the length of spreading code according to the system load adaptively, which is proven to be able to reduce the inter-chip interference, inter-symbol interference and multiple-access interference and thus improve the system performance in terms of bit error rate. The proposed system is simulated in multipath channels with a high level of multipath fragments, modelled by the standard IEEE 802.15.3a indoor channel model 1 (CM1) and channel model 2 (CM2). Numerical results demonstrate that when RAKE receivers are employed, the proposed VLSS DS-UWB system outperforms the conventional system by appropriately allocating the spreading sequences.

Cloning, characterization and promoter analysis of common carphairy/Enhancer-of-split-related gene,her6

Some members of hairy/Enhancer-of-split-related gene (HES) family have important effects on axial mesoderm segmentation and the establishment and maintenance of the somite fringe. In fishes, the her6 gene, a member of the HES family, is the homologue of hes1 in mammals and chicken. In this study, the her6 gene and its full-length cDNA from the common carp (Cyprinus carpio) were isolated and characterized. The genomic sequence of common carp her6 is approximately 1.7 kb, with four exons and three introns, and the full-length cDNA of 1314 bp encodes a putative polypeptide of 271 amino acids. To analyse the promoter sequence of common carp her6, sequences of various lengths upstream from the transcription initiation site of her6 were fused to enhanced green fluorescent protein gene (eGFP) and introduced into zebrafish embryos by microinjection to generate transgenic embryos. Our results show that the upstream sequence of 500 bp can direct highly efficient and tissue-specific expression of eGFP in zebrafish embryos, whereas a fragment of 200 bp containing the TATA box and a partial suppressor of hairless paired site sequence (SPS) is not sufficient to drive eGFP expression in zebrafish embryos.

Tandem mass spectrometry of amidated peptides

The behavior of C-terminal amidated and carboxylated peptides upon low-energy collision-induced dissociation (CID) was investigated. Two sets of 76 sequences of variable amino acid compositions and lengths were synthesized as model compounds. In most cases, C-terminal amidated peptides were found to produce, upon CID, an abundant loss of ammonia from the protonated molecules. To validate such MS/MS signatures, the studied peptides contained amino acids that can potentially release ammonia from their side chains, such as asparagine, glutamine, tryptophan, lysine and arginine. Arginine, and to a lesser extent lysine, was shown to induce a competitive fragmentation leading to the loss of ammonia from their side chains, thus interfering with the targeted backbone neutral release. However, when arginine or lysine was located at the C-terminal position mimicking a tryptic digest, losses of ammonia from the arginine side chain and from the peptide backbone were completely suppressed. Such results were discussed in the frame of peptidomic or proteomic studies in an attempt to reveal the presence of C-terminal amidated peptides or proteins.

Species-Specific Identification of Leptospiraceae by 16S rRNA Gene Sequencing

The genus Leptospira is classified into 13 named species and 4 genomospecies based upon DNA-DNA reassociation studies. Phenotypic tests are unable to distinguish between species of Leptospira, and there is a need for a simplified molecular approach to the identification of leptospires. 16S rRNA gene sequences are potentially useful for species identification of Leptospira, but there are a large number of sequences of various lengths and quality in the public databases. 16S rRNA gene sequences of near full length and bidirectional high redundancy were determined for all type strains of the species of the Leptospiraceae. Three clades were identified within the genus Leptospira, composed of pathogenic species, nonpathogenic species, and another clade of undetermined pathogenicity with intermediate 16S rRNA gene sequence relatedness. All type strains could be identified by 16S rRNA gene sequences, but within both pathogenic and nonpathogenic clades as few as two or three base pairs separated some species. Sequences within the nonpathogenic clade were more similar, and in most cases < or =10 bp distinguished these species. These sequences provide a reference standard for identification of Leptospira species and confirm previously established relationships within the genus. 16S rRNA gene sequencing is a powerful method for identification in the clinical laboratory and offers a simplified approach to the identification of Leptospira species.

A Glycine- and Glutamate-Rich Protein Is Female Salivary Gland-Specific and Abundant in the Malaria Vector Anopheles dirus B (Diptera: Culicidae)

Before transmission, malaria parasites reside in the salivary glands of their female mosquito hosts. Saliva proteins assist in blood feeding and also may influence the ability of mosquitoes to transmit malaria. We attempted to identify and isolate cDNAs encoding proteins expressed at a high level in the salivary glands of a malaria vector, Anopheles dirus B Peyton and Harrison (= An. cracens) (Diptera: Culicidae). A major protein with an estimated molecular mass of 35 kDa and an isoelectric point (pI) of approximately 4 was detected on a two-dimensional (2D) gel. Internal peptide sequences of the protein spot showed high similarity to sequences present in the conserved C-terminal domain of glycine- and glutamate (GE)-rich proteins. A full-length cDNA encoding this protein was isolated from a salivary gland cDNA library of female An. dirus B. The cDNA encoded a 256-residue protein with a calculated molecular mass of 25.4 kDa and a pI of 3.9. BLAST analysis confirmed that it is a member of the GE-rich family. Compositional and sequence analysis of this and other family members revealed a highly acidic N-terminal region of variable length and low sequence conservation and a well conserved C-terminal domain containing 10 identical residues across the 13 known members of the gene family in mosquitoes. The An. dirus B GE-rich transcript was detected by reverse transcription-polymerase chain reaction (PCR) only in the female salivary glands, indicating that this protein is female saliva-specific. The GE-rich proteins may function as a salivary lubricant to facilitate blood feeding.

Dynamic homology and the likelihood criterion.

The use of likelihood as an optimality criterion is explored in the context of dynamic homology. Simple models and procedures are described to allow the analysis of large variable length sequence data sets, alone and in combination with qualitative information (such as morphology). Several approaches are discussed that have different likelihood interpretations in terms of maximum parsimony likelihood and maximum average likelihood. Implementation is discussed and an example in arthropod systematics presented. Topological congruence comparisons with parsimony are made.

Generation and upper bound of one-coincidence sequences in frequency-hopping multiple-access systems

In frequency-hopping one-coincidence sequences, a specified distance between two adjacent symbols is usually employed to reduce the interference among sequences. We propose a new algorithm to construct one-coincidence sequence sets satisfying distance constraints. Compared to the conventional method, our approach provides a more flexible and general way to generate one-coincidence sequences of variable lengths. We also obtain a tighter upper bound for one-coincidence sequences of odd length, which is independent of code lengths.

Minding the gap: frequency of indels in mtDNA control region sequence data and influence on population genetic analyses

Insertions and deletions (indels) result in sequences of various lengths when homologous gene regions are compared among individuals or species. Although indels are typically phylogenetically informative, occurrence and incorporation of these characters as gaps in intraspecific population genetic data sets are rarely discussed. Moreover, the impact of gaps on estimates of fixation indices, such as F(ST), has not been reviewed. Here, I summarize the occurrence and population genetic signal of indels among 60 published studies that involved alignments of multiple sequences from the mitochondrial DNA (mtDNA) control region of vertebrate taxa. Among 30 studies observing indels, an average of 12% of both variable and parsimony-informative sites were composed of these sites. There was no consistent trend between levels of population differentiation and the number of gap characters in a data block. Across all studies, the average influence on estimates of PhiST was small, explaining only an additional 1.8% of among population variance (range 0.0-8.0%). Studies most likely to observe an increase in PhiST with the inclusion of gap characters were those with < 20 variable sites, but a near equal number of studies with few variable sites did not show an increase. In contrast to studies at interspecific levels, the influence of indels for intraspecific population genetic analyses of control region DNA appears small, dependent upon total number of variable sites in the data block, and related to species-specific characteristics and the spatial distribution of mtDNA lineages that contain indels.

Crystal Structure of the First KH Domain of Human Poly(C)-binding Protein-2 in Complex with a C-rich Strand of Human Telomeric DNA at 1.7 Å

Recognition of poly(C) DNA and RNA sequences in mammalian cells is achieved by a subfamily of the KH (hnRNP K homology) domain-containing proteins known as poly(C)-binding proteins (PCBPs). To reveal the molecular basis of poly(C) sequence recognition, we have determined the crystal structure, at 1.7-A resolution, of PCBP2 KH1 in complex with a 7-nucleotide DNA sequence (5'-AACCCTA-3') corresponding to one repeat of the human C-rich strand telomeric DNA. The protein-DNA interaction is mediated by the combination of several stabilizing forces including hydrogen bonding, electrostatic interactions, van der Waals contacts, and shape complementarities. Specific recognition of the three cytosine residues is realized by a dense network of hydrogen bonds involving the side chains of two conserved lysines and one glutamic acid. The co-crystal structure also reveals a protein-protein dimerization interface of PCBP2 KH1 located on the opposite side of the protein from the DNA binding groove. Numerous stabilizing protein-protein interactions, including hydrophobic contacts, stacking of aromatic side chains, and a large number of hydrogen bonds, indicate that the protein-protein interaction interface is most likely genuine. Interaction of PCBP2 KH1 with the C-rich strand of human telomeric DNA suggests that PCBPs may participate in mechanisms involved in the regulation of telomere/telomerase functions.