Gene conversion and concerted evolution in bacterial genomes

Gene conversion is defined as the non-reciprocal transfer of information between homologous sequences. Despite methodological problems to establish non-reciprocity, gene conversion has been demonstrated in a wide variety of bacteria. Besides examples of high-frequency reversion of mutations in repeated genes, gene conversion in bacterial genomes has been implicated in concerted evolution of multigene families. Gene conversion also has a prime importance in the generation of antigenic variation, an interesting mechanism whereby some bacterial pathogens are able to avoid the host immune system. In this review, we analyze examples of bacterial gene conversion (some of them spawned from the current genomic revolution), as well as the molecular models that explain gene conversion and its association with crossovers.     

Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit 'bizarre' secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading 'pseudogenes', even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.     

Gene rearrangements in gekkonid mitochondrial genomes with shuffling,loss, and reassignment of tRNA genes

Background

Vertebrate mitochondrial genomes (mitogenomes) are 16–18 kbp double-stranded circular DNAs that encode a set of 37 genes. The arrangement of these genes and the major noncoding region is relatively conserved through evolution although gene rearrangements have been described for diverse lineages. The tandem duplication-random loss model has been invoked to explain the mechanisms of most mitochondrial gene rearrangements. Previously reported mitogenomic sequences for geckos rarely included gene rearrangements, which we explore in the present study.

Results

We determined seven new mitogenomic sequences from Gekkonidae using a high-throughput sequencing method. The Tropiocolotes tripolitanus mitogenome involves a tandem duplication of the gene block: tRNA^Arg, NADH dehydrogenase subunit 4L, and NADH dehydrogenase subunit 4. One of the duplicate copies for each protein-coding gene may be pseudogenized. A duplicate copy of the tRNA^Arg gene appears to have been converted to a tRNA^Gln gene by a C to T base substitution at the second anticodon position, although this gene may not be fully functional in protein synthesis. The Stenodactylus petrii mitogenome includes several tandem duplications of tRNA^Leu genes, as well as a translocation of the tRNA^Ala gene and a putative origin of light-strand replication within a tRNA gene cluster. Finally, the Uroplatus fimbriatus and U. ebenaui mitogenomes feature the apparent loss of the tRNA^Glu gene from its original position. Uroplatus fimbriatus appears to retain a translocated tRNA^Glu gene adjacent to the 5’ end of the major noncoding region.

Conclusions

The present study describes several new mitochondrial gene rearrangements from Gekkonidae. The loss and reassignment of tRNA genes is not very common in vertebrate mitogenomes and our findings raise new questions as to how missing tRNAs are supplied and if the reassigned tRNA gene is fully functional. These new examples of mitochondrial gene rearrangements in geckos should broaden our understanding of the evolution of mitochondrial gene arrangements.

Electronic supplementary material

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

Two complete mitochondrial genomes in Scolopendra and a comparative analysis of tRNA rearrangements in centipedes

Background

Centipedes are one of the oldest terrestrial arthropods belonging to the sub phylum Myriapoda. With the expansion of our understanding of the application of the two centipedes Scolopendra morsitans and Scolopendra hainanum, belonging to the order Scolopendromorpha, an exhaustive classification was required. Although consensus has been reached on the phylogeny of Chilopoda based on morphological traits, recent analyses based on molecular data exhibited differences in results.

Methods and results

The mitochondrial genome sequences of S. morsitans and S. hainanum were obtained by next-generation sequencing. S. morsitans contains 13 PCGs, two rRNAs, 11 tRNAs, and one CR. whereas S. hainanum contains 12 PCGs, of which ATP8 remains unpredicted, two rRNAs, 14 tRNAs, and one CR. An obvious tRNA rearrangement was found in the genus Scolopendra. S. morsitans exhibited a loss of trnW, trnC, trnI, trnK, trnD, trnA, trnN, trnQ, trnF, trnT, trnS, trnL, and trnV, and a repeat of trnR and trnL. S. hainanum exhibited a loss of trnQ, trnC, trnW, trnI, trnD, trnQ, trnP, and trnV. Phylogenetic analyses of centipedes based on 12 PCGs supported the sister relationship between the orders Geophilomorpha and Lithobiomorpha and a close relationship between Scolopendra dehaani and S. hainanum.

Conclusions

The new mitogenomes determined in this study provide new genomic resources for gene rearrangements and contribute to the understanding of the evolution of gene rearrangement in Chilopoda.

     

Gene conversion as a mechanism for divergence of a chloroplast tRNA gene inserted in the mitochondrial genome of Brassica oleracea.

We have characterized a 1.7 kb sequence, containing a tRNA Leu2 gene shared by the ct and mt genomes of Brassica oleracea. The two sequences are completely homologous except in two short regions where two distinct gene conversion events have occurred between two sets of direct repeats leading to the insertion of 5 bp in the T loop of the mt copy of the ct gene. This is the first evidence that gene conversion represents the initial evolutionary step in inactivation of transferred ct genes in the mt genome. We also indicate that organelle DNA transfer by organelle fusion is an ongoing process which could be useful in genetic engineering.     

Short repetitive sequences in green algal mitochondrial genomes: potential roles in mitochondrial genome evolution

Current data on green algal mitochondrial genomes suggest an unexpected dichotomy within the group with respect to genome structure, organization, and sequence affiliations. The present study suggests that there is a correlation between this dichotomy on one hand and the differences in the abundance, base composition, and distribution of short repetitive sequences we observed among green algal mitochondrial genomes on the other. It is conceivable that the accumulation of GC- rich short repeated sequences in the Chlamydomonas-like but not Prototheca-like mitochondrial genomes might have triggered evolutionary events responsible for the distinct series of evolutionary changes undergone by the two green algal mitochondrial lineages. The similarity in base composition, nucleotide sequence, abundance, and mode of organization we observed between the short repetitive sequences present in Chlamydomonas-like mitochondrial genomes on one hand and fungal and vertebrate homologs on the other might extend to some of the roles that the short repetitive sequences have been shown to have in the latter. Potential involvements we propose for the short repetitive sequences in the evolution of Chlamydomonas-like mitochondrial genomes include fragmentation and scrambling of the ribosomal-RNA-coding regions, extensive gene rearrangements, coding-region deletions, surrogate origins of replication, and chromosomal linearization.      

Tracing plant Mitochondrial DNA evolution: rearrangements of the ancient mitochondrial gene cluster trnA-trnT-nad7 in liverwort phylogeny

Whereas frequent recombination characterizes flowering plant mitochondrial genomes, some mitochondrial gene arrangements may, in contrast, be conserved between streptophyte algae and early land plant clades (bryophytes). Here we explore the evolutionary fate of the mitochondrial gene arrangement trnA-trnT-nad7, which is conserved among the alga Chara, the moss Physcomitrella, and the liverwort Marchantia, although trnT is inverted in orientation in the latter. Surprisingly, we now find that the Chara-type gene arrangement is generally conserved in mosses, but that trnT is lacking between trnA and nad7 in all simple-thalloid and leafy (jungermanniid) liverworts. The ancient gene continuity trnA-trnT-nad7 is, however, conserved in Blasia, representing the sister lineage to all other complex-thalloid (marchantiid) liverworts. The recombinogenic insertion of short sequence stretches, including nad5 and rps7 pseudogene fragments copied from elsewhere in the liverwort mtDNA, likely mediated a subsequent inversion of trnT and flanking sequences in a basal grade of marchantiid liverworts, which was then followed by an independent secondary loss of trnT in derived marchantiid taxa later in evolution. In contrast to the previously observed extreme degree of coding sequence conservation and the assumed absence of active recombination in Marchantia mtDNA, this now reveals a surprisingly dynamic evolution of marchantiid liverwort mitochondrial genomes.     

Analysis of a drosophila tRNA gene cluster: two tRNALeu genes contain intervening sequences

A recombinant DNA phage containing a cluster of Drosophila melanogaster tRNA genes has been isolated and analyzed. The insert of this phage has been mapped by in situ hybridization to chromosomal region 50AB, a known tRNA site. Nucleotide sequencing of the entire Drosophila tRNA coding region reveals seven tRNA genes spanning 2.5 kb of chromosomal DNA. This cluster is separated from other tRNA regions on the chromosome by at least 2.7 kb on one side, and 9.6 kb on the other. Two tRNA genes are nearly identical and contain intervening sequences of length 38 and 45 bases, respectively, in the anticodon loop. These two genes are assigned to be tRNALeu genes because of significant sequence homology with yeast tRNA3Leu, and secondary structure homology with yeast tRNA3Leu intervening sequence. In addition, an 8 base sequence (AAAAUCUU) is conserved in the same location in the intervening sequences of Drosophila tRNALeu genes and a yeast tRNA3Leu gene. Similar sequenes occur in all other tRNAs containing intervening sequences. The remaining five genes are identical tRNAIle genes, which are also identical to a tRNAIle gene from chromosomal region 42A. The 5' flanking regions are only weakly homologous, but each set of isoacceptors contains short regions of strong homology approximately 20 nucleotides preceding the tRNA coding sequences: GCNTTTTG preceding tRNAIle genes; and GANTTTGG preceding tRNALeu genes. The genes are irregularly distributed on both DNA strands; spacing regions are divergent in sequence and length.     

The complete chloroplast and mitochondrial genome sequences of Boea hygrometrica: insights into the evolution of plant organellar genomes

The complete nucleotide sequences of the chloroplast (cp) and mitochondrial (mt) genomes of resurrection plant Boea hygrometrica (Bh, Gesneriaceae) have been determined with the lengths of 153,493 bp and 510,519 bp, respectively. The smaller chloroplast genome contains more genes (147) with a 72% coding sequence, and the larger mitochondrial genome have less genes (65) with a coding faction of 12%. Similar to other seed plants, the Bh cp genome has a typical quadripartite organization with a conserved gene in each region. The Bh mt genome has three recombinant sequence repeats of 222 bp, 843 bp, and 1474 bp in length, which divide the genome into a single master circle (MC) and four isomeric molecules. Compared to other angiosperms, one remarkable feature of the Bh mt genome is the frequent transfer of genetic material from the cp genome during recent Bh evolution. We also analyzed organellar genome evolution in general regarding genome features as well as compositional dynamics of sequence and gene structure/organization, providing clues for the understanding of the evolution of organellar genomes in plants. The cp-derived sequences including tRNAs found in angiosperm mt genomes support the conclusion that frequent gene transfer events may have begun early in the land plant lineage.     

Homologous recombination between highly diverged mitochondrial sequences: examples from maternally and paternally transmitted genomes

Homologous recombination is restricted to sequences of low divergence. This is attributed to the mismatch repairing system (MMR), which does not allow recombination between sequences that are highly divergent. This acts as a safeguard against recombination between nonhomologous sequences that could result in genome imbalance. Here, we report recombination between maternal and paternal mitochondrial genomes of the sea mussel, whose sequences differ by >20%. We propose that the strict maternal inheritance of the animal mitochondrial DNA and the ensuing homoplasmy has relieved the MMR system of the animal mitochondrion from the pressure to tolerate recombination only among sequences with a high degree of similarity.     

Complete DNA sequences of the mitochondrial genomes of the pathogenic yeasts Candida orthopsilosis and Candida metapsilosis: insight into the evolution of linear DNA genomes from mitochondrial telomere mutants

We determined complete mitochondrial DNA sequences of the two yeast species, Candida orthopsilosis and Candida metapsilosis, and compared them with the linear mitochondrial genome of their close relative, C.parapsilosis. Mitochondria of all the three species harbor compact genomes encoding the same set of genes arranged in the identical order. Differences in the length of these genomes result mainly from the presence/absence of introns. Multiple alterations were identified also in the sequences of the ribosomal and transfer RNAs, and proteins. However, the most striking feature of C.orthopsilosis and C.metapsilosis is the existence of strains differing in the molecular form of the mitochondrial genome (circular-mapping versus linear). Their analysis opens a unique window for understanding the role of mitochondrial telomeres in the stability and evolution of molecular architecture of the genome. Our results indicate that the circular-mapping mitochondrial genome derived from the linear form by intramolecular end-to-end fusions. Moreover, we suggest that the linear mitochondrial genome evolved from a circular-mapping form present in a common ancestor of the three species and, at the same time, the emergence of mitochondrial telomeres enabled the formation of linear monomeric DNA forms. In addition, comparison of isogenic C.metapsilosis strains differing in the form of the organellar genome suggests a possibility that, under some circumstances, the linearity and/or the presence of telomeres provide a competitive advantage over a circular-mapping mitochondrial genome.     

Complete sequences of mitochondria genomes of Aedes aegypti and Culex quinquefasciatus and comparative analysis of mitochondrial DNA fragments inserted in the nuclear genomes

We present complete sequences of the mitochondrial genomes for two important mosquitoes, Aedes aegypti and Culex quinquefasciatus, that are major vectors of dengue virus and lymphatic filariasis, respectively. The A. aegypti mitochondrial genome is 16,655 bp in length and that of C. quinquefasciatus is 15,587 bp, yet both contain 13 protein coding genes, 22 transfer RNA (tRNA) genes, one 12S ribosomal RNA (rRNA) gene, one 16S rRNA gene and a control region (CR) in the same order. The difference in the genome size is due to the difference in the length of the control region. We also analyzed insertions of nuclear copies of mtDNA-like sequences (NUMTs) in a comparative manner between the two mosquitoes. The NUMT sequences occupy ~0.008% of the A. aegypti genome and ~0.001% of the C. quinquefasciatus genome. Several NUMTs were found localized in the introns of predicted protein coding genes in both genomes (32 genes in A. aegypti but only four in C. quinquefasciatus). None of these NUMT-containing genes had an ortholog between the two species or had paralogous copies within a genome that was also NUMT-containing. It was further observed that the NUMT-containing genes were relatively longer but had lower GC content compared to the NUMT-less paralogous copies. Moreover, stretches of homologies are present among the genic and non-genic NUMTs that may play important roles in genomic rearrangement of NUMTs in these genomes. Our study provides new insights on understanding the roles of nuclear mtDNA sequences in genome complexities of these mosquitoes.     

Multiple rearrangements in mitochondrial genomes of Isopoda and phylogenetic implications

In this study, we analyse the evolutionary dynamics and phylogenetic implications of gene order rearrangements in five newly sequenced mitochondrial (mt) genomes and four published mt genomes of isopod crustaceans. The sequence coverage is nearly complete for four of the five newly sequenced species, with only the control region and some tRNA genes missing, while in Janira maculosa only two thirds of the genome could be determined. Mitochondrial gene order in isopods seems to be more plastic than that in other crustacean lineages, making all nine known mt gene orders different. Especially the asellote Janira is characterized by many autapomorphies. The following inferred ancestral isopod mt gene order exists slightly modified in modern isopods: nad1, tnrL1, rrnS, control region, trnS1, cob, trnT, nad5, trnF. We consider the inferred gene translocation events leading to gene rearrangements as valuable characters in phylogenetic analyses. In this first study covering major isopod lineages, potential apomorphies were identified, e.g., a shared relative position of trnR in Valvifera. We also report one of the first findings of homoplasy in mitochondrial gene order, namely a shared relative position of trnV in unrelated isopod lineages. In addition to increased taxon sampling secondary structure, modification in tRNAs and GC-skew inversion may be potentially fruitful subjects for future mt genome studies in a phylogenetic context.     

The major resistance gene cluster in lettuce is highly duplicated and spans several megabases.

At least 10 Dm genes conferring resistance to the oomycete downy mildew fungus Bremia lactucae map to the major resistance cluster in lettuce. We investigated the structure of this cluster in the lettuce cultivar Diana, which contains Dm3. A deletion breakpoint map of the chromosomal region flanking Dm3 was saturated with a variety of molecular markers. Several of these markers are components of a family of resistance gene candidates (RGC2) that encode a nucleotide binding site and a leucine-rich repeat region. These motifs are characteristic of plant disease resistance genes. Bacterial artificial chromosome clones were identified by using duplicated restriction fragment length polymorphism markers from the region, including the nucleotide binding site-encoding region of RGC2. Twenty-two distinct members of the RGC2 family were characterized from the bacterial artificial chromosomes; at least two additional family members exist. The RGC2 family is highly divergent; the nucleotide identity was as low as 53% between the most distantly related copies. These RGC2 genes span at least 3.5 Mb. Eighteen members were mapped on the deletion breakpoint map. A comparison between the phylogenetic and physical relationships of these sequences demonstrated that closely related copies are physically separated from one another and indicated that complex rearrangements have shaped this region. Analysis of low-copy genomic sequences detected no genes, including RGC2, in the Dm3 region, other than sequences related to retrotransposons and transposable elements. The related but divergent family of RGC2 genes may act as a resource for the generation of new resistance phenotypes through infrequent recombination or unequal crossing over.     

Large and small rearrangements in the evolution of prokaryotic genomes

Relative frequencies of large and small genome rearrangements (inversions and transpositions) in the evolution of prokaryotic genomes can be evaluated using the ratio between the index S (the ratio of the number of identical pairs of neighboring genes in two genomes to the total number of genes in the sample of interest) and 1–6L/n, where L is the mean difference in intergenic distances and n is the number of genes in the sample. The S value uniformly decreases with the fixation of genome rearrangements, while the decrease rate of 1–6L/n is determined by the rearrangement size. Specifically, large inversions and transpositions lead to a dramatic decrease in the index value, while small rearrangements result in an insignificant decrease. The ratio between these indices was computed for twenty pairs of closely related species belonging to different groups of bacteria and archaea. The pairs examined strongly differed in the relative frequency of large and small rearrangements. However, computer simulation showed that the total variation can be reproduced with the same input parameters of the model. This means that the differences observed can be stochastic and can be interpreted without assuming different mechanisms and factors of genome rearrangements for different groups of prokaryotes. Relative frequencies of large and small rearrangements displayed no noticeable correlations with taxonomic position, total rate of rearrangement fixation, habitation conditions, and the abundance of transposons and repetitive sequences. It is suggested that, in some cases, phage activity increases the frequency of large genome rearrangements.     

Large and small rearrangements in the evolution of prokaryotic genomes

Relative frequencies of large and small genome rearrangements (inversions and transpositions) in the evolution of prokaryotic genomes can be evaluated using the ratio between the index S (the ratio of the number of identical pairs of neighboring genes in two genomes to the total number of genes in the sample of interest) and 1 - 6 x L/n, where L is the mean difference in intergenic distances and n is the number of genes in the sample. The S value uniformly decreases with the fixation of genome rearrangements, while the decrease rate of I - 6 x L/n is determined by the rearrangement size. Specifically, large inversions and transpositions lead to a dramatic decrease in the index value, while small rearrangements result in an insignificant decrease. The ratio between these indices was computed for twenty pairs of closely related species belonging to different groups of bacteria and archaea. The pairs examined strongly differed in the relative frequency of large and small rearrangements. However, computer simulation showed that the total variation can be reproduced with the same input parameters of the model. This means that the differences observed can be stochastic and can be interpreted without assuming different mechanisms and factors of genome rearrangements for different groups of prokaryotes. Relative frequencies of large and small rearrangements displayed no noticeable correlations with taxonomic position, total rate of rearrangement fixation, habitation conditions, and the abundance of transposons and repetitive sequences. It is suggested that, in some cases, phage activity increases the frequency of large genome rearrangements.     

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.