Effects of nucleotide sequence alignment on phylogeny estimation: a case study of 18S rDNAs of apicomplexa

The reconstruction of phylogenetic history is predicated on being able to accurately establish hypotheses of character homology, which involves sequence alignment for studies based on molecular sequence data. In an empirical study investigating nucleotide sequence alignment, we inferred phylogenetic trees for 43 species of the Apicomplexa and 3 of Dinozoa based on complete small-subunit rDNA sequences, using six different multiple-alignment procedures: manual alignment based on the secondary structure of the 18S rRNA molecule, and automated similarity-based alignment algorithms using the PileUp, ClustalW, TreeAlign, MALIGN, and SAM computer programs. Trees were constructed using neighboring-joining, weighted-parsimony, and maximum- likelihood methods. All of the multiple sequence alignment procedures yielded the same basic structure for the estimate of the phylogenetic relationship among the taxa, which presumably represents the underlying phylogenetic signal. However, the placement of many of the taxa was sensitive to the alignment procedure used; and the different alignments produced trees that were on average more dissimilar from each other than did the different tree-building methods used. The multiple alignments from the different procedures varied greatly in length, but aligned sequence length was not a good predictor of the similarity of the resulting phylogenetic trees. We also systematically varied the gap weights (the relative cost of inserting a new gap into a sequence or extending an already-existing gap) for the ClustalW program, and this produced alignments that were at least as different from each other as those produced by the different alignment algorithms. Furthermore, there was no combination of gap weights that produced the same tree as that from the structure alignment, in spite of the fact that many of the alignments were similar in length to the structure alignment. We also investigated the phylogenetic information content of the helical and nonhelical regions of the rDNA, and conclude that the helical regions are the most informative. We therefore conclude that many of the literature disagreements concerning the phylogeny of the Apicomplexa are probably based on differences in sequence alignment strategies rather than differences in data or tree-building methods.      

Cross-hybridizing snake satellite, Drosophila, and mouse DNA sequences may have arisen independently

Previous reports have interpreted hybridization between snake satellite DNA and DNA clones from a variety of distant taxonomic groups as evidence for evolutionary conservation, which implies common ancestry (homology) and/or convergence (analogy) to produce the cross- hybridizing sequences. We have isolated 11 clones from a genomic library of Drosophila melanogaster, using a cloned 2.5-kb snake satellite probe of known nucleotide sequence. We have also analysed published sequence data from snakes, mice, and Drosophila. These data show that (1) all of the cross-hybridization between the snake, fly, and mouse clones can be accounted for by the presence of either of two tandem repeats, [GATA]n and [GACA]n and (2) these tandem repeats are organized differently among the different species. We find no evidence that these sequences are homologous apart from the existence of the simple repeat itself, although their divergence from a common ancestral sequence cannot be ruled out. The sequences contain a variety of homogeneous clusters of tandem repeats of CATA, GA, TA, and CA, as well as GATA and GACA. We suggest that these motifs may have arisen by a self-accelerating process involving slipped-strand mispairing of DNA. Homogeneity of the clusters might simply be the result of a rate of accumulation of tandem repeats that exceeds that of other mutations.      

Structural comparison of fibroblast growth factor-specific heparan sulfates derived from a growing or differentiating neuroepithelial cell line

Heparan sulfate (HS) glycosaminoglycans are essential modulators of fibroblast growth factor (FGF) activity both in vivo and in vitro, and appear to act by cross-linking particular forms of FGF to appropriate FGF receptors. We have recently isolated and characterized two separate HS pools derived from immortalized embryonic day 10 mouse neuroepithelial 2.3D cells: one from cells in log growth phase, which greatly potentiates the activity of FGF-2, and the other from cells undergoing contact-inhibition and differentiation, which preferentially activates FGF-1. These two pools of HS have very similar functional activities to those species isolated from primary neuroepithelial cells at corresponding stages of active proliferation or differentiation. We present here a structural comparison between these cell line HS species to establish the nature of the changes that occur in the biosynthesis of HS. A combination of chemical and enzymatic cleavage, low pressure chromatography and strong anion-exchange HPLC were used to generate full chain models of each species. Overall, the HS pools synthesized in the dividing cell line pools possessed less complex sulfation than those derived from more differentiated, growth arrested cells.      

Kinetic studies of the reverse reaction catalysed by adenosine triphosphate–creatine phosphotransferase. The inhibition by magnesium ions and adenosine diphosphate

1. Kinetic investigations of the reaction catalysed by ATP–creatine phosphotransferase have been carried out. 2. No firm conclusions could be reached about the reaction of Mg²⁺ at the nucleotide-binding site of the enzyme. The value of the kinetic constant for this reaction depends on the value used for the apparent stability constant of the metal ion–nucleotide complex and, to a smaller extent, on the method of plotting the results. 3. At higher concentrations Mg²⁺ is a non-competitive inhibitor of the enzyme with respect to both MgADP⁻ and phosphocreatine. 4. ADP³⁻ is a competitive inhibitor of the enzyme with respect to MgADP⁻ and a non-competitive inhibitor with respect to phosphocreatine. 5. The concentration of phosphocreatine has little, if any, effect on the kinetic constants for the nucleotide reactants.     

Complex phylogeographic history of central African forest elephants and its implications for taxonomy

Background

The mitochondrial genomes of snakes are characterized by an overall evolutionary rate that appears to be one of the most accelerated among vertebrates. They also possess other unusual features, including short tRNAs and other genes, and a duplicated control region that has been stably maintained since it originated more than 70 million years ago. Here, we provide a detailed analysis of evolutionary dynamics in snake mitochondrial genomes to better understand the basis of these extreme characteristics, and to explore the relationship between mitochondrial genome molecular evolution, genome architecture, and molecular function. We sequenced complete mitochondrial genomes from Slowinski's corn snake (Pantherophis slowinskii) and two cottonmouths (Agkistrodon piscivorus) to complement previously existing mitochondrial genomes, and to provide an improved comparative view of how genome architecture affects molecular evolution at contrasting levels of divergence.

Results

We present a Bayesian genetic approach that suggests that the duplicated control region can function as an additional origin of heavy strand replication. The two control regions also appear to have different intra-specific versus inter-specific evolutionary dynamics that may be associated with complex modes of concerted evolution. We find that different genomic regions have experienced substantial accelerated evolution along early branches in snakes, with different genes having experienced dramatic accelerations along specific branches. Some of these accelerations appear to coincide with, or subsequent to, the shortening of various mitochondrial genes and the duplication of the control region and flanking tRNAs.

Conclusion

Fluctuations in the strength and pattern of selection during snake evolution have had widely varying gene-specific effects on substitution rates, and these rate accelerations may have been functionally related to unusual changes in genomic architecture. The among-lineage and among-gene variation in rate dynamics observed in snakes is the most extreme thus far observed in animal genomes, and provides an important study system for further evaluating the biochemical and physiological basis of evolutionary pressures in vertebrate mitochondria.     

Retroviral Splicing Suppressor Requires Three Nonconsensus Uridines in a 5′ Splice Site-Like Sequence

Rous sarcoma virus RNA contains a negative regulator of splicing (NRS) element that aids in maintenance of unspliced RNA. The NRS binds U1 snRNA at a sequence that deviates from the 5' splice site consensus by substitution of U's for A's at three positions: -2, +3, and +4. All three of these U's are important for NRS-mediated splicing suppression. Substitution of a single nonconsensus C or G at any of these sites diminished NRS activity, whereas substitution of a single A generated a preferred 5' splice site within the NRS.     

Structural requirements for efficient prion protein conversion: Cofactors may promote a conversion-competent structure for PrPC

To understand why cross-species infection of prion disease often results in inefficient transmission and reduced protein conversion, most research has focused on defining the effect of variations in PrP primary structures, including sequence compatibility of substrate and seed. By contrast, little research has been aimed at investigating structural differences between different variants of PrP^C and secondary structural requirements for efficient conversion. This is despite a clear role for molecular chaperones in formation of prions in non-mammalian systems, indicating the importance of secondary/tertiary structure during the conversion process. Recent data from our laboratory on the cellular location of disease-specific prion cofactors supports the critical role of specific secondary structural motifs and the stability of these motifs in determining the efficiency of disease-specific prion protein conversion. In this paper we summarize our recent results and build on the hypothesis previously suggested by Wuthrich and colleagues, that stability of certain regions of the prion protein is crucial for protein conversion to abnormal isoforms in vivo. It is suggested that one role for molecular cofactors in the conversion process is to stabilize PrP^C structure in a form that is amenable for conversion to PrP^Sc.Key words: cofactor, structure, cell-free conversion assay, fibrillization, stability, loop region     

A comprehensive sensitivity analysis of microarray breast cancer classification under feature variability

Background  

Large discrepancies in signature composition and outcome concordance have been observed between different microarray breast cancer expression profiling studies. This is often ascribed to differences in array platform as well as biological variability. We conjecture that other reasons for the observed discrepancies are the measurement error associated with each feature and the choice of preprocessing method. Microarray data are known to be subject to technical variation and the confidence intervals around individual point estimates of expression levels can be wide. Furthermore, the estimated expression values also vary depending on the selected preprocessing scheme. In microarray breast cancer classification studies, however, these two forms of feature variability are almost always ignored and hence their exact role is unclear.     

Maps of the Sri Lanka malaria situation preceding the tsunami and key aspects to be considered in the emergency phase and beyond

In Vietnam, a large proportion of all malaria cases and deaths occurs in the central mountainous and forested part of the country. Indeed, forest malaria, despite intensive control activities, is still a major problem which raises several questions about its dynamics. A large-scale malaria morbidity survey to measure malaria endemicity and identify important risk factors was carried out in 43 villages situated in a forested area of Ninh Thuan province, south central Vietnam. Four thousand three hundred and six randomly selected individuals, aged 10–60 years, participated in the survey. Rag Lays (86%), traditionally living in the forest and practising "slash and burn" cultivation represented the most common ethnic group. The overall parasite rate was 13.3% (range [0–42.3] while Plasmodium falciparum seroprevalence was 25.5% (range [2.1–75.6]). Mapping of these two variables showed a patchy distribution, suggesting that risk factors other than remoteness and forest proximity modulated the human-vector interactions. This was confirmed by the results of the multivariate-adjusted analysis, showing that forest work was a significant risk factor for malaria infection, further increased by staying in the forest overnight (OR= 2.86; 95%CI [1.62; 5.07]). Rag Lays had a higher risk of malaria infection, which inversely related to education level and socio-economic status. Women were less at risk than men (OR = 0.71; 95%CI [0.59; 0.86]), a possible consequence of different behaviour. This study confirms that malaria endemicity is still relatively high in this area and that the dynamics of transmission is constantly modulated by the behaviour of both humans and vectors. A well-targeted intervention reducing the "vector/forest worker" interaction, based on long-lasting insecticidal material, could be appropriate in this environment.     

The conversion of cephalosporins to 7 alpha-methoxycephalosporins by cell-free extracts of Streptomyces clavuligerus.

In the presence of S-adenosylmethionine, 2-oxoglutarate, Fe2+ and a reducing agent, cell-free extracts of Streptomyces clavuligerus convert cephalosporin C and O-carbamoyldeacetylcephalosporin C into 7 alpha-methoxy derivatives. No synthesis of a 7 alpha-methoxy derivative of deacetylcephalosporin C was detected in the system used, and the 7 alpha-methoxy derivative of deacetoxycephalosporin C was produced only in relatively small amounts. It appears that the 7 alpha-methoxy group is introduced after the cephalosporin ring system has been formed and that its introduction may represent the final step in a biosynthetic pathway.