Ecological speciation in postglacial European whitefish: rapid adaptive radiations into the littoral,pelagic, and profundal lake habitats

Understanding how a monophyletic lineage of a species diverges into several adaptive forms has received increased attention in recent years, but the underlying mechanisms in this process are still under debate. Postglacial fishes are excellent model organisms for exploring this process, especially the initial stages of ecological speciation, as postglacial lakes represent replicated discrete environments with variation in available niches. Here, we combine data of niche utilization, trophic morphology, and 17 microsatellite loci to investigate the diversification process of three sympatric European whitefish morphs from three northern Fennoscandian lakes. The morphological divergence in the gill raker number among the whitefish morphs was related to the utilization of different trophic niches and was associated with reproductive isolation within and across lakes. The intralacustrine comparison of whitefish morphs showed that these systems represent two levels of adaptive divergence: (1) a consistent littoral–pelagic resource axis; and (2) a more variable littoral–profundal resource axis. The results also indicate that the profundal whitefish morph has diverged repeatedly from the ancestral littoral whitefish morph in sympatry in two different watercourses. In contrast, all the analyses performed revealed clustering of the pelagic whitefish morphs across lakes suggesting parallel postglacial immigration with the littoral whitefish morph into each lake. Finally, the analyses strongly suggested that the trophic adaptive trait, number of gill rakers, was under diversifying selection in the different whitefish morphs. Together, the results support a complex evolutionary scenario where ecological speciation acts, but where both allopatric (colonization history) and sympatric (within watercourse divergence) processes are involved.     

Generation of a neutral FST baseline for testing local adaptation on gill raker number within and between European whitefish ecotypes in the Baltic Sea basin

Divergent selection at ecologically important traits is thought to be a major factor driving phenotypic differentiation between populations. To elucidate the role of different evolutionary processes shaping the variation in gill raker number of European whitefish (Coregonus lavaretus sensu lato) in the Baltic Sea basin, we assessed the relationships between genetic and phenotypic variation among and within three whitefish ecotypes (sea spawners, river spawners and lake spawners). To generate expected neutral distribution of F_ST and to evaluate whether highly variable microsatellite loci resulted in deflated F_ST estimates compared to less variable markers, we performed population genetic simulations under finite island and hierarchical island models. The genetic divergence observed among (F_CT = 0.010) and within (F_ST = 0.014–0.041) ecotypes was rather low. The divergence in gill raker number, however, was substantially higher between sea and river spawners compared to observed microsatellite data and simulated neutral baseline (P_CT > F_CT). This suggests that the differences in gill raker number between sea and river spawners are likely driven by divergent natural selection. We also found strong support for divergent selection on gill raker number among different populations of sea spawners (P_ST > F_ST), most likely caused by highly variable habitat use and diverse diet. The putative role of divergent selection within lake spawners initially inferred from empirical microsatellite data was not supported by simulated F_ST distributions. This work provides a first formal test of divergent selection on gill raker number in Baltic whitefish, and demonstrates the usefulness of population genetic simulations to generate informative neutral baselines for P_ST–F_ST analyses helping to disentangle the effects of stochastic evolutionary processes from natural selection.     

Variable extent of parallelism in respiratory,circulatory, and neurological traits across lake whitefish species pairs

Parallel adaptive radiation events provide a powerful framework for investigations of ecology's contribution to phenotypic diversification. Ecologically driven divergence has been invoked to explain the repeated evolution of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis) species in multiple lakes in eastern North America. Nevertheless, links between most putatively adaptive traits and ecological variation remain poorly defined within and among whitefish species pairs. Here, we examine four species pairs for variation in gill, heart, and brain size; three traits predicted to show strong phenotypic responses to ecological divergence. In each of the species pairs, normals exhibited larger body size standardized gills compared to dwarfs – a pattern that is suggestive of a common ecological driver of gill size divergence. Within lakes, the seasonal hypoxia experienced in the benthic environment is a likely factor leading to the requirement for larger gills in normals. Interestingly, the morphological pathways used to achieve larger gills varied between species pairs from Québec and Maine, which may imply subtle non‐parallelism in gill size divergence related to differences in genetic background. There was also a non‐significant trend toward larger hearts in dwarfs, the more active species of the two, whereas brain size varied exclusively among the lake populations. Taken together, our results suggest that the diversification of whitefish has been driven by parallel and non‐parallel ecological conditions across lakes. Furthermore, the phenotypic response to ecological variation may depend on genetic background of each population.     

Analysis of ROH patterns in the Noriker horse breed reveals signatures of selection for coat color and body size

Overlapping runs of homozygosity (ROH islands) shared by the majority of a population are hypothesized to be the result of selection around a target locus. In this study we investigated the impact of selection for coat color within the Noriker horse on autozygosity and ROH patterns. We analyzed overlapping homozygous regions (ROH islands) for gene content in fragments shared by more than 50% of horses. Long‐term assortative mating of chestnut horses and the small effective population size of leopard spotted and tobiano horses resulted in higher mean genome‐wide ROH coverage (S_ROH) within the range of 237.4–284.2 Mb, whereas for bay, black and roan horses, where rotation mating is commonly applied, lower autozygosity (S_ROH from 176.4–180.0 Mb) was determined. We identified seven common ROH islands considering all Noriker horses from our dataset. Specific islands were documented for chestnut, leopard spotted, roan and bay horses. The ROH islands contained, among others, genes associated with body size (ZFAT, LASP1 and LCORL/NCAPG), coat color (MC1R in chestnut and the factor PATN1 in leopard spotted horses) and morphogenesis (HOXB cluster in all color strains except leopard spotted horses). This study demonstrates that within a closed population sharing the same founders and ancestors, selection on a single phenotypic trait, in this case coat color, can result in genetic fragmentation affecting levels of autozygosity and distribution of ROH islands and enclosed gene content.     

Adaptive divergence in body size overrides the effects of plasticity across natural habitats in the brown trout

The evolution of life-history traits is characterized by trade-offs between different selection pressures, as well as plasticity across environmental conditions. Yet, studies on local adaptation are often performed under artificial conditions, leaving two issues unexplored: (i) how consistent are laboratory inferred local adaptations under natural conditions and (ii) how much phenotypic variation is attributed to phenotypic plasticity and to adaptive evolution, respectively, across environmental conditions? We reared fish from six locally adapted (domesticated and wild) populations of anadromous brown trout (Salmo trutta) in one semi-natural and three natural streams and recorded a key life-history trait (body size at the end of first growth season). We found that population-specific reaction norms were close to parallel across different streams and Q_ST was similar – and larger than F_ST – within all streams, indicating a consistency of local adaptation in body size across natural environments. The amount of variation explained by population origin exceeded the variation across stream environments, indicating that genetic effects derived from adaptive processes have a stronger effect on phenotypic variation than plasticity induced by environmental conditions. These results suggest that plasticity does not “swamp” the phenotypic variation, and that selection may thus be efficient in generating genetic change.     

Seasonal variations in body melanism and size of the wolf spider Pardosa astrigera (Araneae: Lycosidae)

Variations in species morphology and life‐history traits strongly correlate with geographic and climatic characteristics. Most studies on morphological variations in animals focus on ectotherms distributed on a large geographic scale across latitudinal and/or altitudinal gradient. However, the morphological variations of spiders living in the same habitats across different seasons have not been reported. In this study, we used the wolf spider, Pardosa astrigera, as a model to determine seasonal differences in adult body size, melanism, fecundity, and egg diameter both in the overwintering and the first generation for 2010 and 2016. The results showed that in 2010, both females and males of the overwintering generation were significantly darker than the first generation. Moreover, the overwintering females were markedly larger and produced more and bigger eggs than the first generation in both 2010 and 2016. Considering the overwintering P. astrigera experiencing low temperature and/or desiccation stress, these results suggest that substantially darker and larger body of the overwintering generation is adaptive to adverse conditions.     

Directional selection in the evolution of elongated upper canines in clouded leopards and sabre‐toothed cats

Extremely developed or specialized traits such as the elongated upper canines of extinct sabre‐toothed cats are often not analogous to those of any extant species, which limits our understanding of their evolutionary cause. However, an extant species may have undergone directional selection for a similar extreme phenotype. Among living felids, the clouded leopard, Neofelis nebulosa, has exceptionally long upper canines for its body size. We hypothesized that directional selection generated the elongated upper canines of clouded leopards in a manner similar to the process in extinct sabre‐toothed cats. To test this, we developed an approach that compared the effect of directional selection among lineages in a phylogeny using a simulation of trait evolution and approximate Bayesian computation. This approach was applied to analyse the evolution of upper canine length in the Felidae phylogeny. Our analyses consistently showed directional selection favouring longer upper canines in the clouded leopard lineage and a lineage leading to the sabre‐toothed cat with the longest upper canines, Smilodon. Most of our analyses detected an effect of directional selection for longer upper canines in the lineage leading to another sabre‐toothed cat, Homotherium, although this selection may have occurred exclusively in the primitive species. In all the analyses, the clouded leopard and Smilodon lineages showed comparable directional selection. This implies that clouded leopards share a selection advantage with sabre‐toothed cats in having elongated upper canines.     

Rapid evolution and the genomic consequences of selection against interspecific mating

While few species introduced into a new environment become invasive, those that do provide critical information on ecological mechanisms that determine invasions success and the evolutionary responses that follow invasion. Aedes albopictus (the Asian tiger mosquito) was introduced into the naturalized range of Aedes aegypti (the yellow fever mosquito) in the United States in the mid‐1980s, resulting in the displacement of A. aegypti in much of the south‐eastern United States. The rapid displacement was likely due to the superior competitive ability of A. albopictus as larvae and asymmetric mating interference competition, in which male A. albopictus mate with and sterilize A. aegypti females, a process called “satyrization.” The goal of this study was to examine the genomic responses of a resident species to an invasive species in which the mechanism of character displacement is understood. We used double‐digest restriction enzyme DNA sequencing (ddRADseq) to analyse outlier loci between selected and control lines of laboratory‐reared A. aegypti females from two populations (Tucson, AZ and Key West, Florida, USA), and individual females classified as either “resisted” or “mated with” A. albopictus males via mating trials of wild‐derived females from four populations in Florida. We found significant outlier loci in comparing selected and control lines and between mated and nonmated A. aegypti females in the laboratory and wild‐derived populations, respectively. We found overlap in specific outlier loci between different source populations that support consistent genomic signatures of selection within A. aegypti. Our results point to regions of the A. aegypti genome and potential candidate genes that may be involved in mating behaviour, and specifically in avoiding interspecific mating choices.     

Low reproductive isolation and highly variable levels of gene flow reveal limited progress towards speciation between European river and brook lampreys

Ecologically based divergent selection is a factor that could drive reproductive isolation even in the presence of gene flow. Population pairs arrayed along a continuum of divergence provide a good opportunity to address this issue. Here, we used a combination of mating trials, experimental crosses and population genetic analyses to investigate the evolution of reproductive isolation between two closely related species of lampreys with distinct life histories. We used microsatellite markers to genotype over 1000 individuals of the migratory parasitic river lamprey (Lampetra fluviatilis) and freshwater‐resident nonparasitic brook lamprey (Lampetra planeri) distributed in 10 sympatric and parapatric population pairs in France. Mating trials, parentage analyses and artificial fertilizations demonstrated a low level of reproductive isolation between species even though size‐assortative mating may contribute to isolation. Most parapatric population pairs were strongly differentiated due to the joint effects of geographic distance and barriers to migration. In contrast, we found variable levels of gene flow between sympatric populations ranging from panmixia to moderate differentiation, which indicates a gradient of divergence with some population pairs that may correspond to alternative morphs or ecotypes of a single species and others that remain partially isolated. Ecologically based divergent selection may explain these variable levels of divergence among sympatric population pairs, but incomplete genome swamping following secondary contact could have also played a role. Overall, this study illustrates how highly differentiated phenotypes can be maintained despite high levels of gene flow that limit the progress towards speciation.     

Parallel effects of the inversion In(3R)Payne on body size across the North American and Australian clines in Drosophila melanogaster

Chromosomal inversions are thought to play a major role in climatic adaptation. In D. melanogaster, the cosmopolitan inversion In(3R)Payne exhibits latitudinal clines on multiple continents. As many fitness traits show similar clines, it is tempting to hypothesize that In(3R)P underlies observed clinal patterns for some of these traits. In support of this idea, previous work in Australian populations has demonstrated that In(3R)P affects body size but not development time or cold resistance. However, similar data from other clines of this inversion are largely lacking; finding parallel effects of In(3R)P across multiple clines would considerably strengthen the case for clinal selection. Here, we have analysed the phenotypic effects of In(3R)P in populations originating from the endpoints of the latitudinal cline along the North American east coast. We measured development time, egg‐to‐adult survival, several size‐related traits (femur and tibia length, wing area and shape), chill coma recovery, oxidative stress resistance and triglyceride content in homokaryon lines carrying In(3R)P or the standard arrangement. Our central finding is that the effects of In(3R)P along the North American cline match those observed in Australia: standard arrangement lines were larger than inverted lines, but the inversion did not influence development time or cold resistance. Similarly, In(3R)P did not affect egg‐to‐adult survival, oxidative stress resistance and lipid content. In(3R)P thus seems to specifically affect size traits in populations from both continents. This parallelism strongly suggests an adaptive pattern, whereby the inversion has captured alleles associated with growth regulation and clinal selection acts on size across both continents.     

Ecological drivers of body size evolution and sexual size dimorphism in short‐horned grasshoppers (Orthoptera: Acrididae)

Sexual size dimorphism (SSD) is widespread and variable in nature. Although female‐biased SSD predominates among insects, the proximate ecological and evolutionary factors promoting this phenomenon remain largely unstudied. Here, we employ modern phylogenetic comparative methods on eight subfamilies of Iberian grasshoppers (85 species) to examine the validity of different models of evolution of body size and SSD and explore how they are shaped by a suite of ecological variables (habitat specialization, substrate use, altitude) and/or constrained by different evolutionary pressures (female fecundity, strength of sexual selection, length of the breeding season). Body size disparity primarily accumulated late in the history of the group and did not follow a Brownian motion pattern, indicating the existence of directional evolution for this trait. We found support for the converse of Rensch's rule (i.e. females are proportionally bigger than males in large species) across all taxa but not within the two most speciose subfamilies (Gomphocerinae and Oedipodinae), which showed an isometric pattern. Our results do not provide support for the fecundity or sexual selection hypotheses, and we did not find evidence for significant effects of habitat use. Contrary to that expected, we found that species with narrower reproductive window are less dimorphic in size than those that exhibit a longer breeding cycle, suggesting that male protandry cannot solely account for the evolution of female‐biased SSD in Orthoptera. Our study highlights the need to consider alternatives to the classical evolutionary hypotheses when trying to explain why in certain insect groups males remain small.     

The interplay between natural and sexual selection in the evolution of sexual size dimorphism in Sceloporus lizards (Squamata: Phrynosomatidae)

Sexual size dimorphism (SSD) evolves because body size is usually related to reproductive success through different pathways in females and males. Female body size is strongly correlated with fecundity, while in males, body size is correlated with mating success. In many lizard species, males are larger than females, whereas in others, females are the larger sex, suggesting that selection on fecundity has been stronger than sexual selection on males. As placental development or egg retention requires more space within the abdominal cavity, it has been suggested that females of viviparous lizards have larger abdomens or body size than their oviparous relatives. Thus, it would be expected that females of viviparous species attain larger sizes than their oviparous relatives, generating more biased patterns of SSD. We test these predictions using lizards of the genus Sceloporus. After controlling for phylogenetic effects, our results confirm a strong relationship between female body size and fecundity, suggesting that selection for higher fecundity has had a main role in the evolution of female body size. However, oviparous and viviparous females exhibit similar sizes and allometric relationships. Even though there is a strong effect of body size on female fecundity, once phylogenetic effects are considered, we find that the slope of male on female body size is significantly larger than one, providing evidence of greater evolutionary divergence of male body size. These results suggest that the relative impact of sexual selection acting on males has been stronger than fecundity selection acting on females within Sceloporus lizards.     

Cunningham's skinks show low genetic connectivity and signatures of divergent selection across its distribution

Establishing corridors of connecting habitat has become a mainstay conservation strategy to maintain gene flow and facilitate climate‐driven range shifts. Yet, little attention has been given to ascertaining the extent to which corridors will benefit philopatric species, which might exhibit localized adaptation. Measures of genetic connectivity and adaptive genetic variation across species’ ranges can help fill this knowledge gap. Here, we characterized the spatial genetic structure of Cunningham's skink (Egernia cunninghami), a philopatric species distributed along Australia's Great Dividing Range, and assessed evidence of localized adaptation. Analysis of 4,274 SNPs from 94 individuals sampled at four localities spanning 500 km and 4° of latitude revealed strong genetic structuring at neutral loci (mean F_ST ± SD = 0.603 ± 0.237) among the localities. Putatively neutral SNPs and those under divergent selection yielded contrasting spatial patterns, with the latter identifying two genetically distinct clusters. Given low genetic connectivity of the four localities, we suggest that the natural movement rate of this species is insufficient to keep pace with spatial shifts to its climate envelope, irrespective of habitat availability. In addition, our finding of localized adaptation highlights the risk of outbreeding depression should the translocation of individuals be adopted as a conservation management strategy.     

Maternal control of seed weight in rapeseed (Brassica napus L.): the causal link between the size of pod (mother,source) and seed (offspring,sink)

Seed size/weight is one of the key traits related to plant domestication and crop improvement. In rapeseed (Brassica napus L.) germplasm, seed weight shows extensive variation, but its regulatory mechanism is poorly understood. To identify the key mechanism of seed weight regulation, a systematic comparative study was performed. Genetic, morphological and cytological evidence showed that seed weight was controlled by maternal genotype, through the regulation of seed size mainly via cell number. The physiological evidence indicated that differences in the pod length might result in differences in pod wall photosynthetic area, carbohydrates and the final seed weight. We also identified two pleiotropic major quantitative trait loci that acted indirectly on seed weight via their effects on pod length. RNA‐seq results showed that genes related to pod development and hormones were significantly differentially expressed in the pod wall; genes related to development, cell division, nutrient reservoir and ribosomal proteins were all up‐regulated in the seeds of the large‐seed pool. Finally, we proposed a potential seed weight regulatory mechanism that is specific to rapeseed and novel in plants. The results demonstrate a causal link between the size of the pod (mother, source) and the seed (offspring, sink) in rapeseed, which provides novel insight into the maternal control of seed weight and will open a new research field in plants.     

The role of geography and ecology in shaping repeated patterns of morphological and genetic differentiation between European minnows (Phoxinus phoxinus) from the Pyrenees and the Alps

Neutral and selective processes can drive repeated patterns of evolution in different groups of populations experiencing similar ecological gradients. In this paper, we used a combination of nuclear and mitochondrial DNA markers, as well as geometric morphometrics, to investigate repeated patterns of morphological and genetic divergence of European minnows in two mountain ranges: the Pyrenees and the Alps. European minnows (Phoxinus phoxinus) are cyprinid fish inhabiting most freshwater bodies in Europe, including those in different mountain ranges that could act as major geographical barriers to gene flow. We explored patterns of P. phoxinus phenotypic and genetic diversification along a gradient of altitude common to the two mountain ranges, and tested for isolation by distance (IBD), isolation by environment (IBE) and isolation by adaptation (IBA). The results indicated that populations from the Pyrenees and the Alps belong to two well differentiated, reciprocally monophyletic mtDNA lineages. Substantial genetic differentiation due to geographical isolation within and between populations from the Pyrenees and the Alps was also found using rapidly evolving AFLPs markers (isolation by distance or IBD), as well as morphological differences between mountain ranges. Also, morphology varied strongly with elevation and so did genetic differentiation to a lower extent. Despite moderate evidence for IBE and IBA, and therefore of repeated evolution, substantial population heterogeneity was found at the genetic level, suggesting that selection and population specific genetic drift act in concert to affect genetic divergence.