The Holothuria leucospilota genome elucidates sacrificial organ expulsion and bioadhesive trap enriched with amyloid-patterned proteins.

Some tropical sea cucumbers of the family Holothuriidae can efficiently repel or even fatally ensnare predators by sacrificially ejecting a bioadhesive matrix termed the Cuvierian organ (CO), so named by the French zoologist Georges Cuvier who first described it in 1831. Still, the precise mechanisms for how adhesiveness genetically arose in CO and how sea cucumbers perceive and transduce danger signals for CO expulsion during defense have remained unclear. Here, we report the first high-quality, chromosome-level genome assembly of Holothuria leucospilota, an ecologically significant sea cucumber with prototypical CO. The H. leucospilota genome reveals characteristic long-repeat signatures in CO-specific outer-layer proteins, analogous to fibrous proteins of disparate species origins, including spider spidroin and silkworm fibroin. Intriguingly, several CO-specific proteins occur with amyloid-like patterns featuring extensive intramolecular cross-ß structures readily stainable by amyloid indicator dyes. Distinct proteins within the CO connective tissue and outer surface cooperate to give the expelled matrix its apparent tenacity and adhesiveness, respectively. Genomic evidence offers further hints that H. leucospilota directly transduces predator-induced mechanical pressure onto the CO surface through mediation by transient receptor potential channels, which culminates in acetylcholine-triggered CO expulsion in part or in entirety. Evolutionarily, innovative events in two distinct regions of the H. leucospilota genome have apparently spurred CO's differentiation from the respiratory tree to a lethal defensive organ against predators.

TransGEM: a molecule generation model based on Transformer with gene expression data.

MOTIVATION: It is difficult to generate new molecules with desirable bioactivity through ligand-based de novo drug design, and receptor-based de novo drug design is constrained by disease target information availability. The combination of artificial intelligence and phenotype-based de novo drug design can generate new bioactive molecules, independent from disease target information. Gene expression profiles can be used to characterize biological phenotypes. The Transformer model can be utilized to capture the associations between gene expression profiles and molecular structures due to its remarkable ability in processing contextual information. RESULTS: We propose TransGEM (Transformer-based model from gene expression to molecules), which is a phenotype-based de novo drug design model. A specialized gene expression encoder is used to embed gene expression difference values between diseased cell lines and their corresponding normal tissue cells into TransGEM model. The results demonstrate that the TransGEM model can generate molecules with desirable evaluation metrics and property distributions. Case studies illustrate that TransGEM model can generate structurally novel molecules with good binding affinity to disease target proteins. The majority of genes with high attention scores obtained from TransGEM model are associated with the onset of the disease, indicating the potential of these genes as disease targets. Therefore, this study provides a new paradigm for de novo drug design, and it will promote phenotype-based drug discovery. AVAILABILITY AND IMPLEMENTATION: The code is available at https://github.com/hzauzqy/TransGEM.

Whole genome sequencing and genome annotation of Dermacoccus abyssi strain HZAU 226 isolated from spoiled eggs.

Dermacoccus abyssi strain HZAU 226 is a spoilage bacterium isolated from eggs. So far, there are still few genomic resources available on the Dermacoccus abyssi. Here, we reported the complete genome sequence of Dermacoccus abyssi strain HZAU 226. High-quality DNA was extracted using the Qiagen kit, then single-molecule sequencing was performed by GridION sequencer. The raw data was quality-controlled and assembled to obtain the final genome, which consisted of a complete genome of 2,992,060 bp circular chromosome and a 64,524 bp plasmid. The structural and functional annotations of the genome were achieved through the analysis of different available databases, including antibiotic resistance genes, secondary metabolite synthesis genes and stress-related genes. Meanwhile, comparative genomic analyses of the strains were also performed. This is the first report on the complete genome of Dermacoccus abyssi, which will provide genomic resources for the study of spoilage bacteria in eggs.

Reconstructed Genome-Scale Metabolic Model Characterizes Adaptive Metabolic Flux Changes in Peripheral Blood Mononuclear Cells in Severe COVID-19 Patients.

Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) poses a mortal threat to human health. The elucidation of the relationship between peripheral immune cells and the development of inflammation is essential for revealing the pathogenic mechanism of COVID-19 and developing related antiviral drugs. The immune cell metabolism-targeting therapies exhibit a desirable anti-inflammatory effect in some treatment cases. In this study, based on differentially expressed gene (DEG) analysis, a genome-scale metabolic model (GSMM) was reconstructed by integrating transcriptome data to characterize the adaptive metabolic changes in peripheral blood mononuclear cells (PBMCs) in severe COVID-19 patients. Differential flux analysis revealed that metabolic changes such as enhanced aerobic glycolysis, impaired oxidative phosphorylation, fluctuating biogenesis of lipids, vitamins (folate and retinol), and nucleotides played important roles in the inflammation adaptation of PBMCs. Moreover, the main metabolic enzymes such as the solute carrier (SLC) family 2 member 3 (SLC2A3) and fatty acid synthase (FASN), responsible for the reactions with large differential fluxes, were identified as potential therapeutic targets. Our results revealed the inflammation regulation potentials of partial metabolic reactions with differential fluxes and their metabolites. This study provides a reference for developing potential PBMC metabolism-targeting therapy strategies against COVID-19.

Global profile of tRNA-derived small RNAs in gastric cancer patient plasma and identification of tRF-33-P4R8YP9LON4VDP as a new tumor suppressor.

Transfer RNA (tRNA)-derived small RNAs (tsRNAs) have been found to play important roles in the occurrence and development of cancers. However, the tsRNA profile in gastric cancer is unknown. In this study, we aimed to identify the global tsRNA profile in plasma from gastric cancer patients and elucidate the role of tRF-33-P4R8YP9LON4VDP in gastric cancer. Differentially expressed tsRNAs in the plasma of gastric cancer patients and healthy controls were investigated using RNA sequencing. The expression levels of tRF-33-P4R8YP9LON4VDP in the plasma of gastric cancer patients, healthy controls and gastric cancer cell lines were first detected by quantitative reverse transcription-polymerase chain reaction. The effects of tRF-33-P4R8YP9LON4VDP overexpression or downregulation in gastric cancer cells on proliferation, migration, apoptosis, and cell cycle were analyzed using the Cell Counting Kit-8, scratch assay, Transwell assay, and flow cytometry, respectively. There were 21 upregulated and 46 downregulated tsRNAs found in plasma from gastric cancer patients. The significantly upregulated tsRNAs included tRF-18-S3M83004, tRF-31-PNR8YP9LON4VD, tRF-19-3L7L73JD, tRF-33-P4R8YP9LON4VDP, tRF-31-PER8YP9LON4VD, tRF-18-MBQ4NKDJ, and tRF-31-PIR8YP9LON4VD. The significantly downregulated tsRNAs included tRF-41-YDLBRY73W0K5KKOVD, tRF-18-07QSNHD2, tRF-28-86J8WPMN1E0J, tRF-29-86V8WPMN1EJ3, tRF-31-6978WPRLXN4VE, tRF-30-MIF91SS2P46I, tRF-26-MI7O3B1NR8E, tRF-30-RRJ89O9NF5W8, tRF-26-XIP2801MK8E, and tRF-35-V0J8O9YEKPRS93, In vitro studies showed that tRF-33-P4R8YP9LON4VDP inhibited proliferation of gastric cancer cells. In conclusion, tsRNAs such as tRF-33-P4R8YP9LON4VDP could serve as a novel diagnostic biomarker and target for gastric cancer therapeutics.

Low expression of hsa_circ_0001811 in gastric cancer and its role in clinical diagnosis.

BACKGROUND: As far as gastric cancer is concerned, there is a lack of specific early diagnostic biomarkers in clinical diagnosis. Circular RNAs (circRNAs) have been found to be stable in gastric cancer tissues and plasma, so they have the potential to become new type of diagnostic biomarkers for gastric cancer. MATERIALS AND METHODS: The hsa_circ_0001811 expressions in gastric cancer tissues and paired non-cancer tissues, preoperative and postoperative plasma of patients with gastric cancer, and plasma of healthy volunteers were detected using quantitative reverse transcription-polymerase chain reaction. The receiver operating characteristic (ROC) curves were drawn; and the combined ROC curves were used to analyze their diagnostic values. The correlation between the plasma or tissue hsa_circ_0001811 levels and the clinicopathological factors of gastric cancer was further analyzed. RESULTS: Hsa_circ_0001811 was found to be lowly expressed in gastric cancer tissues and plasma from patients with gastric cancer. The area under the ROC curve (AUC) in the gastric cancer tissues and the plasma of patients with gastric cancer was 0.658 and 0.747, respectively. The AUC increased to 0.824 after combination of both. The expression level of hsa_circ_0001811 in gastric cancer tissue correlated with carcinoembryonic antigen (CEA) (P = .0347), tissue differentiation (P = .0138), and lymph node metastasis (P = .0234), while plasma hsa_circ_0001811 level was related to carbohydrate antigen (CA19-9) (P = .0278), lymph node metastasis (P = .0469), distant metastasis (P = .0384), and age (P = .0085). CONCLUSIONS: The above results indicate that hsa_circ_0001811 may become a new biomarker for clinical diagnosis of gastric cancer.

Transcriptomic analysis of sea cucumber (Holothuria leucospilota) coelomocytes revealed the echinoderm cytokine response during immune challenge.

BACKGROUND: The sea cucumber Holothuria leucospilota belongs to echinoderm, which is evolutionally the most primitive group of deuterostomes. Sea cucumber has a cavity between its digestive tract and the body wall that is filled with fluid and suspended coelomic cells similar to blood cells. The humoral immune response of the sea cucumber is based on the secretion of various immune factors from coelomocytes into the coelomic cavity. The aim of this study is to lay out a foundation for the immune mechanisms in echinoderms and their origins in chordates by using RNA-seq. RESULTS: Sea cucumber primary coelomocytes were isolated from healthy H. leucospilota and incubated with lipopolysaccharide (LPS, 10 µg/ml), polyinosinic-polycytidylic acid [Poly (I:C), 10 µg/ml] and heat-inactived Vibrio harveyi (107 cell/ml) for 24 h, respectively. After high-throughput mRNA sequencing on an Illumina HiSeq2500, a de novo transcriptome was assembled and the Unigenes were annotated. Thirteen differentially expressed genes (DEGs) were selected randomly from our data and subsequently verified by using RT-qPCR. The results of RT-qPCR were consistent with those of the RNA-seq (R2 = 0.61). The top 10 significantly enriched signaling pathways and immune-related pathways of the common and unique DEGs were screened from the transcriptome data. Twenty-one cytokine candidate DEGs were identified, which belong to 4 cytokine families, namely, BCL/CLL, EPRF1, IL-17 and TSP/TPO. Gene expression in response to LPS dose-increased treatment (0, 10, 20 and 50 µg/ml) showed that IL-17 family cytokines were significantly upregulated after 10 µg/ml LPS challenge for 24 h. CONCLUSION: A de novo transcriptome was sequenced and assembled to generate the gene expression profiling across the sea cucumber coelomocytes treated with LPS, Poly (I:C) and V. harveyi. The cytokine genes identified in DEGs could be classified into 4 cytokine families, in which the expression of IL-17 family cytokines was most significantly induced after 10 µg/ml LPS challenge for 24 h. Our findings have laid the foundation not only for the research of molecular mechanisms related to the immune response in echinoderms but also for their origins in chordates, particularly in higher vertebrates.

Diastereoselective synthesis and conformational analysis of 4,5-difluoropipecolic acids.

Stereoselectively-fluorinated analogs of pipecolic acid have been investigated through a combined theoretical and experimental approach. Three of the four possible diastereoisomers of 4,5-difluoropipecolic acid were successfully synthesized via deoxyfluorination chemistry, navigating a complex reaction network that included neighboring group participation, rearrangement, and elimination pathways. A DFT-based conformational study, supported by NMR J-based analysis, revealed that the different diastereoisomers of 4,5-difluoropipecolic acid preferentially adopt different puckers of the six-membered ring. These findings could have future relevance for the conformational control of biologically active peptides.

The biogenesis and biological functions of circular RNAs and their molecular diagnostic values in cancers.

BACKGROUND: In addition to non-coding RNAs (lncRNAs) and microRNAs (miRNAs), circular RNAs (circRNAs) are endogenous RNAs with various functions, which have recently become a research hotspot. CircRNAs are a kind of closed circular RNA molecule widely existing in transcriptomes. Due to lack of free ends, they are not easily cleaved by RNase R, thus avoiding degradation. They are more stable than linear RNAs. METHODS: Data were collected through PubMed. The following search terms were used: "circular RNA," "circRNA," "cancer," "mechanism," "biogenesis," "biomarker," "diagnosis." Only articles published in English were included. RESULTS: Most circRNAs express tissue/developmental stage specificity. Moreover, circRNAs are involved in the regulation of a variety of biological activities. In this review, we discuss the formation, classification, and biological functions of circRNAs, especially their molecular diagnostic values in common cancers, including gastric cancer (hsa_circ_002059, circ_LARP4, hsa_circ_0000190, hsa_circ_0000096, circ-SFMBT2, and circ_PVT1), hepatocellular carcinoma (circ_104075, circRNA_100338, circ_MTO1, and circZKSCAN1), colorectal cancer (hsa_circ_0136666 and hsa_circ_0000523), lung cancer (hsa_circ_0006427, circ_100876, and circ_ABCB10), breast cancer (hsa_circ_0089105, circAGFG1, and circEPSTI1), bladder cancer (circFNDC3B and circTFRC), and esophageal squamous cell carcinoma (circ_100876 and circ-DLG1). CONCLUSION: CircRNAs not only play important roles in tumorigenesis, but also may become new diagnostic biomarkers.

An analysis of the mechanism underlying photocatalytic disinfection based on integrated metabolic networks and transcriptional data.

Photocatalytic disinfection has long been used to combat pathogenic bacteria. However, the specific mechanism underlying photocatalytic disinfection and its corresponding targets remain unclear. In this study, an analysis of the potential mechanism underlying photocatalytic disinfection was performed based on integrated metabolic networks and transcriptional data. Two sets of RNA-seq data (wild type and a photocatalysis-resistant mutant mediated by titanium dioxide (TiO2)) were processed to constrain the genome scale metabolic models (GSMM) of E. coli. By analyzing the metabolic network, the differential metabolic flux of every reaction was computed in constrained GSMM, and several significantly differential metabolic fluxes in reactions were extracted and analyzed. Most of these reactions were involved in the transmembrane transport of substances and occurred on the inner membrane or were an important component of the cell membrane. These results, which are consistent with the reported information, validated our analysis process. In addition, our work also identified other new and valuable metabolic pathways, such as the reaction ALCD2x, which has a great effect on the energy production process under bacterial anaerobic conditions. The DHAK reaction is also related to the metabolic process of ATP. These reactions with large differential metabolic fluxes merit further research. Additionally, to provide a strategy to address photocatalysis-resistant mutant bacteria, a metabolic compensation analysis was also performed. The metabolic compensation analysis results provided suggestions for a combined method that can effectively combat resistant bacteria. This method could also be used to explore the mechanisms of drug resistance in other microorganisms.

Tumor-suppressive circular RNAs: Mechanisms underlying their suppression of tumor occurrence and use as therapeutic targets.

Circular RNAs (circRNAs) have a covalently closed circular conformation and are structurally stable. Those circRNAs with tumor-suppressive properties play an important role in tumorigenesis and metastasis and thus may be used as therapeutic targets of cancers. Herein, we review the current understanding of the classification of circRNAs and summarize the functions and mechanisms of circRNAs that have tumor-suppressive roles in various cancers, including liver cancer (circARSP91, circADAMTS13, circADAMTS14, circMTO1, hsa_circ_0079299, and circC3P1), bladder cancer (circFNDC3B, circITCH, circHIPK3, circRNA-3, cdrlas, and circLPAR1), gastric cancer (circLARP4, circYAP1, hsa_cric_0000096, hsa_circ_0000993, and circPSMC3), breast cancer (circ_000911, hsa_circ_0072309, and circASS1), lung cancer (hsa_circ_0000977, circPTK2, circ_0001649, hsa_circ_100395, and circ_0006916), glioma (circ_0001946, circSHPRH, and circFBXW7), and colorectal cancer (circITGA7 and hsa_circ_0014717). Thanks to their structural stability, these tumor-suppressive circRNAs may be used as potential and potent therapeutic targets. Moreover, we propose a new method for the classification of circRNAs. Based on whether they can be translated, circRNAs can be divided into noncoding circRNAs and coding circRNAs.

Extracellular production of recombinant N-glycosylated anti-VEGFR2 monobody in leaky Escherichia coli strain.

OBJECTIVE: To improve the production yield of N-glycosylated anti-VEGFR2 (vascular endothelial growth factor receptor 2) monobody (FN3VEGFR2-Gly) in lpp knockout Escherichia coli cells harboring Campylobacter jejuni N-glycosylation pathway. RESULTS: The leaky CLM37-Δlpp strain efficiently secreted FN3VEGFR2-Gly into culture medium. The extracellular levels of glycosylated FN3VEGFR2-Gly in CLM37-Δlpp culture medium were approximately 11 and 15 times higher compared to those in CLM37 cells via IPTG and auto-induction, respectively. In addition, the highest level of total glycosylated FN3VEGFR2-Gly (70 ± 3.4 mg/L) was found in culture medium via auto-induction. Furthermore, glycosylated FN3VEGFR2-Gly was more stable than unglycosylated FN3VEGFR2-Gly in this expression system, but their bioactivities were relatively similar. CONCLUSIONS: Lpp knockout leaky E. coli strain combined with auto-induction method can enhance the extracellular production of homogenous N-glycosylated FN3VEGFR2-Gly, and facilitate the downstream protein purification. The findings of this study may provide practical implications for the large-scale production and cost-effective harvesting of N-glycosylation proteins.

Increased glycosylation efficiency of recombinant proteins in Escherichia coli by auto-induction.

Escherichia coli cells have been considered as promising hosts for producing N-glycosylated proteins since the successful production of N-glycosylated protein in E. coli with the pgl (N-linked protein glycosylation) locus from Campylobacter jejuni. However, one hurdle in producing N-glycosylated proteins in large scale using E. coli is inefficient glycan glycosylation. In this study, we developed a strategy for the production of N-glycosylated proteins with high efficiency via an optimized auto-induction method. The 10th human fibronectin type III domain (FN3) was engineered with native glycosylation sequon DFNRSK and optimized DQNAT sequon in C-terminus with flexible linker as acceptor protein models. The resulting glycosylation efficiencies were confirmed by Western blots with anti-FLAG M1 antibody. Increased efficiency of glycosylation was obtained by changing the conventional IPTG induction to auto-induction method, which increased the glycosylation efficiencies from 60% and 75% up to 90% and 100% respectively. Moreover, in the condition of inserting the glycosylation sequon in the loop of FN3 (the acceptor sequon with local structural conformation), the glycosylation efficiency was increased from 35% to 80% by our optimized auto-induction procedures. To justify the potential for general application of the optimized auto-induction method, the reconstituted lsg locus from Haemophilus influenzae and PglB from C. jejuni were utilized, and this led to 100% glycosylation efficiency. Our studies provided quantitative evidence that the optimized auto-induction method will facilitate the large-scale production of pure exogenous N-glycosylation proteins in E. coli cells.

Malaria in China, 2011-2015: an observational study.

OBJECTIVE: To ascertain the trends and burden of malaria in China and the costs of interventions for 2011-2015. METHODS: We analysed the spatiotemporal and demographic features of locally transmitted and imported malaria cases using disaggregated surveillance data on malaria from 2011 to 2015, covering the range of dominant malaria vectors in China. The total and mean costs for malaria elimination were calculated by funding sources, interventions and population at risk. FINDINGS: A total of 17 745 malaria cases, including 123 deaths (0.7%), were reported in mainland China, with 15 840 (89%) being imported cases, mainly from Africa and south-east Asia. Almost all counties of China (2855/2858) had achieved their elimination goals by 2015, and locally transmitted cases dropped from 1469 cases in 2011 to 43 cases in 2015, mainly occurring in the regions bordering Myanmar where Anopheles minimus and An. dirus are the dominant vector species. A total of United States dollars (US$) 134.6 million was spent in efforts to eliminate malaria during 2011-2015, with US$ 57.2 million (43%) from the Global Fund to Fight AIDS, Tuberculosis and Malaria and US$ 77.3 million (57%) from the Chinese central government. The mean annual investment (US$ 27 million) per person at risk (574 million) was US$ 0.05 (standard deviation: 0.03). CONCLUSION: The locally transmitted malaria burden in China has decreased. The key challenge is to address the remaining local transmission, as well as to reduce imported cases from Africa and south-east Asia. Continued efforts and appropriate levels of investment are needed in the 2016-2020 period to achieve elimination.

[Analysis of schistosomiasis prevalence in endemic areas using retrospective space-time permutation scan statistics].

Objective: To analyze the characteristics of schistosomasis prevalence by using the spatial epidemiological method, and test the application of retrospective space-time permutation scan statistics in determining mountainous and lake-type endemic areas of schistosomiasis. Methods: The data of schistosomasis in humans, cattle and snails in Jiangxi Province during 2009-2014 and in Yunnan Province during 2004-2013 were collected and analyzed. The temporal and spatial distribution of schistosomiasis endemic areas in the two provinces was analyzed with retrospective space-time permutation scan statistics. Results: The prevalence of schistosomiasis in residents and Oncomelania snails showed a trend of decline in Jiangxi, from 0.21% and 0.03% in 2009 to 0.01% and zero in 2014. A similar trend was found in cattle, from 1.25% in 2012 to 0.12% in 2014. The average annual percentage change (APC) in residents was-47.36%(P < 0.05). The space-time permutation clustering analysis revealed a temporal and spatial clustering of schistosomiasis prevalence from 2009 to 2014 in residents, cattle, and snails, with 3,2 and 1 clustering areas, respectively, all distributed in Poyang Lake Region. A similar declining trend of schistosomiasis prevalence was found in residents, snails and cattle in Yunnan during 2004-2013, from 2.49%,0.70% and 3.76% in 2004 to no infection in residents and snails and 0.02% in cattle in 2013. The APC in residents was-49.17%(P < 0.05). There was a temporal and spatial clustering of schistosomiasis prevalence during 2004-2013 in residents, cattle, and snails, with 2,2 and 6 clustering areas, respectively. Conclusion: A declining trend of schistosomiasis prevalence is shown in lake-type endemic areas in Jiangxi during 2009-2014 and in mountainous endemic areas in Yunnan during 2004-2013. The retrospective space-time permutation scan statistics reveal a clustering of schistosomiasis in humans, cattle, and snails, suggesting its applicability in analyzing the temporal and spatial distribution of schistosomiasis.

[Design and implementation of command system for emergency events of parasitic diseases].

Based on the requirement analysis and functional design of the command system for parasitic disease outbreaks, the system was constructed by workflow technique, function modules and technical architecture. The command system was a multi-platform system, could achieve multiple functions, such as monitoring and early warning of parasitic diseases, emergency video communication, emergency dispatcher, and emergency management. The system can meet the needs in emergency events of parasitic diseases, and increase preparedness level.

PCMD: A multilevel comparison database of intra- and cross-species metabolic profiling in 530 plant species.

Comparative metabolomics plays a crucial role in investigating gene function, exploring metabolite evolution, and accelerating crop genetic improvement. However, a systematic platform for comparing intra- and cross-species metabolites is currently lacking. Here, we report the Plant Comparative Metabolome Database (PCMD; http://yanglab.hzau.edu.cn/PCMD), a multilevel comparison database based on predicted metabolic profiles in 530 plant species. PCMD serves as a platform for comparing metabolite characteristics at various levels, including species, metabolites, pathways, and biological taxonomy. The database also provides a series of user-friendly online tools, such as Species-comparison, Metabolites-enrichment, and ID conversion, enabling users to perform comparisons and enrichment analyses of metabolites across different species. In addition, PCMD establishes a unified system based on existing metabolite-related databases by standardizing metabolite numbering. PCMD is the most species-rich comparative plant metabolomics database currently available, and a case study demonstrated its capability to provide new insights into understanding plant metabolic diversity.

Metabolites of pathogenic microorganisms database (MPMdb) and its seed metabolite applications.

Seed metabolites are the combination of essential compounds required by an organism across various potential environmental conditions. The seed metabolites screening framework based on the network topology approach can capture important biological information of species. This study aims to identify comprehensively the relationship between seed metabolites and pathogenic bacteria. A large-scale data set was compiled, describing the seed metabolite sets and metabolite sets of 124,192 pathogenic strains from 34 genera, by constructing genome-scale metabolic models. The enrichment analysis method was used to screen the specific seed metabolites of each species/genus of pathogenic bacteria. The metabolites of pathogenic microorganisms database (MPMdb) (http://qyzhanglab.hzau.edu.cn/MPMdb/) was established for browsing, searching, predicting, or downloading metabolites and seed metabolites of pathogenic microorganisms. Based on the MPMdb, taxonomic and phylogenetic analyses of pathogenic bacteria were performed according to the function of seed metabolites and metabolites. The results showed that the seed metabolites could be used as a feature for microorganism chemotaxonomy, and they could mirror the phylogeny of pathogenic bacteria. In addition, our screened specific seed metabolites of pathogenic bacteria can be used not only for further tapping the nutritional resources and identifying auxotrophies of pathogenic bacteria but also for designing targeted bactericidal compounds by combining with existing antimicrobial agents.IMPORTANCEMetabolites serve as key communication links between pathogenic microorganisms and hosts, with seed metabolites being crucial for microbial growth, reproduction, external communication, and host infection. However, the large-scale screening of metabolites and the identification of seed metabolites have always been the main technical bottleneck due to the low throughput and costly analysis. Genome-scale metabolic models have become a recognized research paradigm to investigate the metabolic characteristics of species. The developed metabolites of pathogenic microorganisms database in this study is committed to systematically predicting and identifying the metabolites and seed metabolites of pathogenic microorganisms, which could provide a powerful resource platform for pathogenic bacteria research.

Transporter proteins knowledge graph construction and its application in drug development.

Transporters are the main determinant for pharmacokinetics characteristics of drugs, such as absorption, distribution, and excretion of drugs in humans. However, it is difficult to perform drug transporter validation and structure analysis of membrane transporter proteins by experimental methods. Many studies have demonstrated that knowledge graphs (KG) could effectively excavate potential association information between different entities. To improve the effectiveness of drug discovery, a transporter-related KG was constructed in this study. Meanwhile, a predictive frame (AutoInt_KG) and a generative frame (MolGPT_KG) were established based on the heterogeneity information obtained from the transporter-related KG by the RESCAL model. Natural product Luteolin with known transporters was selected to verify the reliability of the AutoInt_KG frame, its ROC-AUC (1:1), ROC-AUC (1:10), PR-AUC (1:1), PR-AUC (1:10) are 0.91, 0.94, 0.91 and 0.78, respectively. Subsequently, the MolGPT_KG frame was constructed to implement efficient drug design based on transporter structure. The evaluation results showed that the MolGPT_KG could generate novel and valid molecules and that these molecules were further confirmed by molecular docking analysis. The docking results showed that they could bind to important amino acids at the active site of the target transporter. Our findings will provide rich information resources and guidance for the further development of the transporter-related drugs.

Pacific white shrimp (Litopenaeus vannamei) vitelline membrane outer layer protein 1 (VMO1) is produced in the hepatopancreas and transported into ovarian oocytes during vitellogenesis.

In birds, vitelline membrane outer layer protein 1 (VMO1) is an exogenous protein that can be absorbed by eggs as a barrier to prevent the mixing of yolk and egg white. However, researches on VMO1 are limited in birds but not other non-avian species until now. In this study, we first identified a novel Vmo1 cDNA (Lv-Vmo1) in Pacific white shrimp (Litopenaeus vannamei), the most important cultured shrimp in the world. We further analyzed its gene organization, phylogenetic relationship and protein structure. The Lv-Vmo1 transcript was specifically expressed in the hepatopancreas without sexual dimorphism. During ovarian development in female, the hepatopancreatic Lv-Vmo1 mRNA levels showed a significant increase. By in situ hybridization, Lv-Vmo1 mRNA was present in three cell types of the hepatopancreas but neither oocytes nor follicle cells of the ovary. In contrast, immunofluorescence revealed that Lv-VMO1 protein was distributed in the cytoplasms of both hepatopancreatic cells and ovarian oocytes. Western blot showed that Lv-VMO1 protein was produced in the hepatopancreas and transported to the ovary via hemolymph circulation. Identification of a species-specific egg-entry guide protein is the key to the receptor-mediated ovarian transduction of cargo, a novel gene editing approach in oviparous animals. This study lays the mechanism for exogenous transport into penaeid shrimp eggs by VMO1, as a foundation for achieving exogenous protein-mediated incorporation into oocytes.