A hotspot of gene order rearrangement by tandem duplication and random loss in the vertebrate mitochondrial genome

Most reported examples of change in vertebrate mitochondrial (mt) gene order could be explained by a tandem duplication followed by random loss of redundant genes (tandem duplication-random loss [TDRL] model). Under this model of evolution, independent loss of genes arising from a single duplication in an ancestral species are predicted, and remnant pseudogenes expected, intermediate states that may remain in rearranged genomes. However, evidence for this is rare and largely scattered across vertebrate lineages. Here, we report new derived mt gene orders in the vertebrate "WANCY" region of four closely related caecilian amphibians. The novel arrangements found in this genomic region (one of them is convergent with the derived arrangement of marsupials), presence of pseudogenes, and positions of intergenic spacers fully satisfy predictions from the TDRL model. Our results, together with comparative data for the available vertebrate complete mt genomes, provide further evidence that the WANCY genomic region is a hotspot for gene order rearrangements and support the view that TDRL is the dominant mechanism of gene order rearrangement in vertebrate mt genomes. Convergent gene rearrangements are not unlikely in hotspots of gene order rearrangement by TDRL.     

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.     

Bird mitochondrial gene order: insight from 3 warbler mitochondrial genomes

Two main gene orders exist in birds: the ancestral gene order and the remnant control region (CR) 2 gene order. These gene orders differ by the presence of 1 or 2 copies of the CR, respectively. Among songbirds, Oscines were thought to follow the ancestral gene order, with the exception of the lyrebird and Phylloscopus warblers. Here, we determined the complete mitochondrial genome sequence of 3 non-Phylloscopus warblers species and found that the blackcap (Sylvia atricapilla) and the reed warbler (Acrocephalus scirpaceus) have 2 almost identical copies of the CR, whereas the eastern orphean warbler (Sylvia crassirostris) follows the remnant CR 2 gene order. Our results contradict previous studies suggesting that Acrocephalus and most sylvioid warblers exhibit the ancestral gene order. We were able to trace this contradiction to a misidentification of gene order from polymerase chain reaction length determination. We thus suggest that passerine gene order evolution needs to be revised.     

The fucosyltransferase gene family: an amazing summary of the underlying mechanisms of gene evolution

The fucosyltransferase gene family encodes enzymes that transfer fucose in 1,2, 1,3/4 and 1,6 linkages on a large variety of glycans. The most ancient genes harbour a split coding sequence, and encode enzyme that transfer fucose at or near O- and N-peptidic sites (serine, threonine or chitobiose unit). Conversely, the more recent genes have a monoexonic coding sequence, and encode enzymes that transfer fucose at the glycan periphery. All basic mechanisms of gene evolution contribute to this amazing scenario: exon shuffling, transposition, point mutations, and duplication. As typical examples: (i) exon shuffling leads to the ancestral organization of the 1,6 fucosyltransferase gene; (ii) the ancestor of 1,2 fucosyltransferase genes is reshaped by retrotransposition at the same locus; (iii) duplication associated to point mutations leads to the most recent 1,3/4 fucosyltransferase genes.     

Phylogeny, recombination, and mechanisms of stepwise mitochondrial genome reorganization in mantellid frogs from Madagascar

In Malagasy frogs of the family Mantellidae, the genus Mantellais known to possess highly reorganized mitochondrial (mt) genomeswith the following characteristics: 1) some rearranged genepositions, 2) 2 distinct genes and a pseudogene correspondingto the transfer RNA gene for methionine (trnM), and 3) 2 controlregions (CRs) with almost identical nucleotide sequences. Theseunique genomic features were observed concentrated between theduplicated CRs surrounding cytochrome b (cob) and nicotinamideadenine dinucleotide dehydrogenase subunit 2 (cnad2) genes.To elucidate the mechanisms and evolutionary pathway that yieldedthe derived genome condition, we surveyed the reorganized genomicportion for all 12 mantellid genera. Our results show that themt genomes of 7 genera retain the ancestral condition. In contrast,adding to Mantella, 4 genera of the subfamily Mantellinae, Blommersia,Guibemantis, Wakea, and Spinomantis, share several derived genomiccharacters. Furthermore, mt genomes of these mantellines showedadditional structural divergences, resulting in different genomeconditions between them. The high frequency of genomic reorganizationdoes not correlate with nucleotide substitution rate. The encounteredmt genomic conditions also suggest the occurrences of stepwisegene duplication and deletion events during the evolution ofmantellines. Simultaneously, the majority of duplication eventsseems to be mediated by general (homologous) or illegitimaterecombination, and general recombination also plays a role inconcerted sequence evolution between multiple CRs. Consideringour observations and recent conditional evidences, the followingoutlines can be expected for recombination processes in mt genomereorganization. 1) The CR is the "hot spot" of recombination;2) highly frequent recombination between CRs may be mediatedby a replication fork barrier lying in the CR; 3) general recombinationhas a potential to cause gene rearrangement in upstream regionsof multiple CRs as the results of gene conversion and unequalcrossing over processes. Our results also suggest that recombinationactivity is not a direct cause of convergent gene rearrangement;rather, homoplasious gene rearrangement seems to be mediatedby persistence of a copied genomic condition through severallineage splits and subsequent parallel deletions.     

基因倍增研究进展

基因倍增是指DNA片段在基因组中复制出一个或更多的拷贝,这种DNA片段可以是一小段基因组序列、整条染色体,甚至是整个基因组。基因倍增是基因组进化最主要的驱动力之一,是产生具有新功能的基因和进化出新物种的主要原因之一。本文综述了脊椎动物、模式植物和酵母在进化过程中基因倍增研究领域的最新进展,并讨论了基因倍增研究的发展方向。     

Linking movement behaviour to dispersal and divergence in plethodontid salamanders

To better understand the evolutionary and ecological effects of dispersal, there is growing emphasis on the need to integrate direct data on movement behaviour into landscape-scale analyses. However, little is known about the general link between movement behaviour and large-scale patterns of dispersal and gene flow. Likewise, although recent studies suggest that nonrandom, directionally biased movement and dispersal can promote evolutionary divergence, the generality of this mechanism is unknown. We test the hypothesis that directionally biased movement and dispersal by plethodontid salamanders interact with the topography of headwater areas to affect genetic and phenotypic divergence. Movements by Gyrinophilus porphyriticus and Eurycea bislineata show contrasting directional biases: upstream bias in G. porphyriticus and downstream bias in E. bislineata. Consistent with predictions of how these biases interact with slope to affect dispersal and gene flow, genetic distance increased with slope in G. porphyriticus and decreased with slope in E. bislineata over a standardized distance of 1 km along six headwater streams. In both species, phenotypic divergence in relative trunk length was positively related to genetic divergence. These results indicate that landscape-scale patterns of dispersal and gene flow are closely related to movement behaviour in G. porphyriticus and E. bislineata, and underscore the value of information on movement behaviour for predicting and interpreting patterns of dispersal and gene flow in complex landscapes. This study also provides new evidence that directionally biased movement and dispersal can be important sources of intra- and interspecific variation in population divergence, and highlights the value of explicit, a priori predictions in landscape genetic studies.     

The evolution of mitochondrial genomes in modern frogs (Neobatrachia): nonadaptive evolution of mitochondrial genome reorganization

Background

Although mitochondrial (mt) gene order is highly conserved among vertebrates, widespread gene rearrangements occur in anurans, especially in neobatrachians. Protein coding genes in the mitogenome experience adaptive or purifying selection, yet the role that selection plays on genomic reorganization remains unclear. We sequence the mitogenomes of three species of Glandirana and hot spots of gene rearrangements of 20 frog species to investigate the diversity of mitogenomic reorganization in the Neobatrachia. By combing these data with other mitogenomes in GenBank, we evaluate if selective pressures or functional constraints act on mitogenomic reorganization in the Neobatrachia. We also look for correlations between tRNA positions and codon usage.

Results

Gene organization in Glandirana was typical of neobatrachian mitogenomes except for the presence of pseudogene trnS (AGY). Surveyed ranids largely exhibited gene arrangements typical of neobatrachian mtDNA although some gene rearrangements occurred. The correlation between codon usage and tRNA positions in neobatrachians was weak, and did not increase after identifying recurrent rearrangements as revealed by basal neobatrachians. Codon usage and tRNA positions were not significantly correlated when considering tRNA gene duplications or losses. Change in number of tRNA gene copies, which was driven by genomic reorganization, did not influence codon usage bias. Nucleotide substitution rates and d_N/d_S ratios were higher in neobatrachian mitogenomes than in archaeobatrachians, but the rates of mitogenomic reorganization and mt nucleotide diversity were not significantly correlated.

Conclusions

No evidence suggests that adaptive selection drove the reorganization of neobatrachian mitogenomes. In contrast, protein-coding genes that function in metabolism showed evidence for purifying selection, and some functional constraints appear to act on the organization of rRNA and tRNA genes. As important nonadaptive forces, genetic drift and mutation pressure may drive the fixation and evolution of mitogenomic reorganizations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-691) contains supplementary material, which is available to authorized users.     

The mitochondrial genome of the lizard Calotes versicolor and a novel gene inversion in South Asian draconine agamids

A complete mitochondrial DNA (mtDNA) sequence was determinedfor the lizard Calotes versicolor (Reptilia; Agamidae). The16,670-bp genome with notable shorter genes for some protein-codingand tRNA genes had the same gene content as that found in othervertebrates. However, a novel gene arrangement was found inwhich the proline tRNA (trnP) gene is located in the light strandinstead of its typical heavy-strand position, providing thefirst known example of gene inversion in vertebrate mtDNAs.A segment of mtDNA encompassing the trnP gene and its flankinggenes and the control region was amplified and sequenced forvarious agamid taxa to investigate timing and mechanism of thegene inversion. The inverted trnP gene organization was sharedby all South Asian draconine agamids examined but by none ofthe other Asian and African agamids. Phylogenetic analyses includingclock-free Bayesian analyses for divergence time estimationsuggested a single occurrence of the gene inversion on a lineageleading to the draconine agamids during the Paleogene period.This gene inversion could not be explained by the tandem duplication/randomloss model for mitochondrial gene rearrangements. Our availablesequence data did not provide evidence for remolding of thetrnP gene by an anticodon switch in a duplicated tRNA gene.Based on results of sequence comparisons and other circumstantialevidence, we hypothesize that inversion of the trnP gene wasoriginally mediated by a homologous DNA recombination and thatthe de novo gene organization that does not disrupt expressionof mitochondrial genes has been maintained in draconine mtDNAsfor such a long period of time.     

Evolution of the trnF(GAA) gene in Arabidopsis relatives and the brassicaceae family: monophyletic origin and subsequent diversification of a plastidic pseudogene

Recently, we used the 5'-trnL(UAA)-trnF(GAA) region of the chloroplast DNA for phylogeographic reconstructions and phylogenetic analysis among the genera Arabidopsis, Boechera, Rorippa, Nasturtium, and Cardamine. Despite the fact that extensive gene duplications are rare among the chloroplast genome of higher plants, within these taxa the anticodon domain of the trnF(GAA) gene exhibit extensive gene duplications with one to eight tandemly repeated copies in close 5' proximity of the functional gene. Interestingly, even in Arabidopsis thaliana we found six putative pseudogenic copies of the functional trnF gene within the 5'-intergenic trnL-trnF spacer. A reexamination of trnL(UAA)-trnF(GAA) regions from numerous published phylogenetic studies among halimolobine, cardaminoid, and other cruciferous taxa revealed not only extensive trnF gene duplications but also favor the hypothesis about a single origin of trnF pseudogene formation during evolution of the Brassicaceae family 16-21 MYA. Conserved sequence motifs from this tandemly repeated region are codistributed nonrandomly throughout the plastome, and we found some similarities with a DNA sequence duplication in the rps7 gene and its adjacent spacer. Our results demonstrate the potential evolutionary dynamics of a plastidic region generally regarded as highly conserved and probably cotranscribed and, as shown here for several genera among cruciferous plants, greatly characterized by parallel gains and losses of duplicated trnF copies.     

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.     

The complete chloroplast genome sequence of Pelargonium x hortorum: organization and evolution of the largest and most highly rearranged chloroplast genome of land plants

The chloroplast genome of Pelargonium x hortorum has been completely sequenced. It maps as a circular molecule of 217,942 bp and is both the largest and most rearranged land plant chloroplast genome yet sequenced. It features 2 copies of a greatly expanded inverted repeat (IR) of 75,741 bp each and, consequently, diminished single-copy regions of 59,710 and 6,750 bp. Despite the increase in size and complexity of the genome, the gene content is similar to that of other angiosperms, with the exceptions of a large number of pseudogenes, the recognition of 2 open reading frames (ORF56 and ORF42) in the trnA intron with similarities to previously identified mitochondrial products (ACRS and pvs-trnA), the losses of accD and trnT-ggu and, in particular, the presence of a highly divergent set of rpoA-like ORFs rather than a single, easily recognized gene for rpoA. The 3-fold expansion of the IR (relative to most angiosperms) accounts for most of the size increase of the genome, but an additional 10% of the size increase is related to the large number of repeats found. The Pelargonium genome contains 35 times as many 31 bp or larger repeats than the unrearranged genome of Spinacia. Most of these repeats occur near the rearrangement hotspots, and 2 different associations of repeats are localized in these regions. These associations are characterized by full or partial duplications of several genes, most of which appear to be nonfunctional copies or pseudogenes. These duplications may also be linked to the disruption of at least 1 but possibly 2 or 3 operons. We propose simple models that account for the major rearrangements with a minimum of 8 IR boundary changes and 12 inversions in addition to several insertions of duplicated sequence.     

Evolution of the deep-sea gulper eel mitochondrial genomes: large-scale gene rearrangements originated within the eels

Recent studies have demonstrated that deviations from the typical vertebrate mitochondrial gene order are more frequent than initially thought. Such deviations, however, are minor, with inversions and/or translocations of a few genes being involved and tandem duplication of the gene regions followed by deletions of genes having been invoked as mechanisms originating in such novel gene order. During the course of molecular phylogenetic studies on the Elopomorpha (eels and their allies), we found that mitochondrial genomes (mitogenomes) from the two deep-sea gulper eels, Eurypharynx pelecanoides (Eurypharyngidae) and Saccopharynx lavenbergi (Saccopharyngidae), exhibit an identical gene order which greatly differs from that of any other vertebrates. Phylogenetic analysis using the mitogenomic data from 59 species of fish not only confirmed a single origin of such a gene order with confidence but also indicated that it had been derived from the typical vertebrate gene order. Detailed comparisons of the gulper eel gene order with that of typical vertebrates suggested that occurrence of a single step, large-scale duplication of gene region extending >12 kb, followed by deletions of genes in a common ancestor of the two species, most parsimoniously accounts for this unusual gene arrangement.     

Multiple sequence rearrangements accompanying the duplication of a tRNAPro gene in wheat mitochondrial DNA

In the course of isolating tRNA genes from wheat mtDNA, we have found the same tRNA^Pro gene in two different Hind III restriction fragments, H-P1 (0.7 kbp) and H-P2 (1.7 kbp). Sequences immediately flanking these duplicate genes are closely related, although not identical; sequence comparisons suggest that multiple rearrangements have occurred in the vicinity of the H-P2 tRNA^Pro gene, relative to the H-P1 version. The chimeric nature of H-P2 is emphasized by the presence of sequences that are also found upstream of the wheat mitochondrial 26S rRNA gene, as well as sequences derived from chloroplast DNA. Comparison of H-P2 with H-P1 plus upstream sequences provides some insight into possible molecular events that might have generated H-P2. In particular, such comparisons suggest a model in which the homologous sequences in H-P2 are seen to be derived from H-P1 plus upstream sequences as a result of an intragenomic, site-specific rearrangement event, followed by amplification of the product, its fixation in the mitochondrial genome, and subsequent sequence divergence (single base changes as well as insertions/deletions of up to 50 nucleotides). The results reported here implicate particular primary sequence motifs in certain of the rearrangements that characterize H-P2.     

Characterization of the wheat mitochondrial orf25 gene

The wheat mitochondrial orf25 nucleotide sequence of 576 pb has been determined. Its derived protein sequence shares 88% and 75% amino acid identity with those of maize and tobacco mitochondria, respectively. The wheat and tobacco orf25 sequences lack four inserts, of 6 bp to 36 bp, that are present in the maize homologue. The wheat orf25 gene is actively transcribed and is preceded by a regulatory sequence block very similar to those located upstream of the wheat coxII and atp6 genes. Our observations support the view that orf25 sequences encode a functional polypeptide in plant mitochondria.     

Plant Mitochondrial Genome Evolution and Cytoplasmic Male Sterility

Mitochondria are responsible for providing energy currency to life processes in the molecular form of ATP and are therefore typically referred to as the power factories of cells. Plant mitochondria are also relevant to the common phenomenon of cytoplasmic male sterility, which is agronomically important in various crop species. Cytoplasmic male sterility (CMS) is a complex trait that may be influenced by patterns of mitochondrial genome evolution, and by intergenomic gene transfer among the organellar and nuclear compartments of plant cells. Here, we review patterns and processes that shape plant mitochondrial genomes, some relevant interactions between organelles, and the general features shared by the majority of cytoplasmic male-sterile genes in plants to further the goal of understanding CMS.     

Codon usage and bias in mitochondrial genomes of parasitic platyhelminthes

Sequences of the complete protein-coding portions of the mitochondrial (mt) genome were analysed for 6 species of cestodes (including hydatid tapeworms and the pork tapeworm) and 5 species of trematodes (blood flukes and liver- and lung-flukes). A near-complete sequence was also available for an additional trematode (the blood fluke Schistosoma malayensis). All of these parasites belong to a large flatworm taxon named the Neodermata. Considerable variation was found in the base composition of the protein-coding genes among these neodermatans. This variation was reflected in statistically-significant differences in numbers of each inferred amino acid between many pairs of species. Both convergence and divergence in nucleotide, and hence amino acid, composition was noted among groups within the Neodermata. Considerable variation in skew (unequal representation of complementary bases on the same strand) was found among the species studied. A pattern is thus emerging of diversity in the mt genome in neodermatans that may cast light on evolution of mt genomes generally.