Review and meta‐analysis of natural selection in mitochondrial complex I in metazoans

Variation in mitochondrial DNA is often assumed to be neutral and is used to construct the genealogical relationships among populations and species. However, if extant variation is the result of episodes of positive selection, these genealogies may be incorrect, although this information itself may provide biologically and evolutionary meaningful information. In fact, positive Darwinian selection has been detected in the mitochondrial‐encoded subunits that comprise complex I from diverse taxa with seemingly dissimilar bioenergetic life histories, but the functional implications of the selected sites are unknown. Complex I produces roughly 40% of the proton flux that is used to synthesize ATP from ADP, and a functional model based on the high‐resolution structure of complex I described a unique biomechanical apparatus for proton translocation. We reported positive selection at sites in this apparatus during the evolution of Pacific salmon, and it appeared this was also the case in published reports from other taxa, but a comparison among studies was difficult because different statistical tests were used to detect selection and oftentimes, specific sites were not reported. Here we review the literature of positive selection in mitochondrial genomes, the statistical tests used to detect selection, and the structural and functional models that are currently available to study the physiological implications of selection. We then search for signatures of positive selection among the coding mitochondrial genomes of 237 species with a common set of tests and verify that the ND5 subunit of complex I is a repeated target of positive Darwinian selection in diverse taxa. We propose a novel hypothesis to explain the results based on their bioenergetic life histories and provide a guide for laboratory and field studies to test this hypothesis.     

Nonsynonymous site heteroplasmy in fish mitochondrial DNA

Heteroplasmic nucleotide polymorphisms are rarely observed in wild animal mitochondrial DNA. The occurrence of such site heteroplasmy is expected to be extremely rare at nonsynonymous sites where the number of nucleotide substitutions per site is low due to functional constraints. This report deals with nonsynonymous mitochondrial heteroplasmy from two wild fish species, chum salmon and Japanese flounder. We detected an A/C nonsynonymous heteroplasmic site corresponding to putative amino acids, Ile or Met, in NADH dehydrogenase subunit-5 (ND5) region of chum salmon. The heteroplasmic site was at the 3rd position of 58th codon. As for Japanese flounder we detected a C/T nonsynonymous heteroplasmic site corresponding to putative amino acids, Leu or Pro, in ND4 region. The heteroplasmic site was at the 2nd position of 450th codon. We also verified heteroplasmy at these sites by sequencing cloned fragments.     

Evolution of the salmonid mitochondrial control region.

To explore the evolutionary nature of the salmonid mitochondrial DNA (mtDNA) control region (D-loop) and its utility for inferring phylogenies, the entire region was sequenced from all eight species of anadromous Pacific salmon, genus Oncorhynchus; the Atlantic salmon, Salmo salar; and the Arctic grayling, Thymallus arcticus. A comparison of aligned sequences demonstrates that the generally conserved sequence elements that have been previously reported for other vertebrates are maintained in these primitive teleost fishes. Results reveal a significantly nonrandom distribution of nucleotide substitutions, insertions, and deletions that suggests that portions of the salmonid D-loop may be under differential selective constraints and that most of the control region of these fishes may evolve at a rate similar to that of the remainder of their mtDNA genomes. Maximum likelihood and Fitch parsimony analyses of 9 kb of aligned salmonid sequence data give evolutionary trees of identical topology. These results are consistent with previous molecular studies of a limited number of salmonid taxa and with more comprehensive, classical analyses of salmonid evolution. Predictions from these data, based on a molecular clock assumption for the mtDNA control region, are also consistent with fossil evidence that suggests that species of Oncorhynchus could be as old as the Middle Pliocene and would have thus given rise to the extant Pacific salmon prior to about 5 or 6 million years ago.     

Identification of the sex chromosome pair in chum salmon (Oncorhynchus keta) and pink salmon (Oncorhynchus gorbuscha)

Fluorescence in situ hybridization (FISH) using a probe to the male-specific GH-Y (growth hormone pseudogene) was used to identify the Y chromosome in the karyotypes of chum salmon (Oncorhynchus keta) and pink salmon (Oncorhynchus gorbuscha). The sex chromosome pair is a small acrocentric chromosome pair in chum salmon and the smallest metacentric chromosome pair in pink salmon. Both of these chromosome pairs are morphologically different from the sex chromosome pairs in chinook salmon (Oncorhynchus tshawytscha) and coho salmon (Oncorhynchus kisutch). The 5S rRNA genes are on multiple chromosome pairs including the sex chromosome pair in chum salmon, but at the centromeres of two autosomal metacentric pairs in pink salmon. The sex chromosome pairs and the chromosomal locations of the 5S rDNA appear to be different in all five of the North American Pacific salmon species and rainbow trout. The implications of these results for evolution of sex chromosomes in salmonids are discussed.     

Putative Diphyllobothrium nihonkaiense acquired from a Pacific salmon (Oncorhynchus keta) eaten in France; genomic identification and case report

We report here a likely case of Diphyllobothrium nihonkaiense contracted in France through the consumption of a Pacific salmon imported from Canada. The species diagnosis was made by molecular analysis of two mitochondrial genes (COI & ND3). This case is rather unusual in that D. nihonkaiense has never been reported along the Pacific coast of North America.     

PCR‐RFLP markers detect 29 mitochondrial haplotypes in Pacific lamprey (Entosphenus tridentatus)

We developed five polymerase chain reaction‐based markers that detect variation in the mitochondrial genome of the Pacific lamprey, Entosphenus tridentatus, across most of its range. Two gene fragments (ND2 and ND5) were amplified and digested with three and two restriction enzymes, respectively, detecting sequence variation at 18 sites (ND2 = 13; ND5 = 5) and yielding 29 composite haplotypes among 1246 lampreys. These sequence‐based markers will be useful in a range of phylogeographical and population genetic studies.     

Conflict Among Individual Mitochondrial Proteins in Resolving the Phylogeny of Eutherian Orders

The phylogenetic relationship among primates, ferungulates (artiodactyls + cetaceans + perissodactyls + carnivores), and rodents was examined using proteins encoded by the H strand of mtDNA, with marsupials and monotremes as the outgroup. Trees estimated from individual proteins were compared in detail with the tree estimated from all 12 proteins (either concatenated or summing up log-likelihood scores for each gene). Although the overall evidence strongly suggests ((primates, ferungulates), rodents), the ND1 data clearly support another tree, ((primates, rodents), ferungulates). To clarify whether this contradiction is due to (1) a stochastic (sampling) error; (2) minor model-based errors (e.g., ignoring site rate variability), or (3) convergent and parallel evolution (specifically between either primates and rodents or ferungulates and the outgroup), the ND1 genes from many additional species of primates, rodents, other eutherian orders, and the outgroup (marsupials + monotremes) were sequenced. The phylogenetic analyses were extensive and aimed to eliminate the following artifacts as possible causes of the aberrant result: base composition biases, unequal site substitution rates, or the cumulative effects of both. Neither more sophisticated evolutionary analyses nor the addition of species changed the previous conclusion. That is, the statistical support for grouping rodents and primates to the exclusion of all other taxa fluctuates upward or downward in quite a tight range centered near 95% confidence. These results and a site-by-site examination of the sequences clearly suggest that convergent or parallel evolution has occurred in ND1 between primates and rodents and/or between ferungulates and the outgroup. While the primate/rodent grouping is strange, ND1 also throws some interesting light on the relationships of some eutherian orders, marsupials, and montremes. In these parts of the tree, ND1 shows no apparent tendency for unexplained convergences. Received: 5 December 1997 / Accepted: 24 February 1998     

Acclimatization of Pink Salmon Oncorhynchus gorbuscha Walbaum in the European North: mtDNA Restriction Data

Pink salmon spawners introduced into the White Sea basin (the Umba River) were compared to the spawners from the basin of the Sea of Okhotsk (the Ola River) using restriction analysis of two fragments of mitochondrial DNA (mtDNA). One of the fragments included genes ND5/ND6, the other, the cytochrome b gene and the D-loop. It was found that mtDNA variation and diversity at the earlier examined nuclear allozyme genes significantly decreased in the odd broodline of pink salmon 8 years after the introduction. The haplotype diversity in the even broodline was considerably lower than in the odd broodline exhibiting virtually no change two generations after the introduction. Based on the results obtained, a possible role of these changes in adaptation of White Sea pink salmon from the odd broodline to the new environment is discussed.     

Acclimatization of salmon Oncorhynchus gorbuscha (Walbaum) in the European north: data from restriction analysis of mtDNA

Pink salmon spawners introduced into the White Sea basin (the Umba River) were compared to the spawners from the basin of the Sea of Okhotsk (the Ola River) using restriction analysis of two fragments of mitochondrial DNA (mtDNA). One of the fragments included genes ND5/ND6, the other, the cytochrome b gene and the control region. It was found that mtDNA variation and diversity at the earlier examined nuclear allozyme genes significantly decreased in the odd broodline of pink salmon 8 years after the introduction. The haplotype diversity in the even broodline was considerably lower than in the odd broodline exhibiting virtually no change two generations after the introduction. Based on the results obtained, a possible role of these changes in adaptation of White Sea pink salmon from the odd broodline to the new environment is discussed.     

Anadromous species and the new international law of the sea

Abstract

This article examines the evolution of international law relating to anadromous species, focusing exclusively on salmon and primarily considering the period since World War II. The discussion concerns major international harvesting in the western North Pacific, eastern North Pacific, and the North Atlantic. Unilateral actions are also described. Special attention is given to the relevant articles of the 1982 Convention on the Law of the Sea and to developments since its conclusion. Evidence is assessed for considering that the customary international law of the sea now recognizes the authority of the state of origin of salmon to prohibit high seas harvesting of salmon.     

Natural selection promotes divergence of transferrin among salmonid species

Transferrin is an iron-binding protein that plays an important role in iron metabolism and resistance to bacterial infection in a variety of organisms. A comparison of transferrin coding sequences from four salmonid species shows that the rate of evolution at nonsynonymous sites is significantly higher than the rate at synonymous sites, suggesting that positive natural selection for new alleles has played an important role in the evolution of transferrin in some salmon species. We hypothesize that the selective agent driving rapid divergence is interactions between host transferrin and the iron-scavenging proteins of pathogenic bacteria.     

Mitochondrial sequence evolution in spiders: intraspecific variation in tRNAs lacking the TPsiC Arm

Analyses of mitochondrial DNA sequences from three species of Habronattus jumping spiders (Chelicerata: Arachnida: Araneae) reveal unusual inferred tRNA secondary structures and gene arrangements, providing new information on tRNA evolution within chelicerate arthropods. Sequences from the protein-coding genes NADH dehydrogenase subunit 1 (ND1), cytochrome oxidase subunit I (COI), and subunit II (COII) were obtained, along with tRNA, tRNA, and large-subunit ribosomal RNA (16S) sequences; these revealed several peculiar features. First, inferred secondary structures of tRNA and, likely, tRNA, lack the TPsiC arm and the variable arm and therefore do not form standard cloverleaf structures. In place of these arms is a 5-6-nt T arm-variable loop (TV) replacement loop such as that originally described from nematode mitochondrial tRNAs. Intraspecific variation occurs in the acceptor stem sequences in both tRNAs. Second, while the proposed secondary structure of the 3' end of 16S is similar to that reported for insects, the sequence at the 5' end is extremely divergent, and the entire gene is truncated about 300 nt with respect to Drosophila yakuba. Third, initiation codons appear to consist of ATY (ATT and ATC) and TTG for ND1 and COII, respectively. Finally, Habronattus shares the same ND1-tRNA-16S gene arrangement as insects and crustaceans, thus illustrating variation in a tRNA gene arrangement previously proposed as a character distinguishing chelicerates from insects and crustaceans.     

Absence of restriction site variation in the mitochondrial ND5 and ND6 genes of Atlantic salmon amplified by the polymerase chain reaction

There was no restriction site variation in Atlantic salmon mitochondrial DNA in an amplified DNA segment beginning in the ND5 gene and terminating in the ND6 gene, although seven haplotypes were observed when the complete molecule was restricted.     

Migratory costs and the evolution of egg size and number in introduced and indigenous salmon populations

The trade-off between reproductive investment and migration should be an important factor shaping the evolution of life-history traits among populations following their radiation into habitats with different migratory costs and benefits. An experimentally induced difference in migratory rigor for families of chinook salmon (Oncorhynchus tshawytscha), of approximately 86 km and 413 m elevation, exacted a cost to somatic energy reserves (approximately 17% reduction in metabolizable mass) and ovarian investment (13.7% reduction in ovarian mass). This cost was associated with a reduction in egg size and paralleled the phenotypic pattern of divergence between two introduced New Zealand populations of common origin, presently breeding at sites with different migration distances. The genetic pattern of divergence of these same populations, detected under common rearing, was consistent with compensation for migratory costs (the population that migrates farther invested more in ovarian mass), but egg number more than egg size was associated with this evolution. These evolutionary patterns are consistent with what is known of the inheritance of these traits and with trade-offs and constraints favoring initial evolution in offspring number over offspring size. Analysis of egg number-size patterns of other Pacific salmon populations in their native range supported the hypothesis that migration strongly influences patterns of reproductive allocation, favoring a higher ratio of egg number to egg size with greater migration distance.     

Trophic ecology of Pacific salmon (<Emphasis Type="Italic">Oncorhynchus</Emphasis> spp.) in the ocean: a synthesis of stable isotope research

Increasing interest in the marine trophic dynamics of Pacific salmon has been motivated by the recognition of their sensitivity to changing climate and to the competitive effects of hatchery fish on wild stocks. It has become more common to use stable isotopes to supplement traditional diet studies of salmon in the ocean; however, there have been no integrated syntheses of these data to determine whether stable isotope analyses support the existing conventional wisdom of feeding strategies of the Pacific salmon. We performed a meta-analysis of stable isotope data to examine the extent of trophic partitioning among five species of Pacific salmon during their marine lives. Pink, sockeye, and chum salmon showed very high overlap in resource use and there was no consistent evidence for chum relying on alternative food webs dominated by gelatinous zooplankton. δ¹⁵N showed that Chinook and coho salmon fed at trophic levels higher than the other three species. In addition, these two species were distinctly enriched in ¹³C, suggesting more extensive use of coastal food webs compared to the more depleted (pelagic) signatures of pink, sockeye, and chum salmon. This paper presents the first synthesis of stable isotope work on Pacific salmon and provides δ¹⁵N and δ¹³C values applicable to research on the fate of the marine derived nutrients these organisms transport to freshwater and riparian ecosystems.     

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

Impact of mutations affecting ND mitochondria-encoded subunits on the activity and assembly of complex I in Chlamydomonas. Implication for the structural organization of the enzyme

The mitochondrial rotenone-sensitive NADH:ubiquinone oxidoreductase (complex I) comprises more than 35 subunits, the majority of which are encoded by the nucleus. In Chlamydomonas reinhardtii, only five components (ND1, ND2, ND4, ND5 and ND6) are coded for by the mitochondrial genome. Here, we characterize two mitochondrial mutants (dum5 and dum17) showing strong reduction or inactivation of complex I activity: dum5 is a 1T deletion in the 3' UTR of nd5 whereas dum17 is a 1T deletion in the coding sequence of nd6. The impact of these mutations and of mutations affecting nd1, nd4 and nd4/nd5 genes on the assembly of complex I is investigated. After separation of the respiratory complexes by blue native (BN)-PAGE or sucrose gradient centrifugation, we demonstrate that the absence of intact ND1 or ND6 subunit prevents the assembly of the 850 kDa whole complex, whereas the loss of ND4 or ND4/ND5 leads to the formation of a subcomplex of 650 kDa present in reduced amount. The implications of our findings for the possible role of these ND subunits on the activity of complex I and for the structural organization of the membrane arm of the enzyme are discussed. In mitochondria from all the strains analyzed, we moreover detected a 160-210 kDa fragment comprising the hydrophilic 49 kDa and 76 kDa subunits of the complex I peripheral arm and showing NADH dehydrogenase activity.     

Adaptive Evolution of cry Genes in Bacillus thuringiensis： Implications for Their Specificity Determination

The cry gene family, produced during the late exponential phase of growth in Bacillus thuringiensis, is a large, still-growing family of homologous genes, in which each gene encodes a protein with strong specific activity against only one or a few insect species. Extensive studies are mostly focusing on the structural and functional relationships of Cry proteins, and have revealed several residues or domains that are important for the target recognition and receptor attachment. In this study, we have employed a maximum likelihood method to detect evidence of adaptive evolution in Cry proteins, and have identified 24 positively selected residues, which are all located in Domain Ⅱ or Ⅲ. Combined with known data from mutagenesis studies, the majority of these residues, at the molecular level, contribute much to the insect specificity determination. We postulate that the potential pressures driving the diversification of Cry proteins may be in an attempt to adapt for the ＂arm race＂ between δ-endotoxins and the targeted insects, or to enlarge their target spectra, hence result in the functional divergence. The sites identified to be under positive selection would provide targets for further structural and functional analyses on Cry proteins.