
==== Front
PNAS Nexus
PNAS Nexus
pnasnexus
PNAS Nexus
2752-6542
Oxford University Press US

39228812
10.1093/pnasnexus/pgae348
pgae348
Biological, Health, and Medical Sciences
AcademicSubjects/MED00010
AcademicSubjects/SCI00010
AcademicSubjects/SOC00010
PNAS_Nexus/evolution
Cell-based assays and comparative genomics revealed the conserved and hidden effects of Wolbachia on insect sex determination
https://orcid.org/0000-0001-9912-3489
Arai Hiroshi Institute of Agrobiological Sciences, National Agriculture and Food Research Organization (NARO), Ibaraki 305-0851, Japan
United Graduate School of Agricultural Science, Tokyo University of Agriculture and Technology, Tokyo 183-8509, Japan

https://orcid.org/0000-0002-7433-5466
Herran Benjamin Institute of Agrobiological Sciences, National Agriculture and Food Research Organization (NARO), Ibaraki 305-0851, Japan

https://orcid.org/0000-0002-2929-0887
Sugimoto Takafumi N Institute of Agrobiological Sciences, National Agriculture and Food Research Organization (NARO), Ibaraki 305-0851, Japan

https://orcid.org/0000-0003-2245-3986
Miyata Mai Faculty of Engineering, University of Fukui, Fukui 910-8507, Japan

https://orcid.org/0000-0003-0929-0991
Sasaki Tetsuhiko Graduate School of Agriculture, Honeybee Science Research Center, Research Institute, Tamagawa University, Tokyo 194-8610, Japan

https://orcid.org/0000-0002-9026-9825
Kageyama Daisuke Institute of Agrobiological Sciences, National Agriculture and Food Research Organization (NARO), Ibaraki 305-0851, Japan

Wilson-Pham Marenda Editor
To whom correspondence should be addressed: Email: dazai39papilio@gmail.com (H.A.), Email: kagymad@affrc.go.jp (D.K.)
Present address: Institute of Infection, Veterinary and Ecological Sciences, University of Liverpool, Liverpool, United Kingdom.

Present address: Department of Biology, University of Nantes, Nantes, France.

Competing Interest: The authors declare no competing interests.

9 2024
22 8 2024
22 8 2024
3 9 pgae34801 3 2024
06 8 2024
03 9 2024
© The Author(s) 2024. Published by Oxford University Press on behalf of National Academy of Sciences.
2024
https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

It is advantageous for maternally transmitted endosymbionts to skew the sex ratio of their hosts toward females. Some endosymbiotic bacteria, such as Wolbachia, cause their insect hosts to exclusively produce female offspring through male killing (MK) or feminization. In some lepidopteran insects, MK is achieved by affecting the sex-determining process in males, and a unique mechanism of MK and its functional link with feminization have been implicated. However, comparative analysis of these phenotypes is often difficult because they have been analyzed in different host–symbiont systems, and transinfection of Wolbachia across different hosts is often challenging. In this study, we demonstrated the effects of nine Wolbachia strains on the splicing of sex-determining genes in Lepidoptera by fixing the host genetic background using a cell culture system. Cell transinfection assays confirmed that three MK-inducing Wolbachia strains and one feminization-inducing Wolbachia strain increased the female-type splicing products of the core sex-determining genes doublesex, masculinizer, and zinc finger protein 2. Regarding Wolbachia strains that do not induce MK/feminization, three had no effect on these sex-determining genes, whereas two strains induced female-type splicing of masculinizer and doublesex but not zinc finger protein 2. Comparative genomics confirmed that homologs of oscar, the Wolbachia gene responsible for MK in Ostrinia, were encoded by four MK/feminizing Wolbachia strains, but not by five non-MK/nonfeminizing strains. These results support the conserved effects underlying MK and feminization induced by oscar-bearing Wolbachia and suggested other potential mechanisms that Wolbachia might employ to manipulate host sex.

cell culture
endosymbiosis
male killing
Oscar
Wolbachia
Japan Society for the Promotion of Science 10.13039/501100001691 20F40719 and 21J00895 20K06084 22K14902 22K20663 23H02229 24H02293 Bio-oriented Technology Research Advancement Institution 10.13039/501100007173 JPJ009237
==== Body
pmcSignificance Statement

Arthropods often carry maternally transmitted symbionts, such as Wolbachia. The lack of paternal transmission frequently led them to evolve host manipulating behaviors in the favor of female hosts, substantiated by male killing and feminization. Although Wolbachia induces these phenotypes in a wide range of insects, the underlying mechanisms, diversity, and commonality remain largely unclear. In this study, a combination of transinfection assays and comparative genomics revealed the conserved effects of male killing and feminizing Wolbachia strains on lepidopteran sex determination. Furthermore, some nonmale killing/nonfeminizing Wolbachia strains were also shown to have an inherent ability to influence sex determination, albeit in a different manner, suggesting the potential for multiple mechanisms to manipulate host sex.

Introduction

Most animal species have distinct male and female forms that are determined by their genomic information during early development in a process called sex determination (1, 2). However, in many cases, sex-determining systems are manipulated or overridden by maternally transmitted cytoplasmic elements (3–6). In arthropods, several bacterial endosymbionts manipulate host reproduction to produce only female progenies (parthenogenesis) and by killing (i.e. male killing [MK]) or feminizing males (feminization [FM]) (7–11). The maternally inherited bacterium Wolbachia (Alphaproteobacteria) is a genus of typical microbes that manipulate host reproduction through various mechanisms (12). For example, Wolbachia species induce parthenogenesis, MK, and FM in diverse insects belonging to phylogenetically distinct taxa (e.g. Diptera, Lepidoptera, Hymenoptera, and Coleoptera) and other arthropods (8, 9, 13–16). Wolbachia species are estimated to inhabit at least 40% of all insect species, making them among the most widespread endosymbionts globally (12, 17).

It has been hypothesized that the mechanisms underlying MK and FM are associated with sex-determination cascades in insects (18). Indeed, some MK Wolbachia strains directly interact with sex-determination systems in lepidopteran insects (19–22). The lepidopteran sex-determination system consists of multiple transcriptional regulators, some of which exhibit sex-linked expression and/or splicing isoforms. Lepidopterans generally have a female heterogametic sex chromosome system (e.g. WZ/ZZ) and employ the dosage compensation system, which equalizes the sex-linked (Z-chromosome-linked) gene dose between males and females (20, 21). The dosage compensation system is regulated by a Z-linked CCCH-type zinc finger motif encoding gene, masculinizer (masc), which is critical for development as a male. The masc gene also regulates the downstream master regulator of sex determination and differentiation, doublesex (dsx), which exhibits sex-dependent splicing isoforms (dsxF in females and dsxM in males) (23, 24). Ostrinia and Homona male moths use these factors in their sex-determination cascades, and MK Wolbachia strains disrupt the expression of the masc gene and induce the “female” isoform of dsx (dsxF), leading to a “mismatch” between genetic sex (male: ZZ sex chromosome constitution) and phenotypic sex (female: based on dsxF) and male death (19–22). In addition to dsx, masc also exhibits splicing variation in the untranslated region between males and females of Ostrinia and Homona moths, and MK Wolbachia suppresses the expression levels of the male-type masc isoforms in male embryos (22, 25). Although the extent to which it is involved in sex-determination remains unclear, Ostrinia moths exhibit male- and female-type splicing variants of the autosomal CCCH-type zinc finger motif encoding gene zinc finger protein 2 (znf2) (25), which is an upstream regulator of dsx in the silkworm Bombyx mori (26). In male embryos, MK Wolbachia (wSca) induces a female-specific transcript variant and reduces the male-type isoforms of znf2. Similarly, the FM Wolbachia strain wFem induces dsxF production in Eurema butterflies (27).

More recently, a Wolbachia protein Oscar (abbreviation of “Osugoroshi [i.e. MK, in Japanese] protein containing CifB C-terminus-like domain and many Ankyrin Repeats”) was shown to recapitulate the Wolbachia-induced MK phenotype in Ostrinia (23). The Oscar protein degrades and interacts Masc, which leads to the failure of dosage compensation and the production of female-type dsx isoforms in Ostrinia male moths (23, 24). Furthermore, an Oscar homolog was found in the prophage region of an MK Wolbachia strain in Homona moths (28). These findings imply a general role for Oscar in MK in lepidopteran insects, which have not been demonstrated. Based on the Wolbachia strains sequenced to date, Oscar is not widely conserved among MK Wolbachia strains (28). In addition, transgenically overexpressed Oscar does not function in Drosophila flies, which do not employ masc in sex determination. Instead, another Wolbachia factor, Wmk, induces male lethality (28, 29). In contrast to Oscar, Wmk homologs are widely conserved among Wolbachia strains regardless of their phenotypes (28, 29). These findings imply that Wolbachia uses multiple mechanisms to induce MK in insects, but their similarity/diversity and functional link with FM remain elusive.

The induction of Wolbachia-induced reproductive manipulation is deeply associated with the host genetic background (30–33). MK and FM can significantly distort the sex ratio of the host population toward females (13, 31, 34–36). Under extremely female-biased sex ratios, male production is strongly selected, and host suppressor mutations have evolved to mask the endosymbiont-induced MK in some insects, such as butterflies, ladybugs, lacewings, and planthoppers (36–40). Therefore, the absence of obvious Wolbachia-induced phenotypes in the natural host does not necessarily indicate the absence of Wolbachia MK capabilities. Thus, investigations of the inherent MK/FM abilities of Wolbachia strains are warranted to elucidate the genetic basis of these phenotypes.

In this study, we reveal a conserved and hidden feminizing effect induced by Wolbachia strains in Lepidoptera. Transinfection assays using a recently established cell culture system (25) revealed that MK/FM Wolbachia strains altered the splicing patterns of genes involved in the moth sex-determination cascade. Genomic analyses confirmed the presence of Oscar in MK/FM Wolbachia strains. Furthermore, an inherent feminizing ability was observed in Oscar-deficient Wolbachia strains that do not induce MK or FM in their native hosts. Based on these findings, we discuss the underlying mechanisms of MK and FM in Lepidoptera and highlight the diversity of causative factors.

Results

Transinfection assays revealed the feminizing abilities of Wolbachia strains in Ostrinia male cells

The NARO-Ossc-M1 (OsM1) cell line derived from a male moth (Ostrinia scapulalis) (25) was transinfected with the native Wolbachia strain wSca and eight nonnative Wolbachia strains from lepidopteran insects (Table 1). Consistent with previous findings (25), transinfection with wSca (native male killer) induced female-type splicing patterns in dsx (OsdsxF), masc (OsmascF), and znf2 (Osznf2F) at 6 weeks posttransinfection (6 wpt) compared with the findings in nontransinfected controls (Steel–Dwass test; Osdsx: P = 0.02; Osmasc: P = 0.02; Osznf2: P = 0.02; Fig. 1). Albeit less or more conspicuously, wCauA, wHm-t, and wFem, having MK/FM abilities in native or nonnative hosts, also increased the expression of OsdsxF (wCauA: P = 0.01, wHm-t: P = 0.04, wFem: P = 0.02), OsmascF (wCauA: P = 0.01, wHm-t: P = 0.04, wFem: P = 0.02), and Osznf2F (wCauA: P = 0.01, wHm-t: P = 0.04, wFem: P = 0.02) at 6 wpt compared with the control levels. In contrast, cells transinfected with either wKue, wCauB, or wCI (inducing cytoplasmic incompatibility but not MK/FM in native hosts) continued to exhibit male-type splicing patterns for dsx (OsdsxM), (wKue: P = 0.48, wCauB: P = 0.99, wCI: P = 0.99), masc (OsmascM) (wKue: P = 0.99, wCauB: P = 0.99, wCI: P = 0.44), and znf2 (Osznf2M) (wKue: P = 0.07, wCauB: P = 0.13, wCI: P = 0.80) similarly as observed in Wolbachia-free control cells. Unexpectedly, wHm-c, which does not induce MK/FM in its host (43, 44), up-regulated OsdsxF (P = 0.02) and OsmascF (P = 0.02), but not Osznf2F (P = 0.07), at 6 wpt compared with the control levels. Similar patterns were also observed for the non-MK/FM wNi1 strain in Trichoplusia ni, including the upregulation of OsdsxF (P = 0.13) and OsmascF (P = 0.13), albeit without significance, whereas no effect was observed on Osznf2F (P = 0.31; Fig. 1).

Fig. 1. Feminizing effects of Wolbachia strains on Ostrinia male cells. Relative expression of OsdsxF/OsdsxM (A), OsmascF/OsmascM (B), and Osznf2F/Osznf2M (C) was quantified in Wolbachia-transinfected and control cells. Different letters above the whisker plots indicate significant differences based on the Steel–Dwass test (P < 0.05). Dot plots indicate replicates.

Table 1. Characteristics of the Wolbachia strains used in this study.

Strain	wKue	wCauA	wCauB	wSca	wCI	wCI + wFem	wNi1	wHm-t	wHm-c	
Phenotype	CI	CI, MKa	CI	MK	CI	FMb	Non-MK,
Non-FM	MK	Non-MK,
Non-CI	
Host	E. kuehniella	C. cautella	C. cautella	O. scapulalis	Eu. mandarina	Eu. mandarina	T. ni	H. magnanima	H. magnanima	
Reference	Sasaki et al. (30)	Sasaki et al. (30)	Sasaki et al. (30)	Sugimoto and Ishikawa (19)	Hiroki et al. (41)	Hiroki et al. (41)	Miyata (pers. obs.)	Arai et al. (42)	Arai et al. (43)	
MK, male killing; NMK, nonmale killing; Fem, feminization; CI, cytoplasmic incompatibility.

a wCauA induces CI in its native host C. cautella, but MK in E. kuehniella (30).

b wFem always coexists with wCI.

Oscar is conserved among MK and FM Wolbachia strains in Lepidoptera

A series of genome assembly and polishing processes were used to construct closed genomes of the Wolbachia strains wSca, wCauA, wCauB, wKue, wCI, wNi1, and wHm-c, which were used for the transactivation assays (Table 2). Comparative genomics using seven newly sequenced strains and the wHm-t strain (28) identified 788 protein clusters shared by the Wolbachia strains, including Cifs (CifB and CifA) and Wmk homologs (Fig. 2A and Table S1). Oscar homologs were the only protein clusters shared by MK Wolbachia (wSca, wCauA, and wHm-t) but are absent in non-MK Wolbachia strains (wCauB, wKue, wCI, wNi1, and wHm-c). As observed in the Oscar proteins of wFur and wSca (23), the Oscar homologs identified in this study encoded ankyrin repeats and the CifB C-terminus-like domains (papain-like cysteine proteinase domains; Fig. 2A and B). An Oscar homolog was also identified in the sequencing data derived from cells coinfected with wFem and wCI. Although we could not reconstruct the genome of the wFem strain due to technical difficulties in isolating wFem from the coinfecting wCI, the Oscar homolog is presumed to be encoded by wFem because the wCI strains used in this study did not carry Oscar homologs in their genomes (Fig. 2C and Table S1). Further diagnostic PCR confirmed that the Oscar homolog was consistently present in wild-caught Eurema mandarina butterflies infected with wFem and wCI (n = 37) but not in those infected with wCI alone (n = 77), suggesting that the Oscar homolog is encoded by wFem. Moreover, wFem-infected females (n = 3) of Eurema hecabe, a close relative of Eu. mandarina, were Oscar positive.

Fig. 2. Phylogenies and structures of the Oscar and Cif proteins. A) Protein clusters conserved among the Wolbachia strains used in this study. The number of clusters is shown on the left, and the number of proteins contained in the clusters is indicated by black ellipses. B) Structure of Oscar and its homologs. aa, amino acids. The structure of Oscar (1,830 aa), derived from the MK wFur strain, is based on Katsuma et al. (23). C) Phylogeny of Oscar and CifB homologs. Accession and gene numbers are given in parentheses.

Table 2. Genomic features of the Wolbachia strains used in this study.

Genome ID	wKue	wCauA	wCauB	wSca	wCI	wNi1	wHm-c	wHm-t	
Contigs	1	1	1	1	1	1	1	1	
Genome size (bp)	1,305,037	1,713,404	1,703,391	1,322,867	1,358,360	1,410,166	1,492,082	1,542,158	
G + C content (%)	35.3	35	34.2	33.7	34.1	34.1	34.0	34.1	
CDS genes	1,276	1,595	1,632	1,258	1,249	1,335	1,347	1,426	
Super group	A	A	B	B	B	B	B	B	
rRNA	3	3	3	3	3	3	3	3	
tRNA	34	34	34	34	33	34	34	34	
Native hosts	E. kuehniella	C. cautella	C. cautella	O. scapulalis	Eu. mandarina	T. ni	H. magnanima	H. magnanima	
Phenotype	CI	CI
MK	CI	MK	CI	Non-MK
Non-FM	Non-MK
Non-FM
Non-CI	MK	
Sources	This study	This study	This study	This study	This study	This study	This study	Arai et al. (28)	
MK, male killing; NMK, non-male killing; CI, cytoplasmic incompatibility.

To assess the phylogenetic associations of the Oscar homologs, we further searched Oscar homologs using the GenBank database, which identified Oscar-like proteins in four Wolbachia strains (wGro in Epagoge grotiana [Tortricidae, Lepidoptera], wSpl in Cydia splendana [Tortricidae, Lepidoptera], wFla in Endotricha flammealis [Pyralidae, Lepidoptera], and wObs in Cardiocondyla obscurior [Formicidae, Hymenoptera]). Since Oscar proteins generally encoded CifB C-terminus-like domains, we constructed a phylogenetic tree of Oscar alongside CifB proteins found in the Wolbachia genomes we sequenced. Phylogenetic analysis classified the Oscar proteins into three clades: type I Oscar (Oscar [wFur and wSca], Ec-Oscar [wCauA], and Oscar-like proteins in wGro [WP_265024289.1], wSpl [WP_264339188.1], and wFla [WP_264954178.1]), type II Oscar (Hm-Oscar [wHm-t] and Em/Eh-Oscar [wFem]), and type III Oscar (wObs [WP_174516118.1]) (Figs. 2C and S1). These Oscar homologs were highly diverse in sequences, and notably, Oscar in the newly sequenced wSca strain was not identical to that in the previously sequenced wSca strain (similarity: 90.4% in 1,797 amino acids; bit score: 3,241, BLASTp) and wFur (similarity: 94.2% in 1,830 amino acids; bit score: 3,432) (45) because of deletions in ankyrin repeats (Table S1). Intriguingly, the type II Oscar homologs Hm-Oscar and Em-Oscar exhibited higher homology than the wSca strains (similarity: 98.5% in 1,181 amino acids; bit score: 2,362). The MLST and genome-based phylogeny of the Wolbachia strains used in our analyses did not match the phylogeny of Oscar, suggesting horizontal transfer and/or rapid divergence of Oscar (Figs. 3 and S1).

Fig. 3. Phylogenies and phenotypes of Wolbachia strains and their effects on the sex-determination system in Ostrinia male cells. wsp and MLST genes were used to construct the phylogenetic tree. The classification of Oscar proteins was based on the phylogeny presented in Fig. 2C. The effects on splicing of each sex-determining gene analyzed in this study are presented as either Yes (feminized) or No (not affected). MK, male killing; NMK, nonmale killing; Fem, feminization; CI, cytoplasmic incompatibility; NA, not assessed.

Oscar is absent in wHm-c and wNi that interact with insect sex-determination systems

The wNi1 and wHm-c strains, which induced female-type splicing of masc (OsmascF) and dsx (OsdsxF) in Ostrinia male cells, did not carry Oscar homologs (Fig. 3 and Table S1). In addition, comparative genomics did not identify any proteins that were exclusively present in Wolbachia strains that affected the host sex-determination system (wHm-t, wSca, wCauA, wNi1, and wHm-c) but not in other strains (wKue, wCauB, and wCI; Fig. 2A). A protein cluster specific to wHm-c and wNi1 (i.e. absent in the other Wolbachia strains) consisted of three hypothetical proteins (63 amino acids: wHmcTK_11260, wNi_10160, and wNi_12630) that did not encode obvious domains and exhibited no homology to known MK genes (Tables S1 and S2). Genes carried by wHm-c or wNi1 that were present in some MK strains (wSca, wCauA, and wHm-t) but absent in non-MK strains (wCauB, wKue, and wCI; Table S2) included a wmk homolog (wNi_11830), whereas others displayed extremely low or no homology to known MK genes (i.e. Oscar (23), Spaid (46), and PVMKp1 (47)).

Discussion

Our cell transinfection assays proved effective in assessing the feminizing ability of Wolbachia strains. Comparative genomics further revealed the presence of Oscar homologs in these strains, which induced female-specific splicing of three sex-determining genes (dsx, masc, and znf2) in male Ostrinia cells. As previously argued (19, 27), the present study highlights a common mechanism underlying Wolbachia-induced MK and FM in lepidopteran insects: Oscar-induced suppression of Masc, leading to female-type sex determination and disruption of the dosage compensation system. Among the Wolbachia strains that feminized sex determination, wSca was the most efficient at inducing the female-type isoforms of the dsx and znf2 genes. This efficiency likely reflects its long-term adaptation to its native host, O. scapulalis, from which the OsM1 cell is derived. Notably, consistent with previous findings (25), the expression ratio of dsx (OsdsxF vs. OsdsxM) was lower than that of masc (OsmascM vs. OsmascF) and znf2 (Osznf2M vs. Osznf2F) in the MK/FM Wolbachia-transinfected cells. These patterns may indicate a rapid effect of Wolbachia on the upstream sex determinants, with a subsequent time lag before influencing the expression levels of downstream factors.

Unexpectedly, our study also revealed that two Oscar-deficient Wolbachia strains (wNi1 and wHm-c) affect the sex-determination system of Ostrinia, albeit likely through a different mechanism from Oscar-bearing Wolbachia strains: Oscar-bearing Wolbachia strains (i.e. wCauA, wSca, wFem, and wHm-t) affected Osdsx, Osmasc, and Osznf2, whereas Oscar-deficient wNi1 and wHm-c affected Osdsx and Osmasc but not Osznf2. These differences probably arose from the different machineries caused by different Wolbachia genes, i.e. oscar and other unknown genes. In B. mori, Znf2 is an upstream factor in the sex-determination cascade that regulates dsx sex-specific splicing (26). However, as the sex-determining gene cascade of Ostrinia (e.g. the hierarchical relationships among Osmasc, Osznf2, and Osdsx) is not fully understood, it remains unclear why only Osznf2 is not affected by wNi1 or wHm-c. One possibility is that Oscar can suppress the functions of Znf2 (or the cascade containing Znf2) as well as Masc, whereas the factors carried by wHm-c and wNi1 cannot act on Znf2 or its upstream components. Our study illustrated that closely related bacteria in the genus Wolbachia can activate different functions to manipulate sex determination in a single insect species. Further investigation of the causes is warranted.

This study provides compelling evidence that Wolbachia strains (wCauA, wNi1, and wHm-c) that do not induce MK or FM in their native hosts (43, 44, 48) retain the inherent ability to manipulate the sex-determination systems of Ostrinia. Despite its ability to alter the sex-determination cascade in Ostrinia cells, wHm-c does not affect the sex-determination cascade in the natural host Homona magnanima (22). These findings suggest that the abilities of Wolbachia to manipulate insect reproduction have been suppressed by the natural hosts but have been maintained through evolution. Wolbachia-induced reproductive manipulation can be influenced by the genetic background of the host (30–33). Under the condition of a female-biased sex ratio, mutations in the host genome that rescue males can be favored by selection. Indeed, the spread of suppressors of symbiont-induced MK has been observed in nature (18, 31, 37–40). For example, the spread of an MK suppressor in the butterfly Hypolimnas bolina was almost complete within 5 years (31, 36). In Cadra cautella, an MK phenotype, presumably induced by the Oscar-bearing wCauA, was recorded in the 1970s (49). In contrast wCauA did not induce MK in C. cautella collected around 2000 (48, 50), although it did induce MK when transferred to the closely related host Ephestia kuehniella (30). These findings suggest that C. cautella may have evolved a suppressor against wCauA-induced MK between the 1970s and 2000.

Since wHm-c relatives are found in many species of tortrix moths (43, 44), it is possible that a suppressor against wHm-c relatives evolved and became fixed in the ancestral species of tortrix moths. This suppressor cannot suppress the feminizing ability of Oscar-bearing wHm-t, and thus, the mechanistic basis of this suppressor would be different from the suppressor of C. cautella against Oscar-bearing wCauA. Consequently, a new suppressor against wHm-t may arise in H. magnanima in the future, leading to an evolutionary arms race between Wolbachia and the host. Recent studies have shown that MK microbes, including Wolbachia, have evolved diverse MK genes/mechanisms (22, 28, 46, 47). In response to these diverse MK mechanisms, host insects may have evolved various suppressors or suppression mechanisms, potentially masking the Wolbachia phenotypes.

In summary, our study demonstrates that the feminizing effect of Wolbachia strains on the sex-determining gene cascade is most likely the mechanistic basis of MK and FM in lepidopteran insects. This study highlights the effectiveness of combining cell culture systems and genomic analyses to uncover the inherent ability of Wolbachia to manipulate sex. These approaches could also prove effective in elucidating the mechanistic interplay between the host and other selfish reproductive manipulators. A molecular understanding of the commonality and diversity of microbial reproductive manipulations will contribute to a better understanding of the evolutionary interactions between selfish elements and their hosts.

Materials and methods

Wolbachia strains and insects

Ephestia kuehniella lines transinfected with wKue, wCauA, and wCauB in each were maintained in the laboratory, as described by Sasaki et al. (30). Females infected with wCauA were mated with uninfected males (50). The insects were reared on a diet consisting of wheat bran, dried yeast, and glycerol (20:1:2) at 25 °C under a 16-h/8-h light/dark photoperiod. Homona magnanima infected with wHm-t (Taoyuan, Taiwan (42)) or wHm-c (Takao, Japan (28)) was maintained using Silkmate 2S (Nosan Co., Yokohama, Japan). The MK host line WT12 was maintained by crossing it with a Wolbachia-free normal sex ratio line, as described by Arai et al. (42). For Eu. mandarina, all-female blood coinfected with wFem and wCI and a normal sex ratio line singly infected with wCI were collected on Tanegashima Island, Japan, and maintained on an artificial diet, as reported by Kageyama et al. (27). Trichoplusia ni harboring wNi1 was collected in Matsudo, Chiba, Japan. The MK wSca strain, maintained in BmM2 cells as described by Herran et al. (25), was used for the following assays.

Transinfection of Wolbachia into cultured cells

Fat bodies septically isolated from Wolbachia-infected insects were placed in the cell lines AeAl2 (wCauA, wCauB, wKue, wNi1, wHm-t, and wHm-c) and BmN4 (wCI and wFem). After confirming that Wolbachia was stably maintained by diagnostic PCR, Wolbachia-positive AeAl2/BmN4 cells were used as donors for Wolbachia transinfection into the O. scapulalis NARO-Ossc-M1 (OsM1) cell line. Wolbachia-infected AeAl2/BmN2 cells were passed through a 5.0-µm filter (cat. no. FJ25ASCCA050PL01, GVS, Via Roma, Italy), and six drops (∼300 µL) of the flow-through were added to the recipient OsM1 cell line.

Reverse-transcription polymerase chain reaction

Total RNA (200–500 ng) extracted from harvested cells using Isogen II (Nippon Gene, Tokyo, Japan) or TRIzol RNA Isolation Reagent (Thermo Fisher Scientific, Waltham, MA, USA) was reverse-transcribed using a PrimeScript II 1st strand cDNA Synthesis Kit (Takara Bio, Shiga, Japan), according to the manufacturer's protocol. qPCR was performed with 5 µL of KOD SYBR (Toyobo), 0.4 µL each of the forward and reverse primers (10 pmol/µL), 2.2 µL of water, and 2.0 µL of the cDNA template. qPCR was performed in a LightCycler 96 System (Roche, Basel, Switzerland) with a temperature profile of 180 s at 95 °C; 40 cycles of 8 s at 98 °C, 10 s at 60 °C, and 10 s at 68 °C; and heating to 90 °C for melting curve analysis. The primers used in this study are listed in Table S3. The relative expression (Osmasc, Osdsx, and Osznf2 vs. the control gene Osef1a) and the ratio of male-to-female splice variants were estimated.

Genomic analyses

wSca (in BmM2 cells, described by Herran et al. (25)), wHm-c, wNi1, wCauA, wKue, wCI, and wFem and wCI (i.e. coinfected into AeAl2 cells) were purified from cells harvested in 150 mL flasks, as described by Iturbe-Ormaetxe et al. (51) with several modifications. Cells were pelleted by centrifugation (900 × g for 5 min) and homogenized with glass beads using Multi-beads Shocker MB3000 at 2,000 rpm for 40 s (Yasui Kikai Co., Osaka, Japan). The lysates were passed through 5.0 and 1.2 µm filters and centrifuged at 15,000 × g for 30 min at 4 °C. Pellet-containing Wolbachia cells were subjected to DNA extraction using a NanoBind Big DNA Tissue kit (Circulomics, Baltimore, MD, USA) according to the manufacturer's protocol. The extracted DNA was sequenced using an Illumina short read (paired-end 150 bp) and MinION with a Rapid sequencing kit (Oxford Nanopore, Oxford, UK) with a Nanopore flange flow cell (R9.4) (Oxford Nanopore). For wCauB, DNA extracted from host insects harboring Wolbachia using a NanoBind Big DNA Tissue kit was sequenced using Illumina (paired-end 150 bp) and Nanopore platforms (ligation sequencing kit with Nanopore MinION flow cell [R10.4]; Oxford Nanopore). The Nanopore data were assembled using Flye v1.6 (52) and polished 3–5 times with Illumina data using Pilon (53) and minimap2 (54). The resulting circular Wolbachia genomes were annotated using DFAST (55). Protein homologies were analyzed using OrthoVenn 3 (https://orthovenn3.bioinfotoolkits.net). Wolbachia cifB, wmk, and oscar (accession numbers are presented in Fig. 2 and Table S1) (23, 29, 56) were used as queries to identify homologs from Wolbachia genomes using local BLASTn and BLASTp searches (default parameters). Motifs in oscar homologs were surveyed using InterPro (https://www.ebi.ac.uk/interpro/) and HHpred (https://toolkit.tuebingen.mpg.de/tools/hhpred). Oscar homologs were further analyzed using Blastn and Blastp searches against GenBank nr database.

Phylogenetic analysis

MLST (gatB, coxA, hcpA, ftsZ, and fbpA genes) and wsp fragments from the Wolbachia strains used in this study were aligned using ClustalW (57) with a Wolbachia endosymbiont in Brugia malayi (wBm) as an out-group. Single-copy genes (n = 638) conserved among Wolbachia strains (wFur, wSca, wHm-t, wHm-c, wNi1, wCauA, wCauB, wCI, wKue, wGro, wSpl, wFla, wObs, and wBm) were obtained using OrthoFinder (58) and concatenated, aligned, and trimmed using SeqKit (59), MAFFT (60), and trimAl (61), respectively. CifB and Oscar proteins were aligned using ClustalW (57), and trimmed using trimAl (61) in the default mode. These alignment files were used to construct a phylogenetic tree based on maximum likelihood with bootstrap resampling of 1,000 replicates using the IQTREE server (http://iqtree.cibiv.univie.ac.at/) and visualized using the iTOL web server (https://itol.embl.de/).

PCR detection of Em/Eh-oscar gene

Eurema mandarina collected from Tanegashima Island (Kagoshima, Japan) (35) and E. hecabe collected from Kohama Island (Okinawa, Japan) were subjected to DNA extraction using a DNeasy kit (Qiagen, Hilden, Germany), according to the manufacturer's protocol. The DNA concentration was adjusted to 10 ng/µL and subjected to PCR using primer sets amplifying the complete sequences of Em/Eh-oscar and Hm-oscar (Table S3). Em/Eh-oscar was amplified using Emerald Amp Max Master Mix (TaKaRa) at 94 °C for 3 min; 35 cycles of 94 °C for 30 s, 62 °C for 30 s, and 72 °C for 3 min; and a final extension at 72 °C for 7 min, as described by Arai et al. (28).

Statistical analysis

For RT-qPCR assays, we used the average cycle threshold value (Ct) for each sample and estimated the relative expression, as described by Herran et al. (25) and Sugimoto et al. (20). Elongation factor a (Ef1a) was used as the control gene. The 2−ΔCt (i.e. CtAve target gene −CtAve ef1a) and 2−ΔCtMal2−ΔCtFem (i.e. relative expression in males vs. females) values were calculated. The expression ratios of the splice variants (i.e. DsxM/DsxF, MascM/MascF, and Znf2M/Znf2F) were analyzed using the Steel–Dwass test in R software v4.0 (62).

Supplementary Material

pgae348_Supplementary_Data

Acknowledgments

The authors thank Kazuyo Watanabe for her helpful advice on cell culture and Yuuka Kimura for technical assistance.

Supplementary Material

Supplementary material including Supplementary Dataset 1 is available at PNAS Nexus online.

Funding

This work was supported by the Japan Society for the Promotion of Science (JSPS) Research Fellowships for Young Scientists (grant numbers 20F40719 and 21J00895), JSPS Grant-in-Aid for Scientific Research (grant numbers 20K06084, 22K14902, 22K20663, 23H02229, and 24H02293), Cabinet Office, Government of Japan, Moonshot Research and Development Program for Agriculture, Forestry, and Fisheries (funding agency: Bio-oriented Technology Research Advancement Institution; no. JPJ009237).

Author Contributions

H.A. collected and maintained the insects; transfected Wolbachia; conducted genome sequencing, qPCR, and data analysis; wrote the original manuscript; and revised the manuscript. B.H. transfected Wolbachia, conducted qPCR, and contributed to the discussion. T.N.S. collected the insects, transfected Wolbachia, conducted qPCR, and contributed to the discussion. M.M. collected the insects and contributed to the discussion. T.S. collected and maintained the insects and contributed to the discussion. D.K. transinfected Wolbachia; designed and conceptualized the experiments; and revised the original manuscript.

Preprints

This manuscript was posted on a preprint: https://doi.org/10.1101/2024.02.12.579973

Data Availability

The sequence read data are publicly available in DDBJ under the accession numbers PRJDB16497 (BioProject), DRA016973 (DDBJ Sequence Read Archive: DRA), and DRA017337 (DRA). Wolbachia genomes (contigs) are available in the DDBJ database under the following accession numbers: wCauA (SAMD00639214, AP028948), wCauB (SAMD00639215, AP028949), wKue (SAMD00639216, AP028950), wCI (SAMD00639217, AP028951), wCI-Fem (SAMD00639218, AP028952), wSca (SAMD00639219, AP028953), wNi1 (SAMD00639220, AP028954), and wHm-c (SAMD00654377, AP028994). Other data generated in this study are included in the manuscript and Supplementary material.
==== Refs
References

1 Werren  JH, Beukeboom  LW. 1998. Sex determination, sex ratios, and genetic conflict. Annu Rev Ecol Syst. 29 :233–261.
2 Mank  JE, Rideout  EJ. 2021. Developmental mechanisms of sex differences: from cells to organisms. Development. 148 :dev199750.34647574
3 Eberhard  WG . 1990. Evolution in bacterial plasmids and levels of selection. Q Rev Biol. 65 :3–22.2186429
4 Murlas Cosmides  L, Tooby  J. 1981. Cytoplasmic inheritance and intragenomic conflict. J Theor Biol. 89 :83–129.7278311
5 Hurst  GDD, Majerus  MEN. 1993. Why do maternally inherited microorganisms kill males?  Heredity. 71 :81–95.
6 Ma  W-J, Vavre  F, Beukeboom  LW. 2014. Manipulation of arthropod sex determination by endosymbionts: diversity and molecular mechanisms. Sex Dev. 8 :59–73.24355929
7 Hurst  GDD, et al  1997. The diversity of inherited parasites of insects: the male-killing agent of the ladybird beetle Coleomegilla maculata is a member of the Flavobacteria. Genet Res. 70 :1–6.
8 Hurst  GDD, et al  1999. Male-killing Wolbachia in two species of insect. Proc Biol Sci.  266 :735–740.
9 Narita  S, Kageyama  D, Nomura  M, Fukatsu  T. 2007. Unexpected mechanism of symbiont-induced reversal of insect sex: feminizing Wolbachia continuously acts on the butterfly Eurema hecabe during larval development. Appl Environ Microbiol. 73 :4332–4341.17496135
10 Ferree  PM, Avery  A, Azpurua  J, Wilkes  T, Werren  JH. 2008. A bacterium targets maternally inherited centrosomes to kill males in Nasonia. Curr Biol. 18 :1409–1414.18804376
11 Hurst  GDD, Frost  CL. 2015. Reproductive parasitism: maternally inherited symbionts in a biparental world. Cold Spring Harb Perspect Biol. 7 :a017699.25934011
12 Werren  JH, Baldo  L, Clark  ME. 2008. Wolbachia: master manipulators of invertebrate biology. Nat Rev Microbiol. 6 :741–751.18794912
13 Dyson  EA, Hurst  GDD. 2004. Persistence of an extreme sex-ratio bias in a natural population. Proc Natl Acad Sci U S A. 101 :6520–6523.15087494
14 Kageyama  D, Traut  W. 2004. Opposite sex–specific effects of Wolbachia and interference with the sex determination of its host Ostrinia scapulalis. Proc Biol Sci.  271 :251–258.15058435
15 Hurst  GD, Johnson  AP, Schulenburg  JH, Fuyama  Y. 2000. Male-killing Wolbachia in Drosophila: a temperature-sensitive trait with a threshold bacterial density. Genetics. 156 :699–709.11014817
16 Herran  B, et al  2021. Feminising Wolbachia disrupt Armadillidium vulgare insulin-like signalling pathway. Cell Microbiol. 23 :e13381.34312980
17 Zug  R, Hammerstein  P. 2012. Still a host of hosts for Wolbachia: analysis of recent data suggests that 40% of terrestrial arthropod species are infected. PLoS One. 7 :e38544.22685581
18 Hornett  EA, Kageyama  D, Hurst  GDD. 2022. Sex determination systems as the interface between male-killing bacteria and their hosts. Proc Biol Sci.  289 :20212781.35414231
19 Sugimoto  TN, Ishikawa  Y. 2012. A male-killing Wolbachia carries a feminizing factor and is associated with degradation of the sex-determining system of its host. Biol Lett. 8 :412–415.22219393
20 Sugimoto  TN, Kayukawa  T, Shinoda  T, Ishikawa  Y, Tsuchida  T. 2015. Misdirection of dosage compensation underlies bidirectional sex-specific death in Wolbachia-infected Ostrinia scapulalis. Insect Biochem Mol Biol. 66 :72–76.26453817
21 Fukui  T, et al  2015. The endosymbiotic bacterium Wolbachia selectively kills male hosts by targeting the masculinizing gene. PLoS Pathog. 11 :e1005048.26172536
22 Arai  H, et al  2023. Diverse molecular mechanisms underlying microbe-inducing male killing in the moth Homona magnanima. Appl Environ Microbiol. 89 :e0209522.37098937
23 Katsuma  S, et al  2022. A Wolbachia factor for male killing in lepidopteran insects. Nat Commun. 13 :6764.36376299
24 Fukui  T, et al  2024. Expression of the Wolbachia male-killing factor Oscar impairs dosage compensation in lepidopteran embryos. FEBS Lett. 598 :331–337.37985236
25 Herran  B, et al  2023. Cell-based analysis reveals that sex-determining gene signals in Ostrinia are pivotally changed by male-killing Wolbachia. PNAS Nexus. 2 :pgac293.36712932
26 Gopinath  G, Arunkumar  KP, Mita  K, Nagaraju  J. 2016. Role of Bmznf-2, a Bombyx mori CCCH zinc finger gene, in masculinisation and differential splicing of Bmtra-2. Insect Biochem Mol Biol. 75 :32–44.27260399
27 Kageyama  D, Yoshimura  K, Sugimoto  TN, Katoh  TK, Watada  M. 2017. Maternally transmitted non-bacterial male killer in Drosophila biauraria. Biol Lett. 13 :20170476.29021316
28 Arai  H, et al  2023. Combined actions of bacteriophage-encoded genes in Wolbachia-induced male lethality. iScience. 26 :106842.37250803
29 Perlmutter  JI, et al  2019. The phage gene wmk is a candidate for male killing by a bacterial endosymbiont. PLoS Pathog. 15 :e1007936.31504075
30 Sasaki  T, Kubo  T, Ishikawa  H. 2002. Interspecific transfer of Wolbachia between two lepidopteran insects expressing cytoplasmic incompatibility: a Wolbachia variant naturally infecting Cadra cautella causes male killing in Ephestia kuehniella. Genetics. 162 :1313–1319.12454075
31 Hornett  EA, et al  2006. Evolution of male-killer suppression in a natural population. PLoS Biol. 4 :e283.16933972
32 Jaenike  J . 2007. Spontaneous emergence of a new Wolbachia phenotype. Evolution. 61 :2244–2252.17767593
33 Veneti  Z, et al  2012. Loss of reproductive parasitism following transfer of male-killing Wolbachia to Drosophila melanogaster and Drosophila simulans. Heredity. 109 :306–312.22892635
34 Owen  DF . 1965. Change in sex ratio in an African butterfly. Nature. 206 :744–744.5833715
35 Kageyama  D, et al  2020. Persistence of a Wolbachia-driven sex ratio bias in an island population of Eurema butterflies. bioRxiv. 10.1101/2020.03.24.005017 , preprint: not peer reviewed.
36 Charlat  S, et al  2007. Extraordinary flux in sex ratio. Science. 317 :214–214.17626876
37 Majerus  TMO, Majerus  MEN. 2010. Intergenomic arms races: detection of a nuclear rescue gene of male-killing in a ladybird. PLoS Pathog. 6 :e1000987.20628578
38 Hayashi  M, Nomura  M, Kageyama  D. 2018. Rapid comeback of males: evolution of male-killer suppression in a green lacewing population. Proc Biol Sci.  285 :20180369.29669904
39 Yoshida  K, Sanada-Morimura  S, Huang  S-H, Tokuda  M. 2021. Silence of the killers: discovery of male-killing suppression in a rearing strain of the small brown planthopper, Laodelphax striatellus. Proc Biol Sci.  288 :20202125.33468006
40 Richardson  KM, et al  2023. A male-killing Wolbachia endosymbiont is concealed by another endosymbiont and a nuclear suppressor. PLoS Biol. 21 :e3001879.36947547
41 Hiroki  M, Kato  Y, Kamito  T, Miura  K. 2002. Feminization of genetic males by a symbiotic bacterium in a butterfly, Eurema hecabe (Lepidoptera: Pieridae). Naturwissenschaften. 89 :167–170.12061400
42 Arai  H, Lin  SR, Nakai  M, Kunimi  Y, Inoue  MN. 2020. Closely related male-killing and nonmale-killing Wolbachia strains in the oriental tea tortrix Homona magnanima. Microb Ecol. 79 :1011–1020.31820073
43 Arai  H, et al  2019. Multiple infection and reproductive manipulations of Wolbachia in Homona magnanima (Lepidoptera: Tortricidae). Microb Ecol. 77 :257–266.29931623
44 Arai  H, et al  2024. Conserved infections and reproductive phenotypes of Wolbachia symbionts in Asian tortrix moths. Environ Microbiol Rep. 16 :e13219.38070178
45 Muro  T, Hikida  H, Fujii  T, Kiuchi  T, Katsuma  S. 2023. Two complete genomes of male-killing Wolbachia infecting Ostrinia moth species illuminate their evolutionary dynamics and association with hosts. Microb Ecol. 86 :1740–1754.36810610
46 Harumoto  T, Lemaitre  B. 2018. Male-killing toxin in a bacterial symbiont of Drosophila. Nature. 557 :252–255.29720654
47 Kageyama  D, et al  2023. A male-killing gene encoded by a symbiotic virus of Drosophila. Nat Commun. 14 :1357.36914655
48 Sasaki  T, Massaki  N, Kubo  T. 2005. Wolbachia variant that induces two distinct reproductive phenotypes in different hosts. Heredity. 95 :389–393.16106260
49 Takahashi  F, Kuwahara  Y. 1970. Studies on sex pheromones of Pyralididae. III. The inheritance of the abnormal sex ratio condition in a strain of the almond moth, Cadra cautella Walker (Phycitinae). Botyu-Kagaku. 35 :11–21.
50 Sasaki  T, Ishikawa  H. 1999. Wolbachia infections and cytoplasmic incompatibility in the almond moth and the Mediterranean flour moth. Zool Sci. 16 :739–744.
51 Iturbe-Ormaetxe  I, Woolfit  M, Rancès  E, Duplouy  A, O’Neill  SL. 2011. A simple protocol to obtain highly pure Wolbachia endosymbiont DNA for genome sequencing. J Microbiol Methods. 84 :134–136.21047535
52 Kolmogorov  M, Yuan  J, Lin  Y, Pevzner  PA. 2019. Assembly of long, error-prone reads using repeat graphs. Nat Biotechnol. 37 :540–546.30936562
53 Walker  BJ, et al  2014. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 9 :e112963.25409509
54 Li  H . 2018. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 34 :3094–3100.29750242
55 Tanizawa  Y, Fujisawa  T, Nakamura  Y. 2018. DFAST: a flexible prokaryotic genome annotation pipeline for faster genome publication. Bioinformatics. 34 :1037–1039.29106469
56 Martinez  J, Klasson  L, Welch  JJ, Jiggins  FM. 2021. Life and death of selfish genes: comparative genomics reveals the dynamic evolution of cytoplasmic incompatibility. Mol Biol Evol. 38 :2–15.32797213
57 Larkin  MA, et al  2007. Clustal W and Clustal X version 2.0. Bioinformatics. 23 :2947–2948.17846036
58 Emms  DM, Kelly  S. 2019. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 20 :238.31727128
59 Shen  W, Le  S, Li  Y, Hu  F. 2016. SeqKit: a cross-platform and ultrafast toolkit for FASTA/Q file manipulation. PLoS One. 11 :e0163962.27706213
60 Katoh  K, Kuma  K, Toh  H, Miyata  T. 2005. MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 33 :511–518.15661851
61 Capella-Gutiérrez  S, Silla-Martínez  JM, Gabaldón  T. 2009. Trimal: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 25 :1972–1973.19505945
62 Ihaka  R, Gentleman  R. 1996. R: a language for data analysis and graphics. J Comp Graph Stat. 5 :299–314.
