Mitochondrial DNA Length Variation and Heteroplasmy in Populations of White Sturgeon (Acipenser Transmontanus)

Southern blot analysis was used to quantify the extent of mtDNA D-loop length variation in two populations of white sturgeon, Acipenser transmontanus. Over 42% of individuals were heteroplasmic for up to six different mtDNA length variants attributable to varying copy numbers of an 82-bp repeat sequence. Chi-square analyses revealed that the frequencies of length genotypes and the incidence of heteroplasmy were significantly different between Fraser and Columbia River sturgeon populations but not between restriction site haplotypes. Heteroplasmic fish have, on average, higher copy number than homoplasmic fish. Forty-five of 101 homoplasmic individuals carry only a single copy of the repeat, while none of the 73 heteroplasmic fish has the single repeat as the predominant variant. On the basis of differences in frequency distributions of copy number within and between fish, we suggest that (1) heteroplasmy is maintained by high recurrent mutation of multiple copy genomes, favoring increased copy number and (2) the mutation pressure toward higher copy number heteroplasmy is partially offset by selection to reduced genome size and segregation to the homoplasmic condition.     

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.     

Pattern of nucleotide substitution and rate heterogeneity in the hypervariable regions I and II of human mtDNA.

This study provides a comprehensive survey of the complex pattern of nucleotide substitution in the control region of human mtDNA, which is of central importance to the studies of human evolution. A total of 1229 different hypervariable region I (HVRI) and 385 different hypervariable region II (HVRII) sequences were analyzed using a complex substitution model. Moreover, we suggest a new method to assign relative rates to each site in the sequence. Estimates are based on maximum-likelihood methods applied to randomly selected subsets of sequences. Our results indicate that the rate of substitution in HVRI is approximately twice as high as in HVRII and that this difference is mainly due to a higher frequency of pyrimidine transitions in HVRI. However, rate heterogeneity is more pronounced in HVRII.     

Application of mtDNA d-loop region for the study of Russian sturgeon population structure from Iranian coastline of the Caspian Sea

The population structure of the Russian sturgeon ( Acipenser gueldenstaedtii ) was investigated using PCR amplification of the mtDNA D-loop region. Seven composite haplotypes were detected and average nucleotide and haplotype diversity over all samples were found to be 0.05 ± 0.00 and 0.75 ± 0.00 (mean ± SE) respectively. Restriction digest of the mtDNA D-loop region detected two genotypes A and B with a relative high frequencies of 0.5 and 0.36 respectively. These two genotypes (A and B) can be considered as potential genetic markers for different biological groups, stocks or seasonal races of Russian sturgeon. A maximum of 10.8% nucleotide diversity was found between haplotypes AAAAB and BBBBF. 100% heteroplasmy was found in Russian sturgeon and restriction digest of the mtDNA D-loop region also exhibited site heteroplasmy.     

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

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

Genetic characterization of horse bone excavated from the Kwakji archaeological site,Jeju, Korea

We determined the nucleotide sequences of the hypervariable D-loop region of mitochondrial DNA (mtDNA) from horse bone (humerus, A.D. 700 to A.D. 800) that was excavated from the Kwakji archaeological site, Jeju, Korea. We compared them with ones from extant horses. We designed three pairs of oligonucleotide primers from the tRNA-Thr and tRNA-Phe gene regions of mtDNA that are highly conserved among many other animal species. We cloned 232, 336, and 644 bp from the horse bone in order to determine the mtDNA D-loop sequence. The sequence was 1,124 bp long; the middle contained 19 tandem repeats of an 8-bp sequence (TGTGCACC) that is specific to equines. The mtDNA D-loop region contained each base (total number, percentage of total) as follows: A (317, 28.20%), C (336, 29.89%), G (169, 15.04%), and T (302, 26.87%). This sequence, like those of other horse populations, was AT rich. Sequence divergence was the lowest (1.71%) between the ancient horse bone and that of the Thoroughbred horse 1. The neighbor-joining and strict consensus tree of three of the most parsimonious trees also suggested that the ancient bone was considerably unrelated to native Jeju horses. The molecular phylogenetic characteristics of the horse bone that was excavated from the Kwakji archaeological site (Jeju, Korea) showed that some horse breeds may have existed on Jeju Island, Korea before Mongolian horses were introduced. The horse bone that was excavated from the Kwakji archaeological site may aid future research on the origin and ancestry of native Jeju horses.     

Phylogeography of the African buffalo based on mitochondrial and Y-chromosomal loci: Pleistocene origin and population expansion of the Cape buffalo subspecies

Population genetics and phylogeography of the African buffalo (Syncerus caffer) are inferred from genetic diversity at mitochondrial D-loop hypervariable region I sequences and a Y-chromosomal microsatellite. Three buffalo subspecies from different parts of Africa are included. Nucleotide diversity of the subspecies Cape buffalo at hypervariable region I is high, with little differentiation between populations. A mutation rate of 13-18% substitutions/million years is estimated for hypervariable region I. The nucleotide diversity indicates an estimated female effective population size of 17 000-32 000 individuals. Both mitochondrial and Y-chromosomal diversity are considerably higher in buffalo from central and southwestern Africa than in Cape buffalo, for which several explanations are hypothesized. There are several indications that there was a late middle to late Pleistocene population expansion in Cape buffalo. This also seems to be the period in which Cape buffalo evolved as a separate subspecies, according to the net sequence divergence with the other subspecies. These two observations are in agreement with the hypothesis of a rapid evolution of Cape buffalo based on fossil data. Additionally, there appears to have been a population expansion from eastern to southern Africa, which may be related to vegetation changes. However, as alternative explanations are also possible, further analyses with autosomal loci are needed.     

Sequence and properties of the human KB cell and mouse L cell D-loop regions of mitochondrial DNA.

The sequences of the displacement-loop (D-loop) regions of mitochondrial DNA (mtDNA) from mouse L cells and human KB cells have been determined and provide physical maps to aid in the identification of sequences involved in the regulation of replication and expression of mammalian mtDNA. Both D-loop regions are bounded by the genes for tRNAPhe and tRNAPro. This region contains the most highly divergent sequences in mtDNA with the exceptions of three small conserved sequence blocks near the 5' ends of D-loop strands, a 225 nucleotide conserved sequence block in the center of the D-loop strand template region, and a short sequence associated with the 3' ends of D-loop strands. A sequence similar to that associated with the 3' termini of D-loop strands overlaps one of the conserved sequence blocks near the 5' ends of D-loop strands. The large, central conserved sequence probably does not code for a protein since no open reading frames are discretely conserved. Numerous symmetric sequences and potential secondary structures exist in these sequences, but none appear to be clearly conserved between species.     

High frequency of mitochondrial DNA D-loop mutations in Ewing’s sarcoma

Somatic mutations and polymorphisms in the noncoding displacement (D)-loop of mitochondrial DNA (mtDNA) are present in a variety of human cancers. To investigate whether Ewing’s sarcoma (EWS) harbors genetic alterations within the D-loop region and their potential association with EWS carcinogenesis, we analyzed and compared the complete mtDNA D-loop sequences from 17 pairs of tumor tissues and corresponding peripheral blood samples using the direct DNA sequencing method. Our results revealed that 12 of the 17 EWS tumor specimens (70.6%) carried 19 somatic mutations in the D-loop of mtDNA, including 11 single-base substitutions, 3 insertions and 5 deletions. Among the tested 17 patients, we screened a total of 40 germline polymorphisms including one novel sequence variant in the D-loop fragment. Most of these identified mutations and germline variations were clustered within two hypervariable segments (HVS1 and HVS2) as well as the homopolymeric C stretch between nucleotide position 303 and 309. In addition, there was no significant correlation between mtDNA D-loop mutations and various clinicopathological factors of EWS. In conclusion, our study reports for the first time that mtDNA D-loop mutations occur at a high frequency in EWS. These data provide evidence of mtDNA alterations’ possible involvement in the initiation and/or progression of this rare malignancy.     

Heteroplasmy of Short Tandem Repeats in Mitochondrial DNA of Atlantic Cod, Gadus Morhua

The mitochondrial DNA of the Atlantic cod (Gadus morhua) contains a tandem array of 40-bp repeats in the D-loop region of the molecule. Variation among molecules in the copy number of these repeats results in mtDNA length variation and heteroplasmy (the presence of more than one form of mtDNA in an individual). In a sample of fish collected from different localities around Iceland and off George's Bank, each individual was heteroplasmic for two or more mtDNAs ranging in repeat copy number from two (common) to six (rare). An earlier report on mtDNA heteroplasmy in sturgeon (Acipenser transmontanus) presented a competitive displacement model for length mutations in mtDNAs containing tandem arrays and the cod data deviate from this model. Depending on the nature of putative secondary structures and the location of D-loop strand termination, additional mechanisms of length mutation may be needed to explain the range of mtDNA length variants maintained in these populations. The balance between genetic drift and mutation in maintaining this length polymorphism is estimated through a hierarchical analysis of diversity of mtDNA length variation in the Iceland samples. Eighty percent of the diversity lies within individuals, 8% among individuals and 12% among localities. An estimate of theta = 2N(eo) mu greater than 1 indicates that this system is characterized by a high mutation rate and is governed primarily by deterministic dynamics. The sequences of repeat arrays from fish collected in Norway, Iceland and George's Bank show no nucleotide variation suggesting that there is very little substructuring to the North Atlantic cod population.     

Stock Structure and Homing Fidelity in Gulf of Mexico Sturgeon (Acipenser Oxyrinchus Desotoi) Based on Restriction Fragment Length Polymorphism and Sequence Analyses of Mitochondrial DNA

Efforts have been proposed worldwide to restore sturgeon populations through the use of hatcheries to supplement natural reproduction and to reintroduce sturgeon where they have become extinct. We examined the population structure and inferred the extent of homing in the anadromous Gulf of Mexico (Gulf) sturgeon (Acipenser oxyrinchus desotoi). Restriction fragment length polymorphism and control region sequence analyses of mitochondrial DNA (mtDNA) were used to identify haplotypes of Gulf sturgeon specimens obtained from eight drainages spanning the subspecies' entire distribution from Louisiana to Florida. Significant differences in haplotype frequencies indicated substantial geographic structuring of populations. A minimum of four regional or river-specific populations were identified (from west to east): (1) Pearl River, LA and Pascagoula River, MS, (2) Escambia and Yellow rivers, FL, (3) Choctawhatchee River, FL, and (4) Apalachicola, Ochlockonee, and Suwannee rivers, FL. Estimates of maternally mediated gene flow between any pair of the four regional or river-specific stocks ranged between 0.15 to 1.2. Tandem repeats in the mtDNA control region of Gulf sturgeon were not perfectly conserved. This result, together with an absence of heteroplasmy and length variation in Gulf sturgeon mtDNA, indicates that the molecular mechanisms of mtDNA control region sequence evolution differ among acipenserids.     

Genetic diversity of Mongolian gazelle Procapra guttorosa Pallas, 1777

The mitochondrial DNA D-loop hypervariable fragment sequence polymorphism was examined in 27 Mongolian gazelles from Mongolia, Russia, and China. Intraspecific polymorphism of the D-loop fragment examined was demonstrated. All haplotypes described were unique. The average nucleotide diversity (pi) for the mtDNA fragment investigated constituted 5.85 +/- 2.92%. A relatively high number of insertions and deletions was observed. In particular, a haplotype with the 77-bp insertion was described. The data obtained point to high genetic diversity of Mongolian populations. There was no correlation between the distribution of haplotypes examined and geographical location of the animal tissue sampling sites.     

MtDNA substitution rate and segregation of heteroplasmy in coding and noncoding regions

The mitochondrial DNA (mtDNA) substitution rate and segregation of heteroplasmy were studied for the non-coding control region (D-loop) and 500 bp of the coding region between nucleotide positions 5550 and 6050, by sequence analysis of blood samples from 194 individuals, representing 33 maternal lineages. No homoplasmic nucleotide substitutions were detected in a total of 292 transmissions. The estimated substitution rate per nucleotide per million years for the control region (micro>0.21, 95% CI 0-0.6) was not significantly different from that for the coding region (micro>0.54, 95% CI 0-1.0). Variation in the length of homopolymeric C streches was observed at three sites in the control region (positions 65, 309 and 16,189), all of which were in the heteroplasmic state. Segregation of heteroplasmic genotypes between generations was observed in several maternal pedigrees. At position 309, a longer poly C tract length was strongly associated with a higher probability for heteroplasmy and rapid segregation between generations. The length heteroplasmy at positions 65 and 16,189 was found at low frequency and was confined to a few families.     

五种鲟鱼线粒体控制区异质性和系统发育分析

利用保守引物得到五种鲟鱼的线粒体DNA（mtDNA）控制区（D-loop）全长,长度在795~813 bp。序列中包括了CBS（conserved sequence block）和TAS（termination-associated sequence）区域。利用最大似然法、最大简约法和贝叶斯法构建了系统发育树,发育树分成两枝,呈现明显的生物地理分布。分析表明,现有的鳇属鱼类不是单系群起源。五种鲟鱼D-loop序列都存在长度和数目不等串联重复序列,长度在78~82 bp之间,重复序列拷贝数在4~6次不等,因此造成了mtDNA广泛的异质性现象。不同种类的重复序列单元十分相似,达氏鳇和史氏鲟重复序列单元相似度为82.93%,西伯利亚鲟和俄罗斯鲟重复序列单元相似度为90.59%。在串联重复序列后是一段不完全重复序列。通过与已有同种的重复序列比对发现不同鲟鱼重复序列相同,不同地理区域相同物种的重复序列可能发生过分子内重组。这些表明重复序列在鲟鱼进化上具有相关意义,推测重复序列可能产生在种分化前,重组发生在种分化后。     

Population variation of human mtDNA control region sequences detected by enzymatic amplification and sequence-specific oligonucleotide probes.

A method for detecting sequence variation of hypervariable segments of the mtDNA control region was developed. The technique uses hybridization of sequence-specific oligonucleotide (SSO) probes to DNA sequences that have been amplified by PCR. The nucleotide sequences of the two hypervariable segments of the mtDNA control region from 52 individuals were determined; these sequences were then used to define nine regions suitable for SSO typing. A total of 23 SSO probes were used to detect sequence variants at these nine regions in 525 individuals from five ethnic groups (African, Asian, Caucasian, Japanese, and Mexican). The SSO typing revealed an enormous amount of variability, with 274 mtDNA types observed among these 525 individuals and with diversity values, for each population, exceeding .95. For each of the nine mtDNA regions significant differences in the frequencies of sequence variants were observed between these five populations. The mtDNA SSO-typing system was successfully applied to a case involving individual identification of skeletal remains; the probability of a random match was approximately 0.7%. The potential useful applications of this mtDNA SSO-typing system thus include the analysis of individual identity as well as population genetic studies.     

Substitution rate variation among sites in mitochondrial hypervariable region I of humans and chimpanzees

Mitochondrial D-loop hypervariable region I (HVI) sequences are widely used in human molecular evolutionary studies, and therefore accurate assessment of rate heterogeneity among sites is essential. We used the maximum-likelihood method to estimate the gamma shape parameter alpha for variable substitution rates among sites for HVI from humans and chimpanzees to provide estimates for future studies. The complete data of 839 humans and 224 chimpanzees, as well as many subsets of these data, were analyzed to examine the effect of sequence sampling. The effects of the genealogical tree and the nucleotide substitution model were also examined. The transition/transversion rate ratio (kappa) is estimated to be about 25, although much larger and biased estimates were also obtained from small data sets at low divergences. Estimates of alpha were 0.28-0.39 for human data sets of different sizes and 0.20-0.39 for data sets including different chimpanzee subspecies. The combined data set of both species gave estimates of 0.42-0.45. While all those estimates suggest highly variable substitution rates among sites, smaller samples tend to give smaller estimates of alpha. Possible causes for this pattern were examined, such as biases in the estimation procedure and shifts in the rate distribution along certain lineages. Computer simulations suggest that the estimation procedure is quite reliable for large trees but can be biased for small samples at low divergences. Thus, an alpha of 0.4 appears suitable for both humans and chimpanzees. Estimates of alpha can be affected by the nucleotide sites included in the data, the overall tree length (the amount of sequence divergence), the number of rate classes used for the estimation, and to a lesser extent, the included sequences. The genealogical tree, the substitution model, and demographic processes such as population expansion do not have much effect.     

Human Mitochondrial DNA Variation and Evolution: Analysis of Nucleotide Sequences from Seven Individuals

We have analyzed nucleotide sequence variation in an approximately 900-base pair region of the human mitochondrial DNA molecule encompassing the heavy strand origin of replication and the D-loop. Our analysis has focused on nucleotide sequences available from seven humans. Average nucleotide diversity among the sequences is 1.7%, several-fold higher than estimates from restriction endonuclease site variation in mtDNA from these individuals and previously reported for other humans. This disparity is consistent with the rapidly evolving nature of this noncoding region. However, several instances of convergent or parallel gain and loss of restriction sites due to multiple substitutions were observed. In addition, other results suggest that restriction site (as well as pairwise sequence) comparisons may underestimate the total number of substitutions that have occurred since the divergence of two mtDNA sequences from a common ancestral sequence, even at low levels of divergence. This emphasizes the importance of recognizing the large standard errors associated with estimates of sequence variability, particularly when constructing phylogenies among closely related sequences. Analysis of the observed number and direction of substitutions revealed several significant biases, most notably a strand dependence of substitution type and a 32-fold bias favoring transitions over transversions. The results also revealed a significantly nonrandom distribution of nucleotide substitutions and sequence length variation. Significantly more multiple substitutions were observed than expected for these closely related sequences under the assumption of uniform rates of substitution. The bias for transitions has resulted in predominantly convergent or parallel changes among the observed multiple substitutions. There is no convincing evidence that recombination has contributed to the mtDNA sequence diversity we have observed.     

口腔鳞状细胞癌线粒体DNA的HVRⅢ区突变的研究（英文）

目的：检测口腔鳞状细胞癌患者线粒体DNA复制控制区（mtDNA D-loop）高变Ⅲ区（hypervariable regionⅢ,HVRⅢ）的突变情况,并探讨其意义。方法：以口腔鳞状细胞癌患者癌旁组织及正常组织作为对照,对7例口腔鳞状细胞癌组织样本的mtDNA D-loop HVRⅢ区进行PCR扩增和测序分析。结果：在7例患者的癌组织、癌旁组织、正常组织样本中共发现72个（56种）核苷酸改变,其中51个（26种）为核苷酸多态性改变;3个肿瘤组织样本中共发现21个突变,其中16个位于HVRⅢ区范围内;癌旁组织及正常组织未发现突变;口腔鳞状细胞癌的mtDNA D-loop HVRⅢ区突变率为42.9%（3/7）。结论：mtDNA D-loop HVRⅢ区的变异可能与口腔鳞状细胞癌的易感性有一定的联系;本研究为寻找新的肿瘤基因诊断和肿瘤遗传易感性的标志物提供了依据。