A dense genetic map of the silkworm, Bombyx mori, covering all chromosomes based on 1018 molecular markers

The fine-scale recombination rate variation and associations with genomic features in a butterfly

Recombination is a key molecular mechanism that has profound implications on both micro- and macroevolutionary processes. However, the determinants of recombination rate variation in holocentric organisms are poorly understood, in particular in Lepidoptera (moths and butterflies). The wood white butterfly (Leptidea sinapis) shows considerable intraspecific variation in chromosome numbers and is a suitable system for studying regional recombination rate variation and its potential molecular underpinnings. Here, we developed a large whole-genome resequencing data set from a population of wood whites to obtain high-resolution recombination maps using linkage disequilibrium information. The analyses revealed that larger chromosomes had a bimodal recombination landscape, potentially caused by interference between simultaneous chiasmata. The recombination rate was significantly lower in subtelomeric regions, with exceptions associated with segregating chromosome rearrangements, showing that fissions and fusions can have considerable effects on the recombination landscape. There was no association between the inferred recombination rate and base composition, supporting a limited influence of GC-biased gene conversion in butterflies. We found significant but variable associations between the recombination rate and the density of different classes of transposable elements, most notably a significant enrichment of short interspersed nucleotide elements in genomic regions with higher recombination rate. Finally, the analyses unveiled significant enrichment of genes involved in farnesyltranstransferase activity in recombination coldspots, potentially indicating that expression of transferases can inhibit formation of chiasmata during meiotic division. Our results provide novel information about recombination rate variation in holocentric organisms and have particular implications for forthcoming research in population genetics, molecular/genome evolution, and speciation.

Linkage mapping and genome annotation give novel insights into gene family expansions and regional recombination rate variation in the painted lady (Vanessa cardui) butterfly

Characterization of gene family expansions and crossing over is crucial for understanding how organisms adapt to the environment. Here, we develop a high-density linkage map and detailed genome annotation of the painted lady butterfly (Vanessa cardui) - a non-diapausing, highly polyphagous species famous for its long-distance migratory behavior and almost cosmopolitan distribution. Our results reveal a complex interplay between regional recombination rate variation, gene duplications and transposable element activity shaping the genome structure of the painted lady. We identify several lineage specific gene family expansions. Their functions are mainly associated with protein and fat metabolism, detoxification, and defense against infection - critical processes for the painted lady's unique life-history. Furthermore, the detailed recombination maps allow us to characterize the regional recombination landscape, data that reveal a strong effect of chromosome size on the recombination rate, a limited impact of GC-biased gene conversion and a positive association between recombination and short interspersed elements.

Gene mapping reveals the association between tyrosine protein kinase Abl1 and the silk yield of Bombyx mori

The Z chromosome of the silkworm contains a major gene that influences silk yield. This major locus on chromosome Z accounts for 35.10% of the phenotypic variance. The location and identification of the gene have been a focus of silkworm genetics research. Unfortunately, identification of this gene has been difficult. We used extreme phenotype subpopulations and selected from a backcross population, BC₁ M, which was obtained using the high-yield strain 872B and the low-yield strain IS-Dazao as parents, for mapping the gene on the chromosome Z. The candidate region was narrowed down to 134 kb at the tip of the chromosome. BmAbl1 in this region correlated with silk gland development by spatiotemporal expression analysis. This gene was differentially expressed in the posterior silk glands of the high- and low-yield strains. In BmAbl1, an insertion-deletion (indel) within the 10th exonic region and an SNP within the 6th intronic region were detected and shown to be associated with cocoon shell weight in 84 Bombyx mori strains with different yields. Nucleotide diversity analysis of BmAbl1 and its 50 kb flanking regions indicated that BmAbl1 has experienced strong artificial selection during silkworm domestication. This study is the first to identify the genes controlling silk yield in the major QTL of the Z chromosome using forward genetics.

The Theory and Applications of Measuring Broad-Range and Chromosome-Wide Recombination Rate from Allele Frequency Decay around a Selected Locus

Recombination is the exchange of genetic material between homologous chromosomes via physical crossovers. High-throughput sequencing approaches detect crossovers genome wide to produce recombination rate maps but are difficult to scale as they require large numbers of recombinants individually sequenced. We present a simple and scalable pooled-sequencing approach to experimentally infer near chromosome-wide recombination rates by taking advantage of non-Mendelian allele frequency generated from a fitness differential at a locus under selection. As more crossovers decouple the selected locus from distal loci, the distorted allele frequency attenuates distally toward Mendelian and can be used to estimate the genetic distance. Here, we use marker selection to generate distorted allele frequency and theoretically derive the mathematical relationships between allele frequency attenuation, genetic distance, and recombination rate in marker-selected pools. We implemented nonlinear curve-fitting methods that robustly estimate the allele frequency decay from batch sequencing of pooled individuals and derive chromosome-wide genetic distance and recombination rates. Empirically, we show that marker-selected pools closely recapitulate genetic distances inferred from scoring recombinants. Using this method, we generated novel recombination rate maps of three wild-derived strains of Drosophila melanogaster, which strongly correlate with previous measurements. Moreover, we show that this approach can be extended to estimate chromosome-wide crossover interference with reciprocal marker selection and discuss how it can be applied in the absence of visible markers. Altogether, we find that our method is a simple and cost-effective approach to generate chromosome-wide recombination rate maps requiring only one or two libraries.

Identification of Genes that Control Silk Yield by RNA Sequencing Analysis of Silkworm (Bombyx mori) Strains of Variable Silk Yield

Silk is an important natural fiber of high economic value, and thus genetic study of the silkworm is a major area of research. Transcriptome analysis can provide guidance for genetic studies of silk yield traits. In this study, we performed a transcriptome comparison using multiple silkworms with different silk yields. A total of 22 common differentially expressed genes (DEGs) were identified in multiple strains and were mainly involved in metabolic pathways. Among these, seven significant common DEGs were verified by quantitative reverse transcription polymerase chain reaction, and the results coincided with the findings generated by RNA sequencing. Association analysis showed that BGIBMGA003330 and BGIBMGA005780 are significantly associated with cocoon shell weight and encode uridine nucleosidase and small heat shock protein, respectively. Functional annotation of these genes suggest that these play a role in silkworm silk gland development or silk protein synthesis. In addition, we performed principal component analysis (PCA) in combination with wild silkworm analysis, which indicates that modern breeding has a stronger selection effect on silk yield traits than domestication, and imply that silkworm breeding induces aggregation of genes related to silk yield.

Genetic variation in ecoraces of tropical tasar silkworm, Antheraea mylitta using SSR markers

The tropical tasar silkworm, Antheraea mylitta, polyphagous sericigenous insect mostly found in the tropical areas of India. It is found in these regions as ecotypes or ecoraces. It feeds primarily on plants, a variety of secondary plants like Terminalia arjuna and T. tomentosa. Tasar culture is a traditional livelihood for lakhs of tribal populace in the areas of Jharkhand, Chhatisgarh, Orissa, Maharashtra, Andhra Pradesh, West Bengal and Uttar Pradesh. In the present study, the genetic diversity of these ecoraces is identified by DNA markers, namely simple sequence repeats (SSRs), most of which produced polymorphic bands.

Targeted Mutagenesis in Bombyx mori Using TALENs

Bombyx mori is a valuable model organism of high economic importance. Its genome sequence is available, as well as basic genetic and molecular genetic tools and markers. The introduction of genome editing methods based on engineered nucleases enables precise manipulations with genomic DNA, including targeted DNA deletions, insertions, or replacements in the genome allowing gene analysis and various applications. We describe here the use of TALENs which have a simple modular design of their DNA-binding domains, are easy to prepare and proved to be efficient in targeting of a wide range of cleavage sites. Our procedure often allows the production of individuals carrying homozygous mutations as early as in the G1 generation.

Allelic Variability in the Intronic Region of the Fibroin Heavy-Chain Gene in Silkworm Bombyx mori L. Strains of Brazilian Germplasm Bank

The domesticated silkworm Bombyx mori L. is currently found only in germplasm banks. Therefore, characterization and conservation of this genetic resource is crucial. Based on previous studies that revealed nucleotide differences in silkworm strains, the intron of the fibroin heavy chain gene (H-fib) can be used for molecular silkworm characterization. The H-fib gene has two exons and a unique intron, and encodes the principal component of the silk fiber, the fibroin heavy chain. Therefore, this study aimed to identify the genetic variability of the unique intron of H-fib gene of 20 silkworm strains maintained at the Universidade Estadual de Maringá Brazilian Germplasm Bank (UBGB) by conformation-sensitive gel electrophoresis (CSGE) and nucleotide sequencing. Genomic DNA extracted from silkworm moths was PCR amplified. CSGE revealed that most of the analyzed silkworm strains had only homoduplex molecules. However, DNA from the Japanese strains B106, B82, and M12-2 had two extra DNA fragments produced by heteroduplex molecules, revealing variation between alleles. Sequencing of the H-fib intron was used to confirm the variation previously detected by CSGE and detected a significant polymorphism characterized by a 17-base pair (bp) deletion, a 2-bp insertion, and eight nucleotide substitutions. Although genetic and allelic variability was detected in some silkworm strains, the intron of the H-fib gene revealed not to be the best molecular marker for the characterization of B. mori strains from UBGB.

Genetic characterization of Bombyx mori (Lepidoptera: Bombycidae) breeding and hybrid lines with different geographic origins

The domesticated silkworm Bombyx mori L. comprises a large number of geographical breeds and hybrid lines. Knowing the genetic structure of those may provide information to improve the conservation of commercial lines by estimating inbreeding over generations and the consequences of excessive use of those lineages. Here, we analyzed the genetic diversity of seven breeds and eight hybrid lines from Eastern Europe and Asia using highly polymorphic microsatellites markers to determine its genetical impact on their use in global breeding programs. No consistent pattern of deviation from Hardy-Weinberg equilibrium was found for most breed and hybrids; and the absence of a linkage disequilibrium also suggests that the strains are in equilibrium. A principal coordinate analysis revealed a clear separation of two silkworm breeds from the rest: one (IBV) originated from India and the other one (RG90) from Romania/Japan. The tendency of the other breeds from different geographic origins to cluster together in a general mix might be due to similar selection pressures (climate and anthropogenic factors) in different geographic locations. Phylogenetic analyses grouped the different silkworm breeds but not the hybrids according to their geographic origin and confirmed the pattern found in the principal coordinate analysis.

Lep-MAP: fast and accurate linkage map construction for large SNP datasets

MOTIVATION

Current high-throughput sequencing technologies allow cost-efficient genotyping of millions of single nucleotide polymorphisms (SNPs) for hundreds of samples. However, the tools that are currently available for constructing linkage maps are not well suited for large datasets. Linkage maps of large datasets would be helpful in de novo genome assembly by facilitating comprehensive genome validation and refinement by enabling chimeric scaffold detection, as well as in family-based linkage and association studies, quantitative trait locus mapping, analysis of genome synteny and other complex genomic data analyses.

RESULTS

We describe a novel tool, called Lepidoptera-MAP (Lep-MAP), for constructing accurate linkage maps with ultradense genome-wide SNP data. Lep-MAP is fast and memory efficient and largely automated, requiring minimal user interaction. It uses simultaneously data on multiple outbred families and can increase linkage map accuracy by taking into account achiasmatic meiosis, a special feature of Lepidoptera and some other taxa with no recombination in one sex (no recombination in females in Lepidoptera). We demonstrate that Lep-MAP outperforms other methods on real and simulated data. We construct a genome-wide linkage map of the Glanville fritillary butterfly (Melitaea cinxia) with over 40 000 SNPs. The data were generated with a novel in-house SOLiD restriction site-associated DNA tag sequencing protocol, which is described in the online supplementary material.

AVAILABILITY AND IMPLEMENTATION

Java source code under GNU general public license with the compiled classes and the datasets are available from http://sourceforge.net/users/lep-map.

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.

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.

A new mapping method for quantitative trait loci of silkworm

BACKGROUND

Silkworm is the basis of sericultural industry and the model organism in insect genetics study. Mapping quantitative trait loci (QTLs) underlying economically important traits of silkworm is of high significance for promoting the silkworm molecular breeding and advancing our knowledge on genetic architecture of the Lepidoptera. Yet, the currently used mapping methods are not well suitable for silkworm, because of ignoring the recombination difference in meiosis between two sexes.

RESULTS

A mixed linear model including QTL main effects, epistatic effects, and QTL × sex interaction effects was proposed for mapping QTLs in an F2 population of silkworm. The number and positions of QTLs were determined by F-test and model selection. The Markov chain Monte Carlo (MCMC) algorithm was employed to estimate and test genetic effects of QTLs and QTL × sex interaction effects. The effectiveness of the model and statistical method was validated by a series of simulations. The results indicate that when markers are distributed sparsely on chromosomes, our method will substantially improve estimation accuracy as compared to the normal chiasmate F2 model. We also found that a sample size of hundreds was sufficiently large to unbiasedly estimate all the four types of epistases (i.e., additive-additive, additive-dominance, dominance-additive, and dominance-dominance) when the paired QTLs reside on different chromosomes in silkworm.

CONCLUSION

The proposed method could accurately estimate not only the additive, dominance and digenic epistatic effects but also their interaction effects with sex, correcting the potential bias and precision loss in the current QTL mapping practice of silkworm and thus representing an important addition to the arsenal of QTL mapping tools.

An ABC transporter mutation is correlated with insect resistance to Bacillus thuringiensis Cry1Ac toxin

Transgenic crops producing insecticidal toxins from Bacillus thuringiensis (Bt) are commercially successful in reducing pest damage, yet knowledge of resistance mechanisms that threaten their sustainability is incomplete. Insect resistance to the pore-forming Cry1Ac toxin is correlated with the loss of high-affinity, irreversible binding to the mid-gut membrane, but the genetic factors responsible for this change have been elusive. Mutations in a 12-cadherin-domain protein confer some Cry1Ac resistance but do not block this toxin binding in in vitro assays. We sought to identify mutations in other genes that might be responsible for the loss of binding. We employed a map-based cloning approach using a series of backcrosses with 1,060 progeny to identify a resistance gene in the cotton pest Heliothis virescens that segregated independently from the cadherin mutation. We found an inactivating mutation of the ABC transporter ABCC2 that is genetically linked to Cry1Ac resistance and is correlated with loss of Cry1Ac binding to membrane vesicles. ABC proteins are integral membrane proteins with many functions, including export of toxic molecules from the cell, but have not been implicated in the mode of action of Bt toxins before. The reduction in toxin binding due to the inactivating mutation suggests that ABCC2 is involved in membrane integration of the toxin pore. Our findings suggest that ABC proteins may play a key role in the mode of action of Bt toxins and that ABC protein mutations can confer high levels of resistance that could threaten the continued utilization of Bt-expressing crops. However, such mutations may impose a physiological cost on resistant insects, by reducing export of other toxins such as plant secondary compounds from the cell. This weakness could be exploited to manage this mechanism of Bt resistance in the field.

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

Lepidoptera, butterflies and moths, is the second largest animal order and includes numerous agricultural pests. To facilitate comparative genomics in Lepidoptera, we isolated BAC clones containing conserved and putative single-copy genes from libraries of three pests, Heliothis virescens, Ostrinia nubilalis, and Plutella xylostella, harboring the haploid chromosome number, n = 31, which are not closely related with each other or with the silkworm, Bombyx mori, (n = 28), the sequenced model lepidopteran. A total of 108–184 clones representing 101–182 conserved genes were isolated for each species. For 79 genes, clones were isolated from more than two species, which will be useful as common markers for analysis using fluorescence in situ hybridization (FISH), as well as for comparison of genome sequence among multiple species. The PCR-based clone isolation method presented here is applicable to species which lack a sequenced genome but have a significant collection of cDNA or EST sequences.

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.

Amplified fragment length polymorphism mapping of quantitative trait loci for economically important traits in the silkworm, Bombyx mori

Cocoon related characteristics are economically important traits in the silkworm, Bombyx mori L. (Lepidoptera: Bombycidae). In this study a genetic linkage map was developed that identified QTL controlling the cocoon weight, cocoon shell weight, and cocoon shell percentage using 161 amplified fragment length polymorphism (AFLP) markers. Twenty PstI/TaqI primer combinations were employed to genotype 78 F(2) progenies derived from a cross between P107 Japanese inbred line and Khorasan Lemon Iranian native strain. Among polymorphic markers, 159 AFLP markers were assigned to 24 linkage groups at the LOD threshold of 2.5 that varied in length from 4 to 299 cM. The total length of the linkage map was 2747 cM, giving an average marker resolution of 19.31 cM. A total of 21 AFLP markers were identified that were distributed over the ten linkage groups linked to the three studied traits using the composite interval mapping method. The explained variation rate by QTL controlling cocoon weight, cocoon shell weight, and cocoon shell percentage ranged from 0.02% to 64.85%, 0.2% to 49.11%, and 0.04% to 84.20%, respectively. These QTL controlled by different actions as well as under dominance, additive, partial dominance, dominance, and over dominance.

Toyama Kametaro and Vernon Kellogg: silkworm inheritance experiments in Japan, Siam, and the United States, 1900-1912

Japanese agricultural scientist Toyama Kametaro's report about the Mendelian inheritance of silkworm cocoon color in Studies on the Hybridology of Insects (1906) spurred changes in Japanese silk production and thrust Toyama and his work into a scholarly exchange with American entomologist Vernon Kellogg. Toyama's work, based on research conducted in Japan and Siam, came under international scrutiny at a time when analyses of inheritance flourished after the "rediscovery" of Mendel's laws of heredity in 1900. The hybrid silkworm studies in Asia attracted the attention of Kellogg, who was concerned with how experimental biology would be used to study the causes of natural selection. He challenged Toyama's conclusions that Mendelism alone could explain the inheritance patterns of silkworm characters such as cocoon color because they had been subject to hundreds of years of artificial selection, or breeding. This examination of the intersection of Japanese sericulture and American entomology probes how practical differences in scientific interests, societal responsibilities, and silkworm materiality were negotiated throughout the processes of legitimating Mendelian genetics on opposite sides of the Pacific. The ways in which Toyama and Kellogg assigned importance to certain silkworm properties show how conflicting intellectual orientations arose in studies of the same organism. Contestation about Mendelism took place not just on a theoretical level, but the debate was fashioned through each scientist's rationale about the categorization of silkworm breeds and races and what counted as "natural". This further mediated the acceptability of the silkworm not as an experimental organism, but as an appropriately "natural" insect with which to demonstrate laws of inheritance. All these shed light on the challenges that came along with the use of agricultural animals to convincingly articulate new biological principles.

Extensive conserved synteny of genes between the karyotypes of Manduca sexta and Bombyx mori revealed by BAC-FISH mapping

Background

Genome sequencing projects have been completed for several species representing four highly diverged holometabolous insect orders, Diptera, Hymenoptera, Coleoptera, and Lepidoptera. The striking evolutionary diversity of insects argues a need for efficient methods to apply genome information from such models to genetically uncharacterized species. Constructing conserved synteny maps plays a crucial role in this task. Here, we demonstrate the use of fluorescence in situ hybridization with bacterial artificial chromosome probes as a powerful tool for physical mapping of genes and comparative genome analysis in Lepidoptera, which have numerous and morphologically uniform holokinetic chromosomes.

Methodology/Principal Findings

We isolated 214 clones containing 159 orthologs of well conserved single-copy genes of a sequenced lepidopteran model, the silkworm, Bombyx mori, from a BAC library of a sphingid with an unexplored genome, the tobacco hornworm, Manduca sexta. We then constructed a BAC-FISH karyotype identifying all 28 chromosomes of M. sexta by mapping 124 loci using the corresponding BAC clones. BAC probes from three M. sexta chromosomes also generated clear signals on the corresponding chromosomes of the convolvulus hawk moth, Agrius convolvuli, which belongs to the same subfamily, Sphinginae, as M. sexta.

Conclusions/Significance

Comparison of the M. sexta BAC physical map with the linkage map and genome sequence of B. mori pointed to extensive conserved synteny including conserved gene order in most chromosomes. Only a few rearrangements, including three inversions, three translocations, and two fission/fusion events were estimated to have occurred after the divergence of Bombycidae and Sphingidae. These results add to accumulating evidence for the stability of lepidopteran genomes. Generating signals on A. convolvuli chromosomes using heterologous M. sexta probes demonstrated that BAC-FISH with orthologous sequences can be used for karyotyping a wide range of related and genetically uncharacterized species, significantly extending the ability to develop synteny maps for comparative and functional genomics.

An integrated genetic linkage map for silkworms with three parental combinations and its application to the mapping of single genes and QTL

Background

Bombyx mori, the domesticated silkworm, is a well-studied model insect with great economic and scientific significance. Although more than 400 mutations have been described in silkworms, most have not been identified, especially those affecting economically-important traits. Simple sequence repeats (SSRs) are effective and economical tools for mapping traits and genetic improvement. The current SSR linkage map is of low density and contains few polymorphisms. The purpose of this work was to develop a dense and informative linkage map that would assist in the preliminary mapping and dissection of quantitative trait loci (QTL) in a variety of silkworm strains.

Results

Through an analysis of > 50,000 genotypes across new mapping populations, we constructed two new linkage maps covering 27 assigned chromosomes and merged the data with previously reported data sets. The integrated consensus map contains 692 unique SSR sites, improving the density from 6.3 cM in the previous map to 4.8 cM. We also developed 497 confirmed neighboring markers for corresponding low-polymorphism sites, with 244 having polymorphisms. Large-scale statistics on the SSR type were suggestive of highly efficient markers, based upon which we searched 16,462 available genomic scaffolds for SSR loci. With the newly constructed map, we mapped single-gene traits, the QTL of filaments, and a number of ribosomal protein genes.

Conclusion

The integrated map produced in this study is a highly efficient genetic tool for the high-throughput mapping of single genes and QTL. Compared to previous maps, the current map offers a greater number of markers and polymorphisms; thus, it may be used as a resource for marker-assisted breeding.

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.

Inheritance and linkage analysis of co-dominant SSR markers on the Z chromosome of the silkworm (Bombyx mori L.)

Microsatellites or simple sequence repeats (SSRs) are co-dominant molecular markers. When we used fluorescent SSR markers to construct a linkage map for the female heterogametic silkworm (Bombyx mori, ZW), we found that some loci did not segregate in a Mendelian ratio of 1:1 in a backcross population. These loci segregated in a 3:1 ratio of single bands compared with double bands. Further examination of band patterns indicated that three types of SSR bands were present: two homozygotes and one heterozygote. In the beginning, we considered to discard these markers. By scoring male and female F1 individuals, we confirmed that these loci were located on the Z chromosome. Using the sex-linked visible mutation sch (K05) and its wild-type (C108), we constructed an F1 male backcross (BC1M) mapping population. The combination of sch backcross and SSR data enabled us to map the SSR markers to the Z chromosome. By adjusting input parameters based on these data, we were able to use Mapmaker software to construct a linkage map. This strategy takes advantage of co-dominant markers for positional cloning of genes on the Z chromosome. We localized sch to the Z chromosome relative to six SSR markers and one PCR marker, covering a total of 76.1 cM. The sch mutation is an important sex-linked visible mutation widely used in breeding of commercial silkworms (e.g. male silkworm selection rearing). Localization of the sch gene may prove helpful in cloning the gene and developing strains for marker-assisted selection in silkworm breeding.

From genome to proteome: great progress in the domesticated silkworm (Bombyx mori L.)

As the only truly domesticated insect, the silkworm not only has great economic value, but it also has value as a model for genetics and molecular biology research. Genomics and proteomics have recently shown vast potential to be essential tools in domesticated silkworm research, especially after the completion of the Bombyx mori genome sequence. This paper reviews the progress of the domesticated silkworm genome, particularly focusing on its genetic map, physical map and functional genome. This review also presents proteomics, the proteomic technique and its application in silkworm research.

A BAC-based integrated linkage map of the silkworm Bombyx mori

An integrated map of the Bombyx mori genome has been constructed using 361.1 Mb of BAC contigs and singletons together with a genetic map containing 1689 independent genes and synteny among Apis, Tribolium, and Bombyx was examined.

Background

In 2004, draft sequences of the model lepidopteran Bombyx mori were reported using whole-genome shotgun sequencing. Because of relatively shallow genome coverage, the silkworm genome remains fragmented, hampering annotation and comparative genome studies. For a more complete genome analysis, we developed extended scaffolds combining physical maps with improved genetic maps.

Results

We mapped 1,755 single nucleotide polymorphism (SNP) markers from bacterial artificial chromosome (BAC) end sequences onto 28 linkage groups using a recombining male backcross population, yielding an average inter-SNP distance of 0.81 cM (about 270 kilobases). We constructed 6,221 contigs by fingerprinting clones from three BAC libraries digested with different restriction enzymes, and assigned a total of 724 single copy genes to them by BLAST (basic local alignment search tool) search of the BAC end sequences and high-density BAC filter hybridization using expressed sequence tags as probes. We assigned 964 additional expressed sequence tags to linkage groups by restriction fragment length polymorphism analysis of a nonrecombining female backcross population. Altogether, 361.1 megabases of BAC contigs and singletons were integrated with a map containing 1,688 independent genes. A test of synteny using Oxford grid analysis with more than 500 silkworm genes revealed six versus 20 silkworm linkage groups containing eight or more orthologs of Apis versus Tribolium, respectively.

Conclusion

The integrated map contains approximately 10% of predicted silkworm genes and has an estimated 76% genome coverage by BACs. This provides a new resource for improved assembly of whole-genome shotgun data, gene annotation and positional cloning, and will serve as a platform for comparative genomics and gene discovery in Lepidoptera and other insects.

Genetic mapping of a food preference gene in the silkworm, Bombyx mori, using restriction fragment length polymorphisms (RFLPs)

The domesticated silkworm, Bombyx mori, has strict food preferences and grows by feeding on mulberry leaves. However, "Sawa-J", an abnormal feeding habit strain selected from the genetic stock, feeds on an artificial diet without mulberry leaf powder. In this study, the food preference gene in Sawa-J was genetically identified using restriction fragment length polymorphisms (RFLPs) of a cDNA clone on each linkage group. Taking advantage of a lack of genetic recombination in females, reciprocal backcrossed F1 (BC1) progenies were independently prepared using a non-feeding strain, C108, as a mating partner of Sawa-J. Our results of linkage analysis and mapping proved that the feeding behavior is primarily controlled by a major recessive gene mapped at 20.2 cM on RFLP linkage group 9 (RFLG9), and clone e73 at a distance of 4.2 cM was found as the first linked molecular marker.

Identification of the female-determining region of the W chromosome in Bombyx mori

The W chromosome of the silkworm Bombyx mori is devoid of functional genes, except for the putative female-determining gene (Fem). To localize Fem, we investigated the presence of W-specific DNA markers on strains in which an autosomal fragment containing dominant marker genes was attached to the W chromosome. We produced new W-chromosomal fragments from the existing Zebra-W strain (T(W;3)Ze chromosome) by X-irradiation, and then carried out deletion mapping of these and sex-limited yellow cocoon strains (T(W;2)Y-Chu, -Abe and -Ban types) from different Japanese stock centers. Of 12 RAPD markers identified in the normal W chromosomes of most silkworm strains in Japan, the newly irradiated W(B-YL-YS)Ze chromosome contained three, the T(W;2)Y-Chu chromosome contained six, and the T(W;2)Y-Abe and -Ban chromosomes contained only one (W-Rikishi). To investigate the ability of the reduced W-chromosome translocation fragments to form heterochromatin bodies, which are found in nuclei of normal adult female sucking stomachs, we examined cells of the normal type p50 strain and the T(W;2)Y-Chu and -Abe strains. A single sex heterochromatin body was found in nuclei of p50 females, whereas we detected only small sex heterochromatin bodies in the T(W;2)Y-Chu strain and no sex heterochromatin body in the T(W;2)Y-Abe strain. Since adult females of all strains were normal and fertile, we conclude that only extremely limited region, containing the W-Rikishi RAPD sequence of the W chromosome, is required to determine femaleness. Based on a comparison of the normal W-chromosome and 7 translocation and W-deletion strains we present a map of Fem relative to the 12 W-specific RAPD markers.

Synteny and chromosome evolution in the lepidoptera: evidence from mapping in Heliconius melpomene

The extent of conservation of synteny and gene order in the Lepidoptera has been investigated previously only by comparing a small subset of linkage groups between the moth Bombyx mori and the butterfly Heliconius melpomene. Here we report the mapping of 64 additional conserved genes in H. melpomene, which contributed 47 markers to a comparative framework of 72 orthologous loci spanning all 21 H. melpomene chromosomes and 27 of the 28 B. mori chromosomes. Comparison of the maps revealed conserved synteny across all chromosomes for the 72 loci, as well as evidence for six cases of chromosome fusion in the Heliconius lineage that contributed to the derived 21-chromosome karyotype. Comparisons of gene order on these fused chromosomes revealed two instances of colinearity between H. melpomene and B. mori, but also one instance of likely chromosomal rearrangement. B. mori is the first lepidopteran species to have its genome sequenced, and the finding that there is conserved synteny and gene order among Lepidoptera indicates that the genomic tools developed in B. mori will be broadly useful in other species.

High recombination frequency creates genotypic diversity in colonies of the leaf-cutting ant Acromyrmex echinatior

Honeybees are known to have genetically diverse colonies because queens mate with many males and the recombination rate is extremely high. Genetic diversity among social insect workers has been hypothesized to improve general performance of large and complex colonies, but this idea has not been tested in other social insects. Here, we present a linkage map and an estimate of the recombination rate for Acromyrmex echinatior, a leaf-cutting ant that resembles the honeybee in having multiple mating of queens and colonies of approximately the same size. A map of 145 AFLP markers in 22 linkage groups yielded a total recombinational size of 2076 cM and an inferred recombination rate of 161 kb cM(-1) (or 6.2 cM Mb(-1)). This estimate is lower than in the honeybee but, as far as the mapping criteria can be compared, higher than in any other insect mapped so far. Earlier studies on A. echinatior have demonstrated that variation in division of labour and pathogen resistance has a genetic component and that genotypic diversity among workers may thus give colonies of this leaf-cutting ant a functional advantage. The present result is therefore consistent with the hypothesis that complex social life can select for an increased recombination rate through effects on genotypic diversity and colony performance.

Primary genetic linkage maps of the ascidian, Ciona intestinalis

For whole-genome analysis in a basal chordate (protochordate), we used F1 pseudo-testcross mapping strategy and amplified fragment length polymorphism (AFLP) markers to construct primary linkage maps of the ascidian tunicate Ciona intestinalis. Two genetic maps consisted of 14 linkage groups, in agreement with the haploid chromosome number, and contained 276 and 125 AFLP loci derived from crosses between British and Neapolitan individuals. The two maps covered 4218.9 and 2086.9 cM, respectively, with an average marker interval of 16.1 and 18.9 cM. We observed a high recombinant ratio, ranging from 25 to 49 kb/cM, which can explain the high degree of polymorphism in this species. Some AFLP markers were converted to sequence tagged sites (STSs) by sequence determination, in order to create anchor markers for the fragmental physical map. Our recombination tools provide basic knowledge of genetic status and whole genome organization, and genetic markers to assist positional cloning in C. intestinalis.

Molecular cloning and chromosomal localization of the Bombyx Sex-lethal gene

We cloned Bm-Sxl, an orthologue of the Drosophila melanogaster Sex-lethal (Sxl) gene from embryos of Bombyx mori. The full-length cDNAs were of 2 sizes, 1528 and 1339 bp, and were named Bm-Sxl-L and Bm-Sxl-S, respectively. Bm-Sxl-L consists of 8 exons and spans more than 20 kb of genomic DNA. The open reading frame (ORF) codes for a protein 336 amino acids in length. Bm-Sxl-S is a splice variant that lacks the second exon. This creates a new translation start 138 nucleotides downstream and an ORF that codes for 46 amino acids fewer at the N-terminus. Linkage analysis using an F2 panel mapped Bm-Sxl to linkage group 16 at 69.8 cM. We isolated 2 BACs that include the Bm-Sxl gene. With BAC-FISH we located Bm-Sxl cytogenetically on the chromosome corresponding to linkage group 16 (LG16) at position >68.8 cM.

Molecular marker systems in insects: current trends and future avenues

Insects comprise the largest species composition in the entire animal kingdom and possess a vast undiscovered genetic diversity and gene pool that can be better explored using molecular marker techniques. Current trends of application of DNA marker techniques in diverse domains of insect ecological studies show that mitochondrial DNA (mtDNA), microsatellites, random amplified polymorphic DNA (RAPD), expressed sequence tags (EST) and amplified fragment length polymorphism (AFLP) markers have contributed significantly for progresses towards understanding genetic basis of insect diversity and for mapping medically and agriculturally important genes and quantitative trait loci in insect pests. Apart from these popular marker systems, other novel approaches including transposon display, sequence-specific amplification polymorphism (S-SAP), repeat-associated polymerase chain reaction (PCR) markers have been identified as alternate marker systems in insect studies. Besides, whole genome microarray and single nucleotide polymorphism (SNP) assays are becoming more popular to screen genome-wide polymorphisms in fast and cost effective manner. However, use of such methodologies has not gained widespread popularity in entomological studies. The current study highlights the recent trends of applications of molecular markers in insect studies and explores the technological advancements in molecular marker tools and modern high throughput genotyping methodologies that may be applied in entomological researches for better understanding of insect ecology at molecular level.

A chicken linkage map based on microsatellite markers genotyped on a Japanese Large Game and White Leghorn cross

A detailed linkage map is necessary for efficient detection of quantitative trait loci (QTL) in chicken resource populations. In this study, microsatellite markers isolated from a (CA)n-enriched library (designated as ABR Markers) were mapped using a population developed from a cross between Japanese Game and White Leghorn chickens. In total, 296 markers including 193 ABR, 43 MCW, 31 ADL, 22 LEI, 3 HUJ, 2 GCT, 1 UMA and 1 ROS were mapped by linkage to chicken chromosomes 1-14, 17-21, 23, 24, 26-28 and Z. In addition, five markers were assigned to the map based on the chicken draft genomic sequence, bringing the total number of markers on the map to 301. The resulting linkage map will contribute to QTL mapping in chicken.

A second-generation integrated map of the silkworm reveals synteny and conserved gene order between lepidopteran insects

A second-generation linkage map was constructed for the silkworm, Bombyx mori, focusing on mapping Bombyx sequences appearing in public nucleotide databases and bacterial artificial chromosome (BAC) contigs. A total of 874 BAC contigs containing 5067 clones (22% of the library) were constructed by PCR-based screening with sequence-tagged sites (STSs) derived from whole-genome shotgun (WGS) sequences. A total of 523 BAC contigs, including 342 independent genes registered in public databases and 85 expressed sequence tags (ESTs), were placed onto the linkage map. We found significant synteny and conserved gene order between B. mori and a nymphalid butterfly, Heliconius melpomene, in four linkage groups (LGs), strongly suggesting that using B. mori as a reference for comparative genomics in Lepidotera is highly feasible.

Construction of a single nucleotide polymorphism linkage map for the silkworm, Bombyx mori, based on bacterial artificial chromosome end sequences

We have developed a linkage map for the silkworm Bombyx mori based on single nucleotide polymorphisms (SNPs) between strains p50T and C108T initially found on regions corresponding to the end sequences of bacterial artificial chromosome (BAC) clones. Using 190 segregants from a backcross of a p50T female x an F1 (p50T x C108T) male, we analyzed segregation patterns of 534 SNPs between p50T and C108T, detected among 3840 PCR amplicons, each associated with a p50T BAC end sequence. This enabled us to construct a linkage map composed of 534 SNP markers spanning 1305 cM in total length distributed over the expected 28 linkage groups. Of the 534 BACs whose ends harbored the SNPs used to construct the linkage map, 89 were associated with 107 different ESTs. Since each of the SNP markers is directly linked to a specific genomic BAC clone and to whole-genome sequence data, and some of them are also linked to EST data, the SNP linkage map will be a powerful tool for investigating silkworm genome properties, mutation mapping, and map-based cloning of genes of industrial and agricultural interest.

Genetic diversity among silkworm (Bombyx mori L., Lep., Bombycidae) germplasms revealed by microsatellites

To determine genetic relationships among strains of silkworm, Bombyx mori L., 31 strains with different origins, number of generations per year, number of molts per generation, and morphological characters were studied using simple sequence repeat (SSR) markers. Twenty-six primer pairs flanking microsatellite sequences in the silkworm genome were assayed. All were polymorphic and unambiguously separated silkworm strains from each other. A total of 188 alleles were detected with a mean value of 7.2 alleles/locus (range 2-17). The average heterozygosity value for each SSR locus ranged from 0 to 0.60, and the highest one was 0.96 (Fl0516 in 4013). The mean polymorphism index content (PIC) was 0.66 (range 0.12-0.89). Unweighted pair group method with arithmetic means (UPGMA) cluster analysis of Nei's genetic distance grouped silkworm strains based on their origin. Seven major ecotypic silkworm groups were analyzed. Principal components analysis (PCA) for SSR data support their UPGMA clustering. The results indicated that SSR markers are an efficient tool for fingerprinting cultivars and conducting genetic-diversity studies in the silkworm.

Localization of Müllerian mimicry genes on a dense linkage map of Heliconius erato

We report a dense genetic linkage map of Heliconius erato, a neotropical butterfly that has undergone a remarkable adaptive radiation in warningly colored mimetic wing patterns. Our study exploited natural variation segregating in a cross between H. erato etylus and H. himera to localize wing color pattern loci on a dense linkage map containing amplified fragment length polymorphisms (AFLP), microsatellites, and single-copy nuclear loci. We unambiguously identified all 20 autosomal linkage groups and the sex chromosome (Z). The map spanned a total of 1430 Haldane cM and linkage groups varied in size from 26.3 to 97.8 cM. The average distance between markers was 5.1 cM. Within this framework, we localized two major color pattern loci to narrow regions of the genome. The first gene, D, responsible for red/orange elements, had a most likely placement in a 6.7-cM region flanked by two AFLP markers on the end of a large 87.5-cM linkage group. The second locus, Sd, affects the melanic pattern on the forewing and was found within a 6.3-cM interval between flanking AFLP loci. This study complements recent linkage analysis of H. erato's comimic, H. melpomene, and forms the basis for marker-assisted physical mapping and for studies into the comparative genetic architecture of wing-pattern mimicry in Heliconius.

DNA fingerprinting using AFLP markers to search for markers associated with yield attributes in the silkworm, Bombyx mori

This study was carried out on 11 Chinese and 12 Japanese silkworm strains maintained by the Center for the Technological Development of Sericulture (CDTS) germplasm bank, located in Pereira, Colombia. The goals were to determine the genetic population structure of the two groups and the association between molecular markers (AFLPs) and important productivity characters. Group analysis showed the separation of the strains according to their geographic origin. The molecular markers and the productivity characters were correlated by multiple variance analysis. The analysis permitted the identification of molecular markers associated with the cocoon weight or the shell weight separately. Some markers were associated with both characters.

Integration of molecular and classical linkage groups of the silkworm, Bombyx mori (n = 28)

Previously published linkage groups (LGs) composed of molecular markers were assigned to classical LGs in the silkworm, Bombyx mori (n = 28). Four markers from the classical linkage map, og, w-1, Lp, and Pfl, were assigned to the molecular linkage maps using sequence tagged sites. In addition, linkage analysis was carried out using BF1 progeny between wild-type and mutant stocks carrying morphological phenotypic markers. As a result, the counterparts for 26 of 28 molecular LGs were identified with their counterparts of the classical LGs. Two visible markers, Sel and Xan, representing different classical LGs, were found to be linked.

Simple sequence repeat-based consensus linkage map of Bombyx mori

We established a genetic linkage map employing 518 simple sequence repeat (SSR, or microsatellite) markers for Bombyx mori (silkworm), the economically and culturally important lepidopteran insect, as part of an international genomics program. A survey of six representative silkworm strains using 2,500 (CA)n- and (CT)n-based SSR markers revealed 17-24% polymorphism, indicating a high degree of homozygosity resulting from a long history of inbreeding. Twenty-nine SSR linkage groups were established in well characterized Dazao and C108 strains based on genotyping of 189 backcross progeny derived from an F(1) male mated with a C108 female. The clustering was further focused to 28 groups by genotyping 22 backcross progeny derived from an F(1) female mated with a C108 male. This set of SSR linkage groups was further assigned to the 28 chromosomes (established linkage groups) of silkworm aided by visible mutations and cleaved amplified polymorphic sequence markers developed from previously mapped genes, cDNA sequences, and cloned random amplified polymorphic DNAs. By integrating a visible mutation p (plain, larval marking) and 29 well conserved genes of insects onto this SSR-based linkage map, a second generation consensus silkworm genetic map with a range of 7-40 markers per linkage group and a total map length of approximately 3431.9 cM was constructed and its high efficiency for genotyping and potential application for synteny studies of Lepidoptera and other insects was demonstrated.

Genetic mapping of Z chromosome and identification of W chromosome-specific markers in the silkworm, Bombyx mori

In the silkworm, Bombyx mori, the female is the heterogametic (ZW) sex and the male is homogametic (ZZ). The female heterogamety is a typical situation in the insect order Lepidoptera. Although the W chromosome in silkworm is strongly female determining, no W-linked gene for a morphological character has been found on it. The Z chromosome carries important traits of economic value as well as genes for various phenotypic traits, but only 2% of molecular information based on its relative size is known. Studies conducted so far indicate that the Z-linked genes are not dosage compensated. In the present study, we constructed a genetic map of randomly amplified polymorphic DNA fragments (RAPD), simple sequence repeats (SSR), and fluorescent intersimple sequence repeat PCR (FISSR) markers for the Z chromosome using a backcross mapping population. A total of 16 Z-linked markers were identified, characterized, and mapped using od, a recessive trait for translucent skin as an anchor marker yielding a total recombination map of 334.5 cM. The linkage distances obtained suggested that the markers were distributed throughout the Z chromosome. Four RAPD and four SSR markers that were linked to W chromosome were also identified. The proposed mapping approach should be useful to identify and map sex-linked traits in the silkworm. The economic and evolutionary significance of Z- and W-linked genes in silkworm, in particular, and lepidopterans, in general, is discussed.

A genetic linkage map of the mimetic butterfly Heliconius melpomene

Heliconius melpomene is a mimetic butterfly that exhibits great geographic variation in color pattern. We present here a genetic linkage map based on analysis of genetic markers in 73 individuals from a single F(2) family, offspring of a cross between H. m. cythera from western Ecuador and H. m. melpomene from French Guiana. A novel "three-step method" is described for the analysis of dominant markers in an F(2) cross, using outbred parental strains and taking advantage of the lack of crossing over in female Lepidoptera. This method is likely to prove useful for future mapping studies in outbred species with crossing over restricted to one sex, such as the Lepidoptera and Drosophila. The resulting linkage map has 21 linkage groups corresponding to the 21 chromosomes of H. melpomene and includes 219 AFLP markers, 23 microsatellites, 19 single-copy nuclear genes, and the color pattern switch genes Yb and Sb. The marker density is high, averaging >1/7 cM. The total map length is 1616 cM and the average chromosome length is 77 cM. The genome size of H. melpomene was estimated to be 292 Mb, giving a relationship of physical-to-map distance of 180 kb/cM. This map forms the basis for future comparative linkage analysis of color pattern evolution in Heliconius.

Genomic tools and cDNA derived markers for butterflies

The Lepidoptera have long been used as examples in the study of evolution, but some questions remain difficult to resolve due to a lack of molecular genetic data. However, as technology improves, genomic tools are becoming increasingly available to tackle unanswered evolutionary questions. Here we have used expressed sequence tags (ESTs) to develop genetic markers for two Müllerian mimic species, Heliconius melpomene and Heliconius erato. In total 1363 ESTs were generated, representing 330 gene objects in H. melpomene and 431 in H. erato. User-friendly bioinformatic tools were used to construct a nonredundant database of these putative genes (available at http://www.heliconius.org), and annotate them with blast similarity searches, InterPro matches and Gene Ontology terms. This database will be continually updated with EST sequences for the Papilionideae as they become publicly available, providing a tool for gene finding in the butterflies. Alignments of the Heliconius sequences with putative homologues derived from Bombyx mori or other public data sets were used to identify conserved PCR priming sites, and develop 55 markers that can be amplified from genomic DNA in both H. erato and H. melpomene. These markers will be used for comparative linkage mapping in Heliconius and will have applications in other phylogenetic and genomic studies in the Lepidoptera.

An AFLP-based interspecific linkage map of sympatric, hybridizing Colias butterflies

Colias eurytheme and C. philodice are sister species with broad sympatry in North America. They hybridize frequently and likely share a significant portion of their genomes through introgression. Both taxa have been ecologically well characterized and exploited to address a broad spectrum of evolutionary issues. Using AFLP markers, we constructed the first linkage map of Colias butterflies. The map is composed of 452 markers spanning 2541.7 cM distributed over 51 linkage groups (40 major groups and 11 small groups with 2-4 markers). Statistical tests indicate that these AFLP markers tend to cluster over the map, with the coefficient of variation of interval sizes being 1.236 (95% C.I. is 1.234-1.240). This nonrandom marker distribution can account for the nonequivalence between the number of linkage groups and the actual haploid chromosome number (N = 31). This study presents the initial step for further marker-assisted research on Colias butterflies, including QTL and introgression analyses. Further investigation of the genomes will help us understand better the roles of introgression and natural selection in the evolution of hybridizing species and devise more appropriate strategies to control these pests.