Identification of a viral gene essential for the genome replication of a domesticated endogenous virus in ichneumonid parasitoid wasps.

Thousands of endoparasitoid wasp species in the families Braconidae and Ichneumonidae harbor "domesticated endogenous viruses" (DEVs) in their genomes. This study focuses on ichneumonid DEVs, named ichnoviruses (IVs). Large quantities of DNA-containing IV virions are produced in ovary calyx cells during the pupal and adult stages of female wasps. Females parasitize host insects by injecting eggs and virions into the body cavity. After injection, virions rapidly infect host cells which is followed by expression of IV genes that promote the successful development of wasp offspring. IV genomes consist of two components: proviral segment loci that serve as templates for circular dsDNAs that are packaged into capsids, and genes from an ancestral virus that produce virions. In this study, we generated a chromosome-scale genome assembly for Hyposoter didymator that harbors H. didymator ichnovirus (HdIV). We identified a total of 67 HdIV loci that are amplified in calyx cells during the wasp pupal stage. We then focused on an HdIV gene, U16, which is transcribed in calyx cells during the initial stages of replication. Sequence analysis indicated that U16 contains a conserved domain in primases from select other viruses. Knockdown of U16 by RNA interference inhibited virion morphogenesis in calyx cells. Genome-wide analysis indicated U16 knockdown also inhibited amplification of HdIV loci in calyx cells. Altogether, our results identified several previously unknown HdIV loci, demonstrated that all HdIV loci are amplified in calyx cells during the pupal stage, and showed that U16 is required for amplification and virion morphogenesis.

Integration of temporal and spatial patterning generates neural diversity.

In the Drosophila optic lobes, 800 retinotopically organized columns in the medulla act as functional units for processing visual information. The medulla contains over 80 types of neuron, which belong to two classes: uni-columnar neurons have a stoichiometry of one per column, while multi-columnar neurons contact multiple columns. Here we show that combinatorial inputs from temporal and spatial axes generate this neuronal diversity: all neuroblasts switch fates over time to produce different neurons; the neuroepithelium that generates neuroblasts is also subdivided into six compartments by the expression of specific factors. Uni-columnar neurons are produced in all spatial compartments independently of spatial input; they innervate the neuropil where they are generated. Multi-columnar neurons are generated in smaller numbers in restricted compartments and require spatial input; the majority of their cell bodies subsequently move to cover the entire medulla. The selective integration of spatial inputs by a fixed temporal neuroblast cascade thus acts as a powerful mechanism for generating neural diversity, regulating stoichiometry and the formation of retinotopy.

H3K9me2 genome-wide distribution in the holocentric insect Spodoptera frugiperda (Lepidoptera: Noctuidae).

BACKGROUND: Eukaryotic genomes are packaged by Histone proteins in a structure called chromatin. There are different chromatin types. Euchromatin is typically associated with decondensed, transcriptionally active regions and heterochromatin to more condensed regions of the chromosomes. Methylation of Lysine 9 of Histone H3 (H3K9me) is a conserved biochemical marker of heterochromatin. In many organisms, heterochromatin is usually localized at telomeric as well as pericentromeric regions but can also be found at interstitial chromosomal loci. This distribution may vary in different species depending on their general chromosomal organization. Holocentric species such as Spodoptera frugiperda (Lepidoptera: Noctuidae) possess dispersed centromeres instead of a monocentric one and thus no observable pericentromeric compartment. To identify the localization of heterochromatin in such species we performed ChIP-Seq experiments and analyzed the distribution of the heterochromatin marker H3K9me2 in the Sf9 cell line and whole 4th instar larvae (L4) in relation to RNA-Seq data. RESULTS: In both samples we measured an enrichment of H3K9me2 at the (sub) telomeres, rDNA loci, and satellite DNA sequences, which could represent dispersed centromeric regions. We also observed that density of H3K9me2 is positively correlated with transposable elements and protein-coding genes. But contrary to most model organisms, H3K9me2 density is not correlated with transcriptional repression. CONCLUSION: This is the first genome-wide ChIP-Seq analysis conducted in S. frugiperda for H3K9me2. Compared to model organisms, this mark is found in expected chromosomal compartments such as rDNA and telomeres. However, it is also localized at numerous dispersed regions, instead of the well described large pericentromeric domains, indicating that H3K9me2 might not represent a classical heterochromatin marker in Lepidoptera. (242 words).

Integrated miRNA and transcriptome profiling to explore the molecular determinism of convergent adaptation to corn in two lepidopteran pests of agriculture.

BACKGROUND: The degree to which adaptation to same environment is determined by similar molecular mechanisms, is a topic of broad interest in evolutionary biology, as an indicator of evolutionary predictability. We wished to address if adaptation to the same host plant in phytophagous insects involved related gene expression patterns. We compared sRNA-Seq and RNA-Seq data between two pairs of taxa of Ostrinia and Spodoptera frugiperda sharing maize as host-plant. For the latter, we had previously carried out a reciprocal transplant experiment by feeding of the larvae of the Corn strain (Sf-C) and the Rice strain (Sf-R) on corn versus rice and characterized the mRNA and miRNA responses. RESULTS: First, we predicted the genes encoding miRNA in Ostrinia nubilalis (On) and O. scapulalis (Os). Respectively 67 and 65 known miRNA genes, as well as 196 and 190 novel ones were predicted with Os genome using sncRNAs extracted from whole larvae feeding on corn or mugwort. In On, a read counts analysis showed that 37 (55.22%) known miRNAs and 19 (9.84%) novel miRNAs were differentially expressed (DE) on mugwort compared to corn (in Os, 25 known miRs (38.46%) and 8 novel ones (4.34%)). Between species on corn, 8 (12.5%) known miRNAs and 8 (6.83%) novel ones were DE while only one novel miRNA showed expression variation between species on mugwort. Gene target prediction led to the identification of 2953 unique target genes in On and 2719 in Os, among which 11.6% (344) were DE when comparing species on corn. 1.8% (54) of On miR targets showed expression variation upon a change of host-plant. We found molecular changes matching convergent phenotype, i.e., a set of nine miRNAs that are regulated either according to the host-plant both in On and Sf-C or between them on the same plant, corn. Among DE miR target genes between taxa, 13.7% shared exactly the same annotation between the two pairs of taxa and had function related to insect host-plant interaction. CONCLUSION: There is some similarity in underlying genetic mechanisms of convergent evolution of two distant Lepidopteran species having adopted corn in their host range, highlighting possible adaptation genes.

A novel reference dated phylogeny for the genus Spodoptera Guenée (Lepidoptera: Noctuidae: Noctuinae): new insights into the evolution of a pest-rich genus.

The noctuid genus Spodoptera currently consists of 31 species with varied host plant breadths, ranging from monophagous and oligophagous non-pest species to polyphagous pests of economic importance. Several of these pest species have become major invaders, colonizing multiple continents outside their native range. Such is the case of the infamous fall armyworm, Spodoptera frugiperda (J.E. Smith), which includes two recognized host strains that have not been treated as separate species. Following its accidental introduction to Africa in 2016, it quickly spread through Africa and Asia to Australia. Given that half the described Spodoptera species cause major crop losses, comparative genomics studies of several Spodoptera species have highlighted major adaptive changes in genetic architecture, possibly relating to their pest status. Several recent population genomics studies conducted on two species enable a more refined understanding of their population structures, migration patterns and invasion processes. Despite growing interest in the genus, the taxonomic status of several Spodoptera species remains unstable and evolutionary studies suffer from the absence of a robust and comprehensive dated phylogenetic framework. We generated mitogenomic data for 14 Spodoptera taxa, which are combined with data from 15 noctuoid outgroups to generate a resolved mitogenomic backbone phylogeny using both concatenation and multi-species coalescent approaches. We combine this backbone with additional mitochondrial and nuclear data to improve our understanding of the evolutionary history of the genus. We also carry out comprehensive dating analyses, which implement three distinct calibration strategies based on either primary or secondary fossil calibrations. Our results provide an updated phylogenetic framework for 28 Spodoptera species, identifying two well-supported ecologically diverse clades that are recovered for the first time. Well-studied larvae in each of these clades are characterized by differences in mandibular shape, with one clade's being more specialized on silica-rich C4 grasses. Interestingly, the inferred timeframe for the genus suggests an earlier origin than previously thought for the genus: about 17-18 million years ago.

Positive selection alone is sufficient for whole genome differentiation at the early stage of speciation process in the fall armyworm.

BACKGROUND: The process of speciation involves differentiation of whole genome sequences between a pair of diverging taxa. In the absence of a geographic barrier and in the presence of gene flow, genomic differentiation may occur when the homogenizing effect of recombination is overcome across the whole genome. The fall armyworm is observed as two sympatric strains with different host-plant preferences across the entire habitat. These two strains exhibit a very low level of genetic differentiation across the whole genome, suggesting that genomic differentiation occurred at an early stage of speciation. In this study, we aim at identifying critical evolutionary forces responsible for genomic differentiation in the fall armyworm. RESULTS: These two strains exhibit a low level of genomic differentiation (FST = 0.0174), while 99.2% of 200 kb windows have genetically differentiated sequences (FST > 0). We found that the combined effect of mild positive selection and genetic linkage to selectively targeted loci are responsible for the genomic differentiation. However, a single event of very strong positive selection appears not to be responsible for genomic differentiation. The contribution of chromosomal inversions or tight genetic linkage among positively selected loci causing reproductive barriers is not supported by our data. Phylogenetic analysis shows that the genomic differentiation occurred by sub-setting of genetic variants in one strain from the other. CONCLUSIONS: From these results, we concluded that genomic differentiation may occur at the early stage of a speciation process in the fall armyworm and that mild positive selection targeting many loci alone is sufficient evolutionary force for generating the pattern of genomic differentiation. This genomic differentiation may provide a condition for accelerated genomic differentiation by synergistic effects among linkage disequilibrium generated by following events of positive selection. Our study highlights genomic differentiation as a key evolutionary factor connecting positive selection to divergent selection.

Control of mitochondrial structure and function by the Yorkie/YAP oncogenic pathway.

Mitochondrial structure and function are highly dynamic, but the potential roles for cell signaling pathways in influencing these properties are not fully understood. Reduced mitochondrial function has been shown to cause cell cycle arrest, and a direct role of signaling pathways in controlling mitochondrial function during development and disease is an active area of investigation. Here, we show that the conserved Yorkie/YAP signaling pathway implicated in the control of organ size also functions in the regulation of mitochondria in Drosophila as well as human cells. In Drosophila, activation of Yorkie causes direct transcriptional up-regulation of genes that regulate mitochondrial fusion, such as opa1-like (opa1) and mitochondria assembly regulatory factor (Marf), and results in fused mitochondria with dramatic reduction in reactive oxygen species (ROS) levels. When mitochondrial fusion is genetically attenuated, the Yorkie-induced cell proliferation and tissue overgrowth are significantly suppressed. The function of Yorkie is conserved across evolution, as activation of YAP2 in human cell lines causes increased mitochondrial fusion. Thus, mitochondrial fusion is an essential and direct target of the Yorkie/YAP pathway in the regulation of organ size control during development and could play a similar role in the genesis of cancer.

Control of nerve cord formation by Engrailed and Gooseberry-Neuro: A multi-step, coordinated process.

One way to better understand the molecular mechanisms involved in the construction of a nervous system is to identify the downstream effectors of major regulatory proteins. We previously showed that Engrailed (EN) and Gooseberry-Neuro (GsbN) transcription factors act in partnership to drive the formation of posterior commissures in the central nervous system of Drosophila. In this report, we identified genes regulated by both EN and GsbN through chromatin immunoprecipitation ("ChIP on chip") and transcriptome experiments, combined to a genetic screen relied to the gene dose titration method. The genomic-scale approaches allowed us to define 175 potential targets of EN-GsbN regulation. We chose a subset of these genes to examine ventral nerve cord (VNC) defects and found that half of the mutated targets show clear VNC phenotypes when doubly heterozygous with en or gsbn mutations, or when homozygous. This strategy revealed new groups of genes never described for their implication in the construction of the nerve cord. Their identification suggests that, to construct the nerve cord, EN-GsbN may act at three levels, in: (i) sequential control of the attractive-repulsive signaling that ensures contralateral projection of the commissural axons, (ii) temporal control of the translation of some mRNAs, (iii) regulation of the capability of glial cells to act as commissural guideposts for developing axons. These results illustrate how an early, coordinated transcriptional control may orchestrate the various mechanisms involved in the formation of stereotyped neuronal networks. They also validate the overall strategy to identify genes that play crucial role in axonal pathfinding.

Characterization and expression profiling of microRNAs in response to plant feeding in two host-plant strains of the lepidopteran pest Spodoptera frugiperda.

BACKGROUND: A change in the environment may impair development or survival of living organisms leading them to adapt to the change. The resulting adaptation trait may reverse, or become fixed in the population leading to evolution of species. Deciphering the molecular basis of adaptive traits can thus give evolutionary clues. In phytophagous insects, a change in host-plant range can lead to emergence of new species. Among them, Spodoptera frugiperda is a major agricultural lepidopteran pest consisting of two host-plant strains having diverged 3 MA, based on mitochondrial markers. In this paper, we address the role of microRNAs, important gene expression regulators, in response to host-plant change and in adaptive evolution. RESULTS: Using small RNA sequencing, we characterized miRNA repertoires of the corn (C) and rice (R) strains of S. frugiperda, expressed during larval development on two different host-plants, corn and rice, in the frame of reciprocal transplant experiments. We provide evidence for 76 and 68 known miRNAs in C and R strains and 139 and 171 novel miRNAs. Based on read counts analysis, 34 of the microRNAs were differentially expressed in the C strain larvae fed on rice as compared to the C strain larvae fed on corn. Twenty one were differentially expressed on rice compared to corn in R strain. Nine were differentially expressed in the R strain compared to C strain when reared on corn. A similar ratio of microRNAs was differentially expressed between strains on rice. We could validate experimentally by QPCR, variation in expression of the most differentially expressed candidates. We used bioinformatics methods to determine the target mRNAs of known microRNAs. Comparison with the mRNA expression profile during similar reciprocal transplant experiment revealed potential mRNA targets of these host-plant regulated miRNAs. CONCLUSIONS: In the current study, we performed the first systematic analysis of miRNAs in Lepidopteran pests feeding on host-plants. We identified a set of the differentially expressed miRNAs that respond to the plant diet, or differ constitutively between the two host plant strains. Among the latter, the ones that are also deregulated in response to host-plant are molecular candidates underlying a complex adaptive trait.

Evolution of H3K27me3-marked chromatin is linked to gene expression evolution and to patterns of gene duplication and diversification.

Histone modifications are critical for the regulation of gene expression, cell type specification, and differentiation. However, evolutionary patterns of key modifications that regulate gene expression in differentiating organisms have not been examined. Here we mapped the genomic locations of the repressive mark histone 3 lysine 27 trimethylation (H3K27me3) in four species of Drosophila, and compared these patterns to those in C. elegans. We found that patterns of H3K27me3 are highly conserved across species, but conservation is substantially weaker among duplicated genes. We further discovered that retropositions are associated with greater evolutionary changes in H3K27me3 and gene expression than tandem duplications, indicating that local chromatin constraints influence duplicated gene evolution. These changes are also associated with concomitant evolution of gene expression. Our findings reveal the strong conservation of genomic architecture governed by an epigenetic mark across distantly related species and the importance of gene duplication in generating novel H3K27me3 profiles.

Diverse patterns of genomic targeting by transcriptional regulators in Drosophila melanogaster.

Annotation of regulatory elements and identification of the transcription-related factors (TRFs) targeting these elements are key steps in understanding how cells interpret their genetic blueprint and their environment during development, and how that process goes awry in the case of disease. One goal of the modENCODE (model organism ENCyclopedia of DNA Elements) Project is to survey a diverse sampling of TRFs, both DNA-binding and non-DNA-binding factors, to provide a framework for the subsequent study of the mechanisms by which transcriptional regulators target the genome. Here we provide an updated map of the Drosophila melanogaster regulatory genome based on the location of 84 TRFs at various stages of development. This regulatory map reveals a variety of genomic targeting patterns, including factors with strong preferences toward proximal promoter binding, factors that target intergenic and intronic DNA, and factors with distinct chromatin state preferences. The data also highlight the stringency of the Polycomb regulatory network, and show association of the Trithorax-like (Trl) protein with hotspots of DNA binding throughout development. Furthermore, the data identify more than 5800 instances in which TRFs target DNA regions with demonstrated enhancer activity. Regions of high TRF co-occupancy are more likely to be associated with open enhancers used across cell types, while lower TRF occupancy regions are associated with complex enhancers that are also regulated at the epigenetic level. Together these data serve as a resource for the research community in the continued effort to dissect transcriptional regulatory mechanisms directing Drosophila development.

A cis-regulatory map of the Drosophila genome.

Systematic annotation of gene regulatory elements is a major challenge in genome science. Direct mapping of chromatin modification marks and transcriptional factor binding sites genome-wide has successfully identified specific subtypes of regulatory elements. In Drosophila several pioneering studies have provided genome-wide identification of Polycomb response elements, chromatin states, transcription factor binding sites, RNA polymerase II regulation and insulator elements; however, comprehensive annotation of the regulatory genome remains a significant challenge. Here we describe results from the modENCODE cis-regulatory annotation project. We produced a map of the Drosophila melanogaster regulatory genome on the basis of more than 300 chromatin immunoprecipitation data sets for eight chromatin features, five histone deacetylases and thirty-eight site-specific transcription factors at different stages of development. Using these data we inferred more than 20,000 candidate regulatory elements and validated a subset of predictions for promoters, enhancers and insulators in vivo. We identified also nearly 2,000 genomic regions of dense transcription factor binding associated with chromatin activity and accessibility. We discovered hundreds of new transcription factor co-binding relationships and defined a transcription factor network with over 800 potential regulatory relationships.

Chromatin insulator factors involved in long-range DNA interactions and their role in the folding of the Drosophila genome.

Chromatin insulators are genetic elements implicated in the organization of chromatin and the regulation of transcription. In Drosophila, different insulator types were characterized by their locus-specific composition of insulator proteins and co-factors. Insulators mediate specific long-range DNA contacts required for the three dimensional organization of the interphase nucleus and for transcription regulation, but the mechanisms underlying the formation of these contacts is currently unknown. Here, we investigate the molecular associations between different components of insulator complexes (BEAF32, CP190 and Chromator) by biochemical and biophysical means, and develop a novel single-molecule assay to determine what factors are necessary and essential for the formation of long-range DNA interactions. We show that BEAF32 is able to bind DNA specifically and with high affinity, but not to bridge long-range interactions (LRI). In contrast, we show that CP190 and Chromator are able to mediate LRI between specifically-bound BEAF32 nucleoprotein complexes in vitro. This ability of CP190 and Chromator to establish LRI requires specific contacts between BEAF32 and their C-terminal domains, and dimerization through their N-terminal domains. In particular, the BTB/POZ domains of CP190 form a strict homodimer, and its C-terminal domain interacts with several insulator binding proteins. We propose a general model for insulator function in which BEAF32/dCTCF/Su(HW) provide DNA specificity (first layer proteins) whereas CP190/Chromator are responsible for the physical interactions required for long-range contacts (second layer). This network of organized, multi-layer interactions could explain the different activities of insulators as chromatin barriers, enhancer blockers, and transcriptional regulators, and suggest a general mechanism for how insulators may shape the organization of higher-order chromatin during cell division.

Divergent transcriptional regulatory logic at the intersection of tissue growth and developmental patterning.

The Yorkie/Yap transcriptional coactivator is a well-known regulator of cellular proliferation in both invertebrates and mammals. As a coactivator, Yorkie (Yki) lacks a DNA binding domain and must partner with sequence-specific DNA binding proteins in the nucleus to regulate gene expression; in Drosophila, the developmental regulators Scalloped (Sd) and Homothorax (Hth) are two such partners. To determine the range of target genes regulated by these three transcription factors, we performed genome-wide chromatin immunoprecipitation experiments for each factor in both the wing and eye-antenna imaginal discs. Strong, tissue-specific binding patterns are observed for Sd and Hth, while Yki binding is remarkably similar across both tissues. Binding events common to the eye and wing are also present for Sd and Hth; these are associated with genes regulating cell proliferation and "housekeeping" functions, and account for the majority of Yki binding. In contrast, tissue-specific binding events for Sd and Hth significantly overlap enhancers that are active in the given tissue, are enriched in Sd and Hth DNA binding sites, respectively, and are associated with genes that are consistent with each factor's previously established tissue-specific functions. Tissue-specific binding events are also significantly associated with Polycomb targeted chromatin domains. To provide mechanistic insights into tissue-specific regulation, we identify and characterize eye and wing enhancers of the Yki-targeted bantam microRNA gene and demonstrate that they are dependent on direct binding by Hth and Sd, respectively. Overall these results suggest that both Sd and Hth use distinct strategies - one shared between tissues and associated with Yki, the other tissue-specific, generally Yki-independent and associated with developmental patterning - to regulate distinct gene sets during development.

Insensitive is a corepressor for Suppressor of Hairless and regulates Notch signalling during neural development.

The Notch intracellular domain functions as a co-activator for the DNA-binding protein Suppressor of Hairless (Su(H)) to mediate myriad cell fate decisions. Notch pathway activity is balanced by transcriptional repression, mediated by Su(H) in concert with its Drosophila corepressor Hairless. We demonstrate that the Drosophila neural BEN-solo protein Insensitive (Insv) is a nuclear factor that inhibits Notch signalling during multiple peripheral nervous system cell fate decisions. Endogenous Insv was particularly critical when repressor activity of Su(H) was compromised. Reciprocally, ectopic Insv generated several Notch loss-of-function phenotypes, repressed most Notch targets in the E(spl)-C, and opposed Notch-mediated activation of an E(spl)m3-luc reporter. A direct role for Insv in transcriptional repression was indicated by binding of Insv to Su(H), and by strong chromatin immunoprecipitation of endogenous Insv to most E(spl)-C loci. Strikingly, ectopic Insv fully rescued sensory organ precursors in Hairless null clones, indicating that Insv can antagonize Notch independently of Hairless. These data shed first light on the in vivo function for a BEN-solo protein as an Su(H) corepressor in the Notch pathway regulating neural development.

Systematic evaluation of factors influencing ChIP-seq fidelity.

We evaluated how variations in sequencing depth and other parameters influence interpretation of chromatin immunoprecipitation-sequencing (ChIP-seq) experiments. Using Drosophila melanogaster S2 cells, we generated ChIP-seq data sets for a site-specific transcription factor (Suppressor of Hairy-wing) and a histone modification (H3K36me3). We detected a chromatin-state bias: open chromatin regions yielded higher coverage, which led to false positives if not corrected. This bias had a greater effect on detection specificity than any base-composition bias. Paired-end sequencing revealed that single-end data underestimated ChIP-library complexity at high coverage. Removal of reads originating at the same base reduced false-positives but had little effect on detection sensitivity. Even at mappable-genome coverage depth of ∼1 read per base pair, ∼1% of the narrow peaks detected on a tiling array were missed by ChIP-seq. Evaluation of widely used ChIP-seq analysis tools suggests that adjustments or algorithm improvements are required to handle data sets with deep coverage.

Adaptive evolution and the birth of CTCF binding sites in the Drosophila genome.

Changes in the physical interaction between cis-regulatory DNA sequences and proteins drive the evolution of gene expression. However, it has proven difficult to accurately quantify evolutionary rates of such binding change or to estimate the relative effects of selection and drift in shaping the binding evolution. Here we examine the genome-wide binding of CTCF in four species of Drosophila separated by between ∼2.5 and 25 million years. CTCF is a highly conserved protein known to be associated with insulator sequences in the genomes of human and Drosophila. Although the binding preference for CTCF is highly conserved, we find that CTCF binding itself is highly evolutionarily dynamic and has adaptively evolved. Between species, binding divergence increased linearly with evolutionary distance, and CTCF binding profiles are diverging rapidly at the rate of 2.22% per million years (Myr). At least 89 new CTCF binding sites have originated in the Drosophila melanogaster genome since the most recent common ancestor with Drosophila simulans. Comparing these data to genome sequence data from 37 different strains of Drosophila melanogaster, we detected signatures of selection in both newly gained and evolutionarily conserved binding sites. Newly evolved CTCF binding sites show a significantly stronger signature for positive selection than older sites. Comparative gene expression profiling revealed that expression divergence of genes adjacent to CTCF binding site is significantly associated with the gain and loss of CTCF binding. Further, the birth of new genes is associated with the birth of new CTCF binding sites. Our data indicate that binding of Drosophila CTCF protein has evolved under natural selection, and CTCF binding evolution has shaped both the evolution of gene expression and genome evolution during the birth of new genes.

The transcription factor ultraspiracle influences honey bee social behavior and behavior-related gene expression.

Behavior is among the most dynamic animal phenotypes, modulated by a variety of internal and external stimuli. Behavioral differences are associated with large-scale changes in gene expression, but little is known about how these changes are regulated. Here we show how a transcription factor (TF), ultraspiracle (usp; the insect homolog of the Retinoid X Receptor), working in complex transcriptional networks, can regulate behavioral plasticity and associated changes in gene expression. We first show that RNAi knockdown of USP in honey bee abdominal fat bodies delayed the transition from working in the hive (primarily "nursing" brood) to foraging outside. We then demonstrate through transcriptomics experiments that USP induced many maturation-related transcriptional changes in the fat bodies by mediating transcriptional responses to juvenile hormone. These maturation-related transcriptional responses to USP occurred without changes in USP's genomic binding sites, as revealed by ChIP-chip. Instead, behaviorally related gene expression is likely determined by combinatorial interactions between USP and other TFs whose cis-regulatory motifs were enriched at USP's binding sites. Many modules of JH- and maturation-related genes were co-regulated in both the fat body and brain, predicting that usp and cofactors influence shared transcriptional networks in both of these maturation-related tissues. Our findings demonstrate how "single gene effects" on behavioral plasticity can involve complex transcriptional networks, in both brain and peripheral tissues.

Genomic data integration for ecological and evolutionary traits in non-model organisms.

Why is it needed to develop system biology initiatives such as ENCODE on non-model organisms?

Establishment and analysis of a reference transcriptome for Spodoptera frugiperda.

BACKGROUND: Spodoptera frugiperda (Noctuidae) is a major agricultural pest throughout the American continent. The highly polyphagous larvae are frequently devastating crops of importance such as corn, sorghum, cotton and grass. In addition, the Sf9 cell line, widely used in biochemistry for in vitro protein production, is derived from S. frugiperda tissues. Many research groups are using S. frugiperda as a model organism to investigate questions such as plant adaptation, pest behavior or resistance to pesticides. RESULTS: In this study, we constructed a reference transcriptome assembly (Sf_TR2012b) of RNA sequences obtained from more than 35 S. frugiperda developmental time-points and tissue samples. We assessed the quality of this reference transcriptome by annotating a ubiquitous gene family--ribosomal proteins--as well as gene families that have a more constrained spatio-temporal expression and are involved in development, immunity and olfaction. We also provide a time-course of expression that we used to characterize the transcriptional regulation of the gene families studied. CONCLUSION: We conclude that the Sf_TR2012b transcriptome is a valid reference transcriptome. While its reliability decreases for the detection and annotation of genes under strong transcriptional constraint we still recover a fair percentage of tissue-specific transcripts. That allowed us to explore the spatial and temporal expression of genes and to observe that some olfactory receptors are expressed in antennae and palps but also in other non related tissues such as fat bodies. Similarly, we observed an interesting interplay of gene families involved in immunity between fat bodies and antennae.