We present a de novo assembly of a eukaryote transcriptome using 454 pyrosequencing data. The Glanville fritillary butterfly (Melitaea cinxia; Lepidoptera: Nymphalidae) is a prominent species in population biology but had no previous genomic data. Sequencing runs using two normalized complementary DNA collections from a genetically diverse pool of larvae, pupae, and adults yielded 608 053 expressed sequence tags (mean length = 110 nucleotides), which assembled into 48 354 contigs (sets of overlapping DNA segments) and 59 943 singletons. BLAST comparisons confirmed the accuracy of the sequencing and assembly, and indicated the presence of c. 9000 unique genes, along with > 6000 additional microarray-confirmed unannotated contigs. Average depth of coverage was 6.5-fold for the longest 4800 contigs (348–2849 bp in length), sufficient for detecting large numbers of single nucleotide polymorphisms. Oligonucleotide microarray probes designed from the assembled sequences showed highly repeatable hybridization intensity and revealed biological differences among individuals. We conclude that 454 sequencing, when performed to provide sufficient coverage depth, allows de novo transcriptome assembly and a fast, cost-effective, and reliable method for development of functional genomic tools for nonmodel species. This development narrows the gap between approaches based on model organisms with rich genetic resources vs. species that are most tractable for ecological and evolutionary studies.

/pdf/rapid-transcriptome-characterization-for-a-nonmodel-organism-2gy5fyc4ip.pdf

Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing.

Although all bilaterian animals have a related set of Hox genes, the genomic organization of this gene complement comes in different flavors In some unrelated species, Hox genes are clustered; in others, they are not This indicates that the bilaterian ancestor had a clustered Hox gene family and that, subsequently, this genomic organization was either maintained or lost Remarkably, the tightest organization is found in vertebrates, raising the embarrassingly finalistic possibility that vertebrates have maintained best this ancestral configuration Alternatively, could they have co-evolved with an increased ;organization' of the Hox clusters, possibly linked to their genomic amplification, which would be at odds with our current perception of evolutionary mechanisms? When discussing the why's and how's of Hox gene clustering, we need to account for three points: the mechanisms of cluster evolution; the underlying biological constraints; and the developmental modes of the animals under consideration By integrating these parameters, general conclusions emerge that can help solve the aforementioned dilemma

The rise and fall of Hox gene clusters

Once called the 'Rosetta stone' of developmental biology, the homeobox continues to fascinate both evolutionary and developmental biologists. The birth of the homeotic, or Hox, gene cluster, and its subsequent evolution, has been crucial in mediating the major transitions in metazoan body plan. Comparative genomics studies indicate that the more recently discovered ParaHox and NK clusters were linked to the Hox cluster early in evolution, and that together they constituted a 'megacluster' of homeobox genes that conspicuously contributed to body-plan evolution.

/pdf/the-genesis-and-evolution-of-homeobox-gene-clusters-si6od80vq1.pdf

The genesis and evolution of homeobox gene clusters.

We review progress in applying molecular genetic and genomic technologies to studies in the domesticated silkworm, Bombyx mori, highlighting its use as a model for Lepidoptera, and in sericulture and biotechnology. Dense molecular linkage maps are being integrated with classical linkage maps for positional cloning and marker-assisted selection. Classical mutations have been identified by a candidate gene approach. Cytogenetic and sequence analyses show that the W chromosome is composed largely of nested full-length long terminal repeat retrotransposons. Z-chromosome-linked sequences show a lack of dosage compensation. The downstream sex differentiation mechanism has been studied via the silkworm homolog of doublesex. Expressed sequence tagged databases have been used to discover Lepidoptera-specific genes, provide evidence for horizontal gene transfer, and construct microarrays. Physical maps using large-fragment bacterial artificial chromosome libraries have been constructed, and whole-genome shotgun sequencing is underway. Germline transformation and transient expression systems are well established and available for functional studies, high-level protein expression, and gene silencing via RNA interference.

/pdf/the-genetics-and-genomics-of-the-silkworm-bombyx-mori-34s0xlx3gu.pdf

The genetics and genomics of the silkworm, bombyx mori*

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

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

We examined chromosomal distribution of major ribosomal DNAs (rDNAs), clustered in the nucleolar organizer regions (NORs), in 18 species of moths and butterflies using fluorescence in situ hybridization with a codling moth (Cydia pomonella) 18S rDNA probe. Most species showed one or two rDNA clusters in their haploid karyotype but exceptions with 4-11 clusters also occurred. Our results in a compilation with previous data revealed dynamic evolution of rDNA distribution in Lepidoptera except Noctuoidea, which showed a highly uniform rDNA pattern. In karyotypes with one NOR, interstitial location of rDNA prevailed, whereas two-NOR karyotypes showed mostly terminally located rDNA clusters. A possible origin of the single interstitial NOR by fusion between two NOR-chromosomes with terminal rDNA clusters lacks support in available data. In some species, spreading of rDNA to new, mostly terminal chromosome regions was found. The multiplication of rDNA clusters without alteration of chromosome numbers rules out chromosome fissions as a major mechanism of rDNA expansion. Based on rDNA dynamics in Lepidoptera and considering the role of ordered nuclear architecture in karyotype evolution, we propose ectopic recombination, i.e., homologous recombination between repetitive sequences of non-homologous chromosomes, as a primary motive force in rDNA repatterning.

/pdf/evolutionary-dynamics-of-rdna-clusters-on-chromosomes-of-1g66to2q9d.pdf

Evolutionary dynamics of rDNA clusters on chromosomes of moths and butterflies (Lepidoptera)

Lepidoptera, i.e. moths and butterflies, have a female heterogametic sex chromosome system, with most females having a WZ constitution while males are ZZ. Besides this predominant WZ/ZZ system, Z/ZZ, WZ1Z2/Z1Z1Z2Z2 and W1W2Z/ZZ systems also occur. Sex is determined by an unknown W-linked gene or genes in Bombyx mori, but by dosage-dependent and equally unknown Z-linked genes in all Z/ZZ species. The female heterogametic sex chromosome system has been conserved for at least 180 MY in the phylogenetic branch that combines Lepidoptera and Trichoptera. The W chromosome, which is present in most lepidopteran species, was incorporated in the sex chromosome system much later, about 90–100 MY ago. The Z chromosomes are highly conserved among Lepidoptera, much like the Z in birds or the X in mammals. The W, on the other hand, is evolving rapidly. It is crammed with repetitive elements which appear to have a high turnover rate but poor in or even devoid of protein-coding genes. It has frequently undergone fusion with autosomes or sporadically lost altogether.

/pdf/sex-chromosome-evolution-in-moths-and-butterflies-1ftz2g1wpi.pdf

Sex chromosome evolution in moths and butterflies

Lepidopteran species have a relatively high number of small holocentric chromosomes (Bombyx mori, 2n = 56). Chromosome identification has long been hampered in this group by the high number and by the absence of suitable markers like centromere position and chromosome bands. In this study, we carried out fluorescence in situ hybridization (FISH) on meiotic chromosome complements using genetically mapped B. mori bacterial artificial chromosomes (BACs) as probes. The combination of two to four either green or red fluorescence-labeled probes per chromosome allowed us to recognize unequivocally each of the 28 bivalents of the B. mori karyotype by its labeling pattern. Each chromosome was assigned one of the already established genetic linkage groups and the correct orientation in the chromosome was defined. This facilitates physical mapping of any other sequence and bears relevance for the ongoing B. mori genome projects. Two-color BAC-FISH karyotyping overcomes the problem of chromosome recognition in organisms where conventional banding techniques are not available.

https://www.genetics.org/content/genetics/170/2/675.full.pdf

The Bombyx mori karyotype and the assignment of linkage groups.

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.

https://www.genetics.org/content/genetics/173/3/1319.full.pdf

A Second-Generation Integrated Map of the Silkworm Reveals Synteny and Conserved Gene Order Between Lepidopteran Insects

We isolated four W chromosome-derived bacterial artificial chromosome (W-BAC) clones from Bombyx mori BAC libraries by the polymerase chain reaction and used them as probes for fluorescence in situ hybridization (FISH) on chromosome preparations from B. mori females. All four W-BAC probes surprisingly highlighted the whole wild-type W sex chromosome and also identified the entire original W-chromosomal region in W chromosome-autosome translocation mutants. This is the first successful identification of a single chromosome by means of BAC-FISH in species with holokinetic chromosomes. Genomic in situ hybridization (GISH) by using female-derived genomic probes highlighted the W chromosome in a similar chromosome-painting manner. Besides the W, hybridization signals of W-BAC probes also occurred in telomeric and/or subtelomeric regions of the autosomes. These signals coincided well with those of female genomic probes except one additional GISH signal that was observed in a large heterochromatin block of one autosome pair. Our results support the opinion that the B. mori W chromosome accumulated transposable elements and other repetitive sequences that also occur, but scattered, elsewhere in the respective genome.

/pdf/w-derived-bac-probes-as-a-new-tool-for-identification-of-the-3kja4f63fo.pdf

Atsuo Yoshido

Papers

Evolutionary dynamics of rDNA clusters on chromosomes of moths and butterflies (Lepidoptera)

Sex chromosome evolution in moths and butterflies

The Bombyx mori karyotype and the assignment of linkage groups.

A Second-Generation Integrated Map of the Silkworm Reveals Synteny and Conserved Gene Order Between Lepidopteran Insects

W-derived BAC probes as a new tool for identification of the W chromosome and its aberrations in Bombyx mori.