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.      

Cloning of 18S and 25S rDNAs from the pathogenic fungus Cryptococcus neoformans

Cryptococcus neoformans is an important pathogenic fungus that has been classified as a basidiomycete. Little is known of the molecular genetics of this fungal pathogen. To begin such studies, we devised a procedure for extraction of DNA from cryptococci; this method involved the use of the cell wall-active enzyme NovoZym 234. Using cloned rDNA of Saccharomyces cerevisiae as a probe, we identified homologous restriction fragments in a Southern blot of digested C. neoformans DNA. An 8.6-kilobase HindIII fragment that hybridized with the yeast rDNA probe was ligated with the vector pBR322 and cloned into Escherichia coli. When the fragment was used as a probe, it hybridized to the 18S and 25S rRNAs of C. neoformans in Northern (RNA) blots of native and denatured RNA. It bound at high stringency only weakly to the rRNAs of the ascomycete S. cerevisiae. The locations of the genes for 5/5.8S, 18S, and 25S subunits in the cloned fragment were identified with labeled rRNA of these different types.     

Unique GH18 chitinase from Euglena gracilis: full-length cDNA cloning and characterization of its catalytic domain

A cDNA of putative chitinase from Euglena gracilis, designated EgChiA, encoded 960 amino acid residues, which is arranged from N-terminus in the order of signal peptide, glycoside hydrolase family 18 (GH18) domain, carbohydrate binding module family 18 (CBM18) domain, GH18 domain, CBM18 domain, and transmembrane helix. It is likely that EgChiA is anchored on the cell surface. The recombinant second GH18 domain of EgChiA, designated as CatD2, displayed optimal catalytic activity at pH 3.0 and 50 °C. The lower the polymerization degree of the chitin oligosaccharides [(GlcNAc)_4–6] used as the substrates, the higher was the rate of degradation by CatD2. CatD2 degraded chitin nanofibers as an insoluble substrate, and it produced only (GlcNAc)₂ and GlcNAc. Therefore, we speculated that EgChiA localizes to the cell surface of E. gracilis and is involved in degradation of chitin polymers into (GlcNAc)₂ or GlcNAc, which are easily taken up by the cells.     

Diversity in methane enrichments from agricultural soil revealed by DGGE separation of PCR amplified 16S rDNA fragments

Denaturing gradient gel electrophoresis (DGGE) profiles of PCR amplified V3 regions of 16S rRNA genes were used to assess the diversity in enrichment cultures with methane as the only carbon and energy source. The enrichments originated from two agricultural soils. One was a sandy soil with low (10%) organic content, the other an organic soil with approximately 50% organic content. DGGE provided a fast evaluation of the distribution of amplifiable sequence types indicating that specific bacterial populations had been enriched from each soil. The DGGE profiles revealed a broader range of amplified V3 fragments in the community derived from organic soil than from sandy soil. Fragments from 19 individual DGGE bands were sequenced and compared with 27 previously published 16S rRNA gene sequences. The sequences confirmed the high diversity with the presence of different methylotrophic populations in each enrichment. No affiliation was found with type I methanotrophs, instead type II methanotroph sequences were found in the enrichments from both soil types. Some of the fragments from the organic soil enrichment were not affiliated with methylotrophs. Most of the sequences clustered distantly on a branch within the α-Proteobacteria. These facts suggested that previously undescribed methylotrophs are abundant in methane enrichments from agricultural soil.     

Nucleotide sequence of a mosquito 18S ribosomal RNA gene.

We have sequenced an 18S ribosomal RNA gene from the mosquito, Aedes albopictus. Computer alignment of the 1950 nucleotide coding region (56% A + T) with 18S rRNA sequences from two insect and three vertebrate species revealed greater sequence divergence among the insects than among the vertebrates. Sequence alignments showed that variable region V4, which has been considered to be the most poorly conserved domain in the 18S rRNA gene, was better conserved among insects and vertebrates than was the V6 domain.     

树鼩CD4全长编码序列的克隆及分子特征分析

树鼩作为多种人类疾病研究模型的可能性已受到广泛关注,但尚缺乏研究其免疫功能的基本标志以及单克隆抗体。该实验首先以树鼩外周血总RNA为材料,通过RT-PCR扩增得到长度为1365bp的树鼩CD4全长编码序列,并确定了数据库中缺失的两个片段,进而通过ClustalW等软件对其序列和分子特征进行分析,发现树鼩CD4氨基酸序列胞外和胞内域保守性较好,且与人类和猴的亲缘关系较近。虽然树鼩和人CD4分子表面大部分区域均带正电荷,但与人CD4胞外域D1相比,树鼩CD4D1结构区域表面带负电荷较多,且多出两个N-糖基化位点。这些差异对抗体的结合可能存在影响。该研究为今后树鼩CD4单克隆抗体制备及功能研究奠定了基础。     

Cloning and sequence analysis of the 18 S ribosomal RNA gene of tomato and a secondary structure model for the 18 S rRNA of angiosperms

Summary The gene of a cytoplasmic 18 S ribosomal RNA (18 S rDNA) of the dicotyledonous plant tomato (ycopersicon esculentum) cv. Rentita has been cloned, and its complete primary structure has been determined. The tomato 18 S rDNA is 1805 by long with a G+C content of 49.6%. Its sequence exhibits 94%–96% positional identity when it is colinearly aligned with the previously reported sequences of the 17–18 S rDNAs of the dicot soybean and the monocots maize and rice. A model of the secondary structure of the 18 S rRNA of angiosperms is presented and its genera-specific structural features are compared with a current eukaryotic 18 S rRNA consensus model.     

树鼩IL-2全长编码序列的克隆及分子特征分析

树鼩作为多种人类疾病模型已受到广泛关注,而免疫因子对于树鼩模型评价至关重要,但目前对其白细胞介素-2(IL-2)的研究鲜有报道。该实验以经ConA(concanavalin)诱导培养的树鼩淋巴细胞总RNA为模板,RT-PCR克隆出465bp的树鼩IL-2全长编码序列,并采用ClustalW软件分析其序列和分子特征。结果表明树鼩IL-2cDNA编码一个由154个氨基酸组成的蛋白质,其cDNA及氨基酸序列与人的同源性分别为93%及80%,且其整体结构与人IL-2相似。MEGA5.0软件构建的进化树表明,树鼩与人及恒河猴的亲缘关系较近。Pymol软件对树鼩和人IL-2氨基酸序列进行的三维结构模建表明,两者的IL-2分子三维空间结构基本相似,表面大部分区域所带电荷相同,但在某些区域差异较大,且树鼩多出一个糖基化位点,这些差异对抗体的结合可能存在影响。该研究为今后树鼩IL-2单克隆抗体的制备及功能研究奠定了基础。     

Leishmania: genus identification based on a specific sequence of the 18S ribosomal RNA sequence

The analysis of PvuII restriction patterns of Leishmania spp. and Trypanosoma spp. genomic DNA showed genus distinctive profiles. A specific PvuII site was detected in the 5' domain of 18S ribosomal DNA of Leishmania. A 20-mer oligonucleotide encompassing this PvuII region was synthesized. This sequence, when utilized as probe, on short exposures of dot tests, detected 10(3) whole promastigotes of all Leishmania species analyzed but did not hybridize with T. cruzi or human nucleic acids. Two other oligonucleotides were synthesized to be used as primers for amplification through polymerase chain reaction of the 18S ribosomal DNA region containing the PvuII site. The probes described may be useful for the detection of Leishmania spp. under clinical and epidemiological trials.     

Cloning of full-length coding sequence of tree shrew CD3E and prediction of its molecular characteristics

The use of tree shrews (Tupaia belangeri) in human disease studies demands essential research tools, in particular cellular markers and their monoclonal antibodies for immunological studies. Here we cloned the full-length cDNAs encoding CD3E from total RNA of the spleen, liver and peripheral blood of tree shrews and analyzed their structural characteristics in comparison with other mammals by Discovery Studio software. The results showed that the open reading frame sequence of tree shrew CD3E was 582 bp, encoding 194 amino acids. The overall structure of tree shrew CD3E protein was similar to its counterparts of other mammals, intracellular and transmembrane domain highly conserved. However, detailed analysis revealed two potential glycosylation sites and different surface charges in the extracellular domain. Availability of the entire open-reading-frame and related sequence information would therefore facilitate the preparation of monoclonal antibodies against tree shrew CD3 and further studies for its function.     

Molecular cloning and sequence analysis of full-length cDNA coding for mouse contrapsin

Four overlapping cDNA clones encoding contrapsin were isolated from a mouse liver cDNA library constructed in the expression vector, lambda gt11. M13 vector sequence analysis revealed that contrapsin cDNA contained an open reading frame of 1,254 bases encoding 418 amino acids. The N-terminal amino acid sequence of the isolated contrapsin matched residues 30 to 48 of the sequence deduced on nucleotide analysis. One clone, which had the longest 3' untranslated region, contained two sets of tandem polyadenylation signals, AATACA and AATAAA, which were located 497 bases apart, while the remaining three clones terminated at the first signal. The entire reading frame sequence of contrapsin cDNA showed 64% homology with that of human alpha-1-antichymotrypsin.     

Nucleotide sequence of the 3''-terminal region of rat 18S ribosomal DNA

Summary The 3-terminal 230 base-pairs (bp) of the gene for 18S rRNA and 40 bp of the adjoining spacer have been sequenced for the Sprague-Dawley rat. This mammalian sequence has been compared with the known sequences of yeast, fruit fly, silkworm, and frog. This study has shown that the nucleotide-sequence differences between rat and frog are the smallest among the five species, probably reflecting their evolutionary closeness and longer maturation time compared to the others. There is little similarity in the nucleotide sequences of the transcribed spacer regions of the five species compared.     

Nucleotide sequence of the 5'-terminal region of rat 18S ribosomal DNA

The 5'-terminal 597 base-pairs (bp) of the Sprague-Dawley rat 18S ribosomal RNA gee and 10 bp of the adjoining transcribed spacer have been sequenced. Previously sequenced 10 large oligonucleotides of rat 18S RNA were located in this region. This mammalian sequence has been compared with the known sequences of yeast and frog 18S rDNA's. The analysis indicates that 534 bp of the 597 bp (89%) are conserved between rat and frog sequences but only 75% of the nucleotides are conserved between rat and yeast in this region. Two large and two small sections have been identified where insertions have been introduced during evolution. Of these 58 bp long inserted sections of the rat rDNA sequence, 50 bp (86%) were G-C base-pairs.     

Isolation and sequence analysis of a full-length cDNA for human M creatine kinase

A full length cDNA for human M creatine kinase has been isolated and sequenced. The cDNA contains 77 bp of 5' untranslated, 338 bp of 3' untranslated sequence and the entire coding region (1146 bp) for human M creatine kinase. The M creatine kinases from different species share considerable sequence homology within the coding region (77-91%) and in amino acid sequence (82-97%). Little or no sequence homology is observed in the 3' untranslated sequence of the mammalian M creatine kinases, although canine and human creatine kinase share overall 80% sequence homology in 5' untranslated sequence. A unique 8 bp sequence was identified in the 5' untranslated regions of mammalian M creatine kinase but is not present in B creatine kinase cDNA. The degree of sequence conservation observed implies an evolutionary constraint on M creatine kinase structure beyond that which would be expected for the maintenance of enzymatic function.