Correction to: Genome-enabled insights into the biology of thrips as crop pests.

An amendment to this paper has been published and can be accessed via the original article.

Genome-enabled insights into the biology of thrips as crop pests.

BACKGROUND: The western flower thrips, Frankliniella occidentalis (Pergande), is a globally invasive pest and plant virus vector on a wide array of food, fiber, and ornamental crops. The underlying genetic mechanisms of the processes governing thrips pest and vector biology, feeding behaviors, ecology, and insecticide resistance are largely unknown. To address this gap, we present the F. occidentalis draft genome assembly and official gene set. RESULTS: We report on the first genome sequence for any member of the insect order Thysanoptera. Benchmarking Universal Single-Copy Ortholog (BUSCO) assessments of the genome assembly (size = 415.8 Mb, scaffold N50 = 948.9 kb) revealed a relatively complete and well-annotated assembly in comparison to other insect genomes. The genome is unusually GC-rich (50%) compared to other insect genomes to date. The official gene set (OGS v1.0) contains 16,859 genes, of which ~ 10% were manually verified and corrected by our consortium. We focused on manual annotation, phylogenetic, and expression evidence analyses for gene sets centered on primary themes in the life histories and activities of plant-colonizing insects. Highlights include the following: (1) divergent clades and large expansions in genes associated with environmental sensing (chemosensory receptors) and detoxification (CYP4, CYP6, and CCE enzymes) of substances encountered in agricultural environments; (2) a comprehensive set of salivary gland genes supported by enriched expression; (3) apparent absence of members of the IMD innate immune defense pathway; and (4) developmental- and sex-specific expression analyses of genes associated with progression from larvae to adulthood through neometaboly, a distinct form of maturation differing from either incomplete or complete metamorphosis in the Insecta. CONCLUSIONS: Analysis of the F. occidentalis genome offers insights into the polyphagous behavior of this insect pest that finds, colonizes, and survives on a widely diverse array of plants. The genomic resources presented here enable a more complete analysis of insect evolution and biology, providing a missing taxon for contemporary insect genomics-based analyses. Our study also offers a genomic benchmark for molecular and evolutionary investigations of other Thysanoptera species.

Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial Hymenoptera.

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO2 detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

Brown marmorated stink bug, Halyomorpha halys (Stål), genome: putative underpinnings of polyphagy, insecticide resistance potential and biology of a top worldwide pest.

BACKGROUND: Halyomorpha halys (Stål), the brown marmorated stink bug, is a highly invasive insect species due in part to its exceptionally high levels of polyphagy. This species is also a nuisance due to overwintering in human-made structures. It has caused significant agricultural losses in recent years along the Atlantic seaboard of North America and in continental Europe. Genomic resources will assist with determining the molecular basis for this species' feeding and habitat traits, defining potential targets for pest management strategies. RESULTS: Analysis of the 1.15-Gb draft genome assembly has identified a wide variety of genetic elements underpinning the biological characteristics of this formidable pest species, encompassing the roles of sensory functions, digestion, immunity, detoxification and development, all of which likely support H. halys' capacity for invasiveness. Many of the genes identified herein have potential for biomolecular pesticide applications. CONCLUSIONS: Availability of the H. halys genome sequence will be useful for the development of environmentally friendly biomolecular pesticides to be applied in concert with more traditional, synthetic chemical-based controls.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome.

BACKGROUND: The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae. RESULTS: The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding. CONCLUSIONS: With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

The genome of the water strider Gerris buenoi reveals expansions of gene repertoires associated with adaptations to life on the water.

BACKGROUND: Having conquered water surfaces worldwide, the semi-aquatic bugs occupy ponds, streams, lakes, mangroves, and even open oceans. The diversity of this group has inspired a range of scientific studies from ecology and evolution to developmental genetics and hydrodynamics of fluid locomotion. However, the lack of a representative water strider genome hinders our ability to more thoroughly investigate the molecular mechanisms underlying the processes of adaptation and diversification within this group. RESULTS: Here we report the sequencing and manual annotation of the Gerris buenoi (G. buenoi) genome; the first water strider genome to be sequenced thus far. The size of the G. buenoi genome is approximately 1,000 Mb, and this sequencing effort has recovered 20,949 predicted protein-coding genes. Manual annotation uncovered a number of local (tandem and proximal) gene duplications and expansions of gene families known for their importance in a variety of processes associated with morphological and physiological adaptations to a water surface lifestyle. These expansions may affect key processes associated with growth, vision, desiccation resistance, detoxification, olfaction and epigenetic regulation. Strikingly, the G. buenoi genome contains three insulin receptors, suggesting key changes in the rewiring and function of the insulin pathway. Other genomic changes affecting with opsin genes may be associated with wavelength sensitivity shifts in opsins, which is likely to be key in facilitating specific adaptations in vision for diverse water habitats. CONCLUSIONS: Our findings suggest that local gene duplications might have played an important role during the evolution of water striders. Along with these findings, the sequencing of the G. buenoi genome now provides us the opportunity to pursue exciting research opportunities to further understand the genomic underpinnings of traits associated with the extreme body plan and life history of water striders.

A model species for agricultural pest genomics: the genome of the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae).

The Colorado potato beetle is one of the most challenging agricultural pests to manage. It has shown a spectacular ability to adapt to a variety of solanaceaeous plants and variable climates during its global invasion, and, notably, to rapidly evolve insecticide resistance. To examine evidence of rapid evolutionary change, and to understand the genetic basis of herbivory and insecticide resistance, we tested for structural and functional genomic changes relative to other arthropod species using genome sequencing, transcriptomics, and community annotation. Two factors that might facilitate rapid evolutionary change include transposable elements, which comprise at least 17% of the genome and are rapidly evolving compared to other Coleoptera, and high levels of nucleotide diversity in rapidly growing pest populations. Adaptations to plant feeding are evident in gene expansions and differential expression of digestive enzymes in gut tissues, as well as expansions of gustatory receptors for bitter tasting. Surprisingly, the suite of genes involved in insecticide resistance is similar to other beetles. Finally, duplications in the RNAi pathway might explain why Leptinotarsa decemlineata has high sensitivity to dsRNA. The L. decemlineata genome provides opportunities to investigate a broad range of phenotypes and to develop sustainable methods to control this widely successful pest.

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.

The whole genome sequence of the Mediterranean fruit fly, Ceratitis capitata (Wiedemann), reveals insights into the biology and adaptive evolution of a highly invasive pest species.

BACKGROUND: The Mediterranean fruit fly (medfly), Ceratitis capitata, is a major destructive insect pest due to its broad host range, which includes hundreds of fruits and vegetables. It exhibits a unique ability to invade and adapt to ecological niches throughout tropical and subtropical regions of the world, though medfly infestations have been prevented and controlled by the sterile insect technique (SIT) as part of integrated pest management programs (IPMs). The genetic analysis and manipulation of medfly has been subject to intensive study in an effort to improve SIT efficacy and other aspects of IPM control. RESULTS: The 479 Mb medfly genome is sequenced from adult flies from lines inbred for 20 generations. A high-quality assembly is achieved having a contig N50 of 45.7 kb and scaffold N50 of 4.06 Mb. In-depth curation of more than 1800 messenger RNAs shows specific gene expansions that can be related to invasiveness and host adaptation, including gene families for chemoreception, toxin and insecticide metabolism, cuticle proteins, opsins, and aquaporins. We identify genes relevant to IPM control, including those required to improve SIT. CONCLUSIONS: The medfly genome sequence provides critical insights into the biology of one of the most serious and widespread agricultural pests. This knowledge should significantly advance the means of controlling the size and invasive potential of medfly populations. Its close relationship to Drosophila, and other insect species important to agriculture and human health, will further comparative functional and structural studies of insect genomes that should broaden our understanding of gene family evolution.

Methuselah/Methuselah-like G protein-coupled receptors constitute an ancient metazoan gene family.

Inconsistent conclusions have been drawn regarding the phylogenetic age of the Methuselah/Methuselah-like (Mth/Mthl) gene family of G protein-coupled receptors, the founding member of which regulates development and lifespan in Drosophila. Here we report the results from a targeted homolog search of 39 holozoan genomes and phylogenetic analysis of the conserved seven transmembrane domain. Our findings reveal that the Mth/Mthl gene family is ancient, has experienced numerous extinction and expansion events during metazoan evolution, and acquired the current definition of the Methuselah ectodomain during its exceptional expansion in arthropods. In addition, our findings identify Mthl1, Mthl5, Mthl14, and Mthl15 as the oldest Mth/Mthl gene family paralogs in Drosophila. Future studies of these genes have the potential to define ancestral functions of the Mth/Mthl gene family.

Recruitment of Saccharomyces cerevisiae Cmr1/Ydl156w to Coding Regions Promotes Transcription Genome Wide.

Cmr1 (changed mutation rate 1) is a largely uncharacterized nuclear protein that has recently emerged in several global genetic interaction and protein localization studies. It clusters with proteins involved in DNA damage and replication stress response, suggesting a role in maintaining genome integrity. Under conditions of proteasome inhibition or replication stress, this protein localizes to distinct sub-nuclear foci termed as intranuclear quality control (INQ) compartments, which sequester proteins for their subsequent degradation. Interestingly, it also interacts with histones, chromatin remodelers and modifiers, as well as with proteins involved in transcription including subunits of RNA Pol I and Pol III, but not with those of Pol II. It is not known whether Cmr1 plays a role in regulating transcription of Pol II target genes. Here, we show that Cmr1 is recruited to the coding regions of transcribed genes of S. cerevisiae. Cmr1 occupancy correlates with the Pol II occupancy genome-wide, indicating that it is recruited to coding sequences in a transcription-dependent manner. Cmr1-enriched genes include Gcn4 targets and ribosomal protein genes. Furthermore, our results show that Cmr1 recruitment to coding sequences is stimulated by Pol II CTD kinase, Kin28, and the histone deacetylases, Rpd3 and Hos2. Finally, our genome-wide analyses implicate Cmr1 in regulating Pol II occupancy at transcribed coding sequences. However, it is dispensable for maintaining co-transcriptional histone occupancy and histone modification (acetylation and methylation). Collectively, our results show that Cmr1 facilitates transcription by directly engaging with transcribed coding regions.

Gene Ages, Nomenclatures, and Functional Diversification of the Methuselah/Methuselah-Like GPCR Family in Drosophila and Tribolium.

Affecting lifespan regulation and oxidative stress resistance, the G-protein coupled receptor (GPCR) gene methuselah (mth) plays important roles in the life history of Drosophila melanogaster. Substantial progress has been made in elucidating the molecular pathways by which mth affects these traits, yet conflicting ideas exist as to how old these genetic interactions are as well as how old the mth gene itself is. Root to these issues is the complex gene family history of the Mth/Mthl GPCR family, which experienced independent expansions in a variety of animal clades, leading to at least six subfamilies in insects. Within insects, drosophilid flies stand out by possessing up to three times more Mth/Mthl receptors due to the expansion of a single subfamily, the mth superclade subfamily, which contains an even younger subfamily introduced here as the melanogaster subgroup subfamily. As a result, most of the 16 Mth/Mthl receptors of D. melanogaster are characterized by n:1 orthology relationships to singleton mth superclade homologs in nondrosophilid species. This challenge is exacerbated by the inconsistent naming of Mth/Mthl orthologs across species. To consolidate this situation, we review established ortholog relationships among insect Mth/Mthl receptors, clarify the gene nomenclatures in two important satellite model species, the fruit fly relative D. virilis and the red flour beetle Tribolium castaneum, and discuss the genetic and functional evolution of the D. melanogaster Mth GPCR.

Histone deacetylases and phosphorylated polymerase II C-terminal domain recruit Spt6 for cotranscriptional histone reassembly.

Spt6 is a multifunctional histone chaperone involved in the maintenance of chromatin structure during elongation by RNA polymerase II (Pol II). Spt6 has a tandem SH2 (tSH2) domain within its C terminus that recognizes Pol II C-terminal domain (CTD) peptides phosphorylated on Ser2, Ser5, or Try1 in vitro. Deleting the tSH2 domain, however, only has a partial effect on Spt6 occupancy in vivo, suggesting that more complex mechanisms are involved in the Spt6 recruitment. Our results show that the Ser2 kinases Bur1 and Ctk1, but not the Ser5 kinase Kin28, cooperate in recruiting Spt6, genome-wide. Interestingly, the Ser2 kinases promote the association of Spt6 in early transcribed regions and not toward the 3' ends of genes, where phosphorylated Ser2 reaches its maximum level. In addition, our results uncover an unexpected role for histone deacetylases (Rpd3 and Hos2) in promoting Spt6 interaction with elongating Pol II. Finally, our data suggest that phosphorylation of the Pol II CTD on Tyr1 promotes the association of Spt6 with the 3' ends of transcribed genes, independently of Ser2 phosphorylation. Collectively, our results show that a complex network of interactions, involving the Spt6 tSH2 domain, CTD phosphorylation, and histone deacetylases, coordinate the recruitment of Spt6 to transcribed genes in vivo.

Dramatic expansion and developmental expression diversification of the methuselah gene family during recent Drosophila evolution.

Functional studies of the methuselah/methuselah-like (mth/mthl) gene family have focused on the founding member mth, but little is known regarding the developmental functions of this receptor or any of its paralogs. We undertook a comprehensive analysis of developmental expression and sequence divergence in the mth/mthl gene family. Using in situ hybridization techniques, we detect expression of six genes (mthl1, 5, 9, 11, 13, and 14) in the embryo during gastrulation and development of the gut, heart, and lymph glands. Four receptors (mthl3, 4, 6, and 8) are expressed in the larval central nervous system, imaginal discs, or both, and two receptors (mthl10 and mth) are expressed in both embryos and larvae. Phylogenetic analysis of all mth/mthl genes in five Drosophila species, mosquito and flour beetle structured the mth/mthl family into several subclades. mthl1, 5, and 14 are present in most species, each forming a separate clade. A newly identified Drosophila mthl gene (CG31720; herein mthl15) formed another ancient clade. The remaining Drosophila receptors, including mth, are members of a large "superclade" that diversified relatively recently during dipteran evolution, in many cases within the melanogaster subgroup. Comparing the expression patterns of the mth/mthl "superclade" paralogs to the embryonic expression of the singleton ortholog in Tribolium suggests both subfunctionalization and acquisition of novel functionalities. Taken together, our findings shed novel light on mth as a young member of an adaptively evolving developmental gene family.