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

Unparalleled mitochondrial heteroplasmy and Wolbachia co-infection in the non-model bee, Amphylaeus morosus

Mitochondrial heteroplasmy is the occurrence of more than one type of mitochondrial DNA within a single individual. Although generally reported to occur in a small subset of individuals within a species, there are some instances of widespread heteroplasmy across entire populations. Amphylaeus morosus is an Australian native bee species in the diverse and cosmopolitan bee family Colletidae. This species has an extensive geographical range along the eastern Australian coast, from southern Queensland to western Victoria, covering approximately 2,000 km. Seventy individuals were collected from five localities across this geographical range and sequenced using Sanger sequencing for the mitochondrial cytochrome c oxidase subunit I (COI) gene. These data indicate that every individual had the same consistent heteroplasmic sites but no other nucleotide variation, suggesting two conserved and widespread heteroplasmic mitogenomes. Ion Torrent shotgun sequencing revealed that heteroplasmy occurred across multiple mitochondrial protein-coding genes and is unlikely explained by transposition of mitochondrial genes into the nuclear genome (NUMTs). DNA sequence data also demonstrated a consistent co-infection of Wolbachia across the A. morosus distribution with every individual infected with both bacterial strains. Our data are consistent with the presence of two mitogenomes within all individuals examined in this species and suggest a major divergence from standard patterns of mitochondrial inheritance. Because the host's mitogenome and the Wolbachia genome are genetically linked through maternal inheritance, we propose three possible hypotheses that could explain maintenance of the widespread and conserved co-occurring bacterial and mitochondrial genomes in this species.

High occurrence of length heteroplasmy in domestic Bactrian camel (Camelus bactrianus)

Heteroplasmy is the presence of more than one mitochondrial DNA (mtDNA) variant within a cell, tissue, or individual. In this study, sequence variation was investigated in the control region (CR) of mitochondrial DNA (mtDNA) from 135 individuals belonging to five primary domestic Bactrian camel breeds in China and Mongolia. Due to variation of the repeat unit G(T/C)(AC)n, length heteroplasmy was detected within each camel by direct sequencing and fragment analysis. A high occurrence of mtDNA heteroplasmy, up to 100 percentages was observed in five camel populations. Our study provides the first evidence of extensive length heteroplasmy in Bactrian camels.

The complete mitochondrial genome sequence of the tubeworm Lamellibrachia satsuma and structural conservation in the mitochondrial genome control regions of Order Sabellida

The control region of the mitochondrial genomes shows high variation in conserved sequence organizations, which follow distinct evolutionary patterns in different species or taxa. In this study, we sequenced the complete mitochondrial genome of Lamellibrachia satsuma from the cold-seep region of Kagoshima Bay, as a part of whole genome study and extensively studied the structural features and patterns of the control region sequences. We obtained 15,037 bp of mitochondrial genome using Illumina sequencing and identified the non-coding AT-rich region or control region (354 bp, AT=83.9%) located between trnH and trnR. We found 7 conserved sequence blocks (CSB), scattered throughout the control region of L. satsuma and other taxa of Annelida. The poly-TA stretches, which commonly form the stem of multiple stem-loop structures, are most conserved in the CSB-I and CSB-II regions. The mitochondrial genome of L. satsuma encodes a unique repetitive sequence in the control region, which forms a unique secondary structure in comparison to Lamellibrachia luymesi. Phylogenetic analyses of all protein-coding genes indicate that L. satsuma forms a monophyletic clade with L. luymesi along with other tubeworms found in cold-seep regions (genera: Lamellibrachia, Escarpia, and Seepiophila). In general, the control region sequences of Annelida could be aligned with certainty within each genus, and to some extent within the family, but with a higher rate of variation in conserved regions.

Comparative Mitogenomics of the Genus Odontobutis (Perciformes: Gobioidei: Odontobutidae) Revealed Conserved Gene Rearrangement and High Sequence Variations

To understand the molecular evolution of mitochondrial genomes (mitogenomes) in the genus Odontobutis, the mitogenome of Odontobutis yaluensis was sequenced and compared with those of another four Odontobutis species. Our results displayed similar mitogenome features among species in genome organization, base composition, codon usage, and gene rearrangement. The identical gene rearrangement of trnS-trnL-trnH tRNA cluster observed in mitogenomes of these five closely related freshwater sleepers suggests that this unique gene order is conserved within Odontobutis. Additionally, the present gene order and the positions of associated intergenic spacers of these Odontobutis mitogenomes indicate that this unusual gene rearrangement results from tandem duplication and random loss of large-scale gene regions. Moreover, these mitogenomes exhibit a high level of sequence variation, mainly due to the differences of corresponding intergenic sequences in gene rearrangement regions and the heterogeneity of tandem repeats in the control regions. Phylogenetic analyses support Odontobutis species with shared gene rearrangement forming a monophyletic group, and the interspecific phylogenetic relationships are associated with structural differences among their mitogenomes. The present study contributes to understanding the evolutionary patterns of Odontobutidae species.

Analysis of the complete mitochondrial genome and characterization of diverse NUMTs of Macaca leonina

As a non-human primate, the pig-tailed macaque has received wide attention because it can be infected by HIV-1. In this study, we determined the complete mtDNA sequence of the northern pig-tailed macaque (Macaca leonina). Unexpectedly, during the amplification of the mtDNA control region (D-loop region) we observed several D-loop-like sequences, which were NUMTs (nuclear mitochondrial sequences) and a total of 14 D-loop-like NUMT haplotypes were later identified in five individuals. The neighbor-joining tree and estimated divergence time based on these D-loop-like NUMT sequences of M. leonina provide some insights into the understanding of the evolutionary history of NUMTs. D-loop-like haplotypes G and H, which also exist in the nuclear genome of mulatta, appear to have been translocated into the nuclear genome before the divergence of M. mulatta and M. leonina. The other D-loop-like NUMT haplotypes were translocated into the nuclear genome of M. leonina after the divergence of the two species. Later sequence conversion was predicted to occur among these 14 D-loop-like NUMT haplotypes. The overall structure of the mtDNA of M. leonina was found to be similar to that seen in other mammalian mitochondrial genomes. Phylogenetic analysis based on the maximum likelihood method shows M. leonina clustered with Macaca silenus among the analyzed mammalian species.

Mitochondrial DNA variability to explore the relationship complexity of Schizothoracine (Teleostei: Cyprinidae)

Despite numerous studies on the taxonomy of a highly complex group of schizothoracine (snow trouts), with over five recognized species from Kashmir, India (Schizothorax niger, Schizothorax esocinus, Schizothorax plagiostomus, Schizothorax curvifrons and Schizothorax labiatus) based on traditional morphological data, the relationships between these species is poorly understood and the taxonomic validity is still under debate. To resolve the evolutionary relationships among these species, we sequenced mitochondrial fragments, including 16Sr RNA, Cytb and the D-loop. Separate analyses of 16S and Cytb showed intermixing of the species and 16S was found more conserved than Cytb. The D-loop was found highly variable and showed length variation between and within species. Length variation was observed in di-nucleotide (TA)n microsatellite repeats with a variable number of repeat units (n = 7-14) that did not show heteroplasmy. Central conserved sequence blocks (CSBs) in D-loop sequences were found comparable to other vertebrate species. All phylogenetic reconstructions recovered the focal taxa as a monophyletic clade within the schizothoracines. Analyses with combined mitochondrial data sets showed close genetic relationships of all the five species. In addition to a close relationship between S. niger and S. curvifrons, two distinct groupings of S. ecoscinus and S. plagiostomus were supported by all the analyses. This study gives an insight into molecular phylogeny of the species and improves our understanding of historical and taxonomic relationships derived from morphological and ecological studies.

Length variation and sequence divergence in mitochondrial control region of Schizothoracine (Teleostei: Cyperinidae) species

Schizothoracine fish commonly called snow trouts inhibit the entire network of snow and spring fed cool waters of Kashmir, India. Over 10 species reported earlier, only five species have been found, these include Schizothorax niger, Schizothorax esocinus, Schizothorax plagiostomus, Schizothorax curvifrons and Schizothorax labiatus. The relationship between these species is contradicting. To understand the evolutionary relation of these species, we examined the sequence information of mitochondrial D-loop of 25 individuals representing five species. Sequence alignment showed D-loop region highly variable and length variation was observed in di-nucleotide (TA)n microsatellite between and within species. Interestingly, all these species have (TA)n microsatellite not associated with longer tandem repeats at the 3' end of the mitochondrial control region and do not show heteroplasmy. Our analysis also indicates the presence of four conserved sequence blocks (CSB), CSB-D, CSB-1, CSB-II and CSB-III, four (Termination Associated Sequence) TAS motifs and 15bp pyrimidine block within the mitochondrial control region, that are highly conserved within genus Schizothorax when compared with other species. The phylogenetic analysis carried by Maximum likelihood (ML), Neighbor Joining (NJ) and Bayesian inference (BI) generated almost identical results. The resultant BI tree showed a close genetic relationship of all the five species and supports two distinct grouping of S. esocinus species. Besides the species relation, the presence of length variation in tandem repeats is attributed to differences in predicting the stability of secondary structures. The role of CSBs and TASs, reported so far as main regulatory signals, would explain the conservation of these elements in evolution.

Heteroplasmy and ancient translocation of mitochondrial DNA to the nucleus in the Chinese Horseshoe Bat (Rhinolophus sinicus) complex

The utility and reliability of mitochondrial DNA sequences in phylogenetic and phylogeographic studies may be compromised by widespread and undetected nuclear mitochondrial copies (numts) as well as heteroplasmy within individuals. Both numts and heteroplasmy are likely to be common across diverse taxa yet few studies have characterised their frequencies and variation at the intra-specific level. Here we report the presence of both numts and heteroplasmy in the mitochondrial control region of the Chinese horseshoe bat Rhinolophus sinicus. In total we generated 123 sequences from 18 bats, which contained two different numt clades (i.e. Numt-1 and Numt-2) and one mtDNA clade. The sequence divergence between Numt-1 and Numt-2 was 16.8% and each numt type was found in all four R. sinicus taxa, suggesting either two ancient translocations of mitochondrial DNA into the nucleus from the same source taxon, or a single translocation from different source taxa that occurred before the split of R. sinicus into different lineages. Within the mtDNA clade, phylogenetic relationships among the four taxa of R. sinicus were similar to those seen in previous results. Based on PCR comparisons, heteroplasmy was inferred between almost all individuals of R. sinicus with respect to sequence variation. Consistent with introgression of mtDNA between Central sinicus and septentrionalis, individuals from these two taxa exhibited similar signatures of repeated sequences in the control region. Our study highlights the importance of testing for the presence of numts and heteroplasmy when applying mtDNA markers to phylogenetic studies.

Lack of Structural Variation but Extensive Length Polymorphisms and Heteroplasmic Length Variations in the Mitochondrial DNA Control Region of Highly Inbred Crested Ibis, Nipponia nippon

The animal mitochondrial DNA (mtDNA) length polymorphism and heteroplasmy are accepted to be universal. Here we report the lack of structural variation but the presence of length polymorphism as well as heteroplasmy in mtDNA control region of an endangered avian species - the Crested Ibis (Nipponia nippon). The complete control region was directly sequenced while the distribution pattern and inheritance of the length variations were examined using both direct sequencing and genotyping of the PCR fragments from captive birds with pedigrees, wild birds and a historical specimen. Our results demonstrated that there was no structural variation in the control region, however, different numbers of short tandem repeats with an identical motif of CA3CA2CA3 at the 3'-end of the control region determined the length polymorphisms among and heteroplasmy within individual birds. There were one to three predominant fragments in every bird; nevertheless multiple minor fragments coexist in all birds. These extremely high polymorphisms were suggested to have derived from the 'replication slippage' of a perfect microsatellite evolution following the step-wise mutational model. The patterns of heteroplasmy were found to be shifted between generations and among siblings but rather stable between blood and feather samples. This study provides the first evidence of a very extensive mtDNA length polymorphism and heteroplasmy in the highly inbred Crested Ibis which carries an mtDNA genome lack of structural genetic diversity. The analysis of pedigreed samples also sheds light on the transmission of mtDNA length heteroplasmy in birds following the genetic bottleneck theory. Further research focusing on the generation and transmission of particular mtDNA heteroplasmy patterns in single germ line of Crested Ibis is encouraged by this study.

Limited phylogenetic distribution of a long tandem-repeat cluster in the mitochondrial control region in Bubo (Aves, Strigidae) and cluster variation in Blakiston's fish owl (Bubo blakistoni)

To investigate the phylogenetic position of Blakiston's fish owl (Bubo blakistoni), we sequenced the mitochondrial (mt) DNA control region and cytochrome b (cyt b) for nine Bubo species. Maximum-likelihood analyses of combined control region and cyt b sequences, and cyt b sequences alone, showed that species formerly placed in genus Ketupa comprise a monophyletic group. Unexpectedly, we discovered a long cluster of 20-25 tandem repeat units 77 or 78bp long in the third control region domain in four of the nine Bubo species for which the control region was sequenced (B. blakistoni, B. flavipes, and B. ketupu in the Ketupa clade; B. lacteus), leading to overall control region lengths of 3.0-3.8kpb estimated from agarose gel electrophoresis. The control region in B. lacteus is the longest (3.8kbp) reported to date in vertebrates. Sequencing of eight repeat units at each end of the cluster in 20 B. blakistoni individuals detected several types of repeat units 77 or 78bp long, and six patterns in the order of unit types. The occurrence of a repeat cluster in all three species examined in the Ketupa clade suggests their common ancestor also had a cluster, whereas a maximum parsimony tree showed repeat-unit types grouping by species, rather than by paralog groups, suggesting independent origins of the clusters. We reconcile these results with a turnover model, in which the range in cluster-length variation and unit types at the 5' end are hypothetically functionally constrained by the protein-binding function of the control region, but otherwise there is a continual turnover of units in evolutionary time, with new unit types arising through mutations, proliferating by duplication of single and double repeat blocks, and being lost through deletion. Estimated free energies for reconstructed secondary structures of single and especially pairs of repeat units were higher than for homologous single-unit blocks in species lacking a repeat cluster, supporting slipped-strand mispairing as the mechanism of cluster turnover.

The complete mitochondrial genome analysis of the tiger (Panthera tigris)

The complete mitochondrial genomes of five tiger samples from three subspecies (P. t. sumatrae, P. t. altica, and P. t. tigris) were successfully obtained by using 26 specifically designed Panthera-specific primer sets. The genome organization and gene arrangement of the five tiger samples were similar to each other; however polymorphic tandem repeat sequences were observed in the control region (CR). This led to a difference in the genome lengths obtained from these five samples with an average size of 16,994 bp for the five tiger mitochondrial genomes. The nucleotide base composition was on average as follows: A, 31.8%; T, 27.0%; C, 26.6%; G, 14.6% and exhibited compositional asymmetry. Most of tiger mitochondrial genome characteristics are similar to those of other common vertebrate species; however, some distinctive features were observed in the CR. First, the repetitive sequence 2 (RS 2) contained two repeat units of 80 bp and the first 15 bp of what would be the third repeat motif. The repetitive sequence 3 (RS 3) contained 47-50 repeat motifs of a shorter 8 bp (ACGTAYAC)(n). Second, length heteroplasmy polycystosine (poly-C) stretches was observed at the end of the HV I locus in all tiger samples.

Comparative analysis of complete mitochondrial DNA control region of four species of Strigiformes

The sequence of the whole mitochondrial (mt) DNA control region (CR) of four species of Strigiformes was obtained. Length of the CR was 3,290 bp, 2,848 bp, 2,444 bp, and 1,771 bp for Asio flammeus, Asio otus, Athene noctua, and Strix aluco, respectively. Interestingly, the length of the control region was maximum in Asio flammeus among all the avian mtDNA control regions sequenced thus far. In addition, the base composition and organization of mtDNA CR of Asio flammeus were identical to those reported for other birds. On the basis of the differential frequencies of base substitutions, the CR may be divided two variable domains, I and III, and a central conserved domain, II. The 3' end of the CR contained many tandem repeats of varying lengths and repeat numbers. In Asio flammeus, the repeated sequences consisted of a 126 bp sequence that was repeated seven times and a 78 bp sequence that was repeated 14 times. In Asio otus, there were also two repeated sequences, namely a 127 bp sequence that was repeated eight times and a 78 bp sequence that was repeated six times. The control region of Athene noctua contained three sets of repeats: a 89 bp sequence that was repeated three times, a 77 bp sequence that was repeated four times, and a 71 bp sequence that was repeated six times. Strix aluco, however, had only one repeated sequence, a 78 bp sequence that was repeated five times. The results of this study seem to indicate that these tandem repeats may have resulted from slipped-strand mispairing during mtDNA replication. Moreover, there are many conserved motifs within the repeated units. These sequences could form stable stem-loop secondary structures, which suggests that these repeated sequences play an important role in regulating transcription and replication of the mitochondrial genome.

The complete mitochondrial genome of the house dust mite Dermatophagoides pteronyssinus (Trouessart): a novel gene arrangement among arthropods

BACKGROUND

The apparent scarcity of available sequence data has greatly impeded evolutionary studies in Acari (mites and ticks). This subclass encompasses over 48,000 species and forms the largest group within the Arachnida. Although mitochondrial genomes are widely utilised for phylogenetic and population genetic studies, only 20 mitochondrial genomes of Acari have been determined, of which only one belongs to the diverse order of the Sarcoptiformes. In this study, we describe the mitochondrial genome of the European house dust mite Dermatophagoides pteronyssinus, the most important member of this largely neglected group.

RESULTS

The mitochondrial genome of D. pteronyssinus is a circular DNA molecule of 14,203 bp. It contains the complete set of 37 genes (13 protein coding genes, 2 rRNA genes and 22 tRNA genes), usually present in metazoan mitochondrial genomes. The mitochondrial gene order differs considerably from that of other Acari mitochondrial genomes. Compared to the mitochondrial genome of Limulus polyphemus, considered as the ancestral arthropod pattern, only 11 of the 38 gene boundaries are conserved. The majority strand has a 72.6% AT-content but a GC-skew of 0.194. This skew is the reverse of that normally observed for typical animal mitochondrial genomes. A microsatellite was detected in a large non-coding region (286 bp), which probably functions as the control region. Almost all tRNA genes lack a T-arm, provoking the formation of canonical cloverleaf tRNA-structures, and both rRNA genes are considerably reduced in size. Finally, the genomic sequence was used to perform a phylogenetic study. Both maximum likelihood and Bayesian inference analysis clustered D. pteronyssinus with Steganacarus magnus, forming a sistergroup of the Trombidiformes.

CONCLUSION

Although the mitochondrial genome of D. pteronyssinus shares different features with previously characterised Acari mitochondrial genomes, it is unique in many ways. Gene order is extremely rearranged and represents a new pattern within the Acari. Both tRNAs and rRNAs are truncated, corroborating the theory of the functional co-evolution of these molecules. Furthermore, the strong and reversed GC- and AT-skews suggest the inversion of the control region as an evolutionary event. Finally, phylogenetic analysis using concatenated mt gene sequences succeeded in recovering Acari relationships concordant with traditional views of phylogeny of Acari.

Microsatellite evolution in the mitochondrial genome of Bechstein's bat (Myotis bechsteinii)

Being highly polymorphic, microsatellites are widely used genetic markers. They are abundant throughout the nuclear genomes of eukaryotes but rare in the mitochondrial genomes (mtDNA) of animals. We describe a short but highly polymorphic AT microsatellite in the mtDNA control region of Bechstein's bat and discuss the role of mutation, genetic drift, and selection in maintaining its variability. As heteroplasmy and hence mutation rate were positively correlated with repeat number, a simple mutation model cannot explain the observed frequency distribution of AT copy numbers. Because of the unimodal distribution of repeat numbers found in heteroplasmic individuals, single step mutations are likely to be the predominant mechanism of copy number alternations. Above a certain copy number (seven repeats), deletions of single dinucleotide repeats seem to be more common than additions, which results in a decrease in frequency of long alleles. Heteroplasmy was inherited from mothers to their offspring and no evidence of paternal inheritance of mitochondria was found. Genetic differences accumulated with more distant ancestry, which suggests that microsatellites can be useful genetic markers in population genetics.

DNA sequence variation in the mitochondrial control region of red-backed voles (Clethrionomys)

The complete mitochondrial DNA (mtDNA) control region was sequenced for 71 individuals from five species of the rodent genus Clethrionomys both to understand patterns of variation and to explore the existence of previously described domains and other elements. Among species, the control region ranged from 942 to 971 bp in length. Our data were compatible with the proposal of three domains (extended terminal associated sequences [ETAS], central, conserved sequence blocks [CSB]) within the control region. The most conserved region in the control region was the central domain (12% of nucleotide positions variable), whereas in the ETAS and CSB domains, 22% and 40% of nucleotide positions were variable, respectively. Tandem repeats were encountered only in the ETAS domain of Clethrionomys rufocanus. This tandem repeat found in C. rufocanus was 24 bp in length and was located at the 5' end of the control region. Only two of the proposed CSB and ETAS elements appeared to be supported by our data; however, a "CSB1-like" element was also documented in the ETAS domain.

The evolutionary history of the coral genus Acropora (Scleractinia, Cnidaria) based on a mitochondrial and a nuclear marker: reticulation, incomplete lineage sorting, or morphological convergence?

This study examines molecular relationships across a wide range of species in the mass spawning scleractinian coral genus Acropora. Molecular phylogenies were obtained for 28 species using DNA sequence analyses of two independent markers, a nuclear intron and the mtDNA putative control region. Although the compositions of the major clades in the phylogenies based on these two markers were similar, there were several important differences. This, in combination with the fact that many species were not monophyletic, suggests either that introgressive hybridization is occurring or that lineage sorting is incomplete. The molecular tree topologies bear little similarity to the results of a recent cladistic analysis based on skeletal morphology and are at odds with the fossil record. We hypothesize that these conflicting results may be due to the same morphology having evolved independently more than once in Acropora and/or the occurrence of extensive interspecific hybridization and introgression in combination with morphology being determined by a small number of genes. Our results indicate that many Acropora species belong to a species complex or syngameon and that morphology has little predictive value with regard to syngameon composition. Morphological species in the genus often do not correspond to genetically distinct evolutionary units. Instead, species that differ in timing of gamete release tend to constitute genetically distinct clades.

Sequence variation in the mitochondrial DNA control region of wild African cheetahs (Acinonyx jubatus)

Five hundred and twenty-five bp of mitochondrial control region were sequenced and analysed for 20 Acinonyx jubatus and one Felis catus. These sequences were compared with published sequences from another domestic cat, 20 ocelots (Leopardus pardalus) and 11 margays (Leopardus weidii). The intraspecific population divergence in cheetahs was found to be less than in the other cats. However variation was present and distinct groups of cheetahs were discernible. The 80 bp RS2 repetitive sequence motif previously described in other felids was found in four copies in cheetah. The repeat units probably have the ability to form secondary structure and may have some function in the regulation of control region replication. The two central repeat units in cheetah show homogenization that may have arisen by convergent evolution.

An empirical genetic assessment of the severity of the northern elephant seal population bottleneck

A bottleneck in population size of a species is often correlated with a sharp reduction in genetic variation. The northern elephant seal (Mirounga angustirostris) has undergone at least one extreme bottleneck, having rebounded from 20-100 individuals a century ago to over 175,000 individuals today. The relative lack of molecular-genetic variation in contemporary populations has been documented, but the extent of variation before the late 19th century remains unknown. We have determined the nucleotide sequence of a 179 base-pair segment of the mitochondrial DNA (mtDNA) control region from seals that lived before, during and after a bottleneck low in 1892. A 'primerless' PCR was used to improve the recovery of information from older samples. Only two mtDNA genotypes were present in all 150+ seals from the 1892 bottleneck on, but we discovered four genotypes in five pre-bottleneck seals. This suggests a much greater amount of mtDNA genotypic variation before this bottleneck, and that the persistence of two genotypes today is a consequence of random lineage sampling. We cannot correlate the loss of mtDNA genotypes with a lowered mean fitness of individuals in the species today. However, we show that the species historically possessed additional genotypes to those present now, and that sampling of ancient DNA could elucidate the genetic consequences of severe reductions in population size.

mtDNA tandem repeats in domestic dogs and wolves: mutation mechanism studied by analysis of the sequence of imperfect repeats

The mitochondrial (mt) DNA control region (CR) of dogs and wolves contains an array of imperfect 10 bp tandem repeats. This region was studied for 14 domestic dogs representing the four major phylogenetic groups of nonrepetitive CR and for 5 wolves. Three repeat types were found among these individuals, distributed so that different sequences of the repeat types were formed in different molecules. This enabled a detailed study of the arrays and of the mutation events that they undergo. Extensive heteroplasmy was observed in all individuals; 85 different array types were found in one individual, and the total number of types was estimated at 384. Among unrelated individuals, no identical molecules were found, indicating a high rate of evolution of the region. By performing a pedigree analysis, array types which had been inherited from mother to offspring and array types which were the result of somatic mutations, respectively, could be identified, showing that about 20% of the molecules within an individual had somatic mutations. By direct pairwise comparison of the mutated and the original array types, the physiognomy of the inserted or deleted elements (indels) and the approximate positions of the mutations could be determined. All mutations could be explained by replication slippage or point mutations. The majority of the indels were 1-5 repeats long, but deletions of up to 17 repeats were found. Mutations were found in all parts of the arrays, but at a higher frequency in the 5' end. Furthermore, the inherited array types within the mother-offspring pair were aligned and compared so that germ line mutations could be studied. The pattern of the germ line mutations was approximately the same as that of the somatic mutations.

Individual identification of cats and dogs using mitochondrial DNA tandem repeats?

Cats and dogs are very common animals in the human environment. In Switzerland, one in five households owns a cat or a dog. Their hairs are very easily transferred and could be used as a frequent trace evidence. Using DNA analysis, identification of these animals should be possible as it is in human identification. However, most of the time, no nuclear DNA can be recovered from the hair. It is therefore necessary to rely on mtDNA. Cats and dogs have tandemly repeated sequences in their mtDNA control region. In this study, the authors show that these tandem repeats are very polymorphic but are also the source of a very high level of heteroplasmy. The authors, therefore, examined if this might prevent their use in forensic identification.

Tandemly repeated sequences in the mitochondrial DNA control region and phylogeography of the Pike-Perches Stizostedion

DNA sequences from the mitochondrial DNA control region are used to test the phylogeographic relationships among the pike-perches, Stizostedion (Teleostei: Percidae) and to examine patterns of variation. Sequences reveal two types of variability: single nucleotide polymorphisms and 6 to 14 copies of 10- to 11-base-pair tandemly repeated sequences. Numbers of copies of the tandem repeats are found to evolve too rapidly to detect phylogenetic signal at any taxonomic level, even among populations. Sequence similarities of the tandem repeats among Stizostedion and other percids suggest concerted evolutionary processes. Predicted folding of the tandem repeats and their proximity to termination-associated sequences indicate that secondary structure mediates slipped-strand mispairing among the d-loop, heavy, and light strands. Neighbor-joining and maximum parsimony analyses of sequences indicate that the genus is divided into clades on the continents of North America and Eurasia. Calibrating genetic distances with divergence times supports the hypothesis that Stizostedion dispersed from Eurasia to North America across a North Pacific Beringial land bridge approximately 4 million years before present, near the beginning of the Pliocene Epoch. The North American S. vitreum and S. canadense appear separated by about 2.75 million years, and the Eurasian S. lucioperca and S. volgensis are diverged by about 1.8 million years, suggesting that speciation occurred during the late Pliocene Epoch.

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 tRNAPhe gene. The potential usefulness of this microsatellite sequence for population-level studies is enhanced by its unique localization in the maternally inherited mitochondrial molecule.

Modelling the secondary structures of slippage-prone hypervariable RNA regions: the example of the tiger beetle 18S rRNA variable region V4

Variable regions within ribosomal RNAs frequently vary in length as a result of incorporating products of slippage. This makes constructing secondary structure models problematic because base homology is difficult or impossible to establish between species. Here, we model such a region by comparing the results of the MFOLD suboptimal folding algorithm for different species to identify conserved structures. Based on the reconstruction of base change on a phylogenetic tree of the species and comparison against null models of character change, we devise a statistical analysis to assess support of these structures from compensatory and semi-compensatory (i.e. G.C to G.U or A.U to G.U) mutations. As a model system we have used variable region V4 from cicindelid (tiger beetle) small subunit ribosomal RNAs (SSU rRNAs). This consists of a mixture of conserved and highly variable subregions and has been subject to extensive comparative analysis in the past. The model that results is similar to a previously described model of this variable region derived from a different set of species and contains a novel structure in the central, highly variable part. The method we describe may be useful in modelling other RNA regions that are subject to slippage.

Mammalian mitochondrial D-loop region structural analysis: identification of new conserved sequences and their functional and evolutionary implications

This paper reports the first comprehensive analysis of Displacement loop (D-loop) region sequences from ten different mammalian orders. It represents a systematic evolutionary study at the molecular level on regulatory homologous regions in organisms belonging to a well defined class, mammalia, which radiated about 150 million years ago (Mya). We have aligned and analyzed 26 complete D-loop region sequences available in the literature and the fat dormouse sequence, recently determined in our laboratory. The novelty of our alignment consists of the extensive manual revision of the preliminary output obtained by computer program to optimize sequence similarity, particularly for the two peripheral domains displaying heterogeneity in length and the presence of repeated sequences. The multialignment is available at the WWW site: http://www.ba.cnr.it/dloop.html. Our comparative study has allowed us to identify new conserved sequence blocks present in all the species under consideration and events of insertion/deletion which have important implications in both functional and evolutionary aspects. In particular we have detected two blocks, about 60 bp long, extended termination associated sequences (ETAS1 and ETAS2) conserved in all the organisms considered. Evaluation against experimental work suggests a possible functional role of ETAS1 and ETAS2 in the regulation of replication and transcription and targeted experimental approaches. The analyses on conserved sequence blocks (CSBs) clearly indicate that CSB1 is the only very essential element, common to all mammalian mt genomes, while CSB2 and CSB3 could be involved in different though related functions, probably species specific, and thus more linked to nuclear mitochondrial coevolutionary processes. Our hypothesis on the different functional implications of the conserved elements, CSBs and TASs, reported so far as main regulatory signals, would explain the different conservation of these elements in evolution. Moreover the intra-order comparison of the D-loop regions highlights peculiar features useful to define the evolutionary dynamics of this region in closely related species.

Extreme sequence heteroplasmy in bat mitochondrial DNA

Mitochondrial heteroplasmy is shown to be extensive when amplification products from the mitochondrial control region are cloned and sequenced from a European bat species. In contrast, a mitochondrial ribosomal RNA gene does not exhibit substantial levels of heteroplasmy when analyzed in an identical way. In the bat, heteroplasmy with respect to length as well as sequence seems to be transmitted from mother to offspring. Thus, the intra-individual sequence diversity seems to accumulate within the female germ line and its extent to be controlled primarily by purifying selection. Similar experiments in humans and a marsupial suggest that heteroplasmy may not be as uncommon among mammals as hitherto thought.

The presence of tandem repeats and the initiation of replication in rabbit mitochondrial DNA

The non-coding region of rabbit mitochondrial DNA (mtDNA) exhibits two sets of tandem repeats between conserved sequence block 1 (CSB1) and the tRNA Phe gene. Both repeated sequences, short repeated (SR) and long repeats (LR), which contain 20 and 153 nucleotides, respectively, are involved in the generation of a high degree of mitochondrial heteroplasmy. Due to the location of these sequences in the regulatory region and their properties in terms of variable conformations, they could affect the initiation of replication of the heavy-strand DNA (H-strand DNA) and subsequently would influence the efficiency of mtDNA replication. The extremities of the displacement loop (D-loop) DNA strand and the 5' ends of RNA primers initiating the H-strand DNA synthesis were characterized in individual rabbits. Mapping at the nucleotide level of the 5' and 3' ends of the D-loop DNA strands indicates that both extremities are heterogeneous. The H-strand replication origin OH is located close to the conserved sequence block CSB1 as in other mammals. In all of the individuals studied so far, DNA molecules with a 5' end 1-2 nucleotides downstream of CSB1 were always present. As H-strand DNA replication is believed to be primed by RNA transcribed from the light-strand promoter (LSP), RNA mapping was carried out to identify the 5' end of H-strand RNA. Neighbouring initiation sites were identified at the nucleotide level in an (A+T)-rich region at nucleotide 1841 and in a stretch of cytosine residues at nucleotides 1849-1852, which are located at the beginning of the first long repeat. A detailed RNA analysis indicates that H-strand RNA molecules are initiated in each long repeat. The amplification of the regulatory region has produced multiple initiation transcription sites and a family of RNA primers of various lengths. These variations in length and the ensuing secondary structures are not critical for their potential function as H-strand DNA replication primers.

Sequence analysis of the complete mitochondrial DNA molecule of the hedgehog, Erinaceus europaeus, and the phylogenetic position of the Lipotyphla

The sequence of the mitochondrial DNA (mtDNA) molecule of the European hedgehog (Erinaceus europaeus) was determined. The length of the sequence presented is 17,442 nucleotides (nt). The molecule is thus the largest eutherian mtDNA molecule so far reported. The organization of the molecule conforms with that of other eutherians, but the control region of the molecule is exceptionally long, 1,988 nt, due to the presence of repeated motifs at two different positions in the 3' part of the control region. The length of the control region is not absolute due to pronounced heteroplasmy caused by variable numbers of the motif TACGCA in one of the repetitive regions. The sequence presented includes 46 repeats of this type. The other repeated region is composed of different AT-rich repeats. This region was identical among four clones studied. Comparison of mitochondrial peptide-coding genes identified a separate position of the hedgehog among several mammalian orders. The concatenated protein sequence of the 13 peptide-coding genes was used in a phylogenetic study using the opossum as outgroup. The position of the hedgehog sequence was basal among the other eutherian sequences included: human, rat, mouse, cow, blue whale, harbor seal, and horse. The analysis did not resolve the relationship among carnivores, perissodactyls, and artiodactyls/cetaceans, suggesting a closer relationship among these orders than acknowledged by classical approaches.

SSCP analysis of pig mitochondrial DNA D-loop region polymorphism

The sequence polymorphism that occurs in the mitochondrial DNA (mtDNA) displacement (D)-loop region is useful as a cytoplasmic DNA marker. We cloned the mtDNA D-loop regions of five breeds of pig by polymerase chain reaction (PCR) and determined their sequences. The sequence diversities in D-loop regions among five breeds of pig were located in the starting area of heavy-strand replication. From these sequences, we designed primers for PCR-mediated single-strand conformation polymorphism (PCR-SSCP) analysis that amplified the most polymorphic 227 bp fragment of the D-loop region. The results of PCR-SSCP analysis clearly showed that four types of polymorphism (A to D) are found in Landrace (A), Large White (A, B), Duroc (A), Göttingen miniature pig (B) and Meishan (C, D). The same polymorphisms were also detected from each porcine embryo by this method. Our results show that PCR-SSCP analysis is useful in detecting polymorphisms in the D-loop region of pigs and pig embryos.

Variable numbers of simple tandem repeats make birds of the order ciconiiformes heteroplasmic in their mitochondrial genomes

We have analyzed a variable domain of the mitochondrial DNA control region of 18 avian species. Intra-individual length variation was identified and characterized in 15 species. The occurrence of heteroplasmy among species is phylogenetically consistent with a current classification of birds. Polymerase chain reaction amplifications, direct sequencing, and Southern analysis of mitochondrial DNA showed that the heteroplasmy is due to variable numbers of direct repeats in a tandem organization, located in the control region close to the tRNAPhe gene. The tandem repeats consist of short sequence motifs that vary in size from 4 to 32 base pairs between species. Sequence complexity of the repeat motifs was low, with almost exclusively Ts and Gs in the heavy-strand. Extensive variation in the copy number of the repeats was seen both intra-specifically and within individuals. This is the first report of mitochondrial heteroplasmy characterized at the sequence level in birds.

The complete mitochondrial DNA sequence of the horse, Equus caballus: extensive heteroplasmy of the control region

The sequence of the mitochondrial (mt) DNA of the horse (Equus caballus) was determined. The length of the sequence presented is 16,660 bp. This figure, however, is not absolute due to pronounced heteroplasmy caused by variable numbers of the motif GTGCACCT in the control region of different molecules. Boundaries of the 13 peptide-coding genes were determined by the presence of start and stop codons, and by analogy with other eutherian mtDNAs. Three genes (COIII, NADH3 and NADH4) were not terminated by a stop codon. Comparison among the peptide-coding genes of the horse and eight other mammals suggests that the boundaries of some mt genes should be redefined. The number of repeats in the control region was determined by sequencing 77 different clones (20 direct plus 57 PCR clones). The number of repeats ranged from 2 to 29. There was a pronounced overrepresentation of clones with many repeats (22-27). Very few clones had a repeat number that was close to the mean number of repeats.

Rapid evolution of a heteroplasmic repetitive sequence in the mitochondrial DNA control region of carnivores

We describe a repetitive DNA region at the 3' end of the mitochondrial DNA (mtDNA) control region and compare it in 21 carnivore species representing eight carnivore families. The sequence and organization of the repetitive motifs can differ extensively between arrays; however, all motifs appear to be derived from the core motif "ACGT." Sequence data and Southern blot analysis demonstrate extensive heteroplasmy. The general form of the array is similar between heteroplasmic variants within an individual and between individuals within a species (varying primarily in the length of the array, though two clones from the northern elephant seal are exceptional). Within certain families, notably ursids, the array structure is also similar between species. Similarity between species was not apparent in other carnivore families, such as the mustelids, suggesting rapid changes in the organization and sequence of some arrays. The pattern of change seen within and between species suggests that a dominant mechanism involved in the evolution of these arrays is DNA slippage. A comparative analysis shows that the motifs that are being reiterated or deleted vary within and between arrays, suggesting a varying rate of DNA turnover. We discuss the evolutionary implications of the observed patterns of variation and extreme levels of heteroplasmy.

Numt, a recent transfer and tandem amplification of mitochondrial DNA to the nuclear genome of the domestic cat

The mitochondrial DNA of plant and animal cells is a transcriptionally active genome that traces its origins to a symbiotic infection of eucaryotic cells by bacterial progenitors. As prescribed by the Serial Endosymbiosis Theory, symbiotic organelles have gradually transferred their genes to the eucaryotic genome, producing a functional interaction of nuclear and mitochondrial genes in organelle function. We report here a recent remarkable transposition of 7.9 kb of a typically 17.0-kb mitochondrial genome to a specific nuclear chromosomal position in the domestic cat. The intergrated segment has subsequently become amplified 38-76 times and now occurs as a tandem repeat macrosatellite with multiple-length alleles resolved by pulse-field gel electrophoresis (PFGE) segregating in cat populations. Sequence determination of the nuclear mitochondrial DNA segment, Numt, revealed a d(CA)-rich 8-bp motif [ACACACGT] repeated imperfectly five times at the deletion junction that is a likely target for recombination. The extent and pattern of sequence divergence of Numt genes from the cytoplasmic mtDNA homologues plus the occurrence of Numt in other species of the family Felidae allowed an estimate for the origins of Numt at 1.8-2.0 million years ago in an ancestor of four modern species in the genus Felis. Numt genes do not function in cats; rather, the locus combines properties of nuclear minisatellites and pseudogenes. These observations provide an empirical glimpse of historic genomic events that may parallel the accommodation of organelles in eucaryotes.

Polymorphic sequence in the D-loop region of equine mitochondrial DNA

The D-loop regions in equine mitochondrial DNA were cloned from three thoroughbred horses by polymerase chain reaction (PCR). The total number of bases in the D-loop region were 1114 bp, 1115 bp and 1146 bp. The equine D-loop region is A/T rich like many other mammalian D-loops. The large central conserved sequence block and small conserved sequence blocks 1, 2 and 3, that are common to other mammals, were observed. Between conserved sequence blocks 1 and 2 there were tandem repeats of an 8 bp equine-specific sequence TGTGCACC, and the number of tandem repeats differed among individual horses. The base composition in the unit of these repeats is G/C rich as are the short repeats in the D-loops of rabbit and pig. Comparing DNA sequences between horse and other mammals, the difference in the D-loop region length is mostly due to the difference in the number of DNA sequences at both extremities. The similarities of the DNA sequences are in the middle part of the D-loop. In comparison of the sequences among three thoroughbred horses, it was determined that the region between tRNA(Pro) and the large central conserved sequence block was the richest in variation. PCR primers in the D-loop region were designed and the expected maternal inheritance was confirmed by PCR-RFLP (restriction fragment length polymorphism).

Evolution of genetic redundancy for advanced players

An ever expanding database on the sequence organization and repetition of genic and non-genic components of nuclear and organelle genomes reveals that the vast majority of sequences are subject to one or other mechanism of DNA turnover (gene conversion, unequal crossing over, slippage, retrotransposition, transposition and others). Detailed studies, using novel methods of experimental detection and analytical procedures, show that such mechanisms can operate one on top of another and that wide variations in their unit lengths, biases, polarities and rates create bizarre and complex patterns of genetic redundancy. The ability of these mechanisms to operate both within and between chromosomes implies that realistic models of the evolutionary dynamics of redundancy, and of the potential interaction with natural selection in a sexual species, need to consider the diffusion of variant repeats across multiple chromosome lineages, in a population context. Recently, important advances in both experimental and analytical approaches have been made along these lines. There is increasing awareness that genetic redundancy and turnover induces a molecular co-evolution between functionally interacting genetic systems in order to maintain essential functions.

Evolution by DNA turnover in the control region of vertebrate mitochondrial DNA

The control region of animal mitochondrial DNA is heterogeneous, including both highly conserved and highly variable sequences. Within the variable regions, variable number tandem repeat sequences have been described for numerous species. Repeats at one location, just upstream of the origin of replication, show an unprecedented level of length variation in somatic tissue. Recent comparison of these sequences in different species indicates a pattern of DNA turnover acting at different rates and over motifs of various sizes.