Phylogenetic relationships of discoglossid frogs (Amphibia:Anura:Discoglossidae) based on complete mitochondrial genomes and nuclear genes

The complete nucleotide sequence of the mitochondrial (mt) genome was determined for three species of discoglossid frogs (Amphibia:Anura:Discoglossidae), representing three of the four recognized genera: Alytes obstetricans, Bombina orientalis, and Discoglossus galganoi. The organization and size of these newly determined mt genomes are similar to those previously reported for other vertebrates. Phylogenetic analyses (maximum likelihood, Bayesian inference, minimum evolution, and maximum parsimony) of mt protein-coding genes at the amino acid level were performed in combination with already published mt genome sequence data of three species of Neobatrachia, one of Pipoidea, and four of Caudata. Phylogenetic analyses based on the deduced amino acid sequences of all mt protein-coding genes arrived at the same topology. The monophyly of Discoglossidae is strongly supported. Within the Discoglossidae, Alytes is consistently recovered as sister group of Discoglossus, to the exclusion of Bombina. The three species representing Neobatrachia exhibited extremely long branches irrespective of the phylogenetic inference method used, and hence their relative position with respect to Discoglossidae and Xenopus may be artefactual due to a severe long branch attraction effect. To further investigate the phylogenetic intrarelationships of discoglossids, nucleotide sequences of four nuclear protein-coding genes (CXCR4, RAG1, RAG2, and Rhodopsin) with sequences available for the three discoglossid genera and Xenopus were retrieved from GenBank, and together with a concatenated nucleotide sequence data set containing all mt protein-coding genes except ND6 were subjected to separate and combined phylogenetic analyses. In all cases, a sister group relationship between Alytes and Discoglossus was recovered with high statistical support.     

ON COMBINING PROTEIN SEQUENCES AND NUCLEIC ACID SEQUENCES IN PHYLOGENETIC ANALYSIS: THE HOMEOBOX PROTEIN CASE

Abstract— Amino acid encoding genes contain character state information that may be useful for phylogenetic analysis on at least two levels. The nucleotide sequence and the translated amino acid sequences have both been employed separately as character states for cladistic studies of various taxa, including studies of the genealogy of genes in multigene families. In essence, amino acid sequences and nucleic acid sequences are two different ways of character coding the information in a gene. Silent positions in the nucleotide sequence (first or third positions in codons that can accrue change without changing the identity of the amino acid that the triplet codes for) may accrue change relatively rapidly and become saturated, losing the pattern of historical divergence. On the other hand, non-silent nucleotide alterations and their accompanying amino acid changes may evolve too slowly to reveal relationships among closely related taxa. In general, the dynamics of sequence change in silent and non-silent positions in protein coding genes result in homoplasy and lack of resolution, respectively. We suggest that the combination of nucleic acid and the translated amino acid coded character states into the same data matrix for phylogenetic analysis addresses some of the problems caused by the rapid change of silent nucleotide positions and overall slow rate of change of non-silent nucleotide positions and slowly changing amino acid positions. One major theoretical problem with this approach is the apparent non-independence of the two sources of characters. However, there are at least three possible outcomes when comparing protein coding nucleic acid sequences with their translated amino acids in a phylogenetic context on a codon by codon basis. First, the two character sets for a codon may be entirely congruent with respect to the information they convey about the relationships of a certain set of taxa. Second, one character set may display no information concerning a phylogenetic hypothesis while the other character set may impart information to a hypothesis. These two possibilities are cases of non-independence, however, we argue that congruence in such cases can be thought of as increasing the weight of the particular phylogenetic hypothesis that is supported by those characters. In the third case, the two sources of character information for a particular codon may be entirely incongruent with respect to phylogenetic hypotheses concerning the taxa examined. In this last case the two character sets are independent in that information from neither can predict the character states of the other. Examples of these possibilities are discussed and the general applicability of combining these two sources of information for protein coding genes is presented using sequences from the homeobox region of 46 homeobox genes fromDrosophila melanogasterto develop a hypothesis of genealogical relationship of these genes in this large multigene family.     

Cloning and characterization of the bovine testicular PH-20 hyaluronidase core domain

The core nucleotide sequence of bovine (Bos taurus) testicular PH-20 hyaluronidase was cloned using one step RT-PCR. The 5' and 3' regions were cloned separately and a sequence overlap of 124 bp facilitated the fusion of these two fragments by overlapping PCR, resulting in a concatenated sequence of 1422 bp. This nucleotide sequence and its deduced amino acid sequence were compared to homologous sequences from eight other mammal species. The bovine sequences were most similar to those of the pig, Sus scrofa (swine Spam1: 79.1% nucleotide and 70.1% amino acid similarity) and least similar to sequences from the Norway rat, Rattus norvegicus (murine Spam1: 61% nucleotide and 53.3% amino acid similarity). A phylogenetic analysis joined the red fox (Vulpes vulpes) sequence as sister to the bull-pig pair. Twelve cysteine residues were conserved among all nine aligned amino acid sequences and five proposed glycosylation sites have been identified. The feasibility of developing an effective, low-cost bovine PH-20 expression system is discussed in light of these new data.     

植物LEAFY 同源基因的研究进展

本文就近10年来LEAFY(简写为LFY)同源基因的研究进展做了综合分析。通过对19种植物中已分离到的LFY同源基因的序列比较分析发现: LFY同源基因编码区核苷酸和氨基酸序列同源性都较高;在双子叶植物基因组中, 拷贝数却有所不同。该基因的表达特性显示其在不同植物中表达的时间和空间有所差异。根据已知序列推导的氨基酸序列构建的系统进化树表明, 单子叶植物与裸子植物的亲缘关系近于双子叶与裸子植物的亲缘关系。上述研究资料为植物成花机理研究提供了重要参考, 且在研究植物系统进化方面也具有重要的意义。     

Cladistic analysis of anuran POMC sequences

Procedures for performing cladistic analyses can provide powerful tools for understanding the evolution of neuropeptide and polypeptide hormone coding genes. These analyses can be done on either amino acid data sets or nucleotide data sets and can utilize several different algorithms that are dependent on distinct sets of operating assumptions and constraints. In some cases, the results of these analyses can be used to gauge phylogenetic relationships between taxa. Selecting the proper cladistic analysis strategy is dependent on the taxonomic level of analysis and the rate of evolution within the orthologous genes being evaluated. For example, previous studies have shown that the amino acid sequence of proopiomelanocortin (POMC), the common precursor for the melanocortins and beta-endorphin, can be used to resolve phylogenetic relationships at the class and order level. This study tested the hypothesis that POMC sequences could be used to resolve phylogenetic relationships at the family taxonomic level. Cladistic analyses were performed on amphibian POMC sequences characterized from the marine toad, Bufo marinus (family Bufonidae; this study), the spadefoot toad, Spea multiplicatus (family Pelobatidae), the African clawed frog, Xenopus laevis (family Pipidae) and the laughing frog, Rana ridibunda (family Ranidae). In these analyses the sequence of Australian lungfish POMC was used as the outgroup. The analyses were done at the amino acid level using the maximum parsimony algorithm and at the nucleotide level using the maximum likelihood algorithm. For the anuran POMC genes, analysis at the nucleotide level using the maximum likelihood algorithm generated a cladogram with higher bootstrap values than the maximum parsimony analysis of the POMC amino acid data set. For anuran POMC sequences, analysis of nucleotide sequences using the maximum likelihood algorithm would appear to be the preferred strategy for resolving phylogenetic relationships at the family taxonomic level.     

Self-incompatibility alleles in Polish wild pear (Pyrus pyraster (L.) Burgsd.): a preliminary analysis

Wild pear (Pyrus pyraster, syn.P. communis var.pyraster) is thought to be one of the species that gave rise to all other members of the genusPyrus, although intraspecific hybridizations with cultivated varieties could cause the disappearance of original species characteristics. S-RNase alleles from 7 different wild pear individuals, collected from various regions of Poland, were cloned on the basis of the PCR method and nucleotide sequence analyses. The hypervariable (HV) region is responsible for allele-specific S-RNase activity in the self-incompatibility mechanism. The high level of polymorphism of its sequences may constitute a source of valuable phylogenetic information. From all individuals, 14 sequences were obtained successfully, and 9 of them were novel alleles. Phylogenetic analysis of these alleles was based on the amino acid sequence interpretation of coding regions and intron nucleotide sequences. The research conducted on a limited pool of availableP. pyraster alleles gives only an initial insight into possible S-RNase allele polymorphisms in wild populations. At this stage, the results do not confirm a strong influence of cultivated pear species on the wild pear.     

Granule-bound starch synthase: structure, function, and phylogenetic utility

Interest in the use of low-copy nuclear genes for phylogenetic analyses of plants has grown rapidly, because highly repetitive genes such as those commonly used are limited in number. Furthermore, because low- copy genes are subject to different evolutionary processes than are plastid genes or highly repetitive nuclear markers, they provide a valuable source of independent phylogenetic evidence. The gene for granule-bound starch synthase (GBSSI or waxy) exists in a single copy in nearly all plants examined so far. Our study of GBSSI had three parts: (1) Amino acid sequences were compared across a broad taxonomic range, including grasses, four dicotyledons, and the microbial homologs of GBSSI. Inferred structural information was used to aid in the alignment of these very divergent sequences. The informed alignments highlight amino acids that are conserved across all sequences, and demonstrate that structural motifs can be highly conserved in spite of marked divergence in amino acid sequence. (2) Maximum-likelihood (ML) analyses were used to examine exon sequence evolution throughout grasses. Differences in probabilities among substitution types and marked among-site rate variation contributed to the observed pattern of variation. Of the parameters examined in our set of likelihood models, the inclusion of among-site rate variation following a gamma distribution caused the greatest improvement in likelihood score. (3) We performed cladistic parsimony analyses of GBSSI sequences throughout grasses, within tribes, and within genera to examine the phylogenetic utility of the gene. Introns provide useful information among very closely related species, but quickly become difficult to align among more divergent taxa. Exons are variable enough to provide extensive resolution within the family, but with low bootstrap support. The combined results of amino acid sequence comparisons, maximum-likelihood analyses, and phylogenetic studies underscore factors that might affect phylogenetic reconstruction. In this case, accommodation of the variable rate of evolution among sites might be the first step in maximizing the phylogenetic utility of GBSSI.      

我国2009年新分离乙脑病毒全基因组序列特征研究

本研究对我国2009年新分离的两株乙脑病毒进行全基因组序列测定和分析,以了解病毒全基因组分子特征。通过RT-PCR和核苷酸序列测定方法获得病毒全基因组序列,采用ClustalX、DNASTAR、MEGA等生物学软件完成核苷酸序列及氨基酸序列分析和系统进化分析等。研究结果显示,新分离两株乙脑病毒YN0911和YN0967株基因组全长均为10 965个核苷酸,编码3 432个氨基酸。这2株乙脑病毒之间核苷酸同源性为98.7%,氨基酸同源性为99.8%。与国际乙脑病毒流行株相比,核苷酸同源性为83.5%～98.9%,氨基酸同源性为94.8%～99.7%。与乙脑病毒疫苗株SA14-14-2相比,在E蛋白有13个氨基酸差异位点,但都位于抗原关键位点之外。这2株病毒在3′UTR区域存在11nt缺失。基于C/PrM区段、E基因、全基因组系统进化分析结果均显示新分离2株乙脑病毒为G I乙脑病毒,并且和越南、四川、贵州、广西以往的分离株遗传进化关系较近。本研究提示我国新分离的2株乙脑病毒均为G I乙脑病毒,决定病毒毒力的关键氨基酸位点未见明显变化。     

Sequence diversity within the reovirus S2 gene: reovirus genes reassort in nature, and their termini are predicted to form a panhandle motif.

To better understand genetic diversity within mammalian reoviruses, we determined S2 nucleotide and deduced sigma 2 amino acid sequences of nine reovirus strains and compared these sequences with those of prototype strains of the three reovirus serotypes. The S2 gene and sigma 2 protein are highly conserved among the four type 1, one type 2, and seven type 3 strains studied. Phylogenetic analyses based on S2 nucleotide sequences of the 12 reovirus strains indicate that diversity within the S2 gene is independent of viral serotype. Additionally, we found marked topological differences between phylogenetic trees generated from S1 and S2 gene nucleotide sequences of the seven type 3 strains. These results demonstrate that reovirus S1 and S2 genes have distinct evolutionary histories, thus providing phylogenetic evidence for lateral transfer of reovirus genes in nature. When variability among the 12 sigma 2-encoding S2 nucleotide sequences was analyzed at synonymous positions, we found that approximately 60 nucleotides at the 5' terminus and 30 nucleotides at the 3' terminus were markedly conserved in comparison with other sigma 2-encoding regions of S2. Predictions of RNA secondary structures indicate that the more conserved S2 sequences participate in the formation of an extended region of duplex RNA interrupted by a pair of stem-loops. Among the 12 deduced sigma 2 amino acid sequences examined, substitutions were observed at only 11% of amino acid positions. This finding suggests that constraints on the structure or function of sigma 2, perhaps in part because of its location in the virion core, have limited sequence diversity within this protein.     

The c-sis gene encodes a precursor of the B chain of platelet-derived growth factor.

The relationship between platelet-derived growth factor (PDGF) and the proto-oncogene c-sis has been determined by amino acid sequence analysis of PDGF and nucleotide sequence analysis of c-sis genomic clones. The nucleotide sequences of five regions of the human c-sis gene which are homologous to sequences of the transforming region (v-sis) of simian sarcoma virus (SSV) were determined. By alignment of the c-sis and v-sis nucleotide sequences the predicted amino acid sequence of a polypeptide homologous to the putative transforming protein p28sis of SSV was deduced. Both predicted sequences use the same termination codon and additional coding sequences may lie 5' to the homologous regions. Amino acid sequence analysis of the PDGF B chain shows identity to the amino acid sequence predicted from the c-sis sequences over 109 amino acid residues. Polymorphism may exist at two amino acid residues. These results suggest that c-sis encodes a polypeptide precursor of the B chain. A partial amino acid sequence of the PDGF A chain is also described. This chain is 60% homologous to the B chain and cannot be encoded by that part of c-sis which has been sequenced but could be encoded by sequences which lie 5' to the five regions of v-sis homology in c-sis, or at a separate locus.     

Phylogenetic network for European mtDNA

The sequence in the first hypervariable segment (HVS-I) of the control region has been used as a source of evolutionary information in most phylogenetic analyses of mtDNA. Population genetic inference would benefit from a better understanding of the variation in the mtDNA coding region, but, thus far, complete mtDNA sequences have been rare. We determined the nucleotide sequence in the coding region of mtDNA from 121 Finns, by conformation-sensitive gel electrophoresis and subsequent sequencing and by direct sequencing of the D loop. Furthermore, 71 sequences from our previous reports were included, so that the samples represented all the mtDNA haplogroups present in the Finnish population. We found a total of 297 variable sites in the coding region, which allowed the compilation of unambiguous phylogenetic networks. The D loop harbored 104 variable sites, and, in most cases, these could be localized within the coding-region networks, without discrepancies. Interestingly, many homoplasies were detected in the coding region. Nucleotide variation in the rRNA and tRNA genes was 6%, and that in the third nucleotide positions of structural genes amounted to 22% of that in the HVS-I. The complete networks enabled the relationships between the mtDNA haplogroups to be analyzed. Phylogenetic networks based on the entire coding-region sequence in mtDNA provide a rich source for further population genetic studies, and complete sequences make it easier to differentiate between disease-causing mutations and rare polymorphisms.     

Chloroplast phylogeny indicates that bryophytes are monophyletic

Opinions on the basal relationship of land plants vary considerably and no phylogenetic tree with significant statistical support has been obtained. Here, we report phylogenetic analyses using 51 genes from the entire chloroplast genome sequences of 20 representative green plant species. The analyses, using translated amino acid sequences, indicated that extant bryophytes (mosses, liverworts, and hornworts) form a monophyletic group with high statistical confidence and that extant bryophytes are likely sisters to extant vascular plants, although the support for monophyletic vascular plants was not strong. Analyses at the nucleotide level could not resolve the basal relationship with statistical confidence. Bryophyte monophyly inferred using amino acid sequences has a good statistical foundation and is not rejected statistically by other data sets. We propose bryophyte monophyly as the currently best hypothesis.     

Beta-subunit of ATP-synthase: a useful marker for studying the phylogenetic relationship of eubacteria

The genes encoding the beta-subunits of ATP-synthases (ATPases) from Bacteroides fragilis DSM 2151, Cytophaga lytica DSM 2039 and 'Taxeobacter ocellatus' were cloned. The nucleotide sequences were determined completely for the genes of the first two organisms and to a major part for that of 'T. ocellatus'. The predicted amino acid sequences were compared with previously published amino acid sequences of beta-subunits. Two characteristic insertions were found in genes from organisms belonging to the so-called bacteroides-cytophaga-flavo-bacterium group. The remaining structure shows a high degree of sequence similarity within this group. These data support the conclusions drawn from comparative 16S rRNA sequence analyses that organisms in this phenotypically heterogeneous group are phylogenetically related. A phylogenetic tree was constructed based on a distance matrix of optimally aligned amino acid sequences of beta-subunits of ATPases of various eubacteria, chloroplasts and mitochondria. It is in good agreement with a tree derived from 16S rRNA sequence analyses.     

A comparison of complete genome sequences of a rabies virus chinese isolate SH06 with the vaccine strains

In this study, we determined the complete nucleotide and deduced amino acid sequence of a primary isolate of rabies virus (SH06) obtained from the brain of a rabid dog. The overall length of the genome was 11 924 nucleotides. Comparison of the genomic sequence showed the homology of SH06 at nucleotide level with full-length genomes of reference vaccine strains ranged from 82.2% with the PV strain to 86.9% with the CTN strain. A full-length genome-based phylogenetic analysis was performed with sequences available from GenBank. Phylogenetic analysis of the complete genome sequences indicated that the SH06 exhibited the highest homology with rabies street virus BD06 and CTN vaccine strain originated from China.     

Evaluation and refinement of tmRNA structure using gene sequences from natural microbial communities

DNA harvested directly from complex natural microbial communities by PCR has been successfully used to predict RNase P RNA structure, and can potentially provide an abundant source of information for structural predictions of other RNAs. In this study, we utilized genetic variation in natural communities to test and refine the secondary and tertiary structural model for the bacterial tmRNA. The variability of proposed tmRNA secondary structures in different organisms and the lack of any predicted tertiary structure suggested that further refinement of the tmRNA could be useful. To increase the phylogenetic representation of tmRNA sequences, and thereby provide additional data for statistical comparative analysis, we amplified, sequenced, and compared tmRNA sequences from natural microbial communities. Using primers designed from gamma proteobacterial sequences, we determined 44 new tmRNA sequences from a variety of environmental DNA samples. Covariation analyses of these sequences, along with sequences from cultured organisms, confirmed most of the proposed tmRNA model but also provided evidence for a new tertiary interaction. This approach of gathering sequence information from natural microbial communities seems generally applicable in RNA structural analysis.     

G-InforBIO: integrated system for microbial genomics

Background  

Genome databases contain diverse kinds of information, including gene annotations and nucleotide and amino acid sequences. It is not easy to integrate such information for genomic study. There are few tools for integrated analyses of genomic data, therefore, we developed software that enables users to handle, manipulate, and analyze genome data with a variety of sequence analysis programs.