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.     

The novel mitochondrial gene arrangement of the cattle tick, Boophilus microplus: fivefold tandem repetition of a coding region.

We sequenced across all of the gene boundaries in the mitochondrial genome of the cattle tick, Boophilus microplus, to determine the arrangement of its genes. The mtDNA of B. microplus has a coding region, composed of tRNA(Glu) and 60 bp of the 3' end of ND1, that is repeated five times. Boophilus microplus is the first coelomate animal known to have more than two copies of a coding sequence. The mitochondrial genome of B. microplus has other unusual features, including (1) reduced T arms in tRNAs, (2) an AT bias in codon use, (3) two control regions that have evolved in concert, (4) three gene rearrangements, and (5) a stem-loop between tRNA(Gln) and tRNA(Phe). The short T arms and small control regions (CRs) of B. microplus and other ticks suggest strong selection for small genomes. Imprecise termination of replication beyond its origin, which can account for the evolution of tandem repeats of coding regions in other mitochondrial genomes, cannot explain the evolution of the fivefold repeated sequence in the mitochondrial genome of B. microplus. Instead, slipped-strand mispairing or recombination are the most plausible explanations for the evolution of these tandem repeats.     

Microarray analysis of defined Mycobacterium tuberculosis populations using RNA amplification strategies

Background

Halibuts are commercially important flatfish species confined to the North Pacific and North Atlantic Oceans. We have determined the complete mitochondrial genome sequences of four specimens each of Atlantic halibut (Hippoglossus hippoglossus), Pacific halibut (Hippoglossus stenolepis) and Greenland halibut (Reinhardtius hippoglossoides), and assessed the nucleotide variability within and between species.

Results

About 100 variable positions were identified within the four specimens in each halibut species, with the control regions as the most variable parts of the genomes (10 times that of the mitochondrial ribosomal DNA). Due to tandem repeat arrays, the control regions have unusually large sizes compared to most vertebrate mtDNAs. The arrays are highly heteroplasmic in size and consist mainly of different variants of a 61-bp motif. Halibut mitochondrial genomes lacking arrays were also detected.

Conclusion

The complexity, distribution, and biological role of the heteroplasmic tandem repeat arrays in halibut mitochondrial control regions are discussed. We conclude that the most plausible explanation for array maintenance includes both the slipped-strand mispairing and DNA recombination mechanisms.     

Complex repeat structures and novel features in the mitochondrial genomes of the diatoms Phaeodactylum tricornutum and Thalassiosira pseudonana

The mitochondrial genome of the raphid pennate diatom Phaeodactylum tricornutum has several novel features compared with the mitochondrial genomes of the centric diatom Thalassiosira pseudonana and the araphid pennate diatom Synedra acus. It is almost double the size (77,356 bp) due to a 35,454 bp sequence block consisting of an elaborate combination of direct repeats, making it the largest stramenopile (heterokont) mitochondrial genome known. In addition, the cox1 gene has a +1 translational frameshift involving Pro codons CCC and CCT, the first translational frameshift to be detected in an algal mitochondrial genome. The nad9 and rps14 genes are fused by the insertion of an in-frame sequence and cotranscribed. The nad11 gene is split into two parts corresponding to the FeS and molybdate-binding domains, but both parts are still on the mitochondrial genome, in contrast to the brown algae where the second domain appears to have been transferred to the nucleus. In contrast to P. tricornutum, the repeat region of T. pseudonana consists of a much smaller 4790 bp string of almost identical double-hairpin elements, evidence of slipped-strand mispairing and active gene conversion. The diatom mitochondrial genomes have undergone considerable gene rearrangement since the three lineages of diatoms diverged, but all three have kept their repeat regions segregated from their relatively compact coding regions.     

Complete mitochondrial genome of the seven-spotted lady beetle, Coccinella septempunctata (Coleoptera: Coccinellidae)

The complete mitochondrial genome (mitogenome) of the seven-spotted lady beetle, Coccinella septempunctata (Coleoptera: Coccinellidae), which is one of the best known insects capable of predation, is described with an emphasis on the noteworthy composition of the A+T-rich region. The C. septempunctata genome consists of 2 rRNAs, 22 tRNAs, 13 protein-coding genes, and 1 control region, designated as the A+T-rich region in insects. Along with an unusually long A+T-rich region (4469 bp), the 18,965-bp long C. septempunctata mitogenome was the largest in Coleoptera. The A+T-rich region is composed of a 2214-bp long non-repeat region composed of 80.17% A/T nucleotides and a 2256-bp long repeat region composed of 65.71% A/T nucleotides. The repeat region harbors 32 identical 70-bp long tandem repeats plus one 15-bp long incomplete first repeat. These repeat sequences may possibly have been caused by slipped-strand mispairing and unequal crossing-over events during DNA replication.     

Mitochondrial DNA sequences of five squamates: phylogenetic affiliation of snakes.

Complete or nearly complete mitochondrial DNA sequences were determined from four lizards (Western fence lizard, Warren's spinytail lizard, Terrestrial arboreal alligator lizard, and Chinese crocodile lizard) and a snake (Texas blind snake). These genomes had a typical gene organization found in those of most mammals and fishes, except for a translocation of the glutamine tRNA gene in the blind snake and a tandem duplication of the threonine and proline tRNA genes in the spinytail lizard. Although previous work showed the existence of duplicate control regions in mitochondrial DNAs of several snakes, the blind snake did not have this characteristic. Phylogenetic analyses based on different tree-building methods consistently supported that the blind snake and a colubrid snake (akamata) make a sister clade relative to all the lizard taxa from six different families. An alternative hypothesis that snakes evolved from a lineage of varanoids was not favored and nearly statistically rejected by the Kishino-Hasegawa test. It is therefore likely that the apparent similarity of the tongue structure between snakes and varanoids independently evolved and that the duplication of the control region occurred on a snake lineage after divergence of the blind snake.     

Molecular instability in the COII-tRNA(Lys) intergenic region of the human mitochondrial genome: multiple origins of the 9-bp deletion and heteroplasmy for expanded repeats.

We have identified two individuals from Glasgow in Scotland who have a deletion of one of two copies of the intergenic 9-bp sequence motif CCCCCTCTA, located between the cytochrome oxidase II (COII) and lysine tRNA (tRNA(Lys)) genes of the human mitochondrial genome. Although this polymorphism is common in Africa and Asia, it has not been reported in Northern Europe. Analysis of the mitochondrial DNA control region sequences of these two individuals suggests that they belong to a lineage that originated independently of the previously characterized African and Asian 9-bp deleted lineages. Among the Scottish population we have also identified a maternal lineage of three generations exhibiting heteroplasmy for two, three and four copies of the CCCCCTCTA motif. Polymerase chain reaction amplification across the COII-tRNA(Lys) intergenic region of these individuals gives different ratios of the three product lengths that are dependent on the concentration of the DNA-binding dye crystal violet. To investigate whether changes in repeat number were generated de novo, we constructed clones containing known numbers of the CCCCCTCTA motif. In the presence of high concentrations of crystal violet we obtained two, three and four copies of this motif when the amplification template contained only four copies. Various DNA-binding drugs are known to stabilize bulged structures in DNA and contribute to the process of slipped-strand mispairing during DNA replication. These results suggest that the COII-tRNA(Lys) intergenic region is unstable owing to slipped-strand mispairing. Although sequences containing four copies of the CCCCCTCTA motif are less stable in vitro, we observed an increase in the proportion of mitochondrial genomes with four repeats between-a mother and a daughter in the heteroplasmic lineage. From this we conclude that drift in the germ-line lineage is a main factor in the maintenance or loss of heteroplasmy.     

鸮形目4种鸟类线粒体调控区全序列的测定与比较研究

利用Long-PCR和Primer Walking的方法对鸮形目的短耳鸮、长耳鸮、纵纹腹小鸮、灰林鸮4种鸟类的线粒体调控区进行了全序列测定。结果表明：短耳鸮的调控区跃度为3290bp;长耳鸮为2848bp;纵纹腹小鸮为2444bp;灰林鸮为1771bp。短耳鸮的调控区长度是4种鸮中最大的,并且是目前已知最大的鸟类线粒体调控区。这4种鸮类调控区的基本结构和其他鸟类相似,按照碱基变化速率的不同可以分为3个区：碱基变化速率较快的外围区域Ⅰ、Ⅲ和保守的中间区域Ⅱ。这4种鸟类调控区的3’端均存在大量的串联重复序列,短耳鸮为126bp单元重复7次和78bp单元重复14次;长耳鸮为127bp单元重复8次和78bp单几重复6次;纵纹腹小鸮有3个重复单元,分别为89bp单元重复3次、77bp单元重复4次和71bp单元重复6次;灰林鸮仅有1个单元的串联重复为78bp重复5次。调控区中串联重复序列可能是由链的滑动错配产生,另外这些重复序列都能形成热力学稳定的多重茎环二级结构,而且在重复序列中还发现一些保守基序,这说明重复序列可能具有一定的生理功能,影响调控区的调重控功能从而影响线粒体基因组的复制和转录。     

The mitochondrial genome of the water spider Argyroneta aquatica (Araneae: Cybaeidae)

We sequenced nearly the entire mitochondrial genome of Argyroneta aquatica, a wholly underwater‐living spider, thereby enhancing the available genomic information for Arachnida. The confirmed sequences contained the complete set of known genes present in other metazoan mitochondrial genomes. However, the mitochondrial gene order of A. aquatica was distinctly different from that of the most distant Chelicerata Limulus polyphemus (Xiphosura), probably because of a series of gene translocations and/or inversions. Comparison of arachnid mitochondrial gene orders for the purpose of phylogenetic inference is only minimally useful, but provides a strong signal in closely related lineages. To test the basal relationships and the evolutionary pattern of tRNA gene rearrangements among Arachnida, phylogenetic analyses using amino acid sequences of the 13 protein‐coding genes were performed. An interesting feature, the five 135‐bp tandem repeats and two 363‐bp tandem repeats, was identified in the putative control region. Although control region tandem repeats have been reported in many other arachnid and metazoan species, this is the first time it has been described in spiders.     

Sequence and organization of the complete mitochondrial genomes of spotted halibut (Verasper variegatus) and barfin flounder (Verasper moseri).

In this work, the mitochondrial genomes for spotted halibut (Verasper variegatus) and barfin flounder (Verasper moseri) were completely sequenced. The entire mitochondrial genome sequences of the spotted halibut and barfin flounder were 17,273 and 17,588 bp in length, respectively. The organization of the two mitochondrial genomes was similar to those reported from other fish mitochondrial genomes containing 37 genes (2 rRNAs, 22 tRNAs and 13 protein-coding genes) and two non-coding regions (control region (CR) and WANCY region). In the CR, the termination associated sequence (ETAS), six central conserved block (CSB-A,B,C,D,E,F), three conserved sequence blocks (CSB1-3) and a region of 61-bp tandem repeat cluster at the end of CSB-3 were identified by similarity comparison with fishes and other vertebrates. The tandem repeat sequences show polymorphism among the different individuals of the two species. The complete mitochondrial genomes of spotted halibut and barfin flounder should be useful for evolutionary studies of flatfishes and other vertebrate species.     

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.     

One-way sequencing of multiple amplicons from tandem repetitive mitochondrial DNA control region

Repetitive DNA sequences not only exist abundantly in eukaryotic nuclear genomes, but also occur as tandem repeats in many animal mitochondrial DNA (mtDNA) control regions. Due to concerted evolution, these repetitive sequences are highly similar or even identical within a genome. When long repetitive regions are the targets of amplification for the purpose of sequencing, multiple amplicons may result if one primer has to be located inside the repeats. Here, we show that, without separating these amplicons by gel purification or cloning, directly sequencing the mitochondrial repeats with the primer outside repetitive region is feasible and efficient. We exemplify it by sequencing the mtDNA control region of the mosquito Aedes albopictus, which harbors typical large tandem DNA repeats. This one-way sequencing strategy is optimal for population surveys.     

A complete mitochondrial genome sequence of Ogura-type male-sterile cytoplasm and its comparative analysis with that of normal cytoplasm in radish (Raphanus sativus L.)

ABSTRACT: BACKGROUND: Plant mitochondrial genome has unique features such as large size, frequent recombination and incorporation of foreign DNA. Cytoplasmic male sterility (CMS) is caused by rearrangement of the mitochondrial genome, and a novel chimeric open reading frame (ORF) created by shuffling of endogenous sequences is often responsible for CMS. The Ogura-type male-sterile cytoplasm is one of the most extensively studied cytoplasms in Brassicaceae. Although the gene orf138 has been isolated as a determinant of Ogura-type CMS, no homologous sequence to orf138 has been found in public databases. Therefore, how orf138 sequence was created is a mystery. In this study, we determined the complete nucleotide sequence of two radish mitochondrial genomes, namely, Ogura- and normal-type genomes, and analyzed them to reveal the origin of the gene orf138. RESULTS: Ogura- and normal-type mitochondrial genomes were assembled to 258,426-bp and 244,036-bp circular sequences, respectively. Normal-type mitochondrial genome contained 33 protein-coding and three rRNA genes, which are well conserved with the reported mitochondrial genome of rapeseed. Ogura-type genomes contained same genes and additional atp9. As for tRNA, normal-type contained 17 tRNAs, while Ogura type contained 17 tRNAs and one additional trnfM. The gene orf138 was specific to Ogura-type mitochondrial genome, and no sequence homologous to it was found in normal-type genome. Comparative analysis of the two genomes revealed that radish mitochondrial genome consists of 11 syntenic regions (length >3kb, similarity >99.9%). It was shown that short repeats and overlapped repeats present in the edge of syntenic regions were involved in recombination events during evolution to interconvert two types of mitochondrial genome. Ogura-type mitochondrial genome has four unique regions (2,803 bp, 1,601 bp, 451 bp and 15,255 bp in size) that are non-syntenic to normal-type genome, and the gene orf138 was found to be located at the edge of the largest unique region. Blast analysis performed to assign the unique regions showed that about 80% of the region was covered by short homologous sequences to the mitochondrial sequences of normal-type radish or other reported Brassicaceae species, although no homology was found for the remaining 20% of sequences. CONCLUSIONS: Ogura-type mitochondrial genome was highly rearranged compared with the normal-type genome by recombination through one large repeat and multiple short repeats. The rearrangement has produced four unique regions in Ogura-type mitochondrial genome, and most of the unique regions are composed of known Brassicaceae mitochondrial sequences. This suggests that the regions unique to the Ogura-type genome were generated by integration and shuffling of pre-existing mitochondrial sequences during the evolution of Brassicaceae, and novel genes such as orf138 could have been created by the shuffling process of mitochondrial genome.     

DNA tertiary structures formed in vitro by misaligned hybridization of multiple tandem repeat sequences.

DNA tertiary structures are shown to be formed by denaturation and reannealing in vitro of molecularly-cloned DNA containing multiple tandem repeat sequences. Electron microscopy of homoduplex DNA molecules containing the human c-Harvey-ras gene revealed knot-like structures which mapped to the position of the 812 bp variable tandem repeat (VTR) sequence. We propose that the structures result from slipped-strand mispairing within the VTR and hybridisation of homologous repetitive sequences in the single-stranded loops so produced. Similar structures were also found in freshly-linearized supercoiled plasmids. More complex knot-like structures were found in homoduplexes of a 4 kb tandem array from the hypervariable region 3' to the human alpha-globin locus. Formation of such DNA tertiary structures in vitro also provides a practical method for identifying and mapping direct tandem repeat arrays that are at least 800 bp long.     

Slipped-strand mispairing: a major mechanism for DNA sequence evolution

Simple repetitive DNA sequences are a widespread and abundant feature of genomic DNA. The following several features characterize such sequences: (1) they typically consist of a variety of repeated motifs of 1-10 bases--but may include much larger repeats as well; (2) larger repeat units often include shorter ones within them; (3) long polypyrimidine and poly-CA tracts are often found; and (4) tandem arrangements of closely related motifs are often found. We propose that slipped-strand mispairing events, in concert with unequal crossing- over, can readily account for all of these features. The frequent occurrence of long tandem repeats of particular motifs (polypyrimidine and poly-CA tracts) appears to result from nonrandom patterns of nucleotide substitution. We argue that the intrahelical process of slipped-strand mispairing is much more likely to be the major factor in the initial expansion of short repeated motifs and that, after initial expansion, simple tandem repeats may be predisposed to further expansion by unequal crossing-over or other interhelical events because of their propensity to mispair. Evidence is presented that single-base repeats (the shortest possible motifs) are represented by longer runs in mammalian introns than would be expected on a random basis, supporting the idea that SSM may be a ubiquitous force in the evolution of the eukaryotic genome. Simple repetitive sequences may therefore represent a natural ground state of DNA unselected for coding functions.      

The complete mitochondrial genomes sequences of Asio flammeus and Asio otus and comparative analysis

Mitochondrial DNA (mtDNA), as a model sys-tem, has been extensively used for molecular phy-logenetic and evolutionary analysis[1]. With the ad-vances in DNA sequencing technology, more andmore researchers prefer to use complete mitochondrialgenome for phylogenetic analysis[2—4]. Since the com-plete sequencing of human mtDNA in 1981 (Andersonet al., 1981)[5], 342 vertebrate mitochondrial genomeshave so far been sequenced. Up to now the completesequences of 29 avian mitochondrial genomes h…