The UniProt-GO Annotation database in 2011

Cucurbitacin B targets STAT3 to induce ferroptosis in non-small cell lung cancer

Cucurbitacin B (CuB) is a compound found in plants like Cucurbitaceae that has shown promise in fighting cancer, particularly in lung cancer. However, the specific impact of CuB on ferroptosis and how it works in lung cancer cells has not been fully understood. Our research has discovered that CuB can effectively slow down the growth of non-small cell lung cancer (NSCLC) cells. Even in small amounts, it was able to inhibit the growth of various NSCLC cell lines. This inhibitory effect was reversed when ferroptosis inhibitors DFO, Lip-1 and Fer-1 were introduced. CuB was found to increase the levels of reactive oxygen species (ROS), lipid ROS, MDA, and ferrous ions within H358 lung cancer cells, leading to a decrease in GSH, mitochondrial membrane potential (MMP) and changes in ferroptosis-related proteins in a dose-dependent manner. These findings were also confirmed in A549 lung cancer cells. In A549 cells, different concentrations of CuB induced the accumulation of intracellular lipid ROS, ferrous ions and changes in ferroptosis-related indicators in a concentration-dependent manner. Meanwhile, the cytotoxic effect induced by CuB in A549 cells was counteracted by ferroptosis inhibitors DFO and Fer-1. Through network pharmacology, we identified potential targets related to ferroptosis in NSCLC cells treated with CuB, with STAT3 targets showing high scores. Further experiments using molecular docking and cell thermal shift assay (CETSA) revealed that CuB interacts with the STAT3 protein. Western blot and immunofluorescence staining demonstrated that CuB inhibits the phosphorylation of STAT3 (P-STAT3) in H358 cells. Silencing STAT3 enhanced CuB-induced accumulation of lipid ROS and iron ions, as well as the expression of ferroptosis-related proteins. On the other hand, overexpression of STAT3 reversed the effects of CuB-induced ferroptosis. The results indicate that CuB has the capability to suppress STAT3 activation, resulting in ferroptosis, and could be a promising treatment choice for NSCLC.

Metatranscriptomic response of deep ocean microbial populations to infusions of oil and/or synthetic chemical dispersant

Oil spills are a frequent perturbation to the marine environment that has rapid and significant impacts on the local microbiome. Previous studies have shown that exposure to synthetic dispersant alone did not enhance heterotrophic microbial activity or oxidation rates of specific hydrocarbon components but increased the abundance of some taxa (e.g., Colwellia). In contrast, exposure to oil, but not dispersants, increased the abundance of other taxa (e.g., Marinobacter) and stimulated hydrocarbon oxidation rates. Here, we advance these findings by interpreting metatranscriptomic data from this experiment to explore how and why specific components of the microbial community responded to distinct organic carbon exposure regimes. Dispersant alone was selected for a unique community and for dominant organisms that reflected treatment- and time-dependent responses. Dispersant amendment also led to diverging functional profiles among the different treatments. Similarly, oil alone was selected for a community that was distinct from treatments amended with dispersants. The presence of oil and dispersants with added nutrients led to substantial differences in microbial responses, likely suggesting increased fitness driven by the presence of additional inorganic nutrients. The oil-only additions led to a marked increase in the expression of phages, prophages, transposable elements, and plasmids (PPTEPs), suggesting that aspects of microbial community response to oil are driven by the "mobilome," potentially through viral-associated regulation of metabolic pathways in ciliates and flagellates that would otherwise throttle the microbial community through grazing.IMPORTANCEMicrocosm experiments simulated the April 2010 Deepwater Horizon oil spill by applying oil and synthetic dispersants (Corexit EC9500A and EC9527A) to deep ocean water samples. The exposure regime revealed severe negative alterations in the treatments' heterotrophic microbial activity and hydrocarbon oxidation rates. We expanded these findings by exploring metatranscriptomic signatures of the microbial communities during the chemical amendments in the microcosm experiments. Here we report how dominant organisms were uniquely associated with treatment- and time-dependent trajectories during the exposure regimes; nutrient availability was a significant factor in driving changes in metatranscriptomic responses. Remarkable signals associated with PPTEPs showed the potential role of mobilome and viral-associated survival responses. These insights underscore the time-dependent environmental perturbations of fragile marine environments under oil and anthropogenic stress.

Exploring the potential mechanisms of Danshen against COVID-19 via network pharmacology analysis and molecular docking

Danshen, a prominent herb in traditional Chinese medicine (TCM), is known for its potential to enhance physiological functions such as blood circulation, immune response, and resolve blood stasis. Despite the effectiveness of COVID-19 vaccination efforts, some individuals still face severe complications post-infection, including pulmonary fibrosis, myocarditis arrhythmias and stroke. This study employs a network pharmacology and molecular docking approach to investigate the potential mechanisms underlying the therapeutic effects of candidate components and targets from Danshen in the treatment of complications in COVID-19. Candidate components and targets from Danshen were extracted from the TCMSP Database, while COVID-19-related targets were obtained from Genecards. Venn diagram analysis identified common targets. A Protein-Protein interaction (PPI) network and gene enrichment analysis elucidated potential therapeutic mechanisms. Molecular docking evaluated interactions between core targets and candidate components, followed by molecular dynamics simulations to assess stability. We identified 59 potential candidate components and 123 targets in Danshen for COVID-19 treatment. PPI analysis revealed 12 core targets, and gene enrichment analysis highlighted modulated pathways. Molecular docking showed favorable interactions, with molecular dynamics simulations indicating high stability of key complexes. Receiver operating characteristic (ROC) curves validated the docking protocol. Our study unveils candidate compounds, core targets, and molecular mechanisms of Danshen in COVID-19 treatment. These findings provide a scientific foundation for further research and potential development of therapeutic drugs.

Phosphorylation of 399S at CsHsp70 of Cymbidium sinense is essential to maintain chlorophyll stability

The Chinese orchids symbolise nobility and gentility in China, and the variation of leaf color makes Cymbidium sinense more diversified and valuable. However, its color variations especially at the protein level still remain largely unexplored. In this study, the proteomics and phosphoproteomics of Cymbidium sinense leaf color variation mutants were studied. A total of 1059 differentially abundant proteins (DAPs) and 1127 differentially abundant phosphorylation sites belonging to 644 phosphoproteins (DAPPs) were identified in the yellow section of leaf variegation mutant of Cymbidium sinense (MY) compared with the green section (MG). Moreover, 349 co-expressing proteins were found in both omics' datasets, while only 26 proteins showed the same expression patterns in the two omics. The interaction network analysis of kinases and phosphatases showed that DAPs and DAPPs in photosynthesis, response to hormones, pigment metabolic process, phosphorylation, glucose metabolic process, and dephosphorylation might contribute to leaf color variation. The abundance of 28 Hsps and 28 phosphorylation sites belonging to 10 Hsps showed significant differences between MG and MY. CsHsp70 was selected to explore the function in Cymbidium sinense leaf variegation. The results showed CsHsp70 is essential for maintaining photosynthetic pigment content and the 399S phosphorylation site is crucial to the function of CsHsp70. Collectively, our findings construct a comprehensive coverage of protein and protein phosphorylation in leaf variegation of C. sinense, providing valuable insights into its formation mechanisms.

Genomic and phenotypic signatures provide insights into the wide adaptation of a global plant invader

Invasive alien species are primary drivers of biodiversity loss and species extinction. Smooth cordgrass (Spartina alterniflora) is one of the most aggressive invasive plants in coastal ecosystems around the world. However, the genomic bases and evolutionary mechanisms underlying its invasion success have remained largely unknown. Here, we assembled a chromosome-level reference genome and performed phenotypic and population genomic analyses between native US and introduced Chinese populations. Our phenotypic comparisons showed that introduced Chinese populations have evolved competitive traits, such as early flowering time and greater plant biomass, during secondary introductions along China's coast. Population genomic and transcriptomic inferences revealed distinct evolutionary trajectories of low- and high-latitude Chinese populations. In particular, genetic mixture among different source populations, together with independent natural selection acting on distinct target genes, may have resulted in high genome dynamics of the introduced Chinese populations. Our study provides novel phenotypic and genomic evidence showing how smooth cordgrass rapidly adapts to variable environmental conditions in its introduced ranges. Moreover, candidate genes related to flowering time, fast growth, and stress tolerance (i.e., salinity and submergence) provide valuable genetic resources for future improvement of cereal crops.

Selection of Catechin Biosynthesis-Related Genes and Functional Analysis from Chromosome-Level Genome Assembly in C. sinensis L. Variety 'Sangmok'

Camellia is an important plant genus that includes well-known species such as C. sinensis, C. oleifera, and C. japonica. The C. sinensis cultivar 'Sangmok', one of Korea's standard types of tea landraces, is a small evergreen tree or shrub. Genome annotation has shown that Korean tea plants have special and unique benefits and superior components, such as catechin. The genome of Camellia sinensis cultivar 'Sangmok' was assembled on the chromosome level, with a length of 2678.62 Mbp and GC content of 38.16%. Further, 15 chromosome-scale scaffolds comprising 82.43% of the assembly (BUSCO completeness, 94.3%) were identified. Analysis of 68,151 protein-coding genes showed an average of 5.003 exons per gene. Among 82,481 coding sequences, the majority (99.06%) were annotated by Uniprot/Swiss-Prot. Further analysis revealed that 'Sangmok' is closely related to C. sinensis, with a divergence time of 60 million years ago. A total of 3336 exclusive gene families in 'Sangmok' were revealed by gene ontology analysis to play roles in auxin transport and cellular response mechanisms. By comparing these exclusive genes with 551 similar catechin genes, 17 'Sangmok'-specific catechin genes were identified by qRT-PCR, including those involved in phytoalexin biosynthesis and related to cytochrome P450. The 'Sangmok' genome exhibited distinctive genes compared to those of related species. This comprehensive genomic investigation enhances our understanding of the genetic architecture of 'Sangmok' and its specialized functions. The findings contribute valuable insights into the evolutionary and functional aspects of this plant species.

Aflibercept Off-Target Effects in Diabetic Macular Edema: An In Silico Modeling Approach

Intravitreal aflibercept injection (IAI) is a treatment for diabetic macular edema (DME), but its mechanism of action (MoA) has not been completely elucidated. Here, we aimed to explore IAI's MoA and its multi-target nature in DME pathophysiology with an in silico (computer simulation) disease model. We used the Therapeutic Performance Mapping System (Anaxomics Biotech property) to generate mathematical models based on the available scientific knowledge at the time of the study, describing the relationship between the modulation of vascular endothelial growth factor receptors (VEGFRs) by IAI and DME pathophysiological processes. We also undertook an enrichment analysis to explore the processes modulated by IAI, visualized the effectors' predicted protein activity, and specifically evaluated the role of VEGFR1 pathway inhibition on DME treatment. The models simulated the potential pathophysiology of DME and the likely IAI's MoA by inhibiting VEGFR1 and VEGFR2 signaling. The action of IAI through both signaling pathways modulated the identified pathophysiological processes associated with DME, with the strongest effects in angiogenesis, blood-retinal barrier alteration and permeability, and inflammation. VEGFR1 inhibition was essential to modulate inflammatory protein effectors. Given the role of VEGFR1 signaling on the modulation of inflammatory-related pathways, IAI may offer therapeutic advantages for DME through sustained VEGFR1 pathway inhibition.

A haplotype-resolved genome provides insight into allele-specific expression in wild walnut (Juglans regia L.)

Wild germplasm resources are crucial for gene mining and molecular breeding because of their special trait performance. Haplotype-resolved genome is an ideal solution for fully understanding the biology of subgenomes in highly heterozygous species. Here, we surveyed the genome of a wild walnut tree from Gongliu County, Xinjiang, China, and generated a haplotype-resolved reference genome of 562.99 Mb (contig N50 = 34.10 Mb) for one haplotype (hap1) and 561.07 Mb (contig N50 = 33.91 Mb) for another haplotype (hap2) using PacBio high-fidelity (HiFi) reads and Hi-C technology. Approximately 527.20 Mb (93.64%) of hap1 and 526.40 Mb (93.82%) of hap2 were assigned to 16 pseudochromosomes. A total of 41039 and 39744 protein-coding gene models were predicted for hap1 and hap2, respectively. Moreover, 123 structural variations (SVs) were identified between the two haplotype genomes. Allele-specific expression genes (ASEGs) that respond to cold stress were ultimately identified. These datasets can be used to study subgenome evolution, for functional elite gene mining and to discover the transcriptional basis of specific traits related to environmental adaptation in wild walnut.

PlantNh-Kcr: a deep learning model for predicting non-histone crotonylation sites in plants

BACKGROUND

Lysine crotonylation (Kcr) is a crucial protein post-translational modification found in histone and non-histone proteins. It plays a pivotal role in regulating diverse biological processes in both animals and plants, including gene transcription and replication, cell metabolism and differentiation, as well as photosynthesis. Despite the significance of Kcr, detection of Kcr sites through biological experiments is often time-consuming, expensive, and only a fraction of crotonylated peptides can be identified. This reality highlights the need for efficient and rapid prediction of Kcr sites through computational methods. Currently, several machine learning models exist for predicting Kcr sites in humans, yet models tailored for plants are rare. Furthermore, no downloadable Kcr site predictors or datasets have been developed specifically for plants. To address this gap, it is imperative to integrate existing Kcr sites detected in plant experiments and establish a dedicated computational model for plants.

RESULTS

Most plant Kcr sites are located on non-histones. In this study, we collected non-histone Kcr sites from five plants, including wheat, tabacum, rice, peanut, and papaya. We then conducted a comprehensive analysis of the amino acid distribution surrounding these sites. To develop a predictive model for plant non-histone Kcr sites, we combined a convolutional neural network (CNN), a bidirectional long short-term memory network (BiLSTM), and attention mechanism to build a deep learning model called PlantNh-Kcr. On both five-fold cross-validation and independent tests, PlantNh-Kcr outperformed multiple conventional machine learning models and other deep learning models. Furthermore, we conducted an analysis of species-specific effect on the PlantNh-Kcr model and found that a general model trained using data from multiple species outperforms species-specific models.

CONCLUSION

PlantNh-Kcr represents a valuable tool for predicting plant non-histone Kcr sites. We expect that this model will aid in addressing key challenges and tasks in the study of plant crotonylation sites.

A chromosome-level genome assembly of Agave hybrid NO.11648 provides insights into the CAM photosynthesis

The subfamily Agavoideae comprises crassulacean acid metabolism (CAM), C3, and C4 plants with a young age of speciation and slower mutation accumulation, making it a model crop for studying CAM evolution. However, the genetic mechanism underlying CAM evolution remains unclear because of lacking genomic information. This study assembled the genome of Agave hybrid NO.11648, a constitutive CAM plant belonging to subfamily Agavoideae, at the chromosome level using data generated from high-throughput chromosome conformation capture, Nanopore, and Illumina techniques, resulting in 30 pseudo-chromosomes with a size of 4.87 Gb and scaffold N50 of 186.42 Mb. The genome annotation revealed 58 841 protein-coding genes and 76.91% repetitive sequences, with the dominant repetitive sequences being the I-type repeats (Copia and Gypsy accounting for 18.34% and 13.5% of the genome, respectively). Our findings also provide support for a whole genome duplication event in the lineage leading to A. hybrid, which occurred after its divergence from subfamily Asparagoideae. Moreover, we identified a gene duplication event in the phosphoenolpyruvate carboxylase kinase (PEPCK) gene family and revealed that three PEPCK genes (PEPCK3, PEPCK5, and PEPCK12) were involved in the CAM pathway. More importantly, we identified transcription factors enriched in the circadian rhythm, MAPK signaling, and plant hormone signal pathway that regulate the PEPCK3 expression by analysing the transcriptome and using yeast one-hybrid assays. Our results shed light on CAM evolution and offer an essential resource for the molecular breeding program of Agave spp.

Agronomic, physiological and transcriptional characteristics provide insights into fatty acid biosynthesis in yellowhorn (Xanthoceras sorbifolium Bunge) during fruit ripening

Yellowhorn (Xanthoceras sorbifolium Bunge) is an oil-bearing tree species in northern China. In this study, we used yellowhorn from Heilongjiang to analyze the morphological and physiological changes of fruit development and conducted transcriptome sequencing. The results showed that the fruit experienced relatively slow growth from fertilization to DAF20 (20 days after flowering). From DAF40 to DAF60, the fruit entered an accelerated development stage, with a rapid increase in both transverse and longitudinal diameters, and the kernel contour developed completely at DAF40. From DAF60 to DAF80, the transverse and vertical diameters of the fruit developed slowly, and the overall measures remained stable until maturity. The soluble sugar, starch, and anthocyanin content gradually accumulated until reaching a peak at DAF80 and then rapidly decreased. RNA-seq analysis revealed differentially expressed genes (DEGs) in the seed coat and kernel, implying that seed components have different metabolite accumulation mechanisms. During the stages of seed kernel development, k-means clustering separated the DEGs into eight sub-classes, indicating gene expression shifts during the fruit ripening process. In subclass 8, the fatty acid biosynthesis pathway was enriched, suggesting that this class was responsible for lipid accumulation in the kernel. WGCNA revealed ten tissue-specific modules for the 12 samples among 20 modules. We identified 54 fatty acid biosynthesis pathway genes across the genome, of which 14 was quantified and confirmed by RT-qPCR. Most genes in the plastid synthesis stage showed high expression during the DAF40-DAF60 period, while genes in the endoplasmic reticulum synthesis stage showed diverse expression patterns. EVM0012847 (KCS) and EVM0002968 (HCD) showed similar high expression in the early stages and low expression in the late stages. EVM0022385 (HCD) exhibited decreased expression from DAF40 to DAF60 and then increased from DAF60 to DAF100. EVM0000575 (KCS) was increasingly expressed from DAF40 to DAF60 and then decreased from DAF60 to DAF100. Finally, we identified transcription factors (TFs) (HB-other, bHLH and ARF) that were predicted to bind to fatty acid biosynthesis pathway genes with significant correlations. These results are conducive to promoting the transcriptional regulation of lipid metabolism and the genetic improvement in terms of high lipid content of yellowhorn.

Acetyltransferase P300 Regulates Glucose Metabolic Reprogramming through Catalyzing Succinylation in Lung Cancer

Aberrant protein post-translational modification is a hallmark of malignant tumors. Lysine succinylation (Ksucc) plays a vital role in cell energy metabolism in various cancers. However, whether succinylation can be catalyzed by acetyltransferase p300 remains unclear. In this study, we unveiled that p300 is a "writer" for succinylation, and p300-mediated Ksucc promotes cell glycometabolism in lung adenocarcinoma (LUAD). Specifically, our succinylome data revealed that EP300 deficiency leads to the systemic reduction of Ksucc, and 79.55% of the p300-succinylated proteins were found in the cytoplasm, which were primarily enriched in the carbohydrate metabolism process. Interestingly, deleting EP300 led to a notable decrease in Ksucc levels on several glycolytic enzymes, especially Phosphoglycerate Kinase 1 (PGK1). Mutation of the succinylated site of PGK1 notably hindered cell glycolysis and lactic acid excretion. Metabolomics in vivo indicated that p300-caused metabolic reprogramming was mainly attributed to the altered carbohydrate metabolism. In addition, 89.35% of LUAD patients exhibited cytoplasmic localization of p300, with higher levels in tumor tissues than adjacent normal tissues. High levels of p300 correlated with advanced tumor stages and poor prognosis of LUAD patients. Briefly, we disclose the activity of p300 to catalyze succinylation, which contributes to cell glucose metabolic reprogramming and malignant progression of lung cancer.

The lysine 2-hydroxyisobutyrylome of Helicobacter pylori: Indicating potential roles of lysine 2-hydroxyisobutyrylation in the bacterial metabolism

Helicobacter pylori (H. pylori) is a pathogen which colonizes the stomach, causing ulcers, chronic gastritis and other related diseases. Protein post-translational modifications (PTMs) in bacteria mainly include glycosylation, ubiquitination, nitrosylation, methylation, phosphorylation and acetylation, all of which have divergent functions in the physiology and pathology of the bacterium. Lysine 2-hydroxyisobutyrylation (Khib) is a newly discovered type of PTM in recent years in some kinds of organisms, and this PTM is involved in the regulation of a variety of metabolic process, such as bacterial glucose metabolism, lipid metabolism and protein synthesis. This study performed the first qualitative lysine 2-hydroxyisobutyrylome in H. pylori, and a total of 4419 Khib sites in 812 proteins were identified. The results show that Khib sites are mainly located in the key functional regions or active domains of proteins involved in nickel-trafficking, energy production, virulence factors, anti-oxidation, metal resistance, and ribosome biosynthesis in H. pylori. The study presented here provides new hints in the metabolism and pathology of H. pylori and the proteins with Khib modification may be potentially promising targets for the further development of antibiotics, especially considering the high occurrence of treatment failure of H. pylori failure due to development of antibiotics-resistance.

Potential Mechanisms Underlying the Therapeutic Roles of Gancao fuzi Decoction in Cold-dampness Obstruction Syndrome-type Knee Osteoarthritis

BACKGROUND

The key active components and potential molecular mechanism of Gancao Fuzi decoction (GFD) in the treatment of cold-dampness obstruction-type knee osteoarthritis (KOA) remain unclear.

OBJECTIVES

To explore the mechanism of GFD in the treatment of cold-dampness obstruction syndrome-type KOA by network pharmacology.

METHODS

The potential active components and targets of the four herbs in GFD (Fuzi, Guizhi, Baizhu, and Gancao) were screened using the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database. The targets of KOA were obtained with the Comparative Toxicogenomics Database (CTD), the GeneCards database, and the DisGeNET database, and the common targets of the drugs and disease were ultimately obtained. Cytoscape (v.3.7.1) was used to draw the active component-target network, and the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) (v.11.0) database was used to construct the protein interaction network. The Database for Annotation, Visualization, and Integrated Discovery (DAVID) was used for the Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the intersecting targets.

RESULTS

A total of 102 potential active components and 208 targets of GFD in the treatment of cold-dampness obstruction syndrome-type KOA were screened. GFD treatment was found to be closely related to many inflammatory signalling pathways in the treatment of KOA.

CONCLUSION

The effect of GFD on cold-dampness obstruction syndrome-type KOA is mediated by multicomponent, multitarget, and multichannel mechanisms, which provides the basis for further experimental study of its pharmacodynamic material basis and mechanism.

Genomic and epigenomic insights into the mechanism of cold response in upland cotton (Gossypium hirsutum)

Functional genome research, including gene transcriptional and posttranslational modifications of histones, can benefit greatly from a high-quality genome assembly. Histone modification plays a significant role in modulating the responses to abiotic stress in plants. However, there are limited reports on the involvement of dynamic changes in histone modification in cold stress response in upland cotton. In this study, the genome of an elite accession, YM11, with considerable cold stress tolerance was de novo assembled, which yielded a genome of 2343.06 Mb with a contig N50 of 88.96 Mb, and a total of 73,821 protein-coding gene models were annotated. Comparisons among YM11 and five Gossypium allopolyploid cotton assemblies highlighted a large amount of structural variations and presence/absence variations. We analyzed transcriptome and metabolome changes in YM11 seedlings subjected to cold stress. Using the CUT&Tag method, genome-wide H3K4me3 and H3K9ac modification patterns and effect of histone changes on gene expression were profiled during cold stress. Significant and consistently changing histone modifications and the gene expressions were screened, of which transcription factors (TFs) were highlighted. Our results suggest a positive correlation between the changes in H3K4me3, H3K9ac modifications and cold stress-responsive gene activation. This genome assembly and comprehensive analysis of genome-wide histone modifications and gene expression provide insights into the genomic variation and epigenetic responses to cold stress in upland cotton.

Advancing drug-target interaction prediction: a comprehensive graph-based approach integrating knowledge graph embedding and ProtBert pretraining

BACKGROUND

The pharmaceutical field faces a significant challenge in validating drug target interactions (DTIs) due to the time and cost involved, leading to only a fraction being experimentally verified. To expedite drug discovery, accurate computational methods are essential for predicting potential interactions. Recently, machine learning techniques, particularly graph-based methods, have gained prominence. These methods utilize networks of drugs and targets, employing knowledge graph embedding (KGE) to represent structured information from knowledge graphs in a continuous vector space. This phenomenon highlights the growing inclination to utilize graph topologies as a means to improve the precision of predicting DTIs, hence addressing the pressing requirement for effective computational methodologies in the field of drug discovery.

RESULTS

The present study presents a novel approach called DTIOG for the prediction of DTIs. The methodology employed in this study involves the utilization of a KGE strategy, together with the incorporation of contextual information obtained from protein sequences. More specifically, the study makes use of Protein Bidirectional Encoder Representations from Transformers (ProtBERT) for this purpose. DTIOG utilizes a two-step process to compute embedding vectors using KGE techniques. Additionally, it employs ProtBERT to determine target-target similarity. Different similarity measures, such as Cosine similarity or Euclidean distance, are utilized in the prediction procedure. In addition to the contextual embedding, the proposed unique approach incorporates local representations obtained from the Simplified Molecular Input Line Entry Specification (SMILES) of drugs and the amino acid sequences of protein targets.

CONCLUSIONS

The effectiveness of the proposed approach was assessed through extensive experimentation on datasets pertaining to Enzymes, Ion Channels, and G-protein-coupled Receptors. The remarkable efficacy of DTIOG was showcased through the utilization of diverse similarity measures in order to calculate the similarities between drugs and targets. The combination of these factors, along with the incorporation of various classifiers, enabled the model to outperform existing algorithms in its ability to predict DTIs. The consistent observation of this advantage across all datasets underlines the robustness and accuracy of DTIOG in the domain of DTIs. Additionally, our case study suggests that the DTIOG can serve as a valuable tool for discovering new DTIs.

The draft genome of the Temminck's tragopan (Tragopan temminckii) with evolutionary implications

BACKGROUND

High-quality genome data of birds play a significant role in the systematic study of their origin and adaptive evolution. The Temminck's tragopan (Tragopan temminckii) (Galliformes, Phasianidae), a larger pheasant, is one of the most abundant and widely distributed species of the genus Tragopan, and was defined as class II of the list of national key protected wild animals in China. The absence of a sequenced genome has restricted previous evolutionary trait studies of this taxa.

RESULTS

The whole genome of the Temminck's tragopan was sequenced using Illumina and PacBio platform, and then de novo assembled and annotated. The genome size was 1.06 Gb, with a contig N50 of 4.17 Mb. A total of 117.22 Mb (11.00%) repeat sequences were identified. 16,414 genes were predicted using three methods, with 16,099 (98.08%) annotated as functional genes based on five databases. In addition, comparative genome analyses were conducted across 12 Galliformes species. The results indicated that T. temminckii was the first species to branch off from the clade containing Lophura nycthemera, Phasianus colchicus, Chrysolophus pictus, Syrmaticus mikado, Perdix hodgsoniae, and Meleagris gallopavo, with a corresponding divergence time of 31.43 million years ago (MYA). Expanded gene families associated with immune response and energy metabolism were identified. Genes and pathways associated with plumage color and feather development, immune response, and energy metabolism were found in the list of positively selected genes (PSGs).

CONCLUSIONS

A genome draft of the Temminck's tragopan was reported, genome feature and comparative genome analysis were described, and genes and pathways related to plumage color and feather development, immune response, and energy metabolism were identified. The genomic data of the Temminck's tragopan considerably contribute to the genome evolution and phylogeny of the genus Tragopan and the whole Galliformes species underlying ecological adaptation strategies.

Two High-Quality Cygnus Genome Assemblies Reveal Genomic Variations Associated with Plumage Color

As an exemplary model for examining molecular mechanisms responsible for extreme phenotypic variations, plumage color has garnered significant interest. The Cygnus genus features two species, Cygnus olor and Cygnus atratus, that exhibit striking disparities in plumage color. However, the molecular foundation for this differentiation has remained elusive. Herein, we present two high-quality genomes for C. olor and C. atratus, procured using the Illumina and Nanopore technologies. The assembled genome of C. olor was 1.12 Gb in size with a contig N50 of 26.82 Mb, while its counterpart was 1.13 Gb in size with a contig N50 of 21.91 Mb. A comparative analysis unveiled three genes (TYR, SLC45A2, and SLC7A11) with structural variants in the melanogenic pathway. Notably, we also identified a novel gene, PWWP domain containing 2A (PWWP2A), that is related to plumage color, for the first time. Using targeted gene modification analysis, we demonstrated the potential genetic effect of the PWWP2A variant on pigment gene expression and melanin production. Finally, our findings offer insight into the intricate pattern of pigmentation and the role of polygenes in birds. Furthermore, these two high-quality genome references provide a comprehensive resource and perspective for comparative functional and genetic studies of evolution within the Cygnus genus.

Discussion on the Antipruritic Mechanism of Qiwei Antipruritic Based on Network Pharmacology and Molecular Docking Technology

Objective

To explore the mechanism of Qiwei antipruritic by using network pharmacology and molecular docking technology.

Methods

The components and related targets of Qiwei antipruritic were screened by using the traditional Chinese medicine system pharmacology database (TCMSP and symmap databases). GeneCards and OMIM databases were used to screen itch-related targets. The protein-protein interaction (PPI) network between active ingredient targets and pruritus disease targets was constructed using STRING database. Cytoscape 3.8.0 software was used to draw the visualization network of "drug-component-target-signaling pathway" and screen the core targets. Gene ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed using R software. AutoDock vina software was used to perform molecular docking of key targets and their corresponding key components.

Results

There were 44 main components of Qiwei antipruritic compound, 118 corresponding targets and 3869 itch-related genes. A total of 82 predicted targets of Qiwei antipruritic in the treatment of pruritus were obtained. Eleven key targets were screened. Among the 23 KEGG enriched pathways, 12 signaling pathways were related to skin pruritus. Molecular docking results showed that the core components of Qiwei antipruritic, including quercetin, kaempferol, β-sitosterol, stigmasterol, luteolin, and preskimmianine, had good binding ability with ESR1, PPARG, IL6, TP53, and EGFR, and the docking scores were all less than -4.

Conclusion

The mechanism of Qiwei antipruritic may be related to histamine activation mechanism, calcium channel mechanism, inhibition of inflammatory signaling pathway, inhibition of neurotransmitters, and regulation of immune pathways. The traditional Chinese medicine compound Qiwei antipruritic can treat clinical pruritus through multiple targets and pathways.

Whole genome structure and resistance genes in carbapenemase-producing multidrug resistant ST378 Klebsiella pneumoniae

BACKGROUND

Carbapenemase-producing Klebsiella pneumoniae (CPKP) is one of the most dangerous multidrug-resistant (MDR) pathogens in human health due to its widespread circulation in the nosocomial environment. CPKP carried by companion dogs, which are close to human beings, should be considered a common threat to public health. However, CPKP dissemination through companion animals is still under consideration of major diagnosis and surveillance systems.

METHODS

Two CPKP isolates which were genotyped to harbor bla NDM-5-encoding IncX3 plasmids, were subjected to the whole-genome study. Whole bacterial DNA was isolated, sequenced, and assembled with Oxford Nanopore long reads and corrected with short reads from the Illumina NovaSeq 6000 platform. The whole-genome structure and positions of antimicrobial resistance (AMR) genes were identified and visualized using CGView. Worldwide datasets were downloaded from the NCBI GenBank database for whole-genome comparative analysis. The whole-genome phylogenetic analysis was constructed using the identified whole-chromosome SNP sites from K. pneumoniae HS11286.

RESULTS

As a result of the whole-genome identification, 4 heterogenous plasmids and a single chromosome were identified, each carrying various AMR genes. Multiple novel structures were identified from the AMR genes, coupled with mobile gene elements (MGE). The comparative whole-genome epidemiology revealed that ST378 K. pneumoniae is a novel type of CPKP, carrying a higher prevalence of AMR genes.

CONCLUSIONS

The characterized whole-genome analysis of this study shows the emergence of a novel type of CPKP strain carrying various AMR genes with variated genomic structures. The presented data in this study show the necessity to develop additional surveillance programs and control measures for a novel type of CPKP strain.

A chromosome-level genome assembly of the Asian giant softshell turtle Pelochelys cantorii

The Asian giant softshell turtle Pelochelys cantorii is one of the largest aquatic turtles in China and has been designated a First Grade Protected Animal in China. To advance conservation research, a combination of Illumina short-read, PacBio long-read, and Hi-C scaffolding technologies was used to develop a high-quality chromosome-level genome assembly for P. cantorii. A total of 262.77 Gb of clean data were produced (121.6 × depth) and then the genome was assembled into 2.16 Gb with a contig N50 of 41.44 Mb and scaffold N50 length of 120.17 Mb, respectively. Moreover, about 99.98% assembly genome sequences were clustered and ordered onto 33 pseudochromosomes. Genome annotation revealed that 21,833 protein-coding genes were predicted, and 96.40% of them were annotated. This new chromosome-level assembly will be an enabling resource for genetic and genomic studies to support fundamental insight into P. cantorii biology.

Multi-omics analysis reveals the evolutionary origin of diterpenoid alkaloid biosynthesis pathways in Aconitum

Diterpenoid alkaloids (DAs) have been often utilized in clinical practice due to their analgesic and anti-inflammatory properties. Natural DAs are prevalent in the family Ranunculaceae, notably in the Aconitum genus. Nevertheless, the evolutionary origin of the biosynthesis pathway responsible for DA production remains unknown. In this study, we successfully assembled a high-quality, pseudochromosome-level genome of the DA-rich species Aconitum vilmorinianum (A. vilmorinianum) (5.76 Gb). An A. vilmorinianum-specific whole-genome duplication event was discovered using comparative genomic analysis, which may aid in the evolution of the DA biosynthesis pathway. We identified several genes involved in DA biosynthesis via integrated genomic, transcriptomic, and metabolomic analyses. These genes included enzymes encoding target ent-kaurene oxidases and aminotransferases, which facilitated the activation of diterpenes and insertion of nitrogen atoms into diterpene skeletons, thereby mediating the transformation of diterpenes into DAs. The divergence periods of these genes in A. vilmorinianum were further assessed, and it was shown that two major types of genes were involved in the establishment of the DA biosynthesis pathway. Our integrated analysis offers fresh insights into the evolutionary origin of DAs in A. vilmorinianum as well as suggestions for engineering the biosynthetic pathways to obtain desired DAs.

A chromosome-level reference genome of the Antarctic blackfin icefish Chaenocephalus aceratus

The blackfin Icefish (Chaenocephalus aceratus) belongs to the family Channichthyidae and the suborder Notothenioidei which lives in the Antarctic. We corrected the mis-scaffolds in the previous linkage map results by Hi-C analysis to obtain improved results for chromosome-level genome assembly. The final assembly analysis resulted in a total of 3,135 scaffolds, a genome size of 1,065.72 Mb, and an N50 of 33.46 Mb. 820.24 Mb, representing 88.88% of the total genome, is anchored to 24 chromosomes. The final gene set of 38,024 genes, including AFGPs, was annotated using RNA evidence, proteins, and ab-initio predictions. The complete percentage of BUSCO analysis is 92.7%. In this study, we aim to contribute to the study of polar fishes by improving the genome sequences of the blackfin icefish with the AFGP genes belonging to the Notothenoidei.

Comparative transcriptomics provides insights into the pathogenic immune response of brown leaf spots in weeping forsythia

Weeping forsythia is an important ornamental, ecological and medicinal plant. Brown leaf spots limit the large-scale production of weeping forsythia as a medicinal crop. Alternaria alternata is a pathogen causing brown leaf spots in weeping forsythia; however, its pathogenesis and the immune response mechanisms of weeping forsythia remain unclear. In this study, we identified two mechanisms based on morphological anatomy, physiological indexes and gene expression analyses. Our results showed that A. alternata induced leaf stomata to open, invaded the mesophyll, dissolved the cell wall, destroyed the cell membrane and decreased the number of chloroplasts by up-regulating the expression of auxin-activated signaling pathway genes. Alternaria alternata also down-regulated iron-ion homeostasis and binding-related genes, which caused an increase in the levels of iron ions and reactive oxygen species in leaves. These processes eventually led to programmed cell death, destroying palisade and spongy tissues and causing the formation of iron rust spots. Alternaria alternata also caused defense and hypersensitive responses in weeping forsythia through signaling pathways mediated by flg22-like and elf18-like polypeptides, ethylene, H2O2 and bacterial secretion systems. Our study provides a theoretical basis for the control of brown leaf spots in weeping forsythia.

Genomic molecular epidemiology of carbapenemase-producing Escherichia coli ST410 isolates by complete genome analysis

The circulation of carbapenemase-producing Escherichia coli (CPEC) in our society is a serious concern for vulnerable patients in nosocomial environments. However, the genomic epidemiology of the circulation of CPEC bacteria among companion animals remains largely unknown. In this study, epidemiological analysis was conducted using complete genome identification of CPEC ST410 isolates obtained from companion animals. To estimate the genomic distance and relatedness of the isolates, a total of 37 whole-genome datasets of E. coli ST410 strains were downloaded and comparatively analysed. As a result of the analysis, the genomic structure of the chromosomes and plasmids was identified, revealing the genomic positions of multiple resistance and virulence genes. The isolates in this study were grouped into the subclade H24/RxC, with fimH24, and substituted quinolone resistance-determining regions (QRDRs) and multiple beta-lactamases, including extended-spectrum β-lactamase (ESBL) and carbapenemase. In addition, the in silico comparison of the whole-genome datasets revealed unidentified ST410 H24/Rx subgroups, including either high pathogenicity islands (HPIs) or H21 serotypes. Considering the genetic variations and resistance gene dissemination of the isolates carried by companion animals, future approaches for preventive measurement must include the "One Health" perspective for public health in our society.

Expression dosage effects of a small number of genes after the artificial doubling of weeping forsythia

Whole genome doubling (WGD) plays a critical role in plant evolution, yet the mechanisms underlying the maintenance of overall equilibrium following an artificial doubling event, as well as its impact on phenotype and adaptability, remain unclear. By comparing the gene expression of naturally occurring weeping forsythia diploids and colchicine-induced autotetraploids under normal growth conditions and cold stress, we identified gene expression dosage responses resulting from ploidy change. Only a small proportion of effectively expressed genes showed dosage effect, and most genes did not exhibit significant expression differences. However, the genes that showed expression dosage effect were largely random. The autotetraploids had slower overall growth rates, possibly resulting from negative gene dosage effects on zeatin synthesis and multiple metabolic delays caused by other negative dosage genes. Our comparative analysis of cold response genes in diploids and autotetraploids revealed that genes related to "response to abscisic acid" and "cold acclimation" were key factors contributing to greater cold tolerance in the autotetraploids. In particular, gene expression related to "cold acclimation" might mitigate the effects of cold stress. Taken together, our findings suggested that overall gene expression equilibrium following WGD of weeping forsythia autotetraploids was achieved through the inactivation of the majority of duplicated genes. Our research provides new insights into the mechanisms regulating expression dosage balance following polyploidization events.

Potential Role of Lysine Acetylation and Autophagy in Brown Film Formation and Postripening of Lentinula edodes Mycelium

Lentinula edodes is one of the most widely cultivated edible mushrooms in the world. When cultivated in sawdust, the surface mycelium of L. edodes needs a long postripening stage wherein it forms a brown film (BF) by secreting and accumulating pigments. BF formation is critical for the high quality and yield of fruiting bodies. Protein lysine acetylation (KAC) is an important post-translational modification that regulates growth and development. Previous studies have shown that deacetylase levels are significantly increased during BF formation in the postripening stage of L. edodes. The aim of this study was to assess the role of protein acetylation during BF formation. To this end, we compared the acetylome of L. edodes mycelia before and after BF formation using anti-acetyl antibody-based label-free quantitative proteomics. We identified 5,613 acetylation sites in 1,991 proteins, and quantitative information was available for 4,848 of these sites in 1,815 proteins. Comparative acetylome analysis showed that the modification of 699 sites increased and that of 562 sites decreased during BF formation. Bioinformatics analysis of the differentially acetylated proteins showed significant enrichment in the tricarboxylic acid (TCA) cycle and proteasome pathways. Furthermore, functional assays showed that BF formation is associated with significant changes in the activities of proteasome, citrate synthase, and isocitrate dehydrogenase. Consistent with this hypothesis, the lysine deacetylase inhibitor trichostatin (TSA) delayed autophagy and BF formation in L. edodes. Taken together, KAC and autophagy play important roles in the mycelial BF formation and postripening stage of L. edodes. IMPORTANCE Mycelial BF formation and postripening of L. edodes affects the quality and quantity of its edible fruiting bodies. In this study, we explored the role of protein KAC in this biological process, with the aim of optimizing the cultivation and yield of L. edodes.

Genomic Resequencing Unravels the Genetic Basis of Domestication, Expansion, and Trait Improvement in Morus Atropurpurea

Mulberry is an economically important plant in the sericulture industry and traditional medicine. However, the genetic and evolutionary history of mulberry remains largely unknown. Here, this work presents the chromosome-level genome assembly of Morus atropurpurea (M. atropurpurea), originating from south China. Population genomic analysis using 425 mulberry accessions reveal that cultivated mulberry is classified into two species, M. atropurpurea and M. alba, which may have originated from two different mulberry progenitors and have independent and parallel domestication in north and south China, respectively. Extensive gene flow is revealed between different mulberry populations, contributing to genetic diversity in modern hybrid cultivars. This work also identifies the genetic architecture of the flowering time and leaf size. In addition, the genomic structure and evolution of sex-determining regions are identified. This study significantly advances the understanding of the genetic basis and domestication history of mulberry in the north and south, and provides valuable molecular markers of desirable traits for mulberry breeding.

Transcriptomic analysis of maize uncovers putative genes involved in metabolic detoxification under four safeners treatment

Isoxadifen-ethyl (IDF) and cyprosulfamide (CSA) can effectively protect maize from nicosulfuron (NIC) injury, while mefenpyr-diethyl (MPR) and fenchlorazole-ethyl (FCO) did not. Their chemical diversity and requirement to use them in combination with the corresponding herbicides suggest that their elicitation of gene expression are complex and whether it is associated with the safening activity remains elusive. In this study, our first objective was to determine whether or not the ability of four safeners to enhance the metabolic rate of nicosulfuron. It was found that nicosulfuron degradation in maize was accelerated by IDF and CSA, but not by FCO and MPR. Transcriptomic analysis showed that the number of genes induced by IDF and CSA were larger than that induced by FCO and MPR. Overall, 34 genes associated with detoxification were identified, including glutathione S-transferase (GST), cytochrome P450 (CYP450), UDP-glucosyltransferase (UGT), transporter and serine. Moreover, 14 detoxification genes were screened for further verification by real-time PCR in two maize inbred lines. Two maize inbred lines exhibited high expression levels of four genes (GST31, GST39, AGXT2 and ADH) after IDF treatment. GST6, GST19, MATE, SCPL18 and UF3GT were specifically up-regulated in telerant maize inbred line under IDF and IDF + NIC treatments. Seven genes, namely GST31, GST6, GST19, UF3GT, MATE, ADH and SCPL18, are induced by IDF and CSA to play a vital role in regulating the detoxification process of NIC. Accordingly, the GST activity in maize was accelerated by IDF and CSA, but not by FCO and MPR. This result is consistent with transcriptome and metabolic data.These results indicate that the mitigation of NIC damage is associated with enhanced herbicide metabolism. IDF and CSA were more effective in protecting maize from NIC injury due to their ability to enhance the detoxification of specific types of herbicides, compared to FCO and MPR. The chemical specificity of four safeners is attributed to the up-regulated genes related to the detoxification pathway.

Active components and molecular mechanisms of Sagacious Confucius' Pillow Elixir to treat cognitive impairment based on systems pharmacology

BACKGROUND

Sagacious Confucius' Pillow Elixir (SCPE) is a common clinical prescription to treat cognitive impairment (CI) in East Asia.

OBJECTIVE

To predict the active components of SCPE, identify the associated signaling pathway, and explore the molecular mechanism using systems pharmacology and an animal study.

METHODS

Systems pharmacology and Python programming language-based molecular docking were used to select and analyze the active components and targets. Senescence-accelerated prone 8 mice were used as a CI model. The molecular mechanism was evaluated using the water maze test, neuropathological observation, cerebrospinal fluid microdialysis, and Western blotting.

RESULTS

Thirty active components were revealed by screening relevant databases and performing topological analysis. Additionally, 376 differentially expressed genes for CI were identified. Pathway enrichment analysis, protein-protein interaction (PPI) network analysis and molecular docking indicated that SCPE played a crucial role in modulating the PI3K/Akt/mTOR signaling pathway, and 23 SCPE components interacted with it. In the CI model, SCPE improved cognitive function, increased the levels of the neurotransmitter 5-hydroxytryptamine (5-HT) and metabolite 5-hydroxyindole acetic acid (5-HIAA), ameliorated pathological damage and regulated the PI3K/AKT/mTOR signaling pathway. SCPE increased the LC3-II/LC3-I, p-PI3K p85/PI3K p85, p-AKT/AKT, and p-mTOR/mTOR protein expression ratios and inhibited P62 expression in the hippocampal tissue of the CI model.

CONCLUSION

Our study revealed that 23 active SCPE components improve CI by increasing the levels of the neurotransmitter 5-HT and metabolite 5-HIAA, suppressing pathological injury and regulating the PI3K/Akt/mTOR signaling pathway to improve cognitive function.

Assembly of Genome and Resequencing Provide Insights into Genetic Differentiation between Parents of Hulong Hybrid Grouper (Epinephelus fuscoguttatus ♀ × E. lanceolatus ♂)

The Hulong hybrid grouper was bred from the brown-marbled grouper (Epinephelus fuscoguttatus) ♀ and the giant grouper (E. lanceolatus) ♂, combining the advantageous traits of both parents. Possessing an excellent performance, this hybrid's cultivation promotes the development of the grouper industry. Its male parent, the giant grouper, possesses the fastest growth and the largest body size among all coral-reef-dwelling fish. This species is not only an economically important species in marine aquaculture, but it is also an ideal male parent in the interspecific crossing of grouper species. In the present study, a high-quality chromosome-level genome of the giant grouper was constructed with a total length of 1.06 Gb, consisting of 24 chromosomes and 69 scaffolds. To analyze the genetic differences between the parents of the Hulong hybrid grouper, the structural variations (SVs) between both parental genomes were detected, and a total of 46,643 SVs were obtained. High-quality SNPs were identified from resequencing data. There were significant differences between the two genomes, and the average FST reached 0.685. A total of 234 highly differentiated regions were detected with an FST > 0.9. The protein-coding genes involved in SVs and highly differentiated regions were significantly enriched in metabolic pathways, including fatty metabolism, carbohydrate metabolism, amino acid metabolism and the TCA cycle. These genes may be related to the differences in feeding preferences and the ability to digest carbohydrates between the two grouper species under natural conditions. In addition, protein-coding genes related to the cell cycle and p53-signaling pathway were also detected. These genes may play important roles in the regulation of body size and growth performance. This research provides genomic resources for further breeding works and evolutionary analyses.

Turnover number predictions for kinetically uncharacterized enzymes using machine and deep learning

The turnover number k_cat, a measure of enzyme efficiency, is central to understanding cellular physiology and resource allocation. As experimental k_cat estimates are unavailable for the vast majority of enzymatic reactions, the development of accurate computational prediction methods is highly desirable. However, existing machine learning models are limited to a single, well-studied organism, or they provide inaccurate predictions except for enzymes that are highly similar to proteins in the training set. Here, we present TurNuP, a general and organism-independent model that successfully predicts turnover numbers for natural reactions of wild-type enzymes. We constructed model inputs by representing complete chemical reactions through differential reaction fingerprints and by representing enzymes through a modified and re-trained Transformer Network model for protein sequences. TurNuP outperforms previous models and generalizes well even to enzymes that are not similar to proteins in the training set. Parameterizing metabolic models with TurNuP-predicted k_cat values leads to improved proteome allocation predictions. To provide a powerful and convenient tool for the study of molecular biochemistry and physiology, we implemented a TurNuP web server.

Identification of ORM1, vWF, SPARC, and PPBP as immune-related proteins involved in immune thrombocytopenia by quantitative LC-MS/MS

BACKGROUND

Immune thrombocytopenia (ITP) is a common autoimmune disease characterized by loss of immune tolerance to platelet autoantigens leading to excessive destruction and insufficient production of platelets.

METHOD

Quantitative liquid chromatography tandem mass spectrometry (LC-MS/MS) was performed to detect the differentially expressed proteins in bone marrow samples from active ITP patients and normal controls.

RESULT

Our bioinformatic analysis identified two upregulated proteins (ORM1 and vWF) and two downregulated proteins (PPBP and SPARC) related to immune function. The four proteins were all found to be related to the tumor necrosis factor (TNF) -α signalling pathway and involved in the pathogenesis of ITP in KEGG pathway analysis.

CONCLUSION

Bioinformatics analysis identified differentially expressed proteins in bone marrow that are involved in the TNF-α signalling pathway and are related to the activation of immune function in ITP patients. These findings could provide new ideas for research on the loss of immune tolerance in ITP patients.

Identification of drought-tolerant hub genes in Iranian KC-2226 genotype of Aegilops tauschii using transcriptomic analysis

Aegilops tauschii, as a donor of D genome to the bread wheat with a valuable source of resistance to different biotic and abiotic stresses, is used to improve the quality of wheat cultivars. Every genotype has a specific genetic content, the investigation of which can lead to the identification of useful genes such as stress tolerance genes, including drought. Therefore, 23 genotypes of Ae. tauschii were selected to evaluate their morphological and physiological traits under greenhouse conditions. Among them, a superior tolerant genotype (KC-2226) was chosen for transcriptomic analysis. Our result showed that 5007 and 3489 genes were deferentially up- and downregulated, respectively. Upregulated genes were involved in photosynthesis, glycolysis/gluconeogenesis, and amino acid biosynthesis whereas downregulated genes were often engaged in DNA synthesis, replication, repair and topological changes. The result of protein-protein interaction network analysis showed that AT1G76550 (1.46), AT1G20950 (1.42), IAR4 (1.19), and PYD2 (1.16) among upregulated genes and THY-1 (44), PCNA1 (41) and TOPII (22) among down-regulated genes had the highest interactions with other genes. In conclusion, Ae. tauschii employs elevated transcription of specific genes involved in photosynthesis, glycolysis and gluconeogenesis and amino acid biosynthesis pathways rather than genes active in DNA synthesis and repair to provide the energy needed for the plant to survive under stress conditions.

Comparative genome and metabolome analyses uncover the evolution and flavonoid biosynthesis between Apocynum venetum and Apocynum hendersonii

Apocynum species have great application prospects in textile and phytoremediation of saline soil, are rich in flavonoids, and possess medicinal significance. Here, we report the draft genomes of Apocynum venetum and Apocynum hendersonii, and elucidate their evolutionary relationship. The high synteny and collinearity between the two suggested that they have experienced the same WGD event. Comparative analysis revealed that flavone 3-hydroxylase (ApF3H) and differentially evolved flavonoid 3-O-glucosyl transferase (ApUFGT) genes are critical for determining natural variation in flavonoid biosynthesis between the species. Overexpression of ApF3H-1 enhanced the total flavonoid content and promoted the antioxidant capacity of transformed plants compared to the wild-type. ApUFGT5 and 6 explained the diversification of flavonoids or their derivatives. These data provide biochemical insight and knowledge on the genetic regulation of flavonoid biosynthesis, supporting the adoption of these genes in breeding programs aimed at the multipurpose utilization of the plants.

A general model to predict small molecule substrates of enzymes based on machine and deep learning

For most proteins annotated as enzymes, it is unknown which primary and/or secondary reactions they catalyze. Experimental characterizations of potential substrates are time-consuming and costly. Machine learning predictions could provide an efficient alternative, but are hampered by a lack of information regarding enzyme non-substrates, as available training data comprises mainly positive examples. Here, we present ESP, a general machine-learning model for the prediction of enzyme-substrate pairs with an accuracy of over 91% on independent and diverse test data. ESP can be applied successfully across widely different enzymes and a broad range of metabolites included in the training data, outperforming models designed for individual, well-studied enzyme families. ESP represents enzymes through a modified transformer model, and is trained on data augmented with randomly sampled small molecules assigned as non-substrates. By facilitating easy in silico testing of potential substrates, the ESP web server may support both basic and applied science.

Genome-Wide Identification and Expression Analysis of the Cytochrome P450 Gene Family in Bemisia tabaci MED and Their Roles in the Insecticide Resistance

The whitefly, Bemisia tabaci MED (Hemiptera: Aleyrodidae), is an omnivorous agricultural pest, which causes huge economic losses to agriculture and is highly resistant to many pesticides. The overexpression of cytochrome P450 may play an important role in host adaptation and insecticide resistance in B. tabaci MED. Therefore, the present study systematically analyzed the cytochrome P450 gene family at the genome-wide level to understand its function in B. tabaci MED. Our analysis identified 58 cytochrome P450 genes in B. tabaci MED, among which 24 were novel. Phylogenetic analysis revealed broad functional and species-specific diversification in B. tabaci MED P450, suggesting the role of multiple P450 genes in detoxifying. Reverse transcription-real time quantitative PCR (RT-qPCR) showed that CYP4CS2, CYP4CS5, CYP4CS6, CYP4CS8, CYP6DW4, CYP6DW5, CYP6DW6, CYP6DZ8, and CYP6EN1 genes increased significantly after two days of exposure to imidacloprid. Interestingly, all nine genes belonged to the CYP4 and CYP6 families. A decrease in the expression of five genes (CYP6DW4, CYP6DW5, CYP6DW6, CYP6DZ8, and CYP4CS6) via RNA interference (RNAi) resulted in a significant increase in the mortalities of whiteflies when exposed to imidacloprid. These results indicate that the overexpression of the P450 genes may play an essential role in imidacloprid tolerance of B. tabaci MED. Thus, the present study provides basic information on P450 genes in B. tabaci MED, which will further help elucidate the insecticide resistance mechanism in the agricultural pest whitefly.

The genome of Orychophragmus violaceus provides genomic insights into the evolution of Brassicaceae polyploidization and its distinct traits

Orychophragmus violaceus, referred to as "eryuelan" (February orchid) in China, is an early-flowering ornamental plant. The high oil content and abundance of unsaturated fatty acids in O. violaceus seeds make it a potential high-quality oilseed crop. Here, we generated a whole-genome assembly for O. violaceus using Nanopore and Hi-C sequencing technologies. The assembled genome of O. violaceus was ∼1.3 Gb in size, with 12 pairs of chromosomes. Through investigation of ancestral genome evolution, we determined that the genome of O. violaceus experienced a tetraploidization event from a diploid progenitor with the translocated proto-Calepineae karyotype. Comparisons between the reconstructed subgenomes of O. violaceus identified indicators of subgenome dominance, indicating that subgenomes likely originated via allotetraploidy. O. violaceus was phylogenetically close to the Brassica genus, and tetraploidy in O. violaceus occurred approximately 8.57 million years ago, close in time to the whole-genome triplication of Brassica that likely arose via an intermediate tetraploid lineage. However, the tetraploidization in Orychophragmus was independent of the hexaploidization in Brassica, as evidenced by the results from detailed phylogenetic analyses and comparisons of the break and fusion points of ancestral genomic blocks. Moreover, identification of multi-copy genes regulating the production of high-quality oil highlighted the contributions of both tetraploidization and tandem duplication to functional innovation in O. violaceus. These findings provide novel insights into the polyploidization evolution of plant species and will promote both functional genomic studies and domestication/breeding efforts in O. violaceus.

A de novo assembled genome of the Tibetan Partridge (Perdix hodgsoniae) and its high-altitude adaptation

The Tibetan Partridge (Perdix hodgsoniae) is an endemic species distributed in high-altitude areas of 3600-5600 m on the Qinghai-Tibet Plateau. To explore how the species is adapted to the high elevation environment, we assembled a draft genome based on both the Illumina and PacBio sequencing platforms with its population genetics and genomics analysis. In total, 134.74 Gb short reads and 30.81 Gb long reads raw data were generated. The 1.05-Gb assembled genome had a contig N50 of 4.56 Mb, with 91.94% complete BUSCOs. The 17 457 genes were annotated, and 11.35% of the genome was composed of repeat sequences. The phylogenetic tree showed that P. hodgsoniae was located at the basal position of the clade, including Golden Pheasant (Chrysolophus pictus), Common Pheasant (Phasianus colchicus), and Mikado Pheasant (Syrmaticus mikado). We found that 1014, 2595, and 2732 of the 6641 one-to-one orthologous genes were under positive selection in P. hodgsoniae, detected using PAML, BUSTED, and aBSREL programs, respectively, of which 965 genes were common under positive selection with 3 different programs. Several positively selected genes and immunity pathways relevant to high-altitude adaptation were detected. Gene family evolution showed that 99 gene families experienced significant expansion events, while 6 gene families were under contraction. The total number of olfactory receptor genes was relatively low in P. hodgsoniae. Genomic data provide an important resource for a further study on the evolutionary history of P. hodgsoniae, which provides a new insight into its high-altitude adaptation mechanisms.

Integrated proteome and acetylome analyses provide novel insights into early somatic embryogenesis of Dimocarpus longan

Longan (Dimocarpus longan) is a precious subtropical fruit with high nutritional value. The somatic embryogenesis (SE) affects the quality and yield of fruit. Apart from clonal propagation, SE has extensive applications in genetic improvement and mutation. Thus, understanding the molecular basis of embryogenesis in longan will help to develop strategies for mass production of quality planting material. Lysine acetylation (Kac) plays an important role in diverse cellular processes, but limited knowledge is available regarding acetylation modifications in plant early SE. In this study, the proteome and acetylome of longan embryogenic callus (ECs) and globular embryos (GEs) were investigated. In total, 7232 proteins and 14,597 Kac sites were identified, and this resulted in the discovery of 1178 differentially expressed proteins and 669 differentially expressed acetylated proteins. KEGG and GO analysis showed that glucose metabolism, carbon metabolism, fatty acid degradation, and oxidative phosphorylation pathways were influenced by Kac modification. Furthermore, sodium butyrate (Sb, a deacetylase inhibitor) led to reduced the proliferation and delayed the differentiation of ECs by regulating the homeostasis of reactive oxygen species (ROS) andindole-3-acetic acid (IAA). Our study provides a comprehensive proteomic and acetylomic analysis to aid in understanding the molecular mechanisms involved in early SE, representing a potential tool for genetic improvement of longan.

Comparative analysis of cold-responsive genes under short-term cold stimulation and cold-adaptive genes under long-term heterogeneous environments reveals a cold adaptation mechanism in weeping forsythia

Identifying cold-related genes can provide insights into the cold adaptation mechanism of weeping forsythia. In this study, we compared the changes in gene expressions and physiological and biochemical indices under short-term cold stimulation with the changes in gene sequences under a long-term heterogeneous environment to investigate the cold adaptation mechanism in weeping forsythia. The data of adaptive gene sequence changes, e.g., single nucleotide polymorphisms, were obtained from previous landscape genomics studies. The physiological and biochemical indicators and transcriptome results showed that weeping forsythia initiated a series of programs, including increasing cell osmotic pressures, scavenging ROS, activating the defense mechanism that crosses with pathogen infection, and upregulating CBF/DREB1 transcription factor 1, to cope with short-term cold stress. A reanalysis of landscape genomic data suggested that weeping forsythia responded to long-term heterogeneous cold stress by the differentiation of genes related to synthesis of aromatic substances and adenosine triphosphate. Our results supported the hypothesis that the adaptation mechanisms of species to short-term environmental stimulation and long-term stress in heterogeneous environments are different. The differences in cold tolerance among populations are not necessarily obtained by changing cold-responsive gene sequences. This study provides new insights into the cold adaptation mechanisms of plants.

Chromosome-scale genomics, metabolomics, and transcriptomics provide insight into the synthesis and regulation of phenols in Vitis adenoclada grapes

Vitis adenoclada is a wild grape unique to China. It exhibits well resistance to heat, humidity, fungal disease, drought, and soil infertility. Here, we report the high-quality, chromosome-level genome assembly of GH6 (V. adenoclada). The 498.27 Mb genome contained 221.78 Mb of transposable elements, 28,660 protein-coding genes, and 481.44 Mb of sequences associated with 19 chromosomes. GH6 shares a common ancestor with PN40024 (Vitis vinifera) from approximately 4.26-9.01 million years ago, whose divergence occurred later than Vitis rotundifolia and Vitis riparia. Widely-targeted metabolome and transcriptome analysis revealed that the profiles and metabolism of phenolic compounds in V. adenoclada varieties significantly were differed from other grape varieties. Specifically, V. adenoclada varieties were rich in phenolic acids and flavonols, whereas the flavan-3-ol and anthocyanin content was lower compared with other varieties that have V. vinifera consanguinity in this study. In addition, ferulic acid and stilbenes content were associated with higher expressions of COMT and STSs in V. adenoclada varieties. Furthermore, MYB2, MYB73-1, and MYB73-2 were presumably responsible for the high expression level of COMT in V. adenoclada berries. MYB12 (MYBF1) was positively correlated with PAL, CHS, FLS and UFGT.Meanwhile, MYB4 and MYBC2-L1 may inhibit the synthesis of flavan-3-ols and anthocyanins in two V. adenoclada varieties (YN2 and GH6). The publication of the V. adenoclada grape genome provides a molecular foundation for further revealing its flavor and quality characteristics, is also important for identifying favorable genes of the East Asian species for future breeding.

Proteomic analysis of protein lysine 2-hydroxyisobutyrylation (Khib) in soybean leaves

BACKGROUND

Protein lysine 2-hydroxyisobutyrylation (K_hib) is a novel post-translational modification (PTM) discovered in cells or tissues of animals, microorganisms and plants in recent years. Proteome-wide identification of K_hib-modified proteins has been performed in several plant species, suggesting that K_hib-modified proteins are involved in a variety of biological processes and metabolic pathways. However, the protein K_hib modification in soybean, a globally important legume crop that provides the rich source of plant protein and oil, remains unclear.

RESULTS

In this study, the K_hib-modified proteins in soybean leaves were identified for the first time using affinity enrichment and high-resolution mass spectrometry-based proteomic techniques, and a systematic bioinformatics analysis of these K_hib-modified proteins was performed. Our results showed that a total of 4251 K_hib sites in 1532 proteins were identified as overlapping in three replicates (the raw mass spectrometry data are available via ProteomeXchange with the identifier of PXD03650). These K_hib-modified proteins are involved in a wide range of cellular processes, particularly enriched in biosynthesis, central carbon metabolism and photosynthesis, and are widely distributed in subcellular locations, mainly in chloroplasts, cytoplasm and nucleus. In addition, a total of 12 sequence motifs were extracted from all identified K_hib peptides, and a basic amino acid residue (K), an acidic amino acid residue (E) and three aliphatic amino acid residues with small side chains (G/A/V) were found to be more preferred around the K_hib site. Furthermore, 16 highly-connected clusters of K_hib proteins were retrieved from the global PPI network, which suggest that K_hib modifications tend to occur in proteins associated with specific functional clusters.

CONCLUSIONS

These findings suggest that K_hib modification is an abundant and conserved PTM in soybean and that this modification may play an important role in regulating physiological processes in soybean leaves. The K_hib proteomic data obtained in this study will help to further elucidate the regulatory mechanisms of K_hib modification in soybean in the future.

annotate_my_genomes: an easy-to-use pipeline to improve genome annotation and uncover neglected genes by hybrid RNA sequencing

BACKGROUND

The advancement of hybrid sequencing technologies is increasingly expanding genome assemblies that are often annotated using hybrid sequencing transcriptomics, leading to improved genome characterization and the identification of novel genes and isoforms in a wide variety of organisms.

RESULTS

We developed an easy-to-use genome-guided transcriptome annotation pipeline that uses assembled transcripts from hybrid sequencing data as input and distinguishes between coding and long non-coding RNAs by integration of several bioinformatic approaches, including gene reconciliation with previous annotations in GTF format. We demonstrated the efficiency of this approach by correctly assembling and annotating all exons from the chicken SCO-spondin gene (containing more than 105 exons), including the identification of missing genes in the chicken reference annotations by homology assignments.

CONCLUSIONS

Our method helps to improve the current transcriptome annotation of the chicken brain. Our pipeline, implemented on Anaconda/Nextflow and Docker is an easy-to-use package that can be applied to a broad range of species, tissues, and research areas helping to improve and reconcile current annotations. The code and datasets are publicly available at https://github.com/cfarkas/annotate_my_genomes.

Quantitative Acetylome Analysis of Soft Wheat Seeds during Artificial Ageing

Lysine acetylation (Kac) is a protein post-translational modification (PTM) widely found in plants that plays vital roles in metabolic pathways. Although seed germination and development are regulated by Kac, its potential function in seed ageing remains to be investigated. Our preliminary study demonstrated that Kac levels were altered during wheat seed artificial ageing. However, its specific role in this process still needs to be elucidated. Here, we performed quantitative acetylation proteomics analysis of soft wheat seeds with different germination rates during artificial ageing. A total of 175 acetylation proteins and 255 acetylation modification sites were remarkably changed. The differentially acetylated proteins were enriched in metabolism; response to harsh intracellular environment, such as ROS; protein storage and processing. Notably, expression, point mutation to mimic Kac by K to Q mutation at K80 and K138, protein purification and enzyme activity detection revealed that the Kac of ROS-scavenging glutathione transferase attenuated its activity, indicating that the defense ability of wheat seeds to stress gradually diminished, and the ageing process was inevitable. Collectively, our data provide a basis for further understanding the roles of Kac in seed ageing and might aid in the development of new techniques to prolong seed viability and food quality.

Identification and Expression Profiling of Two Saudi Arabia Catalase Genes from Wheat and Barley in Response to Abiotic and Hormonal Stresses

Catalase is a crucial enzyme in antioxidant defense systems protecting eukaryotes from oxidative stress. These proteins are present in almost all living organisms and play important roles in controlling plant responses to biotic and abiotic stresses by catalyzing the decomposition of H₂O₂. Despite their importance, little is known about their expression in the majority of monocotyledonous species. Here, we isolated and characterized two novel catalase genes from Triticum turgidum and Hordeum vulgare, designated as TtCAT1 and HvCAT1, respectively. Phylogenetic analysis revealed that TtCAT1 and HvCAT1 presented 492 aa and shared an important identity with other catalase proteins belonging to subfamily 1. Using bioinformatic analysis, we predicted the 3D structure models of TtCAT1 and HvCAT1. Interestingly, analysis showed that the novel catalases harbor a peroxisomal targeting signal (PTS1) located at their C-terminus portion, as shown for other catalase proteins. In addition, this motif is responsible for the in silico peroxisomal localization of both proteins. Finally, RT-qPCR analysis showed that TtCAT1 and HvCAT1 are highly expressed in leaves in normal conditions but faintly in roots. Moreover, both genes are upregulated after the application of different stresses such as salt, osmotic, cold, heavy metal, and hormonal stresses. The positive responses of TtCAT1 and HvCAT1 to the various stimuli suggested that these proteins can help to protect both species against environmental stresses.

The Capparis spinosa var. herbacea genome provides the first genomic instrument for a diversity and evolution study of the Capparaceae family

Background

The caper bush Capparis spinosa L., one of the most economically important species of Capparaceae, is a xerophytic shrub that is well adapted to drought and harsh environments. However, genetic studies on this species are limited because of the lack of its reference genome.

Findings

We sequenced and assembled the Capparis spinosa var. herbacea (Willd.) genome using data obtained from the combination of PacBio circular consensus sequencing and high-throughput chromosome conformation capture. The final genome assembly was approximately 274.53 Mb (contig N50 length of 9.36 Mb, scaffold N50 of 15.15 Mb), 99.23% of which was assigned to 21 chromosomes. In the whole-genome sequence, tandem repeats accounted for 19.28%, and transposable element sequences accounted for 43.98%. The proportion of tandem repeats in the C. spinosa var. herbacea genome was much higher than the average of 8.55% in plant genomes. A total of 21,577 protein-coding genes were predicted, with 98.82% being functionally annotated. The result of species divergence times showed that C. spinosa var. herbacea and Tarenaya hassleriana separated from a common ancestor 43.31 million years ago.

Conclusions

This study reported a high-quality reference genome assembly and genome features for the Capparaceae family. The assembled C. spinosa var. herbacea genome might provide a system for studying the diversity, speciation, and evolution of this family and serve as an important resource for understanding the mechanism of drought and high-temperature resistance.

Cross-Talk of Protein Expression and Lysine Acetylation in Response to TMV Infection in Nicotiana benthamiana

Lysine acetylation (K^ac), a reversible PTM, plays an essential role in various biological processes, including those involving metabolic pathways, pathogen resistance, and transcription, in both prokaryotes and eukaryotes. TMV, the major factor that causes the poor quality of Solanaceae crops worldwide, directly alters many metabolic processes in tobacco. However, the extent and function of K^ac during TMV infection have not been determined. The validation test to detect K^ac level and viral expression after TMV infection and Nicotinamide (NAM) treatment clarified that acetylation was involved in TMV infection. Furthermore, we comprehensively analyzed the changes in the proteome and acetylome of TMV-infected tobacco (Nicotiana benthamiana) seedlings via LC-MS/MS in conjunction with highly sensitive immune-affinity purification. In total, 2082 lysine-acetylated sites on 1319 proteins differentially expressed in response to TMV infection were identified. Extensive bioinformatic studies disclosed changes in acetylation of proteins engaged in cellular metabolism and biological processes. The vital influence of K^ac in fatty acid degradation and alpha-linolenic acid metabolism was also revealed in TMV-infected seedlings. This study first revealed K^ac information in N. benthamiana under TMV infection and expanded upon the existing landscape of acetylation in pathogen infection.

Proteomics and bioinformatics analysis of follicular fluid from patients with polycystic ovary syndrome

Objectives: Polycystic ovary syndrome (PCOS) is a common endocrine and metabolic disorder with heterogeneous manifestations and complex etiology. We used quantitative proteomics analysis based on mass spectrometry to identify the differences in proteomics profiles for follicular fluid obtained from patients with or without PCOS and explore possible mechanisms underlying PCOS.

Methods: Follicular fluid samples were collected from infertile patients with (n = 9) or without (n = 9) PCOS. Total protein was extracted, quantitatively labeled with a tandem mass tag (TMT), and analyzed using liquid chromatography-mass spectrometry (LC‐MS). TMT-based proteomics and bioinformatics analysis were used to determine the differentially expressed proteins (DEPs) and understand the protein networks. The analysis included protein annotation, unsupervised hierarchical clustering, functional classification, functional enrichment and clustering, and protein-protein interaction analysis. Selected DEPs were confirmed by ELISA, and correlation analysis was performed between these DEPs and the clinical characteristics.

Results: In this study, we have identified 1,216 proteins, including 70 DEPs (32 upregulated proteins, 38 downregulated proteins). Bioinformatics analysis revealed that the inflammatory response, complement and coagulation cascades, activation of the immune response, lipid transport, and regulation of protein metabolic processes were co-enriched in patients with PCOS. Based on ELISA results, insulin-like growth factor binding protein 1 (IGFBP1) and apolipoprotein C2 (APOC2) were differentially expressed between patients with and without PCOS. Follicular IGFBP1 showed a positive correlation with the serum levels of high-density lipoprotein cholesterol (HDL-C) (r = 0.3046, p = 0.0419), but negatively correlated with the serum levels of anti-Müllerian hormone (AMH) (r = –0.2924, p = 0.0354) and triglycerides (r = –0.3177, p = 0.0246). Follicular APOC2 was negatively correlated with the serum apolipoprotein A1 (APOA1) levels (r = 0.4509, p = 0.0002).

Conclusion: Our study identified DEPs in the follicular fluid of patients with PCOS. Inflammatory response, complement and coagulation cascades, activation of the immune response, lipid transport, and regulation of protein metabolic process were deregulated in PCOS, which may play essential roles in the pathogenesis of PCOS.