Comparative analysis of the silk gland transcriptomes between the domestic and wild silkworms

Identification of genes associated with the silk gland size using multi-omics in silkworm (Bombyx mori)

Silk gland size in silkworms (Bombyx mori) affects silk output. However, the molecular mechanisms by which genes regulate silk gland size remain unclear. In this study, silk glands from three pure silkworm strains (A798, A306 and XH) with different silk gland weight phenotypes were compared using transcriptomics and proteomics to identify differentially expressed genes (DEGs) and proteins (DEPs). When comparing A798 to A306 and A798 to XH, 830 and 469 DEGs were up-regulated, respectively. These genes were related to the gene ontology terms, metabolic process, transport activity and biosynthesis process. In addition, 372 and 302 up-regulated differentially expressed proteins were detected in A798 to A306 and A798 to XH, respectively, related to the gene ontology terms, ribosome and protein export, ribosome and polypeptide biosynthesis processes. Moreover, combined transcriptomics, proteomics and weighted correlation network analyses showed that five genes (BGIBMGA002524, BGIBMGA002629, BGIBMGA005659, BGIBMGA005711 and BGIBMGA010889) were significantly associated with the silk gland weight. Reverse Transcription-quantitative real-time Polymerase Chain Reaction (RT-qPCR) and Enzyme linked immunosorbent assay (ELISA) were used to verify the mRNA and protein expression of five genes in the silk glands and tissues of 18 silkworm strains. The results showed that four genes have higher expression levels in heavier silk glands. These genes are associated with glycogen metabolism, fatty acid synthesis and branched chain amino acid metabolism, thus potentially promoting growth and silk protein synthesis. These findings provide valuable insights into the molecular mechanisms underlying the relationship between silk gland weight and silk yield in silkworms.

Influence of temperature variation on gene expression and cocoon production in Bombyx mori Linnaeus, 1758 (Lepidoptera: Bombycidae)

Silkworms (Bombyx mori) are lepidopterans of economic importance for global silk production. However, factors that directly affect the yield and quality of silkworm cocoon production, such as diseases and temperature fluctuations, cause great economic losses. Knowing how they respond to rearing temperature during the most critical stage of their life cycle (i.e., fifth instar) could provide information on their adaptation and improve silk production. In the current work, we analyzed transcriptional data from two groups of B. mori that were reared at 26 °C and 34 °C throughout the fifth instar. The silkworms and cocoons were weighed. In total, 3115 transcripts were differentially expressed (DE; including 1696 down-regulated and 1419 up-regulated) among the 29,157 sequences found by transcriptome assembly. We emphasize the genes associated with immunological response, transcription factors, silk biosynthesis, and heat shock proteins, among the DE transcripts in response to the temperature conditions. Silkworms reared at 34 °C presented a reduced mean body weight (-0.944 g in comparison to the 26 °C group), which had a direct impact on the weight of cocoons formed and the silk conversion rate. These changes were statistically significant when compared to silkworms reared at 26 °C. Mortality rates (6 and 9 %, at 26 °C and 34 °C, respectively) were similar to those obtained in breeding fields. The findings provide information on the biological processes involved in the temperature response mechanism of silkworms, as well as information that may be used in future climatization processes at rearing facilities and in breeding for improved thermotolerance.

Gene expression differentials driven by mass rearing and artificial selection in black soldier fly colonies

The micro-evolutionary forces that shape genetic diversity during domestication have been assessed in many plant and animal systems. However, the impact of these processes on gene expression, and consequent functional adaptation to artificial environments, remains under-investigated. In this study, whole-transcriptome dynamics associated with the early stages of domestication of the black soldier fly (BSF), Hermetia illucens, were assessed. Differential gene expression (DGE) was evaluated in relation to (i) generational time within the cultured environment (F2 vs. F3), and (ii) two selection strategies [no artificial selective pressure (NS); and selection for greater larval mass (SEL)]. RNA-seq was conducted on 5th instar BSF larvae (n = 36), representing equal proportions of the NS (F2 = 9; F3 = 9) and SEL (F2 = 9; F3 = 9) groups. A multidimensional scaling plot revealed greater gene expression variability within the NS and F2 subgroups, while the SEL group clustered separately with lower levels of variation. Comparisons between generations revealed 898 differentially expressed genes (DEGs; FDR-corrected p < 0.05), while between selection strategies, 213 DEGs were observed (FDR-corrected p < 0.05). Enrichment analyses revealed that metabolic, developmental, and defence response processes were over-expressed in the comparison between F2 and F3 larvae, while metabolic processes were the main differentiating factor between NS and SEL lines. This illustrates the functional adaptations that occur in BSF colonies across generations due to mass rearing; as well as highlighting genic dynamics associated with artificial selection for production traits that might inform future selective breeding strategies.

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 Silk Transcriptomics Illuminates Distinctive Impact of Artificial Selection in Silkworm Modern Breeding

Early domestication and the following improvement are two important processes in the cocoon silk evolution of silkworms. In contrast to early domestication, understanding of the improvement process is still fuzzy. By systematically comparing the larval silk gland transcriptomes of the wild, early domestic, and improved silkworms, we highlighted a novel landscape of transcriptome in the silk glands of improved ones. We first clarified that silk cocoon protein genes were up-regulated in modern breeding but not in early domestication. Furthermore, we found that differentially expressed genes (DEGs) between improved and early domestic silkworms (2711), as well as between improved and wild silkworms (2264), were obviously more than those between the early domestic and wild silkworms (158), with 1671 DEGs specific in the improved silkworm (IS-DEGs). Hierarchical clustering of all the DEGs consistently indicated that improved silkworms were significantly diverged from the early domestic and wild silkworms, suggesting that modern breeding might cause prompt and drastic dynamic changes of gene expression in the silk gland. We further paid attention to these 1671 IS-DEGs and were surprised to find that down-regulated genes were enriched in basic organonitrogen compound biosynthesis, RNA biosynthesis, and ribosome biogenesis processes, which are generally universally expressed, whereas those up-regulated genes were enriched in organonitrogen compound catabolic processes and functions involving in the dynamic regulation of protein post-translation of modification. We finally highlighted one candidate improvement gene among these up-regulated IS-DEGs, i.e., GDAP2, which may play roles in silk behavior and the overall robustness of the improved silkworm. The findings strongly suggest that modern breeding may facilitate effective control of the basic consumption of nitrogen and a stronger switch of nitrogen resources from other tissues to the silk glands, for an efficient supply for silk production, and implies the importance of brain behavior and robustness in silk yield improvement of modern breeding.

Molecular Characterization of the Functional Genes Associated with Silk Assembly, Transport, and Protection in the Silk Glands of Popular Multivoltine Breeds of Silkworm Bombyx mori. L

Bombyx mori is an agriculturally important insect used extensively for silk production. India, especially the eastern regions, is mostly dependent on the multivoltine breeds of silkworm Bombyx mori and their hybrids/crossbreeds. The multivoltine breeds are indigenous and superior in survival and hardiness but are relatively inferior in terms of qualitative traits, typically the silk quality. Therefore, it is highly relevant to understand the mechanism of silk production in the multivoltine breeds to decipher the reasons for the inferior quality of silk produced by the multivoltine breeds and thus gain leads to improve the quality of silk production in multivoltine breeds. With this background, study was carried to identify differential expression of the major genes associated with silk proteins in the silk gland region of the popular multivoltine breeds. Our results indicated that although fib-L, fib-H, Sericins, and P25 are the major genes associated with silk filament, a few other genes associated with silk assembly, transport, and protection in the silk glands are the ones that largely contribute towards efficient silk production. The differential expression of these genes had a major effect on the movement of silk proteins within the silk gland and the efficiency of silk production as well. The Pearson correlation revealed a positive correlation amongst the genes dealt with in this study, indicating that the concurrent increase in expression of both the types of genes in the silk glands, significantly improves the silk production.

Global Profiling of Genes Expressed in the Silk Glands of Philippine-Reared Mulberry Silkworms (Bombyx mori)

RNA sequencing was used to assemble transcriptome data for Philippine-reared silkworm and compare gene expression profiles of strains reared in high- and low-temperature environments. RNA was isolated from the silk glands of fifth instar larvae and mRNA-enriched libraries were sequenced using Illumina NextSeq 500. Transcriptome reads were assembled using reference-based and de novo assemblers, and assemblies were evaluated using different metrics for transcriptome quality, including the read mapping rate, E90N50, RSEM-eval, and the presence of single-copy orthologs. All transcriptome assemblies were able to reconstruct >40,000 transcripts. Differential expression analysis found 476 differentially expressed genes (DEGs; 222 upregulated, 254 downregulated) in strains reared in different temperatures. Among the top DEGs were myrosinase, heat shock proteins, serine protease inhibitors, dehydrogenases, and regulators of the juvenile hormone. Validation of some of the top DEGs using qPCR supported the findings of the in silico analysis. GO term enrichment analysis reveals an overrepresentation of GO terms related to nucleotide metabolism and biosynthesis, lipid and carbohydrate metabolic processes, regulation of transcription, nucleotide binding, protein binding, metal binding, catalytic activity, oxidoreductase activity, and hydrolase activity. The data provided here will serve as a resource for improving local strains and increasing silk production of Philippine-reared B. mori strains.

Bombyx mori Vps13d is a key gene affecting silk yield

Bombyx mori is an important economic insect, its economic value mainly reflected in the silk yield. The major functional genes affecting the silk yield of B. mori have not been determined yet. Bombyx mori vacuolar protein sorting-associated protein 13d (BmVps13d) has been identified, but its function is not reported. In this study, BmVps13d protein shared 30.84% and 34.35% identity with that of in Drosophila melanogaster and Homo. sapiens, respectively. The expressions of BmVps13d were significantly higher in the midgut and silk gland of JS (high silk yield) than in that of L10 (low silk yield). An insertion of 9 bp nucleotides and two deficiencies of adenine ribonucleotides in the putative promoter region of BmVps13d gene in L10 resulted in the decline of promoter activity was confirmed using dual luciferase assay. Finally, the functions of BmVps13d in B. mori were studied using the CRISPR/Cas9 system, and the mutation of BmVps13d resulted in a 24.7% decline in weight of larvae, as well as a 27.1% (female) decline and a 11.8% (male) decline in the silk yield. This study provides a foundation for studying the molecular mechanism of silk yield and breeding the silkworm with high silk yield.

Comparative analysis of testis transcriptome between a genetic male sterile line (GMS) and its wild-type 898WB in silkworm, Bombyx mori

The silkworm, Bombyx mori, is an important model organism of lepidopteran insects, and its testis is a main male reproductive organ and spermatogenesis place. Studying the testis helps to understand the mechanisms of genetic sterility of lepidopteran insects and to achieve sterile insect technique (SIT) for pest control. Herein, we performed a comparative transcriptome analysis of testes between three biological replicates of the GMS mutant and wild strain 898WB, respectively. In total, 1872 up-regulated genes and 1823 down-regulated genes were identified in the testis of the GMS mutant. Several genes contribute significantly to spermatogenesis and testis development, such as "serine/threonine protein kinase", "organic cation transporter protein", "tyrosine protein kinase", "lncRNAs" and "immune-associated genes". The KEGG pathway analysis shows that the DEGs were annotated to 123 pathways, and 10 pathways were significantly enriched, such as "metabolic pathway", "biosynthesis of amino acids", and "phagosome-lysosome pathway", which are associated with testis development and spermatogenesis. The results of the qPCR expression were consistent with the RNA-seq data, which shows that the RNA-seq results were accurate. The DEGs of the testes between GMS mutant and 898WB were screened by RNA-Seq technology, which provides a reliable reference to understand the molecule mechanism of male sterility of the GMS mutant.

The transcriptomic responses of blunt snout bream (Megalobrama amblycephala) to acute hypoxia stress alone, and in combination with bortezomib

BACKGROUND

Blunt snout bream (Megalobrama amblycephala) is sensitive to hypoxia. A new blunt snout bream strain, "Pujiang No.2", was developed to overcome this shortcoming. As a proteasome inhibitor, bortezomib (PS-341) has been shown to affect the adaptation of cells to a hypoxic environment. In the present study, bortezomib was used to explore the hypoxia adaptation mechanism of "Pujiang No.2". We examined how acute hypoxia alone (hypoxia-treated, HN: 1.0 mg·L^- 1), and in combination with bortezomib (hypoxia-bortezomib-treated, HB: Use 1 mg bortezomib for 1 kg fish), impacted the hepatic ultrastructure and transcriptome expression compared to control fish (normoxia-treated, NN).

RESULTS

Hypoxia tolerance was significantly decreased in the bortezomib-treated group (LOE_crit, loss of equilibrium, 1.11 mg·L^- 1 and 1.32 mg·L^- 1) compared to the control group (LOE_crit, 0.73 mg·L^- 1 and 0.85 mg·L^- 1). The HB group had more severe liver injury than the HN group. Specifically, the activities of alanine aminotransferase (ALT) and aspartate aminotransferase (AST) in the HB group (52.16 U/gprot, 32 U/gprot) were significantly (p < 0.01) higher than those in the HN group (32.85 U/gprot, 21. 68 U/gprot). In addition, more severe liver damage such as vacuoles, nuclear atrophy, and nuclear lysis were observed in the HB group. RNA-seq was performed on livers from the HN, HB and NN groups. KEGG pathway analysis disclosed that many DEGs (differently expressed genes) were enriched in the HIF-1, FOXO, MAPK, PI3K-Akt and AMPK signaling pathway and their downstream.

CONCLUSION

We explored the adaptation mechanism of "Pujiang No.2" to hypoxia stress by using bortezomib, and combined with transcriptome analysis, accurately captured the genes related to hypoxia tolerance advantage.

Co-Expression Network and Time-Course Expression Analyses to Identify Silk Protein Regulatory Factors in Bombyx mori

Simple Summary

Previous studies have reported how the silk production ability of Bombyx mori can be enhanced, but the mechanism that regulates silk protein genes remains unclear. We performed co-expression network analysis using networkz, an in-house program, which led to the identification of 91 transcription factors were co-expressed with silk protein genes. Of them, 13 transcripts were identified to be novel regulatory factors by time-course expression analysis during the fifth instar larvae stage. Their expression patterns were highly relevant to those of silk protein genes. Our results suggest that the two-step expression screening was robust and highly sensitive to screen relative genes, and a complex mechanism regulates silk protein production in B. mori. The novel candidates that were identified herein can serve as key genes to develop methods to enhance the silk protein production ability of B. mori.

Abstract

Bombyx mori is an important economic insect and an animal model in pharmacomedical research. Although its physiology has been studied for many years, the mechanism via which silk protein genes are regulated remains unclear. In this study, we performed two-step expression screening, namely co-expression network and time-course expression analyses to screen silk protein regulation factors. A co-expression network analysis using RNA-seq data that were obtained from various tissues, including the silk glands of B. mori, was performed to identify novel silk protein regulatory factors. Overall, 91 transcription factors, including some known ones, were found to be co-expressed with silk protein genes. Furthermore, time-course expression analysis during the fifth instar larvae stage revealed that the expression pattern of 13 novel transcription factors was highly relevant to that of silk protein genes and their known regulatory factor genes. In particular, the expression peak of several transcription factors (TFs) was detected before the expression of silk protein genes peak. These results indicated that a larger number of genes than expected may be involved in silk protein regulation in B. mori. Functional analyses of function-unknown transcription factors should enhance our understanding of this system.

Transcriptional analysis of Bemisia tabaci MEAM1 cryptic species under the selection pressure of neonicotinoids imidacloprid, acetamiprid and thiamethoxam

BACKGROUND

Neonicotinoids are widely applied in the control of the destructive agricultural pest Bemisia tabaci, and resistance against these chemicals has become a common, severe problem in the control of whiteflies. To investigate the molecular mechanism underlying resistance against nenonicotinoids in whiteflies, RNA-seq technology was applied, and the variation in the transcriptomic profiles of susceptible whiteflies and whiteflies selected by imidacloprid, acetamiprid and thiamethoxam treatment was characterized.

RESULTS

A total of 90.86 GB of clean sequence data were obtained from the 4 transcriptomes. Among the 16,069 assembled genes, 584, 110 and 147 genes were upregulated in the imidacloprid-selected strain (IMI), acetamiprid-selected strain (ACE), and thiamethoxam (THI)-selected strain, respectively, relative to the susceptible strain. Detoxification-related genes including P450s, cuticle protein genes, GSTs, UGTs and molecular chaperone HSP70s were overexpressed in the selected resistant strains, especially in the IMI strain. Five genes were downregulated in all three selected resistant strains, including 2 UDP-glucuronosyltransferase 2B18-like genes (LOC 109030370 and LOC 109032577).

CONCLUSIONS

Ten generations of selection with the three neonicotinoids induced different resistance levels and gene expression profiles, mainly involving cuticle protein and P450 genes, in the three selected resistant whitefly strains. The results provide a reference for research on resistance and cross-resistance against neonicotinoids in B. tabaci.

Comparative Genome Analysis of Genes Regulating Compound Inflorescences in Tomato

Inflorescences are the main factor affecting fruit yield. The quantity and quality of inflorescences are closely related to fruit quality and yield. The presence of compound inflorescences in cherry tomatoes is well established, and it has been discovered by chance that compound racemes also exist in tomatoes. To explore the formation of compound inflorescences in tomato, transcriptome sequencing was performed on Moneymaker (MM) and Compound Inflorescence (CI) plants. In-florescences were collected in three periods (early, middle and late) in three replicates, for a total of 18 samples. Data analysis showed that the DEGs were most enriched in metabolic pathways and plant hormone signal transduction pathways. The DEGs were also enriched in the cell cycle pathway, photosynthesis pathway, carbon metabolism pathway and circadian rhythm pathway. We found that the FALSIFLORA (FA), COMPOUND INFLORESCENCE (S) and ANANTHA (AN) genes were involved in compound inflorescence development, not only revealing novel genes but also providing a rich theoretical basis for compound inflorescence development.

Comparative transcriptome provides insights into the selection adaptation between wild and farmed foxes

The silver fox and blue fox are economically important fur species and were domesticated by humans from their wild counterparts, the arctic fox and red fox, respectively. Farmed foxes show obvious differences from their wild counterparts, including differences in physiology, body size, energy metabolism, and immunity. However, the molecular mechanisms underlying these differences are presently unclear. In this study, we built transcriptome libraries from multiple pooled tissues for each species of farmed fox, used RNA-seq to obtain a comprehensive dataset, and performed selection analysis and sequence-level analyses of orthologous genes to identify the genes that may be influenced by human domestication. More than 153.3, 248.0, 81.6, and 65.8 million clean reads were obtained and assembled into a total of 118,577, 401,520, 79,900, and 186,988 unigenes with an average length range from 521 to 667 bp for AF, BF, RF, and SF, respectively. Selective pressure analysis showed that 11 and 14 positively selected genes were identified, respectively, in the two groups (AF vs. BF and RF vs. SF). Several of these genes were associated with natural immunity (CFI and LRRFIP1), protein synthesis (GOLGA4, CEP19 and SLC35A2), and DNA damage repair (MDC1). Further functional enrichment analyses demonstrated that two positively selected genes (ACO1 and ACAD10) were involved in metabolic process (GO:0008152, p-value = .032), representing a significant enrichment. Sequence analysis of 117 orthologous genes shared by the two groups showed that the LEMD2, RRBP1, and IGBP1 genes might be affected by artificial selection in farmed foxes, with mutation sites located within sequences that are otherwise highly conserved across most mammals. Our results provide a valuable transcriptomic resource for future genetic studies and improvement in the assisted breeding of foxes and other farmed animals.

Gene expression profiling reveals candidate genes for defining spider silk gland types

Molecular studies of the secretory glands involved in spider silk production have revealed candidate genes for silk synthesis and a complicated history of spider silk gene evolution. However, differential gene expression profiles of the multiple silk gland types within an individual orb-web weaving spider are lacking. Each of these gland types produces a functionally distinct silk type. Comparison of gene expression among spider silk gland types would provide insight into the genes that define silk glands generally from non-silk gland tissues, and the genes that define silk glands from each other. Here, we perform 3' tag digital gene expression profiling of the seven silk gland types of the silver garden orb weaver Argiope argentata. Five of these gland types produce silks that are non-adhesive fibers, one silk includes both fibers and glue-like adhesives, and one silk is exclusively glue-like. We identify 1275 highly expressed, significantly upregulated, and tissue specific silk gland specific transcripts (SSTs). These SSTs include seven types of spider silk protein encoding genes known as spidroin genes. We find that the fiber-producing major ampullate and minor ampullate silk glands have more similar expression profiles than any other pair of glands. We also find that a subset of the SSTs is enriched for transmembrane transport and oxidoreductases, and that these transcripts highlight differences and similarities among the major ampullate, minor ampullate, and aggregate silk glands. Furthermore, we show that the wet glue-producing aggregate glands have the most unique SSTs, but still share some SSTs with fiber producing glands. Aciniform glands were the only gland type to share a majority of SSTs with other silk gland types, supporting previous hypotheses that duplication of aciniform glands and subsequent divergence of the duplicates gave rise to the multiple silk gland types within an individual spider.

Transcriptomic analysis of the bagworm moth silk gland reveals a number of silk genes conserved within Lepidoptera

Lepidopteran insects produce cocoons with unique properties. The cocoons are made of silk produced in the larval tissue silk gland and our understanding of the silk genes is still very limited. Here, we investigated silk genes in the bagworm moth Eumeta variegata, a species that has recently been found to produce extraordinarily strong and tough silk. Using short-read transcriptomic analysis, we identified a partial sequence of the fibroin heavy chain gene and its product was found to have a C-terminal structure that is conserved within nonsaturniid species. This is in accordance with the presence of fibroin light chain/fibrohexamerin genes and it is suggested that the bagworm moth is producing silk composed of fibroin ternary complex. This indicates that the fibroin structure has been evolutionarily conserved longer than previously thought. Other than fibroins we identified candidates for sericin genes, expressed strongly in the middle region of the silk gland and encoding serine-rich proteins, and other silk genes, that are structurally conserved with other lepidopteran homologues. The bagworm moth is thus considered to be producing conventional lepidopteran type of silk. We further found a number of genes expressed in a specific region of the silk gland and some genes showed conserved expression with Bombyx mori counterparts. This is the first study allowing comprehensive silk gene identification and expression analysis in the lepidopteran Psychidae family and should contribute to the understanding of silk gene evolution as well as to the development of novel types of silk.

Morus alba: Host reaction for Meloidogyne javanica, biological nematicides assessment and study of these relationships with yield and quality of leaves, cocoon and health of the silkworm

Root-knot nematodes cause damage to several crops and the importance of each species can vary according with the crop and the agricultural region. In Brazil, Meloidogyne javanica is one of the most important nematode species parasitizing mulberry. To define management strategies, it is important to know if the crop species is damaged by the parasitism of the nematode and the best choices for control, as the use of nematicides. Biological nematicides have been extensively used in Brazil, but no information regarding its efficiency to control M. javanica in mulberry is available. Besides, it is not known if biological nematicides could improve the quality of leaves or if they alter the nutrient composition of leaves, which could interfere in the development of the silkworms that are feed with these leaves or in the quality of the silk produced. With the aim to address these questions, we propose a study that will start in the phenotyping of the main Brazilian mulberry cultivars to Meloidogyne species, passing through the test of efficiency of biological nematicides in the control of M. javanica in mulberry cultivar Miura, evaluation of the amount and quality of leaves produced and, using these leaves to feed silkworms, in the analyzes of the impact of these diet in the health of silkworms, and in the production and quality of the silk.

Gene mapping reveals the association between tyrosine protein kinase Abl1 and the silk yield of Bombyx mori

The Z chromosome of the silkworm contains a major gene that influences silk yield. This major locus on chromosome Z accounts for 35.10% of the phenotypic variance. The location and identification of the gene have been a focus of silkworm genetics research. Unfortunately, identification of this gene has been difficult. We used extreme phenotype subpopulations and selected from a backcross population, BC₁ M, which was obtained using the high-yield strain 872B and the low-yield strain IS-Dazao as parents, for mapping the gene on the chromosome Z. The candidate region was narrowed down to 134 kb at the tip of the chromosome. BmAbl1 in this region correlated with silk gland development by spatiotemporal expression analysis. This gene was differentially expressed in the posterior silk glands of the high- and low-yield strains. In BmAbl1, an insertion-deletion (indel) within the 10th exonic region and an SNP within the 6th intronic region were detected and shown to be associated with cocoon shell weight in 84 Bombyx mori strains with different yields. Nucleotide diversity analysis of BmAbl1 and its 50 kb flanking regions indicated that BmAbl1 has experienced strong artificial selection during silkworm domestication. This study is the first to identify the genes controlling silk yield in the major QTL of the Z chromosome using forward genetics.

Silk of the common clothes moth, Tineola bisselliella, a cosmopolitan pest belonging to the basal ditrysian moth line

Many lepidopteran larvae produce silk secretions to build feeding tubes and cocoons that play important protective roles in their lives. Recent research on the silk of bombycoid and pyralid moths has shown that it contains several highly abundant silk components with remarkable mechanical properties. It was also found to contain a number of other proteins of which the functions have yet to be identified. To gain an overview of the silk composition in more primitive lepidopteran species and to identify the core silk components common to most species, we analyzed the cocoon proteins of Tineola bisselliella, which belongs to the basal ditrysian moth line. Using de novo transcriptome sequencing combined with mass spectrometry (MS)-based proteomics, we detected more than 100 secretory proteins in the silk cocoons. Fibroin, sericins, and protease inhibitors were found to be the most abundant proteins, along with several novel candidate silk components. We also verified the tissue and developmental stage specificity of the silk protein expression and characterized the morphology of both the silk glands and silk in T. bisselliella. Our study provides a detailed analysis of silk in the primitive moth, expands the known set of silk-specific genes in Lepidoptera, and helps to elucidate their evolutionary relationships.

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.

Identification of a nomogram based on an 8-lncRNA signature as a novel diagnostic biomarker for head and neck squamous cell carcinoma

Long noncoding RNAs (lncRNAs) have been proposed as diagnostic or prognostic biomarkers of head and neck squamous carcinoma (HNSCC). The current study aimed to develop a lncRNA-based prognostic nomogram for HNSCC. LncRNA expression profiles were downloaded from The Cancer Genome Atlas (TCGA) database. After the reannotation of lncRNAs, the differential analysis identified 253 significantly differentially expressed lncRNAs in training set TCGA-HNSC (n = 300). The prognostic value of each lncRNA was first estimated in univariate Cox analysis, and 41 lncRNAs with P < 0.05 were selected as seed lncRNAs for Cox LASSO regression, which identified 11 lncRNAs. Multivariate Cox analysis was used to establish an 8-lncRNA signature with prognostic value. Patients in the high-signature score group exhibited a significantly worse overall survival (OS) than those in the low-signature score group, and the area under the receiver operating characteristic (ROC) curve for 3-year survival was 0.74. Multivariable Cox regression analysis among the clinical characteristics and signature scores suggested that the signature is an independent prognostic factor. The internal validation cohort, external validation cohort, and 102 HNSCC specimens quantified by qRT-PCR successfully validate the robustness of our nomogram.

Whole-genome resequencing reveals loci under selection during silkworm improvement

Breeding or genetic improvement refers to the process of artificial selection following domestication; as such, it has had a major influence on modern agriculture and animal production. Improvement generally focuses on traits that greatly affect the economic performance. Therefore, understanding the genetic basis underlying improvement will contribute to the identification of genes controlling economic traits and will facilitate future crop and animal breeding. However, genome-wide study of the molecular basis underlying improvement remains rare. The silkworm is a unique, entirely domesticated economically important invertebrate; genetic improvement has had a huge effect on the silkworm regarding silk-related traits. Herein, we performed whole-genomic sequencing on local and genetically improved silkworm lines to identify the genomic regions under strong selection in silkworm breeding/improvement. By genomic-wide selective sweeping analysis, we identified 24 genomic regions with strong selection signals, eight of which contained 13 candidate genes underlying silkworm breeding. Interestingly, six of these genes were annotated with functions related to neural signal response. Among the six genes, BGIBMGA004050 encodes silkworm CREB-regulated_transcription_coactivator_1 (BmCRTC1), which was reported to be involved in energy-sensing pathways. These results suggested that improvement may have affected the nervous system of the silkworm. This research will provide new insights into the genetic basis underlying the genetic improvement of silkworms and possibly of other species.

The beta-1, 4-N-acetylglucosaminidase 1 gene, selected by domestication and breeding, is involved in cocoon construction of Bombyx mori

Holometabolous insects have distinct larval, pupal, and adult stages. The pupal stage is typically immobile and can be subject to predation, but cocoon offers pupal protection for many insect species. The cocoon provides a space in which the pupa to adult metamorphosis occurs. It also protects the pupa from weather, predators and parasitoids. Silk protein is a precursor of the silk used in cocoon construction. We used the silkworm as a model species to identify genes affecting silk protein synthesis and cocoon construction. We used quantitative genetic analysis to demonstrate that β-1,4-N-acetylglucosaminidase 1 (BmGlcNase1) is associated with synthesis of sericin, the main composite of cocoon. BmGlcNase1 has an expression pattern coupled with silk gland development and cocoon shell weight (CSW) variation, and CSW is an index of the ability to synthesize silk protein. Up-regulated expression of BmGlcNase1 increased sericin content by 13.9% and 22.5% while down-regulation reduced sericin content by 41.2% and 27.3% in the cocoons of females and males, respectively. Genomic sequencing revealed that sequence variation upstream of the BmGlcNase1 transcriptional start site (TSS) is associated with the expression of BmGlcNase1 and CSW. Selective pressure analysis showed that GlcNase1 was differentially selected in insects with and without cocoons (ω₁ = 0.044 vs. ω₂ = 0.154). This indicates that this gene has a conserved function in the cocooning process of insects. BmGlcNase1 appears to be involved in sericin synthesis and silkworm cocooning.

The cocoon provides a protected space for the metamorphosis of many insect species. Silk protein is a precursor of the fiber used for cocoon construction. Deciphering the genetic basis underlying silk protein synthesis will improve our understanding of cocoon construction and the adaptations of species that construct cocoons. We used the silkworm (Bombyx mori) as a model to identify genes affecting silk protein synthesis and cocoon construction. Quantitative genetic analysis was used to show that β-1,4-N-acetylglucosaminidase 1 (BmGlcNase1), a gene selected during silkworm domestication and breeding, is associated with sericin synthesis. Transgenic-based functional validation confirmed that BmGlcNase1 positively regulates sericin content in the silkworm cocoon. The selective pressure of GlcNase1 in the evolution of insects with cocoons is higher than those without cocoons. This indicates that it has a conserved function in the cocooning process. These results reveal aspects of the genetic basis of silk protein synthesis and the cocoon construction of insects.

Metabolic Characterisation of the Midgut of Bombyx mori Varieties after BmNPV Infection Using GC-MS-Based Metabolite Profiling

Bombyx mori nucleopolyhedrovirus (BmNPV) is a silkworm disease that is especially harmful to cocoon production and seriously restricts sericultural development. Our laboratory successfully cultivated a new highly BmNPV-resistant silkworm variety, Huakang 2; however, its mechanism of BmNPV resistance remains unclear. To understand its resistance mechanism, we conducted a metabolomic and transcriptomic study of the midgut of silkworm varieties, Baiyu N and Baiyu after BmNPV infection. We identified 451 differential metabolites, which were mostly comprised of small molecules, such as saccharides, acids, amines, alcohols, and glycosides. We found that the primary differences in disease resistance between the silkworm varieties are metabolic-pathways, tryptophan metabolism, oxidative phosphorylation, ABC-transporters, beta-alanine metabolism, and phenylalanine metabolism. Combined analysis with transcriptomic data suggested that tryptophan metabolism and oxidative phosphorylation are closely related to the silkworms' BmNPV resistance. We hypothesize that the roles of the two metabolic pathways in the BmNPV resistance mechanism might be the following: Oxidative phosphorylation generates a large amount of adenosine triphosphate (ATP) in response to BmNPV infection to provide silkworms the energy required for establishing BmNPV resistance. Tryptophan metabolism then activates the aryl hydrocarbon receptor (AhR) through the exogenous virus BmNPV, which activates the silkworm's immune system to defeat BmNPV infections.

A transcriptome atlas of silkworm silk glands revealed by PacBio single-molecule long-read sequencing

The silk gland of the silkworm Bombyx mori is a specialized organ where silk proteins are efficiently synthesized under precise regulation that largely determines the properties of silk fibers. To understand the genes involved in the regulation of silk protein synthesis, considerable research has focused on the transcripts expressed in silk glands; however, the complete transcriptome profile of this organ has yet to be elucidated. Here, we report a full-length silk gland transcriptome obtained by PacBio single-molecule long-read sequencing technology. In total, 11,697 non-redundant transcripts were identified in mixed samples of silk glands dissected from larvae at five developmental stages. When compared with the published reference, the full-length transcripts optimized the structures of 3002 known genes, and a total of 9061 novel transcripts with an average length of 2171 bp were detected. Among these, 1403 (15.5%) novel transcripts were computationally revealed to be lncRNAs, 8135 (89.8%) novel transcripts were annotated to different protein and nucleotide databases, and 5655 (62.4%) novel transcripts were predicted to have complete ORFs. Furthermore, we found 1867 alternative splicing events, 2529 alternative polyadenylation events, 784 fusion events and 6596 SSRs. This study provides a comprehensive set of reference transcripts and greatly revises and expands the available silkworm transcript data. In addition, these data will be very useful for studying the regulatory mechanisms of silk protein synthesis.

Selection of suitable reference genes for qRT-PCR expression analysis of Codonopsis pilosula under different experimental conditions

Codonopsis pilosula is a well-known medicinal plant. Although its transcriptome sequence has been published, suitable reference genes have not been systematically identified for conducting expression analyses via quantitative real-time polymerase chain reaction (qRT-PCR). To screen appropriate genes for use with this species, we applied four different methods-GeNorm, NormFinder, BestKeeper, and RefFinder-to evaluate the stability of 13 candidates: CpiEF1Bb, CpiCACS, CpiF-Box, Cpiβ-Tubulin, CpiGAPDH, CpiActin2, CpiAPT1, CpiActin7, CpiActin8, CpiRPL6, CpiHAF1, CpiTubulin6, and CpiUBQ12. Expression was examined by qRT-PCR for various tissue types, chemical treatments, and developmental stages. For all tested samples, CpiGAPDH proved to be the most stable. Comprehensive analysis indicated that the most stable internal reference genes were CpiF-Box and CpiCACS in different tissues and at different developmental stages, respectively. Under NaCl stress, CpiAPT1 was the best internal reference gene. For methyl jasmonate and abscisic acid treatments, CpiGAPDH and CpiF-Box, respectively, presented the highest degree of expression stability. Based on these findings, we chose CpiSPL9 as the target gene for validating the suitability of these selected reference genes. All of these results provide a foundation for accurate quantification of expression levels by genes of interest in C. pilosula.

Genetic and genomic analysis for cocoon yield traits in silkworm

Domestic species provides a powerful model for examining genetic mechanisms in the evolution of yield traits. The domestic silkworm (Bombyx mori) is an important livestock species in sericulture. While the mechanisms controlling cocoon yield are largely unknown. Here, using B. mori and its wild relative B. mandarina as intercross parents, 100 BC₁ individuals were sequenced by restriction site-associated DNA sequencing (RAD-Seq). The linkage map contained 9,632 markers was constructed. We performed high-resolution quantitative trait locus (QTL) mapping for four cocoon yield traits. A total of 11 QTLs were identified, including one yield-enhancing QTL from wild silkworm. By integrating population genomics and transcriptomic analysis with QTLs, some favourable genes were revealed, including 14 domestication-related genes and 71 differentially expressed genes (DEGs) in the fifth-instar larval silk gland transcriptome between B. mori and B. mandarina. The relationships between the expression of two important candidate genes (KWMTBOMO04917 and KWMTBOMO12906) and cocoon yield were supported by quantitative real-time PCR (qPCR). Our results provide some new insights into the molecular mechanisms of complex yield traits in silkworm. The combined method might be an efficient approach for identifying putative causal genes in domestic livestock and wild relatives.

A Comparison of Co-expression Networks in Silk Gland Reveals the Causes of Silk Yield Increase During Silkworm Domestication

Long-term domestication and selective breeding have increased the silk yield of the domestic silkworm (Bombyx mori) by several times the amount of the silk yield of its wild ancestor (Bombyx mandarina). However, little is known about the molecular mechanisms behind the increase in silk yield during domestication. Based on dynamic patterns of functional divergence in the silk gland between domestic and wild silkworms, we found that at early and intermediate stages of silk gland development, the up-regulated genes of the domestic silkworm were mainly involved in DNA integration, nucleic acid binding, and transporter activity, which are related to the division and growth of cells. This has led to the posterior silk gland (PSG) of the domestic silkworm having significantly more cells (“factories” of fibroin protein synthesis) than that of the wild silkworm. At the late stage of silk gland development, the up-regulated genes in the domestic silkworm was enriched in protein processing and ribosome pathways, suggesting protein synthesis efficiency is greatly improved during silkworm domestication. While there was an increase in fibroin protein synthesis, the production of sericin protein was simultaneously reduced in the silk gland of the domestic silkworm. This reflects that domestic and wild silkworms have been under different selection pressures. Importantly, we found that the network co-expressed with the silk-coding genes of the domestic silkworm was larger than that of the wild silkworm. Furthermore, many more genes co-expressed with silk-coding genes in the domestic silkworm were subjected to artificial selection than those in the wild silkworm. Our results revealed that the increase of silk yield during silkworm domestication is involved in improvement of a biological system which includes not only expansion of “factories” (cells of PSG) of protein synthesis, but also a high expression of silk-coding genes and silk production-related genes such as biological energy, transport, and ribosome pathway genes.

Comparative transcriptome analysis reveals the response mechanism of Cf-16-mediated resistance to Cladosporium fulvum infection in tomato

BACKGROUND

Leaf mold disease caused by Cladosporium fulvum is a serious threat affecting the global production of tomato. Cf genes are associated with leaf mold resistance, including Cf-16, which confers effective resistance to leaf mold in tomato. However, the molecular mechanism of the Cf-16-mediated resistance response is largely unknown.

RESULTS

We performed a comparative transcriptome analysis of C. fulvum-resistant (cv. Ontario7816) and C. fulvum-susceptible (cv. Moneymaker) tomato cultivars to identify differentially expressed genes (DEGs) at 4 and 8 days post inoculation (dpi) with C. fulvum. In total, 1588 and 939 more DEGs were found in Cf-16 tomato than in Moneymaker at 4 and 8 dpi, respectively. Additionally, 1350 DEGs were shared between the 4- and 8-dpi Cf-16 groups, suggesting the existence of common core DEGs in response to C. fulvum infection. The up-regulated DEGs in Cf-16 tomato were primarily associated with defense processes and phytohormone signaling, including salicylic acid (SA) and jasmonic acid (JA). Moreover, SA and JA levels were significantly increased in Cf-16 tomato at the early stages of C. fulvum infection. Contrary to the previous study, the number of up-regulated genes in Cf-16 compared to Cf-10 and Cf-12 tomatoes was significantly higher at the early stages of C. fulvum infection.

CONCLUSION

Our results provide new insight into the Cf-mediated mechanism of resistance to C. fulvum, especially the unique characteristics of Cf-16 tomato in response to this fungus.

Genetic Manipulation of MicroRNAs in the Silk Gland of Silkworm, Bombyx Mori

Background

MicroRNAs (miRNAs) are a class of non-coding RNAs with important post-transcriptional regulatory functions. To reveal the function of miRNAs in vivo, the critical step is to change their expression levels in the tissues or organs. In this work, we explored the application of several important genetic techniques in altering the expression of silk gland-specific miR-274 of silkworm (Bombyx mori).

Results

Injection of synthesized microRNA mimics and antagomirs exerted no effect on the expression of miR-274 in the silk gland, miR-274 sponge specifically absorbed miR-274 and down-regulated its expression, transgenic overexpression of miR-274 precursor significantly up-regulated miR-274, and finally tissue-specific CRISPR/Cas9 system achieved deletion of miR-274.

Conclusions

A practical technical system was established for studying the functions of miRNAs in silk gland of Bombyx mori. Our research provides methodological support for the functional study of miRNAs and other noncoding RNAs in the silk gland and more organs in other species.

Synergism of open chromatin regions involved in regulating genes in Bombyx mori

The dynamic variability of transcription factors (TFs) and their binding sites makes it challenging to conduct genome-wide transcription regulation research. The silkworm Bombyx mori, which produces silk, is one of the most valuable model insects in the order Lepidoptera. The "opening" and "closing" of chromatin in different silk yield strains is associated with changes in silk production, making this insect a good model for studying the transcriptional regulation of genes. However, few studies have examined the open chromatin regions (OCRs) of silkworms, and studying OCR synergism and their function in silk production remains challenging. Here, we performed formaldehyde-assisted isolation of regulatory elements (FAIRE) to isolate OCRs from the silk glands of fifth-instar larvae of the DaZao and D872 strains. In total, 128,908 high confidence OCRs were identified and approximately 80% of OCRs were located in non-coding regions. OCRs upregulated adjacent genes and showed signal-dependent vulnerability to single-nucleotide polymorphisms. Mid- and low-signal OCRs were more likely to have single-nucleotide polymorphisms (SNP). Further, OCRs interacted with each other within a distance of 5 kb. We named the OCR interaction complex as the "cluster of related regions" (COREs). The functions of the CORE and its harbored OCRs showed some differences. Additionally, COREs enriched many silk protein synthesis-associated genes, some of which were upregulated. This study identified numerous high confidence regulation sites and synergistic regulatory modes of OCRs that affect adjacent genes. These results provide insight into silkworm transcriptional regulation and improve our understanding of cis-element cooperation.

Decrease of gene expression diversity during domestication of animals and plants

BACKGROUND

The genetic mechanisms underlying the domestication of animals and plants have been of great interest to biologists since Darwin. To date, little is known about the global pattern of gene expression changes during domestication.

RESULTS

We generated and collected transcriptome data for seven pairs of domestic animals and plants including dog, silkworm, chicken, rice, cotton, soybean and maize and their wild progenitors and compared the expression profiles between the domestic and wild species. Intriguingly, although the number of expressed genes varied little, the domestic species generally exhibited lower gene expression diversity than did the wild species, and this lower diversity was observed for both domestic plants and different kinds of domestic animals including insect, bird and mammal in the whole-genome gene set (WGGS), candidate selected gene set (CSGS) and non-CSGS, with CSGS exhibiting a higher degree of decreased expression diversity. Moreover, different from previous reports which found 2 to 4% of genes were selected by human, we identified 6892 candidate selected genes accounting for 7.57% of the whole-genome genes in rice and revealed that fewer than 8% of the whole-genome genes had been affected by domestication.

CONCLUSIONS

Our results showed that domestication affected the pattern of variation in gene expression throughout the genome and generally decreased the expression diversity across species, and this decrease may have been associated with decreased genetic diversity. This pattern might have profound effects on the phenotypic and physiological changes of domestic animals and plants and provide insights into the genetic mechanisms at the transcriptome level other than decreased genetic diversity and increased linkage disequilibrium underpinning artificial selection.

Evaluation of Reference Genes for Normalization of RT-qPCR Gene Expression Data for Trichoplusia ni Cells During Antheraea pernyi (Lepidoptera: Saturniidae) Multicapsid Nucleopolyhedrovirus (AnpeNPV) Infection

Baculovirus infection impacts global gene expression in the host cell, including the expression of housekeeping genes. Evaluation of candidate reference genes during a viral infection will inform the selection of appropriate reference gene(s) for the normalization of expression data generated by Reverse Transcription Quantitative Real-timePolymerase Chain Reaction (RT-qPCR). Antheraea pernyi multicapsid nucleopolyhedrovirus (AnpeNPV) is able to infect the High Five cells (Tn-Hi5). In the present study, 10 candidate reference genes were evaluated in AnpeNPV-infected Tn-Hi5 cells. Gene expression data were analyzed using geNorm, NormFinder, BestKeeper, and RefFinder. The candidate genes were further validated as reliable reference genes for RT-qPCR by analyzing the expression of three target genes. The results of data analysis using four statistical methods showed that RPS18 was the least stable among the candidate reference genes tested. 18S rRNA and 28S rRNA were not suitable as reference genes for RT-qPCR analysis in AnpeNPV-infected Tn-Hi5 cells. Comprehensive ranking of the 10 candidate reference genes by RefFinder analysis indicated that Ann B, c45128_g1, and ACT were the top three genes. Normalization of the expression of three target genes using the candidate reference genes indicated the same expression pattern when Ann B and c45128_g1 were used as reference genes, with slight differences in the relative expression at each infection time point. Overall, Ann B and c45128_g1 were recommended to be more suitable than the most commonly used reference genes, such as ACT, GAPDH, and TUB, for RT-qPCR data normalization in AnpeNPV-infected Tn-Hi5 cells up to 48 hpi.

Identification of Genes that Control Silk Yield by RNA Sequencing Analysis of Silkworm (Bombyx mori) Strains of Variable Silk Yield

Silk is an important natural fiber of high economic value, and thus genetic study of the silkworm is a major area of research. Transcriptome analysis can provide guidance for genetic studies of silk yield traits. In this study, we performed a transcriptome comparison using multiple silkworms with different silk yields. A total of 22 common differentially expressed genes (DEGs) were identified in multiple strains and were mainly involved in metabolic pathways. Among these, seven significant common DEGs were verified by quantitative reverse transcription polymerase chain reaction, and the results coincided with the findings generated by RNA sequencing. Association analysis showed that BGIBMGA003330 and BGIBMGA005780 are significantly associated with cocoon shell weight and encode uridine nucleosidase and small heat shock protein, respectively. Functional annotation of these genes suggest that these play a role in silkworm silk gland development or silk protein synthesis. In addition, we performed principal component analysis (PCA) in combination with wild silkworm analysis, which indicates that modern breeding has a stronger selection effect on silk yield traits than domestication, and imply that silkworm breeding induces aggregation of genes related to silk yield.

Prognostic value of sorting nexin 10 weak expression in stomach adenocarcinoma revealed by weighted gene co-expression network analysis

AIM

To detect significant clusters of co-expressed genes associated with tumorigenesis that might help to predict stomach adenocarcinoma (SA) prognosis.

METHODS

The Cancer Genome Atlas database was used to obtain RNA sequences as well as complete clinical data of SA and adjacent normal tissues from patients. Weighted gene co-expression network analysis (WGCNA) was used to investigate the meaningful module along with hub genes. Expression of hub genes was analyzed in 362 paraffin-embedded SA biopsy tissues by immunohistochemical staining. Patients were classified into two groups (according to expression of hub genes): Weak expression and over-expression groups. Correlation of biomarkers with clinicopathological factors indicated patient survival.

RESULTS

Whole genome expression level screening identified 6,231 differentially expressed genes. Twenty-four co-expressed gene modules were identified using WGCNA. Pearson’s correlation analysis showed that the tan module was the most relevant to tumor stage (r = 0.24, P = 7 × 10^-6). In addition, we detected sorting nexin (SNX)10 as the hub gene of the tan module. SNX10 expression was linked to T category (P = 0.042, χ² = 8.708), N category (P = 0.000, χ² = 18.778), TNM stage (P = 0.001, χ² = 16.744) as well as tumor differentiation (P = 0.000, χ² = 251.930). Patients with high SNX10 expression tended to have longer disease-free survival (DFS; 44.97 mo vs 33.85 mo, P = 0.000) as well as overall survival (OS; 49.95 vs 40.84 mo, P = 0.000) in univariate analysis. Multivariate analysis showed that dismal prognosis could be precisely predicted clinicopathologically using SNX10 [DFS: P = 0.014, hazard ratio (HR) = 0.698, 95% confidence interval (CI): 0.524-0.930, OS: P = 0.017, HR = 0.704, 95%CI: 0.528-0.940].

CONCLUSION

This study provides a new technique for screening prognostic biomarkers of SA. Weak expression of SNX10 is linked to poor prognosis, and is a suitable prognostic biomarker of SA.

Evidence of peripheral olfactory impairment in the domestic silkworms: insight from the comparative transcriptome and population genetics

BACKGROUND

The insect olfactory system is a highly specific and sensitive chemical detector, which plays important roles in feeding, mating and finding an appropriate oviposition site. The ecological niche of Bombyx mori has changed greatly since domestication from B. mandarina, and its olfactory response to environmental odorants clearly decreased. However, the mechanisms that result in the olfactory impairment are largely unknown.

RESULTS

The antennal transcriptomes were compared between the domestic and wild silkworms. Comparison of the same sex between the domestic and wild silkworms revealed 1410 and 1173 differentially expressed genes (DEGs) in males and females, respectively. To understand the olfactory impairment, we mainly focused on the olfactory-related genes. In total, 30 olfactory genes and 19 odorant-degrading enzymes (ODEs) showed differential expression in the two comparisons, in which 19 and 14 were down-regulated in the domestic silkworm, respectively. Based on population genomic data, the down-regulated odorant receptors (ORs) showed a higher ratio of unique non-synonymous polymorphisms to synonymous polymorphisms (N/S ratio) in the domestic populations than that in the wild silkworms. Furthermore, one deleterious mutation was found in OR30 of the domestic population, which was located in transmembrane helix 6 (TM6).

CONCLUSIONS

Our results suggested that down-regulation of the olfactory-related genes and relaxed selection might be the major reasons for olfactory impairment of the domestic silkworm reared completely indoor environment. Reversely, wild silkworm may increase expression and remove deleterious polymorphisms of olfactory-related genes to retain sensitive olfaction.

Identification and comparison of long non-coding RNAs in the silk gland between domestic and wild silkworms

Under long-term artificial selection, the domestic silkworm (Bombyx mori) has increased its silk yield tremendously in comparison with its wild progenitor, Bombyx mandarina. However, the molecular mechanism of silk yield increase is still unknown. Comparative analysis of long non-coding RNAs (lncRNAs) may provide some insights into understanding this phenotypic variation. In this study, using RNA sequencing technology data of silk gland in domestic and wild silkworms, we identified 599 lncRNAs in the silk gland of the silkworm. Compared with protein-coding genes, the silk gland lncRNA genes tend to have fewer exon numbers, shorter transcript length and lower GC-content. Moreover, we found that three lncRNA genes are significantly and differentially expressed between domestic and wild silkworms. The potential targets of two differentially expressed lncRNAs (DELs) (dw4sg_0040 and dw4sg_0483) and the expression-correlated genes with the two DELs are mainly enriched in the related processes of silk protein translation. This implies that these DELs may affect the phenotypic variation in silk yield between the domestic and wild silkworms through the post-transcriptional regulation of silk protein.

Biomaterials for the Delivery of Growth Factors and Other Therapeutic Agents in Tissue Engineering Approaches to Bone Regeneration

Bone fracture followed by delayed or non-union typically requires bone graft intervention. Autologous bone grafts remain the clinical "gold standard". Recently, synthetic bone grafts such as Medtronic's Infuse Bone Graft have opened the possibility to pharmacological and tissue engineering strategies to bone repair following fracture. This clinically-available strategy uses an absorbable collagen sponge as a carrier material for recombinant human bone morphogenetic protein 2 (rhBMP-2) and a similar strategy has been employed by Stryker with BMP-7, also known as osteogenic protein-1 (OP-1). A key advantage to this approach is its "off-the-shelf" nature, but there are clear drawbacks to these products such as edema, inflammation, and ectopic bone growth. While there are clinical challenges associated with a lack of controlled release of rhBMP-2 and OP-1, these are among the first clinical examples to wed understanding of biological principles with biochemical production of proteins and pharmacological principles to promote tissue regeneration (known as regenerative pharmacology). After considering the clinical challenges with such synthetic bone grafts, this review considers the various biomaterial carriers under investigation to promote bone regeneration. This is followed by a survey of the literature where various pharmacological approaches and molecular targets are considered as future strategies to promote more rapid and mature bone regeneration. From the review, it should be clear that pharmacological understanding is a key aspect to developing these strategies.

Intrinsic antimicrobial properties of silk spun by genetically modified silkworm strains

The domesticated silkworm, Bombyx mori, is a fundamental insect for silk industry. Silk is obtained from cocoons, protective envelopes produced during pupation and composed of single raw silk filaments secreted by the insect silk glands. Currently, silk is used as a textile fibre and to produce new materials for technical and biomedical applications. To enhance the use of both fabrics and silk-based materials, great efforts to obtain silk with antimicrobial properties have been made. In particular, a convincing approach is represented by the enrichment of the textile fibre with antimicrobial peptides, the main effectors of the innate immunity. To this aim, silkworm-based transgenic techniques appear to be cost-effective strategies to obtain cocoons in which antimicrobial peptides are integrated among the silk proteins. Recently, cocoons transgenic for a recombinant silk protein conjugated to the silkworm Cecropin B antimicrobial peptide were obtained and showed enhanced antibacterial properties (Li et al. in Mol Biol Rep 42:19-25, https://doi.org/10.1007/s11033-014-3735-z , 2015a). In this work we used the piggyBac-mediated germline transformation to generate several transgenic B. mori lines able to overexpress Cecropin B or Moricin antimicrobial peptides at the level of the silk gland. The derived cocoons were characterised by increased antimicrobial properties and the resulting silk fibre was able to inhibit the bacterial growth of the Gram-negative Escherichia coli. Our results suggest that the generation of silkworm overexpressing unconjugated antimicrobial peptides in the silk gland might represent an additional strategy to obtain antimicrobial peptide-enriched silk, for the production of new silk-based materials.

Molecular Network-Based Identification of Competing Endogenous RNAs in Thyroid Carcinoma

RNAs may act as competing endogenous RNAs (ceRNAs), a critical mechanism in determining gene expression regulations in many cancers. However, the roles of ceRNAs in thyroid carcinoma remains elusive. In this study, we have developed a novel pipeline called Molecular Network-based Identification of ceRNA (MNIceRNA) to identify ceRNAs in thyroid carcinoma. MNIceRNA first constructs micro RNA (miRNA)-messenger RNA (mRNA)long non-coding RNA (lncRNA) networks from miRcode database and weighted correlation network analysis (WGCNA), based on which to identify key drivers of differentially expressed RNAs between normal and tumor samples. It then infers ceRNAs of the identified key drivers using the long non-coding competing endogenous database (lnCeDB). We applied the pipeline into The Cancer Genome Atlas (TCGA) thyroid carcinoma data. As a result, 598 lncRNAs, 1025 mRNAs, and 90 microRNA (miRNAs) were inferred to be differentially expressed between normal and thyroid cancer samples. We then obtained eight key driver miRNAs, among which hsa-mir-221 and hsa-mir-222 were key driver RNAs identified by both miRNA-mRNA-lncRNA and WGCNA network. In addition, hsa-mir-375 was inferred to be significant for patients' survival with 34 associated ceRNAs, among which RUNX2, DUSP6 and SEMA3D are known oncogenes regulating cellular proliferation and differentiation in thyroid cancer. These ceRNAs are critical in revealing the secrets behind thyroid cancer progression and may serve as future therapeutic biomarkers.

Integrated study on comparative transcriptome and skeletal muscle function in aged rats

The present study aimed to reveal aging-related changes in the skeletal muscle of SD rats by comparing transcriptome analysis, integrated with muscle physiological parameters. Ten rats aged 25 months were set as the old group (OG) and ten rats aged 6 months were set as the young group (YG). After 6 weeks of feeding, the body mass, grip strength, and gastrocnemius muscle mass were determined, and the differentially expressed genes were analyzed by transcriptome sequencing, followed by GO enrichment analysis and KEGG analysis. The results showed that the muscle index and the relative grip strength were lower in OG rats than YG rats. The expressions of AMPK, UCP3, IGF-1, several ion channel associated genes and collagen family genes were down-regulated in OG rats. MGMT, one of the strength determining genes and CHRNa1, a subunit of the acetylcholine receptor were up-regulated in OG rats. The present results supply the global transcriptomic information involved in aging related skeletal muscle dysfunction in rats. The reduced expressions of AMPK, IGF-1, and CASK can explain the losses of muscle mass and function in the aged rats. In addition, the up-regulation of MGMT and CHRNa1 also contribute to muscle wasting and weakness during aging.

QTL analysis of cocoon shell weight identifies BmRPL18 associated with silk protein synthesis in silkworm by pooling sequencing

Mechanisms that regulate silk protein synthesis provide the basis for silkworm variety breeding and silk gland bioreactor optimization. Here, using the pooling sequencing-based methodology, we deciphered the genetic basis for the varied silk production in different silkworm strains. We identified 8 SNPs, with 6 on chromosome 11 and 1 each on chromosomes 22 and 23, that were linked with silk production. After conducting an association analysis between gene expression pattern, silk gland development and cocoon shell weight (CSW), BMGN011620 was found to be regulating silk production. BMGN011620 encodes the 60S ribosomal protein, L18, which is an indispensable component of the 60S ribosomal subunit; therefore we named it BmRPL18. Moreover, the clustering of linked SNPs on chromosome 11 and the analysis of differentially expressed genes reported in previous Omics studies indicated that the genes regulating silk protein synthesis may exhibit a clustering distribution in the silkworm genome. These results collectively advance our understanding of the regulation of silk production, including the role of ribosomal proteins and the clustered distribution of genes involved in silk protein synthesis.

Effects of salinity on the cellular physiological responses of Natrinema sp. J7-2

The halophilic archaea (haloarchaea) live in hyersaline environments such as salt lakes, salt ponds and marine salterns. To cope with the salt stress conditions, haloarchaea have developed two fundamentally different strategies: the "salt-in" strategy and the "compatible-solute" strategy. Although investigation of the molecular mechanisms underlying the tolerance to high salt concentrations has made outstanding achievements, experimental study from the aspect of transcription is rare. In the present study, we monitored cellular physiology of Natrinema sp. J7-2 cells incubated in different salinity media (15%, 25% and 30% NaCl) from several aspects, such as cellular morphology, growth, global transcriptome and the content of intracellular free amino acids. The results showed that the cells were polymorphic and fragile at a low salt concentration (15% NaCl) but had a long, slender rod shape at high salt concentrations (25% and 30% NaCl). The cells grew best in 25% NaCl, mediocre in 30% NaCl and struggled in 15% NaCl. An RNA-seq analysis revealed differentially expressed genes (DEGs) in various salinity media. A total of 1,148 genes were differentially expressed, consisting of 719 DEGs (348 up-regulated and 371 down-regulated genes) between cells in 15% vs 25% NaCl, and 733 DEGs (521 up-regulated and 212 down-regulated genes) between cells in 25% vs 30% NaCl. Moreover, 304 genes were commonly differentially expressed in both 15% vs 25% and 25% vs30% NaCl. The DEGs were enriched in different KEGG metabolic pathways, such as amino acids, glycerolipid, ribosome, nitrogen, protoporphyrin, porphyrin and porhiniods. The intracellular predominant free amino acids consisted of the glutamate family (Glu, Arg and Pro), aspartate family (Asp) and aromatic amino acids (Phe and Trp), especially Glu and Asp.

Transcriptome Analysis of the Sm-Mediated Hypersensitive Response to Stemphylium lycopersici in Tomato

Gray leaf spot disease caused by Stemphylium lycopersici is a major disease in cultivated tomato plants and threatens tomato-growing areas worldwide. Sm is a single dominant gene that confers resistance to tomato gray leaf spot disease agent. However, the underlying molecular mechanism remains unclear. Here, resistant (cv. Motelle, containing the Sm gene) and susceptible (cv. Moneymaker) plants were inoculated with virulent Stemphylium lycopersici isolate at a time point at which both cultivars showed a strong response to S. lycopersici infection. Transcriptome analyses were performed in both cultivars using RNA-seq. The number of differentially expressed genes (DEGs) was higher in Motelle than Moneymaker. Functional classification revealed that most DEGs were involved in plant-pathogen interactions, plant hormone signal transduction, regulation of autophagy, glycerophospholipid metabolism, and α-linolenic acid metabolism. Moreover, the genes that were significantly up-regulated in Sm tomatoes were involved in plant-pathogen interaction pathways. A total of 26 genes were selected for confirmation of differentially expressed levels by quantitative real-time PCR. This knowledge will yield new insights into the molecular mechanism of Sm responses to S. lycopersici infection.

Segmental duplications: evolution and impact among the current Lepidoptera genomes

Background

Structural variation among genomes is now viewed to be as important as single nucleoid polymorphisms in influencing the phenotype and evolution of a species. Segmental duplication (SD) is defined as segments of DNA with homologous sequence.

Results

Here, we performed a systematic analysis of segmental duplications (SDs) among five lepidopteran reference genomes (Plutella xylostella, Danaus plexippus, Bombyx mori, Manduca sexta and Heliconius melpomene) to understand their potential impact on the evolution of these species. We find that the SDs content differed substantially among species, ranging from 1.2% of the genome in B. mori to 15.2% in H. melpomene. Most SDs formed very high identity (similarity higher than 90%) blocks but had very few large blocks. Comparative analysis showed that most of the SDs arose after the divergence of each linage and we found that P. xylostella and H. melpomene showed more duplications than other species, suggesting they might be able to tolerate extensive levels of variation in their genomes. Conserved ancestral and species specific SD events were assessed, revealing multiple examples of the gain, loss or maintenance of SDs over time. SDs content analysis showed that most of the genes embedded in SDs regions belonged to species-specific SDs (“Unique” SDs). Functional analysis of these genes suggested their potential roles in the lineage-specific evolution. SDs and flanking regions often contained transposable elements (TEs) and this association suggested some involvement in SDs formation. Further studies on comparison of gene expression level between SDs and non-SDs showed that the expression level of genes embedded in SDs was significantly lower, suggesting that structure changes in the genomes are involved in gene expression differences in species.

Conclusions

The results showed that most of the SDs were “unique SDs”, which originated after species formation. Functional analysis suggested that SDs might play different roles in different species. Our results provide a valuable resource beyond the genetic mutation to explore the genome structure for future Lepidoptera research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1007-y) contains supplementary material, which is available to authorized users.

Comparative Analysis of Transcriptomes of Macrophage Revealing the Mechanism of the Immunoregulatory Activities of a Novel Polysaccharide Isolated from Boletus speciosus Frost

BACKGROUND

The mechanism of the immunoregulatory activities of polysaccharide is still not clear.

MATERIALS AND METHODS

Here, we performed the B-cell, T-cell, and macrophage cell proliferation, the cell cycle analysis of macrophage cells, sequenced the transcriptomes of control group macrophages, and Boletus speciosus Frost polysaccharide (BSF-1) group macrophages using Illumina sequencing technology to identify differentially expressed genes (DEGs) to determine the molecular mechanisms of immunomodulatory activity of BSF-1 in macrophages.

RESULTS

These results suggested that BSF-1 could promote the proliferation of B-cell, T-cell, and macrophages, promote the proliferation of macrophage cells by abolishing cell cycle arrests in the G0/G1 phases, and promote cell cycle progression in S-phase and G2/M phase, which might induce cell division. A total of 12,498,414 and 11,840,624 bp paired-end reads were obtained for the control group and BSF-1 group, respectively, and they corresponded to a total size of 12.5 G bp and 11.8 G bp, respectively, after the low-quality reads and adapter sequences were removed. Approximately 81.83% of the total number of genes (8,257) were expressed reads per kilobase per million mapped reads (RPKM ≥1) and more than 1366 genes were highly expressed (RPKM >60) in the BSF-1 group. A gene ontology-enrichment analysis generated 13,042 assignments to cellular components, 13,094 assignments to biological processes, and 13,135 assignments to molecular functions. A Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis showed that the mitogen-activated protein kinase (MAPK) signaling pathways are significantly enriched for DEGs between the two cell groups.

CONCLUSION

An analysis of transcriptome resources enabled us to examine gene expression profiles, verify differential gene expression, and select candidate signaling pathways as the mechanisms of the immunomodulatory activity of BSF-1. Based on the experimental data, we believe that the significant antitumor activities of BSF-1 in vivo mainly involve the MAPK signaling pathways.

SUMMARY

Boletus speciosus Frost-1 (BSF-1) could promote the proliferation of B-cell, T-cell, and macrophages, promote the proliferation of macrophage cells by abolishing cell cycle arrests in the G0/G1 phases, and promote cell cycle progression in S-phase and G2/M phase, which might induce cell divisionApproximately 81.83% of the total number of genes (8257) were expressed (reads per kilobase per million mapped reads [RPKM] =1) and more than 1366 genes were highly expressed (RPKM >60) in the BSF-1 groupA gene ontology-enrichment analysis generated 13,042 assignments to cellular components, 13,094 assignments to biological processes, and 13,135 assignments to molecular functionsA Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis showed that the mitogen-activated protein kinase signaling pathways are significantly enriched for DEGs between the two cell groups. Abbreviations used: BSF-1: Boletus speciosus Frost polysaccharide.

Transcriptome Analysis of the Cf-12-Mediated Resistance Response to Cladosporium fulvum in Tomato

Cf-12 is an effective gene for resisting tomato leaf mold disease caused by Cladosporium fulvum (C. fulvum). Unlike many other Cf genes such as Cf-2, Cf-4, Cf-5, and Cf-9, no physiological races of C. fulvum that are virulent to Cf-12 carrying plant lines have been identified. In order to better understand the molecular mechanism of Cf-12 gene resistance response, RNA-Seq was used to analyze the transcriptome changes at three different stages of C. fulvum infection (0, 4, and 8 days post infection [dpi]). A total of 9100 differentially expressed genes (DEGs) between 4 and 0 dpi, 8643 DEGs between 8 and 0 dpi and 2547 DEGs between 8 and 4 dpi were identified. In addition, we found that 736 DEGs shared among the above three groups, suggesting the presence of a common core of DEGs in response to C. fulvum infection. These DEGs were significantly enriched in defense-signaling pathways such as the calcium dependent protein kinases pathway and the jasmonic acid signaling pathway. Additionally, we found that many transcription factor genes were among the DEGs, indicating that transcription factors play an important role in C. fulvum defense response. Our study provides new insight on the molecular mechanism of Cf resistance to C. fulvum, especially the unique features of Cf-12 in responding to C. fulvum infection.

Transcriptome response to elevated atmospheric CO2 concentration in the Formosan subterranean termite, Coptotermes formosanus Shiraki (Isoptera: Rhinotermitidae)

Background

Carbon dioxide (CO₂) is a pervasive chemical stimulus that plays a critical role in insect life, eliciting behavioral and physiological responses across different species. High CO₂ concentration is a major feature of termite nests, which may be used as a cue for locating their nests. Termites also survive under an elevated CO₂ concentration. However, the mechanism by which elevated CO₂ concentration influences gene expression in termites is poorly understood.

Methods

To gain a better understanding of the molecular basis involved in the adaptation to CO₂ concentration, a transcriptome of Coptotermes formosanus Shiraki was constructed to assemble the reference genes, followed by comparative transcriptomic analyses across different CO₂ concentration (0.04%, 0.4%, 4% and 40%) treatments.

Results

(1) Based on a high throughput sequencing platform, we obtained approximately 20 GB of clean data and revealed 189,421 unigenes, with a mean length and an N50 length of 629 bp and 974 bp, respectively. (2) The transcriptomic response of C. formosanus to elevated CO₂ levels presented discontinuous changes. Comparative analysis of the transcriptomes revealed 2,936 genes regulated among 0.04%, 0.4%, 4% and 40% CO₂ concentration treatments, 909 genes derived from termites and 2,027 from gut symbionts. Genes derived from termites appears selectively activated under 4% CO₂ level. In 40% CO₂ level, most of the down-regulated genes were derived from symbionts. (3) Through similarity searches to data from other species, a number of protein sequences putatively involved in chemosensory reception were identified and characterized in C. formosanus, including odorant receptors, gustatory receptors, ionotropic receptors, odorant binding proteins, and chemosensory proteins.

Discussion

We found that most genes associated with carbohydrate metabolism, energy metabolism, and genetic information processing were regulated under different CO₂ concentrations. Results suggested that termites adapt to ∼4% CO₂ level and their gut symbionts may be killed under high CO₂ level. We anticipate that our findings provide insights into the transcriptome dynamics of CO₂ responses in termites and form the basis to gain a better understanding of regulatory networks.

BmncRNAdb: a comprehensive database of non-coding RNAs in the silkworm, Bombyx mori

BACKGROUND

Long non-coding RNAs (lncRNAs) may play critical roles in a wide range of developmental processes of higher organisms. Recently, lncRNAs have been widely identified across eukaryotes and many databases of lncRNAs have been developed for human, mouse, fruit fly, etc. However, there is rare information about them in the only completely domesticated insect, silkworm (Bombyx mori).

DESCRIPTION

In this study, we systematically scanned lncRNAs using the available silkworm RNA-seq data and public unigenes. Finally, we identified and collected 6281 lncRNAs in the silkworm. Besides, we also collected 1986 microRNAs (miRNAs) from previous studies. Then, we organized them into a comprehensive and web-based database, BmncRNAdb. This database offers a user-friendly interface for data browse and online analysis as well as the three online tools for users to predict the target genes of lncRNA or miRNA.

CONCLUSIONS

We have systematically identified and collected the silkworm lncRNAs and constructed a comprehensive database of the silkworm lncRNAs and miRNAs. This work gives a glimpse into lncRNAs of the silkworm and lays foundations for the ncRNAs study of the silkworm and other insects in the future. The BmncRNAdb is freely available at http://gene.cqu.edu.cn/BmncRNAdb/index.php .