InsectBase: a resource for insect genomes and transcriptomes

Knockdown of the aminopeptidase N genes decreases susceptibility of Chilo suppressalis larvae to Cry1Ab/Cry1Ac and Cry1Ca

Bacillus thuringiensis (Bt) insecticide is currently the most widely used bioinsecticide. Bt expressing cry genes are some of the most successful foreign-genome-inserting genes used in transgenic insect-resistant crop development. Cry toxins are resistant to lepidopteran pests, such as Chilo suppressalis, a major insect pest of rice worldwide. Since Cry toxins exert their activity by binding to specific receptors in the midgut of target insects, identification of functional Cry toxin receptors in the midgut of C. suppressalis larvae is crucial to evaluate potential resistance mechanisms and develop effective strategies for inhibiting insect resistance. In this study, we isolated two aminopeptidase N genes (APN6 and APN8) from C. suppressalis and determined that they were expressed in the foregut. APN6 was highly expressed at the fourth instar, and APN8 was highly expressed in adult and pupa. Knockdown of CsAPN6 and CsAPN8 by RNA interference resulted in significantly decreased susceptibility of larvae to Bt rice varieties TT51 (expressing cry1Ac/cry1Ab fusion genes) and T1C-19 (expressing cry1Ca), but not T2A-1 (expressing cry2Aa). These findings suggest that both APN6 and APN8 are involved in the toxicity of Cry1Ac/Cry1Ab and Cry1Ca toxins.

Insect genomes: progress and challenges

In the wake of constant improvements in sequencing technologies, numerous insect genomes have been sequenced. Currently, 1219 insect genome-sequencing projects have been registered with the National Center for Biotechnology Information, including 401 that have genome assemblies and 155 with an official gene set of annotated protein-coding genes. Comparative genomics analysis showed that the expansion or contraction of gene families was associated with well-studied physiological traits such as immune system, metabolic detoxification, parasitism and polyphagy in insects. Here, we summarize the progress of insect genome sequencing, with an emphasis on how this impacts research on pest control. We begin with a brief introduction to the basic concepts of genome assembly, annotation and metrics for evaluating the quality of draft assemblies. We then provide an overview of genome information for numerous insect species, highlighting examples from prominent model organisms, agricultural pests and disease vectors. We also introduce the major insect genome databases. The increasing availability of insect genomic resources is beneficial for developing alternative pest control methods. However, many opportunities remain for developing data-mining tools that make maximal use of the available insect genome resources. Although rapid progress has been achieved, many challenges remain in the field of insect genomics.

The genomes of two parasitic wasps that parasitize the diamondback moth

BACKGROUND

Parasitic insects are well-known biological control agents for arthropod pests worldwide. They are capable of regulating their host's physiology, development and behaviour. However, many of the molecular mechanisms involved in host-parasitoid interaction remain unknown.

RESULTS

We sequenced the genomes of two parasitic wasps (Cotesia vestalis, and Diadromus collaris) that parasitize the diamondback moth Plutella xylostella using Illumina and Pacbio sequencing platforms. Genome assembly using SOAPdenovo produced a 178 Mb draft genome for C. vestalis and a 399 Mb draft genome for D. collaris. A total set that contained 11,278 and 15,328 protein-coding genes for C. vestalis and D. collaris, respectively, were predicted using evidence (homology-based and transcriptome-based) and de novo prediction methodology. Phylogenetic analysis showed that the braconid C. vestalis and the ichneumonid D. collaris diverged approximately 124 million years ago. These two wasps exhibit gene gains and losses that in some cases reflect their shared life history as parasitic wasps and in other cases are unique to particular species. Gene families with functions in development, nutrient acquisition from hosts, and metabolism have expanded in each wasp species, while genes required for biosynthesis of some amino acids and steroids have been lost, since these nutrients can be directly obtained from the host. Both wasp species encode a relative higher number of neprilysins (NEPs) thus far reported in arthropod genomes while several genes encoding immune-related proteins and detoxification enzymes were lost in both wasp genomes.

CONCLUSIONS

We present the annotated genome sequence of two parasitic wasps C. vestalis and D. collaris, which parasitize a common host, the diamondback moth, P. xylostella. These data will provide a fundamental source for studying the mechanism of host control and will be used in parasitoid comparative genomics to study the origin and diversification of the parasitic lifestyle.

A Novel Software and Method for the Efficient Development of Polymorphic SSR Loci Based on Transcriptome Data

Traditional methods for developing polymorphic microsatellite loci without reference sequences are time-consuming and labor-intensive, and the polymorphisms of simple sequence repeat (SSR) loci developed from expressed sequence tag (EST) databases are generally poor. To address this issue, in this study, we developed a new software (PSSRdt) and established an effective method for directly obtaining polymorphism details of SSR loci by analyzing diverse transcriptome data. The new method includes three steps, raw data processing, PSSRdt application, and loci extraction and verification. To test the practicality of the method, we successfully obtained 1940 potential polymorphic SSRs from the transcript dataset combined with 44 pea aphid transcriptomes. Fifty-two SSR loci obtained by the new method were selected for validating the polymorphic characteristics by genotyping in pea aphid individuals. The results showed that over 92% of SSR loci were polymorphic and 73.1% of loci were highly polymorphic. Our new software and method provide an innovative approach to microsatellite development based on RNA-seq data, and open a new path for the rapid mining of numerous loci with polymorphism to add to the body of research on microsatellites.

Birth-and-Death Evolution of the Fatty Acyl-CoA Reductase (FAR) Gene Family and Diversification of Cuticular Hydrocarbon Synthesis in Drosophila

The birth-and-death evolutionary model proposes that some members of a multigene family are phylogenetically stable and persist as a single copy over time, whereas other members are phylogenetically unstable and undergo frequent duplication and loss. Functional studies suggest that stable genes are likely to encode essential functions, whereas rapidly evolving genes reflect phenotypic differences in traits that diverge rapidly among species. One such class of rapidly diverging traits are insect cuticular hydrocarbons (CHCs), which play dual roles in chemical communications as short-range recognition pheromones as well as protecting the insect from desiccation. Insect CHCs diverge rapidly between related species leading to ecological adaptation and/or reproductive isolation. Because the CHC and essential fatty acid biosynthetic pathways share common genes, we hypothesized that genes involved in the synthesis of CHCs would be evolutionary unstable, whereas those involved in fatty acid-associated essential functions would be evolutionary stable. To test this hypothesis, we investigated the evolutionary history of the fatty acyl-CoA reductases (FARs) gene family that encodes enzymes in CHC synthesis. We compiled a unique data set of 200 FAR proteins across 12 Drosophila species. We uncovered a broad diversity in FAR content which is generated by gene duplications, subsequent gene losses, and alternative splicing. We also show that FARs expressed in oenocytes and presumably involved in CHC synthesis are more unstable than FARs from other tissues. Taken together, our study provides empirical evidence that a comparative approach investigating the birth-and-death evolution of gene families can identify candidate genes involved in rapidly diverging traits between species.

Identification of candidate ATP-binding cassette transporter gene family members in Diaphorina citri (Hemiptera: Psyllidae) via adult tissues transcriptome analysis

The ATP-binding cassette (ABC) transporters exist in all living organisms and play major roles in various biological functions by transporting a wide variety of substrates across membranes. The functions of ABC transporters in drug resistance have been extensively studied in vertebrates; however, they are rarely characterized in agricultural pests. The Asian citrus psyllid, Diaphorina citri, is one of the most damaging pests of the Citrus genus because of its transmission of Huanglongbing, also known as Yellow Dragon disease. In this study, the next-generation sequencing technique was applied to research the ABC transporters of D. citri. Fifty-three ABC transporter genes were found in the RNA-Seq data, and among these ABC transporters, 4, 4, 5, 2, 1, 4, 18 and 15 ABC proteins belonged to the ABCA-ABCH subfamilies, respectively. Different expression profiles of 52 genes between imidacloprid-resistant and imidacloprid-susceptible strains were studied by qRT-PCR; 5 ABCGs and 4 ABCHs were significantly upregulated in the imidacloprid-resistant strain. In addition, five of the nine upregulated genes were widely expressed in adult tissues in spatial expression analysis. The results suggest that these genes may play key roles in this phenotype. In general, this study contributed to our current understanding of D. citri resistance to insecticides.

The functional difference of eight chitinase genes between male and female of the cotton mealybug, Phenacoccus solenopsis

The cotton mealybug Phenacoccus solenopsis Tinsley (Hemiptera: Pseudococcidae) is a polyphagous insect that attacks tens of plant and causes substantial economic loss. Insect chitinases are required to remove the old cuticle to allow for continued growth and development. Though insect chitinases have been well studied in tens of insects, their functions in mealybug are still not addressed. Here, we sequenced the transcriptomes of adult males and females, from which eight chitinase genes were identified. We then used the method of rapid amplification of cDNA ends to amplify their full length. Phylogenetic analysis indicated that these genes clustered into five subgroups. Among which, group II PsCht2 had the longest transcript and was highly expressed at second instar nymph. PsCht10, PsCht3-3 and PsIDGF were highly expressed in the adult females, whereas PsCht4 and PsCht4-1 were significantly expressed at the male pupa and adult male. Next, we knocked down all eight chitinase genes by feeding the double-stranded RNA. Knockdown of PsCht4 or PsCht4-1 led to the failure of moult and, silencing PsCht5 resulted in pupation defect, while silencing PsCht10 led to small body size, suggesting these genes have essential roles in development and can be used as a potential target for pest control.

microRNA-14 as an efficient suppressor to switch off ecdysone production after ecdysis in insects

The precise increase and decrease of hormone ecdysone are critical for accurate development in insects. Most previous works focus on transcriptional activation of ecdysone production; however, little is known about the mechanism of switching off ecdysone biosynthesis after ecdysis. Here, we showed that the precursor microRNA-14 (pre-miR-14) encodes two mature miRNAs in silkworm; both of these two mature miRNAs regulate various genes in the ecdysone-signalling pathway. Bmo-miR-14-5p targets on nine genes whereas Bmo-miR-14-3p targets on two genes in the same pathway. These two mature miRNAs increased immediately after the ecdysis, efficiently suppressing the 20-hydroxyecdysone (20E) biosynthesis, the upstream regulation, and the downstream response genes. Knocking down either of two mature miRNAs or both of them delays moult development, impairing development synchrony in antagomir-treated groups. In addition, overexpressing Bmo-miR-14-5p but not Bmo-miR-14-3p significantly affected the 20E titer and increased the moulting time variation, suggesting that Bmo-miR-14-5p, though it is less abundant, has more potent effects in development regulation than Bmo-miR-14-3p. In summary, we present evidence that a pre-miRNA encodes two mature miRNAs targeting on the same pathway, which significantly improves miRNA regulation efficiencies to programmatically switch off ecdysone biosynthesis.

Robust reference gene design and validation for expression studies in the large milkweed bug, Oncopeltus fasciatus, upon cardiac glycoside stress

Despite its history as a developmental and evolutionary model organism, gene expression analysis in the large milkweed bug, Oncopeltus fasciatus, has rarely been explored using quantitative real-time PCR. The strength of this method depends greatly on the endogenous controls used for normalization, which are lacking for the milkweed bug system. Here, to fill in this gap in our knowledge, we validated the stability of a set of ten candidate reference genes identified from the O. fasciatus transcriptome, and did so upon exposure to a dietary toxin, a cardiac glycoside, and across four different exposure periods. To increase robustness against gDNA contaminants, genome resources were used to design intron-bridging primers. A comprehensive stability validation by the Bestkeeper, Normfinder, geNorm and comparative ΔCt methods identified ef1a and tubulin as the most stable genes across treatments and time points, whereas 18S rRNA was the most unstable. However, accounting for the temporal scale indicated that time point confined normalizers might enable higher quantification accuracy for treatment comparison. Overall this study demonstrates: (i) a robust RT-qPCR primer design approach is possible for non-model organisms where genome annotation is often incomplete, and (ii) the importance of detailed reference gene stability exploration in multifactorial experimental designs.

Species-Specific Discrimination of Insect Meals for Aquafeeds by Direct Comparison of Tandem Mass Spectra

Insect protein has the potential to become a sustainable feed ingredient for the rapidly growing aquaculture industry. In the European Union, insect derived protein is placed under the same legislation as processed animal proteins (PAP). It is therefore of interest to develop methods for regulatory use, which unambiguously identify the species origin of insect-based ingredients. We performed (i) total protein quantification of insect samples using the traditional nitrogen-to-protein conversion factor of 6.25 and the sum of anhydrous amino acids, (ii) quantitative amino acid profiling and (iii) high-throughput tandem mass spectrometry to describe and differentiate 18 different commercial-grade insect meal samples derived from Hermetia illucens (8), Tenebrio molitor (5), Alphitobius diaperinus (3) and Acheta domesticus (2). In addition, we investigated and compared different protein extraction and digestion protocols for proteomic analysis. We found that irrespective of sample preparation, shotgun proteomics in combination with direct spectral comparison were able to differentiate insect meal according to their taxonomic classification. The insect specific spectral libraries created in the present work can in future be used to develop more sensitive targeted methods of insect PAP identification and quantification in commercial feed mixtures.

Biodiversity seen through the perspective of insects: 10 simple rules on methodological choices and experimental design for genomic studies

Massively parallel DNA sequencing opens up opportunities for bridging multiple temporal and spatial dimensions in biodiversity research, thanks to its efficiency to recover millions of nucleotide polymorphisms. Here, we identify the current status, discuss the main challenges, and look into future perspectives on biodiversity genomics focusing on insects, which arguably constitute the most diverse and ecologically important group among all animals. We suggest 10 simple rules that provide a succinct step-by-step guide and best-practices to anyone interested in biodiversity research through the study of insect genomics. To this end, we review relevant literature on biodiversity and evolutionary research in the field of entomology. Our compilation is targeted at researchers and students who may not yet be specialists in entomology or molecular biology. We foresee that the genomic revolution and its application to the study of non-model insect lineages will represent a major leap to our understanding of insect diversity.

Rhodnius prolixus: Identification of missing components of the IMD immune signaling pathway and functional characterization of its role in eliminating bacteria

The innate immune system in insects is regulated by specific signalling pathways. Most immune related pathways were identified and characterized in holometabolous insects such as Drosophila melanogaster, and it was assumed they would be highly conserved in all insects. The hemimetabolous insect, Rhodnius prolixus, has served as a model to study basic insect physiology, but also is a major vector of the human parasite, Trypanosoma cruzi, that causes 10,000 deaths annually. The publication of the R. prolixus genome revealed that one of the main immune pathways, the Immune-deficiency pathway (IMD), was incomplete and probably non-functional, an observation shared with other hemimetabolous insects including the pea aphid (Acyrthosiphon pisum) and the bedbug (Cimex lectularius). It was proposed that the IMD pathway is inactive in R. prolixus as an adaptation to prevent eliminating beneficial symbiont gut bacteria. We used bioinformatic analyses based on reciprocal BLAST and HMM-profile searches to find orthologs for most of the "missing" elements of the IMD pathway and provide data that these are regulated in response to infection with Gram-negative bacteria. We used RNAi strategies to demonstrate the role of the IMD pathway in regulating the expression of specific antimicrobial peptides (AMPs) in the fat body of R. prolixus. The data indicate that the IMD pathway is present and active in R. prolixus, which opens up new avenues of research on R. prolixus-T. cruzi interactions.

UDP-Glycosyltransferase Genes in the Striped Rice Stem Borer, Chilo suppressalis (Walker), and Their Contribution to Chlorantraniliprole Resistance

Uridine diphosphate glycosyltransferases (UGTs) are multifunctional detoxification enzymes, which are involved in metabolizing various chemicals and contribute to the development of insecticide resistance. However, the possible roles of UGTs in chlorantraniliprole resistance in Chilo suppressalis have rarely been studied in detail. Based on genome data, 24 UGT genes in C. suppressalis belonging to 11 families were identified, which were designated by the UGT nomenclature committee. Synergism assay data suggested that UGTs are potentially involved in chlorantraniliprole resistance in C. suppressalis.CsUGT40AL1 and CsUGT33AG3 were significantly overexpressed in the chlorantraniliprole resistant strain (12.36- and 5.34-fold, respectively). The two UGTs were highly expressed in the larval Malpighian tubules, fat body, and midgut; however, expression was lowest in the head. Injection of individual dsRNAs reduced the expression of the two target genes (by 69.34% and 48.74%, respectively) and caused significant higher larval mortality (81.33% and 54.67%, respectively). Overexpression of CsUGT40AL1 and CsUGT33AG3 was potentially involved in chlorantraniliprole resistance in C. suppressalis, as confirmed by the RNAi assay. Our findings suggest that overexpression of UGTs may contribute to chlorantraniliprole resistance in C. suppressalis.

Gene Selection and Evolutionary Modeling Affect Phylogenomic Inference of Neuropterida Based on Transcriptome Data

Neuropterida is a super order of Holometabola that consists of the orders Megaloptera (dobsonflies, fishflies, and alderflies), Neuroptera (lacewings) and Raphidioptera (snakeflies). Several proposed higher-level relationships within Neuropterida, such as the relationships between the orders or between the families, have been extensively debated. To further understand the evolutionary history of Neuropterida, we conducted phylogenomic analyses of all 13 published transcriptomes of the neuropterid species, as well as of a new transcriptome of the fishfly species Ctenochauliodes similis of Liu and Yang, 2006 (Megaloptera: Corydalidae: Chauliodinae) that we sequenced. Our phylogenomic data matrix contained 1392 ortholog genes from 22 holometabolan species representing six families from Neuroptera, two families from Raphidioptera, and two families from Megaloptera as the ingroup taxa, and nine orders of Holometabola as outgroups. Phylogenetic reconstruction was performed using both concatenation and coalescent-based approaches under a site-homogeneous model as well as under a site-heterogeneous model. Surprisingly, analyses using the site-homogeneous model strongly supported a paraphyletic Neuroptera, with Coniopterygidae assigned as the sister group of all other Neuropterida. In contrast, analyses using the site-heterogeneous model recovered Neuroptera as monophyletic. The monophyly of Neuroptera was also recovered in concatenation and coalescent-based analyses using genes with stronger phylogenetic signals [i.e., higher average bootstrap support (ABS) values and higher relative tree certainty including all conflicting bipartitions (RTCA) values] under the site-homogeneous model. The present study illustrated how both data selection and model selection influence phylogenomic analyses of large-scale data matrices comprehensively.

The Evolutionary History and Functional Divergence of Trehalase (treh) Genes in Insects

Trehalases (treh) have been found in different organisms, such as bacteria, fungi, yeast, nematodes, insects, vertebrates, and plants. Their biochemical properties are extremely variable and not yet fully understood. Gene expression patterns have shown differences among insect species suggesting a potential functional diversification of trehalase enzymes during their evolution. A second gene family encoding for enzymes with hypothetical trehalase activity has been repeatedly annotated in insect genome as acid trehalases/acid trehalase-like (ath), but its functional role is still not clear. The currently available large amount of genomic data from many insect species may enable a better understanding of the evolutionary history, phylogenetic relationships and possible roles of trehalase encoding genes in this taxon. The aim of the present study is to infer the evolutionary history of trehalases and acid trehalase genes in insects and analyze the trehalase functional divergence during their evolution, combining phylogenetic and genomic synteny/colinearity analyses.

Gilbert MT, Bafna V, Mirarab S. Skmer: assembly-free and alignment-free sample identification using genome skims

The ability to inexpensively describe taxonomic diversity is critical in this era of rapid climate and biodiversity changes. The recent genome-skimming approach extends current barcoding practices beyond short markers by applying low-pass sequencing and recovering whole organelle genomes computationally. This approach discards the nuclear DNA, which constitutes the vast majority of the data. In contrast, we suggest using all unassembled reads. We introduce an assembly-free and alignment-free tool, Skmer, to compute genomic distances between the query and reference genome skims. Skmer shows excellent accuracy in estimating distances and identifying the closest match in reference datasets.

Mechanistic insights into the evolution of DUF26-containing proteins in land plants

Large protein families are a prominent feature of plant genomes and their size variation is a key element for adaptation. However, gene and genome duplications pose difficulties for functional characterization and translational research. Here we infer the evolutionary history of the DOMAIN OF UNKNOWN FUNCTION (DUF) 26-containing proteins. The DUF26 emerged in secreted proteins. Domain duplications and rearrangements led to the appearance of CYSTEINE-RICH RECEPTOR-LIKE PROTEIN KINASES (CRKs) and PLASMODESMATA-LOCALIZED PROTEINS (PDLPs). The DUF26 is land plant-specific but structural analyses of PDLP ectodomains revealed strong similarity to fungal lectins and thus may constitute a group of plant carbohydrate-binding proteins. CRKs expanded through tandem duplications and preferential retention of duplicates following whole genome duplications, whereas PDLPs evolved according to the dosage balance hypothesis. We propose that new gene families mainly expand through small-scale duplications, while fractionation and genetic drift after whole genome multiplications drive families towards dosage balance.

Sex loss in insects: causes of asexuality and consequences for genomes

Boasting a staggering diversity of reproductive strategies, insects provide attractive models for the comparative study of the causes and consequences of transitions to asexuality. We provide an overview of some contemporary studies of reproductive systems in insects and compile an initial database of asexual insect genome resources. Insect systems have already yielded some important insights into various mechanisms by which sex is lost, including genetic, endosymbiont-mediated, and hybridization. Studies of mutation and substitution after loss of sex provide the strongest empirical support for hypothesized effects of asexuality, whereas there is mixed evidence for ecological hypotheses such as increased parasite load and altered niche breadth in asexuals. Most hypotheses have been explored in a select few taxa (e.g. stick insects, aphids), such that much of the great taxonomic breadth of insects remain understudied. Given the variation in the proximate causes of asexuality in insects, we argue for expanding the taxonomic breadth of study systems. Despite some challenges for investigating sex in insects, the increasing cost-effectiveness of genomic sequencing makes data generation for closely-related asexual and sexual lineages increasingly feasible.

Genome-wide identification and characterization of HSP gene superfamily in whitefly (Bemisia tabaci) and expression profiling analysis under temperature stress

Heat shock proteins (HSP) are essential molecular chaperones that play important roles in the stress stimulation of insects. Bemisia tabaci, a phloem feeder and invasive species, can cause extensive crop damage through direct feeding and transmission of plant viruses. Here we employed comprehensive genomics approaches to identity HSP superfamily members in the Middle East Asia Minor 1 whitefly genome. In total, we identified 26 Hsp genes, including three Hsp90, 17 Hsp70, one Hsp60 and five sHSP (small heat shock protein) genes. The HSP gene superfamily of whitefly is expanded compared with the other five insects surveyed here. The gene structures among the same families are relatively conserved. Meanwhile, the motif compositions and secondary structures of BtHsp proteins were predicted. In addition, quantitative polymerase chain reaction analysis showed that the expression patterns of BtHsp gene superfamily were diverse across different tissues of whiteflies. Most Hsp genes were induced or repressed by thermal stress (40°C) and cold treatment (4°C) in whitefly. Silencing the expression of BtHsp70-6 significantly decreased the survival rate of whitefly under 45°C. All the results showed the Hsps conferred thermo-tolerance or cold-tolerance to whiteflies that protect them from being affected by detrimental temperature conditions. Our observations highlighted the molecular evolutionary properties and the response mechanism to temperature assaults of Hsp genes in whitefly.

Evolution of salivary glue genes in Drosophila species

Background

At the very end of the larval stage Drosophila expectorate a glue secreted by their salivary glands to attach themselves to a substrate while pupariating. The glue is a mixture of apparently unrelated proteins, some of which are highly glycosylated and possess internal repeats. Because species adhere to distinct substrates (i.e. leaves, wood, rotten fruits), glue genes are expected to evolve rapidly.

Results

We used available genome sequences and PCR-sequencing of regions of interest to investigate the glue genes in 20 Drosophila species. We discovered a new gene in addition to the seven glue genes annotated in D. melanogaster. We also identified a phase 1 intron at a conserved position present in five of the eight glue genes of D. melanogaster, suggesting a common origin for those glue genes. A slightly significant rate of gene turnover was inferred. Both the number of repeats and the repeat sequence were found to diverge rapidly, even between closely related species. We also detected high repeat number variation at the intrapopulation level in D. melanogaster.

Conclusion

Most conspicuous signs of accelerated evolution are found in the repeat regions of several glue genes.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1364-9) contains supplementary material, which is available to authorized users.

The functions of caspase in whitefly Bemisia tabaci apoptosis in response to ultraviolet irradiation

Whiteflies (Bemisia tabaci) are phloem feeders, and some invasive species are composed of cryptic species complexes that cause extensive crop damage, particularly via the direct transmission of plant viruses. Apoptosis is a type of programmed cell death essential for organismal development and tissue homeostasis. The caspases belong to a family of cysteine proteases that play a central role in the initiation of apoptosis in many organisms. Here, we employed a comprehensive genomics approach to identity caspases in B. tabaci Middle East Asia Minor 1 (MEAM1), an invasive whitefly that carries a cryptic species complex that is devastating to crops. Four caspase genes were identified, and their motif compositions were predicted. Structures were relatively conserved in both putative effector and initiator caspases. Expression patterns of caspase genes differed across insect developmental stages. Three caspase genes were induced immediately after ultraviolet (UV) treatment. Expression levels of Bt-caspase-1 and Bt-caspase-3b increased in the midgut and salivary glands during apoptosis induced by UV treatments, whereas silencing of both genes reduced UV-triggered apoptosis. Our study demonstrates that Bt-caspase-1 and Bt-caspase-3b, respectively, act as putative initiator and effector apoptotic caspases in the MEAM1 whitefly.

The achaete-scute complex contains a single gene that controls bristle development in the semi-aquatic bugs

The semi-aquatic bugs (Heteroptera, Gerromorpha) conquered water surfaces worldwide and diversified to occupy puddles, ponds, streams, lakes, mangroves and even oceans. Critical to this lifestyle is the evolution of sets of hairs that allow these insects to maintain their body weight on the water surface and protect the animals against wetting and drowning. In addition, the legs of these insects are equipped with various grooming combs that are important for cleaning and tidying the hair layers for optimal functional efficiency. Here we show that the hairs covering the legs of water striders represent innervated bristles. Genomic and transcriptomic analyses revealed that in water striders the achaete–scute complex, known to control bristle development in flies, contains only the achaete–scute homologue (ASH) gene owing to the loss of the gene asense. Using RNA interference, we show that ASH plays a pivotal role in the development of both bristles and grooming combs in water striders. Our data suggest that the ASH locus may have contributed to the adaptation to water surface lifestyle through shaping the hydrophobic bristles that prevent water striders from wetting and allow them to exploit water surface tension.

Long-read sequence assembly of the firefly Pyrocoelia pectoralis genome

Background

Fireflies are a family of insects within the beetle order Coleoptera, or winged beetles, and they are one of the most well-known and loved insect species because of their bioluminescence. However, the firefly is in danger of extinction because of the massive destruction of its living environment. In order to improve the understanding of fireflies and protect them effectively, we sequenced the whole genome of the terrestrial firefly Pyrocoelia pectoralis.

Findings

Here, we developed a highly reliable genome resource for the terrestrial firefly Pyrocoelia pectoralis (E. Oliv., 1883; Coleoptera: Lampyridae) using single molecule real time (SMRT) sequencing on the PacBio Sequel platform. In total, 57.8 Gb of long reads were generated and assembled into a 760.4-Mb genome, which is close to the estimated genome size and covered 98.7% complete and 0.7% partial insect Benchmarking Universal Single-Copy Orthologs. The k-mer analysis showed that this genome is highly heterozygous. However, our long-read assembly demonstrates continuousness with a contig N50 length of 3.04 Mb and the longest contig length of 13.69 Mb. Furthermore, 135 589 SSRs and 341 Mb of repeat sequences were detected. A total of 23 092 genes were predicted; 88.44% of genes were annotated with one or more related functions.

Conclusions

We assembled a high-quality firefly genome, which will not only provide insights into the conservation and biodiversity of fireflies, but also provide a wealth of information to study the mechanisms of their sexual communication, bio-luminescence, and evolution.

The genomic features of parasitism, Polyembryony and immune evasion in the endoparasitic wasp Macrocentrus cingulum

BACKGROUND

Parasitoid wasps are well-known natural enemies of major agricultural pests and arthropod borne diseases. The parasitoid wasp Macrocentrus cingulum (Hymenoptera: Braconidae) has been widely used to control the notorious insect pests Ostrinia furnacalis (Asian Corn Borer) and O. nubilalis (European corn borer). One striking phenomenon exhibited by M. cingulum is polyembryony, the formation of multiple genetically identical offspring from a single zygote. Moreover, M. cingulum employs a passive parasitic strategy by preventing the host's immune system from recognizing the embryo as a foreign body. Thus, the embryos evade the host's immune system and are not encapsulated by host hemocytes. Unfortunately, the mechanism of both polyembryony and immune evasion remains largely unknown.

RESULTS

We report the genome of the parasitoid wasp M. cingulum. Comparative genomics analysis of M. cingulum and other 11 insects were conducted, finding some gene families with apparent expansion or contraction which might be linked to the parasitic behaviors or polyembryony of M. cingulum. Moreover, we present the evidence that the microRNA miR-14b regulates the polyembryonic development of M. cingulum by targeting the c-Myc Promoter-binding Protein 1 (MBP-1), histone-lysine N-methyltransferase 2E (KMT2E) and segmentation protein Runt. In addition, Hemomucin, an O-glycosylated transmembrane protein, protects the endoparasitoid wasp larvae from being encapsulated by host hemocytes. Motif and domain analysis showed that only the hemomucin in two endoparasitoids, M. cingulum and Venturia canescens, possessing the ability of passive immune evasion has intact mucin domain and similar O-glycosylation patterns, indicating that the hemomucin is a key factor modulating the immune evasion.

CONCLUSIONS

The microRNA miR-14b participates in the regulation of polyembryonic development, and the O-glycosylation of the mucin domain in the hemomucin confers the passive immune evasion in this wasp. These key findings provide new insights into the polyembryony and immune evasion.

An Integrated Molecular Database on Indian Insects

MOlecular Database on Indian Insects (MODII) is an online database linking several databases like Insect Pest Info, Insect Barcode Information System (IBIn), Insect Whole Genome sequence, Other Genomic Resources of National Bureau of Agricultural Insect Resources (NBAIR), Whole Genome sequencing of Honey bee viruses, Insecticide resistance gene database and Genomic tools. This database was developed with a holistic approach for collecting information about phenomic and genomic information of agriculturally important insects. This insect resource database is available online for free at http://cib.res.in.

AVAILABILITY

http://cib.res.in/.

Lepidoptera genomes: current knowledge, gaps and future directions

Butterflies and moths (Lepidoptera) are one of the most ecologically diverse and speciose insect orders. With recent advances in genomics, new Lepidoptera genomes are regularly being sequenced, and many of them are playing principal roles in genomics studies, particularly in the fields of phylo-genomics and functional genomics. Thus far, assembled genomes are only available for <10 of the 43 Lepidoptera superfamilies. Nearly all are model species, found in the speciose clade Ditrysia. Community support for Lepidoptera genomics is growing with successful management and dissemination of data and analytical tools in centralized databases. With genomic studies quickly becoming integrated with ecological and evolutionary research, the Lepidoptera community will unquestionably benefit from new high-quality reference genomes that are more evenly distributed throughout the order.

Genomes of Diptera

Diptera (true flies) are among the most diverse holometabolan insect orders and were the first eukaryotic order to have a representative genome fully sequenced. 110 fly species have publically available genome assemblies and many hundreds of population-level genomes have been generated in the model organisms Drosophila melanogaster and the malaria mosquito Anopheles gambiae. Comparative genomics carried out in a phylogenetic context is illuminating many aspects of fly biology, providing unprecedented insight into variability in genome structure, gene content, genetic mechanisms, and rates and patterns of evolution in genes, populations, and species. Despite the rich availability of genomic resources in flies, there remain many fly lineages to which new genome sequencing efforts should be directed. Such efforts would be most valuable in fly families or clades that exhibit multiple origins of key fly behaviors such as blood feeding, phytophagy, parasitism, pollination, and mycophagy.

Downregulation of Chilo suppressalis alkaline phosphatase genes associated with resistance to three transgenic Bacillus thuringiensis rice lines

Insecticidal crystal (Cry) proteins produced by the bacterium Bacillus thuringiensis are highly toxic to lepidopteran pests. Strains of transgenic rice expressing cry genes have been developed that are resistant to rice pests. Understanding the mode of action of Cry toxins in rice pests will improve our ability to use them effectively as insecticides. In this study, we tested the hypothesis that alkaline phosphatases (ALPs) are involved in Cry1A, Cry2Aa and Cry1Ca toxicity in Chilo suppressalis, an important insect pest of rice crops in China. We first cloned three novel C. suppressalis alps (Csalps) from the larval midgut of C. suppressalis. RNA interference knockdown of six different Csalp genes (Csalp1, Csalp2, Csalp3, Csalp4, Csalp5 and Csalp6) showed that knockdown of three of these, Csalp1, Csalp2 and Csalp4, reduced larval mortality to the transgenic rice strain TT51, which expresses a fusion protein of Cry1Ab and Cry1Ac, whereas suppression of Csalp1, Csalp2, Csalp3, Csalp4 and Csalp6 transcripts decreased the susceptibility of larvae to the transgenic rice strain T2A-1, which expresses cry2Aa. Moreover, downregulation of Csalp1, Csalp2, Csalp3, Csalp4 and Csalp5 transcripts conferred significant tolerance to the transgenic rice strain T1C-19, which expresses cry1Ca. These results suggest that these ALPs play a key role in the toxicity of Cry1A, Cry2A and Cry1C to C. suppressalis.

WaspBase: a genomic resource for the interactions among parasitic wasps, insect hosts and plants

Insect pests reduce yield and cause economic losses, which are major problems in agriculture. Parasitic wasps are the natural enemies of many agricultural pests and thus have been widely used as biological control agents. Plants, phytophagous insects and parasitic wasps form a tritrophic food chain. Understanding the interactions in this tritrophic system should be helpful for developing parasitic wasps for pest control and deciphering the mechanisms of parasitism. However, the genomic resources for this tritrophic system are not well organized. Here, we describe the WaspBase, a new database that contains 573 transcriptomes of 35 parasitic wasps and the genomes of 12 parasitic wasps, 5 insect hosts and 8 plants. In addition, we identified long non-coding RNA, untranslated regions and 25 widely studied gene families from the genome and transcriptome data of these species. WaspBase provides conventional web services such as Basic Local Alignment Search Tool, search and download, together with several widely used tools such as profile hidden Markov model, Multiple Alignment using Fast Fourier Transform, automated alignment trimming and JBrowse. We also present a collection of active researchers in the field of parasitic wasps, which should be useful for constructing scientific networks in this field.

MGH: a genome hub for the medicinal plant maca (Lepidium meyenii)

Maca (Lepidium meyenii), a Brassicaceae herb plant originated from Andean mountains, has attracted wide interests due to its unique health benefits in reproduction and fertility. Because of its adaptation to the 4000 m high-altitude harsh environment, maca is attracting more and more attention from both crop breeders and basic biologists. After our previous release of the maca genome sequence, there's a growing need to store, query, analyze and integrate various maca resources efficiently. Here, we created Maca Genome Hub (MGH), a genomics and genetics database of maca. Currently, the MGH V1.0 harbors the genome sequence, predicted coding sequences and protein sequences, various annotations, markers and expression data. For the maca research community, we also provided the publications, researchers and related news. MGH is designed to enable users' easy access to analyze, retrieve and visualize the genomic or genetic information through a series of online tools, including the Basic Local Alignment Search Tool, the JBrowse, the query system, the synteny tool and the data downloads. These integrated heterogeneous data, tools and interfaces in MGH allow efficient mining of the latest genomics and genetics data. We hope that MGH will accelerate the research and development in maca.

Navigating the i5k Workspace@NAL: A Resource for Arthropod Genomes

The i5k Workspace@NAL is a genome database tailored toward newly sequenced arthropod genomes and their research communities. With 56 arthropod genomes and counting, the i5k Workspace strives to facilitate public data access, visualization, and community curation across arthropod species. Any researcher with an arthropod genome project who would like to take advantage of the i5k Workspace facilities is encouraged to submit their data. In this chapter, we explain how to use the i5k Workspace@NAL to submit, find, and improve arthropod genomics data.

Nucleases as a barrier to gene silencing in the cotton boll weevil, Anthonomus grandis

RNA interference (RNAi) approaches have been applied as a biotechnological tool for controlling plant insect pests via selective gene down regulation. However, the inefficiency of RNAi mechanism in insects is associated with several barriers, including dsRNA delivery and uptake by the cell, dsRNA interaction with the cellular membrane receptor and dsRNA exposure to insect gut nucleases during feeding. The cotton boll weevil (Anthonomus grandis) is a coleopteran in which RNAi-mediated gene silencing does not function efficiently through dsRNA feeding, and the factors involved in the mechanism remain unknown. Herein, we identified three nucleases in the cotton boll weevil transcriptome denoted AgraNuc1, AgraNuc2, and AgraNuc3, and the influences of these nucleases on the gene silencing of A. grandis chitin synthase II (AgraChSII) were evaluated through oral dsRNA feeding trials. A phylogenetic analysis showed that all three nucleases share high similarity with the DNA/RNA non-specific endonuclease family of other insects. These nucleases were found to be mainly expressed in the posterior midgut region of the insect. Two days after nuclease RNAi-mediated gene silencing, dsRNA degradation by the gut juice was substantially reduced. Notably, after nucleases gene silencing, the orally delivered dsRNA against the AgraChSII gene resulted in improved gene silencing efficiency when compared to the control (non-silenced nucleases). The data presented here demonstrates that A. grandis midgut nucleases are effectively one of the main barriers to dsRNA delivery and emphasize the need to develop novel RNAi delivery strategies focusing on protecting the dsRNA from gut nucleases and enhancing its oral delivery and uptake to crop insect pests.

Multiple miRNAs jointly regulate the biosynthesis of ecdysteroid in the holometabolous insects, Chilo suppressalis

The accurate rise and fall of active hormones is important for insect development. The ecdysteroids must be cleared in a timely manner. However, the mechanism of suppressing the ecdysteroid biosynthesis at the right time remains unclear. Here, we sequenced a small RNA library of Chilo suppressalis and identified 300 miRNAs in this notorious rice insect pest. Microarray analysis yielded 54 differentially expressed miRNAs during metamorphosis development. Target prediction and in vitro dual-luciferase assays confirmed that seven miRNAs (two conserved and five novel miRNAs) jointly targeted three Halloween genes in the ecdysteroid biosynthesis pathway. Overexpression of these seven miRNAs reduced the titer of 20-hydroxyecdysone (20E), induced mortality, and retarded development, which could be rescued by treatment with 20E. Comparative analysis indicated that the miRNA regulation of metamorphosis development is a conserved process but that the miRNAs involved are highly divergent. In all, we present evidence that both conserved and lineage-specific miRNAs have crucial roles in regulating development in insects by controlling ecdysteroid biosynthesis, which is important for ensuring developmental convergence and evolutionary diversity.

Genome-wide identification of long non-coding RNA genes and their association with insecticide resistance and metamorphosis in diamondback moth, Plutella xylostella

Long non-coding RNA (lncRNA) is a class of noncoding RNA >200 bp in length that has essential roles in regulating a variety of biological processes. Here, we constructed a computational pipeline to identify lncRNA genes in the diamondback moth (Plutella xylostella), a major insect pest of cruciferous vegetables. In total, 3,324 lncRNAs corresponding to 2,475 loci were identified from 13 RNA-Seq datasets, including samples from parasitized, insecticide-resistant strains and different developmental stages. The identified P. xylostella lncRNAs had shorter transcripts and fewer exons than protein-coding genes. Seven out of nine randomly selected lncRNAs were validated by strand-specific RT-PCR. In total, 54-172 lncRNAs were specifically expressed in the insecticide resistant strains, among which one lncRNA was located adjacent to the sodium channel gene. In addition, 63-135 lncRNAs were specifically expressed in different developmental stages, among which three lncRNAs overlapped or were located adjacent to the metamorphosis-associated genes. These lncRNAs were either strongly or weakly co-expressed with their overlapping or neighboring mRNA genes. In summary, we identified thousands of lncRNAs and presented evidence that lncRNAs might have key roles in conferring insecticide resistance and regulating the metamorphosis development in P. xylostella.

Bugs Are Not to Be Silenced: Small RNA Pathways and Antiviral Responses in Insects

Like every other organism on Earth, insects are infected with viruses, and they rely on RNA interference (RNAi) mechanisms to circumvent viral infections. A remarkable characteristic of RNAi is that it is both broadly acting, because it is triggered by double-stranded RNA molecules derived from virtually any virus, and extremely specific, because it targets only the particular viral sequence that initiated the process. Reviews covering the different facets of the RNAi antiviral immune response in insects have been published elsewhere. In this review, we build a framework to guide future investigation. We focus on the remaining questions and avenues of research that need to be addressed to move the field forward, including issues such as the activity of viral suppressors of RNAi, comparative genomics, the development of detailed maps of the subcellular localization of viral replication complexes with the RNAi machinery, and the regulation of the antiviral RNAi response.

Large-scale analysis reveals that the genome features of simple sequence repeats are generally conserved at the family level in insects

BACKGROUND

Simple sequence repeats (SSR), also called microsatellites, have been widely used as genetic markers, and have been extensively studied in some model insects. At present, the genomes of more than 100 insect species are available. However, the features of SSRs in most insect genomes remain largely unknown.

RESULTS

We identified 15.01 million SSRs across 136 insect genomes. The number of identified SSRs was positively associated with genome size in insects, but the frequency and density per megabase of genomes were not. Most insect SSRs (56.2-93.1%) were perfect (no mismatch). Imperfect (at least one mismatch) SSRs (average length 22-73 bp) were longer than perfect SSRs (16-30 bp). The most abundant insect SSRs were the di- and trinucleotide types, which accounted for 27.2% and 22.0% of all SSRs, respectively. On average, 59.1%, 36.8%, and 3.7% of insect SSRs were located in intergenic, intronic, and exonic regions, respectively. The percentages of various types of SSRs were similar among insects from the same family. However, they were dissimilar among insects from different families within orders. We carried out a phylogenetic analysis using the SSR frequencies. Species from the same family were generally clustered together in the evolutionary tree. However, insects from the same order but not in the same family did not cluster together. These results indicated that although SSRs undergo rapid expansions and contractions in different populations of the same species, the general genomic features of insect SSRs remain conserved at the family level.

CONCLUSION

Millions of insect SSRs were identified and their genome features were analyzed. Most insect SSRs were perfect and were located in intergenic regions. We presented evidence that the variance of insect SSRs accumulated after the differentiation of insect families.

Insights into the molecular evolution of peptidase inhibitors in arthropods

Peptidase inhibitors are key proteins involved in the control of peptidases. In arthropods, peptidase inhibitors modulate the activity of peptidases involved in endogenous physiological processes and peptidases of the organisms with which they interact. Exploring available arthropod genomic sequences is a powerful way to obtain the repertoire of peptidase inhibitors in every arthropod species and to understand the evolutionary mechanisms involved in the diversification of this kind of proteins. A genomic comparative analysis of peptidase inhibitors in species belonging to different arthropod taxonomic groups was performed. The results point out: i) species or clade-specific presence is shown for several families of peptidase inhibitors; ii) multidomain peptidase inhibitors are commonly found in many peptidase inhibitor families; iii) several families have a wide range of members in different arthropod species; iv) several peptidase inhibitor families show species-specific (or clade-specific) gene family expansions; v) functional divergence may be assumed for particular clades; vi) passive expansions may be used by natural selection to fix adaptations. In conclusion, conservation and divergence of duplicated genes and the potential recruitment as peptidase inhibitors of proteins from other families are the main mechanisms used by arthropods to fix diversity. This diversity would be associated to the control of target peptidases and, as consequence, to adapt to specific environments.

Silencing of juvenile hormone epoxide hydrolase gene (Nljheh) enhances short wing formation in a macropterous strain of the brown planthopper, Nilaparvata lugens

The rice brown planthopper, Nilaparvata lugens, is an important migratory pest in many rice planting areas of Asia. The typical wing dimorphism of N. lugens gives them flexibility to adapt to different environmental cues. As an important hormone in the insect's endocrine regulation, juvenile hormone (JH) has previously been shown to participate in the wing morph determination of N. lugens. In this paper, we investigated the possible wing morph determination roles of two JH metabolic enzymes, JH esterase (JHE) and JH epoxide hydrolase (JHEH). A 1957-bp full-length cDNA sequence encoding JHEH in N. lugens (NlJHEH) was first cloned from a hemipteran insect. Except for an uncertain transmembrane segment prediction, the deduced 454-amino-acid sequence of Nljheh has all of the conserved domains of JHEHs such as the H¹⁴⁷GWP¹⁵⁰, Tyr293 and Tyr368 motif corresponding to the oxyanion hole and the residues Asp222, Glu398, and His425 in the catalytic triad. qRT-PCR results showed that both Nljhe and Nljheh had different expression timeframes between a predominantly brachypterous strain (BS) and a macropterous strain (MS) of N. lugens, indicating that these two enzymes may participate in wing dimorphism regulation in brown planthopper. Silencing Nljheh expression by dsRNA injection enhanced short wing formation in the macropterous strain of N. lugens, while the brachypterizing individuals were mainly females. Compared to the dsgfp injection control, silencing Nljhe had no brachypterizing effect. Our results indicated that NlJHEH plays an important role in the wing morph determination of N. lugens.

Perspectives on Gene Regulatory Network Evolution

Animal development proceeds through the activity of genes and their cis-regulatory modules (CRMs) working together in sets of gene regulatory networks (GRNs). The emergence of species-specific traits and novel structures results from evolutionary changes in GRNs. Recent work in a wide variety of animal models, and particularly in insects, has started to reveal the modes and mechanisms of GRN evolution. I discuss here various aspects of GRN evolution and argue that developmental system drift (DSD), in which conserved phenotype is nevertheless a result of changed genetic interactions, should regularly be viewed from the perspective of GRN evolution. Advances in methods to discover related CRMs in diverse insect species, a critical requirement for detailed GRN characterization, are also described.

Large-scale identification of differentially expressed genes during pupa development reveals solute carrier gene is essential for pupal pigmentation in Chilo suppressalis

Insects undergo metamorphosis, involving an abrupt change in body structure through cell growth and differentiation. Rice stem stripped borer (SSB), Chilo suppressalis, is one of the most destructive rice pests. However, little is known about the regulation mechanism of metamorphosis development in this notorious insect pest. Here, we studied the expression of 22,197 SSB genes at seven time points during pupa development with a customized microarray, identifying 622 differentially expressed genes (DEG) during pupa development. Gene ontology (GO) analysis of these DEGs indicated that the genes related to substance metabolism were highly expressed in the early pupa, which participate in the physiological processes of larval tissue disintegration at these stages. In comparison, highly expressed genes in the late pupal stages were mainly associated with substance biosynthesis, consistent with adult organ formation at these stages. There were 27 solute carrier (SLC) genes that were highly expressed during pupa development. We knocked down SLC22A3 at the prepupal stage, demonstrating that silencing SLC22A3 induced a deficiency in pupa stiffness and pigmentation. The RNAi-treated individuals had white and soft pupa, suggesting that this gene has an essential role in pupal development.

Deep sequencing of the prothoracic gland transcriptome reveals new players in insect ecdysteroidogenesis

Ecdysteroids are steroid hormones that induce molting and determine developmental timing in arthropods. In insect larva, the prothoracic gland (PG) is a major organ for ecdysone synthesis and release. Released ecdysone is converted into the active form, 20-hydroxyecdysone (20E) in the peripheral tissues. All processes from ecdysone synthesis and release from the PG to its conversion to 20E are called ecdysteroidogenesis and are under the regulation of numerous factors expressed in the PG and peripheral tissues. Classical genetic approaches and recent transcriptomic screening in the PG identified several genes responsible for ecdysone synthesis and release, whereas the regulatory mechanism remains largely unknown. We analyzed RNA-seq data of the silkworm Bombyx mori PG and employed the fruit fly Drosophila melanogaster GAL4/UAS binary RNAi system to comprehensively screen for genes involved in ecdysone synthesis and/or release. We found that the genes encoding δ-aminolevulinic acid synthase (CG3017/alas) and putative NAD kinase (CG33156) were highly expressed in the PG of both B. mori and D. melanogaster. Neither alas nor CG33156 RNAi-induced larvae could enter into the pupal stage, and they had a lower abundance of the active form ecdysteroids in their prolonged larval stage. These results demonstrated that alas and CG33156 are indispensable for ecdysteroidogenesis.

Distinct Roles of Met and Interacting Proteins on the Expressions of takeout Family Genes in Brown Planthopper

The takeout family genes encode relatively small proteins that are related to olfaction and are regulated by juvenile hormone (JH). The takeout genes modulate various physiological processes, such as behavioral plasticity in the migratory locust Locusta migraloria and feeding and courtship behaviors in Drosophila. Therefore, to understand the regulatory mechanism of these physiological processes, it is important to study the expressions of the takeout genes that are regulated by JH signaling. We used quantitative real-time PCR (qRTPCR) to study the role of JH signaling in the regulation of the takeout family genes in the brown planthopper Nilaparvata lugens (N. lugens) through the application of Juvenile hormone III (JHIII) and the down-regulation of key genes in the JH signaling pathway. The topical application of JHIII induced the expressions of most of the takeout family genes, and their expressions decreased 2 and 3 days after the JHIII application. Down-regulating the brown planthopper JH receptor NlMethoprene-tolerant (NlMet) and its interacting partners, NlTaiman (NlTai) and Nlß-Ftz-F1 (Nlß-Ftz), through RNAi, exhibited distinct effects on the expressions of the takeout family genes. The down-regulation of NlMet and NlKrüppel-homolog 1 (NlKr-h1) increased the expressions of the takeout family genes, while the down-regulation of the Met interacting partners NlTai and Nlß-Ftz decreased the expressions of most of the takeout family genes. This work advanced our understanding of the molecular function and the regulatory mechanism of JH signaling.

Genome Sizes of Nine Insect Species Determined by Flow Cytometry and k-mer Analysis

The flow cytometry method was used to estimate the genome sizes of nine agriculturally important insects, including two coleopterans, five Hemipterans, and two hymenopterans. Among which, the coleopteran Lissorhoptrus oryzophilus (Kuschel) had the largest genome of 981 Mb. The average genome size was 504 Mb, suggesting that insects have a moderate-size genome. Compared with the insects in other orders, hymenopterans had small genomes, which were averagely about ~200 Mb. We found that the genome sizes of four insect species were different between male and female, showing the organismal complexity of insects. The largest difference occurred in the coconut leaf beetle Brontispa longissima (Gestro). The male coconut leaf beetle had a 111 Mb larger genome than females, which might be due to the chromosome number difference between the sexes. The results indicated that insect invasiveness was not related to genome size. We also determined the genome sizes of the small brown planthopper Laodelphax striatellus (Fallén) and the parasitic wasp Macrocentrus cingulum (Brischke) using k-mer analysis with Illunima Solexa sequencing data. There were slight differences in the results from the two methods. k-mer analysis indicated that the genome size of L. striatellus was 500–700 Mb and that of M. cingulum was ~150 Mb. In all, the genome sizes information presented here should be helpful for designing the genome sequencing strategy when necessary.

Folded gastrulation and T48 drive the evolution of coordinated mesoderm internalization in flies

Gastrulation constitutes a fundamental yet diverse morphogenetic process of metazoan development. Modes of gastrulation range from stochastic translocation of individual cells to coordinated infolding of an epithelial sheet. How such morphogenetic differences are genetically encoded and whether they have provided specific developmental advantages is unclear. Here we identify two genes, folded gastrulation and t48, which in the evolution of fly gastrulation acted as a likely switch from an ingression of individual cells to the invagination of the blastoderm epithelium. Both genes are expressed and required for mesoderm invagination in the fruit fly Drosophila melanogaster but do not appear during mesoderm ingression of the midge Chironomus riparius. We demonstrate that early expression of either or both of these genes in C.riparius is sufficient to invoke mesoderm invagination similar to D.melanogaster. The possible genetic simplicity and a measurable increase in developmental robustness might explain repeated evolution of similar transitions in animal gastrulation.

DOI:http://dx.doi.org/10.7554/eLife.18318.001

In animals, gastrulation is a period of time in early development during which a sphere of cells is reorganized into an embryo with cells arranged into three distinct layers (called germ layers). The process has changed substantially during the course of evolution and thus provides a great experimental system to investigate the genetic basis for differences in animal form and shape.

As an example, true flies use at least two different mechanisms to make the middle germ layer (the mesoderm). In both cases, the mesoderm is made up of cells that move inwards from the boundary of the outer germ layer. In midges and some other flies, these cells migrate individually, while in others including fruit flies, the cells move together as a sheet. Fruit flies and midges shared their last common ancestor 250 million years ago and although the genes that make the mesoderm in fruit flies are well understood, little is known about how the mesoderm forms in midges.

Urbansky, González Avalos et al. investigated which genes were responsible for the evolutionary transition between the different types of cell migration seen in flies. The experiments identified two genes – called folded gastrulation and t48 – that seem to operate as a simple switch between the two ways that mesoderm cells migrate. Both of these genes are active in fruit fly embryos and are required for the group migration of mesoderm cells. However, the genes do not appear to play a major role in the movement of individual mesoderm cells in midges.

Further experiments demonstrate that switching on these genes in midge embryos is sufficient to invoke group mesoderm cell migrations similar to those seen in fruit flies. These findings show that it is possible to identify genetic changes that underlie substantial differences in animal form and shape over hundred million of years. The ease by which Urbansky, González Avalos et al. were able to switch between the two types of mesoderm migration might explain why similar transitions in gastrulation have evolved repeatedly in animals. The next step is to test this hypothesis in other animals.

DOI:http://dx.doi.org/10.7554/eLife.18318.002

Vector-Borne Bacterial Plant Pathogens: Interactions with Hemipteran Insects and Plants

Hemipteran insects are devastating pests of crops due to their wide host range, rapid reproduction, and ability to transmit numerous plant-infecting pathogens as vectors. While the field of plant-virus-vector interactions has flourished in recent years, plant-bacteria-vector interactions remain poorly understood. Leafhoppers and psyllids are by far the most important vectors of bacterial pathogens, yet there are still significant gaps in our understanding of their feeding behavior, salivary secretions, and plant responses as compared to important viral vectors, such as whiteflies and aphids. Even with an incomplete understanding of plant-bacteria-vector interactions, some common themes have emerged: (1) all known vector-borne bacteria share the ability to propagate in the plant and insect host; (2) particular hemipteran families appear to be incapable of transmitting vector-borne bacteria; (3) all known vector-borne bacteria have highly reduced genomes and coding capacity, resulting in host-dependence; and (4) vector-borne bacteria encode proteins that are essential for colonization of specific hosts, though only a few types of proteins have been investigated. Here, we review the current knowledge on important vector-borne bacterial pathogens, including Xylella fastidiosa, Spiroplasma spp., Liberibacter spp., and 'Candidatus Phytoplasma spp.'. We then highlight recent approaches used in the study of vector-borne bacteria. Finally, we discuss the application of this knowledge for control and future directions that will need to be addressed in the field of vector-plant-bacteria interactions.

Genome-wide analysis of chitinase genes and their varied functions in larval moult, pupation and eclosion in the rice striped stem borer, Chilo suppressalis

Some insect chitinases are required to degrade chitin and ensure successful metamorphosis. Although chitinase genes have been well characterized in several model insects, no reports exist for the rice striped stem borer, Chilo suppressalis, a highly destructive pest that causes huge yield losses in rice production. Here, we conducted a genome-level analysis of chitinase genes in C. suppressalis. After amplification of full-length transcripts with rapid amplification of cDNA ends, we identified 12 chitinase genes in C. suppressalis. All these genes had the conserved domains and motifs of glycoside hydrolase family 18 and grouped phylogenetically into five subgroups. C. suppressalis chitinase 1 (CsCht1) was highly expressed in late pupae, whereas CsCht3 was abundant in early pupae. Both CsCht2 and CsCht4 were highly expressed in larvae. CsCht2 was abundant specifically in the third-instar larvae and CsCht4 showed periodic high expression in 2- to 5-day-old larvae in each instar. Tissue specific expression analysis indicated that CsCht1 and CsCht3 were highly expressed in epidermis whereas CsCht2 and CsCht4 were specifically abundant in the midgut. Knockdown of CsCht1 resulted in adults with curled wings, indicating that CsCht1 might have an important role in wing expansion. Silencing of CsCht2 or CsCht4 arrested moulting, suggesting essential roles in larval development. When the expression of CsCht3 was interfered, defects in pupation occurred. Overall, we provide here the first catalogue of chitinase genes in the rice striped stem borer and have elucidated the functions of four chitinases in metamorphosis.

The 2016 database issue of Nucleic Acids Research and an updated molecular biology database collection

The 2016 Database Issue of Nucleic Acids Research starts with overviews of the resources provided by three major bioinformatics centers, the U.S. National Center for Biotechnology Information (NCBI), the European Bioinformatics Institute (EMBL-EBI) and Swiss Institute for Bioinformatics (SIB). Also included are descriptions of 62 new databases and updates on 95 databases that have been previously featured in NAR plus 17 previously described elsewhere. A number of papers in this issue deal with resources on nucleic acids, including various kinds of non-coding RNAs and their interactions, molecular dynamics simulations of nucleic acid structure, and two databases of super-enhancers. The protein database section features important updates on the EBI's Pfam, PDBe and PRIDE databases, as well as a variety of resources on pathways, metabolomics and metabolic modeling. This issue also includes updates on popular metagenomics resources, such as MG-RAST, EBI Metagenomics, and probeBASE, as well as a newly compiled Human Pan-Microbe Communities database. A significant fraction of the new and updated databases are dedicated to the genetic basis of disease, primarily cancer, and various aspects of drug research, including resources for patented drugs, their side effects, withdrawn drugs, and potential drug targets. A further six papers present updated databases of various antimicrobial and anticancer peptides. The entire Database Issue is freely available online on the Nucleic Acids Research website (http://nar.oxfordjournals.org/). The NAR online Molecular Biology Database Collection, http://www.oxfordjournals.org/nar/database/c/, has been updated with the addition of 88 new resources and removal of 23 obsolete websites, which brought the current listing to 1685 databases.

Small RNAs from Bemisia tabaci Are Transferred to Solanum lycopersicum Phloem during Feeding

The phloem-feeding whitefly Bemisia tabaci is a serious pest to a broad range of host plants, including many economically important crops such as tomato. These insects serve as a vector for various devastating plant viruses. It is known that whiteflies are capable of manipulating host-defense responses, potentially mediated by effector molecules in the whitefly saliva. We hypothesized that, beside putative effector proteins, small RNAs (sRNA) are delivered by B. tabaci into the phloem, where they may play a role in manipulating host plant defenses. There is already evidence to suggest that sRNAs can mediate the host-pathogen dialogue. It has been shown that Botrytis cinerea, the causal agent of gray mold disease, takes advantage of the plant sRNA machinery to selectively silence host genes involved in defense signaling. Here we identified sRNAs originating from B. tabaci in the phloem of tomato plants on which they are feeding. sRNAs were isolated and sequenced from tomato phloem of whitefly-infested and control plants as well as from the nymphs themselves, control leaflets, and from the infested leaflets. Using stem-loop RT-PCR, three whitefly sRNAs have been verified to be present in whitefly-infested leaflets that were also present in the whitefly-infested phloem sample. Our results show that whitefly sRNAs are indeed present in tomato tissues upon feeding, and they appear to be mobile in the phloem. Their role in the host-insect interaction can now be investigated.