FAIR Header Reference genome: a TRUSTworthy standard.

The lack of interoperable data standards among reference genome data-sharing platforms inhibits cross-platform analysis while increasing the risk of data provenance loss. Here, we describe the FAIR bioHeaders Reference genome (FHR), a metadata standard guided by the principles of Findability, Accessibility, Interoperability and Reuse (FAIR) in addition to the principles of Transparency, Responsibility, User focus, Sustainability and Technology. The objective of FHR is to provide an extensive set of data serialisation methods and minimum data field requirements while still maintaining extensibility, flexibility and expressivity in an increasingly decentralised genomic data ecosystem. The effort needed to implement FHR is low; FHR's design philosophy ensures easy implementation while retaining the benefits gained from recording both machine and human-readable provenance.

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.

Interspecific gene flow obscures phylogenetic relationships in an important insect pest species complex.

As genomic data proliferates, the prevalence of post-speciation gene flow is making species boundaries and relationships increasingly ambiguous. Although current approaches inferring fully bifurcating phylogenies based on concatenated datasets provide simple and robust answers to many species relationships, they may be inaccurate because the models ignore inter-specific gene flow and incomplete lineage sorting. To examine the potential error resulting from ignoring gene flow, we generated both a RAD-seq and a 500 protein-coding loci highly multiplexed amplicon (HiMAP) dataset for a monophyletic group of 12 species defined as the Bactrocera dorsalis sensu lato clade. With some of the world's worst agricultural pests, the taxonomy of the B. dorsalis s.l. clade is important for trade and quarantines. However, taxonomic confusion confounds resolution due to intra- and interspecific phenotypic variation and convergence, mitochondrial introgression across half of the species, and viable hybrids. We compared the topological convergence of our datasets using concatenated phylogenetic and various multispecies coalescent approaches, some of which account for gene flow. All analyses agreed on species delimitation, but there was incongruence between species relationships. Under concatenation, both datasets suggest identical species relationships with mostly high statistical support. However, multispecies coalescent and multispecies network approaches suggest markedly different hypotheses and detected significant gene flow. We suggest that the network approaches are likely more accurate because gene flow violates the assumptions of the concatenated phylogenetic analyses, but the data-reductive requirements of network approaches resulted in reduced statistical support and could not unambiguously resolve gene flow directions. Our study highlights the importance of testing for gene flow, particularly with phylogenomic datasets, even when concatenated approaches receive high statistical support.

Whole-genome resequencing data support a single introduction of the invasive white pine sawfly, Diprion similis.

Biological introductions are unintended "natural experiments" that provide unique insights into evolutionary processes. Invasive phytophagous insects are of particular interest to evolutionary biologists studying adaptation, as introductions often require rapid adaptation to novel host plants. However, adaptive potential of invasive populations may be limited by reduced genetic diversity-a problem known as the "genetic paradox of invasions." One potential solution to this paradox is if there are multiple invasive waves that bolster genetic variation in invasive populations. Evaluating this hypothesis requires characterizing genetic variation and population structure in the invaded range. To this end, we assemble a reference genome and describe patterns of genetic variation in the introduced white pine sawfly, Diprion similis. This species was introduced to North America in 1914, where it has rapidly colonized the thin-needled eastern white pine (Pinus strobus), making it an ideal invasion system for studying adaptation to novel environments. To evaluate evidence of multiple introductions, we generated whole-genome resequencing data for 64 D. similis females sampled across the North American range. Both model-based and model-free clustering analyses supported a single population for North American D. similis. Within this population, we found evidence of isolation-by-distance and a pattern of declining heterozygosity with distance from the hypothesized introduction site. Together, these results support a single-introduction event. We consider implications of these findings for the genetic paradox of invasion and discuss priorities for future research in D. similis, a promising model system for invasion biology.

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.

Incongruence between molecules and morphology: A seven-gene phylogeny of Dacini fruit flies paves the way for reclassification (Diptera: Tephritidae).

Molecular and morphological research often suggest conflicting results. Selective pressure on certain morphologies can confound understanding of evolutionary relationships. Dacini is one of the most diverse tribes of tephritid flies and contains many economically important pest species. Their black and yellow patterned body markings are presumed to act as wasp mimicry, and the characters separating species and groups are limited and in some cases phenotypically plastic. The traditional taxonomy of the tribe is controversial because groupings are based on unique combinations of morphological characters without the use of cladistic methods, though recent phylogenetic and taxonomic analyses have resulted in significant changes to their taxonomy. The monophyly of the three largest genera in the tribe has been tested with only small numbers of representatives per genus and a limited number of genes. To further understand the taxonomy and evolution of Dacini we sequenced seven genes from 167 Dacini species and five dipteran outgroups to construct a robust phylogeny and test phylogenetic relationships between genera, subgenera, and species complexes. Our phylogeny confirms the monophyly of Dacus, Bactrocera, and Zeugodacus. However, most groups below the genus level are not monophyletic, and only through further revision will we be able to understand their evolution and clarify the taxonomy within this tribe.

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.