When does the female bias arise? Insights from the sex determination cascade of a flea beetle with a strongly skewed sex ratio.

Reproduction-manipulating bacteria like Wolbachia can shift sex ratios in insects towards females, but skewed sex ratios may also arise from genetic conflicts. The flea beetle Altica lythri harbors three main mtDNA strains that are coupled to three different Wolbachia infections. Depending on the mtDNA types, the females produce either offspring with a balanced sex ratio or exclusively daughters. To obtain markers that can monitor when sex bias arises in the beetle's ontogeny, we elucidated the sex determination cascade of A. lythri. We established a RT-PCR method based on length variants of dsx (doublesex) transcripts to determine the sex of morphologically indistinguishable eggs and larvae. In females of one mtDNA type (HT1/HT1*) known to produce only daughters, male offspring were already missing at the egg stage while for females of another type (HT2), the dsx splice variants revealed a balanced sex ratio among eggs and larvae. Our data suggest that the sex determination cascade in A. lythri is initiated by maternally transmitted female-specific tra (transformer) mRNA as primary signal. This tra mRNA seems to be involved in a positive feedback loop that maintains the production of the female splice variant, as known for female offspring in Tribolium castaneum. The translation of the maternally transmitted female tra mRNA must be inhibited in male offspring, but the underlying primary genetic signal remains to be identified. We discuss which differences between the mtDNA types can influence sex determination and lead to the skewed sex ratio of HT1.

Parallel metatranscriptome analysis reveals degradation of plant secondary metabolites by beetles and their gut symbionts.

Switching to a new host plant is a driving force for divergence and speciation in herbivorous insects. This process of incorporating a novel host plant into the diet may require a number of adaptations in the insect herbivores that allow them to consume host plant tissue that may contain toxic secondary chemicals. As a result, herbivorous insects are predicted to have evolved efficient ways to detoxify major plant defences and increase fitness by either relying on their own genomes or by recruiting other organisms such as microbial gut symbionts. In the present study we used parallel metatranscriptomic analyses of Altica flea beetles and their gut symbionts to explore the contributions of beetle detoxification mechanisms versus detoxification by their gut consortium. We compared the gut meta-transcriptomes of two sympatric Altica species that feed exclusively on different host plant species as well as their F1 hybrids that were fed one of the two host plant species. These comparisons revealed that gene expression patterns of Altica are dependent on both beetle species identity and diet. The community structure of gut symbionts was also dependent on the identity of the beetle species, and the gene expression patterns of the gut symbionts were significantly correlated with beetle species and plant diet. Some of the enriched genes identified in the beetles and gut symbionts are involved in the degradation of secondary metabolites produced by plants, suggesting that Altica flea beetles may use their gut microbiota to help them feed on and adapt to their host plants.

The draft genome of the specialist flea beetle Altica viridicyanea (Coleoptera: Chrysomelidae).

BACKGROUND: Altica (Coleoptera: Chrysomelidae) is a highly diverse and taxonomically challenging flea beetle genus that has been used to address questions related to host plant specialization, reproductive isolation, and ecological speciation. To further evolutionary studies in this interesting group, here we present a draft genome of a representative specialist, Altica viridicyanea, the first Alticinae genome reported thus far. RESULTS: The genome is 864.8 Mb and consists of 4490 scaffolds with a N50 size of 557 kb, which covered 98.6% complete and 0.4% partial insect Benchmarking Universal Single-Copy Orthologs. Repetitive sequences accounted for 62.9% of the assembly, and a total of 17,730 protein-coding gene models and 2462 non-coding RNA models were predicted. To provide insight into host plant specialization of this monophagous species, we examined the key gene families involved in chemosensation, detoxification of plant secondary chemistry, and plant cell wall-degradation. CONCLUSIONS: The genome assembled in this work provides an important resource for further studies on host plant adaptation and functionally affiliated genes. Moreover, this work also opens the way for comparative genomics studies among closely related Altica species, which may provide insight into the molecular evolutionary processes that occur during ecological speciation.

Ultrastructural and molecular characterization of Nosema alticae sp. nov. (Microsporidia: Nosematidae), pathogen of the flea beetle, Altica hampei Allard, 1867 (Coleoptera: Chrysomelidae).

In this study, the first microsporidian pathogen from Altica hampei (Coleoptera: Chrysomelidae) is described based on light microscopy, ultrastructural characteristics and comparative 16S SSU rDNA analysis. All developmental stages of the microsporidium are diplokaryotic and in direct contact with the host cell cytoplasm. Giemsa-stained mature spores are oval in shape and measured 3.82 ± 0.35 µm in length and 2.54 ± 0.27 µm in width. The polar filament of the binucleate spores is isofilar with 12-14 coils. Coils are 140.28 ± 4.88 nm (135.59-147.06; n = 36) in diameter and consist of six concentric layers of different electron density and thickness. The spores have a relatively thick (161.72 ± 29.19 nm) trilaminar spore wall. Morphological, ultrastructural and molecular features indicate that the described microsporidium belongs to the genus Nosema and is named Nosema alticae sp. nov.

Long-chain alkanes and fatty acids from Ludwigia octovalvis weed leaf surface waxes as short-range attractant and ovipositional stimulant to Altica cyanea (Weber) (Coleoptera: Chrysomelidae).

The importance of leaf surface wax compounds from the rice-field weed Ludwigia octovalvis (Jacq.) Raven (Onagraceae) was determined in the flea beetle Altica cyanea (Weber) (Coleoptera: Chrysomelidae). Extraction, thin layer chromatography and GC-MS and GC-FID analyses of surface waxes of young, mature and senescent leaves revealed 20, 19 and 19 n-alkanes between n-C15 and n-C35, respectively; whereas 14, 14 and 12 free fatty acids between C12:0 and C22:0 fatty acids were identified in young, mature and senescent leaves, respectively. Tricosane was predominant n-alkane in young and mature leaves, whilst eicosane predominated in senescent leaves. Heneicosanoic acid, palmitic acid and docosanoic acid were the most abundant free fatty acids in young, mature and senescent leaves, respectively. A. cyanea females showed attraction to 0.25 mature leaf equivalent surface waxes compared with young or senescent leaves in a short glass Y-tube olfactometer bioassay. The insects were attracted to a synthetic blend of 0.90, 1.86, 1.83, 1.95, 0.50 and 0.18 µg ml-1 petroleum ether of hexadecane, octadecane, eicosane, tricosane, palmitic acid and alpha-linolenic acid, respectively, comparable with the proportions as present in 0.25 mature leaf equivalent surface waxes. A. cyanea also laid eggs on a filter paper moistened with 0.25 mature leaf equivalent surface waxes or a synthetic blend of 0.90, 1.86, 1.83, 1.95, 0.50 and 0.18 µg ml-1 petroleum ether of hexadecane, octadecane, eicosane, tricosane, palmitic acid and alpha-linolenic acid, respectively. This finding could provide a basis for monitoring of the potential biocontrol agent in the field.

Evidence for feminized genetic males in a flea beetle using newly identified X-linked markers.

The equilibrium of sex ratios in sexually reproducing species is often disrupted by various environmental and genetic factors, including endosymbionts like Wolbachia. In this study, we explore the highly female-biased sex ratio observed in the flea beetle, Altica lythri, and its underlying mechanisms. Ancient hybridization events between Altica species have led to mitochondrial DNA introgression, resulting in distinct mitochondrial haplotypes that go along with different Wolbachia infections (HT1-wLytA1, HT1*- uninfected, HT2-wLytA2, and HT3-wLytB). Notably, beetles with some haplotypes exclusively produce female offspring, suggesting potential Wolbachia-induced phenomena such as feminization of genetic males. However, the observed female bias could also be a consequence of the ancient hybridization resulting in nuclear-cytoplasmic conflicts between introgressed mtDNA and nuclear genes. Through transcriptomic analysis and the program SEX-DETector, we established markers for genotypic sex differentiation for A. lythri, enabling genetic sexing via qPCR. Our findings suggest that feminization of genetic males is contributing to the skewed sex ratios, highlighting the intricate dynamics of sex determination and reproductive strategies in this flea beetle. This study provides valuable insights into the dynamics of genetic conflicts, endosymbionts, and sex ratios, revealing the novel phenomenon of genetic male feminization in the flea beetle A. lythri.

Integrated analysis of miRNA profiles and gut bacterial changes in Altica viridicyanea following antibiotic treatment.

The gut bacteria involves in insect homeostasis by playing essential roles in host physiology, metabolism, innate immunity, and so forth. microRNAs (miRNAs) are endogenous small noncoding RNAs that posttranscriptionally regulate gene expression to affect immune or metabolic processes in insects. For several non-model insects, the available knowledge on the relationship between changes in the gut bacteria and miRNA profiles is limited. In this study, we investigated the gut bacterial diversity, composition, and function from Altica viridicyanea feeding on normal- and antibiotic-treated host plants using 16S rRNA amplicon sequencing; antibiotics have been shown to affect the body weight and development time in A. viridicyanea, suggesting that the gut bacteria of the normal sample were more diverse and abundant than those of the antibiotic-fed group, and most of them were involved in various physical functions by enrichment analysis. Furthermore, we executed small RNA transcriptome sequencing using the two experimental groups to obtain numerous sRNAs, such as piRNAs, siRNAs, and known and novel miRNAs, by data mapping and quality control, and furthermore, a total of 224 miRNAs were identified as significantly differentially expressed miRNAs, of which some DEMs and their target genes participated in immune- and metabolism-related pathways based on GO and KEGG annotation. Besides, regarding the regulatory roles of miRNA and target genes, a interaction network of DEM-target gene pairs from eight immune- or metabolism-related signaling pathways were constructed. Finally, we discovered that DEMs from above pathways were significantly positively or negatively correlated with gut bacterial alterations following antibiotic treatment. Collectively, the observations of this study expand our understanding of how the disturbance of gut bacteria affects miRNA profiles in A. viridicyanea and provide new valuable resources from extreme ranges for future studies on the adaptive evolution in insects.

The lncRNA-mediated ceRNA network of Altica viridicyanea is involved in the regulation of the Toll/Imd signaling pathway under antibiotic treatment.

Long noncoding RNAs (lncRNAs) play significant roles in the regulation of mRNA expression or in shaping the competing endogenous RNA (ceRNA) network by targeting miRNA. The insect gut is one of the most important tissues due to direct contact with external pathogens and functions in the immune defense against pathogen infection through the innate immune system and symbionts, but there are limited observations on the role of the lncRNA-involved ceRNA network of the Toll/Imd pathway and correlation analysis between this network and bacterial microbiota in the Altica viridicyanea gut. In this research, we constructed and sequenced six RNA sequencing libraries using normal and antibiotic-reared samples, generating a total of 17,193 lncRNAs and 26,361 mRNAs from massive clean data by quality control and bioinformatic analysis. Furthermore, a set of 8,539 differentially expressed lncRNAs (DELs) and 13,263 differentially expressed mRNAs (DEMs), of which related to various immune signaling pathways, such as the Toll/Imd, JAK/STAT, NF-κB, and PI3K-Akt signaling pathways, were obtained between the two experimental groups in A. viridicyanea. In addition, numerous GO and KEGG enrichment analyses were used to annotate the DELs and their target genes. Moreover, six Toll family members and nineteen signal genes from the Toll/Imd signaling pathway were identified and characterized using online tools, and phylogenetic analyses of the above genes proved their classification. Next, a lncRNA-miRNA-mRNA network of the Toll/Imd pathway was built, and it contained different numbers of DEMs in this pathway and related DELs based on prediction and annotation. In addition, qRT-PCR validation and sequencing data were conducted to show the expression patterns of the above DELs and DEMs related to the Toll/Imd signaling pathway. Finally, the correlated investigations between DELs or DEMs of the Toll/Imd signaling pathway and most changes in the gut bacterial microbiota revealed significantly positive or negative relationships between them. The present findings provide essential evidence for innate immune ceRNAs in the beetle gut and uncover new potential relationships between innate immune pathways and the gut bacterial microbiota in insects.