Complexity of avian evolution revealed by family-level genomes.

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

Site-directed genotype screening for elimination of antinutritional saponins in quinoa seeds identifies TSARL1 as a master controller of saponin biosynthesis selectively in seeds.

Climate change may result in a drier climate and increased salinization, threatening agricultural productivity worldwide. Quinoa (Chenopodium quinoa) produces highly nutritious seeds and tolerates abiotic stresses such as drought and high salinity, making it a promising future food source. However, the presence of antinutritional saponins in their seeds is an undesirable trait. We mapped genes controlling seed saponin content to a genomic region that includes TSARL1. We isolated desired genetic variation in this gene by producing a large mutant library of a commercial quinoa cultivar and screening the library for specific nucleotide substitutions using droplet digital PCR. We were able to rapidly isolate two independent tsarl1 mutants, which retained saponins in the leaves and roots for defence, but saponins were undetectable in the seed coat. We further could show that TSARL1 specifically controls seed saponin biosynthesis in the committed step after 2,3-oxidosqualene. Our work provides new important knowledge on the function of TSARL1 and represents a breakthrough for quinoa breeding.

A new nuclear phylogeny of the tea family (Theaceae) unravels rapid radiations in genus Camellia.

Molecular analyses of rapidly radiating groups often reveal incongruence between gene trees. This mainly results from incomplete lineage sorting, introgression, and gene tree estimation error, which complicate the estimation of phylogenetic relationships. In this study, we reconstruct the phylogeny of Theaceae using 348 nuclear loci from 68 individuals and two outgroup taxa. Sequence data were obtained by target enrichment using the recently released Angiosperm 353 universal probe set applied to herbarium specimens. The robustness of the topologies to variation in data quality was established under a range of different filtering schemes, using both coalescent and concatenation approaches. Our results confirmed most of the previously hypothesized relationships among tribes and genera, while clarifying additional interspecific relationships within the rapidly radiating genus Camellia. We recovered a remarkably high degree of gene tree heterogeneity indicative of rapid radiation in the group and observed cytonuclear conflicts, especially within Camellia. This was especially pronounced around short branches, which we primarily associate with gene tree estimation error. Our analysis also indicates that incomplete lineage sorting (ILS) contributed to gene-tree conflicts and accounted for approximately 14 % of the explained variation, whereas inferred introgression levels were low. Our study advances the understanding of the evolution of this important plant family and provides guidance on the application of target capture methods and the evaluation of key processes that influence phylogenetic discordances.

Response of an Afro-Palearctic bird migrant to glaciation cycles.

Migration allows animals to exploit spatially separated and seasonally available resources at a continental to global scale. However, responding to global climatic changes might prove challenging, especially for long-distance intercontinental migrants. During glacial periods, when conditions became too harsh for breeding in the north, avian migrants have been hypothesized to retract their distribution to reside within small refugial areas. Here, we present data showing that an Afro-Palearctic migrant continued seasonal migration, largely within Africa, during previous glacial-interglacial cycles with no obvious impact on population size. Using individual migratory track data to hindcast monthly bioclimatic habitat availability maps through the last 120,000 y, we show altered seasonal use of suitable areas through time. Independently derived effective population sizes indicate a growing population through the last 40,000 y. We conclude that the migratory lifestyle enabled adaptation to shifting climate conditions. This indicates that populations of resource-tracking, long-distance migratory species could expand successfully during warming periods in the past, which could also be the case under future climate scenarios.

Comparative genomics provides insights into the aquatic adaptations of mammals.

The ancestors of marine mammals once roamed the land and independently committed to an aquatic lifestyle. These macroevolutionary transitions have intrigued scientists for centuries. Here, we generated high-quality genome assemblies of 17 marine mammals (11 cetaceans and six pinnipeds), including eight assemblies at the chromosome level. Incorporating previously published data, we reconstructed the marine mammal phylogeny and population histories and identified numerous idiosyncratic and convergent genomic variations that possibly contributed to the transition from land to water in marine mammal lineages. Genes associated with the formation of blubber (NFIA), vascular development (SEMA3E), and heat production by brown adipose tissue (UCP1) had unique changes that may contribute to marine mammal thermoregulation. We also observed many lineage-specific changes in the marine mammals, including genes associated with deep diving and navigation. Our study advances understanding of the timing, pattern, and molecular changes associated with the evolution of mammalian lineages adapting to aquatic life.

Identifying loci under positive selection in complex population histories.

Detailed modeling of a species' history is of prime importance for understanding how natural selection operates over time. Most methods designed to detect positive selection along sequenced genomes, however, use simplified representations of past histories as null models of genetic drift. Here, we present the first method that can detect signatures of strong local adaptation across the genome using arbitrarily complex admixture graphs, which are typically used to describe the history of past divergence and admixture events among any number of populations. The method-called graph-aware retrieval of selective sweeps (GRoSS)-has good power to detect loci in the genome with strong evidence for past selective sweeps and can also identify which branch of the graph was most affected by the sweep. As evidence of its utility, we apply the method to bovine, codfish, and human population genomic data containing panels of multiple populations related in complex ways. We find new candidate genes for important adaptive functions, including immunity and metabolism in understudied human populations, as well as muscle mass, milk production, and tameness in specific bovine breeds. We are also able to pinpoint the emergence of large regions of differentiation owing to inversions in the history of Atlantic codfish.

Using ultraconserved elements to track the influence of sea-level change on leafy seadragon populations.

During the Last Glacial Maximum (LGM), global sea levels were 120-130 m lower than today, resulting in the emergence of most continental shelves and extirpation of subtidal organisms from these areas. During the interglacial periods, rapid inundation of shelf regions created a dynamic environment for coastal organisms, such as the charismatic leafy seadragon (Phycodurus eques, Syngnathidae), a brooder with low dispersal ability inhabiting kelp beds in temperate Australia. Reconstructions of the palaeoshoreline revealed that the increase of shallow areas since the LGM was not uniform across the species' range and we investigated the effects of these asymmetries on genetic diversity and structuring. Using targeted capture of 857 variable ultraconserved elements (UCEs, 2,845 single nucleotide polymorphisms) in 68 individuals, we found that the regionally different shelf topographies were paralleled by contrasting population genetic patterns. In the west, populations may not have persisted through sea-level lows because shallow seabed was very limited. Shallow genetic structure, weak expansion signals and a westward cline in genetic diversity indicate a postglacial recolonization of the western part of the range from a more eastern location following sea-level rise. In the east, shallow seabed persisted during the LGM and increased considerably after the flooding of large bays, which resulted in strong demographic expansions, deeper genetic structure and higher genetic diversity. This study suggests that postglacial flooding with rising sea levels produced locally variable signatures in colonizing populations.

Complete Inactivation of Sebum-Producing Genes Parallels the Loss of Sebaceous Glands in Cetacea.

Genomes are dynamic biological units, with processes of gene duplication and loss triggering evolutionary novelty. The mammalian skin provides a remarkable case study on the occurrence of adaptive morphological innovations. Skin sebaceous glands (SGs), for instance, emerged in the ancestor of mammals serving pivotal roles, such as lubrication, waterproofing, immunity, and thermoregulation, through the secretion of sebum, a complex mixture of various neutral lipids such as triacylglycerol, free fatty acids, wax esters, cholesterol, and squalene. Remarkably, SGs are absent in a few mammalian lineages, including the iconic Cetacea. We investigated the evolution of the key molecular components responsible for skin sebum production: Dgat2l6, Awat1, Awat2, Elovl3, Mogat3, and Fabp9. We show that all analyzed genes have been rendered nonfunctional in Cetacea species (toothed and baleen whales). Transcriptomic analysis, including a novel skin transcriptome from blue whale, supports gene inactivation. The conserved mutational pattern found in most analyzed genes, indicates that pseudogenization events took place prior to the diversification of modern Cetacea lineages. Genome and skin transcriptome analysis of the common hippopotamus highlighted the convergent loss of a subset of sebum-producing genes, notably Awat1 and Mogat3. Partial loss profiles were also detected in non-Cetacea aquatic mammals, such as the Florida manatee, and in terrestrial mammals displaying specialized skin phenotypes such as the African elephant, white rhinoceros and pig. Our findings reveal a unique landscape of "gene vestiges" in the Cetacea sebum-producing compartment, with limited gene loss observed in other mammalian lineages: suggestive of specific adaptations or specializations of skin lipids.

Cartilaginous fishes offer unique insights into the evolution of the nuclear receptor gene repertoire in gnathostomes.

Nuclear receptors (NRs) are key transcription factors that originated in the common ancestor of metazoans. The vast majority of NRs are triggered by binding to either endogenous (e.g. retinoic acid) or exogenous (e.g. xenobiotics) ligands, and their evolution and expansion is tightly linked to the function of endocrine systems. Importantly, they represent classic targets of physiological exploitation by endocrine disrupting chemicals. The NR gene repertoire in different lineages has been shaped by gene loss, duplication and mutation, denoting a dynamic evolutionary route. As the earliest diverging class of gnathostomes (jawed vertebrates), cartilaginous fishes offer an exceptional opportunity to address the early diversification of NR gene families and the evolution of the endocrine system in jawed vertebrates. Here we provide an exhaustive analysis into the NR gene composition in five elasmobranch (sharks and rays) and two holocephalan (chimaeras) species. For this purpose, we generated also a low coverage draft genome assembly of the chimaera small-eyed rabbitfish, Hydrolagus affinis. We show that cartilaginous fish retain an archetypal NR gene repertoire, similar to that of mammals and coincident with the two rounds of whole genome duplication that occurred in the gnathostome ancestor. Furthermore, novel gene members of the non-canonical NR0B receptors were found in the genomes of this lineage. Our findings provide an essential view into the early diversification of NRs in gnathostomes, paving the way for functional studies.

Consequences of breed formation on patterns of genomic diversity and differentiation: the case of highly diverse peripheral Iberian cattle.

BACKGROUND: Iberian primitive breeds exhibit a remarkable phenotypic diversity over a very limited geographical space. While genomic data are accumulating for most commercial cattle, it is still lacking for these primitive breeds. Whole genome data is key to understand the consequences of historic breed formation and the putative role of earlier admixture events in the observed diversity patterns. RESULTS: We sequenced 48 genomes belonging to eight Iberian native breeds and found that the individual breeds are genetically very distinct with FST values ranging from 4 to 16% and have levels of nucleotide diversity similar or larger than those of their European counterparts, namely Jersey and Holstein. All eight breeds display significant gene flow or admixture from African taurine cattle and include mtDNA and Y-chromosome haplotypes from multiple origins. Furthermore, we detected a very low differentiation of chromosome X relative to autosomes within all analyzed taurine breeds, potentially reflecting male-biased gene flow. CONCLUSIONS: Our results show that an overall complex history of admixture resulted in unexpectedly high levels of genomic diversity for breeds with seemingly limited geographic ranges that are distantly located from the main domestication center for taurine cattle in the Near East. This is likely to result from a combination of trading traditions and breeding practices in Mediterranean countries. We also found that the levels of differentiation of autosomes vs sex chromosomes across all studied taurine and indicine breeds are likely to have been affected by widespread breeding practices associated with male-biased gene flow.

31° South: The physiology of adaptation to arid conditions in a passerine bird.

Arid environments provide ideal ground for investigating the mechanisms of adaptive evolution. High temperatures and low water availability are relentless stressors for many endotherms, including birds; yet birds persist in deserts. While physiological adaptation probably involves metabolic phenotypes, the underlying mechanisms (plasticity, genetics) are largely uncharacterized. To explore this, we took an intraspecific approach that focused on a species that is resident over a mesic to arid gradient, the Karoo scrub-robin (Cercotrichas coryphaeus). Specifically, we integrated environmental (climatic and primary productivity), physiological (metabolic rates: a measure of energy expenditure), genotypic (genetic variation underlying the machinery of energy production) and microbiome (involved in processing food from where energy is retrieved) data, to infer the mechanism of physiological adaptation. We that found the variation in energetic physiology phenotypes and gut microbiome composition are associated with environmental features as well as with variation in genes underlying energy metabolic pathways. Specifically, we identified a small list of candidate adaptive genes, some of them with known ties to relevant physiology phenotypes. Together our results suggest that selective pressures on energetic physiology mediated by genes related to energy homeostasis and possibly microbiota composition may facilitate adaptation to local conditions and provide an explanation to the high avian intraspecific divergence observed in harsh environments.

Genome sequence and effectorome of Moniliophthora perniciosa and Moniliophthora roreri subpopulations.

BACKGROUND: The hemibiotrophic pathogens Moniliophthora perniciosa (witches' broom disease) and Moniliophthora roreri (frosty pod rot disease) are among the most important pathogens of cacao. Moniliophthora perniciosa has a broad host range and infects a variety of meristematic tissues in cacao plants, whereas M. roreri infects only pods of Theobroma and Herrania genera. Comparative pathogenomics of these fungi is essential to understand Moniliophthora infection strategies, therefore the detection and in silico functional characterization of effector candidates are important steps to gain insight on their pathogenicity. RESULTS: Candidate secreted effector proteins repertoire were predicted using the genomes of five representative isolates of M. perniciosa subpopulations (three from cacao and two from solanaceous hosts), and one representative isolate of M. roreri from Peru. Many putative effectors candidates were identified in M. perniciosa: 157 and 134 in cacao isolates from Bahia, Brazil; 109 in cacao isolate from Ecuador, 92 and 80 in wild solanaceous isolates from Minas Gerais (Lobeira) and Bahia (Caiçara), Brazil; respectively. Moniliophthora roreri showed the highest number of effector candidates, a total of 243. A set of eight core effectors were shared among all Moniliophthora isolates, while others were shared either between the wild solanaceous isolates or among cacao isolates. Mostly, candidate effectors of M. perniciosa were shared among the isolates, whereas in M. roreri nearly 50% were exclusive to the specie. In addition, a large number of cell wall-degrading enzymes characteristic of hemibiotrophic fungi were found. From these, we highlighted the proteins involved in cell wall modification, an enzymatic arsenal that allows the plant pathogens to inhabit environments with oxidative stress, which promotes degradation of plant compounds and facilitates infection. CONCLUSIONS: The present work reports six genomes and provides a database of the putative effectorome of Moniliophthora, a first step towards the understanding of the functional basis of fungal pathogenicity.

Molecular clocks indicate turnover and diversification of modern coleoid cephalopods during the Mesozoic Marine Revolution.

Coleoid cephalopod molluscs comprise squid, cuttlefish and octopuses, and represent nearly the entire diversity of modern cephalopods. Sophisticated adaptations such as the use of colour for camouflage and communication, jet propulsion and the ink sac highlight the unique nature of the group. Despite these striking adaptations, there are clear parallels in ecology between coleoids and bony fishes. The coleoid fossil record is limited, however, hindering confident analysis of the tempo and pattern of their evolution. Here we use a molecular dataset (180 genes, approx. 36 000 amino acids) of 26 cephalopod species to explore the phylogeny and timing of cephalopod evolution. We show that crown cephalopods diverged in the Silurian-Devonian, while crown coleoids had origins in the latest Palaeozoic. While the deep-sea vampire squid and dumbo octopuses have ancient origins extending to the Early Mesozoic Era, 242 ± 38 Ma, incirrate octopuses and the decabrachian coleoids (10-armed squid) diversified in the Jurassic Period. These divergence estimates highlight the modern diversity of coleoid cephalopods emerging in the Mesozoic Marine Revolution, a period that also witnessed the radiation of most ray-finned fish groups in addition to several other marine vertebrates. This suggests that that the origin of modern cephalopod biodiversity was contingent on ecological competition with marine vertebrates.

Combined Transcriptomic and Proteomic Analysis of the Posterior Salivary Gland from the Southern Blue-Ringed Octopus and the Southern Sand Octopus.

This study provides comprehensive proteomic profiles from the venom producing posterior salivary glands of octopus (superorder Octopodiformes) species. A combined transcriptomic and proteomic approach was used to identify 1703 proteins from the posterior salivary gland of the southern blue-ringed octopus, Hapalochlaena maculosa and 1300 proteins from the posterior salivary gland of the southern sand octopus, Octopus kaurna. The two proteomes were broadly similar; clustering of proteins into orthogroups revealed 937 that were shared between species. Serine proteases were particularly diverse and abundant in both species. Other abundant proteins included a large number of secreted proteins, many of which had no known conserved domains, or homology to proteins with known function. On the basis of homology to known venom proteins, 23 putative toxins were identified in H. maculosa and 24 in O. kaurna. These toxins span nine protein families: CAP (cysteine rich secretory proteins, antigen 5, parthenogenesis related), chitinase, carboxylesterase, DNase, hyaluronidase, metalloprotease, phospholipase, serine protease and tachykinin. Serine proteases were responsible for 70.9% and 86.3% of putative toxin expression in H. maculosa and O. kaurna, respectively, as determined using intensity based absolute quantification (iBAQ) measurements. Phylogenetic analysis of the putative toxin serine proteases revealed a similar suite of diverse proteins present in both species. Posterior salivary gland composition of H. maculosa and O. kaurna differ in several key aspects. While O. kaurna expressed the proteinaceous neurotoxin, tachykinin, this was absent from H. maculosa, perhaps reflecting the acquisition of a potent nonproteinaceous neurotoxin, tetrodotoxin (TTX) produced by bacteria in the salivary glands of that species. The dispersal factor, hyaluronidase was particularly abundant in H. maculosa. Chitinase was abundant in both species and is believed to facilitate envenomation in chitinous prey such as crustaceans. Cephalopods represent a largely unexplored source of novel proteins distinct from all other venomous taxa and are of interest for further inquiry, as novel proteinaceous toxins derived from venoms may contribute to pharmaceutical design.

Comparative analysis of complete mitochondrial genomes suggests that relaxed purifying selection is driving high nonsynonymous evolutionary rate of the NADH2 gene in whitefish (Coregonus ssp.).

Several studies have recently reported evidence for positive selection acting on the mitochondrial genome (mitogenome), emphasizing its potential role in adaptive divergence and speciation. In this study we searched 107 full mitogenomes of recently diverged species and lineages of whitefish (Coregonus ssp.) for signals of positive selection. These salmonids show several distinct morphological and ecological differences that may be associated with energetics and therefore potentially positive selection at the mitogenome level. We found that purifying selection and genetic drift were the predominant evolutionary forces acting on the analyzed mitogenomes. However, the NADH dehydrogenase 2 gene (ND2) showed a highly elevated dN/dS ratio compared to the other mitochondrial genes, which was significantly higher in whitefish compared to other salmonids. We therefore further examined nonsynonymous evolution in ND2 by (i) mapping amino acid changes to a protein model structure which showed that they were located away from key functional residues of the protein, (ii) locating them in the sequences of other species of fish (Salmonidae, Anguillidae, Scombridae and Percidae) only to find pronounced overlap of nonsynonymous regions. We thus conclude that relaxed purifying selection is driving the evolution of ND2 by affecting mostly regions that have lower functional relevance.

Shotgun analysis of the marine mussel Mytilus edulis hemolymph proteome and mapping the innate immunity elements.

The marine mussel innate immunity provides protection to pathogen invasion and inflammation. In this regard, the mussel hemolymph takes a main role in the animal innate response. Despite the importance of this body fluid in determining the physiological condition of the animal, little is known about the molecular mechanisms underlying the cellular and humoral responses. In this work, we have applied a MS (nano-LC-MS/MS) strategy integrating genomic and transcriptomic data with the aim to: (i) identify the main protein functional groups that characterize hemolymph and (ii) to map the elements of innate immunity in the marine mussel Mytilus edulis hemolymph proteome. After sample analysis and first protein identification based on MS/MS data comparison, proteins with unknown functions were annotated with blast using public database (nrNCBI) information. Overall 595 hemolymph proteins were identified with high confidence and annotated. These proteins encompass primary cellular metabolic processes: energy production and metabolism of biomolecules, as well as processes related to oxidative stress defence, xenobiotic detoxification, drug metabolism, and immune response. A group of proteins was identified with putative immune effector, receptor, and signaling functions in M. edulis. Data are available via ProteomeXchange with identifier PXD001951 (http://proteomecentral.proteomexchange.org/dataset/PXD001951).

A refined model of the genomic basis for phenotypic variation in vertebrate hemostasis.

BACKGROUND: Hemostasis is a defense mechanism that enhances an organism's survival by minimizing blood loss upon vascular injury. In vertebrates, hemostasis has been evolving with the cardio-vascular and hemodynamic systems over the last 450 million years. Birds and mammals have very similar vascular and hemodynamic systems, thus the mechanism that blocks ruptures in the vasculature is expected to be the same. However, the speed of the process varies across vertebrates, and is particularly slow for birds. Understanding the differences in the hemostasis pathway between birds and mammals, and placing them in perspective to other vertebrates may provide clues to the genetic contribution to variation in blood clotting phenotype in vertebrates. We compiled genomic data corresponding to key elements involved in hemostasis across vertebrates to investigate its genetic basis and understand how it affects fitness. RESULTS: We found that: i) fewer genes are involved in hemostasis in birds compared to mammals; and ii) the largest differences concern platelet membrane receptors and components from the kallikrein-kinin system. We propose that lack of the cytoplasmic domain of the GPIb receptor subunit alpha could be a strong contributor to the prolonged bleeding phenotype in birds. Combined analysis of laboratory assessments of avian hemostasis with the first avian phylogeny based on genomic-scale data revealed that differences in hemostasis within birds are not explained by phylogenetic relationships, but more so by genetic variation underlying components of the hemostatic process, suggestive of natural selection. CONCLUSIONS: This work adds to our understanding of the evolution of hemostasis in vertebrates. The overlap with the inflammation, complement and renin-angiotensin (blood pressure regulation) pathways is a potential driver of rapid molecular evolution in the hemostasis network. Comparisons between avian species and mammals allowed us to hypothesize that the observed mammalian innovations might have contributed to the diversification of mammals that give birth to live young.

Proteomic profiling of cytosolic glutathione transferases from three bivalve species: Corbicula fluminea, Mytilus galloprovincialis and Anodonta cygnea.

Suspension-feeding bivalves are considered efficient toxin vectors with a relative insensitivity to toxicants compared to other aquatic organisms. This fact highlights the potential role of detoxification enzymes, such as glutathione transferases (GSTs), in this bivalve resistance. Nevertheless, the GST system has not been extensively described in these organisms. In the present study, cytosolic GSTs isoforms (cGST) were surveyed in three bivalves with different habitats and life strategies: Corbicula fluminea, Anodonta cygnea and Mytilus galloprovincialis. GSTs were purified by glutathione-agarose affinity chromatography, and the collection of expressed cGST classes of each bivalve were identified using a proteomic approach. All the purified extracts were also characterized kinetically. Results reveal variations in cGST subunits collection (diversity and properties) between the three tested bivalves. Using proteomics, four pi-class and two sigma-class GST subunits were identified in M. galloprovincialis. C. fluminea also yielded four pi-class and one sigma-class GST subunits. For A. cygnea, two mu-class and one pi-class GST subunits were identified, these being the first record of GSTs from these freshwater mussels. The affinity purified extracts also show differences regarding enzymatic behavior among species. The variations found in cGST collection and kinetics might justify diverse selective advantages for each bivalve organism.

Population Genomics Reveals the Underlying Structure of the Small Pelagic European Sardine and Suggests Low Connectivity within Macaronesia.

The European sardine (Sardina pilchardus, Walbaum 1792) is indisputably a commercially important species. Previous studies using uneven sampling or a limited number of makers have presented sometimes conflicting evidence of the genetic structure of S. pilchardus populations. Here, we show that whole genome data from 108 individuals from 16 sampling areas across 5000 km of the species' distribution range (from the Eastern Mediterranean to the archipelago of Azores) support at least three genetic clusters. One includes individuals from Azores and Madeira, with evidence of substructure separating these two archipelagos in the Atlantic. Another cluster broadly corresponds to the center of the distribution, including the sampling sites around Iberia, separated by the Almeria-Oran front from the third cluster that includes all of the Mediterranean samples, except those from the Alboran Sea. Individuals from the Canary Islands appear to belong to the Mediterranean cluster. This suggests at least two important geographical barriers to gene flow, even though these do not seem complete, with many individuals from around Iberia and the Mediterranean showing some patterns compatible with admixture with other genetic clusters. Genomic regions corresponding to the top outliers of genetic differentiation are located in areas of low recombination indicative that genetic architecture also has a role in shaping population structure. These regions include genes related to otolith formation, a calcium carbonate structure in the inner ear previously used to distinguish S. pilchardus populations. Our results provide a baseline for further characterization of physical and genetic barriers that divide European sardine populations, and information for transnational stock management of this highly exploited species towards sustainable fisheries.

Introgression and disruption of migration routes have shaped the genetic integrity of wildebeest populations.

The blue wildebeest (Connochaetes taurinus) is a keystone species in savanna ecosystems from southern to eastern Africa, and is well known for its spectacular migrations and locally extreme abundance. In contrast, the black wildebeest (C. gnou) is endemic to southern Africa, barely escaped extinction in the 1900s and is feared to be in danger of genetic swamping from the blue wildebeest. Despite the ecological importance of the wildebeest, there is a lack of understanding of how its unique migratory ecology has affected its gene flow, genetic structure and phylogeography. Here, we analyze whole genomes from 121 blue and 22 black wildebeest across the genus' range. We find discrete genetic structure consistent with the morphologically defined subspecies. Unexpectedly, our analyses reveal no signs of recent interspecific admixture, but rather a late Pleistocene introgression of black wildebeest into the southern blue wildebeest populations. Finally, we find that migratory blue wildebeest populations exhibit a combination of long-range panmixia, higher genetic diversity and lower inbreeding levels compared to neighboring populations whose migration has recently been disrupted. These findings provide crucial insights into the evolutionary history of the wildebeest, and tangible genetic evidence for the negative effects of anthropogenic activities on highly migratory ungulates.