Multiple origins and rapid evolution of duplicated mitochondrial genes in parthenogenetic geckos (Heteronotia binoei; Squamata, Gekkonidae)

Accumulating evidence for alternative gene orders demonstrates that vertebrate mitochondrial genomes are more evolutionarily dynamic than previously thought. Several lineages of parthenogenetic lizards contain large, tandem duplications that include rRNA, tRNA, and protein-coding genes, as well as the control region. Such duplications are hypothesized as intermediate stages in gene rearrangement, but the early stages of their evolution have not been previously studied. To better understand the evolutionary dynamics of duplicated segments of mitochondrial DNA, we sequenced 10 mitochondrial genomes from recently formed ( approximately 300,000 years ago) hybrid parthenogenetic geckos of the Heteronotia binoei complex and 1 from a sexual form. These genomes included some with an arrangement typical of vertebrates and others with tandem duplications varying in size from 5.7 to 9.4 kb, each with different gene contents and duplication endpoints. These results, together with phylogenetic analyses, indicate independent and frequent origins of the duplications. Small, direct repeats at the duplication endpoints imply slipped-strand error as a mechanism generating the duplications as opposed to a false initiation/termination of DNA replication mechanism that has been invoked to explain duplications in other lizard mitochondrial systems. Despite their recent origin, there is evidence for nonfunctionalization of genes due primarily to deletions, and the observed pattern of gene disruption supports the duplication-deletion model for rearrangement of mtDNA gene order. Conversely, the accumulation of mutations between these recent duplicates provides no evidence for gene conversion, as has been reported in some other systems. These results demonstrate that, despite their long-term stasis in gene content and arrangement in some lineages, vertebrate mitochondrial genomes can be evolutionary dynamic even at short timescales.     

Evolution of craniofacial novelty in parrots through developmental modularity and heterochrony

Parrots (order Psittaciformes) have developed novel cranial morphology. At the same time, they show considerable morphological diversity in the cranial musculoskeletal system, which includes two novel structures: the suborbital arch and the musculus (M.) pseudomasseter. To understand comprehensively the evolutionary pattern and process of novel cranial morphology in parrots, phylogenetic and developmental studies were conducted. Firstly, we undertook phylogenetic analyses based on mitochondrial ribosomal RNA gene sequences to obtain a robust phylogeny among parrots, and secondly we surveyed the cranial morphology of parrots extensively to add new information on the character states. Character mapping onto molecular phylogenies indicated strongly the repeated evolution of both the suborbital arch and the well-developed M. pseudomasseter within parrots. These results also suggested that the direction of evolutionary change is not always identical in the two characters, implying that these characters are relatively independent or decoupled structures behaving as separate modules. Finally, we compared the developmental pattern of jaw muscles among bird species and found a difference in the timing of M. pseudomasseter differentiation between the cockatiel Nymphicus hollandicus (representative of a well-developed condition) and the peach-faced lovebird Agapornis roseicollis (representative of an underdeveloped condition). On the basis of this study, we suggest that in the development of novel traits, modularity and heterochrony facilitate the diversification of parrot cranial morphology.     

Molecular mechanisms of extensive mitochondrial gene rearrangement in plethodontid salamanders

Extensive gene rearrangement is reported in the mitochondrial genomes of lungless salamanders (Plethodontidae). In each genome with a novel gene order, there is evidence that the rearrangement was mediated by duplication of part of the mitochondrial genome, including the presence of both pseudogenes and additional, presumably functional, copies of duplicated genes. All rearrangement-mediating duplications include either the origin of light-strand replication and the nearby tRNA genes or the regions flanking the origin of heavy-strand replication. The latter regions comprise nad6, trnE, cob, trnT, an intergenic spacer between trnT and trnP and, in some genomes, trnP, the control region, trnF, rrnS, trnV, rrnL, trnL1, and nad1. In some cases, two copies of duplicated genes, presumptive regulatory regions, and/or sequences with no assignable function have been retained in the genome following the initial duplication; in other genomes, only one of the duplicated copies has been retained. Both tandem and nontandem duplications are present in these genomes, suggesting different duplication mechanisms. In some of these mitochondrial DNAs, up to 25% of the total length is composed of tandem duplications of noncoding sequence that includes putative regulatory regions and/or pseudogenes of tRNAs and protein-coding genes along with the otherwise unassignable sequences. These data indicate that imprecise initiation and termination of replication, slipped-strand mispairing, and intramolecular recombination may all have played a role in generating repeats during the evolutionary history of plethodontid mitochondrial genomes.     

Complex phylogeographic history of central African forest elephants and its implications for taxonomy

Background

The mitochondrial genomes of snakes are characterized by an overall evolutionary rate that appears to be one of the most accelerated among vertebrates. They also possess other unusual features, including short tRNAs and other genes, and a duplicated control region that has been stably maintained since it originated more than 70 million years ago. Here, we provide a detailed analysis of evolutionary dynamics in snake mitochondrial genomes to better understand the basis of these extreme characteristics, and to explore the relationship between mitochondrial genome molecular evolution, genome architecture, and molecular function. We sequenced complete mitochondrial genomes from Slowinski's corn snake (Pantherophis slowinskii) and two cottonmouths (Agkistrodon piscivorus) to complement previously existing mitochondrial genomes, and to provide an improved comparative view of how genome architecture affects molecular evolution at contrasting levels of divergence.

Results

We present a Bayesian genetic approach that suggests that the duplicated control region can function as an additional origin of heavy strand replication. The two control regions also appear to have different intra-specific versus inter-specific evolutionary dynamics that may be associated with complex modes of concerted evolution. We find that different genomic regions have experienced substantial accelerated evolution along early branches in snakes, with different genes having experienced dramatic accelerations along specific branches. Some of these accelerations appear to coincide with, or subsequent to, the shortening of various mitochondrial genes and the duplication of the control region and flanking tRNAs.

Conclusion

Fluctuations in the strength and pattern of selection during snake evolution have had widely varying gene-specific effects on substitution rates, and these rate accelerations may have been functionally related to unusual changes in genomic architecture. The among-lineage and among-gene variation in rate dynamics observed in snakes is the most extreme thus far observed in animal genomes, and provides an important study system for further evaluating the biochemical and physiological basis of evolutionary pressures in vertebrate mitochondria.     

The evolutionary history of cockatoos (Aves: Psittaciformes: Cacatuidae)

Cockatoos are the distinctive family Cacatuidae, a major lineage of the order of parrots (Psittaciformes) and distributed throughout the Australasian region of the world. However, the evolutionary history of cockatoos is not well understood. We investigated the phylogeny of cockatoos based on three mitochondrial and three nuclear DNA genes obtained from 16 of 21 species of Cacatuidae. In addition, five novel mitochondrial genomes were used to estimate time of divergence and our estimates indicate Cacatuidae diverged from Psittacidae approximately 40.7 million years ago (95% CI 51.6-30.3 Ma) during the Eocene. Our data shows Cacatuidae began to diversify approximately 27.9 Ma (95% CI 38.1-18.3 Ma) during the Oligocene. The early to middle Miocene (20-10 Ma) was a significant period in the evolution of modern Australian environments and vegetation, in which a transformation from mainly mesic to xeric habitats (e.g., fire-adapted sclerophyll vegetation and grasslands) occurred. We hypothesize that this environmental transformation was a driving force behind the diversification of cockatoos. A detailed multi-locus molecular phylogeny enabled us to resolve the phylogenetic placements of the Palm Cockatoo (Probosciger aterrimus), Galah (Eolophus roseicapillus), Gang-gang Cockatoo (Callocephalon fimbriatum) and Cockatiel (Nymphicus hollandicus), which have historically been difficult to place within Cacatuidae. When the molecular evidence is analysed in concert with morphology, it is clear that many of the cockatoo species' diagnostic phenotypic traits such as plumage colour, body size, wing shape and bill morphology have evolved in parallel or convergently across lineages.     

Mitochondrial genomes and avian phylogeny: complex characters and resolvability without explosive radiations

We improve the taxon sampling for avian phylogeny by analyzing 7 new mitochondrial genomes (a toucan, woodpecker, osprey, forest falcon, American kestrel, heron, and a pelican). This improves inference of the avian tree, and it supports 3 major conclusions. The first is that some birds (including a parrot, a toucan, and an osprey) exhibit a complete duplication of the control region (CR) meaning that there are at least 4 distinct gene orders within birds. However, it appears that there are regions of continued gene conversion between the duplicate CRs, resulting in duplications that can be stable for long evolutionary periods. Because of this stable duplicated state, gene order can eventually either revert to the original order or change to the new gene order. The existence of this stable duplicate state explains how an apparently unlikely event (finding the same novel gene order) can arise multiple times. Although rare genomic changes have theoretical advantages for tree reconstruction, they can be compromised if these apparently rare events have a stable intermediate state. Secondly, the toucan and woodpecker improve the resolution of the 6-way split within Neoaves that has been called an "explosive radiation." An explosive radiation implies that normal microevolutionary events are insufficient to explain the observed macroevolution. By showing the avian tree is, in principle, resolvable, we demonstrate that the radiation of birds is amenable to standard evolutionary analysis. Thirdly, and as expected from theory, additional taxa breaking up long branches stabilize the position of some problematic taxa (like the falcon). In addition, we report that within the birds of prey and allies, we did not find evidence pairing New World vultures with storks or accipitrids (hawks, eagles, and osprey) with Falconids.     

Birth and death of duplicated genes in completely sequenced eukaryotes

Gene and genome duplications are commonly regarded as being of major evolutionary significance. But how often does gene duplication occur? And, once duplicated, what are the fates of duplicated genes? How do they contribute to evolution? In a recent article, Lynch and Conery analyze divergence between duplicate genes from six eukaryotic genomes. They estimate the rate of gene duplication, the rate of gene loss after duplication and the strength of selection experienced by duplicate genes. They conclude that although the rate of gene duplications is high, so is the rate of gene loss, and they argue that gene duplications could be a major factor in speciation.     

Tandem direct duplications of mitochondrial DNA in mitochondrial myopathy: analysis of nucleotide sequence and tissue distribution.

Two patients with direct tandem duplications of mitochondrial DNA (mtDNA) and mitochondrial myopathy are described. The breakpoint regions between duplicated segments were amplified using the polymerase chain reaction (PCR), cloned and sequenced. The distribution of normal and abnormal genomes in different tissues was investigated using Southern hybridisation, and in different cells within the same tissue using PCR. In each case the gene for cytochrome oxidase subunit I (MTCOX1) was interrupted, creating reading frames which if transcribed and translated would result in truncated versions of this peptide. Heteroplasmy and mosaicism for the abnormal mtDNA population was apparent.     

Evolutionary Dynamics of Recently Duplicated Genes: Selective Constraints on Diverging Paralogs in the <Emphasis Type="Italic">Drosophila pseudoobscura</Emphasis> Genome

Duplicated genes produce genetic variation that can influence the evolution of genomes and phenotypes. In most cases, for a duplicated gene to contribute to evolutionary novelty it must survive the early stages of divergence from its paralog without becoming a pseudogene. I examined the evolutionary dynamics of recently duplicated genes in the Drosophila pseudoobscura genome to understand the factors affecting these early stages of evolution. Paralogs located in closer proximity have higher sequence identity. This suggests that gene conversion occurs more often between duplications in close proximity or that there is more genetic independence between distant paralogs. Partially duplicated genes have a higher likelihood of pseudogenization than completely duplicated genes, but no single factor significantly contributes to the selective constraints on a completely duplicated gene. However, DNA-based duplications and duplications within chromosome arms tend to produce longer duplication tracts than retroposed and inter-arm duplications, and longer duplication tracts are more likely to contain a completely duplicated gene. Therefore, the relative position of paralogs and the mechanism of duplication indirectly affect whether a duplicated gene is retained or pseudogenized. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The largest unassigned regions of the male- and female-transmitted mitochondrial DNAs in Musculista senhousia (Bivalvia Mytilidae)

Musculista senhousia is a marine mussel with doubly uniparental inheritance (DUI) of mitochondria. In this study we analyzed the largest unassigned region (LUR) of its female- and male-transmitted mitochondrial genomes, described their fine characteristics and searched for shared features. Our results suggest that both LURs contain the control region of their respective mitochondrial genomes. The female-transmitted control region is duplicated in tandem, with the two copies evolving in concert. This makes the F-mtDNA of M. senhousia the first Bivalve mitochondrial genome with this feature. We also compared M. senhousia control regions to that of other Mytilidae, and demonstrated that signals for basic mtDNA functions are retained over evolutionary times even among the fast-evolving mitochondrial genomes of DUI species. Finally, we discussed how similarities between female and male LURs may be explained in the context of DUI evolution and if the duplicated female control region might have influenced the DUI system in this species.     

Automatic genome-wide reconstruction of phylogenetic gene trees

Gene duplication and divergence is a major evolutionary force. Despite the growing number of fully sequenced genomes, methods for investigating these events on a genome-wide scale are still in their infancy. Here, we present SYNERGY, a novel and scalable algorithm that uses sequence similarity and a given species phylogeny to reconstruct the underlying evolutionary history of all genes in a large group of species. In doing so, SYNERGY resolves homology relations and accurately distinguishes orthologs from paralogs. We applied our approach to a set of nine fully sequenced fungal genomes spanning 150 million years, generating a genome-wide catalog of orthologous groups and corresponding gene trees. Our results are highly accurate when compared to a manually curated gold standard, and are robust to the quality of input according to a novel jackknife confidence scoring. The reconstructed gene trees provide a comprehensive view of gene evolution on a genomic scale. Our approach can be applied to any set of sequenced eukaryotic species with a known phylogeny, and opens the way to systematic studies of the evolution of individual genes, molecular systems and whole genomes. Supplementary information: Supplementary data are available at Bioinformatics online.     

Sequence repeats in a polyoma virus DNA region important for gene expression.

The sequenced prototype strains (A2 and A3) of polyoma virus lack sequence duplications characteristic of other papovaviruses. However, we found that five polyoma virus strains (P16, Toronto large plaque, MV, Ts 48, and NG59R) contain tandemly duplicated sequences in a region near the late RNA leader. Although the duplications vary in size (31 to 84 base pairs) and location (between nucleotide [nt] 5068 and nt 5185), the sequence between nt 5114 and nt 5137 is contained within all five duplicated segments. This region is known to be important in polyoma virus early gene expression, and it contains sequences capable of enhancing the expression of nonviral genes. Inspection of the sequences at and around the ends of the repeats indicated that the duplications do not arise by homologous recombination, and there was no indication that a sequence-specific mechanism results in their formation. However, the variation in the structure of the repeats among different polyoma virus strains suggests that these sequence duplications are a recent evolutionary occurrence. The potential biological significance of this variation is discussed.     

Divergent origins and concerted expansion of two segmental duplications on chromosome 16

An unexpected finding of the human genome was the large fraction of the genome organized as blocks of interspersed duplicated sequence. We provide a comparative and phylogenetic analysis of a highly duplicated region of 16p12.2, which is composed of at least four different segmental duplications spanning in excess of 160 kb. We contrast the dispersal of two different segmental duplications (LCR16a and LCR16u). LCR16a, a 20 kb low-copy repeat sequence A from chromosome 16, was shown previously to contain a rapidly evolving novel hominoid gene family (morpheus) that had expanded within the last 10 million years of great ape/human evolution. We compare the dispersal of this genomic segment with a second adjacent duplication called LCR16u. The duplication contains a second putative gene family (KIAA0220/SMG1) that is represented approximately eight times within the human genome. A high degree of sequence identity (approximately 98%) was observed among the various copies of LCR16u. Comparative analyses with Old World monkey species show that LCR16a and LCR16u originated from two distinct ancestral loci. Within the human genome, at least 70% of the LCR16u copies were duplicated in concert with the LCR16a duplication. In contrast, only 30% of the chimpanzee loci show an association between LCR16a and LCR16u duplications. The data suggest that the two copies of genomic sequence were brought together during the chimpanzee/human divergence and were subsequently duplicated as a larger cassette specifically within the human lineage. The evolutionary history of these two chromosome-specific duplications supports a model of rapid expansion and evolutionary turnover among the genomes of man and the great apes.     

Hemiascomycetous yeasts at the forefront of comparative genomics

With more than a dozen species fully sequenced, as many as this partially sequenced, and more in development, yeasts are now used to explore the frontlines of comparative genomics of eukaryotes. Innovative procedures have been developed to compare and annotate genomes at various evolutionary distances, to identify short cis-acting regulatory elements, to map duplications, or to align syntenic blocks. Human and plant pathogens, in addition to yeasts that show a variety of interesting physiological properties, are included in this multidimensional comparative survey, which encompasses a very broad evolutionary range. As major steps of the evolutionary history of hemiascomycetous genomes emerge, precise questions on the general mechanisms of their evolution can be addressed, using both experimental and in silico methods.     

Probing recalcitrant problems in polyclad evolution and systematics with novel mitochondrial genome resources

For their apparent morphological simplicity, the Platyhelminthes or “flatworms” are a diverse clade found in a broad range of habitats. Their body plans have however made them difficult to robustly classify. Molecular evidence is only beginning to uncover the true evolutionary history of this clade. Here we present nine novel mitochondrial genomes from the still undersampled orders Polycladida and Rhabdocoela, assembled from short Illumina reads. In particular we present for the first time in the literature the mitochondrial sequence of a Rhabdocoel, Bothromesostoma personatum (Typhloplanidae, Mesostominae). The novel mitochondrial genomes examined generally contained the 36 genes expected in the Platyhelminthes, with all possessing 12 of the 13 protein-coding genes normally found in metazoan mitochondrial genomes (ATP8 being absent from all Platyhelminth mtDNA sequenced to date), along with two ribosomal RNA genes. The majority presented possess 22 transfer RNA genes, and a single tRNA gene was absent from two of the nine assembled genomes. By comparison of mitochondrial gene order and phylogenetic analysis of the protein coding and ribosomal RNA genes contained within these sequences with those of previously sequenced species we are able to gain a firm molecular phylogeny for the inter-relationships within this clade.Our phylogenetic reconstructions, using both nucleotide and amino acid sequences under several models and both Bayesian and Maximum Likelihood methods, strongly support the monophyly of Polycladida, and the monophyly of Acotylea and Cotylea within that clade. They also allow us to speculate on the early emergence of Macrostomida, the monophyly of a “Turbellarian-like” clade, the placement of Rhabditophora, and that of Platyhelminthes relative to the Lophotrochozoa (=Spiralia). The data presented here therefore represent a significant advance in our understanding of platyhelminth phylogeny, and will form the basis of a range of future research in the still-disputed classifications within this taxon.     

Genetic diversity,molecular phylogeny and selection evidence of the silkworm mitochondria implicated by complete resequencing of 41 genomes

Background  

Mitochondria are a valuable resource for studying the evolutionary process and deducing phylogeny. A few mitochondria genomes have been sequenced, but a comprehensive picture of the domestication event for silkworm mitochondria remains to be established. In this study, we integrate the extant data, and perform a whole genome resequencing of Japanese wild silkworm to obtain breakthrough results in silkworm mitochondrial (mt) population, and finally use these to deduce a more comprehensive phylogeny of the Bombycidae.     

The complete mitochondrial genomes of two common shrimps (Litopenaeus vannamei and Fenneropenaeus chinensis) and their phylogenomic considerations

Complete mitochondrial genomes have proven extremely valuable in helping to understand the evolutionary relationships among metazoans. However, uneven taxon sampling may lead to unclear or even erroneous phylogenetic topologies. The decapod crustaceans are relatively well-sampled, but sampling is still uneven within this group. We have sequenced the mitochondrial genomes of two shrimps Litopenaeus vannamei and Fenneropenaeus chinensis. As seen in other metazoans, the genomes contain a standard set of 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and an AT-rich non-coding region. The gene arrangements are consistent with the pancrustacean ground pattern. Both the pattern of gene rearrangements and phylogenomic analyses using concatenated nucleic acid and amino acid sequences of the 13 mitochondrial protein-coding genes strengthened the support that Caridea and Palinura are primitive members of Pleocyemata. These sequences, in combination with two previously published penaeid mitochondrial genomes, suggest that genera within the family Penaeidae have the following relationship: (((Penaeus+Fenneropenaeus)+Litopenaeus)+Marsupenaeus). The analyses of nucleic acid and amino acid sequences of the mitochondrial genomes also strongly support the monophyly of Penaeidae, Brachyura and Pleocyemata. In addition, the analyses of the average Ka/Ks in the 13 mitochondrial protein-coding genes of penaeid shrimps indicated a strong purifying selection within this group.     

Genomic architectural variation of plant mitochondria—A review of multichromosomal structuring

Since the endosymbiont origin from α-Proteobacteria, mitochondrial genomes have undergone extremely divergent evolutionary trajectories among eukaryotic lineages. Compared with the relatively compact and conserved animal mitochondrial genomes, plant mitochondrial genomes have many unique features, especially their large and complex genomic arrangements. The sizes of fully sequenced plant mitochondrial genomes span over a 100-fold range from 66 kb in Viscum scurruloideum to 11 000 kb in Silene conica. In addition to the typical circular structure, some species of plants also possess linear, and even multichromosomal, architectures. In contrast with the thousands of fully sequenced animal mitochondrial genomes and plant plastid genomes, only around 200 fully sequenced land plant mitochondrial genomes have been published, with many being only draft assemblies. In this review, we summarize some of the known novel characteristics found in plant mitochondrial genomes, with special emphasis on multichromosomal structures described in recent publications. Finally, we discuss the future prospects for studying the inheritance patterns of multichromosomal plant mitochondria and examining architectural variation at different levels of taxonomic organization—including at the population level.     

Long-term conservation of six duplicated structural genes in cephalopod mitochondrial genomes

The complete nucleotide sequences of the mitochondrial (mt) genomes of three cephalopods, Octopus vulgaris (Octopodiformes, Octopoda, Incirrata), Todarodes pacificus (Decapodiformes, Oegopsida, Ommastrephidae), and Watasenia scintillans (Decapodiformes, Oegopsida, Enoploteuthidae), were determined. These three mt genomes encode the standard set of metazoan mt genes. However, W. scintillans and T. pacificus mt genomes share duplications of the longest noncoding region, three cytochrome oxidase subunit genes and two ATP synthase subunit genes, and the tRNA(Asp) gene. Southern hybridization analysis of the W. scintillans mt genome shows that this single genome carries both duplicated regions. The near-identical sequence of the duplicates suggests that there are certain concerted evolutionary mechanisms, at least in cephalopod mitochondria. Molecular phylogenetic analyses of mt protein genes are suggestive, although not statistically significantly so, of a monophyletic relationship between W. scintillans and T. pacificus.     

两种毛蚜线粒体基因组比较分析

【目的】线粒体基因组分析已被应用于昆虫系统发育研究。本研究以蚜科Aphididae重要类群毛蚜亚科物种为代表,测定并比较分析了该类蚜虫的线粒体基因组特征,探讨了基于线粒体基因组信息的蚜虫系统发育关系重建。【方法】以毛蚜亚科三角枫多态毛蚜Periphyllus acerihabitans Zhang和针茅小毛蚜Chaetosiphella stipae Hille Ris Lambers,1947为研究对象,利用长短PCR相结合的方法测定线粒体基因组的序列,分析了基因组的基本特征;基于在线t RNAscan-SE Search Server搜索方法预测了t RNA的二级结构;基于12个物种(本研究获得的2个物种和10个Gen Bank上下载的物种数据)的蛋白编码基因(PCGs)序列,利用最大似然法和贝叶斯法重建了蚜科的系统发育关系。【结果】两种毛蚜均获得了约94%的线粒体基因组数据,P.acerihabitans获得了14 908 bp,控制区为1 205 bp;C.stipae获得了13 893 bp,控制区为609 bp。两种毛蚜同时获得33个基因,包含接近完整的13个蛋白编码基因(PCGs)(nad5不完整),18个tRNA,2个rRNA基因;ka/ks值表明,C.stipae的进化速率更快。从基因组组成、基因排列顺序、核苷酸组成分析、密码子使用情况、t RNA二级结构等特征来分析,两种蚜虫线粒体基因组基本特征相似。系统发育重建结果表明毛蚜亚科、蚜亚科的单系性得到了支持,毛蚜亚科位于蚜科的基部位置。【结论】两种毛蚜线粒体基因组的基本特征相似,符合蚜虫线粒体基因组的一般特征,两种线粒体基因组的长度差异主要来自控制区长度的不同;系统发育重建支持毛蚜亚科与蚜亚科的单系性,毛蚜亚科位于蚜科较为基部的位置。研究结果为蚜虫类系统发育重建提供了参考。