Site heteroplasmy in the mitochondrial cytochrome b gene of the sterlet sturgeon Acipenser ruthenus

Sturgeons are fish species with a complex biology. They are also characterized by complex aspects including polyploidization and easiness of hybridization. As with most of the Ponto-Caspian sturgeons, the populations of Acipenser ruthenus from the Danube have declined drastically during the last decades. This is the first report on mitochondrial point heteroplasmy in the cytochrome b gene of this species. The 1141 bp sequence of the cytb gene in wild sterlet sturgeon individuals from the Lower Danube was determined, and site heteroplasmy evidenced in three of the 30 specimens collected. Two nucleotide sequences were identified in these heteroplasmic individuals. The majority of the heteroplasmic sites are synonymous and do not modify the sequence of amino acids in cytochrome B protein. To date, several cases of point heteroplasmy have been reported in animals, mostly due to paternal leakage of mtDNA. The presence of specific point heteroplasmic sites might be interesting for a possible correlation with genetically distinct groups in the Danube River.     

Structures of two genes coding for melanin-concentrating hormone of chum salmon

Two MCH genes coding for melanin-concentrating hormone (MCH) were isolated from a chum salmon liver DNA library and characterized. They were shown to be intronless genes with 0.63-kb exons, each of which commonly consisted of an about 80-bp 5'-untranslated region, a region coding for 132 amino acids (aa) MCH precursor protein and an approx. 160-bp 3'-untranslated region. About 20 bp upstream from the putative cap site, sequences were found corresponding to the TATA box. The two genes were 86% identical at the nucleotide sequence level. Sequences homologous to the chum salmon MCH genes were present in the genomes of other fish such as catfish, carp and Chinese grass carp, whereas no highly homologous sequence could be detected in other vertebrate genomes.     

Reevaluation of amino acid variability of the human immunodeficiency virus type 1 gp120 envelope glycoprotein and prediction of new discontinuous epitopes

To elucidate the evolutionary mechanisms of the human immunodeficiency virus type 1 gp120 envelope glycoprotein at the single-site level, the degree of amino acid variation and the numbers of synonymous and nonsynonymous substitutions were examined in 186 nucleotide sequences for gp120 (subtype B). Analyses of amino acid variabilities showed that the level of variability was very different from site to site in both conserved (C1 to C5) and variable (V1 to V5) regions previously assigned. To examine the relative importance of positive and negative selection for each amino acid position, the numbers of synonymous and nonsynonymous substitutions that occurred at each codon position were estimated by taking phylogenetic relationships into account. Among the 414 codon positions examined, we identified 33 positions where nonsynonymous substitutions were significantly predominant. These positions where positive selection may be operating, which we call putative positive selection (PS) sites, were found not only in the variable loops but also in the conserved regions (C1 to C4). In particular, we found seven PS sites at the surface positions of the alpha-helix (positions 335 to 347 in the C3 region) in the opposite face for CD4 binding. Furthermore, two PS sites in the C2 region and four PS sites in the C4 region were detected in the same face of the protein. The PS sites found in the C2, C3, and C4 regions were separated in the amino acid sequence but close together in the three-dimensional structure. This observation suggests the existence of discontinuous epitopes in the protein's surface including this alpha-helix, although the antigenicity of this area has not been reported yet.     

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.     

A sensitive denaturing gradient-Gel electrophoresis assay reveals a high frequency of heteroplasmy in hypervariable region 1 of the human mtDNA control region

A population study of heteroplasmy in the hypervariable region 1 (HV1) portion of the human mtDNA control region was performed. Blood samples from 253 randomly chosen individuals were examined using a sensitive denaturing gradient-gel electrophoresis (DGGE) system. This method is capable of detecting heteroplasmic proportions as low as 1% and virtually all heteroplasmy where the minor component is > or = 5%. Heteroplasmy was observed in 35 individuals (13.8%; 95% confidence interval [CI] 9.6-18.0). Of these individuals, 33 were heteroplasmic at one nucleotide position, whereas 2 were heteroplasmic at two different positions (a condition known as "triplasmy"). Although heteroplasmy occurred at a total of 16 different positions throughout HV1, it was most frequently observed at positions 16093 (n=13) and 16129 (n=6). In addition, the majority of heteroplasmic variants occurred at low proportions and could not be detected by direct sequencing of PCR products. This study indicates that low-level heteroplasmy in HV1 is relatively common and that it occurs at a broad spectrum of sites. Our results corroborate those of other recent reports indicating that heteroplasmy in the control region is more common than was previously believed-a finding that is of potential importance to evolutionary studies and forensic applications that are based on mtDNA variation.     

A novel mtDNA C11777A mutation in Leigh syndrome

A novel mitochondrial DNA point mutation, a C-to-A mutation at nucleotide position (np) 11,777, was identified in two unrelated patients out of 100 with Leigh syndrome. This mutation converted a highly evolutionary conserved arginine to a serine at codon 340 in ND4 gene. This codon was also converted by a G-to-A mutation at np 11,778, the most common mutation associated with Leber's hereditary optic neuropathy (LHON), but the amino acid replacement was different (R340S vs. R340H). Cybrid study revealed that the percentage of heteroplasmy was correlated with complex I function and that the novel mutation caused a much more deleterious effect than the np 11,778 LHON mutation in complex I activity.     

Models of amino acid substitution and applications to mitochondrial protein evolution

Models of amino acid substitution were developed and compared using maximum likelihood. Two kinds of models are considered. "Empirical" models do not explicitly consider factors that shape protein evolution, but attempt to summarize the substitution pattern from large quantities of real data. "Mechanistic" models are formulated at the codon level and separate mutational biases at the nucleotide level from selective constraints at the amino acid level. They account for features of sequence evolution, such as transition-transversion bias and base or codon frequency biases, and make use of physicochemical distances between amino acids to specify nonsynonymous substitution rates. A general approach is presented that transforms a Markov model of codon substitution into a model of amino acid replacement. Protein sequences from the entire mitochondrial genomes of 20 mammalian species were analyzed using different models. The mechanistic models were found to fit the data better than empirical models derived from large databases. Both the mutational distance between amino acids (determined by the genetic code and mutational biases such as the transition-transversion bias) and the physicochemical distance are found to have strong effects on amino acid substitution rates. A significant proportion of amino acid substitutions appeared to have involved more than one codon position, indicating that nucleotide substitutions at neighboring sites may be correlated. Rates of amino acid substitution were found to be highly variable among sites.      

A twin study of mitochondrial DNA polymorphisms shows that heteroplasmy at multiple sites is associated with mtDNA variant 16093 but not with zygosity

The mitochondrial theory of ageing proposes that damage to mitochondria and diminished mitochondrial DNA (mtDNA) repair are major contributors to cellular dysfunction and age-related diseases. We investigate the prevalence of heteroplasmy in the mtDNA control region in buccal swab and blood derived samples for 178 women from the TwinsUK cohort (41 DZ pair 39 MZ pairs, 18 singletons, mean age 57.5 range 28–82) and its relationship to age, BMI and fasting insulin and glucose serum levels. The overall estimated prevalence of heteroplasmy for both tissues in the control region measured for 37 sites was 17%. The prevalence of heteroplasmy was higher among the older half of the study subjects than in the younger half (23% vs 10% p<0.03), primarily reflecting the increase in the prevalence of a heteroplasmic dinucleotide CA repeat in variable region II (VRII) with age. The VRII 523–524 heteroplasmic site (heteroplasmic in 25 subjects) was also associated with a decrease in BMI. In addition, concordance rates for common heteroplasmy were observed to be near complete for both dizygotic (DZ = 94%) and monozygotic twin pairs (MZ = 100%), consistent with previous reports that suggest variation in heteroplasmy rates between generations are determined by bottlenecks in maternal transmission of mitochondria. Differences in the prevalence of heteroplasmy were observed overall between samples derived from buccal swabs (19%) and blood (15%, p<0.04). These were particularly marked at position 16093 of hypervariable region I (HVI, 7% vs 0%, respectively, p<4×10⁻¹¹). The presence of the C allele at position 16093 in blood was associated with the presence of heteroplasmy in buccal swabs at this position (p = 3.5×10⁻¹⁴) and also at VRII (p = 2×10⁻⁴) suggesting a possible predisposing role for this site in the accumulation of heteroplasmy. Our data indicate that BMI is potentially associated with control region heteroplasmy.     

Restriction site heteroplasmy in the mitochondrial DNA of the marine fish Sciaenops ocellatus (L.)

Restriction site heteroplasmy involving the enzymes NcoI and XbaI was detected in the mitochondrial DNAs of two individuals of the marine fish Sciaenops ocellatus. This represents only the sixth documented example of mitochondrial DNA restriction site heteroplasmy in animals. Two heteroplasmic individuals were found in a survey of nearly 750 individuals, suggesting that in most studies the incidence of mitochondrial DNA site heteroplasmy may be too low to be routinely detected.     

Detecting Heteroplasmy from High-Throughput Sequencing of Complete Human Mitochondrial DNA Genomes

Heteroplasmy, the existence of multiple mtDNA types within an individual, has been previously detected by using mostly indirect methods and focusing largely on just the hypervariable segments of the control region. Next-generation sequencing technologies should enable studies of heteroplasmy across the entire mtDNA genome at much higher resolution, because many independent reads are generated for each position. However, the higher error rate associated with these technologies must be taken into consideration to avoid false detection of heteroplasmy. We used simulations and phiX174 sequence data to design criteria for accurate detection of heteroplasmy with the Illumina Genome Analyzer platform, and we used artificial mixtures and replicate data to test and refine the criteria. We then applied these criteria to mtDNA sequence reads for 131 individuals from five Eurasian populations that had been generated via a parallel tagged approach. We identified 37 heteroplasmies at 10% frequency or higher at 34 sites in 32 individuals. The mutational spectrum does not differ between heteroplasmic mutations and polymorphisms in the same individuals, but the relative mutation rate at heteroplasmic mutations is significantly higher than that estimated for all mutable sites in the human mtDNA genome. Moreover, there is also a significant excess of nonsynonymous mutations observed among heteroplasmies, compared to polymorphism data from the same individuals. Both mutation-drift and negative selection influence the fate of heteroplasmies to determine the polymorphism spectrum in humans. With appropriate criteria for avoiding false positives due to sequencing errors, next-generation technologies can provide novel insights into genome-wide aspects of mtDNA heteroplasmy.     

An analysis of nucleotide and amino acid variability in the barcode region of cytochrome c oxidase I (cox1) in fishes

The 655 bp cytochrome c oxidase subunit I barcode region of single specimens of 388 species of fishes (four Holocephali, 61 Elasmobranchii and 323 Actinopterygii) was examined. All but two (Urolophus cruciatus and Urolophus sufflavus) showed different cox1 nucleotide sequences (99.5% species discrimination); the two that could not be resolved are suspected to hybridize. Most of the power of cox1 nucleotide sequence analysis for species identification comes from the degenerate nature of the genetic code and the highly variable nature of the third codon position of amino acids. Variation at the third codon position is bimodally distributed, and the more variable mode is dominated by amino acids with four or six codons, while the less variable mode is dominated by amino acids with two codons. The ratio of nonsynonymous to synomymous changes is much less than one, indicating that this gene is subject to strong purifying selection. Consequently, cox1 amino acid sequence diversity is much less than nucleotide sequence diversity and has very poor species resolution power. Fourteen of the 16 amino acid residues recognized as having important functions in the region of cox1 sequenced were completely conserved over all 388 species (and the bovine cox1 sequence), with one fish species varying at one of these sites, and three fish at another site. No significant differences in amino acid conservation were observed between residues in helices, strands and turns. Patterns of nucleotide and amino acid variability were very similar between elasmobranchs and actinopterygians.     

Investigation of Heteroplasmy in the Human Mitochondrial DNA Control Region: A Synthesis of Observations from More Than 5000 Global Population Samples

Instances of point and length heteroplasmy in the mitochondrial DNA control region were compiled and analyzed from over 5,000 global human population samples. These data represent observations from a large and broad population sample, representing nearly 20 global populations. As expected, length heteroplasmy was frequently observed in the HVI, HVII and HVIII C-stretches. Length heteroplasmy was also observed in the AC dinucleotide repeat region, as well as other locations. Point heteroplasmy was detected in approximately 6% of all samples, and while the vast majority of heteroplasmic samples comprised two molecules differing at a single position, samples exhibiting two and three mixed positions were also observed in this data set. In general, the sites at which heteroplasmy was most commonly observed correlated with reported control region mutational hotspots. However, for some sites, observations of heteroplasmy did not mirror established mutation rate data, suggesting the action of other mechanisms, both selective and neutral. Interestingly, these data indicate that the frequency of heteroplasmy differs between particular populations, perhaps reflecting variable mutation rates among different mtDNA lineages and/or artifacts of particular population groups. The results presented here contribute to our general understanding of mitochondrial DNA control region heteroplasmy and provide additional empirical information on the mechanisms contributing to mtDNA control region mutation and evolution. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Complete nucleotide sequence of Japanese flounder (Paralichthys olivaceus) mitochondrial genome: structural properties and cue for resolving teleostean relationships

We cloned and sequenced the complete mitochondrial genome of Japanese flounder (Paralichthys olivaceus). A circular 17,090 bp mitochondrial genome from the flounder contains 37 structural genes as in other vertebrates so far reported. This is the first report of the complete mitochondrial sequence from a higher teleostean fish (Acanthopterygii). The organization including gene order is quite similar to that of other teleostean fishes as well as placental mammals. The putative control region of the Japanese flounder mitochondrial genome contains a length variable region of about a 74 bp tandem repeat cluster. As a preliminary study we adopted the maximum likelihood and neighbor-joining inference methods to examine phylogenetic relationships among teleostean and related fishes. Comparisons of amino acid sequences of protein-coding genes and nucleotide sequences of tRNA genes resolved some middle to deep branches among some teleostean fishes. The flounder mitochondrial genome does not show an indication of evolutionary rate difference among teleosts leading to difficulty in phylogenetic analyses, and our data is useful for future evolutionary studies dealing with higher teleostean fishes.     

Analysis of heteroplasmy in the major noncoding region of mitochondrial DNA in the blood and atherosclerotic plaques of carotid arteries

For identification of somatic mitochondrial DNA (mtDNA) mutations, the mtDNA major noncoding region (D-loop) sequence in blood samples and carotid atherosclerosis plaques from patients with atherosclerosis was analyzed. Five point heteroplasmic positions were observed in 4 of 23 individuals (17%). Only in two cases could heteroplasmy have resulted from somatic mutation, whereas three heteroplasmic positions were found in both vascular tissue and blood. In addition, length heteroplasmy in a polycytosine stretches was registered at nucleotide positions 303–315 in 16 individuals, and also in the 16184–16193 region in four patients. The results suggest that somatic mtDNA mutations can occur during atherosclerosis, but some heteroplasmic mutations may appear in all tissues, possibly being inherited.     

Sequence Heteroplasmy of D-Loop and rRNA Coding Regions in Mitochondrial DNA from Holstein Cows of Independent Maternal Lineages

A mitochondrial DNA (mtDNA) fragment containing the D-loop, phenylalanine tRNA, valine tRNA, and 12S and 16 rRNA genes was cloned and sequenced from 36 cows of 18 maternal lineages to identify the polymorphic sites within those regions and to detect the existence of heteroplasmic mtDNA in cows. Seventeen variable sites were observed within the D-loop and rRNA coding regions of bovine mtDNA within a 2.5-kb span. The hypervariable sites in the D-loop and rRNA coding regions were identified at nucleotide positions 169, 216, and 1594. Heteroplasmic mtDNA (variable mtDNA within a tissue) existed extensively in cows and was detected within the above regions in 11 of 36 cows sequenced. The insertion, deletion, and nucleotide transversion polymorphisms were found only in homopolymer regions. Heteroplasmy was observed frequently and seemingly is persistent in cattle. Though heteroplasmy was demonstrated, most lineages and mtDNA sites showed no heteroplasmy.     

Uncorrected nucleotide bias in mtDNA can mimic the effects of positive Darwinian selection

The relative rates of nucleotide substitution at synonymous and nonsynonymous sites within protein-coding regions have been widely used to infer the action of natural selection from comparative sequence data. It is known, however, that mutational and repair biases can affect rates of evolution at both synonymous and nonsynonymous sites. More importantly, it is also known that synonymous sites are particularly prone to the effects of nucleotide bias. This means that nucleotide biases may affect the calculated ratio of substitution rates at synonymous and nonsynonymous sites. Using a large data set of animal mitochondrial sequences, we demonstrate that this is, in fact, the case. Highly biased nucleotide sequences are characterized by significantly elevated dN/dS ratios, but only when the nucleotide frequencies are not taken into account. When the analysis is repeated taking the nucleotide frequencies at each codon position into account, such elevated ratios disappear. These results suggest that the recently reported differences in dN/dS ratios between vertebrate and invertebrate mitochondrial sequences could be explained by variations in mitochondrial nucleotide frequencies rather than the effects of positive Darwinian selection.     

Quantification of point mutations associated with Leber hereditary optic neuroretinopathy by solid-phase minisequencing

About two-thirds of patients with Leber hereditary optic neuroretinopathy (LHON) harbor mutations in mitochondrial DNA at positions 11778 (ND4) or 3460 (ND1). Thus, the clinical diagnosis of LHON can often be confirmed with mutation analysis. Detection of pathogenic mutations and quantification of heteroplasmy has mainly relied on PCR and restriction site analysis and densitometric scanning. We applied the recently developed solid-phase minisequencing method, based on primerguided nucleotide incorporation, to the simultaneous detection and quantitation of the ND4/11778 and ND1/3460 mutations. The method was highly sensitive, heteroplasmy as low as 1.5% being easily detected. Rapid, reproducible, and accurate results prove solid-phase minisequencing to be the method of choice for quantitative analysis of LHON mutations.     

Rapid shift in genotype of human mitochondrial DNA in a family with Leber's hereditary optic neuropathy

Mitochondrial DNA isolated from white blood cells was investigated in families suffering from Leber's hereditary optic neuropathy. A recently described mutation at nucleotide position 11778 was present in 5 out of 12 families and heteroplasmic mitochondrial DNA was observed in 2 of these 5 families. A rapid shift in genotype was found in one of the families with heteroplasmy: the grandmother had 60 percent mitochondrial DNA mutated at nucleotide position 11778, the mother 55 percent, and the two sons at least 95 percent. These data indicate that the number of mitochondrial DNA molecules transmitted to the progeny passes a developmental bottleneck, as previously proposed to occur in bovine oogenesis.     

Cloning and characterisation of a cDNA encoding Japanese flounder Paralichthys olivaceus IgD

A cDNA containing the gene for Japanese flounder IgD consisted of 3240 bp encoding 998 amino acid residues. The amino acid sequence of the constant region of Japanese flounder IgD shares 38-80% identity with the sequences of previously reported teleost IgDs. The structure of the constant region of Japanese flounder IgD, which contains the micro1, delta1, delta2, delta3, delta4, delta5, delta6, delta7, and TM regions, is similar to the structures of the constant regions of the IgDs of channel catfish and Atlantic salmon. Southern blot hybridisation showed that the Japanese flounder IgD gene exists as a single locus. The Japanese flounder IgD gene was mainly detected in peripheral blood leucocytes (PBLs) and small amounts were detected in the spleen, head and trunk kidney, although IgM mRNA was detected in similar amounts in PBLs, the head kidney, and spleen. The copy number of IgM mRNA in Japanese flounder PBL was 56-fold higher than that of IgD.     

Assessment of concordance among genealogical reconstructions from various mtDNA segments in three species of Pacific salmon (genus Oncorhynchus)

Seven segments of mitochondrial DNA (mtDNA), comprising 97% of the mitochondrial genome, were amplified by polymerase chain reaction (PCR) and examined for restriction site variation using 13 restriction endonucleases in three species of Pacific salmon: pink (Oncorhynchus gorbuscha), chum (O. keta) and sockeye (O. nerka) salmon. The distribution of variability across the seven mtDNA segments differed substantially among species. Little similarity in the distribution of variable restriction sites was found even between the mitochondrial genomes of the even- and odd-year broodlines of pink salmon. Significantly different levels of nucleotide diversity were detected among three groups of genes: six NADH-dehydrogenase genes had the highest; two rRNA genes had the lowest; and a group that included genes for ATPase and cytochrome oxidase subunits, the cytochrome b gene, and the control region had intermediate levels of nucleotide diversity. Genealogies of mtDNA haplotypes were reconstructed for each species, based on the variation in all mtDNA segments. The contributions of variation within different segments to resolution of the genealogical trees were compared within each species. With the exception of sockeye salmon, restriction site data from different genome segments tended to produce rather different trees (and hence rather different genealogies). In the majority of cases, genealogical information in different segments of mitochondrial genome was additive rather than congruent. This finding has a relevance to phylogeographic studies of other organisms and emphasizes the importance of not relying on a limited segment of the mtDNA genome to derive a phylogeographic structure.