Search

Genetic analyses of host‐specific parasites can elucidate the evolutionary histories and biological features of their hosts. Here, we used population‐genomic analyses of ectoparasitic seal lice (Echinophthirius horridus) to shed light on the postglacial history of seals in the Arctic Ocean and the Baltic Sea region. One key question was the enigmatic origin of relict landlocked ringed seal populations in lakes Saimaa and Ladoga in northern Europe. We found that that lice of four postglacially diverged subspecies of the ringed seal (Pusa hispida) and Baltic gray seal (Halichoerus grypus), like their hosts, form genetically differentiated entities. Using coalescent‐based demographic inference, we show that the sequence of divergences of the louse populations is consistent with the geological history of lake formation. In addition, local effective population sizes of the lice are generally proportional to the census sizes of their respective seal host populations. Genome‐based reconstructions of long‐term effective population sizes revealed clear differences among louse populations associated with gray versus ringed seals, with apparent links to Pleistocene and Holocene climatic variation as well as to the isolation histories of ringed seal subspecies. Interestingly, our analyses also revealed ancient gene flow between the lice of Baltic gray and ringed seals, suggesting that the distributions of Baltic seals overlapped to a greater extent in the past than is the case today. Taken together, our results demonstrate how genomic information from specialized parasites with higher mutation and substitution rates than their hosts can potentially illuminate finer scale population genetic patterns than similar data from their hosts. [ABSTRACT FROM AUTHOR]

Hybridization between divergent lineages can result in losses of distinct evolutionary taxa. Alternatively, hybridization can lead to increased genetic variability that may fuel local adaptation and the generation of novel traits and/or taxa. Here, we examined single‐nucleotide polymorphisms generated using genotyping‐by‐sequencing in a population of Dolly Varden char (Pisces: Salmonidae) that is highly admixed within a contact zone between two subspecies (Salvelinus malma malma, Northern Dolly Varden [NDV] and S. m. lordi, Southern Dolly Varden [SDV]) in southwestern Alaska to assess the spatial distribution of hybrids and to test hypotheses on the origin of the admixed population. Ancestry analysis revealed that this admixed population is composed of advanced generation hybrids between NDV and SDV or advanced backcrosses to SDV; no F1 hybrids were detected. Coalescent‐based demographic modelling supported the origin of this population about 55,000 years ago by secondary contact between NDV and SDV with low levels of contemporary gene flow. Ancestry in NDV and SDV varies within the watershed and ancestry in NDV was positively associated with distance upstream from the sea, contingent on habitat‐type sampled, and negatively associated with the number of migrations that individual fish made to the sea. Our results suggest that divergence between subspecies over hundreds of thousands of years may not be associated with significant reproductive isolation, but that elevated diversity owing to hybridization may have contributed to adaptive divergence in habitat use and life history. [ABSTRACT FROM AUTHOR]

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non‐genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central–peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation‐with‐gene‐flow models generate the heterogeneous genome landscapes. These elements weight differently between central‐island and other central–peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low‐recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non‐uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana. [ABSTRACT FROM AUTHOR]

The harbour seal Phoca vitulina is a ubiquitous pinniped species found throughout coastal waters of the Northern Hemisphere. Harbour seal impacts on ecosystem dynamics may be significant due to their high abundance and food web position. Two subspecies exist in North America, P. v. richardii in the Pacific Ocean and P. v. vitulina in the Atlantic. Strong natal philopatry of harbour seals can result in fine‐scale genetic structure and isolation by distance. Management of harbour seals is expected to benefit from improved resolution of seal population structure and dynamics. Here, we use genotyping‐by‐sequencing to genotype 146 harbour seals from the eastern Pacific Ocean (i.e. British Columbia (BC), Oregon and California) and the western Atlantic Ocean (i.e. Québec, Newfoundland and Labrador). Using 12,742 identified variants, we confirm the recently identified elevated genetic diversity in the eastern Pacific relative to the western Atlantic and greatest differentiation between the subspecies. Further, we demonstrate that this is independent of reference genome bias or other potential technical artefacts. Coast‐specific analyses with 8933 and 3828 variants in Pacific and Atlantic subspecies, respectively, identify divergence between BC and Oregon–California, and between Québec and Newfoundland–Labrador. Unexpected PCA outlier clusters were observed in two populations due to cryptic relatedness of individuals; subsequently, closely related samples were removed. Admixture analysis indicates an isolation‐by‐distance signature where Oregon seals contained some of the BC signature, whereas California did not. Additional sampling is needed in the central and north coast of BC to determine whether a discrete separation of populations exists within the region. [ABSTRACT FROM AUTHOR]

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present‐day individuals representing the current native range of both muskox subspecies, the white‐faced and the barren‐ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000‐year‐old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post‐glacial re‐colonizations. In particular, the white‐faced muskox has the lowest genome‐wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long‐term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds. [ABSTRACT FROM AUTHOR]

Phylogeographic studies uncover hidden pathways of divergence and inform conservation. Brown bears (Ursus arctos) have one of the broadest distributions of all land mammals, ranging from Eurasia to North America, and are an important model for evolutionary studies. Although several whole genomes were available for individuals from North America, Europe and Asia, limited whole‐genome data were available from Central Asia, including the highly imperilled brown bears in the Gobi Desert. To fill this knowledge gap, we sequenced whole genomes from nine Asian brown bears from the Gobi Desert of Mongolia, Northern Mongolia and the Himalayas of Pakistan. We combined these data with published brown bear sequences from Europe, Asia and North America, as well as other bear species. Our goals were to determine the evolutionary relationships among brown bear populations worldwide, their genetic diversity and their historical demography. Our analyses revealed five major lineages of brown bears based on a filtered set of 684,081 single nucleotide polymorphisms. We found distinct evolutionary lineages of brown bears in the Gobi, Himalayas, northern Mongolia, Europe and North America. The lowest level of genetic diversity and the highest level of inbreeding were found in Pakistan, the Gobi Desert and Central Italy. Furthermore, the effective population size (Ne) for all brown bears decreased over the last 70,000 years. Our results confirm the genetic distinctiveness and ancient lineage of brown bear subspecies in the Gobi Desert of Mongolia and the Himalayas of Pakistan and highlight their importance for conservation. [ABSTRACT FROM AUTHOR]

Early lineage diversification is central to understand what mutational events drive species divergence. Particularly, gene misregulation in interspecific hybrids can inform about what genes and pathways underlie hybrid dysfunction. In Drosophila hybrids, how regulatory evolution impacts different reproductive tissues remains understudied. Here, we generate a new genome assembly and annotation in Drosophila willistoni and analyse the patterns of transcriptome divergence between two allopatrically evolved D. willistoni subspecies, their male sterile and female fertile hybrid progeny across testis, male accessory gland, and ovary. Patterns of transcriptome divergence and modes of regulatory evolution were tissue‐specific. Despite no indication for cell‐type differences in hybrid testis, this tissue exhibited the largest magnitude of expression differentiation between subspecies and between parentals and hybrids. No evidence for anomalous dosage compensation in hybrid male tissues was detected nor was a differential role for the neo‐ and the ancestral arms of the D. willistoni X chromosome. Compared to the autosomes, the X chromosome appeared enriched for transgressively expressed genes in testis despite being the least differentiated in expression between subspecies. Evidence for fine genome clustering of transgressively expressed genes suggests a role of chromatin structure on hybrid gene misregulation. Lastly, transgressively expressed genes in the testis of the sterile male progeny were enriched for GO terms not typically associated with sperm function, instead hinting at anomalous development of the reproductive tissue. Our thorough tissue‐level portrait of transcriptome differentiation between recently diverged D. willistoni subspecies and their hybrids provides a more nuanced view of early regulatory changes during speciation. [ABSTRACT FROM AUTHOR]

After establishing secondary contact, recently diverged populations may remain reproductively isolated or may hybridize to a varying extent depending on factors such as hybrid fitness and the strength of assortative mating. Here, we used genomic and phenotypic data from three independent contact zones between subspecies of the variable seedeater (Sporophila corvina) to examine how coloration and genetic divergence shape patterns of hybridization. We found that differences in plumage coloration are probably maintained by divergent selection across contact zones; however, the degree of plumage differentiation does not match overall patterns of hybridization. Across two parallel contact zones between populations with divergent phenotypes (entirely black vs. pied plumage), populations hybridized extensively across one contact zone but not the other, suggesting that plumage divergence is not sufficient to maintain reproductive isolation. Where subspecies hybridized, hybrid zones were wide and formed by later‐generation hybrids, suggesting frequent reproduction and high survivorship for hybrid individuals. Moreover, contemporary gene flow has played an important role in shaping patterns of genetic structure between populations. Replicated contact zones between hybridizing taxa offer a unique opportunity to explore how different factors interact to shape patterns of hybridization. Overall, our results demonstrate that divergence in plumage coloration is important in reducing gene flow but insufficient in maintaining reproductive isolation in this clade, and that other factors such as divergence in song and time since secondary contact may also play an important role in driving patterns of reduced hybridization and gene flow. [ABSTRACT FROM AUTHOR]

Oceanic islands are characterized by conditions that favour diversification into endemic lineages that can be very different from their mainland counterparts. This can be the result of fast phenotypic divergence due to drift or the result of slower adaptation to local conditions. This uniqueness can obscure their evolutionary history. Here we used morphological, stable isotope, genetic and genomic data to characterize common quails (Coturnix coturnix) in the Azores archipelago and assess the divergence from neighbouring common quail populations. Historical documents suggested that these quails could have a recent origin associated with the arrival of humans in the last centuries. Our results show that Azorean quails constitute a well‐differentiated lineage with small size and dark throat pigmentation that has lost the migratory ability and that diverged from mainland quail lineages more than 0.8 mya, contrary to the notion of a recent human‐mediated arrival. Even though some Azorean quails carry an inversion that affects 115 Mbp of chromosome 1 and that has been associated with the loss of the migratory behaviour in other common quail populations, half of the analysed individuals do not have that inversion and still do not migrate. The long coexistence and evolution in isolation in the Azores of two chromosomal variants (with and without the inversion) is best explained by balancing selection. Thus, a unique and long evolutionary history led to the island endemic that we know today, C. c. conturbans. [ABSTRACT FROM AUTHOR]

In the context of the current extinction crisis, identifying new conservation units is pivotal to the development of sound conservation measures, especially in highly threatened taxa such as felids. Corsican wildcats are known by Corsican people since a very long time but have been little studied. Meaningful information about their phylogenetic position is lacking. We used ddRADseq to genotype phenotypically homogenous Corsican wildcats at 3671 genome‐wide SNPs and reported for the first time their genetic identity. We compared this genomic information to domestic cats Felis silvestris catus from Corsica and mainland France, European wildcats F. s. silvestris and Sardinian wildcats F. s. lybica. Our premise was that if the Corsican wildcat, as a phenotypic entity, also represents a genetic entity, it deserves conservation measures and to be recognized as a conservation unit. Corsican wildcats appeared highly genetically differentiated from European wildcats and genetically closer to Sardinian wildcats than to domestic cats. Domestic cats from Corsica and mainland France were closer to each other and Sardinian wildcats were intermediate between Corsican wildcats and domestic cats. This suggested that Corsican wildcats do not belong to the F. s. silvestris or catus lineages. The inclusion of more high‐quality Sardinian samples and Near‐Eastern mainland F. s. lybica would constitute the next step toward assessing the status of Corsican wildcat as a subspecies and/or evolutionarily significant unit and tracing back wildcat introduction history of in Corsica. [ABSTRACT FROM AUTHOR]

Present‐day ecology and population structure are the legacies of past climate and habitat perturbations, and this is particularly true for species that are widely distributed at high latitudes. The red knot, Calidris canutus, is an arctic‐breeding, long‐distance migratory shorebird with six recognized subspecies defined by differences in morphology, migration behavior, and annual cycle phenology, in a global distribution thought to have arisen just since the last glacial maximum (LGM). We used nextRAD sequencing of 10,881 single‐nucleotide polymorphisms (SNPs) to assess the neutral genetic structure and phylogeographic history of 172 red knots representing all known global breeding populations. Using population genetics approaches, including model‐based scenario‐testing in an approximate Bayesian computation (ABC) framework, we infer that red knots derive from two main lineages that diverged ca. 34,000 years ago, and thus most probably persisted at the LGM in both Palearctic and Nearctic refugia, followed by at least two instances of secondary contact and admixture. Within two Beringian subspecies (C. c. roselaari and rogersi), we detected previously unknown genetic structure among sub‐populations sharing a migratory flyway, reflecting additional complexity in the phylogeographic history of the region. Conversely, we found very weak genetic differentiation between two Nearctic populations (rufa and islandica) with clearly divergent migratory phenotypes and little or no apparent contact throughout the annual cycle. Together, these results suggest that relative gene flow among migratory populations reflects a complex interplay of historical, geographical, and ecological factors. [ABSTRACT FROM AUTHOR]

Hybrids between species are often sterile or inviable. Hybrid unfitness usually evolves first in the heterogametic sex—a pattern known as Haldane's rule. The genetics of Haldane's rule have been extensively studied in species where the male is the heterogametic (XX/XY) sex, but its basis in taxa where the female is heterogametic (ZW/ZZ), such as Lepidoptera and birds, is largely unknown. Here, we analyse a new case of female hybrid sterility between geographic subspecies of Heliconius pardalinus. The two subspecies mate freely in captivity, but female F1 hybrids in both directions of cross are sterile. Sterility is due to arrested development of oocytes after they become differentiated from nurse cells, but before yolk deposition. We backcrossed fertile male F1 hybrids to parental females and mapped quantitative trait loci (QTLs) for female sterility. We also identified genes differentially expressed in the ovary as a function of oocyte development. The Z chromosome has a major effect, similar to the 'large X effect' in Drosophila, with strong epistatic interactions between loci at either end of the Z chromosome, and between the Z chromosome and autosomal loci on chromosomes 8 and 20. By intersecting the list of genes within these QTLs with those differentially expressed in sterile and fertile hybrids, we identified three candidate genes with relevant phenotypes. This study is the first to characterize hybrid sterility using genome mapping in the Lepidoptera and shows that it is produced by multiple complex epistatic interactions often involving the sex chromosome, as predicted by the dominance theory of Haldane's rule. [ABSTRACT FROM AUTHOR]

In 1859, Charles Darwin proposed that species are not fundamentally different from subspecies or the varieties from which they evolve. A century later, Dobzhansky (1958) suggested that many such lineages are ephemeral and are likely to reverse differentiation through introgression (Figure 1a); only a few evolve complete reproductive isolation and persist in sympatry. In this issue of Molecular Ecology, Bouzid et al. (2021) showed how new analytical methods, when applied to genomic data, allow us to more precisely determine whether or not species formation follows the paths outlined by Darwin and Dobzhansky (Figure 1b). The authors studied the diversification of the lizard Sceloporus occidentalis, finding a continuum of genetic interactions between the preservation of genetic identity to genetic merger, analogous to what is exemplified by ring species. In doing so, they teach us two tales about species formation: that lineages are fractal byproducts of evolutionary processes such as genetic drift and selection, and that lineages are often ephemeral and do not always progress into fully reproductively isolated taxa. Studying ephemeral lineages like those in S. occidentalis allows us to capture divergence at its earliest stages, and potentially to determine the factors that allow lineages to remain distinct despite pervasive gene flow. These lineages thus serve as a natural laboratory to address long standing hypotheses about species formation. [ABSTRACT FROM AUTHOR]

The origins of geographic races in wide‐ranging species are poorly understood. In Texas, the texanus subspecies of Helianthusannuus has long been thought to have acquired its defining phenotypic traits via introgression from a local congener, H. debilis, but previous tests of this hypothesis were inconclusive. Here, we explore the origins of H. a. texanus using whole genome sequencing data from across the entire range of H. annuus and possible donor species, as well as phenotypic data from a common garden study. We found that although it is morphologically convergent with H. debilis, H. a. texanus has conflicting signals of introgression. Genome wide tests (Patterson's D and TreeMix) only found evidence of introgression from H. argophyllus (sister species to H. annuus and also sympatric), but not H. debilis, with the exception of one individual of 109 analysed. We further scanned the genome for localized signals of introgression using PCAdmix and found minimal but nonzero introgression from H. debilis and significant introgression from H. argophyllus in some populations. Given the paucity of introgression from H. debilis, we argue that the morphological convergence observed in Texas is probably from standing genetic variation. We also found that genomic differentiation in H. a. texanus is mostly driven by large segregating inversions, several of which have signatures of natural selection based on haplotype frequencies. [ABSTRACT FROM AUTHOR]

When populations colonize new areas, both strong selection and strong drift can be experienced due to novel environments and small founding populations, respectively. Empirical studies have predominantly focused on the phenotype when assessing the role of selection, and limited neutral‐loci when assessing founder‐induced loss of diversity. Consequently, the extent to which processes interact to influence evolutionary trajectories is difficult to assess. Genomic‐level approaches provide the opportunity to simultaneously consider these processes. Here, we examine the roles of selection and drift in shaping genomic diversity and divergence in historically documented sequential island colonizations by the silvereye (Zosterops lateralis). We provide the first empirical demonstration of the rapid appearance of highly diverged genomic regions following population founding, the position of which are highly idiosyncratic. As these regions rarely contained loci putatively under selection, it is most likely that these differences arise via the stochastic nature of the founding process. However, selection is required to explain rapid evolution of larger body size in insular silvereyes. Reconciling our genomic data with these phenotypic patterns suggests there may be many genomic routes to the island phenotype, which vary across populations. Finally, we show that accelerated divergence associated with multiple founding steps is the product of genome‐wide rather than localized differences, and that diversity erodes due to loss of rare alleles. However, even multiple founder events do not result in divergence and diversity levels seen in evolutionary older subspecies, and therefore do not provide a shortcut to speciation as proposed by founder‐effect speciation models. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Grant's gazelles have recently been proposed to be a species complex comprising three highly divergent mtDNA lineages (Nanger granti, N. notata and N. petersii). The three lineages have nonoverlapping distributions in East Africa, but without any obvious geographical divisions, making them an interesting model for studying the early‐stage evolutionary dynamics of allopatric speciation in detail. Here, we use genomic data obtained by restriction site‐associated (RAD) sequencing of 106 gazelle individuals to shed light on the evolutionary processes underlying Grant's gazelle divergence, to characterize their genetic structure and to assess the presence of gene flow between the main lineages in the species complex. We date the species divergence to 134,000 years ago, which is recent in evolutionary terms. We find population subdivision within N. granti, which coincides with the previously suggested two subspecies, N. g. granti and N. g. robertsii. Moreover, these two lineages seem to have hybridized in Masai Mara. Perhaps more surprisingly given their extreme genetic differentiation, N. granti and N. petersii also show signs of prolonged admixture in Mkomazi, which we identified as a hybrid population most likely founded by allopatric lineages coming into secondary contact. Despite the admixed composition of this population, elevated X chromosomal differentiation suggests that selection may be shaping the outcome of hybridization in this population. Our results therefore provide detailed insights into the processes of allopatric speciation and secondary contact in a recently radiated species complex. [ABSTRACT FROM AUTHOR]

Morphological traits have served generations of biologists as a taxonomic indicator, and have been the main basis for defining and classifying species diversity for centuries. A quantitative integration of behavioural characters, such as vocalizations, in studies on biotic differentiation has arisen more recently, and the relative importance of these different traits in the diversification process remains poorly understood. To provide a framework within which to interpret the evolutionary interplay between morphological and behavioural traits, we generated a draft genome of a cryptic Southeast Asian songbird, the limestone wren‐babbler Napothera crispifrons. We resequenced whole genomes of multiple individuals of all three traditional subspecies and of a distinct leucistic population. We demonstrate strong genomic and mitochondrial divergence among all three taxa, pointing to the existence of three species‐level lineages. Despite its great phenotypic distinctness, the leucistic population was characterized by shallow genomic differentiation from its neighbour, with only a few localized regions emerging as highly diverged. Quantitative bioacoustic analysis across multiple traits revealed deep differences especially between the two taxa characterized by limited plumage differentiation. Our study demonstrates that differentiation in these furtive songbirds has resulted in a complex mosaic of colour‐based and bioacoustic differences among populations. Extreme colour differences can be anchored in few genomic loci and may therefore arise and subside rapidly. [ABSTRACT FROM AUTHOR]

Phenotypic divergence among natural populations can be explained by natural selection or by neutral processes such as drift. Many examples in the literature compare putatively neutral (FST) and quantitative genetic (QST) differentiation in multiple populations to assess their evolutionary signature and identify candidate traits involved with local adaptation. Investigating these signatures in closely related or recently diversified species has the potential to shed light on the divergence processes acting at the interspecific level. Here, we conducted this comparison in two subspecies of snapdragon plants (eight populations of Antirrhinum majus pseudomajus and five populations of A. m. striatum) in a common garden experiment. We also tested whether altitude was involved with population phenotypic divergence. Our results identified candidate phenological and morphological traits involved with local adaptation. Most of these traits were identified in one subspecies but not the other. Phenotypic divergence increased with altitude for a few biomass‐related traits, but only in A. m. striatum. These traits therefore potentially reflect A. m. striatum adaptation to altitude. Our findings imply that adaptive processes potentially differ at the scale of A. majus subspecies. [ABSTRACT FROM AUTHOR]

Parallel evolution can occur through selection on novel mutations, standing genetic variation or adaptive introgression. Uncovering parallelism and introgressed populations can complicate management of threatened species as parallelism may have influenced conservation unit designations and admixed populations are not generally considered under legislations. We examined high coverage whole‐genome sequences of 30 caribou (Rangifer tarandus) from across North America and Greenland, representing divergent intraspecific lineages, to investigate parallelism and levels of introgression contributing to the formation of ecotypes. Caribou are split into four subspecies and 11 extant conservation units, known as designatable units (DUs), in Canada. Using genomes from all four subspecies and six DUs, we undertake demographic reconstruction and confirm two previously inferred instances of parallel evolution in the woodland subspecies and uncover an additional instance of parallelism of the eastern migratory ecotype. Detailed investigations reveal introgression in the woodland subspecies, with introgressed regions found spread throughout the genomes encompassing both neutral and functional sites. Our investigations using whole genomes highlight the difficulties in unequivocally demonstrating parallelism through adaptive introgression in nonmodel species with complex demographic histories, with standing variation and introgression both potentially involved. Additionally, the impact of parallelism and introgression on conservation policy for management units needs to be considered in general, and the caribou designations will need amending in light of our results. Uncovering and decoupling parallelism and differential patterns of introgression will become prevalent with the availability of comprehensive genomic data from nonmodel species, and we highlight the need to incorporate this into conservation unit designations. [ABSTRACT FROM AUTHOR]

As populations diverge many processes can shape genomic patterns of differentiation. Regions of high differentiation can arise due to divergent selection acting on selected loci, genetic hitchhiking of nearby loci, or through repeated selection against deleterious alleles (linked background selection); this divergence may then be further elevated in regions of reduced recombination. Atlantic salmon (Salmo salar) from Europe and North America diverged >600,000 years ago and despite some evidence of secondary contact, the majority of genetic data indicate substantial divergence between lineages. This deep divergence with potential gene flow provides an opportunity to investigate the role of different mechanisms that shape the genomic landscape during early speciation. Here, using 184,295 single nucleotide polymorphisms (SNPs) and 80 populations, we investigate the genomic landscape of differentiation across the Atlantic Ocean with a focus on highly differentiated regions and the processes shaping them. We found evidence of high (mean FST = 0.26) and heterogeneous genomic differentiation between continents. Genomic regions associated with high trans‐Atlantic differentiation ranged in size from single loci (SNPs) within important genes to large regions (1–3 Mbp) on four chromosomes (Ssa06, Ssa13, Ssa16 and Ssa19). These regions showed signatures consistent with selection, including high linkage disequilibrium, despite no significant reduction in recombination. Genes and functional enrichment of processes associated with differentiated regions may highlight continental differences in ocean navigation and parasite resistance. Our results provide insight into potential mechanisms underlying differences between continents, and evidence of near‐fixed and potentially adaptive trans‐Atlantic differences concurrent with a background of high genome‐wide differentiation supports subspecies designation in Atlantic salmon. [ABSTRACT FROM AUTHOR]

The Kinship Theory of Genomic Imprinting (KTGI) posits that, in species where females mate with multiple males, there is selection for a male to enhance the reproductive success of his offspring at the expense of other males and his mating partner. Reciprocal crosses between honey bee subspecies show parent‐of‐origin effects for reproductive traits, suggesting that males modify the expression of genes related to female function in their female offspring. This effect is likely to be greater in the Cape honey bee (Apis mellifera capensis), because a male's daughters have the unique ability to produce female offspring that can develop into reproductive workers or the next queen without mating. We generated reciprocal crosses between Capensis and another subspecies and used RNA‐seq to identify transcripts that are over‐ or underexpressed in the embryos, depending on the parental origin of the gene. As predicted, 21 genes showed expression bias towards the Capensis father's allele in colonies with a Capensis father, with no such bias in the reciprocal cross. A further six genes showed a consistent bias towards expression of the father's allele across all eight colonies examined, regardless of the direction of the cross. Consistent with predictions of the KTGI, six of the 21 genes are associated with female reproduction. No gene consistently showed overexpression of the maternal allele. [ABSTRACT FROM AUTHOR]

Disentangling the contribution of long‐term evolutionary processes and recent anthropogenic impacts to current genetic patterns of wildlife species is key to assessing genetic risks and designing conservation strategies. Here, we used 80 whole nuclear genomes and 96 mitogenomes from populations of the Eurasian lynx covering a range of conservation statuses, climatic zones and subspecies across Eurasia to infer the demographic history, reconstruct genetic patterns, and discuss the influence of long‐term isolation and/or more recent human‐driven changes. Our results show that Eurasian lynx populations shared a common history until 100,000 years ago, when Asian and European populations started to diverge and both entered a period of continuous and widespread decline, with western populations, except Kirov, maintaining lower effective sizes than eastern populations. Population declines and increased isolation in more recent times probably drove the genetic differentiation between geographically and ecologically close westernmost European populations. By contrast, and despite the wide range of habitats covered, populations are quite homogeneous genetically across the Asian range, showing a pattern of isolation by distance and providing little genetic support for the several proposed subspecies. Mitogenomic and nuclear divergences and population declines starting during the Late Pleistocene can be mostly attributed to climatic fluctuations and early human influence, but the widespread and sustained decline since the Holocene is more probably the consequence of anthropogenic impacts which intensified in recent centuries, especially in western Europe. Genetic erosion in isolated European populations and lack of evidence for long‐term isolation argue for the restoration of lost population connectivity. [ABSTRACT FROM AUTHOR]

Variation in the metabolic costs associated with organismal maintenance may play a key role in determining fitness, and thus these differences among individuals are likely to be subject to natural selection. Although the evolvability of maintenance metabolism depends on its underlying genetic architecture, relatively little is known about the nature of genetic variation that underlies this trait. To address this, we measured variation in routine metabolic rate (ṀO2routine), an index of maintenance metabolism, within and among three populations of Atlantic killifish, Fundulus heteroclitus, including a population from a region of genetic admixture between two subspecies. Polygenic association tests among individuals from the admixed population identified 54 single nucleotide polymorphisms (SNPs) that were associated with ṀO2routine, and these SNPs accounted for 43% of interindividual variation in this trait. However, genetic associations with ṀO2routine involved different SNPs if females and males were analysed separately, and there was a sex‐dependent effect of mitochondrial genotype on variation in routine metabolism. These results imply that there are sex‐specific genetic mechanisms, and potential mitonuclear interactions, that underlie variation in ṀO2routine. Additionally, there was evidence for epistatic interactions between 17% of the possible pairs of trait‐associated SNPs, suggesting that epistatic effects on ṀO2routine are common. These data demonstrate not only that phenotypic variation in this ecologically important trait has a polygenic basis with considerable epistasis among loci, but also that these underlying genetic mechanisms, and particularly the role of mitochondrial genotype, may be sex‐specific. [ABSTRACT FROM AUTHOR]

Genomic phylogeography plays an important role in describing evolutionary processes and their geographic, ecological, or cultural drivers. These drivers are often poorly understood in marine environments, which have fewer obvious barriers to mixing than terrestrial environments. Taxonomic uncertainty of some taxa (e.g., cetaceans), due to the difficulty in obtaining morphological data, can hamper our understanding of these processes. One such taxon, the short‐finned pilot whale, is recognized as a single global species but includes at least two distinct morphological forms described from stranding and drive hunting in Japan, the "Naisa" and "Shiho" forms. Using samples (n = 735) collected throughout their global range, we examine phylogeographic patterns of divergence by comparing mitogenomes and nuclear SNP loci. Our results suggest three types within the species: an Atlantic Ocean type, a western/central Pacific and Indian Ocean (Naisa) type, and an eastern Pacific Ocean and northern Japan (Shiho) type. mtDNA control region differentiation indicates these three types form two subspecies, separated by the East Pacific Barrier: Shiho short‐finned pilot whale, in the eastern Pacific Ocean and northern Japan, and Naisa short‐finned pilot whale, throughout the remainder of the species' distribution. Our data further indicate two diverging populations within the Naisa subspecies, in the Atlantic Ocean and western/central Pacific and Indian Oceans, separated by the Benguela Barrier off South Africa. This study reveals a process of divergence and speciation within a globally‐distributed, mobile marine predator, and indicates the importance of the East Pacific Barrier to this evolutionary process. [ABSTRACT FROM AUTHOR]

Patterns of genomic divergence between hybridizing taxa can be heterogeneous along the genome. Both differential introgression and local adaptation may contribute to this pattern. Here, we analysed two teosinte subspecies, Zea mays ssp. parviglumis and ssp. mexicana, to test whether their divergence has occurred in the face of gene flow and to infer which environmental variables have been important drivers of their ecological differentiation. We generated 9,780 DArTseqTM SNPs for 47 populations, and used an additional data set containing 33,454 MaizeSNP50 SNPs for 49 populations. With these data, we inferred features of demographic history and performed genome wide scans to determine the number of outlier SNPs associated with climate and soil variables. The two data sets indicate that divergence has occurred or been maintained despite continuous gene flow and/or secondary contact. Most of the significant SNP associations were to temperature and to phosphorus concentration in the soil. A large proportion of these candidate SNPs were located in regions of high differentiation that had been identified previously as putative inversions. We therefore propose that genomic differentiation in teosintes has occurred by a process of adaptive divergence, with putative inversions contributing to reduced gene flow between locally adapted populations. [ABSTRACT FROM AUTHOR]

Recently evolved species typically share genetic variation across their genomes due to incomplete lineage sorting and/or ongoing gene flow. Given only subtle allele frequency differences at most loci and the expectation that divergent selection may affect only a tiny fraction of the genome, distinguishing closely related species based on multi‐locus data requires substantial genomic coverage. In this study, we used ddRAD‐seq to sample the genomes of five recently diverged, New World "mallards" (Anas spp.), a group of dabbling duck species characterized by diagnosable phenotypic differences but minimal genetic differentiation. With increased genomic sampling, we aimed to characterize population structure within this group and identify genomic regions that may have experienced divergent selection during speciation. We analyzed 3,017 autosomal ddRAD‐seq loci and 177 loci from the Z‐chromosome. In contrast to previous studies, the ddRAD‐seq data were sufficient to assign individuals to their respective species or subspecies and to generate estimates of gene flow in a phylogenetic framework. We find limited evidence of contemporary gene flow between the dichromatic mallard and several monochromatic taxa, but find evidence for historical gene flow between some monochromatic species pairs. We conclude that the overall genetic similarity of these taxa likely reflects retained ancestral polymorphism rather than recent and extensive gene flow. Thus, despite recurring cases of hybridization in this group, our results challenge the current dogma predicting the genetic extinction of the New World monochromatic dabbling ducks via introgressive hybridization with mallards. Moreover, ddRAD‐seq data were sufficient to identify previously unknown outlier regions across the Z‐chromosome and several autosomal chromosomes that may have been involved in the diversification of species in this recent radiation. [ABSTRACT FROM AUTHOR]

Divergence in sexual signals may drive reproductive isolation between lineages, but behavioural barriers can weaken in contact zones. Here, we investigate the role of song as a behavioural and genetic barrier in a contact zone between two subspecies of white-crowned sparrows ( Zonotrichia leucophrys). We employed a reduced genomic data set to assess population structure and infer the history underlying divergence, gene flow and hybridization. We also measured divergence in song and tested behavioural responses to song using playback experiments within and outside the contact zone. We found that the subspecies form distinct genetic clusters, and demographic inference supported a model of secondary contact. Song phenotype, particularly length of the first note (a whistle), was a significant predictor of genetic subspecies identity and genetic distance along the hybrid zone, suggesting a close link between song and genetic divergence in this system. Individuals from both parental and admixed localities responded significantly more strongly to their own song than to the other subspecies song, supporting song as a behavioural barrier. Putative parental and admixed individuals were not significantly different in their strength of discrimination between own and other songs; however, individuals from admixed localities tended to discriminate less strongly, and this difference in discrimination strength was explained by song dissimilarity as well as genetic distance. Therefore, we find that song acts as a reproductive isolating mechanism that is potentially weakening in a contact zone between the subspecies. Our findings also support the hypothesis that intraspecific song variation can reduce gene flow between populations. [ABSTRACT FROM AUTHOR]

Understanding the processes that drive divergence within and among species is a long-standing goal in evolutionary biology. Traditional approaches to assessing differentiation rely on phenotypes to identify intra- and interspecific variation, but many species express subtle morphological gradients in which boundaries among forms are unclear. This intraspecific variation may be driven by differential adaptation to local conditions and may thereby reflect the evolutionary potential within a species. Here, we combine genetic and morphological data to evaluate intraspecific variation within the Nelson's ( Ammodramus nelsoni) and salt marsh ( Ammodramus caudacutus) sparrow complex, a group with populations that span considerable geographic distributions and a habitat gradient. We evaluated genetic structure among and within five putative subspecies of A. nelsoni and A. caudacutus using a reduced-representation sequencing approach to generate a panel of 1929 SNPs among 69 individuals. Although we detected morphological differences among some groups, individuals sorted along a continuous phenotypic gradient. In contrast, the genetic data identified three distinct clusters corresponding to populations that inhabit coastal salt marsh, interior freshwater marsh and coastal brackish-water marsh habitats. These patterns support the current species-level recognition but do not match the subspecies-level taxonomy within each species-a finding which may have important conservation implications. We identified loci exhibiting patterns of elevated divergence among and within these species, indicating a role for local selective pressures in driving patterns of differentiation across the complex. We conclude that this evidence for adaptive variation among subspecies warrants the consideration of evolutionary potential and genetic novelty when identifying conservation units for this group. [ABSTRACT FROM AUTHOR]

Genetic diversity is positively linked to the viability and evolutionary potential of species but is often compromised in threatened taxa. Genetic rescue by gene flow from a more diverse or differentiated source population of the same species can be an effective strategy for alleviating inbreeding depression and boosting evolutionary potential. The helmeted honeyeater Lichenostomus melanops cassidix is a critically endangered subspecies of the common yellow-tufted honeyeater. Cassidix has declined to a single wild population of ~130 birds, despite being subject to intensive population management over recent decades. We assessed changes in microsatellite diversity in cassidix over the last four decades and used population viability analysis to explore whether genetic rescue through hybridization with the neighbouring Lichenostomus melanops gippslandicus subspecies constitutes a viable conservation strategy. The contemporary cassidix population is characterized by low genetic diversity and effective population size ( Ne < 50), suggesting it is vulnerable to inbreeding depression and will have limited capacity to evolve to changing environments. We find that gene flow from gippslandicus to cassidix has declined substantially relative to pre-1990 levels and argue that natural levels of gene flow between the two subspecies should be restored. Allowing gene flow (~4 migrants per generation) from gippslandicus into cassidix (i.e. genetic rescue), in combination with continued annual release of captive-bred cassidix (i.e. demographic rescue) , should lead to positive demographic and genetic outcomes. Although we consider the risk of outbreeding depression to be low, we recommend that genetic rescue be managed within the context of the captive breeding programme, with monitoring of outcomes. [ABSTRACT FROM AUTHOR]

The identification of the genes involved in morphological variation in nature is still a major challenge. Here, we explore a new approach: we combine 178 samples from a natural hybrid zone between two subspecies of the house mouse ( Mus musculus domesticus and Mus musculus musculus), and high coverage of the genome (~ 145K SNPs) to identify loci underlying craniofacial shape variation. Due to the long history of recombination in the hybrid zone, high mapping resolution is anticipated. The combination of genomes from subspecies allows the mapping of both, variation within subspecies and inter-subspecific differences, thereby increasing the overall amount of causal genetic variation that can be detected. Skull and mandible shape were measured using 3D landmarks and geometric morphometrics. Using principal component axes as phenotypes, and a linear mixed model accounting for genetic relatedness in the mapping populations, we identified nine genomic regions associated with skull shape and 10 with mandible shape. High mapping resolution (median size of significant regions = 148 kb) enabled identification of single or few candidate genes in most cases. Some of the genes act as regulators or modifiers of signalling pathways relevant for morphological development and bone formation, including several with known craniofacial phenotypes in mice and humans. The significant associations combined explain 13% and 7% of the skull and mandible shape variation, respectively. In addition, a positive correlation was found between chromosomal length and proportion of variation explained. Our results suggest a complex genetic architecture for shape traits and support a polygenic model. [ABSTRACT FROM AUTHOR]

Recent molecular studies have shown that highly mobile species with continuous distributions can exhibit fine-scale population structure. In this context, we assessed genetic structure within a marine species with high dispersal potential, the Atlantic spotted dolphin ( Stenella frontalis). Using 19 microsatellite loci and mitochondrial control region sequences, population structure was investigated in the western North Atlantic, the Gulf of Mexico and the Azores Islands. Analyses of the microsatellite data identified four distinct genetic clusters, which were supported by the control region sequences. The highest level of divergence was seen between two clusters corresponding to previously described morphotypes that inhabit oceanic and shelf waters. The combined morphological and genetic evidence suggests these two lineages are on distinct evolutionary trajectories and could be considered distinct subspecies despite their parapatry. Further analysis of the continental shelf cluster resulted in three groups: animals inhabiting shelf waters in the western North Atlantic, the eastern Gulf of Mexico and the western Gulf of Mexico. Analyses of environmental data indicate the four genetic clusters inhabit distinct habitats in terms of depth and sea surface temperature. Contemporary dispersal rate estimates suggest all of these populations should be considered as distinct management units. Conversely, no significant genetic differentiation was observed between S. frontalis from offshore waters of the western North Atlantic and the Azores, which are separated by approximately 4500 km. Overall, the hierarchical structure observed within the Atlantic spotted dolphin shows that the biogeography of the species is complex because it is not shaped solely by geographic distance. [ABSTRACT FROM AUTHOR]

Understanding the degree of genetic exchange between subspecies and populations is vital for the appropriate management of endangered species. Blue whales ( Balaenoptera musculus) have two recognized Southern Hemisphere subspecies that show differences in geographic distribution, morphology, vocalizations and genetics. During the austral summer feeding season, the Antarctic blue whale ( B. m. intermedia) is found in polar waters and the pygmy blue whale ( B. m. brevicauda) in temperate waters. Here, we genetically analyzed samples collected during the feeding season to report on several cases of hybridization between the two recognized blue whale Southern Hemisphere subspecies in a previously unconfirmed sympatric area off Antarctica. This means the pygmy blue whales using waters off Antarctica may migrate and then breed during the austral winter with the Antarctic subspecies. Alternatively, the subspecies may interbreed off Antarctica outside the expected austral winter breeding season. The genetically estimated recent migration rates from the pygmy to Antarctic subspecies were greater than estimates of evolutionary migration rates and previous estimates based on morphology of whaling catches. This discrepancy may be due to differences in the methods or an increase in the proportion of pygmy blue whales off Antarctica within the last four decades. Potential causes for the latter are whaling, anthropogenic climate change or a combination of these and may have led to hybridization between the subspecies. Our findings challenge the current knowledge about the breeding behaviour of the world's largest animal and provide key information that can be incorporated into management and conservation practices for this endangered species. [ABSTRACT FROM AUTHOR]

The role of Beringia as a refugium and route for trans-continental exchange of fauna during glacial cycles of the past 2 million years are well documented; less apparent is its contribution as a significant reservoir of genetic diversity. Using mitochondrial DNA sequences and 14 microsatellite loci, we investigate the phylogeographic history of caribou ( Rangifer tarandus) in western North America. Patterns of genetic diversity reveal two distinct groups of caribou. Caribou classified as a Northern group, of Beringian origin, exhibited greater number and variability in mtDNA haplotypes compared to a Southern group originating from refugia south of glacial ice. Results indicate that subspecies R. t. granti of Alaska and R. t. groenlandicus of northern Canada do not constitute distinguishable units at mtDNA or microsatellites, belying their current status as separate subspecies. Additionally, the Northern Mountain ecotype of woodland caribou (presently R. t. caribou) has closer kinship to caribou classified as granti or groenlandicus. Comparisons of mtDNA and microsatellite data suggest that behavioural and ecological specialization is a more recently derived life history characteristic. Notably, microsatellite differentiation among Southern herds is significantly greater, most likely as a result of human-induced landscape fragmentation and genetic drift due to smaller population sizes. These results not only provide important insight into the evolutionary history of northern species such as caribou, but also are important indicators for managers evaluating conservation measures for this threatened species. [ABSTRACT FROM AUTHOR]

The cheetah ( Acinonyx jubatus) has been described as a species with low levels of genetic variation. This has been suggested to be the consequence of a demographic bottleneck 10 000-12 000 years ago (ya) and also led to the assumption that only small genetic differences exist between the described subspecies. However, analysing mitochondrial DNA and microsatellites in cheetah samples from most of the historic range of the species we found relatively deep phylogeographic breaks between some of the investigated populations, and most of the methods assessed divergence time estimates predating the postulated bottleneck. Mitochondrial DNA monophyly and overall levels of genetic differentiation support the distinctiveness of Northern-East African cheetahs ( Acinonyx jubatus soemmeringii). Moreover, combining archaeozoological and contemporary samples, we show that Asiatic cheetahs ( Acinonyx jubatus venaticus) are unambiguously separated from African subspecies. Divergence time estimates from mitochondrial and nuclear data place the split between Asiatic and Southern African cheetahs ( Acinonyx jubatus jubatus) at 32 000-67 000 ya using an average mammalian microsatellite mutation rate and at 4700-44 000 ya employing human microsatellite mutation rates. Cheetahs are vulnerable to extinction globally and critically endangered in their Asiatic range, where the last 70-110 individuals survive only in Iran. We demonstrate that these extant Iranian cheetahs are an autochthonous monophyletic population and the last representatives of the Asiatic subspecies A. j. venaticus. We advocate that conservation strategies should consider the uncovered independent evolutionary histories of Asiatic and African cheetahs, as well as among some African subspecies. This would facilitate the dual conservation priorities of maintaining locally adapted ecotypes and genetic diversity. [ABSTRACT FROM AUTHOR]

Instances of hybridization between mammalian taxa in the wild are rarely documented. To test for introgression between sibling species of horseshoe bat ( Rhinolophus yunanensis and R. pearsoni) and two subspecies of the latter ( R. p. pearsoni and R. p. chinensis), we sequenced two mtDNA and two ncDNA markers in individuals sampled from multiple localities within their overlapping ranges. The interspecific mtDNA gene tree corresponded to the expected taxonomic divisions, and coalescent-based analyses suggested divergence occurred around 4 MYA. However, these relationships strongly conflicted with those recovered from two independent nuclear gene trees, in which R. yunanensis clustered with R. p. pearsoni to the exclusion of R. p. chinensis. This geographically widespread discordance is best explained by large-scale historical introgression of ncDNA from R. yunanensis to R. pearsoni by male-mediated exchange in mixed species colonies during Pleistocene glacial periods, when ranges may have contracted and overlapped more than at present. Further species tree–gene tree conflicts were detected between R. p. pearsoni and R. p. chinensis, also indicating past and/or current introgression in their overlapping regions. However, here the patterns point to asymmetric mtDNA introgression without ncDNA introgression. Analyses of coalescence times indicate this exchange has occurred subsequent to the divergence of these subspecies from their common ancestor. Our work highlights the importance of using multiple data sets for reconstructing phylogeographic histories and resolving taxonomic relationships. [ABSTRACT FROM AUTHOR]

Although new empirical evidence shows that sympatric speciation has occurred in some species, there are few indisputable model organisms for this process of speciation. The two subspecies ( Gymnocypris eckloni eckloni and G. e. scoliostomus) of the schizothoracine Gymnocypris fish species complex from a small glacier lake in the Tibetan Plateau, Lake Sunmcuo, fit several of the key characteristics of the sympatric speciation model. We used combined mitochondrial control region sequences and the cytochrome b gene (1894 bp) to address the phylogenetics and population genetics of 232 specimens of G. e. eckloni and G. e. scoliostomus, as well as all of its closely related sister species. We found that: (i) a total of four old lineages were uncovered in the widespread G. e. eckloni, of which only one was shown to be shared with all G. e. scoliostomus individuals and (ii) the new subspecies ( G. e. scoliostomus) evolved in Lake Sunmcuo from the ancestral G. e. eckloni population within approximately 0.057 Ma. These two taxa of the species complex are morphologically distinct, and reproductive isolation is further suggested. Ecological disruptive selection based on morphological traits (e.g. mouth cleft characters) and food utilization may be a mechanism of incipient speciation of two sympatric populations within Lake Sunmcuo. This study provides the first genetic evidence for sympatric speciation in the schizothoracine fish. [ABSTRACT FROM AUTHOR]

The red squirrel ( Sciurus vulgaris) is a well-known forest animal distributed all over Europe. Still, we are far from having a firm knowledge of the species’ phylogeography. This study investigates the genetic differentiation of S. vulgaris across the species’ Eurasian range, using sequence data from the mitochondrial DNA gene (D-loop, 252 base pairs, cytochrome b, 359 base pairs), and eight variable autosomal microsatellite loci genotyped for 236 individuals. The results reveal the presence of two main mitochondrial phylogroups. The first clade comprises the individuals from the region of Calabria in southern Italy, belonging to the subspecies S. v. meridionalis, while the second clade contains the remainder of the studied individuals. Bayesian analysis of microsatellite genotypes resulted in three main clusterings corresponding to the three S. vulgaris subspecies: infuscatus, meridionalis and fuscoater. Geographical distribution of mtDNA haplotypes and mismatch analysis suggest a common refugium for the red squirrel across most of its present range from which expansion happened rather rapidly. The genotype mixing of italicus with northern populations could be a residual of postglacial expansion. The lack of mixing between the Calabrian lineage and the rest of European red squirrel haplotypes can be seen as evidence for distinct histories throughout the Pleistocene. Calabrian mtDNA probably diverged in an ice age contraction and remained isolated from the neighbouring squirrel populations until very recent times. [ABSTRACT FROM AUTHOR]

An understanding of the role of factors intrinsic to a species’ life history in structuring contemporary genetic variation is a fundamental, but understudied, aspect of evolutionary biology. Here, we assessed the influence of the propensity to outbreak in shaping worldwide genetic variation in Locusta migratoria, a cosmopolitan pest well known for its expression of density-dependent phase polyphenism. We scored 14 microsatellites in nine subspecies from 25 populations distributed over most of the species’ range in regions that vary in the historical frequency and extent of their outbreaks. We rejected the hypothesis that L. migratoria consists of two genetically distinct clusters adapted to habitats either rarely (nonoutbreaking) or cyclically (outbreaking) favourable to increases in population density. We also invalidated the current subspecific taxonomic classification based on morphometrics. Bayesian inferences indicated evidence of a homogenizing effect of outbreaks on L. migratoria population structure. Geographical and ecological barriers to gene flow in conjunction with historical events can also explain the observed patterns. By systematically assessing the effects of null alleles using computer simulations, we also provide a template for the analysis of microsatellite data sets characterized by a high prevalence of null alleles. [ABSTRACT FROM AUTHOR]