Migration-selection balance and local adaptation of mitochondrial haplotypes in rufous-collared sparrows (Zonotrichia capensis) along an elevational gradient

Mitogenomics clarifies the position of the Nearctic magpies (Pica hudsonia and Pica nuttalli) within the Holarctic magpie radiation

Partial separation of a peripheral population may lead to its divergence and, potentially, speciation due to genetic drift followed by selection and geographic isolation. This process may cause taxonomic uncertainty because reproductive isolation in allopatry cannot be verified directly. The two Nearctic allopatric species of magpies (Aves, Corvidae: Pica) serve as a good example of these problems. The Black-billed magpie Pica hudsonia is widely distributed in North America, whereas the Yellow-billed Magpie Pica nuttalli is endemic to a restricted range in California. Their relationships with Palearctic species have been little studied. We obtained complete mitochondrial genomes of both Nearctic magpie species, along with the Eurasian Magpie (Pica pica) and the Oriental Magpie (Pica serica), 20 mitogenomes in total. Phylogenetic analysis reveals a basal position of P. serica, and P. pica as a sister clade to the two Nearctic species. P. hudsonia and P. nuttalli form reciprocal monophyletic subclades, showing recent divergence between and within them. Our data show that the Nearctic magpie lineage diverged from the common ancestor with P. pica, with a single migration wave via the Beringia. Within the Nearctic, we hypothesize a peripatric mode of speciation among Pica taxa due to the divergence and separation of the small marginal population in California below the Sierra-Nevada mountains. Diversifying amino acid substitutions in ND4-ND5-ND6 genes along the branch leading to the New World clade may indicate selection for heat-tolerance. Considering the clear phenotypic differences between P. hudsonia and P. nuttalli, our data, showing their reciprocal monophylies and genetic distinctness, is consistent with the two-species taxonomy.

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.

Population genetic structure and hybrid zone analyses for species delimitation in the Japanese toad (Bufo japonicus)

Hybridization following secondary contact may produce different outcomes depending on the extent to which genetic diversity and reproductive barriers have accumulated during isolation. The Japanese toad, Bufo japonicus, is distributed on the main islands of Japan. In the present study, we applied multiplexed inter-simple sequence repeat genotyping by sequencing to achieve the fine-scale resolution of the genetic cluster in B. j. japonicus and B. j. formosus. We also elucidated hybridization patterns and gene flow degrees across contact zones between the clusters identified. Using SNP data, we found four genetic clusters in B. j. japonicus and B. j. formosus and three contact zones of the cluster pairs among these four clusters. The two oldest diverged lineages, B. j. japonicus and B. j. formosus, formed a narrow contact zone consistent with species distinctiveness. Therefore, we recommend that these two subspecies be elevated to the species level. In contrast, the less diverged pairs of two clusters in B. j. japonicus and B. j. formosus, respectively, admixed over a hundred kilometers, suggesting that they have not yet developed strong reproductive isolation and need to be treated as conspecifics. These results will contribute to resolving taxonomic confusion in Japanese toads.

Disentangling Positive Selection from Relaxed Selection in Animal Mitochondrial Genomes

AbstractDisentangling different types of selection is a common goal in molecular evolution. Elevated dN/dS ratios (the ratio of nonsynonymous to synonymous substitution rates) in focal lineages are often interpreted as signs of positive selection. Paradoxically, relaxed purifying selection can also result in elevated dN/dS ratios, but tests to distinguish these two causes are seldomly implemented. Here, we reevaluated seven case studies describing elevated dN/dS ratios in animal mitochondrial DNA (mtDNA) and their accompanying hypotheses regarding selection. They included flightless lineages versus flighted lineages in birds, bats, and insects and physiological adaptations in snakes, two groups of electric fishes, and primates. We found that elevated dN/dS ratios were often not caused by the predicted mechanism, and we sometimes found strong support for the opposite mechanism. We discuss reasons why energetic hypotheses may be confounded by other selective forces acting on mtDNA and caution against overinterpreting singular molecular signals, including elevated dN/dS ratios.

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.

Widespread genomic signatures of reproductive isolation and sex-specific selection in the Eastern Yellow Robin, Eopsaltria australis

How new species evolve is one of the most fundamental questions in biology. Population divergence, which may lead to speciation, may be occurring in the Eastern Yellow Robin, a common passerine that lives along the eastern coast of Australia. This species is composed of 2 parapatric lineages that have highly divergent mitochondrial DNA; however, similar levels of divergence have not been observed in the nuclear genome. Here we re-examine the nuclear genomes of these mitolineages to test potential mechanisms underlying the discordance between nuclear and mitochondrial divergence. We find that nuclear admixture occurs in a narrow hybrid zone, although the majority of markers across the genome show evidence of reproductive isolation between populations of opposing mitolineages. There is an 8 MB section of a previously identified putative neo-sex chromosome that is highly diverged between allopatric but not parapatric populations, which may be the result of a chromosomal inversion. The neo-sex chromosomal nature of this region, as well as the geographic patterns in which it exhibits divergence, suggest it is unlikely to be contributing to reproductive isolation through mitonuclear incompatibilities as reported in earlier studies. In addition, there are sex differences in the number of markers that are differentiated between populations of opposite mitolineages, with greater differentiation occurring in females, which are heterozygous, than males. These results suggest that, despite the absence of previously observed assortative mating, mitolineages of Eastern Yellow Robin experience at least some postzygotic isolation from each other, in a pattern consistent with Haldane’s Rule.

Avian introgression patterns are consistent with Haldane's Rule

According to Haldane’s Rule, the heterogametic sex will show the greatest fitness reduction in a hybrid cross. In birds, where sex is determined by a ZW system, female hybrids are expected to experience lower fitness compared to male hybrids. This pattern has indeed been observed in several bird groups, but it is unknown whether the generality of Haldane’s Rule also extends to the molecular level. First, given the lower fitness of female hybrids, we can expect maternally inherited loci (i.e., mitochondrial and W-linked loci) to show lower introgression rates than biparentally inherited loci (i.e., autosomal loci) in females. Second, the faster evolution of Z-linked loci compared to autosomal loci and the hemizygosity of the Z-chromosome in females might speed up the accumulation of incompatible alleles on this sex chromosome, resulting in lower introgression rates for Z-linked loci than for autosomal loci. I tested these expectations by conducting a literature review which focused on studies that directly quantified introgression rates for autosomal, sex-linked, and mitochondrial loci. Although most studies reported introgression rates in line with Haldane’s Rule, it remains important to validate these genetic patterns with estimates of hybrid fitness and supporting field observations to rule out alternative explanations. Genomic data provide exciting opportunities to obtain a more fine-grained picture of introgression rates across the genome, which can consequently be linked to ecological and behavioral observations, potentially leading to novel insights into the genetic mechanisms underpinning Haldane’s Rule.

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.

Selection for biparental inheritance of mitochondria under hybridization and mitonuclear fitness interactions

Uniparental inheritance (UPI) of mitochondria predominates over biparental inheritance (BPI) in most eukaryotes. However, examples of BPI of mitochondria, or paternal leakage, are becoming increasingly prevalent. Most reported cases of BPI occur in hybrids of distantly related sub-populations. It is thought that BPI in these cases is maladaptive; caused by a failure of female or zygotic autophagy machinery to recognize divergent male-mitochondrial DNA ‘tags’. Yet recent theory has put forward examples in which BPI can evolve under adaptive selection, and empirical studies across numerous metazoan taxa have demonstrated outbreeding depression in hybrids attributable to disruption of population-specific mitochondrial and nuclear genotypes (mitonuclear mismatch). Based on these developments, we hypothesize that BPI may be favoured by selection in hybridizing populations when fitness is shaped by mitonuclear interactions. We test this idea using a deterministic, simulation-based population genetic model and demonstrate that BPI is favoured over strict UPI under moderate levels of gene flow typical of hybridizing populations. Our model suggests that BPI may be stable, rather than a transient phenomenon, in hybridizing populations.

Patterns of morphological differentiation within Manacus manacus (Aves: Pipridae) in Colombia: revisiting hypotheses of isolation and secondary contact

Determining the factors responsible for phenotypic variation within species is a fundamental task in evolutionary ecology. Isolation by distance, isolation by environment and clines along secondary contact zones between formerly isolated populations are common patterns of morphological variation. In this study, we evaluated whether patterns of morphological variation exhibit association with isolation by distance, isolation by environment or secondary contact in populations of Manacus manacus with contrasting plumage colour. We used field and museum measurements of five morphological characters from 311 individuals from 88 localities distributed across a highly heterogeneous environmental range within Colombia. Climatic variables explained a higher proportion of the variation than geographic distance among localities. We found differences in wing and culmen length between white and yellow phenotypes of M. manacus in Colombia. Overall morphology did not change in a clinal fashion in relation to the contact zone and the few traits that exhibit clinal variation suggest a distant locality for the contact zone, not in line with our expectations by secondary contact between differentiated populations, suggesting alternative mechanisms for differentiation between subspecies. Our results suggest that this species exhibits marked morphological variation associated with changes in temperature and precipitation. These associations are consistent throughout the geographic range analysed, leading us to postulate that dispersal and local adaptation jointly shape the distribution of phenotypic variation.

Habitat-linked genetic structure for white-crowned sparrow (Zonotrichia leucophrys): Local factors shape population genetic structure

Ecological, environmental, and geographic factors all influence genetic structure. Species with broad distributions are ideal systems because they cover a range of ecological and environmental conditions allowing us to test which components predict genetic structure. This study presents a novel, broad geographic approach using molecular markers, morphology, and habitat modeling to investigate rangewide and local barriers causing contemporary genetic differentiation within the geographical range of three white-crowned sparrow (Zonotrichia leucophrys) subspecies: Z. l. gambelii, Z. l. oriantha, and Z. l. pugetensis. Three types of genetic markers showed geographic distance between sampling sites, elevation, and ecosystem type are key factors contributing to population genetic structure. Microsatellite markers revealed white-crowned sparrows do not group by subspecies, but instead indicated four groupings at a rangewide scale and two groupings based on coniferous and deciduous ecosystems at a local scale. Our analyses of morphological variation also revealed habitat differences; sparrows from deciduous ecosystems are larger than individuals from coniferous ecosystems based on principal component analyses. Habitat modeling showed isolation by distance was prevalent in describing genetic structure, but isolation by resistance also had a small but significant influence. Not only do these findings have implications concerning the accuracy of subspecies delineations, they also highlight the critical role of local factors such as habitat in shaping contemporary population genetic structure of species with high dispersal ability.

Mitogenomic analysis of extant condor species provides insight into the molecular evolution of vultures

The evolution of large vultures linked to mountainous habitats was accompanied by extreme physiological and behavioral specializations for energetically efficient flights. However, little is known on the genetic traits associated with the evolution of these obligate soaring scavengers. Mitochondrial DNA plays a vital role in regulating oxidative stress and energy production, and hence may be an important target of selection for flight performance. Herein, we characterized the first mitogenomes of the Andean and California condors, the world's heaviest flying birds and the only living representatives of the Vultur and Gymnogyps genus. We reconstructed the phylogenetic relationships and evaluated possible footprints of convergent evolution associated to the life-history traits and distributional range of vultures. Our phylogenomic analyses supported the independent evolution of vultures, with the origin of Cathartidae in the early Paleogene (~ 61 Mya), and estimated the radiation of extant condors during the late Miocene (~ 11 Mya). Selection analyses indicated that vultures exhibit signals of relaxation of purifying selection relative to other accipitrimorph raptors, possibly indicating the degeneration of flapping flight ability. Overall, our results suggest that the extreme specialization of vultures for efficient soaring flight has compensated the evolution of large body sizes mitigating the selection pressure on mtDNA.

Phenotypic Variation in Mitochondria-Related Performance Traits Across New Zealand Snail Populations

Mitochondrial function is critical for energy homeostasis and should shape how genetic variation in metabolism is transmitted through levels of biological organization to generate stability in organismal performance. Mitochondrial function is encoded by genes in two distinct and separately inherited genomes-the mitochondrial genome and the nuclear genome-and selection is expected to maintain functional mito-nuclear interactions. The documented high levels of polymorphism in genes involved in these mito-nuclear interactions and wide variation for mitochondrial function demands an explanation for how and why variability in such a fundamental trait is maintained. Potamopyrgus antipodarum is a New Zealand freshwater snail with coexisting sexual and asexual individuals and, accordingly, contrasting systems of separate vs. co-inheritance of nuclear and mitochondrial genomes. As such, this snail provides a powerful means to dissect the evolutionary and functional consequences of mito-nuclear variation. The lakes inhabited by P. antipodarum span wide environmental gradients, with substantial across-lake genetic structure and mito-nuclear discordance. This situation allows us to use comparisons across reproductive modes and lakes to partition variation in cellular respiration across genetic and environmental axes. Here, we integrated cellular, physiological, and behavioral approaches to quantify variation in mitochondrial function across a diverse set of wild P. antipodarum lineages. We found extensive across-lake variation in organismal oxygen consumption and behavioral response to heat stress and differences across sexes in mitochondrial membrane potential but few global effects of reproductive mode. Taken together, our data set the stage for applying this important model system for sexual reproduction and polyploidy to dissecting the complex relationships between mito-nuclear variation, performance, plasticity, and fitness in natural populations.

Pervasive Genomic Signatures of Local Adaptation to Altitude Across Highland Specialist Andean Hummingbird Populations

Populations along steep environmental gradients are subject to differentiating selection that can result in local adaptation, despite countervailing gene flow, and genetic drift. In montane systems, where species are often restricted to narrow ranges of elevation, it is unclear whether the selection is strong enough to influence functional differentiation of subpopulations differing by a few hundred meters in elevation. We used targeted capture of 12 501 exons from across the genome, including 271 genes previously implicated in altitude adaptation, to test for adaptation to local elevations for 2 highland hummingbird species, Coeligena violifer (n = 62) and Colibri coruscans (n = 101). For each species, we described population genetic structure across the complex geography of the Peruvian Andes and, while accounting for this structure, we tested whether elevational allele frequency clines in single nucleotide polymorphisms (SNPs) showed evidence for local adaptation to elevation. Although the 2 species exhibited contrasting population genetic structures, we found signatures of clinal genetic variation with shifts in elevation in both. The genes with SNP-elevation associations included candidate genes previously discovered for high-elevation adaptation as well as others not previously identified, with cellular functions related to hypoxia response, energy metabolism, and immune function, among others. Despite the homogenizing effects of gene flow and genetic drift, natural selection on parts of the genome evidently optimizes elevation-specific cellular function even within elevation range-restricted montane populations. Consequently, our results suggest local adaptation occurring in narrow elevation bands in tropical mountains, such as the Andes, may effectively make them "taller" biogeographic barriers.

The role of selection in the evolution of marine turtles mitogenomes

Sea turtles are the only extant chelonian representatives that inhabit the marine environment. One key to successful colonization of this habitat is the adaptation to different energetic demands. Such energetic requirement is intrinsically related to the mitochondrial ability to generate energy through oxidative phosphorylation (OXPHOS) process. Here, we estimated Testudines phylogenetic relationships from 90 complete chelonian mitochondrial genomes and tested the adaptive evolution of 13 mitochondrial protein-coding genes of sea turtles to determine how natural selection shaped mitochondrial genes of the Chelonioidea clade. Complete mitogenomes showed strong support and resolution, differing at the position of the Chelonioidea clade in comparison to the turtle phylogeny based on nuclear genomic data. Codon models retrieved a relatively increased dN/dS (ω) on three OXPHOS genes for sea turtle lineages. Also, we found evidence of positive selection on at least three codon positions, encoded by NADH dehydrogenase genes (ND4 and ND5). The accelerated evolutionary rates found for sea turtles on COX2, ND1 and CYTB and the molecular footprints of positive selection found on ND4 and ND5 genes may be related to mitochondrial molecular adaptation to stress likely resulted from a more active lifestyle in sea turtles. Our study provides insight into the adaptive evolution of the mtDNA genome in sea turtles and its implications for the molecular mechanism of oxidative phosphorylation.

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.

Little genetic structure in a Bornean endemic small mammal across a steep ecological gradient

Janzen's influential "mountain passes are higher in the tropics" hypothesis predicts restricted gene flow and genetic isolation among populations spanning elevational gradients in the tropics. Few studies have tested this prediction, and studies that focus on population genetic structure in Southeast Asia are particularly underrepresented in the literature. Here, we test the hypothesis that mountain treeshrews (Tupaia montana) exhibit limited dispersal across their broad elevational range which spans ~2,300 m on two peaks in Kinabalu National Park (KNP) in Borneo: Mt Tambuyukon (MT) and Mt Kinabalu (MK). We sampled 83 individuals across elevations on both peaks and performed population genomics analyses on mitogenomes and single nucleotide polymorphisms from 4,106 ultraconserved element loci. We detected weak genetic structure and infer gene flow both across elevations and between peaks. We found higher genetic differentiation on MT than MK despite its lower elevation and associated environmental variation. This implies that, contrary to our hypothesis, genetic structure in this system is not primarily shaped by elevation. We propose that this pattern may instead be the result of historical processes and limited upslope gene flow on MT. Importantly, our results serve as a foundational estimate of genetic diversity and population structure from which to track potential future effects of climate change on mountain treeshrews in KNP, an important conservation stronghold for the mountain treeshrew and other montane species.

Dispersal Predicts Hybrid Zone Widths across Animal Diversity: Implications for Species Borders under Incomplete Reproductive Isolation

Hybrid zones occur as range boundaries for many animal taxa. One model for how hybrid zones form and stabilize is the tension zone model, a version of which predicts that hybrid zone widths are determined by a balance between random dispersal into hybrid zones and selection against hybrids. Here, we examine whether random dispersal and proxies for selection against hybrids (genetic distances between hybridizing pairs) can explain variation in hybrid zone widths across 131 hybridizing pairs of animals. We show that these factors alone can explain ∼40% of the variation in zone width among animal hybrid zones, with dispersal explaining far more of the variation than genetic distances. Patterns within clades were idiosyncratic. Genetic distances predicted hybrid zone widths particularly well for reptiles, while this relationship was opposite tension zone predictions in birds. Last, the data suggest that dispersal and molecular divergence set lower bounds on hybrid zone widths in animals, indicating that there are geographic restrictions on hybrid zone formation. Overall, our analyses reinforce the fundamental importance of dispersal in hybrid zone formation and more generally in the ecology of range boundaries.

Panmixia across elevation in thermally sensitive Andean dung beetles

Janzen's seasonality hypothesis predicts that organisms inhabiting environments with limited climatic variability will evolve a reduced thermal tolerance breadth compared with organisms experiencing greater climatic variability. In turn, narrow tolerance breadth may select against dispersal across strong temperature gradients, such as those found across elevation. This can result in narrow elevational ranges and generate a pattern of isolation by environment or neutral genetic differentiation correlated with environmental variables that are independent of geographic distance. We tested for signatures of isolation by environment across elevation using genome-wide SNP data from five species of Andean dung beetles (subfamily Scarabaeinae) with well-characterized, narrow thermal physiologies, and narrow elevational distributions. Contrary to our expectations, we found no evidence of population genetic structure associated with elevation and little signal of isolation by environment. Further, elevational ranges for four of five species appear to be at equilibrium and show no decay of genetic diversity at range limits. Taken together, these results suggest physiological constraints on dispersal may primarily operate outside of a stable realized niche and point to a lower bound on the spatial scale of local adaptation.

Mitonuclear mismatch alters performance and reproductive success in naturally introgressed populations of a montane leaf beetle

Coordination between nuclear and mitochondrial genomes is critical to metabolic processes underlying animals' ability to adapt to local environments, yet consequences of mitonuclear interactions have rarely been investigated in populations where individuals with divergent mitochondrial and nuclear genomes naturally interbreed. Genetic variation in the leaf beetle Chrysomela aeneicollis was assessed along a latitudinal thermal gradient in California's Sierra Nevada. Variation at mitochondrial cytochrome oxidase II (COII) and the nuclear gene phosphoglucose isomerase (PGI) shows concordance and was significantly greater along a 65 km transect than 10 other loci. STRUCTURE analyses using neutral loci identified a southern and northern subpopulation, which interbreed in the central drainage Bishop Creek. COII and PGI were used as indicators of mitochondrial and nuclear genetic variation in field and laboratory experiments conducted on beetles from this admixed population. Fecundity, larval development rate, running speed and male mating frequency were higher for beetles with geographically "matched" than "mismatched" mitonuclear genotypes. Effects of mitonuclear mismatch were largest for individuals with northern nuclear genotypes possessing southern mitochondria and were most pronounced after heat treatment or at high elevation. These findings suggest that mitonuclear incompatibility diminishes performance and reproductive success in nature, effects that could intensify at environmental extremes.

Cryptic and extensive hybridization between ancient lineages of American crows

Most species and therefore most hybrid zones have historically been defined using phenotypic characters. However, both speciation and hybridization can occur with negligible morphological differentiation. Recently developed genomic tools provide the means to better understand cryptic speciation and hybridization. The Northwestern Crow (Corvus caurinus) and American Crow (Corvus brachyrhynchos) are continuously distributed sister taxa that lack reliable traditional characters for identification. In this first population genomic study of Northwestern and American crows, we use genomic SNPs (nuDNA) and mtDNA to investigate the degree of genetic differentiation between these crows and the extent to which they may hybridize. Our results indicate that American and Northwestern crows have distinct evolutionary histories, supported by two nuDNA ancestry clusters and two 1.1%-divergent mtDNA clades dating to the late Pleistocene, when glacial advances may have isolated crow populations in separate refugia. We document extensive hybridization, with geographic overlap of mtDNA clades and admixture of nuDNA across >900 km of western Washington and western British Columbia. This broad hybrid zone consists of late-generation hybrids and backcrosses, but not recent (e.g., F1) hybrids. Nuclear DNA and mtDNA clines had concordant widths and were both centred in southwestern British Columbia, farther north than previously postulated. Overall, our results suggest a history of reticulate evolution in American and Northwestern crows, perhaps due to recurring neutral expansion(s) from Pleistocene glacial refugia followed by lineage fusion(s). However, we do not rule out a contributing role for more recent potential drivers of hybridization, such as expansion into human-modified habitats.

Strong phenotypic divergence in spite of low genetic structure in the endemic Mangrove Warbler subspecies (Setophaga petechia xanthotera) of Costa Rica

Despite the enormous advances in genetics, links between phenotypes and genotypes have been made for only a few nonmodel organisms. However, such links can be essential to understand mechanisms of ecological speciation. The Costa Rican endemic Mangrove Warbler subspecies provides an excellent subject to study differentiation with gene flow, as it is distributed along a strong precipitation gradient on the Pacific coast with no strong geographic barriers to isolate populations. Mangrove Warbler populations could be subject to divergent selection driven by precipitation, which influences soil salinity levels, which in turn influences forest structure and food resources. We used single nucleotide polymorphisms (SNPs) and morphological traits to examine the balance between neutral genetic and phenotypic divergence to determine whether selection has acted on traits and genes with functions related to specific environmental variables. We present evidence showing: (a) associations between environmental variables and SNPs, identifying candidate genes related to bill morphology (BMP) and osmoregulation, (b) absence of population genetic structure in neutrally evolving markers, (c) divergence in bill size across the precipitation gradient, and (d) strong phenotypic differentiation (P ST) which largely exceeds neutral genetic differentiation (F ST) in bill size. Our results indicate an important role for salinity, forest structure, and resource availability in maintaining phenotypic divergence of Mangrove Warblers through natural selection. Our findings add to the growing body of literature identifying the processes involved in phenotypic differentiation along environmental gradients in the face of gene flow.

Molecular genetic studies of Chilean avifauna: an overview about current progress

The biodiversity of Chilean avifauna has been shaped by geological events, geographical factors and extreme climatic regimens. Molecular genetic studies developed to date have begun to provide valuable information to elucidate the effect of these factors on the evolutionary history of Chilean birds. We conducted a search in two bibliographic databases (PubMed and Scopus) over the period between 1998 and 2018 to identify the number and class of articles published in this research field. The quantification of specific research categories was also performed to assess progress in different molecular genetic issues of Chilean birds. Representative studies of the more active research areas: population structure, phylogeography, systematics and taxonomy and genetic diversity, were briefly commented on. Studies on genetic divergence analysis of cryptic species and the underlying genetics at molecular level of adaptive phenotypes in populations, were very scarce. Several types of mitochondrial and nuclear markers have been used, with COI plus D-loop and microsatellite loci, respectively, being the most widely-used markers. However, high-throughput analyses based on complete mitochondrial genomes have also been published. More than one hundred species have been analyzed; the most studied groups being penguins, followed by New World swallows. Although, to date, molecular genetic studies of Chilean birds still require further analysis, the studies available on some species or groups have provided important data to improve our understanding of the origin and evolution of this avifauna.

Parallel Molecular Evolution in Pathways, Genes, and Sites in High-Elevation Hummingbirds Revealed by Comparative Transcriptomics

High-elevation organisms experience shared environmental challenges that include low oxygen availability, cold temperatures, and intense ultraviolet radiation. Consequently, repeated evolution of the same genetic mechanisms may occur across high-elevation taxa. To test this prediction, we investigated the extent to which the same biochemical pathways, genes, or sites were subject to parallel molecular evolution for 12 Andean hummingbird species (family: Trochilidae) representing several independent transitions to high elevation across the phylogeny. Across high-elevation species, we discovered parallel evolution for several pathways and genes with evidence of positive selection. In particular, positively selected genes were frequently part of cellular respiration, metabolism, or cell death pathways. To further examine the role of elevation in our analyses, we compared results for low- and high-elevation species and tested different thresholds for defining elevation categories. In analyses with different elevation thresholds, positively selected genes reflected similar functions and pathways, even though there were almost no specific genes in common. For example, EPAS1 (HIF2α), which has been implicated in high-elevation adaptation in other vertebrates, shows a signature of positive selection when high-elevation is defined broadly (>1,500 m), but not when defined narrowly (>2,500 m). Although a few biochemical pathways and genes change predictably as part of hummingbird adaptation to high-elevation conditions, independent lineages have rarely adapted via the same substitutions.

A continent-scale study of the social structure and phylogeography of the bent-wing bat, Miniopterus schreibersii (Mammalia: Chiroptera), using new microsatellite data

Miniopterus schreibersii is a cave-dwelling bat species with a wide distribution in the western Palearctic spanning southern and central Europe, North Africa, Anatolia, the Middle East, and the Caucasus. We investigated the social structure and its effects on the genetic makeup of this species, using 10 nuclear microsatellite markers and a partial fragment of the mitochondrial cytochrome b gene. Samples were examined from the species' entire circum-Mediterranean range. Local structuring that was previously detected among populations of M. schreibersii using mitochondrial markers was not observed for microsatellite markers, indicating male-biased dispersal for the species. Some support was found for postglacial expansions in Europe, with Anatolia potentially acting as the primary refugium during the Last Glacial Maximum (LGM). However, support for this hypothesis is not as strong as that previously detected using mitochondrial DNA markers. This is likely due to the diminishing effect of male-mediated dispersal, replenishing the nuclear diversity faster than the mitochondrial diversity in regions that are relatively far from the glacial refugia.

Genetic Diversity and Wolbachia Infection Patterns in a Globally Distributed Invasive Ant

Understanding the phylogeographic history of an invasive species may facilitate reconstructing the history and routes of its invasion. The longhorn crazy ant, Paratrechina longicornis, is a ubiquitous agricultural and household pest throughout much of the tropics and subtropics, but little is known about the history of its spread. Here, we examine worldwide genetic variation in P. longicornis and its associated Wolbachia bacterial symbionts. Analyses of mtDNA sequences of 248 P. longicornis workers (one per colony) from 13 geographic regions reveal two highly diverged mtDNA clades that co-occur in most of the geographic regions. These two mtDNA clades are associated with different Wolbachia infection patterns, but are not congruent with patterns of nDNA (microsatellite) variation. Multilocus sequence typing reveals two distinct Wolbachia strains in P. longicornis, namely, wLonA and wLonF. The evolutionary histories of these two strains differ; wLonA appears to be primarily transmitted maternally, and patterns of mtDNA and nDNA variation and wLonA infection status are consistent with a relatively recent Wolbachia-induced selective sweep. In contrast, the observed patterns of mtDNA variation and wLonF infections suggest frequent horizontal transfer and losses of wLonF infections. The lack of nDNA structure among sampled geographic regions coupled with the finding that numerous mtDNA haplotypes are shared among regions implies that inadvertent long-distance movement through human commerce is common in P. longicornis and has shaped the genetic structure of this invasive ant worldwide.

A mirage of cryptic species: Genomics uncover striking mitonuclear discordance in the butterfly Thymelicus sylvestris

Mitochondrial DNA (mtDNA) sequencing has led to an unprecedented rise in the identification of cryptic species. However, it is widely acknowledged that nuclear DNA (nuDNA) sequence data are also necessary to properly define species boundaries. Next generation sequencing techniques provide a wealth of nuclear genomic data, which can be used to ascertain both the evolutionary history and taxonomic status of putative cryptic species. Here, we focus on the intriguing case of the butterfly Thymelicus sylvestris (Lepidoptera: Hesperiidae). We identified six deeply diverged mitochondrial lineages; three distributed all across Europe and found in sympatry, suggesting a potential case of cryptic species. We then sequenced these six lineages using double-digest restriction-site associated DNA sequencing (ddRADseq). Nuclear genomic loci contradicted mtDNA patterns and genotypes generally clustered according to geography, i.e., a pattern expected under the assumption of postglacial recolonization from different refugia. Further analyses indicated that this strong mtDNA/nuDNA discrepancy cannot be explained by incomplete lineage sorting, sex-biased asymmetries, NUMTs, natural selection, introgression or Wolbachia-mediated genetic sweeps. We suggest that this mitonuclear discordance was caused by long periods of geographic isolation followed by range expansions, homogenizing the nuclear but not the mitochondrial genome. These results highlight T. sylvestris as a potential case of multiple despeciation and/or lineage fusion events. We finally argue, since mtDNA and nuDNA do not necessarily follow the same mechanisms of evolution, their respective evolutionary history reflects complementary aspects of past demographic and biogeographic events.

Divergent Fine-Scale Recombination Landscapes between a Freshwater and Marine Population of Threespine Stickleback Fish

Meiotic recombination is a highly conserved process that has profound effects on genome evolution. At a fine-scale, recombination rates can vary drastically across genomes, often localized into small recombination "hotspots" with highly elevated rates, surrounded by regions with little recombination. In most species studied, the location of hotspots within genomes is highly conserved across broad evolutionary timescales. The main exception to this pattern is in mammals, where hotspot location can evolve rapidly among closely related species and even among populations within a species. Hotspot position in mammals is controlled by the gene, Prdm9, whereas in species with conserved hotspots, a functional Prdm9 is typically absent. Due to a limited number of species where recombination rates have been estimated at a fine-scale, it remains unclear whether hotspot conservation is always associated with the absence of a functional Prdm9. Threespine stickleback fish (Gasterosteus aculeatus) are an excellent model to examine the evolution of recombination over short evolutionary timescales. Using a linkage disequilibrium-based approach, we found recombination rates indeed varied at a fine-scale across the genome, with many regions organized into narrow hotspots. Hotspots had highly divergent landscapes between stickleback populations, where only ∼15% of these hotspots were shared. Our results indicate that fine-scale recombination rates may be diverging between closely related populations of threespine stickleback fish. Interestingly, we found only a weak association of a PRDM9 binding motif within hotspots, which suggests that threespine stickleback fish may possess a novel mechanism for targeting recombination hotspots at a fine-scale.

Assessing the fitness consequences of mitonuclear interactions in natural populations

Metazoans exist only with a continuous and rich supply of chemical energy from oxidative phosphorylation in mitochondria. The oxidative phosphorylation machinery that mediates energy conservation is encoded by both mitochondrial and nuclear genes, and hence the products of these two genomes must interact closely to achieve coordinated function of core respiratory processes. It follows that selection for efficient respiration will lead to selection for compatible combinations of mitochondrial and nuclear genotypes, and this should facilitate coadaptation between mitochondrial and nuclear genomes (mitonuclear coadaptation). Herein, we outline the modes by which mitochondrial and nuclear genomes may coevolve within natural populations, and we discuss the implications of mitonuclear coadaptation for diverse fields of study in the biological sciences. We identify five themes in the study of mitonuclear interactions that provide a roadmap for both ecological and biomedical studies seeking to measure the contribution of intergenomic coadaptation to the evolution of natural populations. We also explore the wider implications of the fitness consequences of mitonuclear interactions, focusing on central debates within the fields of ecology and biomedicine.

Landscape effects on the contemporary genetic structure of Ruffed Grouse (Bonasa umbellus) populations

The amount of dispersal that occurs among populations can be limited by landscape heterogeneity, which is often due to both natural processes and anthropogenic activity leading to habitat loss or fragmentation. Understanding how populations are structured and mapping existing dispersal corridors among populations is imperative to both determining contemporary forces mediating population connectivity, and informing proper management of species with fragmented populations. Furthermore, the contemporary processes mediating gene flow across heterogeneous landscapes on a large scale are understudied, particularly with respect to widespread species. This study focuses on a widespread game bird, the Ruffed Grouse (Bonasa umbellus), for which we analyzed samples from the western extent of the range. Using three types of genetic markers, we uncovered multiple factors acting in concert that are responsible for mediating contemporary population connectivity in this species. Multiple genetically distinct groups were detected; microsatellite markers revealed six groups, and a mitochondrial marker revealed four. Many populations of Ruffed Grouse are genetically isolated, likely by macrogeographic barriers. Furthermore, the addition of landscape genetic methods not only corroborated genetic structure results, but also uncovered compelling evidence that dispersal resistance created by areas of unsuitable habitat is the most important factor mediating population connectivity among the sampled populations. This research has important implications for both our study species and other inhabitants of the early successional forest habitat preferred by Ruffed Grouse. Moreover, it adds to a growing body of evidence that isolation by resistance is more prevalent in shaping population structure of widespread species than previously thought.

Landscape genetic structure of Scirpus mariqueter reveals a putatively adaptive differentiation under strong gene flow in estuaries

Estuarine organisms grow in highly heterogeneous habitats, and their genetic differentiation is driven by selective and neutral processes as well as population colonization history. However, the relative importance of the processes that underlie genetic structure is still puzzling. Scirpus mariqueter is a perennial grass almost limited in the Changjiang River estuary and its adjacent Qiantang River estuary. Here, using amplified fragment length polymorphism (AFLP), a moderate-high level of genetic differentiation among populations (range F ST: 0.0310-0.3325) was showed despite large ongoing dispersal. FLOCK assigned all individuals to 13 clusters and revealed a complex genetic structure. Some genetic clusters were limited in peripheries compared with very mixing constitution in center populations, suggesting local adaptation was more likely to occur in peripheral populations. 21 candidate outliers under positive selection were detected, and further, the differentiation patterns correlated with geographic distance, salinity difference, and colonization history were analyzed with or without the outliers. Combined results of AMOVA and IBD based on different dataset, it was found that the effects of geographic distance and population colonization history on isolation seemed to be promoted by divergent selection. However, none-liner IBE pattern indicates the effects of salinity were overwhelmed by spatial distance or other ecological processes in certain areas and also suggests that salinity was not the only selective factor driving population differentiation. These results together indicate that geographic distance, salinity difference, and colonization history co-contributed in shaping the genetic structure of S. mariqueter and that their relative importance was correlated with spatial scale and environment gradient.

New lineages and old species: Lineage diversity and regional distribution of Moina (Crustacea: Cladocera) in China

The distribution and genetic diversity of freshwater zooplankton is understudied in the Eastern Palearctic. Here, we explored the lineage diversity and regional distribution of the genus Moina in China. Members of this genus are often keystone components of freshwater ecosystems and have been frequently subjected to toxicological and physiological studies. Four species of Moina were identified, based on morphology, in 50 of 113 Chinese water bodies examined, and their phylogenetic position was analyzed using both a mitochondrial (mitochondrial cytochrome c oxidase subunit I; COI) and a nuclear marker (the nuclear internal transcribed spacer; ITS-1). Both molecular markers identified four clades corresponding broadly to the morphological species. Mitochondrial DNA analysis showed the presence of four species complexes with eleven lineages across China, five of which were new. However, some lineages (and even individual haplotypes) were widespread in Eurasia, suggesting an ability to disperse over long distances. In contrast, a few lineages exhibited restricted distributions. The nuclear phylogeny also recognized four species of Moina within China and seven very distinct clades. Interestingly, one specimen possessing Moina cf. micrura mtDNA had ITS-1 alleles of the M. cf. brachiata clade. This discordance between mtDNA and nuclear ITS-1 phylogenies is indicative of interspecific introgression and hybridization. Additionally, our COI phylogeny showed apparent paraphyly in two Moina species groups, suggesting introgression of their mitochondrial genomes. Our data shows the regional distribution/diversity of the Moina species complex in a Eurasian context.

The Phytogeographic History of Common Walnut in China

Common walnut (Juglans regia L.) is an economically important hardwood tree species cultivated worldwide for its high quality wood and edible nuts. It is generally accepted that after the last glaciation J. regia survived and grew in almost completely isolated stands in Asia, and that ancient humans dispersed walnuts across Asia and into new habitats via trade and cultural expansion. The history of common walnut in China is a matter of debate, however. We estimated the genetic diversity and spatial genetic structure of 31 walnut populations sampled across its Chinese range using 22 microsatellite markers (13 neutral and 9 non-neutral). Using historical data and population genetic analysis, including approximate Bayesian analysis (ABC), we reconstructed the demographic history of J. regia in China. The genetic data indicated the likely presence of J. regia in glacial refugia in the Xinjiang province (Northwest China), Northeastern China (Beijing, Shandong, and Changbai Mountains), Central China (Qinling and Baishan Mountains and Xi'an), and Southwestern China (Tibet, Yunnan, Guizhou, and Sichuan provinces). Based on DIY-ABC analysis, we identified three ancient lineages of J. regia in China. Two lineages (subpopulation A and subpopulation B+C) diverged about 2.79 Mya, while Southwestern China, and Qinling and Baishan Mountains lineages diverged during the Quaternary glaciations (about 1.13 Mya). Remnants of these once-distinct genetic clusters of J. regia may warrant ecological management if they are to be retained as in situ resources. A population size expansion in Northeastern China was detected in the last five centuries. The present distribution of walnut in China resulted from the combined effects of expansion/contraction from multiple refugia after the Last Glacial Maximum and later human exploitation.

High-altitude shorebird migration in the absence of topographical barriers: avoiding high air temperatures and searching for profitable winds

Nearly 20% of all bird species migrate between breeding and nonbreeding sites annually. Their migrations include storied feats of endurance and physiology, from non-stop trans-Pacific crossings to flights at the cruising altitudes of jetliners. Despite intense interest in these performances, there remains great uncertainty about which factors most directly influence bird behaviour during migratory flights. We used GPS trackers that measure an individual's altitude and wingbeat frequency to track the migration of black-tailed godwits (Limosa limosa) and identify the abiotic factors influencing their in-flight migratory behaviour. We found that godwits flew at altitudes above 5000 m during 21% of all migratory flights, and reached maximum flight altitudes of nearly 6000 m. The partial pressure of oxygen at these altitudes is less than 50% of that at sea level, yet these extremely high flights occurred in the absence of topographical barriers. Instead, they were associated with high air temperatures at lower altitudes and increasing wind support at higher altitudes. Our results therefore suggest that wind, temperature and topography all play a role in determining migratory behaviour, but that their relative importance is context dependent. Extremely high-altitude flights may thus not be especially rare, but they may only occur in very specific environmental contexts.

Selection and environmental adaptation along a path to speciation in the Tibetan frog Nanorana parkeri

Tibetan frogs, Nanorana parkeri, are differentiated genetically but not morphologically along geographical and elevational gradients in a challenging environment, presenting a unique opportunity to investigate processes leading to speciation. Analyses of whole genomes of 63 frogs reveal population structuring and historical demography, characterized by highly restricted gene flow in a narrow geographic zone lying between matrilines West (W) and East (E). A population found only along a single tributary of the Yalu Zangbu River has the mitogenome only of E, whereas nuclear genes of W comprise 89-95% of the nuclear genome. Selection accounts for 579 broadly scattered, highly divergent regions (HDRs) of the genome, which involve 365 genes. These genes fall into 51 gene ontology (GO) functional classes, 14 of which are likely to be important in driving reproductive isolation. GO enrichment analyses of E reveal many overrepresented functional categories associated with adaptation to high elevations, including blood circulation, response to hypoxia, and UV radiation. Four genes, including DNAJC8 in the brain, TNNC1 and ADORA1 in the heart, and LAMB3 in the lung, differ in levels of expression between low- and high-elevation populations. High-altitude adaptation plays an important role in maintaining and driving continuing divergence and reproductive isolation. Use of total genomes enabled recognition of selection and adaptation in and between populations, as well as documentation of evolution along a stepped cline toward speciation.

Cryptic genetic diversity in the mottled rabbitfish Siganus fuscescens with mitochondrial introgression at a contact zone in the South China Sea

The taxonomy of the mottled rabbitfish Siganus fuscescens species complex has long been challenging. In this study, we analyzed microsatellite genotypes, mitochondrial lineages, and morphometric data from 373 S. fuscescens individuals sampled from the northern Philippines and Hong Kong (South China Sea, Philippine Sea and Sulu Sea basins), to examine putative species boundaries in samples comprising three co-occurring mitochondrial lineages previously reported to characterize S. fuscescens (Clade A and Clade B) or S. canaliculatus (Clade C). We report the existence of two cryptic species within S. fuscescens in the northeast region of the South China Sea and northern Philippine Sea, supported by genetic and morphological differences. Individual-based assignment methods recovered concordant groupings of individuals into two nuclear genotype clusters (Cluster 1, Cluster 2) with (1) limited gene flow, if any, between them (FST = 0.241; P < 0.001); (2) low frequency of later-generation hybrids; (3) significant association with mitochondrial Clade A and Clade B, respectively; and (4) subtle yet significant body shape differences as inferred from geometric morphometric analysis. The divergence between mitochondrial Clade C and the two other clades was not matched by genetic differences at microsatellite marker loci. The occurrence of discordant mitonuclear combinations (20.5% of the total number of individuals) is thought to result from mitochondrial introgression, consistent with a scenario of demographic, and presumably spatial, post-Pleistocene expansion of populations from northern regions into a secondary contact zone in the South China Sea. Mitonuclear discordance due to introgression obscures phylogenetic relationships for recently-diverged lineages, and cautions against the use of mitochondrial markers alone for species identification within the mottled rabbitfish species complex in the South China Sea region.

Intraspecific variation in exploratory behavior and elevational affinity in a widely distributed songbird

Populations of the same species can vary substantially in their behavioral and morphometric traits when they are subject to different environmental pressures, which may lead to the development of different adaptive strategies. We quantified variation in exploratory behavior and morphometric traits among two rufous-collared sparrow populations that occur at low and high elevations in central Chile. Moreover, we used census and δ²H values of feather and blood to evaluate migration. We found that individual sparrows inhabiting high elevations were larger and showed more intense exploratory behavior in comparison with those that were captured at lower elevation. Moreover, we observed a steady decline in sparrow abundance during the winter and similar δ²H values for blood collected in the winter and summer at this site, which were significantly lower than blood δ²H values observed at low elevation. This pattern suggests that individuals do not move long distances during winter, and likely they remain at similar elevations in refuge habitats. As predicted, our results support the existent of different adaptive strategies among populations of the same species, and suggest that the combination of behavioral, morphometric, and stable isotope data is a novel and robust integrative approach to assess differences in adaptation across environmental gradients.

Landscape genomics: natural selection drives the evolution of mitogenome in penguins

Background

Mitochondria play a key role in the balance of energy and heat production, and therefore the mitochondrial genome is under natural selection by environmental temperature and food availability, since starvation can generate more efficient coupling of energy production. However, selection over mitochondrial DNA (mtDNA) genes has usually been evaluated at the population level. We sequenced by NGS 12 mitogenomes and with four published genomes, assessed genetic variation in ten penguin species distributed from the equator to Antarctica. Signatures of selection of 13 mitochondrial protein-coding genes were evaluated by comparing among species within and among genera (Spheniscus, Pygoscelis, Eudyptula, Eudyptes and Aptenodytes). The genetic data were correlated with environmental data obtained through remote sensing (sea surface temperature [SST], chlorophyll levels [Chl] and a combination of SST and Chl [COM]) through the distribution of these species.

Results

We identified the complete mtDNA genomes of several penguin species, including ND6 and 8 tRNAs on the light strand and 12 protein coding genes, 14 tRNAs and two rRNAs positioned on the heavy strand. The highest diversity was found in NADH dehydrogenase genes and the lowest in COX genes. The lowest evolutionary divergence among species was between Humboldt (Spheniscus humboldti) and Galapagos (S. mendiculus) penguins (0.004), while the highest was observed between little penguin (Eudyptula minor) and Adélie penguin (Pygoscelis adeliae) (0.097). We identified a signature of purifying selection (Ka/Ks < 1) across the mitochondrial genome, which is consistent with the hypothesis that purifying selection is constraining mitogenome evolution to maintain Oxidative phosphorylation (OXPHOS) proteins and functionality. Pairwise species maximum-likelihood analyses of selection at codon sites suggest positive selection has occurred on ATP8 (Fixed-Effects Likelihood, FEL) and ND4 (Single Likelihood Ancestral Counting, SLAC) in all penguins. In contrast, COX1 had a signature of strong negative selection. ND4 Ka/Ks ratios were highly correlated with SST (Mantel, p-value: 0.0001; GLM, p-value: 0.00001) and thus may be related to climate adaptation throughout penguin speciation.

Conclusions

These results identify mtDNA candidate genes under selection which could be involved in broad-scale adaptations of penguins to their environment. Such knowledge may be particularly useful for developing predictive models of how these species may respond to severe climatic changes in the future.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4424-9) contains supplementary material, which is available to authorized users.

Divergent clades or cryptic species? Mito-nuclear discordance in a Daphnia species complex

BACKGROUND

Genetically divergent cryptic species are frequently detected by molecular methods. These discoveries are often a byproduct of molecular barcoding studies in which fragments of a selected marker are used for species identification. Highly divergent mitochondrial lineages and putative cryptic species are even detected in intensively studied animal taxa, such as the crustacean genus Daphnia. Recently, eleven such lineages, exhibiting genetic distances comparable to levels observed among well-defined species, were recorded in the D. longispina species complex, a group that contains several key taxa of freshwater ecosystems. We tested if three of those lineages represent indeed distinct species, by analyzing patterns of variation of ten nuclear microsatellite markers in six populations.

RESULTS

We observed a discordant pattern between mitochondrial and nuclear DNA, as all individuals carrying one of the divergent mitochondrial lineages grouped at the nuclear level with widespread, well-recognized species coexisting at the same localities (Daphnia galeata, D. longispina, and D. cucullata).

CONCLUSIONS

A likely explanation for this pattern is the introgression of the mitochondrial genome of undescribed taxa into the common species, either in the distant past or after long-distance dispersal. The occurrence of highly divergent but rare mtDNA lineages in the gene pool of widespread species would suggest that hybridization and introgression in the D. longispina species complex is frequent even across large phylogenetic distances, and that discoveries of such distinct clades must be interpreted with caution. However, maintenance of ancient polymorphisms through selection is another plausible alternative that may cause the observed discordance and cannot be entirely excluded.

Signatures of mito-nuclear discordance in Schistosoma turkestanicum indicate a complex evolutionary history of emergence in Europe

High levels of molecular diversity were identified in mitochondrial cytochrome c oxidase (cox1) gene sequences of Schistosoma turkestanicum from Hungary. These cox1 sequences were all specific to Hungary which contrasted with the low levels of diversity seen in the nuclear internal transcribed spacer region (ITS) sequences, the majority of which were shared between China and Iran isolates. Measures of within and between host molecular variation within S. turkestanicum showed there to be substantial differences in molecular diversity, with cox1 being significantly more diverse than the ITS. Measures of haplotype frequencies revealed that each host contained its own subpopulation of genetically unique parasites with significant levels of differentiation. Pairwise mismatch analysis of cox1 sequences indicated S. turkestanicum populations to have a bimodal pairwise difference distribution and to be stable unlike the ITS sequences, which appeared to have undergone a recent population expansion event. Positive selection was also detected in the cox1 sequences, and biochemical modelling of the resulting protein illustrated significant mutational events causing an alteration to the isoelectric point of the cox1 protein, potentially altering metabolism. The evolutionary signature from the cox1 indicates local adaptation and long establishment of S. turkestanicum in Hungary with continual introgression of nuclear genes from Asian isolates. These processes have led to the occurrence of mito-nuclear discordance in a schistosome population.

Genetic population structure in an equatorial sparrow: roles for culture and geography

Female preference for local cultural traits has been proposed as a barrier to breeding among animal populations. As such, several studies have found correlations between male bird song dialects and population genetics over relatively large distances. To investigate whether female choice for local dialects could act as a barrier to breeding between nearby and contiguous populations, we tested whether variation in male song dialects explains genetic structure among eight populations of rufous-collared sparrows (Zonotrichia capensis) in Ecuador. Our study sites lay along a transect, and adjacent study sites were separated by approximately 25 km, an order of magnitude less than previously examined for this and most other species. This transect crossed an Andean ridge and through the Quijos River Valley, both of which may be barriers to gene flow. Using a variance partitioning approach, we show that song dialect is important in explaining population genetics, independent of the geographic variables: distance, the river valley and the Andean Ridge. This result is consistent with the hypothesis that song acts as a barrier to breeding among populations in close proximity. In addition, songs of contiguous populations differed by the same degree or more than between two populations previously shown to exhibit female preference for local dialect, suggesting that birds from these populations would also breed preferentially with locals. As expected, all geographic variables (distance, the river valley and the Andean Ridge) also predicted population genetic structure. Our results have important implications for the understanding whether, and at what spatial scale, culture can affect population divergence.

Digging up the roots of an insular hotspot of genetic diversity: decoupled mito-nuclear histories in the evolution of the Corsican-Sardinian endemic lizard Podarcis tiliguerta

BACKGROUND

Mediterranean islands host a disproportionately high level of biodiversity and endemisms. Growing phylogeographic evidence on island endemics has unveiled unexpectedly complex patterns of intra-island diversification, which originated at diverse spatial and temporal scales. We investigated multilocus genetic variation of the Corsican-Sardinian endemic lizard Podarcis tiliguerta with the aim of shedding more light on the evolutionary processes underlying the origin of Mediterranean island biodiversity.

RESULTS

We analysed DNA sequences of mitochondrial (12S and nd4) and nuclear (acm4 and mc1r) gene fragments in 174 individuals of P. tiliguerta from 81 localities across the full range of the species in a geographic and genealogical framework. We found surprisingly high genetic diversity both at mitochondrial and nuclear loci. Seventeen reciprocally monophyletic allopatric mitochondrial haplogroups were sharply divided into four main mitochondrial lineages (two in Corsica and two in Sardinia) of Miocene origin. In contrast, shallow divergence and shared diversity within and between islands was observed at the nuclear loci. We evaluated alternative biogeographic and evolutionary scenarios to explain such profound discordance in spatial and phylogenetic patterning between mitochondrial and nuclear genomes. While neutral models provided unparsimonious explanations for the observed pattern, the hypothesis of environmental selection driving mitochondrial divergence in the presence of nuclear gene flow is favoured.

CONCLUSIONS

Our study on the genetic variation of P. tiliguerta reveals surprising levels of diversity underlining a complex phylogeographic pattern with a striking example of mito-nuclear discordance. These findings have profound implications, not only for the taxonomy and conservation of P. tiliguerta. Growing evidence on deep mitochondrial breaks in absence of geographic barriers and of climatic factors associated to genetic variation of Corsican-Sardinian endemics warrants additional investigation on the potential role of environmental selection driving the evolution of diversity hotspots within Mediterranean islands.