Phylogeny of locusts and grasshoppers reveals complex evolution of density-dependent phenotypic plasticity

Genomic and transcriptomic perspectives on the origin and evolution of NUMTs in Orthoptera

Nuclear mitochondrial pseudogenes (NUMTs) result from the transfer of mitochondrial DNA (mtDNA) to the nuclear genome. NUMTs, as "frozen" snapshots of mitochondria, can provide insights into diversification patterns. In this study, we analyzed the origins and insertion frequency of NUMTs using genome assembly data from ten species in Orthoptera. We found divergences between NUMTs and contemporary mtDNA in Orthoptera ranging from 0 % to 23.78 %. The results showed that the number of NUMT insertions was significantly positively correlated with the content of transposable elements in the genome. We found that 39.09 %-68.65 % of the NUMTs flanking regions (2,000 bp) contained retrotransposons, and more NUMTs originated from mitochondrial rDNA regions. Based on the analysis of the mitochondrial transcriptome, we found a potential mechanism of NUMT integration: mitochondrial transcripts are reverse transcribed into double-stranded DNA and then integrated into the genome. The probability of this mechanism occurring accounts for 0.30 %-1.02 % of total mitochondrial nuclear transfer events. Finally, based on the phylogenetic tree constructed using NUMTs and contemporary mtDNA, we provide insights into ancient evolutionary events such as species-specific "autaponumts" and "synaponumts" shared among different species, as well as post-integration duplication events.

Influence of climate on desert locust (Schistocerca gregaria Forskål, 1775) Plague and migration prediction in tropics

The outbreak of the desert locust Schistocerca gregaria Forskål, 1775, which originated from the Horn of Africa in 2019-2020 created an episodic plague under bio-geographical settings in the arid and semi-arid areas of South and Southwest Asia. In India, it happened after twenty-seven years due to the persistence of a few favourable conditions caused by its plague, resulting in hundreds of crores in crop damage. Keeping this in mind, the study aims to assess the suitability and likelihood of the desert locust epidemic occurring in India, utilizing two widely recognized statistical models: Weight-of-Evidence (WoE) and Frequency Ratio (FR). This work evaluated nine critical climatic factors for the study considering western and central parts of India. The 'Projected Locust Suitability' (PLS) was calculated by analyzing the correlation of the considered variables and the occurrence of locust swarms and bands. The significance (importance) of each variable on PLS was determined using Principal Component Analysis (PCA) and Random Forest (RF) algorithms. The PLS maps clearly show that 42.7-52.8% of the areas fall under high and very high locust suitability zones. The result suggests that the Ajmer-Gwalior-Allahabad tract is highly prone to future locust occurrences, while the Aligarh-Bareilly-Lakhimpur tract is moderately susceptible. The effectiveness of both modelled PLS maps was determined with the help of the ROC curve. The AUC results indicate that both the WoE (0.92) and the RF (0.90) models worked remarkably well in precisely predicting PLS. The RF-based IncNodePurity analysis indicates that low to moderate temperatures in the presence of cloud cover significantly impact locust occurrence and migration. The present findings are projected to direct the development of sustainable locust management strategies utilizing proper land use policies in the tropical climate.

An Overview of Orthoptera Mass Occurrences in Croatia from 1900 to 2023

During the last century, well-known locust species, such as Calliptamus italicus and Dociostaurus maroccanus, have produced outbreaks of varying degrees in the Balkans. The literature data on outbreaks in the region are scarce, and Croatia is not an exception. This study summarized the data on 23 Orthoptera mass occurrences in Croatia from 1900 to 2023 from 28 localities, representing 12 species. This is a low level of outbreak activity compared with other locust and pest grasshopper species in other parts of the world. The species with the most reporting is C. italicus with altogether six mass occurrences, while second is Barbitistes ocskayi and Miramella irena with three records, and in the third, place D. maroccanus and Gryllotalpa sp., each with two mass occurrences having been reported. One of the most notable swarms is that of Anacridium aegyptium which occurred around Šibenik in 1998, and this paper provides the first account of it, 25 years after it took place. The most recent outbreaks took place in 2022, and the most notable one was that of D. maroccanus swarm in Štikovo. The 2022 and 2023 reports were brief and muted, despite the affected agriculturists claiming significant damages.

The time course of behavioural phase change in the Central American locust Schistocerca piceifrons

Locusts exhibit an extreme form of phenotypic plasticity and can exist as two alternative phenotypes, known as solitarious and gregarious phases. These phases, which can transform from one to another depending on local population density, show distinctly different behavioural characteristics. The proximate mechanisms of behavioural phase polyphenism have been well studied in the desert locust Schistocerca gregaria and the migratory locust Locusta migratoria, and what is known in these species is often treated as a general feature of locusts. However, this approach might be flawed, given that there are approximately 20 locust species that have independently evolved phase polyphenism. Using the Central American locust Schistocerca piceifrons as a study system, we characterised the time course of behavioural phase change using standard locust behavioural assays, using both a logistic regression-based model and analyses of separate behavioural variables. We found that for nymphs of S. piceifrons, solitarisation was a relatively fast, two-step process, but that gregarisation was a much slower process. Additionally, the density of the gregarisation treatment seemed to have no effect on the rate of phase change. These data are at odds with what we know about the time course of behavioural phase change in S. gregaria, suggesting that the mechanisms of locust phase polyphenism in these two species are different and may not be phylogenetically constrained. Our study represents the most in-depth study of behavioural gregarisation and solitarisation in locusts to date.

Sexual repurposing of juvenile aposematism in locusts

Adaptive plasticity requires an integrated suite of functional responses to environmental variation, which can include social communication across life stages. Desert locusts (Schistocerca gregaria) exhibit an extreme example of phenotypic plasticity called phase polyphenism, in which a suite of behavioral and morphological traits differ according to local population density. Male and female juveniles developing at low population densities exhibit green- or sand-colored background-matching camouflage, while at high densities they show contrasting yellow and black aposematic patterning that deters predators. The predominant background colors of these phenotypes (green/sand/yellow) all depend on expression of the carotenoid-binding "Yellow Protein" (YP). Gregarious (high-density) adults of both sexes are initially pinkish, before a YP-mediated yellowing reoccurs upon sexual maturation. Yellow color is especially prominent in gregarious males, but the reason for this difference has been unknown since phase polyphenism was first described in 1921. Here, we use RNA interference to show that gregarious male yellowing acts as an intrasexual warning signal, which forms a multimodal signal with the antiaphrodisiac pheromone phenylacetonitrile (PAN) to prevent mistaken sexual harassment from other males during scramble mating in a swarm. Socially mediated reexpression of YP thus adaptively repurposes a juvenile signal that deters predators into an adult signal that deters undesirable mates. These findings reveal a previously underappreciated sexual dimension to locust phase polyphenism, and promote locusts as a model for investigating the relative contributions of natural versus sexual selection in the evolution of phenotypic plasticity.

A journey towards an integrated understanding of behavioural phase change in locusts

Behavioural phase change initiates and functionally couples the suite of traits that comprise density-dependent polyphenism in locusts. Here I provide a semi-expurgated account of my 25-year research journey studying behavioural phase transition in the desert locust. The journey spans continents, involves a cast of extraordinary colleagues, and travels across levels of biological organisation from deep within the nervous system of individual locusts to mass migration and the evolution and population dynamics of swarming.

Cytogenomic analysis unveils mixed molecular evolution and recurrent chromosomal rearrangements shaping the multigene families on Schistocerca grasshopper genomes

Multigene families are essential components of eukaryotic genomes and play key roles either structurally and functionally. Their modes of evolution remain elusive even in the era of genomics, because multiple multigene family sequences coexist in genomes, particularly in large repetitive genomes. Here, we investigate how the multigene families 18S rDNA, U2 snDNA, and H3 histone evolved in 10 species of Schistocerca grasshoppers with very large and repeat-enriched genomes. Using sequenced genomes and fluorescence in situ hybridization mapping, we find substantial differences between species, including the number of chromosomal clusters, changes in sequence abundance and nucleotide composition, pseudogenization, and association with transposable elements (TEs). The intragenomic analysis of Schistocerca gregaria using long-read sequencing and genome assembly unveils conservation for H3 histone and recurrent pseudogenization for 18S rDNA and U2 snDNA, likely promoted by association with TEs and sequence truncation. Remarkably, TEs were frequently associated with truncated copies, were also among the most abundant in the genome, and revealed signatures of recent activity. Our findings suggest a combined effect of concerted and birth-and-death models driving the evolution of multigene families in Schistocerca over the last 8 million years, and the occurrence of intra- and interchromosomal rearrangements shaping their chromosomal distribution. Despite the conserved karyotype in Schistocerca, our analysis highlights the extensive reorganization of repetitive DNAs in Schistocerca, contributing to the advance of comparative genomics for this important grasshopper genus.

Comparative analysis of phenotypic plasticity sheds light on the evolution and molecular underpinnings of locust phase polyphenism

Locusts exhibit one of nature's most spectacular examples of complex phenotypic plasticity, in which changes in density cause solitary and cryptic individuals to transform into gregarious and conspicuous locusts forming large migrating swarms. We investigated how these coordinated alternative phenotypes might have evolved by studying the Central American locust and three closely related non-swarming grasshoppers in a comparative framework. By experimentally isolating and crowding during nymphal development, we induced density-dependent phenotypic plasticity and quantified the resulting behavioural, morphological, and molecular reaction norms. All four species exhibited clear plasticity, but the individual reaction norms varied among species and showed different magnitudes. Transcriptomic responses were species-specific, but density-responsive genes were functionally similar across species. There were modules of co-expressed genes that were highly correlated with plastic reaction norms, revealing a potential molecular basis of density-dependent phenotypic plasticity. These findings collectively highlight the importance of studying multiple reaction norms from a comparative perspective.

The role of the neuropeptide [His7]-corazonin on phase-related characteristics in the Central American locust

Density-dependent phase polyphenism in locusts is one of the most extreme forms of phenotypic plasticity. Locusts exist along the continuum between two density-dependent phenotypes that differ in nymphal coloration, behavior, morphology, physiology, and reproduction among others. Nymphs of the solitarious phase, found in low population densities, are usually green, relatively inactive, and avoid each other, while gregarious nymphs, found in high density, exhibit a very obvious yellow/orange background with black patterning, and are highly active and attracted to each other. The multifunctional neuropeptide [His⁷]-corazonin has been shown to strongly affect black coloration and several other phase-related characteristics in at least two locust species, even though no effect on phase-related behavioral traits has been found. In this study, we investigate the role of [His⁷]-corazonin in the Central American locust Schistocerca piceifrons (Walker), which evolved density-dependent phase polyphenism independently from the two previously studied locust species. After successfully knocking down the transcript encoding [His⁷]-corazonin (CRZ) using RNA interference, we show that such a knockdown influences both color and morphometrics in this species, but does not influence phase-related behavioral traits. Our results suggest that the role of [His⁷]-corazonin is conserved in different locust species. Finally, our study represents the first controlled study of behavioral solitarization in S. piceifrons.

Comparative analysis of morabine grasshopper genomes reveals highly abundant transposable elements and rapidly proliferating satellite DNA repeats

BACKGROUND

Repetitive DNA sequences, including transposable elements (TEs) and tandemly repeated satellite DNA (satDNAs), collectively called the "repeatome", are found in high proportion in organisms across the Tree of Life. Grasshoppers have large genomes, averaging 9 Gb, that contain a high proportion of repetitive DNA, which has hampered progress in assembling reference genomes. Here we combined linked-read genomics with transcriptomics to assemble, characterize, and compare the structure of repetitive DNA sequences in four chromosomal races of the morabine grasshopper Vandiemenella viatica species complex and determine their contribution to genome evolution.

RESULTS

We obtained linked-read genome assemblies of 2.73-3.27 Gb from estimated genome sizes of 4.26-5.07 Gb DNA per haploid genome of the four chromosomal races of V. viatica. These constitute the third largest insect genomes assembled so far. Combining complementary annotation tools and manual curation, we found a large diversity of TEs and satDNAs, constituting 66 to 75% per genome assembly. A comparison of sequence divergence within the TE classes revealed massive accumulation of recent TEs in all four races (314-463 Mb per assembly), indicating that their large genome sizes are likely due to similar rates of TE accumulation. Transcriptome sequencing showed more biased TE expression in reproductive tissues than somatic tissues, implying permissive transcription in gametogenesis. Out of 129 satDNA families, 102 satDNA families were shared among the four chromosomal races, which likely represent a diversity of satDNA families in the ancestor of the V. viatica chromosomal races. Notably, 50 of these shared satDNA families underwent differential proliferation since the recent diversification of the V. viatica species complex.

CONCLUSION

This in-depth annotation of the repeatome in morabine grasshoppers provided new insights into the genome evolution of Orthoptera. Our TEs analysis revealed a massive recent accumulation of TEs equivalent to the size of entire Drosophila genomes, which likely explains the large genome sizes in grasshoppers. Despite an overall high similarity of the TE and satDNA diversity between races, the patterns of TE expression and satDNA proliferation suggest rapid evolution of grasshopper genomes on recent timescales.

Ontogenetic shift from aposematism and gregariousness to crypsis in a Romaleid grasshopper

Traits of chemically-defended animals can change as an individual grows and matures, and both theoretical and empirical evidence favour a direction of change from crypsis to aposematism. This study examines the suite of traits involved in an unusual opposite shift from aposematism to crypsis in a neotropical toxic-plant-feeding Romaleid grasshopper, Chromacris psittacus (Gerstaecker, 1873). Field surveys, behavioural observations and a rearing experiment compare host plant choice, aggregation, locomotion and thermoregulation between life history stages. Results showed that both nymphs and adults fed exclusively on a narrow range of Solanaceae plants, suggesting that the shift in defensive syndrome is not due to a change in chemical defense. Instead, nymphal aposematism appears linked to aggregation in response to plant-based selection pressures. Slow nymphal development suggests a cost to feeding on toxic plant compounds, and grouping could mitigate this cost. Grouping also increases conspicuousness, and hence can favour warning colourating in chemically-defended insects. The role of diet breadth in aposematism is poorly understood, and these results suggest how constraints imposed by feeding on toxic plants can generate bottom-up selection pressures shaping the adaptive suites of traits of chemically-defended animals.

There is no magic bullet: the importance of testing reference gene stability in RT-qPCR experiments across multiple closely related species

Reverse Transcriptase quantitative Polymerase Chain Reaction (RT-qPCR) is the current gold standard tool for the study of gene expression. This technique is highly dependent on the validation of reference genes, which exhibit stable expression levels among experimental conditions. Often, reference genes are assumed to be stable a priori without a rigorous test of gene stability. However, such an oversight can easily lead to misinterpreting expression levels of target genes if the references genes are in fact not stable across experimental conditions. Even though most gene expression studies focus on just one species, comparative studies of gene expression among closely related species can be very informative from an evolutionary perspective. In our study, we have attempted to find stable reference genes for four closely related species of grasshoppers (Orthoptera: Acrididae) that together exhibit a spectrum of density-dependent phenotypic plasticity. Gene stability was assessed for eight reference genes in two tissues, two experimental conditions and all four species. We observed clear differences in the stability ranking of these reference genes, both between tissues and between species. Additionally, the choice of reference genes clearly influenced the results of a gene expression experiment. We offer suggestions for the use of reference genes in further studies using these four species, which should be taken as a cautionary tale for future studies involving RT-qPCR in a comparative framework.

Extra-genomic instructive influences in morphogenesis: A review of external signals that regulate growth and form

Embryonic development and regeneration accomplish a remarkable feat: individual cells work together to create or repair complex anatomical structures. What is the source of the instructive signals that specify these invariant and robust organ-level outcomes? The most frequently studied source of morphogenetic control is the host genome and its transcriptional circuits. However, it is now apparent that significant information affecting patterning also arrives from outside of the body. Both biotic and physical factors, including temperature and various molecular signals emanating from pathogens, commensals, and conspecific organisms, affect developmental outcomes. Here, we review examples in which anatomical patterning decisions are strongly impacted by lateral signals that originate from outside of the zygotic genome. The endogenous pathways targeted by these influences often show transgenerational effects, enabling them to shape the evolution of anatomies even faster than traditional Baldwin-type assimilation. We also discuss recent advances in the biophysics of morphogenetic controls and speculate on additional sources of important patterning information which could be exploited to better understand the evolution of bodies and to design novel approaches for regenerative medicine.

Eight Million Years of Satellite DNA Evolution in Grasshoppers of the Genus Schistocerca Illuminate the Ins and Outs of the Library Hypothesis

Satellite DNA (satDNA) is an abundant class of tandemly repeated noncoding sequences, showing high rate of change in sequence, abundance, and physical location. However, the mechanisms promoting these changes are still controversial. The library model was put forward to explain the conservation of some satDNAs for long periods, predicting that related species share a common collection of satDNAs, which mostly experience quantitative changes. Here, we tested the library model by analyzing three satDNAs in ten species of Schistocerca grasshoppers. This group represents a valuable material because it diversified during the last 7.9 Myr across the American continent from the African desert locust (Schistocerca gregaria), and this thus illuminates the direction of evolutionary changes. By combining bioinformatic and cytogenetic, we tested whether these three satDNA families found in S. gregaria are also present in nine American species, and whether differential gains and/or losses have occurred in the lineages. We found that the three satDNAs are present in all species but display remarkable interspecies differences in their abundance and sequences while being highly consistent with genus phylogeny. The number of chromosomal loci where satDNA is present was also consistent with phylogeny for two satDNA families but not for the other. Our results suggest eminently chance events for satDNA evolution. Several evolutionary trends clearly imply either massive amplifications or contractions, thus closely fitting the library model prediction that changes are mostly quantitative. Finally, we found that satDNA amplifications or contractions may influence the evolution of monomer consensus sequences and by chance playing a major role in driftlike dynamics.

Gene functions in adult cuticle pigmentation of the yellow mealworm, Tenebrio molitor

In many arthropod species including insects, the cuticle tanning pathway for both pigmentation and sclerotization begins with tyrosine and is responsible for production of both melanin- and quinoid-type pigments, some of which are major pigments for body coloration. In this study we identified and cloned cDNAs of the yellow mealworm, Tenebrio molitor, encoding seven key enzymes involved in this pathway including tyrosine hydroxylase (TmTH), DOPA decarboxylase (TmDDC), laccase 2 (TmLac2), Yellow-y (TmY-y), arylalkylamine N-acetyltransferase (TmAANAT1), aspartate 1-decarboxylase (TmADC) and N-β-alanyldopamine synthase (Tmebony). Expression profiles of these genes during development were analyzed by real-time PCR, revealing development-specific patterns of expression. Loss of function mediated by RNAi of either 1) TmTH or TmLac2, 2) TmDDC or TmY-y, and 3) TmAANAT1, TmADC or Tmebony resulted in pale/white, light yellow/brown and dark/black adult body coloration, respectively. In addition, there are three distinct layer/regional pigmentation differences in rigid types of adult cuticle, a brownish outer exocuticle (EX), a dark pigmented middle mesocuticle (ME) and a transparent inner endocuticle (EN). Decreases in pigmentation of the EX and/or ME layers were observed after RNAi of TmDDC or TmY-y. In TmADC- or Tmebony-deficient adults, a darker pigmented EX layer was observed. In TmAANAT1-deficient adults, trabeculae formed between the dorsal and ventral elytral cuticles as well as the transparent EN layer became highly pigmented. These results demonstrate that knocking down the level of gene expression of specific enzymes of this tyrosine metabolic pathway leads to abnormal pigmentation in individual layers and substructure of the rigid adult exoskeleton of T. molitor.

Revealing hidden density-dependent phenotypic plasticity in sedentary grasshoppers in the genus Schistocerca Stål (Orthoptera: Acrididae: Cyrtacanthacridinae)

Comparative quantification of reaction norms across closely related species in a clade is rare, but such a study can reveal valuable insights into understanding how reaction norms evolve along phylogeny. The grasshopper genus Schistocerca Stål (Orthoptera: Acrididae: Cyrtacanthacridinae) is an ideal group to study the evolution of density-dependent phenotypic plasticity because it includes both swarming locusts and non-swarming sedentary grasshoppers, which show varying degrees of plastic reaction norms in many traits. The swarming locusts exhibit locust phase polyphenism in which cryptically colored and solitary individuals can transform into conspicuously colored and highly gregarious individuals in response to increases in population density. The sedentary grasshoppers do not swarm in nature, and thus it has been assumed that they have little or no expression of plastic reaction norms in many traits, except for color, which has been shown to be a phylogenetically conserved trait. In this study, we have quantified density-dependent reaction norms in behavior, color, body size, and morphometric ratio in the nymphs of four sedentary species within Schistocerca by conducting explicit rearing experiments to induce potential phenotypic changes in response to isolation and crowding. In contrast to our previous assumption, we find that all four species show a certain level of density-dependent plastic reaction norms, which implies that these sedentary species have hidden reaction norms that can only be induced experimentally, some components of which must be phylogenetically conserved. Furthermore, we demonstrate that rearing density differentially affects the expression of reaction norms in different species, suggesting that different reaction norms must have followed independent evolutionary trajectories.

The foraging Gene and Its Behavioral Effects: Pleiotropy and Plasticity

The Drosophila melanogaster foraging (for) gene is a well-established example of a gene with major effects on behavior and natural variation. This gene is best known for underlying the behavioral strategies of rover and sitter foraging larvae, having been mapped and named for this phenotype. Nevertheless, in the last three decades an extensive array of studies describing for's role as a modifier of behavior in a wide range of phenotypes, in both Drosophila and other organisms, has emerged. Furthermore, recent work reveals new insights into the genetic and molecular underpinnings of how for affects these phenotypes. In this article, we discuss the history of the for gene and its role in natural variation in behavior, plasticity, and behavioral pleiotropy, with special attention to recent findings on the molecular structure and transcriptional regulation of this gene.

Genomics of Developmental Plasticity in Animals

Developmental plasticity refers to the property by which the same genotype produces distinct phenotypes depending on the environmental conditions under which development takes place. By allowing organisms to produce phenotypes adjusted to the conditions that adults will experience, developmental plasticity can provide the means to cope with environmental heterogeneity. Developmental plasticity can be adaptive and its evolution can be shaped by natural selection. It has also been suggested that developmental plasticity can facilitate adaptation and promote diversification. Here, we summarize current knowledge on the evolution of plasticity and on the impact of plasticity on adaptive evolution, and we identify recent advances and important open questions about the genomics of developmental plasticity in animals. We give special attention to studies using transcriptomics to identify genes whose expression changes across developmental environments and studies using genetic mapping to identify loci that contribute to variation in plasticity and can fuel its evolution.

Satellite DNAs Unveil Clues about the Ancestry and Composition of B Chromosomes in Three Grasshopper Species

Supernumerary (B) chromosomes are dispensable genomic elements occurring frequently among grasshoppers. Most B chromosomes are enriched with repetitive DNAs, including satellite DNAs (satDNAs) that could be implicated in their evolution. Although studied in some species, the specific ancestry of B chromosomes is difficult to ascertain and it was determined in only a few examples. Here we used bioinformatics and cytogenetics to characterize the composition and putative ancestry of B chromosomes in three grasshopper species, Rhammatocerus brasiliensis, Schistocerca rubiginosa, and Xyleus discoideus angulatus. Using the RepeatExplorer pipeline we searched for the most abundant satDNAs in Illumina sequenced reads, and then we generated probes used in fluorescent in situ hybridization (FISH) to determine chromosomal position. We used this information to infer ancestry and the events that likely occurred at the origin of B chromosomes. We found twelve, nine, and eighteen satDNA families in the genomes of R. brasiliensis, S. rubiginosa, and X. d. angulatus, respectively. Some satDNAs revealed clustered organization on A and B chromosomes varying in number of sites and position along chromosomes. We did not find specific satDNA occurring in the B chromosome. The satDNAs shared among A and B chromosomes support the idea of putative intraspecific ancestry from small autosomes in the three species, i.e., pair S11 in R. brasiliensis, pair S9 in S. rubiginosa, and pair S10 in X. d. angulatus. The possibility of involvement of other chromosomal pairs in B chromosome origin is also hypothesized. Finally, we discussed particular aspects in composition, origin, and evolution of the B chromosome for each species.

Epigenetics and developmental plasticity in orthopteroid insects

Developmental plasticity is a key driver of the extraordinary ecological success of insects. Epigenetic mechanisms provide an important link between the external stimuli that initiate polyphenisms, and the stable changes in gene expression that govern alternative insect morphs. We review the epigenetics of orthopteroid insects, focussing on recent research on locusts and termites, two groups which display high levels of phenotypic plasticity, and for which genome sequences have become available in recent years. We examine research on the potential role of DNA methylation, histone modifications, and non-coding RNAs in the regulation of gene expression in these insects. DNA methylation patterns in orthopteroids share a number of characteristics with those of hymenopteran insects, although methylation levels are much higher, and extend to introns and repeat elements. Future examinations of epigenetic mechanisms in these insects will benefit from comparison of tissues from aged-matched individuals from alternative morphs, and adequate biological replication.