Search

Abstract Background Insect-based food and feed are increasingly attracting attention. As a domesticated insect, the silkworm (Bombyx mori) has a highly nutritious pupa that can be easily raised in large quantities through large-scale farming, making it a highly promising source of food. The ratio of pupa to cocoon (RPC) refers to the proportion of the weight of the cocoon that is attributed to pupae, and is of significant value for edible utilization, as a higher RPC means a higher ratio of conversion of mulberry leaves to pupa. In silkworm production, there is a trade-off between RPC and cocoon shell ratiao(CSR), which refers the ratio of silk protein to the entire cocoon, during metamorphosis process. Understanding the genetic basis of this balance is crucial for breeding edible strains with a high RPC and further advancing its use as feed. Results Using QTL-seq, we identified a quantitative trait locus (QTL) for the balance between RPC and CSR that is located on chromosome 11 and covers a 9,773,115-bp region. This locus is an artificial selection hot spot that contains ten non-overlapping genomic regions under selection that were involved in the domestication and genetic breeding processes. These regions include 17 genes, nine of which are highly expressed in the silk gland, which is a vital component in the trade-off between RPC and CSR. These genes are annotate with function related with epigenetic modifications and the regulation of DNA replication et al. We identified one and two single nucleotide polymorphisms (SNPs) in the exons of teh KWMTBOMO06541 and KWMTBOMO06485 genes that result in amino acid changes in the protein domains. These SNPs have been strongly selected for during the domestication process. The KWMTBOMO06485 gene encodes the Bombyx mori (Bm) tRNA methyltransferase (BmDnmt2) and its knockout results in a significant change in the trade-off between CSR and RPC in both sexes. Conclusions Taken together, our results contribute to a better understanding of the genetic basis of RPC and CSR. The identified QTL and genes that affect RPC can be used for marker-assisted and genomic selection of silkworm strains with a high RPC. This will further enhance the production efficiency of silkworms and of closely-related insects for edible and feed purposes.

Abstract FOXO transcription factors modulate aging-related pathways and influence longevity in multiple species, but the transcriptional targets that mediate these effects remain largely unknown. Here, we identify an evolutionarily conserved FOXO target gene, Oxidative stress-responsive serine-rich protein 1 (OSER1), whose overexpression extends lifespan in silkworms, nematodes, and flies, while its depletion correspondingly shortens lifespan. In flies, overexpression of OSER1 increases resistance to oxidative stress, starvation, and heat shock, while OSER1-depleted flies are more vulnerable to these stressors. In silkworms, hydrogen peroxide both induces and is scavenged by OSER1 in vitro and in vivo. Knockdown of OSER1 in Caenorhabditis elegans leads to increased ROS production and shorter lifespan, mitochondrial fragmentation, decreased ATP production, and altered transcription of mitochondrial genes. Human proteomic analysis suggests that OSER1 plays roles in oxidative stress response, cellular senescence, and reproduction, which is consistent with the data and suggests that OSER1 could play a role in fertility in silkworms and nematodes. Human studies demonstrate that polymorphic variants in OSER1 are associated with human longevity. In summary, OSER1 is an evolutionarily conserved FOXO-regulated protein that improves resistance to oxidative stress, maintains mitochondrial functional integrity, and increases lifespan in multiple species. Additional studies will clarify the role of OSER1 as a critical effector of healthy aging.

Cuproptosis is a newly identified form of copper-dependent cell death that differs from other known pathways. This discovery provides a new way to explore copper-based nanomaterial applications in cancer therapy. This study used a layer-by-layer self-assembling method to load Cu2-xS nanoparticle (NP) cores with the siRNA of the PD-L1 immune escape-related gene and wrap a silk fibroin (SF) shell to form a multifunctional copper-based SF nanoplatform, denoted as CuS-PEI-siRNA-SFNs. CuS-PEI-siRNA-SFNs induced cuproptosis and exerted an antitumor effect via multiple mechanisms, including photothermal therapy (PTT), chemodynamic therapy (CDT), and immune activation. The presence of significant dihydrolipoamide S-acetyltransferase (DLAT) oligomers in 4T1 cells treated with CuS-PEI-siRNA-SFNs indicated the triggering of cuproptosis. Furthermore, in vivo experimental results showed that CuS-PEI-siRNA-SFNs efficiently accumulated in the tumor tissues of 4T1 tumor-bearing mice inhibited primary tumor and lung metastasis, and displayed excellent biosafety and antitumor activity. This study demonstrated that the synergistic effect of cuproptosis, PTT, CDT, and immune activation showed promise for treating metastatic breast cancer.

The color pattern of insects is one of the most diverse adaptive evolutionary phenotypes. However, the molecular regulation of this color pattern is not fully understood. In this study, we found that the transcription factor Bm-mamo is responsible for black dilute (bd) allele mutations in the silkworm. Bm-mamo belongs to the BTB zinc finger family and is orthologous to mamo in Drosophila melanogaster. This gene has a conserved function in gamete production in Drosophila and silkworms and has evolved a pleiotropic function in the regulation of color patterns in caterpillars. Using RNAi and clustered regularly interspaced short palindromic repeats (CRISPR) technology, we showed that Bm-mamo is a repressor of dark melanin patterns in the larval epidermis. Using in vitro binding assays and gene expression profiling in wild-type and mutant larvae, we also showed that Bm-mamo likely regulates the expression of related pigment synthesis and cuticular protein genes in a coordinated manner to mediate its role in color pattern formation. This mechanism is consistent with the dual role of this transcription factor in regulating both the structure and shape of the cuticle and the pigments that are embedded within it. This study provides new insight into the regulation of color patterns as well as into the construction of more complex epidermal features in some insects.

Abstract Alpha klotho (α-Klotho) is an anti-aging molecule associated with aging and several diseases. Previous studies have reported inconsistent levels of α-Klotho in smokers. This study aimed to demonstrate serum α-Klotho levels in smokers among the US population. This cross-sectional study recruited 11,559 participants (aged 40–79 years; 48.2% males). All data were collected from the 2007–2016 National Health and Nutrition Examination Survey. The study comprised adults with reliable Klotho and smoking questionnaire results. The relationship between smoking and serum α-klotho levels was assessed using multivariate linear regression models after adjusting for potential confounders. We also performed a stratified analysis of clinically important variables. The mean serum α-klotho level among the 11,559 participants was 843.85 pg/mL. After full adjustment, habitual smoking was significantly associated with decreased serum levels of α-klotho level (β = − 34.89; 95% CI − 54.97, − 14.81; P = 0.0013) in the total study population. Furthermore, the stratified analysis indicated that the association was insignificant in the 60–79 age group. Quitting smoking was not significantly associated with serum levels of α-klotho as expected (P = 0.1148) in the total study population. However, stratified analyses showed a significant inversed association in the male, those with chronic kidney disease, or those with cancer who quit smoking (all P

Abstract Background Parkinson’s disease (PD) is the second most common neurodegenerative disease in middle-aged and elderly populations, whereas there is no cure for PD so far. Novel animal models and medications await development to elucidate the aetiology of PD and attenuate the symptoms, respectively. Methods A neurotoxin, 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP), was used in the current study to establish a PD pathologic model in silkworms. The time required to complete specific behaviours was recorded. Dopamine content was detected by ultra-performance liquid chromatography (UPLC). The activity of insect tyrosine hydroxylase (TH) was determined using a double-antibody sandwich method. Oxidative stress was assessed by changes in antioxidant enzyme activity and the content of oxidative products. Results MPTP-treated silkworms were characterized by impaired motor ability, reduced dopamine content, and elevated oxidative stress level. The expression of TH, a dopamine biosynthetic enzyme within dopaminergic neurons in the brain, was significantly reduced, indicating that dopaminergic neurons were damaged. Moreover, MPTP-induced motility impairment and reduced dopamine level in the silkworm PD model could be rescued after feeding a combination of levodopa (L-dopa [LD]) and carbidopa (CD). MPTP-induced oxidative damage was also alleviated, in ways consistent with other PD animal models. Interestingly, administration of Lycium barbarum polysaccharide (LBP) improved the motor ability, dopamine level, and TH activity, and the oxidative damage was concomitantly reduced in the silkworm PD model. Conclusions This study provides a promising animal model for elucidating the pathogenesis of PD, as well as a relevant preliminary drug screening (e.g., LBP) and evaluation.

Tong et al. describe a super pangenome assembled from long-read sequences of 545 wild and domesticated silkworms. Naturally selected (diapause, aposemantic coloration) or artificially selected (silk yield and fineness) sets of genes are delineated.

Adjuvant gefitinib improves outcomes in non-small cell lung cancer (NSCLC) patients compared to chemotherapy, but not in all cases. Here, the authors find genomic biomarkers of response to gefitinib in NSCLC patients from the ADJUVANT trial, and propose a score to stratify them by potential benefit from the treatment.

Clinical wound management is always a relatively urgent problem. Moreover, wounds, especially severe wounds with excessive tension or excessive movement are prone to tissue infection, necrosis, and other negative effects during healing. Therefore, research has aimed to develop low-cost complementary treatments to address the urgent need for an innovative low-cost dressing that can adapt to high mechanical requirements and complex wound conditions. At present, tissue engineering to produce artificial skin with a structure similar to that of normal skin is one effective method to solve this challenge in the regeneration and repair of serious wounds. The present study hot pressed flat silk cocoons (FSC) with carboxymethyl chitosan (CMCS) to generate a cross-linked binding without enzymes or cross-linking agents that simulated the 3D structural composites of the skin cuticle. This hybrid membrane showed potential to reduce inflammatory cells and promote neovascularization in skin wound repair. After hot pressing at 130°C and 20 Mpa, the FSC/CMCS composite material was denser than FSC, showed strong light transmission, and could be arbitrarily cut. Simulating the normal skin tissue structure, the hybrid membrane overcame the poor mechanical properties of traditional support materials. Moreover, the combination of protein and polysaccharide simulated the extracellular matrix, thus providing better biocompatibility. The results of this study also demonstrated the excellent mechanical properties of the FSC/CMCS composite support material, which also provided a low-cost and environmentally friendly process for making dressings. In addition, the results of this study preliminarily reveal the mechanism by which the scaffolds promoted the healing of full-thickness skin defects on the back of SD rats. In vivo experiments using a full-thickness skin defect model showed that the FSC/CMCS membranes significantly promoted the rate of wound healing and also showed good effects on blood vessel formation and reduced inflammatory reactions. This bionic support structure, with excellent repair efficacy on deep skin defect wounds, showed potential to further improve the available biomaterial systems, such as skin and other soft tissues.

Fiber optic distributed acoustic sensing (DAS) technology is widely used in security surveillance and geophysical survey applications. The response of the DAS system to external vibrations varies with different types of fiber optic cable connections. The mechanism of mutual influence between the cable’s characteristics and DAS measurement results remains unclear. This study proposed a dynamic model of the interaction between the optical cable and the soil, analyzed the impact of the dynamic parameters of the optical cable and soil on the sensitivity of the DAS system, and validated the theoretical analysis through experiments. The findings suggest that augmenting the cable’s bending stiffness 5.5-fold and increasing its unit mass 4.2-fold result in a discernible reduction of the system’s response to roughly 0.15 times of its initial magnitude. Cables with lower unit mass and bending stiffness are more sensitive to vibration signals. This research provides a foundation for optimizing vibration-enhanced fiber optic cables and broadening the potential usage scenarios for DAS systems.

Bioprospecting of more novel probiotic strains has attained continuous interest. This study aimed to investigate the beneficial effects of Lactobacillus paracasei strain L14, an isolate from a traditional Chinese dairy product, on type 2 diabetes mellitus (T2DM) rats. Preventive supplementation of strain L14 showed excellent anti-diabetic effects on high-fat diet/low-dose streptozotocin (HFD/STZ)-induced T2DM rats. It significantly reduced hyperglycemia, protected pancreatic β-cell and liver function, and ameliorated oxidative stress while considerably improving dyslipidemia and inflammation. Furthermore, the strain modulated the gut microbiota to alleviate gut dysbiosis. Interestingly, most of these biochemical parameters could even restore to normal levels by the intervention of strain L14. The whole-genome sequencing of L14 was performed to provide a critical molecular basis for its probiotic activities. Genes related to antioxidant systems and other beneficial microbial metabolites like exopolysaccharides (EPS) biosynthesis were found. This study demonstrates that probiotic L. paracasei L14 has good potential for applications in functional food and pharmaceutical industries.

A number of research has shown that the plant polyphenol resveratrol, one of the most prominent small molecules, has beneficial protective effects in multiple organisms, including worms, flies, and killifish. To understand the effects of resveratrol on lifespan, we evaluated its effects in the silkworm Bombyx mori. In this study, we found that lifespan was significantly prolonged in both female and male silkworms treated with resveratrol. Silkworm larval weight was significantly increased from day 3 of the 5th larval instar (L5D3) to day 7 of the 5th larval instar (L5D7). However, the weight of the pupa, cocoon, and total cocoon was not significantly different in female silkworms with resveratrol treatment than that in controls. Meanwhile, resveratrol significantly improved the thermotolerance of the silkworms, which enhanced their survival rate. Moreover, antioxidant activity was increased by resveratrol in both female and male silkworms. Furthermore, an antioxidant-related signalling pathway, SIRT7-FoxO-GST, was activated in silkworms with resveratrol treatment. Collectively, these results help us to understand the molecular pathways underlying resveratrol induced pro-longevity effects and indicate that silkworm is a promising animal model for evaluating the effects of lifespan-extending drugs.

Liu et al. analyse the protein composition of developing wing scales in the domestic silk moth. This molecular study of scales provides fundamental information about how such a fine microstructure is constructed and insights into the potential application of His-rich cuticular proteins as new biomaterials.

Bombyx mandarina silk fibers, as a precious resource of wild silk fibers, have some different characteristics to Bombyx mori silk fibers. To make a better use of this precious green material in various fields, it is necessary to investigate the characteristics of the fibers. In this research, the structure and properties of Bombyx mori silk fiber, Hybrid silk fiber, and Bombyx mandarina silk fiber are studied and analyzed. In terms of secondary structure, β-sheet structures are dominant in all three fibers while random coil and α-helix structures account for smaller portions of the fibers. Bombyx mandarina silk fiber exhibits a characteristic of amino acid distribution where the content of glycine is higher than that of alanine. The further analysis indicates the breaking strength of B. mandarina silk fibers is superior to that of B. mori silk fiber although Bombyx mandarina silk fiber is finer than Bombyx mori silk fiber. Our investigations provide a significant guide in the development and utilization of these precious silk resources, and Bombyx mandarina and Hybrid silk fibers show good potential for application in the same fields as Bombyx mori silk fibers.

ROS1 gene rearrangements have been reported in diverse cancer types including non‐small‐cell lung cancer (NSCLC), and with a notably higher prevalence in lung adenocarcinoma. The tyrosine kinase inhibitors, crizotinib, lorlatinib, and entrectinib, have demonstrated favorable efficacy in treating ROS1‐rearranged NSCLCs. Herein, we retrospectively reviewed 17 158 NSCLC patients whose tumor specimen and/or circulating cell‐free DNA underwent comprehensive genomic profiling. A total of 258 unique patients were identified with ROS1 rearrangements, representing an overall prevalence of approximately 1.5% of ROS1 fusions in newly diagnosed and relapsed NSCLC patients. CD74 (38%) was the most common fusion partner of ROS1, followed by EZR (13%), SDC4 (13%), SLC34A2 (10%), and other recurrent fusion partners with lower frequencies, including TPM3, MYH9, and CCDC6. Variant breakpoints occurred in ROS1 introns 33 (37%), 31 (25%), 32 (17%), and 34 (11%) with no obvious hotspots. CD74 (63%) and EZR (50%) were more frequently fused to ROS1 intron 33 than other introns, while ROS1 intron 31 was most frequently fused with SDC4 (79%) and SLC34A2 (81%). Crizotinib progression‐free survival (PFS) was not significantly different between fusion variants involving breakpoints in different ROS1 introns, nor was there a significant difference in PFS between CD74‐ROS1 and non‐CD74‐ROS1 groups of patients. Furthermore, TP53 was most frequently mutated in patients who progressed on crizotinib, and TP53 mutations were significantly associated with shorter crizotinib PFS. ROS1 mutations, including G2032R, were observed in approximately 33% of post‐crizotinib samples. Collectively, we report the prevalence of ROS1 fusions in a large‐scale NSCLC population and the efficacy of crizotinib in treating patients with ROS1‐rearranged NSCLC.

EGFR exon 20 insertions (EGFR e20ins) account for up to 10% of EGFR mutations in lung cancer; however, tumors with EGFR e20ins had poor response rates to EGFR tyrosine kinase inhibitors (TKIs) including gefitinib, erlotinib, afatinib, and osimertinib, and the heterogeneity of EGFR e20ins further complicates the clinical studies. Here, we retrospectively screened next‐generation sequencing (NGS) data from 24 468 lung cancer patients, and a total of 85 unique EGFR e20ins variants were identified in 547 cases (2.24%), with p.A767_V769dup (25.1%) and p.S768_D770dup (17.6%) being the most prevalent ones. Comprehensive genomic profiling revealed that TP53 mutations frequently coexisted with p.H773dup (77.8%, P = 0.0558) and p.A767_V769dup (62.8%, P = 0.0325), while RB1 mutations usually co‐occurred with p.H773_V774insAH (33.3%, P = 0.0551), implying that different EGFR e20ins variants might require distinct genomic context for tumorigenesis and/or maintenance. Despite that treatment regimens were highly diverse for EGFR e20ins‐positive patients, we observed an overall response rate of 14% and a disease control rate (DCR) of 38.4% in 65 patients who received at least one EGFR TKI. The progression‐free survival (PFS) differs significantly in six representative EGFR e20ins variants (P = 0.017), and EGFR p.A763_Y764insFQEA was associated with better PFS than other EGFR e20ins when treating with various EGFR TKIs. Some EGFR e20ins variants showed at least partial response to first‐generation EGFR TKIs, including p.A767_V769dup, p.S768_D770dup, p.N771_H773dup, p.A763_Y764insFQEA, and p.D770_N771insG. Poziotinib achieved higher DCR for p.S768_D770dup than for p.A767_V769dup, whereas osimertinib showed limited effects for these two insertions when used as the first‐line treatment. Overall, our results demonstrated that EGFR e20ins were highly diversified in terms of insertion patterns and co‐occurring mutations and these EGFR e20ins variants showed different clinical responses to various EGFR TKIs, suggesting the clinical importance of selecting proper EGFR TKI treatment based on the specific EGFR e20ins type.

Abstract Background Definitive chemoradiation therapy (dCRT) is the standard treatment for patients with nonsurgical esophageal squamous cell carcinoma (ESCC), yet patients have demonstrated great variations in their responses to dCRT and inevitably progressed following treatment. Methods To identify prognostic biomarkers, we performed targeted next‐generation sequencing of 416 cancer‐related genes on primary tumors from 47 nonsurgical ESCC patients prior to dCRT treatment. The association between genetic alterations and patients' local recurrence‐free survival (LRFS), progression‐free survival (PFS), and overall survival (OS) was analyzed. Results TP53 (78% of patients), NOTCH1 (32%), ARID1A (13%), FAT1 (13%), and CDKN2A (13%) were commonly mutated in ESCC patients, while gene amplifications frequently occurred in MCL1 (36%), FGF19 (34%), MYC (32%), CCND1 (27%), ZNF217 (15%), CDKN2A (13%), and YAP1 (11%). Univariate and multivariate analyses of clinical factors and genetic alterations indicated that sex is an independent prognostic factor, with males tending to have better LRFS (hazard ratio [HR], 0.25; 95%CI, 0.08‐0.77, P = .015) and progression‐free survival (PFS) (HR, 0.35; 95%CI, 0.13‐0.93, P = .030) following dCRT. Meanwhile, YAP1 amplification (n = 7) was an adverse prognostic factor, and patients with this alteration demonstrated a tendency toward worse outcomes with shorter LRFS (HR, 4.06; 95%CI, 1.26‐13.14, P = .019) and OS (HR, 2.78; 95%CI, 0.95‐8.17, P = .062). In a subgroup analysis, while sex and M‐stage were controlled, a much stronger negative effect of YAP1 amplification vs wild‐type in LRFS was observed (log‐rank P = .0067). Conclusion The results suggested that YAP1 amplification is a potentially useful biomarker for predicting treatment outcomes and identifying patients with a high risk of relapse who should be closely monitored.

The domesticated silkworm, Bombyx mori, is an important model system for the order Lepidoptera. Currently, based on third-generation sequencing, the chromosome-level genome of Bombyx mori has been released. However, its transcripts were mainly assembled by using short reads of second-generation sequencing and expressed sequence tags which cannot explain the transcript profile accurately. Here, we used PacBio Iso-Seq technology to investigate the transcripts from 45 developmental stages of Bombyx mori. We obtained 25,970 non-redundant high-quality consensus isoforms capturing ∼60% of previous reported RNAs, 15,431 (∼47%) novel transcripts, and identified 7,253 long non-coding RNA (lncRNA) with a large proportion of novel lncRNA (∼56%). In addition, we found that transposable elements (TEs) exonization account for 11,671 (∼45%) transcripts including 5,980 protein-coding transcripts (∼32%) and 5,691 lncRNAs (∼79%). Overall, our results expand the silkworm transcripts and have general implications to understand the interaction between TEs and their host genes. These transcripts resource will promote functional studies of genes and lncRNAs as well as TEs in the silkworm.

Gastrointestinal tract cancers have high incidence and mortality in China, but their molecular characteristics have not been fully investigated. We sequenced 432 tumor samples from the colorectum, stomach, pancreas, gallbladder, and biliary tract to investigate cancer-related mutations and detail the landscape of microsatellite instability (MSI), tumor mutation burden (TMB), and chromosomal instability (CIN). We observed the highest TMB in colorectal and gastric cancers and the lowest TMB in gastrointestinal stromal tumors (GISTs). Twenty-four hyper-mutated tumors were identified only in colorectal and gastric cancers, with a significant enrichment of mutations in the polymerase genes (POLE, POLD1, and POLH) and mismatch repair (MMR) genes. Additionally, CIN preferentially occurred in colorectal and gastric cancers, while pancreatic, gallbladder, and biliary duct cancers had a much lower CIN. High CIN was correlated with a higher prevalence of malfunctions in chromosome segregation and cell cycle genes, including the copy number loss of WRN, NAT1, NF2, and BUB1B, and the copy number gain of MYC, ERBB2, EGFR, and CDK6. In addition, TP53 mutations were more abundant in high-CIN tumors, while PIK3CA mutations were more frequent in low-CIN tumors. In colorectal and gastric cancers, tumors with MSI demonstrated much fewer copy number changes than microsatellite stable (MSS) tumors. In colorectal and gastric cancers, the molecular characteristics of tumors revealed the mutational diversity between the different anatomical origins of tumors. This study provides novel insights into the molecular landscape of Chinese gastrointestinal cancers and the genetic differences between tumor locations, which could be useful for future clinical patient stratification and targeted interventions.

Wearable self-powered electronics receive considerable attention, but development is hindered by bulky designs, energy storage systems, and safety issues. Here the authors demonstrate a flexible solar-charging integrated unit based on printed magnesium ion aqueous asymmetric supercapacitors.

Non–small cell lung cancer (NSCLC) with activating EGFR mutations in exon 19 and 21 typically responds to EGFR tyrosine kinase inhibitors (TKI); however, for some patients, responses last only a few months. The underlying mechanisms of such short responses have not been fully elucidated. Here, we sequenced the genomes of 16 short-term responders (SR) that had progression-free survival (PFS) of less than 6 months on the first-generation EGFR TKI and compared them to 12 long-term responders (LR) that had more than 24 months of PFS. All patients were diagnosed with advanced lung adenocarcinoma and harbored EGFR 19del or L858R mutations before treatment. Paired tumor samples collected before treatment and after relapse (or at the last follow-up) were subjected to targeted next-generation sequencing of 416 cancer-related genes. SR patients were significantly younger than LR patients (P < .001). Collectively, 88% of SR patients had TP53 variations compared to 13% of LR patients (P < .001). Additionally, 37.5% of SR patients carried EGFR amplifications compared to 8% of LR patients. Other potential primary resistance factors were also identified in the pretreatment samples of 12 SR patients (75%), including PTEN loss; BIM deletion polymorphism; and amplifications of EGFR, ERBB2, MET, HRAS, and AKT2. Comparatively, only three LR patients (25%) were detected with EGFR or AKT1 amplifications that could possibly exert resistance. The diverse preexisting resistance mechanisms in SR patients revealed the complexity of defining treatment strategies even for EGFR-sensitive mutations.

Insect cuticle is critical for the environmental adaptability and insecticide resistance of insects. However, there is no clear understanding of the structure and protein components of the cuticle during each developmental stage of holometabolous insects, and knowledge about the protein components within each layer is vague. We conducted serial sectioning, cuticular structure analysis, and transcriptome sequencing of the larval, pupal, and adult cuticles of Bombyx mori. The deposition processes of epicuticle, exocuticle, and endocuticle during larval, pupal, and adult cuticle formation were similar. Transcriptome analysis showed that these cuticle formations share 74% of the expressed cuticular protein (CP) genes and 20 other structural protein genes, such as larval serum protein and prisilkin. There are seven, six, and eleven stage-specific expressed CP genes in larval, pupal, and adult cuticles, respectively. The types and levels of CP genes may be the key determinants of the properties of each cuticular layer. For example, the CPs of the RR-2 protein family with high contents of histidine (His) are more essential for the exocuticle. Functional analysis suggested that BmorCPAP1-H is involved in cuticle formation. This study not only offers an in-depth understanding of cuticle morphology and protein components but also facilitates the elucidation of molecular mechanisms underlying cuticle formation in future studies.

The insect glycoside hydrolase family 20 β-N-acetylhexosaminidases (HEXs) are key enzymes involved in chitin degradation. In this study, nine HEX genes in Bombyx mori were identified by genome-wide analysis. Bioinformatic analysis based on the transcriptome database indicated that each gene had a distinct expression pattern. qRT-PCR was performed to detect the expression pattern of the chitooligosaccharidolytic β-N-acetylglucosaminidase (BmChiNAG). BmChiNAG was highly expressed in chitin-rich tissues, such as the epidermis. In the wing disc and epidermis, BmChiNAG has the highest expression level during the wandering stage. CRISPR/Cas9-mediated BmChiNAG deletion was used to study the function. In the BmChiNAG-knockout line, 39.2% of female heterozygotes had small and curly wings. The ultrastructure of a cross-section showed that the lack of BmChiNAG affected the stratification of the wing membrane and the formation of the correct wing vein structure. The molting process of the homozygotes was severely hindered during the larva to pupa transition. Epidermal sections showed that the endocuticle of the pupa was not degraded in the mutant. These results indicate that BmChiNAG is involved in chitin catabolism and plays an important role in the molting and wing development of the silkworm, which highlights the potential of BmChiNAG as a pest control target.

Silk fibroin has the merits of biocompatibility, biodegradability, ease of processing, and feasibility of modification, which present it as a promising drug delivery material. This review focuses on the structures of silk fibroin, the controlled transformation of secondary structures, and the formation mechanism of silk fibroin-based nanoparticles (SFNPs). We also discuss the intrinsic multi-responsive, surface functionalization, and transgenic modification of SFNPs for drug delivery.

Many insects spin cocoons to protect the pupae from unfavorable environments and predators. After emerging from the pupa, the moths must escape from the sealed cocoons. Previous works identified cocoonase as the active enzyme loosening the cocoon to form an escape-hatch. Here, using bioinformatics tools, we show that cocoonase is specific to Lepidoptera and that it probably existed before the occurrence of lepidopteran insects spinning cocoons. Despite differences in cocooning behavior, we further show that cocoonase evolved by purification selection in Lepidoptera and that the selection is more intense in lepidopteran insects spinning sealed cocoons. Experimentally, we applied gene editing techniques to the silkworm Bombyx mori, which spins a dense and sealed cocoon, as a model of lepidopteran insects spinning sealed cocoons. We knocked out cocoonase using the CRISPR/Cas9 system. The adults of homozygous knock-out mutants were completely formed and viable but stayed trapped and died naturally in the cocoon. This is the first experimental and phenotypic evidence that cocoonase is the determining factor for breaking the cocoon. This work led to a novel silkworm strain yielding permanently intact cocoons and provides a new strategy for controlling the pests that form cocoons.

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 (ω1 = 0.044 vs. ω2 = 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.

Human sepiapterin reductase (SR) deficiency is an inherited disease caused by SPR gene mutations and is a monoamine neurotransmitter disorder. Here, we investigated whether the silkworm lemon mutant could serve as a model of SR deficiency. A point mutation in the BmSPR gene led to a five amino acid deletion at the carboxyl terminus in the lemon mutant. In addition, classical phenotypes seen in SR deficient patients were observed in the lemon mutant, including a normal phenylalanine level, a decreased dopamine and serotonin content, and an increased neopterin level. A recovery test showed that the replenishment of l-dopa significantly increased the dopamine level in the lemon mutant. The silkworm lemon mutant also showed negative behavioural abilities. These results suggest that the silkworm lemon mutant has an appropriate genetic basis and meets the biochemical requirements to be a model of SR deficiency. Thus, the silkworm lemon mutant can serve as a candidate animal model of SR deficiency, which may be helpful in facilitating accurate diagnosis and effective treatment options of SR deficiency.

Abstract Background Esophageal sarcomatoid carcinoma (ESC) is a rare disease with a mixture of both carcinomatous and sarcomatous components in the tumor. Its genetic background and mechanisms of oncogenesis remain largely unknown. Methods Here we performed targeted next generation sequencing (NGS) on a pan-cancer gene panel in 15 ESC tumors to explore their genetic alterations, and aimed to identify clinically actionable mutations for future treatment instructions. Results TP53 alterations were identified in all patients. Alterations in receptor tyrosine kinases (RTK) were identified in 10 out of 15 patients. Members of downstream RAS and PI3-kinase pathways are also mutated in 10 patients, and PIK3CA is the top mutated gene in these pathways. In addition, we identified mutations on histone modification genes in 5 patients, including histone acetyltransferase gene EP300 and its homologue CREBBP, lysine methyltransferase genes KMT2A and KMT2B, and lysine demethylase gene KDM5A. Finally, mismatch repair (MMR) genes and proofreading gene POLE all together were mutated in one third of the ESC patients. Conclusions This is the first study to unravel the mutational profile of ESC tumors. Our findings could match 9 patients to the targeted therapies currently available in clinical practice or in active clinical trials, suggesting the potential utility of targeted therapies for this rare disease in the future.

Abstract 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.

Diapause is one of the survival strategies of insects for confronting adverse environmental conditions. Bombyx mori displays typical embryonic diapause, and offspring diapause depends on the incubation environment of the maternal embryo in the bivoltine strains of the silkworm. However, the molecular mechanisms of the diapause induction process are still poorly understood. In this study, we compared the differentially expressed miRNAs (DEmiRs) in bivoltine silkworm embryos incubated at diapause- (25 °C) and non-diapause (15 °C)-inducing temperatures during the blastokinesis (BK) and head pigmentation (HP) phases using transcriptome sequencing. There were 411 known miRNAs and 71 novel miRNAs identified during the two phases. Among those miRNAs, there were 108 and 74 DEmiRs in the BK and HP groups, respectively. By the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the predicted target genes of the DEmiRs, we found that aside from metabolism, the targets were also enriched in phototransduction-fly and insect hormone biosynthesis in the BK group and the HP group, respectively. Dual luciferase reporter assay illustrated that bmo-miR-6497-3p directly regulated Bmcycle and subsequently regulated the expression of circadian genes. These results imply that microRNAs, as vitally important regulators, respond to different temperatures and participate in the diapause induction process across species.

Abstract Cancer is a disease of complex genetic alterations, and comprehensive genetic diagnosis is beneficial to match each patient to appropriate therapy. However, acquisition of representative tumor samples is invasive and sometimes impossible. Circulating tumor DNA (ctDNA) is a promising tool to use as a non-invasive biomarker for cancer mutation profiling. Here we implemented targeted next generation sequencing (NGS) with a customized gene panel of 382 cancer-relevant genes on 605 ctDNA samples in multiple cancer types. Overall, tumor-specific mutations were identified in 87% of ctDNA samples, with mutation spectra highly concordant with their matched tumor tissues. 71% of patients had at least one clinically-actionable mutation, 76% of which have suggested drugs approved or in clinical trials. In particular, our study reveals a unique mutation spectrum in Chinese lung cancer patients which could be used to guide treatment decisions and monitor drug-resistant mutations. Taken together, our study demonstrated the feasibility of clinically-useful targeted NGS-based ctDNA mutation profiling to guide treatment decisions in cancer.

Cys2-His2 zinc finger (C2H2-ZF) proteins comprise the largest class of putative eukaryotic transcription factors. The zinc finger motif array is highly divergent, indicating that most proteins will have distinctive binding sites and perform different functions. However, the binding sites and functions of the majority of C2H2-ZF proteins remain unknown. In this study, we identified 327 C2H2-ZF protein genes in the silkworm, 290 in the monarch butterfly, 243 in the fruit fly, 107 in elegans, 673 in mouse, and 1,082 in human. The C2H2-ZF protein genes of the silkworm were classified into three main grouping clades according to a phylogenetic classification, and 312 of these genes could be mapped onto 27 chromosomes. Most silkworm C2H2-ZF protein genes exhibited specific expression in larval tissues. Furthermore, several C2H2-ZF protein genes had sex-specific expression during metamorphosis. In addition, we found that some C2H2-ZF protein genes are involved in metamorphosis and female reproduction by using expression clustering and gene annotation analysis. Among them, five genes were selected, BGIBMGA002091 (CTCF), BGIBMGA006492 (fru), BGIBMGA006230 (wor), BGIBMGA004640 (lola), and BIGBMGA004569, for quantitative real-time PCR analysis from larvae to adult ovaries. The results showed that the five genes had different expression patterns in ovaries, among which BGIBMGA002091 (CTCF) gene expression level was the highest, and its expression level increased rapidly in late pupae and adult stages. These findings provide a basis for further investigation of the functions of C2H2-ZF protein genes in the silkworm, and the results offer clues for further research into the development of metamorphosis and female reproduction in the silkworm.

Flight loss has occurred in many winged insect taxa. The flightless silkmoth Bombyx mori, is domesticated from the wild silkmoth, Bombyx mandarina, which can fly. In this paper, we studied morphological characteristics attributed to flightlessness in silkmoths. Three domestic flightless B. mori strains and one B. mandarina population were used to compare morphological components of the flight apparatus, including wing characteristics (shape, forewing area, loading, and stiffness), flight muscle (weight, ratio, and microscopic detail) and body mass. Compared with B. mandarina, B. mori strains have a larger body, greater wing loading, more flexible wings and a lower flight muscle ratio. The arrangement in microscopy of dorsal longitudinal flight muscles (DLFMs) of B. mori was irregular. Comparative analysis of the sexes suggests that degeneration of flight muscles and reduction of wing mechanical properties (stiffness) are associated with silkmoth flightlessness. The findings provide important clues for further research of the molecular mechanisms of B. mori flight loss.

Previous studies have shown that regeneration gene 3 (reg3) is significantly expressed in gastric mucosa tissues with Helicobacter pylori (HP) cytotoxin-associated gene A (CagA)-positive (HP-CagA+). CagA-positive HP increases the risk of gastric cancer. The purpose of this study was to investigate the correlation between reg3 and HP-CagA+ and explore the effects of reg3 on the proliferation of gastric cancer cells and the development of tissues and organs. We analyzed the expression of reg3 in human tissues and organs. The results showed that reg3 expression in gastric tissues was significantly higher than that in other tissues and organs. In addition, reg3 influenced the prognosis of gastric, lung, and ovarian cancers. Immunohistochemical analysis indicated that the expression of reg3 and CagA in cancerous tissues was higher than that in adjacent tissues. HP-CagA+ infection of gastric cancer cells promotes reg3 expression, suggesting that reg3 may be a target gene of CagA in gastric cancer, which together affects the formation and development of gastric cancer. reg3 and CagA promote cell proliferation, and then affect the development of mouse tissues and organs by regulating G1/S phase transition of the cell cycle via the formation of the cell cycle-dependent complex CDK4/CyclinD1. This is the first study that shows the influence of CagA on the cell cycle and induction of cell proliferation by promoting reg3 expression.

In insects, abdominal appendages, also called prolegs, vary due to adaptive evolution. Mutations on prolegs within species provide insights to better understand the mechanisms underlying appendage development and diversity. In silkworm Bombyx mori, extra-crescents and degenerated abdominal legs (Edl) mutant, belonging to the E pseudoallele group, is a spontaneous mutation that adds crescents and degenerates prolegs on the third abdominal segment (A3). This mutation may be a homeotic transformation of A3 to A2. In this study, the Edl locus was mapped within approximately a 211 Kb region that is 10 Kb upstream of Bmabdominal-A (Bmabd-A). RT-quantitative PCR (RT-qPCR) and Western blot analysis of Bmabd-A expression showed a slight but significant decrease, while the expression of BmUltrabithorax (BmUbx) was up-regulated in the Edl mutant compared to wildtype (Dazao). Moreover, we also found that BmDistal-less (BmDll), which regulated the development of distal proleg structures, was missing at the tips of the A3 prolegs in the Edl mutant compared to BmDll expression in normally developed prolegs in both the wildtype and mutant. Collectively, we identified approximately a 211 Kb region in the Edl locus that regulates BmUbx and Bmabd-A expression and found that changes in BmUbx and Bmabd-A expression may lead to the loss of distal proleg structures in B. mori.

Connexin40 (Cx40)-containing gap-junction channels are expressed in the atrial myocardium and provide a low-resistance passage for rapid impulse propagation. A germline mutation in the GJA5 gene, which encodes Cx40, resulting in a truncated Cx40 (Q49X) was identified in a large Chinese family with lone (idiopathic) atrial fibrillation (AF). This mutation co-segregated with seven AF probands in an autosomal-dominant way over generations. To test the hypothesis that this Cx40 mutant affects the distribution and function of atrial gap junctions, we studied the Q49X mutant in gap-junction-deficient HeLa and N2A cells. The Q49X mutant, unlike wild-type Cx40, was typically localized in the cytoplasm and failed to form gap-junction plaques at cell-cell interfaces. When the Q49X mutant was co-expressed with Cx40 or Cx43, the mutant substantially reduced the gap-junction plaque formation of Cx40 and Cx43. Electrophysiological studies revealed no electrical coupling of cell pairs expressing the mutant alone and a significant decrease in the coupling conductance when the mutant was co-expressed with Cx40 or Cx43. Further colocalization experiments with the organelle residential proteins indicate that Q49X was retained in the endoplasmic reticulum. These findings provide evidence that the Q49X mutant is capable of impairing gap-junction distribution and function of key atrial connexins, which might play a role in the predisposition to and onset of AF.

Wnt is a family of conserved glycoproteins that participate in a variety of important biological processes including embryo development, cell proliferation and differentiation, and tissue regeneration. The Wnt family is a metazoan novelty found in all animal phyla. Studies have revealed that the number of Wnt genes varies among species, presumably due to reproduction and loss of genes during evolution. However, a comprehensive inventory of Wnt genes in Lepidoptera is lacking. In this study, we identified the repertoire of Wnt genes in the silkworm and seven other species of Lepidoptera and obtained eight Wnt genes (Wnt1, Wnt5–Wnt7, Wnt9–Wnt11, and WntA) in each species. Four of these Wnt genes are clustered in two orientations (5′-Wnt9-Wnt1-Wnt6-Wnt10-3′ and 5′-Wnt10-Wnt6-Wnt1-Wnt9-3′) in both moths and butterflies. Transcript analysis of Wnt in silkworm embryonic stages showed that each BmWnt gene had a unique expression pattern during embryological development. Analysis of a larval stage revealed differential expression of Wnt family members in diverse tissues. Our study provides an overview of the Wnt family in Lepidoptera and will inspire further functional study of the Wnt genes in the silkworm.

In this study, eight types of materials including nanoparticles (Cu and CaCO3), metallic ions (Ca2+ and Cu2+), and amino acid substances (serine, tyrosine, sericin amino acid, and fibroin amino acid) were used as additives in silkworm diets to obtain in-situ modified silk fiber composites. The results indicate that tyrosine and fibroin amino acids significantly increase potassium content in silk fibers and induce the transformation of α-helices and random coils to β-sheet structures, resulting in higher crystallinities and better mechanical properties. However, the other additives-modified silk fibers show a decrease in β-sheet contents and a slight increase or even decrease in tensile strengths. This finding provides a green and effective approach to produce mechanically enhanced silk fibers with high crystallinity on a large scale. Moreover, the modification mechanisms of these additives were discussed in this study, which could offer new insights into the design and regulation of modified fibers or composites with desirable properties and functions.

As a source of insect polymorphism, melanism plays an important role in ecological adaption and usually endows advantageous phenotypic-effects on insects. However, due to the mechanistic diversity, there are knowledge gaps in the molecular mechanisms underlying melanism and the related phenotypes. In silk moths, a recessive melanic mutant (sex-controlled melanism, sml) strain exhibits extended adult longevity. We took a transcriptome approach to perform a comparative analysis between this sml strain and a wild-type strain (Dazao). Our analysis resulted in the identification of 59 unique differentially expressed genes in the melanic mutant. Two key genes (laccase2 and yellow) involved in melanin formation were significantly up-regulated in melanic individuals. The laccase2 B-type isoform (BGIBMGA006746) was found to likely participate in the silkworm cuticular melanism process at late pupal stage. Moreover, we discovered 22 cuticular protein encoding genes with the possible function in melanin transport and/or maintenance. Based on our findings, we presume that the longer survival of the melanic sml male moths might be associated with the enhanced antioxidant defense systems and a reduction in the insulin/IGF-1 signaling pathway (IIS). These findings will facilitate the understanding of the molecular basis underlying melanism and the derived phenotypic-effects.

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.

In insects, the integument provides mechanical support for the whole body and protects them from infections, physical and chemical injuries, and dehydration. Diversity in integument properties is often related to body shape, behavior, and survival rate. The stick (sk) silkworm is a spontaneous mutant with a stick-like larval body that is firm to the touch and, thus, less flexible. Analysis of the mechanical properties of the cuticles at day 3 of the fifth instar (L5D3) of sk larvae revealed higher storage modulus and lower loss tangent. Transcriptome sequencing identified a total of 19,969 transcripts that were expressed between wild-type Dazao and the sk mutant at L5D2, of which 11,596 transcripts were novel and detected in the integument. Differential expression analyses identified 710 upregulated genes and 1009 downregulated genes in the sk mutant. Gene Ontology (GO) enrichment analysis indicated that four chitin-binding peritrophin A domain genes and a chitinase gene were upregulated, whereas another four chitin-binding peritrophin A domain genes, a trehalase, and nine antimicrobial peptides were downregulated. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that two functional pathways, namely, fructose and mannose metabolism and tyrosine metabolism, were significantly enriched with differentially-expressed transcripts. This study provides a foundation for understanding the molecular mechanisms underlying the development of the stiff exoskeleton in the sk mutant.

In this study, using round paper tubes (PTs) and rectangular cardboard boxes (CBs) as external constraints to control the size of the cocooning space, we fabricated a series of modified silkworm cocoons (PT cocoons and CB cocoons). Their microstructures, morphologies, compositions, and mechanical properties were characterized and compared with normal silkworm cocoons. These two kinds of modified silkworm cocoons exhibit dense and homogeneous layer structures. Tensile test results indicate that above a size limit of cocooning space, their tensile strengths, Young’s moduli, and strain energy densities increase with the decrease in cocooning space. Especially in comparison with the normal cocoons, the tensile strength and Young’s modulus of the PT-14 cocoon increase by 44% and 100%, respectively. Meanwhile, PT cocoons and CB cocoons, except PT-12, also possess better peeling resistance than normal cocoons. Owing to the dense structure and low porosity, the modified cocoons form robust fiber networks that result in high strength and toughness. This study provides a green and efficient method to fabricate mechanically enhanced silkworm cocoons with special shapes and dense layer structures. The method can be easily subjected to further modification processes and has potential applications in the production of high-performance green cocoon composites and biomimetic materials.

The tyrosine kinases (TKs) are important parts of metazoan signaling pathways and play significant roles in cell growth, development, apoptosis and disease. Genome-wide characterization of TKs has been conducted in many metazoans, however, systematic information about this family in Lepidoptera is still lacking. We retrieved 33 TK-encoding genes in silkworm and classified them into 25 subfamilies by sequence analysis, without members in AXL, FRK, PDGFR, STYK1 and TIE subfamilies. Although domain sequences in each subfamily are conserved, TKs in vertebrates tend to be remarkably conserved and stable. Our results of phylogenetic analysis supported the previous conclusion for the second major expansion of TK family. Gene-Ontology (GO) analysis revealed that a higher proportion of BmTKs played roles in binding, catalysis, signal transduction, metabolism, biological regulation and response to stimulus, compared to all silkworm genes annotated in GO. Moreover, the expression profile analysis of BmTKs among multiple tissues and developmental stages demonstrated that many genes exhibited stage-specific and/or sex-related expression during embryogenesis, molting and metamorphosis, and that 8 BmTKs presented tissue-specific high expression. Our study provides systematic description of silkworm tyrosine kinases, and may also provide further insights into metazoan TKs and assist future studies addressing their functions.

Phenylketonuria (PKU) is an inborn error of metabolism caused by mutations in the phenylalanine hydroxylase (PAH) gene or by defects in the tetrahydrobiopterin (BH4) synthesis pathway. Here, by positional cloning, we report that the 6-pyruvoyl-tetrahydropterin synthase (PTPS) gene, encoding a key enzyme of BH4 biosynthesis, is responsible for the alc (albino C) mutation that displays pale body color, head shaking, and eventually lethality after the first molting in silkworm. Compared to wild type, the alc mutant produced more substrates (phenylalanine (Phe) and tyrosine (Tyr)) and generated less DOPA and dopamine. Application of 2,4-diamino-6-hydroxypyrimidine (DAHP) to block BH4 synthesis in the wild type effectively produced the alc-like phenotype, while BH4 supplementation rescued the defective body color and lethal phenotype in both alc and DAHP-treated individuals. The detection of gene expressions and metabolic substances after drugs treatments in alc and normal individuals imply that silkworms and humans have a high similarity in the drugs metabolic features and the gene pathway related to BH4 and the dopamine biosynthesis. We propose that the alc mutant could be used as an animal model for drug evaluation for BH4-deficient PKU.

WD40 proteins are scaffolding molecules in protein-protein interactions and play crucial roles in fundamental biological processes. Genome-wide characterization of WD40 proteins in animals has been conducted solely in humans. We retrieved 172 WD40 protein genes in silkworm (BmWD40s) and identified these genes in 7 other insects, 9 vertebrates and 5 nematodes. Comparative analysis revealed that the WD40 protein gene family underwent lineage-specific expansions during animal evolution, but did not undergo significant expansion during insect evolution. The BmWD40s were categorized into five clusters and 12 classes according to the phylogenetic classification and their domain architectures, respectively. Sequence analyses indicated that tandem and segmental duplication played minor roles in producing the current number of BmWD40s, and domain recombination events of multi-domain BmWD40s might have occurred mainly after gene duplication events. Gene Ontology (GO) analysis revealed that a higher proportion of BmWD40s was involved in processes, such as binding, transcription-regulation and cellular component biogenesis, compared to all silkworm genes annotated in GO. Microarray-based analysis demonstrated that many BmWD40s had tissue-specific expression and exhibited high and/or sex-related expression during metamorphosis. These findings contribute to a better understanding of the evolution of the animal WD40 protein family and assist the study of the functions of BmWD40s.

Bodies are often made of repeated units, or serial homologs, that develop using the same core gene regulatory network. Local inputs and modifications to this network allow serial homologs to evolve different morphologies, but currently we do not understand which modifications allow these repeated traits to evolve different levels of phenotypic plasticity. Here we describe variation in phenotypic plasticity across serial homologous eyespots of the butterfly Bicyclus anynana, hypothesized to be under selection for similar or different functions in the wet and dry seasonal forms. Specifically, we document the presence of eyespot size and scale brightness plasticity in hindwing eyespots hypothesized to vary in function across seasons, and reduced size plasticity and absence of brightness plasticity in forewing eyespots hypothesized to have the same function across seasons. By exploring the molecular and physiological causes of this variation in plasticity across fore and hindwing serial homologs we discover that: 1) temperature experienced during the wandering stages of larval development alters titers of an ecdysteroid hormone, 20-hydroxyecdysone (20E), in the hemolymph of wet and dry seasonal forms at that stage; 2) the 20E receptor (EcR) is differentially expressed in the forewing and hindwing eyespot centers of both seasonal forms during this critical developmental stage; and 3) manipulations of EcR signaling disproportionately affected hindwing eyespots relative to forewing eyespots. We propose that differential EcR expression across forewing and hindwing eyespots at a critical stage of development explains the variation in levels of phenotypic plasticity across these serial homologues. This finding provides a novel signaling pathway, 20E, and a novel molecular candidate, EcR, for the regulation of levels of phenotypic plasticity across body parts or serial homologs.

Gap junction (GJ) channels provide direct passage for ions and small molecules to be exchanged between neighbouring cells and are crucial for many physiological processes. GJ channels can be gated by transjunctional voltage (known as Vj-gating) and display a wide range of unitary channel conductance (γj), yet the domains responsible for Vj-gating and γj are not fully clear. The first extracellular domain (E1) of several connexins has been shown to line part of their GJ channel pore and play important roles in Vj-gating properties and/or ion permeation selectivity. To test roles of the E1 of Cx50 GJ channels, we generated a chimera, Cx50Cx36E1, where the E1 domain of Cx50 was replaced with that of Cx36, a connexin showing quite distinct Vj-gating and γj from those of Cx50. Detailed characterizations of the chimera and three point mutants in E1 revealed that, although the E1 domain is important in determining γj, the E1 domain of Cx36 is able to effectively function within the context of the Cx50 channel with minor changes in Vj-gating properties, indicating that sequence differences between the E1 domains in Cx36 and Cx50 cannot account for their drastic differences in Vj-gating and γj. Our homology models of the chimera and the E1 mutants revealed that electrostatic properties of the pore-lining residues and their contribution to the electric field in the pore are important factors for the rate of ion permeation of Cx50 and possibly other GJ channels.

Mutations in GJA5 encoding the gap junction protein connexin40 (Cx40) have been linked to lone atrial fibrillation. Some of these mutants result in impaired gap junction function due to either abnormal connexin localization or impaired gap junction channels, which may play a role in promoting atrial fibrillation. However, the effects of the atrial fibrillation-linked Cx40 mutants on hemichannel function have not been studied. Here we investigated two atrial fibrillation-linked germline Cx40 mutants, V85I and L221I. These two mutants formed putative gap junction plaques at cell-cell interfaces, with similar gap junction coupling conductance as that of wild-type Cx40. Connexin deficient HeLa cells expressing either one of these two mutants displayed prominent propidium iodide-uptake distinct from cells expressing wild-type Cx40 or other atrial fibrillation-linked Cx40 mutants, I75F, L229M, and Q49X. Propidium iodide-uptake was sensitive to [Ca2+]o and the hemichannel blockers, carbenoxolone, flufenamic acid and mefloquine, but was not affected by the pannexin 1 channel blocking agent, probenecid, indicating that uptake is most likely mediated via connexin hemichannels. A gain-of-hemichannel function in these two atrial fibrillation-linked Cx40 mutants may provide a novel mechanism underlying the etiology of atrial fibrillation.

Butterfly eyespots may have evolved from the recruitment of pre-existent gene circuits or regulatory networks into novel locations on the wing. Gene expression data suggests one such circuit, the Hedgehog (Hh) signaling pathway and its target gene engrailed (en), was recruited from a role in patterning the anterior-posterior insect wing axis to a role patterning butterfly eyespots. However, while Junonia coenia expresses hh and en both in the posterior compartment of the wing and in eyespot centers, Bicyclus anynana lacks hh eyespot-specific expression. This suggests that Hh signaling may not be functioning in eyespot development in either species or that it functions in J. coenia but not in B. anynana. In order to test these hypotheses, we performed functional tests of Hh signaling in these species. We investigated the effects of Hh protein sequestration during the larval stage on en expression levels, and on wing size and eyespot size in adults. Hh sequestration led to significantly reduced en expression and to significantly smaller wings and eyespots in both species. But while eyespot size in B. anynana was reduced proportionately to wing size, in J. coenia, eyespots were reduced disproportionately, indicating an independent role of Hh signaling in eyespot development in J. coenia. We conclude that while Hh signaling retains a conserved role in promoting wing growth across nymphalid butterflies, it plays an additional role in eyespot development in some, but not all, lineages of nymphalid butterflies. We discuss our findings in the context of alternative evolutionary scenarios that led to the differential expression of hh and other Hh pathway signaling members across nymphalid species.

Catecholamine metabolism plays an important role in the determination of insect body color and cuticle sclerotization. To date, limited research has focused on these processes in silkworm. In the current study, we analyzed the interactions between catecholamines and melanin genes and their effects on the pigmentation patterns and physical properties of sclerotized regions in silkworm, using the melanic mutant melanism (mln) silkworm strain as a model. Injection of β-alanine into mln mutant silkworm induced a change in catecholamine metabolism and turned its body color yellow. Further investigation of the catecholamine content and expression levels of the corresponding melanin genes from different developmental stages of Dazao-mln (mutant) and Dazao (wild-type) silkworm revealed that at the larval and adult stages, the expression patterns of melanin genes precipitated dopamine accumulation corresponding to functional loss of Bm-iAANAT, a repressive effect of excess NBAD on ebony, and upregulation of tan in the Dazao-mln strain. During the early pupal stage, dopamine did not accumulate in Dazao-mln, since upregulation of ebony and black genes led to conversion of high amounts of dopamine into NBAD, resulting in deep yellow cuticles. Scanning electron microscope analysis of a cross-section of adult dorsal plates from both wild-type and mutant silkworm disclosed the formation of different layers in Dazao-mln owing to lack of NADA, compared to even and dense layers in Dazao. Analysis of the mechanical properties of the anterior wings revealed higher storage modulus and lower loss tangent in Dazao-mln, which was closely associated with the altered catecholamine metabolism in the mutant strain. Based on these findings, we conclude that catecholamine metabolism is crucial for the color pattern and physical properties of cuticles in silkworm. Our results should provide a significant contribution to Lepidoptera cuticle tanning research.