Genome and Ontogenetic-Based Transcriptomic Analyses of the Flesh Fly, Sarcophaga bullata.

The flesh fly, Sarcophaga bullata, is a widely-used model for examining the physiology of insect diapause, development, stress tolerance, neurobiology, and host-parasitoid interactions. Flies in this taxon are implicated in myiasis (larval infection of vertebrates) and feed on carrion, aspects that are important in forensic studies. Here we present the genome of S. bullata, along with developmental- and reproduction-based RNA-Seq analyses. We predict 15,768 protein coding genes, identify orthology in relation to closely related flies, and establish sex and developmental-specific gene sets based on our RNA-Seq analyses. Genomic sequences, predicted genes, and sequencing data sets have been deposited at the National Center for Biotechnology Information. Our results provide groundwork for genomic studies that will expand the flesh fly's utility as a model system.

Changes in microRNA abundance may regulate diapause in the flesh fly, Sarcophaga bullata.

Diapause, an alternative developmental pathway characterized by changes in developmental timing and metabolism, is coordinated by molecular mechanisms that are not completely understood. MicroRNA (miRNA) mediated gene silencing is emerging as a key component of animal development and may have a significant role in initiating, maintaining, and terminating insect diapause. In the present study, we test this possibility by using high-throughput sequencing and qRT-PCR to discover diapause-related shifts in miRNA abundance in the flesh fly, Sarcophaga bullata. We identified ten evolutionarily conserved miRNAs that were differentially expressed in diapausing pupae compared to their nondiapausing counterparts. miR-289-5p and miR-1-3p were overexpressed in diapausing pupae and may be responsible for silencing expression of candidate genes during diapause. miR-9c-5p, miR-13b-3p, miR-31a-5p, miR-92b-3p, miR-275-3p, miR-276a-3p, miR-277-3p, and miR-305-5p were underexpressed in diapausing pupae and may contribute to increased expression of heat shock proteins and other factors necessary for the enhanced environmental stress-response that is a feature of diapause. In S. bullata, a maternal effect blocks the programming of diapause in progeny of females that have experienced pupal diapause, and in this study we report that several miRNAs, including miR-263a-5p, miR-100-5p, miR-125-5p, and let-7-5p were significantly overexpressed in such nondiapausing flies and may prevent entry into diapause. Together these miRNAs appear to be integral to the molecular processes that mediate entry into diapause.

Continuous activity and no cycling of clock genes in the Antarctic midge during the polar summer.

The extreme seasonal shifts of day length in polar regions, ranging from constant light in the summer to constant darkness in the winter, pose an intriguing environment for probing activity rhythms and the functioning of circadian clocks. Here, we monitor locomotor activity during the summer on the Antarctic Peninsula and under laboratory conditions, as well as the accompanying patterns of clock gene expression in the Antarctic midge, the only insect endemic to Antarctica. Larvae and adults are most active during the warmest portion of the day, but at a constant temperature they remain continuously active regardless of the photoregime, and activity also persists in constant darkness. The canonical clock genes period, timeless, Clock, and vrille are expressed in the head but we detected no cycling of expression in either the field or under diverse photoregimes in the laboratory. The timekeeping function of the clock has possibly been lost, enabling the midge to opportunistically exploit the unpredictable availability of permissive thermal conditions for growth, development, and reproduction during the short summer in Antarctica.

Compact genome of the Antarctic midge is likely an adaptation to an extreme environment.

The midge, Belgica antarctica, is the only insect endemic to Antarctica, and thus it offers a powerful model for probing responses to extreme temperatures, freeze tolerance, dehydration, osmotic stress, ultraviolet radiation and other forms of environmental stress. Here we present the first genome assembly of an extremophile, the first dipteran in the family Chironomidae, and the first Antarctic eukaryote to be sequenced. At 99 megabases, B. antarctica has the smallest insect genome sequenced thus far. Although it has a similar number of genes as other Diptera, the midge genome has very low repeat density and a reduction in intron length. Environmental extremes appear to constrain genome architecture, not gene content. The few transposable elements present are mainly ancient, inactive retroelements. An abundance of genes associated with development, regulation of metabolism and responses to external stimuli may reflect adaptations for surviving in this harsh environment.

Gene expression changes governing extreme dehydration tolerance in an Antarctic insect.

Among terrestrial organisms, arthropods are especially susceptible to dehydration, given their small body size and high surface area to volume ratio. This challenge is particularly acute for polar arthropods that face near-constant desiccating conditions, as water is frozen and thus unavailable for much of the year. The molecular mechanisms that govern extreme dehydration tolerance in insects remain largely undefined. In this study, we used RNA sequencing to quantify transcriptional mechanisms of extreme dehydration tolerance in the Antarctic midge, Belgica antarctica, the world's southernmost insect and only insect endemic to Antarctica. Larvae of B. antarctica are remarkably tolerant of dehydration, surviving losses up to 70% of their body water. Gene expression changes in response to dehydration indicated up-regulation of cellular recycling pathways including the ubiquitin-mediated proteasome and autophagy, with concurrent down-regulation of genes involved in general metabolism and ATP production. Metabolomics results revealed shifts in metabolite pools that correlated closely with changes in gene expression, indicating that coordinated changes in gene expression and metabolism are a critical component of the dehydration response. Finally, using comparative genomics, we compared our gene expression results with a transcriptomic dataset for the Arctic collembolan, Megaphorura arctica. Although B. antarctica and M. arctica are adapted to similar environments, our analysis indicated very little overlap in expression profiles between these two arthropods. Whereas several orthologous genes showed similar expression patterns, transcriptional changes were largely species specific, indicating these polar arthropods have developed distinct transcriptional mechanisms to cope with similar desiccating conditions.

Combined transcriptomic and metabolomic approach uncovers molecular mechanisms of cold tolerance in a temperate flesh fly.

The ability to respond rapidly to changes in temperature is critical for insects and other ectotherms living in variable environments. In a physiological process termed rapid cold-hardening (RCH), exposure to nonlethal low temperature allows many insects to significantly increase their cold tolerance in a matter of minutes to hours. Additionally, there are rapid changes in gene expression and cell physiology during recovery from cold injury, and we hypothesize that RCH may modulate some of these processes during recovery. In this study, we used a combination of transcriptomics and metabolomics to examine the molecular mechanisms of RCH and cold shock recovery in the flesh fly, Sarcophaga bullata. Surprisingly, out of ∼15,000 expressed sequence tags (ESTs) measured, no transcripts were upregulated during RCH, and likewise RCH had a minimal effect on the transcript signature during recovery from cold shock. However, during recovery from cold shock, we observed differential expression of ∼1,400 ESTs, including a number of heat shock proteins, cytoskeletal components, and genes from several cell signaling pathways. In the metabolome, RCH had a slight yet significant effect on several metabolic pathways, while cold shock resulted in dramatic increases in gluconeogenesis, amino acid synthesis, and cryoprotective polyol synthesis. Several biochemical pathways showed congruence at both the transcript and metabolite levels, indicating that coordinated changes in gene expression and metabolism contribute to recovery from cold shock. Thus, while RCH had very minor effects on gene expression, recovery from cold shock elicits sweeping changes in gene expression and metabolism along numerous cell signaling and biochemical pathways.

Evolutionary convergence to ideal free dispersal strategies and coexistence.

We study a two species competition model in which the species have the same population dynamics but different dispersal strategies and show how these dispersal strategies evolve. We introduce a general dispersal strategy which can result in the ideal free distributions of both competing species at equilibrium and generalize the result of Averill et al. (2011). We further investigate the convergent stability of this ideal free dispersal strategy by varying random dispersal rates, advection rates, or both of these two parameters simultaneously. For monotone resource functions, our analysis reveals that among two similar dispersal strategies, selection generally prefers the strategy which is closer to the ideal free dispersal strategy. For nonmonotone resource functions, our findings suggest that there may exist some dispersal strategies which are not ideal free, but could be locally evolutionarily stable and/or convergent stable, and allow for the coexistence of more than one species.

Heat shock response to hypoxia and its attenuation during recovery in the flesh fly, Sarcophaga crassipalpis.

In this study, pharate adults of the flesh fly Sarcophaga crassipalpis were exposed to two, four, seven, or ten days of severe hypoxia (3% oxygen) to evaluate its impact on emergence and the expression of genes encoding heat shock proteins (Hsps) and heat shock regulatory elements. A four-day exposure to hypoxia significantly reduced survival, but more than seven days was required to reach the LD(50). Eight genes encoding Hsps, at least one from each major family of Hsps (Hsp90, Hsp70, Hsp60, Hsp40, and sHsps) and two genes encoding proteins involved in Hsp regulation (heat shock factor, hsf, and sirtuin) were cloned, and expression levels were assessed during and after hypoxia using qRT-PCR. Most, but not all hsps studied, were significantly up-regulated during hypoxia, and expression levels for most of the hsps reverted to control levels a few hours after return to normoxia. Hsp70 was the most responsive to hypoxia, increasing expression several hundred fold. By contrast, hsp90 and hsp27 showed little response to hypoxia but did respond to recovery. Neither hsf nor sirtuin were elevated by hypoxia, an observation consistent with their assumed post-transcriptional regulatory roles. These data demonstrate a strong Hsp response to hypoxia, suggesting an important role for Hsps in responding to low oxygen environments.