Mitochondrial genome of the Komodo dragon: efficient sequencing method with reptile-oriented primers and novel gene rearrangements.

The mitochondrial genome of the Komodo dragon (Varanus komodoensis) was nearly completely sequenced, except for two highly repetitive noncoding regions. An efficient sequencing method for squamate mitochondrial genomes was established by combining the long polymerase chain reaction (PCR) technology and a set of reptile-oriented primers designed for nested PCR amplifications. It was found that the mitochondrial genome had novel gene arrangements in which genes from NADH dehydrogenase subunit 6 to proline tRNA were extensively shuffled with duplicate control regions. These control regions had 99% sequence similarity over 700 bp. Although snake mitochondrial genomes are also known to possess duplicate control regions with nearly identical sequences, the location of the second control region suggested independent occurrence of the duplication on lineages leading to snakes and the Komodo dragon. Another feature of the mitochondrial genome of the Komodo dragon was the considerable number of tandem repeats, including sequences with a strong secondary structure, as a possible site for the slipped-strand mispairing in replication. These observations are consistent with hypotheses that tandem duplications via the slipped-strand mispairing may induce mitochondrial gene rearrangements and may serve to maintain similar copies of the control region.     

Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews

The complete mitochondrial DNA (mtDNA) control region was amplified and directly sequenced in two species of shrew, Crocidura russula and Sorex araneus (Insectivora, Mammalia). The general organization is similar to that found in other mammals: a central conserved region surrounded by two more variable domains. However, we have found in shrews the simultaneous presence of arrays of tandem repeats in potential locations where repeats tend to occur separately in other mammalian species. These locations correspond to regions which are associated with a possible interruption of the replication processes, either at the end of the three-stranded D-loop structure or toward the end of the heavy-strand replication. In the left domain the repeated sequences (R1 repeats) are 78 bp long, whereas in the right domain the repeats are 12 bp long in C. russula and 14 bp long in S. araneus (R2 repeats). Variation in the copy number of these repeated sequences results in mtDNA control region length differences. Southern blot analysis indicates that level of heteroplasmy (more than one mtDNA form within an individual) differs between species. A comparative study of the R2 repeats in 12 additional species representing three shrew subfamilies provides useful indications for the understanding of the origin and the evolution of these homologous tandemly repeated sequences. An asymmetry in the distribution of variants within the arrays, as well as the constant occurrence of shorter repeated sequences flanking only one side of the R2 arrays, could be related to asymmetry in the replication of each strand of the mtDNA molecule. The pattern of sequence and length variation within and between species, together with the capability of the arrays to form stable secondary structures, suggests that the dominant mechanism involved in the evolution of these arrays in unidirectional replication slippage.      

Conserved PCR primer set designing for closely-related species to complete mitochondrial genome sequencing using a sliding window-based PSO algorithm

Background

Complete mitochondrial (mt) genome sequencing is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. For long template sequencing, i.e., like the entire mtDNA, it is essential to design primers for Polymerase Chain Reaction (PCR) amplicons which are partly overlapping each other. The presented chromosome walking strategy provides the overlapping design to solve the problem for unreliable sequencing data at the 5′ end and provides the effective sequencing. However, current algorithms and tools are mostly focused on the primer design for a local region in the genomic sequence. Accordingly, it is still challenging to provide the primer sets for the entire mtDNA.

Methodology/Principal Findings

The purpose of this study is to develop an integrated primer design algorithm for entire mt genome in general, and for the common primer sets for closely-related species in particular. We introduce ClustalW to generate the multiple sequence alignment needed to find the conserved sequences in closely-related species. These conserved sequences are suitable for designing the common primers for the entire mtDNA. Using a heuristic algorithm particle swarm optimization (PSO), all the designed primers were computationally validated to fit the common primer design constraints, such as the melting temperature, primer length and GC content, PCR product length, secondary structure, specificity, and terminal limitation. The overlap requirement for PCR amplicons in the entire mtDNA is satisfied by defining the overlapping region with the sliding window technology. Finally, primer sets were designed within the overlapping region. The primer sets for the entire mtDNA sequences were successfully demonstrated in the example of two closely-related fish species. The pseudo code for the primer design algorithm is provided.

Conclusions/Significance

In conclusion, it can be said that our proposed sliding window-based PSO algorithm provides the necessary primer sets for the entire mt genome amplification and sequencing.     

Sequence organization of repetitive sequences enriched in small polydisperse circular DNAs from HeLa cells

A total of 36 clones were randomly selected from a recombinant DNA library of small polydisperse circular DNA (spcDNA) molecules from HeLa cells and were shown to contain repetitive sequences of different reiteration frequencies that ranged from several hundred to several hundred thousand per genome. Sequencing of representative clones revealed tandem repeats of alphoid (alpha) satellite DNA, clustered repeats of the Alu family, KpnI family sequences, tandem repeats of an alpha satellite DNA specific to the X chromosome (alpha X), and A + T-rich segments carrying short stretches of poly(A) or poly(T). DNA rearrangement was frequently found in the repetitive sequences enriched in these spcDNA clones. Short regions of homology that were patchy and inverted were often found, especially at the novel joint where spcDNA sequences are circularized. The presence of these inverted repeats suggests that HeLa spcDNAs are formed by a mechanism that involves looping out of the spcDNA region and joining of the flanking DNA by illegitimate recombination.     

中华鲟(Acipenser sinensis)及相关种类的mtDNA控制区串联重复序列及其进化意义

采用 PCR技术和 DNA测序技术 ,发现了我国一级珍稀保护动物中华鲟 ( Acipensersinensis)线粒体 DNA( mt DNA)的控制区 ( D- loop)存在数目不等的串联重复序列 ,该重复序列造成了中华鲟广泛的异质性现象 .从分子水平进行了不同类型重复序列变化规律的研究 ,同时还初探了重复序列在我国其它几种鲟鱼类的存在情况 ,发现在白鲟 ( Psephurus gladius)、达氏鲟 ( A.dabryanus)和史氏鲟 ( A.schrenckii)均存在类似的重复序列结构 .序列比较分析表明 ,不同鲟鱼类重复序列在鲟鱼类进化过程中扮演着一定的角色 ,很有可能碱基差异大小与它们的亲缘关系的远近呈正相关 .     

Hypervariability of the D-loop region in mitochondrial DNA of Russian sturgeon <Emphasis Type="Italic">Acipenser gueldenstaedtii</Emphasis> (Acipenseriformes,Acipenseridae)

The structure of the D-loop region in mitochondrial DNA (mtDNA) of Russian sturgeon Acipenser gueldenstaedtii from the Azov Sea population was studied with the method of direct sequencing. Interindividual heteroplasmy of the length of mtDNA in the region of D-loop realized by the presence of a different number of tandem repeats (82 pairs of bases) was found. Analysis of tandem repeats in the D-loop region in mtDNA in the studied sample (28 individuals) revealed eight mitotypes differed in the pattern of nucleotide substitution and in the number of tandem repeats (2, 3, and 4 repeats). Revealed mitotypes can be considered as potential genetic markers for different biological groups, schools, or seasonal races of A. gueldenstaedtii.     

Cytological characterization of the tandem repetitive sequences and their methylation status in the Antirrhinum majus genome

Tandem repetitive sequences are DNA motifs common in the genomes of eukaryotic species and are often embedded in heterochromatic regions. In most eukaryotes, ribosomal genes, as well as centromeres and telomeres or subtelomeres, are associated with abundant tandem arrays of repetitive sequences and typically represent the final barriers to completion of whole-genome sequencing. The nature of these repeats makes it difficult to estimate their actual sizes. In this study, combining the two cytological techniques DNA fiber-FISH and pachytene chromosome FISH allowed us to characterize the tandem repeats distributed genome wide in Antirrhinum majus and identify four types of tandem repeats, 45S rDNA, 5S rDNA, CentA1, and CentA2, representing the major tandem repetitive components, which were estimated to have a total length of 18.50 Mb and account for 3.59% of the A. majus genome. FISH examination revealed that all the tandem repeats correspond to heterochromatic knobs along the pachytene chromosomes. Moreover, the methylation status of the tandem repeats was investigated in both somatic cells and pollen mother cells from anther tissues using an antibody against 5-methylcytosine combined with sequential FISH analyses. Our results showed that these repeats were hypomethylated in anther tissues, especially in the pollen mother cells at pachytene stage.     

Structure and evolution of the mitochondrial control region of the pollen beetle Meligethes thalassophilus (Coleoptera: Nitidulidae).

The organization of the mitochondrial DNA (mtDNA) control region (CR) of the pollen beetle Meligethes thalassophilus is described. This mtDNA CR represents the longest sequenced for beetles so far, since the entire nucleotide sequence ranges from approximately 5000 to approximately 5500 bp. The CR of M. thalassophilus is organized in three distinct domains: a conserved domain near the tRNAIle gene, a variable domain flanking the 12S rRNA gene, and a relatively large central tandem array made up of a variable number of approximately 170 bp repeats that is responsible for the intraspecific length variation observed. Like other CRs found in insects, the M. thalassophilus CR contains two long homopolymeric runs that may be involved in mtDNA replication. Furthermore, conserved stem-and-loop structures in the repetitive domain were identified and their possible role in generating length variation is examined. Intraspecific comparison of the tandem repeat elements of M. thalassophilus suggests mechanisms of concerted evolution leading to homogenization of the repetitive region. The utility of such an array of tandem repeats as a genetic marker for assessing population-level variability and evolutionary relationships among populations is discussed. Finally, the technical difficulties found in isolating the mtDNA CR in beetles are remarked upon.     

Repetitive sequences in Eurasian lynx (Lynx lynx L.) mitochondrial DNA control region

Mitochondrial DNA (mtDNA) control region (CR) of numerous species is known to include up to five different repetitive sequences (RS1-RS5) that are found at various locations, involving motifs of different length and extensive length heteroplasmy. Two repetitive sequences (RS2 and RS3) on opposite sides of mtDNA central conserved region have been described in domestic cat (Felis catus) and some other felid species. However, the presence of repetitive sequence RS3 has not been detected in Eurasian lynx (Lynx lynx) yet. We analyzed mtDNA CR of 35 Eurasian lynx (L. lynx L.) samples to characterize repetitive sequences and to compare them with those found in other felid species. We confirmed the presence of 80 base pairs (bp) repetitive sequence (RS2) at the 5' end of the Eurasian lynx mtDNA CR L strand and for the first time we described RS3 repetitive sequence at its 3' end, consisting of an array of tandem repeats five to ten bp long. We found that felid species share similar RS3 repetitive pattern and fundamental repeat motif TACAC.     

Detection of species-specific long VNTRs in mitochondrial control region and their application to identifying sympatric Hong Kong grouper (Epinephelus akaara) and yellow grouper (Epinephelus awoara)

Length variation and heteroplasmy were observed in PCR products of the first half of mtDNA control region of both Hong Kong grouper (Epinephelus akaara) and yellow grouper (Epinephelus awoara). DNA sequencing unveiled the phenomena were caused by the presence of species-specific long variable number tandem repeats (VNTRs). This is the first report on the mtDNA VNTRs and their heteroplasmy in groupers. Moreover, these VNTRs are also the longest such structure found in teleost fish. Thereafter, we designed two species-specific PCR reverse primers according to the 3' end sequences of the VNTRs and successfully established assays for the identification of these two sympatric grouper species.     

鸮形目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次。调控区中串联重复序列可能是由链的滑动错配产生,另外这些重复序列都能形成热力学稳定的多重茎环二级结构,而且在重复序列中还发现一些保守基序,这说明重复序列可能具有一定的生理功能,影响调控区的调重控功能从而影响线粒体基因组的复制和转录。     

Ancient DNA analysis of human neolithic remains found in northeastern Siberia

We successfully extracted DNA from a bone sample of a Neolithic skeleton (dated 3,600 +/- 60 years BP) excavated in northeastern Yakutia (east Siberia). Ancient DNA was analyzed by autosomal STRs (short tandem repeats) and by sequencing of the hypervariable region I (HV1) of the mitochondrial DNA (mtDNA) control region. The STR profile, the mitochondrial haplotype, and the haplogroup determined were compared with those of modern Eurasian and Native American populations. The results showed the affinity of this ancient skeleton with both east Siberian/Asian and Native American populations.     

A topoisomerase II cleavage site is associated with a novel mitochondrial DNA deletion

Mitochondrial myopathies and encephalopathies can be caused by nucleotide substitutions, deletions or duplications of the mitochondrial DNA (mtDNA). In one such disorder, Kearns-Sayre Syndrome (KSS), large-scale hetero-plasmic mtDNA deletions are often found. We describe a 14-year-old boy with clinical features of KSS, plus some additional features. Analysis of the entire mitochondrial genome by the polymerase chain reaction and Southern blotting revealed a 7864-bp mtDNA deletion, heteroplasmic in its tissue distribution. DNA sequencing established that the deletion was between nucleotides 6238 and 14103, and flanked by a 4-bp (TCCT) direct repeat sequence. Deletions between direct repeats have been hypothesised to occur by a slipped-mismatching or illegitimate recombination event, or following the DNA cleavage action of topoisomerase II. Analysis of the gene sequence in the region surrounding the mtDNA deletion breakpoint in this patient revealed the presence of putative vertebrate topoisomerase II sites. We suggest that direct repeat sequences, together with putative topoisomerase II sites, may predispose certain regions of the mitochondrial genome to deletions.     

Unique features in the mitochondrial D-loop region of the European seabass Dicentrarchus labrax

We have cloned and sequenced the displacement-loop (D-loop) region of the mitochondrial DNA (mtDNA) from the European seabass Dicentrarchus labrax (Dl). This sequencing revealed the presence of four tandemly repeated elements (R1, R2, R3 and R4); the individual variation in mtDNA total length is entirely accounted for by their variable number. The individuals examined also possessed an imperfect copy of one of the tandem repeats (ΨR2). At least one termination-associated sequence (TAS) is present in each of the repeats and in two copies 5′ upstream from the tandem array as well. The alignment of the Dl D-loop region with D-loop sequences from four other Teleosts and one Chondrosteus showed the Dl sequence to be larger than that of other fish. The extraordinary length of the D1 D-loop sequence is also due to the 5′ and 3′ regions that are flanking the tandem array, the largest ones to date analyzed in fish. In this study, we also report the unique organization and localization of putative TAS and conserved-sequence block (CSB) elements, and the presence of a conserved 218-bp sequence in the D1 D-loop region.     

Small, repetitive DNAs contribute significantly to the expanded mitochondrial genome of cucumber

Closely related cucurbit species possess eightfold differences in the sizes of their mitochondrial genomes. We cloned mitochondrial DNA (mtDNA) fragments showing strong hybridization signals to cucumber mtDNA and little or no signal to watermelon mtDNA. The cucumber mtDNA clones carried short (30-53 bp), repetitive DNA motifs that were often degenerate, overlapping, and showed no homology to any sequences currently in the databases. On the basis of dot-blot hybridizations, seven repetitive DNA motifs accounted for >13% (194 kb) of the cucumber mitochondrial genome, equaling >50% of the size of the Arabidopsis mitochondrial genome. Sequence analysis of 136 kb of cucumber mtDNA revealed only 11.2% with significant homology to previously characterized mitochondrial sequences, 2.4% to chloroplast DNA, and 15% to the seven repetitive DNA motifs. The remaining 71.4% of the sequence was unique to the cucumber mitochondrial genome. There was <4% sequence colinearity surrounding the watermelon and cucumber atp9 coding regions, and the much smaller watermelon mitochondrial genome possessed no significant amounts of cucumber repetitive DNAs. Our results demonstrate that the expanded cucumber mitochondrial genome is in part due to extensive duplication of short repetitive sequences, possibly by recombination and/or replication slippage.     

Heteroplasmy and evidence for recombination in the mitochondrial control region of the flatfish Platichthys flesus

The general assumption that mitochondrial DNA (mtDNA) does not undergo recombination has been challenged recently in invertebrates. Here we present the first direct evidence for recombination in the mtDNA of a vertebrate, the flounder Platichthys flesus. The control region in the mtDNA of this flatfish is characterized by the presence of a variable number of tandem repeats and a high level of heteroplasmy. Two types of repeats were recognized, differing by two C-T point mutations. Most individuals carry a pure "C" or a pure "T" array, but one individual showed a compound "CT" array. Such a compound array is evidence for recombination in the mtDNA control region from the flounder.     

Generation of VNTRs and heteroplasmy by sequence turnover in the mitochondrial control region of two elephant seal species

We describe an unusual repetitive DNA region located in the 3′ end of the light (L)-strand in the mitochondrial control region of two elephant seal species. The array of tandem repeats shows both VNTR (variable-number tandem repeat) and sequence variation and is absent from 12 compared mammalian species, except for the occurrence in the same location of a distinct repetitive region in rabbit mtDNA and a similar repeat in the harbor seal. The sequence composition and arrangement of the repeats differ considerably between the northern elephant seal (Mirounga angustirostris) and the southern species (M. leonina) despite an estimated divergence time of 1 MY (based on an mtDNA-RNA gene and the nonrepetitive control region). Analysis of repeat sequence relationships within and between species indicate that divergence in sequence and structure of repeats has involved both slippagelike and unequal crossingover processes of turnover, generating very high levels of divergence and heteroplasmy. Presented at the NATO Advanced Research Workshop onGenome Organization and Evolution, Spetsai, Greece, 16–22 September 1992     

暗纹东方鲀线粒体基因组核苷酸全序列测定与分析

以暗纹东方鲀(Takifugu fasciatus)肝的线粒体DNA为模板,参照红鳍东方鲀(T.rubripes)等近源鱼类的线粒体基因组DNA序列,设计合成14对特异引物,进行PCR扩增并测序,首次获得了暗纹东方鲀线粒体基因组全序列。结果表明,暗纹东方鲀线粒体基因组序列全长16 444 bp(GenBank登录号为GQ409967),A+T含量为55.8%,其mtDNA结构与其他脊椎动物相似,由22个tRNA基因、2个rRNA基因、13个蛋白质编码基因和1段819 bp非编码的控制区(D-loop)所组成。蛋白质基因除COⅠ和ND6的起始密码子为GTG、CCT以外,均为典型的起始密码子ATG。ND1、ATPase8、COⅢ、ND4L、ND5、Cyt b使用典型的终止密码子TAA,其他的使用不完全终止密码子。除ND6和tRNAGln、tRNAAla、tRNAAsn、tRNACys、tRNATyr、tRNASer、tRNAGlu、tRNAPro在L-链上编码之外,其余基因均在H-链编码。基因排列顺序与已测定的鲀类一致,这显示了鲀类线粒体基因排列顺序上的保守性。tRNA基因核苷酸长度为64～73nt,预测了22个tRNA基因的二级结构,均呈较为典型的三叶草状。基于19种鲀类mtDNA全序列构建的进化树表明,暗纹东方鲀与红鳍东方鲀、中华东方鲀(T.chinensis)聚成一个姊妹群。结果还支持东方鲀属鱼类为一单系类群。     

The Owen mitochondrial genome in sugar beet (Beta vulgaris L.): possible mechanisms of extensive rearrangements and the origin of the mitotype-unique regions

The mitochondrial genomes of normal fertile and male-sterile (Owen CMS) cytoplasms of sugar beet are highly rearranged relative to each other and dozens of inversional recombinations and other reshuffling events must be postulated to interconvert the two genomes. In this paper, a comparative analysis of the entire nucleotide sequences of the two genomes revealed that most of the inversional recombinations involved short repeats present at their endpoints. Attention was also focused on the origin of the Owen CMS-unique mtDNA regions, which occupy 13.6% of the Owen genome and are absent from the normal mtDNA. BLAST search was performed to assign the sequences, and as a result, 7.6% of the unique regions showed significant homology to previously determined mitochondrial sequences, 17.9% to nuclear DNA, 4.6% to mitochondrial episomes, and 0.1% to plastid DNA. Southern blot analysis revealed that additional sequences of nuclear origin may be included within the unique regions. We also found that the copies of many short repeat families are scattered throughout the unique regions. This suggests that, in addition to the incorporation of foreign DNAs, extensive duplication of short repetitive sequences and continued scrambling of mtDNA sequences may be implicated in the generation of the Owen CMS-unique regions.