Comparison of the phylogenetic performance of neodermatan mitochondrial protein-coding genes

We examined the ability of all protein-coding genes of 11 neodermatan mitochondrial genomes to recover the phylogeny of their combination. We subsampled randomly selected lengths (1–6000 bp) of contiguous sequence to consider the effect of length (100–6000 bp) on performance. All genes performed as well as randomly selected regions except ND2 (better) and ND5 (worse). All expected nodes were recovered (> 90% bootstrap) when 4.0 Kb had been sampled; 40 samples of 100 performed as well as four samples of 1000.     

Comparative analysis of sequences preceding protein-coding mitochondrial genes in flowering plants

We examined the nucleotide sequences preceding 23 mitochondrial protein-coding genes held in common by maize, rice, wheat, sugar beet, tobacco, Arabidopsis, and Brassica to look for features related to translation initiation and to assess the degree of conservation in mitochondrial mRNA leaders among these plants. We observed broad variation in sequence similarity as illustrated by dot plot analysis, ranging from a level rivaling that of coding sequences to complete absence of homology due to lineage-specific DNA rearrangements. Genes encoding ATP synthase subunits predominated in the latter category, whereas ones encoding cytochrome c biogenesis proteins and NADH dehydrogenase subunits were primarily of the highly conserved type. Within the region immediately preceding initiation codons, in most cases we did not observe motifs consistent with a bacterial-type Shine-Dalgarno interaction to assist in ribosome binding, nor was any other consensus sequence evident. In fact, indels in the form of tandem repeats were seen among homologues from different plants. We did, however, observe a bias for high adenosine and low cytosine in the proximal approximately 30 nt compared with further upstream. Duplicates of some sequences in our data set were found to be associated with more than one gene within a genome. Indeed, 3 such families of upstream cassettes were identified, and they exhibit a lineage-specific distribution among plants. Moreover, the presence of related sequences at genomic sites distant from known genes raises the possibility of future recruitment as regulatory elements. Our observations point to a dynamic nature in the makeup of the 5' leaders of plant mitochondrial mRNAs and an apparent plasticity in translational control elements.     

Evaluating nuclear protein-coding genes for phylogenetic utility in beetles

Although nuclear protein-coding genes have proven broadly useful for phylogenetic inference, relatively few such genes are regularly employed in studies of Coleoptera, the most diverse insect order. We increase the number of loci available for beetle systematics by developing protocols for three genes previously unused in beetles (alpha-spectrin, RNA polymerase II and topoisomerase I) and by refining protocols for five genes already in use (arginine kinase, CAD, enolase, PEPCK and wingless). We evaluate the phylogenetic performance of each gene in a Bayesian framework against a presumably known test phylogeny. The test phylogeny covers 31 beetle specimens and two outgroup taxa of varying age, including three of the four extant beetle suborders and a denser sampling in Adephaga and in the carabid genus Bembidion. All eight genes perform well for Cenozoic divergences and accurately separate closely related species within Bembidion, but individual genes differ markedly in accuracy over the older Mesozoic and Permian divergences. The concatenated data reconstruct the test phylogeny with high support in both Bayesian and parsimony analyses, indicating that combining data from multiple nuclear loci will be a fruitful approach for assembling the beetle tree of life.     

Housekeeping genes for phylogenetic analysis of eutherian relationships

The molecular relationship of placental mammals has attracted great interest in recent years. However, 2 crucial and conflicting hypotheses remain, one with respect to the position of the root of the eutherian tree and the other the relationship between the orders Rodentia, Lagomorpha (rabbits, hares), and Primates. Although most mitochondrial (mt) analyses have suggested that rodents have a basal position in the eutherian tree, some nuclear data in combination with mt-rRNA genes have placed the root on the so-called African clade or on a branch that includes this clade and the Xenarthra (e.g., anteater and armadillo). In order to generate a new and independent set of molecular data for phylogenetic analysis, we have established cDNA sequences from different tissues of various mammalian species. With this in mind, we have identified and sequenced 8 housekeeping genes with moderately fast rate of evolution from 22 placental mammals, representing 11 orders. In order to determine the root of the eutherian tree, the same genes were also sequenced for 3 marsupial species, which were used as outgroup. Inconsistent with the analyses of nuclear + mt-rRNA gene data, the current data set did not favor a basal position of the African clade or Xenarthra in the eutherian tree. Similarly, by joining rodents and lagomorphs on the same basal branch (Glires hypothesis), the data set is also inconsistent with the tree commonly favored in mtDNA analyses. The analyses of the currently established sequences have helped examination of problematic parts in the eutherian tree at the same time as they caution against suggestions that have claimed that basal eutherian relationships have been conclusively settled.     

A phylogeny of the Caniformia (order Carnivora) based on 12 complete protein-coding mitochondrial genes

Evolutionary relationships of the order Carnivora have been extensively studied. However, phylogenetic studies based on different types of data, species samples, and methods of analysis provide contradictory results. Consequently, phylogenetic relationships of Carnivora remain contentious. Here, the sequence of 12 mitochondrial genes (10,842 nucleotides) from a total of 38 carnivore species was used to investigate the phylogeny of the caniform (dog-like) carnivores. An analysis using maximum parsimony, maximum likelihood, and Bayesian approaches provided a unique and well-supported solution to most contentious relationships within Caniformia. The clade Arctoidea was shown to consist of three major monophyletic groups: Pinnipedia, Ursidae, and Musteloidea. Within Pinnipedia, the families Otariidae and Odobenidae formed a clade, sister to Phocidae. Within Musteloidea, there was a sister relationship between true mustelids (i.e., excluding the skunks) and procyonids, and between ailurids and mephitids (skunks). Despite a high level of confidence obtained at most nodes, uncertainty remained about the relative position of the three major arctoid clades.     

A molecular phylogenetic analysis of the Octocorallia (Cnidaria: Anthozoa) based on mitochondrial protein-coding sequences

Despite their abundance and ecological importance in a wide variety of shallow and deep water marine communities, octocorals (soft corals, sea fans, and sea pens) are a group whose taxonomy and phylogenetic relationships remain poorly known and little studied. The group is currently divided into three orders (O: Alcyonacea, Pennatulacea, and Helioporacea); the large O. Alcyonacea (soft corals and sea fans) is further subdivided into six sub-ordinal groups on the basis of skeletal composition and colony growth form. We used 1429bp of two mitochondrial protein-coding genes, ND2 and msh1, to construct a phylogeny for 103 octocoral genera representing 28 families. In agreement with a previous 18S rDNA phylogeny, our results support a division of Octocorallia into two major clades plus a third, minor clade. We found one large clade (Holaxonia-Alcyoniina) comprising the sea fan sub-order Holaxonia and the majority of soft corals, and a second clade (Calcaxonia-Pennatulacea) comprising sea pens (O. Pennatulacea) and the sea fan sub-order Calcaxonia. Taxa belonging to the sea fan group Scleraxonia and the soft coral family Alcyoniidae were divided among the Holaxonia-Alcyoniina clade and a third, small clade (Anthomastus-Corallium) whose relationship to the two major clades was unresolved. In contrast to the previous studies, we found sea pens to be monophyletic but nested within Calcaxonia; our analyses support the sea fan family Ellisellidae as the sister taxon to the sea pens. We are unable to reject the hypothesis that the calcaxonian and holaxonian skeletal axes each arose once and suggest that the skeletal axis of sea pens is derived from that of Calcaxonia. Topology tests rejected the monophyly of sub-ordinal groups Alcyoniina, Scleraxonia, and Stolonifera, as well as 9 of 14 families for which we sampled multiple genera. The much broader taxon sampling and better phylogenetic resolution afforded by our study relative to the previous efforts greatly clarify the relationships among families and sub-ordinal groups within each of the major clades. The failure of these mitochondrial genes as well as previous 18S rDNA studies to resolve many of the deeper nodes within the tree (including its root) suggest that octocorals underwent a rapid radiation and that large amounts of sequence data will be required in order to resolve the basal relationships within the clade.     

A phylogenetic analysis of the emberizid sparrows based on three mitochondrial genes

Previous molecular phylogenetic studies have examined the taxonomic relationships among a number of typical emberizid sparrow genera. To help clarify these relationships, we sequenced a 1673 base pair fragment for the complete sequence of three mitochondrial genes: adenosine triphosphatase (Atp8 and Atp6) and cytochrome oxidase subunit III (COIII) for 38 sparrow species, along with Passerina amoena (Cardinalidae) and Piranga ludoviciana (Thraupidae) which were selected as the outgroups. Our analysis confirms the monophyly of traditional genera such as Junco, Melospiza, and Zonotrichia. Although Calcarius and Plectrophenax are often thought to be putative emberizids, all our analyses placed these genera basal to all other sparrows examined. As observed with Calcarius, Spizella did not form a monophyletic group, with S. arborea being the sister-taxon to Passerella iliaca. Our analyses also suggest that Aimophila ruficeps is probably more closely related to the "brown towhees" (Pipilo aberti, P. crissalis, and P. fuscus) than its putative congeners. The genus Ammodramus was also not monophyletic, since it appears that Passerculus sandwichensis is more closely related to A. henslowii and A. leconteii then either one is related to its congener A. savannarum. Finally, our analyses exhibited other unsuspected associations, such as the sister-taxon relationships between Amphispiza bilineata and the Chondestes grammacus/Calamospiza melanocorys clade, and Amphispiza belli and Pooecetes gramineus.     

普通野生稻线粒体蛋白质编码基因密码子使用偏好性的分析

以普通野生稻（Oryza rufipogon Griff.）线粒体基因组为对象,分析其蛋白质编码基因的密码子使用特征及与亚洲栽培稻（O.sativa L.）的差异,探讨其密码子偏性形成的影响因素和进化过程。结果显示：普通野生稻线粒体基因组编码序列第1、第2和第3位碱基的GC含量依次为49.18%、42.67%和40.86%;有效密码子数（Nc）分布于45.32~61.00之间,其密码子偏性较弱;Nc值仅与GC₃呈显著相关,密码子第3位的碱基组成对密码子偏性影响较大;第1向量轴上显示9.91%的差异,其与GC_3s、Nc、密码子偏好指数（CBI）和最优密码子使用频率（Fop）的相关性均达到显著水平;而GC₃和GC₁₂的相关性未达到显著水平。因此,普通野生稻线粒体基因组密码子的使用偏性主要受自然选择压力影响而形成。本研究确定了21个普通野生稻线粒体基因组的最优密码子,大多以A或T结尾,与叶绿体密码子具有趋同进化,但是与核基因组具有不同的偏好性。同义密码子相对使用度（RSCU）、PR2偏倚分析和中性绘图分析显示,普通野生稻线粒体基因功能和其密码子使用密切相关,且线粒体密码子使用在普通野生稻、粳稻（O.sativa L.subsp.japonica Kato）和籼稻（O.sativa L.subsp.indica Kato）内具有同质性。     

Choosing appropriate substitution models for the phylogenetic analysis of protein-coding sequences

Although phylogenetic inference of protein-coding sequences continues to dominate the literature, few analyses incorporate evolutionary models that consider the genetic code. This problem is exacerbated by the exclusion of codon-based models from commonly employed model selection techniques, presumably due to the computational cost associated with codon models. We investigated an efficient alternative to standard nucleotide substitution models, in which codon position (CP) is incorporated into the model. We determined the most appropriate model for alignments of 177 RNA virus genes and 106 yeast genes, using 11 substitution models including one codon model and four CP models. The majority of analyzed gene alignments are best described by CP substitution models, rather than by standard nucleotide models, and without the computational cost of full codon models. These results have significant implications for phylogenetic inference of coding sequences as they make it clear that substitution models incorporating CPs not only are a computationally realistic alternative to standard models but may also frequently be statistically superior.     

普通野生稻线粒体蛋白质编码基因密码子使用偏好性的分析

以普通野生稻(Oryza rufipogon Griff.)线粒体基因组为对象,分析其蛋白质编码基因的密码子使用特征及与亚洲栽培稻(O. sativa L.)的差异,探讨其密码子偏性形成的影响因素和进化过程。结果显示:普通野生稻线粒体基因组编码序列第1、第2和第3位碱基的GC含量依次为49.18%、42.67%和40.86%;有效密码子数(Nc)分布于45.32～61.00之间,其密码子偏性较弱; Nc值仅与GC_3呈显著相关,密码子第3位的碱基组成对密码子偏性影响较大;第1向量轴上显示9.91%的差异,其与GC_3s、Nc、密码子偏好指数(CBI)和最优密码子使用频率(Fop)的相关性均达到显著水平;而GC_3和GC₁₂的相关性未达到显著水平。因此,普通野生稻线粒体基因组密码子的使用偏性主要受自然选择压力影响而形成。本研究确定了21个普通野生稻线粒体基因组的最优密码子,大多以A或T结尾,与叶绿体密码子具有趋同进化,但是与核基因组具有不同的偏好性。同义密码子相对使用度(RSCU)、PR2偏倚分析和中性绘图分析显示,普通野生稻线粒体基因功能和其密码子使用密切相关,且线粒体密码子使用在普通野生稻、粳稻(O. sativa L. subsp. japonica Kato)和籼稻(O. sativa L. subsp.indica Kato)内具有同质性。     

Evaluating phylogenetic informativeness and data-type usage for new protein-coding genes across Vertebrata

As a resource for vertebrate phylogenetics, we developed 75 new protein-coding genes using a combination of expressed sequence tags (ESTs) available in Genbank, and targeted amplification of complementary DNA (cDNA). In addition, we performed three additional analyses in order to assess the utility of our approach. First, we profiled the phylogenetic informativeness of these new markers using the online program PhyDesign. Next, we compared the utility of four different data-types used in phylogenetics: nucleotides (NUCL), amino acids (AA), 1st and 2nd codon positions only (N12), and modified sequences to account for codon degeneracy (DEGEN1; Regier et al., 2010). Lastly, we use these new markers to construct a vertebrate phylogeny and address the uncertain relationship between higher-level mammal groups: monotremes, marsupials, and placentals. Our results show that phylogenetic informativeness of the 75 new markers varies, both in the amount of phylogenetic signal and optimal timescale. When comparing the four data-types, we find that the NUCL data-type, due to the high level of phylogenetic signal, performs the best across all divergence times. The remaining three data-types (AA, N12, DEGEN1) are less subject to homoplasy, but have greatly reduced levels of phylogenetic signal relative to NUCL. Our phylogenetic inference supports the Theria hypothesis of mammalian relationships, with marsupials and placentals being sister groups.     

Molecular phylogenetic relationships of pond frogs distributed in the Palearctic region inferred from DNA sequences of mitochondrial 12S ribosomal RNA and cytochrome b genes

The evolutionary relationships of pond frogs distributed in the Far East and Europe were investigated by analyses of nucleotide sequences of mitochondrial 12S ribosomal RNA (12S rRNA) and cytochrome b (cyt b) genes. The nucleotide sequences of a 412-bp segment of the 12S rRNA gene and a 534-bp segment of the cyt b gene were determined by the PCR-direct sequencing method using 19 frogs belonging to six species and one subspecies distributed in the Palearctic region. Phylogenetic trees were constructed by the neighbor-joining and maximum-likelihood methods using Rana catesbeiana or Xenopus laevis as an outgroup. The 412-bp segment of the 12S rRNA gene contained 65 variable sites including gap sites, and the 534-bp segment of the cyt b gene contained 160 variable sites. The nucleotide sequence divergences of the 12S rRNA gene were 0.25-4.83% within the Far Eastern frogs, 0.25-6.22% within the European frogs, and 8.74-11.24% between the Far Eastern and the European frogs, whereas those of the cyt b gene were 3.64-14.73% within the Far Eastern frogs, 0.38-14.42% within the European frogs, and 16.53-23.58% between the Far Eastern and the European frogs. Although most nucleotide substitutions were at the third codon position of the cyt b gene and were silent mutations, 4 amino acid replacements occurred within the Far Eastern frogs, 4 within the European frogs, and 11 between the Far Eastern and the European frogs. The phylogenetic trees constructed from the nucleotide sequence divergences showed slightly different topologies for the 12S rRNA and cyt b genes. R. esculenta from Ukraine was closely related to R. lessonae from Luxembourg in both the 12S rRNA and the cyt b gene sequences.     

Neogastropod phylogenetic relationships based on entire mitochondrial genomes

Background  

The Neogastropoda is a highly diversified group of predatory marine snails (Gastropoda: Caenogastropoda). Traditionally, its monophyly has been widely accepted based on several morphological synapomorphies mostly related with the digestive system. However, recent molecular phylogenetic studies challenged the monophyly of Neogastropoda due to the inclusion of representatives of other caenogastropod lineages (e.g. Littorinimorpha) within the group. Neogastropoda has been classified into up to six superfamilies including Buccinoidea, Muricoidea, Olivoidea, Pseudolivoidea, Conoidea, and Cancellarioidea. Phylogenetic relationships among neogastropod superfamilies remain unresolved.     

Classification analysis of triplet periodicity in protein-coding regions of genes

We introduce a new concept of triplet periodicity class (TPC) and a measure of similarity between such classes. We performed classification of 472288 triplet periodicity (TP) regions found in 578868 genes from 29th release of KEGG databank. Totally 2520 classes were obtained. They contain 94% of 472288 found cases of TP. For 92% of TP regions contained in classes the same linkage of TP to open reading frame (ORF) is observed. For 8% of TP cases we revealed a shift between ORF of a gene and ORF common for majority of genes contained in a TPC. For these 8% of periodic regions the hypothetical amino acid sequences corresponding to ORF built by TPC were made. BLAST program has shown that 2679 hypothetical amino acid sequences have statistically significant similarity with proteins from UniProt databank. We suppose that 8% of TP regions contained in classes possess a mutation originating from ORF shift. Obtained TPCs can be used for identification of genes' coding regions as well as for searching for mutations arisen arising from ORF shift.     

How do insect nuclear ribosomal genes compare to protein-coding genes in phylogenetic utility and nucleotide substitution patterns?

Abstract. The expanding data set on insect molecular systematics allows examination of phylogenetic performance and molecular evolution of different types of gene. Studies combining more than one gene in the same analysis allow examination of the relative contribution and performance of each gene partition and can help inform gene choice for resolving deep and/or problematic divergences. We compared results obtained from analyses of twelve insect data sets in which authors combined one or more nuclear ribosomal genes (28S and/or 18S) with one or more protein-coding genes [elongation factor-1α (EF-1α), histone H3, carbamoylphosphate synthetase domain (CPS domain of CAD, or rudimentary), long-wavelength rhodopsin (LW opsin), glucose-6-phosphate dehydrogenase (G₆pd), phosphoenolpyruvate carboxykinase (PEPCK), arginine kinase, and white]. Data sets examined spanned eight orders of insects (Odonata, Ephemeroptera, Hemiptera, Coleoptera, Trichoptera, Lepidoptera, Diptera and Hymenoptera), providing a broad range of divergence times and taxonomic levels. We estimated the phylogenetic utility of the individual genes (using parsimony methods) and characterized the nucleotide substitution patterns (using Bayesian methods) to ask which type of data is preferable for phylogenetic analysis in insects. Nuclear ribosomal and protein coding genes differed little in our measures of phylogenetic performance and patterns of nucleotide substitution. We recommend combining nuclear ribosomal gene data with nuclear protein-coding gene data because each data set has distinct advantages. We do not recommend using mitochondrial genes for higher-level studies of insect phylogeny because reviewed studies demonstrate substitution patterns that lead to high levels of homoplasy.     

Evolutionary rates, divergence dates, and the performance of mitochondrial genes in Bayesian phylogenetic analysis

The mitochondrial genome is one of the most frequently used loci in phylogenetic and phylogeographic analyses, and it is becoming increasingly possible to sequence and analyze this genome in its entirety from diverse taxa. However, sequencing the entire genome is not always desirable or feasible. Which genes should be selected to best infer the evolutionary history of the mitochondria within a group of organisms, and what properties of a gene determine its phylogenetic performance? The current study addresses these questions in a Bayesian phylogenetic framework with reference to a phylogeny of plethodontid and related salamanders derived from 27 complete mitochondrial genomes; this topology is corroborated by nuclear DNA and morphological data. Evolutionary rates for each mitochondrial gene and divergence dates for all nodes in the plethodontid mitochondrial genome phylogeny were estimated in both Bayesian and maximum likelihood frameworks using multiple fossil calibrations, multiple data partitions, and a clock-independent approach. Bayesian analyses of individual genes were performed, and the resulting trees compared against the reference topology. Ordinal logistic regression analysis of molecular evolution rate, gene length, and the G-shape parameter a demonstrated that slower rate of evolution and longer gene length both increased the probability that a gene would perform well phylogenetically. Estimated rates of molecular evolution vary 84-fold among different mitochondrial genes and different salamander lineages, and mean rates among genes vary 15-fold. Despite having conserved amino acid sequences, cox1, cox2, cox3, and cob have the fastest mean rates of nucleotide substitution, and the greatest variation in rates, whereas rrnS and rrnL have the slowest rates. Reasons underlying this rate variation are discussed, as is the extensive rate variation in cox1 in light of its proposed role in DNA barcoding.     

Phylogenetic analysis of Myriapoda using three nuclear protein-coding genes

We assessed the ability of three nuclear protein-encoding genes-elongation factor-1alpha (EF-1alpha), RNA polymerase II (Pol II), and elongation factor-2 (EF-2)-from 59 myriapod and 12 non-myriapod species to resolve phylogenetic relationships among myriapod classes and orders. In a previous study using EF-1alpha and Pol II (2134 nt combined) from 34 myriapod taxa, Regier and Shultz recovered widely accepted classes, orders, and families but failed to resolve interclass and interordinal relationships. The result was attributed to heterogenous rates of cladogenesis (specifically, the inability of the slowly evolving sequences to capture phylogenetic signal during rapid phylogenetic diversification) but the possibility of inadequate taxon sampling or limited sequence information could not be excluded. In the present study, the myriapod taxon sample was increased by 25 taxa (73%) and sequence length per taxon was effectively doubled through addition of EF-2 (4318 nt combined). Parsimony and Bayesian analyses of the expanded data set recovered a monophyletic Myriapoda, all four myriapod classes and all multiply sampled orders, often with high node support. However, except for three diplopod clades (Colobognatha, Helminothomorpha, and a subgroup of Pentazonia), few interordinal relationships and no interclass relationships were well supported. These results are similar to those of the earlier study by Regier and Shultz, which indicates that taxon sample and sequence length alone do not readily explain the weakly supported resolution in the earlier study. We review recent paleontological evidence to further develop our proposal that heterogeneity in phylogenetic signal provided by our slowly evolving sequences is due to heterogeneity in the temporal structure of myriapod diversification.     

Phylogenetic network and physicochemical properties of nonsynonymous mutations in the protein-coding genes of human mitochondrial DNA

Theories on molecular evolution predict that phylogenetically recent nonsynonymous mutations should contain more non-neutral amino acid replacements than ancient mutations. We analyzed 840 complete coding-region human mitochondrial DNA (mtDNA) sequences for nonsynonymous mutations and evaluated the mutations in terms of the physicochemical properties of the amino acids involved. We identified 465 distinct missense and 6 nonsense mutations. 48% of the amino acid replacements changed polarity, 26% size, 8% charge, 32% aliphaticity, 13% aromaticity, and 44% hydropathy. The reduced-median networks of the amino acid changes revealed relatively few differences between the major continent-specific haplogroups, but a high variation and highly starlike phylogenies within the haplogroups. Some 56% of the mutations were private, and 25% were homoplasic. Nonconservative changes were more common than expected among the private mutations but less common among the homoplasic mutations. The asymptotic maximum of the number of nonsynonymous mutations in European mtDNA was estimated to be 1,081. The results suggested that amino acid replacements in the periphery of phylogenetic networks are more deleterious than those in the central parts, indicating that purifying selection prevents the fixation of some alleles.