The Spatio-Temporal Expression Profiles of Silkworm Pseudogenes Provide Valuable Insights into Their Biological Roles.

Background: Pseudogenes are sequences that have lost the ability to transcribe RNA molecules or encode truncated but possibly functional proteins. While they were once considered to be meaningless remnants of evolution, recent researches have shown that pseudogenes play important roles in various biological processes. However, the studies of pseudogenes in the silkworm, an important model organism, are limited and have focused on single or only a few specific genes. Objective: To fill these gaps, we present a systematic genome-wide studies of pseudogenes in the silkworm. Methods: We identified the pseudogenes in the silkworm using the silkworm genome assemblies, transcriptome, protein sequences from silkworm and its related species. Then we used transcriptome datasets from 832 RNA-seq analyses to construct spatio-temporal expression profiles for these pseudogenes. Additionally, we identified tissue-specifically expressed and differentially expressed pseudogenes to further understand their characteristics. Finally, the functional roles of pseudogenes as lncRNAs were systematically analyzed. Results: We identified a total of 4410 pseudogenes, which were grouped into 4 groups, including duplications (DUPs), unitary pseudogenes (Unitary), processed pseudogenes (retropseudogenes, RETs), and fragments (FRAGs). The most of pseudogenes in the domestic silkworm were generated before the divergence of wild and domestic silkworm, however, the domestication may also involve in the accumulation of pseudogenes. These pseudogenes were clearly divided into 2 cluster, a highly expressed and a lowly expressed, and the posterior silk gland was the tissue with the most tissue-specific pseudogenes (199), implying these pseudogenes may be involved in the development and function of silkgland. We identified 3299 lncRNAs in these pseudogenes, and the target genes of these lncRNAs in silkworm pseudogenes were enriched in the egg formation and olfactory function. Conclusions: This study replenishes the genome annotations for silkworm, provide valuable insights into the biological roles of pseudogenes. It will also contribute to our understanding of the complex gene regulatory networks in the silkworm and will potentially have implications for other organisms as well.

Systematic identification of smORFs in domestic silkworm (Bombyx mori).

The silkworm (Bombyx mori) is not only an excellent model species, but also an important agricultural economic insect. Taking it as the research object, its advantages of low maintenance cost and no biohazard risks are considered. Small open reading frames (smORFs) are an important class of genomic elements that can produce bioactive peptides. However, the smORFs in silkworm had been poorly identified and studied. To further study the smORFs in silkworm, systematic genome-wide identification is essential. Here, we identified and analyzed smORFs in the silkworm using comprehensive methods. Our results showed that at least 738 highly reliable smORFs were found in B. mori and that 34,401 possible smORFs were partially supported. We also identified some differentially expressed and tissue-specific-expressed smORFs, which may be closely related to the characteristics and functions of the tissues. This article provides a basis for subsequent research on smORFs in silkworm, and also hopes to provide a reference point for future research methods for smORFs in other species.

Comparative transcriptome and genome analysis unravels the response of Tatary buckwheat root to nitrogen deficiency.

Tartary buckwheat (Fagopyrum tataricum Garetn.), a dicotyledonous herbaceous crop, has good adaptation to low nitrogen (LN) condition. The plasticity of roots drives the adaption of Tartary buckwheat under LN, but the detailed mechanism behind the response of TB roots to LN remains unclear. In this study, the molecular mechanism of two Tartary buckwheat genotypes' roots with contrasting sensitivity in response to LN was investigated by integrating physiological, transcriptome and whole-genome re-sequencing analysis. LN improved primary and lateral root growth of LN-sensitive genotype, whereas the roots of LN-insensitive genotype showed no response to LN. 2, 661 LN-responsive differentially expressed genes (DEGs) were identified by transcriptome analysis. Of these genes, 17 N transport and assimilation-related and 29 hormone biosynthesis and signaling genes showed response to LN, and they may play important role in Tartary buckwheat root development under LN. The flavonoid biosynthetic genes' expression was improved by LN, and their transcriptional regulations mediated by MYB and bHLH were analyzed. 78 transcription factors, 124 small secreted peptides and 38 receptor-like protein kinases encoding genes involved in LN response. 438 genes were differentially expressed between LN-sensitive and LN-insensitive genotypes by comparing their transcriptome, including 176 LN-responsive DEGs. Furthermore, nine key LN-responsive genes with sequence variation were identified, including FtNRT2.4, FtNPF2.6 and FtMYB1R1. This paper provided useful information on the response and adaptation of Tartary buckwheat root to LN, and the candidate genes for breeding Tartary buckwheat with high N use efficiency were identified.

Cytochrome P450 monooxygenase genes in the wild silkworm, Bombyx mandarina.

Wild (Bombyx mandarina) and domestic silkworms (B. mori) are good models for investigating insect domestication, as 5000 years of artificial breeding and selection have resulted in significant differences between B. mandarina and B. mori. In this study, we improved the genome assemblies to the chromosome level and updated the protein-coding gene annotations for B. mandarina. Based on this updated genome, we identified 68 cytochrome P450 genes in B. mandarina. The cytochrome P450 repository in B. mandarina is smaller than in B. mori. Certain currently unknown key genes, rather than gene number, are critical for insecticide resistance in B. mandarina, which shows greater resistance to insecticides than B. mori. Based on the physical maps of B. mandarina, we located 66 cytochrome P450s on 18 different chromosomes, and 27 of the cytochrome P450 genes were concentrated into seven clusters. KEGG enrichment analysis of the P450 genes revealed the involvement of cytochrome P450 genes in hormone biosynthesis. Analyses of the silk gland transcriptome identified candidate cytochrome P450 genes (CYP306A) involved in ecdysteroidogenesis and insecticide metabolism in B. mandarina.

Comparative genome-wide DNA methylation analysis reveals epigenomic differences in response to heat-humidity stress in Bombyx mori.

DNA methylation is an important epigenetic modification and has been shown to be involved in the response to abiotic stress. However, there are few studies on DNA methylation in insect response to environmental signals. In this study, we conducted a comprehensive comparative analysis of DNA methylation profiles between two silkworm strains with significantly different resistance to heat and humidity by whole-genome bisulfite sequencing (WGBS). We identified, in total, 2934 differentially methylated regions (DMRs) between RT_48h (resistant strain with high-temperature/humidity treatment for 48 h) and ST_48h (sensitive strain with high-temperature/humidity treatment for 48 h) under cytosine context (CG), which corresponded to 1230 DMR-related genes (DMGs), and the DMRs were primarily located in the gene body (exon and intron) region. Gene ontology (GO) and KEGG analysis showed that these DMGs were most significantly enriched in binding, cellular metabolic process, and RNA transport pathways. Moreover, 10 DMGs have been revealed to be involved in the heat-humidity stress response in the silkworm. The results of this study indicated that DNA methylation plays crucial roles in silkworm response to environmental stressors and provides important clues to identify key resistance genes in silkworm under high-temperature/humidity stress response.