Levels of genetic polymorphism: marker loci versus quantitative traits

Species are the units used to measure ecological diversity and alleles are the units of genetic diversity. Genetic variation within and among species has been documented most extensively using allozyme electrophoresis. This reveals wide differences in genetic variability within, and genetic distances among, species, demonstrating that species are not equivalent units of diversity. The extent to which the pattern observed for allozymes can be used to infer patterns of genetic variation in quantitative traits depends on the forces generating and maintaining variability. Allozyme variation is probably not strictly neutral but, nevertheless, heterozygosity is expected to be influenced by population size and genetic distance will be affected by time since divergence. The same is true for quantitative traits influenced by many genes and under weak stabilizing selection. However, the limited data available suggest that allozyme variability is a poor predictor of genetic variation in quantitative traits within populations. It is a better predictor of general phenotypic divergence and of postzygotic isolation between populations or species, but is only weakly correlated with prezygotic isolation. Studies of grasshopper and planthopper mating signal variation and assortative mating illustrate how these characters evolve independently of general genetic and morphological variation. The role of such traits in prezygotic isolation, and hence speciation, means that they will contribute significantly to the diversity of levels of genetic variation within and among species

A Linkage Map and QTL Analysis for Pyrethroid Resistance in the Bed Bug Cimex lectularius.

The rapid evolution of insecticide resistance remains one of the biggest challenges in the control of medically and economically important pests. Insects have evolved a diverse range of mechanisms to reduce the efficacy of the commonly used classes of insecticides and finding the genetic basis of resistance is a major aid to management. In a previously unstudied population, we performed an F2 resistance mapping cross for the common bed bug, Cimex lectularius, for which insecticide resistance is increasingly widespread. Using 334 SNP markers obtained through RAD-sequencing, we constructed the first linkage map for the species, consisting of 14 putative linkage groups (LG), with a length of 407 cM and an average marker spacing of 1.3 cM. The linkage map was used to reassemble the recently published reference genome, facilitating refinement and validation of the current genome assembly. We detected a major QTL on LG12 associated with insecticide resistance, occurring in close proximity (1.2 Mb) to a carboxylesterase encoding candidate gene for pyrethroid resistance. This provides another example of this candidate gene playing a major role in determining survival in a bed bug population following pesticide resistance evolution. The recent availability of the bed bug genome, complete with a full list of potential candidate genes related to insecticide resistance, in addition to the linkage map generated here, provides an excellent resource for future research on the development and spread of insecticide resistance in this resurging pest species

A novel method for predicting the response variability of friction-damped gas turbine blades

Predicting the response of gas turbine blades with underplatform friction dampers is challenging due to the combination of frictional nonlinearity and system uncertainty: a traditional Monte Carlo approach to predicting response distributions requires a large number of nonlinear simulations which is computationally expensive. This paper presents a new approach based on the principle of Maximum Entropy that provides an estimate of the response distribution that is approximately two orders of magnitude faster than Monte Carlo Harmonic Balance Method simulations. The premise is to include the concept of `computational uncertainty': incorporating lack of knowledge of the solution as part of the uncertainty, on the basis that there are diminishing returns in computing precise solutions to an uncertain system. To achieve this, the method uses a describing function approximation of the friction-damped part of the system; chooses an ignorance prior probability density function for the complex value of the describing function based on Coulombs friction law; updates the distribution using an estimate of the mean solution, the admissible domain of solutions, and the principle of Maximum Entropy; then carries out a linear Monte Carlo simulation to estimate the response distribution. The approach is validated by comparison with HBM simulations and experimental tests, using an idealised academic system consisting of a periodic array of beams (with controllable uncertainty) coupled by single-point friction dampers. Comparisons with two- and eight-blade systems show generally good agreement. Predicting the response statistics of the maximum blade amplitude reveals specific well-understood circumstances when the method is less effective. Predictions of the overall blade response statistics agree with Monte Carlo HBM extremely well across a wide range of excitation amplitudes and uncertainty levels. Critically, experimental comparisons reveal the care that is needed in accurately characterising uncertainty in order to obtain agreement of response percentiles. The new method allowed fast iteration of uncertainty parameters and correlations to achieve good agreement, which would not have been possible using traditional methods.Mitsubishi Heavy Industrie

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

Predicting response bounds for friction-damped gas turbine blades with uncertain friction coupling

Friction dampers are often used to reduce high amplitude vibration within gas turbines: they are a robust solution that are able to withstand extreme operating environments. Although the turbine blades are manufactured to tight tolerances, there can be significant variability in the overall response of the assembly. Uncertainties associated with the frictional contact properties are a major factor contributing to this variability. This paper applies a recently developed method for predicting response bounds to friction-damped gas turbines when the characteristics of the friction dampers are unknown, including uncertainty regarding the underlying functional form of the friction law. The approach taken is to represent the frictional contact using a describing function, and formulate an optimisation problem to seek upper and lower bounds on a chosen response metric, such as displacement amplitude. Constraints are chosen that describe known properties of the frictional nonlinearity, without needing to specify a particular constitutive law. The method was validated by comparison with numerical and experimental results from an idealised test system. The experimental test rig consisted of an array of eight beams coupled by pin-contact friction dampers. A modal description of this test rig formed the basis of a numerical model, which uses the Harmonic Balance Method (HBM) for nonlinear simulations. A set of Monte Carlo tests was carried out numerically and experimentally for both a two-beam sub-assembly as well as for the full eight-beam assembly. Comparisons with numerical results showed excellent agreement providing confident verification of the implementation, and comparisons with experimental results revealed that the bounds became less conservative as the system complexity increased. Overall the results are promising: upper and lower response bounds for an array of friction-damped systems can be computed at similar cost to a single HBM simulation, giving reliable bounds that are valid for both parametric and model uncertainties associated with the friction couplings.Mitsubishi Heavy Industrie

Phylogenomics of the adaptive radiation of Triturus newts supports gradual ecological niche expansion towards an incrementally aquatic lifestyle

Newts of the genus Triturus (marbled and crested newts) exhibit substantial variation in the number of trunk vertebrae (NTV) and a higher NTV corresponds to a longer annual aquatic period. Because the Triturus phylogeny has thwarted resolution to date, the evolutionary history of NTV, annual aquatic period, and their potential coevolution has remained unclear. To resolve the phylogeny of Triturus, we generated a c. 6,000 transcriptome-derived marker data set using a custom target enrichment probe set, and conducted phylogenetic analyses using: 1) data concatenation with RAxML, 2) gene-tree summary with ASTRAL, and 3) species-tree estimation with SNAPP. All analyses produce the same, highly supported topology, despite cladogenesis having occurred over a short timeframe, resulting in short internal branch lengths. Our new phylogenetic hypothesis is consistent with the minimal number of inferred changes in NTV count necessary to explain the diversity in NTV observed today. Although a causal relationship between NTV, body form, and aquatic ecology has yet to be experimentally established, our phylogeny indicates that these features have evolved together, and suggest that they may underlie the adaptive radiation that characterizes Triturus

A developmentally descriptive method for quantifying shape in gastropod shells

The growth of snail shells can be described by simple mathematical rules. Variation in a few parameters can explain much of the diversity of shell shapes seen in nature. However, empirical studies of gastropod shell shape variation typically use geometric morphometric approaches, which do not capture this growth pattern. We have developed a way to infer a set of developmentally descriptive shape parameters based on three-dimensional logarithmic helicospiral growth and using landmarks from two-dimensional shell images as input. We demonstrate the utility of this approach, and compare it to the geometric morphometric approach, using a large set of Littorina saxatilis shells in which locally adapted populations differ in shape. Our method can be modified easily to make it applicable to a wide range of shell forms, which would allow for investigations of the similarities and differences between and within many different species of gastropods

Hybridization patterns between two marine snails, Littorina fabalis and L. obtusata

Characterizing the patterns of hybridization between closely related species is crucial to understand the role of gene flow in speciation. In particular, systems comprising multiple contacts between sister species offer an outstanding opportunity to investigate how reproductive isolation varies with environmental conditions, demography and geographic contexts of divergence. The flat periwinkles, Littorina obtusata and L. fabalis (Gastropoda), are two intertidal sister species with marked ecological differences compatible with late stages of speciation. Although hybridization between the two was previously suggested, its extent across the Atlantic shores of Europe remained largely unknown. Here, we combined genetic (microsatellites and mtDNA) and morphological data (shell and male genital morphology) from multiple populations of flat periwinkles in north‐western Iberia to assess the extent of current and past hybridization between L. obtusata and L. fabalis under two contrasting geographic settings of divergence (sympatry and allopatry). Hybridization signatures based on both mtDNA and microsatellites were stronger in sympatric sites, although evidence for recent extensive admixture was found in a single location. Misidentification of individuals into species based on shell morphology was higher in sympatric than in allopatric sites. However, despite hybridization, species distinctiveness based on this phenotypic trait together with male genital morphology remained relatively high. The observed variation in the extent of hybridization among locations provides a rare opportunity for future studies on the consequences of different levels of gene flow for reinforcement, thus informing about the mechanisms underlying the completion of speciation

Genomic architecture of parallel ecological divergence : beyond a single environmental contrast

The study of parallel ecological divergence provides important clues to the operation of natural selection. Parallel divergence often occurs in heterogeneous environments with different kinds of environmental gradients in different locations, but the genomic basis underlying this process is unknown. We investigated the genomics of rapid parallel adaptation in the marine snail Littorina saxatilis in response to two independent environmental axes (crab-predation versus wave-action and low-shore versus high-shore). Using pooled whole-genome resequencing, we show that sharing of genomic regions of high differentiation between environments is generally low but increases at smaller spatial scales. We identify different shared genomic regions of divergence for each environmental axis and show that most of these regions overlap with candidate chromosomal inversions. Several inversion regions are divergent and polymorphic across many localities. We argue that chromosomal inversions could store shared variation that fuels rapid parallel adaptation to heterogeneous environments, possibly as balanced polymorphism shared by adaptive gene flow

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