Recombination and Evolution of Duplicate Control Regions in the Mitochondrial Genome of the Asian Big-Headed Turtle,Platysternon megacephalum

Complete mitochondrial (mt) genome sequences with duplicate control regions (CRs) have been detected in various animal species. In Testudines, duplicate mtCRs have been reported in the mtDNA of the Asian big-headed turtle, Platysternon megacephalum, which has three living subspecies. However, the evolutionary pattern of these CRs remains unclear. In this study, we report the completed sequences of duplicate CRs from 20 individuals belonging to three subspecies of this turtle and discuss the micro-evolutionary analysis of the evolution of duplicate CRs. Genetic distances calculated with MEGA 4.1 using the complete duplicate CR sequences revealed that within turtle subspecies, genetic distances between orthologous copies from different individuals were 0.63% for CR1 and 1.2% for CR2app:addword:respectively, and the average distance between paralogous copies of CR1 and CR2 was 4.8%. Phylogenetic relationships were reconstructed from the CR sequences, excluding the variable number of tandem repeats (VNTRs) at the 3′ end using three methods: neighbor-joining, maximum likelihood algorithm, and Bayesian inference. These data show that any two CRs within individuals were more genetically distant from orthologous genes in different individuals within the same subspecies. This suggests independent evolution of the two mtCRs within each P. megacephalum subspecies. Reconstruction of separate phylogenetic trees using different CR components (TAS, CD, CSB, and VNTRs) suggested the role of recombination in the evolution of duplicate CRs. Consequently, recombination events were detected using RDP software with break points at ≈290 bp and ≈1,080 bp. Based on these results, we hypothesize that duplicate CRs in P. megacephalum originated from heterological ancestral recombination of mtDNA. Subsequent recombination could have resulted in homogenization during independent evolutionary events, thus maintaining the functions of duplicate CRs in the mtDNA of P. megacephalum.     

Phenotypic Effects of Isolated pRiA4 TL-DNA rol Genes in the Presence of Intact TR-DNA in Transgenic Plants of Solanum dulcamara L.

The effect of the rol genes, together with the TR-DNA of pRiA4on the phenotype of Solanum dulcamara plants, was analysed.Plants transformed by Agrobacterium strain BN1010: :rolA (rolA^7plus;TR⁺)exhibited severe leaf wrinkling, whereas plants transformedby strain BN1010: :rolC (ro/C⁺TR⁺) had a typical ‘hairyroot’ phenotype. Leaf discs excised from these latterplants produced roots on hormone-free medium. BN1010: :rolABC(rolABC⁺TR⁺) transformed plants had an exaggerated transformedphenotype. Some of the BN1010: :rolABC transformants had positivelygeotropic root growth which correlated with the presence ofmultiple copies of the TR-DNA. S. dulcamara plants, transformedby the TR-DNA region only, exhibited epinasty. Scanning electronmicroscopy of plants containing various regions of agropineRi T-DNA revealed that transformation causes changes in basicplant siructure.     

Structure,DNA sequence variation and phylogenetic implications of the mitochondrial control region in horseshoe bats

There have been few studies of the structural and evolutionary characteristics of the mitochondrial control region (CR) in rhinolophids, yet this could have important consequences for the interpretation of phylogenetic relationships within this group. Here we sequenced and analyzed the CR of 37 individuals from 12 Rhinolophus species, including 2 species from GenBank. The length of the CR ranged from 1335 to 1514 bp, and the base composition was very similar among species. The CR of horseshoe bats, like that of other mammals, could be subdivided into a central conserved domain (CD) and two flanking variable domains, extended termination associated sequences (ETAS), and conserved sequence blocks (CSB). Besides the common conserved blocks (ETAS1, ETAS2, F-B boxes, CSB1, CSB2 and CSB3) found in 3 domains, an ETAS2-like and a CSB1-like element were also detected in the ETAS and CSB domains, respectively, in all individuals. Notwithstanding a short tandem repeat (11 or 13 bp) between CSB1 and CSB2 in all specimens, the base composition, copy number and arrays are all variable. A long tandem repeat (79 bp) was only identified in the ETAS domain in one individual of R. pusillus. Phylogenetic reconstructions based on the CR sequences indicated that the molecular phylogenetic relationships among some Rhinolophus species were inconsistent with the results of phenetic analyses, but similar to phylogenetic constructions using cytochrome b. An unidentified species R. sp and 3 species from the philippinensis-group that were clearly morphologically different comprised a monophyletic group, which could have resulted from morphological independent evolution.     

Complete mitochondrial genome of the Eurasian flying squirrel Pteromys volans (Sciuromorpha, Sciuridae) and revision of rodent phylogeny

In this study, the complete mitochondrial genome of the Eurasian flying squirrel Pteromys volans (Rodentia, Sciuromorpha, Sciuridae) was sequenced and characterized in detail. The entire mitochondrial genome of P. volans consisted of 16,513 bp and contained 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and two non-coding regions. Its gene arrangement pattern was consistent with the mammalian ground pattern. The overall base composition and AT contents were similar to those of other rodent mitochondrial genomes. The light-strand origin generally identified between tRNA ^Asn and tRNA ^Cys consisted of a secondary structure with an 11-bp stem and an 11-bp loop. The large control region was constructed of three characteristic domains, ETAS, CD, and CSB without any repeat sequences. Each domain contained ETAS1, subsequences A, B, and C, and CSB1, respectively. In order to examine phylogenetic contentious issues of the monophyly of rodents and phylogenetic relationships among five rodent suborders, here, phylogenetic analyses based on nucleotide sequence data of the 35 rodent and 3 lagomorph mitochondrial genomes were performed using the Bayesian inference and maximum likelihood method. The result strongly supported the rodent monophyly with high node confidence values (BP 100 % in ML and BPP 1.00 in BI) and also monophylies of four rodent suborders (BP 85–100 % in ML and BPP 1.00 in BI), except for Anomalumorpha in which only one species was examined here. Also, phylogenetic relationships among the five rodent suborders were suggested and discussed in detail.     

Genes Encoding the Group I Intron-containing tRNALeu and Subunit L of NADH Dehydrogenase from the Cyanobacterium Synechococcus PCC 6301

A part of the tRNA^Leu (UAA) gene containing a 240-nucleotidegroup I intron was amplified by PCR from cyanobacterium SynechococcusPCC 6301 genomic DNA. The pre-tRNA synthesized from the clonedPCR product was efficiently self-spliced in vitro under physiologicalconditions. The gene encoding the tRNA^Leu (UAA), trnL-UAA, wasisolated from a Synechococcus PCC 6301 genomic library and thenucleotide sequence of a 2,167-bp portion was determined. ThetrnL-UAA consists of a 34-bp 5' exon, a 240-bp group I intronand a 50-bp 3' exon. In addition, three open reading frames(ORF1, ORF2 and ORF3) were found in the 5' and 3' flanking regionsof trnL-UAA. The predicted protein sequence of ORF3, which islocated 74-bp upstream from trnL-UAA on the opposite strand,shows 66.2% amino acid identity to that of the SynechocystisPCC 6803 gene encoding subunit L of NADH dehydrogenase (ndhL).     

Evidence for multiple functional copies of the male sex-determining locus, sry, in African murine rodents

Southern hybridization data suggest that the male sex-determining locus, Sry, is often duplicated in rodents. Here we explore DNA sequence evolution of orthologous and paralogous copies of Sry isolated from six species of African murines. PCR amplification followed by direct sequencing revealed from two to four copies of Sry per species. All copies include a long open reading frame, with a stop codon that coincides closely with the stop codon of the house mouse, Mus musculus, a species known to have a single copy of Sry. A phylogenetic analysis suggests that there are at least seven paralogous copies of Sry in this group of rodents. Putative orthologues are identical; sequence divergence among putative paralogues ranges from 1 to 8% (excluding the CAG repeat), with much lower levels of divergence in the high-mobility group (HMG-box) region than in the C-terminal region. A high proportion of nucleotide substitutions in both regions result in amino-acid replacement. The long open reading frame, conserved HMG-box, and pattern of evolution of the putative paralogues suggest that they are functional. Received: 4 October 1996 / Accepted: 17 January 1997     

Sequence evolution in and around the mitochondrial control region in birds

By cloning and sequencing 3.4 kilobases of snow goose mtDNA we found that the ND5 gene is followed by the genes for cytochrome b, tRNA^Thr, tRNA^Pro, ND6, tRNA^Glu, the control region, tRNA^Phe, and srRNA. This order is identical to that of chicken, quail, and duck mtDNA but differs from that of mammals and a frog (Xenopus). The mean extent of difference due to base substitution between goose and chicken is generally closer to the same comparison between rat and mouse but less than that between human and cow. For one of the nine regions compared (tRNA^Glu), the bird differences appear to be anomalous, possibly implicating altered functional constraints. Within the control region, several short sequences common to mammals are also conserved in the birds. Comparison of the goose control region with that of quail and chicken suggests that a sequence element with similarity to CSB-1 duplicated once prior to the divergence of goose and chicken and again on the lineage leading to chicken. Between goose (or duck) and chicken there are four times more transversions at the third positions of fourfold-degenerate codons in mitochondrial than in nuclear genes.Abbreviations CSB conserved sequence block - cytb cytochrome b - ND NADH dehydrogenase - srRNA small-subunit ribosomal RNA Deceased July 21, 1991 Correspondence to: T.W. Quinn at the University of Denver     

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

Background

The 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.

Results

The complete mitochondrial genome sequences of the sixteen ardeid birds exhibited four distinct mitochondrial gene orders in which two of them, named as “duplicate tRNA^Glu–CR” and “duplicate tRNA^Thr–tRNA^Pro 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.

Conclusions

The 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 material

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

Complete mtDNA of Meretrix lusoria (Bivalvia: Veneridae) reveals the presence of an atp8 gene, length variation and heteroplasmy in the control region

The complete nucleotide sequence of the mitochondrial genome of the clam Meretrix lusoria (Bivalvia: Veneridae) was determined. It comprises 20,268 base pairs (bp) and contains 13 protein-coding genes, including ATPase subunit 8 (atp8), two ribosomal RNAs, 22 transfer RNAs, and a non-coding control region. The atp8 encodes a protein of 39 amino acids. All genes are encoded on the same strand. A putative control region (CR or D-loop) was identified in the major non-coding region (NCR) between the tRNA^Gly and tRNA^Gln. A 1087 bp tandem repeat fragment was identified that comprises nearly 11 copies of a 101 bp motif and accounts for approximately 41% of the NCR. The 101 bp tandem repeat motif of the NCR can be folded into a stem–loop secondary structure. Samples of eight individuals from Hainan and Fujian provinces were collected and their NCR regions were successfully amplified and sequenced. The data revealed a highly polymorphic VNTR (variable number of tandem repeats) associated with high levels of heteroplasmy in the D-loop region. The size of the CR ranged from 1942 to 3354 bp depending upon the copy number of the repeat sequence.     

Transformation and Characterization of Transgenic Plants of Solanum dulcamara L.--Incidence of Transgene Silencing

A transformation system is described for Solanum dulcamara usingthe supervirulentAgrobacterium tumefaciens strain 1065, carryingboth the ß-glucuronidase (gus) and neomycin phosphotransferaseII (npt II) genes adjacent to the right and left T-DNA borders,respectively. Leaf explants were more efficient for the productionof transformed plants compared to stem explants on medium containing50 mg l^-1of kanamycin sulphate. A 1:10 (v:v) dilution of anovernight culture ofAgrobacterium gave optimal transformationin terms of transgenic plant regeneration. From a total of 174kanamycin-resistant plants selected by their antibiotic resistance,16 failed to exhibit GUS activity. Southern analysis revealedthat these GUS-negative transformants originated from threeindependently transformed cell lines. Restriction enzyme analysesshowed that the GUS-negative plants had both the gus and nptII genes integrated into their genome (one plant had a singlecopy of each gene; the other two plants had multiple copies),with major rearrangement of the gus gene occurring in plantswith several copies of the transgene. GUS-negative plants showedleaf malformations, delayed flowering and a reduction in flower,fruit and seed production compared to GUS-positive and non-transformed(control) plants. Although gene silencing of the gus gene occurred,albeit at a low frequency (9.2%), the transformation systemdescribed generates large numbers of phenotypically normal,stably transformed plants. Copyright 2000 Annals of Botany Company Agrobacterium -mediated transformation, gene silencing, Solanum dulcamara L. (Bittersweet, Woody Nightshade), T-DNA truncation, transgene expression     

云斑车蝗线粒体基因组全序列测定与分析

采用长距 PCR 扩增及保守引物步移法并结合克隆测序测定并注释了云斑车蝗 Gastrimargus marmoratus （Thunberg）的线粒体基因组全序列。结果表明：云斑车蝗线粒体基因组全序列为15 904 bp（GenBank登录号为EU527334）,A+T含量略高于非洲飞蝗Locusta migratoria,为76.04%,包括13个蛋白质编码基因,22个tRNA 基因,2个rRNA基因和一段1 057 bp的A+T富集区。蛋白质基因的起始密码子中,除COⅠ和ND5为TTG以外,均为昆虫典型的起始密码子ATN。ND5基因使用了不完全终止密码子T,其余基因均为典型的TAA或TAG。预测了22个tRNA基因的二级结构,发现tRNA^Ser（AGN）缺少DHU臂, tRNA^Ser（UGY）的反密码子环上有9个碱基。预测了云斑车蝗12S和16S rRNA二级结构,分别包括3个结构域30个茎环和6个结构域44个茎环。A+T富集区含有3个串联重复序列。     

Complete Mitochondrial DNA Sequence of Conger myriaster (Teleostei: Anguilliformes): Novel Gene Order for Vertebrate Mitochondrial Genomes and the Phylogenetic Implications for Anguilliform Families

The complete nucleotide sequence of the mitochondrial genome was determined for a conger eel, Conger myriaster (Elopomorpha: Anguilliformes), using a PCR-based approach that employs a long PCR technique and many fish-versatile primers. Although the genome [18,705 base pairs (bp)] contained the same set of 37 mitochondrial genes [two ribosomal RNA (rRNA), 22 transfer RNA (tRNA), and 13 protein-coding genes] as found in other vertebrates, the gene order differed from that recorded for any other vertebrates. In typical vertebrates, the ND6, tRNA^Glu, and tRNA^Pro genes are located between the ND5 gene and the control region, whereas the former three genes, in C. myriaster, have been translocated to a position between the control region and the tRNA^Phe gene that are contiguously located at the 5′ end of the 12S rRNA gene in typical vertebrates. This gene order is similar to the recently reported gene order in four lineages of birds in that the latter lack the ND6, tRNA^Glu, and tRNA^Pro genes between the ND5 gene and the control region; however, the relative position of the tRNA^Pro to the ND6–tRNA^Glu genes in C. myriaster was different from that in the four birds, which presumably resulted from different patterns of tandem duplication of gene regions followed by gene deletions in two distantly related groups of organisms. Sequencing of the ND5–cyt b region in 11 other anguilliform species, representing 11 families, plus one outgroup species, revealed that the same gene order as C. myriaster was shared by another 4 families, belonging to the suborder Congroidei. Although the novel gene orders of four lineages of birds were indicated to have multiple independent origins, phylogenetic analyses using nucleotide sequences from the mitochondrial 12S rRNA and cyt b genes suggested that the novel gene orders of the five anguilliform families had originated in a single ancestral species. Received: 13 July 2000 / Accepted: 30 November 2000     

白鱼线粒体DNA控制区结构和种群遗传多样性分析

用特异性引物对白鱼（Anabarilius grahami）DNA进行PCR扩增,获得了白鱼线粒体DNA控制区基因全序列（930 bp）。控制区T、C、A和G碱基组成为29.8%、22.5%、33.0%和14.7%。对照其他已报道的鱼类控制区结构,对白鱼控制区结构进行了分析,识别了其终止序列区、中央保守区和保守序列区,找到了终止相关的序列TAS以及保守序列（CSB-F、CSB-D、CSB-1、CSB-2、CSB-3）。同时运用DNA分析软件对白鱼一个驯养种群（中国科学院昆明动物研究所珍稀鱼类繁育中心）及两个自然地理种群（江川县明星鱼洞、江川县牛摩村）进行了遗传多样性分析。结果显示：两个自然种群存在较强基因交流,未出现遗传分化;人工驯养种群遗传多样性最高,种群复壮程度较好。     

李氏大足蝗线粒体全基因组序列分析

采用长距PCR 扩增及保守引物步移法测定并注释了李氏大足蝗( Aeropus licenti Chang)的线粒体基因组全序列。结果表明,李氏大足蝗的线粒体基因组全长15 597 bp,A+T 含量为74.8%,37个基因位置与飞蝗的一致,基因间隔序列共计17处105 bp,间隔长度从1～21 bp不等;有10对基因间存在58 bp重叠,重叠碱基数在1～17 bp之间。13个蛋白质编码基因中找到4种可能的起始密码子;有12个基因在基因3'端找到完全的TAA或TAG 终止密码子,只有ND5基因终止密码子为不完整的T。除tRNA^Ser（AGN）外,其余21个tRNA基因的二级结构均属典型的三叶草结构。tRNA^Ser（AGN）的DHU臂缺失,在相应的位置上只形成一个环。预测的lrRNA二级结构总共有6个结构域（结构域Ⅲ缺失）,47个茎环结构;预测的srRNA的二级结构包含3个结构域,31个茎环结构。A+T 丰富区长度为712 bp。     

Complete mitochondrial genome of Odontobutis haifengensis (Perciformes,Odontobutiae): A unique rearrangement of tRNAs and additional non-coding regions identified in the genus Odontobutis

Herein, the complete mitochondrial genome of Odontobutis haifengensis was sequenced for the first time. The O. haifengensis mitogenome was 17,016 bp in length and included 13 protein-coding genes, 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a control region (CR). The genome organization, base composition, codon usage, and gene rearrangement was similar to other Odontobutis species. Furthermore, a tRNA gene rearrangement within the SLH cluster was found to be identical to other Odontobutis species. Moreover, the gene order and the positions of additional intergenic non-coding regions suggests that the observed unique gene rearrangement resulted from a tandem duplication and random loss of large-scale gene regions. Additionally, phylogenetic analysis showed that Odontobutis species form a monophyletic clade due to the conserved mitochondrial gene rearrangement. This study provides useful information that aids in a better understanding of mitogenomic diversity and evolutionary patterns of Odontobutidae species.     

Ribosomal RNA genes from Prevotella bryantii: organization and heterogeneity

FiveP. bryantii B₁4 16S rRNA gene copies and their flanking regions were cloned and analyzed. A genomic library was constructed and screened with oligonucleotide DNA probe specific for 16S rRNA gene ofP. bryantii. Five out of six different copies of 16S RNA gene were recovered and sequenced. Only minor differences (0.3–1.2%) between copies were detected within the 1541 bp long sequence. The impact of the sequence variability of 16S rRNA gene copies on phylogenetic positioning ofP. bryantii was determined. All five sequences from clonedP. bryantii B₁4 16S rRNA genes were placed in the same operational taxonomy unit. Control regions of all five analyzed rRNA operatons were almost identical and three candidate for promoter sequences were identified by Neutral Network Promoter Prediction. Spacer regions between 16S-rRNA and 23S rRNA genes in all five cloned copies were 543 bp long and genes for tRNA^lle and tRNA^Ala were identified inside this regions.     

Cloning and characterization of a microsatellite in the mitochondrial control region of the African side-necked turtle,Pelomedusa subrufa

The nucleotide sequence of the African side-necked turtle mitochondrial control region and its flanking tRNA genes was determined. This 73% A+T-rich region is 1194 bp long. Several conserved motifs involved in the regulation of the mitochondrial genome replication process, including one conserved sequence block (CSB1), and three termination-associated sequences were identified. The most remarkable feature found in this control region was the presence of six microsatellite-containing tandem repeats between the CSB1 motif and the tRNA^Phe gene. The potential usefulness of this microsatellite sequence for population-level studies is enhanced by its unique localization in the maternally inherited mitochondrial molecule.     

Complete Mitochondrial DNA Sequences of Six Snakes: Phylogenetic Relationships and Molecular Evolution of Genomic Features

Complete mitochondrial DNA (mtDNA) sequences were determined for representative species from six snake families: the acrochordid little file snake, the bold boa constrictor, the cylindrophiid red pipe snake, the viperid himehabu, the pythonid ball python, and the xenopeltid sunbeam snake. Thirteen protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 2 control regions were identified in these mtDNAs. Duplication of the control region and translocation of the tRNA^Leu gene were two notable features of the snake mtDNAs. The duplicate control regions had nearly identical nucleotide sequences within species but they were divergent among species, suggesting concerted sequence evolution of the two control regions. In addition, the duplicate control regions appear to have facilitated an interchange of some flanking tRNA genes in the viperid lineage. Phylogenetic analyses were conducted using a large number of sites (9570 sites in total) derived from the complete mtDNA sequences. Our data strongly suggested a new phylogenetic relationship among the major families of snakes: ((((Viperidae, Colubridae), Acrochordidae), (((Pythonidae, Xenopeltidae), Cylindrophiidae), Boidae)), Leptotyphlopidae). This conclusion was distinct from a widely accepted view based on morphological characters in denying the sister-group relationship of boids and pythonids, as well as the basal divergence of nonmacrostomatan cylindrophiids. These results imply the significance to reconstruct the snake phylogeny with ample molecular data, such as those from complete mtDNA sequences.[Reviewing Editor: Dr. Bill Ballard]     

Cloning and Characterization of New Repetitive Sequences in Field Bean (Vicia fabaL.)

Five new repetitive sequences have been isolated from theViciafabagenome, by cloning bands visible on agarose gel electrophoresisafter digestion of genomic DNA with various restriction enzymes.The sequences were 109 to 584 bp long, their abundance rangingfrom 5x10⁴to 5x10⁵copies per haploid genome. Southern blot andinsituhybridization revealed that four of five newly isolatedrepeats were dispersed in theV. fabagenome. One of the repeats(TIII15) showed tandem organization with several major hybridizationspots on mitotic chromosomesin situ.These sites were distributedin euchromatic as well as in heterochromatic chromosomal regions,and in several loci they were simultaneously localized withpreviously describedFokI repeated elements. The sequence ofTIII15 comprises four 26–27 bp subrepeats, but sharesno homology toFokI elements which have similar sequence organization.All newly described sequences were highly specific forV. faba,withlittle or no hybridization to DNA of otherViciaspecies, andno hybridization to DNA of other legumes tested.Copyright 1999Annals of Botany Company Vicia faba, field bean, repeated DNA sequences, FISH, PRINS, genome organization, copy number.