Sexy ways: the methodical approaches to study plant sex chromosomes.

Sex chromosomes have evolved in many plant species with separate sexes. Current plant research is shifting from examining the structure of sex chromosomes to exploring their functional aspects. New studies are progressively unveiling the specific genetic and epigenetic mechanisms responsible for shaping distinct sexes in plants. While the fundamental methods of molecular biology and genomics are generally employed for the analysis of sex chromosomes, it is often necessary to modify classical procedures not only to simplify and expedite analyses but sometimes to make them possible at all. In this review, we demonstrate how, at the level of structural and functional genetics, cytogenetics, and bioinformatics, it is essential to adapt established procedures for sex chromosome analysis.

Transcriptome analysis of thermomorphogenesis in ovules and during early seed development in Brassica napus.

BACKGROUND: Plant sexual reproduction is highly sensitive to elevated ambient temperatures, impacting seed development and production. We previously phenotyped this effect on three rapeseed cultivars (DH12075, Topas DH4079, and Westar). This work describes the transcriptional response associated with the phenotypic changes induced by heat stress during early seed development in Brassica napus. RESULTS: We compared the differential transcriptional response in unfertilized ovules and seeds bearing embryos at 8-cell and globular developmental stages of the three cultivars exposed to high temperatures. We identified that all tissues and cultivars shared a common transcriptional response with the upregulation of genes linked to heat stress, protein folding and binding to heat shock proteins, and the downregulation of cell metabolism. The comparative analysis identified an enrichment for a response to reactive oxygen species (ROS) in the heat-tolerant cultivar Topas, correlating with the phenotypic changes. The highest heat-induced transcriptional response in Topas seeds was detected for genes encoding various peroxidases, temperature-induced lipocalin (TIL1), or protein SAG21/LEA5. On the contrary, the transcriptional response in the two heat-sensitive cultivars, DH12075 and Westar, was characterized by heat-induced cellular damages with the upregulation of genes involved in the photosynthesis and plant hormone signaling pathways. Particularly, the TIFY/JAZ genes involved in jasmonate signaling were induced by stress, specifically in ovules of heat-sensitive cultivars. Using a weighted gene co-expression network analysis (WGCNA), we identified key modules and hub genes involved in the heat stress response in studied tissues of either heat-tolerant or sensitive cultivars. CONCLUSIONS: Our transcriptional analysis complements a previous phenotyping analysis by characterizing the growth response to elevated temperatures during early seed development and reveals the molecular mechanisms underlying the phenotypic response. The results demonstrated that response to ROS, seed photosynthesis, and hormonal regulation might be the critical factors for stress tolerance in oilseed rape.

Telomeric retrotransposons show propensity to form G-quadruplexes in various eukaryotic species.

BACKGROUND: Canonical telomeres (telomerase-synthetised) are readily forming G-quadruplexes (G4) on the G-rich strand. However, there are examples of non-canonical telomeres among eukaryotes where telomeric tandem repeats are invaded by specific retrotransposons. Drosophila melanogaster represents an extreme example with telomeres composed solely by three retrotransposons-Het-A, TAHRE and TART (HTT). Even though non-canonical telomeres often show strand biased G-distribution, the evidence for the G4-forming potential is limited. RESULTS: Using circular dichroism spectroscopy and UV absorption melting assay we have verified in vitro G4-formation in the HTT elements of D. melanogaster. Namely 3 in Het-A, 8 in TART and 2 in TAHRE. All the G4s are asymmetrically distributed as in canonical telomeres. Bioinformatic analysis showed that asymmetric distribution of potential quadruplex sequences (PQS) is common in telomeric retrotransposons in other Drosophila species. Most of the PQS are located in the gag gene where PQS density correlates with higher DNA sequence conservation and codon selection favoring G4-forming potential. The importance of G4s in non-canonical telomeres is further supported by analysis of telomere-associated retrotransposons from various eukaryotic species including green algae, Diplomonadida, fungi, insects and vertebrates. Virtually all analyzed telomere-associated retrotransposons contained PQS, frequently with asymmetric strand distribution. Comparison with non-telomeric elements showed independent selection of PQS-rich elements from four distinct LINE clades. CONCLUSION: Our findings of strand-biased G4-forming motifs in telomere-associated retrotransposons from various eukaryotic species support the G4-formation as one of the prerequisites for the recruitment of specific retrotransposons to chromosome ends and call for further experimental studies.

Unique peptidic agonists of a juvenile hormone receptor with species-specific effects on insect development and reproduction.

Juvenile hormones (JHs) control insect metamorphosis and reproduction. JHs act through a receptor complex consisting of methoprene-tolerant (Met) and taiman (Tai) proteins to induce transcription of specific genes. Among chemically diverse synthetic JH mimics (juvenoids), some of which serve as insecticides, unique peptidic juvenoids stand out as being highly potent yet exquisitely selective to a specific family of true bugs. Their mode of action is unknown. Here we demonstrate that, like established JH receptor agonists, peptidic juvenoids act upon the JHR Met to halt metamorphosis in larvae of the linden bug, Pyrrhocoris apterus. Peptidic juvenoids induced ligand-dependent dimerization between Met and Tai proteins from P. apterus but, consistent with their selectivity, not from other insects. A cell-based split-luciferase system revealed that the Met-Tai complex assembled within minutes of agonist presence. To explore the potential of juvenoid peptides, we synthesized 120 new derivatives and tested them in Met-Tai interaction assays. While many substituents led to loss of activity, improved derivatives active at sub-nanomolar range outperformed hitherto existing peptidic and classical juvenoids including fenoxycarb. Their potency in inducing Met-Tai interaction corresponded with the capacity to block metamorphosis in P. apterus larvae and to stimulate oogenesis in reproductively arrested adult females. Molecular modeling demonstrated that the high potency correlates with high affinity. This is a result of malleability of the ligand-binding pocket of P. apterus Met that allows larger peptidic ligands to maximize their contact surface. Our data establish peptidic juvenoids as highly potent and species-selective novel JHR agonists.

HiC-TE: a computational pipeline for Hi-C data analysis to study the role of repeat family interactions in the genome 3D organization.

MOTIVATION: The role of repetitive DNA in the 3D organization of the interphase nucleus is a subject of intensive study. In studies of 3D nucleus organization, mutual contacts of various loci can be identified by Hi-C sequencing. Typical analyses use binning of read pairs by location to reduce noise. We use binning by repeat families instead to make similar conclusions about repeat regions. RESULTS: To achieve this, we combined Hi-C data, reference genome data and tools for repeat analysis into a Nextflow pipeline identifying and quantifying the contacts of specific repeat families. As an output, our pipeline produces heatmaps showing contact frequency and circular diagrams visualizing repeat contact localization. Using our pipeline with tomato data, we revealed the preferential homotypic interactions of ribosomal DNA, centromeric satellites and some LTR retrotransposon families and, as expected, little contact between organellar and nuclear DNA elements. While the pipeline can be applied to any eukaryotic genome, results in plants provide better coverage, since the built-in TE-greedy-nester software only detects tandems and LTR retrotransposons. Other repeats can be fed via GFF3 files. This pipeline represents a novel and reproducible way to analyze the role of repetitive elements in the 3D organization of genomes. AVAILABILITY AND IMPLEMENTATION: https://gitlab.fi.muni.cz/lexa/hic-te/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Nucleic acids movement and its relation to genome dynamics of repetitive DNA: Is cellular and intercellular movement of DNA and RNA molecules related to the evolutionary dynamic genome components?: Is cellular and intercellular movement of DNA and RNA molecules related to the evolutionary dynamic genome components?

There is growing evidence of evolutionary genome plasticity. The evolution of repetitive DNA elements, the major components of most eukaryotic genomes, involves the amplification of various classes of mobile genetic elements, the expansion of satellite DNA, the transfer of fragments or entire organellar genomes and may have connections with viruses. In addition to various repetitive DNA elements, a plethora of large and small RNAs migrate within and between cells during individual development as well as during evolution and contribute to changes of genome structure and function. Such migration of DNA and RNA molecules often results in horizontal gene transfer, thus shaping the whole genomic network of interconnected species. Here, we propose that a high evolutionary dynamism of repetitive genome components is often related to the migration/movement of DNA or RNA molecules. We speculate that the cytoplasm is probably an ideal compartment for such evolutionary experiments.

Quadruplex-Forming Motif Inserted into 3'UTR of Ty1his3-AI Retrotransposon Inhibits Retrotransposition in Yeast.

Guanine quadruplexes (G4s) serve as regulators of replication, recombination and gene expression. G4 motifs have been recently identified in LTR retrotransposons, but their role in the retrotransposon life-cycle is yet to be understood. Therefore, we inserted G4s into the 3'UTR of Ty1his3-AI retrotransposon and measured the frequency of retrotransposition in yeast strains BY4741, Y00509 (without Pif1 helicase) and with G4-stabilization by N-methyl mesoporphyrin IX (NMM) treatment. We evaluated the impact of G4s on mRNA levels by RT-qPCR and products of reverse transcription by Southern blot analysis. We found that the presence of G4 inhibited Ty1his3-AI retrotransposition. The effect was stronger when G4s were on a transcription template strand which leads to reverse transcription interruption. Both NMM and Pif1p deficiency reduced the retrotransposition irrespective of the presence of a G4 motif in the Ty1his3-AI element. Quantity of mRNA and products of reverse transcription did not fully explain the impact of G4s on Ty1his3-AI retrotransposition indicating that G4s probably affect some other steps of the retrotransposon life-cycle (e.g., translation, VLP formation, integration). Our results suggest that G4 DNA conformation can tune the activity of mobile genetic elements that in turn contribute to shaping the eukaryotic genomes.

Long-lived termite kings and queens activate telomerase in somatic organs.

Kings and queens of termites, like queens of other advanced eusocial insects, are endowed with admirable longevity, which dramatically exceeds the life expectancies of their non-reproducing nest-mates and related solitary insects. In the quest to find the mechanisms underlying the longevity of termite reproductives, we focused on somatic maintenance mediated by telomerase. This ribonucleoprotein is well established for pro-longevity functions in vertebrates, thanks primarily to its ability of telomere extension. However, its participation in lifespan regulation of insects, including the eusocial taxa, remains understudied. Here, we report a conspicuous increase of telomerase abundance and catalytic activity in the somatic organs of primary and secondary reproductives of the termite Prorhinotermes simplex and confirm a similar pattern in two other termite species. These observations stand in contrast with the telomerase downregulation characteristic for most adult somatic tissues in vertebrates and also in solitary insects and non-reproducing castes of termites. At the same time, we did not observe caste-specific differences in telomere lengths that might explain the differential longevity of termite castes. We conclude that although the telomerase activation in termite reproductives is in line with the broadly assumed association between telomerase and longevity, its direct phenotypic impact remains to be elucidated.

TE-greedy-nester: structure-based detection of LTR retrotransposons and their nesting.

MOTIVATION: Transposable elements (TEs) in eukaryotes often get inserted into one another, forming sequences that become a complex mixture of full-length elements and their fragments. The reconstruction of full-length elements and the order in which they have been inserted is important for genome and transposon evolution studies. However, the accumulation of mutations and genome rearrangements over evolutionary time makes this process error-prone and decreases the efficiency of software aiming to recover all nested full-length TEs. RESULTS: We created software that uses a greedy recursive algorithm to mine increasingly fragmented copies of full-length LTR retrotransposons in assembled genomes and other sequence data. The software called TE-greedy-nester considers not only sequence similarity but also the structure of elements. This new tool was tested on a set of natural and synthetic sequences and its accuracy was compared to similar software. We found TE-greedy-nester to be superior in a number of parameters, namely computation time and full-length TE recovery in highly nested regions. AVAILABILITY AND IMPLEMENTATION: http://gitlab.fi.muni.cz/lexa/nested. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

What Can Long Terminal Repeats Tell Us About the Age of LTR Retrotransposons, Gene Conversion and Ectopic Recombination?

LTR retrotransposons constitute a significant part of plant genomes and their evolutionary dynamics play an important role in genome size changes. Current methods of LTR retrotransposon age estimation are based only on LTR (long terminal repeat) divergence. This has prompted us to analyze sequence similarity of LTRs in 25,144 LTR retrotransposons from fifteen plant species as well as formation of solo LTRs. We found that approximately one fourth of nested retrotransposons showed a higher LTR divergence than the pre-existing retrotransposons into which they had been inserted. Moreover, LTR similarity was correlated with LTR length. We propose that gene conversion can contribute to this phenomenon. Gene conversion prediction in LTRs showed potential converted regions in 25% of LTR pairs. Gene conversion was higher in species with smaller genomes while the proportion of solo LTRs did not change with genome size in analyzed species. The negative correlation between the extent of gene conversion and the abundance of solo LTRs suggests interference between gene conversion and ectopic recombination. Since such phenomena limit the traditional methods of LTR retrotransposon age estimation, we recommend an improved approach based on the exclusion of regions affected by gene conversion.

Complex Evolution of Insect Insulin Receptors and Homologous Decoy Receptors, and Functional Significance of Their Multiplicity.

Evidence accumulates that the functional plasticity of insulin and insulin-like growth factor signaling in insects could spring, among others, from the multiplicity of insulin receptors (InRs). Their multiple variants may be implemented in the control of insect polyphenism, such as wing or caste polyphenism. Here, we present a comprehensive phylogenetic analysis of insect InR sequences in 118 species from 23 orders and investigate the role of three InRs identified in the linden bug, Pyrrhocoris apterus, in wing polymorphism control. We identified two gene clusters (Clusters I and II) resulting from an ancestral duplication in a late ancestor of winged insects, which remained conserved in most lineages, only in some of them being subject to further duplications or losses. One remarkable yet neglected feature of InR evolution is the loss of the tyrosine kinase catalytic domain, giving rise to decoys of InR in both clusters. Within the Cluster I, we confirmed the presence of the secreted decoy of insulin receptor in all studied Muscomorpha. More importantly, we described a new tyrosine kinase-less gene (DR2) in the Cluster II, conserved in apical Holometabola for ∼300 My. We differentially silenced the three P. apterus InRs and confirmed their participation in wing polymorphism control. We observed a pattern of Cluster I and Cluster II InRs impact on wing development, which differed from that postulated in planthoppers, suggesting an independent establishment of insulin/insulin-like growth factor signaling control over wing development, leading to idiosyncrasies in the co-option of multiple InRs in polyphenism control in different taxa.

Correction to: Nested plant LTR retrotransposons target specific regions of other elements, while all LTR retrotransposons often target palindromes and nucleosome-occupied regions: in silico study.

[This corrects the article DOI: 10.1186/s13100-019-0186-z.].

Nested plant LTR retrotransposons target specific regions of other elements, while all LTR retrotransposons often target palindromes and nucleosome-occupied regions: in silico study.

BACKGROUND: Nesting is common in LTR retrotransposons, especially in large genomes containing a high number of elements. RESULTS: We analyzed 12 plant genomes and obtained 1491 pairs of nested and original (pre-existing) LTR retrotransposons. We systematically analyzed mutual nesting of individual LTR retrotransposons and found that certain families, more often belonging to the Ty3/gypsy than Ty1/copia superfamilies, showed a higher nesting frequency as well as a higher preference for older copies of the same family ("autoinsertions"). Nested LTR retrotransposons were preferentially located in the 3'UTR of other LTR retrotransposons, while coding and regulatory regions (LTRs) are not commonly targeted. Insertions displayed a weak preference for palindromes and were associated with a strong positional pattern of higher predicted nucleosome occupancy. Deviation from randomness in target site choice was also found in 13,983 non-nested plant LTR retrotransposons. CONCLUSIONS: We reveal that nesting of LTR retrotransposons is not random. Integration is correlated with sequence composition, secondary structure and the chromatin environment. Insertion into retrotransposon positions with a low negative impact on family fitness supports the concept of the genome being viewed as an ecosystem of various elements.

Changes in vitellogenin expression caused by nematodal and fungal infections in insects.

This study examined the expression and role of vitellogenin (Vg) in the body of the firebug Pyrrhocoris apterus (Heteroptera, Insecta) during infection elicited by two entomopathogenic organisms, the nematode Steinernema carpocapsae and the fungus Isaria fumosorosea Infection by S. carpocapsae significantly upregulated Vg mRNA expression in the male body. The corresponding increase in Vg protein expression was also confirmed by electrophoretic and immunoblotting analyses. Remarkably, in females, the opposite tendency was noted. Nematodal infection significantly reduced both Vg mRNA and Vg protein expression levels in fat body and hemolymph, respectively. We speculate that infection of reproductive females reduces Vg expression to a level that is still sufficient for defense, but is insufficient for reproduction. This circumstance reduces energy expenditure and helps the individual to cope with the infection. Importantly, purified Vg significantly inhibited growth of Xenorhabdus spp., an entomotoxic bacteria isolated from S. carpocapsae. However, the effect of Vg against I. fumosorosea was not so obvious. The fungus significantly stimulated Vg gene expression in males; however, a similar increase was not recapitulated at the protein level. Nevertheless, in females, both mRNA and protein Vg levels were significantly reduced after the fungal infection. The obtained data demonstrate that Vg is probably an important defense protein, possibly with a specific activity. This considerably expands the known spectrum of Vg functions, as its primary role was thought to be limited to regulating egg development in the female body.

Exquisite ligand stereoselectivity of a Drosophila juvenile hormone receptor contrasts with its broad agonist repertoire.

The sesquiterpenoid juvenile hormone (JH) is vital to insect development and reproduction. Intracellular JH receptors have recently been established as basic helix-loop-helix transcription factor (bHLH)/PAS proteins in Drosophila melanogaster known as germ cell-expressed (Gce) and its duplicate paralog, methoprene-tolerant (Met). Upon binding JH, Gce/Met activates its target genes. Insects possess multiple native JH homologs whose molecular activities remain unexplored, and diverse synthetic compounds including insecticides exert JH-like effects. How the JH receptor recognizes its ligands is unknown. To determine which structural features define an active JH receptor agonist, we tested several native JHs and their nonnative geometric and optical isomers for the ability to bind the Drosophila JH receptor Gce, to induce Gce-dependent transcription, and to affect the development of the fly. Our results revealed high ligand stereoselectivity of the receptor. The geometry of the JH skeleton, dictated by two stereogenic double bonds, was the most critical feature followed by the presence of an epoxide moiety at a terminal position. The optical isomerism at carbon C11 proved less important even though Gce preferentially bound a natural JH enantiomer. The results of receptor-ligand-binding and cell-based gene activation assays tightly correlated with the ability of different geometric JH isomers to induce gene expression and morphogenetic effects in the developing insects. Molecular modeling supported the requirement for the proper double-bond geometry of JH, which appears to be its major selective mechanism. The strict stereoselectivity of Gce toward the natural hormone contrasts with the high potency of synthetic Gce agonists of disparate chemistries.

Metabolomic and transcriptomic data on major metabolic/biosynthetic pathways in workers and soldiers of the termite Prorhinotermes simplex (Isoptera: Rhinotermitidae) and chemical synthesis of intermediates of defensive (E)-nitropentadec-1-ene biosynthesis.

Production of nitro compounds has only seldom been recorded in arthropods. The aliphatic nitroalkene (E)-nitropentadec-1-ene (NPD), identified in soldiers of the termite genus Prorhinotermes, was the first case documented in insects in early seventies. Yet, the biosynthetic origin of NPD has long remained unknown. We previously proposed that NPD arises through the condensation of amino acids glycine and/or l-serine with tetradecanoic acid along a biosynthetic pathway analogous to the formation of sphingolipids. Here, we provide a metabolomics and transcriptomic data of the Prorhinotermes simplex termite workers and soldiers. Data are related to NPD biosynthesis in P. simplex soldiers. Original metabolomics data were deposited in MetaboLights metabolomics database and are become publicly available after publishing the original article. Additionally, chemical synthesis of biosynthetic intermediates of NPD in nonlabeled and stable labeled forms are reported. Data extend our poor knowledge of arthropod metabolome and transcriptome and would be useful for comparative study in termites or other arthropods. The data were used for de-replication of NPD biosynthesis and published separately (Jirosová et al., 2017) [1].

Co-option of the sphingolipid metabolism for the production of nitroalkene defensive chemicals in termite soldiers.

The aliphatic nitroalkene (E)-1-nitropentadec-1-ene (NPD), reported in early seventies in soldiers of the termite genus Prorhinotermes, was the first documented nitro compound produced by insects. Yet, its biosynthetic origin has long remained unknown. Here, we investigated in detail the biosynthesis of NPD in P. simplex soldiers. First, we track the dynamics in major metabolic pathways during soldier ontogeny, with emphasis on likely NPD precursors and intermediates. Second, we propose a hypothesis of NPD formation and verify its individual steps using in vivo incubations of putative precursors and intermediates. Third, we use a de novo assembled RNA-Seq profiles of workers and soldiers to identify putative enzymes underlying NPD formation. And fourth, we describe the caste- and age-specific expression dynamics of candidate initial genes of the proposed biosynthetic pathway. Our observations provide a strong support to the following biosynthetic scenario of NPD formation, representing an analogy of the sphingolipid pathway starting with the condensation of tetradecanoic acid with l-serine and leading to the formation of a C16 sphinganine. The C16 sphinganine is then oxidized at the terminal carbon to give rise to 2-amino-3-hydroxyhexadecanoic acid, further oxidized to 2-amino-3-oxohexadecanoic acid. Subsequent decarboxylation yields 1-aminopentadecan-2-one, which then proceeds through six-electron oxidation of the amino moiety to give rise to 1-nitropentadecan-2-one. Keto group reduction and hydroxyl moiety elimination lead to NPD. The proposed biosynthetic sequence has been constructed from age-related quantitative dynamics of individual intermediates and confirmed by the detection of labeled products downstream of the administered labeled intermediates. Comparative RNA-Seq analyses followed by qRT-PCR validation identified orthologs of serine palmitoyltransferase and 3-ketodihydrosphingosine reductase genes as highly expressed in the NPD production site, i.e. the frontal gland of soldiers. A dramatic onset of expression of the two genes in the first days of soldier's life coincides with the start of NPD biosynthesis, giving further support to the proposed biosynthetic hypothesis.

Adipokinetic hormones control amylase activity in the cockroach (Periplaneta americana) gut.

This study examined the biochemical characteristics of α-amylase and hormonal (adipokinetic hormone: AKH) stimulation of α-amylase activity in the cockroach (Periplaneta americana) midgut. We applied two AKHs in vivo and in vitro, then measured resultant amylase activity and gene expression, as well as the expression of AKH receptor (AKHR). The results revealed that optimal amylase activity is characterized by the following: pH: 5.7, temperature: 38.4 °C, Km (Michaelis-Menten constant): 2.54 mg starch/mL, and Vmax (maximum reaction velocity): 0.185 µmol maltose/mL/min. In vivo application of AKHs resulted in significant increase of amylase activity: by two-fold in the gastric caeca and 4-7 fold in the rest of the midgut. In vitro experiments supported results seen in vivo: a 24-h incubation with the hormones resulted in the increase of amylase activity by 1.4 times in the caeca and 4-9 times in the midgut. Further, gene expression analyses reveal that AKHR is expressed in both the caeca and the rest of the midgut, although expression levels in the former were 23 times higher than levels in the latter. A similar pattern was found for the amylase (AMY) gene. Hormonal treatment did not affect the expression of either gene. This study is the first to provide evidence indicating direct AKH stimulation of digestive enzyme activity in the insect midgut, supported by specific AKHR gene expression in this organ.

Gene Expression Dynamics in Major Endocrine Regulatory Pathways along the Transition from Solitary to Social Life in a Bumblebee, Bombus terrestris.

Understanding the social evolution leading to insect eusociality requires, among other, a detailed insight into endocrine regulatory mechanisms that have been co-opted from solitary ancestors to play new roles in the complex life histories of eusocial species. Bumblebees represent well-suited models of a relatively primitive social organization standing on the mid-way to highly advanced eusociality and their queens undergo both, a solitary and a social phase, separated by winter diapause. In the present paper, we characterize the gene expression levels of major endocrine regulatory pathways across tissues, sexes, and life-stages of the buff-tailed bumblebee, Bombus terrestris, with special emphasis on critical stages of the queen's transition from solitary to social life. We focused on fundamental genes of three pathways: (1) Forkhead box protein O and insulin/insulin-like signaling, (2) Juvenile hormone (JH) signaling, and (3) Adipokinetic hormone signaling. Virgin queens were distinguished by higher expression of forkhead box protein O and downregulated insulin-like peptides and JH signaling, indicated by low expression of methyl farnesoate epoxidase (MFE) and transcription factor Krüppel homolog 1 (Kr-h1). Diapausing queens showed the expected downregulation of JH signaling in terms of low MFE and vitellogenin (Vg) expressions, but an unexpectedly high expression of Kr-h1. By contrast, reproducing queens revealed an upregulation of MFE and Vg together with insulin signaling. Surprisingly, the insulin growth factor 1 (IGF-1) turned out to be a queen-specific hormone. Workers exhibited an expression pattern of MFE and Vg similar to that of reproducing queens. Males were characterized by high Kr-h1 expression and low Vg level. The tissue comparison unveiled an unexpected resemblance between the fat body and hypopharyngeal glands across all investigated genes, sexes, and life stages.

Characterization and expression analysis of adipokinetic hormone and its receptor in eusocial aphid Pseudoregma bambucicola.

Aphids display an extraordinary phenotypic plasticity ranging from widespread reproductive and wing polyphenisms to the occurrence of sterile or subfertile soldier morphs restricted to eusocial species of the subfamilies Eriosomatinae and Hormaphidinae. Individual morphs are specialized by their behavior, anatomy, and physiology to perform different roles in aphid societies at different stages of the life cycle. The capacity of the insects to cope with environmental stressors is under the control of a group of neuropeptides of the adipokinetic hormone/red pigment-concentrating hormone family (AKH/RPCH) that bind to a specific receptor (AKHR). Here, we describe the molecular characteristics of AKH and AKHR in the eusocial aphid Pseudoregma bambucicola. The sequence of the bioactive AKH decapeptide and the intron position in P. bambucicola AKH preprohormone were found to be identical to those in a phylogenetically distant aphid Dreyfusia spp. (Adelgidae). We detected four transcript variants of AKHR that are translated into three protein isoforms. Further, we analyzed AKH/AKHR expression in different tissues and insects of different castes. In wingless females, a remarkable amount of AKH mRNA was only expressed in the heads. In contrast, AKHR transcript levels increased in the order gut