Dynamic transcriptome landscape of Asian domestic honeybee (Apis cerana) embryonic development revealed by high-quality RNA sequencing

Microplastics enhanced the allelopathy of pyrogallol on toxic Microcystis with additional risks: Microcystins release and greenhouse gases emissions

Cyanobacteria blooms (CBs) caused by eutrophication pose a global concern, especially Microcystis aeruginosa (M. aeruginosa), which could release harmful microcystins (MCs). The impact of microplastics (MPs) on allelopathy in freshwater environments is not well understood. This study examined the joint effect of adding polystyrene (PS-MPs) as representative MPs and two concentrations (2 and 8 mg/L) of pyrogallol (PYR) on the allelopathy of M. aeruginosa. The results showed that the addition of PS-MPs intensified the inhibitory effect of 8 mg/L PYR on the growth and photosynthesis of M. aeruginosa. After a 7-day incubation period, the cell density decreased to 69.7 %, and the chl-a content decreased to 48 % compared to the condition without PS-MPs (p < 0.05). Although the growth and photosynthesis of toxic Microcystis decreased with the addition of PS-MPs, the addition of PS-MPs significantly resulted in a 3.49-fold increase in intracellular MCs and a 1.10-fold increase in extracellular MCs (p < 0.05). Additionally, the emission rates of greenhouse gases (GHGs) (carbon dioxide, nitrous oxide and methane) increased by 2.66, 2.23 and 2.17-fold, respectively (p < 0.05). In addition, transcriptomic analysis showed that the addition of PS-MPs led to the dysregulation of gene expression related to DNA synthesis, membrane function, enzyme activity, stimulus detection, MCs release and GHGs emissions in M. aeruginosa. PYR and PS-MPs triggered ROS-induced membrane damage and disrupted photosynthesis in algae, leading to increased MCs and GHG emissions. PS-MPs accumulation exacerbated this issue by impeding light absorption and membrane function, further heightening the release of MCs and GHGs emissions. Therefore, PS-MPs exhibited a synergistic effect with PYR in inhibiting the growth and photosynthesis of M. aeruginosa, resulting in additional risks such as MCs release and GHGs emissions. These results provide valuable insights for the ecological risk assessment and control of algae bloom in freshwater ecosystems.

Expression Profile, Regulatory Network, and Putative Role of microRNAs in the Developmental Process of Asian Honey Bee Larval Guts

MiRNAs, as a kind of key regulators in gene expression, play vital roles in numerous life activities from cellular proliferation and differentiation to development and immunity. However, little is known about the regulatory manner of miRNAs in the development of Asian honey bee (Apis cerana) guts. Here, on basis of our previously gained high-quality transcriptome data, transcriptome-wide identification of miRNAs in the larval guts of Apis cerana cerana was conducted, followed by investigation of the miRNAs' differential expression profile during the gut development. In addition to the regulatory network, the potential function of differentially expressed miRNAs (DEmiRNAs) was further analyzed. In total, 330, 351, and 321 miRNAs were identified in the 4-, 5-, and 6-day-old larval guts, respectively; among these, 257 miRNAs were shared, while 38, 51, and 36 ones were specifically expressed. Sequences of six miRNAs were confirmed by stem-loop RT-PCR and Sanger sequencing. Additionally, in the "Ac4 vs. Ac5" comparison group, there were seven up-regulated and eight down-regulated miRNAs; these DEmiRNAs could target 5041 mRNAs, involving a series of GO terms and KEGG pathways associated with growth and development, such as cellular process, cell part, Wnt, and Hippo. Comparatively, four up-regulated and six down-regulated miRNAs detected in the "Ac5 vs. Ac6" comparison group and the targets were associated with diverse development-related terms and pathways, including cell, organelle, Notch and Wnt. Intriguingly, it was noticed that miR-6001-y presented a continuous up-regulation trend across the developmental process of larval guts, implying that miR-6001-y may be a potential essential modulator in the development process of larval guts. Further investigation indicated that 43 targets in the "Ac4 vs. Ac5" comparison group and 31 targets in the "Ac5 vs. Ac6" comparison group were engaged in several crucial development-associated signaling pathways such as Wnt, Hippo, and Notch. Ultimately, the expression trends of five randomly selected DEmiRNAs were verified using RT-qPCR. These results demonstrated that dynamic expression and structural alteration of miRNAs were accompanied by the development of A. c. cerana larval guts, and DEmiRNAs were likely to participate in the modulation of growth as well as development of larval guts by affecting several critical pathways via regulation of the expression of target genes. Our data offer a basis for elucidating the developmental mechanism underlying Asian honey bee larval guts.

A comparative transcriptome analysis of the head of 1 and 9 days old worker honeybees (Apis mellifera)

The role of bees in the environment, economic, biodiversity and pharmaceutical industries is due to its social behavior, which is oriented from the brain and hypopharyngeal gland that is the center of royal jelly (RJ) production. Limited studies have been performed on the head gene expression profile at the RJ production stage. The aim of this study was to compare the gene expressions in 9 and 1-day-old (DO) honeybee workers in order to achieve better understanding about head gene expression pattern. After sequencing of RNAs, transcriptome and their networks were compared. The head expression profile undergoes various changes. 1662 gene transcripts had differential expressions which 1125 and 537 were up and down regulated, respectively, in 9_DO compared with 1_DO honey bees. The day 1th had more significant role in the expression of genes related to RJ production as major RJ protein 1, 2, 3, 5, 6 and 9 encoding genes, but their maximum secretion occurred at day 9th. All process related to hypopharyngeal glands activities as CYP450 gene, fatty acid synthase gene, vitamin B6 metabolism and some of genes involved in fatty acid elongation and degradation process had an upward trend from 1_DO and were age-dependent. By increasing the age, the activity of pathways related to immune system increased for keeping the health of bees against the chemical compound. The expression of aromatic amino acid genes involved in Phenylalanine, tyrosine and tryptophan biosynthesis pathway are essential for early stage of life. In 9_DO honeybees, the energy supplying, reducing stress, protein production and export pathways have a crucial role for support the body development and the social duties. It can be stated that the activity of honeybee head is focused on energy supply instead of storage, while actively trying to improve the level of cell dynamics for increasing the immunity and reducing stress. Results of current study identified key genes of certain behaviors of honeybee workers. Deeper considering of some pathways will be evaluated in future studies.

Dynamics of maternal gene expression in Rhodnius prolixus

The study of developmental processes in Rhodnius prolixus has recently advanced with the sequencing of the genome. In this work, we analyze the maternal gene expression driving oogenesis and early embryogenesis in R. prolixus. We examined the transcriptional profile of mRNAs to establish the genes expressed across the ovary, unfertilized eggs and different embryonic stages of R. prolixus until the formation of the germ band anlage (0, 12, 24, and 48 h post egg laying). We identified 81 putative maternal and ovary-related genes and validated their expression by qRT-PCR. We validate the function of the ortholog gene Bicaudal-D (Rp-BicD) by in situ hybridization and parental RNAi. Consistent with a role in oogenesis and early development of R. prolixus, we show that lack of Rp-BicD does not significantly affect oogenesis but impairs the formation of the blastoderm. Based on our findings, we propose three times of action for maternal genes during oogenesis and embryogenesis in R. prolixus.

Transcriptional Profiles of Diploid Mutant Apis mellifera Embryos after Knockout of csd by CRISPR/Cas9

Simple Summary

In honey bees, males are haploid while females are diploid, leading to a fundamental difference in genetic materials between the sexes. In order to better control the comparison of gene expression between males and females, diploid mutant males were generated by knocking out the sex-determining gene, complementary sex determiner (csd), in fertilized embryos. The diploid mutant drones had male external morphological features, as well as male gonads. RNA sequencing was performed on the diploid mutant embryos and one-day-old larvae. The transcriptome analysis showed that several female-biased genes, such as worker-enriched antennal (Wat), vitellogenin (Vg), and some venom-related genes, were down-regulated in the diploid mutant males. In contrast, some male-biased genes, like takeout and apolipophorin-III-like protein (A4), were up-regulated. Moreover, the co-expression gene networks suggested that csd might interact very closely with fruitless (fru), feminizer (fem) might have connections with hexamerin 70c (hex70c), and transformer-2 (tra2) might play roles with troponin T (TpnT). Foundational information about the differences in the gene expression caused by sex differentiation was provided in this study. It is believed that this study will pave the ground for further research on the different mechanisms between males and females in honey bees.

Abstract

In honey bees, complementary sex determiner (csd) is the primary signal of sex determination. Its allelic composition is heterozygous in females, and hemizygous or homozygous in males. To explore the transcriptome differences after sex differentiation between males and females, with genetic differences excluded, csd in fertilized embryos was knocked out by CRISPR/Cas9. The diploid mutant males at 24 h, 48 h, 72 h, and 96 h after egg laying (AEL) and the mock-treated females derived from the same fertilized queen were investigated through RNA-seq. Mutations were detected in the target sequence in diploid mutants. The diploid mutant drones had typical male morphological characteristics and gonads. Transcriptome analysis showed that several female-biased genes, such as worker-enriched antennal (Wat), vitellogenin (Vg), and some venom-related genes, were down-regulated in the diploid mutant males. In contrast, some male-biased genes, such as takeout and apolipophorin-III-like protein (A4), had higher expressions in the diploid mutant males. Weighted gene co-expression network analysis (WGCNA) indicated that there might be interactions between csd and fruitless (fru), feminizer (fem) and hexamerin 70c (hex70c), transformer-2 (tra2) and troponin T (TpnT). The information provided by this study will benefit further research on the sex dimorphism and development of honey bees and other insects in Hymenoptera.

Egg transcriptome profile responds to maternal virus infection in honey bees, Apis mellifera

Trans-generational disease effects include vertical pathogen transmission but also immune priming to enhance offspring immunity. Accordingly, the survival consequences of maternal virus infection can vary and its molecular consequences during early development are poorly understood. The honey bee queen is long-lived and represents the central hub for vertical virus transmission as the sole reproductive individual in her colony. Even though virus symptoms in queens are mild, viral infection may have severe consequences for the offspring. Thus, transcriptome patterns during early developmental are predicted to respond to maternal virus infection. To test this hypothesis, gene expression patterns were compared among pooled honey bee eggs laid by queens that were either infected with Deformed wing virus (DWV¹), Sacbrood virus (SBV²), both viruses (DWV and SBV), or no virus. Whole transcriptome analyses revealed significant expression differences of a few genes, some of which have hitherto no known function. Despite the paucity of single gene effects, functional enrichment analyses revealed numerous biological processes in the embryos to be affected by virus infection. Effects on several regulatory pathways were consistent with maternal responses to virus infection and correlated with responses to DWV and SBV in honey bee larvae and pupae. Overall, effects on egg transcriptome patterns were specific to each virus and the results of dual-infection samples suggested synergistic effects of DWV and SBV. We interpret our results as consequences of maternal infections. Thus, this first study to document and characterize virus-associated changes in the transcriptome of honey bee eggs represents an important contribution to understanding trans-generational virus effects, although more in-depth studies are needed to understand the detailed mechanisms of how viruses affect honey bee embryos.

An Overview of Embryogenesis: External Morphology and Transcriptome Profiling in the Hemipteran Insect Nilaparvata lugens

During embryogenesis of insects, the morphological and transcriptional changes are important signatures to obtain a better understanding of insect patterning and evolution. The brown planthopper Nilaparvata lugens is a serious insect pest of rice plants, but its embryogenesis has not uncovered. Here, we described embryonic development process of the pest and found it belongs to an intermediate-germ mode. The RNA-seq data from different times (6, 30, 96, and 150 h, after egg laying) of embryogenesis were then analyzed, and a total of 10,895 genes were determined as differentially expressed genes (DEGs) based on pairwise comparisons. Afterward, 1,898 genes, differentially expressed in at least two comparisons of adjacent embryonic stages were divided into 10 clusters using K means cluster analysis (KMCA). Eight-gene modules were established using a weighted gene co-expression network analysis (WGCNA). Gene expression patterns in the different embryonic stages were identified by combining the functional enrichments of the stage-specific clusters and modules, which displayed the expression level and reprogramming of multiple developmental genes during embryogenesis. The "hub" genes at each embryonic stage with possible crucial roles were identified. Notably, we found a "center" set of genes that were related to overall membrane functions and might play important roles in the embryogenesis process. After parental RNAi of the MSTRG.3372, the hub gene, the embryo was observed as abnormal. Furthermore, some homologous genes in classic embryonic development processes and signaling pathways were also involved in embryogenesis of this insect. An improved comprehensive finding of embryogenesis within the N. lugens reveals better information on genetic and genomic studies of embryonic development and might be a potential target for RNAi-based control of this insect pest.

Egg-size plasticity in Apis mellifera: Honey bee queens alter egg size in response to both genetic and environmental factors

Social evolution has led to distinct life-history patterns in social insects, but many colony-level and individual traits, such as egg size, are not sufficiently understood. Thus, a series of experiments was performed to study the effects of genotypes, colony size and colony nutrition on variation in egg size produced by honey bee (Apis mellifera) queens. Queens from different genetic stocks produced significantly different egg sizes under similar environmental conditions, indicating standing genetic variation for egg size that allows for adaptive evolutionary change. Further investigations revealed that eggs produced by queens in large colonies were consistently smaller than eggs produced in small colonies, and queens dynamically adjusted egg size in relation to colony size. Similarly, queens increased egg size in response to food deprivation. These results could not be solely explained by different numbers of eggs produced in the different circumstances but instead seem to reflect an active adjustment of resource allocation by the queen in response to colony conditions. As a result, larger eggs experienced higher subsequent survival than smaller eggs, suggesting that honey bee queens might increase egg size under unfavourable conditions to enhance brood survival and to minimize costly brood care of eggs that fail to successfully develop, and thus conserve energy at the colony level. The extensive plasticity and genetic variation of egg size in honey bees has important implications for understanding life-history evolution in a social context and implies this neglected life-history stage in honey bees may have trans-generational effects.

Comparative Identification of MicroRNAs in Apis cerana cerana Workers' Midguts in Responseto Nosema ceranae Invasion

Here, the expression profiles and differentially expressed miRNAs (DEmiRNAs) in the midguts of Apis cerana cerana workers at 7 d and 10 d post-inoculation (dpi) with N. ceranae were investigated via small RNA sequencing and bioinformatics. Five hundred and twenty nine (529) known miRNAs and 25 novel miRNAs were identified in this study, and the expression of 16 predicted miRNAs was confirmed by Stem-loop RT-PCR. A total of 14 DEmiRNAs were detected in the midgut at 7 dpi, including eight up-regulated and six down-regulated miRNAs, while 12 DEmiRNAs were observed in the midgut at 10 dpi, including nine up-regulated and three down-regulated ones. Additionally, five DEmiRNAs were shared, while nine and seven DEmiRNAs were specifically expressed in midguts at 7 dpi and 10 dpi. Gene ontology analysis suggested some DEmiRNAs and corresponding target mRNAs were involved in various functions including immune system processes and response to stimulus. KEGG pathway analysis shed light on the potential functions of some DEmiRNAs in regulating target mRNAs engaged in material and energy metabolisms, cellular immunity and the humoral immune system. Further investigation demonstrated a complex regulation network between DEmiRNAs and their target mRNAs, with miR-598-y, miR-252-y, miR-92-x and miR-3654-y at the center. Our results can facilitate future exploration of the regulatory roles of miRNAs in host responses to N. ceranae, and provide potential candidates for further investigation of the molecular mechanisms underlying eastern honeybee-microsporidian interactions.

High-Efficiency CRISPR/Cas9-Mediated Gene Editing in Honeybee (Apis mellifera) Embryos

The honeybee (Apis mellifera) is an important insect pollinator of wild flowers and crops, playing critical roles in the global ecosystem. Additionally, the honeybee serves as an ideal social insect model. Therefore, functional studies on honeybee genes are of great interest. However, until now, effective gene manipulation methods have not been available in honeybees. Here, we reported an improved CRISPR/Cas9 gene-editing method by microinjecting sgRNA and Cas9 protein into the region of zygote formation within 2 hr after queen oviposition, which allows one-step generation of biallelic knockout mutants in honeybee with high efficiency. We first targeted the Mrjp1 gene. Two batches of honeybee embryos were collected and injected with Mrjp1 sgRNA and Cas9 protein at the ventral cephalic side and the dorsal posterior side of the embryos, respectively. The gene-editing rate at the ventral cephalic side was 93.3%, which was much higher than that (11.8%) of the dorsal-posterior-side injection. To validate the high efficiency of our honeybee gene-editing system, we targeted another gene, Pax6, and injected Pax6 sgRNA and Cas9 protein at the ventral cephalic side in the third batch. A 100% editing rate was obtained. Sanger sequencing of the TA clones showed that 73.3% (for Mrjp1) and 76.9% (for Pax6) of the edited current-generation embryos were biallelic knockout mutants. These results suggest that the CRISPR/Cas9 method we established permits one-step biallelic knockout of target genes in honeybee embryos, thereby demonstrating an efficient application to functional studies of honeybee genes. It also provides a useful reference to gene editing in other insects with elongated eggs.