Molecular evolution and expression of the CRAL_TRIO protein family in insects

Expression of proteins supporting visual function in heterobranch gastropods

To sense light, animals often utilize mechanisms that rely on visual pigments composed of opsin and retinal. The photon-induced isomerization of 11-cis-retinal to the all-trans configuration triggers phototransduction cascades, resulting in a change in the membrane potential of the photoreceptor. In mollusks, the most abundant opsin in the eye is Gq-coupled rhodopsin (Gq-rhodopsin). The Gq-rhodopsin-based visual pigment is bistable, with the regeneration of 11-cis-retinal occurring in a light-dependent manner without leaving the opsin moiety. 11-cis-retinal is also regenerated by the action of retinochrome in the cell bodies. Retinal binding protein (RALBP) mediates retinal transport between Gq-rhodopsin and retinochrome in the cytoplasm. However, recent studies have identified additional bistable opsins in mollusks, including Opn5 and xenopsin. It is unknown whether these bistable opsins require RALBP and retinochrome for the continuous regeneration of 11-cis-retinal. In the present study, we examined the expression of RALBP and retinochrome in the photoreceptors expressing Opn5 or Xenopsin in the heterobranch gastropods Limax and Peronia. Our findings revealed that retinochrome, but not RALBP, was present in some of the Opn5A-positive brain photosensory neurons of Limax. The ciliary cells in the dorsal eye of Peronia, which express Xenopsin2, lacked both retinochrome and RALBP. Therefore, bistable opsins do not necessarily depend on the RALBP-retinochrome system in a cell. We also examined the expression of other proteins that support visual function, such as β-arrestin, Gq, and Go, in all types of photoreceptors in these animals, and uncovered differences in the molecular composition among the photoreceptors.

Changes in the proteome of Apis mellifera acutely exposed to sublethal dosage of glyphosate and imidacloprid

Uncontrolled use of pesticides has caused a dramatic reduction in the number of pollinators, including bees. Studies on the effects of pesticides on bees have reported effects on both metabolic and neurological levels under chronic exposure. In this study, variations in the differential expression of head and thorax-abdomen proteins in Africanized A. mellifera bees treated acutely with sublethal doses of glyphosate and imidacloprid were studied using a proteomic approach. A total of 92 proteins were detected, 49 of which were differentially expressed compared to those in the control group (47 downregulated and 2 upregulated). Protein interaction networks with differential protein expression ratios suggested that acute exposure of A. mellifera to sublethal doses of glyphosate could cause head damage, which is mainly associated with behavior and metabolism. Simultaneously, imidacloprid can cause damage associated with metabolism as well as, neuronal damage, cellular stress, and impairment of the detoxification system. Regarding the thorax-abdomen fractions, glyphosate could lead to cytoskeleton reorganization and a reduction in defense mechanisms, whereas imidacloprid could affect the coordination and impairment of the oxidative stress response.

Comprehensive transcriptome sequencing of silkworm Midguts: Uncovering extensive isoform diversity and alternative splicing in BmNPV-Sensitive and BmNPV-resistant strains

The silkworm, Bombyx mori, stands out as one of the few economically valuable insects within the realm of model organisms. However, Bombyx mori nucleopolyhedrovirus (BmNPV) poses a significant threat, decreasing the quality and quantity of silkworm cocoons. Over the past few decades, a multitude of researchers has delved into the mechanisms that underlie silkworm resistance to BmNPV, employing diverse methodologies and approaching the problem from various angles. Despite this extensive research, the role of alternative splicing (AS) in the silkworm's response to BmNPV infection has been largely unexplored. This study leveraged both third-generation (Oxford Nanopore Technologies) and second-generation (Illumina) high-throughput sequencing technologies to meticulously identify and analyze AS patterns in the context of BmNPV response, utilizing two distinct silkworm strains-the susceptible strain 306 and the resistant strain NB. Consequently, we identified five crucial genes (Dsclp, LOC692903, LOC101743583, LOC101742498, LOC101743809) that are linked to the response to BmNPV infection through AS and differential expression. Additionally, a thorough comparative analysis was conducted on their diverse transcriptomic expression profiles, including alternative polyadenylation, simple sequence repeats, and transcription factors.

Changes in Gene Expression of Whiteflies, Bemisia tabaci MED Feeding on Tomato Plants Infected by One of the Criniviruses, Tomato Chlorosis Virus through Transcriptome Analysis

Tomato chlorosis virus (ToCV), transmitted by the whitefly, Bemisia tabaci (Gennadius; Hemiptera: Aleyrodidae) has been continuously emerging on tomato plants and causing a significant economic loss throughout China. In the current study, RNA-Seq analysis was used to explore the gene expression profiles of B. tabaci Mediterranean (MED) that fed on both ToCV-infected and -uninfected tomato plants for 6, 12, 24, and 48 hours, respectively. The results revealed that dynamic changes occurred in the gene expressions of whiteflies at different time intervals after they acquired the virus. A total of 1709, 461, 4548, and 1748 differentially expressed genes (DEGs) were identified after a 6, 12, 24, and 48 hours feeding interval for the viral acquisition, respectively. The least number of expressed genes appeared in whiteflies with the 12 hours feeding treatment, and the largest numbers of those found in those with 24 hours feeding treatment. Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that B. tabaci MED responded to ToCV acquisition through altering its nerve system development, fertility, detoxification, glucose metabolism, and immune function before it lost its ability to transmit the virus. The number of DEGs, degree of differential gene expressions, expression level of the same gene, involved biological processes, and metabolic functions in whiteflies post the 12 hours feeding, and viral acquisition were different from those from other three feeding treatments, which could be a significant finding suggesting an effective control of B. tabaci MED should be done less than 12 hours after whiteflies started feeding on ToCV-infected tomatoes. Our results further provided a clarified understanding in how B. tabaci was protected from viral acquisitions through comparison of the differential profile of gene expressions in whiteflies feeding on plants that were infected by semipersistent viruses.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07905-7.

A plant CitPITP1 protein-coding exon sequence serves as a promoter in bacteria

Genetic manipulation of plant genes in prokaryotes has been widely used in molecular biology, but the function of a DNA sequence is far from being fully known. Here, we discovered that a plant protein-coding gene containing the CRAL_TRIO domain serves as a promoter in bacteria. We firstly characterized CitPITP1 from Citrus, which contains the CRAL_TRIO domain, and identified a 64-bp sequence (key64) that is critical for prokaryotic promoter activity. In vitro experiments indicated that the bacterial RNA polymerase subunit RpoD specifically binds to key64. We then expanded our research to fungi, plant and animal species to identify key64-like sequences. Five such prokaryotic promoters were isolated from Amborella, Rice, Arabidopsis and Citrus. Two conserved motifs were identified, and mutation analysis indicated that the nucleotides at positions 7, 29 and 30 are crucial for key64-like transcription activity. We detected full-length recombinant CitPITP1 from E. coli, and visualized a CitPITP1-GFP fusion protein in plant cells, supporting the idea that CitPITP1 encodes a protein. However, although exon 4 of CitPITP1 contained key64, it did not demonstrate promoter activity in plants. Our study describes a new basal promoter, provides evidence for neofunction of gene elements across different kingdoms, and provides new knowledge for the modular design of promoters.

The evolution of insect metallothioneins

Metallothioneins (MTs) are a family of cysteine-rich metal-binding proteins that are important in the chelating and detoxification of toxic heavy metals. Until now, the short length and the low sequence complexity of MTs have hindered the inference of robust phylogenies, hampering the study of their evolution. To address this longstanding question, we applied an iterative BLAST search pipeline that allowed us to build a unique dataset of more than 300 MT sequences in insects. By combining phylogenetics and synteny analysis, we reconstructed the evolutionary history of MTs in insects. We show that the MT content in insects has been shaped by lineage-specific tandem duplications from a single ancestral MT. Strikingly, we also uncovered a sixth MT, MtnF, in the model organism Drosophila melanogaster. MtnF evolves faster than other MTs and is characterized by a non-canonical length and higher cysteine content. Our methodological framework not only paves the way for future studies on heavy metal detoxification but can also allow us to identify other previously unidentified genes and other low complexity genomic features.

Evolution of Phototransduction Genes in Lepidoptera

Vision is underpinned by phototransduction, a signaling cascade that converts light energy into an electrical signal. Among insects, phototransduction is best understood in Drosophila melanogaster. Comparison of D. melanogaster against three insect species found several phototransduction gene gains and losses, however, lepidopterans were not examined. Diurnal butterflies and nocturnal moths occupy different light environments and have distinct eye morphologies, which might impact the expression of their phototransduction genes. Here we investigated: 1) how phototransduction genes vary in gene gain or loss between D. melanogaster and Lepidoptera, and 2) variations in phototransduction genes between moths and butterflies. To test our prediction of phototransduction differences due to distinct visual ecologies, we used insect reference genomes, phylogenetics, and moth and butterfly head RNA-Seq and transcriptome data. As expected, most phototransduction genes were conserved between D. melanogaster and Lepidoptera, with some exceptions. Notably, we found two lepidopteran opsins lacking a D. melanogaster ortholog. Using antibodies we found that one of these opsins, a candidate retinochrome, which we refer to as unclassified opsin (UnRh), is expressed in the crystalline cone cells and the pigment cells of the butterfly, Heliconius melpomene. Our results also show that butterflies express similar amounts of trp and trpl channel mRNAs, whereas moths express ∼50× less trp, a potential adaptation to darkness. Our findings suggest that while many single-copy D. melanogaster phototransduction genes are conserved in lepidopterans, phototransduction gene expression differences exist between moths and butterflies that may be linked to their visual light environment.

Immediate Transcriptional Response to a Temperature Pulse under a Fluctuating Thermal Regime

The response of ectotherms to temperature stress is complex, non-linear, and is influenced by life stage and previous thermal exposure. Mortality is higher under constant low temperatures than under a fluctuating thermal regime (FTR) that maintains the same low temperature but adds a brief, daily pulse of increased temperature. Long term exposure to FTR has been shown to increase transcription of genes involved in oxidative stress, immune function, and metabolic pathways, which may aid in recovery from chill injury and oxidative damage. Previous research suggests the transcriptional response that protects against sub-lethal damage occurs rapidly under exposure to fluctuating temperatures. However, existing studies have only examined gene expression after a week or over many months. Here we characterize gene expression during a single temperature cycle under FTR. Development of pupating alfalfa leafcutting bees (Megachile rotundata) was interrupted at the red-eye stage and were transferred to 6°C with a 1-h pulse to 20°C and returned to 6°C. RNA was collected before, during, and after the temperature pulse and compared to pupae maintained at a static 6°C. The warm pulse is sufficient to cause expression of transcripts that repair cell membrane damage, modify membrane composition, produce antifreeze proteins, restore ion homeostasis, and respond to oxidative stress. This pattern of expression indicates that even brief exposure to warm temperatures has significant protective effects on insects exposed to stressful cold temperatures that persist beyond the warm pulse. Megachile rotundata's sensitivity to temperature fluctuations indicates that short exposures to temperature changes affect development and physiology. Genes associated with developmental patterning are expressed after the warm pulse, suggesting that 1 h at 20°C was enough to resume development in the pupae. The greatest difference in gene expression occurred between pupae collected after the warm pulse and at constant low temperatures. Although both were collected at the same time and temperature, the transcriptional response to one FTR cycle included multiple transcripts previously identified under long-term FTR exposure associated with recovery from chill injury, indicating that the effects of FTR occur rapidly and are persistent.

Copy Number Variation and Expression Analysis Reveals a Nonorthologous Pinta Gene Family Member Involved in Butterfly Vision

Vertebrate (cellular retinaldehyde-binding protein) and Drosophila (prolonged depolarization afterpotential is not apparent [PINTA]) proteins with a CRAL-TRIO domain transport retinal-based chromophores that bind to opsin proteins and are necessary for phototransduction. The CRAL-TRIO domain gene family is composed of genes that encode proteins with a common N-terminal structural domain. Although there is an expansion of this gene family in Lepidoptera, there is no lepidopteran ortholog of pinta. Further, the function of these genes in lepidopterans has not yet been established. Here, we explored the molecular evolution and expression of CRAL-TRIO domain genes in the butterfly Heliconius melpomene in order to identify a member of this gene family as a candidate chromophore transporter. We generated and searched a four tissue transcriptome and searched a reference genome for CRAL-TRIO domain genes. We expanded an insect CRAL-TRIO domain gene phylogeny to include H. melpomene and used 18 genomes from 4 subspecies to assess copy number variation. A transcriptome-wide differential expression analysis comparing four tissue types identified a CRAL-TRIO domain gene, Hme CTD31, upregulated in heads suggesting a potential role in vision for this CRAL-TRIO domain gene. RT-PCR and immunohistochemistry confirmed that Hme CTD31 and its protein product are expressed in the retina, specifically in primary and secondary pigment cells and in tracheal cells. Sequencing of eye protein extracts that fluoresce in the ultraviolet identified Hme CTD31 as a possible chromophore binding protein. Although we found several recent duplications and numerous copy number variants in CRAL-TRIO domain genes, we identified a single copy pinta paralog that likely binds the chromophore in butterflies.

Coping with living in the soil: the genome of the parthenogenetic springtail Folsomia candida

Background

Folsomia candida is a model in soil biology, belonging to the family of Isotomidae, subclass Collembola. It reproduces parthenogenetically in the presence of Wolbachia, and exhibits remarkable physiological adaptations to stress. To better understand these features and adaptations to life in the soil, we studied its genome in the context of its parthenogenetic lifestyle.

Results

We applied Pacific Bioscience sequencing and assembly to generate a reference genome for F. candida of 221.7 Mbp, comprising only 162 scaffolds. The complete genome of its endosymbiont Wolbachia, was also assembled and turned out to be the largest strain identified so far. Substantial gene family expansions and lineage-specific gene clusters were linked to stress response. A large number of genes (809) were acquired by horizontal gene transfer. A substantial fraction of these genes are involved in lignocellulose degradation. Also, the presence of genes involved in antibiotic biosynthesis was confirmed. Intra-genomic rearrangements of collinear gene clusters were observed, of which 11 were organized as palindromes. The Hox gene cluster of F. candida showed major rearrangements compared to arthropod consensus cluster, resulting in a disorganized cluster.

Conclusions

The expansion of stress response gene families suggests that stress defense was important to facilitate colonization of soils. The large number of HGT genes related to lignocellulose degradation could be beneficial to unlock carbohydrate sources in soil, especially those contained in decaying plant and fungal organic matter. Intra- as well as inter-scaffold duplications of gene clusters may be a consequence of its parthenogenetic lifestyle. This high quality genome will be instrumental for evolutionary biologists investigating deep phylogenetic lineages among arthropods and will provide the basis for a more mechanistic understanding in soil ecology and ecotoxicology.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3852-x) contains supplementary material, which is available to authorized users.

Gene Duplication and Gene Expression Changes Play a Role in the Evolution of Candidate Pollen Feeding Genes in Heliconius Butterflies

Heliconius possess a unique ability among butterflies to feed on pollen. Pollen feeding significantly extends their lifespan, and is thought to have been important to the diversification of the genus. We used RNA sequencing to examine feeding-related gene expression in the mouthparts of four species of Heliconius and one nonpollen feeding species, Eueides isabella. We hypothesized that genes involved in morphology and protein metabolism might be upregulated in Heliconius because they have longer proboscides than Eueides, and because pollen contains more protein than nectar. Using de novo transcriptome assemblies, we tested these hypotheses by comparing gene expression in mouthparts against antennae and legs. We first looked for genes upregulated in mouthparts across all five species and discovered several hundred genes, many of which had functional annotations involving metabolism of proteins (cocoonase), lipids, and carbohydrates. We then looked specifically within Heliconius where we found eleven common upregulated genes with roles in morphology (CPR cuticle proteins), behavior (takeout-like), and metabolism (luciferase-like). Closer examination of these candidates revealed that cocoonase underwent several duplications along the lineage leading to heliconiine butterflies, including two Heliconius-specific duplications. Luciferase-like genes also underwent duplication within lepidopterans, and upregulation in Heliconius mouthparts. Reverse-transcription PCR confirmed that three cocoonases, a peptidase, and one luciferase-like gene are expressed in the proboscis with little to no expression in labial palps and salivary glands. Our results suggest pollen feeding, like other dietary specializations, was likely facilitated by adaptive expansions of preexisting genes—and that the butterfly proboscis is involved in digestive enzyme production.

Multifaceted biological insights from a draft genome sequence of the tobacco hornworm moth, Manduca sexta

Manduca sexta, known as the tobacco hornworm or Carolina sphinx moth, is a lepidopteran insect that is used extensively as a model system for research in insect biochemistry, physiology, neurobiology, development, and immunity. One important benefit of this species as an experimental model is its extremely large size, reaching more than 10 g in the larval stage. M. sexta larvae feed on solanaceous plants and thus must tolerate a substantial challenge from plant allelochemicals, including nicotine. We report the sequence and annotation of the M. sexta genome, and a survey of gene expression in various tissues and developmental stages. The Msex_1.0 genome assembly resulted in a total genome size of 419.4 Mbp. Repetitive sequences accounted for 25.8% of the assembled genome. The official gene set is comprised of 15,451 protein-coding genes, of which 2498 were manually curated. Extensive RNA-seq data from many tissues and developmental stages were used to improve gene models and for insights into gene expression patterns. Genome wide synteny analysis indicated a high level of macrosynteny in the Lepidoptera. Annotation and analyses were carried out for gene families involved in a wide spectrum of biological processes, including apoptosis, vacuole sorting, growth and development, structures of exoskeleton, egg shells, and muscle, vision, chemosensation, ion channels, signal transduction, neuropeptide signaling, neurotransmitter synthesis and transport, nicotine tolerance, lipid metabolism, and immunity. This genome sequence, annotation, and analysis provide an important new resource from a well-studied model insect species and will facilitate further biochemical and mechanistic experimental studies of many biological systems in insects.

Transcriptome-Wide Differential Gene Expression in Bicyclus anynana Butterflies: Female Vision-Related Genes Are More Plastic

Vision is energetically costly to maintain. Consequently, over time many cave-adapted species downregulate the expression of vision genes or even lose their eyes and associated eye genes entirely. Alternatively, organisms that live in fluctuating environments, with different requirements for vision at different times, may evolve phenotypic plasticity for expression of vision genes. Here, we use a global transcriptomic and candidate gene approach to compare gene expression in the heads of a polyphenic butterfly. Bicyclus anynana have two seasonal forms that display sexual dimorphism and plasticity in eye morphology, and female-specific plasticity in opsin gene expression. Nonchoosy dry season females downregulate opsin expression, consistent with the high physiological cost of vision. To identify other genes associated with sexually dimorphic and seasonally plastic differences in vision, we analyzed RNA-sequencing data from whole head tissues. We identified two eye development genes (klarsicht and warts homologs) and an eye pigment biosynthesis gene (henna) differentially expressed between seasonal forms. By comparing sex-specific expression across seasonal forms, we found that klarsicht, warts, henna, and another eye development gene (domeless) were plastic in a female-specific manner. In a male-only analysis, white (w) was differentially expressed between seasonal forms. Reverse transcription polymerase chain reaction confirmed that warts and white are expressed in eyes only, whereas klarsicht, henna and domeless are expressed in both eyes and brain. We find that differential expression of eye development and eye pigment genes is associated with divergent eye phenotypes in B. anynana seasonal forms, and that there is a larger effect of season on female vision-related genes.