Recent secondary contact, genome-wide admixture, and asymmetric introgression of neo-sex chromosomes between two Pacific island bird species

Secondary contact between closely related taxa represents a "moment of truth" for speciation-an opportunity to test the efficacy of reproductive isolation that evolved in allopatry and to identify the genetic, behavioral, and/or ecological barriers that separate species in sympatry. Sex chromosomes are known to rapidly accumulate differences between species, an effect that may be exacerbated for neo-sex chromosomes that are transitioning from autosomal to sex-specific inheritance. Here we report that, in the Solomon Islands, two closely related bird species in the honeyeater family-Myzomela cardinalis and Myzomela tristrami-carry neo-sex chromosomes and have come into recent secondary contact after ~1.1 my of geographic isolation. Hybrids of the two species were first observed in sympatry ~100 years ago. To determine the genetic consequences of hybridization, we use population genomic analyses of individuals sampled in allopatry and in sympatry to characterize gene flow in the contact zone. Using genome-wide estimates of diversity, differentiation, and divergence, we find that the degree and direction of introgression varies dramatically across the genome. For sympatric birds, autosomal introgression is bidirectional, with phenotypic hybrids and phenotypic parentals of both species showing admixed ancestry. In other regions of the genome, however, the story is different. While introgression on the Z/neo-Z-linked sequence is limited, introgression of W/neo-W regions and mitochondrial sequence (mtDNA) is highly asymmetric, moving only from the invading M. cardinalis to the resident M. tristrami. The recent hybridization between these species has thus enabled gene flow in some genomic regions but the interaction of admixture, asymmetric mate choice, and/or natural selection has led to the variation in the amount and direction of gene flow at sex-linked regions of the genome.

Flexibility underlies differences in mitochondrial respiratory performance between migratory and non-migratory White-crowned Sparrows (Zonotrichia leucophrys)

Migration is one of the most energy-demanding behaviors observed in birds. Mitochondria are the primary source of energy used to support these long-distance movements, yet how mitochondria meet the energetic demands of migration is scarcely studied. We quantified changes in mitochondrial respiratory performance in the White-crowned Sparrow (Zonotrichia leucophrys), which has a migratory and non-migratory subspecies. We hypothesized that the long-distance migratory Gambel's subspecies (Z. l. gambelii) would show higher mitochondrial respiratory performance compared to the non-migratory Nuttall's subspecies (Z. l. nuttalli). We sampled Gambel's individuals during spring pre-migration, active fall migration, and a period with no migration or breeding (winter). We sampled Nuttall's individuals during periods coinciding with fall migration and the winter period of Gambel's annual cycle. Overall, Gambel's individuals had higher citrate synthase, a proxy for mitochondrial volume, than Nuttall's individuals. This was most pronounced prior to and during migration. We found that both OXPHOS capacity (state 3) and basal respiration (state 4) of mitochondria exhibit high seasonal flexibility within Gambel's individuals, with values highest during active migration. These values in Nuttall's individuals were most similar to Gambel's individuals in winter. Our observations indicate that seasonal changes in mitochondrial respiration play a vital role in migration energetics.

Migration increases mitochondrial oxidative capacity without increasing reactive oxygen species emission in a songbird

Birds remodel their flight muscle metabolism prior to migration to meet the physiological demands of migratory flight, including increases in both oxidative capacity and defence against reactive oxygen species. The degree of plasticity mediated by changes in these mitochondrial properties is poorly understood but may be explained by two non-mutually exclusive hypotheses: variation in mitochondrial quantity or in individual mitochondrial function. We tested these hypotheses using yellow-rumped warblers (Setophaga coronata), a Nearctic songbird which biannually migrates 2000-5000 km. We predicted higher flight muscle mitochondrial abundance and substrate oxidative capacity, and decreased reactive oxygen species emission in migratory warblers captured during autumn migration compared with a short-day photoperiod-induced non-migratory phenotype. We assessed mitochondrial abundance via citrate synthase activity and assessed isolated mitochondrial function using high-resolution fluororespirometry. We found 60% higher tissue citrate synthase activity in the migratory phenotype, indicating higher mitochondrial abundance. We also found 70% higher State 3 respiration (expressed per unit citrate synthase) in mitochondria from migratory warblers when oxidizing palmitoylcarnitine, but similar H2O2 emission rates between phenotypes. By contrast, non-phosphorylating respiration was higher and H2O2 emission rates were lower in the migratory phenotype. However, flux through electron transport system complexes I-IV, II-IV and IV was similar between phenotypes. In support of our hypotheses, these data suggest that flight muscle mitochondrial abundance and function are seasonally remodelled in migratory songbirds to increase tissue oxidative capacity without increasing reactive oxygen species formation.

Conservation Mitonuclear Replacement: Facilitated mitochondrial adaptation for a changing world

Most species will not be able to migrate fast enough to cope with climate change, nor evolve quickly enough with current levels of genetic variation. Exacerbating the problem are anthropogenic influences on adaptive potential, including the prevention of gene flow through habitat fragmentation and the erosion of genetic diversity in small, bottlenecked populations. Facilitated adaptation, or assisted evolution, offers a way to augment adaptive genetic variation via artificial selection, induced hybridization, or genetic engineering. One key source of genetic variation, particularly for climatic adaptation, are the core metabolic genes encoded by the mitochondrial genome. These genes influence environmental tolerance to heat, drought, and hypoxia, but must interact intimately and co-evolve with a suite of important nuclear genes. These coadapted mitonuclear genes form some of the important reproductive barriers between species. Mitochondrial genomes can and do introgress between species in an adaptive manner, and they may co-introgress with nuclear genes important for maintaining mitonuclear compatibility. Managers should consider the relevance of mitonuclear genetic variability in conservation decision-making, including as a tool for facilitating adaptation. I propose a novel technique dubbed Conservation Mitonuclear Replacement (CmNR), which entails replacing the core metabolic machinery of a threatened species-the mitochondrial genome and key nuclear loci-with those from a closely related species or a divergent population, which may be better-adapted to climatic changes or carry a lower genetic load. The most feasible route to CmNR is to combine CRISPR-based nuclear genetic editing with mitochondrial replacement and assisted reproductive technologies. This method preserves much of an organism's phenotype and could allow populations to persist in the wild when no other suitable conservation options exist. The technique could be particularly important on mountaintops, where rising temperatures threaten an alarming number of species with almost certain extinction in the next century.

High-altitude adaptation is accompanied by strong signatures of purifying selection in the mitochondrial genomes of three Andean waterfowl

Evidence from a variety of organisms points to convergent evolution on the mitochondria associated with a physiological response to oxygen deprivation or temperature stress, including mechanisms for high-altitude adaptation. Here, we examine whether demography and/or selection explains standing mitogenome nucleotide diversity in high-altitude adapted populations of three Andean waterfowl species: yellow-billed pintail (Anas georgica), speckled teal (Anas flavirostris), and cinnamon teal (Spatula cyanoptera). We compared a total of 60 mitogenomes from each of these three duck species (n = 20 per species) across low and high altitudes and tested whether part(s) or all of the mitogenome exhibited expected signatures of purifying selection within the high-altitude populations of these species. Historical effective population sizes (Ne) were inferred to be similar between high- and low-altitude populations of each species, suggesting that selection rather than genetic drift best explains the reduced genetic variation found in mitochondrial genes of high-altitude populations compared to low-altitude populations of the same species. Specifically, we provide evidence that establishment of these three Andean waterfowl species in the high-altitude environment, coincided at least in part with a persistent pattern of negative purifying selection acting on oxidative phosphorylation (OXPHOS) function of the mitochondria. Our results further reveal that the extent of gene-specific purifying selection has been greatest in the speckled teal, the species with the longest history of high-altitude occupancy.

No Signs of Adaptations for High Flight Intensity in the Mitochondrial Genome of Birds

Mitochondrial genomes are expected to show adaptations for optimizing aerobic respiration in birds that make intense use of flight. However, there is limited empirical evidence of such a relationship. We here examine correlates of several mitochondrial genome characteristics and flight use across a diverse sample of 597 bird species. We developed an index of flight use intensity that ranged from 0 in flightless species to 9 in migratory hummingbirds and examined its association with nucleobase composition, amino acid class composition, and amino acid site allelic variation using phylogenetic comparative methods. We found no evidence of mitochondrial genome adaptations to flight intensity. Neither nucleotide composition nor amino acid properties showed consistent patterns related to flight use. While specific sites in mitochondrial genes exhibited variation associated with flight intensity, there was limited association between specific amino acid residues and flight intensity levels. Our findings suggest a complex genetic architecture for aerobic performance traits, where multiple genes in both mitochondria and the nucleus may contribute to overall performance. Other factors, such as gene expression regulation and anatomical adaptations, may play a more significant role in influencing flight performance than changes in mitochondrial genes. These findings highlight the need for comprehensive genomic analyses to unravel the intricate relationship between genetic variants and aerobic performance in birds.

Long-duration wind tunnel flights reveal exponential declines in protein catabolism over time in short- and long-distance migratory warblers

During migration, long-distance migratory songbirds may fly nonstop for days, whereas shorter-distance migrants complete flights of 6 to 10 h. Fat is the primary fuel source, but protein is also assumed to provide a low, consistent amount of energy for flight. However, little is known about how the use of these fuel sources differs among bird species and in response to flight duration. Current models predict that birds can fly until fat stores are exhausted, with little consideration of protein's limits on flight range or duration. We captured two related migratory species-ultra long-distance blackpoll warblers (Setophaga striata) and short-distance yellow-rumped warblers (Setophaga coronata)-during fall migration and flew them in a wind tunnel to examine differences in energy expenditure, overall fuel use, and fuel mixture. We measured fat and fat-free body mass before and after flight using quantitative magnetic resonance and calculated energy expenditure from body composition changes and doubly labeled water. Three blackpolls flew voluntarily for up to 28 h-the longest wind tunnel flight to date-and ended flights with substantial fat reserves but concave flight muscle, indicating that protein loss, rather than fat, may actually limit flight duration. Interestingly, while blackpolls had significantly lower mass-specific metabolic power in flight than that of yellow-rumped warblers and fuel use was remarkably similar in both species, with consistent fat use but exceptionally high rates of protein loss at the start of flight that declined exponentially over time. This suggests that protein may be a critical, dynamic, and often overlooked fuel for long-distance migratory birds.

Short- and long-distance avian migrants differ in exercise endurance but not aerobic capacity

Background

Migratory birds differ markedly in their migration strategies, particularly those performing short- versus long-distance migrations. In preparation for migration, all birds undergo physiological and morphological modifications including enlargement of fat stores and pectoral muscles to fuel and power their flights, as well as cardiovascular and biochemical adjustments that improve lipid and oxygen delivery and uptake by flight muscles. While the magnitude of these changes varies in relation to migration strategy, the consequence of these variations on aerobic performance is unknown. We tested whether the aerobic performance of four Old-world flycatcher species (Muscicapidae) varied according to migration strategy by comparing minimum resting metabolic rates (RMRmin), exercise-induced maximum metabolic rates (MMR), and exercise endurance times of short-distance and long-distance migratory birds.

Results

As expected, RMRmin did not vary between short-distance and long-distance migrants but differed between the species within a migration strategy and between sexes. Unexpectedly, MMR did not vary with migration strategy, but MMR and blood haemoglobin content were positively related among the birds tested. Exercise endurance times differed substantially between migration strategies with long-distance migrants sustaining exercise for > 60% longer than short-distance migrants. Blood haemoglobin content had a significant positive effect on endurance among all birds examined.

Conclusions

The lack of difference in RMRmin and MMR between long- and short-distance migrants during this stage of migration suggests that the attributes favouring the greater aerobic endurance of long-distance migrants did not come at the expense of increased maintenance costs or require greater aerobic capacity.

Do habitat and elevation promote hybridization during secondary contact between three genetically distinct groups of warbling vireo (Vireo gilvus)?

Following postglacial expansion, secondary contact can occur between genetically distinct lineages. These genetic lineages may be associated with specific habitat or environmental variables and therefore, their distributions in secondary contact could reflect such conditions within these areas. Here we used mtDNA, microsatellite, and morphological data to study three genetically distinct groups of warbling vireo (Vireo gilvus) and investigate the role that elevation and habitat play in their distributions. We studied two main contact zones and within each contact zone, we examined two separate transects. Across the Great Plains contact zone, we found that hybridization between eastern and western groups occurs along a habitat and elevational gradient, whereas hybridization across the Rocky Mountain contact zone was not as closely associated with habitat or elevation. Hybrids in the Great Plains contact zone were more common in transitional areas between deciduous and mixed-wood forests, and at lower elevations (<1000 m). Hybridization patterns were similar along both Great Plains transects indicating that habitat and elevation play a role in hybridization between distinct eastern and western genetic groups. The observed patterns suggest adaptation to different habitats, perhaps originating during isolation in multiple Pleistocene refugia, is facilitating hybridization in areas where habitat types overlap.

Fat Stores and Antioxidant Capacity Affect Stopover Decisions in Three of Four Species of Migratory Passerines With Different Migration Strategies: An Experimental Approach

During migratory stopovers, birds must make decisions about when and where to travel and these decisions are likely contingent on their fuel stores, food availability, and antioxidant capacity as well as seasonal changes in key environmental factors. We conducted a field experiment on an offshore stopover site (Block Island, Rhode Island, United States: 41°130N, 71°330W) during autumn migration to test the hypothesis that birds with greater fuel stores and non-enzymatic antioxidant capacity have shorter stopovers than lean birds with low antioxidant capacity, and to determine the extent to which this depends on migration strategy. We used a 2 × 2 factorial field experiment (two levels each of available food and dietary polyphenols) with four species of songbirds kept in captivity for 3–5 days to produce experimental groups with different fuel stores and antioxidant capacity. We attached digital VHF transmitters to assess stopover duration and departure direction using automated telemetry. Non-enzymatic antioxidant capacity increased during refueling for Red-eyed Vireos (Vireo olivaceus) and Blackpoll Warblers (Setophaga striata) fed ad lib diets, and for ad lib fed Hermit Thrushes (Catharus guttatus) supplemented with polyphenols, but not for Yellow-rumped Warblers (Setophaga coronata coronata). Glutathione peroxidase (GPx) decreased during captivity and was influenced by dietary treatment only in Red-eyed Vireos. Oxidative damage decreased during captivity for all species except Yellow-rumped Warblers. Stopover duration was shorter for Vireos and Blackpolls fed ad lib as compared to those fed maintenance. Ad lib fed Hermit Thrushes supplemented with polyphenols had shorter stopovers than those fed ad lib, as did thrushes fed at maintenance and supplemented with polyphenols compared with those fed at maintenance alone. There was no influence of condition on stopover duration for Yellow-rumped Warblers. Departure direction was not strongly related to condition, and birds primarily reoriented north when departing Block Island. Thus, fat stores and oxidative status interacted to influence the time passerines spent on stopover, and condition-dependent departure decisions were related to a bird’s migration strategy. Therefore, seasonal variation in macro- and micro-nutrient resources available for refueling at stopover sites can affect body condition and antioxidant capacity and in turn influence the timing and success of migration.

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.

Cytonuclear discordance in the crowned-sparrows, Zonotrichia atricapilla and Zonotrichia leucophrys

The golden-crowned (Zonotrichia atricapilla) and white-crowned (Z. leucophrys) sparrows have been presented as a compelling case for rapid speciation. They display divergence in song and plumage with overlap in their breeding ranges implying reproductive isolation, but have almost identical mitochondrial genomes. Previous research proposed hybridization and subsequent mitochondrial introgression as an alternate explanation, but lacked robust nuclear gene trees to distinguish between introgression and incomplete lineage sorting. We test for signatures of these processes between Z. atricapilla and Z. leucophrys, and investigate the relationships among Z. leucophrys subspecies, using mitochondrial sequencing and a reduced representation nuclear genomic dataset. Contrary to the paraphyly evident in mitochondrial gene trees, we confirmed the reciprocal monophyly of Z. atricapilla and Z. leucophrys using large panels of single nucleotide polymorphisms (SNPs). The pattern of cytonuclear discordance is consistent with limited, historical hybridization and mitochondrial introgression, rather than a recent origin and incomplete lineage sorting between recent sister species. We found evidence of nuclear phylogeographic structure within Z. leucophrys with two distinct clades. Altogether, our results indicate deeper divergences between Z. atricapilla and Z. leucophrys than inferred using mitochondrial markers. Our results demonstrate the limitations of relying solely on mitochondrial DNA for taxonomy, and raise questions about the possibility of selection on the mitochondrial genome during temperature oscillations (e.g. during the Pleistocene). Historical mitochondrial introgression facilitated by past environmental changes could cause erroneous dating of lineage splitting in other taxa when based on mitochondrial DNA alone.

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.

Comparing divergence landscapes from reduced-representation and whole genome resequencing in the yellow-rumped warbler (Setophaga coronata) species complex

Researchers seeking to generate genomic data for non-model organisms are faced with a number of trade-offs when deciding which method to use. The selection of reduced representation approaches versus whole genome resequencing will ultimately affect the marker density, sequencing depth, and the number of individuals that can multiplexed. These factors can affect researchers' ability to accurately characterize certain genomic features, such as landscapes of divergence-how F_ST varies across the genomes. To provide insight into the effect of sequencing method on the estimation of divergence landscapes, we applied an identical bioinformatic pipeline to three generations of sequencing data (GBS, ddRAD, and WGS) produced for the same system, the yellow-rumped warbler species complex. We compare divergence landscapes generated using each method for the myrtle warbler (Setophaga coronata coronata) and the Audubon's warbler (S. c. auduboni), and for Audubon's warblers with deeply divergent mtDNA resulting from mitochondrial introgression. We found that most high-F_ST peaks were not detected in the ddRAD data set, and that while both GBS and WGS were able to identify the presence of large peaks, WGS was superior at a finer scale. Comparing Audubon's warblers with divergent mitochondrial haplotypes, only WGS allowed us to identify small (10-20 kb) regions of elevated differentiation, one of which contained the nuclear-encoded mitochondrial gene NDUFAF3. We calculated the cost per base pair for each method and found it was comparable between GBS and WGS, but significantly higher for ddRAD. These comparisons highlight the advantages of WGS over reduced representation methods when characterizing landscapes of divergence.

Geographical disjunction and environmental conditions drive intraspecific differentiation in the chalk-hill blue butterfly

Drivers of evolution are often related to geographical isolation and/or diverging environmental conditions. Spatial variation in neutral genetic markers mostly reflects past geographical isolation, i.e. long-lasting allopatry, whereas morphology is often driven by local environmental conditions, resulting in more rapid evolution. In Europe, most thermophilic species persisted during the past glacial periods in geographically disjunct refugia, representing long-lasting isolates, frequently with diverging environmental conditions. This situation has driven the evolution of intraspecific signatures in species. Here, we analysed wing shape and wing pigmentation of the chalk-hill blue butterfly, Polyommatus coridon, across its entire distribution range restricted to the western Palaearctic. In addition, we compiled abiotic environmental parameters for each sampling site. Wing colour patterns differentiated a western and an eastern lineage. These lineages might represent two main Pleistocene refugia and differentiation centres, one located on the Italian Peninsula and the other in the Balkan region. The two lineages showed evidence of hybridization across Central Europe, from the Alps and across Germany. The intraspecific differentiation was strongest in the width of the brown band on the outer margin of the wings. The morphological structures obtained are in line with genetic signatures found in previous studies, but the latter are more fine-grained. Current environmental conditions, such as mean temperatures, were only marginally correlated with colour patterns. Our study underlines that Pleistocene range shifts, often resulting in allopatric isolation, shape intraspecific phenotypic structures within species; that pigmentation responds in a more sensitive manner to spatial disjunction than wing shape; and that morphometric and genetic structures in P. coridon provide concordant patterns and thus support identical biogeographical conclusions.

Complex histories of gene flow and a mitochondrial capture event in a nonsister pair of birds

Hybridization, introgression, and reciprocal gene flow during speciation, specifically the generation of mitonuclear discordance, are increasingly observed as parts of the speciation process. Genomic approaches provide insight into where, when, and how adaptation operates during and after speciation and can measure historical and modern introgression. Whether adaptive or neutral in origin, hybridization can cause mitonuclear discordance by placing the mitochondrial genome of one species (or population) in the nuclear background of another species. The latter, introgressed species may eventually have its own mtDNA replaced or “captured” by other species across its entire geographical range. Intermediate stages in the capture process should be observable. Two nonsister species of Australasian monarch‐flycatchers, Spectacled Monarch (Symposiachrus trivirgatus) mostly of Australia and Indonesia and Spot‐winged Monarch (S. guttula) of New Guinea, present an opportunity to observe this process. We analysed thousands of single nucleotide polymorphisms (SNPs) derived from ultraconserved elements of all subspecies of both species. Mitochondrial DNA sequences of Australian populations of S. trivirgatus form two paraphyletic clades, one being sister to and presumably introgressed by S. guttula despite little nuclear signal of introgression. Population genetic analyses (e.g., tests for modern and historical gene flow and selection) support at least one historical gene flow event between S. guttula and Australian S. trivirgatus. We also uncovered introgression from the Maluku Islands subspecies of S. trivirgatus into an island population of S. guttula, resulting in apparent nuclear paraphyly. We find that neutral demographic processes, not adaptive introgression, are the most likely cause of these complex population histories. We suggest that a Pleistocene extinction of S. guttula from mainland Australia resulted from range expansion by S. trivirgatus.

Asymmetric introgression reveals the genetic architecture of a plumage trait

Genome-wide variation in introgression rates across hybrid zones offers a powerful opportunity for studying population differentiation. One poorly understood pattern of introgression is the geographic displacement of a trait implicated in lineage divergence from genome-wide population boundaries. While difficult to interpret, this pattern can facilitate the dissection of trait genetic architecture because traits become uncoupled from their ancestral genomic background. We studied an example of trait displacement generated by the introgression of head plumage coloration from personata to alba subspecies of the white wagtail. A previous study of their hybrid zone in Siberia revealed that the geographic transition in this sexual signal that mediates assortative mating was offset from other traits and genetic markers. Here we show that head plumage is associated with two small genetic regions. Despite having a simple genetic architecture, head plumage inheritance is consistent with partial dominance and epistasis, which could contribute to its asymmetric introgression.

Hybrid zones are windows into the evolutionary process. Semenov et al. find that the head plumage differences between white wagtail subspecies have a simple genetic basis involving two small genetic regions, in which partially dominant and epistatic interactions help to explain how this sexual signal has become decoupled from other plumage traits.

Complete mitochondrial genome of Chroicocephalus brunnicephalus from India: phylogeny with other Larids

The complete mitogenome sequence of the brown-headed gull, Chroicocephalus brunnicephalus was determined in this study. The 16,771 bp genome consists of 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes, and a control region (CR). The decoded mitogenome was AT-rich (54.77%) with nine overlapping and 17 intergenic spacer regions. Most of the PCGs were started by a typical ATG initiation codon except for cox1 and nad3. Further, the usual termination codons (AGG, TAG, TAA, and AGA) were used by 11 PCGs except for cox3 and nad4. The concatenated PCGs based Bayesian phylogeny clearly discriminates all the Laridae species and reflects the sister relationship of C. brunnicephalus with C. ridibundus. The present mitogenome-based phylogeny was congruent with the earlier hypothesis and confirmed the evolutionary position of the brown-headed gull as masked species. The generated mitogenome of C. brunnicephalus is almost identical to the previously generated mitogenome from China except for two base pairs in CR. To visualize the population structure of this migratory species, we propose more sampling from different geographical locations and the generation of additional molecular data to clarify the reality.

Elevation as a selective force on mitochondrial respiratory chain complexes of the Phrynocephalus lizards in the Tibetan plateau

Animals living in extremely high elevations have to adapt to low temperatures and low oxygen availability (hypoxia), but the underlying genetic mechanisms associated with these adaptations are still unclear. The mitochondrial respiratory chain can provide >95% of the ATP in animal cells, and its efficiency is influenced by temperature and oxygen availability. Therefore, the respiratory chain complexes (RCCs) could be important molecular targets for positive selection associated with respiratory adaptation in high-altitude environments. Here, we investigated positive selection in 5 RCCs and their assembly factors by analyzing sequences of 106 genes obtained through RNA-seq of all 15 Chinese Phrynocephalus lizard species, which are distributed from lowlands to the Tibetan plateau (average elevation >4,500 m). Our results indicate that evidence of positive selection on RCC genes is not significantly different from assembly factors, and we found no difference in selective pressures among the 5 complexes. We specifically looked for positive selection in lineages where changes in habitat elevation happened. The group of lineages evolving from low to high altitude show stronger signals of positive selection than lineages evolving from high to low elevations. Lineages evolving from low to high elevation also have more shared codons under positive selection, though the changes are not equivalent at the amino acid level. This study advances our understanding of the genetic basis of animal respiratory metabolism evolution in extreme high environments and provides candidate genes for further confirmation with functional analyses.

Introgression obscures lineage boundaries and phylogeographic history in the western banded gecko, Coleonyx variegatus (Squamata: Eublepharidae)

The geomorphological formation of the Baja California peninsula and the Gulf of California is a principal driver of diversification for the reptiles of North America’s warm deserts. The western banded gecko, Coleonyx variegatus, is distributed throughout the Mojave, Sonoran and Peninsular deserts. In this study we use multilocus sequence data to address deep phylogeographic structure within C. variegatus. Analyses of mtDNA data recover six divergent clades throughout the range of C. variegatus. Topology of the mtDNA gene tree suggests separate origins of peninsular populations with an older lineage in the south and a younger one in the north. In contrast, analyses of multilocus nuclear data provide support for four lineages, corresponding to the subspecies C. v. abbotti, C. v. peninsularis, C. v. sonoriensis and C. v. variegatus. Phylogenetic analyses of the nuclear data recover C. v. abbotti and C. v. peninsularis as a clade, indicating a single origin of the peninsular populations. Discordance between the nuclear and mtDNA data is largely the result of repeated episodes of mtDNA introgression that have obscured both lineage boundaries and biogeographic history. Dating analyses of the combined nuclear and mtDNA data suggest that the peninsular clade diverged from the continental group in the Late Miocene.

Robbery in progress: Historical museum collections bring to light a mitochondrial capture within a bird species widespread across southern Australia, the Copperback Quail-thrush Cinclosoma clarum

We surveyed mitochondrial, autosomal, and Z chromosome diversity within and between the Copperback Quail-thrush Cinclosoma clarum and Chestnut Quail-thrush C. castanotum, which together span the arid and semi-arid zones of southern Australia, and primarily from specimens held in museum collections. We affirm the recent taxonomic separation of the two species and then focus on diversity within the more widespread of the two species, C. clarum. To guide further study of the system and what it offers to understanding the genomics of the differentiation and speciation processes, we develop and present a hypothesis to explain mitonuclear discordance that emerged in ourdata. Following a period of historical allopatry, secondary contact has resulted in an eastern mitochondrial genome replacing the western mitochondrial genome in western populations. This is predicted under a population-level invasion in the opposite direction, that of the western population invading the range of the eastern one. Mitochondrial captures can be driven by neutral, demographic processes, or adaptive mechanisms, and we favor the hypothesized capture being driven by neutral means. We cannot fully reject the adaptive process but suggest how these alternatives may be further tested. We acknowledge an alternative hypothesis, which finds some support in phenotypic data published elsewhere, namely that outcomes of secondary contact have been more complex than our current genomic data suggest. Discriminating and reconciling these two alternative hypotheses, which may not be mutually exclusive, could be tested with closer sampling at levels of population, individual, and nucleotide than has so far been possible. This would be further aided by knowledge of the genetic basis to phenotypic variation described elsewhere.

Genome-wide data reveal discordant mitonuclear introgression in the intermediate horseshoe bat (Rhinolophus affinis)

Closely related taxa often exhibit mitonuclear discordance attributed to introgression of mitochondrial DNA (mtDNA), yet few studies have considered the underlying causes of mtDNA introgression. Here we test for demographic versus adaptive processes as explanations for mtDNA introgression in three subspecies of the intermediate horseshoe bat (Rhinolophus affinis). We generated sequences of 1692 nuclear genes and 13 mitochondrial protein-coding genes for 48 individuals. Phylogenetic reconstructions based on 320 exon sequences and 2217 single nucleotide polymorphisms (SNPs) both revealed conflicts between the species tree and mtDNA tree. These results, together with geographic patterns of mitonuclear discordance, and shared identical or near-identical mtDNA sequences, suggest extensive introgression of mtDNA between the two parapatric mainland subspecies. Under demographic hypotheses, we would also expect to uncover traces of ncDNA introgression, however, population structure and gene flow analyses revealed little nuclear admixture. Furthermore, we found inconsistent estimates of the timing of population expansion and that of the most recent common ancestor for the clade containing introgressed haplotypes. Without a clear demographic explanation, we also examined whether introgression likely arises from adaptation. We found that five mtDNA genes contained fixed amino acid differences between introgressed and non-introgressed individuals, including putative positive selection found in one codon, although this did not show introgression. While our evidence for rejecting demographic hypotheses is arguably stronger than that for rejecting adaptation, we find no definitive support for either explanation. Future efforts will focus on larger-scale resequencing to decipher the underlying causes of discordant mitonuclear introgression in this system.

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.

How long is 3 km for a butterfly? Ecological constraints and functional traits explain high mitochondrial genetic diversity between Sicily and the Italian Peninsula

Populations inhabiting Mediterranean islands often show contrasting genetic lineages, even on islands that were connected to the mainland during glacial maxima. This pattern is generated by forces acting in historical and contemporary times. Understanding these phenomena requires comparative studies integrating genetic structure, functional traits and dispersal constraints. Using as a model the butterfly species living across the Messina strait (3 km wide) separating Sicily from the Italian Peninsula, we aimed to unravel the mechanisms limiting the dispersal of matrilines and generating genetic differentiation across a narrow sea strait. We analysed the mitochondrial COI gene of 84 butterfly species out of 90 documented in Sicily and compared them with populations from the neighbouring southern Italian Peninsula (1,398 sequences) and from the entire Palaearctic region (8,093 sequences). For each species, we regressed 13 functional traits and 2 ecological constraints to dispersal (winds experienced at the strait and climatic suitability) against genetic differentiation between Sicily and Italian Peninsula to understand the factors limiting dispersal. More than a third of the species showed different haplogroups across the strait and most of them also represented endemic haplogroups for this island. One fifth of Sicilian populations (and 32.3% of endemic lineages) had their closest relatives in distant areas, instead of the neighbouring Italian Peninsula, which suggests high relictuality. Haplotype diversity was significantly explained by the length of the flight period, an intrinsic phenology trait, while genetic differentiation was explained by both intrinsic traits (wingspan and degree of generalism) and contemporary local constraints (winds experienced at the strait and climatic suitability). A relatively narrow sea strait can produce considerable differentiation among butterfly matrilines and this phenomenon showed a largely deterministic fingerprint. Because of unfavourable winds, populations of the less dispersive Sicilian butterflies tended to differentiate into endemic variants or to maintain relict populations. Understanding these phenomena required the integration of DNA sequences, species traits and physical constraints for a large taxon at continental scale. Future studies may reveal if the patterns here shown for mitochondrial DNA are also reflected in the nuclear genome or, alternatively, are the product of limited female dispersal.

Effects of mitonuclear combination and thermal acclimation on the energetic phenotype

Activity of the oxidative phosphorylation complexes rely on intimately associated subunits encoded by the mitochondrial and nuclear genomes. Given the key role of this system in adenosine triphosphate production, genes from both genomes must coevolve. A combination of northern redbelly dace (Chrosomus eos) or finescale dace (C. neogaeus) mitochondrial genome with a C. eos nuclear genome allows for a close examination of a naturally occurring disruption of mitonuclear coevolution. We, therefore, investigated the combined effect of mitonuclear genotypes, acclimation, and temperature on the activity of enzymes linked with the energy metabolism in a sympatric population of wild type and cybrid. As expected, the activity of the nuclear-encoded citrate synthase was only influenced by temperature while the cytochrome c oxidase (composed of nuclear and mitochondrial subunits from wild type and cybrid individuals) responded differently to temperature. This study provides clear evidence of the extent by which mitonuclear coadaptation could influence aerobic metabolism.

Lack of gene flow between Phytophthora infestans populations of two neighboring countries with the largest potato production

Gene flow is an important evolutionary force that enables adaptive responses of plant pathogens in response to changes in the environment and plant disease management strategies. In this study, we made a direct inference concerning gene flow in the Irish famine pathogen Phytophthora infestans between two of its hosts (potato and tomato) as well as between China and India. This was done by comparing sequence characteristics of the eukaryotic translation elongation factor 1 alpha (eEF-1α) gene, generated from 245 P. infestans isolates sampled from two countries and hosts. Consistent with previous results, we found that eEF-1α gene was highly conserved and point mutation was the only mechanism generating any sequence variation. Higher genetic variation was found in the eEF-1α sequences in the P. infestans populations sampled from tomato compared to those sampled from potato. We also found the P. infestans population from India displayed a higher genetic variation in the eEF-1α sequences compared to China. No gene flow was detected between the pathogen populations from the two countries, which is possibly attributed to the geographic barrier caused by Himalaya Plateau and the minimum cross-border trade of potato and tomato products. The implications of these results for a sustainable management of late blight diseases are discussed.

Plastic but repeatable: rapid adjustments of mitochondrial function and density during reproduction in a wild bird species

Most of the energy fluxes supporting animal performance flow through mitochondria. Hence, inter-individual differences in performance might be rooted in inter-individual variations in mitochondrial function and density. Furthermore, because the energy required by an individual often changes across life stages, mitochondrial function and density are also expected to show within-individual variation (i.e. plasticity). No study so far has repeatedly measured mitochondrial function and density in the same individuals to simultaneously test for within-individual repeatability and plasticity of mitochondrial traits. Here, we repeatedly measured mitochondrial DNA copy number (a proxy of density) and respiration rates from blood cells of female pied flycatchers (Ficedula hypoleuca) at the incubation and chick-rearing stages. Mitochondrial density and respiration rates were all repeatable (R = [0.45; 0.80]), indicating high within-individual consistency in mitochondrial traits across life-history stages. Mitochondrial traits were also plastic, showing a quick (i.e. 10 days) downregulation from incubation to chick-rearing in mitochondrial density, respiratory activity, and cellular regulation by endogenous substrates and/or ATP demand. These downregulations were partially compensated by an increase in mitochondrial efficiency at the chick-rearing stage. Therefore, our study provides clear evidence for both short-term plasticity and high within-individual consistency in mitochondrial function and density during reproduction in a wild bird species.

When barcoding fails: development of diagnostic nuclear markers for the sibling caddisfly species Sericostoma personatum (Spence in Kirby & Spence, 1826) and Sericostoma flavicorne Schneider, 1845

The larval stages of the central European sibling caddisfly species Sericostomapersonatum (Spence in Kirby and Spence, 1826) and S.flavicorne Schneider, 1845 are morphologically similar and can only be distinguished by differences in coloration in late larval instars. Identification using the mitochondrial barcoding gene, i.e., the Cytochrome c Oxidase 1, is impossible, as both species share the same highly differentiated haplotypes due to introgression. Nuclear gene markers obtained through double digest restriction site associate sequencing (ddRAD seq), however, can reliably distinguish both species, yet the method is expensive as well as time-consuming and therefore not practicable for species determination. To facilitate accurate species identification without sequencing genome-wide markers, we developed nine diagnostic nuclear RFLP markers based on ddRAD seq data. The markers were successfully tested on geographically distinct populations of the two Sericostoma species in western Germany, on known hybrids, and on another sericostomatid caddisfly species, Oecismusmonedula (Hagen, 1859) that sometimes shares the habitat and can be morphologically confounded with Sericostoma. We describe a simple and fast protocol for reliable species identification of S.personatum and S.flavicorne independent of the life cycle stage of the specimens.

Mitochondrial variation in small brown planthoppers linked to multiple traits and probably reflecting a complex evolutionary trajectory

While it has been proposed in several taxa that the mitochondrial genome is associated with adaptive evolution to different climatic conditions, making links between mitochondrial haplotypes and organismal phenotypes remains a challenge. Mitonuclear discordance occurs in the small brown planthopper (SBPH), Laodelphax striatellus, with one mitochondrial haplogroup (HGI) more common in the cold climate region of China relative to another form (HGII) despite strong nuclear gene flow, providing a promising model to investigate climatic adaptation of mitochondrial genomes. We hypothesized that cold adaptation through HGI may be involved, and considered mitogenome evolution, population genetic analyses, and bioassays to test this hypothesis. In contrast to our hypothesis, chill-coma recovery tests and population genetic tests of selection both pointed to HGII being involved in cold adaptation. Phylogenetic analyses revealed that HGII is nested within HGI, and has three nonsynonymous changes in ND2, ND5 and CYTB in comparison to HGI. These molecular changes likely increased mtDNA copy number, cold tolerance and fecundity of SBPH, particularly through a function-altering amino acid change involving M114T in ND2. Nuclear background also influenced fecundity and chill recovery (i.e., mitonuclear epistasis) and protein modelling indicates possible nuclear interactions for the two nonsynonymous changes in ND2 and CYTB. The high occurrence frequency of HGI in the cold climate region of China remains unexplained, but several possible reasons are discussed. Overall, our study points to a link between mtDNA variation and organismal-level evolution and suggests a possible role of mitonuclear interactions in maintaining mtDNA diversity.

Divergent mitochondrial lineages arose within a large, panmictic population of the Savannah sparrow (Passerculus sandwichensis)

Unusual patterns of mtDNA diversity can reveal interesting aspects of a species' biology. However, making such inferences requires discerning among the many alternative scenarios that could underlie any given mtDNA pattern. Next-generation sequencing methods provide large, multilocus data sets with increased power to resolve unusual mtDNA patterns. A mtDNA-based phylogeography of the Savannah sparrow (Passerculus sandwichensis) previously identified two sympatric, but divergent (~2%) clades within the nominate subspecies group and a third clade that consisted of birds sampled from northwest Mexico. We revisited the phylogeography of this species using a population genomic data set to resolve the processes leading to the evolution of sympatric and divergent mtDNA lineages. We identified two genetic clusters in the genomic data set corresponding to (a) the nominate subspecies group and (b) northwestern Mexico birds. Following divergence, the nominate clade maintained a large, stable population, indicating that divergent mitochondrial lineages arose within a panmictic population. Simulations based on parameter estimates from this model further confirmed that this demographic history could produce observed levels of mtDNA diversity. Patterns of divergent, sympatric mtDNA lineages are frequently interpreted as admixture of historically isolated lineages. Our analyses reject this interpretation for Savannah sparrows and underscore the need for genomic data sets to resolve the evolutionary mechanisms behind anomalous, locus-specific patterns.

Contrasting phylogeographic histories between the continent and islands of East Asia: Massive mitochondrial introgression and long-term isolation of hares (Lagomorpha: Lepus)

Hares of the genus Lepus are distributed worldwide, and introgressive hybridization is thought to be pervasive among species, leading to reticulate evolution and taxonomic confusion. Here, we performed phylogeographic analyses of the following species of hare across East Asia: L. timidus, L. mandshuricus, L. coreanus, and L. brachyurus collected from far-eastern Russia, South Korea, and Japan. Nucleotide sequences of one mitochondrial DNA and eight nuclear gene loci were examined, adding sequences of hares in China from databases. All nuclear DNA analyses supported the clear separation of three phylogroups: L. timidus, L. brachyurus, and the L. mandshuricus complex containing L. coreanus. On the other hand, massive mitochondrial introgression from two L. timidus lineages to the L. mandshuricus complex was suggested in continental East Asia. The northern population of the L. mandshuricus complex was mainly associated with introgression from the continental lineage of L. timidus, possibly since the last glacial period, whereas the southern population of the L. mandshuricus complex experienced introgression from another L. timidus lineage related to the Hokkaido population, possibly before the last glacial period. In contrast to continental hares, no evidence of introgression was found in L. brachyurus in the Japanese Archipelago, which showed the oldest divergence amongst East Asian hare lineages. Our findings suggest that glacial-interglacial climate changes in the circum-Japan Sea region promoted distribution shifts and introgressive hybridization among continental hare species, while the geographic structure of the region contributed to long-term isolation of hares on the islands, preventing inter-species gene flow.

Exploring the hybrid speciation continuum in birds

Hybridization is increasingly recognized as a creative evolutionary force contributing to adaptation and speciation. Homoploid hybrid speciation-the process in which hybridization results in a stable, fertile, and reproductively isolated hybrid lineage where there is no change in ploidy-has been documented in several taxa. Hybridization can directly contribute to reproductive isolation or reinforce it at a later stage. Alternatively, hybridization might not be related to the evolution of reproductive isolation. To account for these different scenarios, I propose to discriminate between two types of hybrid speciation: type I where reproductive isolation is a direct consequence of hybridization and type II where it is the by-product of other processes. I illustrate the applicability of this classification scheme with avian examples. To my knowledge, seven hybrid bird species have been proposed: Italian sparrow, Audubon's warbler, Genovesa mockingbird, Hawaiian duck, red-breasted goose, golden-crowned manakin, and a recent lineage of Darwin's finches on the island of Daphne Major ("Big Bird"). All studies provide convincing evidence for hybridization, but do not always confidently discriminate between scenarios of hybrid speciation and recurrent introgressive hybridization. The build-up of reproductive isolation between the hybrid species and their parental taxa is mainly driven by premating isolation mechanisms and comparable to classical speciation events. One hybrid species can be classified as type I ("Big Bird") while three species constitute type II hybrid species (Italian sparrow, Audubon's warbler, and golden-crowned manakin). The diversity in hybrid bird species across a range of divergence times also provides an excellent opportunity to study the evolution of hybrid genomes in terms of genome stabilization and adaptation.

Genetic determination of migration strategies in large soaring birds: evidence from hybrid eagles

The relative contributions of genetic and social factors in shaping the living world are a crucial question in ecology. The annual migration of birds to their wintering grounds and back provides significant knowledge in this field of research. Migratory movements are predominantly genetically determined in passerine birds, while in large soaring birds, it is presumed that social (cultural) factors play the largest role. In this study, we show that genetic factors in soaring birds are more important than previously assumed. We used global positioning system (GPS)-telemetry to compare the autumn journeys and wintering ranges of two closely related large raptorial bird species, the greater spotted eagle Clanga clanga and the lesser spotted eagle Clanga pomarina, and hybrids between them. The timing of migration in hybrids was similar to that of one parental species, but the wintering distributions and home range sizes were similar to those of the other. Tracking data were supported by habitat suitability modelling, based on GPS fixes and ring recoveries. These results suggest a strong genetic influence on migration strategy via a trait-dependent dominance effect, although we cannot rule out the contribution of social interactions.

Gene flow prevents mitonuclear co-adaptation: A comparative portrait of sympatric wild types and cybrids in the fish Chrosomus eos

Allospecific mtDNA can occasionally be beneficial for the fitness of populations. It is, however, difficult to assess the effect of mtDNA in natural conditions due to genetic and/or environmental interactions. In the fish Chrosomus eos, the transfer of C. neogaeus mitochondria occurs in a single generation and results in natural cybrids. For a few lakes in Quebec, C. eos can harbor either a C. eos mtDNA (wild types) or a C. neogaeus mtDNA (cybrids). Moreover, mtDNA of cybrids originated either from Mississippian or Atlantic glacial refuges. Such diversity provides a useful system for in situ assessment of allospecific mtDNA effects. We determined genetic, epigenetic and transcriptomic variation as well as mitochondrial enzymatic activity (complex IV) changes among wild types and cybrids either in sympatry or allopatry. Wild types and cybrids did not segregate spatially within a lake. Moreover, no significant genetic differentiation was detected among wild types and cybrids indicating sustained gene flow. Mitochondrial complex IV activity was higher for cybrids in both sympatry and allopatry while no difference was detected among cybrid haplotypes. Epigenetic and transcriptomic analyses revealed only subtle differences between sympatric wild types and cybrids compared to differences between sites. Altogether, these results indicate a limited influence of allospecific mtDNA in nuclear gene expression when controlling for genetic and environmental effects. The absence of a reproductive barrier between wild types and cybrids results in random association of either C. eos or C. neogaeus mtDNA with C. eos nDNA at each generation, and prevents mitonuclear co-adaptation in sympatry.

Extensive phenotypic diversification coexists with little genetic divergence and a lack of population structure in the White Wagtail subspecies complex (Motacilla alba)

Geographically clustered phenotypes often demonstrate consistent patterns in molecular markers, particularly mitochondrial DNA (mtDNA) traditionally used in phylogeographic studies. However, distinct evolutionary trajectories among traits and markers can lead to their discordance. First, geographic structure in phenotypic traits and nuclear molecular markers can be co-aligned but inconsistent with mtDNA (mito-nuclear discordance). Alternatively, phenotypic variation can have little to do with patterns in neither mtDNA nor nuclear markers. Disentangling between these distinct patterns can provide insight into the role of selection, demography and gene flow in population divergence. Here, we examined a previously reported case of strong inconsistency between geographic structure in mtDNA and plumage traits in a widespread polytypic bird species, the White Wagtail (Motacilla alba). We tested whether this pattern is due to mito-nuclear discordance or discrepancy between morphological evolution and both nuclear and mtDNA markers. We analysed population differentiation and structure across six out of nine commonly recognized subspecies using 17 microsatellite loci and a combination of microsatellites and plumage indices in a comprehensively sampled region of a contact between two subspecies. We did not find support for the mito-nuclear discordance hypothesis: nuclear markers indicated a subtle signal of genetic clustering only partially consistent with plumage groups, similar to previous findings that relied on mtDNA. We discuss evolutionary factors that could have shaped the intricate patterns of phenotypic diversification in the White wagtail and the role that repeated selection on plumage 'hotspots' and hybridization may have played.

Sex chromosomes and speciation in birds and other ZW systems

Theory and empirical patterns suggest a disproportionate role for sex chromosomes in evolution and speciation. Focusing on ZW sex determination (females ZW, males ZZ; the system in birds, many snakes, and lepidopterans), I review how evolutionary dynamics are expected to differ between the Z, W and the autosomes, discuss how these differences may lead to a greater role of the sex chromosomes in speciation and use data from birds to compare relative evolutionary rates of sex chromosomes and autosomes. Neutral mutations, partially or completely recessive beneficial mutations, and deleterious mutations under many conditions are expected to accumulate faster on the Z than on autosomes. Sexually antagonistic polymorphisms are expected to arise on the Z, raising the possibility of the spread of preference alleles. The faster accumulation of many types of mutations and the potential for complex evolutionary dynamics of sexually antagonistic traits and preferences contribute to a role for the Z chromosome in speciation. A quantitative comparison among a wide variety of bird species shows that the Z tends to have less within-population diversity and greater between-species differentiation than the autosomes, likely due to both adaptive evolution and a greater rate of fixation of deleterious alleles. The W chromosome also shows strong potential to be involved in speciation, in part because of its co-inheritance with the mitochondrial genome. While theory and empirical evidence suggest a disproportionate role for sex chromosomes in speciation, the importance of sex chromosomes is moderated by their small size compared to the whole genome.

The effects of allospecific mitochondrial genome on the fitness of northern redbelly dace (Chrosomus eos)

Instantaneous mitochondrial introgression events allow the disentangling of the effects of hybridization from those of allospecific mtDNA. Such process frequently occurred in the fish Chrosomus eos, resulting in cybrid individuals composed of a C. eos nuclear genome but with a C. neogaeus mtDNA. This provides a valuable model to address the fundamental question: How well do introgressed individuals perform in their native environment? We infer where de novo production of cybrids occurred to discriminate native environments from those colonized by cybrids in 25 sites from two regions (West-Qc and East-Qc) in Quebec (Canada). We then compared the relative abundance of wild types and cybrids as a measure integrating both fitness and de novo production of cybrids. According to mtDNA variation, 12 introgression events are required to explain the diversity of cybrids. Five cybrid lineages could not be associated with in situ introgression events. This includes one haplotype carried by 93% of the cybrids expected to have colonized West-Qc. These cybrids also displayed a nearly complete allopatric distribution with wild types. We still inferred de novo production of cybrids at seven sites, that accounted for 70% of the cybrids in East-Qc. Wild-type and cybrid individuals coexist in all East-Qc sites while cybrids were less abundant. Allopatry of cybrids restricted to the postglacial expansion suggests the existence of higher fitness for cybrids in specific conditions, allowing for the colonization of different environments and expanding the species' range. However, allospecific mtDNA does not provide a higher fitness to cybrids in their native environment compared to wild types, making the success of an introgressed lineage uncertain.

Hypoxia induces selective modifications to the acetylome in the brain of zebrafish (Danio rerio)

Reversible protein acetylation is an important regulatory mechanism for modulating protein function. The cellular protein acetylome is in large part dictated by the cellular redox balance, and in particular [NAD⁺]. While the relationship between hypoxia, redox balance, energy charge and resulting mitochondrial dysfunction has been examined in the context of hypoxia-linked pathologies, little is known about the direct effects of decreases in environmental oxygen on reversible lysine acetylation, and the resulting modifications to mitochondrial metabolism. To address this knowledge gap, we exposed zebrafish (Danio rerio) to 16 h of hypoxia (2.21 kPa) and quantified acetylation levels of 1220 proteins using whole-cell proteomics in samples of brain taken from normoxic and hypoxic zebrafish. In addition, we examined the effects of hypoxia on cytoplasmic and mitochondrial redox status, whole-cell energetics, the activity of the mitochondrial NAD⁺-dependent deacetylase SIRT3, and electron transport chain complex activities to determine if there is an association between hypoxia-induced metabolic disturbances, protein acetylation, and mitochondrial function. Our results (1) reveal several key changes in the acetylation status of proteins in the brain, primarily within the mitochondria; (2) show significant fluctuations in cytoplasmic and mitochondrial redox status within the brain during hypoxia exposure; and (3) provide evidence that lysine acetylation may be related to large changes in electron transport and ATP-synthase complex activities and adenylate status in zebrafish exposed to hypoxic stress. Together, these data provide new insights into the role of protein modifications in mitochondrial metabolism during hypoxia.

Deciphering the origin of mito-nuclear discordance in two sibling caddisfly species

An increasing number of phylogenetic studies have reported discordances among nuclear and mitochondrial markers. These discrepancies are highly relevant to widely used biodiversity assessment approaches, such as DNA barcoding, that rely almost exclusively on mitochondrial markers. Although the theoretical causes of mito-nuclear discordances are well understood, it is often extremely challenging to determine the principal underlying factor in a given study system. In this study, we uncovered significant mito-nuclear discordances in a pair of sibling caddisfly species. Application of genome sequencing, ddRAD and DNA barcoding revealed ongoing hybridization, as well as historical hybridization in Pleistocene refugia, leading us to identify introgression as the ultimate cause of the observed discordance pattern. Our novel genomic data, the discovery of a European-wide hybrid zone and the availability of established techniques for laboratory breeding make this species pair an ideal model system for studying species boundaries with ongoing gene flow.

Migratory orientation in a narrow avian hybrid zone

BACKGROUND

Zones of contact between closely related taxa with divergent migratory routes, termed migratory divides, have been suggested as areas where hybrid offspring may have intermediate and inferior migratory routes, resulting in low fitness of hybrids and thereby promoting speciation. In the Rocky Mountains of Canada there is a narrow hybrid zone between Audubon's and myrtle warblers that is likely maintained by selection against hybrids. Band recoveries and isotopic studies indicate that this hybrid zone broadly corresponds to the location of a possible migratory divide, with Audubon's warblers migrating south-southwest and myrtle warblers migrating southeast. We tested a key prediction of the migratory divide hypothesis: that genetic background would be predictive of migratory orientation among warblers in the center of the hybrid zone.

METHODS

We recorded fall migratory orientation of wild-caught migrating warblers in the center of the hybrid zone as measured by video-based monitoring of migratory restlessness in circular orientation chambers. We then tested whether there was a relationship between migratory orientation and genetic background, as measured using a set of species-specific diagnostic genetic markers.

RESULTS

We did not detect a significant association between orientation and genetic background. There was large variation among individuals in orientation direction. Mean orientation was towards the NE, surprising for birds on fall migration, but aligned with the mountain valley in which the study took place.

CONCLUSIONS

Only one other study has directly analyzed migratory orientation among naturally-produced hybrids in a migratory divide. While the other study showed an association between genetic background and orientation, we did not observe such an association in yellow-rumped warblers. We discuss possible reasons, including the possibility of a lack of a strong migratory divide in this hybrid zone and/or methodological limitations that may have prevented accurate measurements of long-distance migratory orientation.

Population structure of Purple Sandpipers (Calidris maritima) as revealed by mitochondrial DNA and microsatellites

The Purple Sandpiper (Calidris maritima) is a medium‐sized shorebird that breeds in the Arctic and winters along northern Atlantic coastlines. Migration routes and affiliations between breeding grounds and wintering grounds are incompletely understood. Some populations appear to be declining, and future management policies for this species will benefit from understanding their migration patterns. This study used two mitochondrial DNA markers and 10 microsatellite loci to analyze current population structure and historical demographic trends. Samples were obtained from breeding locations in Nunavut (Canada), Iceland, and Svalbard (Norway) and from wintering locations along the coast of Maine (USA), Nova Scotia, New Brunswick, and Newfoundland (Canada), and Scotland (UK). Mitochondrial haplotypes displayed low genetic diversity, and a shallow phylogeny indicating recent divergence. With the exception of the two Canadian breeding populations from Nunavut, there was significant genetic differentiation among samples from all breeding locations; however, none of the breeding populations was a monophyletic group. We also found differentiation between both Iceland and Svalbard breeding populations and North American wintering populations. This pattern of divergence is consistent with a previously proposed migratory pathway between Canadian breeding locations and wintering grounds in the United Kingdom, but argues against migration between breeding grounds in Iceland and Svalbard and wintering grounds in North America. Breeding birds from Svalbard also showed a genetic signature intermediate between Canadian breeders and Icelandic breeders. Our results extend current knowledge of Purple Sandpiper population genetic structure and present new information regarding migration routes to wintering grounds in North America.

Integrative Approaches for Studying Mitochondrial and Nuclear Genome Co-evolution in Oxidative Phosphorylation

In animals, interactions among gene products of mitochondrial and nuclear genomes (mitonuclear interactions) are of profound fitness, evolutionary, and ecological significance. Most fundamentally, the oxidative phosphorylation (OXPHOS) complexes responsible for cellular bioenergetics are formed by the direct interactions of 13 mitochondrial-encoded and ∼80 nuclear-encoded protein subunits in most animals. It is expected that organisms will develop genomic architecture that facilitates co-adaptation of these mitonuclear interactions and enhances biochemical efficiency of OXPHOS complexes. In this perspective, we present principles and approaches to understanding the co-evolution of these interactions, with a novel focus on how genomic architecture might facilitate it. We advocate that recent interdisciplinary advances assist in the consolidation of links between genotype and phenotype. For example, advances in genomics allow us to unravel signatures of selection in mitochondrial and nuclear OXPHOS genes at population-relevant scales, while newly published complete atomic-resolution structures of the OXPHOS machinery enable more robust predictions of how these genes interact epistatically and co-evolutionarily. We use three case studies to show how integrative approaches have improved the understanding of mitonuclear interactions in OXPHOS, namely those driving high-altitude adaptation in bar-headed geese, allopatric population divergence in Tigriopus californicus copepods, and the genome architecture of nuclear genes coding for mitochondrial functions in the eastern yellow robin.

Recovering the Genetic Identity of an Extinct-in-the-Wild Species: The Puzzling Case of the Alagoas Curassow

The conservation of many endangered taxa relies on hybrid identification, and when hybrids become morphologically indistinguishable from the parental species, the use of molecular markers can assign individual admixture levels. Here, we present the puzzling case of the extinct in the wild Alagoas Curassow (Pauxi mitu), whose captive population descends from only three individuals. Hybridization with the Razor-billed Curassow (P. tuberosa) began more than eight generations ago, and admixture uncertainty affects the whole population. We applied an analysis framework that combined morphological diagnostic traits, Bayesian clustering analyses using 14 microsatellite loci, and mtDNA haplotypes to assess the ancestry of all individuals that were alive from 2008 to 2012. Simulated data revealed that our microsatellites could accurately assign an individual a hybrid origin until the second backcross generation, which permitted us to identify a pure group among the older, but still reproductive animals. No wild species has ever survived such a severe bottleneck, followed by hybridization, and studying the recovery capability of the selected pure Alagoas Curassow group might provide valuable insights into biological conservation theory.