Advances in insect phylogeny at the dawn of the postgenomic era

A global overview of insect-fern interactions and its ecological trends

Historically, ferns have been described as underutilized by insects. However, studies have shown a diversity of insects interacting with ferns, although the evolutionary and ecological drivers of these interactions are still to be untangled. To fill these gaps, we compiled more than 100 yr of global data on insect-fern interactions from the literature comprising 374 fern and 649 insect species. With this database we assessed how fern trophic specialization, phylogenetic relationships and climate have shaped their interactions with insects. Our findings showed that interactions between ferns and insects can be explained by the phylogenetic relations among them. We observed that insect orders part of the Endopterygota clade tend to interact with similar fern species, which might be a result of the inheritance of Endopterygota ancestors probably due to phylogenetic niche conservationism. Under an ecological context, fern specialization increased with temperature, precipitation, and climatic stability. Our results show that climate might be one of the main factors explaining the spatial variation of insect-fern interactions, postulate also supported by the observed phylogenetic clustering of the studied ferns species. Our study highlights the intricate and multifaceted nature of insect-fern interactions, where evolutionary history and ecological factors converge to shape these relationships.

Revisiting the four Hexapoda classes: Protura as the sister group to all other hexapods

Insects represent the most diverse animal group, yet previous phylogenetic analyses based on morphological and molecular data have failed to agree on the evolutionary relationships of early insects and their six-legged relatives (together constituting the clade Hexapoda). In particular, the phylogenetic positions of the three early-diverging hexapod lineages-the coneheads (Protura), springtails (Collembola), and two-pronged bristletails (Diplura)-have been debated for over a century, with alternative topologies implying drastically different scenarios of the evolution of the insect body plan and hexapod terrestrialization. We addressed this issue by sampling all hexapod orders and experimenting with a broad range of across-site compositional heterogeneous models designed to tackle ancient divergences. Our analyses support Protura as the earliest-diverging hexapod lineage ("Protura-sister") and Collembola as a sister group to Diplura, a clade corresponding to the original composition of Entognatha, and characterized by the shared possession of internal muscles in the antennal flagellum. The previously recognized 'Elliplura' hypothesis is recovered only under the site-homogeneous substitution models with partial supermatrices. Our cross-validation analysis shows that the site-heterogeneous CAT-GTR model, which recovers "Protura-sister," fits significantly better than homogeneous models. Furthermore, the morphologically unusual Protura are also supported as the earliest-diverging hexapod lineage by other lines of evidence, such as mitogenomes, comparative embryology, and sperm morphology, which produced results similar to those in this study. Our backbone phylogeny of hexapods will facilitate the exploration of the underpinnings of hexapod terrestrialization and megadiversity.

Evolutionary and phylogenetic insights from the mitochondrial genomic analysis of Diceraeus melacanthus and D. furcatus (Hemiptera: Pentatomidae)

The mitochondrial genomes of D. melacanthus and D. furcatus were sequenced and used to investigate the phylogenetic relationships with 54 species of Pentatomidae. Their mitogenomes are 17,197 and 15,444 bp-long, respectively, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22/21 transfer RNA genes, with conserved gene arrangement. Leu, Lys, and Ser were the most common amino acids in their PCGs. PCGs evolutionary analysis indicated their mitogenomes are under purifying selection, and the most conserved genes are from the cytochrome complex, reinforcing their suitability as markers for molecular taxonomy. We identified 490 mtSSRs in 56 Pentatomidae species, with large variation and a positive correlation between mtSSR number and genome size. Three mtSSRs were identified in each Diceraeus species. Only the mtSSR in the nad6 (D. melacanthus) and nad4 (D. furcatus) appear to have application as molecular markers for species characterization. Phylogenetic analysis confirmed the monophyly of Pentatomidae. However, our analysis challenged the monophyly of Pentatominae and Podopinae. We also detected unexpected relationships among some tribes and genera, highlighting the complexity of the internal taxonomic structure of Pentatomidae. Both Diceraeus species were grouped in the same clade with the remaining Carpocorini analyzed.

Clogmia albipunctata (Williston, 1893) midgut physiology: pH control and functional relationship with Lower Diptera (nematoceran) especially with hematophagous species

Clogmia albipunctata (Williston, 1893) is a non-hematophagous insect belonging to the order Diptera, suborder Nematocera (Lower Diptera) and family Psychodidae. In the present work, we investigated how C. albipunctata control their midgut pH under different physiological conditions, comparing their midgut physiology with some nematoceran hematophagous species. The C. albipunctata midgut pH was measured after ingestion of sugar, protein and under the effect of the alkalinizing hormone released in the hemolymph of the hematophagous sand fly Lutzomyia longipalpis obtained just after a blood meal. The midgut pH of unfed or sugar-fed C. albipunctata is 5.5-6, and its midgut underwent alkalinization after protein ingestion or under treatment with hemolymph collected from blood fed L. longipalpis. These results suggested that in nematocerans, mechanisms for pH control seem shared between hematophagous and non-hematophagous species. This kind of pH control is convenient for successful blood digestion. The independent evolution of many hematophagous groups from the Lower Diptera suggests that characteristics involved in midgut pH control were already present in non-hematophagous species and represent a readiness for adaptation to this feeding mode.

Characterization of the complete mitochondrial genome of Paecilomyces variotii and comparative evolutionary mitochondriomics of 36 fungi

BACKGROUD

Paecilomyces variotii has important economic value in stimulating crop growth, biodegradation, and other aspects. Up to now, there are no research reports on its mitochondrial genome.

METHODS AND RESULTS

The mitochondrial genome of Paecilomyces variotii was determined with the next-generation sequencing method (Illumina, NovaSeq), and its characteristics were analyzed using various bioinformatics approaches. The length of complete mitochondrial genome sequence of P. variotii is 40,965 bp and consists of 14 protein-coding genes, 2 ribosomal RNA genes, 1 ribosomal protein S3 gene, 26 transport RNA genes. The results of phylogenetics analysis using Bayesian inference and Maximum likelihood methods showed that P. variotii belongs to the Eurotiales order in the Thermoascaceae family, and 9 genera within the Eurotiomycetes class were effectively distinguished with high support rates (bootstrap value > 92% and posterior probabilities > 99%). The analysis of synonymous substitution rates and nonsynonymous substitution rates indicated that the Ka/Ks values of the 14 PCGs in the mitochondrial genomes of the two orders in the Eurotiomycetes class ranged from 0 to 0.4333.

CONCLUSIONS

This study revealed the structural and sequence information characteristics of the mitochondrial genome of P. variotii, and the phylogenetic results strongly support its classification within the family Thermoascaceae, consistent with traditional morphological taxonomy studies. The 14 PCGs in the mitochondrial genomes of the two orders in the Eurotiomycetes class are subject to strong purifying (negative) selection. The results of this research provides an important molecular basis for the development of genomics, evolutionary genetics and molecular markers of P. variotii in the future.

Sensory pathway in aquatic basal polyneoptera: Antennal sensilla and brain morphology in stoneflies

Aquatic insects represent a great portion of Arthropod diversity and the major fauna in inland waters. The sensory biology and neuroanatomy of these insects are, however, poorly investigated. This research aims to describe the antennal sensilla of nymphs of the stonefly Dinocras cephalotes using scanning electron microscopy and comparing them with the adult sensilla. Besides, central antennal pathways in nymphs and adults are investigated by neuron mass-tracing with tetramethylrhodamine, and their brain structures are visualized with an anti-synapsin antibody. No dramatic changes occur in the antennal sensilla during nymphal development, while antennal sensilla profoundly change from nymphs to adults when switching from an aquatic to an aerial lifestyle. However, similar brain structures are used in nymphs and adults to process diverging sensory information, perceived through different sensilla in water and air. These data provide valuable insights into the evolution of aquatic heterometabolous insects, maintaining a functional sensory system throughout development, including a distinct adaptation of the peripheral olfactory systems during the transition from detection of water-soluble chemicals to volatile compounds in the air. From a conservation biology perspective, the present data contribute to a better knowledge of the biology of stoneflies, which are very important bioindicators in rivers.

Evolutionary genomics of camouflage innovation in the orchid mantis

The orchid mantises achieve camouflage with morphological modifications in body color and pattern, providing an interesting model for understanding phenotypic innovation. However, a reference genome is lacking for the order Mantodea. To unveil the mechanisms of plant-mimicking body coloration and patterns, we performed de novo assembly of two chromosome-level genomes of the orchid mantis and its close relative, the dead leaf mantis. Comparative genomic analysis revealed that the Scarlet gene plays an important role in the synthesis of xanthommatin, an important pigment for mantis camouflage coloration. Combining developmental transcriptomic analysis and genetic engineering experiments, we found that the cuticle was an essential component of the 'petal-like' enlargement, and specific expression in the ventral femur was controlled by Wnt signaling. The prolonged expression of Ultrabithorax (Ubx) accompanied by femoral expansion suggested that Ubx determines leg remodeling in the early developmental stage. We also found evidence of evolution of the Trypsin gene family for insectivory adaptation and ecdysone-dependent sexual dimorphism in body size. Overall, our study presents new genome catalogs and reveals the genetic and evolutionary mechanisms underlying the unique camouflage of the praying mantis, providing evolutionary developmental insights into phenotypic innovation and adaptation.

Evolutionary history of Nasutitermes kemneri (Termitidae, Nasutitermitinae), a termite from the South American diagonal of open formations

Little is known about the phylogeography of termites in the Neotropical region. Here, we explored the genetic patterns and phylogeographical processes in the evolutionary history of Nasutitermes kemneri, an endemic termite of the South American diagonal of open formations (DOF) formed by the Chaco, Cerrado, and Caatinga phytogeographic domains. We sampled 60 individuals across the three domains of the DOF, and using the mitochondrial genes 16S, COI, and COII, as well as the nuclear gene ITS, evaluated the genetic diversity and divergence time of the populations, along with their genetic structure. The results show a strong genetic and spatial structure within the samples, evidencing the existence of two well-differentiated genetic groups: the Northeastern and the Southwestern populations, which diverged about 2.5 Mya, during the Pliocene-Pleistocene boundary. The Northeastern population, which encompasses Caatinga and northern portions of Cerrado, has an intricate structure and seems to have suffered repetitive retraction-expansion events due to climactic fluctuations during the Quaternary. The Southwestern population, which ranges from central-south Cerrado to the northeast peripherical portions of the Chaco, displays a star-shaped haplotype structure, indicating that this region may have acted as a refugia during interglacial periods.

Characterizing Hox genes in mayflies (Ephemeroptera), with Hexagenia limbata as a new mayfly model

BACKGROUND

Hox genes are key regulators of appendage development in the insect body plan. The body plan of mayfly (Ephemeroptera) nymphs differs due to the presence of abdominal appendages called gills. Despite mayflies' phylogenetic position in Paleoptera and novel morphology amongst insects, little is known of their developmental genetics, such as the appendage-regulating Hox genes. To address this issue we present an annotated, early instar transcriptome and embryonic expression profiles for Antennapedia, Ultrabithorax, and Abdominal A proteins in the mayfly Hexagenia limbata, identify putative Hox protein sequences in the mayflies H. limbata, Cloeon dipterum, and Ephemera danica, and describe the genomic organization of the Hox gene cluster in E. danica.

RESULTS

Transcriptomic sequencing of early instar H. limbata nymphs yielded a high-quality assembly of 83,795 contigs, of which 22,975 were annotated against Folsomia candida, Nilaparvata lugens, Zootermopsis nevadensis and UniRef90 protein databases. Homeodomain protein phylogeny and peptide annotations identified coding sequences for eight of the ten canonical Hox genes (excluding zerknüllt/Hox3 and fushi tarazu) in H. limbata and C. dipterum, and all ten in E. danica. Mayfly Hox protein sequences and embryonic expression patterns of Antp, Ubx, and Abd-A appear highly conserved with those seen in other non-holometabolan insects. Similarly, the genomic organization of the Hox cluster in E. danica resembles that seen in most insects.

CONCLUSIONS

We present evidence that mayfly Hox peptide sequences and the embryonic expression patterns for Antp, Ubx, and Abd-A are extensively conserved with other insects, as is organization of the mayfly Hox gene cluster. The protein data suggest mayfly Antp, Ubx, and Abd-A play appendage promoting and repressing roles during embryogenesis in the thorax and abdomen, respectively, as in other insects. The identified expression of eight Hox genes, including Ubx and abd-A, in early instar nymphs further indicates a post-embryonic role, possibly in gill development. These data provide a basis for H. limbata as a complementary Ephemeridae model to the growing repertoire of mayfly model species and molecular techniques.

Multiscale analysis of the randomization limits of the chromosomal gene organization between Lepidoptera and Diptera

How chromosome gene organization and gene content evolve among distantly related and structurally malleable genomes remains unresolved. This is particularly the case when considering different insect orders. We have compared the highly contiguous genome assemblies of the lepidopteran Danaus plexippus and the dipteran Drosophila melanogaster, which shared a common ancestor around 290 Ma. The gene content of 23 out of 30 D. plexippus chromosomes was significantly associated with one or two of the six chromosomal elements of the Drosophila genome, denoting common ancestry. Despite the phylogenetic distance, 9.6% of the 1-to-1 orthologues still reside within the same ancestral genome neighbourhood. Furthermore, the comparison D. plexippus-Bombyx mori indicated that the rates of chromosome repatterning are lower in Lepidoptera than in Diptera, although still within the same order of magnitude. Concordantly, 14 developmental gene clusters showed a higher tendency to retain full or partial clustering in D. plexippus, further supporting that the physical association between the SuperHox and NK clusters existed in the ancestral bilaterian. Our results illuminate the scope and limits of the evolution of the gene organization and content of the ancestral chromosomes to the Lepidoptera and Diptera while helping reconstruct portions of the genome in their most recent common ancestor.

Complete mitochondrial genome sequencing of Oxycarenus laetus (Hemiptera: Lygaeidae) from two geographically distinct regions of India

Oxycarenus laetus is a seed-sap sucking pest affecting a variety of crops, including cotton plants. Rising incidence and pesticide resistance by O. laetus have been reported from India and neighbouring countries. In this study, O. laetus samples were collected from Bhatinda and Coimbatore (India). Pure mtDNA was isolated and sequenced using Illumina MiSeq. Both the samples were found to be identical species (99.9%), and the complete genome was circular (15,672 bp), consisting of 13 PCGs, 2 rRNA, 23 tRNA genes, and a 962 bp control region. The mitogenome is 74.1% AT-rich, 0.11 AT, and - 0.19 GC skewed. All the genes had ATN as the start codon except cox1 (TTG), and an additional trnT was predicted. Nearly all tRNAs folded into the clover-leaf structure, except trnS1 and trnV. The intergenic space between trnH and nad4, considered as a synapomorphy of Lygaeoidea, was displaced. Two 5 bp motifs AATGA and ACCTA, two tandem repeats, and a few microsatellite sequences, were also found. The phylogenetic tree was constructed using 36 mitogenomes from 7 super-families of Hemiptera by employing rigorous bootstrapping and ML. Ours is the first study to sequence the complete mitogenome of O. laetus or any Oxycarenus species. The findings from this study would further help in the evolutionary studies of Lygaeidae.

The evolution of insect biodiversity

Insects comprise over half of all described animal species. Together with the Protura (coneheads), Collembola (springtails) and Diplura (two-pronged bristletails), insects form the Hexapoda, a terrestrial arthropod lineage characterised by possessing six legs. Exponential growth of genome-scale data for the hexapods has substantially altered our understanding of the origin and evolution of insect biodiversity. Phylogenomics has provided a new framework for reconstructing insect evolutionary history, resolving their position among the arthropods and some long-standing internal controversies such as the placement of the termites, twisted-winged insects, lice and fleas. However, despite the greatly increased size of phylogenomic datasets, contentious relationships among key insect clades remain unresolved. Further advances in insect phylogeny cannot rely on increased depth and breadth of genome and taxon sequencing. Improved modelling of the substitution process is fundamental to countering tree-reconstruction artefacts, while gene content, modelling of duplications and deletions, and comparative morphology all provide complementary lines of evidence to test hypotheses emerging from the analysis of sequence data. Finally, the integration of molecular and morphological data is key to the incorporation of fossil species within insect phylogeny. The emerging integrated framework of insect evolution will help explain the origins of insect megadiversity in terms of the evolution of their body plan, species diversity and ecology. Future studies of insect phylogeny should build upon an experimental, hypothesis-driven approach where the robustness of hypotheses generated is tested against increasingly realistic evolutionary models as well as complementary sources of phylogenetic evidence.

Effects of Functional Depletion of Doublesex on Male Development in the Sawfly, Athalia rosae

Simple Summary

The sawfly, Athalia rosae, exploits a haplodiploid mode of reproduction, in which fertilized eggs develop into diploid females, whereas unfertilized eggs parthenogenetically develop into haploid males. The doublesex (dsx) gene is a well-conserved transcription factor that regulates sexual differentiation in insects. In the present study, we knocked down the A. rosae ortholog of dsx (Ardsx) during several developmental stages with repeated double-stranded RNA (dsRNA) injections. As a result, knockdown of Ardsx in haploid males caused almost complete male-to-female sex reversal, but the resulting eggs were infertile. The same knockdown approach using diploid males caused complete male-to-female sex reversal; they were able to produce fertile eggs and exhibited female behaviors. The same RNAi treatment did not affect female differentiation. These results demonstrated that dsx in the sawfly is essential for male development and its depletion caused complete male-to-female sex reversal. This is the first demonstration of functional depletion of dsx not causing intersexuality but inducing total sex reversal in males instead.

Abstract

The doublesex (dsx) gene, which encodes a transcription factor, regulates sexual differentiation in insects. Sex-specific splicing of dsx occurs to yield male- and female-specific isoforms, which promote male and female development, respectively. Thus, functional disruption of dsx leads to an intersexual phenotype in both sexes. We previously identified a dsx ortholog in the sawfly, Athalia rosae. Similar to dsx in other insects, dsx in the sawfly yields different isoforms in males and females as a result of alternative splicing. The sawfly exploits a haplodiploid mode of reproduction, in which fertilized eggs develop into diploid females, whereas unfertilized eggs parthenogenetically develop into haploid males. In the present study, we knocked down the A. rosae ortholog of dsx (Ardsx) during several developmental stages with repeated double-stranded RNA (dsRNA) injections. Knockdown of Ardsx via parental RNA interference (RNAi), which enables knockdown of genes in offspring embryos, led to a lack of internal and external genitalia in haploid male progeny. Additional injection of dsRNA targeting Ardsx in these animals caused almost complete male-to-female sex reversal, but the resulting eggs were infertile. Notably, the same knockdown approach using diploid males obtained by sib-crossing caused complete male-to-female sex reversal; they were morphologically and behaviorally females. The same RNAi treatment did not affect female differentiation. These results indicate that dsx in the sawfly is essential for male development and its depletion caused complete male-to-female sex reversal. This is the first demonstration of functional depletion of dsx not causing intersexuality but inducing total sex reversal in males instead.

Middle Jurassic origin in India: a new look at evolution of Vermileonidae and time-scaled relationships of lower brachyceran flies

Vermileonidae (wormlions) comprises 61 described species in 12 genera and is one of the oddest and most rarely collected dipteran families. Larvae of Vermileonidae are famous for their pitfall prey behaviour, whereas the adults are rarely seen. Here we report, for the first time, the complete mitochondrial genome from members of Vermileonidae. Phylogenetic reconstruction based on a representative sampling of the order reveals new insights into relationships between the Vermileonidae and other members of lower brachyceran flies. A sister-group relationship between Vermileonidae and Xylophagidae is supported, and the higher-level clade relationships are Tabanomorpha + (Muscomorpha + (Xylophagomorpha + Stratiomyomorpha)). Combining mitochondrial genome data with a morphological phylogeny, geographical distribution and geological history, we propose that the Vermileonidae originated in India during the Middle Jurassic, spreading to Africa via land bridges during the Late Jurassic, to the Palaearctic after the collision of India with Laurasia in the Late Oligocene, and to the Nearctic in the Early Miocene, via either North Atlantic or Beringian land bridges. Wet forested regions have proved to be barriers to their dispersal, explaining their absence from Central Africa, South America and eastern North America.

Evolution of Toll, Spatzle and MyD88 in insects: the problem of the Diptera bias

BACKGROUND

Arthropoda, the most numerous and diverse metazoan phylum, has species in many habitats where they encounter various microorganisms and, as a result, mechanisms for pathogen recognition and elimination have evolved. The Toll pathway, involved in the innate immune system, was first described as part of the developmental pathway for dorsal-ventral differentiation in Drosophila. Its later discovery in vertebrates suggested that this system was extremely conserved. However, there is variation in presence/absence, copy number and sequence divergence in various genes along the pathway. As most studies have only focused on Diptera, for a comprehensive and accurate homology-based approach it is important to understand gene function in a number of different species and, in a group as diverse as insects, the use of species belonging to different taxonomic groups is essential.

RESULTS

We evaluated the diversity of Toll pathway gene families in 39 Arthropod genomes, encompassing 13 different Insect Orders. Through computational methods, we shed some light into the evolution and functional annotation of protein families involved in the Toll pathway innate immune response. Our data indicates that: 1) intracellular proteins of the Toll pathway show mostly species-specific expansions; 2) the different Toll subfamilies seem to have distinct evolutionary backgrounds; 3) patterns of gene expansion observed in the Toll phylogenetic tree indicate that homology based methods of functional inference might not be accurate for some subfamilies; 4) Spatzle subfamilies are highly divergent and also pose a problem for homology based inference; 5) Spatzle subfamilies should not be analyzed together in the same phylogenetic framework; 6) network analyses seem to be a good first step in inferring functional groups in these cases. We specifically show that understanding Drosophila's Toll functions might not indicate the same function in other species.

CONCLUSIONS

Our results show the importance of using species representing the different orders to better understand insect gene content, origin and evolution. More specifically, in intracellular Toll pathway gene families the presence of orthologues has important implications for homology based functional inference. Also, the different evolutionary backgrounds of Toll gene subfamilies should be taken into consideration when functional studies are performed, especially for TOLL9, TOLL, TOLL2_7, and the new TOLL10 clade. The presence of Diptera specific clades or the ones lacking Diptera species show the importance of overcoming the Diptera bias when performing functional characterization of Toll pathways.

A review of the hexapod tracheal system with a focus on the apterygote groups

Respiratory systems are key innovations for the radiation of terrestrial arthropods. It is therefore surprising that there is still a considerable lack of knowledge. In this review of the available information on tracheal systems of hexapods (with a focus on the apterygote lineages Protura, Collembola, Diplura, Archaeognatha and Zygentoma), we summarize available data on the spiracles (number, position and morphology), the shape and variability of tracheal branching patterns including anastomoses, the tracheal fine structure and the respiratory proteins. The available data are strongly fragmented, and information for most subgroups is missing. In various cases, individual observations for one species account for the knowledge of the entire order. The available data show that there are strong differences between but also within apterygote orders. We conclude that the available data are insufficient to derive detailed conclusions on the hexapod ground plan and outline the possible evolutionary scenarios for the tracheal system in this group.

How are the mitochondrial genomes reorganized in Hexapoda? Differential evolution and the first report of convergences within Hexapoda

The evolution of the Hexapoda mitochondrial genome has been the focus of several genetic and evolutionary studies over the last decades. However, they have concentrated on certain taxonomic orders of economic or health importance. The recent increase of mitochondrial genomes sequencing of diverse taxonomic orders generates an important opportunity to clarify the evolution of this group of organisms. However, there is no comparative study that investigates the evolution of the Hexapoda mitochondrial genome. In order to verify the level of rearrangement and the mitochondrial genome evolution, we performed a comparative genomic analysis of the Hexapoda mitochondrial genome available in the NCBI database. Using a combination of bioinformatics methods to carefully examine the mitochondrial gene rearrangements in 1198 Hexapoda species belonging to 32 taxonomic orders, we determined that there is a great variation in the rate of rearrangement by gene and by taxonomic order. A higher rate of genetic reassortment is observed in Phthiraptera, Thysanoptera, Protura, and Hymenoptera; compared to other taxonomic orders. Twenty-four events of convergence in the genetic order between different taxonomic orders were determined, most of them not previously reported; which proves the great evolutionary dynamics within Hexapoda.

The complete mitochondrial genome of Trifida elongate and comparative analysis of 43 leafhoppers

In this study, the mitochondrial genome of Trifida elongate was sequenced, and comparative analysis of T. elongate and other 43 leafhoppers was performed based on the mitochondrial genome. The mitochondrial genome sequence length of T. elongate was 14,924 bp. It comprised 13 protein-coding genes, 2 ribosomal RNA (rRNA) genes, 22 transport RNA (tRNA) genes, and 1 non-coding control region. The control region is located between the rrnS and trnI genes, is characterized by two tandem repeats and three simple sequence repeats. Phylogenetic analysis showed that T. elongate is closely related with Bolanusoides shaanxiensis and Limassolla lingchuanensis (bootstrap value = 92% and posterior probabilities = 1). Analysis of synonymous and non-synonymous nucleotide substitutions showed that Ka/Ks value of the 13 protein-coding genes of 8 subfamily leafhoppers were less than 1 ranging from 0.0315 to 0.9928. atp8 had the highest Ka/Ks value whereas cox1 had the lowest Ka/Ks value. This study provides information on the structure and sequence characteristics of the mitochondrial genome of T. elongata. Typhlocybinae is clustered with (Cicadellinae+(Idiocerinae+(Mileewinae+(Nirvaninae+(Evacanthinae+Ledrinae))))).

Inferring Evolutionary Timescales without Independent Timing Information: An Assessment of "Universal" Insect Rates to Calibrate a Collembola (Hexapoda) Molecular Clock

Previous estimates of nucleotide substitution rates are routinely applied as secondary or “universal” molecular clock calibrations for estimating evolutionary timescales in groups that lack independent timing information. A major limitation of this approach is that rates can vary considerably among taxonomic groups, but the assumption of rate constancy is rarely evaluated prior to using secondary rate calibrations. Here I evaluate whether an insect mitochondrial DNA clock is appropriate for estimating timescales in Collembola—a group of insect-like arthropods characterized by high levels of cryptic diversity. Relative rates of substitution in cytochrome oxidase subunit 1 (COI) were inferred via Bayesian analysis across a topologically constrained Hexapod phylogeny using a relaxed molecular clock model. Rates for Collembola did not differ significantly from the average rate or from the rates estimated for most other groups (25 of 30), suggesting that (1) their apparent cryptic diversity cannot be explained by accelerated rates of molecular evolution and (2) clocks calibrated using “universal” insect rates may be appropriate for estimating evolutionary timescales in this group. However, of the 31 groups investigated, 10 had rates that deviated significantly from the average (6 higher, 4 lower), underscoring the need for caution and careful consideration when applying secondary insect rate calibrations. Lastly, this study exemplifies a relatively simple approach for evaluating rate constancy within a taxonomic group to determine whether the use of secondary rates are appropriate for molecular clock calibrations.

Mitochondrial genomes of stick insects (Phasmatodea) and phylogenetic considerations

Phasmatodea represents an order of hemimetabolous insects. This group includes species with extreme forms of masquerade crypsis, whereby they imitate twigs, bark, lichen, moss, and leaves. In this study, we sequenced and annotated three mitochondrial genomes (mitogenomes) from Phasmatodea. The lengths of the novel mitogenomes range from 14,162 bp to 15,879 bp. The gene content and organization correspond to those inferred for the ancestral insect. We conducted phylogenetic analyses together with the existing mitogenomes of polyneopterans and mayflies. In most cases, the Phasmatodea was non-monophyletic, with Embioptera and Zoraptera nested inside. The mitogenome sequences from Embioptera and Zoraptera suffered from high substitution rates and displayed very long branches in phylogenetic trees. The monophyletic Phasmatodea was recovered only when the analysis employed the site-heterogeneous CAT-GTR model in PhyloBayes and used the nucleotide dataset PCG_nt. The Euphasmatodea was well established by various data types and inference methods. In addition, the clade Heteropterygidae and the subfamilies Lonchodinae and Necrosciinae were strongly supported. The Australasian clade Lanceocercata was recovered across analyses. However, the Clitumninae was non-monophyletic.

Suppressing tawny crazy ant (Nylanderia fulva) by RNAi technology

The tawny crazy ant (Nylanderia fulva) is a new invasive pest in the United States. At present, its management mainly relies on the use of synthetic insecticides, which are generally ineffective at producing lasting control of the pest, necessitating alternative environmentally friendly measures. In this study, we evaluated the feasibility of gene silencing to control this ant species. Six housekeeping genes encoding actin (NfActin), coatomer subunit β (NfCOPβ), arginine kinase (NfArgK), and V-type proton ATPase subunits A (NfvATPaseA), B (NfvATPaseB) and E (NfvATPaseE) were cloned. Phylogenetic analysis revealed high sequence similarity to homologs from other ant species, particularly the Florida carpenter ant (Camponotus floridanus). To silence these genes, vector L4440 was used to generate six specific RNAi constructs for bacterial expression. Heat-inactivated, dsRNA-expressing Escherichia coli were incorporated into artificial diet. Worker ants exhibited reduced endogenous gene expression after feeding on such diet for 9 d. However, only ingestion of dsRNAs of NfCOPβ (a gene involved in protein trafficking) and NfArgK (a cellular energy reserve regulatory gene in invertebrates) caused modest but significantly higher ant mortality than the control. These results suggest that bacterially expressed dsRNA can be orally delivered to ant cells as a mean to target its vulnerabilities. Improved efficacy is necessary for the RNAi-based approach to be useful in tawny crazy ant management.

Putting the genome in insect phylogenomics

Next-generation sequencing technologies provide a substantial increase in the size of molecular phylogenetic datasets that can be obtained for studies of insect systematics. Several new genome reduction approaches are leveraging these technologies to generate large phylogenomic datasets: targeted amplicon sequencing, target capture, and transcriptome sequencing. Although cost effective, these approaches provide limited data for questions outside of phylogenetics. For many groups of insects, sequencing the entire genome at modest coverage is feasible. Using these genomic reads, an automated Target Restricted Assembly Method (aTRAM) can use the results of blast searches to assemble thousands of single copy ortholog genes across a group of interest. These locally assembled genes can then be compiled into very large phylogenomic datasets. These genomic libraries have the advantage in that they also contain reads from the mitochondrial genome and symbiont genomes, as well the entire insect genome, and can be leveraged for additional studies beyond phylogenetics.

The mitochondrial genomes of palaeopteran insects and insights into the early insect relationships

Phylogenetic relationships of basal insects remain a matter of discussion. In particular, the relationships among Ephemeroptera, Odonata and Neoptera are the focus of debate. In this study, we used a next-generation sequencing approach to reconstruct new mitochondrial genomes (mitogenomes) from 18 species of basal insects, including six representatives of Ephemeroptera and 11 of Odonata, plus one species belonging to Zygentoma. We then compared the structures of the newly sequenced mitogenomes. A tRNA gene cluster of IMQM was found in three ephemeropteran species, which may serve as a potential synapomorphy for the family Heptageniidae. Combined with published insect mitogenome sequences, we constructed a data matrix with all 37 mitochondrial genes of 85 taxa, which had a sampling concentrating on the palaeopteran lineages. Phylogenetic analyses were performed based on various data coding schemes, using maximum likelihood and Bayesian inferences under different models of sequence evolution. Our results generally recovered Zygentoma as a monophyletic group, which formed a sister group to Pterygota. This confirmed the relatively primitive position of Zygentoma to Ephemeroptera, Odonata and Neoptera. Analyses using site-heterogeneous CAT-GTR model strongly supported the Palaeoptera clade, with the monophyletic Ephemeroptera being sister to the monophyletic Odonata. In addition, a sister group relationship between Palaeoptera and Neoptera was supported by the current mitogenomic data.

Comparative transcriptomic analysis of a wing-dimorphic stonefly reveals candidate wing loss genes

Background

The genetic basis of wing development has been well characterised for model insect species, but remains poorly understood in phylogenetically divergent, non-model taxa. Wing-polymorphic insect species potentially provide ideal systems for unravelling the genetic basis of secondary wing reduction. Stoneflies (Plecoptera) represent an anciently derived insect assemblage for which the genetic basis of wing polymorphism remains unclear. We undertake quantitative RNA-seq of sympatric full-winged versus vestigial-winged nymphs of a widespread wing-dimorphic New Zealand stonefly, Zelandoperla fenestrata, to identify genes potentially involved in wing development and secondary wing loss.

Results

Our analysis reveals substantial differential expression of wing-development genes between full-winged versus vestigial-winged stonefly ecotypes. Specifically, of 23 clusters showing significant similarity to Drosophila wing development-related genes and their pea aphid orthologues, nine were significantly upregulated in full-winged stonefly ecotypes, whereas only one cluster (teashirt) was substantially upregulated in the vestigial-winged ecotype.

Conclusions

These findings suggest remarkable conservation of key wing-development pathways throughout 400 Ma of insect evolution. The finding that two Juvenile Hormone pathway clusters were significantly upregulated in vestigial-winged Zelandoperla supports the hypothesis that Juvenile Hormone may play a key role in modulating insect wing polymorphism, as has previously been suggested for other insect lineages.

CRISPR Disruption of BmOvo Resulted in the Failure of Emergence and Affected the Wing and Gonad Development in the Silkworm Bombyx mori

The domesticated silkworm is an economically important insect that is widely used as a lepidopteran insect model. Although somatic sex determination in the silkworm is well characterized, germline sex determination is not. Here, we used the transgenic-based CRISPR/Cas9 genome editing system to study the function of the Ovo gene in Bombyx mori. BmOvo is the homolog of a factor important in germline sex determination in Drosophila melanogaster. BmOvo mutants had abnormally shaped eggs that were disordered in the ovarioles, and gonad development was abnormal. Interestingly, wing discs and wings did not develop properly, and most of the mutants failed to eclose. Gene expression analyses by qRT-PCR showed that BmOvo gene was highly expressed in the wing disc and epidermis. Genes involved in the WNT signaling pathway and wing development genes BmWCP10 and BmE74 were downregulated in the BmOvo mutants when compared with wild-type animals. These results demonstrate that the BmOvo gene product plays an important role in wing metamorphosis. Thus, this study provides new insights into the multiple functions of BmOvo beyond germline sex determination.

Mutualistic interactions reshuffle the effects of climate change on plants across the tree of life

Climatically induced local species extinctions may trigger coextinction cascades, thus driving many more species to extinction than originally predicted by species distribution models. Using seven pollination networks across Europe that include the phylogeny and life history traits of plants, we show a substantial variability across networks in climatically predicted plant extinction-and particularly the subsequent coextinction-rates, with much higher values in Mediterranean than Eurosiberian networks. While geographic location best predicts the probability of a plant species to be driven to extinction by climate change, subsequent coextinctions are best predicted by the local network of interactions. These coextinctions not only increase the total number of plant species being driven to extinction but also add a bias in the way the major taxonomic and functional groups are pruned.

The phylogeny and evolutionary timescale of stoneflies (Insecta: Plecoptera) inferred from mitochondrial genomes

Phylogenetic analysis based on mitochondrial genomic data from 25 stonefly species recovered a well-supported tree resolving higher-level relationships within Plecoptera (stoneflies). The monophyly of both currently recognized suborders was strongly supported, concordant with previous molecular analyses of Plecoptera. The southern hemisphere suborder Antarctoperlaria formed two clades: Eustheniidae + Diamphipnoidae and Austroperlidae + Gripopterygidae; consistent with relationships proposed based on morphology. The largely northern hemisphere suborder Arctoperlaria also divided into two groups, Euholognatha and Systellognatha, each composed of the five families traditionally assigned to each infraorder (the placement Scopuridae by mt genome data remains untested at this time). Within Euholognatha, strong support for the clade Nemouridae + Notonemouridae confirmed the northern origin of the currently southern hemisphere restricted Notonemouridae. Other family level relationships within the Arctoperlaria differ from those recovered by previous morphology and molecular based analyses. A fossil-calibrated divergence estimation suggests the formation of two suborders dates back to the Jurassic (181 Ma), with subsequent diversification of most stonefly families during the Cretaceous. This result confirms the hypothesis that initial divergence between the suborders was driven by the breakup of the supercontinent Pangaea into Laurasia and Gondwanaland (commencing 200 Ma and complete by 150 Ma).

Phthiraptera of Canada

There are approximately 463 species of parasitic lice recorded in Canada, in three suborders: Amblycera, six families; Ischnocera, two families; Anoplura, eight families. At least an additional 361 species may eventually be recorded based on presence of suitable hosts and proximity to known distributions. Approximately 41 species are introduced non-native species. Only about 54% of the expected chewing louse fauna has been recorded, and considerable collecting effort is needed, especially for lice infesting passerine birds, shorebirds, and seabirds. The sucking louse fauna is well known, with approximately 88% of the expected fauna recorded. Investigations into ecology of lice and the nature of relationships with their hosts are badly needed. Barcode Index Numbers are available for only 13 species of parasitic lice in Canada.

Multispecies coalescent analysis confirms standing phylogenetic instability in Hexapoda

The multispecies coalescent (MSC) has been increasingly used in phylogenomic analyses due to the accommodation of gene tree topological heterogeneity by taking into account population-level processes, such as incomplete lineage sorting. In this sense, the phylogeny of insect species, which are characterized by their large effective population sizes, is suitable for a coalescent-based analysis. Furthermore, studies so far recovered short internal branches at early divergences of the insect tree of life, indicating fast evolutionary radiations that increase the probability of incomplete lineage sorting in deep time. Here, we investigated the performance of the MSC for a phylogenomic data set of hexapods compiled by Misof et al. (2014, Science 346:763). Our analysis recovered the monophyly of most insect orders, and major phylogenetic relationships were in agreement with current insect systematics. We identified, however, some evolutionary associations that were consistently problematic. Most noticeable, Hexapod monophyly was disrupted by the sister group relationship between the remiped crustacean and Insecta. Additionally, the interordinal relationships within Polyneoptera and Neuropteroidea were found to be phylogenetically unstable. We show that these controversial phylogenetic arrangements were also poorly supported by previous analyses, and therefore, we evaluated their robustness to stochastic errors from sampling sites and terminals, confirming standing problems in hexapod phylogeny in the genomics age.

Characterizing the embryonic development of B. hygida (Diptera: Sciaridae) following enzymatic treatment to permeabilize the serosal cuticle

Understanding the evolution of the developmental programs active during dipteran embryogenesis depends on comparative studies. As a counterpoint to the intensively investigated and highly derived cyclorrhaphan flies that include the model organism Drosophila melanogaster, we are studying the basal Diptera Bradysia hygida, a member of the Sciaridae family that is amenable to laboratory cultivation. Here we describe the B. hygida embryogenesis, which lasts 9 days at 22 °C. The use of standard fixation D. melanogaster protocols resulted in embryos refractory to DAPI staining and to overcome this, a new enzyme-based method was developed. Calcofluor-White staining of enzimatically-treated embryos revealed that this method removes chitin from the serosal cuticle surrounding the B. hygida embryo. Chitin is one of the main components of serosal cuticles and searches in a B. hygida embryonic transcriptome database revealed conservation of the chitin synthesis pathway, further supporting the occurrence of chitin biosynthesis in B. hygida embryos. Combining the enzymatic treatment protocol with the use of both DIC and fluorescence microscopy allowed the first complete description of the B. hygida embryogenesis. Our results constitute an important step towards the understanding of early development of a basal Diptera and pave the way for future evo-devo studies.

Fragmentary Gene Sequences Negatively Impact Gene Tree and Species Tree Reconstruction

Species tree reconstruction from genome-wide data is increasingly being attempted, in most cases using a two-step approach of first estimating individual gene trees and then summarizing them to obtain a species tree. The accuracy of this approach, which promises to account for gene tree discordance, depends on the quality of the inferred gene trees. At the same time, phylogenomic and phylotranscriptomic analyses typically use involved bioinformatics pipelines for data preparation. Errors and shortcomings resulting from these preprocessing steps may impact the species tree analyses at the other end of the pipeline. In this article, we first show that the presence of fragmentary data for some species in a gene alignment, as often seen on real data, can result in substantial deterioration of gene trees, and as a result, the species tree. We then investigate a simple filtering strategy where individual fragmentary sequences are removed from individual genes but the rest of the gene is retained. Both in simulations and by reanalyzing a large insect phylotranscriptomic data set, we show the effectiveness of this simple filtering strategy.

Topological support and data quality can only be assessed through multiple tests in reviewing Blattodea phylogeny

Assessing support for molecular phylogenies is difficult because the data is heterogeneous in quality and overwhelming in quantity. Traditionally, node support values (bootstrap frequency, Bayesian posterior probability) are used to assess confidence in tree topologies. Other analyses to assess the quality of phylogenetic data (e.g. Lento plots, saturation plots, trait consistency) and the resulting phylogenetic trees (e.g. internode certainty, parameter permutation tests, topological tests) exist but are rarely applied. Here we argue that a single qualitative analysis is insufficient to assess support of a phylogenetic hypothesis and relate data quality to tree quality. We use six molecular markers to infer the phylogeny of Blattodea and apply various tests to assess relationship support, locus quality, and the relationship between the two. We use internode-certainty calculations in conjunction with bootstrap scores, alignment permutations, and an approximately unbiased (AU) test to assess if the molecular data unambiguously support the phylogenetic relationships found. Our results show higher support for the position of Lamproblattidae, high support for the termite phylogeny, and low support for the position of Anaplectidae, Corydioidea and phylogeny of Blaberoidea. We use Lento plots in conjunction with mutation-saturation plots, calculations of locus homoplasy to assess locus quality, identify long branch attraction, and decide if the tree's relationships are the result of data biases. We conclude that multiple tests and metrics need to be taken into account to assess tree support and data robustness.

The first two mitochondrial genomes of the family Aphelinidae with novel gene orders and phylogenetic implications

Chalcidoidea is one of the most diverse group in Hymenoptera by possessing striking mitochondrial gene arrangement. By using next generation sequencing method, the first two nearly complete mitochondrial genomes in the family Aphelinidae (Insecta, Hymenopetra, Chalcidoidea) were obtained in this study. Almost all previously sequenced mitochondrial genome of Chalcidoidea species have a large inversion including six genes (atp6-atp8-trnD-trnK-cox2-trnL2-cox1) as compared with ancestral mitochondrial genome, but these two Encarsia mitochondrial genomes had a large inversion including nine genes (nad3-trnG-atp6-atp8-trnD-trnK-cox2-trnL2-cox1), which was only congruent with the species in the genus Nasonia. Moreover, we found that one shuffling changes (trnD and trnK) happened in the species E. obtusiclava but not in another species E. formosa within the same genus, of which such shuffling within the same genus at this region was only detected in Polisters within Insecta. Phylogenetic analysis displayed that different data matrix (13PCG+ 2 rRNA or 13 PCG) and inference methods (BI or ML) indicate the identical topology with high nodal supports that Aphelinidae formed a sister group with (Trichogrammatidae + Aganoidae) and the monophyly of Pteramalidae. Our results also indicated the validity of assembling and feasibility of next-generation technology to obtain the mitochondrial genomes of parasitic Hymenoptera.

DNA barcoding of microgastrine parasitoid wasps (Hymenoptera: Braconidae) using high-throughput methods more than doubles the number of species known for Australia

The Microgastrinae are a hugely diverse subfamily of endoparasitoid wasps of lepidopteran caterpillars. They are important in agriculture as biological control agents and play a significant ecological role in the regulation of caterpillar populations. Whilst the group has been the focus of intensive rearing and DNA barcoding studies in the Northern Hemisphere, the Australian fauna has received little attention. In total, 99 species have been described from or have been introduced into Australia, but the real species diversity for the region is clearly much larger than this. In this study, museum ethanol samples and recent field collections were mined for hundreds of specimens of microgastrine wasps, which were then barcoded for the COI region, ITS2 ribosomal spacer and the wingless nuclear genes, using a pooled sequencing approach on an Illumina Miseq system. Full COI sequences were obtained for 525 individuals which, when combined with 162 publicly available sequences, represented 417 haplotypes, and a total of 236 species were delimited using a consensus approach. By more than doubling the number of known microgastrine wasp species in Australia, our study highlights the value of DNA barcoding in the context of employing high-throughput sequencing methods of bulk ethanol museum collections for biodiversity assessment.

Construction of a Species-Level Tree of Life for the Insects and Utility in Taxonomic Profiling

Although comprehensive phylogenies have proven an invaluable tool in ecology and evolution, their construction is made increasingly challenging both by the scale and structure of publically available sequences. The distinct partition between gene-rich (genomic) and species-rich (DNA barcode) data is a feature of data that has been largely overlooked, yet presents a key obstacle to scaling supermatrix analysis. I present a phyloinformatics framework for draft construction of a species-level phylogeny of insects (Class Insecta). Matrix-building requires separately optimized pipelines for nuclear transcriptomic, mitochondrial genomic, and species-rich markers, whereas tree-building requires hierarchical inference in order to capture species-breadth while retaining deep-level resolution. The phylogeny of insects contains 49,358 species, 13,865 genera, 760 families. Deep-level splits largely reflected previous findings for sections of the tree that are data rich or unambiguous, such as inter-ordinal Endopterygota and Dictyoptera, the recently evolved and relatively homogeneous Lepidoptera, Hymenoptera, Brachycera (Diptera), and Cucujiformia (Coleoptera). However, analysis of bias, matrix construction and gene-tree variation suggests confidence in some relationships (such as in Polyneoptera) is less than has been indicated by the matrix bootstrap method. To assess the utility of the insect tree as a tool in query profiling several tree-based taxonomic assignment methods are compared. Using test data sets with existing taxonomic annotations, a tendency is observed for greater accuracy of species-level assignments where using a fixed comprehensive tree of life in contrast to methods generating smaller de novo reference trees. Described herein is a solution to the discrepancy in the way data are fit into supermatrices. The resulting tree facilitates wider studies of insect diversification and application of advanced descriptions of diversity in community studies, among other presumed applications. [Data integration; data mining; insects; phylogenomics; phyloinformatics; tree of life.].

Entomological Collections in the Age of Big Data

With a million described species and more than half a billion preserved specimens, the large scale of insect collections is unequaled by those of any other group. Advances in genomics, collection digitization, and imaging have begun to more fully harness the power that such large data stores can provide. These new approaches and technologies have transformed how entomological collections are managed and utilized. While genomic research has fundamentally changed the way many specimens are collected and curated, advances in technology have shown promise for extracting sequence data from the vast holdings already in museums. Efforts to mainstream specimen digitization have taken root and have accelerated traditional taxonomic studies as well as distribution modeling and global change research. Emerging imaging technologies such as microcomputed tomography and confocal laser scanning microscopy are changing how morphology can be investigated. This review provides an overview of how the realization of big data has transformed our field and what may lie in store.

Egg structure and outline of embryonic development of the basal mantodean, Metallyticus splendidus Westwood, 1835 (Insecta, Mantodea, Metallyticidae)

The egg structure and outline of the embryonic development of Metallyticus splendidus of one of the basal Mantodea representatives, Metallyticidae, were described in the present study. The results obtained were compared with those from the previous studies, to reconstruct and discuss the groundplan of Mantodea and Dictyoptera. In M. splendidus, the egg is spheroidal, it has a convex ventral side at the center in which numerous micropyles are grouped, and it possesses a conspicuous hatching line in its anterior half. These are the groundplan features of mantodean eggs and the "grouped micropyles in the ventral side of the egg" are regarded as an apomorphic groundplan feature of Dictyoptera. A small circular embryo is formed by a simple concentration of blastoderm cells, which then undergoes embryogenesis of the typical short germ band type. Blastokinesis is of the "non-reversion type" and the embryo keeps its original superficial position and original orientation throughout embryonic development. During the middle stages of development, the embryo undergoes rotation around the egg's anteroposterior axis. These features are a part of the groundplan of Mantodea. It is uncertain whether sharing of the "non-reversion type" of blastokinesis by Mantodea and blaberoidean Blattodea can be regarded as homology or homoplasy.

Costs of immune responses are related to host body size and lifespan

A central assumption in ecological immunology is that immune responses are costly, with costs manifesting directly (e.g., increases in metabolic rate and increased amino acid usage) or as tradeoffs with other life processes (e.g., reduced growth and reproductive success). Across taxa, host longevity, timing of maturity, and reproductive effort affect the organization of immune systems. It is reasonable, therefore, to expect that these and related factors should also affect immune activation costs. Specifically, species that spread their breeding efforts over a long lifetime should experience lower immune costs than those that mature and breed quickly and die comparatively early. Likewise, body mass should affect immune costs, as body size affects the extent to which hosts are exposed to parasites as well as how hosts can combat infections (via its effects on metabolic rates and other factors). Here, we used phylogenetic meta-regression to reveal that, in general, animals incur costs of immune activation, but small species that are relatively long-lived incur the largest costs. These patterns probably arise because of the relative need for defense when infection risk is comparatively high and fitness can only be realized over a comparatively long period. However, given the diversity of species considered here and the overall modest effects of body mass and life history on immune costs, much more research is necessary before generalizations are appropriate.

Embryonic development of Eucorydia yasumatsui Asahina, with special reference to external morphology (Insecta: Blattodea, Corydiidae)

As the first step in the comparative embryological study of Blattodea, with the aim of reconstructing the groundplan and phylogeny of Dictyoptera and Polyneoptera, the embryonic development of a corydiid was examined and described in detail using Eucorydia yasumatsui. Ten to fifteen micropyles are localized on the ventral side of the egg, and aggregated symbiont bacterial "mycetomes" are found in the egg. The embryo is formed by the fusion of paired blastodermal regions, with higher cellular density on the ventral side of the egg. This type of embryo formation, regarded as one of the embryological autapomorphies of Polyneoptera, was first demonstrated for "Blattaria" in the present study. The embryo undergoes embryogenesis of the short germ band type, and elongates to its full length on the ventral side of the egg. The embryo undergoes katatrepsis and dorsal closure, and then finally, it acquires its definitive form, keeping its original position on the ventral side of the egg, with its anteroposterior axis never reversed throughout development. The information obtained was compared with that of previous studies on other insects. "Micropyles grouped on the ventral side of the egg" is thought to be a part of the groundplan of Dictyoptera, and "possession of bacteria in the form of mycetomes" to be an apomorphic groundplan of Blattodea. Corydiid embryos were revealed to perform blastokinesis of the "non-reversion type (N)", as reported in blaberoid cockroaches other than Corydiidae ("Ectobiidae," Blaberidae, etc.) and in Mantodea; the embryos of blattoid cockroaches (Blattidae and Cryptocercidae) and Isoptera undergo blastokinesis of the "reversion type (R)," in which the anteroposterior axis of the embryo is reversed during blastokinesis. Dictyopteran blastokinesis types can be summarized as "Mantodea (N) + Blattodea [= Blaberoidea (N) + Blattoidea (R) + Isoptera (R)]".

Hypoxia Inducible Factor (HIF) transcription factor family expansion, diversification, divergence and selection in eukaryotes

Hypoxia inducible factor (HIF) transcription factors are crucial for regulating a variety of cellular activities in response to oxygen stress (hypoxia). In this study, we determine the evolutionary history of HIF genes and their associated transactivation domains, as well as perform selection and functional divergence analyses across their four characteristic domains. Here we show that the HIF genes are restricted to metazoans: At least one HIF-α homolog is found within the genomes of non-bilaterians and bilaterian invertebrates, while most vertebrate genomes contain between two and six HIF-α genes. We also find widespread purifying selection across all four characteristic domain types, bHLH, PAS, NTAD, CTAD, in HIF-α genes, and evidence for Type I functional divergence between HIF-1α, HIF-2α /EPAS, and invertebrate HIF genes. Overall, we describe the evolutionary histories of the HIF transcription factor gene family and its associated transactivation domains in eukaryotes. We show that the NTAD and CTAD domains appear de novo, without any appearance outside of the HIF-α subunits. Although they both appear in invertebrates as well as vertebrate HIF- α sequences, there seems to have been a substantial loss across invertebrates or were convergently acquired in these few lineages. We reaffirm that HIF-1α is phylogenetically conserved among most metazoans, whereas HIF-2α appeared later. Overall, our findings can be attributed to the substantial integration of this transcription factor family into the critical tasks associated with maintenance of oxygen homeostasis and vascularization, particularly in the vertebrate lineage.

Identification and functional characterization of the sex-determining gene doublesex in the sawfly, Athalia rosae (Hymenoptera: Tenthredinidae)

Sexual fate of the sawfly, Athalia rosae (Hymenoptera: Tenthredinidae) is determined by the complementary sex determination (CSD) mechanism as is the case in honeybees. However, to date, genes involved in sex determination have not been identified in this species. In this study, we attempted to identify orthologs of complementary sex-determiner (csd), feminizer (fem), and doublesex (dsx) from the A. rosae genome, all of which are crucial components of the sex determination cascade in the honeybee. As a result, we identified a sawfly ortholog of dsx (designated as Ardsx). Rapid amplification of cDNA ends (RACE) using total RNA extracted from male and female larvae identified three male-specific variants and three female-specific variants. Comparison between the full-length Ardsx cDNAs and the genomic sequence revealed that exon 5 was differentially spliced between the male- and female-specific variants. RT-PCR analysis demonstrated that Ardsx pre-mRNA was spliced alternatively in a sex-dependent manner at almost all the developmental stages. RNAi-mediated knockdown of Ardsx in males caused severe defects in the reproductive organs and, notably, induced development of the ovipository apparatus containing the dorsal pair of blades and the sheath. These males also showed abnormalities in testes and seminal vesicles and lacked mature sperm. The present study provides the first direct evidence that dsx is essential for sexual development in hymenopteran species.

Predictors of elevational biodiversity gradients change from single taxa to the multi-taxa community level

The factors determining gradients of biodiversity are a fundamental yet unresolved topic in ecology. While diversity gradients have been analysed for numerous single taxa, progress towards general explanatory models has been hampered by limitations in the phylogenetic coverage of past studies. By parallel sampling of 25 major plant and animal taxa along a 3.7 km elevational gradient on Mt. Kilimanjaro, we quantify cross-taxon consensus in diversity gradients and evaluate predictors of diversity from single taxa to a multi-taxa community level. While single taxa show complex distribution patterns and respond to different environmental factors, scaling up diversity to the community level leads to an unambiguous support for temperature as the main predictor of species richness in both plants and animals. Our findings illuminate the influence of taxonomic coverage for models of diversity gradients and point to the importance of temperature for diversification and species coexistence in plant and animal communities.

Explaining species richness patterns is a key question in ecology. Peters et al. sample diverse plant and animal groups across elevation on Mt. Kilimanjaro to show that, while disparate factors drive distributions of individual taxa, diversity overall decreases with elevation, mostly driven by effects of temperature.