An algorithm to predict 3 ' intron splice sites in Plasmodium falciparum genomic sequences

A new algorithm, PfAGSS, for predicting 3' splice sites in Plasmodium falciparum genomic sequences is described. Application of this program to the published P. falciparum chromosome 2 and 3 data suggests that existing programs result in a high error rate in assigning 3' intron boundaries. (C) 2001 Elsevier Science B.V. All rights reserved.

A simulated annealing algorithm for finding consensus sequences

Motivation: A consensus sequence for a family of related sequences is, as the name suggests, a sequence that captures the features common to most members of the family. Consensus sequences are important in various DNA sequencing applications and are a convenient way to characterize a family of molecules. Results: This paper describes a new algorithm for finding a consensus sequence, using the popular optimization method known as simulated annealing. Unlike the conventional approach of finding a consensus sequence by first forming a multiple sequence alignment, this algorithm searches for a sequence that minimises the sum of pairwise distances to each of the input sequences. The resulting consensus sequence can then be used to induce a multiple sequence alignment. The time required by the algorithm scales linearly with the number of input sequences and quadratically with the length of the consensus sequence. We present results demonstrating the high quality of the consensus sequences and alignments produced by the new algorithm. For comparison, we also present similar results obtained using ClustalW. The new algorithm outperforms ClustalW in many cases.

Molecular phylogenetics of Lentibulariaceae inferred from plastid rps16 Intron and trnL-F DNA sequences: implications for character evolution and biogeography

Phylogenetic relationships among 75 species of Lentibulariaceae, representing the three recognized genera, were assessed by cladistic analysis of DNA sequences from the plastid rps16 intron and the trnL-F region. Sequence data from the two loci were analyzed both separately and in combination. Consensus trees from all analyses are congruent, and parsimony jackknife results demonstrate strong support for relationships both between and within each of the three demonstrably monophyletic genera. The genus Pinguicula is sister to a Genlisea-Utricularia clade, the phylogenetic structure within this clade closely follows Taylor's recent sectional delimitations based on morphology. Three principal clades are shown within Utricularia, with the basal sections Polypoinpholyx and Pleiochasia together forming the sister lineage of the remaining Utricularia species. Of the fundamental morphological specializations, the stoloniferous growth form apparently arose independently within Genlisea and Utricularia three times, and within Utricularia itself, perhaps more than once. The epiphytic habit has evolved independently at least three times, in Pinguicula, in Utricularia section Phyllaria, and within the two sections Orchidioides and Iperua (in the latter as bromeliad tank-epiphytes). The suspended aquatic habit may have evolved independently within sections Utricularia and Vesiculina. Biogeographic optimization on the phylogeny demonstrates patterns commonly associated with the boreotropics hypothesis and limits the spatial origin of Lentibulariaceae to temperate Eurasia or tropical America.

Inferring genome trees by using a filter to eliminate phylogenetically discordant sequences and a distance matrix based on mean normalized BLASTP scores

Darwin's paradigm holds that the diversity of present-day organisms has arisen via a process of genetic descent with modification, as on a bifurcating tree. Evidence is accumulating that genes are sometimes transferred not along lineages but rather across lineages. To the extent that this is so, Darwin's paradigm can apply only imperfectly to genomes, potentially complicating or perhaps undermining attempts to reconstruct historical relationships among genomes (i.e., a genome tree). Whether most genes in a genome have arisen via treelike (vertical) descent or by lateral transfer across lineages can be tested if enough complete genome sequences are used. We define a phylogenetically discordant sequence (PDS) as an open reading frame (ORF) that exhibits patterns of similarity relationships statistically distinguishable from those of most other ORFs in the same genome. PDSs represent between 6.0 and 16.8% (mean, 10.8%) of the analyzable ORFs in the genomes of 28 bacteria, eight archaea, and one eukaryote (Saccharomyces cerevisiae). In this study we developed and assessed a distance-based approach, based on mean pairwise sequence similarity, for generating genome trees. Exclusion of PDSs improved bootstrap support for basal nodes but altered few topological features, indicating that there is little systematic bias among PDSs. Many but not all features of the genome tree from which PDSs were excluded are consistent with the 16S rRNA tree.

C-terminal sequences in R-Ras are involved in integrin regulation and in plasma membrane microdomain distribution

The small GTPases R-Ras and H-Ras are highly homologous proteins with contrasting biological properties, for example, they differentially modulate integrin affinity: H-Ras suppresses integrin activation in fibroblasts whereas R-Ras can reverse this effect of H-Ras. To gain insight into the sequences directing this divergent phenotype, we investigated a panel of H-Ras/R-Ras chimeras and found that sequences in the R-Ras hypervariable C-terminal region including amino acids 175-203 are required for the R-Ras ability to increase integrin activation in CHO cells; however, the proline-rich site in this region, previously reported to bind the adaptor protein Nck, was not essential for this effect. In addition, we found that the GTPase TC21 behaved similarly to R-Ras. Because the C-termini of Ras proteins can control their subcellular localization, we compared the localization of H-Ras and R-Ras. In contrast to H-Ras, which migrates out of lipid rafts upon activation, we found that activated R-Ras remained localized to lipid rafts. However, functionally distinct H-Ras/R-Ras chimeras containing different C-terminal R-Ras segments localized to lipid rafts irrespective of their integrin phenotype. (C) 2003 Elsevier Inc. All rights reserved.

Ex vivo analysis of T-cell responses to Epstein-Barr virus-encoded oncogene latent membrane protein 1 reveals highly conserved epitope sequences in virus isolates from diverse geographic regions

Epstein-Barr virus (EBV)-encoded oncogene latent membrane protein (LMP) 1, which is consistently expressed in multiple EBV-associated malignancies, has been proposed as a potential target antigen for any future vaccine designed to control these malignancies. However, the high degree of genetic variation in the LMP1 sequence has been considered a major impediment for its use as a potential immunotherapeutic target for the treatment of EBV-associated malignancies. In the present study, we have employed a highly efficient strategy, based on ex vivo functional assays, to conduct an extensive sequence-wide analysis of LMP1-specific T-cell responses in a large panel of healthy virus carriers of diverse ethnic origin and nasopharyngeal carcinoma patients. By comparing the frequencies of T cells specific for overlapping peptides spanning LMP1, we mapped a number of novel HLA class I- and class II-restricted LMP1 T-cell epitopes, including an epitope with dual HLA class I restriction. More importantly, extensive sequence analysis of LMP1 revealed that the majority of the T-cell epitopes were highly conserved in EBV isolates from Caucasian, Papua New Guinean, African, and Southeast Asian populations, while unique geographically constrained genetic variation was observed within one HLA A2 supertype-restricted epitope. These findings indicate that conserved LMP1 epitopes should be considered in designing epitope-based immunotherapeutic strategies against EBV-associated malignancies in different ethnic populations.

The value of idiosyncratic markers and changes to conserved tRNA sequences from the mitochondrial genome of hard ticks (Acari: Ixodida: Ixodidae) for phylogenetic inference

Idiosyncratic markers are features of genes and genomes that are so unusual that it is unlikely that they evolved more than once in a lineage of organisms. Here we explore further the potential of idiosyncratic markers and changes to typically conserved tRNA sequences for phylogenetic inference. Hard ticks were chosen as the model group because their phylogeny has been studied extensively. Fifty-eight candidate markers from hard ticks ( family Ixodidae) and 22 markers from the subfamily Rhipicephalinae sensu lato were mapped onto phylogenies of these groups. Two of the most interesting markers, features of the secondary structure of two different tRNAs, gave strong support to the hypothesis that species of the Prostriata ( Ixodes spp.) are monophyletic. Previous analyses of genes and morphology did not strongly support this relationship, instead suggesting that the Prostriata is paraphyletic with respect to the Metastriata ( the rest of the hard ticks). Parallel or convergent evolution was not found in the arrangements of mitochondrial genes in ticks nor were there any reversals to the ancestral arthropod character state. Many of the markers identified were phylogenetically informative, whereas others should be informative with study of additional taxa. Idiosyncratic markers and changes to typically conserved nucleotides in tRNAs that are phylogenetically informative were common in this data set, and thus these types of markers might be found in other organisms.

Effects of additional sequences directly downstream from the AUG on the expression of GFP gene

We have studied the expression of the green fluorescent protein (GFP) gene to gain more understanding of the effects of additional nucleotide triplets (codons) downstream from the initiation codon on the translation of the GFP mRNA in CHO and Cos1 cells. A leader sequence of six consecutive identical codons (GUG, CUC, AGU or UCA) was introduced into a humanized GFP (hm gfp) gene downstream from the AUG to produce four GFP gene variants. Northern blot and RT-PCR analysis indicated that mRNA transcription from the GFP gene was not significantly affected by any of the additional sequences. However, immunoblotting and FACS analysis revealed that AGU and UCA GFP variants produced GFP at a mean level per cell 3.5-fold higher than the other two GFP variants and the hm gfp gene. [35S]-Methionine labeling and immunoprecipitation demonstrate that GFP synthesis was very active in UCA variant transfected-cells, but not in GUG variant and hm gfp transfected-cells. Moreover, proteasome inhibitor MG-132 treatment indicated that the GFPs encoded by each of the GFP variants and the hm gfp were equally stable, and this together with the comparable mRNA levels observed for each construct suggested that the different steady-state GFP concentrations observed reflected different translation efficiencies of the various GFP genes. In addition, the CUC GFP variant, when transiently transfected into CHO or COS-1 cells, did not produce any GFP expressing cells (fully green cells), and the GUG variant produced GFP expressing cells less than 10%, while AGU and UCA GFP variants up to 30–35% in a time course study from 8 to 36 h posttransfection. Analysis of the potential secondary structure of the GFP variant mRNAs especially in the translation initiation region suggested that the secondary structure of the GFP mRNAs was unlikely to explain the different translation efficiencies of the GFP variants. The present findings indicate that a change of the initiation context of the GFP gene by addition of extra coding sequence can alter the translation efficiency of GFP mRNA, providing a means of more efficient expression of GFP in eukaryotic cells.

Phylogeny of the Sporobolus indicus complex, based on internal transcribed spacer (ITS) sequences

The entire internal transcribed spacer ( ITS) region, including the 5.8S subunit of the nuclear ribosomal DNA ( rDNA), was sequenced by direct double-stranded sequencing of polymerase chain reaction (PCR) amplified fragments. The study included 40 Sporobolus ( Family Poaceae, subfamily Chloridoideae) seed collections from 14 putative species ( all 11 species from the S. indicus complex and three Australian native species). These sequences, along with those from two out-group species [ Pennisetum alopecuroides ( L.) Spreng. and Heteropogon contortus ( L.) P. Beauv. ex Roemer & Schultes, Poaceae, subfamily Panicoideae], were analysed by the parsimony method (PAUP; version 4.0b4a) to infer phylogenetic relationships among these species. The length of the ITS1, 5.8S subunit and ITS2 region were 222, 164 and 218 base pairs ( bp), respectively, in all species of the S. indicus complex, except for the ITS2 region of S. diandrus P. Beauv. individuals, which was 217 bp long. Of the 624 characters included in the analysis, 245 ( 39.3%) of the 330 variable sites contained potential phylogenetic information. Differences in sequences among the members of the S. pyramidalis P. Beauv., S. natalensis (Steud.) Dur & Schinz and S. jacquemontii Kunth. collections were 0%, while differences ranged from 0 to 2% between these and other species of the complex. Similarly, differences in sequences among collections of S. laxus B. K. Simon, S. sessilis B. K. Simon, S. elongatus R. Br. and S. creber De Nardi were 0%, compared with differences of 1-2% between these four species and the rest of the complex. When comparing S. fertilis ( Steud.) Clayton and S. africanus (Poir.) Robyns & Tourney, differences between collections ranged from 0 to 1%. Parsimony analysis grouped all 11 species of the S. indicus complex together, indicating a monophyletic origin. For the entire data set, pair-wise distances among members of the S. indicus complex varied from 0.00 to 1.58%, compared with a range of 20.08-21.44% among species in the complex and the Australian native species studied. A strict consensus phylogenetic tree separated 11 species of the S. indicus complex into five major clades. The phylogeny, based on ITS sequences, was found to be congruent with an earlier study on the taxonomic relationship of the weedy Sporobolus grasses revealed from random amplified polymorphic DNA ( RAPD). However, this cladistic analysis of the complex was not in agreement with that created on past morphological analyses and therefore gives a new insight into the phylogeny of the S. indicus complex.

Role of flanking sequences and phosphorylation in the recognition of the simian-virus-40 large T-antigen nuclear localization sequences by importin-a

The nuclear import of simian-virus-40 large T-antigen (tumour antigen) is enhanced via phosphorylation by the protein kinase CK2 at Ser(112) in the vicinity of the NLS (nuclear localization sequence). To determine the structural basis of the effect of the sequences flanking the basic cluster KKKRK, and the effect of phosphorylation on the recognition of the NLS by the nuclear import factor importin-alpha (Impalpha), we co-crystallized non-autoinhibited Impalpha with peptides corresponding to the phosphorylated and non-phosphorylated forms of the NLS, and determined the crystal structures of the complexes. The structures show that the amino acids N-terminally flanking the basic cluster make specific contacts with the receptor that are distinct from the interactions between bipartite NLSs and Impalpha. We confirm the important role of flanking sequences using binding assays. Unexpectedly, the regions of the peptides containing the phosphorylation site do not make specific contacts with the receptor. Binding assays confirm that phosphorylation does not increase the affinity of the T-antigen NLS to Impalpha. We conclude that the sequences flanking the basic clusters in NLSs play a crucial role in nuclear import by modulating the recognition of the NLS by Impalpha, whereas phosphorylation of the T-antigen enhances nuclear import by a mechanism that does not involve a direct interaction of the phosphorylated residue with Impalpha.