Whether coding or regulatory sequence change is more important to the evolution of phenotypic novelty is one of biology's major unresolved questions. The field of evo-devo has shown that in early development changes to regulatory regions are the dominant mode of genetic change, but whether this extends to the evolution of novel phenotypes in the adult organism is unclear. Here, we conduct ten RNA-Seq experiments across both novel and conserved tissues in the honey bee to determine to what extent postdevelopmental novelty is based on changes to the coding regions of genes. We make several discoveries. First, we show that with respect to novel physiological functions in the adult animal, positively selected tissue-specific genes of high expression underlie novelty by conferring specialized cellular functions. Such genes are often, but not always taxonomically restricted genes (TRGs). We further show that positively selected genes, whether TRGs or conserved genes, are the least connected genes within gene expression networks. Overall, this work suggests that the evo-devo paradigm is limited, and that the evolution of novelty, postdevelopment, follows additional rules. Specifically, evo-devo stresses that high network connectedness (repeated use of the same gene in many contexts) constrains coding sequence change as it would lead to negative pleiotropic effects. Here, we show that in the adult animal, the converse is true: Genes with low network connectedness (TRGs and tissue-specific conserved genes) underlie novel phenotypes by rapidly changing coding sequence to perform new-specialized functions.

BackgroundStick and leaf insects (Phasmatodea) are an exclusively leaf-feeding order of insects with no record of omnivory, unlike other “herbivorous” Polyneoptera. They represent an ideal system for investigating the adaptations necessary for obligate folivory, including plant cell wall degrading enzymes (PCWDEs). However, their physiology and internal anatomy is poorly understood, with limited genomic resources available.ResultsWe de novo assembled transcriptomes for the anterior and posterior midguts of six diverse Phasmatodea species, with RNA-Seq on one exemplar species, Peruphasma schultei. The latter’s assembly yielded >100,000 transcripts, with over 4000 transcripts uniquely or more highly expressed in specific midgut sections. Two to three dozen PCWDE encoding gene families, including cellulases and pectinases, were differentially expressed in the anterior midgut. These genes were also found in genomic DNA from phasmid brain tissue, suggesting endogenous production. Sequence alignments revealed catalytic sites on most PCWDE transcripts. While most phasmid PCWDE genes showed homology with those of other insects, the pectinases were homologous to bacterial genes.ConclusionsWe identified a large and diverse PCWDE repertoire endogenous to the phasmids. If these expressed genes are translated into active enzymes, then phasmids can theoretically break plant cell walls into their monomer components independently of microbial symbionts. The differential gene expression between the two midgut sections provides the first molecular hints as to their function in living phasmids. Our work expands the resources available for industrial applications of animal-derived PCWDEs, and facilitates evolutionary analysis of lower Polyneopteran digestive enzymes, including the pectinases whose origin in Phasmatodea may have been a horizontal transfer event from bacteria.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-917) contains supplementary material, which is available to authorized users.

No abstract

BackgroundDNA methylation is a common regulator of gene expression, including acting as a regulator of developmental events and behavioral changes in adults. Using the unique system of genetic caste determination in Pogonomyrmex barbatus, we were able to document changes in DNA methylation during development, and also across both ancient and contemporary hybridization events.Methodology/Principal FindingsSodium bisulfite sequencing demonstrated in vivo methylation of symmetric CG dinucleotides in P. barbatus. We also found methylation of non-CpG sequences. This validated two bioinformatics methods for predicting gene methylation, the bias in observed to expected ratio of CpG dinucleotides and the density of CpG/TpG single nucleotide polymorphisms (SNP). Frequencies of genomic DNA methylation were determined for different developmental stages and castes using ms-AFLP assays. The genetic caste determination system (GCD) is probably the product of an ancestral hybridization event between P. barbatus and P. rugosus. Two lineages obligately co-occur within a GCD population, and queens are derived from intra-lineage matings whereas workers are produced from inter-lineage matings. Relative DNA methylation levels of queens and workers from GCD lineages (contemporary hybrids) were not significantly different until adulthood. Virgin queens had significantly higher relative levels of DNA methylation compared to workers. Worker DNA methylation did not vary among developmental stages within each lineage, but was significantly different between the currently hybridizing lineages. Finally, workers of the two genetic caste determination lineages had half as many methylated cytosines as workers from the putative parental species, which have environmental caste determination.Conclusions/SignificanceThese results suggest that DNA methylation may be a conserved regulatory mechanism moderating division of labor in both bees and ants. Current and historic hybridization appear to have altered genomic methylation levels suggesting a possible link between changes in overall DNA methylation and the origin and regulation of genetic caste determination in P. barbatus.

No abstract

No abstract

Honey bee queens undergo dramatic behavioral (e.g., reduced sexual receptivity), physiological (e.g., ovary activation, ovulation, and modulation of pheromone production) and transcriptional changes after they complete mating. To elucidate how queen post-mating changes are influenced by seminal fluid, the non-spermatozoa-containing component of semen, we injected queens with semen or seminal fluid alone. We assessed queen sexual receptivity (as measured by likelihood to take mating flights), ovary activation, worker retinue response (which is influenced by queen pheromone production), and transcriptional changes in queen abdominal fat body and brain tissues. Injection with either seminal fluid or semen resulted in decreased sexual receptivity, increased attractiveness of queens to workers, and altered expression of several genes that are also regulated by natural mating in queens. The post-mating and transcriptional changes of queens receiving seminal fluid were not significantly different from queens injected with semen, suggesting that components in seminal fluid, such as seminal fluid proteins, are largely responsible for stimulating post-mating changes in queens. Managed honey bee (Apis mellifera) colonies provide an estimated $175 billion in pollination services for a large variety of crops worldwide with an estimated $14.6 billion value in the United States alone 1-3. While the number of honey bee colonies has increased between 1961 and 2007 4 , the U.S. and some regions of Europe have exhibited high annual colony losses, with the U.S. experiencing an average annual loss of 33% recorded since 2006 5-13. Premature queen loss is one of the top causes of colony loss reported by beekeepers both in the U.S. and in Europe 6,7,13. Therefore, improving colony health and survivorship in large part will depend on understanding factors affecting queen health and productivity as well as investigating the role of drones (haploid male bees) on queen health. A honey bee colony consists of a single egg-laying queen, hundreds to thousands of seasonal drones, and tens of thousands of sterile female workers 14. Approximately 1 week after emergence as an adult, the queen will take one or more mating flights to large congregations of males called drone congregation areas (DCAs) 14,15. Honey bee queens are polyandrous, which means they mate with multiple drones; queens have recently been reported to mate as many as 34-77 drones 16. During copulation, drone semen is transferred to the queen's oviduct where only 3-5% of each drone's spermatozoa are stored in a specialized organ called the spermatheca 14,15,17,18. The process of mating initiates dramatic behavioral, physiological, and molecular changes in queens 19-24. Behaviorally, mated queens become sexually unreceptive, as evidenced by the cessation of mating flights and their reduced attraction to light 14,25. The ovaries of newly-mated queens begin activating and their ovarioles thicken and start to produce eggs, upon which the queens initiate egg-laying and continue to do so...

The sequencing depth necessary for documenting differential gene expression using RNA-Seq has been little explored outside of model systems. In particular, the depth required to analyze large-scale patterns of differential transcription factor expression is not known. The goal of the present study is to explore the effectiveness of shallow (relatively low read depth) RNA-Seq. We focus on two tissues in the honey bee: the sting gland and the digestive tract. The sting gland is an experimentally well-understood tissue that we use to benchmark the utility of this approach. We use the digestive tract to test the results obtained with the sting gland, and to conduct RNA-Seq between tissue types. Using a list of experimentally verified genes conferring tissue-specific functions in the sting gland, we show that relatively little read depth is necessary to identify them. We argue that this result should be broadly applicable, since genes important for tissue-specific functions often have robust expression patterns, and because we obtained similar results in our analysis of the digestive tract. Furthermore, we demonstrate that the differential expression of transcription factors, which are transcribed at low levels compared to other genes, can nevertheless often be determined using shallow RNA-Seq. Overall, we find over 150 differentially expressed transcription factors in our tissues at a read depth of only 12 million. This work shows the utility of low-depth sequencing for identifying genes important for tissue-specific functions. It also verifies the often-held belief that transcription factors show low levels of expression, while demonstrating that, in spite of this, they are frequently amenable to shallow RNA-Seq. Our findings should be of benefit to researchers using RNA-Seq in many different biological systems.