Huizhen Guo

Expression map of a complete set of gustatory receptor genes in chemosensory organs of Bombyx mori

Most lepidopteran species are herbivores, and interaction with host plants affects their gene expression and behavior as well as their genome evolution. Gustatory receptors (Grs) are expected to mediate host plant selection, feeding, oviposition and courtship behavior. However, due to their high diversity, sequence divergence and extremely low level of expression it has been difficult to identify precisely a complete set of Grs in Lepidoptera. By manual annotation and BAC sequencing, we improved annotation of 43 gene sequences compared with previously reported Grs in the most studied lepidopteran model, the silkworm, Bombyx mori, and identified 7 new tandem copies of BmGr30 on chromosome 7, bringing the total number of BmGrs to 76. Among these, we mapped 68 genes to chromosomes in a newly constructed chromosome distribution map and 8 genes to scaffolds; we also found new evidence for large clusters of BmGrs, especially from the bitter receptor family. RNA-seq analysis of diverse BmGr expression patterns in chemosensory organs of larvae and adults enabled us to draw a precise organ specific map of BmGr expression. Interestingly, most of the clustered genes were expressed in the same tissues and more than half of the genes were expressed in larval maxillae, larval thoracic legs and adult legs. For example, BmGr63 showed high expression levels in all organs in both larval and adult stages. By contrast, some genes showed expression limited to specific developmental stages or organs and tissues. BmGr19 was highly expressed in larval chemosensory organs (especially antennae and thoracic legs), the single exon genes BmGr53 and BmGr67 were expressed exclusively in larval tissues, the BmGr27-BmGr31 gene cluster on chr7 displayed a high expression level limited to adult legs and the candidate CO2 receptor BmGr2 was highly expressed in adult antennae, where few other Grs were expressed. Transcriptional analysis of the Grs in B. mori provides a valuable new reference for finding genes involved in plant-insect interactions in Lepidoptera and establishing correlations between these genes and vital insect behaviors like host plant selection and courtship for mating.

Transcriptome analysis of interactions between silkworm and cytoplasmic polyhedrosis virus.

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) specifically infects silkworm midgut (MG) and multiplication occurs mainly in posterior midgut (PM). In this study, MG and fat body (FB) were extracted at 0, 3, 24, and 72 h after BmCPV infection. The total sequence reads of each sample were more than 1510000, and the mapping ratio exceeded 95.3%. Upregulated transcripts increased in MG during the infection process. Gene ontology (GO) categories showed that antioxidants were all upregulated in FB but not in MG. BGI001299, BGI014434, BGI012068, and BGI009201 were MG-specific genes with transmembrane transport function, the expression of which were induced by BmCPV. BGI001299, BGI014434, and BGI012068 expressed in entire MG and may be involved in BmCPV invasion. BGI009201 expressed only in PM and may be necessary for BmCPV proliferation. BmPGRP-S2 and BGI012452 (a putative serine protease) were induced by BmCPV and may be involved in immune defense against BmCPV. The expression level of BmCPV S1, S2, S3, S6, and S7 was high and there was no expression of S9 in MG 72 h, implying that the expression time of structural protein coding genes is earlier. These results provide insights into the mechanism of BmCPV infection and host defense.

A Comprehensive Analysis of the Chorion Locus in Silkmoth

Despite more than 40 years of intense study, essential features of the silkmoth chorion (eggshell) are still not fully understood. To determine the precise structure of the chorion locus, we performed extensive EST analysis, constructed a bacterial artificial chromosome (BAC) contig, and obtained a continuous genomic sequence of 871,711 base pairs. We annotated 127 chorion genes in two segments interrupted by a 164 kb region with 5 non-chorion genes, orthologs of which were on chorion bearing scaffolds in 4 ditrysian families. Detailed transcriptome analysis revealed expression throughout choriogenesis of most chorion genes originally categorized as “middle”, and evidence for diverse regulatory mechanisms including cis-elements, alternative splicing and promoter utilization, and antisense RNA. Phylogenetic analysis revealed multigene family associations and faster evolution of early chorion genes and transcriptionally active pseudogenes. Proteomics analysis identified 99 chorion proteins in the eggshell and micropyle localization of 1 early and 6 Hc chorion proteins.

Selection of reference genes for analysis of stress-responsive genes after challenge with viruses and temperature changes in the silkworm Bombyx mori

Viruses and high temperature (HT) are the primary threats to silkworms. Changes in the expression of stress-response genes can be measured using quantitative polymerase chain reaction (qPCR) after exposure to viruses or HT. However, appropriate reference genes (RGs) for qPCR data normalization have not been established in this organism. In this study, we summarized the RGs used in the previous silkworm studies after infection with Bombyx mori nucleopolyhedrovirus (BmNPV), B. mori cytoplasmic polyhedrosis virus (BmCPV), or B. mori densovirus (BmDNV) or after HT treatment. The expression levels of these RGs were extracted from silkworm transcriptome data to screen for candidate RGs that were unaffected by the experimental conditions. Actin-1 (A1), actin-3 (A3), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and translation initiation factor 4a (TIF-4A) were selected for further qPCR verification. The results of RNA-seq and qPCR showed that GAPDH and TIF-4A were suitable RGs after BmNPV challenge or HT stress, whereas TIF-4A was an appropriate RG for BmCPV or BmDNV-Z challenge in silkworms. These results suggested that TIF-4A may be the most appropriate RG for gene expression analysis after challenge with viruses or HT in silkworms.

Enhancement of antiviral capacity of transgenic silkworms against cytoplasmic polyhedrosis virus via knockdown of multiple viral genes

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV), a major pathogen of silkworms, causes serious economic losses in sericulture. The BmCPV genome contains 10 discrete dsRNA segments; among these, S1, S2, S3, S4, S6, and S7 encode virus structural proteins, whereas S5, S8, S9, and S10 encode nonstructural proteins. In an attempt to create an anti-BmCPV silkworm strain, we constructed transgenic RNAi vector pb-CNS for knockdown of S5, S8, S9, and S10, and pb-SNS targeting S1, S2, S4, S5, and S8. Transgenic silkworm line CNS and SNS were generated via microinjection of the practical diapause silkworm strain Furong. Following infection via the oral administration of a high dose of BmCPV, the mortality rates of the nontransgenic control, CNS, and SNS were 91%, 37%, and 41%, respectively. qPCR showed that the viral mRNA content in CNS and SNS was significantly lower than that in the nontransgenic line. The economic traits of CNS and SNS were not affected. These results suggest that the knockdown of multiple BmCPV genes significantly enhances the antiviral capacity of the silkworm.

Overexpression of host plant urease in transgenic silkworms

Bombyx mori and mulberry constitute a model of insect–host plant interactions. Urease hydrolyzes urea to ammonia and is important for the nitrogen metabolism of silkworms because ammonia is assimilated into silk protein. Silkworms do not synthesize urease and acquire it from mulberry leaves. We synthesized the artificial DNA sequence ure-as using the codon bias of B. mori to encode the signal peptide and mulberry urease protein. A transgenic vector that overexpresses ure-as under control of the silkworm midgut-specific P2 promoter was constructed. Transgenic silkworms were created via embryo microinjection. RT-PCR results showed that urease was expressed during the larval stage and qPCR revealed the expression only in the midgut of transgenic lines. Urea concentration in the midgut and hemolymph of transgenic silkworms was significantly lower than in a nontransgenic line when silkworms were fed an artificial diet. Analysis of the daily body weight and food conversion efficiency of the fourth and fifth instar larvae and economic characteristics indicated no differences between transgenic silkworms and the nontransgenic line. These results suggested that overexpression of host plant urease promoted nitrogen metabolism in silkworms.

Construction, complete sequence, and annotation of a BAC contig covering the silkworm chorion locus.

The silkmoth chorion was studied extensively by F.C. Kafatos’ group for almost 40 years. However, the complete structure of the chorion locus was not obtained in the genome sequence of Bombyx mori published in 2008 due to repetitive sequences, resulting in gaps and an incomplete view of the locus. To obtain the complete sequence of the chorion locus, expressed sequence tags (ESTs) derived from follicular epithelium cells were used as probes to screen a bacterial artificial chromosome (BAC) library. Seven BACs were selected to construct a contig which covered the whole chorion locus. By Sanger sequencing, we successfully obtained complete sequences of the chorion locus spanning 871,711 base pairs on chromosome 2, where we annotated 127 chorion genes. The dataset reported here will recruit more researchers to revisit one of the oldest model systems which has been used to study developmentally regulated gene expression. It also provides insights into egg development and fertilization mechanisms and is relevant to applications related to improvements in breeding procedures and transgenesis.

A new approach for comprehensively describing heterogametic sex chromosomes

Abstract Notwithstanding the rapid developments in sequencing techniques, Y and W sex chromosomes have still been mostly excluded from whole genome sequencing projects due to their high repetitive DNA content. Therefore, Y and W chromosomes are poorly described in most species despite their biological importance. Several methods were developed for identifying Y or W-linked sequences among unmapped scaffolds. However, it is not enough to discover functional regions from short unmapped scaffolds. Here, we provide a new and simple strategy based on k-mer comparison for comprehensive analysis of the W chromosome in Bombyx mori. Using this novel method, we effectively assembled de novo 1281 W-derived genome contigs (totaling 1.9 Mbp), and identified 156 W-linked transcript RNAs and 345 W-linked small RNAs. This method will help in the elucidation of mechanisms of sexual development and exploration of W chromosome biological functions, and provide insights into the evolution of sex chromosomes. Moreover, we showed this method can be employed in identifying heterogametic sex chromosomes (W and Y chromosomes) in many other species where genomic information is still scarce.

Enhanced antiviral immunity against Bombyx mori cytoplasmic polyhedrosis virus via overexpression of peptidoglycan recognition protein S2 in transgenic silkworms

ABSTRACT In insect innate immunity, peptidoglycan recognition proteins act as pattern recognition receptors, helping hosts combat invasive microorganisms. Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) is the main silkworm pathogen that invades the midgut columnar cell layer. We previously reported that B. mori peptidoglycan recognition protein S2 (BmPGRP‐S2) was upregulated in silkworm larvae after BmCPV infection. Here, we constructed a transgenic vector overexpressing BmPGRP‐S2 under the control of a midgut‐specific promoter. Transgenic silkworm lines (PGRPS2‐1 and PGRPS2‐2) were generated via embryonic microinjection. BmPGRP‐S2 was successfully overexpressed in transgenic silkworms and BmE cells. After oral inoculation with BmCPV, the mortality of PGRPS2‐1 and PGRPS2‐2 decreased by approximately 36% and 32%, respectively, compared with that of the non‐transgenic line, and BmCPV mRNA contents were significantly lower. In the PGRPS2‐1 line, imd, relish, and the antimicrobial peptide (AMP) genes attacin2, gloverin2, and moricin showed increased expression after viral infection; however, the Toll pathway was not activated. These results indicate that BmPGRP‐S2 overexpression can activate the Imd pathway and induce AMP upregulation, enhancing silkworm antiviral resistance. HighlightsBombyx mori cytoplasmic polyhedrosis virus (BmCPV) is a major silkworm pathogen.BmCPV infection increases peptidoglycan recognition protein S2 (BmPGRP‐S2) levels.High BmPGRP‐S2 induces Imd pathway component and antimicrobial peptide expression.Transgenic lines with high BmPGRP‐S2 were more resistant, with lower mortality.

Increased antiviral capacity of transgenic silkworm via knockdown of multiple genes on Bombyx mori bidensovirus

ABSTRACT Bombyx mori bidensovirus (BmBDV) causes fatal flacherie disease leading to severe economic losses in sericultures. The BmDNV‐Z genome contains two single‐stranded DNA molecules, VD1 and VD2. For generating silkworm lines with antiviral properties, two transgenic RNA interference (RNAi) vectors were constructed. Open reading frames (ORFs) 1–4 of VD1 were knockdown by vector pb‐BDV1 while ORF1a, ORF1b, and ORF3 of VD2 were knockdown by vector pb‐BDV2. Transgenic silkworm lines BDV1‐I and BDV2‐I were generated via RNAi microinjection. Mortality rates of BDV1‐I and BDV2‐I were reduced by 45% and 39%, respectively, and quantitative PCR showed that VD1 and VD2 contents in BDV1‐I and BDV2‐I were significantly lower than in the non‐transgenic line. However, economic traits showed no obvious differences. Thus, knockdown of multiple BmDNV‐Z genes provides strong resistance to BDV1‐I and BDV2‐I lines, and these can be used in sericulture without hampering silk production. HighlightsTransgenic silkworm lines were generated by silencing VD1 and VD2 genes of BmDNV‐Z.Replication of VD1 and VD2 was suppressed in BDV1‐I and BDV2‐I, respectively.Mortalities of BDV1‐I and BDV2‐I were reduced by 45% and 39%, respectively.Economic characteristics of transgenic BDV1‐I and BDV2‐I were unaffected.