Active genic machinery for epigenetic RNA modifications in bees.

Epitranscriptomics is an emerging field of investigation dedicated to the study of post-transcriptional RNA modifications. RNA methylations regulate RNA metabolism and processing, including changes in response to environmental cues. Although RNA modifications are conserved from bacteria to eukaryotes, there is little evidence of an epitranscriptomic pathway in insects. Here we identified genes related to RNA m6 A (N6-methyladenine) and m5 C (5-methylcytosine) methylation machinery in seven bee genomes (Apis mellifera, Melipona quadrifasciata, Frieseomelitta varia, Eufriesea mexicana, Bombus terrestris, Megachile rotundata and Dufourea novaeangliae). In A. mellifera, we validated the expression of methyltransferase genes and found that the global levels of m6 A and m5 C measured in the fat body and brain of adult workers differ significantly. Also, m6 A levels were differed significantly mainly between the fourth larval instar of queens and workers. Moreover, we found a conserved m5 C site in the honeybee 28S rRNA. Taken together, we confirm the existence of epitranscriptomic machinery acting in bees and open avenues for future investigations on RNA epigenetics in a wide spectrum of hymenopteran species.

Expression profiles of neotropical termites reveal microbiota-associated, caste-biased genes and biotechnological targets.

Termites are well recognized by their complex development trajectories, involving dynamic differentiation process between non-reproductive castes, workers and soldiers. These insects are associated with endosymbiotic microorganisms, which help in lignocellulose digestion and nitrogen metabolism. Aiming to identify genes harbouring biotechnological potential, we analyzed workers and soldiers RNA-Seq data of three neotropical termites: Heterotermes tenuis (Isoptera: Rhinotermitidae), Velocitermes heteropterus (Isoptera: Termitidae) and Cornitermes cumulans (Isoptera: Termitidae). We observed differences in the microbiota associated with each termite family, and found protists' genes in both Termitidae species. We found an opposite pattern of caste-biased gene expression between H. tenuis and the termitids studied. Moreover, the two termitids are considerably different concerning the number of differentially expressed genes (DEGs). Functional annotation indicated considerable differences in caste-biased gene content between V. heteropterus and C. cumulans, even though they share similar diet and biological niche. Among the most DEGs, we highlighted those involved in caste differentiation and cellulose digestion, which are attractive targets for studying more efficient technologies for termite control, biomass digestion and other biotechnological applications.

Hormonal control and target genes of ftz-f1 expression in the honeybee Apis mellifera: a positive loop linking juvenile hormone, ftz-f1, and vitellogenin.

Ftz-f1 is an orphan member of the nuclear hormone receptor superfamily. A 20-hydroxyecdysone pulse allows ftz-f1 gene expression, which then regulates the activity of downstream genes involved in major developmental progression events. In honeybees, the expression of genes like vitellogenin (vg), prophenoloxidase and juvenile hormone-esterase during late pharate-adult development is known to be hormonally controlled in both queens and workers by increasing juvenile hormone (JH) titres in the presence of declining levels of ecdysteroids. Since Ftz-f1 is known for mediating intracellular JH signalling, we hypothesized that ftz-f1 could mediate JH action during the pharate-adult development of honeybees, thus controlling the expression of these genes. Here, we show that ftz-f1 has caste-specific transcription profiles during this developmental period, with a peak coinciding with the increase in JH titre, and that its expression is upregulated by JH and downregulated by ecdysteroids. RNAi-mediated knock down of ftz-f1 showed that the expression of genes essential for adult development (e.g. vg and cuticular genes) depends on ftz-f1 expression. Finally, a double-repressor hypothesis-inspired vg gene knock-down experiment suggests the existence of a positive molecular loop between JH, ftz-f1 and vg.

MicroRNA signatures characterizing caste-independent ovarian activity in queen and worker honeybees (Apis mellifera L.).

Queen and worker honeybees differ profoundly in reproductive capacity. The queen of this complex society, with 200 highly active ovarioles in each ovary, is the fertile caste, whereas the workers have approximately 20 ovarioles as a result of receiving a different diet during larval development. In a regular queenright colony, the workers have inactive ovaries and do not reproduce. However, if the queen is sensed to be absent, some of the workers activate their ovaries, producing viable haploid eggs that develop into males. Here, a deep-sequenced ovary transcriptome library of reproductive workers was used as supporting data to assess the dynamic expression of the regulatory molecules and microRNAs (miRNAs) of reproductive and nonreproductive honeybee females. In this library, most of the differentially expressed miRNAs are related to ovary physiology or oogenesis. When we quantified the dynamic expression of 19 miRNAs in the active and inactive worker ovaries and compared their expression in the ovaries of virgin and mated queens, we noted that some miRNAs (miR-1, miR-31a, miR-13b, miR-125, let-7 RNA, miR-100, miR-276, miR-12, miR-263a, miR-306, miR-317, miR-92a and miR-9a) could be used to identify reproductive and nonreproductive statuses independent of caste. Furthermore, integrative gene networks suggested that some candidate miRNAs function in the process of ovary activation in worker bees.

Organization, evolution and transcriptional profile of hexamerin genes of the parasitic wasp Nasonia vitripennis (Hymenoptera: Pteromalidae).

Hexamerins and prophenoloxidases (PPOs) proteins are members of the arthropod-haemocyanin superfamily. In contrast to haemocyanin and PPO, hexamerins do not bind oxygen, but mainly play a role as storage proteins that supply amino acids for insect metamorphosis. We identified seven genes encoding hexamerins, three encoding PPOs, and one hexamerin pseudogene in the genome of the parasitoid wasp Nasonia vitripennis. A phylogenetic analysis of hexamerins and PPOs from this wasp and related proteins from other insect orders suggests an essentially order-specific radiation of hexamerins. Temporal and spatial transcriptional profiles of N. vitripennis hexamerins suggest that they have physiological functions other than metamorphosis, which are arguably coupled with its lifestyle.

Genetic polymorphisms in TLR4, CR1 and Duffy genes are not associated with malaria resistance in patients from Baixo Amazonas region, Brazil.

The main purpose of this research was to analyze the relation of the genetic polymorphisms frequently expressed by antigen-presenting cells, erythrocytes and malaria susceptibility/resistance with the human malaria infection cases. The sample used consisted of 23 Plasmodium vivax (Pv)- and P. falciparum (Pf)-infected patients, and 21 healthy individuals as a control group, from the Baixo Amazonas population in Pará, Brazil. The Asp299Gly polymorphisms in the Toll-like receptor 4 (TLR4), and Gly42Asp, Arg89Cys, Ala100Thr, and T-33C in the Duffy gene (FY) were analyzed by restriction fragment length polymorphism-polymerase chain reaction. The Lys1590Glu and Arg1601Gly polymorphisms in the complement receptor type 1 (CR1) were analyzed by DNA sequencing. According to the results obtained and statistical analysis considering a significance level or alpha = 0.01, we conclude that the low heterozygote frequency (2.27%) for the Asp299Gly mutation, detected in the TLR4 gene, is not related to the Pv and Pf infections in the patients analyzed. Also, the promoter region GATA-1 analysis of the FY gene in the Pv-infected patients showed that the heterozygote frequency for the T-33C mutation (11.36% of the infected patients and 20.45% of the control patients) is not related to infection resistance. Regarding the CR1 gene, the observed heterozygote frequency (9.09%) for the Arg1601Gly mutation in Pf-infected patients when compared to heterozygote frequency in the control group (18.18%) suggests that there is no correlation with infection resistance.

Caste development and reproduction: a genome-wide analysis of hallmarks of insect eusociality.

The honey bee queen and worker castes are a model system for developmental plasticity. We used established expressed sequence tag information for a Gene Ontology based annotation of genes that are differentially expressed during caste development. Metabolic regulation emerged as a major theme, with a caste-specific difference in the expression of oxidoreductases vs. hydrolases. Motif searches in upstream regions revealed group-specific motifs, providing an entry point to cis-regulatory network studies on caste genes. For genes putatively involved in reproduction, meiosis-associated factors came out as highly conserved, whereas some determinants of embryonic axes either do not have clear orthologs (bag of marbles, gurken, torso), or appear to be lacking (trunk) in the bee genome. Our results are the outcome of a first genome-based initiative to provide an annotated framework for trends in gene regulation during female caste differentiation (representing developmental plasticity) and reproduction.