Retrotransposable elements on the W chromosome of the silkworm, Bombyx mori

W chromosome sequences of two bombycid moths provide an insight into the origin of Fem

Fem is a W-linked gene that encodes a piRNA precursor, and its product, Fem piRNA, is a master factor of female determination in Bombyx mori. Fem has low similarity to any known sequences, and the origin of Fem remains unclear. So far, two hypotheses have been proposed for the origin of Fem: The first hypothesis is that Fem is an allele of Masc, which assumes that the W chromosome was originally a homologous chromosome of the Z chromosome. The second hypothesis is that Fem arose by the transposition of Masc to the W chromosome. To explore the origin of Fem, we determined the W chromosome sequences of B. mori and, as a comparison, a closely relative bombycid species of Trilocha varians with a Fem-independent sex determination system. To our surprise, although the sequences of W and Z chromosomes show no homology to each other, a few pairs of homologues are shared by W and Z chromosomes, indicating the W chromosome of both species originated from Z chromosome. In addition, the W chromosome of T. varians lacks Fem, while the W chromosome of B. mori has over 100 copies of Fem. The high-quality assembly of the W chromosome of B. mori arose the third hypothesis about the origin of Fem: Fem is a chimeric sequence of multiple transposons. More than half of one transcriptional unit of Fem shows a significant homology to RTE-BovB. Moreover, the Fem piRNA-producing region could correspond to the boundary of the two transposons, gypsy and satellite DNA.

Multiple independent origins of the female W chromosome in moths and butterflies

Lepidoptera, the most diverse group of insects, exhibit female heterogamy (Z0 or ZW), which is different from most other insects (male heterogamy, XY). Previous studies suggest a single origin of the Z chromosome. However, the origin of the lepidopteran W chromosome remains poorly understood. Here, we assemble the genome from females down to the chromosome level of a model insect (Bombyx mori) and identify a W chromosome of approximately 10.1 megabase using a newly developed tool. In addition, we identify 3593 genes that were not previously annotated in the genomes of B. mori. Comparisons of 21 lepidopteran species (including 17 ZW and four Z0 systems) and three trichopteran species (Z0 system) reveal that the formation of Ditrysia W involves multiple mechanisms, including previously proposed canonical and noncanonical models, as well as a newly proposed mechanism called single-Z turnover. We conclude that there are multiple independent origins of the W chromosome in the Ditrysia (most moths and all butterflies) of Lepidoptera.

Gene gain and loss from the Asian corn borer W chromosome

BACKGROUND

Sex-limited chromosomes Y and W share some characteristics, including the degeneration of protein-coding genes, enrichment of repetitive elements, and heterochromatin. However, although many studies have suggested that Y chromosomes retain genes related to male function, far less is known about W chromosomes and whether they retain genes related to female-specific function.

RESULTS

Here, we built a chromosome-level genome assembly of the Asian corn borer, Ostrinia furnacalis Guenée (Lepidoptera: Crambidae, Pyraloidea), an economically important pest in corn, from a female, including both the Z and W chromosome. Despite deep conservation of the Z chromosome across Lepidoptera, our chromosome-level W assembly reveals little conservation with available W chromosome sequence in related species or with the Z chromosome, consistent with a non-canonical origin of the W chromosome. The W chromosome has accumulated significant repetitive elements and experienced rapid gene gain from the remainder of the genome, with most genes exhibiting pseudogenization after duplication to the W. The genes that retain significant expression are largely enriched for functions in DNA recombination, the nucleosome, chromatin, and DNA binding, likely related to meiotic and mitotic processes within the female gonad.

CONCLUSIONS

Overall, our chromosome-level genome assembly supports the non-canonical origin of the W chromosome in O. furnacalis, which experienced rapid gene gain and loss, with the retention of genes related to female-specific function.

Accumulation of retrotransposons contributes to W chromosome differentiation in the willow beauty Peribatodes rhomboidaria (Lepidoptera: Geometridae)

The W chromosome of Lepidoptera is typically gene-poor, repeat-rich and composed of heterochromatin. Pioneering studies investigating this chromosome reported an abundance of mobile elements. However, the actual composition of the W chromosome varies greatly between species, as repeatedly demonstrated by comparative genomic hybridization (CGH) or fluorescence in situ hybridization (FISH). Here we present an analysis of repeats on the W chromosome in the willow beauty, Peribatodes rhomboidaria (Geometridae), a species in which CGH predicted an abundance of W-enriched or W-specific sequences. Indeed, comparative analysis of male and female genomes using RepeatExplorer identified ten putative W chromosome-enriched repeats, most of which are LTR or LINE mobile elements. We analysed the two most abundant: PRW LINE-like and PRW Bel-Pao. The results of FISH mapping and bioinformatic analysis confirmed their enrichment on the W chromosome, supporting the hypothesis that mobile elements are the driving force of W chromosome differentiation in Lepidoptera. As the W chromosome is highly underrepresented in chromosome-level genome assemblies of Lepidoptera, this recently introduced approach, combining bioinformatic comparative genome analysis with molecular cytogenetics, provides an elegant tool for studying this elusive and rapidly evolving part of the genome.

A step forward in the genome characterization of the sugarcane borer, Diatraea saccharalis: karyotype analysis, sex chromosome system and repetitive DNAs through a cytogenomic approach

Moths of the family Crambidae include a number of pests that cause economic losses to agricultural crops. Despite their economic importance, little is known about their genome architecture and chromosome evolution. Here, we characterized the chromosomes and repetitive DNA of the sugarcane borer Diatraea saccharalis using a combination of low-pass genome sequencing, bioinformatics, and cytogenetic methods, focusing on the sex chromosomes. Diploid chromosome numbers differed between the sexes, i.e., 2n = 33 in females and 2n = 34 in males. This difference was caused by the occurrence of a WZ₁Z₂ trivalent in female meiosis, indicating a multiple sex-chromosome system WZ₁Z₂/Z₁Z₁Z₂Z₂. A strong interstitial telomeric signal was observed on the W chromosome, indicating a fusion of the ancestral W chromosome with an autosome. Among repetitive DNAs, transposable elements (TEs) accounted for 39.18% (males) to 41.35% (females), while satDNAs accounted for only 0.214% (males) and 0.215% (females) of the genome. FISH mapping revealed different chromosomal organization of satDNAs, such as single localized clusters, spread repeats, and non-clustered repeats. Two TEs mapped by FISH were scattered. Although we found a slight enrichment of some satDNAs in the female genome, they were not differentially enriched on the W chromosome. However, we found enriched FISH signals for TEs on the W chromosome, suggesting their involvement in W chromosome degeneration and differentiation. These data shed light on karyotype and repetitive DNA dynamics due to multiple chromosome fusions in D. saccharalis, contribute to the understanding of genome structure in Lepidoptera and are important for future genomic studies.

Identification and Characterization of the Masculinizing Function of the Helicoverpa armigera Masc Gene

The Masculinizer (Masc) gene has been known to control sex development and dosage compensation in lepidopterans. However, it remains unclear whether its ortholog exists and plays the same roles in distantly related lepidopterans such as Helicoverpa armigera. To address this question, we cloned Masc from H. armigera (HaMasc), which contains all essential functional domains of BmMasc, albeit with less than 30% amino acid sequence identity with BmMasc. Genomic PCR and qPCR analyses showed that HaMasc is a Z chromosome-linked gene since its genomic content in males (ZZ) was two times greater than that in females (ZW). RT-PCR and RT-qPCR analyses revealed that HaMasc expression was sex- and stage-biased, with significantly more transcripts in males and eggs than in females and other stages. Transfection of a mixture of three siRNAs of HaMasc into a male embryonic cell line of H. armigera led to the appearance of female-specific mRNA splicing isoforms of H. armigeradoublesex (Hadsx), a downstream target gene of HaMasc in the H. armigera sex determination pathway. The knockdown of HaMasc, starting from the third instar larvae resulted in a shift of Hadsx splicing from male to female isoforms, smaller male pupa and testes, fewer but larger/longer spermatocytes and sperm bundles, delayed pupation and internal fusion of the testes and follicles. These data demonstrate that HaMasc functions as a masculinizing gene in the H. armigera sex-determination cascade.

The Dryas iulia Genome Supports Multiple Gains of a W Chromosome from a B Chromosome in Butterflies

In butterflies and moths, which exhibit highly variable sex determination mechanisms, the homogametic Z chromosome is deeply conserved and is featured in many genome assemblies. The evolution and origin of the female W sex chromosome, however, remains mostly unknown. Previous studies have proposed that a ZZ/Z0 sex determination system is ancestral to Lepidoptera, and that W chromosomes may originate from sex-linked B chromosomes. Here, we sequence and assemble the female Dryas iulia genome into 32 highly contiguous ordered and oriented chromosomes, including the Z and W sex chromosomes. We then use sex-specific Hi-C, ATAC-seq, PRO-seq, and whole-genome DNA sequence data sets to test if features of the D. iulia W chromosome are consistent with a hypothesized B chromosome origin. We show that the putative W chromosome displays female-associated DNA sequence, gene expression, and chromatin accessibility to confirm the sex-linked function of the W sequence. In contrast with expectations from studies of homologous sex chromosomes, highly repetitive DNA content on the W chromosome, the sole presence of domesticated repetitive elements in functional DNA, and lack of sequence homology with the Z chromosome or autosomes is most consistent with a B chromosome origin for the W, although it remains challenging to rule out extensive sequence divergence. Synteny analysis of the D. iulia W chromosome with other female lepidopteran genome assemblies shows no homology between W chromosomes and suggests multiple, independent origins of the W chromosome from a B chromosome likely occurred in butterflies.

Chromosome-level genome reference and genome editing of the tea geometrid

The tea geometrid is a destructive insect pest on tea plants, which seriously affects tea production in terms of both yield and quality and causes severe economic losses. The tea geometrid also provides an important study system to address the ecological adaptive mechanisms underlying its unique host plant adaptation and protective resemblance. In this study, we fully sequenced and de novo assembled the reference genome of the tea geometrid, Ectropis grisescens, using long sequencing reads. We presented a highly continuous, near-complete genome reference (787.4 Mb; scaffold N50: 26.9 Mb), along with the annotation of 18,746 protein-coding genes and 53.3% repeat contents. Importantly, we successfully placed 97.8% of the assembly in 31 chromosomes based on Hi-C interactions and characterized the sex chromosome based on sex-biased sequencing coverage. Multiple quality-control assays and chromosome-scale synteny with the model species all supported the high quality of the presented genome reference. We focused biological annotations on gene families related to the host plant adaptation and camouflage in the tea geometrid and performed comparisons with other representative lepidopteran species. Important findings include the E. grisescens-specific expansion of CYP6 P450 genes that might be involved in metabolism of tea defensive chemicals and unexpected massive expansion of gustatory receptor gene families that suggests potential polyphagy for this tea pest. Furthermore, we developed an efficient genome editing system based on CRISPR/Cas9 technology and successfully implement mutagenesis of a Hox gene in the tea geometrid. Our study provides key genomic resources both for exploring unique mechanisms underlying the ecological adaptation of tea geometrids and for developing environment-friendly strategies for tea pest management.

The Role of Satellite DNAs in Genome Architecture and Sex Chromosome Evolution in Crambidae Moths

Tandem repeats are important parts of eukaryotic genomes being crucial e.g., for centromere and telomere function and chromatin modulation. In Lepidoptera, knowledge of tandem repeats is very limited despite the growing number of sequenced genomes. Here we introduce seven new satellite DNAs (satDNAs), which more than doubles the number of currently known lepidopteran satDNAs. The satDNAs were identified in genomes of three species of Crambidae moths, namely Ostrinia nubilalis, Cydalima perspectalis, and Diatraea postlineella, using graph-based computational pipeline RepeatExplorer. These repeats varied in their abundance and showed high variability within and between species, although some degree of conservation was noted. The satDNAs showed a scattered distribution, often on both autosomes and sex chromosomes, with the exception of both satellites in D. postlineella, in which the satDNAs were located at a single autosomal locus. Three satDNAs were abundant on the W chromosomes of O. nubilalis and C. perspectalis, thus contributing to their differentiation from the Z chromosomes. To provide background for the in situ localization of the satDNAs, we performed a detailed cytogenetic analysis of the karyotypes of all three species. This comparative analysis revealed differences in chromosome number, number and location of rDNA clusters, and molecular differentiation of sex chromosomes.

The Sex Determination Cascade in the Silkworm

In insects, sex determination pathways involve three levels of master regulators: primary signals, which determine the sex; executors, which control sex-specific differentiation of tissues and organs; and transducers, which link the primary signals to the executors. The primary signals differ widely among insect species. In Diptera alone, several unrelated primary sex determiners have been identified. However, the doublesex (dsx) gene is highly conserved as the executor component across multiple insect orders. The transducer level shows an intermediate level of conservation. In many, but not all examined insects, a key transducer role is performed by transformer (tra), which controls sex-specific splicing of dsx. In Lepidoptera, studies of sex determination have focused on the lepidopteran model species Bombyx mori (the silkworm). In B. mori, the primary signal of sex determination cascade starts from Fem, a female-specific PIWI-interacting RNA, and its targeting gene Masc, which is apparently specific to and conserved among Lepidoptera. Tra has not been found in Lepidoptera. Instead, the B. mori PSI protein binds directly to dsx pre-mRNA and regulates its alternative splicing to produce male- and female-specific transcripts. Despite this basic understanding of the molecular mechanisms underlying sex determination, the links among the primary signals, transducers and executors remain largely unknown in Lepidoptera. In this review, we focus on the latest findings regarding the functions and working mechanisms of genes involved in feminization and masculinization in Lepidoptera and discuss directions for future research of sex determination in the silkworm.

A chromosome-level genome assembly of rice leaffolder, Cnaphalocrocis medinalis

The rice leaffolder, Cnaphalocrocis medinalis Guenée (Crambidae, Lepidoptera), is an important agricultural pest that causes serious losses to rice production in rice-growing regions with high humidity and temperature. However, a lack of genomic resources limits in-depth understanding of its biological characteristics and ecological adaptation. Here, we sequenced the genome of rice leaffolder using the Illumina and PacBio platforms, yielding a genome assembly of 528.3 Mb with a contig N50 of 524.6 kb. A high percentage (96.4%) of Benchmarking Universal Single-Copy Orthologs (BUSCOs) were successfully detected, suggesting high-level completeness of the genome assembly. In total, 39.5% of the genome consists of repeat sequences and 15,045 protein-coding genes were annotated. Comparative phylogenomic analysis showed that some gene families associated with hormone biosynthesis expanded in rice leaffolder. Next, we used the Hi-C technique to produce a chromosome-level genome assembly with a scaffold N50 of 16.1 Mb by anchoring 3,248 scaffolds to 31 chromosomes. The rice leaffolder genome showed high chromosomal synteny with the genome of four other lepidopteran insects. By comparing coverage ratios from the genome resequencing of male and female pupae, we identified near intact Z and W chromosomes. The W chromosome is estimated as 20.75 Mb, which is the most complete known W chromosome in Lepidoptera. The protein-coding genes on the W chromosome were significantly enriched in metabolic pathways. In all, the high-quality genome assembly and the near-intact W chromosome of rice leaffolder should be a useful resource for the fields of insect migration, chromosome evolution and pest control.

An Efficient Duplex PCR Method for Sex Identification of the European Grapevine Moth Lobesia botrana (Lepidoptera: Tortricidae) at Any Developmental Stage

Many genetic studies in insects require sex identification of individuals in all developmental stages. The most common sex chromosome system in lepidopterans is WZ/ZZ; the W chromosome is present only in females. Based on two W chromosome-specific short sequences (CpW2 and CpW5) described in Cydia pomonella (L.) (Lepidoptera: Tortricidae), we identified homologous female-specific sequences in Lobesia botrana Den. & Schiff, a polyphagous and very harmful species present in Chile since 2008. From this starting point, we extended the sequence information using the inverse PCR method, identifying the first W-specific sequences described up to now for the moth. Finally, we developed a duplex PCR method for rapid and sensitive determination of sex in L. botrana from larva to adult. The method showed a detection limit of 1 pg of genomic DNA; a blind panel of samples exhibited exact correspondence with the morphological identification. These results will be very useful for studies requiring sex-specific analyses at any developmental stage, contributing also to the understanding of gene expression in the insect, as well as to the eventual development of control protocols against the moth, such as the development of genetic sexing strains for the implementation of the sterile insect technique.

Characterization of the First W-Specific Protein-Coding Gene for Sex Identification in Helicoverpa armigera

Helicoverpa armigera is a globally-important crop pest with a WZ (female)/ZZ (male) sex chromosome system. The absence of discernible sexual dimorphism in its egg and larval stages makes it impossible to address any sex-related theoretical and applied questions before pupation unless a W-specific sequence marker is available for sex diagnosis. To this end, we used one pair of morphologically pre-sexed pupae to PCR-screen 17 non-transposon transcripts selected from 4855 W-linked candidate reads identified by mapping a publicly available egg transcriptome of both sexes to the male genome of this species and detected the read SRR1015458.67499 only in the female pupa. Subsequent PCR screenings of this read and the previously reported female-specific RAPD (random amplified polymorphic DNA) marker AF18 with ten more pairs of pre-sexed pupae and different annealing positions and/or temperatures as well as its co-occurrence with the female-specific transcript splicing isoforms of doublesex gene of H. armigera (Hadsx) and amplification and sequencing of their 5′ unknown flanking sequences in three additional pairs of pre-sexed pupae verified that SRR1015458.67499 is a single copy protein-coding gene unique to W chromosome (named GUW1) while AF18 is a multicopy MITE transposon located on various chromosomes. Test application of GUW1 as a marker to sex 30 neonates of H. armigera yielded a female/male ratio of 1.14: 1.00. Both GUW1 and Hadsx splicing isoforms assays revealed that the H. armigera embryo cell line QB-Ha-E-1 is a male cell line. Taken together, GUW1 is not only a reliable DNA marker for sexing all stages of H. armigera and its cell lines, but also represents the first W-specific protein-coding gene in lepidopterans.

Evolution of multiple sex-chromosomes associated with dynamic genome reshuffling in Leptidea wood-white butterflies

Sex-chromosome systems tend to be highly conserved and knowledge about their evolution typically comes from macroevolutionary inference. Rapidly evolving complex sex-chromosome systems represent a rare opportunity to study the mechanisms of sex-chromosome evolution at unprecedented resolution. Three cryptic species of wood-white butterflies—Leptidea juvernica, L. sinapis and L. reali—have each a unique set of multiple sex-chromosomes with 3–4 W and 3–4 Z chromosomes. Using a transcriptome-based microarray for comparative genomic hybridisation (CGH) and a library of bacterial artificial chromosome (BAC) clones, both developed in L. juvernica, we identified Z-linked Leptidea orthologs of Bombyx mori genes and mapped them by fluorescence in situ hybridisation (FISH) with BAC probes on multiple Z chromosomes. In all three species, we determined synteny blocks of autosomal origin and reconstructed the evolution of multiple sex-chromosomes. In addition, we identified W homologues of Z-linked orthologs and characterised their molecular differentiation. Our results suggest that the multiple sex-chromosome system evolved in a common ancestor as a result of dynamic genome reshuffling through repeated rearrangements between the sex chromosomes and autosomes, including translocations, fusions and fissions. Thus, the initial formation of neo-sex chromosomes could not have played a role in reproductive isolation between these Leptidea species. However, the subsequent species-specific fissions of several neo-sex chromosomes could have contributed to their reproductive isolation. Then, significantly increased numbers of Z-linked genes and independent neo-W chromosome degeneration could accelerate the accumulation of genetic incompatibilities between populations and promote their divergence resulting in speciation.

Activity and inactivity of moth sex chromosomes in somatic and meiotic cells

Moths and butterflies (Lepidoptera) are the most species-rich group of animals with female heterogamety, females mostly having a WZ, males a ZZ sex chromosome constitution. We studied chromatin conformation, activity, and inactivity of the sex chromosomes in the flour moth Ephestia kuehniella and the silkworm Bombyx mori, using immunostaining with anti-H3K9me2/3, anti-RNA polymerase II, and fluoro-uridine (FU) labelling of nascent transcripts, with conventional widefield fluorescence microscopy and 'spatial structured illumination microscopy' (3D-SIM). The Z chromosome is euchromatic in somatic cells and throughout meiosis. It is transcriptionally active in somatic cells and in the postpachaytene stage of meiosis. The W chromosome in contrast is heterochromatic in somatic cells as well as in meiotic cells at pachytene, but euchromatic and transcriptionally active like all other chromosomes at postpachytene. As the W chromosomes are apparently devoid of protein-coding genes, their transcripts must be non-coding. We found no indication of 'meiotic sex chromosome inactivation' (MSCI) in the two species.

Molecular sexing of Lepidoptera

We developed a universal method of Lepidoptera molecular sexing. The method is based on comparing the number of copies of the same gene in different sexes. Males of the majority of lepidopteran species have two Z chromosomes, whereas females have only one Z chromosome. Correspondingly, the number of copies of each gene located on this chromosome differs by two times between males and females. For quantitative estimation, we used qPCR. Via multiple alignment of the kettin (a Z chromosome gene) nucleotide sequences, we detected the most conserved fragment and designed primers with broad interspecies specificity for Lepidoptera. Using these primers, we successfully determined the sex of three lepidopteran species belonging to different superfamilies. The developed method is a simple, cost-effective and high-throughput technique for routine sexing. The sex of lepidopteran individuals can be examined at any developmental stage.

Genetic Diversity on the Sex Chromosomes

Levels and patterns of genetic diversity can provide insights into a population’s history. In species with sex chromosomes, differences between genomic regions with unique inheritance patterns can be used to distinguish between different sets of possible demographic and selective events. This review introduces the differences in population history for sex chromosomes and autosomes, provides the expectations for genetic diversity across the genome under different evolutionary scenarios, and gives an introductory description for how deviations in these expectations are calculated and can be interpreted. Predominantly, diversity on the sex chromosomes has been used to explore and address three research areas: 1) Mating patterns and sex-biased variance in reproductive success, 2) signatures of selection, and 3) evidence for modes of speciation and introgression. After introducing the theory, this review catalogs recent studies of genetic diversity on the sex chromosomes across species within the major research areas that sex chromosomes are typically applied to, arguing that there are broad similarities not only between male-heterogametic (XX/XY) and female-heterogametic (ZZ/ZW) sex determination systems but also any mating system with reduced recombination in a sex-determining region. Further, general patterns of reduced diversity in nonrecombining regions are shared across plants and animals. There are unique patterns across populations with vastly different patterns of mating and speciation, but these do not tend to cluster by taxa or sex determination system.

Silkworm genetic sexing through W chromosome-linked, targeted gene integration

Sex separation methods are critical for genetic sexing systems in commercial insect production and sterile insect techniques. Integration of selectable marker genes into a sex chromosome is particularly useful in insects with a heterogametic sex determination system. Here, we describe targeted gene integration of fluorescent marker expression cassettes into a randomly amplified polymorphic DNA (RAPD) marker region in the W chromosome of the lepidopteran model insect Bombyx mori using transcriptional activator-like effector nuclease (TALEN)-mediated genome editing. This silkworm strain shows ubiquitous female-specific red or green fluorescence from the embryonic to adult stages. Furthermore, we developed a binary, female-specific, embryonic lethality system combining the TALEN and the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) technology. This system includes one strain with TALEN-mediated, W-specific Cas9 expression driven by the silkworm germ cell-specific nanos (nos) promoter and another strain with U6-derived single-guide RNA (sgRNA) expression targeting transformer 2 (tra2), an essential gene for silkworm embryonic development. Filial 1 (F1) hybrids exhibit complete female-specific lethality during embryonic stages. Our study provides a promising approach for B. mori genetic sexing and sheds light on developing sterile insect techniques in other insect species, especially in lepidopteran pests with WZ/ZZ sex chromosome systems.

Genome-Wide Characterization of Endogenous Retroviruses in Bombyx mori Reveals the Relatives and Activity of env Genes

Endogenous retroviruses (ERVs) are retroviral sequences that remain fixed in the host genome, where they could play an important role. Some ERVs have been identified in insects and proven to have infectious properties. However, no information is available regarding Bombyx mori ERVs (BmERVs) to date. Here, we systematically identified 256 potential BmERVs in the silkworm genome via a whole-genome approach. BmERVs were relatively evenly distributed across each of the chromosomes and accounted for about 25% of the silkworm genome. All BmERVs were classified as young ERVs, with insertion times estimated to be less than 10 million years. Seven BmERVs possessing the env genes were identified. With the exception of the Orf133 Helicoverpa armigera nuclear polyhedrosis virus, the env sequences of BmERVs were distantly related to genes encoding F (Fa and Fb) and GP64 proteins from Group I and Group II NPVs. In addition, only the amino acid sequence of the BmERV-21 envelope protein shared a similar putative furin-like cleavage site and fusion peptide with Group II baculoviruses. All of the env genes in the seven BmERVs were verified to exist in the genome and be expressed in the midgut and fat bodies, which suggest that BmERVs might play an important role in the host biology.

Insights into the Structure of the Spruce Budworm (Choristoneura fumiferana) Genome, as Revealed by Molecular Cytogenetic Analyses and a High-Density Linkage Map

Genome structure characterization can contribute to a better understanding of processes such as adaptation, speciation, and karyotype evolution, and can provide useful information for refining genome assemblies. We studied the genome of an important North American boreal forest pest, the spruce budworm, Choristoneura fumiferana, through a combination of molecular cytogenetic analyses and construction of a high-density linkage map based on single nucleotide polymorphism (SNP) markers obtained through a genotyping-by-sequencing (GBS) approach. Cytogenetic analyses using fluorescence in situ hybridization methods confirmed the haploid chromosome number of n = 30 in both sexes of C. fumiferana and showed, for the first time, that this species has a WZ/ZZ sex chromosome system. Synteny analysis based on a comparison of the Bombyx mori genome and the C. fumiferana linkage map revealed the presence of a neo-Z chromosome in the latter species, as previously reported for other tortricid moths. In this neo-Z chromosome, we detected an ABC transporter C2 (ABCC2) gene that has been associated with insecticide resistance. Sex-linkage of the ABCC2 gene provides a genomic context favorable to selection and rapid spread of resistance against Bacillus thuringiensis serotype kurstaki (Btk), the main insecticide used in Canada to control spruce budworm populations. Ultimately, the linkage map we developed, which comprises 3586 SNP markers distributed over 30 linkage groups for a total length of 1720.41 cM, will be a valuable tool for refining our draft assembly of the spruce budworm genome.

Sex Chromosomes of the Iconic Moth Abraxas grossulariata (Lepidoptera, Geometridae) and Its Congener A. sylvata

The magpie moth, Abraxas grossulariata, is an iconic species in which female heterogamety was discovered at the beginning of the 20th century. However, the sex chromosomes of this species have not yet been cytologically identified. We describe the sex chromosomes of A. grossulariata and its congener, A. sylvata. Although these species split only around 9.5 million years ago, and both species have the expected WZ/ZZ chromosomal system of sex determination and their sex chromosomes share the major ribosomal DNA (rDNA) representing the nucleolar organizer region (NOR), we found major differences between their karyotypes, including between their sex chromosomes. The species differ in chromosome number, which is 2n = 56 in A. grossularita and 2n = 58 in A. sylvata. In addition, A. grossularita autosomes exhibit massive autosomal blocks of heterochromatin, which is a very rare phenomenon in Lepidoptera, whereas the autosomes of A. sylvata are completely devoid of distinct heterochromatin. Their W chromosomes differ greatly. Although they are largely composed of female-specific DNA sequences, as shown by comparative genomic hybridization, cross-species W-chromosome painting revealed considerable sequence differences between them. The results suggest a relatively rapid molecular divergence of Abraxas W chromosomes by the independent spreading of female-specific repetitive sequences.

The genome of the Hi5 germ cell line from Trichoplusia ni, an agricultural pest and novel model for small RNA biology

We report a draft assembly of the genome of Hi5 cells from the lepidopteran insect pest, Trichoplusia ni, assigning 90.6% of bases to one of 28 chromosomes and predicting 14,037 protein-coding genes. Chemoreception and detoxification gene families reveal T. ni-specific gene expansions that may explain its widespread distribution and rapid adaptation to insecticides. Transcriptome and small RNA data from thorax, ovary, testis, and the germline-derived Hi5 cell line show distinct expression profiles for 295 microRNA- and >393 piRNA-producing loci, as well as 39 genes encoding small RNA pathway proteins. Nearly all of the W chromosome is devoted to piRNA production, and T. ni siRNAs are not 2´-O-methylated. To enable use of Hi5 cells as a model system, we have established genome editing and single-cell cloning protocols. The T. ni genome provides insights into pest control and allows Hi5 cells to become a new tool for studying small RNAs ex vivo.

A common moth called the cabbage looper is becoming increasingly relevant to the scientific community. Its caterpillars are a serious threat to cabbage, broccoli and cauliflower crops, and they have started to resist the pesticides normally used to control them. Moreover, the insect’s germline cells – the ones that will produce sperm and eggs – are used in laboratories as ‘factories’ to artificially produce proteins of interest.

The germline cells also host a group of genetic mechanisms called RNA silencing. One of these processes is known as piRNA, and it protects the genome against ‘jumping genes’. These genetic elements can cause mutations by moving from place to place in the DNA: in germline cells, piRNA suppresses them before the genetic information is transmitted to the next generation. Not all germline cells grow equally well under experimental conditions, or are easy to use to examine piRNA mechanisms in a laboratory. The germline cells from the cabbage looper, on the other hand, have certain characteristics that would make them ideal to study piRNA in insects.

However, the genome of the moth had not yet been fully resolved. This hinders research on new ways of controlling the pest, on how to use the germline cells to produce more useful proteins, or on piRNA.

Decoding a genome requires several steps. First, the entire genetic information is broken in short sections that can then be deciphered. Next, these segments need to be ‘assembled’ – put together, and in the right order, to reconstitute the entire genome. Certain portions of the genome, which are formed of repeats of the same sections, can be difficult to assemble. Finally, the genome must be annotated: the different regions – such as the genes – need to be identified and labeled.

Here, Fu et al. assembled and annotated the genome of the cabbage looper, and in the process developed strategies that could be used for other species with a lot of repeated sequences in their genomes. Having access to the looper’s full genetic information makes it possible to use their germline cells to produce new types of proteins, for example for pharmaceutical purposes. Fu et al. went on to make working with these cells even easier by refining protocols so that modern research techniques, such as the gene-editing technology CRISPR-Cas9, can be used on the looper germline cells.

The mapping of the genome also revealed that the genes involved in removing toxins from the insects’ bodies are rapidly evolving, which may explain why the moths readily become resistant to insecticides. This knowledge could help finding new ways of controlling the pest.

Finally, the genes involved in RNA silencing were labeled: results show that an entire chromosome is the source of piRNAs. Combined with the new protocols developed by Fu et al., this could make cabbage looper germline cells the default option for any research into the piRNA mechanism. How piRNA works in the moth could inform work on human piRNA, as these processes are highly similar across the animal kingdom.

W-enriched satellite sequence in the Indian meal moth, Plodia interpunctella (Lepidoptera, Pyralidae)

The W chromosome of most lepidopteran species represents the largest heterochromatin entity in the female genome. Although satellite DNA is a typical component of constitutive heterochromatin, there are only a few known satellite DNAs (satDNAs) located on the W chromosome in moths and butterflies. In this study, we isolated and characterized new satDNA (PiSAT1) from microdissected W chromosomes of the Indian meal moth, Plodia interpunctella. Even though the PiSAT1 is mainly localized near the female-specific segment of the W chromosome, short arrays of this satDNA also occur on autosomes and/or the Z chromosome. Probably due to the predominant location in the non-recombining part of the genome, PiSAT1 exhibits a relatively large nucleotide variability in its monomers. However, at least a part of all predicted functional motifs is located in conserved regions. Moreover, we detected polyadenylated transcripts of PiSAT1 in all developmental stages and in both sexes (female and male larvae, pupae and adults). Our results suggest a potential structural and functional role of PiSAT1 in the P. interpunctella genome, which is consistent with accumulating evidence for the important role of satDNAs in eukaryotic genomes.

CpSAT-1, a transcribed satellite sequence from the codling moth, Cydia pomonella

Satellite DNA (satDNA) is a non-coding component of eukaryotic genomes, located mainly in heterochromatic regions. Relevance of satDNA began to emerge with accumulating evidence of its potential yet hardly comprehensible role that it can play in the genome of many organisms. We isolated the first satDNA of the codling moth (Cydia pomonella, Tortricidae, Lepidoptera), a species with holokinetic chromosomes and a single large heterochromatic element, the W chromosome in females. The satDNA, called CpSAT-1, is located on all chromosomes of the complement, although in different amounts. Surprisingly, the satellite is almost missing in the heterochromatic W chromosome. Additionally, we isolated mRNA from all developmental stages (1st-5th instar larva, pupa, adult), both sexes (adult male and female) and several tissues (Malpighian tubules, gut, heart, testes, and ovaries) of the codling moth and showed the CpSAT-1 sequence was transcribed in all tested samples. Using CpSAT-1 specific primers we amplified, cloned and sequenced 40 monomers from cDNA and gDNA, respectively. The sequence analysis revealed a high mutation rate and the presence of potentially functional motifs, mainly in non-conserved regions of the monomers. Both the chromosomal distribution and the sequence analysis suggest that CPSAT-1 has no function in the C. pomonella genome.

Guardian small RNAs and sex determination

The W chromosome of the silkworm Bombyx mori has been known to determine femaleness for more than 80 years. However, the feminizing gene has not been molecularly identified, because the B. mori W chromosome is almost fully occupied by a large number of transposable elements. The W chromosome-derived feminizing factor of B. mori was recently shown to be a female-specific PIWI-interacting RNA (piRNA). piRNAs are small RNAs that potentially repress invading “non-self” elements (e.g., transposons and virus-like elements) by associating with PIWI proteins. Our results revealed that female-specific piRNA precursors, which we named Fem, are transcribed from the sex-determining region of the W chromosome at the early embryonic stage and are processed into a single mature piRNA (Fem piRNA). Fem piRNA forms a complex with Siwi (silkworm Piwi), which cleaves a protein-coding mRNA transcribed from the Z chromosome. RNA interference of this Z-linked gene, which we named Masc, revealed that this gene encodes a protein required for masculinization and dosage compensation. Fem and Masc both participate in the ping-pong cycle of the piRNA amplification loop by associating with the 2 B. mori PIWI proteins Siwi and BmAgo3 (silkworm Ago3), respectively, indicating that the piRNA-mediated interaction between the 2 sex chromosomes is the primary signal for the B. mori sex determination cascade. Fem is a non-transposable repetitive sequence on the W chromosome, whereas Masc is a single-copy protein-coding gene. It is of great interest how the piRNA system recognizes “self ”Masc mRNA as “non-self” RNA.

Two Conserved Cysteine Residues Are Required for the Masculinizing Activity of the Silkworm Masc Protein

We have recently discovered that the Masculinizer (Masc) gene encodes a CCCH tandem zinc finger protein, which controls both masculinization and dosage compensation in the silkworm Bombyx mori. In this study, we attempted to identify functional regions or residues that are required for the masculinizing activity of the Masc protein. We constructed a series of plasmids that expressed the Masc derivatives and transfected them into a B. mori ovary-derived cell line, BmN-4. To assess the masculinizing activity of the Masc derivatives, we investigated the splicing patterns of B. mori doublesex (Bmdsx) and the expression levels of B. mori IGF-II mRNA-binding protein, a splicing regulator of Bmdsx, in Masc cDNA-transfected BmN-4 cells. We found that two zinc finger domains are not required for the masculinizing activity. We also identified that the C-terminal 288 amino acid residues are sufficient for the masculinizing activity of the Masc protein. Further detailed analyses revealed that two cysteine residues, Cys-301 and Cys-304, in the highly conserved region among lepidopteran Masc proteins are essential for the masculinizing activity in BmN-4 cells. Finally, we showed that Masc is a nuclear protein, but its nuclear localization is not tightly associated with the masculinizing activity.

Positional cloning of the sex-linked giant egg (Ge) locus in the silkworm, Bombyx mori

The giant egg (Ge) locus is a Z-linked mutation that leads to the production of large eggs. Cytological observations suggest that an unusual translocation of a large fragment of the W chromosome bearing a putative egg size-determining gene, Esd, gave rise to giant egg mutants. However, there is currently no molecular evidence confirming either a W-Z translocation or the presence of Esd on the W chromosome. To elucidate the origin of giant egg mutants, we performed positional cloning. We observed that the Bombyx mori. orthologue of the human Phytanoyl-CoA dioxygenase domain containing 1 gene (PHYHD1) is disrupted in giant egg mutants. PHYHD1 is highly conserved in eukaryotes and is predicted to be a Fe(II) and 2-oxoglutarate-dependent oxygenase. Exon skipping in one of the two available Ge mutants is probably caused by the insertion of a non-long terminal repeat transposon into intron 4 in the vicinity of the 5' splice site. Segmental duplication in Ge(2) , an independent allele, was caused by unequal recombination between short interspersed elements inserted into introns 3 and 5. Our results indicate that (1) Bombyx PHYHD1 is responsible for the Ge mutants and that (2) the Ge locus is unrelated to the W-linked putative Esd. To our knowledge, this is the first report describing the phenotypic defects caused by mutations in PHYHD1 orthologues.

The Glanville fritillary genome retains an ancient karyotype and reveals selective chromosomal fusions in Lepidoptera

Previous studies have reported that chromosome synteny in Lepidoptera has been well conserved, yet the number of haploid chromosomes varies widely from 5 to 223. Here we report the genome (393 Mb) of the Glanville fritillary butterfly (Melitaea cinxia; Nymphalidae), a widely recognized model species in metapopulation biology and eco-evolutionary research, which has the putative ancestral karyotype of n=31. Using a phylogenetic analyses of Nymphalidae and of other Lepidoptera, combined with orthologue-level comparisons of chromosomes, we conclude that the ancestral lepidopteran karyotype has been n=31 for at least 140 My. We show that fusion chromosomes have retained the ancestral chromosome segments and very few rearrangements have occurred across the fusion sites. The same, shortest ancestral chromosomes have independently participated in fusion events in species with smaller karyotypes. The short chromosomes have higher rearrangement rate than long ones. These characteristics highlight distinctive features of the evolutionary dynamics of butterflies and moths.

Butterflies and moths (Lepidoptera) vary in chromosome number. Here, the authors sequence the genome of the Glanville fritillary butterfly, Melitaea cinxia, show it has the ancestral lepidopteran karyotype and provide insight into how chromosomal fusions have shaped karyotype evolution in butterflies and moths.

Flipping the doublesex switch with a piRNA

Recent work in the silkworm Bombyx mori has uncovered a novel Piwi-interacting RNA regulator of the sex determination switch doublesex.

A single female-specific piRNA is the primary determiner of sex in the silkworm

The silkworm Bombyx mori uses a WZ sex determination system that is analogous to the one found in birds and some reptiles. In this system, males have two Z sex chromosomes, whereas females have Z and W sex chromosomes. The silkworm W chromosome has a dominant role in female determination, suggesting the existence of a dominant feminizing gene in this chromosome. However, the W chromosome is almost fully occupied by transposable element sequences, and no functional protein-coding gene has been identified so far. Female-enriched PIWI-interacting RNAs (piRNAs) are the only known transcripts that are produced from the sex-determining region of the W chromosome, but the function(s) of these piRNAs are unknown. Here we show that a W-chromosome-derived, female-specific piRNA is the feminizing factor of B. mori. This piRNA is produced from a piRNA precursor which we named Fem. Fem sequences were arranged in tandem in the sex-determining region of the W chromosome. Inhibition of Fem-derived piRNA-mediated signalling in female embryos led to the production of the male-specific splice variants of B. mori doublesex (Bmdsx), a gene which acts at the downstream end of the sex differentiation cascade. A target gene of Fem-derived piRNA was identified on the Z chromosome of B. mori. This gene, which we named Masc, encoded a CCCH-type zinc finger protein. We show that the silencing of Masc messenger RNA by Fem piRNA is required for the production of female-specific isoforms of Bmdsx in female embryos, and that Masc protein controls both dosage compensation and masculinization in male embryos. Our study characterizes a single small RNA that is responsible for primary sex determination in the WZ sex determination system.

Anatomical and functional analysis of domestication effects on the olfactory system of the silkmoth Bombyx mori

The silkmoth Bombyx mori is the main producer of silk worldwide and has furthermore become a model organism in biological research, especially concerning chemical communication. However, the impact domestication might have had on the silkmoth's olfactory sense has not yet been investigated. Here, we show that the pheromone detection system in B. mori males when compared with their wild ancestors Bombyx mandarina seems to have been preserved, while the perception of environmental odorants in both sexes of domesticated silkmoths has been degraded. In females, this physiological impairment was mirrored by a clear reduction in olfactory sensillum numbers. Neurophysiological experiments with hybrids between wild and domesticated silkmoths suggest that the female W sex chromosome, so far known to have the sole function of determining femaleness, might be involved in the detection of environmental odorants. Moreover, the coding of odorants in the brain, which is usually similar among closely related moths, differs strikingly between B. mori and B. mandarina females. These results indicate that domestication has had a strong impact on odour detection and processing in the olfactory model species B. mori.

Large scale full-length cDNA sequencing reveals a unique genomic landscape in a lepidopteran model insect, Bombyx mori

The establishment of a complete genomic sequence of silkworm, the model species of Lepidoptera, laid a foundation for its functional genomics. A more complete annotation of the genome will benefit functional and comparative studies and accelerate extensive industrial applications for this insect. To realize these goals, we embarked upon a large-scale full-length cDNA collection from 21 full-length cDNA libraries derived from 14 tissues of the domesticated silkworm and performed full sequencing by primer walking for 11,104 full-length cDNAs. The large average intron size was 1904 bp, resulting from a high accumulation of transposons. Using gene models predicted by GLEAN and published mRNAs, we identified 16,823 gene loci on the silkworm genome assembly. Orthology analysis of 153 species, including 11 insects, revealed that among three Lepidoptera including Monarch and Heliconius butterflies, the 403 largest silkworm-specific genes were composed mainly of protective immunity, hormone-related, and characteristic structural proteins. Analysis of testis-/ovary-specific genes revealed distinctive features of sexual dimorphism, including depletion of ovary-specific genes on the Z chromosome in contrast to an enrichment of testis-specific genes. More than 40% of genes expressed in specific tissues mapped in tissue-specific chromosomal clusters. The newly obtained FL-cDNA sequences enabled us to annotate the genome of this lepidopteran model insect more accurately, enhancing genomic and functional studies of Lepidoptera and comparative analyses with other insect orders, and yielding new insights into the evolution and organization of lepidopteran-specific genes.

High-throughput sequencing of a single chromosome: a moth W chromosome

Y and W chromosomes have mostly been excluded from whole genome sequencing projects. Due to the high amount of repetitive sequences they are 'difficult' to assemble and therefore need special treatment in the form of, e.g. adapted assembly programs, a range of different libraries, and accurate maps, if possible. A minimum requirement for these approaches is pure template DNA. We therefore microdissected the W chromatin of highly polyploid cells from the flour moth, Ephestia kuehniella, and used Roche/454 and Sanger sequencing to generate 72.6 Mbp of DNA sequence. Nominal coverage was 4.3× of the 16.7 Mbp of W chromosomal DNA. We used these data to assess the genetic content of the W chromosome. This approach allowed us to determine constituent families of transposable elements, microsatellites, and recent insertion sites of mitochondrial DNA. However, no conventional protein-coding gene has yet been found. The sequence collection is a rich source for the definition of W-specific PCR markers and the reconstruction of W chromosome loci, as a step towards full reconstruction of the chromosome.

Recent gene-capture on the UV sex chromosomes of the moss Ceratodon purpureus

Sex chromosomes evolve from ordinary autosomes through the expansion and subsequent degeneration of a region of suppressed recombination that is inherited through one sex. Here we investigate the relative timing of these processes in the UV sex chromosomes of the moss Ceratodon purpureus using molecular population genetic analyses of eight newly discovered sex-linked loci. In this system, recombination is suppressed on both the female-transmitted (U) sex chromosome and the male-transmitted (V) chromosome. Genes on both chromosomes therefore should show the deleterious effects of suppressed recombination and sex-limited transmission, while purifying selection should maintain homologs of genes essential for both sexes on both sex chromosomes. Based on analyses of eight sex-linked loci, we show that the nonrecombining portions of the U and V chromosomes expanded in at least two events (~0.6-1.3 MYA and ~2.8-3.5 MYA), after the divergence of C. purpureus from its dioecious sister species, Trichodon cylindricus and Cheilothela chloropus. Both U- and V-linked copies showed reduced nucleotide diversity and limited population structure, compared to autosomal loci, suggesting that the sex chromosomes experienced more recent selective sweeps that the autosomes. Collectively these results highlight the dynamic nature of gene composition and molecular evolution on nonrecombining portions of the U and V sex chromosomes.

Rapid turnover of the W chromosome in geographical populations of wild silkmoths, Samia cynthia ssp

Our previous studies revealed a considerably high level of chromosomal polymorphism in wild silkmoths, Samia cynthia ssp. (Lepidoptera: Saturniidae). Geographical populations of this species complex differ in chromosome numbers and show derived sex chromosome systems including Z0/ZZ in S. cynthia ricini (2n = 27/28; Vietnam), neo-Wneo-Z/neo-Zneo-Z in S. cynthia walkeri (2n = 26/26; Sapporo, Hokkaido) and neo-WZ1Z2/Z1Z1Z2Z2 in S. cynthia subsp. indet. (2n = 25/26; Nagano, Honshu). In this study, we collected specimens of S. cynthia pryeri in Japanese islands Kyushu, Shikoku and Honshu, with an ancestral-like karyotype of 2n = 28 in both sexes and a WZ/ZZ sex chromosome system, except for one population, in which females have lost the W chromosome. However, the S. cynthia pryeri W chromosome showed a very unusual morphology: It was composed of a highly heterochromatic body, which remained condensed throughout the whole cell cycle and of a euchromatin-like "tail." We examined molecular composition of the W and neo-W chromosomes in S. cynthia subspecies by comparative genomic hybridisation and fluorescence in situ hybridisation with W chromosome painting probes prepared from laser-microdissected W chromatin of S. cynthia pryeri. These methods revealed that the molecular composition of highly heterochromatic part of the S. cynthia pryeri W chromosome is very different and lacks homology in the genomes of other subspecies, whereas the euchromatin-like part of the W chromosome corresponds to a heterochromatic part of the neo-W chromosomes in S. cynthia walkeri and S. cynthia subsp. indet. Our findings suggest that the curious WZ system of S. cynthia pryeri may represent an ancestral state of the Samia species complex but do not exclude an alternative hypothesis of its derived origin.

Linkage map of the peppered moth, Biston betularia (Lepidoptera, Geometridae): a model of industrial melanism

We have constructed a linkage map for the peppered moth (Biston betularia), the classical ecological genetics model of industrial melanism, aimed both at localizing the network of loci controlling melanism and making inferences about chromosome dynamics. The linkage map, which is based primarily on amplified fragment length polymorphisms (AFLPs) and genes, consists of 31 linkage groups (LGs; consistent with the karyotype). Comparison with the evolutionarily distant Bombyx mori suggests that the gene content of chromosomes is highly conserved. Gene order is conserved on the autosomes, but noticeably less so on the Z chromosome, as confirmed by physical mapping using bacterial artificial chromosome fluorescence in situ hybridization (BAC-FISH). Synteny mapping identified three pairs of B. betularia LGs (11/29, 23/30 and 24/31) as being orthologous to three B. mori chromosomes (11, 23 and 24, respectively). A similar finding in an outgroup moth (Plutella xylostella) indicates that the B. mori karyotype (n=28) is a phylogenetically derived state resulting from three chromosome fusions. As with other Lepidoptera, the B. betularia W chromosome consists largely of repetitive sequence, but exceptionally we found a W homolog of a Z-linked gene (laminin A), possibly resulting from ectopic recombination between the sex chromosomes. The B. betularia linkage map, featuring the network of known melanization genes, serves as a resource for melanism research in Lepidoptera. Moreover, its close resemblance to the ancestral lepidopteran karyotype (n=31) makes it a useful reference point for reconstructing chromosome dynamic events and ancestral genome architectures. Our study highlights the unusual evolutionary stability of lepidopteran autosomes; in contrast, higher rates of intrachromosomal rearrangements support a special role of the Z chromosome in adaptive evolution and speciation.

Altered expression of testis-specific genes, piRNAs, and transposons in the silkworm ovary masculinized by a W chromosome mutation

BACKGROUND

In the silkworm, Bombyx mori, femaleness is strongly controlled by the female-specific W chromosome. Originally, it was presumed that the W chromosome encodes female-determining gene(s), accordingly called Fem. However, to date, neither Fem nor any protein-coding gene has been identified from the W chromosome. Instead, the W chromosome is occupied with numerous transposon-related sequences. Interestingly, the silkworm W chromosome is a source of female-enriched PIWI-interacting RNAs (piRNAs). piRNAs are small RNAs of 23-30 nucleotides in length, which are required for controlling transposon activity in animal gonads. A recent study has identified a novel mutant silkworm line called KG, whose mutation in the W chromosome causes severe female masculinization. However, the molecular nature of KG line has not been well characterized yet.

RESULTS

Here we molecularly characterize the KG line. Genomic PCR analyses using currently available W chromosome-specific PCR markers indicated that no large deletion existed in the KG W chromosome. Genetic analyses demonstrated that sib-crosses within the KG line suppressed masculinization. Masculinization reactivated when crossing KG females with wild type males. Importantly, the KG ovaries exhibited a significantly abnormal transcriptome. First, the KG ovaries misexpressed testis-specific genes. Second, a set of female-enriched piRNAs was downregulated in the KG ovaries. Third, several transposons were overexpressed in the KG ovaries.

CONCLUSIONS

Collectively, the mutation in the KG W chromosome causes broadly altered expression of testis-specific genes, piRNAs, and transposons. To our knowledge, this is the first study that describes a W chromosome mutant with such an intriguing phenotype.

The silkworm W chromosome is a source of female-enriched piRNAs

In the silkworm, Bombyx mori, the W chromosome plays a dominant role in female determination. However, neither protein-coding genes nor transcripts have so far been isolated from the W chromosome. Instead, a large amount of functional transposable elements and their remnants are accumulated on the W chromosome. PIWI-interacting RNAs (piRNAs) are 23-30-nt-long small RNAs that potentially act as sequence-specific guides for PIWI proteins to silence transposon activity in animal gonads. In this study, by comparing ovary- and testis-derived piRNAs, we identified numerous female-enriched piRNAs. Our data indicated that female-enriched piRNAs are derived from the W chromosome. Moreover, comparative analyses on piRNA profiles from a series of W chromosome mutant strains revealed a striking enrichment of a specific set of transposon-derived piRNAs in the putative sex-determining region. Collectively, we revealed the nature of the silkworm W chromosome as a source of piRNAs.

Identification and characterization of the fusion transcript, composed of the apterous homolog and a putative protein phosphatase gene, generated by 1.5-Mb interstitial deletion in the vestigial (Vg) mutant of Bombyx mori

The vestigial (Vg) mutant is a Z-linked mutant that causes vestigial wings in the silkworm, Bombyx mori. We have previously reported a 1.5-Mb interstitial deletion on the Z chromosome bearing the Vg mutation (Z(Vg) chromosome). In this study, we found that exons 3-8 of a gene named Bmptp-Z encoding a putative tyrosine-specific protein phosphatase are deleted by the 1.5-Mb interstitial deletion. We found that a gene encoding the Bombyx homolog of Drosophila Apterous (BmAp-A) protein is located 4.5 kb downstream of the distal breakpoint of the 1.5-Mb interstitial deletion. Moreover, an in-frame fusion transcript composed of the 5' part of Bmptp-Z and the 3' part of Bmap-A is generated specific to the Z(Vg) chromosome. Effects of the in-frame fusion transcript on the vestigial phenotype are discussed.

Evolution of sex chromosomes in insects

Sex chromosomes have many unusual features relative to autosomes. Y (or W) chromosomes lack genetic recombination, are male- (female-) limited, and show an abundance of genetically inert heterochromatic DNA but contain few functional genes. X (or Z) chromosomes also show sex-biased transmission (i.e., X chromosomes show female-biased and Z-chromosomes show male-biased inheritance) and are hemizygous in the heterogametic sex. Their unusual ploidy level and pattern of inheritance imply that sex chromosomes play a unique role in many biological processes and phenomena, including sex determination, epigenetic chromosome-wide regulation of gene expression, the distribution of genes in the genome, genomic conflict, local adaptation, and speciation. The vast diversity of sex chromosome systems in insects--ranging from the classical male heterogametic XY system in Drosophila to ZW systems in Lepidoptera or mobile genes determining sex as found in house flies--implies that insects can serve as unique model systems to study various functional and evolutionary aspects of these different processes.

Repression of tyrosine hydroxylase is responsible for the sex-linked chocolate mutation of the silkworm, Bombyx mori

Pigmentation patterning has long interested biologists, integrating topics in ecology, development, genetics, and physiology. Wild-type neonatal larvae of the silkworm, Bombyx mori, are completely black. By contrast, the epidermis and head of larvae of the homozygous recessive sex-linked chocolate (sch) mutant are reddish brown. When incubated at 30 degrees C, mutants with the sch allele fail to hatch; moreover, homozygous mutants carrying the allele sch lethal (sch(l)) do not hatch even at room temperature (25 degrees C). By positional cloning, we narrowed a region containing sch to 239,622 bp on chromosome 1 using 4,501 backcross (BC1) individuals. Based on expression analyses, the best sch candidate gene was shown to be tyrosine hydroxylase (BmTh). BmTh coding sequences were identical among sch, sch(l), and wild-type. However, in sch the approximately 70-kb sequence was replaced with approximately 4.6 kb of a Tc1-mariner type transposon located approximately 6 kb upstream of BmTh, and in sch(l), a large fragment of an L1Bm retrotransposon was inserted just in front of the transcription start site of BmTh. In both cases, we observed a drastic reduction of BmTh expression. Use of RNAi with BmTh prevented pigmentation and hatching, and feeding of a tyrosine hydroxylase inhibitor also suppressed larval pigmentation in the wild-type strain, pnd(+) and in a pS (black-striped) heterozygote. Feeding L-dopa to sch neonate larvae rescued the mutant phenotype from chocolate to black. Our results indicate the BmTh gene is responsible for the sch mutation, which plays an important role in melanin synthesis producing neonatal larval color.

Identification of 2chromosome region translocated onto the W chromosome by RFLP with EST-cDNA clones in the Gensei-kouken strains of the mulberry silkworm, Bombyx mori L

In silkworms, sex-limited strains are either obtained spontaneously or induced by X-rays or gamma rays. When a fragment of an autosome carrying a dominant allele of those genes responsible for certain characters is translocated onto a W chromosome, the female of the successive generations will express these phenotypic characters and sex discrimination can be facilitated. Gensei-kouken strains are sex-limited strains of silkworms developed by irradiating the pupae with gamma rays, by which a portion of the second chromosome is translocated onto the W chromosome. In these improved strains, the females are yellow-blooded and spin yellow cocoons. By using the EST-cDNA clones mapped on the Z chromosome, we identified the sex according to the polymorphic banding pattern or intensity of the signals. Furthermore, by using the clones on the second chromosome, the region of the second chromosome translocated onto the W chromosome was also defined. In both the A95 and A 96 strains selected for the present study, only the mid-portion of the second chromosome was translocated. The differences in length of the fragments translocated in these strains are discussed.

Dosage analysis of Z chromosome genes using microarray in silkworm, Bombyx mori

In many organisms, dosage compensation is needed to equalize sex-chromosome gene expression in males and females. Several genes on silkworm Z chromosome were previously detected to show a higher expression level in males and lacked dosage compensation. Whether silkworm lacks global dosage compensation still remains poorly known. Here, we analyzed male:female (M:F) ratios of expression of chromosome-wide Z-linked genes in the silkworm using microarray data. The expression levels of genes on Z chromosome in each tissue were significantly higher in males compared to females, which indicates no global dosage compensation in silkworm. Interestingly, we also found some genes with no bias (M:F ratio: 0.8-1.2) on the Z chromosome. Comparison of male-biased (M:F ratio more than 1.5) and unbiased genes indicated that the two sets of the genes have functional differences. Analysis of gene expression by sex showed that M:F ratios were, to some extent, associated with their expression levels. These results provide useful clues to further understanding roles of dosage of Z chromosome and some Z-linked sexual differences in silkworms.

Cytogenetic characterization and AFLP-based genetic linkage mapping for the butterfly Bicyclus anynana, covering all 28 karyotyped chromosomes

Background

The chromosome characteristics of the butterfly Bicyclus anynana, have received little attention, despite the scientific importance of this species. This study presents the characterization of chromosomes in this species by means of cytogenetic analysis and linkage mapping.

Methodology/Principal Findings

Physical genomic features in the butterfly B. anynana were examined by karyotype analysis and construction of a linkage map. Lepidoptera possess a female heterogametic W-Z sex chromosome system. The WZ-bivalent in pachytene oocytes of B. anynana consists of an abnormally small, heterochromatic W-chromosome with the Z-chromosome wrapped around it. Accordingly, the W-body in interphase nuclei is much smaller than usual in Lepidoptera. This suggests an intermediate stage in the process of secondary loss of the W-chromosome to a ZZ/Z sex determination system. Two nucleoli are present in the pachytene stage associated with an autosome and the WZ-bivalent respectively. Chromosome counts confirmed a haploid number of n = 28. Linkage mapping had to take account of absence of crossing-over in females, and of our use of a full-sib crossing design. We developed a new method to determine and exclude the non-recombinant uninformative female inherited component in offspring. The linkage map was constructed using a novel approach that uses exclusively JOINMAP-software for Lepidoptera linkage mapping. This approach simplifies the mapping procedure, avoids over-estimation of mapping distance and increases the reliability of relative marker positions. A total of 347 AFLP markers, 9 microsatellites and one single-copy nuclear gene covered all 28 chromosomes, with a mapping distance of 1354 cM. Conserved synteny of Tpi on the Z-chromosome in Lepidoptera was confirmed for B. anynana. The results are discussed in relation to other mapping studies in Lepidoptera.

Conclusions/Significance

This study adds to the knowledge of chromosome structure and evolution of an intensively studied organism. On a broader scale it provides an insight in Lepidoptera sex chromosome evolution and it proposes a simpler and more reliable method of linkage mapping than used for Lepidoptera to date.

Recent spread of a retrotransposon in the Silene latifolia genome, apart from the Y chromosome

Transposable elements often accumulate in nonrecombining regions, such as Y chromosomes. Contrary to this trend, a new Silene retrotransposon described here, has spread recently all over the genome of plant Silene latifolia, except its Y chromosome. This coincided with the latest steps of sex chromosome evolution in this species.

Silkworm nucleotide databases--current trends and future prospects

The domesticated silkworm, Bombyx mori serves as an ideal representative of lepidopteran species for a variety of scientific studies. As a result, databases have been created to organize information pertaining to the silkworm genome that is subject to constant updating. Of these, four main databases are important for store nucleotide information in the form of genomic data, ESTs and microsatelites. These databases also store data related to other lepidoptera and important insects, which help in insect biological research. Though a considerable amount of nucleotide data is currently available, there is a paucity of data related to silkworm and other lepidopteran proteins. Hence, the focus of this article is to present the current status of nucleotide databases of silkworm, avenues for improvement and possibilities for databases that could be created in the future.