Diverse pathways to speciation revealed by marine snails.

Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems.

Chromosomal inversion polymorphisms are widespread across the species ranges of rough periwinkles (Littorina saxatilis and L. arcana).

Inversions are thought to play a key role in adaptation and speciation, suppressing recombination between diverging populations. Genes influencing adaptive traits cluster in inversions, and changes in inversion frequencies are associated with environmental differences. However, in many organisms, it is unclear if inversions are geographically and taxonomically widespread. The intertidal snail, Littorina saxatilis, is one such example. Strong associations between putative polymorphic inversions and phenotypic differences have been demonstrated between two ecotypes of L. saxatilis in Sweden and inferred elsewhere, but no direct evidence for inversion polymorphism currently exists across the species range. Using whole genome data from 107 snails, most inversion polymorphisms were found to be widespread across the species range. The frequencies of some inversion arrangements were significantly different among ecotypes, suggesting a parallel adaptive role. Many inversions were also polymorphic in the sister species, L. arcana, hinting at an ancient origin.

How chromosomal inversions reorient the evolutionary process.

Inversions are structural mutations that reverse the sequence of a chromosome segment and reduce the effective rate of recombination in the heterozygous state. They play a major role in adaptation, as well as in other evolutionary processes such as speciation. Although inversions have been studied since the 1920s, they remain difficult to investigate because the reduced recombination conferred by them strengthens the effects of drift and hitchhiking, which in turn can obscure signatures of selection. Nonetheless, numerous inversions have been found to be under selection. Given recent advances in population genetic theory and empirical study, here we review how different mechanisms of selection affect the evolution of inversions. A key difference between inversions and other mutations, such as single nucleotide variants, is that the fitness of an inversion may be affected by a larger number of frequently interacting processes. This considerably complicates the analysis of the causes underlying the evolution of inversions. We discuss the extent to which these mechanisms can be disentangled, and by which approach.

Using replicate hybrid zones to understand the genomic basis of adaptive divergence.

Combining hybrid zone analysis with genomic data is a promising approach to understanding the genomic basis of adaptive divergence. It allows for the identification of genomic regions underlying barriers to gene flow. It also provides insights into spatial patterns of allele frequency change, informing about the interplay between environmental factors, dispersal and selection. However, when only a single hybrid zone is analysed, it is difficult to separate patterns generated by selection from those resulting from chance. Therefore, it is beneficial to look for repeatable patterns across replicate hybrid zones in the same system. We applied this approach to the marine snail Littorina saxatilis, which contains two ecotypes, adapted to wave-exposed rocks vs. high-predation boulder fields. The existence of numerous hybrid zones between ecotypes offered the opportunity to test for the repeatability of genomic architectures and spatial patterns of divergence. We sampled and phenotyped snails from seven replicate hybrid zones on the Swedish west coast and genotyped them for thousands of single nucleotide polymorphisms. Shell shape and size showed parallel clines across all zones. Many genomic regions showing steep clines and/or high differentiation were shared among hybrid zones, consistent with a common evolutionary history and extensive gene flow between zones, and supporting the importance of these regions for divergence. In particular, we found that several large putative inversions contribute to divergence in all locations. Additionally, we found evidence for consistent displacement of clines from the boulder-rock transition. Our results demonstrate patterns of spatial variation that would not be accessible without continuous spatial sampling, a large genomic data set and replicate hybrid zones.

Speciation in marine environments: Diving under the surface.

Marine environments are inhabited by a broad representation of the tree of life, yet our understanding of speciation in marine ecosystems is extremely limited compared with terrestrial and freshwater environments. Developing a more comprehensive picture of speciation in marine environments requires that we 'dive under the surface' by studying a wider range of taxa and ecosystems is necessary for a more comprehensive picture of speciation. Although studying marine evolutionary processes is often challenging, recent technological advances in different fields, from maritime engineering to genomics, are making it increasingly possible to study speciation of marine life forms across diverse ecosystems and taxa. Motivated by recent research in the field, including the 14 contributions in this issue, we highlight and discuss six axes of research that we think will deepen our understanding of speciation in the marine realm: (a) study a broader range of marine environments and organisms; (b) identify the reproductive barriers driving speciation between marine taxa; (c) understand the role of different genomic architectures underlying reproductive isolation; (d) infer the evolutionary history of divergence using model-based approaches; (e) study patterns of hybridization and introgression between marine taxa; and (f) implement highly interdisciplinary, collaborative research programmes. In outlining these goals, we hope to inspire researchers to continue filling this critical knowledge gap surrounding the origins of marine biodiversity.

Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe.

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat-related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model-based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.

Perfume and Flavor Engineering: A Chemical Engineering Perspective.

In the last two decades, scientific methodologies for the prediction of the design, performance and classification of fragrance mixtures have been developed at the Laboratory of Separation and Reaction Engineering. This review intends to give an overview of such developments. It all started with the question: what do we smell? The Perfumery Ternary Diagram enables us to determine the dominant odor for each perfume composition. Evaporation and 1D diffusion model is analyzed based on vapor-liquid equilibrium and Fick's law for diffusion giving access to perfume performance parameters. The effect of matrix and skin is addressed and the trail of perfumes analyzed. Classification of perfumes with the perfumery radar is discussed. The methodology is extended to flavor and taste engineering. Finally, future research directions are suggested.

Membrane-Supported Layered Coordination Polymer as an Advanced Sustainable Catalyst for Desulfurization.

The application of a catalytic membrane in the oxidative desulfurization of a multicomponent model diesel formed by most refractory sulfur compounds present in fuel is reported here for the first time. The catalytic membrane was prepared by the impregnation of the active lamellar [Gd(H4nmp)(H2O)2]Cl·2H2O (UAV-59) coordination polymer (CP) into a polymethyl methacrylate (PMMA, acrylic glass) supporting membrane. The use of the catalytic membrane in the liquid-liquid system instead of a powder catalyst arises as an enormous advantage associated with the facility of catalyst handling while avoiding catalyst mass loss. The optimization of various parameters allowed to achieve a near complete desulfurization after 3 h under sustainable conditions, i.e., using an aqueous H2O2 as oxidant and an ionic liquid as extraction solvent ([BMIM]PF6, 1:0.5 ratio diesel:[BMIM]PF6). The performance of the catalytic membrane and of the powdered UAV-59 catalyst was comparable, with the advantage that the former could be recycled successfully for a higher number of desulfurization cycles without the need of washing and drying procedures between reaction cycles, turning the catalytic membrane process more cost-efficient and suitable for future industrial application.

Phylogeographic history of flat periwinkles, Littorina fabalis and L. obtusata.

BACKGROUND: The flat periwinkles, Littorina fabalis and L. obtusata, are two sister species widely distributed throughout the Northern Atlantic shores with high potential to inform us about the process of ecological speciation in the intertidal. However, whether gene flow has occurred during their divergence is still a matter of debate. A comprehensive assessment of the genetic diversity of these species is also lacking and their main glacial refugia and dispersal barriers remain largely unknown. In order to fill these gaps, we sequenced two mitochondrial genes and two nuclear fragments to perform a phylogeographic analysis of flat periwinkles across their distribution range. RESULTS: We identified two main clades largely composed by species-specific haplotypes corresponding to L. obtusata and L. fabalis, with moderate to strong support, respectively. Importantly, a model of divergence with gene flow between the two species (from L. obtusata to L. fabalis) was better supported, both in Iberia and in northern-central Europe. Three mitochondrial clades were detected within L. fabalis and two within L. obtusata, with strong divergence between Iberia and the remaining populations. The largest component of the genetic variance within each species was explained by differences between geographic regions associated with these clades. Our data suggests that overall intraspecific genetic diversity is similar between the two flat periwinkle species and that populations from Iberia tend to be less diverse than populations from northern-central Europe. CONCLUSIONS: The phylogeographic analysis of this sister-species pair supports divergence with gene flow. This system thus provides us with the opportunity to study the contribution of gene flow and natural selection during diversification. The distribution of the different clades suggests the existence of glacial refugia in Iberia and northern-central Europe for both species, with a main phylogeographic break between these regions. Although the genetic diversity results are not fully conclusive, the lower diversity observed in Iberia could reflect marginal conditions at the southern limit of their distribution range during the current interglacial period.

Evolutionary history of two cryptic species of northern African jerboas.

BACKGROUND: Climatic variation and geologic change both play significant roles in shaping species distributions, thus affecting their evolutionary history. In Sahara-Sahel, climatic oscillations shifted the desert extent during the Pliocene-Pleistocene interval, triggering the diversification of several species. Here, we investigated how these biogeographical and ecological events have shaped patterns of genetic diversity and divergence in African Jerboas, desert specialist rodents. We focused on two sister and cryptic species, Jaculus jaculus and J. hirtipes, where we (1) evaluated their genetic differentiation, (2) reconstructed their evolutionary and demographic history; (3) tested the level of gene flow between them, and (4) assessed their ecological niche divergence. RESULTS: The analyses based on 231 individuals sampled throughout North Africa, 8 sequence fragments (one mitochondrial and seven single copy nuclear DNA, including two candidate genes for fur coloration: MC1R and Agouti), 6 microsatellite markers and ecological modelling revealed: (1) two distinct genetic lineages with overlapping distributions, in agreement with their classification as different species, J. jaculus and J. hirtipes, with (2) low levels of gene flow and strong species divergence, (3) high haplotypic diversity without evident geographic structure within species, and (4) a low level of large-scale ecological divergence between the two taxa, suggesting species micro-habitat specialization. CONCLUSIONS: Overall, our results suggest a speciation event that occurred during the Pliocene-Pleistocene transition. The contemporary distribution of genetic variation suggests ongoing population expansions. Despite the largely overlapping distributions at a macrogeographic scale, our genetic results suggest that the two species remain reproductively isolated, as only negligible levels of gene flow were observed. The overlapping ecological preferences at a macro-geographic scale and the ecological divergence at the micro-habitat scale suggest that local adaptation may have played a crucial role in the speciation process of these species.

Predictors of Fabry disease in patients with hypertrophic cardiomyopathy: How to guide the diagnostic strategy?

BACKGROUND: Fabry disease (FD) is a treatable cause of hypertrophic cardiomyopathy (HCM). We aimed to determine the independent predictors of FD and to define a clinically useful strategy to discriminate FD among HCM. METHODS: Multicenter study including 780 patients with the ESC definition of HCM. FD screening was performed by enzymatic assay in males and genetic testing in females. Multivariate regression analysis identified independent predictors of FD in HCM. A discriminant function analysis defined a score based on the weighted combination of these predictors. RESULTS: FD was found in 37 of 780 patients with HCM (4.7%): 31 with p.F113L mutation due to a founder effect; and 6 with other variants (p.C94S; p.M96V; p.G183V; p.E203X; p.M290I; p.R356Q/p.G360R). FD prevalence in HCM adjusted for the founder effect was 0.9%. Symmetric HCM (OR 3.464, CI95% 1.151-10.430), basal inferolateral late gadolinium enhancement (LGE) (OR 10.677, CI95% 3.633-31.380), bifascicular block (OR 10.909, CI95% 2.377-50.059) and ST-segment depression (OR 4.401, CI95% 1.431-13.533) were independent predictors of FD in HCM. The score ID FABRY-HCM [-0.729 + (2.781xBifascicular block) + (0.590xST depression) + (0.831xSymmetric HCM) + (2.130xbasal inferolateral LGE)] had a negative predictive value of 95.8% for FD, with a cut-off of 1.0, meaning that, in the absence of both bifascicular block and basal inferolateral LGE, FD is a less probable cause of HCM, being more appropriate to perform HCM gene panel than targeted FD screening. CONCLUSION: FD prevalence in HCM was 0.9%. Bifascicular block and basal inferolateral LGE were the most powerful predictors of FD in HCM. In their absence, HCM gene panel is the most appropriate step in etiological study of HCM.

Founder effect of Fabry disease due to p.F113L mutation: Clinical profile of a late-onset phenotype.

BACKGROUND: Knowledge on clinical profiles of late-onset phenotypes of Fabry disease (FD) is essential to better define their natural history. Our study aims to demonstrate a founder effect of FD due to the GLA gene mutation c.337T>C (p.F113L) in the Portuguese region of Guimarães; and to characterize the clinical profile of this late-onset phenotype in a large cohort of genetically related adult patients, living in the same region. METHODS AND RESULTS: FD screening was performed in 150 adult patients with hypertrophic cardiomyopathy (HCM) and found 25 Fabry patients (16.6%). The p.F113L mutation was found in 21 of them, leading to a genealogy study and haplotype analysis of the p.F113L patients. Genealogy research revealed a 12-generation family tree with a common ancestor to p.F113L patients, suggesting a founder effect that was supported by haplotype findings. Pedigree analysis was performed and 120 consecutive p.F113L patients underwent a predefined diagnostic evaluation of FD multiorgan involvement. This late-onset phenotype was characterized by common and/or potentially severe cardiac manifestations (left ventricular hypertrophy 40.8%, atrial fibrillation 5%, non-sustained ventricular tachycardia 12.5%, atrioventricular block 18.3%, bifascicular block 13.4%). Extracardiac manifestations included albuminuria>30 mg/24 h 36.1%, chronic kidney disease≥G3 7.6%, brain white matter lesions 54.4%, stroke 3.3%, sensorineural deafness 44.5%, cornea verticillata 13.9%. Plasma lyso-GB3 was undetectable in females, regardless of clinical manifestations. CONCLUSION: A founder effect of FD due to p.F113L mutation was documented by genealogy and genetics in a Portuguese region. In this late-onset phenotype, although cardiac manifestations carry the highest prognostic impact, extracardiac involvement is common.

Is embryo abortion a post-zygotic barrier to gene flow between Littorina ecotypes?

Genetic incompatibilities contribute to reproductive isolation between many diverging populations, but it is still unclear to what extent they play a role if divergence happens with gene flow. In contact zones between the "Crab" and "Wave" ecotypes of the snail Littorina saxatilis, divergent selection forms strong barriers to gene flow, while the role of post-zygotic barriers due to selection against hybrids remains unclear. High embryo abortion rates in this species could indicate the presence of such barriers. Post-zygotic barriers might include genetic incompatibilities (e.g. Dobzhansky-Muller incompatibilities) but also maladaptation, both expected to be most pronounced in contact zones. In addition, embryo abortion might reflect physiological stress on females and embryos independent of any genetic stress. We examined all embryos of >500 females sampled outside and inside contact zones of three populations in Sweden. Females' clutch size ranged from 0 to 1,011 embryos (mean 130 ± 123), and abortion rates varied between 0% and 100% (mean 12%). We described female genotypes by using a hybrid index based on hundreds of SNPs differentiated between ecotypes with which we characterized female genotypes. We also calculated female SNP heterozygosity and inversion karyotype. Clutch size did not vary with female hybrid index, and abortion rates were only weakly related to hybrid index in two sites but not at all in a third site. No additional variation in abortion rate was explained by female SNP heterozygosity, but increased female inversion heterozygosity added slightly to increased abortion. Our results show only weak and probably biologically insignificant post-zygotic barriers contributing to ecotype divergence, and the high and variable abortion rates were marginally, if at all, explained by hybrid index of females.

Signatures of Selection on Standing Genetic Variation Underlie Athletic and Navigational Performance in Racing Pigeons.

Racing pigeons have been selectively bred to find their way home quickly over what are often extremely long distances. This breed is of substantial commercial value and is also an excellent avian model to gain empirical insights into the evolution of traits associated with flying performance and spatial orientation. Here, we investigate the molecular basis of the superior athletic and navigational capabilities of racing pigeons using whole-genome and RNA sequencing data. We inferred multiple signatures of positive selection distributed across the genome of racing pigeons. The strongest signature overlapped the CASK gene, a gene implicated in the formation of neuromuscular junctions. However, no diagnostic alleles were found between racing pigeons and other breeds, and only a small proportion of highly differentiated variants were exclusively detected in racing pigeons. We can thus conclude that very few individual genetic changes, if any, are either strictly necessary or sufficient for superior athletics and navigation. Gene expression analysis between racing and nonracing breeds revealed modest differences in muscle (213) and brain (29). These transcripts, however, showed only slightly elevated levels of genetic differentiation between the two groups, suggesting that most differential expression is not causative but likely a consequence of alterations in regulatory networks. Our results show that the unique suite of traits that enable fast flight, long endurance, and accurate navigation in racing pigeons, do not result from few loci acting as master switches but likely from a polygenic architecture that leveraged standing genetic variation available at the onset of the breed formation.

Multiple chromosomal rearrangements in a hybrid zone between Littorina saxatilis ecotypes.

Both classical and recent studies suggest that chromosomal inversion polymorphisms are important in adaptation and speciation. However, biases in discovery and reporting of inversions make it difficult to assess their prevalence and biological importance. Here, we use an approach based on linkage disequilibrium among markers genotyped for samples collected across a transect between contrasting habitats to detect chromosomal rearrangements de novo. We report 17 polymorphic rearrangements in a single locality for the coastal marine snail, Littorina saxatilis. Patterns of diversity in the field and of recombination in controlled crosses provide strong evidence that at least the majority of these rearrangements are inversions. Most show clinal changes in frequency between habitats, suggestive of divergent selection, but only one appears to be fixed for different arrangements in the two habitats. Consistent with widespread evidence for balancing selection on inversion polymorphisms, we argue that a combination of heterosis and divergent selection can explain the observed patterns and should be considered in other systems spanning environmental gradients.

Comparative and demographic analysis of orang-utan genomes.

'Orang-utan' is derived from a Malay term meaning 'man of the forest' and aptly describes the southeast Asian great apes native to Sumatra and Borneo. The orang-utan species, Pongo abelii (Sumatran) and Pongo pygmaeus (Bornean), are the most phylogenetically distant great apes from humans, thereby providing an informative perspective on hominid evolution. Here we present a Sumatran orang-utan draft genome assembly and short read sequence data from five Sumatran and five Bornean orang-utan genomes. Our analyses reveal that, compared to other primates, the orang-utan genome has many unique features. Structural evolution of the orang-utan genome has proceeded much more slowly than other great apes, evidenced by fewer rearrangements, less segmental duplication, a lower rate of gene family turnover and surprisingly quiescent Alu repeats, which have played a major role in restructuring other primate genomes. We also describe a primate polymorphic neocentromere, found in both Pongo species, emphasizing the gradual evolution of orang-utan genome structure. Orang-utans have extremely low energy usage for a eutherian mammal, far lower than their hominid relatives. Adding their genome to the repertoire of sequenced primates illuminates new signals of positive selection in several pathways including glycolipid metabolism. From the population perspective, both Pongo species are deeply diverse; however, Sumatran individuals possess greater diversity than their Bornean counterparts, and more species-specific variation. Our estimate of Bornean/Sumatran speciation time, 400,000 years ago, is more recent than most previous studies and underscores the complexity of the orang-utan speciation process. Despite a smaller modern census population size, the Sumatran effective population size (N(e)) expanded exponentially relative to the ancestral N(e) after the split, while Bornean N(e) declined over the same period. Overall, the resources and analyses presented here offer new opportunities in evolutionary genomics, insights into hominid biology, and an extensive database of variation for conservation efforts.

A prominent role of KRAB-ZNF transcription factors in mammalian speciation?

The mechanisms of speciation have been one of the most debated topics in evolutionary biology. Among all reproductive barriers, postzygotic reproductive isolation is perhaps the one that has attracted the most attention from geneticists. Despite remarkable advances in the identification of loci involved in Drosophila speciation, little is known about the genes, functions, and biochemical interactions of the molecules underlying hybrid sterility and inviability in mammals. Here, we discuss the main evolutionary and molecular features that make transcription factors (TFs), especially the family of zinc finger proteins with a Krüppel-associated box domain (KRAB-ZNF), strong candidates to play an important role in postzygotic reproductive isolation. Motivated by the recent identification of the gene encoding PR domain zinc finger protein 9 (Prdm9; a KRAB-ZNF gene) as the first hybrid sterility gene identified in mammals, we further propose integrative approaches to study KRAB-ZNF genes with the main goal of characterizing the molecular pathways and interactions involved in hybrid incompatibilities.

Magadi tilapia ecological specialization: filling the early gap in the speciation continuum.

Cichlid fish are well known for their high speciation rates, which are usually accompanied by spectacular and rapid diversification in eco-morphological and secondary sexual traits. This is best illustrated by the famous repeated explosive radiations in the African Great Lakes Tanganyika, Malawi and Victoria, each lake harbouring several hundreds of mostly endemic species. Correspondingly, cichlids diversified very rapidly in many other lakes across their range. Although the larger radiations, unparalleled in vertebrates, are certainly the most intriguing, they are also the most intricate and difficult to address because of their complex nature. This is where smaller, simpler systems may prove to be the most useful. In this issue of Molecular Ecology, Kavembe et al. (2016) report very recent genetic diversification accompanied by ecological specialization in cichlids of the small and ecologically extreme Lake Magadi, in Kenya. Combining geometric morphometrics, stable isotope analysis, population genomics using RADSeq data and coalescent-based modelling techniques, the authors characterize the eco-morphological differences between genetically distinct populations of Magadi tilapia (Alcolapia grahami), which are consistent with the different environmental conditions they experience, and infer their history of divergence. The simplicity of the focal system and the use of a multidisciplinary approach make this work particularly important for our understanding of the early stages of speciation, in both cichlids and other organisms.

Pool and conquer: new tricks for (c)old problems.

The early period of genetics is closely associated with the study of chromosomal inversions. For almost a century, evolutionary biologists found evidence for the role of inversions in central processes such as adaptation and speciation. In spite of huge efforts, many questions remain about the evolutionary forces underlying the distribution and dynamics of inversions in natural populations. Fortunately, old problems can be solved with new tools. In this issue, Kapun et al. (2014) present a remarkable combination of resources and techniques, including publicly available data, karyotyping, statistical estimation of haplotypes, Pool-Seq data and experimental evolution, setting the ground for exciting developments in the field.