Divergent decapentaplegic expression patterns in compound eye development and the evolution of insect metamorphosis

In the fruit fly Drosophila, the patterning genes decapentaplegic and wingless contribute to the spatial control of retina development in an antagonistic manner. We examined the expression patterns of these genes in the developing visual system of the hemimetabolous grasshopper Schistocerca americana and the primitive holometabolous beetle species Tribolium castaneum. The pattern of wingless expression was strongly conserved as a pair of lateral domains at the anterior margins of both the developing retina and the developing optic lobes. The expression of decapentaplegic, on the other hand, is different. Unlike in Drosophila, no decapentaplegic expression was detected before the onset of photoreceptor differentiation in the retinal precursor tissue of either grasshopper or beetle. Moreover, the subsequent expression of decapentaplegic in the latter species was not concentrated in the moving front of retina differentiation, as in Drosophila, but observed in anterior and posterior regions. Our results indicate that Drosophila eye development contains elements of both ancestral and derived regulatory gene functions. The requirement for decapentaplegic as an antagonist of wingless during the early development of the Drosophila retina might have originated during the evolution of insect metamorphosis.     

Molecular evolution and expression of the CRAL_TRIO protein family in insects

CRAL_TRIO domain proteins are known to bind small lipophilic molecules such as retinal, inositol and Vitamin E and include such gene family members as PINTA, α-tocopherol transfer (ATT) proteins, retinoid binding proteins, and clavesins. In insects, very little is known about either the molecular evolution of this family of proteins or their ligand specificity. Here we characterize insect CRAL_TRIO domain proteins and present the first insect CRAL_TRIO protein phylogeny constructed by performing reciprocal BLAST searches of the reference genomes of Drosophila melanogaster, Anopheles gambiae, Apis mellifera, Tribolium castaneum, Bombyx mori, Manduca sexta and Danaus plexippus. We find several highly conserved amino acid residues in the CRAL_TRIO domain-containing genes across insects and a gene expansion resulting in more than twice as many gene family members in lepidopterans than in other surveyed insect species, but no lepidopteran homolog of the PINTA gene in Drosophila. In addition, we examined the expression pattern of CRAL_TRIO domain genes in Manduca sexta heads using RNA-Seq data. Of the 42 gene family members found in the M. sexta reference genome, we found 30 expressed in the head tissue with similar expression profiles between males and females. Our results suggest this gene family underwent a large expansion in Lepidoptera, making the lepidopteran CRAL_TRIO domain family distinct from other holometabolous insect lineages.     

Concerted evolution within the Drosophila dumpy gene

We have determined by reverse Southern analysis and direct sequence comparisons that most of the dumpy gene has evolved in the dipteran and other insect orders by purifying selection acting on amino acid replacements. One region, however, is evolving rapidly due to unequal crossing over and/or gene conversion. This region, called "PIGSFEAST," or PF, encodes in D. melanogaster 30-47 repeats of 102 amino acids rich in serines, threonines, and prolines. We show that the processes of concerted evolution have been operating on all species of Drosophila examined to date, but that an adjacent region has expanded in Anopheles gambiae, Aedes aegypti, and Tribolium castaneum, while the PF repeats are reduced in size and number. In addition, processes of concerted evolution have radically altered the codon usage patterns in D. melanogaster, D. pseudoobscura, and D. virilis compared with the rest of the dumpy gene. We show also that the dumpy gene is expressed on the inner surface of the micropyle of the mature oocyte and propose that, as in the abalone system, concerted evolution may be involved in adaptive changes affecting Dumpy's possible role in sperm-egg recognition.     

Identifying concerted evolution and gene conversion in mammalian gene pairs lasting over 100 million years

Background  

Concerted evolution occurs in multigene families and is characterized by stretches of homogeneity and higher sequence similarity between paralogues than between orthologues. Here we identify human gene pairs that have undergone concerted evolution, caused by ongoing gene conversion, since at least the human-mouse divergence. Our strategy involved the identification of duplicated genes with greater similarity within a species than between species. These genes were required to be present in multiple mammalian genomes, suggesting duplication early in mammalian divergence. To eliminate genes that have been conserved due to strong purifying selection, our analysis also required at least one intron to have retained high sequence similarity between paralogues.     

Diversity in insect axis formation: two orthodenticle genes and hunchback act in anterior patterning and influence dorsoventral organization in the honeybee (Apis mellifera)

Axis formation is a key step in development, but studies indicate that genes involved in insect axis formation are relatively fast evolving. Orthodenticle genes have conserved roles, often with hunchback, in maternal anterior patterning in several insect species. We show that two orthodenticle genes, otd1 and otd2, and hunchback act as maternal anterior patterning genes in the honeybee (Apis mellifera) but, unlike other insects, act to pattern the majority of the anteroposterior axis. These genes regulate the expression domains of anterior, central and posterior gap genes and may directly regulate the anterior gap gene giant. We show otd1 and hunchback also influence dorsoventral patterning by regulating zerknült (zen) as they do in Tribolium, but that zen does not regulate the expression of honeybee gap genes. This suggests that interactions between anteroposterior and dorsal-ventral patterning are ancestral in holometabolous insects. Honeybee axis formation, and the function of the conserved anterior patterning gene orthodenticle, displays unique characters that indicate that, even when conserved genes pattern the axis, their regulatory interactions differ within orders of insects, consistent with relatively fast evolution in axis formation pathways.     

Pax group III genes and the evolution of insect pair-rule patterning

Pair-rule genes were identified and named for their role in segmentation in embryos of the long germ insect Drosophila. Among short germ insects these genes exhibit variable expression patterns during segmentation and thus are likely to play divergent roles in this process. Understanding the details of this variation should shed light on the evolution of the genetic hierarchy responsible for segmentation in Drosophila and other insects. We have investigated the expression of homologs of the Drosophila Pax group III genes paired, gooseberry and gooseberry-neuro in short germ flour beetles and grasshoppers. During Drosophila embryogenesis, paired acts as one of several pair-rule genes that define the boundaries of future parasegments and segments, via the regulation of segment polarity genes such as gooseberry, which in turn regulates gooseberry-neuro, a gene expressed later in the developing nervous system. Using a crossreactive antibody, we show that the embryonic expression of Pax group III genes in both the flour beetle Tribolium and the grasshopper Schistocerca is remarkably similar to the pattern in Drosophila. We also show that two Pax group III genes, pairberry1 and pairberry2, are responsible for the observed protein pattern in grasshopper embryos. Both pairberry1 and pairberry2 are expressed in coincident stripes of a one-segment periodicity, in a manner reminiscent of Drosophila gooseberry and gooseberry-neuro. pairberry1, however, is also expressed in stripes of a two-segment periodicity before maturing into its segmental pattern. This early expression of pairberry1 is reminiscent of Drosophila paired and represents the first evidence for pair-rule patterning in short germ grasshoppers or any hemimetabolous insect.     

Comparative genomic analysis of the <Emphasis Type="Italic">Tribolium</Emphasis> immune system

Background  

Tribolium castaneum is a species of Coleoptera, the largest and most diverse order of all eukaryotes. Components of the innate immune system are hardly known in this insect, which is in a key phylogenetic position to inform us about genetic innovations accompanying the evolution of holometabolous insects. We have annotated immunity-related genes and compared them with homologous molecules from other species.     

Isolation of an abdominal-A gene from the locust Schistocerca gregaria and its expression during early embryogenesis

Using sequence homology to Drosophila homeobox-containing genes, we have cloned a homologue of abdominal-A from the locust Schistocerca gregaria. The Schistocerca clone encodes a stretch of 78 amino acids including the homeodomain and its flanking regions identical to the corresponding region of abdominal-A. We have shown by in situ hybridization that this gene is transcribed and have used an antibody raised against its protein product to examine the expression of abdominal-A during early Schistocerca embryogenesis. Schistocerca is a short germ insect. Although the segmented body plan is very similar to that of Drosophila, the segments are generated sequentially by a process of growth, not simultaneously by subdivision of a syncytial blastoderm. In both organisms, abdominal-A is expressed throughout the abdomen from a sharp anterior boundary located within the first abdominal segment (A1). The initial activation of the genes in the two species differs. Schistocerca initiates expression in a small group of cells in the anterior of A2, shortly after this segment is defined by the appearance of engrailed protein. This contrasts with the appearance of abdominal-A expression in Drosophila, which appears simultaneously throughout the entire abdomen.     

The expression and function of the achaete-scute genes in Tribolium castaneum reveals conservation and variation in neural pattern formation and cell fate specification

The study of achaete-scute (ac/sc) genes has recently become a paradigm to understand the evolution and development of the arthropod nervous system. We describe the identification and characterization of the ac/sc genes in the coleopteran insect species Tribolium castaneum. We have identified two Tribolium ac/sc genes - achaete-scute homolog (Tc-ASH) a proneural gene and asense (Tc-ase) a neural precursor gene that reside in a gene complex. Focusing on the embryonic central nervous system we find that Tc-ASH is expressed in all neural precursors and the proneural clusters from which they segregate. Through RNAi and misexpression studies we show that Tc-ASH is necessary for neural precursor formation in Tribolium and sufficient for neural precursor formation in Drosophila. Comparison of the function of the Drosophila and Tribolium proneural ac/sc genes suggests that in the Drosophila lineage these genes have maintained their ancestral function in neural precursor formation and have acquired a new role in the fate specification of individual neural precursors. Furthermore, we find that Tc-ase is expressed in all neural precursors suggesting an important and conserved role for asense genes in insect nervous system development. Our analysis of the Tribolium ac/sc genes indicates significant plasticity in gene number, expression and function, and implicates these modifications in the evolution of arthropod neural development.     

The role of Broad in the development of Tribolium castaneum: implications for the evolution of the holometabolous insect pupa

The evolution of complete metamorphosis in insects is a key innovation that has led to the successful diversification of holometabolous insects, yet the origin of the pupa remains an enigma. Here, we analyzed the expression of the pupal specifier gene broad (br), and the effect on br of isoform-specific, double-stranded RNA-mediated silencing, in a basal holometabolous insect, the beetle Tribolium castaneum. All five isoforms are weakly expressed during the penultimate instar and highly expressed during the prepupal period of the final instar. Application of hydroprene, a juvenile hormone analog, during the penultimate instar caused a repeat of the penultimate br expression patterns, and the formation of supernumerary larvae. Use of dsRNA against the br core region, or against a pair of either the br-Z2 or br-Z3 isoform with the br-Z1 or br-Z4 isoform, produced mobile animals with well-differentiated adult-like appendages, but which retained larval-like urogomphi and epidermis. Disruption of either the br-Z2 or the br-Z3 isoform caused the formation of shorter wings. Disruption of both br-Z1 and br-Z4 caused the appearance of pupal traits in the adults, but disruption of br-Z5 had no morphological effect. Our findings show that the br isoform functions are broadly conserved within the Holometabola and suggest that evolution of br isoform expression may have played an important role in the evolution of the pupa in holometabolous insects.     

The evolution of the adenylate-forming protein family in beetles: multiple luciferase gene paralogues in fireflies and glow-worms

Bioluminescence in beetles is dependent upon the enzyme luciferase. It has been hypothesised luciferase evolved from a fatty acyl-CoA synthetase gene deriving a novel bioluminescent function (neofunctionalization) after a gene duplication event. We evaluated this hypothesis within a phylogenetic framework using independent evidence obtained from the genome of Tribolium castaneum, published luciferase genes and novel luciferase and luciferase-like sequences. This phylogenetic study provides evidence for a large gene family of luciferase and luciferase-like paralogues in bioluminescent and non-bioluminescent beetles. All luciferase sequences formed a clade supporting a protoluciferase existing prior to the divergence of the Lampyridae, Elateridae and Phengodidae (Elateroidea). Multiple luciferase genes were identified from members of the Photurinae and the Luciolinae indicating complex gene duplication events within lampyrid genomes. The majority of luciferase residues were identified to be under purifying selection as opposed to positive selection. We conclude that beetle luciferase may have arisen from a process of subfunctionalization as opposed to neofunctionalization early on in the evolution of the Elateroidea.     

Sex-specific splicing of the honeybee doublesex gene reveals 300 million years of evolution at the bottom of the insect sex-determination pathway

Sex-determination mechanisms vary greatly among taxa. It has been proposed that genetic sex-determination pathways evolve in reverse order from the final step in the pathway to the first step. Consistent with this hypothesis, doublesex (dsx), the most downstream gene in the Drosophila sex-determination cascade that determines most sexual phenotypes also determines sex in other dipterans and the silk moth, while the upstream genes vary among these species. However, it is unknown when dsx was recruited to the sex-determination pathway during insect evolution. Furthermore, sex-specific splicing of dsx, by which dsx determines sex, is different in pattern and mechanism between the moth and the fly, raising an interesting question of how these insects have kept the executor of sex determination while allowing flexibility in the means of execution. To address these questions, here we study the dsx gene of the honeybee Apis mellifera, a member of the most basal lineage of holometabolous insects. We report that honeybee dsx is sex-specifically spliced and that it produces both the fly-type and moth-type splicing forms, indicating that the use of different splicing forms of Dsx in controlling sexual differentiation was present in the common ancestor of holometabolous insects. Our data suggest that in ancestral holometabolous insects the female Dsx form is the default and the male form is generated by suppressing the splicing of the female form. Thus, it is likely that the dsx splicing activator system in flies, where the male form is the default, arose during early dipteran evolution.     

Purifying selection is a prevailing motif in the evolution of ketoacyl synthase domains of polyketide synthases from lichenized fungi

We analysed ketoacyl synthase domains of type I polyketide synthase (PKS) gene fragments of 163 lichenized and 51 non-lichenized fungi in a Bayesian phylogenetic framework. Lichenized taxa from several unrelated taxonomic groups, some of which produce identical secondary metabolites, were included. We found 12 clades of non-reducing PKS genes, which represent monophyletic PKS paralogues. PAML and SELECTON analyses indicated that purifying selection is the prevailing selective force in the evolution of the keto synthase domain of these paralogues. We detected no unambiguous correlation between PKS clades and the distribution of lichen substances. Together with the strong evidence for purifying selection, the wide distribution of certain paralogues in ascomycetes suggested early gene duplication events in the evolutionary history of this gene family in the Ascomycota.     

Accelerated evolution of mitochondrial but not nuclear genomes of Hymenoptera: new evidence from crabronid wasps

Mitochondrial genes in animals are especially useful as molecular markers for the reconstruction of phylogenies among closely related taxa, due to the generally high substitution rates. Several insect orders, notably Hymenoptera and Phthiraptera, show exceptionally high rates of mitochondrial molecular evolution, which has been attributed to the parasitic lifestyle of current or ancestral members of these taxa. Parasitism has been hypothesized to entail frequent population bottlenecks that increase rates of molecular evolution by reducing the efficiency of purifying selection. This effect should result in elevated substitution rates of both nuclear and mitochondrial genes, but to date no extensive comparative study has tested this hypothesis in insects. Here we report the mitochondrial genome of a crabronid wasp, the European beewolf (Philanthus triangulum, Hymenoptera, Crabronidae), and we use it to compare evolutionary rates among the four largest holometabolous insect orders (Coleoptera, Diptera, Hymenoptera, Lepidoptera) based on phylogenies reconstructed with whole mitochondrial genomes as well as four single-copy nuclear genes (18S rRNA, arginine kinase, wingless, phosphoenolpyruvate carboxykinase). The mt-genome of P. triangulum is 16,029 bp in size with a mean A+T content of 83.6%, and it encodes the 37 genes typically found in arthropod mt genomes (13 protein-coding, 22 tRNA, and two rRNA genes). Five translocations of tRNA genes were discovered relative to the putative ancestral genome arrangement in insects, and the unusual start codon TTG was predicted for cox2. Phylogenetic analyses revealed significantly longer branches leading to the apocritan Hymenoptera as well as the Orussoidea, to a lesser extent the Cephoidea, and, possibly, the Tenthredinoidea than any of the other holometabolous insect orders for all mitochondrial but none of the four nuclear genes tested. Thus, our results suggest that the ancestral parasitic lifestyle of Apocrita is unlikely to be the major cause for the elevated substitution rates observed in hymenopteran mitochondrial genomes.     

Drosophila fushi tarazu. a gene on the border of homeotic function

BACKGROUND: Hox genes specify cell fate and regional identity during animal development. These genes are present in evolutionarily conserved clusters thought to have arisen by gene duplication and divergence. Most members of the Drosophila Hox complex (HOM-C) have homeotic functions. However, a small number of HOM-C genes, such as the segmentation gene fushi tarazu (ftz), have nonhomeotic functions. If these genes arose from a homeotic ancestor, their functional properties must have changed significantly during the evolution of modern Drosophila. RESULTS: Here, we have asked how Drosophila ftz evolved from an ancestral homeotic gene to obtain a novel function in segmentation. We expressed Ftz proteins at various developmental stages to assess their potential to regulate segmentation and to generate homeotic transformations. Drosophila Ftz protein has lost the inherent ability to mediate homeosis and functions exclusively in segmentation pathways. In contrast, Ftz from the primitive insect Tribolium (Tc-Ftz) has retained homeotic potential, generating homeotic transformations in larvae and adults and retaining the ability to repress homothorax, a hallmark of homeotic genes. Similarly, Schistocerca Ftz (Sg-Ftz) caused homeotic transformations of antenna toward leg. Primitive Ftz orthologs have moderate segmentation potential, reflected by weak interactions with the segmentation-specific cofactor Ftz-F1. Thus, Ftz orthologs represent evolutionary intermediates that have weak segmentation potential but retain the ability to act as homeotic genes. CONCLUSIONS: ftz evolved from an ancestral homeotic gene as a result of changes in both regulation of expression and specific alterations in the protein-coding region. Studies of ftz orthologs from primitive insects have provided a "snap-shot" view of the progressive evolution of a Hox protein as it took on segmentation function and lost homeotic potential. We propose that the specialization of Drosophila Ftz for segmentation resulted from loss and gain of specific domains that mediate interactions with distinct cofactors.     

A genome-wide analysis of genes and gene families involved in innate immunity of Bombyx mori

A genome-wide analysis of innate immunity-related genes and gene families was conducted using the silkworm, Bombyx mori. We identified orthologs for a large number of genes involved in insect immunity that have been reported from Drosophila melanogaster (Diptera), Anopheles gambiae (Diptera), Apis mellifera (Hymenoptera) and Tribolium castaneum (Coleoptera). B. mori has a unique recognition gene and antimicrobial peptide genes that are not present in the Drosophila, Anopheles, Apis and Tribolium genomes, suggesting a lineage-specific gene evolution for lepidopteran insects. The comparative analysis of the insect immune repertoires indicated a dynamic and flexible gene expansion in recognition, modulation and effector mechanisms due to different selection pressures. Differential gene regulation by different bacterial species was found in PGRP and Serpin genes, suggesting that Bombyx has a highly selective gene regulation system depending on bacterial species.     

Elongation factor-1 alpha occurs as two copies in bees: implications for phylogenetic analysis of EF-1 alpha sequences in insects

We report the complete sequence of a paralogous copy of elongation factor-1 alpha (EF-1 alpha) in the honeybee, Apis mellifera (Hymenoptera: Apidae). This copy differs from a previously described copy in the positions of five introns and in 25% of the nucleotide sites in the coding regions. The existence of two paralogous copies of EF-1 alpha in Drosophila and Apis suggests that two copies of EF-1 alpha may be widespread in the holometabolous insect orders. To distinguish between a single, ancient gene duplication and parallel, independent fly and bee gene duplications, we performed a phylogenetic analysis of hexapod EF-1 alpha sequences. Unweighted parsimony analysis of nucleotide sequences suggests an ancient gene duplication event, whereas weighted parsimony analysis of nucleotides and unweighted parsimony analysis of amino acids suggests the contrary: that EF-1 alpha underwent parallel gene duplications in the Diptera and the Hymenoptera. The hypothesis of parallel gene duplication is supported both by congruence among nucleotide and amino acid data sets and by topology-dependent permutation tail probability (T-PTP) tests. The resulting tree topologies are also congruent with current views on the relationships among the holometabolous orders included in this study (Diptera, Hymenoptera, and Lepidoptera). More sequences, from diverse orders of holometabolous insects, will be needed to more accurately assess the historical patterns of gene duplication in EF-1 alpha.