Environmental concentrations of microplastic-induced gut microbiota and metabolite disruption in silkworm, Bombyx mori

Microplastics (MPs) existing extensively in various ecosystems can be ingested by marine organisms and enter the food chain, resulting the health risks from the presence of MPs in aquatic and terrestrial ecosystems. In the present study, an ideal model for Lepidoptera, the silkworm, Bombyx mori, was exposed to environmental concentrations (0.125 μg, 0.25 μg or 0.5 μg/diet) of MPs for 5 days, and the global changes in gut microbes and metabolites were subsequently examined via 16S rDNA sequencing and GC‒MS-based metabolomics. The results showed that MPs exposure did not seriously threaten survival but may regulate signaling pathways involved in development and cocoon production. MPs exposure induced gut microbiota perturbation according to the indices of α-diversity and β-diversity, and the functional prediction of the altered microbiome and associated metabolites demonstrated the potential roles of the altered microbiome following MPs exposure in the metabolic and physiological states of silkworm. The metabolites markedly altered following MPs exposure may play vital biological roles in energy metabolism, lipid metabolism, xenobiotic detoxification and the immune system by directly or indirectly affecting the physiological state of silkworms. These findings contribute to assessing the health risks of MPs exposure in model insects and provide novel insight into the toxicity mechanism of MPs.

Rapid production of chimeric silkworm/spider silk with improved mechanical properties by infection of nonpermissive Bombyx mori with recombinant AcMNPV harboring native-size of spidroin genes

Spider silks with excellent mechanical properties attract more attention from scientists worldwide, and the dragline silk that serves as the framework of the spider's web is considered one of the strongest fibers. However, it is unfeasible for large-scale production of spider silk due to its highly territorial, cannibalistic, predatory, and solitary behavior. Herein, to alleviate some of these problems and explore aneasy way to produce spider fibers, we constructed recombinant baculovirus Autographa californica multiple nucleopolyhedrovirus (AcMNPV) simultaneously expressing Trichonephila clavipes native ampullate spidroin 2 (MaSp-G) and spidroin 1 (MaSp-C) driven by the promoters of silkworm fibroin genes, to infect the nonpermissive Bombyx mori larvae at the fifth instar. MaSp-G and MaSp-C were co-expressed in the posterior silk glands (PSGs) of infected silkworms and successfully secreted into the lumen of the silk gland for fibroin globule assembly. The integration of MaSp-G and MaSp-C into silkworm silk fibers significantly improved the mechanical properties of these chimeric silk fibers, especially the strength and extensibility, which may be caused by the increment of β-sheet in the chimeric silkworm/spider silk fiber. These results demonstrated that silkworms could be developed as the nonpermissive heterologous host for the mass production of chimeric silkworm/spider silk fibers via the recombinant baculovirus AcMNPV.

Identification and characterization of extrachromosomal circular DNA in the silk gland of Bombyx mori

The silk gland cells of silkworm are special cells which only replicate DNA in the nucleus without cell division throughout the larval stage. The extrachromosomal circular DNAs (eccDNAs) have not yet been reported in the silk gland of silkworms. Herein, we have explored the characterization of eccDNAs in the posterior silk gland of silkworms. A total of 35 346 eccDNAs were identified with sizes ranging from 30 to 13 569 549 bp. Motif analysis revealed that dual direct repeats are flanking the 5' and 3' breaking points of eccDNA. The sequences exceeding 1 kb length in eccDNAs present palindromic sequence characteristics flanking the 5' and 3' breaking points of the eccDNA. These motifs might support possible models for eccDNA generation. Genomic annotation of the eccDNA population revealed that most eccDNAs (58.6%) were derived from intergenic regions, whereas full or partial genes were carried by 41.4% of eccDNAs. It was found that silk protein genes fib-H, fib-L, and P25, as well as the transcription factors SGF and sage, which play an important regulatory role in silk protein synthesis, could be carried by eccDNAs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses showed that the genes carried by eccDNAs were mainly associated with the development and metabolism-related signaling pathways. Moreover, it was found that eccDNAfib-L could promote the transcription of fib-L gene. Overall, the results of the present study not only provide a novel perspective on the mechanism of silk gland development and silk protein synthesis but also complement previously reported genome-scale eccDNA data supporting that eccDNAs are common in eukaryotes.

A novel chimeric RNA originating from BmCPV S4 and Bombyx mori HDAC11 transcripts regulates virus proliferation

Polymerases encoded by segmented negative-strand RNA viruses cleave 5'-m7G-capped host transcripts to prime viral mRNA synthesis ("cap-snatching") to generate chimeric RNA, and trans-splicing occurs between viral and cellular transcripts. Bombyx mori cytoplasmic polyhedrosis virus (BmCPV), an RNA virus belonging to Reoviridae, is a major pathogen of silkworm (B. mori). The genome of BmCPV consists of 10 segmented double-stranded RNAs (S1-S10) from which viral RNAs encoding a protein are transcribed. In this study, chimeric silkworm-BmCPV RNAs, in which the sequence derived from the silkworm transcript could fuse with both the 5' end and the 3' end of viral RNA, were identified in the midgut of BmCPV-infected silkworms by RNA_seq and further confirmed by RT-PCR and Sanger sequencing. A novel chimeric RNA, HDAC11-S4 RNA 4, derived from silkworm histone deacetylase 11 (HDAC11) and the BmCPV S4 transcript encoding viral structural protein 4 (VP4), was selected for validation by in situ hybridization and Northern blotting. Interestingly, our results indicated that HDAC11-S4 RNA 4 was generated in a BmCPV RNA-dependent RNA polymerase (RdRp)-independent manner and could be translated into a truncated BmCPV VP4 with a silkworm HDAC11-derived N-terminal extension. Moreover, it was confirmed that HDAC11-S4 RNA 4 inhibited BmCPV proliferation, decreased the level of H3K9me3 and increased the level of H3K9ac. These results indicated that during infection with BmCPV, a novel mechanism, different from that described in previous reports, allows the genesis of chimeric silkworm-BmCPV RNAs with biological functions.

Interaction between Bombyx mori Cytoplasmic Polyhedrosis Virus NSP8 and BmAgo2 Inhibits RNA Interference and Enhances Virus Proliferation

Some insect viruses encode suppressors of RNA interference (RNAi) to counteract the antiviral RNAi pathway. However, it is unknown whether Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) encodes an RNAi suppressor. In this study, the presence of viral small interfering RNA (vsiRNA) in BmN cells infected with BmCPV was confirmed by small RNA sequencing. The Dual-Luciferase reporter test demonstrated that BmCPV infection may prevent firefly luciferase (Luc) gene silencing caused by particular short RNA. It was also established that the inhibition relied on the nonstructural protein NSP8, which suggests that NSP8 was a possible RNAi suppressor. In cultured BmN cells, the expressions of viral structural protein 1 (vp1) and NSP9 were triggered by overexpression of nsp8, suggesting that BmCPV proliferation was enhanced by NSP8. A pulldown assay was conducted with BmCPV genomic double-stranded RNA (dsRNA) labeled with biotin. The mass spectral detection of NSP8 in the pulldown complex suggests that NSP8 is capable of direct binding to BmCPV genomic dsRNA. The colocalization of NSP8 and B. mori Argonaute 2 (BmAgo2) was detected by an immunofluorescence assay, leading to the hypothesis that NSP8 interacts with BmAgo2. Coimmunoprecipitation further supported the present investigation. Moreover, vasa intronic protein, a component of RNA-induced silencing complex (RISC), could be detected in the coprecipitation complex of NSP8 by mass spectrum analysis. NSP8 and the mRNA decapping protein (Dcp2) were also discovered to colocalize to processing bodies (P bodies) for RNAi-mediated gene silencing in Saccharomyces cerevisiae. These findings revealed that by interacting with BmAgo2 and suppressing RNAi, NSP8 promoted BmCPV growth. IMPORTANCE It has been reported that the RNAi pathway is inhibited by binding RNAi suppressors encoded by some insect-specific viruses belonging to Dicistroviridae, Nodaviridae, or Birnaviridae to dsRNAs to protect dsRNAs from being cut by Dicer-2. However, it is unknown whether BmCPV, belonging to Spinareoviridae, encodes an RNAi suppressor. In this study, we found that nonstructural protein NSP8 encoded by BmCPV inhibits small interfering RNA (siRNA)-induced RNAi and that NSP8, as an RNAi suppressor, can bind to viral dsRNAs and interact with BmAgo2. Moreover, vasa intronic protein, a component of RISC, was found to interact with NSP8. Heterologously expressed NSP8 and Dcp2 were colocalized to P bodies in yeast. These results indicated that NSP8 promoted BmCPV proliferation by binding itself to BmCPV genomic dsRNAs and interacting with BmAgo2 through suppression of siRNA-induced RNAi. Our findings deepen our understanding of the game between BmCPV and silkworm in regulating viral infection.

[Comparison of aspirin treatment strategies for primary prevention of cardiovascular diseases: A decision-analytic Markov modelling study]

OBJECTIVE

To compare the expected population impact of benefit and risk of aspirin treatment strategies for the primary prevention of cardiovascular diseases recommended by different guidelines in the Chinese Electronic Health Records Research in Yinzhou (CHERRY) study.

METHODS

A decision-analytic Markov model was used to simulate and compare different strategies of aspirin treatment, including: Strategy ①: Aspirin treatment for Chinese adults aged 40-69 years with a high 10-year cardiovascular risk, recommended by the 2020 Chinese Guideline on the Primary Prevention of Cardiovascular Diseases; Strategy ②: Aspirin treatment for Chinese adults aged 40-59 years with a high 10-year cardiovascular risk, recommended by the 2022 United States Preventive Services Task Force Recommendation Statement on Aspirin Use to Prevent Cardiovascular Disease; Strategy ③: Aspirin treatment for Chinese adults aged 40-69 years with a high 10-year cardiovascular risk and blood pressure well-controlled (< 150/90 mmHg), recommended by the 2019 Guideline on the Assessment and Management of Cardio-vascular Risk in China. The high 10-year cardiovascular risk was defined as the 10-year predicted risk over 10% based on the 2019 World Health Organization non-laboratory model. The Markov model simulated different strategies for ten years (cycles) with parameters mainly from the CHERRY study or published literature. Quality-adjusted life year (QALY) and the number needed to treat (NNT) for each ischemic event (including myocardial infarction and ischemic stroke) were calculated to assess the effectiveness of the different strategies. The number needed to harm (NNH) for each bleeding event (including hemorrhagic stroke and gastrointestinal bleeding) was calculated to assess the safety. The NNT for each net benefit (i.e., the difference of the number of ischemic events could be prevented and the number of bleeding events would be added) was also calculated. One-way sensitivity analysis on the uncertainty of the incidence rate of cardiovascular diseases and probabilistic sensitivity analysis on the uncertainty of hazard ratios of interventions were conducted.

RESULTS

A total of 212 153 Chinese adults, were included in this study. The number of people who were recommended for aspirin treatment Strategies ①-③ was 34 235, 2 813, and 25 111, respectively. The Strategy ③ could gain the most QALY of 403 [95% uncertainty interval (UI): 222-511] years. Compared with Strategy ①, Strategy ③ had similar efficiency but better safety, with the extra NNT of 4 (95%UI: 3-4) and NNH of 39 (95%UI: 19-132). The NNT per net benefit was 131 (95%UI: 102-239) for Strategy ①, 256 (95%UI: 181-737) for Strategy ②, and 132 (95%UI: 104-232) for Strategy ③, making Strategy ③ the most favorable option with a better QALY and safety, along with similar efficiency in terms of net benefit. The results were consistent in the sensitivity analyses.

CONCLUSION

The aspirin treatment strategies recommended by the updated guidelines on the primary prevention of cardiovascular diseases showed a net benefit for high-risk Chinese adults from developed areas. However, to balance effectiveness and safety, aspirin is suggested to be used for primary prevention of cardiovascular diseases with consideration for blood pressure control, resulting in better intervention efficiency.

BmCPV replication is suppressed by the activation of the NF-κB/autophagy pathway through the interaction of vsp21 translated by vcircRNA_000048 with ubiquitin carboxyl-terminal hydrolase

Bombyx mori cypovirus (BmCPV), a typical double-stranded RNA virus, was demonstrated to generate a viral circRNA, vcircRNA_000048, which encodes a vsp21 with 21 amino acid residues to suppress viral replication. However, the regulatory mechanism of vsp21 on virus infection remained unclear. This study discovered that vsp21 induces reactive oxygen species (ROS) generation, activates autophagy, and attenuates virus replication by inducing autophagy. Then we confirmed that the effect of vsp21-induced autophagy on viral replication was attributed to the activation of the NF-κB signaling pathway. Furthermore, we clarified that vsp21 interacted with ubiquitin carboxyl-terminal hydrolase (UCH) and that ubiquitination and degradation of phospho-IκB-α were enhanced by vsp21 via competitive binding to UCH. Finally, we validated that vsp21 activates the NF-κB/autophagy pathway to suppress viral replication by interacting with UCH. These findings provided new insights into regulating viral multiplication and reovirus-host interaction.

[Application of discrete event simulation model in analysis on cost-effectiveness of epidemiology screening]

Discrete event simulation (DES) model is based on individual data, by which discrete events over time are simulated to reflect disease progression. The effects of individual characteristics on disease progression could be considered in the DES model. Moreover, unlike state-transition models, DES model without setting of fixed cycle can contribute to more accurate estimation of event time, especially in the evaluation of the long-term effectiveness of screening strategies for complex diseases in which time dimension needs to be considered. This article introduces the general principles, construction steps, analytic methods and other relevant issues of the DES model. Based on a research case of estimating the cost-effectiveness of screening for abdominal aortic aneurysms in women aged 65 years and above in the United Kingdom, key points in applications of the DES model in analysis on effectiveness of complex disease screening are discussed in detail, including model construction and analysis and interpretation of the results. DES model can predict occurring time of discrete events accurately by establishing the distribution function of their occurring time and is increasingly used to evaluate the screening strategies for complex diseases in which time dimension needs to be considered. In the construction of DES model, it is necessary to pay close attention to the clear presentation of model structure and simulation process and follow the relevant reporting specification to conduct cost-effectiveness analysis to ensure the transparency and repeatability of the research.

BmNPV circular RNA-encoded peptide VSP39 promotes viral replication

CircRNAs are covalently closed single-stranded circular RNA molecules, which are not easily degraded by endonucleases and play vital roles in many biological processes. Currently, most studies on circRNAs focus on endogenous circRNAs in cells, and there are few studies on virus-encoded circRNAs. In this study, a viral circRNA (circRNA-000010) derived from the region (-/bp: 114514-115,319) of the complementary strand of Bombyx mori Nucleopolyhedrovirus (BmNPV) genome was identified with the circRNA-sequencing. The authenticity of viral circRNA-000010 was further confirmed by reverse transcription PCR, reverse transcription-rolling circle amplification (TCA), in situ hybridization, immunofluorescent staining, and Northern blotting. The results of overexpression and knockdown experiments showed that circRNA-000010 promoted viral replication. Furthermore, a viral small peptide VSP39 with 39 amino acid residues translated by circRNA-000010 but not its linear molecule was confirmed. Finally, VSP39 was found to promote viral replication. Our findings indicated that a viral circRNA encoded by BmNPV promoted viral replication. These findings will provide new clues for further understanding coding information of the BmNPV genome and open a new insight for investigating host-virus interactions.

A recombinant baculovirus vector vaccine (BacMCP) against the infectious spleen and kidney necrosis virus (ISKNV)

The infectious spleen and kidney necrosis virus (ISKNV) is a highly lethal virus, which has brought significant losses to aquaculture. Therefore, a new vaccine against ISKNV with high efficiency, safety and convenience must be developed. While baculoviruses are more commonly used as protein expression systems for vaccine antigen production, this paper used baculovirus technology to develop a live-vector vaccine, BacMCP, which contains the coding sequence of the major capsid protein (MCP) (GenBank accession no. AF371960) of ISKNV and is driven by a CMV promoter. Real-time PCR and immunofluorescence showed that the MCP gene was successfully delivered to and expressed in fish cells and tissues inoculated with BacMCP. Immune-related gene (IgM, TGF-β, IL-1, IL-8, TNF-α) expression was induced in BacMCP-treated groups of largemouth bass compared with control groups. Specific antibodies could be detected in the serum of BacMCP injection-vaccinated largemouth bass by ELISA. After injection or immersion vaccination with BacMCP for 21 days, largemouth bass were infected with ISKNV. The immune effect of the injected immunization on fish in different sizes was evaluated. The vaccine efficacy of injection-vaccinated bass was 100% in small bass and 85.7% in large bass. The vaccine efficacy of immersion-vaccinated small bass was 77.3%. This study suggested that BacMCP can be used as a vector-based vaccine candidate to prevent the diseases caused by ISKNV infection.

Prediction using T2-weighted magnetic resonance imaging-based radiomics of residual uterine myoma regrowth after high-intensity focused ultrasound ablation

OBJECTIVES

To develop and evaluate magnetic resonance imaging (MRI)-based radiomics models for predicting residual myoma regrowth within 1 year after high-intensity focused ultrasound (HIFU) ablation of uterine myomas.

METHODS

A retrospective analysis of residual myoma regrowth within 1 year was performed on 428 myomas in 339 patients who were diagnosed with uterine myoma and treated with HIFU ablation in two hospital centers. In total, 851 radiomics features were extracted from T2-weighted images (T2WI) obtained 1 day after HIFU ablation, and the least absolute shrinkage and selection operator in the training cohort (n = 243) was employed to select radiomics features. Support vector machines were adopted to develop radiomics, clinicoradiological and combined radiomics-clinical models to predict residual myoma regrowth, defined as an increase in residual myoma volume of > 10% between that at day 1 post HIFU and that at follow-up MRI within 1 year. These models were validated in both internal (n = 81) and external (n = 104) test cohorts. The predictive performance and clinical application of these models were assessed using receiver-operating-characteristics-curve analysis, the area under the curve (AUC) and decision-curve analysis.

RESULTS

The AUCs of the T2WI-based radiomics prediction model in the internal and external test cohorts were 0.834 (95% CI, 0.747-0.920) and 0.801 (95% CI, 0.712-0.889), respectively, and those of the clinicoradiological model were 0.888 (95% CI, 0.816-0.960) and 0.912 (95% CI, 0.851-0.973), respectively. The combined model had better predictive performance than either the radiomics or the clinicoradiological model, with AUC values of 0.922 (95% CI, 0.857-0.987) and 0.930 (95% CI, 0.880-0.980) in the internal and external test cohorts, respectively. Decision-curve analysis also indicated that application of the combined model has clinical value, this model achieving more net benefits than the other two models.

CONCLUSION

T2WI-based radiomics features can predict effectively the occurrence of residual myoma regrowth within 1 year after HIFU ablation of uterine myomas, which serves as an accurate and convenient reference for clinical decision-making. © 2022 The Authors. Ultrasound in Obstetrics & Gynecology published by John Wiley & Sons Ltd on behalf of International Society of Ultrasound in Obstetrics and Gynecology.

Grass Carp Reovirus triggers autophagy enhancing virus replication via the Akt/mTOR pathway

Autophagy impacts the replication cycle of many viruses. Grass Carp Reovirus (GCRV) is an agent that seriously affects the development of the grass carp aquaculture industry. The role of autophagy in GCRV infection is not clearly understood. In this study, we identified that GCRV infection triggered autophagy in CIK cells, which was demonstrated by transmission electron microscopy, the conversion of LC3B I to LC3B II and the level of autophagy substrate p62. Furthermore, we found that GCRV infection activated Akt-mTOR signaling pathway, and the conversion of LC3B I to LC3B II was increased by inhibiting mTOR with rapamycin (Rap) but decreased by activating Akt with insulin. We then assessed the effects of autophagy on GCRV replication. We found that inducing autophagy with Rap promoted GCRV proliferation but inhibiting autophagy with 3 MA or CQ inhibited GCRV replication in CIK cells. Moreover, it was found that enhancing Akt-mTOR activity by insulin, GCRV VP7 protein and viral titers of GCRV were decreased. Collectively, these results indicated that GCRV infection induced autophagy involved in GCRV replication via the Akt-mTOR signal pathway. Thus, new insights into GCRV pathogenesis and antiviral treatment strategies are provided.

Circ-Udg Derived from Cyprinid Herpesvirus 2 Promotes Viral Replication

Cyprinid herpesvirus 2 (CyHV-2) has caused great losses to the gibel carp (Carassius auratus gibelio) industry. Previous studies showed that certain DNA viruses can encode circular RNAs (circRNAs) to regulate virus infection, which provides new clues for the treatment of viral disease. Whether CyHV-2 can encode circRNAs is still unknown. Here, 10 CyHV-2-derived circRNAs were identified, and the function of circ-udg, a circRNA derived from the CyHV-2 uracil DNA glycosylase (udg) gene, was studied. Although the expression level of circ-udg was lower than that of the parental gene, udg, its expression level was elevated in tandem with the proliferation of CyHV-2 and udg. In vitro experiments confirmed that circ-udg could promote the proliferation of CyHV-2. Moreover, circ-udg could encode a truncated UDG protein consisting of 147-amino-acid residues (termed circ-udg-P147). Both UDG and circ-udg-P147 were found to promote CyHV-2 proliferation, but the promoting effect of circ-udg on CyHV-2 proliferation was attenuated after circ-udg lost the ability to encode circ-udg-P147. Also, circ-udg-P147 could not change the transcription level of the udg gene. Interestingly, the UDG protein level was increased by circ-udg-P147. These results deepen the understanding of the genetic information carried by the genome of CyHV-2 and provide a new target for the treatment of gibel carp bleeding disease caused by CyHV-2. IMPORTANCE The outbreak of C. auratus gibelio gill hemorrhagic disease caused by CyHV-2 brought great losses to the gibel carp industry. Therefore, exploring the interaction between CyHV-2 and host and the molecular mechanism of viral infection is of great significance in preventing and treating the gibel carp gill hemorrhagic disease. Although some progress has been made in the study of CyHV-2, the mechanism of interaction between CyHV-2 and crucian carp is still unclear. In this study, we found that CyHV-2 can encode circRNA to regulate virus replication. Our study provides novel information on CyHV-2 functional genomics, a reference for research into the circRNA of other viruses, and theoretical guidance for preventing and treating gibel carp bleeding disease.

[Analysis on the epidemic characteristics and genetic characteristics of varicella in Beijing from 2019 to 2021]

The reported incidence of varicella in Beijing from 2019 to 2021 were 63.8/100 000, 32.2/100 000 and 38.6/100 000, respectively. There were two VZV epidemics in Beijing each year, one peaked in May and the other in November. However, the first VZV epidemic almost disappeared in 2020. Among the cases involved in the varicella outbreaks in school, the proportion of the students with no history of vaccine immunization, 1 dose of immunization and 2 doses of immunization were 33.12%, 44.79% and 22.08%, respectively. The major body of VZV breakthrough cases was children aged 6-14 years (523/755, 69.27%). The proportion of moderate- or severe-rash were 55.32%, 39.06%, 29.96% in the three groups of cases with no immunization history, 1 dose of immunization and 2 doses of immunization, respectively (P<0.001). A total of 1 089 varicella samples were collected, and 837 (76.86%) were confirmed to be PCR-positive for VZV and were identified as VZV wild strains. 311 VZV strains were sequenced successfully, and 307 strains were clade 2 (98.72%), 1 clade 3 (0.32%) and 3 Clade 5 (0.96%). Compared with the representative strains, the nucleotide similarities of ORF22 fragments were between 99.4% and 100%, and amino acid similarities were between 99.4% and 100%.

[Effectiveness of statin treatment strategies for primary prevention of cardiovascular diseases in a community-based Chinese population: A decision-analytic Markov model]

OBJECTIVE

To evaluate the effectiveness of statin treatment strategies based on risk assessment for the primary prevention of cardiovascular diseases by the Western guidelines in a community-based Chinese population from economically developed areas using data from the Chinese electronic health records research in Yinzhou (CHERRY) study.

METHODS

A Markov model was used to evaluate the effectiveness of the following statin treatment strategies, including: (1) usual care without cardiovascular risk assessment(Strategy 0); (2) using the World Health Organization (WHO) non-laboratory-based risk charts with statin treatment for high-risk group (risk ≥ 20%) (Strategy 1); (3) using the WHO laboratory-based risk charts with statin treatment for high-risk group (risk ≥ 20%) (Strategy 2); and (4) using the Prediction for Atherosclerotic cardiovascular disease Risk in China (China-PAR) model with statin treatment for high-risk group (risk ≥ 10%, Strategy 3). According to the guidelines, adults in the medium-risk group received lifestyle intervention, and adults in the high-risk group received life-style intervention and statin treatment under these strategies. The Markov model simulated different strategies for ten years (cycles) using parameters from the CHERRY study, published data, meta-analyses and systematic reviews for Chinese. The number of cardiovascular events or deaths, as well as the number need to treat (NNT) with statin per cardiovascular event or death prevented, were calculated to compare the effectiveness of different strategies. One-way sensitivity analysis on the uncertainty of incidence rate of cardiovascular diseases, and probabilistic sensitivity analysis on the uncertainty of hazard ratios of interventions were conducted.

RESULTS

Totally 225 811 Chinese adults aged 40-79 years without cardiovascular diseases at baseline were enrolled. In contrast to the usual care without risk assessment-based statin treatment strategy, Strategy 1 using the WHO non-laboratory-based risk charts could prevent 3 482 [95% uncertainty interval (UI): 2 110-4 661] cardiovascular events, Strategy 2 using the WHO laboratory-based risk charts could prevent 3 685 (95%UI: 2 255-4 912) events, and Strategy 3 using the China-PAR model could prevent 3 895 (95%UI: 2 396-5 181) events. NNTs with statin per cardiovascular event prevented were 22 (95%UI: 14-54), 21 (95%UI: 14-52), and 27 (95%UI: 17-67), respectively. Strategy 3 could prevent more cardiovascular events, while Strategies 1 and 2 required fewer numbers need to treat with statin per cardiovascular event prevented. The results were consistent in the sensitivity analyses.

CONCLUSION

The statin treatment strategies based on risk assessment for the primary prevention of cardiovascular diseases recommended by the Western guidelines could achieve substantive health benefits in adults from developed areas of China. Using the China-PAR model for cardiovascular risk assessment could prevent more cardiovascular diseases while using the WHO risk charts seems more efficient.

Viral Circular RNAs and Their Possible Roles in Virus-Host Interaction

Circular RNAs (circRNAs) as novel regulatory molecules have been recognized in diverse species, including viruses. The virus-derived circRNAs play various roles in the host biological process and the life cycle of the viruses. This review summarized the circRNAs from the DNA and RNA viruses and discussed the biogenesis of viral and host circRNAs, the potential roles of viral circRNAs, and their future perspective. This review will elaborate on new insights gained on viruses encoded circRNAs during virus infection.

[Applications of microsimulation model for cost-effectiveness analysis on screening in epidemiology]

Microsimulation model simulates individuals and estimates transition probabilities within the population using individual participant data. This approach could deal with the heterogeneous characteristics among the people or personal history of diseases and may be relevant in addressing cost-effectiveness problems of screening for complex conditions in epidemiology. This paper introduces the general principles, basic steps involved in implementation, analytic methods, and other related issues of the microsimulation model. Based on a practical research case of estimating the cost-effectiveness of microalbuminuria screening for chronic kidney disease in the United States, critical points in applications of the microsimulation model for cost-effectiveness analysis of screening were discussed in detail, including model development, model analysis, and the interpretation of the results. The microsimulation model considers the dynamic nature of complex diseases by estimating a broad range of individual characteristics and increasingly used to provide insights into complex problems that the Markov model does not efficiently address. For better supporting evidence-informed decision-making in public health, future studies should be aware of the accuracy of parameters in the decision-analytic model and the transparency of the models and results, as well as complying with the relevant reporting standards.

AB0894 Efficacy and Safety of Guselkumab in Biologic-Naïve Patients With Active Axial Psoriatic Arthritis: Study Design of a Phase 4, Randomized, Double-Blind, Placebo-Controlled Trial

Background

Established criteria for classifying axial psoriatic arthritis (PsA) are lacking, and assessments of disease activity often rely on measures developed for ankylosing spondylitis (AS). There is an unmet need to systematically identify and measure efficacy of treatments for axial PsA patients (pts). Guselkumab (GUS), a selective interleukin (IL)-23p19 inhibitor, was efficacious in improving signs and symptoms of active PsA in 2 phase 3, randomized, placebo (PBO)-controlled studies: DISCOVER-1 and DISCOVER-2. In a post-hoc pooled analysis of DISCOVER-1&2 pts with investigator-confirmed sacroiliitis, GUS-treated pts had greater improvements in axial symptoms compared with PBO.1 Imaging in DISCOVER-1&2 was restricted to the sacroiliac (SI) joints, occurring prior to/at screening as confirmed by the investigator, and locally read.

Objectives

To design a new, dedicated study to evaluate the effects of GUS on axial PsA prospectively.

Methods

Cumulative evidence from DISCOVER-1&2, including exposure–response relationship, covariate adjustment for modest baseline imbalances across treatment groups, subgroup analyses, and comparisons within and across these studies, was considered in designing a new trial. Data from the pivotal registrational studies suggest similar efficacy with GUS every-4-weeks (Q4W) and Q8W regimens in treating PsA signs and symptoms, including symptoms of axial involvement. Power calculations were based on mean changes in Bath AS Disease Activity Index (BASDAI) scores in DISCOVER-1&2.

Results

The phase 4, randomized, PBO-controlled STAR study is specifically designed to prospectively assess efficacy outcomes in PsA pts with magnetic resonance imaging (MRI)-confirmed axial inflammation. Based on observed mean (SD) changes from baseline in BASDAI score from DISCOVER-1&2 (Table 1), 405 pts, randomized (1:1:1) to GUS Q4W, GUS at W0, W4, then Q8W, or PBO →GUS Q8W at W24, are planned for enrollment (Figure 1). STAR eligibility criteria include PsA ≥6 months and active disease (≥3 swollen & ≥3 tender joints, C-reactive protein [CRP] ≥0.3 mg/dL) despite prior non-biologic disease-modifying antirheumatic drugs, apremilast, and/or non-steroidal anti-inflammatory drugs. Pts will be naïve to biologics and Janus kinase inhibitors and have BASDAI ≥4, spinal pain score (visual analog scale [VAS]) ≥4, and screening MRI-confirmed axial disease (positive spine and/or SI joints defined as centrally read Spondyloarthritis Research Consortium of Canada [SPARCC] score ≥3). Follow-up MRIs of spine and SI joints will be obtained at W0, W24, and W52 and also centrally read, with readers blinded to treatment group and timepoint. Spinal/SI joint inflammation will be scored using the SPARCC method, with the former also assessed using the CAN-DEN method. The primary endpoint is mean change in BASDAI at W24; controlled (hierarchical) secondary endpoints, all at W24, include AS Disease Activity Score (ASDAS-CRP), Disease Activity Index for PsA (DAPSA), ≥40% improvement in Assessment in AS criteria (ASAS40), and mean changes in spine/SI joint SPARCC scores.

Table 1.

Power calculations for the primary endpoint in the Phase 4 STAR study.

Historical trial data*Observed mean (SD) change in BASDAI from W0-24Effect sizePower(N=135; α=0.05)**PBO-1.28 (2.24)GUS 100 mg Q4W-2.51 (2.00)1.23>99%GUS 100 mg Q8W-2.61 (2.47)1.33>99%

* From the pooled DISCOVER-1&2 trials

**Power calculations based on N=135 per study group (1:1:1 randomization) and 2-sided significance of 0.05 using a 2-sample T-test assuming equal variances

BASDAI, Bath Ankylosing Spondylitis Disease Activity Index; GUS, guselkumab; PBO, placebo; Q4W, every 4 weeks; Q8W, every 8 weeks; SD, standard deviation; W, week

Conclusion

The phase 4 STAR study will allow for an in-depth, prospective evaluation of the effects of selectively inhibiting the IL-23p19 subunit with GUS in pts with MRI-confirmed axial PsA.

References

[1]Mease, et al. Lancet Rheum. 2021;3(10):e715-e723.

Disclosure of Interests

Dafna D Gladman Consultant of: AbbVie, Amgen, BMS, Celgene, Eli Lilly, Galapagos, Gilead, Janssen, Novartis, Pfizer, and UCB, Grant/research support from: Abbvie, Amgen, BMS, Celgene, Eli Lilly, Janssen, Novartis, Pfizer and UCB, Philip J Mease Speakers bureau: AbbVie, Aclaris, Amgen, BMS, Celgene, Eli Lilly, Galapagos, Gilead, GlaxoSmithKline, Inmagene, Janssen, Novartis, Pfizer, SUN Pharma, and UCB, Consultant of: AbbVie, Aclaris, Amgen, BMS, Celgene, Eli Lilly, Galapagos, Gilead, GlaxoSmithKline, Inmagene, Janssen, Novartis, Pfizer, SUN Pharma, and UCB, Grant/research support from: AbbVie, Aclaris, Amgen, BMS, Celgene, Eli Lilly, Galapagos, Gilead, GlaxoSmithKline, Inmagene, Janssen, Novartis, Pfizer, SUN Pharma, and UCB, Paul Bird Speakers bureau: AbbVie, Eli Lilly, Gilead, Janssen, MSD, Pfizer, and UCB, Consultant of: Eli Lilly, Gilead, Janssen, Novartis, and Pfizer, Enrique Soriano Speakers bureau: AbbVie, Amgen, Bristol Myers Squibb, Eli Lilly, Janssen, Novartis, Pfizer, Roche, and UCB, Consultant of: AbbVie, Janssen, Novartis, and Roche, Grant/research support from: AbbVie, Janssen, Novartis, Pfizer, Roche, and UCB, Soumya D Chakravarty Shareholder of: Johnson & Johnson, Employee of: Janssen Scientific Affairs, LLC, May Shawi Shareholder of: Johnson & Johnson, Employee of: Janssen Global Services, LLC, Stephen Xu Shareholder of: Johnson & Johnson, Employee of: Janssen Research & Development, LLC, Sean Quinn Shareholder of: Johnson & Johnson, Employee of: Janssen Scientific Affairs, LLC, Cinty Gong Shareholder of: Johnson & Johnson, Employee of: Janssen Scientific Affairs, LLC, Evan Leibowitz Shareholder of: Johnson & Johnson, Employee of: Janssen Scientific Affairs, LLC, Lai-Shan Tam Consultant of: Janssen, Pfizer, Sanofi, AbbVie, Boehringer Ingelheim, and Lilly, Grant/research support from: Amgen, Boehringer Ingelheim, Janssen, GSK, Novartis, and Pfizer, Philip Helliwell Speakers bureau: AbbVie, Janssen, and Novartis, Consultant of: Galapagos and Janssen, Grant/research support from: AbbVie, Janssen, and Pfizer, Arthur Kavanaugh Consultant of: Abbvie, Amgen, Bristol Myers Squibb, Eli Lilly, Gilead, Janssen, Novartis, Pfizer, and UCB, Atul Deodhar Speakers bureau: AbbVie, Eli Lilly, Janssen, Novartis, Pfizer, and UCB, Consultant of: AbbVie, Amgen, Boehringer Ingelheim, Bristol Myers Squibb, Celgene, Eli Lilly, Galapagos, Glaxo Smith & Kline, Janssen, Novartis, Pfizer, and UCB, Grant/research support from: AbbVie, Eli Lilly, Glaxo Smith & Kline, Novartis, Pfizer, and UCB, Mikkel Østergaard Speakers bureau: AbbVie, Boehringer-Ingelheim, Bristol Myers Squibb, Celgene, Eli-Lilly, Hospira, Janssen, Merck, Novartis, Novo, Orion, Pfizer, Regeneron, Roche, Sandoz, Sanofi, and UCB, Consultant of: AbbVie, Boehringer-Ingelheim, Bristol Myers Squibb, Celgene, Eli-Lilly, Hospira, Janssen, Merck, Novartis, Novo, Orion, Pfizer, Regeneron, Roche, Sandoz, Sanofi, and UCB, Grant/research support from: AbbVie, Bristol Myers Squibb, Celgene, and Novartis, Xenofon Baraliakos Speakers bureau: AbbVie, Biocad, Chugai, Eli Lilly, Janssen, MSD, Novartis, Pfizer, Roche, and UCB, Consultant of: AbbVie, Biocad, Chugai, Eli Lilly, Janssen, MSD, Novartis, Pfizer, Roche, and UCB, Grant/research support from: AbbVie, Biocad, Chugai, Eli Lilly, Janssen, MSD, Novartis, Pfizer, Roche, and UCB

Micropeptide vsp21 translated by Reovirus circular RNA 000048 attenuates viral replication

To date, some DNA viruses and single-stranded RNA viruses have been found to generate circRNAs. However, the reports on circRNAs produced by double-stranded RNA viruses are very limited. In this study, Bombyx mori cypovirus (BmCPV), a typical double-stranded RNA virus belonging to the Reoviridae, was demonstrated to generate viral circRNAs (vcircRNAs) and a vcircRNA_000048 whose sequence corresponds with the region 164-1245 nt on the BmCPV genomic dsRNA S5 segment (GQ294468.1) was validated by PCR, Sanger sequencing, reverse transcription-rolling circle amplification, and Northern blotting. Furthermore, we verified that vcircRNA_000048 translates a micropeptide vsp21 with 21 amino acid residues in an IRES-dependent manner, and vsp21 attenuates the viral replication. These findings provided a novel clue to understanding the regulation of viral multiplication and interaction of reovirus with the host.

β-Arrestin 2 acts an adaptor protein that facilitates viral replication in silkworm

β-Arrestin 2 is known to be a widely distributed adaptor protein in mammals but its function has never been reported in Lepidoptera insects. Herein, the β-Arrestin 2 (BmArrestin 2) gene from silkworm was cloned and characterized. The spatiotemporal expression level of BmArrestin 2 was highest in the gonads at the 3rd day of 5th instar, whereas the highest and lowest abundance of BmArrestin 2 were identified in the tracheal and testis, respectively. BmArrestin 2 is mainly distributed in the cytoplasm. Furthermore, in BmN cells，overexpression of BmArrestin 2 promoted Bombyx mori nucleopolyhedrovirus (BmNPV) and B. mori cytoplasmic polyhedrosis virus (BmCPV) replication as the increment of the concentration of plasmid transfection, whereas silencing the gene with specific siRNA inhibited viral replication. Replication of BmNPV and BmCPV also was weakened using BmArrestin 2 antiserum as the increment of the concentration. Immunofluorescent staining revealed the invasion of recombinant BmNPV or BmCPV was decreased after blocking endogenous BmArrestin 2. On the other hand, BmArrestin 2 co-localizes with recombinant BmNPV and BmCPV virions in BmN cells. These results suggest that BmArrestin 2 may represent a novel target for antiviral strategies, as it is an adaptor protein that plays a key role in virus replication.

Identification and Characterization of BmNPV m6A Sites and Their Possible Roles During Viral Infection

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the most serious pathogens and causes serious economic losses in sericulture. At present, there is no epigenetic modification of BmNPV transcripts, especially of m6A, and this modification mediates diverse cellular and viral functions. This study showed that m6A modifications are widespread in BmNPV transcripts in virally infected cells and the identified m6A peaks with a conserved RRACH sequence. m6A sites predominantly appear in the coding sequences (CDS) and the 3'-end of CDS. About 37% of viral genes with m6A sites deleted from the viral genome did not produce any infectious virions in KOV-transfected cells. Among the viral genes related to replication and proliferation, ie-1 mRNA was identified with a higher m6A level than other viral genes. The m6A sites in the ie-1 mRNA may be negatively related to the protein expression. Viral replication was markedly inhibited in cells overexpressed with BmYTHDF3 in a dose-dependent manner, and a contrary effect was found in si-BmYTHDF3-transfected cells. Collectively, the identification of putative m6A modification in BmNPV transcripts provides a foundation for comprehensively understanding the viral infection, replication, and pathobiology in silkworms.

BmCPV-Derived Circular DNA vcDNA-S7 Mediated by Bombyx mori Reverse Transcriptase (RT) Regulates BmCPV Infection

Circular DNAs derived from single-stranded RNA viruses play important roles in counteracting viral infection. However, whether double-stranded RNA viruses generate functional circular DNAs is still unknown. Using circDNA sequencing, divergent PCR, DNA in situ hybridization and rolling circular amplification, we presently confirmed that in silkworm, Bombyx mori cytoplasmic polyhedrosis virus (BmCPV), a double-stranded RNA virus belonging to cypovirus, is prone to produce a BmCPV-derived circular DNA termed as vcDNA-S7. We have also found that vcDNA-S7 formation is mediated by endogenous reverse transcriptase (RT), and the proliferation of BmCPV can be inhibited by vcDNA-S7 in vitro and in vivo. Moreover, we have discovered that the silkworm RNAi immune pathway is activated by vcDNA-S7, while viral small interfering RNAs (vsiRNAs) derived from transcribed RNA by vcDNA-S7 can be detected by small RNA deep sequencing. These results suggest that BmCPV-derived vcDNA-S7, mediated by RT, can serve as a template for the biogenesis of antiviral siRNAs, which may lead to the repression of BmCPV infection. To our knowledge, this is the first demonstration that a circular DNA, produced by double stranded RNA viruses, is capable of regulating virus infection.

Bombyx mori cypovirus (BmCPV) induces PINK1-Parkin mediated mitophagy via interaction of VP4 with host Tom40

The mechanism by which infection by Bombyx mori cytoplasmic nucleopolyhedrosis virus (BmCPV) causes autophagy has not been studied in detail. Herein we have observed by electron microscopy that infection with BmCPV causes autophagosome and mitochondrial structure damage in Bombyx mori midgut. In BmN cells infected with BmCPV and expressing eGFP-LC3, fluorescence spots and LC3-II levels increased, suggesting that BmCPV infection causes autophagy. Autophagy inducer rapamycin (Rap) and autophagy inhibitor 3-methyladenine (3-MA) were used to monitor the effects of mitophagy on BmCPV proliferation. It was found BmCPV proliferation to be promoted by mitophagy. Transient transfection experiments in cultured BmN cells showed that mitophagy can be triggered by expression of BmCPV structural protein VP4. Moreover, VP4 caused upregulation of p-Drp1, PINK1 and Parkin proteins in the mitophagy pathway and downregulation of mitochondrial membrane protein Tom20. Furthermore, interaction between VP4 with Tom40 was confirmed by Co-IP, western blot and colocalization experiment, and overexpression of Tom40 reduce the level of mitochondrial autophagy induced by VP4. These results suggested that VP4 induced PINK1-Parkin-mediated mitophagy interacting with Tom40. These findings deepen our understanding of the interaction between BmCPV and silkworm and also provide a molecular target for screening anti-BmCPV drugs.

Bombyx mori Akirin hijacks a viral peptide vSP27 encoded by BmCPV circRNA and activates the ROS-NF-κB pathway against viral infection

Bombyx mori cypovirus (BmCPV), a member of the family Reoviridae, is a model of Cypovirus, has a 10 segmented double-stranded RNA genome. However, so far, only one viral small peptide vSP27 with negative regulation on viral infection was identified; the mechanisms underlying host-BmCPV interaction are still unknown. Here, we identified that vSP27 was translated from a BmCPV derived circular RNA (circRNA-vSP27). Subsequently, results showed that vSP27 induced generation of ROS activated the NF-κB signaling pathway, induced the expression of antimicrobial peptides, and suppressed BmCPV infection. On the other hand, we identified a nuclear protein Akirin that could hijack vSP27, positively regulate the NF-κB pathway, and lead to inhibiting the viral infection. Altogether, our data suggested that BmCPV derived circRNA-vSP27 with small peptide translation activity may be employed by the host immunity in defense against the BmCPV infection.

Response to Bombyx mori nucleopolyhedrovirus infection in silkworm: Gut metabolites and microbiota

The diversity of microbiota and metabolites in the digestive tract gut is important in physiology and homeostasis, nutrient uptake and virus infection. In lepidopteran insect model silkworms, little is known about how microbiota and metabolites are altered after oral infection with BmNPV. Herein, we used 16S rDNA sequencing and metabolomics to show that the gut microbiota and metabolites of silkworm midgut are significantly altered after BmNPV infection. Kyoto Encyclopedia of Genes and Genomes analysis revealed enrichment of flavone and flavonol biosynthesis, glycosyltransferases, and electron transfer carriers signaling pathways, suggesting potential roles in viral infection. Infection also changed the abundance of metabolites in the digestive tract gut, where most pathways were related to metabolism of amino acids, fatty acids and other pathways (e.g., citrate cycle). In addition, a correlation was observed between digestive tract gut microbiota and metabolites. These results shed light on the interaction between digestive tract gut microbiota, metabolites and host-virus interaction, and enhance our understanding of viral infection links to midgut microbiota and metabolic activity in the silkworm.

Hepatocellular carcinoma progression mediated by hepatitis B virus-encoded circRNA HBV_circ_1 through interaction with CDK1

Hepatitis B virus (HBV) produces circular RNA (circRNA), whose functions have not yet been clearly elucidated. In this study, a novel circRNA HBV_circ_1 produced by HBV was identified in HBV-positive HepG2.2.15 cells and HBV-related hepatocellular carcinoma (HCC) tissue (HCCT). Microarray analysis of 68 HCCT samples showed that HBV_circ_1 abundance was significantly higher than that in paracancerous tissues. In addition, survival rate of HBV_circ_1-positive patients was significantly lower compared with HBV_circ_1-negative patients. Transient expression indicated that HBV_circ_1 enhanced cell proliferation, migration, and invasion and inhibited apoptosis in vitro. Furthermore, ectopical HBV_circ_1 expression increased tumor size in vivo. HBV_circ_1 was confirmed to interact with cyclin-dependent kinase 1 (CDK1) to regulate cell proliferation. These results suggest that HCC progression may be promoted by interaction of HBV_circ_1 with CDK1. Our data not only showed a novel clue to understand carcinogenesis and progress of HBV-related HCC but also provided a new target for the development of therapeutic drugs.

[Epidemiological characteristics of Chlamydia pneumoniae in cases with acute respiratory infection in Beijing, 2015-2019]

Objective: To analyze the epidemiological characteristics of Chlamydia pneumoniae infection among patients with acute respiratory infection in Beijing from 2015 to 2019. Methods: The epidemiological data of acute respiratory infection patients from 35 sentinel hospitals in Beijing were collected by the respiratory pathogen surveillance system in Beijing. The clinical samples were collected to detect Chlamydia pneumoniae, and the sequence of the VD4 region of the ompA gene in positive samples was analyzed. Results: From January 2015 to December 2019, the overall positive rate of Chlamydia pneumoniae among patients with acute respiratory infection in Beijing was 0.34% (129/37 460). The positive rate of Chlamydia pneumoniae generally increased in March, reaching the peak in May, and started to drop in July, with a duration of about 5-8 months. The epidemic season in different years fluctuated by 1-2 months. The positive monthly rate of Chlamydia pneumoniae was no less than 0.30% in every epidemic season. The positive rate of Chlamydia pneumoniae was the highest in the 5-44 years old group and the highest in 10-14 year-olds. The risk of Chlamydia pneumoniae infection increased with age in patients younger than 25 years old and decreased in those older one aged than 25 years of age. The positive rates in male and female patients were 0.33% (68/20 830) and 0.37% (61/16 528), respectively, and there was no significant difference between the two groups (χ²=0.486, P=0.486). The positive rate of Chlamydia pneumoniae in patients with common pneumonia was higher than that in patients with upper pneumonia and severe pneumonia (χ²=36.797, P<0.01). Other respiratory pathogens were also detected in the Chlamydia pneumoniae samples, and the top four pathogens appeared as Haemophilus influenzae (15 cases), Streptococcus pneumoniae (13 cases), Rhinovirus (8 cases), and Stenotrophomonas maltophilia (7 cases). 101 strains of 129 Chlamydia pneumoniae positive samples were identified as type A by sequencing. Conclusions: The annual epidemic pattern of Chlamydia pneumoniae in Beijing, is unimodal, and the epidemic season generally appears from March to July. The seasonal characteristics of Chlamydia pneumoniae in Beijing can be used for the differential diagnosis of Chlamydia pneumoniae from other respiratory pathogens. Chlamydia pneumoniae is most common in people aged 5-44 years, and the primary genotype is type A. People aged 10-44 years old suffer the highest incidence. If the nucleic acid positive rate of Chlamydia pneumoniae exceeds 0.30% for two consecutive months, the high prevalence period of Chlamydia pneumoniae can be preliminarily expected. Chlamydia pneumoniae infection has a higher probability of progressing to severe pneumonia from general pneumonia.

[Clinical characteristics of patients infected with Chlamydia pneumoniae in Beijing from 2015 to 2019]

Objective: To aralyze the clinical characteristics of Chlamydia pneumoniae infection in Beijing. Methods: Based on Beijing Respiratory Pathogen Surveillance System, acute respiratory infection patients were recruited from 35 different types of sentinel hospitals in Beijing. Their epidemiological and clinical data were systematically collected and clinical specimens were also obtained. Nuclear acid testing was performed for 30 types of respiratory pathogens (including Chlamydia pneumoniae). The identified patients of Chlamydia pneumoniae infection were divided into two groups, the acute upper respiratory tract infection (AURI) group and pneumoniae group. The differences in clinical characteristics, laboratory examination and prognosis were compared by using independent sample t test, Mann-Whitney U test, χ² test or Fisher's exact probability test. Results: A total of 119 patients of Chlamydia pneumoniae infection were included, including 12 patients in the AURI group and 107 patients in pneumoniae group. Chlamydia pneumoniae infection mainly occurred in people aged from 5 to 44 years, accounting for 81.5% (97/119). The three most common clinical symptoms were cough (92.4%, 110/119), fever (88.8%, 95/107), and sputum production (76.5%, 91/119). White blood cell counts increased in 39.3% (46/117) of patients. Neutrophile granulocyte proportion increased in 39.7% (46/116) of patients. Platelet count increased in 36.9% (41/111) of patients. An increase of the creatine kinase MB isoenzyme (CKMB) was observed in 12 pneumonia patients (24.5%, 12/49). Radiological examination showed that 90.6% (87/105) of patients in the pneumoniae group had pulmonary parenchymal changes; the lesion occurred most commonly in the lower right lung lobe (34.3%, 36/105) and the lower left lung lobe (27.6%, 29/105). Although 73.8% (79/107) of patients in the pneumoniae group were hospitalized, no case received intensive care unit or mechanical ventilation. As to outcomes, one patient developed respiratory failure and 6 patients suffered myocardial injury. No death was observed in this study. The median days of hospitalization and course of illness for pneumonia patients M(P₂₅,P₇₅) were 10.0 (7.0, 13.0) days and 18.0 (13.5, 22.0) days, respectively. Conclusion: Generally, Chlamydia pneumoniae infections in Beijing from 2015 to 2019 were mild, and the main clinical manifestations were cough, fever and sputum. However, most patients in the pneumoniae group caused by Chlamydia pneumoniae still required hospitalization but with a better outcome.

Proteomic analysis of the exosomes secreted from Ctenopharyngodon idellus kidney cells infected with grass carp reovirus reveals their involvement in the cellular responses to viral infection

Exosomes are small membrane-enclosed vesicles secreted by various types of cells. Exosomes not only participate in different physiological processes in cells, but also involve in the cellular responses to viral infection. Grass carp reovirus (GCRV) is a non-enveloped virus with segmented, double-stranded RNA genome. Nowadays, the exact role of exosomes in regulating the life cycle of GCRV infection is still unclear. In this study, the exosomes secreted from Ctenopharyngodon idellus kidney (CIK) cells infected or uninfected with GCRV were isolated, and the differential protein expression profiles were analyzed by proteomic technologies. A total of 1297 proteins were identified in the isolated exosomes. The differentially abundant proteins were further analyzed with functional categories, and numerous important pathways were regulated by exosomes in GCRV-infected CIK cells. These exosomal proteins were estimated to interact with the genes (proteins) of the top 10 most enriched signaling pathways. Furthermore, GW4869 exosome inhibitor suppressed the expression level of VP7 in GCRV-infected cells, suggesting that exosomes play a crucial role in the life cycle of GCRV infection. These findings could shed new lights on understanding the functional roles of exosomes in the cellular responses to GCRV infection.

Metatranscriptomic Analysis Reveals an Imbalance of Hepatopancreatic Flora of Chinese Mitten Crab Eriocheir sinensis with Hepatopancreatic Necrosis Disease

Simple Summary

The cause of Chinese mitten crab Eriocheir sinensis hepatopancreas necrosis disease (HPND) remains a mystery. In this study, metatranscriptomics sequencing was conducted to characterize the changes in the structure and gene expression of hepatopancreatic flora of crabs with and without typical symptoms of HPND; an imbalance of hepatopancreatic flora can be found in the crab with HPND, and the detected microbial taxa decreased, whereas the prevalence of Spiroplasma eriocheiris significantly increased in the hepatopancreatic flora of crabs with typical symptoms of HPND, and the relative abundances of the virus and microsporidia in crabs with HPND were very low and did not increase with disease progression. The differentially-expressed genes (DEGs) in hepatopancreatic flora between crabs with and without HPND were enriched ribosome, retinol metabolism, and biosynthesis of unsaturated fatty acid KEGG pathways. These results suggested that an imbalance of hepatopancreatic flora was associated with crab HPND, and the enriched pathways of DEGs were associated with the pathological mechanism of HPND.

Abstract

Hepatopancreas necrosis disease (HPND) of the Chinese mitten crab Eriocheir sinensis causes huge economic loss in China. However, the pathogenic factors and pathogenesis are still a matter of dissension. To search for potential pathogens, the hepatopancreatic flora of diseased crabs with mild symptoms, diseased crabs with severe symptoms, and crabs without visible symptoms were investigated using metatranscriptomics sequencing. The prevalence of Absidia glauca and Candidatus Synechococcus spongiarum decreased, whereas the prevalence of Spiroplasma eriocheiris increased in the hepatopancreatic flora of crabs with HPND. Homologous sequences of 34 viral species and 4 Microsporidian species were found in the crab hepatopancreas without any significant differences between crabs with and without HPND. Moreover, DEGs in the hepatopancreatic flora between crabs with severe symptoms and without visible symptoms were enriched in the ribosome, retinol metabolism, metabolism of xenobiotics by cytochrome P450, drug metabolism—cytochrome P450, biosynthesis of unsaturated fatty acids, and other glycan degradation. Moreover, the relative abundance of functions of DEDs in the hepatopancreatic flora changed with the pathogenesis process. These results suggested that imbalance of hepatopancreatic flora was associated with crab HPND. The identified DEGs were perhaps involved in the pathological mechanism of HPND; nonetheless, HPND did not occur due to virus or microsporidia infection.

Incidence of Carassius auratus Gibelio Gill Hemorrhagic Disease Caused by CyHV-2 Infection Can Be Reduced by Vaccination with Polyhedra Incorporating Antigens

Encapsulation of antigens within protein microcrystals (polyhedra) is a promising approach for the stable delivery of vaccines. In this study, a vaccine was encapsulated into polyhedra against cyprinid herpesvirus II (CyHV-2). CyHV-2 typically infects gibel carp, Carassius auratus gibelio, causing gill hemorrhagic disease. The vaccine was constructed using a codon-optimized sequence, D4ORF, comprising the ORF72 (region 1–186 nt), ORF66 (region 993–1197 nt), ORF81 (region 603–783 nt), and ORF82 (region 85–186 nt) genes of CyHV-2. The H1-D4ORF and D4ORF-VP3 sequences were, respectively, obtained by fusing the H1-helix sequence (region 1–90 nt) ofBombyx mori cypovirus(BmCPV) polyhedrin to the 5′ terminal end of D4ORF and by fusing a partial sequence (1–279 nt) of the BmCPV VP3 gene to the 3′ terminal end of D4ORF. Furthermore, BmNPV-H1-D4ORF-polh and BmNPV-D4ORF-VP3-polh recombinant B. mori nucleopolyhedroviruses (BmNPVs), belonging to the family Baculoviridae, and co-expressing BmCPV polyhedrin and H1-D4ORF or D4ORF-VP3, were constructed. H1-D4ORF and D4ORF-VP3 fusion proteins were confirmed to be encapsulated into recombinant cytoplasmic polyhedra by Western blotting. Degradation of vaccine proteins was assessed by SDS-PAGE, and the results showed that the encapsulated vaccine proteins in polyhedra could be protected from degradation. Furthermore, when gibel carp were vaccinated with the purified polyhedra from BmNPV-H1-D4ORF-polh and BmNPV-D4ORF-VP3-polh via injection, the antibody titers in the serum of the vaccinated fish reached 1:6400–1:12,800 at 3 weeks post-vaccination. Therelative percentage of survival of immunized gibel carp reached 64.71% and 58.82%, respectively, following challenge with CyHV-2. These results suggest that incorporating vaccine protein into BmCPV polyhedra may be a novel approach for developing aquaculture microencapsulated vaccines.

Characterization of the lipidomic profile of BmN cells in response to Bombyx mori cytoplasmic polyhedrosis virus infection

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV)that belongs to the genus Cypovirus in the family of Reoviridae is one of the problematic pathogens in sericulture. In our previous study, we have found that lipid-related constituents in the host cellular membrane are associated with the BmCPV life cycle. It is important to note that the lipids not only affect the cellular biological processes, they also impact the virus life cycle. However, the intracellular lipid homeostasis in BmN cells after BmCPV infection remains unclear. Here, the lipid metabolism in BmCPV-infected BmN cells was studied by lipidomics analysis. Our results revealed that the intracellular lipid homeostasis was disturbed in BmN cells upon BmCPV infection. Major lipids constituents in cellular membrane were found to be significantly induced upon BmCPV infection, which included triglycerides, phosphatidylcholine, phosphatidylethanolamine, sphingomyelin, phospholipids, glucoside ceramide, monoetherphosphatidylcholin, ceramide, ceramide phosphoethanolamine and cardiolipin. Further analysis of the pathways related to these altered lipids (such as PE and PC) showed that glycerophospholipid metabolism was one of the most enriched pathways. These results suggested that BmCPV may manipulate the lipid metabolism of cells for their own interest. The findings may facilitate a better understanding of the roles of lipid metabolic changes during virus infection in future studies.

circEgg regulates histone H3K9me3 by sponging bmo-miR-3391-5p and encoding circEgg-P122 protein in the silkworm, Bombyx mori

A large number of circular RNAs (circRNAs) have been found in different organisms; however, their function in the regulation of histone modification remains unknown. In this study, we found that the circRNA circEgg, cyclized by the 9th-13th exon of Bombyx mori histone-lysine N-methyltransferase eggless (BmEgg) gene, mainly distributes in the cytoplasm, its expression levels changed with silkworm developmental stages, and the linear transcript level of the BmEgg gene was decreased when circEgg was overexpressed. Moreover, circEgg was found to repress histone H3 lysine 9 methylation (H3K9me3), promote histone H3 lysine 9 acetylation (H3K9ac), and positively regulate histone deacetylase (HDAC) Rpd3 (BmHDAC Rpd3) gene expression by sponging the microRNA bmo-miR-3391-5p. Furthermore, circEgg encodes a circEgg-P122 protein which appears to inhibit H3K9me3. These results suggest that circEgg regulates histone modification by sponging bmo-miR-3391-5p and encoding circEgg-P122 protein. To our knowledge, this is the first report showing that a circRNA produced by BmEgg plays an important role in histone epigenetic modification.

AIV polyantigen epitope expressed by recombinant baculovirus induces a systemic immune response in chicken and mouse models

BACKGROUND

The protective efficacy of avian influenza virus (AIV) vaccines is unsatisfactory due to the presence of various serotypes generated by genetic reassortment. Thus, immunization with a polyantigen chimeric epitope vaccine may be an effective strategy for protecting poultry from infection with different AIV subtypes.

METHODS

Baculovirus has recently emerged as a novel and attractive gene delivery vehicle for animal cells. In the present study, a recombinant baculovirus BmNPV-CMV/THB-P10/CTLT containing a fused codon-optimized sequence (CTLT) of T lymphocyte epitopes from H1HA, H9HA, and H7HA AIV subtypes, and another fused codon-optimized sequence (THB) of Th and B cell epitopes from H1HA, H9HA, and H7HA AIV subtypes, driven by a baculovirus P10 promoter and cytomegalovirus CMV promoter, respectively, was constructed.

RESULTS

Western blotting and cellular immunofluorescence demonstrated that the CTLT (THB) can be expressed in rBac-CMV/THB-P10/CTLT-infected silkworm cells (mammalian HEK293T cells). Furthermore, the recombinant virus, rBac-CMV-THB-CTLT, was used to immunize both chickens and mice.

CONCLUSIONS

The results of an indirect ELISA, immunohistochemistry, and T lymphocyte proliferation assay indicated that specific humoral and cellular responses were detected in both chicken and mice. These results suggest that rBac-CMV/THB-P10/CTLT can be developed as a potential vaccine against different AIV subtypes.

Interleukin-17 suppresses grass carp reovirus infection in Ctenopharyngodon idellus kidney cells by activating NF-κB signaling

The grass carp accounts for a large proportion of aquacultural production in China, but the hemorrhagic disease caused by grass carp reovirus (GCRV) infection often causes huge economic losses to the industry. Interleukin 17 (IL-17) is an important cytokine that plays a critical role in the inflammatory and immune responses. Although IL-17 family members have been extensively studied in mammals, our knowledge of the activity of IL-17 proteins in teleosts in response to viral infection is still limited. In this study, the role of IL-17 in GCRV infection and its mechanism were investigated. The expression levels of IL-17AF1, IL-17AF2, and IL-17AF3 in Ctenopharyngodon idella kidney (CIK) cells gradually increased from 6 h after infection with GCRV. The nuclear translocation of p65, which acts in the NF-κB signaling pathway, was also increased by GCRV infection. The overexpression of IL-17AF1, IL-17AF2, or IL-17AF3 also promoted the nuclear translocation of p65 and the levels of phospho-IκBα in CIK cells, and reduced the expression of the viral structural protein VP7. An NF-κB signal inhibitor abolished the inhibition of GCRV infection by IL-17 proteins. These results suggested that the NF-κB signaling pathway was activated by the overexpression of IL-17 proteins, resulting in the inhibition of viral infection. In conclusion, in this study, we demonstrated that IL-17AF1, IL-17AF2, and IL-17AF3 acted as immune cytokines, exerting an antiviral effect by activating the NF-κB signaling pathway.

N 6-Methyladenosine Level in Silkworm Midgut/Ovary Cell Line Is Associated With Bombyx mori Nucleopolyhedrovirus Infection

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the most serious pathogens in sericulture and causes huge economic loss annually. The roles of N6-methyladenosine (m6A) modification in silkworms following BmNPV infection are currently unclear. Here, methylated RNA immunoprecipitation with next-generation sequencing were applied to investigate the m6A profiles in silkworm midgut following BmNPV infection. A total of 9144 and 7384 m6A peaks were identified from the BmNPV-infected (TEST) and uninfected silkworm midguts (CON), respectively, which were distributed predominantly near stop codons. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of common m6A peaks in nuclear genes revealed that these m6A-related transcripts were associated with crucial signaling pathways. Comparative transcriptome analysis showed that 1221 differential expressed m6A peaks were identified between TEST and CON, indicating that m6A modification is regulated following BmNPV infection. GO and KEGG pathway analysis of the differentially expressed m6A peaks showed their association with signal transduction, translation, and degradation. To understand further the effect of the m6A machinery on virus infection, expression levels of m6A-related genes were altered in silencing and overexpression experiments. Expression of viral structural protein VP39 was increased in BmN cells by siRNA-mediated depletion of methyltransferase-like (METTL) enzyme genes (BmMETTL3, BmMETTL14) and cytoplasmic YTH-domain family 3 (BmYTHDF3), while the reverse results were found after overexpression of the m6A-related enzymes in BmN cells. Overall, m6A modification might be a novel epigenetic mechanism that regulation BmNPV infection and interference with this mechanism may provide a novel antiviral strategy for preventing BmNPV disease.

Identification and characterization of novel type of RNAs, circRNAs in crucian carp Carassius auratus gibelio

Circular RNAs (circRNAs) with regulatory potency activity was identified from varieties of species. Crucian carp (Carassius auratus gibelio) is one of the most freshwater aquaculture species in China. Every year, huge economic damage to the farming was caused by the virus and bacterial infection. Until now, there is any information about circRNA reported from the Crucian carp. In this study, the expression pattern of circRNA in Crucian carp was investigated with transcriptomic analysis. The results showed that only 37 circRNAs were identified from the Crucian carp, and these circRNAs biogenesis was formed with canonical GU-AG splicing mechanism with unevenly distributed on the chromosomes. Wherein, most of the circRNAs were derived from the sense overlapping strategy. Reverse transcript PCR and Sanger sequencing data indicated that these circRNAs were existed authenticity in Crucian carp. The bioinformatics analysis indicated that circRNAs identified from the Crucian carp with potential miRNA sponge regulate the expression level of mRNAs. GO annotation and KEGG pathway analysis of these circRNAs showed that more than 20% circRNAs were related with catalytic activity and binding in the category of molecular function, and these circRNAs were enriched in 9 signaling pathways, such as, Wnt signaling pathway, MAPK signaling pathway, Ubiquitin mediated proteolysis et al. 220 mRNAs would be regulated by the circRNAs via miRNAs mediation. These target mRNAs were further analyzed with functional annotation and KEGG analysis. GO annotation analysis showed that several genes were related with function of nucleotide binding, transcription regulatory activity. KEGG pathway analysis showed that 5 genes were enriched in the pathway of Endocytosis. The circRNA-miRNA-mRNA regulation network indicated that one miRNA can link one or more circRNA and one or more mRNA. Overall, these results will not only help us to further understand the novel RNA transcripts in Crucian carp, but also provide the novel clue to investigate the interaction between host and pathogens from this novel circRNA molecule.

Recombinant baculovirus BacCarassius-D4ORFs has potential as a live vector vaccine against CyHV-2

Cyprinid herpesvirus II (CyHV-2) is highly contagious and pathogenic to Carassius auratus gibelio (gibel carp), causing enormous economic losses in aquaculture in Yancheng city, Jiangsu province, China; however, to date, there is no effective way to protect C. auratus gibelio from CyHV-2 infection. In this study, a recombinant baculovirus vector vaccine, BacCarassius-D4ORFs, containing a fused codon-optimized sequence D4ORFs comprising the ORF72 (region 1-186 nt), ORF66 (region 993-1197 nt), ORF81 (region 603-783 nt) and ORF82 (region 85-186 nt) genes of CyHV-2, driven by a Megalobrama amblycephala β-actin promoter, was constructed. Then, qPCR, Western blotting and immunofluorescence assays showed that the fused gene D4ORFs was successfully delivered and expressed in fish cells or tissues by transduction with BacCarassius-D4ORFs. The fused gene D4ORFs could not be detected by PCR in the C. auratus gibelio injected with BacCarassius-D4ORFs after 7 weeks. Specific antibody against ORF72 could be detected in the serum of vaccinated C. auratus gibelio by injection with BacCarassius-D4ORFs. Furthermore, when C. auratus gibelio were vaccinated with BacCarassius-D4ORFs via the oral or injection route, followed by challenge with CyHV-2, the relative survival rate of immunized C. auratus gibelio reached 59.3% and 80.01%, respectively. These results suggested that BacCarassius-D4ORFs has the potential to be used as a vector-based vaccine for the prevention and treatment of disease caused by CyHV-2 infection.

Probiotics maintain intestinal secretory immunoglobulin A levels in healthy formula-fed infants: a randomised, double-blind, placebo-controlled study

Formula-fed infants are more susceptible to infectious diseases because they lack the maternal immune factors transferred from breast milk, while their own immune system is still immature. As timely probiotic administration was suggested to promote immune system development in formula-fed infants, this study aimed at assessing the safety and the effects of a probiotic supplement (Bifidobacterium infantis R0033, Bifidobacterium bifidum R0071, and Lactobacillus helveticus R0052) on mucosal immune competence and digestive function in formula-fed infants. Healthy infants (3.5-6 months old) were randomised to receive either probiotic- (n=66) or placebo-supplemented (n=66) formula once a day for four weeks. In the probiotics group, faecal secretory immunoglobulin A (SIgA) levels remained similar between visit 2 (baseline; V2) and visit 3 (end-of-treatment; V3), but decreased in the placebo group. Changes in SIgA levels following treatment (log₁₀ΔV3-V2 [95%CI]) between the probiotic and placebo groups were statistically significant (23 ng/dl [-57;102] and -137 ng/dl [-212;-62], respectively (P=0.0044; ANCOVA)). While log₁₀ΔV3-V2 [95%CI] for salivary SIgA levels increased in both groups, this trend was more pronounced in the probiotics than in the placebo group with an increase of 123 ng/dl [9;236] and 37 ng/dL [-72;147], respectively (P=0.2829; ANCOVA). The weekly average number of stools/day was significantly higher in the probiotics group compared to placebo during the last week of treatment for the per protocol population. There was no difference in microbiota composition or anthropometric parameters between groups. No serious adverse event was reported, and all adverse events were mild and unrelated to the product or study. Our results show that formula-fed infants receiving probiotics maintained higher faecal SIgA levels at the end of the four-week treatment period, suggesting a positive effect of probiotics on SIgA production. This study demonstrates the safety of this probiotic formulation in infants. Formula-fed infants may benefit from probiotics supplementation to sustain the development of mucosal immunity.

Bombyx mori cypovirus encoded small peptide inhibits viral multiplication

Bombyx mori cypovirus (BmCPV) is one of the most infectious pathogen in sericulture and a member of the family Reoviridae. It specifically infects the midgut of silkworm. The BmCPV genome consists of 10 dsRNAs segments (S1-S10), which have generally been assumed to be monocistronic. In this study, a small open reading frame encoding the peptide S5-sORF, containing 27 amino acid residues, was predicted in a region of the negative (-) strand of BmCPV segment S5. An immunofluorescence assay detected S5-sORF in the cytoplasm and nuclei of BmCPV-infected cells, and it was also detected in the virion with western blotting, suggesting that S5-sORF may be assembled into the BmCPV virion. Viral gene expression was inhibited by overexpressed S5-sORF, and viral multiplication was dose-dependently suppressed by the S5-sORF peptide. A viable recombinant virus, BmCPV-S5-sORF^mut, in which the start codon (ATG) of S5-sORF was mutated to a stop codon (TGA), was generated with reverse genetics. The proliferation of BmCPV was increased by the abolition of S5-sORF expression. Furthermore, the RNA transcript of S5-sORF and small peptide of S5-sORF were involved in BmCPV replication. The expression of genes related to the innate immune pathways and apoptosis in the silkworm were not significantly affected by S5-sORF overexpression. Our results suggest that a viral nucleotide sequence is utilized by the host to generate an antiviral peptide, which may be a novel strategy protecting the host from viral infection.

CRISPR/Cas9 Initiated Transgenic Silkworms as a Natural Spinner of Spider Silk

Using transgenic silkworms with their natural spinning apparatus has proven to be a promising way to spin spider silk-like fibers. The challenges are incorporating native-size spider silk proteins and achieving an inheritable transgenic silkworm strain. In this study, a CRISPR/Cas9 initiated fixed-point strategy was used to successfully incorporate spider silk protein genes into the Bombyx mori genome. Native-size spider silk genes (up to 10 kb) were inserted into an intron of the fibroin heavy or light chain (FibH or FibL) ensuring that any sequence changes induced by the CRISPR/Cas9 would not impact protein production. The resulting fibers are as strong as native spider silks (1.2 GPa tensile strength). The transgenic silkworms have been tracked for several generations with normal inheritance of the transgenes. This strategy demonstrates the feasibility of using silkworms as a natural spider silk spinner for industrial production of high-performance fibers.

Functional characterization of BmOVOs in silkworm, Bombyx mori

Background

In our previous study, we identified four isoforms of the Bmovo gene, Bmovo-1, Bmovo-2, Bmovo-3 and Bmovo-4 from the silkworm ovary and verified that ovarian development was regulated by the BmOVO proteins. Results: To understand the regulatory mechanisms of ovarian development, the regulation of four BmOVO isoforms on the B. mori ovarian tumor (Bmotu) promoter activity was investigated with luciferase reporter assays. The results showed the Bmotu promoter activity was positively regulated by BmOVO-1, BmOVO-2, BmOVO-3 and BmOVO-4 in a dose-dependent manner, of which BmOVO-2 had the highest transcriptional activation. However, the first (A1) and third acidic domains (A3) at the N-terminus of BmOVO-1 are transcriptional repression domains, while the fourth (A4) and fifth acidic domains (A5) are transcriptional activation domains. A recombinant BmOVO zinc-finger domain was found to bind to the GTACCGTTGTA sequence located at the Bmotu promoter. Furthermore, the Bmotu promoter activity was negatively regulated by ‘Tal-like’ peptide, which can trigger BmOVO-1 degradation at the N-terminus.

Conclusions

These results will help us to further understand the regulatory mechanisms of BmOVO isoforms on Bmotu promoter activity and ovarian development in the silkworm.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5697-y) contains supplementary material, which is available to authorized users.

Cultured cells and wing disc size of silkworm can be controlled by the Hippo pathway

Hippo signalling represents a cell proliferation and organ-size control pathway. Yorki (Yki), a component of the Hippo pathway, induces the transcription of a number of targets that promote cell proliferation and survival. The functions of Yki have been characterized in Drosophila and mammals, while there are few reports on silkworm, Bombyx mori. In the present study, we found that BmYki3 facilitates cell migration and cell division, and enlarges the cultured cell and wing disc size. Co-immunoprecipitation results indicated that BmYki3 may interact with thymosin, E3 ubiquitin-protein ligase, protein kinase ASK1, dedicator of cytokinesis protein 1, calcium-independent phospholipase A2 and beta-spectrin. RNA-seq results indicated that 4444 genes were upregulated and 10 291 genes were downregulated after BmYki3 was overexpressed in the cultured cells. GO annotation indicated that the up/downregulated genes were enriched in 268/382 GO terms (p < 0.01); KEGG analysis showed that the up/downregulated genes were enriched in 49/101 pathways. These findings provided novel information to understand the functions of BmYki3 in a cell proliferation and organ-size control pathway.