Find research datasets worth reusing

Data from: A dedicated target capture approach reveals variable genetic markers across micro- and macro-evolutionary time scales in palms

Understanding the genetics of biological diversification across micro- and macro-evolutionary time scales is a vibrant field of research for molecular ecologists as rapid advances in sequencing technologies promise to overcome former limitations. In palms, an emblematic, economically and ecologically important plant family with high diversity in the tropics, studies of diversification at the population and species levels are still hampered by a lack of genomic markers suitable for the genotyping of large numbers of recently diverged taxa. To fill this gap, we used a whole genome sequencing approach to develop target sequencing for molecular markers in 4,184 genome regions, including 4,051 genes and 133 non-genic putatively neutral regions. These markers were chosen to cover a wide range of evolutionary rates allowing future studies at the family, genus, species and population levels. Special emphasis was given to the avoidance of copy number variation during marker selection. In addition, a set of 149 well-known sequence regions previously used as phylogenetic markers by the palm biological research community were included in the target regions, to open the possibility to combine and jointly analyse already available data sets with genomic data to be produced with this new toolkit. The bait set was effective for species belonging to all three palm subfamilies tested (Arecoideae, Ceroxyloideae and Coryphoideae), with high mapping rates, specificity and efficiency. The number of high quality Single Nucleotide Polymorphisms (SNPs) detected at both the subfamily and population levels facilitates efficient analyses of genomic diversity across micro- and macro-evolutionary time scales.

Data from: Genome-wide evolutionary dynamics of influenza B viruses on a global scale

The global-scale epidemiology and genome-wide evolutionary dynamics of influenza B remain poorly understood compared with influenza A viruses. We compiled a spatio-temporally comprehensive dataset of influenza B viruses, comprising over 2,500 genomes sampled worldwide between 1987 and 2015, including 382 newly-sequenced genomes that fill substantial gaps in previous molecular surveillance studies. Our contributed data increase the number of available influenza B virus genomes in Europe, Africa and Central Asia, improving the global context to study influenza B viruses. We reveal Yamagata-lineage diversity results from co-circulation of two antigenically-distinct groups that also segregate genetically across the entire genome, without evidence of intra-lineage reassortment. In contrast, Victoria-lineage diversity stems from geographic segregation of different genetic clades, with variability in the degree of geographic spread among clades. Differences between the lineages are reflected in their antigenic dynamics, as Yamagata-lineage viruses show alternating dominance between antigenic groups, while Victoria-lineage viruses show antigenic drift of a single lineage. Structural mapping of amino acid substitutions on trunk branches of influenza B gene phylogenies further supports these antigenic differences and highlights two potential mechanisms of adaptation for polymerase activity. Our study provides new insights into the epidemiological and molecular processes shaping influenza B virus evolution globally.

opencc-zeroDec 2011View details →

Data from: Direct genetic evidence for reproductive philopatry and associated fine-scale migrations in female blacktip reef sharks (Carcharhinus melanopterus) in French Polynesia.

Conservation of top predators has been emphasised as essential in an ecosystem due to their role in trophic chain regulation. Optimization of conservation strategies for these endangered marine top predators requires direct estimates of breeding patterns and connectivity, which are essential to understand their population dynamics. While there have been some attempts to investigate breeding patterns of reef sharks from litter reconstruction using molecular analyses, direct fine-scale migrations of female sharks for parturition as well as connectivity at a medium scale, such as between islands, remain mostly unknown. We used microsatellite DNA markers and a likelihood-based parentage analysis to determine breeding patterns of female blacktip reef sharks in Moorea (Society Islands, French Polynesia). Such an approach demonstrates that while most females gave birth at their home island, some migrated to specific nursery areas outside their home range, sometimes going to another island (up to 50 km across deep ocean). Analysis also revealed that females migrated to the same nursery on every birthing event. Many offspring showed a high level of inbreeding indicating an overall reduced population size, restricted movements and dispersal, or specific mating behaviour. Females represent the vectors that transport the genes at nursery grounds, and their fidelity should thus define reproductive units. As females seem to be philopatric, males could be the ones dispersing genes between metapopulations. These results highlight the need for conservation of coastal zones where female reef sharks seems to exhibit philopatry during the breeding season.

Data from: Fine-scale analysis of an assassin bug's behaviour: predatory strategies to bypass the sensory systems of prey

Some predators sidestep environments that render them conspicuous to the sensory systems of prey. However, these challenging environments are unavoidable for certain predators. Stenolemus giraffa is an assassin bug that feeds on web-building spiders; the web is the environment in which this predator finds its prey, but it also forms part of its preys' sophisticated sensory apparatus, blurring the distinction between environment and sensory systems. Stenolemus giraffa needs to break threads in the web that obstruct its path to the spiders, and such vibrations can alert the spiders. Using laser vibrometry, this study demonstrates how S. giraffa avoids alerting the spiders during its approach. When breaking threads, S. giraffa attenuates the vibrations produced by holding on to the loose ends of the broken thread and causing them to sag prior to release. In addition, S. giraffa releases the loose ends of a broken thread one at a time (after several seconds or minutes) and in this way spaces out the production of vibrations in time. Furthermore, S. giraffa was found to maximally reduce the amplitude of vibrations when breaking threads that are prone to produce louder vibrations. Finally, S. giraffa preferred to break threads in the presence of wind, suggesting that this araneophagic insect exploits environmental noise that temporarily impairs the spiders' ability to detect vibrations. The predatory behaviour of S. giraffa seems to be adaptated in intricate manner for bypassing the sophisticated sensory systems of web-building spiders. These findings illustrate how the physical characteristics of the environment, along with the sensory systems of prey can shape the predatory strategies of animals.

opencc-zeroOct 2019View details →

At what spatial scale(s) do mammals respond to urbanization?

Spatial scale is fundamental in understanding species-landscape relationships because species' responses to landscape characteristics typically vary across scales. Nonetheless, such scales are often unidentified or unreliably predicted by theory. Many landscapes worldwide are urbanizing, yet the spatial scaling of species' responses to urbanization is poorly understood. We investigated the spatial scaling of urbanization effects on a community of 15 mammal species using ~ 60,000 wildlife detections collected from a constellation of 207 camera traps across an extensive urban park system. We embedded a bivariate Gaussian kernel in hierarchical multi-species models to determine two scales of effect (a scale of maximal effect and a broader scale of cumulative landscape effect) for two biological responses (occupancy and site visit frequency) across two seasons (winter and summer) for each species. We then assessed whether scales of effect varied according to theoretical predictions associated with biological responses and species traits (body size and mobility). Scales of effect ranged from < 50 m to > 9,000 m and varied among species, but not as predicted by theory. Species' occupancy generally showed a weak response to urbanization and the scale of this effect was both highly uncertain and consistent across species. We did not detect any relationship between scales of effect and species' body size or mobility, nor was there any evident pattern of scaling across biological response or seasons. These results imply that 1) urbanization effects on mammals manifest across a very broad spectrum of spatial scales, and 2) current theories that a priori predict the scale at which urbanization affects mammals may be of limited use within a given system. Overall, this study suggests that developing general theory regarding the scaling of species-landscape relationships requires additional empirical work conducted across multiple species, systems, and timescales.

opencc-zeroDec 2011View details →

Data from: Identifying biases at different spatial and temporal scales of diversification: a case study in the Neotropical parrotlet genus Forpus

The temporal origins of the extraordinary biodiversity of the Neotropical region are highly debated. Recent empirical work has found support for alternative models on the tempo of speciation in Neotropical species further fuelling the debate. However, relationships within many Neotropical lineages are poorly understood and it is unclear how this uncertainty impacts inferences on the evolution of taxa in the region. We examined the robustness of diversification patterns in the avian genus Forpus by testing whether the use of different units of biodiversity (i.e. biological species and statistically inferred species) impacted diversification rates and inferences regarding important biogeographic breaks in the genus. We found that the best-fit model of diversification for the biological species data set was a declining rate of diversification; whereas a model of constant diversification was the best-fit model for statistically inferred species or subspecies. Moreover, the relative importance of different landscape features in delimiting genetic structure across the landscape varied across data sets with differing units of biodiversity. Patterns based on divergence times among biological species indicated old speciation events across major geographic and river barriers, in contrast, data sets more inclusive of the diversity in Forpus illustrate the role of both old divergence across major landscape features as well as more recent divergences that are possibly attributed to Pleistocene climatic changes. Overall, these results indicate that conflicting models on the temporal origins of Neotropical birds may be attributable to sampling biases.

Data from: Multiple scaling behavior and nonlinear traits in music scores

We present a statistical analysis of music scores from different composers using detrended fluctuation analysis. We find different fluctuation profiles that correspond to distinct auto-correlation structures of the musical pieces. Further, we reveal evidence for the presence of nonlinear auto-correlations by estimating the detrended fluctuation analysis of the magnitude series, a result validated by a corresponding study of appropriate surrogate data. The amount and the character of nonlinear correlations vary from one composer to another. Finally, we performed a simple experiment in order to evaluate the pleasantness of the musical surrogate pieces in comparison with the original music and find that nonlinear correlations could play an important role in the aesthetic perception of a musical piece.

Data from: The impact of variable degrees of freedom and scale parameters in Bayesian methods for genomic prediction in Chinese Simmental beef cattle

Three conventional Bayesian approaches (BayesA, BayesB and BayesCπ) have been demonstrated to be powerful in predicting genomic merit for complex traits in livestock. A priori, these Bayesian models assume that the non-zero SNP effects (marginally) follow a t-distribution depending on two fixed hyperparameters, degrees of freedom and scale parameters. In this study, we performed genomic prediction in Chinese Simmental beef cattle and treated degrees of freedom and scale parameters as unknown with inappropriate priors. Furthermore, we compared the modified methods (BayesFA, BayesFB and BayesFCπ) with their corresponding counterparts using simulation datasets. We found that the modified methods with distribution assumed to the two hyperparameters were beneficial for improving the predictive accuracy. Our results showed that the predictive accuracies of the modified methods were slightly higher than those of their counterparts especially for traits with low heritability and a small number of QTLs. Moreover, cross-validation analysis for three traits, namely carcass weight, live weight and tenderloin weight, in 1136 Simmental beef cattle suggested that predictive accuracy of BayesFCπ noticeably outperformed BayesCπ with the highest increase (3.8%) for live weight using the cohort masking cross-validation.

Data from: Lateralized feeding behavior is associated with asymmetrical neuroanatomy and lateralized gene expressions in the brain in scale-eating cichlid fish

Lateralized behavior ('handedness') is unusual, but consistently found across diverse animal lineages, including humans. It is thought to reflect brain anatomical and/or functional asymmetries, but its neuro-molecular mechanisms remain largely unknown. Lake Tanganyika scale-eating cichlid fish, Perissodus microlepis show pronounced asymmetry in their jaw morphology as well as handedness in feeding behavior – feeding scales preferentially only from one or the other side of their victims. This makes them an ideal model in which to investigate potential laterality in neuroanatomy and transcription in the brain in relation to behavioral handedness. After determining behavioral handedness in P. microlepis (preferred attack side), we estimated the volume of the hemispheres of brain regions and captured their gene expression profiles. Our analyses revealed that the degree of behavioral handedness is mirrored at the level of neuroanatomical asymmetry, particularly in the tectum opticum. Transcriptome analyses showed that different brain regions (tectum opticum, telencephalon, hypothalamus and cerebellum) display distinct expression patterns, potentially reflecting their developmental interrelationships. For numerous genes in each brain region, their extent of expression differences between hemispheres was found to be correlated with the degree of behavioral lateralization. Interestingly, the tectum opticum and telencephalon showed divergent biases on the direction of up- or down-regulation of the laterality candidate genes (e.g., grm2) in the hemispheres, highlighting the connection of handedness with gene expression profiles and the different roles of these brain regions. Hence, handedness in predation behavior may be caused by asymmetric size of brain hemispheres and also by lateralized gene expressions in the brain.

Data from: Spatiotemporal scaling of plant species richness and functional diversity in a temperate semi-natural grassland

The accumulation of biodiversity in space and time has been modelled extensively using the species-area relationship and the species-time relationship, respectively. Recently, these models have been combined into time-area curves in order to investigate spatiotemporal scaling of species richness. This study expands on previous research by applying these spatiotemporal models to functional diversity. Understanding spatiotemporal dynamics of ecological traits is important due to their crucial role in ecosystem functioning and mediating species responses to environmental change. We present a new function based on the semi-logarithmic species-area relationship, which was applied with a power function to vegetation survey data from Scottish machair grassland for both species richness and two measures of functional diversity. When taking a whole-study approach using non-linear mixed effects models, the semi-logarithmic function used here shows a positive time-area interaction for species richness, contrasting with the negative interaction of the power law found in previous investigations. Although there was a negative time-area interaction for functional diversity measures at the whole-study scale, parameter estimates were inconsistent at the individual site level. Overall, the results reveal differing spatiotemporal dynamics of species and their traits and suggest that the appropriate scale for space-for-time substitutions depends on the aspect of biodiversity being investigated. The new model developed in this study, and the novel application to functional diversity, opens up future possible research into spatiotemporal dynamics of biodiversity.

Data from: Meta‐analysis of chromosome‐scale crossover rate variation in eukaryotes and its significance to evolutionary genomics

Understanding the distribution of crossovers along chromosomes is crucial to evolutionary genomics because the crossover rate determines how strongly a genome region is influenced by natural selection. Nevertheless, generalities in the chromosome-scale distribution of crossovers have not been investigated formally. We fill this gap by synthesizing joint information on genetic and physical maps across 62 animal, plant, and fungal species. Our quantitative analysis reveals a strong and taxonomically wide-spread reduction of the crossover rate in the center of chromosomes relative to their peripheries. We demonstrate that this pattern is poorly explained by the position of the centromere, but find that the magnitude of the relative reduction in the crossover rate in chromosome centers increases with chromosome length. That is, long chromosomes often display a dramatically low crossover rate in their center whereas short chromosomes exhibit a relatively homogeneous crossover rate. This observation is compatible with a model in which crossovers are initiated from the chromosome tips, an idea with preliminary support from mechanistic investigations of meiotic recombination. Consequently, we show that organisms achieve a higher genome-wide crossover rate by evolving smaller chromosomes. Summarizing theory and providing empirical examples, we finally highlight that taxonomically wide-spread and systematic heterogeneity in crossover rate along chromosomes generates predictable broad-scale trends in genetic diversity and population differentiation by modifying the impact of natural selection among regions within a genome. We conclude by emphasizing that chromosome-scale heterogeneity in crossover rate should urgently be incorporated into analytical tools in evolutionary genomics, and in the interpretation of emerging patterns.

opencc-zeroDec 2014View details →

Data from: Microform-scale variations in peatland permeability and their ecohydrological implications

1. The acrotelm-catotelm model of peatland hydrological and biogeochemical processes posits that the permeability of raised bogs is largely homogenous laterally but varies strongly with depth through the soil profile; uppermost peat layers are highly permeable while deeper layers are, effectively, impermeable. 2. We measured down-core changes in peat permeability, plant macrofossil assemblages, dry bulk density and degree of humification beneath two types of characteristic peatland microform – ridges and hollows – at a raised bog in Wales. Six 14C dates were also collected for one hollow and an adjacent ridge. 3. Contrary to the acrotelm-catotelm model, we found that deeper peat can be as highly permeable as near-surface peat and that its permeability can vary by more than an order of magnitude between microforms over horizontal distances of 1-5 metres. 4. Our palaeo-ecological data paint a complicated picture of microform persistence. Some microforms can remain in the same position on a bog for millennia, growing vertically upwards as the bog grows. However, adjacent areas on the bog (< 10 m distant) show switches between microform type over time, indicating a lack of persistence. 5. Synthesis. We suggest that the acrotelm-catotelm model should be used cautiously; spatial variations in peatland permeability do not fit the simple patterns suggested by the model. To understand how peatlands as a whole function both hydrologically and ecologically it is necessary to understand how patterns of peat physical properties and peatland vegetation develop and persist.

opencc-zeroDec 2014View details →

Data from: Shape up or ship out: migratory behaviour predicts morphology across spatial scale in a freshwater fish

1. Migration is a widespread phenomenon, with powerful ecological and evolutionary consequences. Morphological adaptations to reduce the energetic costs associated with migratory transport are commonly documented for migratory species. However, few studies have investigated whether variation in body morphology can be explained by variation in migratory strategy within a species. 2. We address this question in roach Rutilus rutilus, a partially migratory freshwater fish that migrates from lakes into streams during winter. We both compare body shape between populations that differ in migratory opportunity (open vs. closed lakes), and between individuals from a single population that vary in migratory propensity (migrants and residents from a partially migratory population). Following hydrodynamic theory, we posit that migrants should have a more shallow body depth, to reduce the costs associated with migrating into streams with higher flow conditions than the lakes the residents occupy all year round. 3. We find evidence both across and within populations to support our prediction, with individuals from open lakes and migrants from the partially migratory population having a more slender, shallow-bodied morphology than fish from closed lakes and all-year residents. 4. Our data suggest that a shallow body morphology is beneficial to migratory individuals and our study is one of the first to link migratory strategy and intraspecific variation in body shape.

opencc-zeroDec 2013View details →

Data from: Target capture and massively parallel sequencing of ultraconserved elements for comparative studies at shallow evolutionary time scales

Comparative genetic studies of non-model organisms are transforming rapidly due to major advances in sequencing technology. A limiting factor in these studies has been the identification and screening of orthologous loci across an evolutionarily distant set of taxa. Here, we evaluate the efficacy of genomic markers targeting ultraconserved DNA elements (UCEs) for analyses at shallow evolutionary timescales. Using sequence capture and massively parallel sequencing to generate UCE data for five co-distributed Neotropical rainforest bird species, we recovered 776–1516 UCE loci across the five species. Across species, 53–77% of the loci were polymorphic, containing between 2.0 and 3.2 variable sites per polymorphic locus, on average. We performed species tree construction, coalescent modeling, and species delimitation, and we found that the five co-distributed species exhibited discordant phylogeographic histories. We also found that species trees and divergence times estimated from UCEs were similar to the parameters obtained from mtDNA. The species that inhabit the understory had older divergence times across barriers, contained a higher number of cryptic species, and exhibited larger effective population sizes relative to the species inhabiting the canopy. Because orthologous UCEs can be obtained from a wide array of taxa, are polymorphic at shallow evolutionary timescales, and can be generated rapidly at low cost, they are an effective genetic marker for studies investigating evolutionary patterns and processes at shallow timescales.

Data from: Stable coexistence of genetically divergent Atlantic cod ecotypes at multiple spatial scales

Coexistence in the same habitat of closely related yet genetically different populations is a phenomenon that challenges our understanding of local population structure and adaptation. Identifying the underlying mechanisms for such coexistence can yield new insight into adaptive evolution, diversification, and the potential for organisms to adapt and persist in response to a changing environment. Recent studies have documented cryptic, sympatric populations of Atlantic cod (Gadus morhua) in coastal areas. We analyzed genetic origin of 6483 individual cod sampled annually over 14 years from 125 locations along the Norwegian Skagerrak coast and document stable coexistence of two genetically divergent Atlantic cod ecotypes throughout the study area and study period. A 'fjord' ecotype dominated in numbers deep inside fjords while a 'North Sea' ecotype was the only type found in offshore North Sea. Both ecotypes coexisted in similar proportions throughout coastal habitats at all spatial scales. The size-at-age of the North Sea ecotype on average exceeded that of the fjord ecotype by 20% in length and 80% in weight across all habitats. Different growth and size among individuals of the two types might be one of several ecologically significant variables that allow for stable coexistence of closely related populations within the same habitat. Management plans, biodiversity initiatives, and other mitigation strategies that do not account for the mixture of species ecotypes are unlikely to meet objectives related to the sustainability of fish and fisheries.

Data from: An evaluation of transcriptome-based exon capture for frog phylogenomics across multiple scales of divergence (Class: Amphibia, Order: Anura)

Custom sequence capture experiments are becoming an efficient approach for gathering large sets of orthologous markers in nonmodel organisms. Transcriptome-based exon capture utilizes transcript sequences to design capture probes, typically using a reference genome to identify intron–exon boundaries to exclude shorter exons (<200 bp). Here, we test directly using transcript sequences for probe design, which are often composed of multiple exons of varying lengths. Using 1260 orthologous transcripts, we conducted sequence captures across multiple phylogenetic scales for frogs, including outgroups ~100 Myr divergent from the ingroup. We recovered a large phylogenomic data set consisting of sequence alignments for 1047 of the 1260 transcriptome-based loci (~561 000 bp) and a large quantity of highly variable regions flanking the exons in transcripts (~70 000 bp), the latter improving substantially by only including ingroup species (~797 000 bp). We recovered both shorter (<100 bp) and longer exons (>200 bp), with no major reduction in coverage towards the ends of exons. We observed significant differences in the performance of blocking oligos for target enrichment and nontarget depletion during captures, and differences in PCR duplication rates resulting from the number of individuals pooled for capture reactions. We explicitly tested the effects of phylogenetic distance on capture sensitivity, specificity, and missing data, and provide a baseline estimate of expectations for these metrics based on a priori knowledge of nuclear pairwise differences among samples. We provide recommendations for transcriptome-based exon capture design based on our results, cost estimates and offer multiple pipelines for data assembly and analysis.

Data from: Large-scale functional networks identified from resting-state EEG using spatial ICA

Several methods have been applied to EEG or MEG signals to detect functional networks. In recent works using MEG/EEG and fMRI data, temporal ICA analysis has been used to extract spatial maps of resting-state networks with or without an atlas-based parcellation of the cortex. Since the links between the fMRI signal and the electromagnetic signals are not fully established, and to avoid any bias, we examined whether EEG alone was able to derive the spatial distribution and temporal characteristics of functional networks. To do so, we propose a two-step original method: 1) An individual multi-frequency data analysis including EEG-based source localisation and spatial independent component analysis, which allowed us to characterize the resting-state networks. 2) A group-level analysis involving a hierarchical clustering procedure to identify reproducible large-scale networks across the population. Compared with large-scale resting-state networks obtained with fMRI, the proposed EEG-based analysis revealed smaller independent networks thanks to the high temporal resolution of EEG, hence hierarchical organization of networks. The comparison showed a substantial overlap between EEG and fMRI networks in motor, premotor, sensory, frontal, and parietal areas. However, there were mismatches between EEG-based and fMRI-based networks in temporal areas, presumably resulting from a poor sensitivity of fMRI in these regions or artefacts in the EEG signals. The proposed method opens the way for studying the high temporal dynamics of networks at the source level thanks to the high temporal resolution of EEG. It would then become possible to study detailed measures of the dynamics of connectivity.

opencc-zeroNov 2019View details →

Data from: Fine-scale appendage structure of the Cambrian trilobitomorph Naraoia spinosa and its ontogenetic and ecological implications

Trilobitomorphs are a species-rich Palaeozoic arthropod assemblage that unites trilobites with several other lineages that share similar appendage structure. Post-embryonic development of the exoskeleton is well documented for some trilobitomorphs, especially trilobites, but little is known of the ontogeny of their soft parts, limiting understanding of their autecology. Here we document appendage structure of the Cambrian naraoiid trilobitomorph Naraoia spinosa by computed microtomography, resulting in three-dimensional reconstructions of appendages at both juvenile and adult stages. The adult has dense, strong spines on the protopods of post-antennal appendages, implying a predatory/scavenging behavior. Absence of such gnathobasic structures but instead tiny protopodal bristles and a number of endopodal setae suggests a detritus-feeding strategy for the juvenile. Our data add strong morphological evidence for ecological-niche shifting by Cambrian arthropods during their life cycles. A conserved number of appendages across the sampled developmental stages demonstrates that Naraoia ceased budding off new appendages by the mid-juvenile stage.

Data from: Correlated paternity measures mate monopolization and scales with the magnitude of sexual selection

Indirect measures of sexual selection have been criticized because they can overestimate the magnitude of selection. In particular, they do not account for the degree to which mating opportunities can be monopolized by individuals of the sex that compete for mates. We introduce a measure of mate monopolization (m) based on the magnitude of correlated paternity and evaluate its ability to track changes in the magnitude of sexual selection. Simulation models were used to compare how well m tracked changes in the selection differential (s) for a trait regulating mating success. We further evaluated the association between m and s using ten replicate mating arrays of Sagittaria latifolia in which plants with contrasting alleles at microsatellite loci competed for siring opportunities. The computer models and mating arrays both demonstrated a positive linear association between m and s, supporting the utility of m as an index of sexual selection. Commonly used measures of sexual selection are not easily applied to organisms, such as the flowering plants, for which mating events are difficult to observe. The measure of mate monopolization introduced here could prove to be a useful addition to studies of sexual selection in these organisms.

Data from: Allometry of wing twist and camber in a flower chafer during free flight: how do wing deformations scale with body size?

Intraspecific variation in adult body mass can be particularly high in some insect species, mandating adjustment of the wing's structural properties to support the weight of the larger body mass in air. Insect wings elastically deform during flapping, dynamically changing the twist and camber of the relatively thin and flat aerofoil. We examined how wing deformations during free flight scale with body mass within a species of rose chafers (Coleoptera: Protaetia cuprea) in which individuals varied more than threefold in body mass (0.38–1.29 g). Beetles taking off voluntarily were filmed using three high-speed cameras and the instantaneous deformation of their wings during the flapping cycle was analysed. Flapping frequency decreased in larger beetles but, otherwise, flapping kinematics remained similar in both small and large beetles. Deflection of the wing chord-wise varied along the span, with average deflections at the proximal trailing edge higher by 0.2 and 0.197 wing lengths compared to the distal trailing edge in the downstroke and the upstroke, respectively. These deflections scaled with wing chord to the power of 1.0, implying a constant twist and camber despite the variations in wing and body size. This suggests that the allometric growth in wing size includes adjustment of the flexural stiffness of the wing structure to preserve wing twist and camber during flapping.