Comparisons of ape and human sequences that regulate mitochondrial DNA transcription and D-loop DNA synthesis

Mitogenome sequences of domestic cats demonstrate lineage expansions and dynamic mutation processes in a mitochondrial minisatellite

BACKGROUND

As a population genetic tool, mitochondrial DNA is commonly divided into the ~ 1-kb control region (CR), in which single nucleotide variant (SNV) diversity is relatively high, and the coding region, in which selective constraint is greater and diversity lower, but which provides an informative phylogeny. In some species, the CR contains variable tandemly repeated sequences that are understudied due to heteroplasmy. Domestic cats (Felis catus) have a recent origin and therefore traditional CR-based analysis of populations yields only a small number of haplotypes.

RESULTS

To increase resolution we used Nanopore sequencing to analyse 119 cat mitogenomes via a long-amplicon approach. This greatly improves discrimination (from 15 to 87 distinct haplotypes in our dataset) and defines a phylogeny showing similar starlike topologies within all major clades (haplogroups), likely reflecting post-domestication expansion. We sequenced RS2, a CR tandem array of 80-bp repeat units, placing RS2 array structures within the phylogeny and increasing overall haplotype diversity. Repeat number varies between 3 and 12 (median: 4) with over 30 different repeat unit types differing largely by SNVs. Five SNVs show evidence of independent recurrence within the phylogeny, and seven are involved in at least 11 instances of rapid spread along repeat arrays within haplogroups.

CONCLUSIONS

In defining mitogenome variation our study provides key information for the forensic genetic analysis of cat hair evidence, and for the first time a phylogenetically informed picture of tandem repeat variation that reveals remarkably dynamic mutation processes at work in the mitochondrion.

Indian Red Jungle fowl reveals a genetic relationship with South East Asian Red Jungle fowl and Indian native chicken breeds as evidenced through whole mitochondrial genome sequences

Background: Native chickens are dispersed in a wide geographical range and have hereditary assets that are kept by farmers for various purposes. Mitochondrial DNA (mtDNA) is a widely utilized marker in molecular studies because of its quick advancement, matrilineal legacy, and simple molecular structure. Method and Results: We performed NGS sequencing to investigate mitochondrial genomes and to evaluate the hereditary connections, diversity, and measure of gene stream estimation in Indian native chicken breeds and Red Jungle fowl. The chicken breeds were genotyped using the D-loop region and 23 haplotypes were identified. When compared to Indian native breeds, more haplotypes were identified in the NADH dehydrogenase subunits, Cytochrome c oxidase, Cytochrome b, ATP synthase subunit 6, and Ribosomal RNA genes. The phylogenetic examination indicated that the analyzed chicken breeds were divided into six significant clades, namely A, B, C, D, E, and F, of which the F clade indicated the domestication of chicken breeds in India. Additionally, our work affirmed that the Indian Red Jungle Fowl is the origin for both reference Red Jungle Fowl as well as all Indian breeds, which is reflected in the dendrogram as well as network analysis based on the whole mtDNA and D-loop region. Indian Red Jungle Fowl is distributed as an outgroup, suggesting that this ancestry was reciprocally monophyletic. Conclusion: The mtDNA sequences of Indian native chickens provided novel insights into adaptation mechanisms and the significance of important mtDNA variations in understanding the maternal lineages of native birds.

Comparative mitogenomics and phylogenetics of the family Carangidae with special emphasis on the mitogenome of the Indian Scad Decapterus russelli

Carangids are abundant and commercially important marine fish that contribute to a significant portion of the fisheries in many parts of the world. In the present study, we characterized the complete mitogenome of the Indian scad, Decapterus russelli and performed a comprehensive comparative mitogenomic analysis of the family Carangidae. The comparative mitogenomics provided valuable insights into the structure, variability, and features of the coding and non-coding regions that evolved across species over millions of years. The structural features of tRNAs revealed changes in the frequency of mismatched and wobble base pairs, which is reflected in the base composition of H and L strands. The highly conserved sequence motif of the mTERF binding site in carangids over the ~ 400 MYA of their divergence demonstrated the functional importance of these sites. The control region of carangids was characterized by the presence of discontinuous repeat units with a high rate of sequence divergence in the form of base substitutions, insertions, and deletions. The maintenance of secondary structures in the control region independent of the rapid evolution of primary structure suggested the effect of selective constraints on their maintenance. Maximum likelihood (ML) and Bayesian inference (BI) phylogeny revealed a similar topology consistent with previous taxonomic studies. The extant carangids diverged through the evolutionary events experienced during the Cretaceous, Paleogene, and Neogene periods.

Human mitochondrial DNA diversity is compatible with the multiregional continuity theory of the origin of Homo sapiens

Confidence intervals for estimates of human mtDNA sequence diversity, chimpanzee-human mtDNA sequence divergence, and the time of splitting of the pongid-hominid lineages are presented. Consistent with all the data used in estimating the coalescence time for human mitochondrial lineages to a common ancestral mitochondrion is a range of dates from less than 79,000 years ago to more than 1,139,000 years ago. Consequently, the hypothesis that a migration of modern humans (Homo sapiens) out of Africa in the range of 140,000 to 280,000 years ago resulted in the complete replacement, without genetic interchange, of earlier Eurasian hominid populations (Homo erectus) is but one of several possible interpretations of the mtDNA data. The data are also compatible with the hypothesis, suggested earlier and supported by fossil evidence, of a single, more ancient expansion of the range of Homo erectus from Africa, followed by a gradual transition to Homo sapiens in Europe, Asia, and Africa.

Mitochondrial genomes of African pangolins and insights into evolutionary patterns and phylogeny of the family Manidae

BACKGROUND

This study used next generation sequencing to generate the mitogenomes of four African pangolin species; Temminck's ground pangolin (Smutsia temminckii), giant ground pangolin (S. gigantea), white-bellied pangolin (Phataginus tricuspis) and black-bellied pangolin (P. tetradactyla).

RESULTS

The results indicate that the mitogenomes of the African pangolins are 16,558 bp for S. temminckii, 16,540 bp for S. gigantea, 16,649 bp for P. tetradactyla and 16,565 bp for P. tricuspis. Phylogenetic comparisons of the African pangolins indicated two lineages with high posterior probabilities providing evidence to support the classification of two genera; Smutsia and Phataginus. The total GC content between African pangolins was observed to be similar between species (36.5% - 37.3%). The most frequent codon was found to be A or C at the 3rd codon position. Significant variations in GC-content and codon usage were observed for several regions between African and Asian pangolin species which may be attributed to mutation pressure and/or natural selection. Lastly, a total of two insertions of 80 bp and 28 bp in size respectively was observed in the control region of the black-bellied pangolin which were absent in the other African pangolin species.

CONCLUSIONS

The current study presents reference mitogenomes of all four African pangolin species and thus expands on the current set of reference genomes available for six of the eight extant pangolin species globally and represents the first phylogenetic analysis with six pangolin species using full mitochondrial genomes. Knowledge of full mitochondrial DNA genomes will assist in providing a better understanding on the evolution of pangolins which will be essential for conservation genetic studies.

Novel duplication pattern of the mitochondrial control region in Cantor's Giant softshell turtle Pelochelys cantorii

Cantor's Giant Softshell Turtle, Pelochelys cantorii has become one of the most critically endangered species in the world. When comparative analyses of the P. cantorii complete mitochondrial genome sequences were conducted, we discovered a duplication of a segment of the control region in the mitochondrial genome of P. cantorii. The duplication is characterized by two copies of conserved sequence box 2 (CSB2) and CSB3 in a single control region. In contrast to previous reports of duplications involving the control regions of other animals, this particular pattern of duplications appears to be unique to P. cantorii. Copies of the CSB2 and CSB3 show many of the conserved sequence features typically found in mitochondrial control regions, and rare differences were found between the paralogous copies. Using the primer design principle of simple sequence repeats (SSR) and the reference sequence of the duplicated CSBs, specific primers were designed to amplify the duplicated CSBs. These primers were validated among different individuals and populations of P. cantorii. This unique duplication structure suggests the two copies of the CSB2 and CSB3 may have arisen through occasional tandem duplication and subsequent concerted evolution.

The mitochondrial genome of the Saunders's gull Chroicocephalus saundersi (Charadriiformes: Laridae) and a higher phylogeny of shorebirds (Charadriiformes)

The complete mitogenome of Chroicocephalus saundersi was characterized and compared with the 6 published Charadriiformes mitogenomes. The mitogenome of C. saundersi is a closed circular molecule 16,739 bp in size, and contains 37 genes and a control region. The AT and GC skews are positive and negative, respectively, and in agreement with those of the other Charadriiformes mitogenomes. The mitogenome of C. saundersi contains 3 start codons (ATG, GTG, and ATT), 4 stop codons (TAA, TAG, AGG, and AGA), and an incomplete stop codon (T-) in 13 PCGs. A codon usage analysis of all available Charadriiformes mitogenomes showed that the ATG (78%) and TAA (50.5%) were the most common start codon and stop codon, respectively. An unusual start codon, ATT, is commonly found in the ND3s of Charadriiformes mitogenomes, whereas the more common start codons, ATC and ATA, are rarely found. In all the Laridae species, one extra cytosine was inserted at position 174 in ND3. The control region of C. saundersi is 1180-bp long, with a nucleotide composition of 30.2% A, 28.6% T, 27.3% C, and 14.0% G. Variable numbers of tandem repeats (VNTRs) with nine copies of the 10 bp repeat sequence (AACAACAAAC) are found within the CSB domain of the control region. The ML/BI analyses, based on the amino acids of the 13 mitochondrial PCGs, strongly support the monophyly of the order Charadriiformes, with the suborder Lari considered sister to the Scolopaci, which is in turn a sister group to the suborder Charadrii.

Molecular Characterization of Sudanese and Southern Sudanese Chicken Breeds Using mtDNA D-Loop

The objective of this study was to assess the genetic relationships and diversity and to estimate the amount of gene flow among the five chicken populations from Sudan and South Sudan and commercial strain of egg line White Leghorn chickens. The chicken populations were genotyped using mtDNA D-loop as a molecular marker. PCR product of the mtDNA D-loop segment was 600 bp and 14 haplotypes were identified. The neighbor-joining phylogenetic tree indicated that the indigenous Sudanese chickens can be grouped into two clades, IV and IIIa only. Median joining networks analysis showed that haplotype LBB49 has the highest frequency. The hierarchal analysis of molecular variance (AMOVA) showed that genetic variation within the population was 88.6% and the differentiation among the population was 11.4%. When the populations was redefined into two geographical zones, rich and poor Savanna, the results were fractioned into three genetic variations: between individuals within population 95.5%, between populations within the group 0.75%, and genetic variation between groups 3.75%. The pair wise F st showed high genetic difference between Betwil populations and the rest with F st ranging from 0.1492 to 0.2447. We found that there is large number of gene exchanges within the Sudanese indigenous chicken (Nm = 4.622).

The complete mitochondrial genome of Aspiorhynchus laticeps and its phylogenetic analysis

The complete nucleotide sequence of the mitochondrial genome (mitogenome) of Aspiorhynchus laticeps was determined. The length of the complete mitochondrial DNA sequence of A. laticeps is 16,591 bp, which consists of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a non-coding region 'D-loop'. Except for the D-loop, another non-coding region named replication origin of L-strand (OL) region was also found. According to the phylogenetic analysis, A. laticeps has a closer relationship with Schizothorax.

Mitochondrial genome sequence of Sibynophis chinensis (Squamata, Colubridae)

The complete mitochondrial (mt) genome of the Chinese many-toothed snake, Sibynophis chinensis, was sequenced and found to be 17,163 bp in length. The arrangement of 13 protein-coding genes, tRNAs and rRNAs was identical to that of other common snake mt genomes. The mt protein-coding genes of S. chinensis utilized ATA, ATG, ATA and GTG as initiation codons and AGA, AGG, TAA, TAG and T as termination codons. Among three tRNA clusters (LQM, WANCY and HSL), LQM was found instead of IQM, which is common in other vertebrates. We also identified two control regions that contained several conserved elements known as conserved sequence blocks and termination-associated sequences related to mt replication and transcription.

U-turn DNA bending by human mitochondrial transcription factor A

Transcription factor A (TFAM) is involved in the transcription regulation, maintenance and compaction of the mitochondrial genome. Recent structural data on TFAM showed its mode of operation and clarified previous biochemical and genetic results. In solution, TFAM is highly dynamic. According to crystal structures of its complex with the cognate light-strand promoter (LSP) binding sequence, it intertwines and dramatically bends DNA, thereby allowing interactions with the transcription initiation machinery. Recent studies have shown TFAM sliding on non-specific DNA, which induces compaction by increasing DNA flexibility. Finally, the structural localization of disease-related TFAM mutations suggests functional impairment at the molecular level.

The mitochondrial genome of Apodemus peninsulae (Rodentia, Muridae)

The complete mitochondrial (mt) genome of the Korean field mouse Apodemus peninsulae was sequenced and found to be 16,266 bp in length. The mt protein-coding genes of A. peninsulae had ATG, GTG, ATC, and ATA as initiation codons and TAA, TAG, TA, and T as termination codons. Two forms each of trnL and trnS and the three tRNA clusters, IQM, WANCY, and HSL were identified, as in the typical Rodentia mt genome. Among tRNAs, abnormal cloverleaf structure of trnS((AGY)) was identified in DHU arm. The l-strand replication origin has the potential to form a stable stem-loop structure and control region has several conserved sequence elements.

Mitochondrial genome phylogeny among Asiatic black bear Ursus thibetanus subspecies and comprehensive analysis of their control regions

The complete mitochondrial genome (16,824 bp) of an Asiatic black bear Ursus thibetanus ussuricus (Mammalia, Carnivora, Ursidae) was newly sequenced and characterized in detail. It is the second mitochondrial genome from this subspecies which has been completely sequenced. The two U. t. ussuricus individuals were compared with each other and then with individuals from the other four U. thibetanus subspecies and the other nine ursid species, focusing especially on the control regions in the 14 mitochondrial genomes. Within these control regions, tandem repeats of basically 10 bp (5'-ACGCACGTGT-3' or its derivatives) were found in Domain II. Plausible secondary structures of the repeat region were compared between the North and South Korean individuals of U. t. ussuricus. According to the maximum likelihood and Bayesian inference trees inferred from the nucleotide sequences of 13 protein-coding and two rRNA genes, the ursine members within the monophyletic ursid clade can be divided into at least three groups: A, B, and C. According to this analysis, U. thibetanus subspecies were found with Ursus americanus and Ursus malayanus within Group A, showing the following relationships with nodal bootstrap values above 91% and Bayesian posterior probabilities of 1.00: ([(U. t. thibetanus, U. t. formosanus), U. t. spp.], U. t. ussuricus), U. t. mupinensis. In addition, we present a hypothetical scenario of the evolution of the major repeat motifs in the control region.

Mosaic gene conversion after a tandem duplication of mtDNA sequence in Diomedeidae (albatrosses)

Although the tandem duplication of mitochondrial (mt) sequences, especially those of the control region (CR), has been detected in metazoan species, few studies have focused on the features of the duplicated sequence itself, such as the gene conversion rate, distribution patterns of the variation, and relative rates of evolution between the copies. To investigate the features of duplicated mt sequences, we partially sequenced the mt genome of 16 Phoebastria albatrosses belonging to three species (P. albatrus, P. nigripes, and P. immutabilis). More than 2,300 base pairs of tandemly-duplicated sequence were shared by all three species. The observed gene arrangement was shared in the three Phoebastria albatrosses and suggests that the duplication event occurred in the common ancestor of the three species. Most of the copies in each individual were identical or nearly identical, and were maintained through frequent gene conversions. By contrast, portions of CR domains I and III had different phylogenetic signals, suggesting that gene conversion had not occurred in those sections after the speciation of the three species. Several lines of data, including the heterogeneity of the rate of molecular evolution, nucleotide differences, and putative secondary structures, suggests that the two sequences in CR domain I are maintained through selection; however, additional studies into the mechanisms of gene conversion and mtDNA synthesis are required to confirm this hypothesis.

Comparative analysis of mitochondrial control region in polyploid hybrids of red crucian carp (Carassius auratus) x blunt snout bream (Megalobrama amblycephala)

The entire sequences of the mitochondrial (mt)DNA control region (CR) and portions of its flanking genes in the red crucian carp (RC) and blunt snout bream (BSB) as well as their polyploid hybrids (3nRB, 4nRB and 5nRB) were determined and subjected to a comparative analysis. The mtDNA-CRs of these five fish species ranged from 923 to 937 bp in length, they had the same flanking gene arrangement as other vertebrates and the pattern of nucleotide substitution bias was also similar to that in other vertebrates. Our data are consistent with the viewpoint of three domains [extended terminal associated sequence (ETAS domain), central conserved sequence block domain and conserved sequence block (CSB) domain] within the mtDNA-CR of mammals. On the basis our comparative analysis of the mtDNA-CRs of these five fish species, we were able to identify the consensus sequences of functional conserved units, including the ETAS, CSB-F, CSB-D, CSB-E, CSB1, CSB2 and CSB3 and putative promoter. The percentage of variable nucleotide positions (41.98%) in the central domain was lower than those in the ETAS and conserved domain (71.70 and 47.12%, respectively), suggesting that the central domain was the most conserved part of the mtDNA-CR. These results provide useful and important information for the further study of mtDNA-CR structure in fish. The sequence similarities of mtDNA-CR among the 3nRB, 4nRB, 5nRB hybrids and their respective female parents were higher than those among the 3nRB, 4nRB, 5nRB hybrids and their respective male parents, providing the direct evidence of stringent maternal inheritance of mtDNA-CR in the 3nRB, 4nRB and 5nRB hybrids.

The structure of the Australian and South American marsupial mitochondrial control region

BACKGROUND AND AIMS

The mitochondrial control region (CR) was studied across five marsupialian orders, in order to give a detailed overview of its features.

RESULTS

The CR is organised into three domains similar to the CR of placental mammals. However, the conservation of different features among the marsupial orders is in general more strict. In the first domain, two conserved blocks extended termination-associated sequences (ETAS 1 and ETAS 2) are present in all marsupial orders. In the third domain, the three conserved sequence blocks (CSB 1, CSB 2 and CSB 3) are present and complete, with CSB 1 being duplicated in four of five marsupial orders.

CONCLUSIONS

The nucleotide frequency and secondary structures of the repeats were typical for marsupial species. The repeats are generally AT-rich except in Dasyuridae and Paucituberculata, which show a significant increase in GC content.

Complete mitochondrial genome of Otis tarda (Gruiformes: Otididae) and phylogeny of Gruiformes inferred from mitochondrial DNA sequences

The complete nucleotide sequence of mitochondrial genome of the Great bustard (Otis tarda) was determined by using polymerase chain reaction (PCR) method. The genome is 16,849 bp in size, containing 13 protein-coding, 2 ribosomal and 22 transfer RNA genes. Sequences of the tRNA genes can be folded into canonical cloverleaf secondary structure except for tRNA-Cys and tRNA-Ser (AGY), which lose "DHU" arm. Sequence analysis showed that the O. tarda mitochondrial control region (mtCR) contained many elements in common with other avian mtCRs. A microsatellite repeat was found in the 3'-peripheral domain of the O. tarda mtCR. Based on the mitochondrial DNA sequences of 12S rRNA, 16S rRNA and tRNA-Val, a phylogenetic study of Gruiformes was performed. The result showed that Otididae was a sister group to "core Gruiformes" and Charadriiformes with strong support (97% posterior probability values) in Bayesian analysis. The taxonomic status of Rhynochetidae, Mesitornithidae, Pedionomidae and Turnicidae that traditionally belonged to Gruiformes was also discussed in this paper.

Evolution of the mitochondrial genome in mammals living at high altitude: new insights from a study of the tribe Caprini (Bovidae, Antilopinae)

Organisms living at high altitude are exposed to severe environmental stress associated with decreased oxygen pressure, cold temperatures, increased levels of UV radiation, steep slopes, and scarce food supplies, which may have imposed important selective constraints on the evolution of the mitochondrial genome. Within mammals, the tribe Caprini is of particular interest for studying the evolutionary effects of life at high altitude, as most species live in mountain regions, where they developed morphological and physiological adaptations for climbing. In this report, we analyzed the complete mitochondrial genome of 24 ruminants, including 20 species of Caprini. The phylogenetic analyses based on 16,117 nucleotides suggested the existence of a new large clade, here named subtribe Caprina, containing all genera, but Pantholops (Pantholopina), Capricornis, Naemorhedus, and Ovibos (Ovibovina). The alignment of the control region showed that all Caprini have between two and four tandem repeats of ~75 bp in the RS2 region, and that several of these copies emerged from recent and independent duplication events. We proposed therefore that the maintenance of at least two RS2 repeats in the control region of Caprini is positively selected, probably for producing a higher number of D-loop strands 3'-ending at different locations. The analyses of base composition at third-codon positions of protein-coding genes revealed that Caprini have the highest percentage of A nucleotide and the lowest percentage of G nucleotide, a pattern which suggests increased rates of cytosine deamination (C-->T transitions) on the H strand of mtDNA. Two nonexclusive hypotheses related to high-altitude life can explain such a mutational pattern: more severe oxidative stress (ROS) and higher metabolic rates. By comparing the relative rates of nonsynonymous and synonymous substitutions in protein-coding genes, we identified that Caprini have higher levels of adaptive variation in the ATPase complex. In addition, we detected several changes in mitochondrial genes that should be tested for their potential role in mountain adaptation.

Characterization of human control region sequences for Spanish individuals in a forensic mtDNA data set

Population data on the hypervariable regions of the mitochondrial DNA (mtDNA) genome are used to convey the relative rarity of mtDNA profiles obtained from evidence samples and of profiles used to identify missing persons. In this study, mtDNA profiles of Spanish individuals (n=312) were analyzed to describe haplogroup distributions and to determine relevant single nucleotide polymorphisms (SNPs) of those haplogroups. All nine common European haplogroups were observed in the sample, and these were divided into subgroups when possible. Haplogroup H was the most common haplogroup. The haplogroups U, J, T, and V were the next most frequent groups, each occurring at a frequency of 6.4% or greater. In addition, African and Asian sequences were present though rare in the samples. The data were compared with and found to be similar to other published data sets. There were 109 SNPs observed in the data set, including 10 positions not previously reported. The most variable sites are consistent with other studies.

Evaluation of Variation in Control Region Sequences for Hispanic Individuals in the SWGDAM mtDNA Data Set

The Scientific Working Group on DNA Analysis Methods (SWGDAM) Hispanic data set was analyzed to determine the diversity, phylogeny, and relevant single nucleotide polymorphisms (SNPs) that describe haplogroup patterns for Hispanic Americans (N=686), and to assess the degree of admixture regarding mitochondrial DNA (mtDNA). The largest component of admixture based on mtDNA analysis derives from the four major haplogroups previously observed in Native American ancestry, including A (29.3%), B (15.7%), C (20.6%), and D (4.8%). European (17.8%) and African (11.8%) haplogroups also were observed within this data set. Hispanic SWGDAM samples from the southwest, compared with other SWGDAM Hispanic samples, were observed to have a greater percent of Native American haplogroups present (79.9%), and fewer African American haplogroups (4.5%). A total of 234 SNPs were observed in the data set, including 36 newly reported variable positions. These SWGDAM Hispanic data set SNPs ranged from having 1 to 31 changes (Length=L) on the phylogenetic tree, with site 16519 being the most variable. On average, there were 3.9 character changes for each variable position on the tree. The most variable sites (with 13 or more changes each listed from fastest to slowest) observed were 16519 (L=31), 16189 (L=23), 152 (L=23), 16311 (L=19), 146 (L=17), 195 (L=17), 16093 (L=15), 16362 (L=14), 16129 (L=13), 150 (L=13), and 153 (L=13). These sites are consistent with other reports on highly variable positions. A total of 27 SNPs were chosen to identify all clusters containing 1% (N=7) or more individuals in the SWGDAM Hispanic data set. The descriptive analyses revealed that the SWGDAM Hispanic data set is similar to published Native American and Hispanic data sets.

Nuclear insertions help and hinder inference of the evolutionary history of gorilla mtDNA

Numts are fragments of mitochondrial DNA (mtDNA) that have been translocated to the nucleus, where they can persist while their mitochondrial counterparts continue to rapidly evolve. Thus, numts represent 'molecular fossils' useful for comparison with mitochondrial variation, and are particularly suited for studies of the fast-evolving hypervariable segment of the mitochondrial control region (HV1). Here we used information from numts found in western gorillas (Gorilla gorilla) and eastern gorillas (Gorilla beringei) to estimate that these two species diverged about 1.3 million years ago (Ma), an estimate similar to recent calculations for the divergence of chimpanzee and bonobo. We also describe the sequence of a gorilla numt still possessing a segment lost from all contemporary gorilla mtDNAs. In contrast to that sequence, many numts of the HV1 are highly similar to authentic mitochondrial organellar sequences, making it difficult to determine whether purported mitochondrial sequences truly derive from that genome. We used all available organellar HV1 and corresponding numt sequences from gorillas in a phylogenetic analysis aimed at distinguishing these two types of sequences. Numts were found in several clades in the tree. This, in combination with the fact that only a limited amount of the extant variation in gorillas has been sampled, suggests that categorization of new sequences by the indirect means of phylogenetic comparison would be prone to uncertainty. We conclude that for taxa such as gorillas that contain numerous numts, direct approaches to the authentication of HV1 sequences, such as amplification strategies relying upon the circularity of the mtDNA molecule, remain necessary.

Characterization of human control region sequences of the African American SWGDAM forensic mtDNA data set

The scientific working group on DNA analysis Methods (SWGDAM) mitochondrial DNA (mtDNA) population data set is used to infer the relative rarity of control region mtDNA profiles obtained from evidence samples and of profiles used for identification of missing persons. In this study, the African American haplogroup patterns in the SWGDAM data were analyzed in a phylogenetic context to determine relevant single nucleotide polymorphisms (SNPs) and to describe haplogroup distributions for Africans observed in these data sets. Over 200 SNPs (n=217) were observed in the African American data set (n=1148). These SNPs ranged from having 1-39 changes in the phylogenetic tree, with sites 152 and 16519 being the most variable. On average there were 5.8 changes for a character on the tree. The most variable sites (with 19 or more changes each) observed included 16093, 16129, 16189, 16311, 16362, 16519, 146, 150, 152, 189, and 195. These rapidly changing sites are consistent with other published analyses. Only 34 SNPs are needed to identify all clusters containing 10 or more individuals in the African American data set. The results show that the African American SWGDAM mtDNA data set contains variation consistent with that described in continental African populations. Thirteen of the 18 haplogroups previously observed in African populations were observed and include: L1a, L1b, L1c, L2a, L2b, L2c, L3b, L3d, L3e1, L3e2, L3e3, L3e4 and L3f. Haplogroup L2a is the most commonly observed cluster (18.8%) in the African American data set. The next most common haplogroups in the African American data set include the clusters L1c (11.0%), L1b (9.1%), L3e2 (9.0%) and L3b (8.1%). Approximately 8% of the haplogroups observed within African Americans were common in European Caucasians or East Asians; these were H (n=32), J (n=4), K (n=5), T (n=2), U5 (n=6), U6 (n=9 also known from North Africa), A (n=12), B (n=7), C (n=4), and M (n=16), respectively. The European Caucasian and East Asian haplogroups are expected due to admixture between individuals with recent ancestry in Western Eurasia and sub-Saharan Africa. The genetic characterization of these relevant data sets is fully consistent with other published mtDNA genetic variation. The sequence diversity observed in this data set makes it a valuable tool for forensic applications.

Mutagenesis by transient misalignment in the human mitochondrial DNA control region

To study spontaneous base substitutions in human mitochondrial DNA (mtDNA), we reconstructed the mutation spectra of the hypervariable segments I and II (HVS I and II) using published data on polymorphisms from various human populations. Classification analysis revealed numerous mutation hotspots in HVS I and II mutation spectra. Statistical analysis suggested that strand dislocation mutagenesis, operating in monotonous runs of nucleotides, plays an important role in generating base substitutions in the mtDNA control region. The frequency of mutations compatible with the primer strand dislocation in the HVS I region was almost twice as high as that for template strand dislocation. Frequencies of mutations compatible with the primer and template strand dislocation models are almost equal in the HVS II region. Further analysis of strand dislocation models suggested that an excess of pyrimidine transitions in mutation spectra, reconstructed on the basis of the L-strand sequence, is caused by an excess of both L-strand pyrimidine transitions and H-strand purine transitions. In general, no significant bias toward parent H-strand-specific dislocation mutagenesis was found in the HVS I and II regions.

Control region sequences for East Asian individuals in the Scientific Working Group on DNA Analysis Methods forensic mtDNA data set

The Scientific Working Group on DNA Analysis Methods (SWGDAM) mitochondrial DNA (mtDNA) population data set is used to infer the relative rarity of mtDNA profiles obtained from evidence samples and of profiles used to identify missing persons. In this study, the East Asian haplogroup patterns in the SWGDAM data sets were analyzed in a phylogenetic context to determine relevant single nucleotide polymorphisms (SNPs) and to describe haplogroup distributions for Asians (n = 753; with a breakdown of individuals from China n = 356, Korea n = 182, Japan n = 163, and Thailand n = 52). We focus on the patterns observed in the SWGDAM Chinese data set and refer to interesting differences in the smaller subgroup data sets for the other East Asian populations (Japanese, Korean, and Thai). A total of 218 SNPs were observed in the data set, including 37 observed positions not previously reported. In the largest of the East Asian SWGDAM data sets (Chinese), these SNPs ranged from having 1 to 29 changes in the phylogenetic tree, with site 16519 being the most variable. On average there were 4.5 changes for a character on the tree. The most variable sites (with 14 or more changes each listed from fastest to slowest) observed were 16519 (L = 29), 16311 (L = 27), 152 (L = 24), 146 (L = 21), 16172 (L = 17), 16189 (L = 17), 195 (L = 16), 16362 (L = 15), 16093 (L = 14), 16129 (L = 14) and 150 (L = 14). These rapidly changing sites are consistent with other published analyses. Only 28 SNPs are needed to identify all clusters containing 1% (n = 7) or more individuals in the East Asian data set. All 36 haplogroups previously observed in East Asian populations were also seen in the SWGDAM data sets and include: A, B, B4, B4a, B4b, B5a, B5b, C, D, D4, D4a, D4b, D5, D5a, F, F1, F1a, F1b, F1c, F2a, G2, G2a, M, M7a1, M7b, M7b1, M7b2, M7c, M8a, M9, M10, N9a, R, R9a, Y, and Z. Haplogroups A, B4a, D4, and F1a were the most commonly observed clusters in the Chinese data set (the largest of the data sets) with each of these occurring in more than 6% of the samples in the data set. The next most common haplogroups in the Chinese data set include the clusters C, M7b1, and N9a with each observed at frequencies greater than or equal to 4%. European Caucasian, and African haplogroups were rarely observed within the East Asian data sets. The various analyses revealed that the data set was similar to published East Asian data sets such as those from Han Chinese.

Unreliable mtDNA data due to nuclear insertions: a cautionary tale from analysis of humans and other great apes

Analysis of mitochondrial DNA sequence variation has been used extensively to study the evolutionary relationships of individuals and populations, both within and across species. So ubiquitous and easily acquired are mtDNA data that it has been suggested that such data could serve as a taxonomic 'barcode' for an objective species classification scheme. However, there are technical pitfalls associated with the acquisition of mtDNA data. One problem is the presence of translocated pieces of mtDNA in the nuclear genome of many taxa that may be mistaken for authentic organellar mtDNA. We assessed the extent to which such 'numt' sequences may pose an overlooked problem in analyses of mtDNA from humans and apes. Using long-range polymerase chain reaction (PCR), we generated necessarily authentic mtDNA sequences for comparison with sequences obtained using typical methods for a segment of the mtDNA control region in humans, chimpanzees, bonobos, gorillas and orangutans. Results revealed that gorillas are notable for having such a variety of numt sequences bearing high similarity to authentic mtDNA that any analysis of mtDNA using standard approaches is rendered impossible. Studies on humans, chimpanzees, bonobos or orangutans are apparently less problematic. One implication is that explicit measures need to be taken to authenticate mtDNA sequences in newly studied taxa or when any irregularities arise. Furthermore, some taxa may not be amenable to analysis of mtDNA variation at all.

Survey of gorillas (Gorilla gorilla gorilla) and chimpanzees (Pan troglodytes troglodytes) in Southwestern Cameroon

A study on the distribution and population density of the western lowland gorilla (Gorilla g. gorilla) and the central chimpanzee (Pan t. troglodytes) was undertaken between December 1997 and August 2000 in the Campo and Ma'an Forests in southwestern Cameroon. The aim of this survey was to estimate the densities of the apes in different parts of the area, to assess the importance of the region for the conservation of these endangered species and to determine the influence of human activities such as logging and hunting. The survey was based on night nest counts on a total of 665.5 km of line transects. The overall density in the Campo Forest was estimated at 0.2 gorillas/km(2) and at 0.63-0.78 chimpanzees/km(2). The overall density of chimpanzees in the Ma'an Forest was estimated at 0.8-1 individuals/km(2). Gorilla density in this area was too low to allow an estimation. The highest gorilla nest density was found in secondary forest. The gorilla density in unlogged forest was significantly lower. Chimpanzees showed a clear preference for less disturbed areas. In unlogged forest, old secondary forests (logging more than 23 years ago) and areas of recent logging with large remaining patches of primary forest, significantly higher densities were calculated than inside the more heavily exploited logging concession. In areas with both logging and high hunting pressure both species were rare or even absent. The Campo Ma'an area is considered a very important area for the conservation of gorillas and chimpanzees. Conservation measures are urgently required to reduce the impact of logging and hunting. The creation of the Campo Ma'an National Park in January 2000 is an important measure to preserve the unique biodiversity in this so far hardly protected area.

FORENSICS ANDMITOCHONDRIALDNA: Applications, Debates, and Foundations

Debate on the validity and reliability of scientific methods often arises in the courtroom. When the government (i.e., the prosecution) is the proponent of evidence, the defense is obliged to challenge its admissibility. Regardless, those who seek to use DNA typing methodologies to analyze forensic biological evidence have a responsibility to understand the technology and its applications so a proper foundation(s) for its use can be laid. Mitochondrial DNA (mtDNA), an extranuclear genome, has certain features that make it desirable for forensics, namely, high copy number, lack of recombination, and matrilineal inheritance. mtDNA typing has become routine in forensic biology and is used to analyze old bones, teeth, hair shafts, and other biological samples where nuclear DNA content is low. To evaluate results obtained by sequencing the two hypervariable regions of the control region of the human mtDNA genome, one must consider the genetically related issues of nomenclature, reference population databases, heteroplasmy, paternal leakage, recombination, and, of course, interpretation of results. We describe the approaches, the impact some issues may have on interpretation of mtDNA analyses, and some issues raised in the courtroom.

Characterization of the red knot (Calidris canutus) mitochondrial control region

We sequenced the complete mitochondrial control regions of 11 red knots (Calidris canutus). The control region is 1168 bp in length and is flanked by tRNA glutamate (glu) and the gene ND6 at its 5' end and tRNA phenylalanine (phe) and the gene 12S on its 3' end. The sequence possesses conserved sequence blocks F, E, D, C, CSB-1, and the bird similarity box (BSB), as expected for a mitochondrial copy. Flanking tRNA regions show correct secondary structure, and a relative rate test indicated no significant difference between substitution rates in the sequence we obtained versus the known mitochondrial sequence of turnstones (Charadriiformes: Scolopacidae). These characteristics indicate that the sequence is mitochondrial in origin. To confirm this, we sequenced the control region of a single individual using both purified mitochondrial DNA and genomic DNA. The sequences were identical using both methods. The sequence and methods presented in this paper may now serve as a reference for future studies using knot and other avian control regions. Furthermore, the discovery of five variable sites in 11 knots towards the 3' end of the control region, and the variability of this region in contrast to the more conserved central domain in the alignment between knots and other Charadriiformes, highlights the importance of this area as a source of variation for future studies in knots and other birds.

The crocodilian mitochondrial control region: general structure, conserved sequences, and evolutionary implications

We present the first comprehensive analysis of the crocodilian control region. We have analyzed sequences from all three families of Crocodylia (Crocodylidae, Gavialidae, Alligatoridae), incorporating all genera except Paleosuchus and Melanosuchus. Within the control region of other vertebrates, several sequence motifs and their order appear to be conserved. Herein, we compare aligned crocodilian D-loop sequences to homologous sequences from other vertebrates ranging from fish to birds. Among other findings, we have discovered that while domain I tends to be shorter than the same region in mammals and birds, it contains sequences similar in structure to both the goose-hairpin and termination associated sequences (TAS). Domain II is highly conservative with regard to size among the taxa examined and contains several of the conserved sequence boxes characterized in other vertebrates. Domain III contains several interesting sequence motifs including tandemly repeated sequences, a long poly-A region in the Crocodylidae, and possible bidirection promoter sequences.

Mitochondrial DNA sequence evolution and phylogeny of the Atlantic Alcidae, including the extinct great auk (Pinguinus impennis)

The Atlantic auk assemblage includes four extant species, razorbill (Alca torda), dovekie (Alle alle), common murre (Uria aalge), and thick-billed murre (U. lomvia), and one recently extinct species, the flightless great auk (Pinguinus impennis). To determine the phylogenetic relationships among the species, a contiguous 4.2-kb region of the mitochondrial genome from the extant species was amplified using PCR. This region included one ribosomal RNA gene, four transfer RNA genes, two protein-coding genes, the control region, and intergenic spacers. Sets of PCR primers for amplifying the same region from great auk were designed from sequences of the extant species. The authenticity of the great auk sequence was ascertained by alternative amplifications, cloning, and separate analyses in an independent laboratory. Phylogenetic analyses of the entire assemblage, made possible by the great auk sequence, fully resolved the phylogenetic relationships and split it into two primary lineages, Uria versus Alle, Alca, and Pinguinus. A sister group relationship was identified between Alca and Pinguinus to the exclusion of ALLE: Phylogenetically, the flightless great auk originated late relative to other divergences within the assemblage. This suggests that three highly divergent species in terms of adaptive specializations, Alca, Alle, and Pinguinus, evolved from a single lineage in the Atlantic Ocean, in a process similar to the initial adaptive radiation of alcids in the Pacific Ocean.

Structure and evolution of the avian mitochondrial control region

The structural and evolutionary characteristics of the mitochondrial control region were studied by using control region sequences of 68 avian species. The distribution of the variable nucleotide positions within the control region was found to be genus specific and not dependant on the level of divergence, as suggested before. Saturation was shown to occur at the level of divergence of 10% in pairwise comparisons of the control region sequences, as has also been reported for the third codon positions in ND2 and cytochrome b genes of mtDNA. The ratio of control region vs cytochrome b divergence in pairwise comparisons of the sequences was shown to vary from 0.13 to 21.65, indicating that the control region is not always the most variable region of the mtDNA, but also that there are differences in the rate of divergence among the lineages. Only two of the conserved sequence blocks localized earlier for other species, D box and CSB-1, were found to show a considerable amount of sequence conservation across the avian and mammalian sequences. Additionally, a novel avian-specific sequence block was found.

Mitochondrial DNA variability of West New Guinea populations

This paper reports human mitochondrial DNA variability in West New Guinea (the least known, western side of the island of New Guinea), not yet described from a molecular perspective. The study was carried out on 202 subjects from 12 ethnic groups, belonging to six different Papuan language families, representative of both mountain and coastal plain areas. Mitochondrial DNA hypervariable region 1 (HVS 1) and the presence of the 9-bp deletion (intergenic region COII-tRNA(Lys)) were investigated. HVS 1 sequencing identified 73 polymorphic sites defining 89 haplotypes; the 9-bp deletion, which is considered a marker of Austronesian migration in the Pacific, was found to be absent in the whole West New Guinea study sample. Statistical analysis applied to the resulting haplotypes reveal high heterogeneity and an intersecting distribution of genetic variability in these populations, despite their cultural and geographic diversity. The results of subsequent phylogenetic approaches subdivide mtDNA diversity in West New Guinea into three main clusters (groups I-III), defined by sets of polymorphisms which are also shared by some individuals from Papua New Guinea. Comparisons with worldwide HVS 1 sequences stored in the MitBASE database show the absence of these patterns outside Oceania and a few Indonesian subjects, who also lack the 9-bp deletion. This finding, which is consistent with the effects of genetic drift and prolonged isolation of West New Guinea populations, lead us to regard these patterns as New Guinea population markers, which may harbor the genetic memory of the earliest human migrations to the island.

The complete sequence of the zebrafish (Danio rerio) mitochondrial genome and evolutionary patterns in vertebrate mitochondrial DNA

We describe the complete sequence of the 16,596-nucleotide mitochondrial genome of the zebrafish (Danio rerio); contained are 13 protein genes, 22 tRNAs, 2 rRNAs, and a noncoding control region. Codon usage in protein genes is generally biased toward the available tRNA species but also reflects strand-specific nucleotide frequencies. For 19 of the 20 amino acids, the most frequently used codon ends in either A or C, with A preferred over C for fourfold degenerate codons (the lone exception was AUG: methionine). We show that rates of sequence evolution vary nearly as much within vertebrate classes as between them, yet nucleotide and amino acid composition show directional evolutionary trends, including marked differences between mammals and all other taxa. Birds showed similar compositional characteristics to the other nonmammalian taxa, indicating that the evolutionary trend in mammals is not solely due to metabolic rate and thermoregulatory factors. Complete mitochondrial genomes provide a large character base for phylogenetic analysis and may provide for robust estimates of phylogeny. Phylogenetic analysis of zebrafish and 35 other taxa based on all protein-coding genes produced trees largely, but not completely, consistent with conventional views of vertebrate evolution. It appears that even with such a large number of nucleotide characters (11,592), limited taxon sampling can lead to problems associated with extensive evolution on long phyletic branches.

Models of sequence evolution for DNA sequences containing gaps

Most evolutionary tree estimation methods for DNA sequences ignore or inefficiently use the phylogenetic information contained within shared patterns of gaps. This is largely due to the computational difficulties in implementing models for insertions and deletions. A simple way to incorporate this information is to treat a gap as a fifth character (with the four nucleotides being the other four) and to incorporate it within a Markov model of nucleotide substitution. This idea has been dismissed in the past, since it treats a multiple-site insertion or deletion as a sequence of independent events rather than a single event. While this is true, we have found that under many circumstances it is better to incorporate gap information inadequately than to ignore it, at least for topology estimation. We propose an extension to a class of nucleotide substitution models to incorporate the gap character and show that, for data sets (both real and simulated) with short and medium gaps, these models do lead to effective use of the information contained within insertions and deletions. We also implement an ad hoc method in which the likelihood at columns containing multiple-site gaps is downweighted in order to avoid giving them undue influence. The precision of the estimated tree, assessed using Markov chain Monte Carlo techniques to find the posterior distribution over tree space, improves under these five-state models compared with standard methods which effectively ignore gaps.

On the estimation of the rate of nucleotide substitution for the control region of human mitochondrial DNA

To apply molecular clock for studying human evolution, the pattern of nucleotide substitution for the control region of human mtDNA was analyzed in detail. It is well known that the rate of nucleotide substitution for the control region is much higher than that for any other part of mtDNA. In this study, the higher substitution rate was attributed to the higher rate of transition-type substitution between pyrimidines within the D-loop part, whereas the rates of other types of substitution were essentially the same over the entire mtDNA molecule. Even within the control region, the rate and pattern of nucleotide substitution were different between the D-loop part and the rest. The rate and pattern for the non-D-loop part were very similar to those for fourfold-degenerate sites in the protein-coding region. In contrast, the D-loop and non-D-loop parts showed similarities in the base composition, whereas the base composition of fourfold-degenerate sites slightly different from that of the both parts of the control region. It is concluded, therefore, that the nucleotide frequencies of the control region should be used to estimate the number of substitutions (d) between the control region sequences. However, a method to verify the accuracy of the estimation of d by means of the transition/transversion (s/v) ratio was theoretically studied. It was suggested that the s/v ratio becomes constant over a wide range of d values only when the estimation of d is unbiased. On the basis of this result, the estimates of d previously obtained between human sequences were evaluated.

Identification of a novel human mitochondrial D-loop RNA species which exhibits upregulated expression following cellular immortalization

We report the identification and characterization of a novel human mitochondrial RNA species approximately 0.47 kb long that is transcribed from the mtDNA L-strand and is derived from the D-loop. Its expression increases when human cells become immortal, a key event in tumorigenesis. The RNA is therefore designated IDL (Immortalization-associated D-Loop). Sequence and hybrid cell analyses suggest that the increased level of IDL RNA in immortal cells is due to a recessive change, possibly in the activity of a trans-acting factor that controls IDL RNA expression.

Peopling of Sahul: mtDNA variation in aboriginal Australian and Papua New Guinean populations

We examined genetic affinities of Aboriginal Australian and New Guinean populations by using nucleotide variation in the two hypervariable segments of the mtDNA control region (CR). A total of 318 individuals from highland Papua New Guinea (PNG), coastal PNG, and Aboriginal Australian populations were typed with a panel of 29 sequence-specific oligonucleotide (SSO) probes. The SSO-probe panel included five new probes that were used to type an additional 1,037 individuals from several Asian populations. The SSO-type data guided the selection of 78 individuals from Australia and east Indonesia for CR sequencing. A gene tree of these CR sequences, combined with published sequences from worldwide populations, contains two previously identified highland PNG clusters that do not include any Aboriginal Australians; the highland PNG clusters have coalescent time estimates of approximately 80,000 and 122,000 years ago, suggesting ancient isolation and genetic drift. SSO-type data indicate that 84% of the sample of PNG highlander mtDNA belong to these two clusters. In contrast, the Aboriginal Australian sequences are intermingled throughout the tree and cluster with sequences from multiple populations. Phylogenetic and multidimensional-scaling analyses of CR sequences and SSO types split PNG highland and Aboriginal Australian populations and link Aboriginal Australian populations with populations from the subcontinent of India. These mtDNA results do not support a close relationship between Aboriginal Australian and PNG populations but instead suggest multiple migrations in the peopling of Sahul.