Phylogenetic relationships of agaric fungi based on nuclear large subunit ribosomal DNA sequences

Phylogenetic relationships of mushrooms and their relatives within the order Agaricales were addressed by using nuclear large subunit ribosomal DNA sequences. Approximately 900 bases of the 5' end of the nucleus-encoded large subunit RNA gene were sequenced for 154 selected taxa representing most families within the Agaricales. Several phylogenetic methods were used, including weighted and equally weighted parsimony (MP), maximum likelihood (ML), and distance methods (NJ). The starting tree for branch swapping in the ML analyses was the tree with the highest ML score among previously produced MP and NJ trees. A high degree of consensus was observed between phylogenetic estimates obtained through MP and ML. NJ trees differed according to the distance model that was used; however, all NJ trees still supported most of the same terminal groupings as the MP and ML trees did. NJ trees were always significantly suboptimal when evaluated against the best MP and ML trees, by both parsimony and likelihood tests. Our analyses suggest that weighted MP and ML provide the best estimates of Agaricales phylogeny. Similar support was observed between bootstrapping and jackknifing methods for evaluation of tree robustness. Phylogenetic analyses revealed many groups of agaricoid fungi that are supported by moderate to high bootstrap or jackknife values or are consistent with morphology-based classification schemes. Analyses also support separate placement of the boletes and russules, which are basal to the main core group of gilled mushrooms (the Agaricineae of Singer). Examples of monophyletic groups include the families Amanitaceae, Coprinaceae (excluding Coprinus comatus and subfamily Panaeolideae), Agaricaceae (excluding the Cystodermateae), and Strophariaceae pro parte (Stropharia, Pholiota, and Hypholoma); the mycorrhizal species of Tricholoma (including Leucopaxillus, also mycorrhizal); Mycena and Resinomycena; Termitomyces, Podabrella, and Lyophyllum; and Pleurotus with Hohenbuehelia. Several groups revealed by these data to be nonmonophyletic include the families Tricholomataceae, Cortinariaceae, and Hygrophoraceae and the genera Clitocybe, Omphalina, and Marasmius. This study provides a framework for future systematics studies in the Agaricales and suggestions for analyzing large molecular data sets.     

Relative efficiencies of the maximum-likelihood, neighbor-joining, and maximum-parsimony methods when substitution rate varies with site

The relative efficiencies of the maximum-likelihood (ML), neighbor- joining (NJ), and maximum-parsimony (MP) methods in obtaining the correct topology and in estimating the branch lengths for the case of four DNA sequences were studied by computer simulation, under the assumption either that there is variation in substitution rate among different nucleotide sites or that there is no variation. For the NJ method, several different distance measures (Jukes-Cantor, Kimura two- parameter, and gamma distances) were used, whereas for the ML method three different transition/transversion ratios (R) were used. For the MP method, both the standard unweighted parsimony and the dynamically weighted parsimony methods were used. The results obtained are as follows: (1) When the R value is high, dynamically weighted parsimony is more efficient than unweighted parsimony in obtaining the correct topology. (2) However, both weighted and unweighted parsimony methods are generally less efficient than the NJ and ML methods even in the case where the MP method gives a consistent tree. (3) When all the assumptions of the ML method are satisfied, this method is slightly more efficient than the NJ method. However, when the assumptions are not satisfied, the NJ method with gamma distances is slightly better in obtaining the correct topology than is the ML method. In general, the two methods show more or less the same performance. The NJ method may give a correct topology even when the distance measures used are not unbiased estimators of nucleotide substitutions. (4) Branch length estimates of a tree with the correct topology are affected more easily than topology by violation of the assumptions of the mathematical model used, for both the ML and the NJ methods. Under certain conditions, branch lengths are seriously overestimated or underestimated. The MP method often gives serious underestimates for certain branches. (5) Distance measures that generate the correct topology, with high probability, do not necessarily give good estimates of branch lengths. (6) The likelihood-ratio test and the confidence-limit test, in Felsenstein's DNAML, for examining the statistical of branch length estimates are quite sensitive to violation of the assumptions and are generally too liberal to be used for actual data. Rzhetsky and Nei's branch length test is less sensitive to violation of the assumptions than is Felsenstein's test. (7) When the extent of sequence divergence is < or = 5% and when > or = 1,000 nucleotides are used, all three methods show essentially the same efficiency in obtaining the correct topology and in estimating branch lengths.(ABSTRACT TRUNCATED AT 400 WORDS)      

呼肠孤病毒科的系统发育分析

运用邻接法(neibor-joining,NJ)、最大似然法(maximum likehood,ML)对国际病毒分类委员会(Intemational Committee on Taxonomy of Viruses,ICTV)第八次报告中所有已报道的呼肠孤病毒科外层衣壳蛋白序列进行了系统发育分析.结果表明:两种方法构建的系统树拓扑结构基本一致,尽管部分分支略有差异,但都能够反映出科内各属的系统发育关系,分析结果支持了该科的分类地位.同时还比较两种建树方法的异同点,并对建树产生的差异原因作了探讨.     

Phylogenetic relationships among six species of Epistylis inferred from 18S-ITS1 sequences

Phylogenetic relationships among six species of Epistylis (i. e. E. plicatilis, E. urceolata, E. chrysemydis, E. hentscheli, E. wenrichi, and E. galea) were investigated using sequences of the first internal transcribed spacer region (ITS-1) of ribosomal DNA (rDNA). Amplified rDNA fragment sequences consisted of 215 or 217 bases of the flanking 18S and 5.8S regions, and the entire ITS-1 region (from 145 to 155 bases). There were more than 33 variable bases between E. galea and the other five species in both the 18S region and the ITS-1 region. The affiliation of them was assessed using Neighbor-joining (NJ), maximum parsimony (MP) and maximum likelihood (ML) analyses. In all the NJ, MP and ML analyses E. galea, whose macronucleic position and shape are distinctly different from those of the other five species, was probably diverged from the ancestor of Epistylis earlier than the other five species. The topology in which E. plicatilis and E. hentscheli formed a strongly supported sister clade to E. urceol     

Phylogenetic relationships among the Lorisoidea as indicated by craniodental morphology and mitochondrial sequence data

The phylogeny of the Afro-Asian Lorisoidea is controversial. While postcranial data attest strongly to the monophyly of the Lorisidae, most molecular analyses portray them as paraphyletic and group the Galagidae alternately with the Asian or African lorisids. One of the problems that has bedevilled phylogenetic analysis of the group in the past is the limited number of taxa sampled for both ingroup families. We present the results of a series of phylogenetic analyses based on 635 base pairs (bp) from two mitochondrial genes (12S and 16S rRNA) with and without 36 craniodental characters, for 11 galagid and five lorisid taxa. The outgroup was the gray mouse lemur (Microcebus murinus). Analyses of the molecular data included maximum parsimony (MP), neighbor joining (NJ), maximum likelihood (ML), and Bayesian methods. The model-based analyses and the combined "molecules+morphology" analyses supported monophyly of the Lorisidae and Galagidae. The lorisids form two geographically defined clades. We find no support for the taxonomy of Galagidae as proposed recently by Groves [Primate Taxonomy, Washington, DC: Smithsonian Institution Press. 350 p, 2001]. The taxonomy of Nash et al. [International Journal of Primatology 10:57-80, 1989] is supported by the combined "molecules+morphology" analysis; however, the model-based analyses suggest that Galagoides may be an assemblage of species united by plesiomorphic craniodental characters.     

介形纲丽足目和速足目及相关类群18S rDNA分子系统发育的研究

文中分析现生介形类 (Ostracoda) 4目 2 1科 2 9属的 18SrDNA部分序列 ,采用最大似然法 (ML)、邻接法 (NJ)和最大简约法 (MP) ,尝试构建介形类的分子系统树 ;结合介形类的形态特征和化石记录 ,主要对速足目(Podocopida)、丽足目 (Myodocopida)及其超科级分类阶元的系统发生关系进行探讨。 3种分析方法均支持形态学上Podocopida ,Myodocopida和海萤超科 (Cypridinacea)的界定 ;但对Podocopida目土菱介超科 (Bairdiacea)的系统地位提出质疑 ,该类群可能不是单系发生的自然类群。上述分析显示 ,Podocopida,Myodocopida,Platycopida和Halo cypridina组成一个单系群 ;介形类在目、超科、科和属的水平上可能发生过多次辐射分化     

REEXAMINATION OF PHYLOGENETIC RELATIONSHIPS WITHIN THE COLONIAL VOLVOCALES (CHLOROPHYTA): AN ANALYSIS OF atpB AND rbcL GENE SEQUENCES

The chloroplast-encoded atp B gene was sequenced from 33 strains representing 28 species of the colonial Volvocales (the Volvocaceae and its relatives) to reexamine phylogenetic relationships as previously deduced by morphological data and rbc L gene sequence data.1128 base pairs in the coding regions of the atp B gene were analyzed by MP, NJ, and ML analyses. Although supported with relatively low bootstrap values (75% and 65% in the NJ and ML analyses, respectively), three anisogamous/oogamous volvocacean genera— Eudorina, Pleodorina, and Volvox, excluding the section Volvox (= Euvolvox, illegitimate name), constituted a large monophyletic group (Eudorina group). Outside the Eudorina group, a robust lineage composed of three species of Volvox sect. Volvox was resolved as in the rbc L gene trees, rejecting the hypothesis of the previous cladistic analysis based on morphological data that the genus Volvox is monophyletic. In addition, the NJ and ML trees suggested that Eudorina is a nonmonophyletic genus as inferred from the morphological data and rbc L gene sequences. Although phylogenetic status of the genus Gonium is ambiguous in the rbc L gene trees and the paraphyly of this genus is resolved in the cladistic analysis based on morphological data, the atp B gene sequence data suggest monophyly of Gonium with relatively low bootstrap values (56–61%) in the NJ and ML trees. On the basis of the combined sequence data (2256 base pairs) from atp B and rbc L genes, Gonium was resolved as a robust monophyletic genus in the NJ and ML trees (with 68–86% bootstrap values), and Eudorina elegans Ehrenberg represented a paraphyletic species positioned most basally within the Eudorina group. However, phylogenetic status and relationships of the families of the colonial Volvocales were still almost ambiguous even in the combined analysis.     

A heuristic approach of maximum likelihood method for inferring phylogenetic tree and an application to the mammalian SOX-3 origin of the testis-determining gene SRY

Applying the tree bisection and reconnection (TBR) algorithm, we have developed a heuristic method (maximum likelihood (ML)-TBR) for inferring the ML tree based on tree topology search. For initial trees from which iterative processes start in ML-TBR, two cases were considered: one is 100 neighbor-joining (NJ) trees based on the bootstrap resampling and the other is 100 randomly generated trees. The same ML tree was obtained in both cases. All different iterative processes started from 100 independent initial trees ultimately converged on one optimum tree with the largest log-likelihood value, suggesting that a limited number of initial trees will be quite enough in ML-TBR. This also suggests that the optimum tree corresponds to the global optimum in tree topology space and thus probably coincides with the ML tree inferred by intact ML analysis. This method has been applied to the inference of phylogenetic tree of the SOX family members. The mammalian testis-determining gene SRY is believed to have evolved from SOX-3, a member of the SOX family, based on several lines of evidence, including their sequence similarity, the location of SOX-3 on the X chromosome and some aspects of their expression. This model should be supported directly from the phylogenetic tree of the SOX family, but no evidence has been provided to date. A recently published NJ tree shows implausibly remote origin of SRY, suggesting that a more sophisticated method is required for understanding this problem. The ML tree inferred by the present method showed that the SRYs of marsupial and placental mammals form a monophyletic cluster which had diverged from the mammalian SOX-3 in the early evolution of mammals.     

MOLECULAR SYSTEMATICS OF RIVER DOLPHINS INFERRED FROM COMPLETE MITOCHONDRIAL CYTOCHROME-B GENE SEQUENCES

1,140 bp of the complete mitochondrial cytochrome- b gene sequences of baiji ( Lipotes vexillifer ), franciscana ( Pontoporia blainvillei ), and Ganges river dolphin ( Platanista gangetica gangetica ) were determined to address the systematic position and phylogeny of extant river dolphins with combination of homologous sequences of other cetaceans. The neighbor-joining (NJ), maximum parsimony (MP), and maximum likelihood (ML) phylogenetic analyses all identified the river dolphins into three lineages, i. e., Platanista, Lipotes , and Inia + Pontoporia . The Lipotes did not have sister relationship with either Platanista or Inia + Pontoporia , which strongly supported the referral of Lipotes to a separate family, i. e. , Lipotidae. There were very high sequence divergences between all river dolphin genera, suggesting a relatively longer period of separation time than those among other odontocete families.     

从核rDNA序列探讨隙蛛亚科Coelotinae的分类地位

隙蛛亚科Coelotinae主要分布于东亚地区,其中我国的已有种类占到全世界种数的一半以上,因此对于我国隙蛛类蜘蛛的研究已经成为世界暗蛛科研究的重点之一。隙蛛亚科属于无筛器类群,于1893年,由Cambridge以隙蛛属为模式属而建立,归属于无筛器的漏斗蛛科。之后,虽然经历了数次修订     

Phylogenetic relationships among six species of Epistylis inferred from 18S-ITS1 sequences

Phylogenetic relationships among six species of Epistylis (i. e. E. plicatilis, E. urceolata, E. chrysemydis, E. hentscheli, E. wenrichi, and E. galea) were investigated using sequences of the first internal transcribed spacer region (ITS-1) of ribosomal DNA (rDNA). Amplified rDNA fragment sequences consisted of 215 or 217 bases of the flanking 18S and 5.8S regions, and the entire ITS-1 region (from 145 to 155 bases). There were more than 33 variable bases between E. galea and the other five species in both the 18S region and the ITS-1 region. The affiliation of them was assessed using Neighbor-joining (NJ), maximum parsimony (MP) and maximum likelihood (ML) analyses. In all the NJ, MP and ML analyses E. galea, whose macronucleic position and shape are distinctly different from those of the other five species, was probably diverged from the ancestor of Epistylis earlier than the other five species. The topology in which E. plicatilis and E. hentscheli formed a strongly supported sister clade to E. urceolata, E. chrysemydis, and E. wenrichi was consistent with variations in the thickness of the peristomial lip. We concluded that the macronucleus and peristomial lip might be the important phylogenetic characteristics within the genus Epistylis.     

Long-Branch Attraction Bias and Inconsistency in Bayesian Phylogenetics

Bayesian inference (BI) of phylogenetic relationships uses the same probabilistic models of evolution as its precursor maximum likelihood (ML), so BI has generally been assumed to share ML''s desirable statistical properties, such as largely unbiased inference of topology given an accurate model and increasingly reliable inferences as the amount of data increases. Here we show that BI, unlike ML, is biased in favor of topologies that group long branches together, even when the true model and prior distributions of evolutionary parameters over a group of phylogenies are known. Using experimental simulation studies and numerical and mathematical analyses, we show that this bias becomes more severe as more data are analyzed, causing BI to infer an incorrect tree as the maximum a posteriori phylogeny with asymptotically high support as sequence length approaches infinity. BI''s long branch attraction bias is relatively weak when the true model is simple but becomes pronounced when sequence sites evolve heterogeneously, even when this complexity is incorporated in the model. This bias—which is apparent under both controlled simulation conditions and in analyses of empirical sequence data—also makes BI less efficient and less robust to the use of an incorrect evolutionary model than ML. Surprisingly, BI''s bias is caused by one of the method''s stated advantages—that it incorporates uncertainty about branch lengths by integrating over a distribution of possible values instead of estimating them from the data, as ML does. Our findings suggest that trees inferred using BI should be interpreted with caution and that ML may be a more reliable framework for modern phylogenetic analysis.     

Monophyly of the family Desmoscolecidae (Nematoda,Demoscolecida) and Its Phylogenetic position inferred from 18S rDNA sequences

To infer the monophyletic origin and phylogenetic relationships of the order Desmoscolecida, a unique and puzzling group of mainly free-living marine nematodes, we newly determined nearly complete 18S rDNA sequences for six marine desmoscolecid nematodes belonging to four genera (Desmoscolex, Greeffiella, Tricoma and Paratricoma). Based on the present data and those of 72 nematode species previously reported, the first molecular phylogenetic analysis focusing on Desmoscolecida was done by using neighbor joining (NJ), maximum parsimony (MP), maximum likelihood (ML) and Bayesian inference (BI) methods. All four resultant trees consistently and strongly supported that the family Desmoscolecidae forms a monophyletic group with very high node confidence values. The monophyletic clade of desmocolecid nematodes was placed as a sister group of the clade including some members of Monhysterida and Araeolaimida, Cyartonema elegans (Cyartonematidae) and Terschellingia longicaudata (Linhomoeidae) in all the analyses. However, the present phylogenetic trees do not show any direct attraction between the families Desmoscolecidae and Cyartonematidae. Within the monophyletic clade of the family Desmoscolecidae in all of the present phylogenetic trees, there were consistently observed two distinct sub-groups which correspond to the subfamilies Desmoscolecinae [Greeffiella sp. + Desmoscolex sp.] and Tricominae [Paratricoma sp. + Tricoma sp].     

Molecular phylogeny of oribatid mites (Oribatida, Acari): evidence for multiple radiations of parthenogenetic lineages

Nucleotide sequences of the D3 expansion segment and its flanking regions of the 28S rDNA gene were used to evaluate phylogenetic relationships among representative sexual and asexual oribatid mites (Oribatida, Acariformes). The aim of this study was to investigate the hypothesis that oribatid mites consist of species-rich clusters of asexual species that may have radiated while being parthenogenetic. Furthermore, the systematic position of the astigmate mites (Astigmata, Acariformes) which have been hypothesised to represent a paedomorphic lineage within the oribatid mites, is investigated. This is the first phylogenetic tree for oribatid mites s.l. (incl. Astigmata) based on nucleotide sequences. Intraspecific genetic variation in the D3 region was very low, confirming the hypothesis that this region is a good species marker. Results from neighbour joining (NJ) and maximum parsimony (MP) algorithms indicate that several species-rich parthenogenetic groups like Camisiidae, Nanhermanniidae and Malaconothridae are monophyletic, consistent with the hypothesis that some oribatid mite groups diversified despite being parthenogenetic. The MP and maximum likelihood (ML) method indicated that the D3 region is a good tool for elucidating the relationship of oribatid mite species on a small scale(genera, families) but is not reliable for large-scale taxonomy, because branches from the NJ algorithm collapsed in the MP and ML tree. In all trees calculated by different algorithms the Astigmata clustered within the oribatid mites, as proposed earlier.     

Efficiencies of different genes and different tree-building methods in recovering a known vertebrate phylogeny

The relative efficiencies of different protein-coding genes of the mitochondrial genome and different tree-building methods in recovering a known vertebrate phylogeny (two whale species, cow, rat, mouse, opossum, chicken, frog, and three bony fish species) was evaluated. The tree-building methods examined were the neighbor joining (NJ), minimum evolution (ME), maximum parsimony (MP), and maximum likelihood (ML), and both nucleotide sequences and deduced amino acid sequences were analyzed. Generally speaking, amino acid sequences were better than nucleotide sequences in obtaining the true tree (topology) or trees close to the true tree. However, when only first and second codon positions data were used, nucleotide sequences produced reasonably good trees. Among the 13 genes examined, Nd5 produced the true tree in all tree-building methods or algorithms for both amino acid and nucleotide sequence data. Genes Cytb and Nd4 also produced the correct tree in most tree-building algorithms when amino acid sequence data were used. By contrast, Co2, Nd1, and Nd41 showed a poor performance. In general, large genes produced better results, and when the entire set of genes was used, all tree-building methods generated the true tree. In each tree-building method, several distance measures or algorithms were used, but all these distance measures or algorithms produced essentially the same results. The ME method, in which many different topologies are examined, was no better than the NJ method, which generates a single final tree. Similarly, an ML method, in which many topologies are examined, was no better than the ML star decomposition algorithm that generates a single final tree. In ML the best substitution model chosen by using the Akaike information criterion produced no better results than simpler substitution models. These results question the utility of the currently used optimization principles in phylogenetic construction. Relatively simple methods such as the NJ and ML star decomposition algorithms seem to produce as good results as those obtained by more sophisticated methods. The efficiencies of the NJ, ME, MP, and ML methods in obtaining the correct tree were nearly the same when amino acid sequence data were used. The most important factor in constructing reliable phylogenetic trees seems to be the number of amino acids or nucleotides used.      

基于COI基因序列的九种真蝽属昆虫的分子系统学研究(半翅目:蝽科)

通过对真蝽属Pentatoma 9种昆虫线粒体COI基因约798bp的序列进行分子进化分析,并以同蝽科宽铗同蝽Acanthosoma labiduroides为外群,采用最大简约法、最大似然法和邻接法构建了分子系统树,来探讨真蝽属的系统发育关系.研究结果支持褐真蝽群P. semiannulata-group的划分,绿角真蝽Pentatoma viridicornuta应划归到褐真蝽群P. Semiannulata-group;红足真蝽群中的角肩真蝽P. angulata与红足真蝽P. rufipes遗传距离较小,它们是否为1个物种值得关注;真蝽属各群之间的系统发育关系以及是否可分为3个属值得进一步研究.     

Relationships within Cornales and circumscription of Cornaceae-matK and rbcL sequence data and effects of outgroups and long branches

Phylogenetic relationships in Cornales were assessed using sequences rbcL and matK. Various combinations of outgroups were assessed for their suitability and the effects of long branches and outgroups on tree topology were examined using RASA 2.4 prior to conducting phylogenetic analyses. RASA identified several potentially problematic taxa having long branches in individual data sets that may have obscured phylogenetic signal, but when data sets were combined RASA no longer detected long branch problems. t(RASA) provides a more conservative measurement for phylogenetic signal than the PTP and skewness tests. The separate matK and rbcL sequence data sets were measured as not containing phylogenetic signal by RASA, but PTP and skewness tests suggested the reverse [corrected]. Nonetheless, the matK and rbcL sequence data sets suggested relationships within Cornales largely congruent with those suggested by the combined matK-rbcL sequence data set that contains significant phylogenetic signal as measured by t(RASA), PTP, and skewness tests. Our analyses also showed that a taxon having a long branch on the tree may not be identified as a "long-branched" taxon by RASA. The long branches identified by RASA had little effect on the arrangement of other taxa in the tree, but the placements of the long-branched taxa themselves were often problematic. Removing the long-branched taxa from analyses generally increased bootstrap support, often substantially. Use of non-optimal outgroups (as identified by RASA) decreased phylogenetic resolution in parsimony analyses and suggested different relationships in maximum likelihood analyses, although usually weakly supported clades (less than 50% support) were impacted. Our results do not recommend using t(RASA) as a sole criterion to discard data or taxa in phylogenetic analyses, but t(RASA) and the taxon variance ratio obtained from RASA may be useful as a guide for improved phylogenetic analyses. Results of parsimony and ML analyses of the sequence data using optimal outgroups suggested by RASA revealed four major clades within Cornales: (1) Curtisia-Grubbia, (2) Cornus-Alangium, (3) Nyssa-Camptotheca-Davidia-Mastixia-Diplopanax, and (4) Hydrangeaceae-Loasaceae, with clades (2) and (3) forming a monophyletic group sister to clade (4) and clade (1) sister to the remainder of Cornales. However, there was not strong bootstrap support for relationships among the major clades. The placement of Hydrostachys could not be reliably determined, although most analyses place the genus within Hydrangeaceae; ML analyses, for example, placed the genus as the sister of Hydrangeeae. Our results supported a Cornales including the systematically problematic Hydrostachys, a Cornaceae consisting of Cornus and Alangium, a Nyssaceae consisting of Nyssa and Camptotheca, a monogeneric Davidiaceae, a Mastixiaceae consisting of Mastixia and Diplopanax, and an expanded Grubbiaceae consisting of Grubbia and Curtisia, and two larger families, Hydrangeaceae and Loasaceae.     

Phylogenetic relationships between the six superoxide dismutase proteins (FeSOD) of Trichomonas vaginalis and FeSOD6 genetic diversity

The parasitic protozoan Trichomonas vaginalis is known to contain several types of Fe-containing superoxide dismutase proteins (FeSOD). Using three different methods of phylogenetic analysis, maximum parsimony (MP), neighbor joining (NJ), and maximum likelihood (ML) methods, we examined the phylogenetic relationships among the six FeSOD (FeSOD1-FeSOD6) based on their amino acid sequences. All the analyses consistently suggested that the six proteins formed a monophyletic group implying that they probably be originated from an ancestral protein form through repeated duplication events. Although MP tree was totally unresolved, the NJ and ML trees revealed that FeSOD6 placed the most basal position and thus emerged earlier than the other five gene types during the evolution of T. vaginalis. Phylogenetic relationships among the five remaining proteins were (FeSOD2, FeSOD3), (FeSOD4, (FeSOD1, FeSOD5)) although weakly supported in terms of bootstrapping values. In addition to this, we newly designed two PCR primer specifically amplifying full-length FeSOD6 gene and examined its genetic diversity among 12 T. vaginalis isolates from five countries and three continents. They had the same nucleotide sequences except those of three Korean isolates which showed one to three different nucleotides.     

Phylogeny and comparative genomic analysis of Pteriomorphia (Mollusca: Bivalvia) based on complete mitochondrial genomes

The subclass Pteriomorphia is a morphologically diverse and economically important group of Mollusca. We retrieved 42 mitochondrial genomes (mtGenomes) of Pteriomorphia and concatenated protein-coding genes, rRNAs and tRNAs to assess phylogenetic relationships and divergence times among the families with maximum likelihood (ML) and Bayesian inference (BI) analyses. Both ML and BI analyses strongly support the same topology except for the position of Atrina pectinata. Our study confirms the monophyly of the families Arcidae, Mytilidae, Pteriidae, Ostreidae and Pectinidae. Within Pteriomorphia, we recovered two clusters, one comprising Mytilidae, Arcidae and Pectinidae, the other consisting of Ostreidae, Pteriidae and Pinnidae, but we did not confirm a basal position for any family. The phylogenetic trees suggest that Ostreidae, Pteriidae and Pinnidae should be grouped as the order Ostreoida. Divergence times of major families are estimated as follows: Arcidae, 315.9 Ma; Pectinidae, 384.4 Ma; Ostreidae, 240.8 Ma; Mytilidae, 390.8 Ma. Comparative analysis indicates a low-level codon usage bias (with an average of 50.29) in mtGenomes of Pteriomorphia. In Mytilidae and Ostreidae, the codon usage bias was under mutation pressure rather than selection. Contrastingly, mutation is not the main factor in defining the codon usage in Pectinidae and Pteriidae. Among Ostreidae, Pectinidae and Mytilidae, Ka/Ks ratios range from 0.00 to 1.22 and most values (89.11%) are less than 0.20, indicating that most genes are under strong negative or purifying selection. The protein-coding gene orders show dramatically different patterns in Pteriomorphia. There is no gene block even consisting of two genes that is shared by five families.     

Selective constraints on amino acids estimated by a mechanistic codon substitution model with multiple nucleotide changes

Background

Empirical substitution matrices represent the average tendencies of substitutions over various protein families by sacrificing gene-level resolution. We develop a codon-based model, in which mutational tendencies of codon, a genetic code, and the strength of selective constraints against amino acid replacements can be tailored to a given gene. First, selective constraints averaged over proteins are estimated by maximizing the likelihood of each 1-PAM matrix of empirical amino acid (JTT, WAG, and LG) and codon (KHG) substitution matrices. Then, selective constraints specific to given proteins are approximated as a linear function of those estimated from the empirical substitution matrices.

Results

Akaike information criterion (AIC) values indicate that a model allowing multiple nucleotide changes fits the empirical substitution matrices significantly better. Also, the ML estimates of transition-transversion bias obtained from these empirical matrices are not so large as previously estimated. The selective constraints are characteristic of proteins rather than species. However, their relative strengths among amino acid pairs can be approximated not to depend very much on protein families but amino acid pairs, because the present model, in which selective constraints are approximated to be a linear function of those estimated from the JTT/WAG/LG/KHG matrices, can provide a good fit to other empirical substitution matrices including cpREV for chloroplast proteins and mtREV for vertebrate mitochondrial proteins.

Conclusions/Significance

The present codon-based model with the ML estimates of selective constraints and with adjustable mutation rates of nucleotide would be useful as a simple substitution model in ML and Bayesian inferences of molecular phylogenetic trees, and enables us to obtain biologically meaningful information at both nucleotide and amino acid levels from codon and protein sequences.