The first complete mitochondrial genome of marigold pest thrips, Neohydatothrips samayunkur (Sericothripinae) and comparative analysis

Complete mitogenomes from the order Thysanoptera are limited to representatives of the subfamily Thripinae. Therefore, in the present study, we sequenced the mitochondrial genome of Neohydatothrips samayunkur (15,295 bp), a member of subfamily Sericothripinae. The genome possesses the canonical 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and two ribosomal RNA genes (rRNAs) as well as two putative control regions (CRs). The majority strand was 77.42% A + T content, and 22.58% G + C with weakly positive AT skew (0.04) and negative GC skew (-0.03). The majority of PCGs start with ATN codons as observed in other insect mitochondrial genomes. The GCG codon (Alanine) was not used in N. samayunkur. Most tRNAs have the typical cloverleaf secondary structure, however the DHU stem and loop were absent in trnV and trnS1, while the TΨC loop was absent in trnR and trnT. The two putative control regions (CR1 and CR2) show 99% sequence similarity indicated a possible duplication, and shared 57 bp repeats were identified. N. samayunkur showed extensive gene rearrangements, with 11 PCGs, 22 tRNAs, and two rRNAs translocated when compared to the ancestral insect. The gene trnL2 was separated from the 'trnL2-cox2' gene block, which is a conserved, ancestral gene order found in all previously sequenced thrips mitogenomes. Both maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees resulted in similar topologies. The phylogenetic position of N. samayunkur indicates that subfamily Sericothripinae is sister to subfamily Thripinae. More molecular data from different taxonomic groups is needed to understand thrips phylogeny and evolution.

Wolbachia Infection Associated with Increased Recombination in Drosophila

Wolbachia is a maternally-transmitted endosymbiotic bacteria that infects a large diversity of arthropod and nematode hosts. Some strains of Wolbachia are parasitic, manipulating host reproduction to benefit themselves, while other strains of Wolbachia exhibit obligate or facultative mutualisms with their host. The effects of Wolbachia on its host are many, though primarily relate to host immune and reproductive function. Here we test the hypothesis that Wolbachia infection alters the frequency of homologous recombination during meiosis. We use D. melanogaster as a model system, and survey recombination in eight wild-derived Wolbachia-infected (strain wMel) and Wolbachia-uninfected strains, controlling for genotype. We measure recombination in two intervals of the genome. Our results indicate that Wolbachia infection is associated with increased recombination in one genomic interval and not the other. The effect of Wolbachia infection on recombination is thus heterogenous across the genome. Our data also indicate a reproductive benefit of Wolbachia infection; infected females show higher fecundity than their uninfected genotypic controls. Given the prevalence of Wolbachia infection in natural populations, our findings suggest that Wolbachia infection is likely to contribute to recombination rate and fecundity variation among individuals in nature.

Transforming insect population control with precision guided sterile males with demonstration in flies

The sterile insect technique (SIT) is an environmentally safe and proven technology to suppress wild populations. To further advance its utility, a novel CRISPR-based technology termed precision guided SIT (pgSIT) is described. PgSIT mechanistically relies on a dominant genetic technology that enables simultaneous sexing and sterilization, facilitating the release of eggs into the environment ensuring only sterile adult males emerge. Importantly, for field applications, the release of eggs will eliminate burdens of manually sexing and sterilizing males, thereby reducing overall effort and increasing scalability. Here, to demonstrate efficacy, we systematically engineer multiple pgSIT systems in Drosophila which consistently give rise to 100% sterile males. Importantly, we demonstrate that pgSIT-generated sterile males are fit and competitive. Using mathematical models, we predict pgSIT will induce substantially greater population suppression than can be achieved by currently-available self-limiting suppression technologies. Taken together, pgSIT offers to potentially transform our ability to control insect agricultural pests and disease vectors.

Spatiotemporally controlled genetic perturbation for efficient large-scale studies of cell non-autonomous effects

Studies in genetic model organisms have revealed much about the development and pathology of complex tissues. Most have focused on cell-intrinsic gene functions and mechanisms. Much less is known about how transformed, or otherwise functionally disrupted, cells interact with healthy ones toward a favorable or pathological outcome. This is largely due to technical limitations. We developed new genetic tools in Drosophila melanogaster that permit efficient multiplexed gain- and loss-of-function genetic perturbations with separable spatial and temporal control. Importantly, our novel tool-set is independent of the commonly used GAL4/UAS system, freeing the latter for additional, non-autonomous, genetic manipulations; and is built into a single strain, allowing one-generation interrogation of non-autonomous effects. Altogether, our design opens up efficient genome-wide screens on any deleterious phenotype, once plasmid or genome engineering is used to place the desired miRNA(s) or ORF(s) into our genotype. Specifically, we developed tools to study extrinsic effects on neural tumor growth but the strategy presented has endless applications within and beyond neurobiology, and in other model organisms.

Genome-wide patterns of local adaptation in Western European Drosophila melanogaster natural populations

Signatures of spatially varying selection have been investigated both at the genomic and transcriptomic level in several organisms. In Drosophila melanogaster, the majority of these studies have analyzed North American and Australian populations, leading to the identification of several loci and traits under selection. However, several studies based mainly in North American populations showed evidence of admixture that likely contributed to the observed population differentiation patterns. Thus, disentangling demography from selection might be challenging when analyzing these populations. European populations could help identify loci under spatially varying selection provided that no recent admixture from African populations would have occurred. In this work, we individually sequence the genome of 42 European strains collected in populations from contrasting environments: Stockholm (Sweden) and Castellana Grotte (Southern Italy). We found low levels of population structure and no evidence of recent African admixture in these two populations. We thus look for patterns of spatially varying selection affecting individual genes and gene sets. Besides single nucleotide polymorphisms, we also investigated the role of transposable elements in local adaptation. We concluded that European populations are a good dataset to identify candidate loci under spatially varying selection. The analysis of the two populations sequenced in this work in the context of all the available D. melanogaster data allowed us to pinpoint genes and biological processes likely to be relevant for local adaptation. Identifying and analyzing populations with low levels of population structure and admixture should help to disentangle selective from non-selective forces underlying patterns of population differentiation in other species as well.

Firefly genomes illuminate parallel origins of bioluminescence in beetles

Fireflies and their luminous courtships have inspired centuries of scientific study. Today firefly luciferase is widely used in biotechnology, but the evolutionary origin of bioluminescence within beetles remains unclear. To shed light on this long-standing question, we sequenced the genomes of two firefly species that diverged over 100 million-years-ago: the North American Photinus pyralis and Japanese Aquatica lateralis. To compare bioluminescent origins, we also sequenced the genome of a related click beetle, the Caribbean Ignelater luminosus, with bioluminescent biochemistry near-identical to fireflies, but anatomically unique light organs, suggesting the intriguing hypothesis of parallel gains of bioluminescence. Our analyses support independent gains of bioluminescence in fireflies and click beetles, and provide new insights into the genes, chemical defenses, and symbionts that evolved alongside their luminous lifestyle.

Using genomic data to study insecticide resistance in the house fly, Musca domestica

The house fly, Musca domestica, is a major pest at livestock facilities throughout the world. Insecticides have been the most common control strategy for flies, but many populations have evolved resistance. The speed by which we are able to identify the mutations responsible for resistance has been a major challenge for the development of high throughput resistance monitoring assays as new insecticides are introduced for control. This is particularly true for mutations that cause trans regulation of a gene, which then results in resistance. In this paper we take advantage of the conserved homology of dipteran chromosomes to assign 3069 genes to chromosomes. Of these, 234 were of toxicological interest (CYPs, esterases/hydrolases, glutathione S-transferases (GSTs) and target sites). The chromosomal location of genes known from linkage analysis studies matched the location predicted by homology mapping in ten out of ten cases, indicating a high reliability of our approach. The CYPs, esterases/hydrolases and GSTs were not randomly distributed throughout the genome. They clustered on chromosomes, but the pattern was different between the CYPs, esterases/hydrolases and GSTs. Examples are provided for how the availability of the house fly genome, combined with an ability to assign genes to chromosomes, will help to accelerate research in house flies.

The origin of the odorant receptor gene family in insects

The origin of the insect odorant receptor (OR) gene family has been hypothesized to have coincided with the evolution of terrestriality in insects. Missbach et al. (2014) suggested that ORs instead evolved with an ancestral OR co-receptor (Orco) after the origin of terrestriality and the OR/Orco system is an adaptation to winged flight in insects. We investigated genomes of the Collembola, Diplura, Archaeognatha, Zygentoma, Odonata, and Ephemeroptera, and find ORs present in all insect genomes but absent from lineages predating the evolution of insects. Orco is absent only in the ancestrally wingless insect lineage Archaeognatha. Our new genome sequence of the zygentoman firebrat Thermobia domestica reveals a full OR/Orco system. We conclude that ORs evolved before winged flight, perhaps as an adaptation to terrestriality, representing a key evolutionary novelty in the ancestor of all insects, and hence a molecular synapomorphy for the Class Insecta.

Ancestral polymorphisms explain the role of chromosomal inversions in speciation

Understanding the role of chromosomal inversions in speciation is a fundamental problem in evolutionary genetics. Here, we perform a comprehensive reconstruction of the evolutionary histories of the chromosomal inversions in Drosophila persimilis and D. pseudoobscura. We provide a solution to the puzzling origins of the selfish Sex-Ratio arrangement in D. persimilis and uncover surprising patterns of phylogenetic discordance on this chromosome. These patterns show that, contrary to widely held views, all fixed chromosomal inversions between D. persimilis and D. pseudoobscura were already present in their ancestral population long before the species split. Our results suggest that patterns of higher genomic divergence and an association of reproductive isolation genes with chromosomal inversions may be a direct consequence of incomplete lineage sorting of ancestral polymorphisms. These findings force a reconsideration of the role of chromosomal inversions in speciation, not as protectors of existing hybrid incompatibilities, but as fertile grounds for their formation.

The Complete Mitochondrial Genome of Ugyops sp. (Hemiptera: Delphacidae)

The complete mitochondrial genome (mitogenome) of Ugyops sp. (Hemiptera: Delphacidae) was sequenced, making it the first determined mitogenome from the subfamily Asiracinae, the basal clade of the family Delphacidae. The mitogenome was 15,259 bp in length with A + T content of 77.65% and contained 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a control region. The gene order was identical with that of the ancestral insect. The nucleotide composition analysis indicated that the whole mitogenome was strongly A-skewed (0.288) and highly C-skewed (-0.270). For PCGs on the J-strand, the AT skew was positive, and the GC skew was negative. All PCGs started with canonical ATN codons, except for cox1 and nad5, which used CTG and GTG as start codon, respectively. All tRNAs could fold into typical cloverleaf secondary structures, with the exception of trnS1 (AGN), in which the dihydrouridine arm was reduced to a simple loop. The control region included a poly-T stretch downstream of the small rRNA gene (rrnS), a subregion of higher A + T content and tandemly repeated sequence near trnI. The mitogenome of Ugyops sp. could be very helpful in exploring the diversity and evolution of mitogenomes in Delphacidae.

Chromatin immunoprecipitation (ChIP) method for non-model fruit flies (Diptera: Tephritidae) and evidence of histone modifications

Interactions between DNA and proteins located in the cell nucleus play an important role in controlling physiological processes by specifying, augmenting and regulating context-specific transcription events. Chromatin immunoprecipitation (ChIP) is a widely used methodology to study DNA-protein interactions and has been successfully used in various cell types for over three decades. More recently, by combining ChIP with genomic screening technologies and Next Generation Sequencing (e.g. ChIP-seq), it has become possible to profile DNA-protein interactions (including covalent histone modifications) across entire genomes. However, the applicability of ChIP-chip and ChIP-seq has rarely been extended to non-model species because of a number of technical challenges. Here we report a method that can be used to identify genome wide covalent histone modifications in a group of non-model fruit fly species (Diptera: Tephritidae). The method was developed by testing and refining protocols that have been used in model organisms, including Drosophila melanogaster. We demonstrate that this method is suitable for a group of economically important pest fruit fly species, viz., Bactrocera dorsalis, Ceratitis capitata, Zeugodacus cucurbitae and Bactrocera tryoni. We also report an example ChIP-seq dataset for B. tryoni, providing evidence for histone modifications in the genome of a tephritid fruit fly for the first time. Since tephritids are major agricultural pests globally, this methodology will be a valuable resource to study taxa-specific evolutionary questions and to assist with pest management. It also provides a basis for researchers working with other non-model species to undertake genome wide DNA-protein interaction studies.

Measuring Sister Chromatid Cohesion Protein Genome Occupancy in Drosophila melanogaster by ChIP-seq

This chapter presents methods to conduct and analyze genome-wide chromatin immunoprecipitation of the cohesin complex and the Nipped-B cohesin loading factor in Drosophila cells using high-throughput DNA sequencing (ChIP-seq). Procedures for isolation of chromatin, immunoprecipitation, and construction of sequencing libraries for the Ion Torrent Proton high throughput sequencer are detailed, and computational methods to calculate occupancy as input-normalized fold-enrichment are described. The results obtained by ChIP-seq are compared to those obtained by ChIP-chip (genomic ChIP using tiling microarrays), and the effects of sequencing depth on the accuracy are analyzed. ChIP-seq provides similar sensitivity and reproducibility as ChIP-chip, and identifies the same broad regions of occupancy. The locations of enrichment peaks, however, can differ between ChIP-chip and ChIP-seq, and low sequencing depth can splinter broad regions of occupancy into distinct peaks.

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.

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.

Comparative Gene Mapping as a Tool to Understand the Evolution of Pest Crop Insect Chromosomes

The extent of the conservation of synteny and gene order in aphids has been previously investigated only by comparing a small subset of linkage groups between the pea aphid Acyrthosiphon pisum and a few other aphid species. Here we compared the localization of eight A. pisum scaffolds (covering more than 5 Mb and 83 genes) in respect to the Drosophila melanogaster Muller elements identifying orthologous loci spanning all the four A. pisum chromosomes. Comparison of the genetic maps revealed a conserved synteny across different loci suggesting that the study of the fruit fly Muller elements could favour the identification of chromosomal markers useful for the study of chromosomal rearrangements in aphids. A. pisum is the first aphid species to have its genome sequenced and the finding that there are several chromosomal regions in synteny between Diptera and Hemiptera indicates that the genomic tools developed in A. pisum will be broadly useful not only for the study of other aphids but also for other insect species.

Patterns of Genome-Wide Diversity and Population Structure in the Drosophila athabasca Species Complex

The Drosophila athabasca species complex contains three recently diverged, prezygotically isolated semispecies (Western-Northern, Eastern-A, and Eastern-B) that are distributed across North America and share zones of sympatry. Inferences based on a handful of loci suggest that this complex might be an ideal system for studying the genetics of incipient speciation and the evolution of prezygotic isolating mechanisms, but patterns of differentiation have not been characterized systematically. Here, we assembled a draft genome for D. athabasca and analyze whole-genome re-sequencing data for 28 individuals from across the species range to characterize genome-wide patterns of diversity and population differentiation among semispecies. Patterns of differentiation on the X-chromosome vs. autosomes vary, with the X-chromosome showing better phylogenetic resolution and increased levels of between semispecies divergence. Despite low levels of overall differentiation and a lack of phylogenetic resolution of the autosomes for the most closely related semispecies, individuals do exhibit distinct genetic clustering. Demographic analyses provide some support for a model of isolation with migration within D. athabasca, with divergence times <20 kya. The young divergence times of the semispecies of D. athabasca, together with strong levels of sexual isolation, makes them a promising system for studying the evolution of prezygotic isolation and speciation.

Genome Size in North American Fireflies: Substantial Variation Likely Driven by Neutral Processes

Eukaryotic genomes show tremendous size variation across taxa. Proximate explanations for genome size variation include differences in ploidy and amounts of noncoding DNA, especially repetitive DNA. Ultimate explanations include selection on physiological correlates of genome size such as cell size, which in turn influence body size, resulting in the often-observed correlation between body size and genome size. In this study, we examined body size and repetitive DNA elements in relationship to the evolution of genome size in North American representatives of a single beetle family, the Lampyridae (fireflies). The 23 species considered represent an excellent study system because of the greater than 5-fold range of genome sizes, documented here using flow cytometry, and the 3-fold range in body size, measured using pronotum width. We also identified common genomic repetitive elements using low-coverage sequencing. We found a positive relationship between genome size and repetitive DNA, particularly retrotransposons. Both genome size and these elements were evolving as expected given phylogenetic relatedness. We also tested whether genome size varied with body size and found no relationship. Together, our results suggest that genome size is evolving neutrally in fireflies.

Genome-wide identification of lncRNAs associated with chlorantraniliprole resistance in diamondback moth Plutella xylostella (L.)

BACKGROUND

Long noncoding RNAs (lncRNAs) are now considered important regulatory factors, with a variety of biological functions in many species including insects. Some lncRNAs have the ability to show rapid responses to diverse stimuli or stress factors and are involved in responses to insecticide. However, there are no reports to date on the characterization of lncRNAs associated with chlorantraniliprole resistance in Plutella xylostella.

RESULTS

Nine RNA libraries constructed from one susceptible (CHS) and two chlorantraniliprole-resistant P. xylostella strains (CHR, ZZ) were sequenced, and 1309 lncRNAs were identified, including 877 intergenic lncRNAs, 190 intronic lncRNAs, 76 anti-sense lncRNAs and 166 sense-overlapping lncRNAs. Of the identified lncRNAs, 1059 were novel. Furthermore, we found that 64 lncRNAs were differentially expressed between CHR and CHS and 83 were differentially expressed between ZZ and CHS, of which 22 were differentially expressed in both CHR and ZZ. Most of the differentially expressed lncRNAs were hypothesized to be associated with chlorantraniliprole resistance in P. xylostella. The targets of lncRNAs via cis- (<10 kb upstream and downstream) or trans- (Pearson's correlation, r > 0.9 or < -0.9, P < 0.05) regulatory effects were also identified; many of the differently expressed lncRNAs were correlated with various important protein-coding genes involved in insecticide resistance, such as the ryanodine receptor, uridine diphosphate glucuronosyltransferase (UGTs), cytochrome P450, esterase and the ATP-binding cassette transporter.

CONCLUSIONS

Our results represent the first global identification of lncRNAs associated with chlorantraniliprole resistance in P. xylostella. These results will facilitate future studies of the regulatory mechanisms of lncRNAs in chlorantraniliprole and other insecticide resistance and in other biological processes in P. xylostella.

Genome sequencing of the winged midge, Parochlus steinenii, from the Antarctic Peninsula

Background

In the Antarctic, only two species of Chironomidae occur naturally-the wingless midge, Belgica antarctica , and the winged midge, Parochlus steinenii . B. antarctica is an extremophile with unusual adaptations. The larvae of B. antarctica are desiccation- and freeze-tolerant and the adults are wingless. Recently, the compact genome of B. antarctica was reported and it is the first Antarctic eukaryote to be sequenced. Although P. steinenii occurs naturally in the Antarctic with B. antarctica , the larvae of P. steinenii are cold-tolerant but not freeze-tolerant and the adults are winged. Differences in adaptations in the Antarctic midges are interesting in terms of evolutionary processes within an extreme environment. Herein, we provide the genome of another Antarctic midge to help elucidate the evolution of these species.

Results

The draft genome of P. steinenii had a total size of 138 Mbp, comprising 9513 contigs with an N50 contig size of 34,110 bp, and a GC content of 32.2%. Overall, 13,468 genes were predicted using the MAKER annotation pipeline, and gene ontology classified 10,801 (80.2%) predicted genes to a function. Compared with the assembled genome architecture of B. antarctica , that of P. steinenii was approximately 50 Mbp longer with 6.2-fold more repeat sequences, whereas gene regions were as similarly compact as in B. antarctica .

Conclusions

We present an annotated draft genome of the Antarctic midge, P. steinenii . The genomes of P. steinenii and B. antarctica will aid in the elucidation of evolution in harsh environments and provide new resources for functional genomic analyses of the order Diptera.

Genomics of the Asian rice gall midge and its interactions with rice

Understanding virulence and manipulative strategies of gall formers will reveal new facets of plant defense and insect counter defense. Among the gall midges, the Asian rice gall midge (AGM) has emerged as a model for studies on plant-insect interactions. Data from several genomics, transcriptomics and metabolomics studies have revealed diverse strategies adopted by AGM to successfully invade the host while overcoming its defense. Adaptive skills of AGM transcend from its genomic and transcriptomic make-up. Information arising from studies on genetics, mitochondrial genome and miRNAs, amongst other parameters, highlights AGM's capacity to maneuver the host defense, reorient host metabolome and redirect its morphogenesis.

Genomics of interaction between the brown planthopper and rice

Rice (Oryza sativa L.) and the brown planthopper (Nilaparvata lugens (Stål)) form a model system for dissection of the mechanism of interaction between insect pest and crop. In this review, we focus on the genomics of BPH-rice interaction. On the side of rice, a number of BPH-resistance genes have been identified genetically. Thirteen of these genes have been cloned which shed a light on the molecular basis of the interaction. On the aspect of BPH, a lot of salivary proteins have been identified using transcriptome and proteome techniques. The genetic loci of virulence were mapped in BPH genome based on the linkage map. The understanding of interaction between BPH and rice will provide novel insights into efficient control of this pest.

Improved annotation of the insect vector of citrus greening disease: biocuration by a diverse genomics community

Database URL

https://citrusgreening.org/.

Genome of the Asian longhorned beetle (Anoplophora glabripennis), a globally significant invasive species, reveals key functional and evolutionary innovations at the beetle-plant interface

BACKGROUND

Relatively little is known about the genomic basis and evolution of wood-feeding in beetles. We undertook genome sequencing and annotation, gene expression assays, studies of plant cell wall degrading enzymes, and other functional and comparative studies of the Asian longhorned beetle, Anoplophora glabripennis, a globally significant invasive species capable of inflicting severe feeding damage on many important tree species. Complementary studies of genes encoding enzymes involved in digestion of woody plant tissues or detoxification of plant allelochemicals were undertaken with the genomes of 14 additional insects, including the newly sequenced emerald ash borer and bull-headed dung beetle.

RESULTS

The Asian longhorned beetle genome encodes a uniquely diverse arsenal of enzymes that can degrade the main polysaccharide networks in plant cell walls, detoxify plant allelochemicals, and otherwise facilitate feeding on woody plants. It has the metabolic plasticity needed to feed on diverse plant species, contributing to its highly invasive nature. Large expansions of chemosensory genes involved in the reception of pheromones and plant kairomones are consistent with the complexity of chemical cues it uses to find host plants and mates.

CONCLUSIONS

Amplification and functional divergence of genes associated with specialized feeding on plants, including genes originally obtained via horizontal gene transfer from fungi and bacteria, contributed to the addition, expansion, and enhancement of the metabolic repertoire of the Asian longhorned beetle, certain other phytophagous beetles, and to a lesser degree, other phytophagous insects. Our results thus begin to establish a genomic basis for the evolutionary success of beetles on plants.

A Mitochondrial Genome of Rhyparochromidae (Hemiptera: Heteroptera) and a Comparative Analysis of Related Mitochondrial Genomes

The Rhyparochromidae, the largest family of Lygaeoidea, encompasses more than 1,850 described species, but no mitochondrial genome has been sequenced to date. Here we describe the first mitochondrial genome for Rhyparochromidae: a complete mitochondrial genome of Panaorus albomaculatus (Scott, 1874). This mitochondrial genome is comprised of 16,345 bp, and contains the expected 37 genes and control region. The majority of the control region is made up of a large tandem-repeat region, which has a novel pattern not previously observed in other insects. The tandem-repeats region of P. albomaculatus consists of 53 tandem duplications (including one partial repeat), which is the largest number of tandem repeats among all the known insect mitochondrial genomes. Slipped-strand mispairing during replication is likely to have generated this novel pattern of tandem repeats. Comparative analysis of tRNA gene families in sequenced Pentatomomorpha and Lygaeoidea species shows that the pattern of nucleotide conservation is markedly higher on the J-strand. Phylogenetic reconstruction based on mitochondrial genomes suggests that Rhyparochromidae is not the sister group to all the remaining Lygaeoidea, and supports the monophyly of Lygaeoidea.

Next-Generation Sequencing of Two Mitochondrial Genomes from Family Pompilidae (Hymenoptera: Vespoidea) Reveal Novel Patterns of Gene Arrangement

Animal mitochondrial genomes have provided large and diverse datasets for evolutionary studies. Here, the first two representative mitochondrial genomes from the family Pompilidae (Hymenoptera: Vespoidea) were determined using next-generation sequencing. The sequenced region of these two mitochondrial genomes from the species Auplopus sp. and Agenioideus sp. was 16,746 bp long with an A + T content of 83.12% and 16,596 bp long with an A + T content of 78.64%, respectively. In both species, all of the 37 typical mitochondrial genes were determined. The secondary structure of tRNA genes and rRNA genes were predicted and compared with those of other insects. Atypical trnS1 using abnormal anticodons TCT and lacking D-stem pairings was identified. There were 49 helices belonging to six domains in rrnL and 30 helices belonging to three domains in rrns present. Compared with the ancestral organization, four and two tRNA genes were rearranged in mitochondrial genomes of Auplopus and Agenioideus, respectively. In both species, trnM was shuffled upstream of the trnI-trnQ-trnM cluster, and trnA was translocated from the cluster trnA-trnR-trnN-trnS1-trnE-trnF to the region between nad1 and trnL1, which is novel to the Vespoidea. In Auplopus, the tRNA cluster trnW-trnC-trnY was shuffled to trnW-trnY-trnC. Phylogenetic analysis within Vespoidea revealed that Pompilidae and Mutillidae formed a sister lineage, and then sistered Formicidae. The genomes presented in this study have enriched the knowledge base of molecular markers, which is valuable in respect to studies about the gene rearrangement mechanism, genomic evolutionary processes and phylogeny of Hymenoptera.