Transcriptome analysis in the silkworm Bombyx mori overexpressing piRNA-resistant Masculinizer gene

Abstract

Dosage compensation is a process that produces a similar expression of sex-linked and autosomal genes. In the silkworm Bombyx mori with a WZ sex-determination system, the expression from the single Z in WZ females matches that of ZZ males due to the suppression of Z-linked genes in males. A primary maleness determinant gene, Masculinizer (Masc), is also required for dosage compensation. In females, P-element induced wimpy testis (PIWI) is complexed with the W chromosome-derived female-specific Feminizer (Fem) PIWI-interacting RNA (piRNA) and cleaves Masc mRNA. When Fem piRNA-resistant Masc cDNA (Masc-R) is overexpressed in both sexes, only female larvae are dead during the larval stage. In this study, transcriptome analysis was performed in neonate larvae to examine the effects of Masc-R overexpression on a global gene expression profile. Z-linked genes were globally repressed in Masc-R-overexpressing females due to force-driven dosage compensation. In contrast, Masc-R overexpression had little effect on the expression of Z-linked genes and the male-specific isoform of B. mori insulin-like growth factor II mRNA-binding protein in males, indicating that excessive Masc expression strengthens neither dosage compensation nor maleness in males. Fourteen genes were differentially expressed between Masc-R-overexpressing and control neonate larvae in both sexes, suggesting Masc functions other than dosage compensation and masculinization.HighlightsTranscriptome analysis was performed in Masc-R-overexpressing neonate Bombyx mori larvae.Z-linked genes were globally suppressed in Masc-R-overexpressing females.Masc-R overexpression had little effect on BmImpM expression in males.Several genes may be controlled by Masc-R regardless of sex.

Masculinizer-induced dosage compensation is achieved by transcriptional downregulation of both copies of Z-linked genes in the silkworm, Bombyx mori

“…The libraries were sequenced by Hiseq 2500 sequencing system (Illumina) with 100 bp paired-end reads. Obtained RNA-seq data were submitted to DDBJ under accession number DRA013569 [19]. Adapter sequences and low-quality reads were removed with T rim G alore !…”