Qingyou Xia

Complete Resequencing of 40 Genomes Reveals Domestication Events and Genes in Silkworm (Bombyx)

The Taming of the Silkworm Silkworms, Bombyx mori, represent one of the few domesticated insects, having been domesticated over 10,000 years ago. Xia et al. (p. 433, published online 27 August) sequenced 29 domestic and 11 wild silkworm lines and identified genes that were most likely to be selected during domestication. These genes represent those that enhance silk production, reproduction, and growth. Furthermore, silkworms were probably only domesticated once from a large progenitor population, rather than on multiple occasions, as has been observed for other domesticated animals. Silkworm genomes show signatures of selection associated with domestication. A single–base pair resolution silkworm genetic variation map was constructed from 40 domesticated and wild silkworms, each sequenced to approximately threefold coverage, representing 99.88% of the genome. We identified ~16 million single-nucleotide polymorphisms, many indels, and structural variations. We find that the domesticated silkworms are clearly genetically differentiated from the wild ones, but they have maintained large levels of genetic variability, suggesting a short domestication event involving a large number of individuals. We also identified signals of selection at 354 candidate genes that may have been important during domestication, some of which have enriched expression in the silk gland, midgut, and testis. These data add to our understanding of the domestication processes and may have applications in devising pest control strategies and advancing the use of silkworms as efficient bioreactors.

Single base-resolution methylome of the silkworm reveals a sparse epigenomic map

Epigenetic regulation in insects may have effects on diverse biological processes. Here we survey the methylome of a model insect, the silkworm Bombyx mori, at single-base resolution using Illumina high-throughput bisulfite sequencing (MethylC-Seq). We conservatively estimate that 0.11% of genomic cytosines are methylcytosines, all of which probably occur in CG dinucleotides. CG methylation is substantially enriched in gene bodies and is positively correlated with gene expression levels, suggesting it has a positive role in gene transcription. We find that transposable elements, promoters and ribosomal DNAs are hypomethylated, but in contrast, genomic loci matching small RNAs in gene bodies are densely methylated. This work contributes to our understanding of epigenetics in insects, and in contrast to previous studies of the highly methylated genomes of Arabidopsis and human, demonstrates a strategy for sequencing the epigenomes of organisms such as insects that have low levels of methylation.

KAIKObase: An integrated silkworm genome database and data mining tool

BackgroundThe silkworm, Bombyx mori, is one of the most economically important insects in many developing countries owing to its large-scale cultivation for silk production. With the development of genomic and biotechnological tools, B. mori has also become an important bioreactor for production of various recombinant proteins of biomedical interest. In 2004, two genome sequencing projects for B. mori were reported independently by Chinese and Japanese teams; however, the datasets were insufficient for building long genomic scaffolds which are essential for unambiguous annotation of the genome. Now, both the datasets have been merged and assembled through a joint collaboration between the two groups.DescriptionIntegration of the two data sets of silkworm whole-genome-shotgun sequencing by the Japanese and Chinese groups together with newly obtained fosmid- and BAC-end sequences produced the best continuity (~3.7 Mb in N50 scaffold size) among the sequenced insect genomes and provided a high degree of nucleotide coverage (88%) of all 28 chromosomes. In addition, a physical map of BAC contigs constructed by fingerprinting BAC clones and a SNP linkage map constructed using BAC-end sequences were available. In parallel, proteomic data from two-dimensional polyacrylamide gel electrophoresis in various tissues and developmental stages were compiled into a silkworm proteome database. Finally, a Bombyx trap database was constructed for documenting insertion positions and expression data of transposon insertion lines.ConclusionFor efficient usage of genome information for functional studies, genomic sequences, physical and genetic map information and EST data were compiled into KAIKObase, an integrated silkworm genome database which consists of 4 map viewers, a gene viewer, and sequence, keyword and position search systems to display results and data at the level of nucleotide sequence, gene, scaffold and chromosome. Integration of the silkworm proteome database and the Bombyx trap database with KAIKObase led to a high-grade, user-friendly, and comprehensive silkworm genome database which is now available from URL: http://sgp.dna.affrc.go.jp/KAIKObase/.

Large Scale Full-Length cDNA Sequencing Reveals a Unique Genomic Landscape in a Lepidopteran Model Insect, Bombyx mori

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

Sequence structure and expression pattern of a novel anionic defensin-like gene from silkworm (Bombyx mori).

A defensin-like gene, BmdefA, was rediscovered in the silkworm genome and expressed sequence tags databases. The open reading frame of BmdefA encodes a prepropeptide consisting of a 22-residue signal peptide, a 34-residue propeptide, and a 36-residue mature peptide with a molecular mass of 4.0xa0kDa. The mature peptide possesses the characteristic six-cysteine motif of insect defensins, and its predicted isoelectric point is 4.12, indicating it is a novel anionic defensin. An intron is present in BmdefA and several cis-regulatory elements are in the regulating region. It is transcribed constitutively at a high level in the hemocyte, silk gland, head, and ovary of the silkworm larvae, and in the fat body of early-stage pupae and moth. BmdefA is also strongly induced by immune challenge. These results suggest that BmdefA plays an important role in both immunity and metamorphosis.

Genetic diversity, molecular phylogeny and selection evidence of the silkworm mitochondria implicated by complete resequencing of 41 genomes

BackgroundMitochondria are a valuable resource for studying the evolutionary process and deducing phylogeny. A few mitochondria genomes have been sequenced, but a comprehensive picture of the domestication event for silkworm mitochondria remains to be established. In this study, we integrate the extant data, and perform a whole genome resequencing of Japanese wild silkworm to obtain breakthrough results in silkworm mitochondrial (mt) population, and finally use these to deduce a more comprehensive phylogeny of the Bombycidae.ResultsWe identified 347 single nucleotide polymorphisms (SNPs) in the mt genome, but found no past recombination event to have occurred in the silkworm progenitor. A phylogeny inferred from these whole genome SNPs resulted in a well-classified tree, confirming that the domesticated silkworm, Bombyx mori, most recently diverged from the Chinese wild silkworm, rather than from the Japanese wild silkworm. We showed that the population sizes of the domesticated and Chinese wild silkworms both experience neither expansion nor contraction. We also discovered that one mt gene, named cytochrome b, shows a strong signal of positive selection in the domesticated clade. This gene is related to energy metabolism, and may have played an important role during silkworm domestication.ConclusionsWe present a comparative analysis on 41 mt genomes of B. mori and B. mandarina from China and Japan. With these, we obtain a much clearer picture of the evolution history of the silkworm. The data and analyses presented here aid our understanding of the silkworm in general, and provide a crucial insight into silkworm phylogeny.

Cloning and evolutionary analysis of tobacco MAPK gene family

The mitogen-activated protein (MAP) kinase cascade is an important signaling module which is involved in biotic and abiotic stress responses as well as plant growth and development. In this study, we identified 17 tobacco MAPKs including 11 novel tobacco MAPK genes that have not been identified before. Comparative analysis with MAPK gene families from other plants, such as Athaliana thaliana, rice and poplar, suggested that tobacco MAPKs (such as NtMPK1, NtMPK3 and NtMPK8) might play similar functions in response to abiotic and biotic stresses. QRT-PCR analysis revealed that a total of 14 NtMPKs were regulated by SA and/or MeJA, suggesting their potential roles involved in plant defense response. In addition, 6 NtMPKs were induced by drought treatment, implying their roles in response to drought stress. Our results indicated that most of tobacco MAPK might be involved in plant defense response, which provides the basis for further analysis on physiological functions of tobacco MAPKs.

Detection and analysis of alternative splicing in the silkworm by aligning expressed sequence tags with the genomic sequence

We identified 277 alternative splice forms in silkworm genes based on aligning expressed sequence tags with genomic sequences, using a transcipt assembly program. A large fraction (74%) of these alternative splices are located in protein‐coding regions and alter protein products, whereas only 26% are in untranslated regions. From the alternative splices located in protein‐coding regions, some (43%) affect protein domains that bind various biological molecules. The vast majority of the detected alternative forms in this study appear to be novel, and potentially affect biologically meaningful control of function in silkworm genes. Our results indicate that alternative splicing in silkworm largely produces protein diversity and functional diversity, and is a widely used mechanism for regulating gene expression.

BmHrp28 is a RNA‐binding protein that binds to the female‐specific exon 4 of Bombyx mori dsx pre‐mRNA

The Bombyx mori sex determination gene Bmdsx is alternatively spliced in the male and female to produce the male‐ and female‐specific proteins. In an effort to better understand the mechanism of the alternative splicing regulation of Bmdsx, we conducted a gel‐shift assay followed by LC‐MS/MS analysis to identify the putative proteins bound to the cis‐element CE1+6 in the exon 4 of Bmdsx. A protein named as BmHrp28 which is homologous to the Drosophila Hrp48, a member of the hnRNPA/B family, was identified and expressed in Escherichia coli for testing RNA‐protein binding in vitro. All of the results showed that BmHrp28 specifically bound to the CE1+6 RNA probe. BmHrp28 has two RNA recognition motifs at the N‐terminal and a glycine‐rich motif at the C‐terminal. It might be one of the factors involved in the male‐specific splicing of Bmdsx.