Search

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 rate of 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, though 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. This article is protected by copyright. All rights reserved.

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.

OBJECTIVES To determine the proportion of solitary rib lesions on pre-treatment 68 Gallium-labelled prostate-specific membrane antigen (PSMA)/computed tomography (CT) scans in men with prostate cancer that are malignant and examine any predictive factors. PATIENTS AND METHODS This retrospective single tertiary referral institution cohort study of men reviewed the results of 68 Ga-PSMA-11 positron emission tomography (PET)/CT scans performed for primary staging prior to treatment of prostate cancer from July 2014 to September 2019. Men with PSMA uptake outside the prostate in only the rib lesion were included. A solitary rib lesion was considered to be malignant if it increased in size on follow-up imaging. A lesion was considered benign if the prostate-specific antigen (PSA) level remained

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.

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.

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.

In the wake of constant improvements in sequencing technologies, numerous insect genomes have been sequenced. Currently, 1219 insect genome-sequencing projects have been registered with the National Center for Biotechnology Information, including 401 that have genome assemblies and 155 with an official gene set of annotated protein-coding genes. Comparative genomics analysis showed that the expansion or contraction of gene families was associated with well-studied physiological traits such as immune system, metabolic detoxification, parasitism and polyphagy in insects. Here, we summarize the progress of insect genome sequencing, with an emphasis on how this impacts research on pest control. We begin with a brief introduction to the basic concepts of genome assembly, annotation and metrics for evaluating the quality of draft assemblies. We then provide an overview of genome information for numerous insect species, highlighting examples from prominent model organisms, agricultural pests and disease vectors. We also introduce the major insect genome databases. The increasing availability of insect genomic resources is beneficial for developing alternative pest control methods. However, many opportunities remain for developing data-mining tools that make maximal use of the available insect genome resources. Although rapid progress has been achieved, many challenges remain in the field of insect genomics.

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.

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 non-fertilizing 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 (Danausplexippus) and Carolina sphinx month (Manduca sexta). In both species, we identified approximatey 700 sperm proteins, with half present in both morphs and the majority of the remainder specific to 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, while 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 exhibit greater evolutionary lability, which may reflect adaptive response to shifting regimes of sexual selection. Additionally, we provide the first (to our knowledge) scanning electron micrographs of lepidopteran sperm.

Secondary spontaneous pneumothorax (SSP) is traditionally managed with an intercostal chest tube attached to an underwater seal. We investigated whether use of a one-way flutter valve shortened patients’ length of stay (LoS).This open-label randomised controlled trial enrolled patients presenting with SSP and randomised to either a chest tube and underwater seal (standard care: SC) or ambulatory care (AC) with a flutter valve. The type of flutter valve used depended on whether at randomisation the patient already had a chest tube in place: in those without a chest tube a pleural vent (PV) was used; in those with a chest tubein situ, an Atrium Pneumostat (AP) valve was attached. The primary end-point was LoS.Between March 2017 and March 2020, 41 patients underwent randomisation: 20 to SC and 21 to AC (13=PV, 8=AP). There was no difference in LoS in the first 30 days following treatment intervention: AC (median=6 days, IQR 14.5) and SC (median=6 days, IQR 13.3). In patients treated with PV there was a high rate of early treatment failure (6/13; 46%), compared to patients receiving SC (3/20; 15%) (p=0.11) Patients treated with AP had no (0/8 0%) early treatment failures and a median LoS of 1.5 days (IQR 23.8).There was no difference in LoS between ambulatory and standard care. Pleural Vents had high rates of treatment failure and should not be used in SSP. Atrium Pneumostats are a safer alternative, with a trend towards lower LoS.

The rate of divergence for Z or X chromosomes is theoretically predicted to be greater than autosomes, but the possible explanations 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, though 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.

Staging/restaging of prostate cancer utilizing Gallium-68 (68Ga) prostate-specific membrane antigen (PSMA) positron emission tomography/computed tomography (PET/CT) in combination with an intravenous urogram allows improved discrimination between radiotracer activity in the renal tract and small pelvic nodes or local recurrences. Within this pictorial essay, we describe the imaging protocol utilized at our institution and present cases which demonstrate the utility of this combined imaging approach.

Rationale, aims and objectives Public health systems are under pressure to meet increasing demand for health care in environments of increasing financial resource constraint. There is therefore a need to maximise health outcomes given limited public healthcare expenditure. This paper aims to establish the extent of literature and approaches to efficiency improvement in public health systems of developed countries. Methods The Rapid Evidence Assessment model was used to address the review question. Two database searches returned a result of 3,526 unique titles, which were individually screened for potential relevance. 144 titles were selected for full review to determine relevance. 73 papers were included in the final review. Results Data on country, study design, key findings and links to efficiency improvement were extracted and synthesized. Synthesis of findings revealed that the literature on this topic is disparate and non-cohesive. A range of isolated approaches were described, and no evidence or consensus on a single best-practice approach to efficiency improvement was identified. Conclusions Combining the factors identified in this review has the potential to inform a framework for supporting efficiency improvement in public health systems. By considering these factors central health system management bodies can support efficiency improvement to deliver both financial and health services benefits.

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 insectBombyx mori(silkworm moth), with particular focus on the genePMFBP1(polyamine modulated factor 1 binding protein 1). In humans and mouse,PMFBP1is essential for spermatogenesis, and mutations of this gene are associated with acephalic spermatozoa, which cause infertility. We identified aB.morigene 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 generatedBmPMFBP1mutants using a transgenic CRISPR/Cas9 system. Mutant males were sterile while the fertility of mutant females was comparable to wildtype females. InB.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 theBmPMFBP1deficiency 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-deficientsex-lethalmutants rescued theΔBmPMFBP1infertility phenotype. These results indicateBmPMFBP1functions only in eupyrene spermatogenesis, and highlight that distinct genes underlie the development of the two sperm morphs commonly found in Lepidoptera. Bioinformatic analyses suggestPMFBP1may have evolved independently in lepidoptera and mammals, and that despite the shared name, are likely not homologous genes.

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. Electronic supplementary material The online version of this article (10.1186/s12864-018-4783-x) contains supplementary material, which is available to authorized users.

Background Malignant pleural effusion affects more than 750,000 persons each year across Europe and the United States. Pleurodesis with the administration of talc in hospitalized patients is the most common treatment, but indwelling pleural catheters placed for drainage offer an ambulatory alternative. We examined whether talc administered through an indwelling pleural catheter was more effective at inducing pleurodesis than the use of an indwelling pleural catheter alone. Methods Over a period of 4 years, we recruited patients with malignant pleural effusion at 18 centers in the United Kingdom. After the insertion of an indwelling pleural catheter, patients underwent drainage regularly on an outpatient basis. If there was no evidence of substantial lung entrapment (nonexpandable lung, in which lung expansion and pleural apposition are not possible because of visceral fibrosis or bronchial obstruction) at 10 days, patients were randomly assigned to receive either 4 g of talc slurry or placebo through the indwelling pleural catheter on an outpatient basis. Talc or placebo was administered on a single-blind basis. Follow-up lasted for 70 days. The primary outcome was successful pleurodesis at day 35 after randomization. Results The target of 154 patients undergoing randomization was reached after 584 patients were approached. At day 35, a total of 30 of 69 patients (43%) in the talc group had successful pleurodesis, as compared with 16 of 70 (23%) in the placebo group (hazard ratio, 2.20; 95% confidence interval, 1.23 to 3.92; P = 0.008). No significant between- group differences in effusion size and complexity, number of inpatient days, mortality, or number of adverse events were identified. No significant excess of blockages of the indwelling pleural catheter was noted in the talc group. Conclusions Among patients without substantial lung entrapment, the outpatient administration of talc through an indwelling pleural catheter for the treatment of malignant pleural effusion resulted in a significantly higher chance of pleurodesis at 35 days than an indwelling catheter alone, with no deleterious effects.

We report the discovery of a neo-sex chromosome in Monarch butterfly, Danaus plexippus, and several of its close relatives. Z-linked scaffolds in the D. plexippus genome assembly were identified via sex-specific differences in Illumina sequencing coverage. Additionally, a majority of the D. plexippus genome assembly was assigned to chromosomes based on counts of 1-to-1 orthologs relative to the butterfly Melitaea cinxia (with replication using two other lepidopteran species), in which genome scaffolds have been mapped to linkage groups. Sequencing-coverage based assessments of Z-linkage combined with homology based chromosomal assignments provided strong evidence for a Z-autosome fusion in the Danaus lineage, involving the autosome homologous to chromosome 21 in M. cinxia. Coverage analysis also identified three notable assembly errors resulting in chimeric Z-autosome scaffolds. Cytogenetic analysis further revealed a large W-chromosome that is partially euchromatic, consistent with being a neo-W chromosome. The discovery of a neo-Z and the provisional assignment of chromosome linkage for >90% of D. plexippus genes lays the foundation for novel insights concerning sex chromosome evolution in this female-heterogametic model species for functional and evolutionary genomics.

Where previously described, patterns of sex chromosome dosage compensation in the Lepidoptera (moths and butterflies) have several unusual characteristics. Other female-heterogametic (ZW/ZZ) species exhibit female Z-linked expression that is reduced compared with autosomal expression and male Z expression. In the Lepidoptera, however, Z expression typically appears balanced between sexes but overall reduced relative to autosomal expression, that is Z ≈ ZZ

Summary 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.

Immune genes presumably rapidly evolve as pathogens exert strong selection pressures on host defense, but the evolution of immune genes is also constrained by trade-offs with other biological functions and shaped by the environmental context. Thus, immune genes may exhibit complex evolutionary patterns, particularly when organisms disperse to or live in variable environments. We examined the evolutionary patterns of the full set of known canonical immune genes within and among populations of monarch butterflies (Danaus plexippus), and relative to a closely related species (D. gilippus). Monarchs represent a system with a known evolutionary history, in which North American monarchs dispersed to form novel populations across the world, providing an opportunity to explore the evolution of immunity in the light of population expansion into novel environments. By analyzing a whole-genome resequencing dataset across populations, we found that immune genes as a whole do not exhibit consistent patterns of selection, differentiation, or genetic variation, but that patterns are specific to functional classes. Species comparisons betweenD. plexippusandD. gilippusand analyses of monarch populations both revealed consistently low levels of genetic variation in signaling genes, suggesting conservation of these genes over evolutionary time. Modulation genes showed the opposite pattern, with signatures of relaxed selection across populations. In contrast, recognition and effector genes exhibited less consistent patterns. When focusing on genes with exceptionally strong signatures of selection or differentiation, we also found population-specific patterns, consistent with the hypothesis that monarch populations do not face uniform selection pressures with respect to immune function.

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.

Herbivorous insects have evolved many mechanisms to overcome plant chemical defenses, 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 favor down-regulation of endogenous immunity by providing an alternative (exogenous) defense against parasitism. However, studies on genome-wide transcriptomic responses to plant defenses and the interplay between host plant toxicity and parasite infection remain rare. Monarch butterflies (Danaus plexippus) are specialist herbivores that feed on 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 medicinal effects to monarchs against 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. Our results demonstrate that monarch larvae differentially express several hundred genes when feeding on A. curassavica and A. incarnata, two species that are similar in nutritional content but 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, a small number of immune genes were down-regulated, consistent with the hypothesis that medicinal plants can reduce reliance on endogenous immunity.

The codling moth Cydia pomonella, a major invasive pest of pome fruit, has spread around the globe in the last half century. We generated a chromosome-level scaffold assembly including the Z chromosome and a portion of the W chromosome. This assembly reveals the duplication of an olfactory receptor gene (OR3), which we demonstrate enhances the ability of C. pomonella to exploit kairomones and pheromones in locating both host plants and mates. Genome-wide association studies contrasting insecticide-resistant and susceptible strains identify hundreds of single nucleotide polymorphisms (SNPs) potentially associated with insecticide resistance, including three SNPs found in the promoter of CYP6B2. RNAi knockdown of CYP6B2 increases C. pomonella sensitivity to two insecticides, deltamethrin and azinphos methyl. The high-quality genome assembly of C. pomonella informs the genetic basis of its invasiveness, suggesting the codling moth has distinctive capabilities and adaptive potential that may explain its worldwide expansion., The codling moth, Cydia pomonella, is one of the major pests of pome fruit (apples and pears) and walnuts. Here, the authors sequence and analyze its genome, providing insights on olfactory and detoxification processes that may underlie its worldwide expansion.

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.

SUMMARYColour morphs associated with different physiological adaptations offer unique opportunities to study the genomic basis of otherwise elusive adaptive traits. These complex balanced polymorphisms are often controlled by groups of tightly linked genes, and understanding how these ‘supergenes’ evolve and are maintained is an active area of research in evolutionary biology (Schwanderet al. 2014, Thompson and Jiggins 2014). Within the Atlantic, the common murre (Uria aalge, a colonial seabird) displays a plumage colour dimorphism (‘bridled’ and ‘unbridled’) that seems to be associated with differences in thermal adaptation (Birkhead 1984; Reiertsenet al. 2012). The genes associated with bridling and how these genes affect thermal adaptation are unknown. Using whole genome resequencing, we investigate the genomic basis of differences in colouration and thermal adaptation between the two morphs, and how the association between the two traits is maintained despite random mating. We identify a 60 kb genomic region of high differentiation laying in the intergenic area amongst three candidate genes for colouration and thermal adaptation:retinoic acid receptor beta(RARB),thyroid hormone receptor beta(THRB), andnuclear receptor subfamily 1 group D member 2(NR1D2orRev-erbβ). Differentiation is due to a complex copy number variant (CNV) that suppresses recombination locally. We show that this CNV acts as a ‘supergene’ and maintain association between regulatory elements likely affecting gene expression of one or more of the identified candidate genes. Our analyses also provide insights into the origin of the dimorphism: while copy number proliferation in the unbridled haplotype was potentially mediated by transposable elements (TEs), the bridled haplotype seems to have introgressed from the more cold-adapted sister species, the thick-billed murre (U. lomvia). Our results highlight the role of copy number variants in adaptation, especially when association among traits is maintained in the face of gene flow. They also shed light into the molecular mechanisms of adaptive thermogenesis in birds, which is poorly understood.HighlightsDifferences in plumage colour in Atlantic common murres are associated with different thermal adaptationsA single region is highly differentiated between bridled and unbridled morphsA complex copy number variant in a non-coding region underlies the dimorphismTransposable elements and adaptive introgression from the thick-billed murre seem to explain the origin of the dimorphism

Sperm are among the most variable cells in nature. Both within and between species, sperm display more diversity in form than one would expect of an ostensibly single-function cell. Morphologically, some of this variation within species has been demonstrated to be non-adaptive, but for many species that consistently produce multiple sperm morphs, the significance of variable sperm remains unknown. Here, we investigate the molecular evolution of dimorphic sperm in Lepidoptera, the butterflies and moths; males of this order produce both fertilizing eupyrene sperm and a secondary apyrene type that lacks DNA. Based on population genetic analyses in two species, the monandrous Carolina sphinx moth and the highly polyandrous monarch butterfly, we see evidence for increased selection in fertilizing sperm, but only in the polyandrous species. This signal comes primarily from a decrease in non-synonymous polymorphism in sperm proteins compared to the rest of the genome, indicative of strong purifying selection. Additionally, rates of non-synonymous divergence are comparable between sperm genes and the rest of the genome, suggesting that many alleles reach fixation owing to positive selection as well. Investigation of the distribution of fitness effects of new non-synonymous mutations in monarch sperm confirms stronger selection on sperm proteins in monarchs, with very few neutral variants and a preponderance of deleterious and beneficial variants. These results mirror findings on sperm evolution in other taxa; increased sperm competition decreases within-population morphological variation. Our results suggest that sperm competition can be a powerful selective force at the sequence level as well.

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.

Many animals with genetic sex determination harbor heteromorphic sex chromosomes, where the heterogametic sex has half the gene dose of the homogametic sex. This imbalance, if reflected in the abundance of transcripts or proteins, has the potential to deleteriously disrupt interactions between X-linked and autosomal loci in the heterogametic sex. Classical theory predicts that molecular mechanisms will evolve to provide dosage compensation that recovers expression levels comparable to ancestral expression prior to sex chromosome divergence. Such dosage compensating mechanisms may also, secondarily, result in balanced sex-linked gene expression between males and females. However, numerous recent studies addressing sex chromosome dosage compensation (SCDC) in a diversity of animals have yielded a surprising array of patterns concerning dosage compensation in the heterogametic sex, as well as dosage balance between sexes. These results substantially contradict longstanding theory, catalyzing both novel perspectives and new approaches in dosage compensation research. In this review, we summarize the theory, analytical approaches, and recent results concerning evolutionary patterns of SCDC in animals. We also discuss methodological challenges and discrepancies encountered in this research, which often underlie conflicting results. Finally, we discuss what outstanding questions and opportunities exist for future research on SCDC.

Most speciation events probably occur gradually, without complete and immediate reproductive isolation, but the full extent of gene flow between diverging species has rarely been characterized on a genome-wide scale. Documenting the extent and timing of admixture between diverging species can clarify the role of geographic isolation in speciation. Here we use new methodology to quantify admixture at different stages of divergence in Heliconius butterflies, based on whole-genome sequences of 31 individuals. Comparisons between sympatric and allopatric populations of H. melpomene, H. cydno, and H. timareta revealed a genome-wide trend of increased shared variation in sympatry, indicative of pervasive interspecific gene flow. Up to 40% of 100-kb genomic windows clustered by geography rather than by species, demonstrating that a very substantial fraction of the genome has been shared between sympatric species. Analyses of genetic variation shared over different time intervals suggested that admixture between these species has continued since early in speciation. Alleles shared between species during recent time intervals displayed higher levels of linkage disequilibrium than those shared over longer time intervals, suggesting that this admixture took place at multiple points during divergence and is probably ongoing. The signal of admixture was significantly reduced around loci controlling divergent wing patterns, as well as throughout the Z chromosome, consistent with strong selection for Müllerian mimicry and with known Z-linked hybrid incompatibility. Overall these results show that species divergence can occur in the face of persistent and genome-wide admixture over long periods of time.

Reproductive proteins often diverge rapidly between species. This pattern is frequently attributed to postmating sexual selection. Heliconius butterflies offer a good opportunity to examine this hypothesis by contrasting patterns of reproductive protein evolution between clades with divergent mating systems. Pupal-mating Heliconius females typically mate only once, limiting opportunity for postmating sexual selection. In contrast, adult-mating females remate throughout life. Reproductive protein evolution is therefore predicted to be slower and show little evidence of positive selection in the pupal-mating clade. We examined this prediction by sequencing 18 seminal fluid protein genes from a dozen Heliconius species and a related outgroup. Two proteins exhibited dN/dS > 1, implicating positive selection in the rapid evolution of at least a few Heliconius seminal fluid proteins. However, contrary to predictions, the average evolutionary rate of seminal fluid proteins was greater among pupal-mating Heliconius. Based on these results, we suggest that positive selection and relaxed constraint can generate conflicting patterns of reproductive protein evolution between mating systems. As predicted, some loci may show elevated evolutionary rates in promiscuous taxa relative to monandrous taxa resulting from adaptations to postmating sexual selection. However, when monandry is derived (as in Heliconius), the opposite pattern may result from relaxed selective constraints.

Seminal fluid proteins (SFPs) directly influence a wide range of reproductive processes, including fertilization, sperm storage, egg production, and immune response. Like many other reproductive proteins, the molecular evolution of SFPs is generally characterized by rapid and frequently adaptive evolution. However, the evolutionary processes underlying this often-documented pattern have not yet been confidently determined. A robust understanding of the processes governing SFP evolution will ultimately require identifying SFPs and characterizing their evolution in many different taxa, often where only limited genomic resources are available. Here, we report the first comprehensive molecular genetic and evolutionary analysis of SFPs conducted in Lepidoptera (moths and butterflies). We have identified 51 novel SFPs from two species of Heliconius butterflies (Heliconius erato and Heliconius melpomene) by combining "indirect" bioinformatic and expression analyses of expressed sequence tags from male accessory gland and wing tissues with "direct" proteomic analyses of spermatophores. Proteomic analyses identified fewer SFPs than the indirect criteria but gave consistent results. Of 51 SFPs, 40 were identified in both species but fewer than half could be functionally annotated via similarity searches (Blast, IPRscan, etc.). The majority of annotated Heliconius SFPs were predicted to be chymotrypsins. Comparisons of Heliconius SFPs with those from fruit fly, mosquito, honeybee, and cricket suggest that gene turnover is high among these proteins and that SFPs are rarely conserved across insect orders. Pairwise estimates of evolutionary rates between SFPs and nonreproductive proteins show that, on average, Heliconius SFPs are evolving rapidly. At least one of these SFPs is evolving adaptively (dN/dS > 1), implicating a role for positive selection in this rapid evolution. This work establishes a strong precedent for future research on the causes and consequences of reproductive protein evolution in the Lepidoptera. Butterflies and moths have an extremely rich history of organismal research, which will provide an informative ecological context for further molecular evolutionary investigations.

Data on patterns of variation within hybrid zones, combined with studies of life history, mate choice, and hybrid performance, allow estimates of the contribution of different pre-zygotic and post-zygotic barriers to reproductive isolation. We examine the role of behavioural barriers to gene exchange in the maintenance of a hybrid zone between North American field crickets Gryllus firmus and Gryllus pennsylvanicus. We consider these barriers in the context of previous studies that documented temporal and ecological isolation and a one-way post-mating incompatibility (i.e. G. firmus females do not produce offspring when they mate only with heterospecific males). Based on no-choice mating experiments in the laboratory, we demonstrate strong behavioural pre-mating barriers between the two species, but no apparent fecundity or fertility costs for G. firmus females when they mate with both conspecific and heterospecific males. Furthermore, we show that G. firmus females do not discriminate between hybrids and conspecifics, whereas G. pennsylvanicus females do. This observation could explain the asymmetric allele introgression observed in the hybrid zone. We also document a failure of heterospecific males to induce normal oviposition in G. firmus females, which may be due to rapid evolution of accessory gland proteins and may serve as an additional barrier to gene exchange. © 2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2009, 97, 390–402.

In contrast to XX/XY species, Z-linked expression is overall reduced in female WZ/ZZ species compared to males or the autosomal expression. This pattern (Z

The evolution of heterogametic sex chromosomes is often—but not always—accompanied by the evolution of dosage compensating mechanisms that mitigate the impact of sex-specific gene dosage on levels of gene expression. One emerging view of this process is that such mechanisms may only evolve in male-heterogametic (XY) species but not in female-heterogametic (ZW) species, which will consequently exhibit “incomplete” sex chromosome dosage compensation. However, recent results suggest that at least some Lepidoptera (moths and butterflies) may prove to be an exception to this prediction. Studies in bombycoid moths indicate the presence of a chromosome-wide epigenetic mechanism that effectively balances Z chromosome gene expression between the sexes by reducing Z-linked expression in males. In contrast, strong sex chromosome dosage effects without any reduction in male Z-linked expression were previously reported in a pyralid moth, suggesting a lack of any such dosage compensating mechanism. Here we report an analysis of sex chromosome dosage compensation in Heliconius butterflies, sampling multiple individuals for several different adult tissues (head, abdomen, leg, mouth, and antennae). Methodologically, we introduce a novel application of linear mixed-effects models to assess dosage compensation, offering a unified statistical framework that can estimate effects specific to chromosome, to sex, and their interactions (i.e., a dosage effect). Our results show substantially reduced Z-linked expression relative to autosomes in both sexes, as previously observed in bombycoid moths. This observation is consistent with an increasing body of evidence that some lepidopteran species possess an epigenetic dosage compensating mechanism that reduces Z chromosome expression in males to levels comparable with females. However, this mechanism appears to be imperfect in Heliconius, resulting in a modest dosage effect that produces an average 5–20% increase in male expression relative to females on the Z chromosome, depending on the tissue. Thus our results in Heliconius reflect a mixture of previous patterns reported for Lepidoptera. In Heliconius, a moderate pattern of incomplete dosage compensation persists apparently despite the presence of an epigenetic dosage compensating mechanism. The chromosomal distributions of sex-biased genes show an excess of male-biased and a dearth of female-biased genes on the Z chromosome relative to autosomes, consistent with predictions of sexually antagonistic evolution.

We describe a simple and straightforward procedure for the purification and separation of apyrene and eupyrene forms of lepidopteran sperm. The procedure is generally applicable to both butterfly and moth species with results varying according to the relative amounts of sperm produced and size of sperm storage organs. The technique relies upon inherent differences between eupyene sperm bundles and free apyrene sperm morphology. These differences allow for separation of the sperm morphs by repeated "panning" of sperm bundles into the center of a plastic dish. The purified eupyrene sperm bundles can then be removed and apyrene sperm collected from the supernatant by centrifugation. Efficacy of the purification process was confirmed by light microscopy and gel electrophoresis of the resulting fractions. Both one- and two-dimensional gel electrophoresis identified significant protein differences between the fractions further suggesting that the panning procedure effectively separated eurpyrene from apyrene sperm. The panning procedure should provide a convenient and accessible technique for further studies of sperm biology in lepidopterans.

(Received 18 November 2003 and in revised form 27 April 2004 ) Summary Microsatellite primers are often developed in one species and used to assess neutral variability in related species. Such analyses may be confounded by ascertainment bias (i.e. a decline in amplification success and allelic variability with increasing genetic distance from the source of the microsatellites). In addition, other factors, such as the size of the microsatellite, whether it consists of perfect or interrupted tandem repeats, and whether it is autosomal or X-linked, can affect variation. To test the relative importance of these factors on microsatellite variation, we examine patterns of amplification and allelic diversity in 52 microsatellite loci amplified from five individuals in each of six populations of Cyrtodiopsis stalk-eyed flies that range from 2 . 2 % to 11 . 2% mitochondrial DNA sequence divergence from the population used for microsatellite development. We find that amplification success and most measures of allelic diversity declined with genetic distance from the source population, in some cases an order of magnitude faster than in birds or mammals. The median and range of the repeat array length did not decline with genetic distance. In addition, for loci on the X chromosome, we find evidence of lower observed heterozygosity compared with loci on autosomes. The differences in variability between X-linked and autosomal loci are not adequately explained by differences in effective population sizes of the chromosomes. We suggest, instead, that periodic selection events associated with X-chromosome meiotic drive, which is present in many of these populations, reduces X-linked variation.

While most adult Lepidoptera use flower nectar as their primary food source, butterflies in the genus Heliconius have evolved the novel ability to acquire amino acids from consuming pollen. Heliconius butterflies collect pollen on their proboscis, moisten the pollen with saliva, and use a combination of mechanical disruption and chemical degradation to release free amino acids that are subsequently re-ingested in the saliva. Little is known about the molecular mechanisms of this complex pollen feeding adaptation. Here we report an initial shotgun proteomic analysis of saliva from Heliconius melpomene. Results from liquid-chromatography tandem mass-spectrometry confidently identified 31 salivary proteins, most of which contained predicted signal peptides, consistent with extracellular secretion. Further bioinformatic annotation of these salivary proteins indicated the presence of four distinct functional classes: proteolysis (10 proteins), carbohydrate hydrolysis (5), immunity (6), and "housekeeping" (4). Additionally, six proteins could not be functionally annotated beyond containing a predicted signal sequence. The presence of several salivary proteases is consistent with previous demonstrations that Heliconius saliva has proteolytic capacity. It is likely that these proteins play a key role in generating free amino acids during pollen digestion. The identification of proteins functioning in carbohydrate hydrolysis is consistent with Heliconius butterflies consuming nectar, like other lepidopterans, as well as pollen. Immune-related proteins in saliva are also expected, given that ingestion of pathogens is a likely route to infection. The few "housekeeping" proteins are likely not true salivary proteins and reflect a modest level of contamination that occurred during saliva collection. Among the unannotated proteins were two sets of paralogs, each seemingly the result of a relatively recent tandem duplication. These results offer a first glimpse into the molecular foundation of Heliconius pollen feeding and provide a substantial advance towards comprehensively understanding this striking evolutionary novelty. publisher: Elsevier articletitle: Pollen feeding proteomics: Salivary proteins of the passion flower butterfly, Heliconius melpomene journaltitle: Insect Biochemistry and Molecular Biology articlelink: http://dx.doi.org/10.1016/j.ibmb.2015.04.004 content_type: article copyright: Copyright © 2015 Elsevier Ltd. All rights reserved. ispartof: Insect Biochemistry and Molecular Biology vol:63 pages:7-13 ispartof: location:England status: published

The application of mass spectrometry based proteomics to sperm biology has greatly accelerated progress in understanding the molecular composition and function of spermatozoa. To date, these approaches have been largely restricted to model organisms, all of which produce a single sperm morph capable of oocyte fertilisation. Here we apply high-throughput mass spectrometry proteomic analysis to characterise sperm composition in Manduca sexta, the tobacco hornworm moth, which produce heteromorphic sperm, including one fertilisation competent (eupyrene) and one incompetent (apyrene) sperm type. This resulted in the high confidence identification of 896 proteins from a co-mixed sample of both sperm types, of which 167 are encoded by genes with strict one-to-one orthology in Drosophila melanogaster. Importantly, over half (55.1%) of these orthologous proteins have previously been identified in the D. melanogaster sperm proteome and exhibit significant conservation in quantitative protein abundance in sperm between the two species. Despite the complex nature of gene expression across spermatogenic stages, a significant correlation was also observed between sperm protein abundance and testis gene expression. Lepidopteran-specific sperm proteins (e.g., proteins with no homology to proteins in non-Lepidopteran taxa) were present in significantly greater abundance on average than those with homology outside the Lepidoptera. Given the disproportionate production of apyrene sperm (96% of all mature sperm in Manduca) relative to eupyrene sperm, these evolutionarily novel and highly abundant proteins are candidates for possessing apyrene-specific functions. Lastly, comparative genomic analyses of testis-expressed, ovary-expressed and sperm genes identified a concentration of novel sperm proteins shared amongst Lepidoptera of potential relevance to the evolutionary origin of heteromorphic spermatogenesis. As the first published Lepidopteran sperm proteome, this whole-cell proteomic characterisation will facilitate future evolutionary genetic and developmental studies of heteromorphic sperm production and parasperm function. Furthermore, the analyses presented here provide useful annotation information regarding sex-biased gene expression, novel Lepidopteran genes and gene function in the male gamete to complement the newly sequenced and annotated Manduca genome.

Secondary plant compounds are strong deterrents of insect oviposition and feeding, but may also be attractants for specialist herbivores. These insect-plant interactions are mediated by insect gustatory receptors (Grs) and olfactory receptors (Ors). An analysis of the reference genome of the butterfly Heliconius melpomene, which feeds on passion-flower vines (Passiflora spp.), together with whole-genome sequencing within the species and across the Heliconius phylogeny has permitted an unprecedented opportunity to study the patterns of gene duplication and copy-number variation (CNV) among these key sensory genes. We report in silico gene predictions of 73 Gr genes in the H. melpomene reference genome, including putative CO2, sugar, sugar alcohol, fructose, and bitter receptors. The majority of these Grs are the result of gene duplications since Heliconius shared a common ancestor with the monarch butterfly or the silkmoth. Among Grs but not Ors, CNVs are more common within species in those gene lineages that have also duplicated over this evolutionary time-scale, suggesting ongoing rapid gene family evolution. Deep sequencing (∼1 billion reads) of transcriptomes from proboscis and labial palps, antennae, and legs of adult H. melpomene males and females indicates that 67 of the predicted 73 Gr genes and 67 of the 70 predicted Or genes are expressed in these three tissues. Intriguingly, we find that one-third of all Grs show female-biased gene expression (n = 26) and nearly all of these (n = 21) are Heliconius-specific Grs. In fact, a significant excess of Grs that are expressed in female legs but not male legs are the result of recent gene duplication. This difference in Gr gene expression diversity between the sexes is accompanied by a striking sexual dimorphism in the abundance of gustatory sensilla on the forelegs of H. melpomene, suggesting that female oviposition behaviour drives the evolution of new gustatory receptors in butterfly genomes., Author Summary Insects and their chemically-defended hostplants engage in a co-evolutionary arms race but the genetic basis by which suitable host plants are identified by insects is poorly understood. Host plant specializations require specialized sensors by the insects to exploit novel ecological niches. Adult male and female Heliconius butterflies feed on nectar and, unusually for butterflies, on pollen from flowers while their larvae feed on the leaves of passion-flower vines. We have discovered–between sub-species of butterflies-fixed differences in copy-number variation among several putative sugar receptor genes that are located on different chromosomes, raising the possibility of local adaptation around the detection of sugars. We also show that the legs of adult female butterflies, which are used by females when selecting a host plant on which to lay their eggs, express more gustatory (taste) receptor genes than those of male butterflies. These female-biased taste receptors show a significantly higher level of gene duplication than a set of taste receptors expressed in both sexes. Sex-limited behaviour may therefore influence the long-term evolution of physiologically important gene families resulting in a strong genomic signature of ecological adaptation.

The evolutionary importance of hybridization and introgression has long been debated. Hybrids are usually rare and unfit, but even infrequent hybridization can aid adaptation by transferring beneficial traits between species. Here we use genomic tools to investigate introgression in Heliconius, a rapidly radiating genus of neotropical butterflies widely used in studies of ecology, behaviour, mimicry and speciation. We sequenced the genome of Heliconius melpomene and compared it with other taxa to investigate chromosomal evolution in Lepidoptera and gene flow among multiple Heliconius species and races. Among 12, 669 predicted genes, biologically important expansions of families of chemosensory and Hox genes are particularly noteworthy. Chromosomal organization has remained broadly conserved since the Cretaceous period, when butterflies split from the Bombyx (silkmoth) lineage. Using genomic resequencing, we show hybrid exchange of genes between three co-mimics, Heliconius melpomene, Heliconius timareta and Heliconius elevatus, especially at two genomic regions that control mimicry pattern. We infer that closely related Heliconius species exchange protective colour-pattern genes promiscuously, implying that hybridization has an important role in adaptive radiation. © 2012 Macmillan Publishers Limited.

Reproductive proteins often diverge rapidly between species. This pattern is frequently attributed to postmating sexual selection. Heliconius butterflies offer a good opportunity to examine this hypothesis by contrasting patterns of reproductive protein evolution between clades with divergent mating systems. Pupal-mating Heliconius females typically mate only once, limiting opportunity for postmating sexual selection. In contrast, adult-mating females remate throughout life. Reproductive protein evolution is therefore predicted to be slower and show little evidence of positive selection in the pupal-mating clade. We examined this prediction by sequencing 18 seminal fluid protein genes from a dozen Heliconius species and a related outgroup. Two proteins exhibited dN/dS1, implicating positive selection in the rapid evolution of at least a few Heliconius seminal fluid proteins. However, contrary to predictions, the average evolutionary rate of seminal fluid proteins was greater among pupal-mating Heliconius. Based on these results, we suggest that positive selection and relaxed constraint can generate conflicting patterns of reproductive protein evolution between mating systems. As predicted, some loci may show elevated evolutionary rates in promiscuous taxa relative to monandrous taxa resulting from adaptations to postmating sexual selection. However, when monandry is derived (as in Heliconius), the opposite pattern may result from relaxed selective constraints.