Exclusion from a social group is an effective way to avoid parasite transmission. This type of social removal has also been proposed as a form of collective defense, or social immunity, in eusocial insect groups. If parasitic modification of host behavior is widespread in social insects, the underlying physiological and neuronal mechanisms remain to be investigated. We studied this phenomenon in honey bees parasitized by the mite Varroa destructor or microsporidia Nosema ceranae, which make bees leave the hive precociously. We characterized the chemical, behavioral and neurogenomic changes in parasitized bees, and compared the effects of both parasites. Analysis of cuticular hydrocarbon (CHC) profiles by gas chromatography coupled with mass spectrophotometry (GC-MS) showed changes in honey bees parasitized by either Nosema ceranae or Varroa destructor after 5 days of infestation. Levels of 10-HDA, an antiseptic important for social immunity, did not change in response to parasitism. Behavioral analysis of N. ceranae- or V. destructor- parasitized bees revealed no significant differences in their behavioral acts or social interactions with nestmates. Digital gene expression (DGE) analysis of parasitized honey bee brains demonstrated that, despite the difference in developmental stage at which the bee is parasitized, Nosema and Varroa-infested bees shared more gene changes with each other than with honey bee brain expression gene sets for forager or nurse castes. Parasitism by Nosema or Varroa induces changes to both the CHC profiles on the surface of the bee and transcriptomic profiles in the brain, but within the social context of the hive, does not result in observable effects on her behavior or behavior towards her. While parasitized bees are reported to leave the hive as foragers, their brain transcription profiles suggest that their behavior is not driven by the same molecular pathways that induce foraging behavior.

Plant life-long organogenesis involves sequential, time and tissue specific expression of developmental genes. This requires activities of Polycomb Group (PcG) and trithorax Group complexes (trxG), respectively responsible for repressive Histone 3 trimethylation at lysine 27 (H3K27me3) and activation-related Histone 3 trimethylation at lysine 4 (H3K4me3). However, the genome-wide dynamics in histone modifications that occur during developmental processes have remained elusive. Here, we report the distributions of H3K27me3 and H3K4me3 along with expression changes, in a developmental series including Arabidopsis thaliana leaf and three stages of flower development. We found that chromatin mark levels are highly dynamic over the time series on nearly half of all Arabidopsis genes. Moreover, during early flower morphogenesis, changes in H3K4me3 prevail over changes in H3K27me3 and quantitatively correlate with expression changes, while H3K27me3 changes occur later. Notably, we found that H3K4me3 increase during the early activation of PcG target genes while H3K27me3 level remain relatively constant at the locus. Our results reveal that H3K4me3 predicts changes in gene expression better than H3K27me3, unveil unexpected chromatin mechanisms at gene activation and underline the relevance of tissue-specific temporal epigenomics.

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. 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. 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.

Sexual differentiation of inflorescences and flowers is important for reproduction and affects crop plant productivity. We report here on a molecular study of the process of sexual differentiation in the immature inflorescence of oil palm (Elaeis guineensis). This species is monoecious and exhibits gender diphasy, producing male and female inflorescences separately on the same plant in alternation. Three main approaches were used: small RNA-seq to characterise and study the expression of miRNA genes; RNA-seq to monitor mRNA accumulation patterns; hormone quantification to assess the role of cytokinins and auxins in inflorescence differentiation. Our study allowed the characterisation of 30 previously unreported palm MIRNA genes. In differential gene and miRNA expression studies, we identified a number of key developmental genes and miRNA-mRNA target modules previously described in relation to their developmental regulatory role in the cereal panicle, notably the miR156/529/535-SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE (SPL) gene regulatory module. Gene enrichment analysis highlighted the importance of hormone-related genes, and this observation was corroborated by the detection of much higher levels of cytokinins in the female inflorescence. Our data illustrate the importance of branching regulation within the developmental window studied, during which the female inflorescence, unlike its male counterpart, produces flower clusters on new successive axes by sympodial growth.

Over the past decade, the spotted wing Drosophila, Drosophila suzukii, has invaded Europe and America and has become a major agricultural pest in these areas, thereby prompting intense research activities to better understand its biology. Two draft genome assemblies based on short-read sequencing were released in 2013 for this species. Although valuable, these resources contain pervasive assembly errors and are highly fragmented, two features limiting their values. Our purpose here was to improve the assembly of theD. suzukiigenome. For this, we generated PacBio long-read sequencing data at 160X sequence coverage and assembled a novel, contiguousD. suzukiigenome. We obtained a high-quality assembly of 270 Mb (with 546 contigs, a N50 of 2.6Mb, a L50 of 15, and a BUSCO score of 95%) that we called WT3-2.0. We found that despite 16 rounds of full-sib crossings theD. suzukiistrain that we sequenced has maintained high levels of polymorphism in some regions of its genome (ca. 19Mb). As a consequence, the quality of the assembly of these regions was reduced. We explored possible origins of this high residual diversity, including the presence of structural variants and a possible heterogeneous admixture pattern of North American and Asian ancestry. Overall, our WT3-2.0 assembly provides a higher quality genomic resource compared to the previous one in terms of general assembly statistics, sequence quality and gene annotation. This newD. suzukiigenome assembly is therefore an improved resource for high-throughput sequencing approaches, as well as manipulative genetic technologies to studyD. suzukii.

Evidence is accumulating that evolutionary changes are not only common during biological invasions but may also contribute directly to invasion success. The genomic basis of such changes is still largely unexplored. Yet, understanding the genomic response to invasion may help to predict the conditions under which invasiveness can be enhanced or suppressed. Here we characterized the genome response of the spotted wing drosophilaDrosophila suzukiiduring the worldwide invasion of this pest insect species, by conducting a genome-wide association study to identify genes involved in adaptive processes during invasion. Genomic data from 22 population samples were analyzed to detect genetic variants associated with the status (invasive versus native) of the sampled populations based on a newly developed statistic, we calledC₂, that contrasts allele frequencies corrected for population structure. This new statistical framework has been implemented in an upgraded version of the program BayPass. We identified a relatively small set of single nucleotide polymorphisms (SNPs) that show a highly significant association with the invasive status of populations. In particular, two genesRhoGEF64Candcpo, the latter contributing to natural variation in several life-history traits (including diapause) inDrosophila melanogaster, contained SNPs significantly associated with the invasive status in the two separate main invasion routes ofD. suzukii. Our methodological approaches can be applied to any other invasive species, and more generally to any evolutionary model for species characterized by non-equilibrium demographic conditions for which binary covariables of interest can be defined at the population level.