Patterns of natural selection acting on the mitochondrial genome of a locally adapted fish species

Differential Mitochondrial Genome Expression of Four Hylid Frog Species under Low-Temperature Stress and Its Relationship with Amphibian Temperature Adaptation

Extreme weather poses huge challenges for animals that must adapt to wide variations in environmental temperature and, in many cases, it can lead to the local extirpation of populations or even the extinction of an entire species. Previous studies have found that one element of amphibian adaptation to environmental stress involves changes in mitochondrial gene expression at low temperatures. However, to date, comparative studies of gene expression in organisms living at extreme temperatures have focused mainly on nuclear genes. This study sequenced the complete mitochondrial genomes of five Asian hylid frog species: Dryophytes japonicus, D. immaculata, Hyla annectans, H. chinensis and H. zhaopingensis. It compared the phylogenetic relationships within the Hylidae family and explored the association between mitochondrial gene expression and evolutionary adaptations to cold stress. The present results showed that in D. immaculata, transcript levels of 12 out of 13 mitochondria genes were significantly reduced under cold exposure (p < 0.05); hence, we put forward the conjecture that D. immaculata adapts by entering a hibernation state at low temperature. In H. annectans, the transcripts of 10 genes (ND1, ND2, ND3, ND4, ND4L, ND5, ND6, COX1, COX2 and ATP8) were significantly reduced in response to cold exposure, and five mitochondrial genes in H. chinensis (ND1, ND2, ND3, ND4L and ATP6) also showed significantly reduced expression and transcript levels under cold conditions. By contrast, transcript levels of ND2 and ATP6 in H. zhaopingensis were significantly increased at low temperatures, possibly related to the narrow distribution of this species primarily at low latitudes. Indeed, H. zhaopingensis has little ability to adapt to low temperature (4 °C), or maybe to enter into hibernation, and it shows metabolic disorder in the cold. The present study demonstrates that the regulatory trend of mitochondrial gene expression in amphibians is correlated with their ability to adapt to variable climates in extreme environments. These results can predict which species are more likely to undergo extirpation or extinction with climate change and, thereby, provide new ideas for the study of species extinction in highly variable winter climates.

Evolutionary characteristics of the mitochondrial NADH dehydrogenase subunit 6 gene in some populations of four sympatric Mustela species (Mustelidae, Mammalia) from central Europe

BACKGROUND

Selection on or reticulate evolution of mtDNA is documented in various mammalian taxa and could lead to misleading phylogenetic conclusions if not recognized. We sequenced the MT-ND6 gene of four sympatric Mustelid species of the genus Mustela from some central European populations. We hypothesised positive selection on MT-ND6, given its functional importance and the different body sizes and life histories of the species, even though climatic differences may be unimportant for adaptation in sympatry.

METHODS AND RESULTS

MT-ND6 genes were sequenced in 187 sympatric specimens of weasels, Mustela nivalis, stoats, M. erminea, polecats, M. putorius, and steppe polecats, M. eversmannii, from eastern Austria and of fourteen allopatric polecats from eastern-central Germany. Median joining networks, neighbour joining and maximum likelihood analyses as well as Bayesian inference grouped all species according to earlier published phylogenetic models. However, polecats and steppe polecats, two very closely related species, shared the same two haplotypes. We found only negative selection within the Mustela sequences, including 131 downloaded ones covering thirteen species. Positive selection was observed on three MT-ND6 codons of other mustelid genera retrieved from GenBank.

CONCLUSIONS

Negative selection for MT-ND6 within the genus Mustela suggests absence of both environmental and species-specific effects of cellular energy metabolism despite large species-specific differences in body size. The presently found shared polymorphism in European polecats and steppe polecats may result from ancestral polymorphism before speciation and historical or recent introgressive hybridization; it may indicate mtDNA capture of steppe polecats by M. putorius in Europe.

Evidence for Selection on Mitochondrial OXPHOS Genes in the Mediterranean Killifish Aphanius fasciatus Valenciennes, 1821

Mitochondrial oxidative phosphorylation (OXPHOS) genes are a system subject to selection under determined environmental constraints despite a neutral evolution model that has long been hypothesized for the mitochondrial genome. In this study, the sequences of ND1, Cytb, and COI OXPHOS genes were analyzed in six populations of the eurythermal and euryhaline killifish A. fasciatus, to detect non-synonymous mutations leading to amino acid changes and to check whether selection acted on them using tests of recombination and selection. The results indicate a high COI and Cytb gene diversity and a high percentage of private haplotypes in all populations. In the Greek population, non-synonymous nucleotide substitutions were observed in the N-terminal region of COI and Cytb. Positively selected sites were also found. The information we obtained from the mitochondrial DNA sequences of A. fasciatus adds to the growing data on selective pressure acting on mitochondrial DNA in non-model species. These results should be explored from the perspective of the local adaptation of eurythermal and euryhaline species and supported using experimental evidence to better understand the interplay between historical climatic events and local adaptation and how each of them contributes to shaping the genetic structure of this species.

Ancient reindeer mitogenomes reveal island-hopping colonisation of the Arctic archipelagos

Climate warming at the end of the last glacial period had profound effects on the distribution of cold-adapted species. As their range shifted towards northern latitudes, they were able to colonise previously glaciated areas, including remote Arctic islands. However, there is still uncertainty about the routes and timing of colonisation. At the end of the last ice age, reindeer/caribou (Rangifer tarandus) expanded to the Holarctic region and colonised the archipelagos of Svalbard and Franz Josef Land. Earlier studies have proposed two possible colonisation routes, either from the Eurasian mainland or from Canada via Greenland. Here, we used 174 ancient, historical and modern mitogenomes to reconstruct the phylogeny of reindeer across its whole range and to infer the colonisation route of the Arctic islands. Our data shows a close affinity among Svalbard, Franz Josef Land and Novaya Zemlya reindeer. We also found tentative evidence for positive selection in the mitochondrial gene ND4, which is possibly associated with increased heat production. Our results thus support a colonisation of the Eurasian Arctic archipelagos from the Eurasian mainland and provide some insights into the evolutionary history and adaptation of the species to its High Arctic habitat.

Relaxed purifying selection pressure drives accelerated and dynamic gene rearrangements in thrips (Insecta: Thysanoptera) mitochondrial genomes

Insect mitochondrial genomes (mitogenome) generally present a typical gene order, which is considered as the ancestral arrangement. All sequenced mitogenomes in the Thysanoptera display high levels of gene rearrangement. Due to limited number of thrips mitogenomes sequenced, how gene rearrangement may be shaped by evolution remain unclear. Here, we analyzed 33 thrips mitogenomes, including 14 newly sequenced. These mitogenomes were diverse in organization, nucleotides substitution and gene arrangements. We found 28 highly rearranged gene orders with the breakpoints of gene rearrangements from 25 to 33. Reconstruction of the ancestors mitochondrial gene arrangements states indicated that Tubulifera have more complex pathways than Terebrantia in the gene order evolution. Molecular calibration estimated that divergence of two suborders occurred in the middle Triassic while the radiation of thrips was associated with the arose and flourish of angiosperm. Our evolutionary hypothesis testing suggests that relaxation of selection pressure enabled the early phase of Thysanoptera evolution, followed by a stronger selective pressure fixed diversification. Our analyses found gene inversion increases the nonsynonymous substitution rates and provide an evolutionary hypothesis driving the diverse gene orders.

Disproportionate role of nuclear-encoded proteins in organismal and mitochondrial thermal performance in a copepod

Determining the mechanisms by which organisms evolve thermal tolerance is crucial to predicting how populations may respond to changes in local temperature regimes. Although evidence of relationships between mitochondrial background and thermal adaptation have been found, the presence of both nuclear-encoded and mitochondrial DNA (mtDNA)-encoded proteins warrants experiments aimed at parsing out the relative role of each genome in thermal adaptation. We investigated the relative role of mtDNA-encoded products in thermal tolerance between two divergent populations of Tigriopus californicus using first-generation (F1) hybrids that vary in maternally inherited mtDNA but are heterozygous for population-specific alleles across nuclear loci. We tested two measures of thermal tolerance, (1) survivorship to acute thermal stress and (2) thermal stability of mitochondrial performance in Complex I-fueled ATP synthesis, both across a range of increasing temperatures. We found that the southern population (San Diego, CA, USA) outperformed the northern population (Strawberry Hill, OR, USA) in survivorship, and that both reciprocal F1 hybrid crosses had intermediate survival. Mitochondria from the San Diego population displayed greater stability in ATP synthesis with increasing temperatures compared with those from Strawberry Hill. Interestingly, hybrids from both cross directions had synthesis profiles that were very similar to that of Strawberry Hill. Taken together, these results suggest that the relative role of the mtDNA in these phenotypes is negligible compared with that of elements encoded by nuclear DNA in this system.

Cryopreservation effect on DNA methylation profile in rainbow trout spermatozoa

Spermatozoa are the cells that are most commonly used for cryopreservation of valuable genetic resources in aquaculture. It is known that fish spermatozoa transmit to the embryo not only their genetic but also their epigenetic profile, especially DNA methylation. Therefore, any alteration of the DNA methylation profile in spermatozoa induces the risk of transmitting epigenetic alterations to the offspring. The aim of this study was to assess the effect of cryopreservation on DNA methylation in rainbow trout spermatozoa. To trigger variable cellular response after freezing-thawing, spermatozoa from mature males were cryopreserved with dimethyl sulfoxide, methanol or glycerol as cryoprotectant. We observed that dimethyl sulfoxide was the best to preserve thawed spermatozoa functions. Methanol only slightly preserved all the cellular parameters, while glycerol failed to protect motility and fertilization ability. The consequences on DNA methylation were assessed using Reduced Representation Bisulfite Sequencing (RRBS). Sperm cryopreservation did not thoroughly impact DNA methylation, although 335-564 differentially methylated cytosines were characterized depending on the cryoprotectant. Very few of them were shared between cryoprotectants, and no correlation with the extent of cellular damage was found. Our study showed that DNA methylation was only slightly altered after sperm cryopreservation, and this may render further analysis of the risk for the progeny very challenging.

A Comprehensive Analysis of the Fowleria variegata (Valenciennes, 1832) Mitochondrial Genome and Its Phylogenetic Implications within the Family Apogonidae

Controversies surrounding the phylogenetic relationships within the family Apogonidae have persisted due to the limited molecular data, obscuring the evolution of these diverse tropical marine fishes. This study presents the first complete mitochondrial genome of Fowleria variegata, a previously unrepresented genus, using high-throughput Illumina sequencing. Through a comparative mitogenomic analysis, F. variegate was shown to exhibit a typical genome architecture and composition, including 13 protein-coding, 22 tRNA and 2 rRNA genes and a control region, consistent with studies of other Apogonidae species. Nearly all protein-coding genes started with ATG, while stop codons TAA/TAG/T were observed, along with evidence of strong functional constraints imposed via purifying selection. Phylogenetic reconstruction based on maximum likelihood and Bayesian approaches provided robust evidence that F. variegata forms a basal lineage closely related to P. trimaculatus within Apogonidae, offering novel perspectives into the molecular evolution of this family. By generating new mitogenomic resources and evolutionary insights, this study makes important headway in elucidating the phylogenetic relationships and mitogenomic characteristics of Apogonidae fishes. The findings provide critical groundwork for future investigations into the drivers of diversification, speciation patterns, and adaptive radiation underlying the extensive ecological diversity and biological success of these marine fishes using phylogenomics and population genomics approaches.

Solving the conundrum of intra-specific variation in metabolic rate: A multidisciplinary conceptual and methodological toolkit: New technical developments are opening the door to an understanding of why metabolic rate varies among individual animals of a species: New technical developments are opening the door to an understanding of why metabolic rate varies among individual animals of a species

Researchers from diverse disciplines, including organismal and cellular physiology, sports science, human nutrition, evolution and ecology, have sought to understand the causes and consequences of the surprising variation in metabolic rate found among and within individual animals of the same species. Research in this area has been hampered by differences in approach, terminology and methodology, and the context in which measurements are made. Recent advances provide important opportunities to identify and address the key questions in the field. By bringing together researchers from different areas of biology and biomedicine, we describe and evaluate these developments and the insights they could yield, highlighting the need for more standardisation across disciplines. We conclude with a list of important questions that can now be addressed by developing a common conceptual and methodological toolkit for studies on metabolic variation in animals.

Dusp1 regulates thermal tolerance limits in zebrafish by maintaining mitochondrial integrity

Temperature tolerance restricts the distribution of a species. However, the molecular and cellular mechanisms that set the thermal tolerance limits of an organism are poorly understood. Here, we report on the function of dual-specificity phosphatase 1 (DUSP1) in thermal tolerance regulation. Notably, we found that dusp1 ^-/- zebrafish grew normally but survived within a narrowed temperature range. The higher susceptibility of these mutant fish to both cold and heat challenges was attributed to accelerated cell death caused by aggravated mitochondrial dysfunction and over-production of reactive oxygen species in the gills. The DUSP1-MAPK-DRP1 axis was identified as a key pathway regulating these processes in both fish and human cells. These observations suggest that DUSP1 may play a role in maintaining mitochondrial integrity and redox homeostasis. We therefore propose that maintenance of cellular redox homeostasis may be a key mechanism for coping with cellular thermal stress and that the interplay between signaling pathways regulating redox homeostasis in the most thermosensitive tissue (i.e., gills) may play an important role in setting the thermal tolerance limit of zebrafish.

Signatures of positive selection in the mitochondrial genome of neotropical freshwater stingrays provide clues about the transition from saltwater to freshwater environment

Neotropical freshwater stingrays (subfamily Potamotrygoninae) are carnivorous bottom feeder batoids widely distributed in most river basins of South America. They represent the unique extant group of elasmobranchs that evolved to live exclusively in freshwater environments. These species are exploited either by commercial fisheries (e.g., for food or ornamental industry) or by indigenous communities allocated along with their natural range. Restrictive life history characteristics coupled with habitat degradation make Potamotrygoninae species highly vulnerable to human impacts and highlight the necessity of studies to inform basic biological aspects, from ecology to genetics, to guide their conservation and clarify the molecular basis of adaptation to the freshwater environment. We used available and newly assembled Potamotrygon spp. mitogenomes to perform a comparative investigation of their molecular evolution. A phylogenetic estimation using the mitogenome of Potamotrygon falkneri and other Elasmobranchii supports monophyly for Potamotrygonidae and indicates a close relationship to Dasyatidae. A synteny analysis comprising 3 Potamotrygon and other 51 batoids revealed a highly conserved mitogenomic context. We detected various amino acid sites under positive selection exclusively in Potamotrygon spp., within the sequences of ND4, ND5, ND6, and COXII genes. Positively selected mutational events in key genes of energetic metabolism may be related to the physiological adaptation of Potamotrygon spp. during the ancient incursion into freshwater. This broad comparative mitogenomic study provides novel insights into the evolutionary history of neotropical freshwater stingrays and their relatives and stands out as a valuable resource to aid in current and future research on elasmobranch molecular evolution.

Comparative Mitogenomics of Two Sympatric Catfishes of Exostoma (Siluriformes: Sisoridae) from the Lower Yarlung Tsangpo River and Its Application for Phylogenetic Consideration

The genus Exostoma is a group of stenotopic and rheophilic glyptosternine catfishes distributed in South and Southeast Asia. So far, comprehensive studies on mitogenomics referring to this genus are very scarce. In this study, we first sequenced and annotated the complete mitochondrial genomes of Exostoma tibetanum and Exostoma tenuicaudatum—two sympatric congeners from the lower Yarlung Tsangpo River, Tibet, China. The mitogenomes of both species contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one light-strand origin of replication, and one control region, with lengths of 16,528 bp and 16,533 bp, respectively. The mitogenome architecture, nucleotide composition, and codon usage of protein-coding genes were almost identical between the two Exostoma species, although some estimated parameters varied. Phylogenetic analysis strongly supported the monophyly of Exostoma in the subfamily Glyptosternae, and Exostoma tibetanum had the closest relationship to Exostoma tenuicaudatum. The divergence time estimation demonstrated that these two species diverged approximately 1.51 Ma during the early Pleistocene, which was speculated to be triggered by the river system changes caused by the uplift of the southeastern Tibetan Plateau. Selection pressure analyses indicated that all protein-coding genes of Exostoma species underwent a strong purifying selection, while minority positive sites from NADH dehydrogenase complex genes were detected. These findings are expected to promote our understanding of the molecular phylogeny of the genus Exostoma and provide valuable mitogenomic resources for the subfamily Glyptosternae

Mitogenome selection in the evolution of key ecological strategies in the ancient hexapod class Collembola

A longstanding question in evolutionary biology is how natural selection and environmental pressures shape the mitochondrial genomic architectures of organisms. Mitochondria play a pivotal role in cellular respiration and aerobic metabolism, making their genomes functionally highly constrained. Evaluating selective pressures on mitochondrial genes can provide functional and ecological insights into the evolution of organisms. Collembola (springtails) are an ancient hexapod group that includes the oldest terrestrial arthropods in the fossil record, and that are closely associated with soil environments. Of interest is the diversity of habitat stratification preferences (life forms) exhibited by different species within the group. To understand whether signals of positive selection are linked to the evolution of life forms, we analysed 32 published Collembola mitogenomes in a phylomitogenomic framework. We found no evidence that signatures of selection are correlated with the evolution of novel life forms, but rather that mutations have accumulated as a function of time. Our results highlight the importance of nuclear-mitochondrial interactions in the evolution of collembolan life forms and that mitochondrial genomic data should be interpreted with caution, as complex selection signals may complicate evolutionary inferences.

Whole-Genome Survey Analyses Provide a New Perspective for the Evolutionary Biology of Shimofuri Goby, Tridentiger bifasciatus

The shimofuri goby (Tridentiger bifasciatus) is a small and highly adaptable goby, distributed along the coasts of China, the Sea of Japan, and the west coastal and estuarine areas of the Northwest Pacific. Next-generation sequencing was used to generate genome-wide survey data to provide essential characterization of the shimofuri goby genome and for the further mining of genomic information. The genome size of the shimofuri goby was estimated to be approximately 887.60 Mb through K-mer analysis, with a heterozygosity ratio and repeat sequence ratio of 0.47% and 32.60%, respectively. The assembled genome was used to identify microsatellite motifs (Simple Sequence Repeats, SSRs), extract single-copy homologous genes and assemble the mitochondrial genome. A total of 288,730 SSRs were identified. The most frequent SSRs were dinucleotide repeats (with a frequency of 61.15%), followed by trinucleotide (29.87%), tetranucleotide (6.19%), pentanucleotide (1.13%), and hexanucleotide repeats (1.66%). The results of the phylogenetic analysis based on single-copy homologous genes showed that the shimofuri goby and Rhinogobius similis can be clustered into one branch. The shimofuri goby was originally thought to be the same as the chameleon goby (Tridentiger trigonocephalus) due to their close morphological resemblance. However, a complete mitochondrial genome was assembled and the results of the phylogenetic analysis support the inclusion of the shimofuri goby as a separate species. PSMC analysis indicated that the shimofuri goby experienced a bottleneck event during the Pleistocene Glacial Epoch, in which its population size decreased massively, and then it began to recover gradually after the Last Glacial Maximum. This study provides a reference for the further assembly of the complete genome map of the shimofuri goby, and is a valuable genomic resource for the study of its evolutionary biology.

Mitochondrial Genome Contributes to the Thermal Adaptation of the Oomycete Phytophthora infestans

As a vital element of climate change, elevated temperatures resulting from global warming present new challenges to natural and agricultural sustainability, such as ecological disease management. Mitochondria regulate the energy production of cells in responding to environmental fluctuation, but studying their contribution to the thermal adaptation of species is limited. This knowledge is needed to predict future disease epidemiology for ecology conservation and food security. Spatial distributions of the mitochondrial genome (mtDNA) in 405 Phytophthora infestans isolates originating from 15 locations were characterized. The contribution of MtDNA to thermal adaptation was evaluated by comparative analysis of mtDNA frequency and intrinsic growth rate, relative population differentiation in nuclear and mtDNA, and associations of mtDNA distribution with local geography climate conditions. Significant variation in frequency, intrinsic growth rate, and spatial distribution was detected in mtDNA. Population differentiation in mtDNA was significantly higher than that in the nuclear genome, and spatial distribution of mtDNA was strongly associated with local climatic conditions and geographic parameters, particularly air temperature, suggesting natural selection caused by a local temperature is the main driver of the adaptation. Dominant mtDNA grew faster than the less frequent mtDNA. Our results provide useful insights into the evolution of pathogens under global warming. Given its important role in biological functions and adaptation to local air temperature, mtDNA intervention has become an increasing necessity for future disease management. To secure ecological integrity and food production under global warming, a synergistic study on the interactive effect of changing temperature on various components of biological and ecological functions of mitochondria in an evolutionary frame is urgently needed.

Adaptive mitochondrial genome functioning in ecologically different farm-impacted natural seedbeds of the endemic blue mussel Mytilus chilensis

We assessed the adaptive contribution of the mitochondrial genes involved with the respiratory chain and oxidative phosphorylation of the blue mussel Mytilus chilensis, a native and heavily exploited species in the inner sea of Chiloé Island, southern Chile. The assembled mitochondrial transcriptome of individuals from two ecologically different farm-impacted natural seedbeds, Cochamó (41°S) and Yaldad (42°S), represented about 4.5% of the whole de novo transcriptome of the species and showed location and tissue (gills, mantle) specific expression differences in 13 protein-coding mitochondrial genes. The RNA-Seq analysis detected differences in the number of up-regulated mitogenes between individuals from Cochamó (7) and Yaldad (11), some being tissue-specific (ND4L and COX2). However, the analysis did not detect transcripts-per-million (TPM = 0) of ND2 and ND5 in gills and ATP6 in mantle samples from Cochamó. Likewise, for ND6 and ATP8 in any sample. Several monomorphic location-specific mitochondrial genetic variants were detected in samples from Cochamó (78) and Yaldad (207), representing standing genetic variability to optimize mitochondrial functioning under local habitats. Overall, these mitochondrial transcriptomic differences reflect the impact of environmental conditions on the mitochondrial genome functioning and offer new markers to assess the effects on mussel fitness of habitat translocations, a routine industry practice. Likewise, these mitochondrial markers should help monitor and maintain adaptive population differences in this keystone and heavily exploited native species.

Positive selection over the mitochondrial genome and its role in the diversification of gentoo penguins in response to adaptation in isolation

Although mitochondrial DNA has been widely used in phylogeography, evidence has emerged that factors such as climate, food availability, and environmental pressures that produce high levels of stress can exert a strong influence on mitochondrial genomes, to the point of promoting the persistence of certain genotypes in order to compensate for the metabolic requirements of the local environment. As recently discovered, the gentoo penguins (Pygoscelis papua) comprise four highly divergent lineages across their distribution spanning the Antarctic and sub-Antarctic regions. Gentoo penguins therefore represent a suitable animal model to study adaptive processes across divergent environments. Based on 62 mitogenomes that we obtained from nine locations spanning all four gentoo penguin lineages, we demonstrated lineage-specific nucleotide substitutions for various genes, but only lineage-specific amino acid replacements for the ND1 and ND5 protein-coding genes. Purifying selection (dN/dS < 1) is the main driving force in the protein-coding genes that shape the diversity of mitogenomes in gentoo penguins. Positive selection (dN/dS > 1) was mostly present in codons of the Complex I (NADH genes), supported by two different codon-based methods at the ND1 and ND4 in the most divergent lineages, the eastern gentoo penguin from Crozet and Marion Islands and the southern gentoo penguin from Antarctica respectively. Additionally, ND5 and ATP6 were under selection in the branches of the phylogeny involving all gentoo penguins except the eastern lineage. Our study suggests that local adaptation of gentoo penguins has emerged as a response to environmental variability promoting the fixation of mitochondrial haplotypes in a non-random manner. Mitogenome adaptation is thus likely to have been associated with gentoo penguin diversification across the Southern Ocean and to have promoted their survival in extreme environments such as Antarctica. Such selective processes on the mitochondrial genome may also be responsible for the discordance detected between nuclear- and mitochondrial-based phylogenies of gentoo penguin lineages.

Comparative mitogenomics of Clupeoid fish provides insights into the adaptive evolution of mitochondrial oxidative phosphorylation (OXPHOS) genes and codon usage in the heterogeneous habitats

Clupeoid fish can be considered excellent candidates to understand the role of mitochondrial DNA in adaptive evolution, as they have colonized different habitats (marine, brackish, freshwater, tropical and temperate regions) over millions of years. Here, we investigate patterns of tRNA location, codon usage bias, and lineage-specific diversifying selection signals to provide novel insights into how evolutionary improvements of mitochondrial metabolic efficiency have allowed clupeids to adapt to different habitats. Based on whole mitogenome data of 70 Clupeoids with a global distribution we find that purifying selection was the dominant force acting and that the mutational deamination pressure in mtDNA was stronger than the codon/amino acid constraints. The codon usage pattern appears evolved to achieve high translational efficiency (codon/amino acid-related constraints), as indicated by the complementarity of most codons to the GT-saturated tRNA anticodon sites (retained by deamination-induced pressure) and usage of the codons of the tRNA genes situated near to the control region (fixed by deamination pressure) where transcription efficiency was high. The observed shift in codon preference patterns between marine and euryhaline/freshwater Clupeoids indicates possible selection for improved translational efficiency in mitochondrial genes while adapting to low-salinity habitats. This mitogenomic plasticity and enhanced efficiency of the metabolic machinery may have contributed to the evolutionary success and abundance of Clupeoid fish.

Temperature-Specific and Sex-Specific Fitness Effects of Sympatric Mitochondrial and Mito-Nuclear Variation in Drosophila obscura

Simple Summary

Does variation in the mitochondrial DNA sequence influence the survival and reproduction of an individual? What is the purpose of genetic variation of the mitochondrial DNA between individuals from the same population? As a simple laboratory model, Drosophila species can give us the answer to this question. Creating experimental lines with different combinations of mitochondrial and nuclear genomic DNA and testing how successful these lines were in surviving in different experimental set-ups enables us to deduce the effect that both genomes have on fitness. This study on D. obscura experimentally validates theoretical models that explain the persistence of mitochondrial DNA variation within populations. Our results shed light on the various mechanisms that maintain this type of variation. Finally, by conducting the experiments on two experimental temperatures, we have shown that environmental variations can support mitochondrial DNA variation within populations.

Abstract

The adaptive significance of sympatric mitochondrial (mtDNA) variation and the role of selective mechanisms that maintain it are debated to this day. Isofemale lines of Drosophila obscura collected from four populations were backcrossed within populations to construct experimental lines, with all combinations of mtDNA Cyt b haplotypes and nuclear genetic backgrounds (nuDNA). Individuals of both sexes from these lines were then subjected to four fitness assays (desiccation resistance, developmental time, egg-to-adult viability and sex ratio) on two experimental temperatures to examine the role of temperature fluctuations and sex-specific selection, as well as the part that interactions between the two genomes play in shaping mtDNA variation. The results varied across populations and fitness components. In the majority of comparisons, they show that sympatric mitochondrial variants affect fitness. However, their effect should be examined in light of interactions with nuDNA, as mito-nuclear genotype was even more influential on fitness across all components. We found both sex-specific and temperature-specific differences in mitochondrial and mito-nuclear genotype ranks in all fitness components. The effect of temperature-specific selection was found to be more prominent, especially in desiccation resistance. From the results of different components tested, we can also infer that temperature-specific mito-nuclear interactions rather than sex-specific selection on mito-nuclear genotypes have a more substantial role in preserving mtDNA variation in this model species.

Complete mitochondrial genomes reveal robust phylogenetic signals and evidence of positive selection in horseshoe bats

BACKGROUND

In genus Rhinolophus, species in the Rhinolophus philippinensis and R. macrotis groups are unique because the horseshoe bats in these group have relatively low echolocation frequencies and flight speeds compared with other horseshoe bats with similar body size. The different characteristics among bat species suggest particular evolutionary processes may have occurred in this genus. To study the adaptive evidence in the mitochondrial genomes (mitogenomes) of rhinolophids, especially the mitogenomes of the species with low echolocation frequencies, we sequenced eight mitogenomes and used them for comparative studies of molecular phylogeny and adaptive evolution.

RESULTS

Phylogenetic analysis using whole mitogenome sequences produced robust results and provided phylogenetic signals that were better than those obtained using single genes. The results supported the recent establishment of the separate macrotis group. The signals of adaptive evolution discovered in the Rhinolophus species were tested for some of the codons in two genes (ND2 and ND6) that encode NADH dehydrogenases in oxidative phosphorylation system complex I. These genes have a background of widespread purifying selection. Signals of relaxed purifying selection and positive selection were found in ND2 and ND6, respectively, based on codon models and physicochemical profiles of amino acid replacements. However, no pronounced overlap was found for non-synonymous sites in the mitogenomes of all the species with low echolocation frequencies. A signal of positive selection for ND5 was found in the branch-site model when R. philippinensis was set as the foreground branch.

CONCLUSIONS

The mitogenomes provided robust phylogenetic signals that were much more informative than the signals obtained using single mitochondrial genes. Two mitochondrial genes that encoding proteins in the oxidative phosphorylation system showed some evidence of adaptive evolution in genus Rhinolophus and the positive selection signals were tested for ND5 in R. philippinensis. These results indicate that mitochondrial protein-coding genes were targets of adaptive evolution during the evolution of Rhinolophus species, which might have contributed to a diverse range of acoustic adaptations in this genus.

Contrasting Host-Parasite Population Structure: Morphology and Mitogenomics of a Parasitic Flatworm on Pelagic Deepwater Cichlid Fishes from Lake Tanganyika

Little phylogeographic structure is presumed for highly mobile species in pelagic zones. Lake Tanganyika is a unique ecosystem with a speciose and largely endemic fauna famous for its remarkable evolutionary history. In bathybatine cichlid fishes, the pattern of lake-wide population differentiation differs among species. We assessed the congruence between the phylogeographic structure of bathybatine cichlids and their parasitic flatworm Cichlidogyrus casuarinus to test the magnifying glass hypothesis. Additionally, we evaluated the use of a PoolSeq approach to study intraspecific variation in dactylogyrid monogeneans. The lake-wide population structure of C. casuarinus ex Hemibates stenosoma was assessed based on a portion of the cox1 gene combined with morphological characterisation. Additionally, intraspecific mitogenomic variation among 80 parasite samples from one spatially constrained metapopulation was assessed using shotgun NGS. While no clear geographic genetic structure was detected in parasites, both geographic and host-related phenotypic variation was apparent. The incongruence with the genetic north-south gradient observed in H. stenosoma may be explained by the broad host range of this flatworm including eupelagic bathybatine host species that form panmictic populations across the lake. In addition, we present the first parasite mitogenome from Lake Tanganyika and propose a methodological framework for studying the intraspecific mitogenomic variation of dactylogyrid monogeneans.

Convergent Adaptation in Mitochondria of Phylogenetically Distant Birds: Does it Exist?

In a wide range of taxa, proteins encoded by mitochondrial genomes are involved in adaptation to lifestyle that requires oxygen starvation or elevation of metabolism rate. It remains poorly understood to what extent adaptation to similar conditions is associated with parallel changes in these proteins. We search for a genetic signal of parallel or convergent evolution in recurrent molecular adaptation to high altitude, migration, diving, wintering, unusual flight abilities, or loss of flight in mitochondrial genomes of birds. Developing on previous work, we design an approach for the detection of recurrent coincident changes in genotype and phenotype, indicative of an association between the two. We describe a number of candidate sites involved in recurrent adaptation in ND genes. However, we find that the majority of convergence events can be explained by random coincidences without invoking adaptation.

The evolutionary history of manatees told by their mitogenomes

The manatee family encompasses three extant congeneric species: Trichechus senegalensis (African manatee), T. inunguis (Amazonian manatee), and T. manatus (West Indian manatee). The fossil record for manatees is scant, and few phylogenetic studies have focused on their evolutionary history. We use full mitogenomes of all extant manatee species to infer the divergence dates and biogeographical histories of these species and the effect of natural selection on their mitogenomes. The complete mitochondrial genomes of T. inunguis (16,851 bp), T. senegalensis (16,882 bp), and T. manatus (16,882 bp), comprise 13 protein-coding genes, 2 ribosomal RNA genes (rRNA - 12S and 16S), and 22 transfer RNA genes (tRNA), and (D-loop/CR). Our analyses show that the first split within Trichechus occurred during the Late Miocene (posterior mean 6.56 Ma and 95% HPD 3.81–10.66 Ma), followed by a diversification event in the Plio-Pleistocene (posterior mean 1.34 Ma, 95% HPD 0.1–4.23) in the clade composed by T. inunguis and T. manatus; T. senegalensis is the sister group of this clade with higher support values (pp > 0.90). The branch-site test identified positive selection on T. inunguis in the 181st position of the ND4 amino acid gene (LRT = 6.06, p = 0.0069, BEB posterior probability = 0.96). The ND4 gene encodes one subunit of the NADH dehydrogenase complex, part of the oxidative phosphorylation machinery. In conclusion, our results provide novel insight into the evolutionary history of the Trichechidae during the Late Miocene, which was influenced by geological events, such as Amazon Basin formation.

Targeted Sequencing of Mitochondrial Genes Reveals Signatures of Molecular Adaptation in a Nearly Panmictic Small Pelagic Fish Species

Ongoing climatic changes, with predictable impacts on marine environmental conditions, are expected to trigger organismal responses. Recent evidence shows that, in some marine species, variation in mitochondrial genes involved in the aerobic conversion of oxygen into ATP at the cellular level correlate with gradients of sea surface temperature and gradients of dissolved oxygen. Here, we investigated the adaptive potential of the European sardine Sardina pilchardus populations offshore the Iberian Peninsula. We performed a seascape genetics approach that consisted of the high throughput sequencing of mitochondria’s ATP6, COI, CYTB and ND5 and five microsatellite loci on 96 individuals coupled with environmental information on sea surface temperature and dissolved oxygen across five sampling locations. Results show that, despite sardines forming a nearly panmictic population around Iberian Peninsula, haplotype frequency distribution can be explained by gradients of minimum sea surface temperature and dissolved oxygen. We further identified that the frequencies of the most common CYTB and ATP6 haplotypes negatively correlate with minimum sea surface temperature across the sampled area, suggestive of a signature of selection. With signatures of selection superimposed on highly connected populations, sardines may be able to follow environmental optima and shift their distribution northwards as a response to the increasing sea surface temperatures.

Integrative genomic phylogeography reveals signs of mitonuclear incompatibility in a natural hybrid goby population

Hybridization between divergent lineages generates new allelic combinations. One mechanism that can hinder the formation of hybrid populations is mitonuclear incompatibility, that is, dysfunctional interactions between proteins encoded in the nuclear and mitochondrial genomes (mitogenomes) of diverged lineages. Theoretically, selective pressure due to mitonuclear incompatibility can affect genotypes in a hybrid population in which nuclear genomes and mitogenomes from divergent lineages admix. To directly and thoroughly observe this key process, we de novo sequenced the 747-Mb genome of the coastal goby, Chaenogobius annularis, and investigated its integrative genomic phylogeographics using RNA-sequencing, RAD-sequencing, genome resequencing, whole mitogenome sequencing, amplicon sequencing, and small RNA-sequencing. Chaenogobius annularis populations have been geographically separated into Pacific Ocean (PO) and Sea of Japan (SJ) lineages by past isolation events around the Japanese archipelago. Despite the divergence history and potential mitonuclear incompatibility between these lineages, the mitogenomes of the PO and SJ lineages have coexisted for generations in a hybrid population on the Sanriku Coast. Our analyses revealed accumulation of nonsynonymous substitutions in the PO-lineage mitogenomes, including two convergent substitutions, as well as signals of mitochondrial lineage-specific selection on mitochondria-related nuclear genes. Finally, our data implied that a microRNA gene was involved in resolving mitonuclear incompatibility. Our integrative genomic phylogeographic approach revealed that mitonuclear incompatibility can affect genome evolution in a natural hybrid population.

Mitochondria and the thermal limits of ectotherms

Temperature is a critical abiotic factor shaping the distribution and abundance of species, but the mechanisms that underpin organismal thermal limits remain poorly understood. One possible mechanism underlying these limits is the failure of mitochondrial processes, as mitochondria play a crucial role in animals as the primary site of ATP production. Conventional measures of mitochondrial performance suggest that these organelles can function at temperatures much higher than those that limit whole-organism function, suggesting that they are unlikely to set organismal thermal limits. However, this conclusion is challenged by recent data connecting sequence variation in mitochondrial genes to whole-organism thermal tolerance. Here, we review the current state of knowledge of mitochondrial responses to thermal extremes and ask whether they are consistent with a role for mitochondrial function in shaping whole-organism thermal limits. The available data are fragmentary, but it is possible to draw some conclusions. There is little evidence that failure of maximal mitochondrial oxidative capacity as assessed in vitro sets thermal limits, but there is some evidence to suggest that temperature effects on ATP synthetic capacity may be important. Several studies suggest that loss of mitochondrial coupling is associated with the thermal limits for organismal growth, although this needs to be rigorously tested. Most studies have utilized isolated mitochondrial preparations to assess the effects of temperature on these organelles, and there remain many untapped opportunities to address these questions using preparations that retain more of their biological context to better connect these subcellular processes with whole-organism thermal limits.

An experimental test of temperature-dependent selection on mitochondrial haplotypes in Callosobruchus maculatus seed beetles

Mitochondrial DNA (mtDNA) consists of few but vital maternally inherited genes that interact closely with nuclear genes to produce cellular energy. How important mtDNA polymorphism is for adaptation is still unclear. The assumption in population genetic studies is often that segregating mtDNA variation is selectively neutral. This contrasts with empirical observations of mtDNA haplotypes affecting fitness-related traits and thermal sensitivity, and latitudinal clines in mtDNA haplotype frequencies. Here, we experimentally test whether ambient temperature affects selection on mtDNA variation, and whether such thermal effects are influenced by intergenomic epistasis due to interactions between mitochondrial and nuclear genes, using replicated experimental evolution in Callosobruchus maculatus seed beetle populations seeded with a mixture of different mtDNA haplotypes. We also test for sex-specific consequences of mtDNA evolution on reproductive success, given that mtDNA mutations can have sexually antagonistic fitness effects. Our results demonstrate natural selection on mtDNA haplotypes, with some support for thermal environment influencing mtDNA evolution through mitonuclear epistasis. The changes in male and female reproductive fitness were both aligned with changes in mtDNA haplotype frequencies, suggesting that natural selection on mtDNA is sexually concordant in stressful thermal environments. We discuss the implications of our findings for the evolution of mtDNA.

An evolutionary, or "Mitocentric" perspective on cellular function and disease

The incidence of common, metabolic diseases (e.g. obesity, cardiovascular disease, diabetes) with complex genetic etiology has been steadily increasing nationally and globally. While identification of a genetic model that explains susceptibility and risk for these diseases has been pursued over several decades, no clear paradigm has yet been found to disentangle the genetic basis of polygenic/complex disease development. Since the evolution of the eukaryotic cell involved a symbiotic interaction between the antecedents of the mitochondrion and nucleus (which itself is a genetic hybrid), we suggest that this history provides a rational basis for investigating whether genetic interaction and co-evolution of these genomes still exists. We propose that both mitochondrial and Mendelian, or "mito-Mendelian" genetics play a significant role in cell function, and thus disease risk. This paradigm contemplates the natural variation and co-evolution of both mitochondrial and nuclear DNA backgrounds on multiple mitochondrial functions that are discussed herein, including energy production, cell signaling and immune response, which collectively can influence disease development. At the nexus of these processes is the economy of mitochondrial metabolism, programmed by both mitochondrial and nuclear genomes.

Signals of selection in the mitogenome provide insights into adaptation mechanisms in heterogeneous habitats in a widely distributed pelagic fish

Oceans are vast, dynamic, and complex ecosystems characterized by fluctuations in environmental parameters like sea surface temperature (SST), salinity, oxygen availability, and productivity. Environmental variability acts as the driver of organismal evolution and speciation as organisms strive to cope with the challenges. We investigated the evolutionary consequences of heterogeneous environmental conditions on the mitogenome of a widely distributed small pelagic fish of Indian ocean, Indian oil sardine, Sardinella longiceps. Sardines were collected from different eco-regions of the Indian Ocean and selection patterns analyzed in coding and non-coding regions. Signals of diversifying selection were observed in key functional regions involved in OXPHOS indicating OXPHOS gene regulation as the critical factor to meet enhanced energetic demands. A characteristic control region with 38–40 bp tandem repeat units under strong selective pressure as evidenced by sequence conservation and low free energy values was also observed. These changes were prevalent in fishes from the South Eastern Arabian Sea (SEAS) followed by the Northern Arabian Sea (NAS) and rare in Bay of Bengal (BoB) populations. Fishes belonging to SEAS exhibited accelerated substitution rate mainly due to the selective pressures to survive in a highly variable oceanic environment characterized by seasonal hypoxia, variable SST, and food availability.

Population-level inferences from environmental DNA-Current status and future perspectives

Environmental DNA (eDNA) extracted from water samples has recently shown potential as a valuable source of population genetic information for aquatic macroorganisms. This approach offers several potential advantages compared with conventional tissue-based methods, including the fact that eDNA sampling is noninvasive and generally more cost-efficient. Currently, eDNA approaches have been limited to single-marker studies of mitochondrial DNA (mtDNA), and the relationship between eDNA haplotype composition and true haplotype composition still needs to be thoroughly verified. This will require testing of bioinformatic and statistical software to correct for erroneous sequences, as well as biases and random variation in relative sequence abundances. However, eDNA-based population genetic methods have far-reaching potential for both basic and applied research. In this paper, we present a brief overview of the achievements of eDNA-based population genetics to date, and outline the prospects for future developments in the field, including the estimation of nuclear DNA (nuDNA) variation and epigenetic information. We discuss the challenges associated with eDNA samples as opposed to those of individual tissue samples and assess whether eDNA might offer additional types of information unobtainable with tissue samples. Lastly, we provide recommendations for determining whether an eDNA approach would be a useful and suitable choice in different research settings. We limit our discussion largely to contemporary aquatic systems, but the advantages, challenges, and perspectives can to a large degree be generalized to eDNA studies with a different spatial and temporal focus.

Positive selection and precipitation effects on the mitochondrial NADH dehydrogenase subunit 6 gene in brown hares (Lepus europaeus) under a phylogeographic perspective

Previous studies in hares and jackrabbits have indicated that positive selection has shaped the genetic diversity of mitochondrial genes involved in oxidative phosphorylation, which may affect cellular energy production and cause regional adaptation to different environmental (climatic) pressures. In the present study, we sequenced the NADH dehydrogenase subunit 6 (MT-ND6) gene of 267 brown hares (L. europaeus) from Europe and Asia Minor and tested for positive selection and adaptations acting on amino acid sequences (protein variants). Molecular diversity indices and spatial clustering were assessed by DnaSP, Network, and Geneland, while the presence of selection signals was tested by codeml in PAML, and by using the Datamonkey Adaptive Evolution web server. The SPSS software was used to run multinomial regression models to test for possible effects of climate parameters on the currently obtained protein variants. Fifty-eight haplotypes were revealed with a haplotype diversity of 0.817, coding for 17 different protein variants. The MT-ND6 phylogeographic pattern as determined by the nucleotide sequences followed the earlier found model based on the neutrally evolving D-loop sequences, and reflected the earlier found phylogeographic Late Pleistocene scenario. Based on several selection tests, only one codon position consistently proved to be under positive selection. It did occur exclusively in the evolutionarily younger hares from Europe and it gave rise to several protein variants from the southeastern and south-central Balkans. The occurrence of several of those variants was significantly favored under certain precipitation conditions, as proved by our multinomial regression models. Possibly, the great altitudinal variation in the Balkans may have lead to bigger changes in precipitation across that region and this may have imposed an evolutionarily novel selective pressure on the protein variants and could have led to regional adaptation.

Preliminary assessment of adaptive evolution of mitochondrial protein coding genes in darters (Percidae: Etheostomatinae)

Background: Mitochondrial DNA of vertebrates contains genes for 13 proteins involved in oxidative phosphorylation. Some of these genes have been shown to undergo adaptive evolution in a variety of species. This study examines all mitochondrial protein coding genes in 11 darter species to determine if any of these genes show evidence of positive selection. Methods: The mitogenome from four darter was sequenced and annotated. Mitogenome sequences for another seven species were obtained from GenBank. Alignments of each of the protein coding genes were subject to codon-based identification of positive selection by Selecton, MEME and FEL. Results: Evidence of positive selection was obtained for six of the genes by at least one of the methods. CYTB was identified as having evolved under positive selection by all three methods at the same codon location. Conclusions: Given the evidence for positive selection of mitochondrial protein coding genes in darters, a more extensive analysis of mitochondrial gene evolution in all the extant darter species is warranted.

Signals of positive selection in mitochondrial protein-coding genes of woolly mammoth: Adaptation to extreme environments?

The mammoths originated in warm and equatorial Africa and later colonized cold and high-latitude environments. Studies on nuclear genes suggest that woolly mammoth had evolved genetic variations involved in processes relevant to cold tolerance, including lipid metabolism and thermogenesis, and adaptation to extremely varied light and darkness cycles. The mitochondria is a major regulator of cellular energy metabolism, thus the mitogenome of mammoths may also exhibit adaptive evolution. However, little is yet known in this regard. In this study, we analyzed mitochondrial protein-coding genes (MPCGs) sequences of 75 broadly distributed woolly mammoths (Mammuthus primigenius) to test for signatures of positive selection. Results showed that a total of eleven amino acid sites in six genes, namely ND1, ND4, ND5, ND6, CYTB, and ATP6, displayed strong evidence of positive selection. Two sites were located in close proximity to proton-translocation channels in mitochondrial complex I. Biochemical and homology protein structure modeling analyses demonstrated that five amino acid substitutions in ND1, ND5, and ND6 might have influenced the performance of protein-protein interaction among subunits of complex I, and three substitutions in CYTB and ATP6 might have influenced the performance of metabolic regulatory chain. These findings suggest metabolic adaptations in the mitogenome of woolly mammoths in relation to extreme environments and provide a basis for further tests on the significance of the variations on other systems.

Rapid divergence, molecular evolution, and morphological diversification of coastal host-parasite systems from southern Brazil

This study assessed the role of historical processes on the geographic isolation, molecular evolution, and morphological diversification of host-parasite populations from the southern Brazilian coast. Adult specimens of Scleromystax barbatus and Scleromystax macropterus were collected from the sub-basin of the Nhundiaquara River and the sub-basin of the Paranaguá Bay, state of Paraná, Brazil. Four species of Gyrodactylus were recovered from the body surface of both host species. Morphometric analysis of Gyrodactylus spp. and Scleromystax spp. indicated that subpopulations of parasites and hosts could be distinguished from different sub-basins and locations, but the degree of morphological differentiation seems to be little related to geographic distance between subpopulations. Phylogenetic relationships based on DNA sequences of Gyrodactylus spp. and Scleromystax spp. allowed distinguishing lineages of parasites and hosts from different sub-basins. However, the level of genetic structuring of parasites was higher in comparison to host species. Evidence of positive selection in mtDNA sequences is likely associated with local adaptation of lineages of parasites and hosts. A historical demographic analysis revealed that populations of Gyrodactylus and Scleromystax have expanded in the last 250 000 years. The genetic variation of parasites and hosts is consistent with population-specific selection, population expansions, and recent evolutionary co-divergence.

An evaluation of alternative explanations for widespread cytonuclear discordance in annual sunflowers (Helianthus)

Cytonuclear discordance is commonly observed in phylogenetic studies, yet few studies have tested whether these patterns reflect incomplete lineage sorting or organellar introgression. Here, we used whole-chloroplast sequence data in combination with over 1000 nuclear single-nucleotide polymorphisms to clarify the extent of cytonuclear discordance in wild annual sunflowers (Helianthus), and to test alternative explanations for such discordance. Our phylogenetic analyses indicate that cytonuclear discordance is widespread within this group, both in terms of the relationships among species and among individuals within species. Simulations of chloroplast evolution show that incomplete lineage sorting cannot explain these patterns in most cases. Instead, most of the observed discordance is better explained by cytoplasmic introgression. Molecular tests of evolution further indicate that selection may have played a role in driving patterns of plastid variation - although additional experimental work is needed to fully evaluate the importance of selection on organellar variants in different parts of the geographic range. Overall, this study represents one of the most comprehensive tests of the drivers of cytonuclear discordance and highlights the potential for gene flow to lead to extensive organellar introgression in hybridizing taxa.

Whole-genome sequences of Malawi cichlids reveal multiple radiations interconnected by gene flow

The hundreds of cichlid fish species in Lake Malawi constitute the most extensive recent vertebrate adaptive radiation. Here we characterize its genomic diversity by sequencing 134 individuals covering 73 species across all major lineages. The average sequence divergence between species pairs is only 0.1-0.25%. These divergence values overlap diversity within species, with 82% of heterozygosity shared between species. Phylogenetic analyses suggest that diversification initially proceeded by serial branching from a generalist Astatotilapia-like ancestor. However, no single species tree adequately represents all species relationships, with evidence for substantial gene flow at multiple times. Common signatures of selection on visual and oxygen transport genes shared by distantly related deep-water species point to both adaptive introgression and independent selection. These findings enhance our understanding of genomic processes underlying rapid species diversification, and provide a platform for future genetic analysis of the Malawi radiation.

Spatio-temporal patterns of genetic variation in Arbacia lixula, a thermophilous sea urchin in expansion in the Mediterranean

The genetic structure of 13 populations of the amphiatlantic sea urchin Arbacia lixula, as well as temporal genetic changes in three of these localities, were assessed using ten hypervariable microsatellite loci. This thermophilous sea urchin is an important engineer species triggering the formation of barren grounds through its grazing activity. Its abundance seems to be increasing in most parts of the Mediterranean, probably favoured by warming conditions. Significant genetic differentiation was found both spatially and temporally. The main break corresponded to the separation of western Atlantic populations from those in eastern Atlantic and the Mediterranean Sea. A less marked, but significant differentiation was also found between Macaronesia (eastern Atlantic) and the Mediterranean. In the latter area, a signal of differentiation between the transitional area (Alboran Sea) and the rest of the Mediterranean was detected. However, no genetic structure is found within the Mediterranean (excluding Alboran) across the Siculo-Tunisian Strait, resulting from either enough gene flow to homogenize distance areas or/and a recent evolutionary history marked by demographic expansion in this basin. Genetic temporal variation at the Alboran Sea is as important as spatial variation, suggesting that temporal changes in hydrological features can affect the genetic composition of the populations. A picture of genetic homogeneity in the Mediterranean emerges, implying that the potential expansion of this keystone species will not be limited by intraspecific genetic features and/or potential impact of postulated barriers to gene flow in the region.

Role of oceanography in shaping the genetic structure in the North Pacific hake Merluccius productus

Determining the relative influence of biotic and abiotic factors on genetic connectivity among populations remains a major challenge in evolutionary biology and in the management and conservation of species. North Pacific hake (Merluccius productus) inhabits upwelling regions in the California Current ecosystem from the Gulf of California to the Gulf of Alaska. In this study, we examined mitochondrial DNA (mtDNA) and microsatellite variation to estimate levels of genetic differentiation of M. productus in relation to the role of oceanographic features as potential barriers to gene flow. Samples were obtained from nine sites spanning a large part of the geographic range of the species, from Puget Sound, Washington to Costa Rica. The microsatellite results revealed three genetically discrete populations: one spanning the eastern Pacific coast, and two apparently resident populations circumscribed to the Puget Sound and the northern Gulf of California (F_ST = 0.032, p = 0.036). Cytochrome b sequence data indicated that isolation between the Puget Sound and northern Gulf of California populations from the coastal Pacific were recent phenomena (18.5 kyr for Puget Sound and 40 kyr for the northern Gulf of California). Oceanographic data obtained from the Gulf of California support the hypothesis that permanent fronts within the region, and strong gradients at the entrance to the Gulf of California act as barriers to gene flow. A seascape genetics approach found significant genetic–environment associations, where the daytime sea surface temperature and chlorophyll concentrations were the best predictive variables for the observed genetic differentiation. Considering the potential causes of genetic isolation among the three populations, e.g. spawning areas in different latitudes associated with upwelling processes, oceanographic barriers, asymmetric migration and specialized diet, oceanographic barriers appear to be a likely mechanism restricting gene flow.

Genetic homogeneity of the invasive lionfish across the Northwestern Atlantic and the Gulf of Mexico based on Single Nucleotide Polymorphisms

Despite the devastating impact of the lionfish (Pterois volitans) invasion on NW Atlantic ecosystems, little genetic information about the invasion process is available. We applied Genotyping by Sequencing techniques to identify 1,220 single nucleotide polymorphic sites (SNPs) from 162 lionfish samples collected between 2013 and 2015 from two areas chronologically identified as the first and last invaded areas in US waters: the east coast of Florida and the Gulf of Mexico. We used population genomic analyses, including phylogenetic reconstruction, Bayesian clustering, genetic distances, Discriminant Analyses of Principal Components, and coalescence simulations for detection of outlier SNPs, to understand genetic trends relevant to the lionfish’s long-term persistence. We found no significant differences in genetic structure or diversity between the two areas (F_STp-values > 0.01, and t-test p-values > 0.05). In fact, our genomic analyses showed genetic homogeneity, with enough gene flow between the east coast of Florida and Gulf of Mexico to erase previous signals of genetic divergence detected between these areas, secondary spreading, and bottlenecks in the Gulf of Mexico. These findings suggest rapid genetic changes over space and time during the invasion, resulting in one panmictic population with no signs of divergence between areas due to local adaptation.

Mitochondrial Recombination Reveals Mito-Mito Epistasis in Yeast

Genetic variation in mitochondrial DNA (mtDNA) provides adaptive potential although the underlying genetic architecture of fitness components within mtDNAs is not known. To dissect functional variation within mtDNAs, we first identified naturally occurring mtDNAs that conferred high or low fitness in Saccharomyces cerevisiae by comparing growth in strains containing identical nuclear genotypes but different mtDNAs. During respiratory growth under temperature and oxidative stress conditions, mitotype effects were largely independent of nuclear genotypes even in the presence of mito-nuclear interactions. Recombinant mtDNAs were generated to determine fitness components within high- and low-fitness mtDNAs. Based on phenotypic distributions of isogenic strains containing recombinant mtDNAs, we found that multiple loci contributed to mitotype fitness differences. These mitochondrial loci interacted in epistatic, nonadditive ways in certain environmental conditions. Mito-mito epistasis (i.e., nonadditive interactions between mitochondrial loci) influenced fitness in progeny from four different crosses, suggesting that mito-mito epistasis is a widespread phenomenon in yeast and other systems with recombining mtDNAs. Furthermore, we found that interruption of coadapted mito-mito interactions produced recombinant mtDNAs with lower fitness. Our results demonstrate that mito-mito epistasis results in functional variation through mitochondrial recombination in fungi, providing modes for adaptive evolution and the generation of mito-mito incompatibilities.

Adaptive Patterns of Mitogenome Evolution Are Associated with the Loss of Shell Scutes in Turtles

The mitochondrial genome encodes several protein components of the oxidative phosphorylation (OXPHOS) pathway and is critical for aerobic respiration. These proteins have evolved adaptively in many taxa, but linking molecular-level patterns with higher-level attributes (e.g., morphology, physiology) remains a challenge. Turtles are a promising system for exploring mitochondrial genome evolution as different species face distinct respiratory challenges and employ multiple strategies for ensuring efficient respiration. One prominent adaptation to a highly aquatic lifestyle in turtles is the secondary loss of keratenized shell scutes (i.e., soft-shells), which is associated with enhanced swimming ability and, in some species, cutaneous respiration. We used codon models to examine patterns of selection on mitochondrial protein-coding genes along the three turtle lineages that independently evolved soft-shells. We found strong evidence for positive selection along the branches leading to the pig-nosed turtle (Carettochelys insculpta) and the softshells clade (Trionychidae), but only weak evidence for the leatherback (Dermochelys coriacea) branch. Positively selected sites were found to be particularly prevalent in OXPHOS Complex I proteins, especially subunit ND2, along both positively selected lineages, consistent with convergent adaptive evolution. Structural analysis showed that many of the identified sites are within key regions or near residues involved in proton transport, indicating that positive selection may have precipitated substantial changes in mitochondrial function. Overall, our study provides evidence that physiological challenges associated with adaptation to a highly aquatic lifestyle have shaped the evolution of the turtle mitochondrial genome in a lineage-specific manner.