Phylogenomic analysis reveals bees and wasps (Hymenoptera) at the base of the radiation of Holometabolous insects

Different phylogenomic approaches to resolve the evolutionary relationships among model fish species

Comparative genomics holds the promise to magnify the information obtained from individual genome sequencing projects, revealing common features conserved across genomes and identifying lineage-specific characteristics. To implement such a comparative approach, a robust phylogenetic framework is required to accurately reconstruct evolution at the genome level. Among vertebrate taxa, teleosts represent the second best characterized group, with high-quality draft genome sequences for five model species (Danio rerio, Gasterosteus aculeatus, Oryzias latipes, Takifugu rubripes, and Tetraodon nigroviridis), and several others are in the finishing lane. However, the relationships among the acanthomorph teleost model fishes remain an unresolved taxonomic issue. Here, a genomic region spanning over 1.2 million base pairs was sequenced in the teleost fish Dicentrarchus labrax. Together with genomic data available for the above fish models, the new sequence was used to identify unique orthologous genomic regions shared across all target taxa. Different strategies were applied to produce robust multiple gene and genomic alignments spanning from 11,802 to 186,474 amino acid/nucleotide positions. Ten data sets were analyzed according to Bayesian inference, maximum likelihood, maximum parsimony, and neighbor joining methods. Extensive analyses were performed to explore the influence of several factors (e.g., alignment methodology, substitution model, data set partitions, and long-branch attraction) on the tree topology. Although a general consensus was observed for a closer relationship between G. aculeatus (Gasterosteidae) and Di. labrax (Moronidae) with the atherinomorph O. latipes (Beloniformes) sister taxon of this clade, with the tetraodontiform group Ta. rubripes and Te. nigroviridis (Tetraodontiformes) representing a more distantly related taxon among acanthomorph model fish species, conflicting results were obtained between data sets and methods, especially with respect to the choice of alignment methodology applied to noncoding parts of the genomic region under study. This may limit the use of intergenic/noncoding sequences in phylogenomics until more robust alignment algorithms are developed.

Comparative mitogenomics of Braconidae (Insecta: Hymenoptera) and the phylogenetic utility of mitochondrial genomes with special reference to Holometabolous insects

BACKGROUND

Animal mitochondrial genomes are potential models for molecular evolution and markers for phylogenetic and population studies. Previous research has shown interesting features in hymenopteran mitochondrial genomes. Here, we conducted a comparative study of mitochondrial genomes of the family Braconidae, one of the largest families of Hymenoptera, and assessed the utility of mitochondrial genomic data for phylogenetic inference at three different hierarchical levels, i.e., Braconidae, Hymenoptera, and Holometabola.

RESULTS

Seven mitochondrial genomes from seven subfamilies of Braconidae were sequenced. Three of the four sequenced A+T-rich regions are shown to be inverted. Furthermore, all species showed reversal of strand asymmetry, suggesting that inversion of the A+T-rich region might be a synapomorphy of the Braconidae. Gene rearrangement events occurred in all braconid species, but gene rearrangement rates were not taxonomically correlated. Most rearranged genes were tRNAs, except those of Cotesia vestalis, in which 13 protein-coding genes and 14 tRNA genes changed positions or/and directions through three kinds of gene rearrangement events. Remote inversion is posited to be the result of two independent recombination events. Evolutionary rates were lower in species of the cyclostome group than those of noncyclostomes. Phylogenetic analyses based on complete mitochondrial genomes and secondary structure of rrnS supported a sister-group relationship between Aphidiinae and cyclostomes. Many well accepted relationships within Hymenoptera, such as paraphyly of Symphyta and Evaniomorpha, a sister-group relationship between Orussoidea and Apocrita, and monophyly of Proctotrupomorpha, Ichneumonoidea and Aculeata were robustly confirmed. New hypotheses, such as a sister-group relationship between Evanioidea and Aculeata, were generated. Among holometabolous insects, Hymenoptera was shown to be the sister to all other orders. Mecoptera was recovered as the sister-group of Diptera. Neuropterida (Neuroptera + Megaloptera), and a sister-group relationship with (Diptera + Mecoptera) were supported across all analyses.

CONCLUSIONS

Our comparative studies indicate that mitochondrial genomes are a useful phylogenetic tool at the ordinal level within Holometabola, at the superfamily within Hymenoptera and at the subfamily level within Braconidae. Variation at all of these hierarchical levels suggests that the utility of mitochondrial genomes is likely to be a valuable tool for systematics in other groups of arthropods.

A phylogenomic approach to resolve the arthropod tree of life

Arthropods were the first animals to conquer land and air. They encompass more than three quarters of all described living species. This extraordinary evolutionary success is based on an astoundingly wide array of highly adaptive body organizations. A lack of robustly resolved phylogenetic relationships, however, currently impedes the reliable reconstruction of the underlying evolutionary processes. Here, we show that phylogenomic data can substantially advance our understanding of arthropod evolution and resolve several conflicts among existing hypotheses. We assembled a data set of 233 taxa and 775 genes from which an optimally informative data set of 117 taxa and 129 genes was finally selected using new heuristics and compared with the unreduced data set. We included novel expressed sequence tag (EST) data for 11 species and all published phylogenomic data augmented by recently published EST data on taxonomically important arthropod taxa. This thorough sampling reduces the chance of obtaining spurious results due to stochastic effects of undersampling taxa and genes. Orthology prediction of genes, alignment masking tools, and selection of most informative genes due to a balanced taxa-gene ratio using new heuristics were established. Our optimized data set robustly resolves major arthropod relationships. We received strong support for a sister group relationship of onychophorans and euarthropods and strong support for a close association of tardigrades and cycloneuralia. Within pancrustaceans, our analyses yielded paraphyletic crustaceans and monophyletic hexapods and robustly resolved monophyletic endopterygote insects. However, our analyses also showed for few deep splits that were recently thought to be resolved, for example, the position of myriapods, a remarkable sensitivity to methods of analyses.

[When Tribolium complements the genetics of Drosophila]

With its recently sequenced genome, the red flour beetle Tribolium castaneum became one of the few model organisms with all the main genetic tools. As a coleoptera, it belongs to the most species-rich order of animals. Tribolium is also a worldwide pest for stored dried foods. Regarding developmental biology, Tribolium offers a complementary model to the highly derived Drosophila. For example, the function of many gap and pair-rule segmentation genes is different in both species. These differences reveal the evolutionary plasticity between two modes of development, with a long germ band in fly and a short one in Tribolium. This beetle allowed the identification of a new type of ecdysone receptor for holometabolous insects. Finally, in the search for the juvenile hormone receptor, a crucial result was obtained with experiments that could be performed only with Tribolium, and not with Drosophila. Tribolium, in association with Drosophila, should help to understand the general rules of development in insects.

Loss of lipid synthesis as an evolutionary consequence of a parasitic lifestyle

Evolutionary loss of traits can result from negative selection on a specific phenotype, or if the trait is selectively neutral, because the phenotype associated with the trait has become redundant. Even essential traits may be lost, however, if the resulting phenotypic deficiencies can be compensated for by the environment or a symbiotic partner. Here we demonstrate that loss of an essential me-tabolic trait in parasitic wasps has evolved through environmental compensation. We tested 24 species for the ability to synthesize lipids de novo and collected additional data from the literature. We found the majority of adult parasitoid species to be incapable of synthesizing lipids, and phylogenetic analyses showed that the evolution of lack of lipogenesis is concurrent with that of parasitism in insects. Exploitive host manipulation, in which the host is forced to synthesize lipids to the benefit of the parasitoid, presumably facilitates loss of lipogenesis through environmental compensation. Lipogenesis re-evolved in a small number of parasitoid species, particularly host generalists. The wide range of host species in which generalists are able to develop may impede effective host manipulation and could have resulted in regaining of lipogenic ability in generalist parasitoids. As trait loss through environmental compensation is unnoticed at the phenotypic level, it may be more common than currently anticipated, especially in species involved in intricate symbiotic relationships with other species.

Evolution of a genomic regulatory domain: the role of gene co-option and gene duplication in the Enhancer of split complex

The Drosophila Enhancer of split complex [E(spl)-C] is a remarkable complex of genes many of which are effectors or modulators of Notch signaling. The complex contains different classes of genes including four bearded genes and seven basic helix-loop-helix (bHLH) genes. We examined the evolution of this unusual complex by identifying bearded and bHLH genes in the genome sequences of Arthropods. We find that a four-gene E(spl)-C, containing three bHLH genes and one bearded gene, is an ancient component of the genomes of Crustacea and Insects. The complex is well conserved in insects but is highly modified in Drosophila, where two of the ancestral genes of the complex are missing, and the remaining two have been duplicated multiple times. Through examining the expression of E(spl)-C genes in honeybees, aphids, and Drosophila, we determined that the complex ancestrally had a role in Notch signaling. The expression patterns of genes found inserted into the complex in some insects, or that of ancestral E(spl)-C genes that have moved out of the complex, imply that the E(spl)-C is a genomic domain regulated as a whole by Notch signaling. We hypothesize that the E(spl)-C is a Notch-regulated genomic domain conserved in Arthropod genomes for around 420 million years. We discuss the consequence of this conserved domain for the recruitment of novel genes into the Notch signaling cascade.

Gene duplication and the origins of morphological complexity in pancrustacean eyes, a genomic approach

BACKGROUND

Duplication and divergence of genes and genetic networks is hypothesized to be a major driver of the evolution of complexity and novel features. Here, we examine the history of genes and genetic networks in the context of eye evolution by using new approaches to understand patterns of gene duplication during the evolution of metazoan genomes. We hypothesize that 1) genes involved in eye development and phototransduction have duplicated and are retained at higher rates in animal clades that possess more distinct types of optical design; and 2) genes with functional relationships were duplicated and lost together, thereby preserving genetic networks. To test these hypotheses, we examine the rates and patterns of gene duplication and loss evident in 19 metazoan genomes, including that of Daphnia pulex - the first completely sequenced crustacean genome. This is of particular interest because the pancrustaceans (hexapods+crustaceans) have more optical designs than any other major clade of animals, allowing us to test specifically whether the high amount of disparity in pancrustacean eyes is correlated with a higher rate of duplication and retention of vision genes.

RESULTS

Using protein predictions from 19 metazoan whole-genome projects, we found all members of 23 gene families known to be involved in eye development or phototransduction and deduced their phylogenetic relationships. This allowed us to estimate the number and timing of gene duplication and loss events in these gene families during animal evolution. When comparing duplication/retention rates of these genes, we found that the rate was significantly higher in pancrustaceans than in either vertebrates or non-pancrustacean protostomes. Comparing patterns of co-duplication across Metazoa showed that while these eye-genes co-duplicate at a significantly higher rate than those within a randomly shuffled matrix, many genes with known functional relationships in model organisms did not co-duplicate more often than expected by chance.

CONCLUSIONS

Overall, and when accounting for factors such as differential rates of whole-genome duplication in different groups, our results are broadly consistent with the hypothesis that genes involved in eye development and phototransduction duplicate at a higher rate in Pancrustacea, the group with the greatest variety of optical designs. The result that these genes have a significantly high number of co-duplications and co-losses could be influenced by shared functions or other unstudied factors such as synteny. Since we did not observe co-duplication/co-loss of genes for all known functional modules (e.g. specific regulatory networks), the interactions among suites of known co-functioning genes (modules) may be plastic at the temporal scale of analysis performed here. Other factors in addition to gene duplication - such as cis-regulation, heterotopy, and co-option - are also likely to be strong factors in the diversification of eye types.

The transcription factor Krüppel homolog 1 is linked to hormone mediated social organization in bees

BACKGROUND

Regulation of worker behavior by dominant queens or workers is a hallmark of insect societies, but the underlying molecular mechanisms and their evolutionary conservation are not well understood. Honey bee and bumble bee colonies consist of a single reproductive queen and facultatively sterile workers. The queens' influences on the workers are mediated largely via inhibition of juvenile hormone titers, which affect division of labor in honey bees and worker reproduction in bumble bees. Studies in honey bees identified a transcription factor, Krüppel-homolog 1 (Kr-h1), whose expression in worker brains is significantly downregulated in the presence of a queen or queen pheromone and higher in forager bees, making this gene an ideal candidate for examining the evolutionary conservation of socially regulated pathways in Hymenoptera.

RESULTS

In contrast to honey bees, bumble bees foragers do not have higher Kr-h1 levels relative to nurses: in one of three colonies levels were similar in nurses and foragers, and in two colonies levels were higher in nurses. Similarly to honey bees, brain Kr-h1 levels were significantly downregulated in the presence versus absence of a queen. Furthermore, in small queenless groups, Kr-h1 levels were downregulated in subordinate workers with undeveloped ovaries relative to dominant individuals with active ovaries. Brain Kr-h1 levels were upregulated by juvenile hormone treatment relative to a vehicle control. Finally, phylogenetic analysis indicates that KR-H1 orthologs are presence across insect orders. Though this protein is highly conserved between honey bees and bumble bees, there are significant differences between orthologs of insects from different orders.

CONCLUSIONS

Our results suggest that Kr-h1 is associated with juvenile hormone mediated regulation of reproduction in bumble bees. The expression of this transcription factor is inhibited by the queen and associated with endocrine mediated regulation of social organization in two species of bees. Thus, KR-H1 may transcriptionally regulate a conserved genetic module that is part of a pathway that has been co-opted to function in social behavior, and adjusts the behavior of workers to their social environmental context.

The origins of species richness in the Hymenoptera: insights from a family-level supertree

BACKGROUND

The order Hymenoptera (bees, ants, wasps, sawflies) contains about eight percent of all described species, but no analytical studies have addressed the origins of this richness at family-level or above. To investigate which major subtaxa experienced significant shifts in diversification, we assembled a family-level phylogeny of the Hymenoptera using supertree methods. We used sister-group species-richness comparisons to infer the phylogenetic position of shifts in diversification.

RESULTS

The supertrees most supported by the underlying input trees are produced using matrix representation with compatibility (MRC) (from an all-in and a compartmentalised analysis). Whilst relationships at the tips of the tree tend to be well supported, those along the backbone of the tree (e.g. between Parasitica superfamilies) are generally not. Ten significant shifts in diversification (six positive and four negative) are found common to both MRC supertrees. The Apocrita (wasps, ants, bees) experienced a positive shift at their origin accounting for approximately 4,000 species. Within Apocrita other positive shifts include the Vespoidea (vespoid wasps/ants containing 24,000 spp.), Anthophila + Sphecidae (bees/thread-waisted wasps; 22,000 spp.), Bethylidae + Chrysididae (bethylid/cuckoo wasps; 5,200 spp.), Dryinidae (dryinid wasps; 1,100 spp.), and Proctotrupidae (proctotrupid wasps; 310 spp.). Four relatively species-poor families (Stenotritidae, Anaxyelidae, Blasticotomidae, Xyelidae) have undergone negative shifts. There are some two-way shifts in diversification where sister taxa have undergone shifts in opposite directions.

CONCLUSIONS

Our results suggest that numerous phylogenetically distinctive radiations contribute to the richness of large clades. They also suggest that evolutionary events restricting the subsequent richness of large clades are common. Problematic phylogenetic issues in the Hymenoptera are identified, relating especially to superfamily validity (e.g. "Proctotrupoidea", "Mymarommatoidea"), and deeper apocritan relationships. Our results should stimulate new functional studies on the causes of the diversification shifts we have identified. Possible drivers highlighted for specific adaptive radiations include key anatomical innovations, the exploitation of rich host groups, and associations with angiosperms. Low richness may have evolved as a result of geographical isolation, specialised ecological niches, and habitat loss or competition.

Tubulin evolution in insects: gene duplication and subfunctionalization provide specialized isoforms in a functionally constrained gene family

Background

The completion of 19 insect genome sequencing projects spanning six insect orders provides the opportunity to investigate the evolution of important gene families, here tubulins. Tubulins are a family of eukaryotic structural genes that form microtubules, fundamental components of the cytoskeleton that mediate cell division, shape, motility, and intracellular trafficking. Previous in vivo studies in Drosophila find a stringent relationship between tubulin structure and function; small, biochemically similar changes in the major alpha 1 or testis-specific beta 2 tubulin protein render each unable to generate a motile spermtail axoneme. This has evolutionary implications, not a single non-synonymous substitution is found in beta 2 among 17 species of Drosophila and Hirtodrosophila flies spanning 60 Myr of evolution. This raises an important question, How do tubulins evolve while maintaining their function? To answer, we use molecular evolutionary analyses to characterize the evolution of insect tubulins.

Results

Sixty-six alpha tubulins and eighty-six beta tubulin gene copies were retrieved and subjected to molecular evolutionary analyses. Four ancient clades of alpha and beta tubulins are found in insects, a major isoform clade (alpha 1, beta 1) and three minor, tissue-specific clades (alpha 2-4, beta 2-4). Based on a Homarus americanus (lobster) outgroup, these were generated through gene duplication events on major beta and alpha tubulin ancestors, followed by subfunctionalization in expression domain. Strong purifying selection acts on all tubulins, yet maximum pairwise amino acid distances between tubulin paralogs are large (0.464 substitutions/site beta tubulins, 0.707 alpha tubulins). Conversely orthologs, with the exception of reproductive tissue isoforms, show little sequence variation except in the last 15 carboxy terminus tail (CTT) residues, which serve as sites for post-translational modifications (PTMs) and interactions with microtubule-associated proteins. CTT residues overwhelming comprise the co-evolving residues between Drosophila alpha 2 and beta 3 tubulin proteins, indicating CTT specializations can be mediated at the level of the tubulin dimer. Gene duplications post-dating separation of the insect orders are unevenly distributed, most often appearing in major alpha 1 and minor beta 2 clades. More than 40 introns are found in tubulins. Their distribution among tubulins reveals that insertion and deletion events are common, surprising given their potential for disrupting tubulin coding sequence. Compensatory evolution is found in Drosophila beta 2 tubulin cis-regulation, and reveals selective pressures acting to maintain testis expression without the use of previously identified testis cis-regulatory elements.

Conclusion

Tubulins have stringent structure/function relationships, indicated by strong purifying selection, the loss of many gene duplication products, alpha-beta co-evolution in the tubulin dimer, and compensatory evolution in beta 2 tubulin cis-regulation. They evolve through gene duplication, subfunctionalization in expression domain and divergence of duplication products, largely in CTT residues that mediate interactions with other proteins. This has resulted in the tissue-specific minor insect isoforms, and in particular the highly diverse α3, α4, and β2 reproductive tissue-specific tubulin isoforms, illustrating that even a highly conserved protein family can participate in the adaptive process and respond to sexual selection.

Ribosomal protein genes of holometabolan insects reject the Halteria, instead revealing a close affinity of Strepsiptera with Coleoptera

The phylogenetic relationships among holometabolan insect orders remain poorly known, despite a wealth of previous studies. In particular, past attempts to clarify the sister-group of the enigmatic order Strepsiptera with rRNA genes have led to intense debate about long-branch attraction (the 'Strepsiptera problem'), without resolving the taxonomic question at hand. Here, we appealed to alternative nuclear sequences of 27 ribosomal proteins (RPs) to generate a data matrix of 10,731 nucleotides for 22 holometabolan taxa, including two strepsipteran species. Phylogenetic relationships among holometabolan insects were analyzed under several nucleotide-coding schemes to explore differences in signal and systematic biases. Saturation and compositional bias particularly affected third positions, which greatly differed in AT content (18-72%). Such confounding factors were best reduced by R-Y coding and removal of third codon positions, resulting in more strongly supported topologies, whereas amino acid coding gave poor resolution. The placement of Strepsiptera with Coleoptera (the Coleopterida) was recovered under most coding schemes and analytical methods, if often with modest support and ambiguity. In contrast, an alternative sister-group with Diptera (the Halteria) was only found in one analysis using parsimony, and weakly supported. The topologies here generally support a Coleoptera+Strepsiptera as sister-group to Mecopterida (Siphonaptera+Mecoptera+Diptera+Lepidoptera+Trichoptera), while Hymenoptera were always recovered as sister-group to the remaining Holometabola.

Antimicrobial peptide-like genes in Nasonia vitripennis: a genomic perspective

BACKGROUND

Antimicrobial peptides (AMPs) are an essential component of innate immunity which can rapidly respond to diverse microbial pathogens. Insects, as a rich source of AMPs, attract great attention of scientists in both understanding of the basic biology of the immune system and searching molecular templates for anti-infective drug design. Despite a large number of AMPs have been identified from different insect species, little information in terms of these peptides is available from parasitic insects.

RESULTS

By using integrated computational approaches to systemically mining the Hymenopteran parasitic wasp Nasonia vitripennis genome, we establish the first AMP repertoire whose members exhibit extensive sequence and structural diversity and can be distinguished into multiple molecular types, including insect and fungal defensin-like peptides (DLPs) with the cysteine-stabilized alpha-helical and beta-sheet (CSalphabeta) fold; Pro- or Gly-rich abaecins and hymenoptaecins; horseshoe crab tachystatin-type AMPs with the inhibitor cystine knot (ICK) fold; and a linear alpha-helical peptide. Inducible expression pattern of seven N. vitripennis AMP genes were verified, and two representative peptides were synthesized and functionally identified to be antibacterial. In comparison with Apis mellifera (Hymenoptera) and several non-Hymenopteran model insects, N. vitripennis has evolved a complex antimicrobial immune system with more genes and larger protein precursors. Three classical strategies that are likely responsible for the complexity increase have been recognized: 1) Gene duplication; 2) Exon duplication; and 3) Exon-shuffling.

CONCLUSION

The present study established the N. vitripennis peptidome associated with antimicrobial immunity by using a combined computational and experimental strategy. As the first AMP repertoire of a parasitic wasp, our results offer a basic platform for further studying the immunological and evolutionary significances of these newly discovered AMP-like genes in this class of insects.

giant is a bona fide gap gene in the intermediate germband insect, Oncopeltus fasciatus

Drosophila undergoes a form of development termed long germ segmentation, where all segments are specified nearly simultaneously so that by the blastoderm stage, the entire body plan has been determined. This mode of segmentation is evolutionarily derived. Most insects undergo short or intermediate germ segmentation, where only anterior segments are specified early, and posterior segments are sequentially specified during germband elongation. These embryological differences imply that anterior and posterior segments might rely upon different molecular mechanisms. In Drosophila, embryos mutant for giant show a gap in the anterior as well fusions of several abdominal segments. In Tribolium, a short germ beetle, giant is required for segmental identity, but not formation, in gnathal segments and also for segmentation of the entire abdomen. This raises the possibility that giant might not act as a gap gene in short and intermediate germ insects. Oncopeltus fasciatus is an intermediate germ insect that is an outgroup to the clade containing Drosophila and Tribolium. We cloned the Oncopeltus homolog of giant and determined its expression and function during segmentation. We find that Oncopeltus giant is a canonical gap gene in the maxillary and labial segments and also plays a gap-like role in the first four abdominal segments. Our results suggest that giant was a bona fide gap gene in the ancestor of these insects with this role being lost in the lineage leading towards Tribolium. This highlights the conservation of anterior patterning and evolutionary plasticity of the genetic regulation controlling posterior segmentation, even in short and intermediate germ insects.

A functional DNA methylation system in the pea aphid, Acyrthosiphon pisum

Methylation of cytosine is one of the main epigenetic mechanisms involved in controlling gene expression. Here we show that the pea aphid (Acyrthosiphon pisum) genome possesses homologues to all the DNA methyltransferases found in vertebrates, and that 0.69% (+/-0.25%) of all cytosines are methylated. Identified methylation sites are predominantly restricted to the coding sequence of genes at CpG sites. We identify twelve methylated genes, including genes that interact with juvenile hormone, a key endocrine signal in insects. Bioinformatic prediction using CpG ratios for all predicted genes suggest that a large proportion of genes are methylated within the pea aphid.

Analysis of a normalised expressed sequence tag (EST) library from a key pollinator, the bumblebee Bombus terrestris

BACKGROUND

The bumblebee, Bombus terrestris (Order Hymenoptera), is of widespread importance. This species is extensively used for commercial pollination in Europe, and along with other Bombus spp. is a key member of natural pollinator assemblages. Furthermore, the species is studied in a wide variety of biological fields. The objective of this project was to create a B. terrestris EST resource that will prove to be valuable in obtaining a deeper understanding of this significant social insect.

RESULTS

A normalised cDNA library was constructed from the thorax and abdomen of B. terrestris workers in order to enhance the discovery of rare genes. A total of 29'428 ESTs were sequenced. Subsequent clustering resulted in 13'333 unique sequences. Of these, 58.8 percent had significant similarities to known proteins, with 54.5 percent having a "best-hit" to existing Hymenoptera sequences. Comparisons with the honeybee and other insects allowed the identification of potential candidates for gene loss, pseudogene evolution, and possible incomplete annotation in the honeybee genome. Further, given the focus of much basic research and the perceived threat of disease to natural and commercial populations, the immune system of bumblebees is a particularly relevant component. Although the library is derived from unchallenged bees, we still uncover transcription of a number of immune genes spanning the principally described insect immune pathways. Additionally, the EST library provides a resource for the discovery of genetic markers that can be used in population level studies. Indeed, initial screens identified 589 simple sequence repeats and 854 potential single nucleotide polymorphisms.

CONCLUSION

The resource that these B. terrestris ESTs represent is valuable for ongoing work. The ESTs provide direct evidence of transcriptionally active regions, but they will also facilitate further functional genomics, gene discovery and future genome annotation. These are important aspects in obtaining a greater understanding of this key pollinator species.

Identification and characterization of Nasonia Pax genes

Pax genes are a group of critical developmental transcriptional regulators in both invertebrates and vertebrates, characterized by the presence of a paired DNA-binding domain. Pax proteins also often contain an octapeptide motif and a C-terminal homeodomain. The genome of Nasonia vitripennis (Hymenoptera) has recently become available, and analysis of this genome alongside Apis mellifera allowed us to contribute to the phylogeny of this gene family in insects. Nasonia, a parasitic wasp, has independently evolved a similar mode of development to that of the well-studied Drosophila, making it an excellent model system for comparative studies of developmental gene networks. We report the characterization of the seven Nasonia Pax genes. We describe their genomic organization, and the embryonic expression of three of them, and uncover wider conservation of the octapeptide motif than previously described.

Inventory and phylogenomic distribution of meiotic genes in Nasonia vitripennis and among diverse arthropods

The parasitoid jewel wasp Nasonia vitripennis reproduces by haplodiploidy (arrhenotokous parthenogenesis). In diploid females, meiosis occurs during oogenesis, but in haploid males spermatogenesis is ameiotic and involves a single equational division. Here we describe the phylogenomic distribution of meiotic genes in N. vitripennis and in 10 additional arthropods. Homologues for 39 meiosis-related genes (including seven meiosis-specific genes) were identified in N. vitripennis. The meiotic genes missing from N. vitripennis are also sporadically absent in other arthropods, suggesting that certain meiotic genes are dispensable for meiosis. Among an additional set of 15 genes thought to be specific for male meiosis in Drosophila, two genes (bol and crl) were identified in N. vitripennis and Apis mellifera (both for which canonical meiosis is absent in males) and in other arthropods. The distribution of meiotic genes across arthropods and the impact of gene duplications and reproductive modes on meiotic gene evolution are discussed.

Contrasting patterns of selective constraints in nuclear-encoded genes of the oxidative phosphorylation pathway in holometabolous insects and their possible role in hybrid breakdown in Nasonia

The principal energy generating system in animals is the oxidative phosphorylation (OXPHOS) pathway, which depends on the tight interaction of nuclear and mitochondrial encoded genes to function properly. Mitochondrial genes accumulate substitutions more quickly than nuclear genes, yet the impact of selection on mitochondrial genes is significantly reduced relative to nuclear genes due to the non-recombining nature of the mitochondrial genome and its predicted smaller effective population size. It has therefore been hypothesized that the nuclear encoded genes of the OXPHOS pathway are under strong selective pressure to compensate for the accumulation of deleterious nucleotide substitutions in mitochondrial encoded OXPHOS genes; a process known as compensatory co-adaptation. We evaluated this hypothesis by analyzing nuclear encoded OXPHOS genes for signatures of positive selection as well as evolutionary constraints at amino acid sites. We considered OXPHOS genes of six holometabolous insects and their orthologs from three Nasonia parasitoid wasps; hybrids of which suffer from an increased mortality rate caused by cytonuclear genic incompatibilities. Although nuclear OXPHOS genes are typically highly conserved, we found significant evidence for elevated amino acid divergence in four of the 59 studied nuclear encoded OXPHOS genes. We also found that three of these four genes, as well as six other OXPHOS genes, contain amino acid substitutions between Nasonia species at evolutionarily constrained sites. It is possible that these genes account for the reported incompatibility in Nasonia hybrids and their characterization may lead to a better understanding of the role of positive selection in the genetics of speciation.

454 pyrosequencing based transcriptome analysis of Zygaena filipendulae with focus on genes involved in biosynthesis of cyanogenic glucosides

Background

An essential driving component in the co-evolution of plants and insects is the ability to produce and handle bioactive compounds. Plants produce bioactive natural products for defense, but some insects detoxify and/or sequester the compounds, opening up for new niches with fewer competitors. To study the molecular mechanism behind the co-adaption in plant-insect interactions, we have investigated the interactions between Lotus corniculatus and Zygaena filipendulae. They both contain cyanogenic glucosides which liberate toxic hydrogen cyanide upon breakdown. Moths belonging to the Zygaena family are the only insects known, able to carry out both de novo biosynthesis and sequestration of the same cyanogenic glucosides as those from their feed plants. The biosynthetic pathway for cyanogenic glucoside biosynthesis in Z. filipendulae proceeds using the same intermediates as in the well known pathway from plants, but none of the enzymes responsible have been identified. A genomics strategy founded on 454 pyrosequencing of the Z. filipendulae transcriptome was undertaken to identify some of these enzymes in Z. filipendulae.

Results

Comparisons of the Z. filipendulae transcriptome with the sequenced genomes of Bombyx mori, Drosophila melanogaster, Tribolium castaneum, Apis mellifera and Anopheles gambiae indicate a high coverage of the Z. filipendulae transcriptome. 11% of the Z. filipendulae transcriptome sequences were assigned to Gene Ontology categories. Candidate genes for enzymes functioning in the biosynthesis of cyanogenic glucosides (cytochrome P450 and family 1 glycosyltransferases) were identified based on sequence length, number of copies and presence/absence of close homologs in D. melanogaster, B. mori and the cyanogenic butterfly Heliconius. Examination of biased codon usage, GC content and selection on gene candidates support the notion of cyanogenesis as an "old" trait within Ditrysia, as well as its origins being convergent between plants and insects.

Conclusion

Pyrosequencing is an attractive approach to gain access to genes in the biosynthesis of bio-active natural products from insects and other organisms, for which the genome sequence is not known. Based on analysis of the Z. filipendulae transcriptome, promising gene candidates for biosynthesis of cyanogenic glucosides was identified, and the suitability of Z. filipendulae as a model system for cyanogenesis in insects is evident.

Evolving enhancer-promoter interactions within the tinman complex of the flour beetle, Tribolium castaneum

Modifications of cis-regulatory DNAs, particularly enhancers, underlie changes in gene expression during animal evolution. Here, we present evidence for a distinct mechanism of regulatory evolution, whereby a novel pattern of gene expression arises from altered gene targeting of a conserved enhancer. The tinman gene complex (Tin-C) controls the patterning of dorsal mesodermal tissues, including the dorsal vessel or heart in Drosophila. Despite broad conservation of Tin-C gene expression patterns in the flour beetle (Tribolium castaneum), the honeybee (Apis mellifera) and the fruit fly (Drosophila melanogaster), the expression of a key pericardial determinant, ladybird, is absent from the dorsal mesoderm of Tribolium embryos. Evidence is presented that this loss in expression is replaced by expression of C15, the neighboring gene in the complex. This switch in expression from ladybird to C15 appears to arise from an inversion within the tinman complex, which redirects a conserved ladybird 3' enhancer to regulate C15. In Drosophila, this enhancer fails to activate C15 expression owing to the activity of an insulator at the intervening ladybird early promoter. By contrast, a chromosomal inversion allows the cardiac enhancer to bypass the ladybird insulator in Tribolium. Given the high frequency of genome rearrangements in insects, it is possible that such enhancer switching might be widely used in the diversification of the arthropods.

Sex determination in honeybees: two separate mechanisms induce and maintain the female pathway

Sex determination in honeybees is realized by the csd and the fem gene that establish and maintain, throughout development, sexual fates via the control of alternative splicing.

Organisms have evolved a bewildering diversity of mechanisms to generate the two sexes. The honeybee (Apis mellifera) employs an interesting system in which sex is determined by heterozygosity at a single locus (the Sex Determination Locus) harbouring the complementary sex determiner (csd) gene. Bees heterozygous at Sex Determination Locus are females, whereas bees homozygous or hemizygous are males. Little is known, however, about the regulation that links sex determination to sexual differentiation. To investigate the control of sexual development in honeybees, we analyzed the functions and the regulatory interactions of genes involved in the sex determination pathway. We show that heterozygous csd is only required to induce the female pathway, while the feminizer (fem) gene maintains this decision throughout development. By RNAi induced knockdown we show that the fem gene is essential for entire female development and that the csd gene exclusively processes the heterozygous state. Fem activity is also required to maintain the female determined pathway throughout development, which we show by mosaic structures in fem-repressed intersexuals. We use expression of Fem protein in males to demonstrate that the female maintenance mechanism is controlled by a positive feedback splicing loop in which Fem proteins mediate their own synthesis by directing female fem mRNA splicing. The csd gene is only necessary to induce this positive feedback loop in early embryogenesis by directing splicing of fem mRNAs. Finally, fem also controls the splicing of Am-doublesex transcripts encoding conserved male- and female-specific transcription factors involved in sexual differentiation. Our findings reveal how the sex determination process is realized in honeybees differing from Drosophila melanogaster.

Sexual differentiation is a fundamental process in the animal kingdom, and different species have evolved a bewildering diversity of mechanisms to generate the two sexes in the proper proportions. Sex determination in honeybees (Apis mellifera) provides an interesting and unusual system to study, as it is governed by heterozygosity of a single locus harbouring the complementary sex determiner gene (csd), in contrast to the well-studied sex chromosome system of Drosophila melanogaster. We show that the female sex determination pathway is exclusively induced by the csd gene in early embryogenesis. Later on and throughout development this inductive signal is maintained via a positive feedback loop of the feminizer (fem) gene, in which the Fem protein mediates its own synthesis. The findings reveal how the sex determination process in honeybees is realized by the regulation and function of two genes differing from Drosophila.

An EST screen from the annelid Pomatoceros lamarckii reveals patterns of gene loss and gain in animals

Background

Since the drastic reorganisation of the phylogeny of the animal kingdom into three major clades of bilaterians; Ecdysozoa, Lophotrochozoa and Deuterostomia, it became glaringly obvious that the selection of model systems with extensive molecular resources was heavily biased towards only two of these three clades, namely the Ecdysozoa and Deuterostomia. Increasing efforts have been put towards redressing this imbalance in recent years, and one of the principal phyla in the vanguard of this endeavour is the Annelida.

Results

In the context of this effort we here report our characterisation of an Expressed Sequence Tag (EST) screen in the serpulid annelid, Pomatoceros lamarckii. We have sequenced over 5,000 ESTs which consolidate into over 2,000 sequences (clusters and singletons). These sequences are used to build phylogenetic trees to estimate relative branch lengths amongst different taxa and, by comparison to genomic data from other animals, patterns of gene retention and loss are deduced.

Conclusion

The molecular phylogenetic trees including the P. lamarckii sequences extend early observations that polychaetes tend to have relatively short branches in such trees, and hence are useful taxa with which to reconstruct gene family evolution. Also, with the availability of lophotrochozoan data such as that of P. lamarckii, it is now possible to make much more accurate reconstructions of the gene complement of the ancestor of the bilaterians than was previously possible from comparisons of ecdysozoan and deuterostome genomes to non-bilaterian outgroups. It is clear that the traditional molecular model systems for protostomes (e.g. Drosophila melanogaster and Caenorhabditis elegans), which are restricted to the Ecdysozoa, have undergone extensive gene loss during evolution. These ecdysozoan systems, in terms of gene content, are thus more derived from the bilaterian ancestral condition than lophotrochozoan systems like the polychaetes, and thus cannot be used as good, general representatives of protostome genomes. Currently sequenced insect and nematode genomes are less suitable models for deducing bilaterian ancestral states than lophotrochozoan genomes, despite the array of powerful genetic and mechanistic manipulation techniques in these ecdysozoans. A distinct category of genes that includes those present in non-bilaterians and lophotrochozoans, but which are absent from ecdysozoans and deuterostomes, highlights the need for further lophotrochozoan data to gain a more complete understanding of the gene complement of the bilaterian ancestor.

Evolutionary conservation and changes in insect TRP channels

Background

TRP (Transient Receptor Potential) channels respond to diverse stimuli and thus function as the primary integrators of varied sensory information. They are also activated by various compounds and secondary messengers to mediate cell-cell interactions as well as to detect changes in the local environment. Their physiological roles have been primarily characterized only in mice and fruit flies, and evolutionary studies are limited. To understand the evolution of insect TRP channels and the mechanisms of integrating sensory inputs in insects, we have identified and compared TRP channel genes in Drosophila melanogaster, Bombyx mori, Tribolium castaneum, Apis mellifera, Nasonia vitripennis, and Pediculus humanus genomes as part of genome sequencing efforts.

Results

All the insects examined have 2 TRPV, 1 TRPN, 1 TRPM, 3 TRPC, and 1 TRPML subfamily members, demonstrating that these channels have the ancient origins in insects. The common pattern also suggests that the mechanisms for detecting mechanical and visual stimuli and maintaining lysosomal functions may be evolutionarily well conserved in insects. However, a TRPP channel, the most ancient TRP channel, is missing in B. mori, A. mellifera, and N. vitripennis. Although P. humanus and D. melanogaster contain 4 TRPA subfamily members, the other insects have 5 TRPA subfamily members. T. castaneum, A. mellifera, and N. vitripennis contain TRPA5 channels, which have been specifically retained or gained in Coleoptera and Hymenoptera. Furthermore, TRPA1, which functions for thermotaxis in Drosophila, is missing in A. mellifera and N. vitripennis; however, they have other Hymenoptera-specific TRPA channels (AmHsTRPA and NvHsTRPA). NvHsTRPA expressed in HEK293 cells is activated by temperature increase, demonstrating that HsTRPAs function as novel thermal sensors in Hymenoptera.

Conclusion

The total number of insect TRP family members is 13-14, approximately half that of mammalian TRP family members. As shown for mammalian TRP channels, this may suggest that single TRP channels are responsible for integrating diverse sensory inputs to maintain the insect sensory systems. The above results demonstrate that there are both evolutionary conservation and changes in insect TRP channels. In particular, the evolutionary processes have been accelerated in the TRPA subfamily, indicating divergence in the mechanisms that insects use to detect environmental temperatures.

Resolving insect phylogeny: The significance of cephalic structures of the Nannomecoptera in understanding endopterygote relationships

The Nannomecoptera are among the most enigmatic and controversial taxa of endopterygote insects, the phylogenetic resolution of which is crucial to understanding the evolution of neopteran insects. Once considered a subordinate lineage within the Mecoptera, renewed interest in nannochoristids has been prompted by evidence that the Nannomecoptera are not admissible to the clade of extant scorpion flies but are more likely to belong to the clade Siphonaptera + Nannomecoptera + Diptera. The overarching purpose of the present account is to provide novel and extensive morphological character traits in addition to those already existing for adult structures. The aim is to determine if these traits support molecular data sets that have been suggested elsewhere to clarify the phyletic position of Nannochoristidae. This account focuses on nannomecopteran larvae, which unlike those of other mecopterans have received little attention. Thus, the thrust of this investigation is to provide detailed anatomical data on nannochoristid larvae for a targeted inquiry into their phylogenetic affinities. The described characters are discussed and presented in a data matrix comprising representatives of all endopterygote orders. While the study is based primarily on the New Zealand species Nannochorista philpotti, it is proposed that all nannomecopteran larvae will prove to be similar to this taxon in most if not all structural features of significance to a higher-level phylogenetic context.

The small heat shock protein (sHSP) genes in the silkworm, Bombyx mori, and comparative analysis with other insect sHSP genes

Background

Small heat shock proteins (sHSPs) are products of heat shock response and of other stress responses, and ubiquitous in all three domains of life, archaea, bacteria, and eukarya. They mainly function as molecular chaperones to protect proteins from being denatured in extreme conditions. Study on insect sHSPs could provide some insights into evolution of insects that have adapted to diverse niches in the world.

Results

Taking advantage of the newly assembled genome sequence, we performed a genome-wide analysis of the candidate sHSP genes in the silkworm, Bombyx mori. Based on known silkworm sHSP sequences, we identified 16 silkworm sHSP genes. Most of them are distributed on two silkworm chromosomes 5 and 27, respectively. 15 of 16 silkworm sHSPs have expression evidence. The comparative analysis of insect sHSPs from B. mori, Drosophila melanogaster, Apis mellifera, Tribolium castaneum, and Anopheles gambiae revealed that there is only one orthologous cluster whereas remaining clusters are species-specific on the phylogenetic tree. This suggested that most of sHSPs might have diverged in function across insects investigated. In addition, the data presented in this study also revealed that sHSPs in the insect orthologous cluster are highly conserved in both sequence and expression pattern. In sum, insect sHSPs show a completely different evolutionary pattern from that found in vertebrate sHSPs.

Conclusion

B. mori has the largest number of insect sHSP genes characterized to date, including 16 genes. The inference that most species-specific sHSPs might have diverged in function across insects investigated will help us understand the adaptability of these insects to diverse environments.

A phylogenomic approach to resolve the basal pterygote divergence

One of the most fascinating Bauplan transitions in the animal kingdom was the invention of insect wings, a change that also contributed to the success and enormous diversity of this animal group. However, the origin of insect flight and the relationships of basal winged insect orders are still controversial. Three hypotheses have been proposed to explain the phylogeny of winged insects: 1) the traditional Palaeoptera hypothesis (Ephemeroptera + Odonata, Neoptera), 2) the Metapterygota hypothesis (Ephemeroptera, Odonata + Neoptera), and 3) the Chiastomyaria hypothesis (Odonata, Ephemeroptera + Neoptera). Neither phylogenetic analyses of single genes nor even multiple marker systems (e.g., molecular markers + morphological characters) have yet been able to conclusively resolve basal pterygote divergences. A possible explanation for the lack of resolution is that the divergences took place in the mid-Devonian within a short period of time and attempts to solve this problem have been confounded by the major challenge of finding molecular markers to accurately track these short ancient internodes. Although phylogenomic data are available for Neoptera and some wingless (apterygote) orders, they are lacking for the crucial Odonata and Ephemeroptera orders. We adopt a multigene approach including data from two new expressed sequence tag projects-from the orders Ephemeroptera (Baetis sp.) and Odonata (Ischnura elegans)-to evaluate the potential of phylogenomic analyses in clarifying this unresolved issue. We analyzed two data sets that differed in represented taxa, genes, and overall sequence lengths: maxspe (15 taxa, 125 genes, and 31,643 amino acid positions) and maxgen (8 taxa, 150 genes, and 42,541 amino acid positions). Maximum likelihood and Bayesian inference analyses both place the Odonata at the base of the winged insects. Furthermore, statistical hypotheses testing rejected both the Palaeoptera and the Metapterygota hypotheses. The comprehensive molecular data set developed here provides conclusive support for odonates as the most basal winged insect order (Chiastomyaria hypothesis). Data quality assessment indicates that proteins involved in cellular processes and signaling harbor the most informative phylogenetic signal.

Single-copy nuclear genes resolve the phylogeny of the holometabolous insects

BACKGROUND

Evolutionary relationships among the 11 extant orders of insects that undergo complete metamorphosis, called Holometabola, remain either unresolved or contentious, but are extremely important as a context for accurate comparative biology of insect model organisms. The most phylogenetically enigmatic holometabolan insects are Strepsiptera or twisted wing parasites, whose evolutionary relationship to any other insect order is unconfirmed. They have been controversially proposed as the closest relatives of the flies, based on rDNA, and a possible homeotic transformation in the common ancestor of both groups that would make the reduced forewings of Strepsiptera homologous to the reduced hindwings of Diptera. Here we present evidence from nucleotide sequences of six single-copy nuclear protein coding genes used to reconstruct phylogenetic relationships and estimate evolutionary divergence times for all holometabolan orders.

RESULTS

Our results strongly support Hymenoptera as the earliest branching holometabolan lineage, the monophyly of the extant orders, including the fleas, and traditionally recognized groupings of Neuropteroidea and Mecopterida. Most significantly, we find strong support for a close relationship between Coleoptera (beetles) and Strepsiptera, a previously proposed, but analytically controversial relationship. Exploratory analyses reveal that this relationship cannot be explained by long-branch attraction or other systematic biases. Bayesian divergence times analysis, with reference to specific fossil constraints, places the origin of Holometabola in the Carboniferous (355 Ma), a date significantly older than previous paleontological and morphological phylogenetic reconstructions. The origin and diversification of most extant insect orders began in the Triassic, but flourished in the Jurassic, with multiple adaptive radiations producing the astounding diversity of insect species for which these groups are so well known.

CONCLUSION

These findings provide the most complete evolutionary framework for future comparative studies on holometabolous model organisms and contribute strong evidence for the resolution of the 'Strepsiptera problem', a long-standing and hotly debated issue in insect phylogenetics.

Episodes in insect evolution

This article derives from a society-wide symposium organized by Timothy Bradley and Adriana Briscoe and presented at the 2009 annual meeting of the Society for Integrative and Comparative Biology in Boston, Massachusetts. David Grimaldi provided the opening presentation in which he outlined the major evolutionary events in the formation and subsequent diversification of the insect clade. This presentation was followed by speakers who detailed the evolutionary history of specific physiological and/or behavioral traits that have caused insects to be both ecologically successful and fascinating as subjects for biological study. These include a review of the evolutionary history of the insects, the origins of flight, osmoregulation, the evolution of tracheal systems, the evolution of color vision, circadian clocks, and the evolution of eusociality. These topics, as covered by the speakers, provide an overview of the pattern and timing of evolutionary diversification and specialization in the group of animals we know as insects.

Evolutionary history of the iroquois/Irx genes in metazoans

BACKGROUND

The iroquois (iro/Irx) genes encode transcriptional regulators that belong to the TALE superclass of homeodomain proteins and have key functions during development in both vertebrates and insects. The Irx genes occur in one or two genomic clusters containing three genes each within the Drosophila and several vertebrate genomes, respectively. The similar genomic organization in Drosophila and vertebrates is widely considered as a result of convergent evolution, due to independent tandem gene duplications. In this study, we investigate the evolutionary history of the Irx genes at the scale of the whole metazoan kingdom.

RESULTS

We identified in silico the putative full complement of Irx genes in the sequenced genomes of 36 different species representative of the main metazoan lineages, including non bilaterian species, several arthropods, non vertebrate chordates, and a basal vertebrate, the sea lamprey. We performed extensive phylogenetic analyses of the identified Irx genes and defined their genomic organizations. We found that, in most species, there are several Irx genes, these genes form two to four gene clusters, and the Irx genes are physically linked to a structurally and functionally unrelated gene known as CG10632 in Drosophila.

CONCLUSION

Three main conclusions can be drawn from our study. First, an Irx cluster composed of two genes, araucan/caupolican and mirror, is ancestral to the crustaceans+insects clade and has been strongly conserved in this clade. Second, three Irx genes were probably present in the last common ancestor of vertebrates and the duplication that has given rise to the six genes organized into two clusters found in most vertebrates, likely occurred in the gnathostome lineage after its separation from sea lampreys. Third, the clustered organization of the Irx genes in various evolutionary lineages may represent an exceptional case of convergent evolution or may point to the existence of an Irx gene cluster ancestral to bilaterians.

The Fruitless gene in Nasonia displays complex sex-specific splicing and contains new zinc finger domains

The transcription factor Fruitless exerts a broad range of functions during Drosophila development, the most apparent of which is the determination of sexual behavior in males. Although fruitless sequences are found in other insect orders, little is known about fruitless structure and function outside Diptera. We have performed a thorough analysis of fruitless transcripts in the haplo-diploid wasp Nasonia vitripennis and found both sex-specific and non-sex-specific transcripts similar to those found in Drosophila. In Nasonia, however, a novel, large fruitless transcript is present in females only. Putative binding sites for sex-specific splicing factors found in Nasonia fruitless and doublesex as well as Apis mellifera doublesex transcripts were sufficient to identify a corresponding female-specific fruitless exon in A. mellifera, suggesting that similar factors in both hymenopteran species could be responsible for sex-specific splicing of both genes. Furthermore, new C(2)H(2) zinc finger domains found in Nasonia fruitless transcripts were also identified in the fruitless locus of major holometabolous insect species but not in drosophilids. Conservation of important domains and sex-specific splicing in Diptera and Hymenoptera support the hypothesis that fruitless is an ancient gene and has conserved functions in insects. Considerable divergences in other parts of the gene are expected to underlie species-specific differences and may help to explain diversity observed in insect sexual behaviors.

Allatostatin C and its paralog allatostatin double C: the arthropod somatostatins

Arthropods do not have one, but two genes encoding an allatostatin C-like peptide. The newly discovered paralog gene was called Ast-CC, and the peptide which it is predicted to make was called allatostatin double C (ASTCC). Genes for both allatostatin C (ASTC) and its paralog were found in the tick Ixodes scapularis as well as dipteran, lepidopteran, coleopteran, aphidoidean and phthirapteran insect species. In addition partial or complete cDNAs derived from Ast-CCs were found in a number of species, including Drosophila melanogaster, Bombyx mori and Rhodnius prolixus. The ASTCC precursors have a second conserved peptide sequence suggesting that they may produce two biologically active peptides. The predicted precursors encoded by the Ast-CCs have some unusual features, particularly in Drosophila, where they lack a signal peptide, and have instead a peptide anchor. These unusual structural features suggest that they are perhaps expressed by cells that are not specialized in neuropeptide synthesis and that in Drosophila ASTCC may be a juxtacrine. Data from the Fly Atlas project show that in Drosophila Ast-CC is little expressed. Nevertheless a P-element insertion in this gene is embryonic lethal, suggesting that it is an essential gene. Similarity between the precursors and receptors of ASTC/ASTCC and somatostatin suggests that ASTC/ASTCC and somatostatin have a common ancestor.

Resolving arthropod phylogeny: exploring phylogenetic signal within 41 kb of protein-coding nuclear gene sequence

This study attempts to resolve relationships among and within the four basal arthropod lineages (Pancrustacea, Myriapoda, Euchelicerata, Pycnogonida) and to assess the widespread expectation that remaining phylogenetic problems will yield to increasing amounts of sequence data. Sixty-eight regions of 62 protein-coding nuclear genes (approximately 41 kilobases (kb)/taxon) were sequenced for 12 taxonomically diverse arthropod taxa and a tardigrade outgroup. Parsimony, likelihood, and Bayesian analyses of total nucleotide data generally strongly supported the monophyly of each of the basal lineages represented by more than one species. Other relationships within the Arthropoda were also supported, with support levels depending on method of analysis and inclusion/exclusion of synonymous changes. Removing third codon positions, where the assumption of base compositional homogeneity was rejected, altered the results. Removing the final class of synonymous mutations--first codon positions encoding leucine and arginine, which were also compositionally heterogeneous--yielded a data set that was consistent with a hypothesis of base compositional homogeneity. Furthermore, under such a data-exclusion regime, all 68 gene regions individually were consistent with base compositional homogeneity. Restricting likelihood analyses to nonsynonymous change recovered trees with strong support for the basal lineages but not for other groups that were variably supported with more inclusive data sets. In a further effort to increase phylogenetic signal, three types of data exploration were undertaken. (1) Individual genes were ranked by their average rate of nonsynonymous change, and three rate categories were assigned--fast, intermediate, and slow. Then, bootstrap analysis of each gene was performed separately to see which taxonomic groups received strong support. Five taxonomic groups were strongly supported independently by two or more genes, and these genes mostly belonged to the slow or intermediate categories, whereas groups supported only by a single gene region tended to be from genes of the fast category, arguing that fast genes provide a less consistent signal. (2) A sensitivity analysis was performed in which increasing numbers of genes were excluded, beginning with the fastest. The number of strongly supported nodes increased up to a point and then decreased slightly. Recovery of Hexapoda required removal of fast genes. Support for Mandibulata (Pancrustacea + Myriapoda) also increased, at times to "strong" levels, with removal of the fastest genes. (3) Concordance selection was evaluated by clustering genes according to their ability to recover Pancrustacea, Euchelicerata, or Myriapoda and analyzing the three clusters separately. All clusters of genes recovered the three concordance clades but were at times inconsistent in the relationships recovered among and within these clades, a result that indicates that the a priori concordance criteria may bias phylogenetic signal in unexpected ways. In a further attempt to increase support of taxonomic relationships, sequence data from 49 additional taxa for three slow genes (i.e., EF-1 alpha, EF-2, and Pol II) were combined with the various 13-taxon data sets. The 62-taxon analyses supported the results of the 13-taxon analyses and provided increased support for additional pancrustacean clades found in an earlier analysis including only EF-1 alpha, EF-2, and Pol II.

Structural and evolutionary innovation of the heterodimerization interface between USP and the ecdysone receptor ECR in insects

Understanding how the variability of protein structure arises during evolution and leads to new structure-function relationships ultimately promoting evolutionary novelties is a major goal of molecular evolution and is critical for interpreting genome sequences. We addressed this issue using the ecdysone receptor (ECR), a major developmental factor that controls development and reproduction of arthropods. The functional ECR is a heterodimer of two nuclear receptors: ECR, which binds ecdysteroids, and its obligatory partner ultraspirade (USP), which is orthologous to the retinoid X receptor of vertebrates. Both genes underwent a dramatic increase of evolutionary rate in Mecopterida, the major insect terminal group containing Dipteras and Lepidopteras. We therefore questioned the implication of this event in terms of coevolution of their dimerization interface. A structural comparison revealed a 30% larger ligand-binding domain (LBD) heterodimerization surface in the Lepidoptera Heliothis when compared with basal insects, associated with a symmetrization of the interface, which is exceptional for nuclear receptors. Reconstruction of ancestral sequences and homology modeling of the ancestral Mecopterida ECR-USP reveal that this enlarged dimerization surface is a synapomorphy for Mecopterida. Furthermore, we show that the residues implicated in the new dimerization surface underwent specific evolutionary constraints in Mecopterida indicative of their new and conserved role in the dimerization interface. Most of all, the novel surface originates from a 15 degrees torsion of a subdomain of USP LBD toward its partner ECR, which is a long-range consequence of the peculiar position of a Mecopterida-specific insertion in loop L1-3, located outside of the interaction surface, in a less crucial domain of the partner protein. These results indicate that the coevolution between ECR and USP occurred through a novel mechanism of intramolecular epistasis that will undoubtedly be generalized for other molecules because it uses flexibility of a less-constrained region of a protein to modify the structure of another, critical part of the molecule.

The myofibrillar protein, projectin, is highly conserved across insect evolution except for its PEVK domain

All striated muscles respond to stretch by a delayed increase in tension. This physiological response, known as stretch activation, is, however, predominantly found in vertebrate cardiac muscle and insect asynchronous flight muscles. Stretch activation relies on an elastic third filament system composed of giant proteins known as titin in vertebrates or kettin and projectin in insects. The projectin insect protein functions jointly as a "scaffold and ruler" system during myofibril assembly and as an elastic protein during stretch activation. An evolutionary analysis of the projectin molecule could potentially provide insight into how distinct protein regions may have evolved in response to different evolutionary constraints. We mined candidate genes in representative insect species from Hemiptera to Diptera, from published and novel genome sequence data, and carried out a detailed molecular and phylogenetic analysis. The general domain organization of projectin is highly conserved, as are the protein sequences of its two repeated regions-the immunoglobulin type C and fibronectin type III domains. The conservation in structure and sequence is consistent with the proposed function of projectin as a scaffold and ruler. In contrast, the amino acid sequences of the elastic PEVK domains are noticeably divergent, although their length and overall unusual amino acid makeup are conserved. These patterns suggest that the PEVK region working as an unstructured domain can still maintain its dynamic, and even its three-dimensional, properties, without the need for strict amino acid conservation. Phylogenetic analysis of the projectin proteins also supports a reclassification of the Hymenoptera in relation to Diptera and Coleoptera.

Evolutionary origin and genomic organisation of runt-domain containing genes in arthropods

BACKGROUND

Gene clusters, such as the Hox gene cluster, are known to have critical roles in development. In eukaryotes gene clusters arise primarily by tandem gene duplication and divergence. Genes within a cluster are often co-regulated, providing selective pressure to maintain the genome organisation, and this co-regulation can result in temporal or spatial co-linearity of gene expression. It has been previously noted that in Drosophila melanogaster, three of the four runt-domain (RD) containing genes are found in a relatively tight cluster on chromosome 1, raising the possibility of a putative functional RD gene cluster in D. melanogaster.

RESULTS

To investigate the possibility of such a gene cluster, orthologues of the Drosophila melanogaster RD genes were identified in several endopterygotan insects, two exopterygotan insects and two non-insect arthropods. In all insect species four RD genes were identified and orthology was assigned to the Drosophila sequences by phylogenetic analyses. Although four RD genes were found in the crustacean D. pulex, orthology could not be assigned to the insect sequences, indicating independent gene duplications from a single ancestor following the split of the hexapod lineage from the crustacean lineage.In insects, two chromosomal arrangements of these genes was observed; the first a semi-dispersed cluster, such as in Drosophila, where lozenge is separated from the core cluster of three RD genes often by megabases of DNA. The second arrangement was a tight cluster of the four RD genes, such as in Apis mellifera.This genomic organisation, particularly of the three core RD genes, raises the possibility of shared regulatory elements. In situ hybridisation of embryonic expression of the four RD genes in Drosophila melanogaster and the honeybee A. mellifera shows no evidence for either spatial or temporal co-linearity of expression during embryogenesis.

CONCLUSION

All fully sequenced insect genomes contain four RD genes and orthology can be assigned to these genes based on similarity to the D. melanogaster protein sequences. Examination of the genomic organisation of these genes provides evidence for a functional RD gene cluster. RD genes from non-insect arthropods are also clustered, however the lack of orthology between these and insect RD genes suggests this cluster is likely to have resulted from a duplication event independent from that which created the insect RD gene cluster. Analysis of embryonic RD gene expression in two endopterygotan insects, A. mellifera and D. melanogaster, did not show evidence for coordinated gene expression, therefore while the functional significance of this gene cluster remains unknown its maintenance during insect evolution implies some functional significance to the cluster.

Heads and tails: evolution of antero-posterior patterning in insects

In spite of their varied appearances, insects share a common body plan whose layout is established by patterning genes during embryogenesis. We understand in great molecular detail how the Drosophila embryo patterns its segments. However, Drosophila has a type of embryogenesis that is highly derived and varies extensively as compared to most insects. Therefore, the study of other insects is invaluable for piecing together how the ancestor of all insects established its segmented body plan, and how this process can be plastic during evolution. In this review, we discuss the evolution of Antero-Posterior (A-P) patterning mechanisms in insects. We first describe two distinct modes of insect development - long and short germ development - and how these two modes of patterning are achieved. We then summarize how A-P patterning occurs in the long-germ Drosophila, where most of our knowledge comes from, and in the well-studied short-germ insect, Tribolium. Finally, using examples from other insects, we highlight differences in patterns of expression, which suggest foci of evolutionary change.

Mitochondrial genome organization and phylogeny of two vespid wasps

We sequenced the entire mitochondrial genome of Abispa ephippium (Hymenoptera: Vespoidea: Vespidae: Eumeninae) and most of the mitochondrial genome of Polistes humilis synoecus (Hymenoptera: Vespoidea: Vespidae: Polistinae). The arrangement of genes differed between the two genomes and also differed slightly from that inferred to be ancestral for the Hymenoptera. The genome organization for both vespids is different from that of all other mitochondrial genomes previously reported. A number of tRNA gene rearrangements were identified that represent potential synapomorphies for a subset of the Vespidae. Analysis of all available hymenopteran mitochondrial genome sequences recovered an uncontroversial phylogeny, one consistent with analyses of other types of data.

Characterization of the multicopper oxidase gene family in Anopheles gambiae

The multicopper oxidase (MCO) family of enzymes includes laccases, which oxidize a broad range of substrates including diphenols, and several oxidases with specific substrates such as iron, copper or ascorbic acid. We have identified five putative MCO genes in the genome of Anopheles gambiae and have cloned cDNAs encompassing the full coding region for each gene. MCO1 mRNA was detected in all developmental stages and in all of the larval and adult tissues tested. We observed an increase in MCO1 transcript abundance in the midguts and Malphighian tubules of adult females following a blood meal and in adult abdominal carcasses in response to an immune challenge. Two alternatively spliced isoforms of MCO2 mRNA were identified. The A isoform of MCO2 was previously detected in larval and pupal cuticle where it probably catalyzes sclerotization reactions (He, N., Botelho, J.M.C., McNall, R.J., Belozerov, V., Dunn, W.A., Mize, T., Orlando, R., Willis, J.H., 2007. Proteomic analysis of cast cuticles from Anopheles gambiae by tandem mass spectrometry. Insect Biochem. Mol. Biol. 37, 135-146). The B isoform was transcriptionally upregulated in ovaries in response to a blood meal. MCO3 mRNA was detected in the adult midgut, Malpighian tubules, and male reproductive tissues; like MCO1, it was upregulated in response to an immune challenge or a blood meal. MCO4 and MCO5 were observed primarily in eggs and in the abdominal carcass of larvae. A phylogenetic analysis of insect MCO genes identified putative orthologs of MCO1 and MCO2 in all of the insect genomes tested, whereas MCO3, MCO4 and MCO5 were found only in the two mosquito species analyzed. MCO2 orthologs have especially high sequence similarity, suggesting that they are under strong purifying selection; the A isoforms are more conserved than the B isoforms. The mosquito specific group shares a common ancestor with MCO2. This initial study of mosquito MCOs suggests that MCO2 may be required for egg development or eggshell tanning in addition to cuticle tanning, while MCO1 and MCO3 may be involved in metal metabolism or immunity.

Conserved cluster organization of insect Runx genes

Runx gene family transcription factors play important regulatory roles in metazoan development. The Drosophila genome contains four Runx genes, two of which are well studied (runt and lozenge) and two of which have not been explored in depth yet (CG42267 and CG34145). In the absence of paralog-specific orthologs in vertebrates, we investigated the evolutionary conservation of the Drosophila Runx genes in the genomes of mosquito, red flour beetle, honeybee, and distantly related Bilateria. All four Drosophila Runx genes are conserved in other endopterygote insects. Strikingly, their genetic linkage and genomic arrangement is also highly conserved despite past recombination events in the same genomic region. The preservation of Runx gene cluster organization for at least 300 million years resembles that of insect Hox and selected Wnt genes. We propose a model for the conserved cluster organization of developmental gene family paralogs based on differential coduplication of regulatory elements that act over long distances.

The evolution of developmental gene networks: lessons from comparative studies on holometabolous insects

Recent comparative studies have revealed significant differences in the developmental gene networks operating in three holometabolous insects: the beetle Tribolium castaneum, the parasitic wasp Nasonia vitripennis and the fruitfly Drosophila melanogaster. I discuss these differences in relation to divergent and convergent changes in cellular embryology. I speculate on how segmentation gene networks have evolved to operate in divergent embryological contexts, and highlight the role that co-option might have played in this process. I argue that insects represent an important example of how diversification in life-history strategies between lineages can lead to divergence in the genetic and cellular mechanisms controlling the development of homologous adult structures.