Isolation of BAC clones containing conserved genes from libraries of three distantly related moths: a useful resource for comparative genomics of Lepidoptera

The Role of Repetitive Sequences in Repatterning of Major Ribosomal DNA Clusters in Lepidoptera

Genes for major ribosomal RNAs (rDNA) are present in multiple copies mainly organized in tandem arrays. The number and position of rDNA loci can change dynamically and their repatterning is presumably driven by other repetitive sequences. We explored a peculiar rDNA organization in several representatives of Lepidoptera with either extremely large or numerous rDNA clusters. We combined molecular cytogenetics with analyses of second- and third-generation sequencing data to show that rDNA spreads as a transcription unit and reveal association between rDNA and various repeats. Furthermore, we performed comparative long read analyses among the species with derived rDNA distribution and moths with a single rDNA locus, which is considered ancestral. Our results suggest that satellite arrays, rather than mobile elements, facilitate homology-mediated spread of rDNA via either integration of extrachromosomal rDNA circles or ectopic recombination. The latter arguably better explains preferential spread of rDNA into terminal regions of lepidopteran chromosomes as efficiency of ectopic recombination depends on the proximity of homologous sequences to telomeres.

Conservation and lineage-specific rearrangements in the GOBP/PBP gene complex of distantly related ditrysian Lepidoptera

General odorant binding proteins (GOBPs) and pheromone binding proteins (PBPs) form a monophyletic subfamily of insect odorant binding proteins (OBPs) specific for Lepidoptera, butterflies and moths. The GOBP/PBP genes include six subgroups (GOBP1–2, PBP-A–D) previously reported to form a complex arrayed in a conserved order in representative moths (superfamily Bombycoidea) and butterflies (Nymphalidae). Although our knowledge of lepidopteran genomes has increased greatly recently, the structure of the GOBP/PBP complex has been studied only for species that represent limited lineages of the highly diverged Ditrysia. To understand the evolution of this functionally important gene complex, we determined 69–149 kb genomic sequences that include GOBP2 and five PBP genes in three Ostrinia moths (Pyraloidea), O. nubilalis, O. furnacalis, and O. latipennis, using bacterial artificial chromosome (BAC) and fosmid clones. The structure of the GOBP2/PBP gene cluster was well conserved despite the different sex pheromone composition utilized by the three moths. Five expressed PBP genes in Ostrinia moths were the result of two duplications of PBP-A genes. Surprisingly, an allele containing a fusion gene between tandemly arrayed PBP-A genes was observed in O. nubilalis. We also revealed duplication and intra-chromosomal translocation of the GOBP1 gene in P. xylostella by fluorescence in situ hybridization (FISH) analysis. Additionally, we compared the structure of the GOBP/PBP gene complex of seventeen species covering six superfamilies and twelve families of the lepidopteran clade, Ditrysia, and found the gene order was basically conserved despite the frequent occurrence of lineage-specific gains, losses, inversions and translocations of these genes, compared with their neighboring genes. Our findings support the hypothesis that the structure of the GOBP/PBP gene complex was already established in the common ancestor of Ditrysia.

A FISH-based chromosome map for the European corn borer yields insights into ancient chromosomal fusions in the silkworm

A significant feature of the genomes of Lepidoptera, butterflies and moths, is the high conservation of chromosome organization. Recent remarkable progress in genome sequencing of Lepidoptera has revealed that syntenic gene order is extensively conserved across phylogenetically distant species. The ancestral karyotype of Lepidoptera is thought to be n=31; however, that of the most well-studied moth, Bombyx mori, is n=28, and diverse studies suggest that three chromosomal fusion events occurred in this lineage. To identify the boundaries between predicted ancient fusions involving B. mori chromosomes 11, 23 and 24, we constructed fluorescence in situ hybridization (FISH)-based chromosome maps of the European corn borer, Ostrinia nubilalis (n=31). We first determined a 511 Mb genomic sequence of the Asian corn borer, O. furnacalis, a congener of O. nubilalis, and isolated bacterial artificial chromosomes and fosmid clones that were expected to localize in candidate regions for the boundaries using these sequences. Combined with FISH and genetic analysis, we narrowed down the candidate regions to 40 kb-1.5 Mb, in strong agreement with a previous estimate based on the genome of a butterfly, Melitaea cinxia. The significant difference in the lengths of the candidate regions where no functional genes were observed may reflect the evolutionary time after fusion events.

EST sequencing and fosmid library construction in a non-model moth, Mamestra brassicae, for comparative mapping

Genome data are useful for both basic and applied research; however, it is difficult to carry out large-scale genome analyses using species with limited genetic or genomic resources. Here, we describe a cost-effective method to analyze the genome of a non-model species, using the cabbage moth, Mamestra brassicae (Lepidoptera: Noctuidae). First, we conducted expression sequence tag (EST) analysis. In this analysis, we performed PCR-based prescreening of a non-normalized embryonic cDNA library to eliminate already sequenced cDNAs from further sequencing, which significantly increased the percentage of unique genes. Next, we constructed a fosmid library of M. brassicae and isolated 120 clones containing 119 putative single copy genes by PCR-based screening with primer sets designed from the ESTs. Finally, we showed that the isolated fosmid clones could be used as probes for multicolor fluorescence in situ hybridization (FISH) analysis against an M. brassicae chromosome and confirmed conserved gene order between M. brassicae and the silkworm, Bombyx mori. Thus, we developed new genomic resources for comparative genome analysis in M. brassicae using robust and relatively low cost methods that can be applied to any non-model organism.

Application of high-resolution DNA melting for genotyping in lepidopteran non-model species: Ostrinia furnacalis (Crambidae)

Development of an ideal marker system facilitates a better understanding of the genetic diversity in lepidopteran non-model organisms, which have abundant species, but relatively limited genomic resources. Single nucleotide polymorphisms (SNPs) discovered within single-copy genes have proved to be desired markers, but SNP genotyping by current techniques remain laborious and expensive. High resolution melting (HRM) curve analysis represents a simple, rapid and inexpensive genotyping method that is primarily confined to clinical and diagnostic studies. In this study, we evaluated the potential of HRM analysis for SNP genotyping in the lepidopteran non-model species Ostrinia furnacalis (Crambidae). Small amplicon and unlabeled probe assays were developed for the SNPs, which were identified in 30 females of O. furnacalis from 3 different populations by our direct sequencing. Both assays were then applied to genotype 90 unknown female DNA by prior mixing with known wild-type DNA. The genotyping results were compared with those that were obtained using bi-directional sequencing analysis. Our results demonstrated the efficiency and reliability of the HRM assays. HRM has the potential to provide simple, cost-effective genotyping assays and facilitates genotyping studies in any non-model lepidopteran species of interest.

Sex-linked pheromone receptor genes of the European corn borer, Ostrinia nubilalis, are in tandem arrays

Background

Tuning of the olfactory system of male moths to conspecific female sex pheromones is crucial for correct species recognition; however, little is known about the genetic changes that drive speciation in this system. Moths of the genus Ostrinia are good models to elucidate this question, since significant differences in pheromone blends are observed within and among species. Odorant receptors (ORs) play a critical role in recognition of female sex pheromones; eight types of OR genes expressed in male antennae were previously reported in Ostrinia moths.

Methodology/Principal Findings

We screened an O. nubilalis bacterial artificial chromosome (BAC) library by PCR, and constructed three contigs from isolated clones containing the reported OR genes. Fluorescence in situ hybridization (FISH) analysis using these clones as probes demonstrated that the largest contig, which contained eight OR genes, was located on the Z chromosome; two others harboring two and one OR genes were found on two autosomes. Sequence determination of BAC clones revealed the Z-linked OR genes were closely related and tandemly arrayed; moreover, four of them shared 181-bp direct repeats spanning exon 7 and intron 7.

Conclusions/Significance

This is the first report of tandemly arrayed sex pheromone receptor genes in Lepidoptera. The localization of an OR gene cluster on the Z chromosome agrees with previous findings for a Z-linked locus responsible for O. nubilalis male behavioral response to sex pheromone. The 181-bp direct repeats might enhance gene duplications by unequal crossovers. An autosomal locus responsible for male response to sex pheromone in Heliothis virescens and H. subflexa was recently reported to contain at least four OR genes. Taken together, these findings support the hypothesis that generation of additional copies of OR genes can increase the potential for male moths to acquire altered specificity for pheromone components, and accordingly, facilitate differentiation of sex pheromones.