Whole Genomes Reveal Evolutionary Relationships and Mechanisms Underlying Gene-Tree Discordance in Neodiprion Sawflies.

Rapidly evolving taxa are excellent models for understanding the mechanisms that give rise to biodiversity. However, developing an accurate historical framework for comparative analysis of such lineages remains a challenge due to ubiquitous incomplete lineage sorting and introgression. Here, we use a whole-genome alignment, multiple locus-sampling strategies, and summary-tree and SNP-based species-tree methods to infer a species tree for eastern North American Neodiprion species, a clade of pine-feeding sawflies (Order: Hymenopteran; Family: Diprionidae). We recovered a well-supported species tree that-except for three uncertain relationships-was robust to different strategies for analyzing whole-genome data. Nevertheless, underlying gene-tree discordance was high. To understand this genealogical variation, we used multiple linear regression to model site concordance factors estimated in 50-kb windows as a function of several genomic predictor variables. We found that site concordance factors tended to be higher in regions of the genome with more parsimony-informative sites, fewer singletons, less missing data, lower GC content, more genes, lower recombination rates, and lower D-statistics (less introgression). Together, these results suggest that incomplete lineage sorting, introgression, and genotyping error all shape the genomic landscape of gene-tree discordance in Neodiprion. More generally, our findings demonstrate how combining phylogenomic analysis with knowledge of local genomic features can reveal mechanisms that produce topological heterogeneity across genomes.

Reference genome for the Mojave poppy bee (Perdita meconis), a specialist pollinator of conservation concern.

The Mojave poppy bee, Perdita meconis Griswold (Hymenoptera: Anthophila: Andrenidae), is a species of conservation concern that is restricted to the eastern Mojave Desert of North America. It is a specialist pollinator of two poppy genera, Arctomecon and Argemone (Papaveraceae), and is being considered for listing under the US Endangered Species Act along with one of its pollinator hosts, the Las Vegas bearpoppy (Arctomecon californica). Here, we present a near chromosome-level genome of the Mojave poppy bee to provide a genomic resource that will aid conservation efforts and future research. We isolated DNA from a single, small (<7 mm), male specimen collected using non-ideal preservation methods then performed whole-genome sequencing using PacBio HiFi technology. After quality and contaminant filtering, the final draft genome assembly is 327 Mb, with an N50 length of 17.5 Mb. Annotated repetitive elements compose 37.3% of the genome, although a large proportion (24.87%) of those are unclassified repeats. Additionally, we annotated 18,245 protein-coding genes and 19,433 transcripts. This genome represents one of only a few genomes from the large bee family Andrenidae and one of only a few genomes for pollinator specialists. We highlight both the potential of this genome as a resource for future research, and how high-quality genomes generated from small, non-ideal (in terms of preservation) specimens could facilitate biodiversity genomics.

HiFiAdapterFilt, a memory efficient read processing pipeline, prevents occurrence of adapter sequence in PacBio HiFi reads and their negative impacts on genome assembly.

BACKGROUND: Pacific Biosciences HiFi read technology is currently the industry standard for high accuracy long-read sequencing that has been widely adopted by large sequencing and assembly initiatives for generation of de novo assemblies in non-model organisms. Though adapter contamination filtering is routine in traditional short-read analysis pipelines, it has not been widely adopted for HiFi workflows. RESULTS: Analysis of 55 publicly available HiFi datasets revealed that a read-sanitation step to remove sequence artifacts derived from PacBio library preparation from read pools is necessary as adapter sequences can be erroneously integrated into assemblies. CONCLUSIONS: Here we describe the nature of adapter contaminated reads, their consequences in assembly, and present HiFiAdapterFilt, a simple and memory efficient solution for removing adapter contaminated reads prior to assembly.

Rapid Viral Symbiogenesis via Changes in Parasitoid Wasp Genome Architecture.

Viral genome integration provides a complex route to biological innovation that has rarely but repeatedly occurred in one of the most diverse lineages of organisms on the planet, parasitoid wasps. We describe a novel endogenous virus in braconid wasps derived from pathogenic alphanudiviruses. Limited to a subset of the genus Fopius, this recent acquisition allows an unprecedented opportunity to examine early endogenization events. Massive amounts of virus-like particles (VLPs) are produced in wasp ovaries. Unlike most endogenous viruses of parasitoid wasps, the VLPs do not contain DNA, translating to major differences in parasitism-promoting strategies. Rapid changes include genomic rearrangement, loss of DNA processing proteins, and wasp control of viral gene expression. These events precede the full development of tissue-specific viral gene expression observed in older associations. These data indicate that viral endogenization can rapidly result in functional and evolutionary changes associated with genomic novelty and adaptation in parasitoids.

Wheat streak mosaic virus alters the transcriptome of its vector, wheat curl mite (Aceria tosichella Keifer), to enhance mite development and population expansion.

Wheat streak mosaic virus (WSMV; genus Tritimovirus; family Potyviridae) is an economically important wheat virus that is transmitted by the wheat curl mite (WCM; Aceria tosichella Keifer) in a persistent manner. Virus-vector coevolution may potentially influence vector gene expression to prolong viral association and thus increase virus transmission efficiency and spread. To understand the transcriptomic responses of WCM to WSMV, RNA sequencing was performed to assemble and analyse transcriptomes of WSMV viruliferous and aviruliferous mites. Among 7291 de novo-assembled unigenes, 1020 were differentially expressed between viruliferous and aviruliferous WCMs using edgeR at a false discovery rate ≤0.05. Differentially expressed unigenes were enriched for 108 gene ontology terms, with the majority of the unigenes showing downregulation in viruliferous mites in comparison to only a few unigenes that were upregulated. Protein family and metabolic pathway enrichment analyses revealed that most downregulated unigenes encoded enzymes and proteins linked to stress response, immunity and development. Mechanistically, these predicted changes in mite physiology induced by viral association could be suggestive of pathways needed for promoting virus-vector interactions. Overall, our data suggest that transcriptional changes in viruliferous mites facilitate prolonged viral association and alter WCM development to expedite population expansion, both of which could enhance viral transmission.

Phylogenomics supports incongruence between ecological specialization and taxonomy in a charismatic clade of buck moths.

Local adaptation can be a fundamental component of speciation, but its dynamics in relation to gene flow are not necessarily straightforward. Herbivorous taxa with localized host plant or habitat specialization across their geographic range are ideal models for investigating the patterns and constraints of local adaptation and its impact on diversification. The charismatic, day-flying moths of the Hemileuca maia species complex (Lepidoptera: Saturniidae) are such taxa, as they are geographically widespread, exhibit considerable ecological and morphological variability and host and habitat specificity, but apparently lack genetic differentiation across their range. Here, we use genomewide single nucleotide polymorphisms to assess relationships and population structure of this group across North America and investigate the scales where genomic divergence correlates with adaptive ecological characteristics. In contrast to previous genetic studies of the group, we find broad- and fine-scale genetic differentiation between lineages, which is at odds with various levels of taxonomic description and recognition of conservation units. Furthermore, ecological specialization only explains some fine-scale genetic differentiation, and across much of the group's range, local adaptation is apparently occurring in the face of strong gene flow. These results provide unprecedented insight into drivers of speciation in this group, the relationship between taxonomy and genomics-informed species boundaries and conservation management of internationally protected entities. Broadly, this system provides a model for understanding how local adaptation in an herbivore can arise and be maintained in the face of apparently strong gene flow, and the importance of geographic isolation in generating genomic divergence, despite a lack of ecological divergence.

Transpacific coalescent pathways of coconut rhinoceros beetle biotypes: Resistance to biological control catalyses resurgence of an old pest.

Biological control agents have several advantages over chemical control for pest management, including the capability to restore ecosystem balance with minimal non-target effects and a lower propensity for targets to develop resistance. These factors are particularly important for invasive species control. The coconut rhinoceros beetle (Oryctes rhinoceros Linnaeus) is a major palm pest that invaded many Pacific islands in the early 20th century through human-mediated dispersal. Application of the Oryctes nudivirus in the 1960s successfully halted the beetle's first invasion wave and made it a textbook example of successful biological control. However, a recently discovered O. rhinoceros biotype that is resistant to the nudivirus appears to be correlated with a new invasion wave. We performed a population genomics analysis of 172 O. rhinoceros from seven regions, including native and invasive populations, to reconstruct invasion pathways and explore correlation between recent invasions and biotypes. With ddRAD sequencing, we generated data sets ranging from 4,000 to 209,000 loci using stacks and ipyrad software pipelines and compared genetic signal in downstream clustering and phylogenetic analyses. Analysis suggests that the O. rhinoceros resurgence is mediated by the nudivirus-resistant biotype. Genomic data have been proven essential to understanding the new O. rhinoceros biotype's invasion patterns and interactions with the original biotype. Such information is crucial to optimization of strategies for quarantine and control of resurgent pests. Our results demonstrate that while invasions are relatively rare events, new introductions can have significant ecological consequences, and quarantine vigilance is required even in previously invaded areas.

Comparative Proteomic Profiling between Each of Two Consecutive Developmental Stages of the Solanum Fruit Fly, Bactrocera latifrons (Hendel).

The Solanum fruit fly, Bactrocera latifrons (Hendel), has a complex life cycle including multiple stages (egg, larva, pupa, and adult). Understanding the details of "what", "when", "where", "why", and "how" many hundred thousand proteins operate in this insect, interact, and express between each two consecutive developmental stages at molecular level not only can expand our knowledge, but also lead to the development of novel fruit fly control techniques. We tried to find what, when, and where in this study. Why and how will be presented in upcoming papers. We conducted a proteome profiling using 2-D gel electrophoresis and mass spectrometry. Samples of 3-day-old eggs, 1- and 10-day-old larvae, 1- and 10-day-old pupae, 1- and 9-day-old females and males of B. latifrons were used. A custom peptide database, derived from the de novo B. latifrons whole genome assembly was used for peptide identification. Differentially expressed proteins (DEPs) with significant fold expression and protein functions between two consecutive developmental stages were identified, annotated, described, and listed in gel images and/or charts. With this foundational information, we are not only providing valuable information, but also any impacts due to the biotic or abiotic environmental factors can be identified and manipulated, and lead to further research on gene editing and biomarker discovery.

De novo metatranscriptome assembly and coral gene expression profile of Montipora capitata with growth anomaly.

BACKGROUND: Scleractinian corals are a vital component of coral reef ecosystems, and of significant cultural and economic value worldwide. As anthropogenic and natural stressors are contributing to a global decline of coral reefs, understanding coral health is critical to help preserve these ecosystems. Growth anomaly (GA) is a coral disease that has significant negative impacts on coral biology, yet our understanding of its etiology and pathology is lacking. In this study we used RNA-seq along with de novo metatranscriptome assembly and homology assignment to identify coral genes that are expressed in three distinct coral tissue types: tissue from healthy corals ("healthy"), GA lesion tissue from diseased corals ("GA-affected") and apparently healthy tissue from diseased corals ("GA-unaffected"). We conducted pairwise comparisons of gene expression among these three tissue types to identify genes and pathways that help us to unravel the molecular pathology of this coral disease. RESULTS: The quality-filtered de novo-assembled metatranscriptome contained 76,063 genes, of which 13,643 were identified as putative coral genes. Overall gene expression profiles of coral genes revealed high similarity between healthy tissue samples, in contrast to high variance among diseased samples. This indicates GA has a variety of genetic effects at the colony level, including on seemingly healthy (GA-unaffected) tissue. A total of 105 unique coral genes were found differentially expressed among tissue types. Pairwise comparisons revealed the greatest number of differentially expressed genes between healthy and GA-affected tissue (93 genes), followed by healthy and GA-unaffected tissue (33 genes), and GA-affected and -unaffected tissue (7 genes). The putative function of these genes suggests GA is associated with changes in the activity of genes involved in developmental processes and activation of the immune system. CONCLUSION: This is one of the first transcriptome-level studies to investigate coral GA, and the first metatranscriptome assembly for the M. capitata holobiont. The gene expression data, metatranscriptome assembly and methodology developed through this study represent a significant addition to the molecular information available to further our understanding of this coral disease.

Prediction of a peptidome for the western tarnished plant bug Lygus hesperus.

Many strategies for controlling insect pests require an understanding of their hormonal signaling agents, peptides being the largest and most diverse single class of these molecules. Lygus hesperus is a pest species of particular concern, as it is responsible for significant damage to a wide variety of commercially important plant crops. At present, little is known about the peptide hormones of L. hesperus. Here, transcriptomic data were used to predict a peptidome for L. hesperus. Fifty-three L. hesperus transcripts encoding peptide precursors were identified, with a subset amplified by PCR for sequence verification. The proteins deduced from these transcripts allowed for the prediction of a 119-sequence peptidome for L. hesperus. The predicted peptides include isoforms of allatostatin A, allatostatin B (AST-B), allatostatin C, allatotropin, bursicon, CCHamide, corazonin, crustacean cardioactive peptide, crustacean hyperglycemic hormone/ion transport peptide, diuretic hormone 31, GSEFLamide, insulin-like peptide, myosuppressin, neuroparsin, neuropeptide F, orcokinin, orcomyotropin, pyrokinin, short neuropeptide F, SIFamide, sulfakinin and tachykinin-related peptide. Of note were several isoforms of AST-B that possess -WX7Wamide carboxyl-termini rather than the stereotypical -WX6Wamide (e.g., KWQDMQNPGWamide), an allatotropin ending in -SARGFamide rather than -TARGFamide (GLKNGPLNSARGFamide), a GSEFLamide ending in -GTEFLamide (TVGTEFLamide), several orcokinins with PMDEIDR- rather than NFDEIDR- amino-termini (e.g., PMDEIDRAGFTHFV), and an eight rather than 12 amino acid long isoform of SIFamide (PPFNGSIFamide). Collectively, the L. hesperus peptidome predicted here provides a resource for initiating physiological investigations of peptidergic signaling in this species, including studies directed at the biological control of this agricultural pest.

Comprehensive annotation of Glossina pallidipes salivary gland hypertrophy virus from Ethiopian tsetse flies: a proteogenomics approach.

Glossina pallidipes salivary gland hypertrophy virus (GpSGHV; family Hytrosaviridae) can establish asymptomatic and symptomatic infection in its tsetse fly host. Here, we present a comprehensive annotation of the genome of an Ethiopian GpSGHV isolate (GpSGHV-Eth) compared with the reference Ugandan GpSGHV isolate (GpSGHV-Uga; GenBank accession number EF568108). GpSGHV-Eth has higher salivary gland hypertrophy syndrome prevalence than GpSGHV-Uga. We show that the GpSGHV-Eth genome has 190 291 nt, a low G+C content (27.9 %) and encodes 174 putative ORFs. Using proteogenomic and transcriptome mapping, 141 and 86 ORFs were mapped by transcripts and peptides, respectively. Furthermore, of the 174 ORFs, 132 had putative transcriptional signals [TATA-like box and poly(A) signals]. Sixty ORFs had both TATA-like box promoter and poly(A) signals, and mapped by both transcripts and peptides, implying that these ORFs encode functional proteins. Of the 60 ORFs, 10 ORFs are homologues to baculovirus and nudivirus core genes, including three per os infectivity factors and four RNA polymerase subunits (LEF4, 5, 8 and 9). Whereas GpSGHV-Eth and GpSGHV-Uga are 98.1 % similar at the nucleotide level, 37 ORFs in the GpSGHV-Eth genome had nucleotide insertions (n = 17) and deletions (n = 20) compared with their homologues in GpSGHV-Uga. Furthermore, compared with the GpSGHV-Uga genome, 11 and 24 GpSGHV ORFs were deleted and novel, respectively. Further, 13 GpSGHV-Eth ORFs were non-canonical; they had either CTG or TTG start codons instead of ATG. Taken together, these data suggest that GpSGHV-Eth and GpSGHV-Uga represent two different lineages of the same virus. Genetic differences combined with host and environmental factors possibly explain the differential GpSGHV pathogenesis observed in different G. pallidipes colonies.

LARVAL X-RAY IRRADIATION INFLUENCES PROTEIN EXPRESSION IN PUPAE OF THE ORIENTAL FRUIT FLY, BACTROCERA DORSALIS.

The sterile insect technique (SIT) was developed to eradicate the new world screwworm from the southern United States and Mexico, and became a component of many area-wide integrated pest management programs, particularly useful in managing tephritid fruit flies. SIT is based on the idea of rearing and sterilizing male pests, originally by ionizing radiation, and then releasing into field, where they compete for and mate with wild females. Mating with sterile males leads to reduced fecundity to lower pest populations. There are concerns with the use and distribution of radioisotopes for SIT programs, which have led to developing X-ray irradiation protocols to sterilize insects. We considered the possibility that X-ray irradiation exerts sublethal impacts aside form sterilizing insects. Such effects may not be directly observable, which led us to the hypothesis that X-ray irradiation in one life stage creates alterations in biological fitness and protein expression in the subsequent stage. We tested our hypothesis by irradiating larvae of Bactrocera dorsalis. There are two major points. One, exposing larvae to X-ray treatments led to reduced adult emergence, fecundity, fertility, and flight capacity from the corresponding pupae and emerged adults. Two, the X-ray treatments led to substantial expression changes in 27 pupal proteins. We assorted the 67 spots representing these proteins into three groups, metabolism, development, and structure. Our interpretation is these X-ray induced changes in biological performance and protein expression indicate their adult counterparts may be disabled in their abilities to successfully compete for and mate wild females in native habitats.

Identification of a carboxylesterase associated with resistance to naled in Bactrocera dorsalis (Hendel).

Compared to other organophosphate-resistant and -susceptible (S) lines of Bactrocera dorsalis, the carboxylesterase (CBE) BdE5 in the naled-resistant (nal-r) line has been found to possess remarkable quantitative elevation. Our study attempts to identify the role of BdE5 in naled resistance, and we discovered several points of interest. Firstly, activity staining on native PAGE revealed that the percentage of flies with intensive BdE5 bands in the nal-r line was substantially higher than in the S line, indicating that the BdE5 band correlates with naled susceptibility. Secondly, in vitro and in vivo inhibition assays showed that BdE5 was inhibited by naled in both lines; under diagnostic doses of naled, the overall extent of inhibition on CBEs was much greater in the S line than in the nal-r line. Thirdly, NanoLC-nanoESi-MS/MS analysis used the NCBI database to identify and annotate BdE5 as an esterase FE4-like (XP_011200445.1) in B. dorsalis. Fourthly, rapid amplification of cDNA ends was used to obtain the 2012-bp full-length BdE5 cDNA, which contained an open reading frame of 1770bp and encoded a putative protein of 590 amino acid residues. Phylogenetic analysis revealed that BdE5 is a secreted ß-esterase (E clade) closely related to CG6414 (NP_570089), a CBE in Drosophila melanogaster. Finally, our relative quantification real-time PCR data showed a significant elevation in transcript levels of the BdE5 gene in nal-r line. Our results confirmed that BdE5 is correlated with naled resistance and provides further understanding about the identification and molecular characteristics of BdE5 in B. dorsalis.

Molecular Markers Detect Cryptic Predation on Coffee Berry Borer (Coleoptera: Curculionidae) by Silvanid and Laemophloeid Flat Bark Beetles (Coleoptera: Silvanidae, Laemophloeidae) in Coffee Beans.

The coffee berry borer, Hypothenemus hampei (Ferrari) (Coleoptera: Curculionidae), is a serious pest of coffee worldwide. It was first detected in Hawai'i in 2010. Two predatory beetles, Cathartus quadricollis (Coleoptera: Silvanidae) and Leptophloeus sp. (Coleoptera: Laemophloeidae), have been observed in H. hampei-infested coffee. Under laboratory conditions, colony-reared C. quadricollis and Leptophloeus sp. prey upon all life stages of H. hampei. However, the H. hampei life cycle occurs almost exclusively within a coffee bean obscured from direct observation. Thus, it is unknown if C. quadricollis and Leptophloeus sp. consume H. hampei as prey in the wild. To demonstrate predation of H. hampei by C. quadricollis and Leptophloeus sp., a molecular assay was developed utilizing species-specific primers targeting short regions of the mitochondrial COI gene to determine species presence. Using these primers, wild C. quadricollis and Leptophloeus sp. were collected and screened for the presence of H. hampei DNA using PCR. Analysis of collections from five coffee farms revealed predation of C. quadricollis and Leptophloeus sp. on H. hampei. Further laboratory testing showed that H. hampei DNA could be detected in predators for as long as 48 h after feeding, indicating the farm-caught predators had preyed on H. hampei within 2 d of sampling. This study demonstrates the utility of molecular markers for the study of the ecology of predators and prey with cryptic behavior, and suggests C. quadricollis and Leptophloeus sp. might be useful biocontrol agents against H. hampei.

MicroRNAs in the oriental fruit fly, Bactrocera dorsalis: extending Drosophilid miRNA conservation to the Tephritidae.

BACKGROUND: The oriental fruit fly, Bactrocera dorsalis, is an important plant pest species in the family Tephritidae. It is a phytophagous species with broad host range, and while not established in the mainland United States, is a species of great concern for introduction. Despite the vast amount of information available from the closely related model organism Drosophila melanogaster, information at the genome and transcriptome level is still very limited for this species. Small RNAs act as regulatory molecules capable of determining transcript levels in the cells. The most studied small RNAs are micro RNAs, which may impact as much as 30 % of all protein coding genes in animals. RESULTS: We have sequenced small RNAs (sRNAs) from the Tephritid fruit fly, B. dorsalis (oriental fruit fly), specifically sRNAs corresponding to the 17 to 28 nucleotides long fraction of total RNA. Sequencing yielded more than 16 million reads in total. Seventy five miRNAs orthologous to known miRNAs were identified, as well as five additional novel miRNAs that might be specific to the genera, or to the Tephritid family. We constructed a gene expression profile for the identified miRNAs, and used comparative analysis with D. melanogaster to support our expression data. In addition, several miRNA clusters were identified in the genome that show conservancy with D. melanogaster. Potential targets for the identified miRNAs were also searched. CONCLUSIONS: The data presented here adds to our growing pool of information concerning the genome structure and characteristics of true fruit flies. It provides a basis for comparative studies with other Dipteran and within Tephritid species, and can be used for applied research such as in the development of new control strategies based on gene silencing and transgenesis.

A genomic perspective to assessing quality of mass-reared SIT flies used in Mediterranean fruit fly (Ceratitis capitata) eradication in California.

BACKGROUND: Temperature sensitive lethal (tsl) mutants of the tephritid C. capitata are used extensively in control programs involving sterile insect technique in California. These flies are artificially reared and treated with ionizing radiation to render males sterile for further release en masse into the field to compete with wild males and disrupt establishment of invasive populations. Recent research suggests establishment of C. capitata in California, despite the fact that over 250 million sterile flies are released weekly as part of the state's preventative program. In this project, genome-level quality assessment was performed, measured as expression differences between the Vienna-7 tsl mutants used in SIT programs and wild flies. RNA-seq was performed to provide a genome-wide map of the messenger RNA populations in C. capitata, and to investigate significant expression changes in Vienna-7 mass reared flies. RESULTS: Flies from the Vienna-7 colony showed a markedly reduced abundance of transcripts related to visual and chemical responses, including light stimuli, neural development and signaling pathways when compared to wild flies. In addition, genes associated with muscle development and locomotion were shown to be reduced. This suggests that the Vienna-7 line may be less competitive in mating and host plant finding where these stimuli are utilized. Irradiated flies showed several transcripts representing stress associated with irradiation. CONCLUSIONS: There are significant changes at the transcriptome level that likely alter the competitiveness of mass reared flies and provide justification for pursuing methods for strain improvement, increasing competitiveness of mass-reared flies, or exploring alternative SIT approaches to increase the efficiency of eradication programs.

Characterizing the developmental transcriptome of the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae) through comparative genomic analysis with Drosophila melanogaster utilizing modENCODE datasets.

BACKGROUND: The oriental fruit fly, Bactrocera dorsalis, is an important pest of fruit and vegetable crops throughout Asia, and is considered a high risk pest for establishment in the mainland United States. It is a member of the family Tephritidae, which are the most agriculturally important family of flies, and can be considered an out-group to well-studied members of the family Drosophilidae. Despite their importance as pests and their relatedness to Drosophila, little information is present on B. dorsalis transcripts and proteins. The objective of this paper is to comprehensively characterize the transcripts present throughout the life history of B. dorsalis and functionally annotate and analyse these transcripts relative to the presence, expression, and function of orthologous sequences present in Drosophila melanogaster. RESULTS: We present a detailed transcriptome assembly of B. dorsalis from egg through adult stages containing 20,666 transcripts across 10,799 unigene components. Utilizing data available through Flybase and the modENCODE project, we compared expression patterns of these transcripts to putative orthologs in D. melanogaster in terms of timing, abundance, and function. In addition, temporal expression patterns in B. dorsalis were characterized between stages, to establish the constitutive or stage-specific expression patterns of particular transcripts. A fully annotated transcriptome assembly is made available through NCBI, in addition to corresponding expression data. CONCLUSIONS: Through characterizing the transcriptome of B. dorsalis through its life history and comparing the transcriptome of B. dorsalis to the model organism D. melanogaster, a database has been developed that can be used as the foundation to functional genomic research in Bactrocera flies and help identify orthologous genes between B. dorsalis and D. melanogaster. This data provides the foundation for future functional genomic research that will focus on improving our understanding of the physiology and biology of this species at the molecular level. This knowledge can also be applied towards developing improved methods for control, survey, and eradication of this important pest.

Functional genomics and microbiome profiling of the Asian longhorned beetle (Anoplophora glabripennis) reveal insights into the digestive physiology and nutritional ecology of wood feeding beetles.

BACKGROUND: Wood-feeding beetles harbor an ecologically rich and taxonomically diverse assemblage of gut microbes that appear to promote survival in woody tissue, which is devoid of nitrogen and essential nutrients. Nevertheless, the contributions of these apparent symbionts to digestive physiology and nutritional ecology remain uncharacterized in most beetle lineages. RESULTS: Through parallel transcriptome profiling of beetle- and microbial- derived mRNAs, we demonstrate that the midgut microbiome of the Asian longhorned beetle (Anoplophora glabripennis), a member of the beetle family Cerambycidae, is enriched in biosynthetic pathways for the synthesis of essential amino acids, vitamins, and sterols. Consequently, the midgut microbiome of A. glabripennis can provide essential nutrients that the beetle cannot obtain from its woody diet or synthesize itself. The beetle gut microbiota also produce their own suite of transcripts that can enhance lignin degradation, degrade hemicellulose, and ferment xylose and wood sugars. An abundance of cellulases from several glycoside hydrolase families are expressed endogenously by A. glabripennis, as well as transcripts that allow the beetle to convert microbe-synthesized essential amino acids into non-essential amino acids. A. glabripennis and its gut microbes likely collaborate to digest carbohydrates and convert released sugars and amino acid intermediates into essential nutrients otherwise lacking from their woody host plants. CONCLUSIONS: The nutritional provisioning capabilities of the A. glabripennis gut microbiome may contribute to the beetles' unusually broad host range. The presence of some of the same microbes in the guts of other Cerambycidae and other wood-feeding beetles suggests that partnerships with microbes may be a facilitator of evolutionary radiations in beetles, as in certain other groups of insects, allowing access to novel food sources through enhanced nutritional provisioning.

Genetic Diversity of Bactrocera dorsalis (Diptera: Tephritidae) on the Hawaiian Islands: Implications for an Introduction Pathway Into California.

Population genetic diversity of the oriental fruit fly, Bactrocera dorsalis (Hendel), on the Hawaiian islands of Oahu, Maui, Kauai, and Hawaii (the Big Island) was estimated using DNA sequences of the mitochondrial cytochrome c oxidase subunit I gene. In total, 932 flies representing 36 sampled sites across the four islands were sequenced for a 1,500-bp fragment of the gene named the C1500 marker. Genetic variation was low on the Hawaiian Islands with >96% of flies having just two haplotypes: C1500-Haplotype 1 (63.2%) or C1500-Haplotype 2 (33.3%). The other 33 flies (3.5%) had haplotypes similar to the two dominant haplotypes. No population structure was detected among the islands or within islands. The two haplotypes were present at similar frequencies at each sample site, suggesting that flies on the various islands can be considered one population. Comparison of the Hawaiian data set to DNA sequences of 165 flies from outbreaks in California between 2006 and 2012 indicates that a single-source introduction pathway of Hawaiian origin cannot explain many of the flies in California. Hawaii, however, could not be excluded as a maternal source for 69 flies. There was no clear geographic association for Hawaiian or non-Hawaiian haplotypes in the Bay Area or Los Angeles Basin over time. This suggests that California experienced multiple, independent introductions from different sources.

Chromosome-scale genome assembly of the Hunt bumble bee, Bombus huntii Greene, 1860, a species of agricultural interest.

The Hunt bumble bee, Bombus huntii, is a widely distributed pollinator in western North America. The species produces large colony sizes in captive rearing conditions, experiences low parasite and pathogen loads, and has been demonstrated to be an effective pollinator of tomatoes grown in controlled environment agriculture systems. These desirable traits have galvanized producer efforts to develop commercial B. huntii colonies for growers to deliver pollination services to crops. To better understand B. huntii biology and support population genetic studies and breeding decisions, we sequenced and assembled the B. huntii genome from a single haploid male. High-fidelity sequencing of the entire genome using PacBio, along with HiC sequencing, led to a comprehensive contig assembly of high continuity. This assembly was further organized into a chromosomal arrangement, successfully identifying 18 chromosomes spread across the 317.4 Mb assembly with a BUSCO score indicating 97.6% completeness. Synteny analysis demonstrates shared chromosome number (n = 18) with B. terrestris, a species belonging to a different subgenus, matching the expectation that presence of 18 haploid chromosomes is an ancestral trait at least between the subgenera Pyrobombus and Bombus sensu stricto. In conclusion, the assembly outcome, alongside the minimal tissue sampled destructively, showcase efficient techniques for producing a comprehensive, highly contiguous genome.