Resolution of phylogenetic relationships among recently evolved species as a function of amount of DNA sequence: an empirical study based on woodpeckers (Aves: Picidae)

Synonymous substitutions in the 13 mitochondrial encoded protein genes form a large pool of characters that should approach the ideal for phylogenetic analysis of being independently and identically distributed. Pooling sequences from multiple mitochondrial protein-coding genes should result in statistically more powerful estimates of relationships among species that diverged sufficiently recently that most nucleotide substitutions are synonymous. Cytochrome oxidase I (COI) was sequenced for woodpecker species for which cytochrome b (cyt b) sequences were available. A pairing-design test based on the normal distribution indicated that cyt b evolves more rapidly than COI when all nucleotides are compared but their rates are equal for synonymous substitutions. Nearly all of the phylogenetically informative substitutions among woodpeckers are synonymous. Statistical support for relationships, as measured by bootstrap proportions, increased as the number of nucleotides increased from 1047 (cyt b) to 1512 (COI) to 2559 nucleotides (aggregate data set). Pseudo-bootstrap replicates showed the same trend and increasing the amount of sequence beyond the actual length of 2559 nucleotides to 5120 (2x) resulted in stronger bootstrap support, even though the amount of phylogenetic information was the same. However, the amount of sequence required to resolve an internode depends on the length of the internode and its depth in the phylogeny.     

我国沿海棱鳀属鱼类的物种鉴定与系统发育

利用DNA条形码技术对中国沿海分布的6种棱鳀属(Thryssa)鱼类样品进行了物种鉴定, 并每种取5尾用于探讨该属系统发育关系。结果显示: 棱鳀属鱼类的主要形态鉴别特征为上颌骨伸达位置和第一鳃耙的下鳃耙数量。在525 bp的目的片段上有175个变异位点, 其中简约信息位点172个, 单一信息位点3个, 无插入缺失现象, 转换数为182, 颠换数为57。A+T含量明显高于G+C含量, 并且表现出明显的反G偏倚。结合GenBank中相关的同源序列进行比较发现, 所有序列明显分为10个组群, 表明已提交的棱鳀属鱼类COI基因序列中仍存在一定的问题。从各组群间的遗传距离和氨基酸遗传差异水平可以看出, 10个组群应为不同的有效种, 但是否存在隐存种还有待于进一步确定。从NJ树上可以看出, 长颌棱鳀(T. setirostris)是最先分化出的物种, 保持着最原始的特征, 而中颌棱鳀(T. mystax)与黄吻棱鳀(T. vitrirostris)聚类到一起, 二者间存在共享单倍型。棱鳀属鱼类最早分化于中新世早期。在今后的研究中仍需要结合更多的分子标记对中颌棱鳀和黄吻棱鳀的分类地位作进一步的探讨。     

Effect of genetic convergence on phylogenetic inference

Phylogenetic reconstructions are a major component of many studies in evolutionary biology, but their accuracy can be reduced under certain conditions. Recent studies showed that the convergent evolution of some phenotypes resulted from recurrent amino acid substitutions in genes belonging to distant lineages. It has been suggested that these convergent substitutions could bias phylogenetic reconstruction toward grouping convergent phenotypes together, but such an effect has never been appropriately tested. We used computer simulations to determine the effect of convergent substitutions on the accuracy of phylogenetic inference. We show that, in some realistic conditions, even a relatively small proportion of convergent codons can strongly bias phylogenetic reconstruction, especially when amino acid sequences are used as characters. The strength of this bias does not depend on the reconstruction method but varies as a function of how much divergence had occurred among the lineages prior to any episodes of convergent substitutions. While the occurrence of this bias is difficult to predict, the risk of spurious groupings is strongly decreased by considering only 3rd codon positions, which are less subject to selection, as long as saturation problems are not present. Therefore, we recommend that, whenever possible, topologies obtained with amino acid sequences and 3rd codon positions be compared to identify potential phylogenetic biases and avoid evolutionarily misleading conclusions.     

A cryptic, intergeneric cytochrome c oxidase I pseudogene in tyrant flycatchers (family: Tyrannidae)

Nuclear mitochondrial pseudogenes, or "numts", are nonfunctional copies of mitochondrial genes that have been translocated to the nuclear genome. Numts have been used to study differences in mutation rates between the nuclear and mitochondrial genomes, but have also been implicated as troublesome for phylogenetic studies and DNA-based species identification (i.e., DNA barcoding). In this study, a suspected numt discovered during a study of mitochondrial cytochrome c oxidase I (COI) diversity in North American birds was targeted and sequenced from tyrant flycatchers (family: Tyrannidae). In total, the numt was found in five taxa representing two genera. Substitution rates were compared between COI and numt sequences. None of the numt sequences harboured stop codons nor frameshift mutations, but phylogenetic analysis revealed they had accumulated more amino acid substitutions than the mitochondrial COI sequences. Mitochondrial COI appeared to be preferentially amplified in most cases, but methods for numt detection are discussed for cases like this where sequences lack obvious features for identification. Because of its persistence across a broad taxonomic lineage, this numt could form a valuable model system for studying evolution in numts. The full size of the numt and its location within the nuclear genome are yet to be determined.     

The distribution and molecular diversity of the Eastern Atlantic and Mediterranean chthamalids (Crustacea, Cirripedia)

The three chthamalids Chthamalus stellatus , C. montagui and Euraphia depressa are common inhabitants of the intertidal zone in the Eastern Atlantic, Mediterranean Sea and Black Sea. In this study, we investigated the occurrence of these barnacles in a wide range of their distribution. Population divergences of these two species have been inferred using three molecular markers — internal transcribed spacer (ITS), elongation factor 1α (EF-1α) and cytochrome oxidase subunit I (COI). ITS sequences of C. stellatus were identical throughout the species range, whereas ITS sequences of C. montagui indicated that the Black Sea and Mediterranean populations are isolated from the Atlantic population. The COI and EF-1α sequences were the most variable and informative. They indicated a high genetic divergence between Atlantic, Mediterranean and Black Sea populations for C. montagui . In addition significant genetic structure was found among the populations of C. stellatus based on EF-1α but not COI. Interestingly, our molecular dating analysis correlated the pattern of diversification in C. montagui to major geological changes that occurred in the Mediterranean during the end of the Messinian and Pleiocene periods. We suggest that palaeohistory shaped the divergences between Chthamalus populations that have probably been maintained by current hydrographic conditions. Finally, COI phylogenetic analysis placed the genus Euraphia within the Chthamalus clade, suggesting the need for a taxonomic revision of Euraphia . This study represents the most detailed phylogeographical analysis of intertidal Mediterranean species to date, and shows that geological events have strongly shaped the current diversity pattern of this fauna.     

Effect of co-evolving amino acid residues on topology of phylogenetic trees

The presence in proteins of amino acid residues that change in concert during evolution is associated with keeping constant the protein spatial structure and functions. As in the case with morphological features, correlated substitutions may become the cause of homoplasies--the independent evolution of identical non-homological adaptations. Our data obtained on model phylogenetic trees and corresponding sets of sequences have shown that the presence of correlated substitutions distorts the results of phylogenetic reconstructions. A method for accounting for co-evolving amino acid residues in phylogenetic analysis is proposed. According to this method, only a single site from the group of correlated amino acid positions should remain, whereas other positions should not be used in further phylogenetic analysis. Simulations performed have shown that replacement on the average of 8% of variable positions in a pair of model sequences by coordinately evolving amino acid residues is able to change the tree topology. The removal of such amino acid residues from sequences before phylogenetic analysis restores the correct topology.     

Phylogenetic analysis of the mitochondrial cytochrome oxidase subunit I gene of five species of the Culicoides imicola species complex

The phylogenetic status of members of the Culicoides imicola Kieffer (Diptera: Ceratopogonidae) species complex of haematophagous midges is unknown, and simple means to identify the members using all life stages are unavailable. In this study, the status of three confirmed (C. imicola s.s., C. bolitinos Meiswinkel and C. loxodontis Meiswinkel) and two provisional (C. tuttifrutti Meiswinkel and C. kwagga Meiswinkel) members of the complex from South Africa was assessed using phylogenetic analysis of partial DNA and amino acid sequences of the mitochondrial cytochrome oxidase subunit I (COI) gene. The four or five individuals of each species analysed contained one or two haplotypes each. Interspecific divergence was significant and characterized by strong A <--> T transversion bias. Phylogenetic trees constructed using neighbour-joining, parsimony and maximum likelihood showed each species to be distinct. Combinations of sites for two restriction enzymes in the COI sequences were species-specific and could form the basis of a diagnostic PCR assay.     

The Phylogenetic Informativeness of Nucleotide and Amino Acid Sequences for Reconstructing the Vertebrate Tree

To aid in future efforts to accurately reconstruct the vertebrate tree, a quantitative measure of phylogenetic informativeness was applied to nucleotide and amino acid sequences for a set of 11 genes. We identified orthologues and assembled published fossil-calibrated divergence times between taxa that had been sequenced for each gene. Rates of molecular evolution for each site were estimated to characterize the molecular evolutionary pattern of genes and to calculate the phylogenetic informativeness. The fast-evolving gene albumin yielded the highest informativeness over the period from 60 million years ago to 500 million years ago. In contrast, calmodulin yielded the lowest informativeness, presumably because functional constraint minimized substitutions in the amino acid sequence. The gene c-myc showed an intermediate level of informativeness. The nucleotide sequence of cytochrome b showed extremely high utility for recent epochs, but low utility for times before 100 million years ago. We ranked nine other genes for their utility during the epochs of the divergence of the muroid rodents, early placental mammals, early vertebrates, and early metazoa, yielding results consistent with, but more precise than, previous studies. Interestingly, DNA sequence always exceeded amino acid sequence in informativeness over all time scales, yet support values were at best moderately higher. For epochs not subject to strong phylogenetic conflict due to convergence, we advocate gleaning the additional power of the threefold increase in number of characters that is present for DNA sequences over resorting to the less noisy but less informative amino acid sequences.     

Genomic structure of the luciferase gene and phylogenetic analysis in the Hotaria-group fireflies

The luminescent fireflies have species specific flash patterns, being recognized as sexual communication. The luciferase gene is the sole enzyme responsible for bioluminescence. We describe here the complete nucleotide sequence and the exon-intron structure of the luciferase gene of the Hotaria-group fireflies, H. unmunsana, H. papariensis and H. tsushimana. The luciferase gene of the Hotaria-group firefly including the known H. parvula spans 1950 bp and consisted of six introns and seven exons coding for 548 amino acid residues, suggesting highly conserved structure among the Hotaria-group fireflies. Although only one luciferase gene was cloned from H. papariensis, each of the two sequences of the gene was found in H. unmunsana (U1 and Uc) and H. tsushimana (T1 and T2). The amino acid sequence divergence among H. unmunsana, H. papariensis, and H. tsushimana only ranged from zero to three amino acid residues, but H. parvula differed by 10-11 amino acid residues from the other Hotaria-group fireflies, suggesting a divergent relationship of this species. Phylogenetic analysis using the deduced amino acid sequences of the luciferase gene resulted in a monophyletic group in the Hotaria excluding H. parvula, suggesting a close relationship among H. unmunsana, H. papariensis and H. tsushimana. Additionally, we also analyzed the mitochondrial cytochrome oxidase I (COI) gene of the Hotaria-group fireflies. The deduced amino acid sequence of the COI gene of H. unmunsana was identical to that of H. papariensis and H. tsushimana, but different by three positions from H. parvula. In terms of nucleotide sequences of the COI gene, intraspecific sequence divergence was sometimes larger than interspecies level, and phylogenetic analysis placed the three species into monophyletic groups unresolved among them, but excluded H. parvula. In conclusion, our results suggest that H. unmunsana, H. papariensis and H. tsushimana are very closely related or might be an identical species, at least based on the luciferase and COI genes.     

Evolution of the primate beta-globin gene region. High rate of variation in CpG dinucleotides and in short repeated sequences between man and chimpanzee

A 5500 base-pair fragment including the beta-globin gene downstream from codon 122 and about 4000 base-pairs of its 5' flanking sequence was cloned from chimpanzee DNA and thoroughly sequenced before being compared with the corresponding human sequence: 88 point differences (83 substitutions and 5 deletions or insertions of 1 base-pair) were detected as well as seven more important deletion/insertion events. These changes occur preferentially in two kinds of structure. First, 40% of the CpG dinucleotides present in either human or chimpanzee sequences are affected by nucleotide variations. This corresponds to a divergence level considerably higher than that expected. Second, most short repeated sequences found in the 5' extragenic sequence are involved in mutational events (amplification or contraction of the number of basic motifs as well as point substitutions or deletions/insertions of 1 base-pair). Considering the very low level of nucleotide sequence divergence between these two closely related species, our data provide direct evidence for CpG and tandem array instability.     

c-mos variation in songbirds: molecular evolution, phylogenetic implications, and comparisons with mitochondrial differentiation

Nucleotide sequences from the c-mos proto-oncogene have previously been used to reconstruct the phylogenetic relationships between distantly related vertebrate taxa. To explore c-mos variation at shallower levels of avian divergence, we compared c-mos sequences from representative passerine taxa that span a range of evolutionary differentiation, from basal passerine lineages to closely allied genera. Phylogenetic reconstructions based on these c-mos sequences recovered topologies congruent with previous DNA-DNA hybridization-based reconstructions, with many nodes receiving high support, as indicated by bootstrap and reliability values. One exception was the relationship of Acanthisitta to the remaining passerines, where the c-mos-based searches indicated a three-way polytomy involving the Acanthisitta lineage and the suboscine and oscine passerine clades. We also compared levels of c-mos and mitochondrial differentiation across eight oscine passerine taxa and found that c-mos nucleotide substitutions accumulate at a rate similar to that of transversion substitutions in mitochondrial protein-coding genes. These comparisons suggest that nuclear-encoded loci such as c-mos provide a temporal window of phylogenetic resolution that overlaps the temporal range where mitochondrial protein-coding sequences have their greatest utility and that c-mos substitutions and mtDNA transversions can serve as complementary, informative, and independent phylogenetic markers for the study of avian relationships.     

线蛱蝶亚科蝶类部分类群线粒体COI基因的分子系统发生分析

以线粒体细胞色素氧化酶I（COI）基因作分子标记,对线蛱蝶亚科蝴蝶进行序列测定．序列分析的结果表明．经比对和处理后的序列总长度是645bp,其中有199个变异位点,147个简约信息位点;所编码的氨基酸序列中有18个变异位点,7个信息位点．A＋T平均含量为69．6％,G＋C平均含量为30．4％,碱基组成出现AT偏斜．以蛱蝶亚科及秀蛱蝶亚科物种为外类群,用NJ、MP及贝叶斯法重建了该亚科的系统发生树,探讨了它们主要类群间的系统发生关系．分子系统树显示,线蛱蝶亚科由以下3大支系：环蛱蝶族＋翠蛱蝶族、线蛱蝶族、丽蛱蝶族构成;其中,环蛱蝶族为单系群（NJ树也支持线蛱蝶族的单系性）;翠蛱蝶族与环蛱蝶族亲缘关系较近：丽蛱蝶族可能是该亚科较早分化出的一支．     

Species status and phylogeography of two closely related Coptolabrus species (Coleoptera: Carabidae) in South Korea inferred from mitochondrial and nuclear gene sequences

We investigated the species status and intraspecific phylogeography in South Korea of two ground beetle species, Coptolabrus jankowskii and Coptolabrus smaragdinus (Coleoptera: Carabidae), using statistical parsimony networks and nested clade analyses based on sequences from the mitochondrial cytochrome oxidase subunit I (COI) and nuclear phosphoenolpyruvate carboxykinase (PepCK) and wingless (Wg) genes. Although traditional parsimony tree construction generally failed to resolve interspecific relationships and construct biologically meaningful genealogies, analysis using statistical parsimony networks yielded statistically significant inter- and intraspecific genealogical structures. We found that although these two species represent a notable case of trans-species polymorphisms in both mitochondrial and nuclear gene sequences, their status as separate species was evidenced by the nonrandom association between species and nested clades at various nesting levels. The exceptional occurrence of shared identical or very similar COI sequences was considered to be the result of introgressive hybridization. In addition, range expansion and fragmentation events across the Korean Peninsula and adjacent islands were inferred from nested clade phylogeographical analyses. The COI gene revealed the geographical divergence of major eastern and western clades and historical biogeographical events within each major clade, whereas the nuclear PepCK gene, which did not reveal corresponding east-west clades, indicated past fragmentation and range expansion across wide areas that may have been the result of older biogeographical events. Thus, phylogeographical inferences drawn from analyses of mitochondrial and nuclear genes can reveal different and potentially complementary information about phylogeographical processes.     

青海湖地区狗獾分类地位和狗獾属进化历史探讨

本研究选择线粒体细胞色素b (cytochrome b, Cyt b)和控制区(control region, CR)片段作为分子标记, 探讨了青海湖地区狗獾(Meles sp.)的系统发育地位和狗獾属分歧时间。研究结果支持目前将狗獾属分为4个种的结论。Cyt b和CR片段序列拼接后总长1,652 bp, 23条序列共定义了21个单倍型。研究结果表明欧亚大陆狗獾分为东西两个支系, 每个支系进一步分为两个种: 东部支系包括亚洲狗獾(M. leucurus)和日本狗獾(M. anakuma); 西部支系包括欧洲狗獾(M. meles)和西南亚狗獾(M. canescens)。贝叶斯树和单倍型网络关系图都支持青海湖地区狗獾属于亚洲狗獾。分歧时间的估算结果与古生物学证据相符, 东部支系和西部支系在2.24 Ma左右产生分歧, 西南亚狗獾在1.27 Ma左右从欧洲狗獾分出, 而日本狗獾和亚洲狗獾的分化时间为0.99 Ma左右。     

Essential factors determining codon usage in ubiquitin genes

Summary Ubiquitin is ubiquitous in all eukaryotes and its amino acid sequence shows extreme conservation. Ubiquitin genes comprise direct repeats of the ubiquitin coding unit with no spacers. The nucleotide sequences coding for 13 ubiquitin genes from 11 species reported so far have been compiled and analyzed. The G+C content of codon third base reveals a positive linear correlation with the genome G+C content of the corresponding species. The slope strongly suggests that the overall G+C content of codons of polyubiquitin genes clearly reflects the genome G+C content by AT/GC substitutions at the codon third position. The G+C content of ubiquitin codon third base also shows a positive linear correlation with the overall G+C content of coding regions of compiled genes, indicating the codon choices among synonymous codons reflect the average codon usage pattern of corresponding species. On the other hand, the monoubiquitin gene, which is different from the polyubiquitin gene in gene organization, gene expression, and function of the encoding protein, shows a different codon usage pattern compared with that of the polyubiquitin gene. From comparisons of the levels of synonymous substitutions among ubiquitin repeats and the homology of the amino acid sequence of the tail of monomeric ubiquitin genes, we propose that the molecular evolution of ubiquitin genes occurred as follows: Plural primitive ubiquitin sequences were dispersed on genome in ancestral eukaryotes. Some of them situated in a particular environment fused with the tail sequence to produce monomeric ubiquitin genes that were maintained across species. After divergence of species, polyubiquitin genes were formed by duplication of the other primitive ubiquitin sequences on different chromosomes. Differences in the environments in which ubiquitin genes are embedded reflect the differences in codon choice and in gene expression pattern between poly- and monomeric ubiquitin genes.     

New molecular data for tardigrade phylogeny, with the erection of Paramacrobiotus gen. nov.

Up to few years ago, the phylogenies of tardigrade taxa have been investigated using morphological data, but relationships within and between many taxa are still unresolved. Our aim has been to verify those relationships adding molecular analysis to morphological analysis, using nearly complete 18S ribosomal DNA gene sequences (five new) of 19 species, as well as cytochrome oxidase subunit 1 (COI) mitochondrial DNA gene sequences (15 new) from 20 species, from a total of seven families. The 18S rDNA tree was calculated by minimum evolution, maximum parsimony (MP) and maximum likelihood (ML) analyses. DNA sequences coding for COI were translated to amino acid sequences and a tree was also calculated by neighbour-joining, MP and ML analyses. For both trees (18S rDNA and COI) posterior probabilities were calculated by MrBayes. Prominent findings are as follows: the molecular data on Echiniscidae (Heterotardigrada) are in line with the phylogenetic relationships identifiable by morphological analysis. Among Eutardigrada, orders Apochela and Parachela are confirmed as sister groups. Ramazzottius (Hypsibiidae) results more related to Macrobiotidae than to the genera here considered of Hypsibiidae. Macrobiotidae and Macrobiotus result not monophyletic and confirm morphological data on the presence of at least two large groups within Macrobiotus. Using 18S rDNA and COI mtDNA genes, a new phylogenetic line has been identified within Macrobiotus , corresponding to the ' richtersi-areolatus group'. Moreover, cryptic species have been identified within the Macrobiotus ' richtersi group' and within Richtersius . Some evolutionary lines of tardigrades are confirmed, but others suggest taxonomic revision. In particular, the new genus Paramacrobiotus gen. n. has been identified, corresponding to the phylogenetic line represented by the ' richtersi-areolatus group'.     

Adaptive intraspecific divergence: an example using the animal cytochrome b gene

The topologies of phylogenetic trees characterized by a high level of intraspecific divergence between the phylogenetic DNA groups (clades) are often explained in terms of the theory of Pleistocene refugia. To elucidate the issue of the adaptive role of intraspecific divergence, the changes in the physicochemical properties of amino acids in the course ofcladogenesis (MM01 model of the TreeSAAP 3.2 package) were analyzed in this work using as an example the nucleotide sequences of the cytochrome b gene in some species of northern animals (lemmings, redbacked voles, chipmunks, flying squirrels, ermines). It was shown that the process of intraspecific divergence was rarely accompanied by radical amino acid substitutions in cytochrome b caused by adaptation (directional selection). In connection with this, the hypothesis is discussed according to which the adaptive variants formed in the species at the peak of cold were lost with climatic warming due to the drift or selection against individuals adapted to cold.     

COI (cytochrome oxidase-I) sequence based studies of Carangid fishes from Kakinada coast,India

Mitochondrial DNA, cytochrome oxidase-1 gene sequences were analyzed for species identification and phylogenetic relationship among the very high food value and commercially important Indian carangid fish species. Sequence analysis of COI gene very clearly indicated that all the 28 fish species fell into five distinct groups, which are genetically distant from each other and exhibited identical phylogenetic reservation. All the COI gene sequences from 28 fishes provide sufficient phylogenetic information and evolutionary relationship to distinguish the carangid species unambiguously. This study proves the utility of mtDNA COI gene sequence based approach in identifying fish species at a faster pace.