Genomic sequence around butterfly wing development genes: annotation and comparative analysis

Background

Analysis of genomic sequence allows characterization of genome content and organization, and access beyond gene-coding regions for identification of functional elements. BAC libraries, where relatively large genomic regions are made readily available, are especially useful for species without a fully sequenced genome and can increase genomic coverage of phylogenetic and biological diversity. For example, no butterfly genome is yet available despite the unique genetic and biological properties of this group, such as diversified wing color patterns. The evolution and development of these patterns is being studied in a few target species, including Bicyclus anynana, where a whole-genome BAC library allows targeted access to large genomic regions.

Methodology/Principal Findings

We characterize ∼1.3 Mb of genomic sequence around 11 selected genes expressed in B. anynana developing wings. Extensive manual curation of in silico predictions, also making use of a large dataset of expressed genes for this species, identified repetitive elements and protein coding sequence, and highlighted an expansion of Alcohol dehydrogenase genes. Comparative analysis with orthologous regions of the lepidopteran reference genome allowed assessment of conservation of fine-scale synteny (with detection of new inversions and translocations) and of DNA sequence (with detection of high levels of conservation of non-coding regions around some, but not all, developmental genes).

Conclusions

The general properties and organization of the available B. anynana genomic sequence are similar to the lepidopteran reference, despite the more than 140 MY divergence. Our results lay the groundwork for further studies of new interesting findings in relation to both coding and non-coding sequence: 1) the Alcohol dehydrogenase expansion with higher similarity between the five tandemly-repeated B. anynana paralogs than with the corresponding B. mori orthologs, and 2) the high conservation of non-coding sequence around the genes wingless and Ecdysone receptor, both involved in multiple developmental processes including wing pattern formation.     

Mapping and recombination analysis of two moth colour mutations,Black moth and Wild wing spot,in the silkworm Bombyx mori

Many lepidopteran insects exhibit body colour variations, where the high phenotypic diversity observed in the wings and bodies of adults provides opportunities for studying adaptive morphological evolution. In the silkworm Bombyx mori, two genes responsible for moth colour mutation, Bm and Ws, have been mapped to 0.0 and 14.7 cM of the B. mori genetic linkage group 17; however, these genes have not been identified at the molecular level. We performed positional cloning of both genes to elucidate the molecular mechanisms that underlie the moth wing- and body-colour patterns in B. mori. We successfully narrowed down Bm and Ws to ~2-Mb-long and 100-kb-long regions on the same scaffold Bm_scaf33. Gene prediction analysis of this region identified 77 candidate genes in the Bm region, whereas there were no candidate genes in the Ws region. Fluorescence in-situ hybridisation analysis in Bm mutant detected chromosome inversion, which explains why there are no recombination in the corresponding region. The comparative genomic analysis demonstrated that the candidate regions of both genes shared synteny with a region associated with wing- and body-colour variations in other lepidopteran species including Biston betularia and Heliconius butterflies. These results suggest that the genes responsible for wing and body colour in B. mori may be associated with similar genes in other Lepidoptera.     

Construction of a bacterial artificial chromosome library of Endoclita excrescens as a tool for comparative gene mapping in Lepidoptera

Karyotype evolution in one of the most diverse and species‐rich group of insects, moths and butterflies (Lepidoptera), has interesting features that remain to be resolved. Recent studies showed that fluorescence in situ hybridization using bacterial artificial chromosome clones (BAC‐FISH) is an efficient cytogenetic method for identification and gene mapping of lepidopteran chromosomes. Using comparative mapping by BAC‐FISH, extensive synteny of genes was revealed between chromosomes of different lepidopteran species based on Bombyx mori genomic information. However, this comparative mapping has been done only in representatives of advanced groups of Lepidoptera. Here we constructed a BAC library of Endoclita excrescens, which belongs to the primitive lepidopteran family Hepialidae. High molecular weight DNA for the library construction was prepared from the pupae by using a rapid nuclear isolation method known in plants. The BAC clones of E. excrescens contain 66.6 kb inserts on average. The successful application of BAC‐FISH showed that the BAC library of E. excrescens is a useful tool for comparative gene mapping on chromosomes of this species.     

Genome-wide SNP identification for the construction of a high-resolution genetic map of Japanese flounder (Paralichthys olivaceus): applications to QTL mapping of Vibrio anguillarum disease resistance and comparative genomic analysis

High-resolution genetic maps are essential for fine mapping of complex traits, genome assembly, and comparative genomic analysis. Single-nucleotide polymorphisms (SNPs) are the primary molecular markers used for genetic map construction. In this study, we identified 13,362 SNPs evenly distributed across the Japanese flounder (Paralichthys olivaceus) genome. Of these SNPs, 12,712 high-confidence SNPs were subjected to high-throughput genotyping and assigned to 24 consensus linkage groups (LGs). The total length of the genetic linkage map was 3,497.29 cM with an average distance of 0.47 cM between loci, thereby representing the densest genetic map currently reported for Japanese flounder. Nine positive quantitative trait loci (QTLs) forming two main clusters for Vibrio anguillarum disease resistance were detected. All QTLs could explain 5.1–8.38% of the total phenotypic variation. Synteny analysis of the QTL regions on the genome assembly revealed 12 immune-related genes, among them 4 genes strongly associated with V. anguillarum disease resistance. In addition, 246 genome assembly scaffolds with an average size of 21.79 Mb were anchored onto the LGs; these scaffolds, comprising 522.99 Mb, represented 95.78% of assembled genomic sequences. The mapped assembly scaffolds in Japanese flounder were used for genome synteny analyses against zebrafish (Danio rerio) and medaka (Oryzias latipes). Flounder and medaka were found to possess almost one-to-one synteny, whereas flounder and zebrafish exhibited a multi-syntenic correspondence. The newly developed high-resolution genetic map, which will facilitate QTL mapping, scaffold assembly, and genome synteny analysis of Japanese flounder, marks a milestone in the ongoing genome project for this species.     

A rearrangement of the Z chromosome topology influences the sex-linked gene display in the European corn borer, <Emphasis Type="Italic">Ostrinia nubilalis</Emphasis>

Males are homogametic (ZZ) and females are heterogametic (WZ) with respect to the sex chromosomes in many species of butterflies and moths (insect order Lepidoptera). Genes on the Z chromosome influence traits involved in larval development, environmental adaptation, and reproductive isolation. To facilitate the investigation of these traits across Lepidoptera, we developed 43 degenerate primer pairs to PCR amplify orthologs of 43 Bombyx mori Z chromosome-linked genes. Of the 34 orthologs that amplified by PCR in Ostrinia nubilalis, 6 co-segregated with the Z chromosome anchor markers kettin (ket) and lactate dehydrogenase (ldh), and produced a consensus genetic linkage map of ~89 cM in combination with 5 AFLP markers. The O. nubilalis and B. mori Z chromosomes are comparatively co-linear, although potential gene inversions alter terminal gene orders and a translocation event disrupted synteny at one chromosome end. Compared to B. mori orthologs, O. nubilalis Z chromosome-linked genes showed conservation of tissue-specific and growth-stage-specific expression, although some genes exhibited species-specific expression across developmental stages or tissues. The O. nubilalis Z chromosome linkage map provides new tools for isolating quantitative trait loci (QTL) involved in sex-linked traits that drive speciation and it exposes genome rearrangements as a possible mechanism for differential gene regulation in Lepidoptera.     

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.     

Molecular Evolution of the Odorant and Gustatory Receptor Genes in Lepidopteran Insects: Implications for Their Adaptation and Speciation

Lepidoptera (comprised of butterflies and moths) is one of the largest groups of insects, including more than 160,000 described species. Chemoreception plays important roles in the adaptation of these species to a wide range of niches, e.g., plant hosts, egg-laying sites, and mates. This study investigated the molecular evolution of the lepidopteran odorant (Or) and gustatory receptor (Gr) genes using recently identified genes from Bombyx mori, Danaus plexippus, Heliconius melpomene, Plutella xylostella, Heliothis virescens, Manduca sexta, Cydia pomonella, and Spodoptera littoralis. A limited number of cases of large lineage-specific gene expansion are observed (except in the P. xylostella lineage), possibly due to selection against tandem gene duplication. There has been strong purifying selection during the evolution of both lepidopteran odorant and gustatory genes, as shown by the low ω values estimated through CodeML analysis, ranging from 0.0093 to 0.3926. However, purifying selection has been relaxed on some amino acid sites in these receptors, leading to sequence divergence, which is a precursor of positive selection on these sequences. Signatures of positive selection were detected only in a few loci from the lineage-specific analysis. Estimation of gene gains and losses suggests that the common ancestor of the Lepidoptera had fewer Or genes compared to extant species and an even more reduced number of Gr genes, particularly within the bitter receptor clade. Multiple gene gains and a few gene losses occurred during the evolution of Lepidoptera. Gene family expansion may be associated with the adaptation of lepidopteran species to plant hosts, especially after angiosperm radiation. Phylogenetic analysis of the moth sex pheromone receptor genes suggested that chromosomal translocations have occurred several times. New sex pheromone receptors have arisen through tandem gene duplication. Positive selection was detected at some amino acid sites predicted to be in the extracellular and transmembrane regions of the newly duplicated genes, which might be associated with the evolution of the new pheromone receptors.     

Hox genes are essential for the development of eyespots in Bicyclus anynana butterflies

The eyespot patterns found on the wings of nymphalid butterflies are novel traits that originated first in hindwings and subsequently in forewings, suggesting that eyespot development might be dependent on Hox genes. Hindwings differ from forewings in the expression of Ultrabithorax (Ubx), but the function of this Hox gene in eyespot development as well as that of another Hox gene Antennapedia (Antp), expressed specifically in eyespots centers on both wings, are still unclear. We used CRISPR-Cas9 to target both genes in Bicyclus anynana butterflies. We show that Antp is essential for eyespot development on the forewings and for the differentiation of white centers and larger eyespots on hindwings, whereas Ubx is essential not only for the development of at least some hindwing eyespots but also for repressing the size of other eyespots. Additionally, Antp is essential for the development of silver scales in male wings. In summary, Antp and Ubx, in addition to their conserved roles in modifying serially homologous segments along the anterior–posterior axis of insects, have acquired a novel role in promoting the development of a new set of serial homologs, the eyespot patterns, in both forewings (Antp) and hindwings (Antp and Ubx) of B. anynana butterflies. We propose that the peculiar pattern of eyespot origins on hindwings first, followed by forewings, could be due to an initial co-option of Ubx into eyespot development followed by a later, partially redundant, co-option of Antp into the same network.     

Samia cynthia versus Bombyx mori: comparative gene mapping between a species with a low-number karyotype and the model species of Lepidoptera

We performed gene-based comparative FISH mapping between a wild silkmoth, Samia cynthia ssp. with a low number of chromosomes (2n = 25-28) and the model species, Bombyx mori (2n = 56), in order to identify the genomic components that make up the chromosomes in a low-number karyotype. Mapping of 64 fosmid probes containing orthologs of B. mori genes revealed that the homologues of either two or four B. mori chromosomes constitute the S. c. ricini (Vietnam population, 2n = 27??/28??, Z0/ZZ) autosomes. Where tested, even the gene order was conserved between S. c. ricini and B. mori. This was also true for the originally autosomal parts of the neo-sex chromosomes in S. c. walkeri (Sapporo population, 2n = 26??/26??, neo-Wneo-Z/neo-Zneo-Z) and S. cynthia subsp. indet. (Nagano population, 2n = 25??/26??, neo-WZ₁Z₂/Z₁Z₁Z₂Z₂). The results are evidence for an internal stability of lepidopteran chromosomes even when all autosomes had undergone fusion processes to form a low-number karyotype.     

Functional endogenous viral elements in the genome of the parasitoid wasp Cotesia congregata: insights into the evolutionary dynamics of bracoviruses

Bracoviruses represent the most complex endogenous viral elements (EVEs) described to date. Nudiviral genes have been hosted within parasitoid wasp genomes since approximately 100 Ma. They play a crucial role in the wasp life cycle as they produce bracovirus particles, which are injected into parasitized lepidopteran hosts during wasp oviposition. Bracovirus particles encapsidate multiple dsDNA circles encoding virulence genes. Their expression in parasitized caterpillars is essential for wasp parasitism success. Here, we report on the genomic organization of the proviral segments (i.e. master sequences used to produce the encapsidated dsDNA circles) present in the Cotesia congregata parasitoid wasp genome. The provirus is composed of a macrolocus, comprising two-thirds of the proviral segments and of seven dispersed loci, each containing one to three segments. Comparative genomic analyses with closely related species gave insights into the evolutionary dynamics of bracovirus genomes. Conserved synteny in the different wasp genomes showed the orthology of the proviral macrolocus across different species. The nudiviral gene odv-e66-like1 is conserved within the macrolocus, suggesting an ancient co-localization of the nudiviral genome and bracovirus proviral segments. By contrast, the evolution of proviral segments within the macrolocus has involved a series of lineage-specific duplications.     

Rearing Temperature Influences Adult Response to Changes in Mating Status

Rearing environment can have an impact on adult behavior, but it is less clear how rearing environment influences adult behavior plasticity. Here we explore the effect of rearing temperature on adult mating behavior plasticity in the butterfly Bicyclus anynana, a species that has evolved two seasonal forms in response to seasonal changes in temperature. These seasonal forms differ in both morphology and behavior. Females are the choosy sex in cohorts reared at warm temperatures (WS butterflies), and males are the choosy sex in cohorts reared at cooler temperatures (DS butterflies). Rearing temperature also influences mating benefits and costs. In DS butterflies, mated females live longer than virgin females, and mated males live shorter than virgin males. No such benefits or costs to mating are present in WS butterflies. Given that choosiness and mating costs are rearing temperature dependent in B. anynana, we hypothesized that temperature may also impact male and female incentives to remate in the event that benefits and costs of second matings are similar to those of first matings. We first examined whether lifespan was affected by number of matings. We found that two matings did not significantly increase lifespan for either WS or DS butterflies relative to single matings. However, both sexes of WS but not DS butterflies experienced decreased longevity when mated to a non-virgin relative to a virgin. We next observed pairs of WS and DS butterflies and documented changes in mating behavior in response to changes in the mating status of their partner. WS but not DS butterflies changed their mating behavior in response to the mating status of their partner. These results suggest that rearing temperature influences adult mating behavior plasticity in B. anynana. This developmentally controlled behavioral plasticity may be adaptive, as lifespan depends on the partner’s mating status in one seasonal form, but not in the other.     

Eyespots deflect predator attack increasing fitness and promoting the evolution of phenotypic plasticity

Some eyespots are thought to deflect attack away from the vulnerable body, yet there is limited empirical evidence for this function and its adaptive advantage. Here, we demonstrate the conspicuous ventral hindwing eyespots found on Bicyclus anynana butterflies protect against invertebrate predators, specifically praying mantids. Wet season (WS) butterflies with larger, brighter eyespots were easier for mantids to detect, but more difficult to capture compared to dry season (DS) butterflies with small, dull eyespots. Mantids attacked the wing eyespots of WS butterflies more frequently resulting in greater butterfly survival and reproductive success. With a reciprocal eyespot transplant, we demonstrated the fitness benefits of eyespots were independent of butterfly behaviour. Regardless of whether the butterfly was WS or DS, large marginal eyespots pasted on the hindwings increased butterfly survival and successful oviposition during predation encounters. In previous studies, DS B. anynana experienced delayed detection by vertebrate predators, but both forms suffered low survival once detected. Our results suggest predator abundance, identity and phenology may all be important selective forces for B. anynana. Thus, reciprocal selection between invertebrate and vertebrate predators across seasons may contribute to the evolution of the B. anynana polyphenism.     

FISH identification of Helicoverpa armigera and Mamestra brassicae chromosomes by BAC and fosmid probes

Since the Bombyx mori genome sequence was published, conserved synteny between B. mori and some other lepidopteran species has been revealed by either FISH (fluorescence in situ hybridization) with BAC (bacterial artificial chromosome) probes or linkage analysis. However, no species belonging to the Noctuidae, the largest lepidopteran family which includes serious polyphagous pests, has been analyzed so far with respect to genome-wide conserved synteny and gene order. For that purpose, we selected the noctuid species Helicoverpa armigera and Mamestra brassicae, both with n = 31 chromosomes. Gene-defined fosmid clones from M. brassicae and BAC clones from a closely related species of H. armigera, Heliothis virescens, were used for a FISH analysis on pachytene chromosomes. We recognized all H. armigera chromosomes from specific cross-hybridization signals of 146 BAC probes. With 100 fosmid clones we identified and characterized all 31 bivalents of M. brassicae. Synteny and gene order were well conserved between the two noctuid species. The comparison with the model species B. mori (n = 28) showed the same phenomenon for 25 of the 28 chromosomes. Three chromosomes (#11, #23 and #24) had two counterparts each in H. armigera and M. brassicae. Since n = 31 is the modal chromosome number in Lepidoptera, the noctuid chromosomes probably represent an ancestral genome organization of Lepidoptera. This is the first identification of a full karyotype in Lepidoptera by means of BAC cross-hybridization between species. The technique shows the potential to expand the range of analyzed species efficiently.     

A homolog of the human Hermansky–Pudluck syndrome-5 (HPS5) gene is responsible for the oa larval translucent mutants in the silkworm, Bombyx mori

Normally, many granules containing uric acid accumulate in the larval integument of the silkworm, Bombyx mori. These uric acid granules cause the wild-type larval integument to be white or opaque, and the absence of these granules results in a translucent integument. Although about 30 B. mori loci governing larval translucency have been mapped, most have not been molecularly identified yet. Here, based on a structural analysis of a deletion of chromosome 14 that included the oa (aojyuku translucent) locus, we concluded that the BmHPS5 encoding a Bombyx homolog of the HPS5 subunit of biogenesis of lysosome-related organelles complex-2 is the candidate for the oa locus. Nucleotide sequence analyses of cDNAs and genomic DNAs in three mutant strains, each of which were homozygous for the respective allele of the oa locus (oa, oa ² , and oa ^v ), revealed that each mutant strain has a frame shift or a premature stop codon (caused by deletion or nonsense mutation, respectively) in the BmHPS5 gene. Our findings indicate that some genes that cause the translucent phenotype in Bombyx, some HPS-associated genes in humans, and some genes that cause mutant eye color phenotypes in Drosophila are homologous and participate in an evolutionarily conserved mechanism that leads to biogenesis of lysosome-related organelles.     

Characterization of an Ancient Lepidopteran Lateral Gene Transfer

Bacteria to eukaryote lateral gene transfers (LGT) are an important potential source of material for the evolution of novel genetic traits. The explosion in the number of newly sequenced genomes provides opportunities to identify and characterize examples of these lateral gene transfer events, and to assess their role in the evolution of new genes. In this paper, we describe an ancient lepidopteran LGT of a glycosyl hydrolase family 31 gene (GH31) from an Enterococcus bacteria. PCR amplification between the LGT and a flanking insect gene confirmed that the GH31 was integrated into the Bombyx mori genome and was not a result of an assembly error. Database searches in combination with degenerate PCR on a panel of 7 lepidopteran families confirmed that the GH31 LGT event occurred deep within the Order approximately 65–145 million years ago. The most basal species in which the LGT was found is Plutella xylostella (superfamily: Yponomeutoidea). Array data from Bombyx mori shows that GH31 is expressed, and low dN/dS ratios indicates the LGT coding sequence is under strong stabilizing selection. These findings provide further support for the proposition that bacterial LGTs are relatively common in insects and likely to be an underappreciated source of adaptive genetic material.     

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.     

The Lepidoptera Odorant Binding Protein gene family: Gene gain and loss within the GOBP/PBP complex of moths and butterflies

Butterflies and moths differ significantly in their daily activities: butterflies are diurnal while moths are largely nocturnal or crepuscular. This life history difference is presumably reflected in their sensory biology, and especially the balance between the use of chemical versus visual signals. Odorant Binding Proteins (OBP) are a class of insect proteins, at least some of which are thought to orchestrate the transfer of odor molecules within an olfactory sensillum (olfactory organ), between the air and odor receptor proteins (ORs) on the olfactory neurons. A Lepidoptera specific subclass of OBPs are the GOBPs and PBPs; these were the first OBPs studied and have well documented associations with olfactory sensilla. We have used the available genomes of two moths, Manduca sexta and Bombyx mori, and two butterflies, Danaus plexippus and Heliconius melpomene, to characterize the GOBP/PBP genes, attempting to identify gene orthologs and document specific gene gain and loss. First, we identified the full repertoire of OBPs in the M. sexta genome, and compared these with the full repertoire of OBPs from the other three lepidopteran genomes, the OBPs of Drosophila melanogaster and select OBPs from other Lepidoptera. We also evaluated the tissue specific expression of the M. sexta OBPs using an available RNAseq databases. In the four lepidopteran species, GOBP2 and all PBPs reside in single gene clusters; in two species GOBP1 is documented to be nearby, about 100 kb from the cluster; all GOBP/PBP genes share a common gene structure indicating a common origin. As such, the GOBP/PBP genes form a gene complex. Our findings suggest that (1) the lepidopteran GOBP/PBP complex is a monophyletic lineage with origins deep within Lepidoptera phylogeny, (2) within this lineage PBP gene evolution is much more dynamic than GOBP gene evolution, and (3) butterflies may have lost a PBP gene that plays an important role in moth pheromone detection, correlating with a shift from olfactory (moth) to visual (butterfly) communication, at least regarding long distance mate recognition. These findings will be clarified by additional lepidopteran genomic data, but the observation that moths and butterflies share most of the PBP/GOBP genes suggests that they also share common chemosensory-based behavioral pathways.     

The first gene-based map of Lupinus angustifolius L.-location of domestication genes and conserved synteny with Medicago truncatula

We report the first gene-based linkage map of Lupinus angustifolius (narrow-leafed lupin) and its comparison to the partially sequenced genome of Medicago truncatula. The map comprises 382 loci in 20 major linkage groups, two triplets, three pairs and 11 unlinked loci and is 1,846 cM in length. The map was generated from the segregation of 163 RFLP markers, 135 gene-based PCR markers, 75 AFLP and 4 AFLP-derived SCAR markers in a mapping population of 93 recombinant inbred lines, derived from a cross between domesticated and wild-type parents. This enabled the mapping of five major genes controlling key domestication traits in L. angustifolius. Using marker sequence data, the L. angustifolius genetic map was compared to the partially completed M. truncatula genome sequence. We found evidence of conserved synteny in some regions of the genome despite the wide evolutionary distance between these legume species. We also found new evidence of widespread duplication within the L. angustifolius genome.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorised users.     

The genome sequence of Samia ricini,a new model species of lepidopteran insect

Samia ricini, a gigantic saturniid moth, has the potential to be a novel lepidopteran model species. Samia ricini is far more resistant to diseases than the current model species Bombyx mori, and therefore can be more easily reared. In addition, genetic resources available for S. ricini rival those for B. mori: at least 26 ecoraces of S. ricini are reported and S. ricini can hybridize with wild Samia species, which are distributed throughout Asian countries, and produce fertile progenies. Physiological traits such as food preference, integument colour and larval spot pattern differ among S. ricini strains and wild Samia species so that those traits can be targeted in forward genetic analyses. To facilitate genetic research in S. ricini, we determined its whole genome sequence. The assembled genome of S. ricini was 458 Mb with 155 scaffolds, and the scaffold N50 length of the assembly was ~ 21 Mb. In total, 16,702 protein coding genes were predicted. While the S. ricini genome was mostly collinear with that of B. mori with some rearrangements and few S. ricini‐specific genes were discovered, chorion genes and fibroin genes seemed to have expanded in the S. ricini lineage. As the first step of genetic analyses, causal genes for “Blue,” “Yellow,” “Spot,” and “Red cocoon” phenotypes were mapped to chromosomes.