Molecular phylogeny of silkmoths reveals the origin of domesticated silkmoth, Bombyx mori from Chinese Bombyx mandarina and paternal inheritance of Antheraea proylei mitochondrial DNA

Molecular phylogeny of some of the economically important silkmoths was derived using three mitochondrial genes, 12S rRNA, 16S rRNA, and COI, and the control region (CR). Maximum likelihood (ML) analyses showed two distinct clades, one consisting of moths from Bombycidae family and the other from Saturniidae family. The mitochondrial CR showed length polymorphisms with indels. The ML analyses for complete mitochondrial genome sequences of Bombyx mori (strains Aojuku, C108, Backokjam, and Xiafang), Japanese and Chinese strains of B. mandarina (Japanese mandarina and Chinese mandarina) and, Antheraea pernyi revealed two distinct clades, one comprising of B. mori strains and the other with B. mandarina, and A. pernyi forming an outgroup. Pairwise distances revealed that all of the strains of B. mori studied are closer to Chinese than to Japanese mandarina. Phylogenetic analyses based on whole mitochondrial genome sequences, the finding of a tandem triplication of a 126bp repeat element only in Japanese mandarina, and chromosome number variation in B. mandarina suggest that B. mori must have shared its recent common ancestor with Chinese mandarina. Another wild species of the Bombycidae family, Theophila religiosa, whose phylogenetic status was not clear, clustered together with the other bombycid moths in the study. Analysis of the interspecific hybrid, A. proylei gave evidence for paternal inheritance of mitochondrial DNA.     

Mitochondrial genome nucleotide substitution pattern between domesticated silkmoth, Bombyx mori, and its wild ancestors, Chinese Bombyx mandarina and Japanese Bombyx mandarina

Bombyx mori and Bombyx mandarina are morphologically and physiologically similar. In this study, we compared the nucleotide variations in the complete mitochondrial (mt) genomes between the domesticated silkmoth, B. mori, and its wild ancestors, Chinese B. mandarina (ChBm) and Japanese B. mandarina (JaBm). The sequence divergence and transition mutation ratio between B. mori and ChBm are significantly smaller than those observed between B. mori and JaBm. The preference of transition by DNA strands between B. mori and ChBm is consistent with that between B. mori and JaBm, however, the regional variation in nucleotide substitution rate shows a different feature. These results suggest that the ChBm mt genome is not undergoing the same evolutionary process as JaBm, providing evidence for selection on mtDNA. Moreover, investigation of the nucleotide sequence divergence in the A+T-rich region of Bombyx mt genomes also provides evidence for the assumption that the A+T-rich region might not be the fastest evolving region of the mtDNA of insects.     

Nucleotide sequence variation in mitochondrial COI gene among 147 silkworm (Bombyx mori) strains from Japanese, Chinese, European and moltinism classes

We characterized the nucleotide sequences of PCR-amplified mitochondrial COI fragments of 147 silkworm (Bombyx mori) strains that have been maintained in the National Institute of Agrobiological Sciences. Coding sequences (714 bp) of the 147 COI fragments were classified into eight haplotypes based on nucleotide differences at eight segregating sites. No length variation was identified in this region. The 5'-noncoding region showed different features, wherein changes in the number of Ts in the T-stretch, together with two base substitutions, were observed. As a result, the 147 COI noncoding sequences were classified into six haplotypes. Combining the coding and noncoding regions, we identified 14 haplotypes. One of the 14 haplotypes, Hap1A was exclusively abundant in the Japanese native strain class, while this haplotype was less frequent in the other three native strain classes. This finding suggests that the Japanese strain class underwent significant genetic differentiation from the Chinese, European, and moltinism classes, when the each class is regarded as a population. Comparison of the nucleotide sequences to those of B. mandarina (which inhabits Japan) revealed changes that are significantly larger than those within either B. mori or B. mandarina. Furthermore, we detected no common haplotypes between them, which suggests the concept of suppressed gene flow between the two species.     

Homologues of fibroin L-chain and P25 of Bombyx mori are present in Dendrolimus spectabilis and Papilio xuthus but not detectable in Antheraea yamamai

Low molecular mass protein components of fibroin, whose electrophoretic patterns before and after the reductive cleavage of disulfide bonds were similar to those of L-chain and P25 of Bombyx mori, were identified in fibroin samples of Bombyx mandarina, Dendrolimus spectabilis and Papilio xuthus but not of Antheraea yamamai. Fibroin of A. yamamai is suggested to form a dimer of H-chain. Full length cDNA sequences were cloned for the homologues of L-chain and P25 from B. mandarina, D. spectabilis and P. xuthus. The deduced sequences of L-chain and P25 of B. mandarina are almost identical to those of B. mori, each containing a single amino acid change. Homologues of L-chain and P25 of D. spectabilis and P. xuthus show about 50% overall identity, respectively, with those of B. mori, but essential structural features; i.e. the three Cys residues in an L-chain and the eight Cys residues and one of the potential N-glycosylation sites in P25, are conserved in both species. These results, together with the published results for Galleria mellonella, suggest that the three-components (H-chain, L-chain and P25) complex of fibroin is rather common among Lepidopteran silk-producing insects, in contrast to the H-H dimer type found in the saturnid silkworm.     

野桑蚕、蓖麻蚕及家蚕基因组的RFLP分析

以家蚕Bombyx mori丝素重链基因、丝胶基因1和胰凝乳蛋白酶抑制因子13基因为探针,对野桑蚕B.mandarina、蓖麻蚕Philosamia cynthia ricini和家蚕B.mori基因组DNA进行限制性片段长度多态性分析。结果发现,在野桑蚕、蓖麻蚕基因组中存在着家蚕丝素重链基因、丝胶基因1的同源序列,而在中日野桑蚕以及蓖麻蚕品种间存在着限制性酶切位点差异;丝胶基因1在中国野桑蚕基因组的EcoRⅠ酶切图谱较日本野桑蚕与家蚕更为一致,表明家蚕与中国野桑蚕亲缘关系更近。此外,在野桑蚕基因组中发现了家蚕胰凝乳蛋白酶抑制因子13基因的同源序列,并且在家蚕品种间以及中日野桑蚕之间也存在着多态性。这些结果表明不同绢丝昆虫在适应生存环境的进化过程中,基因组发生了结构改变。     

Characterization of mitochondrial genome of Chinese wild mulberry silkworm, Bomyx mandarina (Lepidoptera: Bombycidae)

The complete mitochondrial genome of Chinese Bombyx mandarina(ChBm) was determined.The cir-cular genome is 15682 bp long, and contains a typical gene complement, order, and arrangement iden-tical to that of Bombyx mori(B.mori) and Japanese Bombyx mandarina(JaBm) except for two addi-tional tRNA-like structures:tRNASer(TGA)-like and tRNAIle(TAT)-like.All protein-coding sequences are initi-ated with a typical ATN codon except for the COI gene, which has a 4-bp TTAG putative initiator codon.Eleven of 13 protein-coding genes(PCGs) have a complete termination codon(all TAA), but the re-maining two genes terminate with incomplete codons.All tRNAs have the typical clover-leaf structures of mitochondrial tRNAs, with the exception of tRNASer(TGA)-like, with a four stem-and-loop structure.The length of the A T-rich region of ChBm is 484 bp, shorter than those of JaBm(747 bp) and B.mori(494―499 bp).Phylogenetic analysis among B.mori, ChBm, JaBm, and Antheraea pernyi(Anpe) showed that B.mori is more closely related to ChBm than JaBm.The earliest divergence time estimate for B.mori-ChBm and B.mori-JaBm is about 1.08±0.18―1.41±0.24 and 1.53±0.20―2.01±0.26 Mya, respec-tively.ChBm and JaBm diverged around 1.11±0.16―1.45±0.21 Mya.     

家蚕P450基因CYP18A1的克隆、序列分析及转录活性

蜕皮激素对昆虫生长、发育和繁殖有重要调控作用,尤其对蜕皮和变态过程。利用GenBank上登录的蜕皮激素C₂₆羟基化酶候选基因CYP18A1的氨基酸序列对家蚕Bombyx mori全基因组数据库进行BLASTP比对,发现了家蚕直向同源基因(ortholog),其完全编码序列经RT-PCR检测和克隆、测序验证后,再以此为信息探针检索家蚕表达序列标签(expressed sequence tags,EST)数据库进行拼接延伸,获得了一条包括5′非翻译区在内的长度为1 737 bp的cDNA序列,验证结果也表明与电子克隆序列完全一致(GenBank登录号为EF421988,P450命名委员会将其命名为CYP18A1)。该基因的开放阅读框为1 623 bp,编码541个氨基酸,含有包括P450s特征结构域在内的所有昆虫P450基因的5个保守结构域,其推定的分子量为61.67 kD,等电点为 8.54。将该基因cDNA序列与家蚕基因组序列进行比对,结果表明该基因具有6个外显子,5个内含子,外显子／内含子边界符合经典的GT-AG规则。同源性分析也发现家蚕CYP18A1与其他昆虫的直向同源基因具有较高相似性。用RT-PCR方法对家蚕主要发育变态时期与组织进行检测,显示出该基因的转录表达不仅具有时空特异性,而且在表达时期上与已报道的蚕体内蜕皮激素含量变化有紧密的一致性。该研究进一步证实了CYP18A1基因与昆虫体内蜕皮激素代谢平衡相关联。     

M chromosome of the wild silkworm, Bombyx mandarina (n = 27), corresponds to two chromosomes in the domesticated silkworm, Bombyx mori (n = 28).

Chromosomes of Bombyx mori (n = 28) and of Bombyx mandarina (n = 27) were studied cytogenetically to resolve the origin of the large M chromosome in the Japaneses type of B. mandarina. In the F1 progeny from the reciprocal cross between B. mandarina and B. mori, the mitotic chromosome number was 2n = 55, and a chromosome configuration of 26 bivalents plus 1 trivalent was observed at metaphase I of germ cells. The trivalent chromosome consisted of the M chromosome from B. mandarina and two chromosomes from B. mori. When males of B. mori were mated to the F1 females, nuclei with two types of chromosome number (2n = 55 and 2n = 56) and two sets of chromosome pairs (26 bivalents plus 1 trivalent versus 28 bivalents) were observed in the metaphase I stage. Linkage analysis showed that the 14th chromosome of B. mori was involved in these two types of chromosome segregation. This result indicates that the M chromosome in B. mandarina arose from a fusion between a chromosome corresponding to the 14th linkage group and another, yet unidentified linkage group.     

Genome Sequence of a Bombyx mori Nucleopolyhedrovirus Strain with Cubic Occlusion Bodies

Bombyx mori nucleopolyhedrovirus (BmNPV) is a typical species of Baculoviridae. The complete genome sequence of a BmNPV strain with cubic occlusion bodies is reported here. The genome of this strain consists of 127,465 nucleotides with a G+C content of 40.36% and is 97.3% and 97.5% identical to those of BmNPV strain T3 and Bombyx mandarina NPV S1, respectively. Despite the abnormal polyhedra it forms, the polyhedrin gene of the BmNPV cubic strain is 100% identical to those of the other two strains. Baculovirus repeated ORFs and homologous repeat regions cause the major differences in genome size of these BmNPV isolates.     

基于amy基因的中国野桑蚕遗传多样性及其与家蚕的系统发育关系

为探索中国野桑蚕Bombyx mandarina的遗传多样性及其与家蚕B. mori的系统发育关系, 采用PCR产物直接测序法（少数样本克隆测序）获得34个家蚕和野桑蚕样本淀粉酶基因amy序列片段（715 bp）。分析发现56个多态性位点, 鉴定出28种单倍型（haplotype）; 核苷酸多样性π=0.01390±0.00103, 单倍型多样度Hd=0.988±0.011。核苷酸不配对分析（mismatch analysis）和Fu’s Fs 检测表明中国野桑蚕曾发生过种群扩张。分子方差分析（AMOVA）表明, 遗野桑蚕传差异主要在种群内, 种群间和地理组群间差异不显著。聚类树上34个样本聚为3枝/3蔟, 野蚕和家蚕都不按地理区域或系统（类型）聚类, A枝由来自不同地区的野蚕和不同类型的家蚕混合构成, 并且进一步分成3个亚枝, 每一亚枝也同时包含家蚕和野蚕, B枝由3个家蚕和1个野蚕混合构成, C枝全部由来自不同地区的野蚕构成。网络分析没有发现“祖先单倍型”和优势单倍型。结果提示, 淀粉酶基因是一个多态性丰富的分子标记, 中国野桑蚕遗传多样性十分丰富, 据此推测家蚕起源于多种生态类型混杂的野桑蚕。     

Proteomic studies of isoforms of the P25 component of Bombyx mori fibroin

It is recognized that P25 is one of three polypeptide components of the fibroin synthesized in the larval silk gland (SG) of silkworm, having two glycosylated isoforms. In the present study, however, eight P25 isoforms were separated by proteomics, including two-dimensional gel electrophoresis of whole SG proteins, and were identified by the peptide mass fingerprinting method. Four of the eight isoforms were identified as Bombyx mandarina P25s, although the SG of Bombyx mori has never been considered to contain the P25 from B. mandarina. It is suggested that this diversity of P25 isoforms depends on phosphorylation modification in addition to glycosylation.     

Isolation and nucleotide sequence of cecropin B cDNA clones from the silkworm, Bombyx mori.

Two cDNA clones encoding cecropin B, an antibacterial protein, were isolated from a fat body cDNA library of the silkworm, Bombyx mori. Amino acid sequences of these clones, deduced from nucleotide sequences, were identical, including signal peptide regions. However, the nucleotide sequences were different at 30 positions. Deduced amino acid sequences of Bombyx mori cecropin B showed higher homology with cecropins from Lepidoptera than with those from Diptera.     

Diversity in copy number and structure of a silkworm morphogenetic gene as a result of domestication

The carotenoid-binding protein (CBP) of the domesticated silkworm, Bombyx mori, a major determinant of cocoon color, is likely to have been substantially influenced by domestication of this species. We analyzed the structure of the CBP gene in multiple strains of B. mori, in multiple individuals of the wild silkworm, B. mandarina (the putative wild ancestor of B. mori), and in a number of other lepidopterans. We found the CBP gene copy number in genomic DNA to vary widely among B. mori strains, ranging from 1 to 20. The copies of CBP are of several types, based on the presence of a retrotransposon or partial deletion of the coding sequence. In contrast to B. mori, B. mandarina was found to possess a single copy of CBP without the retrotransposon insertion, regardless of habitat. Several other lepidopterans were found to contain sequences homologous to CBP, revealing that this gene is evolutionarily conserved in the lepidopteran lineage. Thus, domestication can generate significant diversity of gene copy number and structure over a relatively short evolutionary time.     

Intrinsic bent DNA colocalizes with the sequence involved in the Nd-sD mutation in the Bombyx mori fibroin light chain gene

Multiple sequence alignments of the Bombyx mori fibroin light chain gene (fib-L) from hybrids and from Chinese and Japanese strains demonstrated that 51.6% of the fib-L third intron is conserved. One of these conserved segments, 41 bp long, contains the sequence CGTTATTATACATATT, which is duplicated in the B. mori Nd-s(D) mutant. In the present work, electrophoretic mobility assays and computational analyses revealed a major peak of intrinsic bent DNA within the segment that undergoes breakage in the previously-described Nd-s(D) mutation. This result suggested that this intrinsically-curved region might mediate DNA cleavage and enhance recombination events in the third intron of the Bombyx mori fib-L gene.     

家蚕核糖体蛋白L7基因cDNA序列的克隆和分析

为了研究家蚕孤雌生殖的调节机制,应用二维凝胶电泳(2DE)技术分离正常生殖的家蚕卵与孤雌生殖家蚕卵的差异蛋白质,在蛋白质水平上筛选与家蚕孤雌生殖过程相关的重要蛋白质.利用MALDI-TOF-TOF MS分析这些差异蛋白,获得了大量小肽的序列特征.BLAST搜索本实验室构建的cDNA文库,获得了1个与家蚕孤雌生殖相关的核糖体蛋白L7基因.根据已有的cDNA文库,采用RACE方法克隆得到该核糖体蛋白基因的全长cDNA.利用生物信息学的方法和工具,对这个基因在核酸水平和蛋白质水平分别作了详细的分析和讨论并进行蛋白结构预测.结果表明,核糖体L7基因的cDNA全长为858 bp,编码区包含6个外显子,共编码268个氨基酸残基,蛋白的疏水性平均值为-0.586,分子量大小为30 kD,极性的最大值为 39.616,最小值为0.451,等电点为10.52.分子系统分析显示,该蛋白与Apis, Lysiphlebus 和Meladema中的核糖体蛋白L7具有较高的同源性.     

家蚕羧酸酯酶基因克隆及差异表达

家蚕浓核病毒 (Bombyx mori densonucleosis virus,BmDNV)是蚕业生产上危害比较严重的一类病毒。用完全抗浓核病中国镇江株(BmDNV-Z)的家蚕品系秋丰、感性品系华八及以华八为轮回亲本回交8代和自交8代构建的近等基因系BC₈为材料,采用mRNA荧光差显技术首次分离克隆了家蚕羧酸酯酶（B. mori carboxylesterase,BmCarE）基因全长cDNA,并用实时荧光定量PCR检测了添毒后12 h、36 h、72 h BmCarE在感、抗BmDNV-Z家蚕品系中肠内的表达差异。结果表明： （1）添毒后12 h不同品系家蚕中肠BmCarE表达差异最大,抗性品系BC₈和秋丰分别是感性品系华八的17.714倍和3.602倍,三者彼此间的差异达到极显著水平;（2）同一品系添毒后12 h与添清水后12 h BmCarE表达也有较大差异, BC₈添毒是BC₈添清水的15.08倍, 秋丰添毒是秋丰添清水的3.39倍, 差异达到极显著水平,而华八添毒和添清水的BmCarE表达量均低,二者差异不显著;（3）同一品系添毒后不同时间BmCarE表达也有较大差异, BC₈和秋丰添毒后12 h BmCarE表达量最高,显著高于各自添毒后36 h和72 h表达水平,而添毒后36 h与72 h表达无显著差异;华八添毒后12 h、36 h和72 h,BmCarE表达无显著差异。上述结果提示羧酸酯酶基因可能与家蚕抗浓核病毒有一定关系。     

Structural features of B family chorion sequences in the silkmoth Bombyx mori, and their evolutionary implications.

Partial protein sequences, and DNA sequences of corresponding cDNA and genomic clones were obtained and analyzed to reveal the primary structural features of major, developmentally middle or late components of the B chorion multigene family in Bombyx mori. Comparisons with other types of sequences confirm and clarify the tripartite domain structure of chorion proteins. Glycine-, leucine- and tyrosine-containing, tandemly repetitive peptides form the bulk of the amino-terminal and carboxy-terminal domains ('arms'). Extensive sequence homologies suggest a common evolutionary origin for the amino-terminal arms of some B. mori B sequences and the corresponding portions of members of a different (A) chorion multigene family in Antheraea polyphemus, a distantly related silkmoth.     

Survey and analysis of microsatellites in the silkworm, Bombyx mori: frequency, distribution, mutations, marker potential and their conservation in heterologous species

We studied microsatellite frequency and distribution in 21.76-Mb random genomic sequences, 0.67-Mb BAC sequences from the Z chromosome, and 6.3-Mb EST sequences of Bombyx mori. We mined microsatellites of >/=15 bases of mononucleotide repeats and >/=5 repeat units of other classes of repeats. We estimated that microsatellites account for 0.31% of the genome of B. mori. Microsatellite tracts of A, AT, and ATT were the most abundant whereas their number drastically decreased as the length of the repeat motif increased. In general, tri- and hexanucleotide repeats were overrepresented in the transcribed sequences except TAA, GTA, and TGA, which were in excess in genomic sequences. The Z chromosome sequences contained shorter repeat types than the rest of the chromosomes in addition to a higher abundance of AT-rich repeats. Our results showed that base composition of the flanking sequence has an influence on the origin and evolution of microsatellites. Transitions/transversions were high in microsatellites of ESTs, whereas the genomic sequence had an equal number of substitutions and indels. The average heterozygosity value for 23 polymorphic microsatellite loci surveyed in 13 diverse silkmoth strains having 2-14 alleles was 0.54. Only 36 (18.2%) of 198 microsatellite loci were polymorphic between the two divergent silkworm populations and 10 (5%) loci revealed null alleles. The microsatellite map generated using these polymorphic markers resulted in 8 linkage groups. B. mori microsatellite loci were the most conserved in its immediate ancestor, B. mandarina, followed by the wild saturniid silkmoth, Antheraea assama.