- Research article
- Open Access
Nymphalid eyespots are co-opted to novel wing locations following a similar pattern in independent lineages
- Sandra R Schachat1, 2Email author,
- Jeffrey C Oliver3 and
- Antónia Monteiro4, 5, 6Email author
https://doi.org/10.1186/s12862-015-0300-x
© Schachat et al.; licensee BioMed Central. 2015
- Received: 10 September 2014
- Accepted: 29 January 2015
- Published: 14 February 2015
Abstract
Background
Variation in the number of repeated traits, or serial homologs, has contributed greatly to animal body plan diversity. Eyespot color patterns of nymphalid butterflies, like arthropod and vertebrate limbs, are an example of serial homologs. These eyespot color patterns originated in a small number of wing sectors on the ventral hindwing surface and later appeared in novel wing sectors, novel wings, and novel wing surfaces. However, the details of how eyespots were co-opted to these novel wing locations are currently unknown.
Results
We used a large data matrix of eyespot/presence absence data, previously assembled from photographs of contemporary species, to perform a phylogenetic investigation of eyespot origins in nine independent nymphalid lineages. To determine how the eyespot gene regulatory network acquired novel positional information, we used phylogenetic correlation analyses to test for non-independence in the origination of eyespots. We found consistent patterns of eyespot gene network redeployment in the nine lineages, where eyespots first redeployed from the ventral hindwing to the ventral forewing, then to new sectors within the ventral wing surface, and finally to the dorsal wing surface. Eyespots that appeared in novel wing sectors modified the positional information of their serial homolog ancestors in one of two ways: by changing the wing or surface identity while retaining sector identity, or by changing the sector identity while retaining wing and surface identity.
Conclusions
Eyespot redeployment to novel sectors, wings, and surfaces happened multiple times in different nymphalid subfamilies following a similar pattern. This indicates that parallel mutations altering expression of the eyespot gene regulatory network led to its co-option to novel wing locations over time.
Keywords
- Serial homology
- Correlation analysis
- Ancestral states
- Phylogeny
- Wing patterns
Background
Animal bodies are often made up of repeated units, or serial homologs, which vary in number, size, shape, and position and largely contribute to body plan diversity. Examples of serial homologs include vertebrate limbs [1] and teeth [2]; insect limbs [3], bristles [4,5], and trichomes [6,7]; and nymphalid butterfly eyespots [8]. In many animal lineages, the development of serial homologs and other complex traits is mediated via gene regulatory networks [9]. Serial homologs develop using a largely shared network [10]. However, in many of the examples above, it is still unclear how the gene regulatory network, responsible for the development of each serial homolog, became deployed at each body location and gained individuation from the other units to generate body plan diversity.
Examples of eyespot diversity in nymphalid butterflies. (A) Hamanumida daedalus (YPM ENT 422238). (B) Hamadryas februa (YPM ENT 815633). (C) Marpesia eleuchea (YPM ENT 815891). (D) Neorina crishna (YPM ENT 407679). (E) Bicyclus anynana (YPM ENT 814497). (F) Megisto cymela (YPM ENT 405705). All photographs show the ventral wing surface of female butterflies; all specimens represent taxa investigated in this study and are housed in the Entomology Collection of the Yale Peabody Museum of Natural History.
Here we take a phylogenetic approach to examine how eyespots evolved at novel wing locations. We used correlation analyses to elucidate patterns of independence and non-independence across pairs of eyespots, and therefore, the patterns of eyespot deployment on the wing. We used ancestral state reconstructions to evaluate an iterative series of eyespot evolution models, allowing tests of potential pathways through which eyespots may be deployed across different wing locations. These models provide a test for whether the origin of an eyespot in one wing location is contingent on the presence of another eyespot in a different wing location. Comparison among models of evolutionary change thus allow us to evaluate support for the possibility that regulatory information coding for the more recent (i.e. “derived”) eyespot evolved from regulatory information that coded for the older (i.e. “ancestral”) eyespot. In contrast, those cases in which eyespot origins are not contingent on the presence of another eyespot on the wings may indicate that independent enhancers, or regulatory sequences, have evolved for each eyespot. By examining the evolutionary history of multiple eyespots, our goal is to understand the relationship between pre-existing positional information for extant eyespots and novel positional information required to deploy eyespots to new locations. Furthermore, by evaluating these models in multiple lineages of nymphalid butterflies we can assess the degree to which there are shared pathways by which novel positional information evolves to generate diversity in butterfly eyespot number and location.
Methods
Data collection and ancestral state reconstructions
We collected character data from 394 nymphalid species based on phylogenetic sampling of previous studies [16,17]. The majority of specimens are from holdings in the Yale Peabody Natural History Museum, the American Museum of Natural History, and the Harvard Museum of Comparative Zoology. For each species, we scored eyespot presence/absence for all wing surfaces (the dorsal and ventral surfaces of the forewing and hindwing), for a grand total of 38 different wing sectors for a single side of the body (Additional file 1). We scored eyespots for each of the sexes, and used two specimens per sex when possible. An eyespot is defined as a circular pattern element containing two or more concentric rings of color [8]. Additional pattern elements, such as patches and single-color spots, were not used in these analyses. All patterns were scored from photographs available at www.lepdata.org. We used the estimate of nymphalid phylogeny published by Wahlberg et al. [17].
Identification of independent clades for analysis
A phylogeny of the Nymphalidae showing all tribes (and subtribes within the tribe Satyrini), modified from [ 17 ]. The clades analyzed in this study are marked by an eyespot; because Parargina, Mycalesina, and Lethina form a clade, they were analyzed together, as indicated by the eyespot at the node shared by these three subtribes. Colored branches represent different subfamilies.
Eyespot models of evolution
For each clade, we restricted our analyses to those pairs of eyespots that were inferred to have a shared history in the clade being investigated. A shared history was inferred when both eyespots had a higher likelihood of being present than absent in at least one shared ancestral node. Ancestral state likelihoods were calculated in Mesquite [18], using an asymmetrical (2-parameter) likelihood model. For each pair of eyespots identified as having a shared history, we categorized one eyespot as “ancestral” (the eyespot that evolved first) and the other eyespot as “derived”. To objectively categorize the eyespots, we again used the likelihoods of presence/absence for each eyespot using an asymmetrical likelihood model in Mesquite [18]. If one eyespot had a greater likelihood of being present at a node ancestral to the “shared node” identified above, while the other eyespot had a greater likelihood of being absent at that ancestral node, the first eyespot was categorized as the ancestral eyespot. When the deepest node where each eyespot was more likely present than absent was the same for each eyespot, the eyespot with the greater likelihood of being present was categorized as the ancestral eyespot.
Four models of eyespot evolution used to identify correlated eyespots. (A) Independent constrained model, in which transition rates in an eyespot are independent of the state of the other eyespot and the origination rates for the two eyespots are set to equal; (B) independent unconstrained model, in which transition rates in an eyespot are independent of the state of the other eyespot; (C) dependent constrained model, in which the transition rates of the derived eyespot are conditional on the state of the ancestral eyespot and in which the derived eyespot may only evolve when ancestral eyespot is present (the transition rate indicated by white arrow is set to zero). (D) dependent unconstrained model, in which the derived eyespot may originate regardless of whether the ancestral eyespot is present or absent. Legend: A−: ancestral eyespot absent; A+: ancestral eyespot present; D−: derived eyespot absent; D+: derived eyespot present. In each panel, colors of arrows indicate unique rate parameters. Models increase in complexity, i.e., in the number of different rates of transition between character states (arrows of different colors) that have to be estimated from the data from left (model A; three rates) to right (model D; 6 rates).
Results
Eyespot models of evolution
The eyespots of extant Nymphalidae have three kinds of positional information: wing (forewing vs. hindwing), wing surface (dorsal vs. ventral), and wing sector. Because the first eyespots arose on only one wing, on only one wing surface, and only in certain wing sectors [13], these three positional identities were present when Nymphalid eyespots first arose. In this study, eyespots are assumed to redeploy to new locations by acquiring only one new kind of positional information at a time. Therefore, the eyespot pairs used for analyses were those that differ only in wing identity (same wing sector, same wing surface), those that differ only in wing surface identity (same wing sector, same wing), and those that differ only in wing sector identity (same wing, same wing surface).
Illustration of all the wing sectors where at least one eyespot was found in species scored for this study, and correspondent eyespot sector nomenclature. The five eyespots identified by [13] as having appeared first are outlined in red. The ventral surface is rendered in light brown, on the left; the dorsal surface is rendered in dark brown, on the right.
Significant correlations between pairs of eyespots identified in this study, organized by derived eyespot. The arrow points from the ancestral eyespot to the derived eyespot, which is colored yellow. Ventral surface is light brown; dorsal surface is dark brown.
Patterns of eyespot redeployment
Different patterns of eyespot redeployment: (A) between wings and surfaces (eyespots M3 and Cu1), (B) within wing surface (eyespots ScR, Cu2, and Pc), and (C) ambiguous (eyespots Rs/R5, M1, and M2). The arrow points from the ancestral eyespots to the derived eyespots. The derived eyespots in each schematic are colored yellow. Ventral surface is light brown; dorsal surface is dark brown.
Eyespots in sectors M3 and Cu1 (Figure 4) belong to the first two categories (Figure 6A). Previous work estimated that the first eyespots appeared only on the hindwing [13], and if we are correct in our assumption that eyespot redeployment occurs only by step-wise changes in positional information, then the M3 and Cu1 eyespots on the ventral forewing must have originated as a modification of the positional information coding for the M3 and Cu1 eyespots on the ventral hindwing. Positive correlation data (Additional file 5: Table S3) also supports the inference that M3 and Cu1 eyespots originated on the dorsal surface from modification of positional information coding for their ventral serial homologs.
Eyespots in sectors ScR, Cu2, and Pc belong to the third category, where eyespots are redeployed to different wing sectors within a single wing surface (Figure 6B): our results indicate that these eyespots originate from modifications of positional information coding for eyespots present in other sectors on the same wing.
The origin of eyespots R5/Rs, M1, and M2 remains ambiguous, as we recovered support for models that fit in the first, second, and third categories. At present we cannot determine whether these eyespots originate as modification of positional information coding for eyespots in other wings, wing surfaces, or from the M3 and Cu1 eyespots within the same wing surface (Figure 6C).
Within wing surfaces, eyespots do not appear sequentially in adjacent wing sectors. For example, the origin of eyespots in the Pc wing sector is never contingent on the presence of eyespots in the adjacent Cu2 wing sector.
None of the significant correlations that we recovered suggest any differences in patterns of ventral eyespot redeployment among nymphalid subfamilies, tribes, and subtribes (Additional file 4: Tables S2, Additional file 5: Table S3, Additional file 6: Table S4). The few cases of non-independence in dorsal eyespot origins are also consistent among all lineages examined in that they show that, for both wings, only eyespots in the sectors Rs, M1, M2, M3, and Cu1 redeployed from the ventral to the dorsal wing surface.
Discussion
Here we document how nymphalid butterfly eyespot number evolved over the past ~90 million years. A recent study showed that eyespots originated in a few wing sectors on the hindwing [13]. Here we extend that study by showing that these eyespots gradually redeployed to new wings (forewing), new surfaces (the dorsal surface), and new wing sectors to produce the extant Nymphalidae eyespot diversity. In addition, these redeployments appear to have occurred in a parallel fashion across independent lineages.
We found that eyespots belonging to different wing sectors (ScR, Cu2, or Pc as opposed to M3 or Cu1) follow different patterns of gene network redeployment. ScR, Cu2, and Pc likely arise via modification of regulatory positional information that was already available for eyespots on the same wing surface, whereas M3 and Cu1 appear to arise via modification of regulatory information available for homologous eyespots on other wing surfaces. Therefore, although eyespots can appear in a continuous row within a single wing surface, they evolved positional information for several of these sectors in different ways.
Similarities in redeployment across the phylogeny
Frequency of eyespot presence in each wing sector for each Nymphalid subfamily. The size of each eyespot corresponds to the percentage of taxa sampled with an eyespot present in females in the given wing sector. Eyespots that are never present in any of the taxa sampled are not shown; no eyespots are present in the families Calinaginae and Libytheinae. The width of each branch is proportional to the number of taxa sampled per subfamily. Ventral surface is light brown; dorsal surface is dark brown.
Developmental hypotheses for eyespot redeployments
The developmental mechanisms that underlie eyespot redeployment are not currently understood. Molecular work in model species has illustrated that alterations to enhancers of key genes are often responsible for the deployment of serial homologs in different places in the body, and to evolution of serial homolog number [20]. Those key genes control the activation of the whole gene network that regulates the development of each serial homolog [20]. Two types of modification to the enhancers have been described in various Drosophila lineages. The first includes elimination of distinct enhancers of a key regulator of trichome development, shavenbaby, leading to loss of trichomes in specific areas of the larval body [6,7]. The second involves modifications to pre-existent single enhancers of scute that extended this gene’s expression domain to more anterior body locations, and increased the number of bristles in those areas [5,7]. These alternative modes for variation in serial homolog number can lead to different patterns of correlations regarding gains and losses of serial homologs. The first mode predicts that gains and losses of serial homologs should not be evolutionary correlated with each other, as the enhancers were completely independent of each other, whereas gains and losses of serial homologs in the second mode predicts some form of correlation across serial homologs, i.e., gains (or losses) of new homologs only happen contingent on the presence of a pre-existent serial homolog. It is possible that similar modification to enhancers of eyespot network genes led to the different patterns of redeployment observed in our study.
When eyespots first originated it is possible that the activating signals for network deployment may have been present across every single wing sector. These signals could have involved morphogens diffusing from the veins and from the wing margin as previously proposed [11]. However, it also possible that negative regulators of the network, present in particular wings, surfaces, and/or sectors may have initially prevented the network from being expressed in those areas. This would be the case if binding sites for these negative regulators were initially present in eyespot network genes, at the time of eyespot origins. Over time, however, the removal of these negative regulatory interactions may have allowed the network to become active in additional wings (forewings), surfaces (dorsal surface), and sectors.
Adjacent sectors on an insect wing look developmentally homogeneous, but in fact, they are not. During wing development in Drosophila, each wing sector expresses different combinations of transcription factors that gives that sector a unique genetic identity [21,22]. The same applies to dorsal and ventral surfaces and fore and hindwings that express different selector genes [23-25]. For example, the transcription factor Spalt, expressed in the center of the Drosophila wing, is also expressed in the butterfly wing, in a sector-specific pattern [26]. Regulatory evolution of the genes in the eyespot network (losing or gaining binding sites for genes such as Spalt, as well as for other wing sector-restricted genes, activators or repressors) may have produced the pattern of successive eyespot origins uncovered in this study. Future studies should examine the regulatory composition of transcription factors expressed in each sector of a nymphalid wing. This information could help us make sense of the particular history of eyespot origins as well as help in the design of experiments to understand the molecular basis of eyespot network deployment.
Comparison with evolution of serial homology in other systems
Butterfly eyespots are a relatively recent example of serial homology in arthropods. Other serial homologous traits with earlier origins, such as arthropod limbs and body segments, vary far less than eyespots: segment number is stable in most high-rank clades of arthropods, including insects [27], and arthropod limb number has largely remained stable during the past 250 million years [28]. Butterfly eyespot number, on the other hand, continues to evolve and to vary among closely related taxa, and this recent variation may yield clues to its genetic basis.
The ancestral condition for eyespot distribution on the wings also contrasts sharply with other, closely studied serially homologous traits such as arthropod limbs and vertebrate teeth, which are more likely to have arisen as multiple repeated units that later become reduced in number and individuated [2,3]. The evolution of eyespot number appears to be most similar to the evolution of the paired vertebrate limbs, where fins (limbs) appear first in the anterior part of the body (pectoral fins) and later appear in a more posterior region (pelvic fins) [1].
Other serially homologous traits that vary among closely related extant insects include the bristles and trichomes in Drosophila. In these species, the genetic basis of variation across closely related species is better understood. Sensory bristles originated in four rows on the thorax in the ancestor of Drosophila and other higher Diptera [29]. The lineage leading to Drosophila lost some anterior bristles but a few members of the clade subsequently re-evolved them [5,30]. This secondary gain of bristles was due to cis-regulatory evolution at the scute locus [5,31]. The recent and independent loss of trichomes, hair-like projections that appear on the larval body in multiple lineages of Drosophila, was also due to cis-regulatory evolution at a single locus, shavenbaby [7,32]. The evolution of trichomes and their genetic underpinnings has been studied in a phylogenetic context, and ancestral patterns of trichome distribution have been proposed [33], but ancestral trichome distributions on different body locations have not yet been reconstructed.
Nymphalid eyespots represent the reverse of the scenario found for Drosophila bristles, in that bristles originally appeared in four complete rows, whereas eyespots arose in a few wing sectors on the ventral hindwing and later originated in new locations, on the forewing and on the dorsal wing surface. Continuing advances in the systematics of higher Diptera [34], stronger hypotheses for ancestral trichome distributions, and identification of the genetic underpinnings of butterfly eyespot development [14] will facilitate the integration of developmental and phylogenetic data to produce a complete understanding of the evolution of serial homology through both loss and gain (redeployment) of individuated traits.
Conclusions
In butterflies of the family Nymphalidae, eyespots first arose in only a few wing sectors on the ventral wing surface. Over tens of millions of years, the eyespot gene network redeployed to new wing locations; this happened independently in multiple lineages. However, our results show that eyespot gene network redeployment followed the same pattern in all subfamilies. The first eyespots, located in five wing sectors on the ventral hindwing, redeployed to the ventral forewing. From there, these ventral eyespots redeployed to novel locations within the same wing surface, and to the dorsal surface.
Availability of supporting data
The data sets supporting the results of this article are included within the article (and its additional files).
Declarations
Acknowledgements
We thank Donald J Harvey, Carla M Penz, and Jeffrey W Shultz for helpful suggestions and discussion, and Lawrence F Gall for assistance in the entomology collection of the Peabody Museum of Natural History. Numerous undergraduates assisted with photography and initial wing-scoring. NSF IOS 0818731 and Ministry of Education of Singapore R-154-000-602-112 grants supported this research. The authors declare no conflicts of interest.
Authors’ Affiliations
References
- Ruvinsky I, Gibson-Brown JJ. Genetic and developmental bases of serial homology in vertebrate limb evolution. Development. 2000;127:5233–44.PubMedGoogle Scholar
- Stock DW, Jackman WR, Trapani J. Developmental genetic mechanisms of evolutionary tooth loss in cypriniform fishes. Development. 2006;133:3127–37.View ArticlePubMedGoogle Scholar
- Gebelein B, Culi J, Ryoo HD, Zhang W, Mann RS. Specificity of Distalless repression and limb primordia development by abdominal Hox proteins. Dev Cell. 2002;3:487–98.View ArticlePubMedGoogle Scholar
- Usui K, Goldstone C, Gibert J-M, Simpson P. Redundant mechanisms mediate bristle patterning on the Drosophila thorax. Proc Natl Acad Sci U S A. 2008;105:20112–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Marcellini S, Simpson P. Two or four bristles: functional evolution of an enhancer of scute in Drosophilidae. PLoS Biol. 2006;4:e386.View ArticlePubMed CentralPubMedGoogle Scholar
- McGregor AP, Orgogozo V, Delon I, Zanet J, Srinivasan DG, Payre F, et al. Morphological evolution through multiple cis-regulatory mutations at a single gene. Nature. 2007;448:587–90.View ArticlePubMedGoogle Scholar
- Sucena É, Delon I, Jones I, Payre F, Stern DL. Regulatory evolution of shavenbaby/ovo underlies multiple cases of morphological parallelism. Nature. 2003;424:935–8.View ArticlePubMedGoogle Scholar
- Monteiro A. Alternative models for the evolution of eyespots and of serial homology on lepidopteran wings. BioEssays. 2008;30:358–66.View ArticlePubMedGoogle Scholar
- Stathopoulos A, Levine M. Genomic regulatory networks and animal development. Dev Cell. 2005;9:449–62.View ArticlePubMedGoogle Scholar
- Gilbert SF, Bolker JA. Homologies of process and modular elements of embryonic construction. J Exp Zool. 2001;291:1–12.View ArticlePubMedGoogle Scholar
- Nijhout HF. Elements of butterfly wing patterns. J Exp Zool. 2001;291:213–25.View ArticlePubMedGoogle Scholar
- Held LI. Rethinking butterfly eyespots. Evol Biol. 2012;40:158–68.View ArticleGoogle Scholar
- Oliver JC, Beaulieu JM, Gall LF, Piel WH, Monteiro A. Nymphalid eyespot serial homologues originate as a few individualized modules. Proc Biol Sci. 2014;281:20133262.View ArticlePubMedGoogle Scholar
- Oliver JC, Tong X, Gall LF, Piel WH, Monteiro A. A single origin for nymphalid butterfly eyespots followed by widespread loss of associated gene expression. PLoS Genet. 2012;8:e1002893.View ArticlePubMed CentralPubMedGoogle Scholar
- Force A, Cresko WA, Pickett FB, Proulx SR, Amemiya C, Lynch M. The origin of Subfunctions and modular gene regulation. Genetics. 2005;170:433–46.View ArticlePubMed CentralPubMedGoogle Scholar
- Oliver J, Beaulieu J, Gall L, Piel W, Monteiro A: Data from: Nymphalid eyespot serial homologs originate as a few individualized modules. Dryad Digit Repos 2014. http://dx.doi.org/10.5061/dryad.f55c5
- Wahlberg N, Leneveu J, Kodandaramaiah U, Peña C, Nylin S, Freitas AVL, et al. Nymphalid butterflies diversify following near demise at the Cretaceous/Tertiary boundary. Proc R Soc B. 2009;276:4295–302.View ArticlePubMed CentralPubMedGoogle Scholar
- Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis, version 2.75. [http://mesquiteproject.org] 2013.
- Pagel M, Meade A: BayesTraits, version 1.0. [http://www.evolution.rdg.ac.uk/BayesTraits.html] 2013.
- Arnosti DN. Analysis and function of transcriptional regulatory elements: insights from Drosophila. Annu Rev Entomol. 2003;48:579–602.View ArticlePubMedGoogle Scholar
- Blair SS. Wing vein patterning in Drosophila and the analysis of intercellular signaling. Annu Rev Cell Dev Biol. 2007;23:293–319.View ArticlePubMedGoogle Scholar
- Monteiro A, Prijs J, Bax M, Hakkaart T, Brakefield PM. Mutants highlight the modular control of butterfly eyespot patterns. Evol Dev. 2003;187:180–7.View ArticleGoogle Scholar
- Carroll SB, Gates J, Keys DN, Paddock SW, Panganiban GE, Selegue JE, et al. Pattern formation and eyespot determination in butterfly wings. Science. 1994;265(80):109–14.View ArticlePubMedGoogle Scholar
- Weatherbee SD, Nijhout HF, Grunert LW, Halder G, Galant R, Selegue JE, et al. Ultrabithorax function in butterfly wings and the evolution of insect wing patterns. Curr Biol. 1999;9:109–15.View ArticlePubMedGoogle Scholar
- Tong X, Hrycaj S, Podlaha O, Popadic A, Monteiro A. Over-expression of Ultrabithorax alters embryonic body plan and wing patterns in the butterfly Bicyclus anynana. Dev Biol. 2014;394:1–10.View ArticleGoogle Scholar
- Monteiro A. The origin, development, and evolution of butterfly eyespots. Annu Rev Entomol. 2015;60:253–71.View ArticlePubMedGoogle Scholar
- Minelli A, Maruzzo D, Fusco G. Multi-scale relationships between numbers and size in the evolution of arthropod body features. Arthropod Struct Dev. 2010;39:468–77.View ArticlePubMedGoogle Scholar
- McShea DW. Perspective: metazoan complexity and evolution: is there a trend? Evol (N Y). 1996;50:477.Google Scholar
- Culi J, Modolell J. Proneural gene self-stimulation in neural precursors: an essential mechanism for sense organ development that is regulated by notch signaling. Genes Dev. 1998;12:2036–47.View ArticlePubMed CentralPubMedGoogle Scholar
- Pistillo D, Skaer N, Simpson P. scute expression in Calliphora vicina reveals an ancestral pattern of longitudinal stripes on the thorax of higher Diptera. Development. 2002;129:563–72.PubMedGoogle Scholar
- Stern DL, Orgogozo V. The loci of evolution: how predictable is genetic evolution? Evol (N Y). 2008;62:2155–77.Google Scholar
- Stern DL, Frankel N. The structure and evolution of cis-regulatory regions: the shavenbaby story. Philos Trans R Soc Lond B Biol Sci. 2013;368:20130028.View ArticlePubMed CentralPubMedGoogle Scholar
- Arif S, Murat S, Almudi I, Nunes MDS, Bortolamiol-Becet D, McGregor NS, et al. Evolution of mir-92a underlies natural morphological variation in Drosophila melanogaster. Curr Biol. 2013;23:523–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Lambkin CL, Sinclair BJ, Pape T, Courtney GW, Skevington JH, Meier R, et al. The phylogenetic relationships among infraorders and superfamilies of Diptera based on morphological evidence. Syst Entomol. 2013;38:164–79.View ArticleGoogle Scholar
Copyright
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.