Tandem Duplication-random Loss Model Research Articles

Gene Rearrangements in the Mitochondrial Genome of Gonatopsis borealis and Onychoteuthis compacta Reveal Their Phylogenetic Implications for Oegopsida.

The complete mitochondrial genome provides crucial information for comprehending gene rearrangement, molecular evolution, and phylogenetic analysis. Here, we have determined the complete mitogenome sequence of Gonatopsis borealis and Onychoteuthis compacta for the first time. Their genome sizes were 20,148 bp and 20,491 bp, respectively, including 18 protein-coding genes, COI-COIII, ATP6, and ATP8 are duplicated, 23 transfer RNA genes, and 2 ribosomal RNA (rRNA) genes (12S and 16S rRNA). Specifically, the overall A+T content is 70.69% and 72.67%. It shows a significant AT bias. The whole mitogenomes indicate positive AT skew (0.070 and 0.062). Furthermore, the gene order has been rearranged within Oegopsida. The tandem duplication random loss model was determined as most likely to explain the observed gene rearrangements. Phylogenetic analysis was performed, and the result tree was found to be consistent with the morphological identification classification. Estimation of divergence time for 35 species showed that the main differentiation of Oegopsida occurred in 140.70Mya. These results will help to better understand the gene rearrangements and evolution of G. borealis and O. compacta and lay a foundation for further phylogeny genetic studies of Oegopsida.

Evolution of the Noncoding Features of Sea Snake Mitochondrial Genomes within Elapidae.

Mitochondrial genomes of four elapid snakes (three marine species [Emydocephalus ijimae, Hydrophis ornatus, and Hydrophis melanocephalus], and one terrestrial species [Sinomicrurus japonicus]) were completely sequenced by a combination of Sanger sequencing, next-generation sequencing and Nanopore sequencing. Nanopore sequencing was especially effective in accurately reading through long tandem repeats in these genomes. This led us to show that major noncoding regions in the mitochondrial genomes of those three sea snakes contain considerably long tandem duplications, unlike the mitochondrial genomes previously reported for same and other sea snake species. We also found a transposition of the light-strand replication origin within a tRNA gene cluster for the three sea snakes. This change can be explained by the Tandem Duplication—Random Loss model, which was further supported by remnant intervening sequences between tRNA genes. Mitochondrial genomes of true snakes (Alethinophidia) have been shown to contain duplicate major noncoding regions, each of which includes the control region necessary for regulating the heavy-strand replication and transcription from both strands. However, the control region completely disappeared from one of the two major noncoding regions for two Hydrophis sea snakes, posing evolutionary questions on the roles of duplicate control regions in snake mitochondrial genomes. The timing and molecular mechanisms for these changes are discussed based on the elapid phylogeny.

The first mitochondrial genome of the genus Exhippolysmata (Decapoda: Caridea: Lysmatidae), with gene rearrangements and phylogenetic associations in Caridea

The complete mitochondrial genome (mitogenome) of animals can provide useful information for evolutionary and phylogenetic analyses. The mitogenome of the genus Exhippolysmata (i.e., Exhippolysmata ensirostris) was sequenced and annotated for the first time, its phylogenetic relationship with selected members from the infraorder Caridea was investigated. The 16,350 bp mitogenome contains the entire set of 37 common genes. The mitogenome composition was highly A + T biased at 64.43% with positive AT skew (0.009) and negative GC skew (− 0.199). All tRNA genes in the E. ensirostris mitogenome had a typical cloverleaf secondary structure, except for trnS1 (AGN), which appeared to lack the dihydrouridine arm. The gene order in the E. ensirostris mitogenome was rearranged compared with those of ancestral decapod taxa, the gene order of trnL2-cox2 changed to cox2-trnL2. The tandem duplication-random loss model is the most likely mechanism for the observed gene rearrangement of E. ensirostris. The ML and BI phylogenetic analyses place all Caridea species into one group with strong bootstrap support. The family Lysmatidae is most closely related to Alpheidae and Palaemonidae. These results will help to better understand the gene rearrangements and evolutionary position of E. ensirostris and lay a foundation for further phylogenetic studies of Caridea.

Unique Duplication of trnN in Odontoptilum angulatum (Lepidoptera: Pyrginae) and Phylogeny within Hesperiidae.

Simple SummaryPyrginae is one of the major groups of the Hesperiidae, and has some particular characteristics. The annotated complete mitogenome from this subfamily is reported here. The gene order of the new mitogenome with the duplication of trnN differs from the typical Lepidoptera-specific arrangement and is unique to Hesperiidae. The presence of a pseudo gene in the mt genome of Odontoptilum angulatum supports the duplication of trnN, following the TDRL model. Comparison of the newly generated mitogenome of Odontoptilum angulatum to all available mitochondrial genomes of other rearranged Pyrginae species revealed that the condition of Odontoptilum angulatum is of independent origin. Therefore, we hypothesize that the gene block trnN–trnS1–trnE is the hot spot of gene rearrangement in the Tagiadini of Pyrginae. Phylogenetic analyses based on 13 protein-coding genes and entire RNA genes of mitogenomes show the monophyly of Pyrginae.To explore the variation and relationship between gene rearrangement and phylogenetic effectiveness of mitogenomes among lineages of the diversification of the tribe Tagiadini in the subfamily Pyrginae, we sequenced the complete mitogenome of Odontoptilum angulatum. The genome is 15,361 bp with the typical 37 genes, a large AT-rich region and an additional trnN (trnN2), which is completely identical to trnN (sequence similarity: 100%). The gene order differs from the typical Lepidoptera-specific arrangement and is unique to Hesperiidae. The presence of a “pseudo-trnS1” in the non-coding region between trnN1 and trnN2 supports the hypothesis that the presence of an extra trnN can be explained by the tandem duplication-random loss (TDRL) model. Regarding the phylogenetic analyses, we found that the dataset comprising all 37 genes produced the highest node support, as well as a monophyly of Pyrginae, indicating that the inclusion of RNAs improves the phylogenetic signal. Relationships among the subfamilies in Hesperiidae were also in general agreement with the results of previous studies. The monophyly of Tagiadini is strongly supported. Our study provides a new orientation for application of compositional and mutational biases of mitogenomes in phylogenetic analysis of Tagiadini and even all Hesperiidae based on larger taxon sampling in the future.

Gene rearrangements in the mitochondrial genome of Chiromantes eulimene (Brachyura: Sesarmidae) and phylogenetic implications for Brachyura

Complete mitochondrial genome (mitogenome) provides important information for better understanding of gene rearrangement, molecular evolution and phylogenetic analysis. Here we determined the complete mitogenome sequence of Chiromantes eulimene (Brachyura: Sesarmidae) for the first time. The total length is 15,894 bp and includes 13 protein-coding genes (PCGs), 22 transfer RNAs, two ribosomal RNAs, as well as a putative control region. The genome composition is highly A + T biased (75.5%), and exhibits a negative AT-skew (−0.017) and GC-skew (−0.206). All of the 13 PCGs are initiated by the start codon ATN, with an exception (GTG) in ND1. The typical stop codon (TAA or TAG) is detected in ten PCGs, whereas the remaining three PCGs (COI, COII and Cyt b) terminate by an incomplete T. The gene order in C. eulimene mitogenome was rearranged compared with that of the ancestor of Decapoda. The gene order of F-ND5-H changed to H-F-ND5. Like other sesarmid crabs, the I-Q-M gene cluster in the pancrustacean ground pattern became Q-I-M order in C. eulimene genome. Tandem duplication-random loss model and slipped-strand mispairing mechanism are determined as most likely to explain the observed gene rearrangements. Phylogenetic analysis places all Sesarmidae species into one group. Almost all families except Xanthidae, Gecarcinidae and Homolidae form a monophyletic clade and the polyphyly of Eriphioidea, Ocypodoidea and Grapsoidea is well supported. These results will help to better understand the gene rearrangements and evolutionary position of C. eulimene and lay a foundation for further phylogenetic studies of Brachyura.

Deep transcriptome profiling sheds light on key players in nucleus implantation induced immune response in the pearl oyster Pinctada martensii

Complete mitochondrial genome (mitogenome) provides important information for better understanding of gene rearrangement, molecular evolution and phylogenetic analysis. Here we determined the complete mitogenome sequence of Chiromantes eulimene (Brachyura: Sesarmidae) for the first time. The total length is 15,894 bp and includes 13 protein-coding genes (PCGs), 22 transfer RNAs, two ribosomal RNAs, as well as a putative control region. The genome composition is highly A + T biased (75.5%), and exhibits a negative AT-skew (−0.017) and GC-skew (−0.206). All of the 13 PCGs are initiated by the start codon ATN, with an exception (GTG) in ND1. The typical stop codon (TAA or TAG) is detected in ten PCGs, whereas the remaining three PCGs (COI, COII and Cyt b) terminate by an incomplete T. The gene order in C. eulimene mitogenome was rearranged compared with that of the ancestor of Decapoda. The gene order of F-ND5-H changed to H-F-ND5. Like other sesarmid crabs, the I-Q-M gene cluster in the pancrustacean ground pattern became Q-I-M order in C. eulimene genome. Tandem duplication-random loss model and slipped-strand mispairing mechanism are determined as most likely to explain the observed gene rearrangements. Phylogenetic analysis places all Sesarmidae species into one group. Almost all families except Xanthidae, Gecarcinidae and Homolidae form a monophyletic clade and the polyphyly of Eriphioidea, Ocypodoidea and Grapsoidea is well supported. These results will help to better understand the gene rearrangements and evolutionary position of C. eulimene and lay a foundation for further phylogenetic studies of Brachyura.

A transfer RNA gene rearrangement in the lepidopteran mitochondrial genome

Gene arrangements in the mitochondrial genomes (mitogenomes) of insects are conserved across the major lineages, but can be rearranged within derived groups and may provide valuable phylogenetic characters. In this study, we sequenced the entire mitogenome of Parasa consocia, a moth of the family Limacodidae (Lepidoptera: Zygaenoidea). Compared with other lepidopterans and ancestral insects, the P. consocia mitogenome features a transfer RNA gene arrangement novel among lepidopterans between the ND3 and ND5 genes: RANSEF (the underline signifies an inverted gene), which differs from the ARNSEF arrangement of ancestral insects. This rearrangement can be explained by the tandem duplication-random loss model. We inferred a phylogenetic hypothesis for the lepidopteran superfamily based on mitochondrial amino-acid sequences using the Bayesian-inference and maximum-likelihood methods. Our results showed that P. consocia belongs to the Zygaenoidea superfamily and supported the following phylogenetic relationship: Yponomeutoidea + (Tortricoidea + Zygaenoidea + (Papilionoidea + (Pyraloidea + (Noctuoidea + (Geometroidea + Bombycoidea)))))). Comparative analyses indicated that mitogenomes are a useful phylogenetic tool at the subfamily level within the order Lepidoptera. Our findings also suggest that mitogenomes are likely to represent a valuable tool for systematics in other groups of lepidopterans.

Complete mitochondrial genomes of two gelechioids, Mesophleps albilinella and Dichomeris ustalella (Lepidoptera: Gelechiidae), with a description of gene rearrangement in Lepidoptera

We sequenced the entire mitochondrial genome (mitogenome) of two gelechioids, Mesophleps albilinella and Dichomeris ustalella, and compared their genome organization and sequence composition to those of available gelechioid mitogenomes for an enhanced understanding of Gelechioidea genomic characteristics. We compared all available lepidopteran mitogenome arrangements, including that of M. albilinella, which is unique in Gelechioidea, to comprehend the extensiveness and mechanisms of gene rearrangement in Lepidoptera. The genomes of M. albilinella and D. ustalella are 15,274 and 15,410bp in size, respectively, with the typical sets of mitochondrial (mt) genes. The COI gene begins with CGA (arginine) in all sequenced gelechioids, including M. albilinella and D. ustalella, reinforcing the feature as a synapomorphic trait, at least in the Gelechioidea. Each 353- and 321-bp long A+T-rich region of M. albilinella and D. ustalella contains one (D. ustalella) or two (M. albilinella) tRNA-like structures. The M. albilinella mitogenome has a unique gene arrangement among the Gelechioidea: ARNESF (the underline signifies an inverted gene) at the ND3 and ND5 junction, as opposed to the ARNSEF that is found in ancestral insects. An extensive search of available lepidopteran mitogenomes, including that of M. albilinella, turned up six rearrangements that differ from those of ancestral insects. Most of the rearrangements can be explained by the tandem duplication-random loss model, but inversion, which requires recombination, is also found in two cases, including M. albilinella. Excluding the MIQ rearrangement at the A+T-rich region and ND2 junction, which is found in nearly all Ditrysia, most of the remaining rearrangements found in Lepidoptera appear to be independently derived in that they are automorphic at several taxonomic scales, although current mitogenomic data are limited, particularly for congeneric data.

The complete mitochondrial genomes of sixteen ardeid birds revealing the evolutionary process of the gene rearrangements.

BackgroundThe animal mitochondrial genome is generally considered to be under selection for both compactness and gene order conservation. As more mitochondrial genomes are sequenced, mitochondrial duplications and gene rearrangements have been frequently identified among diverse animal groups. Although several mechanisms of gene rearrangement have been proposed thus far, more observational evidence from major taxa is needed to validate specific mechanisms. In the current study, the complete mitochondrial DNA of sixteen bird species from the family Ardeidae was sequenced and the evolution of mitochondrial gene rearrangements was investigated. The mitochondrial genomes were then used to review the phylogenies of these ardeid birds.ResultsThe complete mitochondrial genome sequences of the sixteen ardeid birds exhibited four distinct mitochondrial gene orders in which two of them, named as “duplicate tRNAGlu–CR” and “duplicate tRNAThr–tRNAPro and CR”, were newly discovered. These gene rearrangements arose from an evolutionary process consistent with the tandem duplication - random loss model (TDRL). Additionally, duplications in these gene orders were near identical in nucleotide sequences within each individual, suggesting that they evolved in concert. Phylogenetic analyses of the sixteen ardeid species supported the idea that Ardea ibis, Ardea modesta and Ardea intermedia should be classified as genus Ardea, and Ixobrychus flavicollis as genus Ixobrychus, and indicated that within the subfamily Ardeinae, Nycticorax nycticorax is closely related to genus Egretta and that Ardeola bacchus and Butorides striatus are closely related to the genus Ardea.ConclusionsThe duplicate tRNAThr–CR gene order is found in most ardeid lineages, suggesting this gene order is the ancestral pattern within these birds and persisted in most lineages via concerted evolution. In two independent lineages, when the concerted evolution stopped in some subsections due to the accumulation of numerous substitutions and deletions, the duplicate tRNAThr–CR gene order was transformed into three other gene orders. The phylogenetic trees produced from concatenated rRNA and protein coding genes have high support values in most nodes, indicating that the mitochondrial genome sequences are promising markers for resolving the phylogenetic issues of ardeid birds when more taxa are added.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-573) contains supplementary material, which is available to authorized users.

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

BackgroundVertebrate 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.ResultsWe 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: tRNAArg, 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 tRNAArg gene appears to have been converted to a tRNAGln 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 tRNALeu genes, as well as a translocation of the tRNAAla 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 tRNAGlu gene from its original position. Uroplatus fimbriatus appears to retain a translocated tRNAGlu gene adjacent to the 5’ end of the major noncoding region.ConclusionsThe 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 materialThe online version of this article (doi:10.1186/1471-2164-15-930) contains supplementary material, which is available to authorized users.

The complete mitochondrial genome of the relict frog Leiopelma archeyi: Insights into the root of the frog Tree of Life

Determining the root of the anuran Tree of Life is still a contentious and open question in frog systematics. Two genera with disjunct distributions have been traditionally considered the most basal among extant frogs: Leiopelma, which is endemic to New Zealand, and Ascaphus, which lives in North America. However, their specific phylogenetic position is rather elusive because each genus shows many autapomorphies, and together they retain many symplesiomorphic characters. Therefore, several alternative hypotheses have been proposed regarding the relative phylogenetic position of both Leiopelma and Ascaphus. In order to distinguish among these competing phylogenetic hypotheses, we sequenced the complete mitochondrial (mt) genome of Leiopelma archeyi and used it along with previously reported frog mt genomes (including that of Ascaphus truei) to infer a robust phylogeny of major anuran lineages. The reconstructed maximum likelihood and Bayesian inference phylogenies recovered identical topology, which supports the sister group relationship of Ascaphus and Leiopelma, and the placement of this clade at the base of the anuran tree. Interestingly, the mt genome of L. archeyi displays a novel gene arrangement in frog mt genomes affecting the relative position of cytochrome b, trnT, NADH dehydrogenase subunit 6, trnE, and trnP genes. The tandem duplication—random loss model of gene order change explains the origin of this novel frog mt genome arrangement, which is convergent with others reported in some fishes and salamanders. These results, together with comparative data for other available vertebrate mt genomes, provide evidence that the 5′ end of the control region is a hot spot for gene order rearrangement.

Evolution of Caenorhabditis Mitochondrial Genome Pseudogenes and Caenorhabditis briggsae Natural Isolates

Although most metazoan mitochondrial genomes are highly streamlined and encode little noncoding DNA outside of the "AT" region, the accumulation of mitochondrial pseudogenes and other types of noncoding DNA has been observed in a growing number of animal groups. The nematode species Caenorhabditis briggsae harbors two mitochondrial DNA (mtDNA) pseudogenes, named Psinad5-1 and Psinad5-2, presumably derived from the nad5 protein-coding gene. Here, we provide an in-depth analysis of mtDNA pseudogene evolution in C. briggsae natural isolates and related Caenorhabditis species. Mapping the observed presence and absence of the pseudogenes onto phylogenies suggests that Psinad5-1 originated in the ancestor to C. briggsae and its recently discovered outcrossing relative species Caenorhabditis sp. 5 and Caenorhabditis sp. 9. However, Psinad5-1 was not detected in Caenorhabditis sp. 9 natural isolates, suggesting a lineage-specific loss of this pseudogene in this species. Our results corroborated the previous finding that Psinad5-2 originated within C. briggsae. The observed pattern of mitochondrial pseudogene gain and loss in Caenorhabditis was inconsistent with predictions of the tandem duplication-random loss model of mitochondrial genome evolution and suggests that intralineage recombination-like mechanisms might play a major role in Caenorhabditis mtDNA evolution. Natural variation was analyzed at the pseudogenes and flanking mtDNA sequences in 141 geographically diverse C. briggsae natural isolates. Although phylogenetic analysis placed the majority of isolates into the three previously established major intraspecific clades of C. briggsae, two new and unexpected haplotypes fell outside of these conventional groupings. Psinad5-2 copy number variation was observed among C. briggsae isolates collected from the same geographic site. Patterns of nucleotide diversity were analyzed in Psinad5-1 and Psinad5-2, and confidence intervals were found to overlap values from synonymous sites in protein-coding genes, consistent with neutral expectations. Our findings provide new insights into the mode and tempo of mitochondrial genome and pseudogene evolution both within and between Caenorhabditis nematode species.

A variant of the tandem duplication — random loss model of genome rearrangement

In [K. Chaudhuri, K. Chen, R. Mihaescu, S. Rao, On the tandem duplication-random loss model of genome rearrangement, in: SODA, 2006, pp. 564–570], Chaudhuri, Chen, Mihaescu and Rao study algorithmic properties of the tandem duplication — random loss model of genome rearrangement, well-known in evolutionary biology. In their model, the cost of one step of duplication-loss of width k is α k for α = 1 or α ≥ 2 . In this paper, we study a variant of this model, where the cost of one step of width k is 1 if k ≤ K and ∞ if k > K , for any value of the parameter K ∈ N ∪ { ∞ } . We first show that permutations obtained after p steps of width K define classes of pattern-avoiding permutations. We also compute the numbers of duplication-loss steps of width K necessary and sufficient to obtain any permutation of S n , in the worst case and on average. In this second part, we may also consider the case K = K ( n ) , a function of the size n of the permutation on which the duplication-loss operations are performed.