Sex-biased gene content is associated with sex chromosome turnover in Danaini butterflies.

Sex chromosomes play an outsized role in adaptation and speciation, and thus deserve particular attention in evolutionary genomics. In particular, fusions between sex chromosomes and autosomes can produce neo-sex chromosomes, which offer important insights into the evolutionary dynamics of sex chromosomes. Here, we investigate the evolutionary origin of the previously reported Danaus neo-sex chromosome within the tribe Danaini. We assembled and annotated genomes of Tirumala septentrionis (subtribe Danaina), Ideopsis similis (Amaurina), Idea leuconoe (Euploeina) and Lycorea halia (Itunina) and identified their Z-linked scaffolds. We found that the Danaus neo-sex chromosome resulting from the fusion between a Z chromosome and an autosome corresponding to the Melitaea cinxia chromosome (McChr) 21 arose in a common ancestor of Danaina, Amaurina and Euploina. We also identified two additional fusions as the W chromosome further fused with the synteny block McChr31 in I. similis and independent fusion occurred between ancestral Z chromosome and McChr12 in L. halia. We further tested a possible role of sexually antagonistic selection in sex chromosome turnover by analysing the genomic distribution of sex-biased genes in I. leuconoe and L. halia. The autosomes corresponding to McChr21 and McChr31 involved in the fusions are significantly enriched in female- and male-biased genes, respectively, which could have hypothetically facilitated fixation of the neo-sex chromosomes. This suggests a role of sexual antagonism in sex chromosome turnover in Lepidoptera. The neo-Z chromosomes of both I. leuconoe and L. halia appear fully compensated in somatic tissues, but the extent of dosage compensation for the ancestral Z varies across tissues and species.

Comparative analysis of animal lifespan.

Comparative studies of aging are a promising approach to identifying general properties of and processes leading to aging. While to date, many comparative studies of aging in animals have focused on relatively narrow species groups, methodological innovations now allow for studies that include evolutionary distant species. However, comparative studies of aging across a wide range of species that have distinct life histories introduce additional challenges in experimental design. Here, we discuss these challenges, highlight the most pressing problems that need to be solved, and provide suggestions based on current approaches to successfully carry out comparative aging studies across the animal kingdom.

Genome sequence of Ophryocystis elektroscirrha, an apicomplexan parasite of monarch butterflies: cryptic diversity and response to host-sequestered plant chemicals.

Apicomplexa are ancient and diverse organisms which have been poorly characterized by modern genomics. To better understand the evolution and diversity of these single-celled eukaryotes, we sequenced the genome of Ophryocystis elektroscirrha, a parasite of monarch butterflies, Danaus plexippus. We contextualize our newly generated resources within apicomplexan genomics before answering longstanding questions specific to this host-parasite system. To start, the genome is miniscule, totaling only 9 million bases and containing fewer than 3,000 genes, half the gene content of two other sequenced invertebrate-infecting apicomplexans, Porospora gigantea and Gregarina niphandrodes. We found that O. elektroscirrha shares different orthologs with each sequenced relative, suggesting the true set of universally conserved apicomplexan genes is very small indeed. Next, we show that sequencing data from other potential host butterflies can be used to diagnose infection status as well as to study diversity of parasite sequences. We recovered a similarly sized parasite genome from another butterfly, Danaus chrysippus, that was highly diverged from the O. elektroscirrha reference, possibly representing a distinct species. Using these two new genomes, we investigated potential evolutionary response by parasites to toxic phytochemicals their hosts ingest and sequester. Monarch butterflies are well-known to tolerate toxic cardenolides thanks to changes in the sequence of their Type II ATPase sodium pumps. We show that Ophryocystis completely lacks Type II or Type 4 sodium pumps, and related proteins PMCA calcium pumps show extreme sequence divergence compared to other Apicomplexa, demonstrating new avenues of research opened by genome sequencing of non-model Apicomplexa.

A molecular cell biology toolkit for the study of meiosis in the silkworm Bombyx mori.

Meiosis is usually described as 4 essential and sequential processes: (1) homolog pairing; (2) synapsis, mediated by the synaptonemal complex; (3) crossing over; and (4) segregation. In this canonical model, the maturation of crossovers into chiasmata plays a vital role in holding homologs together and ensuring their segregation at the first meiotic division. However, Lepidoptera (moths and butterflies) undergo 3 distinct meiotic processes, only one of which is canonical. Lepidoptera males utilize 2 meiotic processes: canonical meiosis that produces nucleated fertile sperm, and a noncanonical meiosis that produces anucleated nonfertile sperm which are nonetheless essential for reproduction. Lepidoptera females, which carry heteromorphic sex chromosomes, undergo a completely achiasmate (lacking crossovers) meiosis, thereby requiring an alternative mechanism to ensure proper homolog segregation. Here, we report that the development of a molecular cell biology toolkit designed to properly analyze features of meiosis, including the synaptonemal complex structure and function, in the silkworm Bombyx mori. In addition to standard homology searches to identify Bombyx orthologs of known synaptonemal complex encoding genes, we developed an ortholog discovery app (Shinyapp) to identify Bombyx orthologs of proteins involved in several meiotic processes. We used this information to clone genes expressed in the testes and then created antibodies against their protein products. We used the antibodies to confirm the localization of these proteins in normal male spermatocytes, as well as using in vitro assays to confirm orthologous interactions. The development of this toolkit will facilitate further study of the unique meiotic processes that characterize meiosis in Lepidoptera.

The Prmt5-Vasa module is essential for spermatogenesis in Bombyx mori.

In lepidopteran insects, dichotomous spermatogenesis produces eupyrene spermatozoa, which are nucleated, and apyrene spermatozoa, which are anucleated. Both sperm morphs are essential for fertilization, as eupyrene sperm fertilize the egg, and apyrene sperm is necessary for the migration of eupyrene sperm. In Drosophila, Prmt5 acts as a type II arginine methyltransferase that catalyzes the symmetrical dimethylation of arginine residues in the RNA helicase Vasa. Prmt5 is critical for the regulation of spermatogenesis, but Vasa is not. To date, functional genetic studies of spermatogenesis in the lepidopteran model Bombyx mori has been limited. In this study, we engineered mutations in BmPrmt5 and BmVasa through CRISPR/Cas9-based gene editing. Both BmPrmt5 and BmVasa loss-of-function mutants had similar male and female sterility phenotypes. Through immunofluorescence staining analysis, we found that the morphs of sperm from both BmPrmt5 and BmVasa mutants have severe defects, indicating essential roles for both BmPrmt5 and BmVasa in the regulation of spermatogenesis. Mass spectrometry results identified that R35, R54, and R56 of BmVasa were dimethylated in WT while unmethylated in BmPrmt5 mutants. RNA-seq analyses indicate that the defects in spermatogenesis in mutants resulted from reduced expression of the spermatogenesis-related genes, including BmSxl, implying that BmSxl acts downstream of BmPrmt5 and BmVasa to regulate apyrene sperm development. These findings indicate that BmPrmt5 and BmVasa constitute an integral regulatory module essential for spermatogenesis in B. mori.

Genome sequence and silkomics of the spindle ermine moth, Yponomeuta cagnagella, representing the early diverging lineage of the ditrysian Lepidoptera.

Many lepidopteran species produce silk, cocoons, feeding tubes, or nests for protection from predators and parasites for caterpillars and pupae. Yet, the number of lepidopteran species whose silk composition has been studied in detail is very small, because the genes encoding the major structural silk proteins tend to be large and repetitive, making their assembly and sequence analysis difficult. Here we have analyzed the silk of Yponomeuta cagnagella, which represents one of the early diverging lineages of the ditrysian Lepidoptera thus improving the coverage of the order. To obtain a comprehensive list of the Y. cagnagella silk genes, we sequenced and assembled a draft genome using Oxford Nanopore and Illumina technologies. We used a silk-gland transcriptome and a silk proteome to identify major silk components and verified the tissue specificity of expression of individual genes. A detailed annotation of the major genes and their putative products, including their complete sequences and exon-intron structures is provided. The morphology of silk glands and fibers are also shown. This study fills an important gap in our growing understanding of the structure, evolution, and function of silk genes and provides genomic resources for future studies of the chemical ecology of Yponomeuta species.

BmPMFBP1 regulates the development of eupyrene sperm in the silkworm, Bombyx mori.

Sperm deliver the male complement of DNA to the ovum, and thus play a key role in sexual reproduction. Accordingly, spermatogenesis has outstanding significance in fields as disparate as infertility treatments and pest-control, making it a broadly interesting and important focus for molecular genetics research in a wide range of species. Here we investigate spermatogenesis in the model lepidopteran insect Bombyx mori (silkworm moth), with particular focus on the gene PMFBP1 (polyamine modulated factor 1 binding protein 1). In humans and mouse, PMFBP1 is essential for spermatogenesis, and mutations of this gene are associated with acephalic spermatozoa, which cause infertility. We identified a B. mori gene labeled as "PMFBP1" in GenBank's RefSeq database and sought to assess its role in spermatogenesis. Like in mammals, the silkworm version of this gene (BmPMFBP1) is specifically expressed in testes. We subsequently generated BmPMFBP1 mutants using a transgenic CRISPR/Cas9 system. Mutant males were sterile while the fertility of mutant females was comparable to wildtype females. In B. mori, spermatogenesis yields two types of sperm, the nucleated fertile eupyrene sperm, and anucleated unfertile apyrene sperm. Mutant males produced abnormal eupyrene sperm bundles but normal apyrene sperm bundles. For eupyrene sperm, nuclei were mislocated and disordered inside the bundles. We also found the BmPMFBP1 deficiency blocked the release of eupyrene sperm bundles from testes to ejaculatory seminalis. We found no obvious abnormalities in the production of apyrene sperm in mutant males, and double-matings with apyrene-deficient sex-lethal mutants rescued the ΔBmPMFBP1 infertility phenotype. These results indicate BmPMFBP1 functions only in eupyrene spermatogenesis, and highlight that distinct genes underlie the development of the two sperm morphs commonly found in Lepidoptera. Bioinformatic analyses suggest PMFBP1 may have evolved independently in lepidoptera and mammals, and that despite the shared name, are likely not homologous genes.

Sex-specific aging in animals: Perspective and future directions.

Sex differences in aging occur in many animal species, and they include sex differences in lifespan, in the onset and progression of age-associated decline, and in physiological and molecular markers of aging. Sex differences in aging vary greatly across the animal kingdom. For example, there are species with longer-lived females, species where males live longer, and species lacking sex differences in lifespan. The underlying causes of sex differences in aging remain mostly unknown. Currently, we do not understand the molecular drivers of sex differences in aging, or whether they are related to the accepted hallmarks or pillars of aging or linked to other well-characterized processes. In particular, understanding the role of sex-determination mechanisms and sex differences in aging is relatively understudied. Here, we take a comparative, interdisciplinary approach to explore various hypotheses about how sex differences in aging arise. We discuss genomic, morphological, and environmental differences between the sexes and how these relate to sex differences in aging. Finally, we present some suggestions for future research in this area and provide recommendations for promising experimental designs.

Support for faster and more adaptive Z chromosome evolution in two divergent lepidopteran lineages.

The rateof divergence for Z or X chromosomes is usually observed to be greater than autosomes, but the proposed evolutionary causes for this pattern vary, as do empirical results from diverse taxa. Even among moths and butterflies (Lepidoptera), which generally share a single-origin Z chromosome, the handful of available studies give mixed support for faster or more adaptive evolution of the Z chromosome, depending on the species assayed. Here, we examine the molecular evolution of Z chromosomes in two additional lepidopteran species: the Carolina sphinx moth and the monarch butterfly, the latter of which possesses a recent chromosomal fusion yielding a segment of newly Z-linked DNA. We find evidence for both faster and more adaptive Z chromosome evolution in both species, although this effect is strongest in the neo-Z portion of the monarch sex chromosome. The neo-Z is less male-biased than expected of a Z chromosome, and unbiased and female-biased genes drive the signal for adaptive evolution here. Together these results suggest that male-biased gene accumulation and haploid selection have opposing effects on long-term rates of adaptation and may help explain the discrepancies in previous findings as well as the repeated evolution of neo-sex chromosomes in Lepidoptera.

Population genomics reveals variable patterns of immune gene evolution in monarch butterflies (Danaus plexippus).

Humoral and cellular immune responses provide animals with major defences against harmful pathogens. While it is often assumed that immune genes undergo rapid diversifying selection, this assumption has not been tested in many species. Moreover, it is likely that different classes of immune genes experience different levels of evolutionary constraint, resulting in varying selection patterns. We examined the evolutionary patterns for a set of 91 canonical immune genes of North American monarch butterflies (Danaus plexippus), using as an outgroup the closely related soldier butterfly (Danaus eresimus). As a comparison to these immune genes, we selected a set of control genes that were paired with each immune for approximate size and genomic location. As a whole, these immune genes had a significant but modest reduction in Tajima's D relative to paired-control genes, but otherwise did not show distinct patterns of population genetic variation or evolutionary rates. When further partitioning these immune genes into four functional classes (recognition, signalling, modulation, and effector), we found distinct differences among these groups. Relative to control genes, recognition genes exhibit increased nonsynonymous diversity and divergence, suggesting reduced constraints on evolution, and supporting the notion that coevolution with pathogens results in diversifying selection. In contrast, signalling genes showed an opposite pattern of reduced diversity and divergence, suggesting evolutionary constraints and conservation. Modulator and effector genes showed no statistical differences from controls. These results are consistent with patterns found in immune genes in fruit flies and Pieris butterflies, suggesting that consistent selective pressures on different classes of immune genes broadly govern the evolution of innate immunity among insects.

Correction: Gtsf1 is essential for proper female sex determination and transposon silencing in the silkworm, Bombyx mori.

[This corrects the article DOI: 10.1371/journal.pgen.1009194.].

Gtsf1 is essential for proper female sex determination and transposon silencing in the silkworm, Bombyx mori.

Sex determination pathways are astoundingly diverse in insects. For instance, the silk moth Bombyx mori uniquely use various components of the piRNA pathway to produce the Fem signal for specification of the female fate. In this study, we identified BmGTSF1 as a novel piRNA factor which participates in B. mori sex determination. We found that BmGtsf1 has a distinct expression pattern compared to Drosophila and mouse. CRISPR/Cas9 induced mutation in BmGtsf1 resulted in partial sex reversal in genotypically female animals by shifting expression of the downstream targets BmMasc and Bmdsx to the male pattern. As levels of Fem piRNAs were substantially reduced in female mutants, we concluded that BmGtsf1 plays a critical role in the biogenesis of the feminizing signal. We also demonstrated that BmGTSF1 physically interacted with BmSIWI, a protein previously reported to be involved in female sex determination, indicating BmGTSF1 function as the cofactor of BmSIWI. BmGtsf1 mutation resulted in piRNA pathway dysregulation, including piRNA biogenesis defects and transposon derepression, suggesting BmGtsf1 is also a piRNA factor in the silkworm. Furthermore, we found that BmGtsf1 mutation leads to gametogenesis defects in both male and female. Our data suggested that BmGtsf1 is a new component involved in the sex determination pathway in B. mori.

A chromosome-level genome assembly provides new insights into paternal genome elimination in the cotton mealybug Phenacoccus solenopsis.

Mealybugs (Hemiptera: Pseudococcidae) are economically important agricultural pests with several compelling biological phenomena including paternal genome elimination (PGE). However, limited high-quality genome assemblies of mealybugs hinder a full understanding of this striking and unusual biological phenomenon. Here, we generated a chromosome-level genome assembly of cotton mealybug, Phenacoccus solenopsis, by combining Illumina short reads, PacBio long reads and Hi-C scaffolding. The assembled genome was 292.54 Mb with a contig N50 of 489.8 kb and a scaffold N50 of 49.0 Mb. Hi-C scaffolding assigned 84.42% of the bases to five chromosomes. A total of 110.75 Mb (37.9%) repeat sequences and 11,880 protein-coding genes were predicted. The completeness of the genome assembly was estimated to be 95.5% based on BUSCO genes. In addition, 27,086 (95.3%) full-length PacBio transcripts were uniquely mapped to the assembled scaffolds, suggesting the high quality of the genome assembly. We showed that cotton mealybugs lack differentiated sex chromosomes by analysing genome resequencing data of males and females. DAPI staining confirmed that one chromosome set in males becomes heterochromatin at an early embryo stage. Chromatin immunoprecipitation assays with sequencing analysis demonstrated that the epigenetic modifications H3K9me3 and H3K27me3 are distributed across the whole genome in males, suggesting that these two modifications might be involved in maintaining heterochromatin status. Both markers were more likely to be distributed in repeat regions, while H3K27me3 had higher overall enrichment. Our results provide a valuable genomic resource and shed new light on the genomic and epigenetic basis of PGE in cotton mealybugs.

Dichotomy of Dosage Compensation along the Neo Z Chromosome of the Monarch Butterfly.

Mechanisms of sex chromosome dosage compensation (SCDC) differ strikingly among animals. In Drosophila flies, chromosome-wide transcription is doubled from the single X chromosome in hemizygous (XY) males, whereas in Caenorhabditis nematodes, expression is halved for both X copies in homozygous (XX) females [1, 2]. Unlike other female-heterogametic (WZ female and ZZ male) animals, moths and butterflies exhibit sex chromosome dosage compensation patterns typically seen only in male-heterogametic species [3]. The monarch butterfly carries a newly derived Z chromosome segment that arose from an autosomal fusion with the ancestral Z [4]. Using a highly contiguous genome assembly, we show that gene expression is balanced between sexes along the entire Z chromosome but with distinct modes of compensation on the two segments. On the ancestral Z segment, depletion of H4K16ac corresponds to nearly halving of biallelic transcription in males, a pattern convergent to nematodes. Conversely, the newly derived Z segment shows a Drosophila-like mode of compensation, with enriched H4K16ac levels corresponding to doubled monoallelic transcription in females. Our work reveals that, contrary to the expectation of co-opting regulatory mechanisms readily in place, the evolution of plural modes of dosage compensation is also possible along a single sex chromosome within a species.

Transcriptomics of monarch butterflies (Danaus plexippus) reveals that toxic host plants alter expression of detoxification genes and down-regulate a small number of immune genes.

Herbivorous insects have evolved many mechanisms to overcome plant chemical defences, including detoxification and sequestration. Herbivores may also use toxic plants to reduce parasite infection. Plant toxins could directly interfere with parasites or could enhance endogenous immunity. Alternatively, plant toxins could favour down-regulation of endogenous immunity by providing an alternative (exogenous) defence against parasitism. However, studies on genomewide transcriptomic responses to plant defences and the interplay between plant toxicity and parasite infection remain rare. Monarch butterflies (Danaus plexippus) are specialist herbivores of milkweeds (Asclepias spp.), which contain toxic cardenolides. Monarchs have adapted to cardenolides through multiple resistance mechanisms and can sequester cardenolides to defend against bird predators. In addition, high-cardenolide milkweeds confer monarch resistance to a specialist protozoan parasite (Ophryocystis elektroscirrha). We used this system to study the interplay between the effects of plant toxicity and parasite infection on global gene expression. We compared transcriptional profiles between parasite-infected and uninfected monarch larvae reared on two milkweed species. Our results demonstrate that monarch differentially express several hundred genes when feeding on A. curassavica and A. incarnata, two species that differ substantially in cardenolide concentrations. These differentially expressed genes include genes within multiple families of canonical insect detoxification genes, suggesting that they play a role in monarch toxin resistance and sequestration. Interestingly, we found little transcriptional response to infection. However, parasite growth was reduced in monarchs reared on A. curassavica, and in these monarchs, several immune genes were down-regulated, consistent with the hypothesis that medicinal plants can reduce reliance on endogenous immunity.

Evolutionary Proteomics Reveals Distinct Patterns of Complexity and Divergence between Lepidopteran Sperm Morphs.

Spermatozoa are one of the most strikingly diverse animal cell types. One poorly understood example of this diversity is sperm heteromorphism, where males produce multiple distinct morphs of sperm in a single ejaculate. Typically, only one morph is capable of fertilization and the function of the nonfertilizing morph, called parasperm, remains to be elucidated. Sperm heteromorphism has multiple independent origins, including Lepidoptera (moths and butterflies), where males produce a fertilizing eupyrene sperm and an apyrene parasperm, which lacks a nucleus and nuclear DNA. Here we report a comparative proteomic analysis of eupyrene and apyrene sperm between two distantly related lepidopteran species, the monarch butterfly (Danaus plexippus) and Carolina sphinx moth (Manduca sexta). In both species, we identified ∼700 sperm proteins, with half present in both morphs and the majority of the remainder observed only in eupyrene sperm. Apyrene sperm thus have a distinctly less complex proteome. Gene ontology (GO) analysis revealed proteins shared between morphs tend to be associated with canonical sperm cell structures (e.g., flagellum) and metabolism (e.g., ATP production). GO terms for morph-specific proteins broadly reflect known structural differences, but also suggest a role for apyrene sperm in modulating female neurobiology. Comparative analysis indicates that proteins shared between morphs are most conserved between species as components of sperm, whereas morph-specific proteins turn over more quickly, especially in apyrene sperm. The rapid divergence of apyrene sperm content is consistent with a relaxation of selective constraints associated with fertilization and karyogamy. On the other hand, parasperm generally exhibit greater evolutionary lability, and our observations may therefore reflect adaptive responses to shifting regimes of sexual selection.

Nonfertilizing sperm in Lepidoptera show little evidence for recurrent positive selection.

Sperm are among the most variable cells in nature. Some of this variation results from nonadaptive errors in spermatogenesis, but many species consistently produce multiple sperm morphs, the adaptive significance of which remains unknown. Here, we investigate the evolution of dimorphic sperm in Lepidoptera, the butterflies and moths. Males of this order produce both fertilizing sperm and a secondary, nonfertilizing type that lacks DNA. Previous organismal studies suggested a role for nonfertilizing sperm in sperm competition, but this hypothesis has never been evaluated from a molecular framework. We combined published data sets with new sequencing in two species, the monandrous Carolina sphinx moth and the highly polyandrous monarch butterfly. Based on population genetic analyses, we see evidence for increased adaptive evolution in fertilizing sperm, but only in the polyandrous species. This signal comes primarily from a decrease in nonsynonymous polymorphism in sperm proteins compared to the rest of the genome, suggesting stronger purifying selection, consistent with selection via sperm competition. Nonfertilizing sperm proteins, in contrast, do not show an effect of mating system and do not appear to evolve differently from the background genome in either species, arguing against the involvement of nonfertilizing sperm in direct sperm competition. Based on our results and previous work, we suggest that nonfertilizing sperm may be used to delay female remating in these insects and decrease the risk of sperm competition rather than directly affect its outcome.

The genomic features of parasitism, Polyembryony and immune evasion in the endoparasitic wasp Macrocentrus cingulum.

BACKGROUND: Parasitoid wasps are well-known natural enemies of major agricultural pests and arthropod borne diseases. The parasitoid wasp Macrocentrus cingulum (Hymenoptera: Braconidae) has been widely used to control the notorious insect pests Ostrinia furnacalis (Asian Corn Borer) and O. nubilalis (European corn borer). One striking phenomenon exhibited by M. cingulum is polyembryony, the formation of multiple genetically identical offspring from a single zygote. Moreover, M. cingulum employs a passive parasitic strategy by preventing the host's immune system from recognizing the embryo as a foreign body. Thus, the embryos evade the host's immune system and are not encapsulated by host hemocytes. Unfortunately, the mechanism of both polyembryony and immune evasion remains largely unknown. RESULTS: We report the genome of the parasitoid wasp M. cingulum. Comparative genomics analysis of M. cingulum and other 11 insects were conducted, finding some gene families with apparent expansion or contraction which might be linked to the parasitic behaviors or polyembryony of M. cingulum. Moreover, we present the evidence that the microRNA miR-14b regulates the polyembryonic development of M. cingulum by targeting the c-Myc Promoter-binding Protein 1 (MBP-1), histone-lysine N-methyltransferase 2E (KMT2E) and segmentation protein Runt. In addition, Hemomucin, an O-glycosylated transmembrane protein, protects the endoparasitoid wasp larvae from being encapsulated by host hemocytes. Motif and domain analysis showed that only the hemomucin in two endoparasitoids, M. cingulum and Venturia canescens, possessing the ability of passive immune evasion has intact mucin domain and similar O-glycosylation patterns, indicating that the hemomucin is a key factor modulating the immune evasion. CONCLUSIONS: The microRNA miR-14b participates in the regulation of polyembryonic development, and the O-glycosylation of the mucin domain in the hemomucin confers the passive immune evasion in this wasp. These key findings provide new insights into the polyembryony and immune evasion.

The Z chromosome is enriched for sperm proteins in two divergent species of Lepidoptera.

Genes that promote sexual conflict, such as those with a sex-limited fitness benefit, are expected to accumulate differentially on sex chromosomes relative to autosomes. Few tests of this hypothesis exist for male homogametic (ZZ) taxa, however, and most use RNA expression data to identify such genes. Here, we employ a different identification method by using proteomic analysis of sperm cells to identify genes with a sex-limited benefit. We tested for a bias in genomic location of sperm protein genes in two species of Lepidoptera. An excess of sperm protein genes was identified on the Z chromosomes of both the Carolina sphinx moth (Manduca sexta) and the monarch butterfly (Danaus plexippus). Taking into consideration a Z-autosome fusion in monarchs, we discover that the ancestrally sex-linked portion of the genome is the source of this enrichment, while the newly sex-linked portion still appears similar to autosomes in relative abundance of sperm protein genes. Together, these results point to an enrichment of male-beneficial genes on the Z chromosome and demonstrate the usefulness of proteomic datasets in sexual conflict research.

Contrasting patterns of evolutionary constraint and novelty revealed by comparative sperm proteomic analysis in Lepidoptera.

BACKGROUND: Rapid evolution is a hallmark of reproductive genetic systems and arises through the combined processes of sequence divergence, gene gain and loss, and changes in gene and protein expression. While studies aiming to disentangle the molecular ramifications of these processes are progressing, we still know little about the genetic basis of evolutionary transitions in reproductive systems. Here we conduct the first comparative analysis of sperm proteomes in Lepidoptera, a group that exhibits dichotomous spermatogenesis, in which males produce a functional fertilization-competent sperm (eupyrene) and an incompetent sperm morph lacking nuclear DNA (apyrene). Through the integrated application of evolutionary proteomics and genomics, we characterize the genomic patterns potentially associated with the origination and evolution of this unique spermatogenic process and assess the importance of genetic novelty in Lepidopteran sperm biology. RESULTS: Comparison of the newly characterized Monarch butterfly (Danaus plexippus) sperm proteome to those of the Carolina sphinx moth (Manduca sexta) and the fruit fly (Drosophila melanogaster) demonstrated conservation at the level of protein abundance and post-translational modification within Lepidoptera. In contrast, comparative genomic analyses across insects reveals significant divergence at two levels that differentiate the genetic architecture of sperm in Lepidoptera from other insects. First, a significant reduction in orthology among Monarch sperm genes relative to the remainder of the genome in non-Lepidopteran insect species was observed. Second, a substantial number of sperm proteins were found to be specific to Lepidoptera, in that they lack detectable homology to the genomes of more distantly related insects. Lastly, the functional importance of Lepidoptera specific sperm proteins is broadly supported by their increased abundance relative to proteins conserved across insects. CONCLUSIONS: Our results identify a burst of genetic novelty amongst sperm proteins that may be associated with the origin of heteromorphic spermatogenesis in ancestral Lepidoptera and/or the subsequent evolution of this system. This pattern of genomic diversification is distinct from the remainder of the genome and thus suggests that this transition has had a marked impact on lepidopteran genome evolution. The identification of abundant sperm proteins unique to Lepidoptera, including proteins distinct between specific lineages, will accelerate future functional studies aiming to understand the developmental origin of dichotomous spermatogenesis and the functional diversification of the fertilization incompetent apyrene sperm morph.