Identification and Characterization of Small RNA Markers of Age in the Blow Fly Cochliomyia macellaria (Fabricius) (Diptera: Calliphoridae).

Blow fly development is important in decomposition ecology, agriculture, and forensics. Much of the impact of these species is from immature samples, thus knowledge of their development is important to enhance or ameliorate their effects. One application of this information is the estimation of immature insect age to provide temporal information for death investigations. While traditional markers of age such as stage and size are generally accurate, they lack precision in later developmental stages. We used miRNA sequencing to measure miRNA expression, throughout development, of the secondary screwworm, Cochliomyia macellaria (Fabricius) (Diptera: Calliphoridae) and identified 217 miRNAs present across the samples. Ten were identified to be significantly differentially expressed in larval samples and seventeen were found to be significantly differentially expressed in intrapuparial samples. Twenty-eight miRNAs were identified to be differentially expressed between sexes. Expression patterns of two miRNAs, miR-92b and bantam, were qPCR-validated in intrapuparial samples; these and likely food-derived miRNAs appear to be stable markers of age in C. macellaria. Our results support the use of miRNAs for developmental markers of age and suggest further investigations across species and under a range of abiotic and biotic conditions.

Evaluation of Development Datasets for Hermetia illucens (L.) (Diptera: Stratiomyidae) for Estimating the Time of Placement of Human and Swine Remains in Texas, USA.

A basic tenet of forensic entomology is development data of an insect can be used to predict the time of colonization (TOC) by insect specimens collected from remains, and this prediction is related to the time of death and/or time of placement (TOP). However, few datasets have been evaluated to determine their accuracy or precision. The black soldier fly, Hermetia illucens (L.) (Diptera: Stratiomyidae) is recognized as an insect of forensic importance. This study examined the accuracy and precision of several development datasets for the black soldier fly by estimating the TOP of five sets of human and three sets of swine remains in San Marcos and College Station, TX, respectively. Data generated from this study indicate only one of these datasets consistently (time-to-prepupae 52%; time-to-eclosion 75%) produced TOP estimations that occurred within a day of the actual TOP of the remains. It is unknown if the precolonization interval (PreCI) of this species is long, but it has been observed that the species can colonize within 6 d after death. This assumption remains untested by validation studies. Accounting for this PreCI improved accuracy for the time-to-prepupae group, but reduced accuracy in the time-to-eclosion group. The findings presented here highlight a need for detailed, forensic-based development data for the black soldier fly that can reliably and accurately be used in casework. Finally, this study outlines the need for a basic understanding of the timing of resource utilization (i.e., duration of the PreCI) for forensically relevant taxa so that reasonable corrections may be made to TOC as related to minimum postmortem interval (mPMI) estimates.

Biogeographic parallels in thermal tolerance and gene expression variation under temperature stress in a widespread bumble bee.

Global temperature changes have emphasized the need to understand how species adapt to thermal stress across their ranges. Genetic mechanisms may contribute to variation in thermal tolerance, providing evidence for how organisms adapt to local environments. We determine physiological thermal limits and characterize genome-wide transcriptional changes at these limits in bumble bees using laboratory-reared Bombus vosnesenskii workers. We analyze bees reared from latitudinal (35.7-45.7°N) and altitudinal (7-2154 m) extremes of the species' range to correlate thermal tolerance and gene expression among populations from different climates. We find that critical thermal minima (CTMIN) exhibit strong associations with local minimums at the location of queen origin, while critical thermal maximum (CTMAX) was invariant among populations. Concordant patterns are apparent in gene expression data, with regional differentiation following cold exposure, and expression shifts invariant among populations under high temperatures. Furthermore, we identify several modules of co-expressed genes that tightly correlate with critical thermal limits and temperature at the region of origin. Our results reveal that local adaptation in thermal limits and gene expression may facilitate cold tolerance across a species range, whereas high temperature responses are likely constrained, both of which may have implications for climate change responses of bumble bees.

The genomes of a monogenic fly: views of primitive sex chromosomes.

The production of male and female offspring is often determined by the presence of specific sex chromosomes which control sex-specific expression, and sex chromosomes evolve through reduced recombination and specialized gene content. Here we present the genomes of Chrysomya rufifacies, a monogenic blow fly (females produce female or male offspring, exclusively) by separately sequencing and assembling each type of female and the male. The genomes (> 25X coverage) do not appear to have any sex-linked Muller F elements (typical for many Diptera) and exhibit little differentiation between groups supporting the morphological assessments of C. rufifacies homomorphic chromosomes. Males in this species are associated with a unimodal coverage distribution while females exhibit bimodal coverage distributions, suggesting a potential difference in genomic architecture. The presence of the individual-sex draft genomes herein provides new clues regarding the origination and evolution of the diverse sex-determining mechanisms observed within Diptera. Additional genomic analysis of sex chromosomes and sex-determining genes of other blow flies will allow a refined evolutionary understanding of how flies with a typical X/Y heterogametic amphogeny (male and female offspring in similar ratios) sex determination systems evolved into one with a dominant factor that results in single sex progeny in a chromosomally monomorphic system.

De Novo Genome Assemblies for Three North American Bumble Bee Species: Bombus bifarius, Bombus vancouverensis, and Bombus vosnesenskii.

Bumble bees are ecologically and economically important insect pollinators. Three abundant and widespread species in western North America, Bombus bifarius, Bombus vancouverensis, and Bombus vosnesenskii, have been the focus of substantial research relating to diverse aspects of bumble bee ecology and evolutionary biology. We present de novo genome assemblies for each of the three species using hybrid assembly of Illumina and Oxford Nanopore Technologies sequences. All three assemblies are of high quality with large N50s (> 2.2 Mb), BUSCO scores indicating > 98% complete genes, and annotations producing 13,325 - 13,687 genes, comparing favorably with other bee genomes. Analysis of synteny against the most complete bumble bee genome, Bombus terrestris, reveals a high degree of collinearity. These genomes should provide a valuable resource for addressing questions relating to functional genomics and evolutionary biology in these species.

Local adaptation across a complex bioclimatic landscape in two montane bumble bee species.

Understanding evolutionary responses to variation in temperature and precipitation across species ranges is of fundamental interest given ongoing climate change. The importance of temperature and precipitation for multiple aspects of bumble bee (Bombus) biology, combined with large geographic ranges that expose populations to diverse environmental pressures, make these insects well-suited for studying local adaptation. Here, we analyzed genome-wide sequence data from two widespread bumble bees, Bombus vosnesenskii and Bombus vancouverensis, using multiple environmental association analysis methods to investigate climate adaptation across latitude and altitude. The strongest signatures of selection were observed in B. vancouverensis, but despite unique responses between species for most loci, we detected several shared responses. Genes relating to neural and neuromuscular function and ion transport were especially evident with respect to temperature variables, while genes relating to cuticle formation, tracheal and respiratory system development, and homeostasis were associated with precipitation variables. Our data thus suggest that adaptive responses for tolerating abiotic variation are likely to be complex, but that several parallels among species can emerge even for these complex traits and landscapes. Results provide the framework for future work into mechanisms of thermal and desiccation tolerance in bumble bees and a set of genomic targets that might be monitored for future conservation efforts.

Gene expression correlates of facultative predation in the blow fly Chrysomya rufifacies (Diptera: Calliphoridae).

Effects of intraguild predation (IGP) on omnivores and detritivores are relatively understudied when compared to work on predator guilds. Functional genetic work in IGP is even more limited, but its application can help answer a range of questions related to ultimate and proximate causes of this behavior. Here, we integrate behavioral assays and transcriptomic analysis of facultative predation in a blow fly (Diptera: Calliphoridae) to evaluate the prevalence, effect, and correlated gene expression of facultative predation by the invasive species Chrysomya rufifacies. Field work observing donated human cadavers indicated facultative predation by C. rufifacies on the native blow fly Cochliomyia macellaria was rare under undisturbed conditions, owing in part to spatial segregation between species. Laboratory assays under conditions of starvation showed predation had a direct fitness benefit (i.e., survival) to the predator. As a genome is not available for C. rufifacies, a de novo transcriptome was developed and annotated using sequence similarity to Drosophila melanogaster. Under a variety of assembly parameters, several genes were identified as being differentially expressed between predators and nonpredators of this species, including genes involved in cell-to-cell signaling, osmotic regulation, starvation responses, and dopamine regulation. Results of this work were integrated to develop a model of the processes and genetic regulation controlling facultative predation.

Distance, elevation and environment as drivers of diversity and divergence in bumble bees across latitude and altitude.

Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site-associated DNA sequencing (RADseq) in two bumble bee species, Bombus vosnesenskii and Bombus bifarius, across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A. Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure while B. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, with B. vosnesenskii exhibiting relatively consistent levels of genetic diversity across its range, while B. bifarius has reduced genetic diversity at low latitudes, where it is restricted to high-elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.

A scalable and memory-efficient algorithm for de novo transcriptome assembly of non-model organisms.

BACKGROUND: With increased availability of de novo assembly algorithms, it is feasible to study entire transcriptomes of non-model organisms. While algorithms are available that are specifically designed for performing transcriptome assembly from high-throughput sequencing data, they are very memory-intensive, limiting their applications to small data sets with few libraries. RESULTS: We develop a transcriptome assembly algorithm that recovers alternatively spliced isoforms and expression levels while utilizing as many RNA-Seq libraries as possible that contain hundreds of gigabases of data. New techniques are developed so that computations can be performed on a computing cluster with moderate amount of physical memory. CONCLUSIONS: Our strategy minimizes memory consumption while simultaneously obtaining comparable or improved accuracy over existing algorithms. It provides support for incremental updates of assemblies when new libraries become available.

Population genomics reveals a candidate gene involved in bumble bee pigmentation.

Variation in bumble bee color patterns is well-documented within and between species. Identifying the genetic mechanisms underlying such variation may be useful in revealing evolutionary forces shaping rapid phenotypic diversification. The widespread North American species Bombus bifarius exhibits regional variation in abdominal color forms, ranging from red-banded to black-banded phenotypes and including geographically and phenotypically intermediate forms. Identifying genomic regions linked to this variation has been complicated by strong, near species level, genome-wide differentiation between red- and black-banded forms. Here, we instead focus on the closely related black-banded and intermediate forms that both belong to the subspecies B. bifarius nearcticus. We analyze an RNA sequencing (RNAseq) data set and identify a cluster of single nucleotide polymorphisms (SNPs) within one gene, Xanthine dehydrogenase/oxidase-like, that exhibit highly unusual differentiation compared to the rest of the sequenced genome. Homologs of this gene contribute to pigmentation in other insects, and results thus represent a strong candidate for investigating the genetic basis of pigment variation in B. bifarius and other bumble bee mimicry complexes.

Structural and Genetic Investigation of the Egg and First-Instar Larva of an Egg-Laying Population of Blaesoxipha plinthopyga (Diptera: Sarcophagidae), a Species of Forensic Importance.

Flies in the family Sarcophagidae incubate their eggs and are known to be ovoviviparous (i.e., ovolarviparous), but a laboratory-maintained colony of Blaesoxipha plinthopyga (Wiedemann) deposited clutches of viable eggs over 10 generations. A description of the egg and first-instar larva of this species is provided along with genetic data (genome size and cytochrome oxidase I sequences). The egg is similar to previously described eggs of other Sarcophagidae but differs in the configuration of the micropyle. In the first-instar larva, the oral ridges are much more developed than has been described for other species. B. plinthopyga has forensic importance, and the present descriptive information is critical for proper case management.

Bacteria mediate oviposition by the black soldier fly, Hermetia illucens (L.), (Diptera: Stratiomyidae).

There can be substantial negative consequences for insects colonizing a resource in the presence of competitors. We hypothesized that bacteria, associated with an oviposition resource and the insect eggs deposited on that resource, serve as a mechanism regulating subsequent insect attraction, colonization, and potentially succession of insect species. We isolated and identified bacterial species associated with insects associated with vertebrate carrion and used these bacteria to measure their influence on the oviposition preference of adult black soldier flies which utilizes animal carcasses and is an important species in waste management and forensics. We also ascertained that utilizing a mixture of bacteria, rather than a single species, differentially influenced behavioral responses of the flies, as did bacterial concentration and the species of fly from which the bacteria originated. These studies provide insight into interkingdom interactions commonly occurring during decomposition, but not commonly studied.

Proteus mirabilis interkingdom swarming signals attract blow flies.

Flies transport specific bacteria with their larvae that provide a wider range of nutrients for those bacteria. Our hypothesis was that this symbiotic interaction may depend on interkingdom signaling. We obtained Proteus mirabilis from the salivary glands of the blow fly Lucilia sericata; this strain swarmed significantly and produced a strong odor that attracts blow flies. To identify the putative interkingdom signals for the bacterium and flies, we reasoned that as swarming is used by this bacterium to cover the food resource and requires bacterial signaling, the same bacterial signals used for swarming may be used to communicate with blow flies. Using transposon mutagenesis, we identified six novel genes for swarming (ureR, fis, hybG, zapB, fadE and PROSTU_03490), then, confirming our hypothesis, we discovered that fly attractants, lactic acid, phenol, NaOH, KOH and ammonia, restore swarming for cells with the swarming mutations. Hence, compounds produced by the bacterium that attract flies also are utilized for swarming. In addition, bacteria with the swarming mutation rfaL attracted fewer blow flies and reduced the number of eggs laid by the flies. Therefore, we have identified several interkingdom signals between P. mirabilis and blow flies.