ReadTools: A universal toolkit for handling sequence data from different sequencing platforms

The genomic distribution of transposable elements is driven by spatially variable purifying selection

It is widely accepted that the genomic distribution of transposable elements (TEs) mainly reflects the outcome of purifying selection and insertion bias (1). Nevertheless, the relative importance of these two evolutionary forces could not be tested thoroughly. Here, we introduce an experimental system, which allows separating purifying selection from TE insertion bias. We used experimental evolution to study the TE insertion patterns in Drosophila simulans founder populations harboring 1040 insertions of an active P-element. After 10 generations at a large population size, we detected strong selection against P-element insertions. The exception were P-element insertions in genomic regions for which a strong insertion bias has been proposed (2-4). Because recurrent P-element insertions cannot explain this pattern, we conclude that purifying selection, with variable strength along the chromosomes, is the major determinant of the genomic distribution of P-elements. Genomic regions with relaxed purifying selection against P-element insertions exhibit normal levels of purifying selection against base substitutions. This suggests that different types of purifying selection operate on base substitutions and P-element insertions. Our results highlight the power of experimental evolution to understand basic evolutionary processes, which are difficult to infer from patterns of natural variation alone.

Human skin-resident CD8+ T cells require RUNX2 and RUNX3 for induction of cytotoxicity and expression of the integrin CD49a

The integrin CD49a marks highly cytotoxic epidermal-tissue-resident memory (T_RM) cells, but their differentiation from circulating populations remains poorly defined. We demonstrate enrichment of RUNT family transcription-factor-binding motifs in human epidermal CD8⁺CD103⁺CD49a⁺ T_RM cells, paralleled by high RUNX2 and RUNX3 protein expression. Sequencing of paired skin and blood samples revealed clonal overlap between epidermal CD8⁺CD103⁺CD49a⁺ T_RM cells and circulating memory CD8⁺CD45RA^-CD62L⁺ T cells. In vitro stimulation of circulating CD8⁺CD45RA^-CD62L⁺ T cells with IL-15 and TGF-β induced CD49a expression and cytotoxic transcriptional profiles in a RUNX2- and RUNX3-dependent manner. We therefore identified a reservoir of circulating cells with cytotoxic T_RM potential. In melanoma patients, high RUNX2, but not RUNX3, transcription correlated with a cytotoxic CD8⁺CD103⁺CD49a⁺ T_RM cell signature and improved patient survival. Together, our results indicate that combined RUNX2 and RUNX3 activity promotes the differentiation of cytotoxic CD8⁺CD103⁺CD49a⁺ T_RM cells, providing immunosurveillance of infected and malignant cells.

Epigenetic regulation of T cell lineages in skin and blood following hematopoietic stem cell transplantation

Allogeneic hematopoietic stem-cell transplantation (HSCT) seeks to reconstitute the host's immune system from donor stem cells. The success of HSCT is threatened by complications including leukemia relapse or graft-versus-host-disease (GvHD). To investigate the underlying regulatory processes in central and peripheral T cell recovery, we performed sequential multi-omics analysis of T cells of the skin and blood during HSCT. We detected rapid effector T cell reconstitution, while emergence of regulatory T cells was delayed. Epigenetic and gene-regulatory programs were associated with recovering T cells and diverged greatly between skin and blood T cells. The BRG1/BRM-associated factor chromatin remodeling complex and histone deacetylases (HDACs) were epigenetic regulators involved in restoration of T cell homeostasis after transplantation. In isolated T cells of patients after HSCT, we observed class I HDAC-inhibitors to modulate their dysbalance. The present study highlights the importance of epigenetic regulation in the recovery of T cells following HSCT.

Disturbances in microbial skin recolonization and cutaneous immune response following allogeneic stem cell transfer

The composition of the gut microbiome influences the clinical course after allogeneic hematopoietic stem cell transplantation (HSCT), but little is known about the relevance of skin microorganisms. In a single-center, observational study, we recruited a cohort of 50 patients before undergoing conditioning treatment and took both stool and skin samples up to one year after HSCT. We could confirm intestinal dysbiosis following HSCT and report that the skin microbiome is likewise perturbed in HSCT-recipients. Overall bacterial colonization of the skin was decreased after conditioning. Particularly patients that developed acute skin graft-versus-host disease (aGVHD) presented with an overabundance of Staphylococcus spp. In addition, a loss in alpha diversity was indicative of aGVHD development already before disease onset and correlated with disease severity. Further, co-localization of CD45⁺ leukocytes and staphylococci was observed in the skin of aGVHD patients even before disease development and paralleled with upregulated genes required for antigen-presentation in mononuclear phagocytes. Overall, our data reveal disturbances of the skin microbiome as well as cutaneous immune response in HSCT recipients with changes associated with cutaneous aGVHD.

CD45RA+CD62L- ILCs in human tissues represent a quiescent local reservoir for the generation of differentiated ILCs

Innate lymphoid cells (ILCs) are highly plastic and predominantly mucosal tissue-resident cells that contribute to both homeostasis and inflammation depending on the microenvironment. The discovery of naïve-like ILCs suggests an ILC differentiation process that is akin to naïve T cell differentiation. Delineating the mechanisms that underlie ILC differentiation in tissues is crucial for understanding ILC biology in health and disease. Here, we showed that tonsillar ILCs expressing CD45RA lacked proliferative activity, indicative of cellular quiescence. CD62L distinguished two subsets of CD45RA⁺ ILCs. CD45RA⁺CD62L⁺ ILCs (CD62L⁺ ILCs) resembled circulating naïve ILCs because they lacked the transcriptional, metabolic, epigenetic, and cytokine production signatures of differentiated ILCs. CD45RA⁺CD62L^- ILCs (CD62L^- ILCs) were epigenetically similar to CD62L⁺ ILCs but showed a transcriptional, metabolic, and cytokine production signature that was more akin to differentiated ILCs. CD62L⁺ and CD62L^- ILCs contained uni- and multipotent precursors of ILC1s/NK cells and ILC3s. Differentiation of CD62L⁺ and CD62L^- ILCs led to metabolic reprogramming including up-regulation of genes associated with glycolysis, which was needed for their effector functions after differentiation. CD62L^- ILCs with preferential differentiation capacity toward IL-22-producing ILC3s accumulated in the inflamed mucosa of patients with inflammatory bowel disease. These data suggested distinct differentiation potential of CD62L⁺ and CD62L^- ILCs between tissue microenvironments and identified that manipulation of these cells is a possible approach to restore tissue-immune homeostasis.

Y-Chromosomal Insights into Breeding History and Sire Line Genealogies of Arabian Horses

The Y chromosome is a valuable genetic marker for studying the origin and influence of paternal lineages in populations. In this study, we conducted Y-chromosomal lineage-tracing in Arabian horses. First, we resolved a Y haplotype phylogeny based on the next generation sequencing data of 157 males from several breeds. Y-chromosomal haplotypes specific for Arabian horses were inferred by genotyping a collection of 145 males representing most Arabian sire lines that are active around the globe. These lines formed three discrete haplogroups, and the same haplogroups were detected in Arabian populations native to the Middle East. The Arabian haplotypes were clearly distinct from the ones detected in Akhal Tekes, Turkoman horses, and the progeny of two Thoroughbred foundation sires. However, a haplotype introduced into the English Thoroughbred by the stallion Byerley Turk (1680), was shared among Arabians, Turkomans, and Akhal Tekes, which opens a discussion about the historic connections between Oriental horse types. Furthermore, we genetically traced Arabian sire line breeding in the Western World over the past 200 years. This confirmed a strong selection for relatively few male lineages and uncovered incongruences to written pedigree records. Overall, we demonstrate how fine-scaled Y-analysis contributes to a better understanding of the historical development of horse breeds.

Delayed antiretroviral therapy in HIV-infected individuals leads to irreversible depletion of skin- and mucosa-resident memory T cells

People living with HIV (PLWH) are at increased risk for developing skin and mucosal malignancies despite systemic reconstitution of CD4⁺ T cells upon antiretroviral therapy (ART). The underlying mechanism of chronic tissue-related immunodeficiency in HIV is unclear. We found that skin CD4⁺ tissue-resident memory T (Trm) cells were depleted after HIV infection and replenished only upon early ART initiation. TCR clonal analysis following early ART suggested a systemic origin for reconstituting CD4⁺ Trm cells. Single-cell RNA sequencing in PLWH that received late ART treatment revealed a loss of CXCR3⁺ Trm cells and a tolerogenic skin immune environment. Human papilloma virus-induced precancerous lesion biopsies showed reduced CXCR3⁺ Trm cell frequencies in the mucosa in PLWH versus HIV^- individuals. These results reveal an irreversible loss of CXCR3⁺ Trm cells confined to skin and mucosa in PLWH who received late ART treatment, which may be a precipitating factor in the development of HPV-related cancer.

Evolution of phenotypic variance in response to a novel hot environment

Shifts in trait means are widely considered as evidence for adaptive responses, but the impact on phenotypic variance remains largely unexplored. Classic quantitative genetics provides a theoretical framework to predict how selection on phenotypic mean affects the variance. In addition to this indirect effect, it is also possible that the variance of the trait is the direct target of selection, but experimentally characterized cases are rare. Here, we studied gene expression variance of Drosophila simulans males before and after 100 generations of adaptation to a novel hot laboratory environment. In each of the two independently evolved populations, the variance of 125 and 97 genes was significantly reduced. We propose that the drastic loss in environmental complexity from nature to the laboratory may have triggered selection for reduced variance. Our observation that selection could drive changes in the variance of gene expression could have important implications for studies of adaptation processes in natural and experimental populations.

Highly Parallel Genomic Selection Response in Replicated Drosophila melanogaster Populations with Reduced Genetic Variation

Many adaptive traits are polygenic and frequently more loci contributing to the phenotype are segregating than needed to express the phenotypic optimum. Experimental evolution with replicated populations adapting to a new controlled environment provides a powerful approach to study polygenic adaptation. Because genetic redundancy often results in nonparallel selection responses among replicates, we propose a modified evolve and resequence (E&R) design that maximizes the similarity among replicates. Rather than starting from many founders, we only use two inbred Drosophila melanogaster strains and expose them to a very extreme, hot temperature environment (29 °C). After 20 generations, we detect many genomic regions with a strong, highly parallel selection response in 10 evolved replicates. The X chromosome has a more pronounced selection response than the autosomes, which may be attributed to dominance effects. Furthermore, we find that the median selection coefficient for all chromosomes is higher in our two-genotype experiment than in classic E&R studies. Because two random genomes harbor sufficient variation for adaptive responses, we propose that this approach is particularly well-suited for the analysis of polygenic adaptation.

The genetic architecture of temperature adaptation is shaped by population ancestry and not by selection regime

Background

Understanding the genetic architecture of temperature adaptation is key for characterizing and predicting the effect of climate change on natural populations. One particularly promising approach is Evolve and Resequence, which combines advantages of experimental evolution such as time series, replicate populations, and controlled environmental conditions, with whole genome sequencing. Recent analysis of replicate populations from two different Drosophila simulans founder populations, which were adapting to the same novel hot environment, uncovered very different architectures—either many selection targets with large heterogeneity among replicates or fewer selection targets with a consistent response among replicates.

Results

Here, we expose the founder population from Portugal to a cold temperature regime. Although almost no selection targets are shared between the hot and cold selection regime, the adaptive architecture was similar. We identify a moderate number of targets under strong selection (19 selection targets, mean selection coefficient = 0.072) and parallel responses in the cold evolved replicates. This similarity across different environments indicates that the adaptive architecture depends more on the ancestry of the founder population than the specific selection regime.

Conclusions

These observations will have broad implications for the correct interpretation of the genomic responses to a changing climate in natural populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02425-9.

The transcription factor Bcl11b promotes both canonical and adaptive NK cell differentiation

Epigenetic landscapes can provide insight into regulation of gene expression and cellular diversity. Here, we examined the transcriptional and epigenetic profiles of seven human blood natural killer (NK) cell populations, including adaptive NK cells. The BCL11B gene, encoding a transcription factor (TF) essential for T cell development and function, was the most extensively regulated, with expression increasing throughout NK cell differentiation. Several Bcl11b-regulated genes associated with T cell signaling were specifically expressed in adaptive NK cell subsets. Regulatory networks revealed reciprocal regulation at distinct stages of NK cell differentiation, with Bcl11b repressing RUNX2 and ZBTB16 in canonical and adaptive NK cells, respectively. A critical role for Bcl11b in driving NK cell differentiation was corroborated in BCL11B-mutated patients and by ectopic Bcl11b expression. Moreover, Bcl11b was required for adaptive NK cell responses in a murine cytomegalovirus model, supporting expansion of these cells. Together, we define the TF regulatory circuitry of human NK cells and uncover a critical role for Bcl11b in promoting NK cell differentiation and function.

Parallel gene expression evolution in natural and laboratory evolved populations

Ecological adaptation is frequently inferred by the comparison of natural populations from different environments. Nevertheless, inference of the selective forces suffers the challenge that many environmental factors covary. With well-controlled environmental conditions, experimental evolution provides a powerful approach to complement the analysis of natural populations. On the other hand, it is apparent that laboratory conditions differ in many ways from natural environments, which raises the question as to what extent selection responses in experimental evolution studies can inform us about adaptation processes in the wild. In this study, we compared the expression profiles of replicated Drosophila melanogaster populations which have been exposed to two distinct temperature regimes (18/28 and 10/20°C) in the laboratory for more than 80 generations. Using gene-wise differential expression analysis and co-expression network analysis, we identified 541 genes and three coregulated gene modules that evolved in the same direction in both temperature regimes, and most of these changes probably reflect an adaptation to the space constraint or diurnal temperature fluctuation that is common in both selection regimes. In total, 203 genes and seven modules evolved temperature-specific expression changes. Remarkably, we detected a significant overlap of these temperature-adaptive genes/modules from experimental evolution with temperature-adaptive genes inferred from natural Drosophila populations covering two different temperature clines. We conclude that well-designed experimental evolution studies are a powerful tool to dissect evolutionary responses.

Population genomics for wildlife conservation and management

Biodiversity is under threat worldwide. Over the past decade, the field of population genomics has developed across nonmodel organisms, and the results of this research have begun to be applied in conservation and management of wildlife species. Genomics tools can provide precise estimates of basic features of wildlife populations, such as effective population size, inbreeding, demographic history and population structure, that are critical for conservation efforts. Moreover, population genomics studies can identify particular genetic loci and variants responsible for inbreeding depression or adaptation to changing environments, allowing for conservation efforts to estimate the capacity of populations to evolve and adapt in response to environmental change and to manage for adaptive variation. While connections from basic research to applied wildlife conservation have been slow to develop, these connections are increasingly strengthening. Here we review the primary areas in which population genomics approaches can be applied to wildlife conservation and management, highlight examples of how they have been used, and provide recommendations for building on the progress that has been made in this field.

Long-term skin-resident memory T cells proliferate in situ and are involved in human graft-versus-host disease

The skin contains a population of tissue-resident memory T cells (Trm) that is thought to contribute to local tissue homeostasis and protection against environmental injuries. Although information about the regulation, survival program, and pathophysiological roles of Trm has been obtained from murine studies, little is known about the biology of human cutaneous Trm Here, we showed that host-derived CD69+ αβ memory T cell clones in the epidermis and dermis remain stable and functionally competent for at least 10 years in patients with allogeneic hematopoietic stem cell transplantation. Single-cell RNA sequencing revealed low expression of genes encoding tissue egress molecules by long-term persisting Trm in the skin, whereas tissue retention molecules and stem cell markers were displayed by Trm The transcription factor RUNX3 and the surface molecule galectin-3 were preferentially expressed by host T cells at the RNA and protein levels, suggesting two new markers for human skin Trm Furthermore, skin lesions from patients developing graft-versus-host disease (GVHD) showed a large number of cytokine-producing host-derived Trm, suggesting a contribution of these cells to the pathogenesis of GVHD. Together, our studies highlighted the relationship between the local human skin environment and long-term persisting Trm, which differs from murine skin. Our results also indicated that local tissue inflammation occurs through host-derived Trm after allogeneic hematopoietic stem cell transplantation.

Fitness effects for Ace insecticide resistance mutations are determined by ambient temperature

BACKGROUND

Insect pest control programs often use periods of insecticide treatment with intermittent breaks, to prevent fixing of mutations conferring insecticide resistance. Such mutations are typically costly in an insecticide-free environment, and their frequency is determined by the balance between insecticide treatment and cost of resistance. Ace, a key gene in neuronal signaling, is a prominent target of many insecticides and across several species, three amino acid replacements (I161V, G265A, and F330Y) provide resistance against several insecticides. Because temperature disturbs neuronal signaling homeostasis, we reasoned that the cost of insecticide resistance could be modulated by ambient temperature.

RESULTS

Experimental evolution of a natural Drosophila simulans population at hot and cold temperature regimes uncovered a surprisingly strong effect of ambient temperature. In the cold temperature regime, the resistance mutations were strongly counter selected (s = - 0.055), but in a hot environment, the fitness costs of resistance mutations were reduced by almost 50% (s = - 0.031). We attribute this unexpected observation to the advantage of the reduced enzymatic activity of resistance mutations in hot environments.

CONCLUSION

We show that fitness costs of insecticide resistance genes are temperature-dependent and suggest that the duration of insecticide-free periods need to be adjusted for different climatic regions to reflect these costs. We suggest that such environment-dependent fitness effects may be more common than previously assumed and pose a major challenge for modeling climate change.

A discrete subset of epigenetically primed human NK cells mediates antigen-specific immune responses

Adaptive features of natural killer (NK) cells have been reported in various species with different underlying mechanisms. It is unclear, however, which NK cell populations are capable of mounting antigen-specific recall responses and how such functions are regulated at the molecular level. Here, we identify and characterize a discrete population of CD49a+CD16- NK cells in the human liver that displays increased epigenetic potential to elicit memory responses and has the functional properties to exert antigen-specific immunity in the skin as an effector site. Integrated chromatin-based epigenetic and transcriptomic profiling revealed unique characteristics of hepatic CD49a+CD16- NK cells when compared with conventional CD49a-CD16+ NK cells, thereby defining active genomic regions and molecules underpinning distinct NK cell reactivity. In contrast to conventional NK cells, our results suggest that adaptive CD49a+CD16- NK cells are able to bypass the KIR receptor-ligand system upon antigen-specific stimulation. Furthermore, these cells were highly migratory toward chemokine gradients expressed in epicutaneous patch test lesions as an effector site of adaptive immune responses in the skin. These results define pathways operative in human antigen-specific adaptive NK cells and provide a roadmap for harnessing this NK cell subset for specific therapeutic or prophylactic vaccine strategies.

Long-Term Dynamics Among Wolbachia Strains During Thermal Adaptation of Their Drosophila melanogaster Hosts

Climate change is a major evolutionary force triggering thermal adaptation in a broad range of species. While the consequences of global warming are being studied for an increasing number of species, limited attention has been given to the evolutionary dynamics of endosymbionts in response to climate change. Here, we address this question by studying the dynamics of Wolbachia, a well-studied endosymbiont of Drosophila melanogaster. D. melanogaster populations infected with 13 different Wolbachia strains were exposed to novel hot and cold laboratory environments for up to 180 generations. The short-term dynamics suggested a temperature-related fitness difference resulting in the increase of clade V strains in the cold environment only. Our long-term analysis now uncovers that clade V dominates in all replicates after generation 60 irrespective of temperature treatment. We propose that adaptation of the Drosophila host to either temperature or Drosophila C virus (DCV) infection are the cause of the replicated, temporally non-concordant Wolbachia dynamics. Our study provides an interesting case demonstrating that even simple, well-controlled experiments can result in complex, but repeatable evolutionary dynamics, thus providing a cautionary note on too simple interpretations on the impact of climate change.

Neuronal Function and Dopamine Signaling Evolve at High Temperature in Drosophila

Neuronal activity is temperature sensitive and affects behavioral traits important for individual fitness, such as locomotion and courtship. Yet, we do not know enough about the evolutionary response of neuronal phenotypes in new temperature environments. Here, we use long-term experimental evolution of Drosophila simulans populations exposed to novel temperature regimes. Here, we demonstrate a direct relationship between thermal selective pressure and the evolution of neuronally expressed molecular and behavioral phenotypes. Several essential neuronal genes evolve lower expression at high temperatures and higher expression at low temperatures, with dopaminergic neurons standing out by displaying the most consistent expression change across independent replicates. We functionally validate the link between evolved gene expression and behavioral changes by pharmacological intervention in the experimentally evolved D. simulans populations as well as by genetically triggered expression changes of key genes in D. melanogaster. As natural temperature clines confirm our results for Drosophila and Anopheles populations, we conclude that neuronal dopamine evolution is a key factor for temperature adaptation.

Anti-Apoptotic Molecule BCL2 Is a Therapeutic Target in Steroid-Refractory Graft-Versus-Host Disease

Graft-versus-host disease (GVHD) is the leading cause of mortality after hematopoietic stem cell transplantation and primarily affects barrier organs such as the skin. One-third of cases are refractory to steroid treatment resulting in poor outcomes and the need for novel therapies. Longitudinal analysis of T-cell transcriptomes in patients before the appearance of GVHD symptoms revealed the upregulation of anti-apoptotic regulator B-cell lymphoma 2 (BCL2) at GVHD initiation. To determine the potential of BCL2 inhibition in active GVHD, we analyzed tissues of 88 patients with acute or chronic GVHD. BCL2 RNA was elevated in multiple organs affected by GVHD and expression correlated with transplant-related mortality and steroid-refractory GVHD. BCL2-expressing lymphocytes were present in skin lesions and peripheral blood of patients with acute and chronic GVHD. Inhibition of BCL2 increased the CD4 to CD8 ratio in allogeneic T cells in vitro and induced apoptosis of T cells from patients with steroid-pretreated chronic GVHD ex vivo. In addition, the higher ratio of regulatory to nonregulatory T cells upon blockage of BCL2 could add to the anti-inflammatory effect of BCL2 blockage. Collectively, our results highlight BCL2 as an important factor for GVHD development and introduce BCL2 inhibition as previously unreported and urgently needed targeted therapy in the treatment of steroid-refractory GVHD.

Rapid sex-specific adaptation to high temperature in Drosophila

The pervasive occurrence of sexual dimorphism demonstrates different adaptive strategies of males and females. While different reproductive strategies of the two sexes are well-characterized, very little is known about differential functional requirements of males and females in their natural habitats. Here, we study the impact environmental change on the selection response in both sexes. Exposing replicated Drosophila populations to a novel temperature regime, we demonstrate sex-specific changes in gene expression, metabolic and behavioral phenotypes in less than 100 generations. This indicates not only different functional requirements of both sexes in the new environment but also rapid sex-specific adaptation. Supported by computer simulations we propose that altered sex-biased gene regulation from standing genetic variation, rather than new mutations, is the driver of rapid sex-specific adaptation. Our discovery of environmentally driven divergent functional requirements of males and females has important implications-possibly even for gender aware medical treatments.

The horse Y chromosome as an informative marker for tracing sire lines

Analysis of the Y chromosome is the best-established way to reconstruct paternal family history in humans. Here, we applied fine-scaled Y-chromosomal haplotyping in horses with biallelic markers and demonstrate the potential of our approach to address the ancestry of sire lines. We de novo assembled a draft reference of the male-specific region of the Y chromosome from Illumina short reads and then screened 5.8 million basepairs for variants in 130 specimens from intensively selected and rural breeds and nine Przewalski's horses. Among domestic horses we confirmed the predominance of a young'crown haplogroup' in Central European and North American breeds. Within the crown, we distinguished 58 haplotypes based on 211 variants, forming three major haplogroups. In addition to two previously characterised haplogroups, one observed in Arabian/Coldblooded and the other in Turkoman/Thoroughbred horses, we uncovered a third haplogroup containing Iberian lines and a North African Barb Horse. In a genealogical showcase, we distinguished the patrilines of the three English Thoroughbred founder stallions and resolved a historic controversy over the parentage of the horse 'Galopin', born in 1872. We observed two nearly instantaneous radiations in the history of Central and Northern European Y-chromosomal lineages that both occurred after domestication 5,500 years ago.

A 24 h Age Difference Causes Twice as Much Gene Expression Divergence as 100 Generations of Adaptation to a Novel Environment

Gene expression profiling is one of the most reliable high-throughput phenotyping methods, allowing researchers to quantify the transcript abundance of expressed genes. Because many biotic and abiotic factors influence gene expression, it is recommended to control them as tightly as possible. Here, we show that a 24 h age difference of Drosophilasimulans females that were subjected to RNA sequencing (RNA-Seq) five and six days after eclosure resulted in more than 2000 differentially expressed genes. This is twice the number of genes that changed expression during 100 generations of evolution in a novel hot laboratory environment. Importantly, most of the genes differing in expression due to age introduce false positives or negatives if an adaptive gene expression analysis is not controlled for age. Our results indicate that tightly controlled experimental conditions, including precise developmental staging, are needed for reliable gene expression analyses, in particular in an evolutionary framework.