Investigation of Inter- and Intra-Day Variability of Tear Fluid Regarding Flow Rate, Protein Concentration as well as Protein Composition

Purpose

The purpose of this study was to present the determination of inter- and intra-day variations in tear flow rate, and tear fluid protein concentration, as well as protein composition regarding their impact for future biomarker studies.

Methods

Tear fluid was collected noninvasively from 18 healthy subjects by performing Schirmer tests at 4 different time points repetitive in a period of 2 days. The tear flow rate on the Schirmer test strips was measured. Proteins were extracted from strips and quantified using amino acid analysis. Protein composition was analyzed by the strips data-independent (DIA) based mass spectrometry. To exclude any impairments to health, volunteers underwent a detailed neurological as well as an ophthalmological examination.

Results

Whether tear fluid was collected from oculus sinister or oculus dexter did not affect the tear flow rate (P ≈ 0.63) or protein concentration (P ≈ 0.97) of individual subjects. Moreover, protein concentration was independent from the tear volume, so that a change in volume may only influence the total protein amount. When the examination days were compared, investigation of tear flow rate (P ≈ 0.001) and protein concentration (P ≈ 0.0003) indicated significant differences. Further, mass spectrometric analysis of tear fluid revealed 11 differentially regulated proteins when comparing both examination days.

Conclusions

Our findings provide evidence of inter-day variation in tear flow rate, tear proteome concentration, and composition in healthy subjects, suggesting that inter-day variation needs to be taken into consideration in biomarker research of tear fluid. Identified proteins were assigned to functions in the immune response, oxidative and reducing processes, as well as mannose metabolism.

DrugRep-HeSiaGraph: when heterogenous siamese neural network meets knowledge graphs for drug repurposing

BACKGROUND

Drug repurposing is an approach that holds promise for identifying new therapeutic uses for existing drugs. Recently, knowledge graphs have emerged as significant tools for addressing the challenges of drug repurposing. However, there are still major issues with constructing and embedding knowledge graphs.

RESULTS

This study proposes a two-step method called DrugRep-HeSiaGraph to address these challenges. The method integrates the drug-disease knowledge graph with the application of a heterogeneous siamese neural network. In the first step, a drug-disease knowledge graph named DDKG-V1 is constructed by defining new relationship types, and then numerical vector representations for the nodes are created using the distributional learning method. In the second step, a heterogeneous siamese neural network called HeSiaNet is applied to enrich the embedding of drugs and diseases by bringing them closer in a new unified latent space. Then, it predicts potential drug candidates for diseases. DrugRep-HeSiaGraph achieves impressive performance metrics, including an AUC-ROC of 91.16%, an AUC-PR of 90.32%, an accuracy of 84.63%, a BS of 0.119, and an MCC of 69.31%.

CONCLUSION

We demonstrate the effectiveness of the proposed method in identifying potential drugs for COVID-19 as a case study. In addition, this study shows the role of dipeptidyl peptidase 4 (DPP-4) as a potential receptor for SARS-CoV-2 and the effectiveness of DPP-4 inhibitors in facing COVID-19. This highlights the practical application of the model in addressing real-world challenges in the field of drug repurposing. The code and data for DrugRep-HeSiaGraph are publicly available at https://github.com/CBRC-lab/DrugRep-HeSiaGraph .

G-quadruplex DNA selective targeting for anticancer therapy: a computational study of a novel PtII monofunctional complex activated by adaptive binding

Targeting of G-quadruplex (G-Q) nucleic acids, which are helical four-stranded structures formed from guanine-rich nucleic acid sequences, has emerged in recent years as an appealing opportunity for drug intervention in anticancer therapy. Small-molecule drugs can stabilize quadruplex structures, promoting selective downregulation of gene expression and telomerase inhibition and also activating DNA damage responses. Thus, rational design of small molecular ligands able to selectively interact with and stabilize G-Q structures is a promising strategy for developing potent anti-cancer drugs with selective toxicity towards cancer cells over normal ones. Here, the outcomes of a thorough computational investigation of a recently synthesized monofunctional PtII complex (Pt1), whose selectivity for G-Q is activated by what is called adaptive binding, are reported. Quantum mechanics and molecular dynamics calculations have been employed for studying the classical key steps of the mechanism of action of PtII complexes, the conversion of the non-charged and non-planar Pt1 complex into a planar and charged PtII (Pt2) complex able to play the role of a G-Q binder and, finally, the interaction of Pt2 with G-Q. The information obtained from such an investigation allows us to rationalize the behavior of the novel PtII complex proposed to be activated by adaptive binding toward selective interaction with G-Q or similar molecules and can be exploited for designing ligands with more effective recognition ability toward G-quadruplex DNA.

MuDCoD: multi-subject community detection in personalized dynamic gene networks from single-cell RNA sequencing

MOTIVATION

With the wide availability of single-cell RNA-seq (scRNA-seq) technology, population-scale scRNA-seq datasets across multiple individuals and time points are emerging. While the initial investigations of these datasets tend to focus on standard analysis of clustering and differential expression, leveraging the power of scRNA-seq data at the personalized dynamic gene co-expression network level has the potential to unlock subject and/or time-specific network-level variation, which is critical for understanding phenotypic differences. Community detection from co-expression networks of multiple time points or conditions has been well-studied; however, none of the existing settings included networks from multiple subjects and multiple time points simultaneously. To address this, we develop Multi-subject Dynamic Community Detection (MuDCoD) for multi-subject community detection in personalized dynamic gene networks from scRNA-seq. MuDCoD builds on the spectral clustering framework and promotes information sharing among the networks of the subjects as well as networks at different time points. It clusters genes in the personalized dynamic gene networks and reveals gene communities that are variable or shared not only across time but also among subjects.

RESULTS

Evaluation and benchmarking of MuDCoD against existing approaches reveal that MuDCoD effectively leverages apparent shared signals among networks of the subjects at individual time points, and performs robustly when there is no or little information sharing among the networks. Applications to population-scale scRNA-seq datasets of human-induced pluripotent stem cells during dopaminergic neuron differentiation and CD4+ T cell activation indicate that MuDCoD enables robust inference for identifying time-varying personalized gene modules. Our results illustrate how personalized dynamic community detection can aid in the exploration of subject-specific biological processes that vary across time.

AVAILABILITY AND IMPLEMENTATION

MuDCoD is publicly available at https://github.com/bo1929/MuDCoD as a Python package. Implementation includes simulation and real-data experiments together with extensive documentation.

Chemoproteomic Approaches to Studying RNA Modification-Associated Proteins

The function of cellular RNA is modulated by a host of post-transcriptional chemical modifications installed by dedicated RNA-modifying enzymes. RNA modifications are widespread in biology, occurring in all kingdoms of life and in all classes of RNA molecules. They regulate RNA structure, folding, and protein-RNA interactions, and have important roles in fundamental gene expression processes involving mRNA, tRNA, rRNA, and other types of RNA species. Our understanding of RNA modifications has advanced considerably; however, there are still many outstanding questions regarding the distribution of modifications across all RNA transcripts and their biological function. One of the major challenges in the study of RNA modifications is the lack of sequencing methods for the transcriptome-wide mapping of different RNA-modification structures. Furthermore, we lack general strategies to characterize RNA-modifying enzymes and RNA-modification reader proteins. Therefore, there is a need for new approaches to enable integrated studies of RNA-modification chemistry and biology.In this Account, we describe our development and application of chemoproteomic strategies for the study of RNA-modification-associated proteins. We present two orthogonal methods based on nucleoside and oligonucleotide chemical probes: 1) RNA-mediated activity-based protein profiling (RNABPP), a metabolic labeling strategy based on reactive modified nucleoside probes to profile RNA-modifying enzymes in cells and 2) photo-cross-linkable diazirine-containing synthetic oligonucleotide probes for identifying RNA-modification reader proteins.We use RNABPP with C5-modified cytidine and uridine nucleosides to capture diverse RNA-pyrimidine-modifying enzymes including methyltransferases, dihydrouridine synthases, and RNA dioxygenase enzymes. Metabolic labeling facilitates the mechanism-based cross-linking of RNA-modifying enzymes with their native RNA substrates in cells. Covalent RNA-protein complexes are then isolated by denaturing oligo(dT) pulldown, and cross-linked proteins are identified by quantitative proteomics. Once suitable modified nucleosides have been identified as mechanism-based proteomic probes, they can be further deployed in transcriptome-wide sequencing experiments to profile the substrates of RNA-modifying enzymes at nucleotide resolution. Using 5-fluorouridine-mediated RNA-protein cross-linking and sequencing, we analyzed the substrates of human dihydrouridine synthase DUS3L. 5-Ethynylcytidine-mediated cross-linking enabled the investigation of ALKBH1 substrates. We also characterized the functions of these RNA-modifying enzymes in human cells by using genetic knockouts and protein translation reporters.We profiled RNA readers for N6-methyladenosine (m6A) and N1-methyladenosine (m1A) using a comparative proteomic workflow based on diazirine-containing modified oligonucleotide probes. Our approach enables quantitative proteome-wide analysis of the preference of RNA-binding proteins for modified nucleotides across a range of affinities. Interestingly, we found that YTH-domain proteins YTHDF1/2 can bind to both m6A and m1A to mediate transcript destabilization. Furthermore, m6A also inhibits stress granule proteins from binding to RNA.Taken together, we demonstrate the application of chemical probing strategies, together with proteomic and transcriptomic workflows, to reveal new insights into the biological roles of RNA modifications and their associated proteins.

Conformational Morphing by a DNA Analogue Featuring 7-Deazapurines and 5-Halogenpyrimidines and the Origins of Adenine-Tract Geometry

Several efforts are currently directed at the creation and cellular implementation of alternative genetic systems composed of pairing components that are orthogonal to the natural dA/dT and dG/dC base pairs. In an alternative approach, Watson-Crick-type pairing is conserved, but one or all of the four letters of the A, C, G, and T alphabet are substituted by modified components. Thus, all four nucleobases were altered to create halogenated deazanucleic acid (DZA): dA was replaced by 7-deaza-2'-deoxyadenosine (dzA), dG by 7-deaza-2'-deoxyguanosine (dzG), dC by 5-fluoro-2'-deoxycytidine (FdC), and dT by 5-chloro-2'-deoxyuridine (CldU). This base-pairing system was previously shown to retain function in Escherichia coli. Here, we analyze the stability, hydration, structure, and dynamics of a DZA Dickerson-Drew Dodecamer (DDD) of sequence 5'-FdC-dzG-FdC-dzG-dzA-dzA-CldU-CldU-FdC-dzG-FdC-dzG-3'. Contrary to similar stabilities of DDD and DZA-DDD, osmotic stressing revealed a dramatic loss of hydration for the DZA-DDD relative to that for the DDD. The parent DDD 5'-d(CGCGAATTCGCG)-3' features an A-tract, a run of adenosines uninterrupted by a TpA step, and exhibits a hallmark narrow minor groove. Crystal structures─in the presence of RNase H─and MD simulations show increased conformational plasticity ("morphing") of DZA-DDD relative to that of the DDD. The narrow dzA-tract minor groove in one structure widens to resemble that in canonical B-DNA in a second structure. These changes reflect an indirect consequence of altered DZA major groove electrostatics (less negatively polarized compared to that in DNA) and hydration (reduced compared to that in DNA). Therefore, chemical modifications outside the minor groove that lead to collapse of major groove electrostatics and hydration can modulate A-tract geometry.

Underlying Mechanisms for the Sex- and Chemical-Specific Hepatotoxicity of Perfluoroalkyl Phosphinic Acids in Common Carp (Cyprinus carpio)

The hepatotoxicities of perfluoroalkyl and polyfluoroalkyl substances (PFASs) have been extensively investigated, while little is known about the sex-specific differences. In this study, common carp were exposed to the emerging perfluoroalkyl phosphinic acids (6:6 and 8:8 PFPiAs) for 14 days to disclose sex-specific hepatotoxicity. Apparent hepatotoxicity, including cell necrosis, apoptosis, and steatosis, was observed in both male and female carp liver. The observed hepatocyte steatosis was predominantly attributed to the dysregulation of hepatic lipid metabolism but was based on sex-specific mechanisms. It was manifested as inhibited oxidative decomposition of fatty acids (FAs) in the female liver, whereas it enhanced the uptake of FAs into the male liver, both of which led to excessive lipid accumulation. Untargeted lipidomics validated that the metabolism pathways of FA, sphingolipid, glycerolipid, and glycerophospholipid were disrupted by both compounds, leading to the generation of reactive oxygen species and oxidative stress. The oxidative stress further evolved into inflammation, manifested as promoted expression of proinflammatory cytokines and repressed expression of anti-inflammatory cytokines. Consistently, all of the changes were more noticeable in male carp, suggesting that male fish were more susceptible to PFPiA disruption. 8:8 PFPiA was less accumulated but caused stronger hepatotoxicity than 6:6 PFPiA, possibly because of the stronger binding capacity of 8:8 PFPiA to nuclear transcription factors mediating lipid metabolism and inflammation. The findings of this study highlight the significance of sex- and chemical-dependent bioaccumulation and the toxicity of PFASs in organisms.

Bisulfite and Nanopore Sequencing for Pseudouridine in RNA

Nucleophilic addition of bisulfite to pyrimidine bases has been known for a half century, and the reaction has been in use for at least a quarter of a century for identifying 5-methylcytidine in DNA. This account focuses on the chemistry of bisulfite with pseudouridine, an isomer of the RNA nucleoside uridine in which the uracil base is connected to C1' of ribose via C5 instead of N1. Pseudouridine, Ψ, is the most common nucleotide modification found in cellular RNA overall, in part due to its abundance in rRNAs and tRNAs. It has a stabilizing influence on RNA structure because N1 is now available for additional hydrogen bonding and because the heterocycle is slightly better at π stacking. The isomerization of U to Ψ in RNA strands is catalyzed by 13 different enzymes in humans and 11 in E. coli; some of these enzymes are implicated in disease states which is testament to the biological importance of pseudouridine in cells. Recently, pseudouridine came into the limelight as the key modification that, after N1 methylation, enables mRNA vaccines to be delivered efficiently into human tissue with minimal generation of a deleterious immunogenic response. Here we describe the bisulfite reaction with pseudouridine which gives rise to a chemical sequencing method to map the modified base in the epitranscriptome. Unlike the reaction with cytidine, the addition of bisulfite to Ψ leads irreversibly to form an adduct that is bypassed during cDNA synthesis by reverse transcriptases yielding a characteristic deletion signature. Although there were hints to the structure of the bisulfite adduct(s) 30 to 50 years ago, it took modern spectroscopic and computational methods to solve the mystery. Raman spectroscopy along with extensive NMR, ECD, and computational work led to the assignment of the major product as the (R) diastereomer of an oxygen adduct at C1' of a ring-opened pseudouridine. Mechanistically, this arose from a succession of conjugate addition, E2 elimination, and a [2,3] sigmatropic rearrangement, all of which are stereodefined reactions. A minor reaction with excess bisulfite led to the (S) isomer of a S-adducted SO3- group. Understanding structure and mechanism aided the design of a Ψ-specific sequencing reaction and guided attempts to improve the utility and specificity of the method. Separately, we have been investigating the use of nanopore direct RNA sequencing, a single-molecule method that directly analyzes RNA strands isolated from cells after end-ligation of adaptor sequences. By combining the electrical current and base-calling data from the nanopore with dwell-time analysis from the helicase employed to deliver RNA to the nanopore, we were able to map Ψ sites in nearly all sequence contexts. This analysis was employed to find Ψ residues in the SARS-CoV-2 vRNA, to analyze the sequence context effects of mRNA vaccine synthesis via in vitro transcription, and to evaluate the impact of stress on chemical modifications in the E. coli ribosome. Most recently, we found that bisulfite treatment of RNA leading to Ψ adducts could modulate the nanopore signal to help in mapping modifications of low occupancy.

Three E2F target-related genes signature for predicting prognosis, immune features, and drug sensitivity in hepatocellular carcinoma

Background: Hepatocellular carcinoma (HCC) is extremely malignant and difficult to treat. The adenoviral early region 2 binding factors (E2Fs) target pathway is thought to have a major role in tumor growth. This study aimed to identify a predictive E2F target signature and facilitate individualized treatment for HCC patients. Methods: We constructed an E2F target-related gene profile using univariate COX and LASSO regression models and proved its predictive efficacy in external cohorts. Furthermore, we characterized the role of the E2F target pathway in pathway enrichment, immune cell infiltration, and drug sensitivity of HCC. Results: Lasso Cox regression created an E2F target-related gene signature of GHR, TRIP13, and CDCA8. HCC patients with high risk were correlated with shorter survival time, immune evasion, tumor stem cell characteristics and high sensitivity to Tipifarnib and Camptothecin drugs. Conclusion: Hepatocellular carcinoma prognosis was predicted by an E2F target signature. This finding establishes the theoretical usefulness of the E2F target route in customized identification and treatment for future research.

Revealing the diagnostic value and immune infiltration of senescence-related genes in endometriosis: a combined single-cell and machine learning analysis

Introduction: Endometriosis is a prevalent and recurrent medical condition associated with symptoms such as pelvic discomfort, dysmenorrhea, and reproductive challenges. Furthermore, it has the potential to progress into a malignant state, significantly impacting the quality of life for affected individuals. Despite its significance, there is currently a lack of precise and non-invasive diagnostic techniques for this condition. Methods: In this study, we leveraged microarray datasets and employed a multifaceted approach. We conducted differential gene analysis, implemented weighted gene co-expression network analysis (WGCNA), and utilized machine learning algorithms, including random forest, support vector machine, and LASSO analysis, to comprehensively explore senescence-related genes (SRGs) associated with endometriosis. Discussion: Our comprehensive analysis, which also encompassed profiling of immune cell infiltration and single-cell analysis, highlights the therapeutic potential of this gene assemblage as promising targets for alleviating endometriosis. Furthermore, the integration of these biomarkers into diagnostic protocols promises to enhance diagnostic precision, offering a more effective diagnostic journey for future endometriosis patients in clinical settings. Results: Our meticulous investigation led to the identification of a cluster of genes, namely BAK1, LMNA, and FLT1, which emerged as potential discerning biomarkers for endometriosis. These biomarkers were subsequently utilized to construct an artificial neural network classifier model and were graphically represented in the form of a Nomogram.

PET and Optical Imaging of Caveolin-1 in Gastric Tumors

Previous studies have suggested tumoral caveolin-1 (CAV1) as a predictive biomarker for the response to anti-HER2 antibody drug therapies in gastric tumors. In this study, radiolabeled and fluorescently labeled anti-CAV1 antibodies were developed and tested as an immunoPET or optical imaging agent to detect CAV1 in HER2-positive/CAV1-high NCIN87 gastric tumors. The expression of CAV1 receptors in NCIN87 gastric tumors and nontumor murine organs was determined by Western blot. Binding assays were performed to validate the anti-CAV1 antibody specificity for CAV1-expressing NCIN87 cancer cells. Subcutaneous and orthotopic NCIN87 xenografts were used for PET imaging and ex vivo biodistribution of the radioimmunoconjugate. Additional HER2-PET and CAV1-optical imaging was also performed to determine CAV1 in the HER2-positive tumors. 89Zr-labeled anti-CAV1 antibody was able to bind to CAV1-expressing NCIN87 cells with a Bmax value of 2.7 × 103 CAV1 receptors/cell in vitro. ImmunoPET images demonstrated the localization of the antibody in subcutaneous NCIN87 xenografts. In the orthotopic model, CAV1 expression was also observed by optical imaging in the HER2-positive tumors previously imaged with HER2-PET. Ex vivo biodistribution analysis further confirmed these imaging results. The preclinical data from this study demonstrate the potential of using CAV1-PET and optical imaging for detecting gastric tumors.

Incidental finding at methylation-specific multiplex ligation-dependent probe amplification (MS-MLPA): how to proceed?

Introduction: Since the advent of new generation sequencing, professionals are aware of the possibility of obtaining findings unrelated to the pathology under study. However, this possibility is usually forgotten in the case of studies aimed at a single gene or region. We report a case of a 16-month-old girl with clinical suspicion of Silver-Russell syndrome (SRS). Methods: Following the international SRS consensus, methylation alterations and copy number variations (CNVs) at 11p15 region and maternal uniparental disomy of chromosome 7 were analysed and discarded by MS-MLPA. Results: Unexpectedly, the 11p15 region MS-MLPA showed a decrease in the signal of a copy number reference probe. Deletions affecting a single probe are inconclusive. So, we faced the ethical dilemma of whether it was appropriate to confirm this alteration with independent techniques and to offer a diagnostic possibility that was in no way related to clinical suspicion. Fortunately, in this particular case, the informed consent had not been specific to a particular pathology but to any disorder associated with growth failure. Performed alternative studies allowed the final diagnosis of 22q deletion syndrome. Conclusion: We demonstrate the importance of informing patients about the possibility of obtaining incidental findings in genetic techniques (not only in next generation sequencing) during pre-test genetic counselling consultations. In addition, we highlight the relevance of including in the informed consent the option of knowing these unexpected incidental findings as in some cases, this will help to elucidate the definitive diagnosis and provide the correct follow-up and treatment.

Characterization of the complete mitochondrial genome and phylogenetic analysis of bean thrips Megalurothrips usitatus (Bagnall, 1913) (Thysanoptera: thripidae)

Megalurothrips usitatus (Bagnall, 1913) (Thysanoptera: Thripidae) is a widely distributed pest in Asia that primarily affects the production of snap beans and cowpea. The complete mitochondrial genome of Megalurothrips usitatus has been sequenced and annotated in this study, which is 17,209 bp long and contains 13 protein-coding genes (PCGs), two rRNAs, and 22 tRNA genes. Most of the protein-coding genes (PCGs) start with ATG except ND4 using TTG. Meanwhile, eight PCGs stop with TAA, four PCGs have an incomplete stop codon, and the gene Cytb ends with TAG. Phylogenetic analysis showed that M. usitatus is closely related to Frankliniella intonsa and F. occidentalis, providing a basis for the study of the mitochondrial evolution of Thripinae.

Mechanisms of metabolic adaptation in the duckweed Lemna gibba: an integrated metabolic, transcriptomic and flux analysis

BACKGROUND

Duckweeds are small, rapidly growing aquatic flowering plants. Due to their ability for biomass production at high rates they represent promising candidates for biofuel feedstocks. Duckweeds are also excellent model organisms because they can be maintained in well-defined liquid media, usually reproduce asexually, and because genomic resources are becoming increasingly available. To demonstrate the utility of duckweed for integrated metabolic studies, we examined the metabolic adaptation of growing Lemna gibba cultures to different nutritional conditions.

RESULTS

To establish a framework for quantitative metabolic research in duckweeds we derived a central carbon metabolism network model of Lemna gibba based on its draft genome. Lemna gibba fronds were grown with nitrate or glutamine as nitrogen source. The two conditions were compared by quantification of growth kinetics, metabolite levels, transcript abundance, as well as by 13C-metabolic flux analysis. While growing with glutamine, the fronds grew 1.4 times faster and accumulated more protein and less cell wall components compared to plants grown on nitrate. Characterization of photomixotrophic growth by 13C-metabolic flux analysis showed that, under both metabolic growth conditions, the Calvin-Benson-Bassham cycle and the oxidative pentose-phosphate pathway are highly active, creating a futile cycle with net ATP consumption. Depending on the nitrogen source, substantial reorganization of fluxes around the tricarboxylic acid cycle took place, leading to differential formation of the biosynthetic precursors of the Asp and Gln families of proteinogenic amino acids. Despite the substantial reorganization of fluxes around the tricarboxylic acid cycle, flux changes could largely not be associated with changes in transcripts.

CONCLUSIONS

Through integrated analysis of growth rate, biomass composition, metabolite levels, and metabolic flux, we show that Lemna gibba is an excellent system for quantitative metabolic studies in plants. Our study showed that Lemna gibba adjusts to different nitrogen sources by reorganizing central metabolism. The observed disconnect between gene expression regulation and metabolism underscores the importance of metabolic flux analysis as a tool in such studies.

PLK1 maintains DNA methylation and cell viability by regulating phosphorylation-dependent UHRF1 protein stability

PLK1 is a key serine/threonine kinase as well as a master mitotic regulator, but it has never been reported that PLK1 regulates DNA methylation. In the present study, we for the first time found that PLK1 inhibition disrupted global DNA methylation and elevated the expression level of tumor suppressor genes. Mechanistically, we found that PLK1 interacts UHRF1 protein to induce its phosphorylation at serine 265. Phosphorylation is required for the maintenance of UHRF1 protein stability by recruiting a deubiquitinase USP7. Conversely, PLK1 inhibition decreases UHRF1 protein interaction with USP7 and activates the ubiquitin-proteasome pathway, thereby accelerating UHRF1 protein degradation. UHRF1 degradation decreases the recruitment of DNMT1 to chromatin, and decreases the level of genome-wide DNA methylation, thereby elevating the expression of tumor suppressor genes and decreasing cell viability. We here presented the first report on the novel role of PLK1 in DNA methylation maintenance through UHRF1-DNMT1 pathway, and revealed a novel anticancer mechanism of PLK1 inhibitors.

Balancing biomass reaction stoichiometry and measured fluxes in flux balance analysis

MOTIVATION

Flux balance analysis (FBA) is widely recognized as an important method for studying metabolic networks. When incorporating flux measurements of certain reactions into an FBA problem, it is possible that the underlying linear program may become infeasible, e.g. due to measurement or modeling inaccuracies. Furthermore, while the biomass reaction is of central importance in FBA models, its stoichiometry is often a rough estimate and a source of high uncertainty.

RESULTS

In this work, we present a method that allows modifications to the biomass reaction stoichiometry as a means to (i) render the FBA problem feasible and (ii) improve the accuracy of the model by corrections in the biomass composition. Optionally, the adjustment of the biomass composition can be used in conjunction with a previously introduced approach for balancing inconsistent fluxes to obtain a feasible FBA system. We demonstrate the value of our approach by analyzing realistic flux measurements of E.coli. In particular, we find that the growth-associated maintenance (GAM) demand of ATP, which is typically integrated with the biomass reaction, is likely overestimated in recent genome-scale models, at least for certain growth conditions. In light of these findings, we discuss issues related to the determination and inclusion of GAM values in constraint-based models. Overall, our method can uncover potential errors and suggest adjustments in the assumed biomass composition in FBA models based on inconsistencies between the model and measured fluxes.

AVAILABILITY AND IMPLEMENTATION

The developed method has been implemented in our software tool CNApy available from https://github.com/cnapy-org/CNApy.

Similarities and differences in gene expression profiles of BRCA1 methylated and mutated epithelial ovarian cancers

Introduction

BRCA1 methylated (BRCA1met) epithelial ovarian cancer (EOC) is a recently defined and not well-investigated subset of neoplasms. To date, no studies have focused on the transcriptional profiles of BRCA1met cases, and, as a matter of fact, we still do not know if this subset of EOCs is similar, and to what extent, to BRCA1 mutated (BRCA1mut) cases.

Methods

We compared a group of 17 BRCA1met cases against 10 BRCA1mut cases using a subset of carefully selected 17 BRCAwt EOCs as a control group.

Results

First, BRCA1met cases showed a downregulation of the relative transcript, while this association was not observed for BRCA1mut EOCs. The BRCA1met group exhibited a general upregulation of homologous recombination (HR)-related genes, as well as BRCA1mut. Overall, BRCA1met had a different gene expression profile, characterized by diffuse downregulation, whereas BRCA1mut showed a general upregulation (p < 0.0001). Both BRCA1-defective groups showed a slightly activated immune response mediated by interferon (IFN) gamma pathways.

Discussion

In conclusion, even if the expression profile of many genes related to DNA damage and repair system is shared between BRCA1mut and BRCA1met EOCs supporting that BRCA1met EOCs may benefit from PARPi therapies, our data demonstrate that BRCA1mut and BRCA1met EOCs show different expression profiles, suggesting a different mechanism of carcinogenesis that can be reflected in different responses to therapies and disease recovery.

SHIP: identifying antimicrobial resistance gene transfer between plasmids

MOTIVATION

Plasmids are carriers for antimicrobial resistance (AMR) genes and can exchange genetic material with other structures, contributing to the spread of AMR. There is no reliable approach to identify the transfer of AMR genes across plasmids. This is mainly due to the absence of a method to assess the phylogenetic distance of plasmids, as they show large DNA sequence variability. Identifying and quantifying such transfer can provide novel insight into the role of small mobile elements and resistant plasmid regions in the spread of AMR.

RESULTS

We developed SHIP, a novel method to quantify plasmid similarity based on the dynamics of plasmid evolution. This allowed us to find conserved fragments containing AMR genes in structurally different and phylogenetically distant plasmids, which is evidence for lateral transfer. Our results show that regions carrying AMR genes are highly mobilizable between plasmids through transposons, integrons, and recombination events, and contribute to the spread of AMR. Identified transferred fragments include a multi-resistant complex class 1 integron in Escherichia coli and Klebsiella pneumoniae, and a region encoding tetracycline resistance transferred through recombination in Enterococcus faecalis.

AVAILABILITY AND IMPLEMENTATION

The code developed in this work is available at https://github.com/AbeelLab/plasmidHGT.

Antiprotozoal activity of auranofin on Trypanosoma cruzi, Leishmania tropica and Toxoplasma gondii: in vitro and ex vivo study

BACKGROUND

Three obligate intracellular protozoan parasite species, which are responsible for significant morbidity and mortality and settle in macrophage cells, affect more than one-half of the world's population, namely, Trypanosoma cruzi, Leishmania tropica and Toxoplasma gondii, which are causative agents of Chagas disease, leishmaniasis and toxoplasmosis, respectively. In the current study, it was aimed to investigate the in vitro and ex vivo antiprotozoal activity of auranofin on T. cruzi, L. tropica and T. gondii.

METHODS

The in vitro drug efficacy (IC50) of auranofin was investigated by haemocytometry and the CellTiter-Glo assay methods and the ex vivo drug efficacy (IC50) by light microscopic examination of Giemsa-stained slides. Also, the cytotoxic activity (CC50) of auranofin was examined by the CellTiter-Glo assay. The selectivity index (SI) was calculated for auranofin.

RESULTS

According to IC50, CC50 and SI data, auranofin did not exhibit cytotoxic activity on Vero cells, but exhibited antiprotozoal activity on epimastigotes and intracellular amastigotes of T. cruzi, promastigotes and intracellular amastigotes of L. tropica and intracellular tachyzoites of T. gondii (p<0.05).

CONCLUSIONS

The detection antiprotozoal activity of auranofin on T. cruzi, L. tropica and T. gondii according to the IC50, CC50 and SI values is considered an important and promising development. This is significant because auranofin may be an effective alternative treatment for Chagas disease, leishmaniasis and toxoplasmosis in the future.

Heterologous Production in the Synechocystis Chassis Suggests the Biosynthetic Pathway of Astaxanthin in Cyanobacteria

Astaxanthin is a carotenoid species with the highest antioxidant capability. Its natural resource is very rare. The biosynthesis of astaxanthin from β-carotene includes a hydroxylation step and a ketolation step, for which the corresponding enzymes have been characterized in a few species. However, the sequence of these two reactions is unclear, and may vary with different organisms. In this study, we aimed to elucidate this sequence in Synechocystis, which is an ideal cyanobacterial synthetic biology chassis. We first silenced the endogenous carotene oxygenase gene SyneCrtO to avoid its possible interference in the carotenoid metabolic network. We then introduced the β-carotene ketolase gene from Haematococcus pluvialis (HpBKT) and the CrtZ-type carotene β-hydroxylase gene from Pantoea agglomerans (PaCrtZ) to this δCrtO strain. Our pigment analysis demonstrated that both the endogenous CrtR-type carotene hydroxylase SyneCrtR and HpBKT have the preference to use β-carotene as their substrate for hydroxylation and ketolation reactions to produce zeaxanthin and canthaxanthin, respectively. However, the endogenous SyneCrtR is not able to further catalyze the 3,3'-hydroxylation of canthaxanthin to generate astaxanthin. From our results, a higher accumulation of canthaxanthin and a much lower level of astaxanthin, as confirmed using liquid chromatography-tandem mass spectrometry analysis, were detected in our transgenic BKT+/CrtZ+/δCrtO cells. Therefore, we proposed that the bottleneck for the heterologous production of astaxanthin in Synechocystis might exist at the hydroxylation step, which requires a comprehensive screening or genetic engineering for the corresponding carotene hydroxylase to enable the industrial production of astaxanthin.

Development of an Integrated Sample Amplification Control for Salivary Point-of-Care Pathogen Testing

Background

The COVID-19 pandemic has led to a rise in point-of-care (POC) and home-based tests, but concerns over usability, accuracy, and effectiveness have arisen. The incorporation of internal amplification controls (IACs), essential control for translational POC diagnostics, could mitigate false-negative and false-positive results due to sample matrix interference or inhibition. Although emerging POC nucleic acid amplification tests (NAATs) for detecting SARS-CoV-2 show impressive analytical sensitivity in the lab, the assessment of clinical accuracy with IACs is often overlooked. In some cases, the IACs were run spatially, complicating assay workflow. Therefore, the multiplex assay for pathogen and IAC is needed.

Results

We developed a one-pot duplex reverse transcriptase loop-mediated isothermal amplification (RT-LAMP) assay for saliva samples, a non-invasive and simple collected specimen for POC NAATs. The ORF1ab gene of SARS-CoV-2 was used as a target and a human 18S ribosomal RNA in human saliva was employed as an IAC to ensure clinical reliability of the RT-LAMP assay. The optimized assay could detect SARS-CoV-2 viral particles down to 100 copies/μL of saliva within 30 minutes without RNA extraction. The duplex RT-LAMP for SARS-CoV-2 and IAC is successfully amplified in the same reaction without cross-reactivity. The valid results were easily visualized in triple-line lateral flow immunoassay, in which two lines (flow control and IAC lines) represent valid negative results and three lines (flow control, IAC, and test line) represent valid positive results. This duplex assay demonstrated a clinical sensitivity of 95%, specificity of 100%, and accuracy of 96% in 30 clinical saliva samples.

Significance

IACs play a crucial role in ensuring user confidence with respect to the accuracy and reliability of at-home and POC molecular diagnostics. We demonstrated the multiplex capability of SARS-COV-2 and human18S ribosomal RNA RT-LAMP without complicating assay design. This generic platform can be extended in a similar manner to include human18S ribosomal RNA IACs into different clinical sample matrices.

stVAE deconvolves cell-type composition in large-scale cellular resolution spatial transcriptomics

MOTIVATION

Recent rapid developments in spatial transcriptomic techniques at cellular resolution have gained increasing attention. However, the unique characteristics of large-scale cellular resolution spatial transcriptomic datasets, such as the limited number of transcripts captured per spot and the vast number of spots, pose significant challenges to current cell-type deconvolution methods.

RESULTS

In this study, we introduce stVAE, a method based on the variational autoencoder framework to deconvolve the cell-type composition of cellular resolution spatial transcriptomic datasets. To assess the performance of stVAE, we apply it to five datasets across three different biological tissues. In the Stereo-seq and Slide-seqV2 datasets of the mouse brain, stVAE accurately reconstructs the laminar structure of the pyramidal cell layers in the cortex, which are mainly organized by the subtypes of telencephalon projecting excitatory neurons. In the Stereo-seq dataset of the E12.5 mouse embryo, stVAE resolves the complex spatial patterns of osteoblast subtypes, which are supported by their marker genes. In Stereo-seq and Pixel-seq datasets of the mouse olfactory bulb, stVAE accurately delineates the spatial distributions of known cell types. In summary, stVAE can accurately identify spatial patterns of cell types and their relative proportions across spots for cellular resolution spatial transcriptomic data. It is instrumental in understanding the heterogeneity of cell populations and their interactions within tissues.

AVAILABILITY AND IMPLEMENTATION

stVAE is available in GitHub (https://github.com/lichen2018/stVAE) and Figshare (https://figshare.com/articles/software/stVAE/23254538).

Orderly Aggregated Catalytic Hairpin Assembly for Synchronous Ultrasensitive Detecting and High-Efficiency Co-Localization Imaging of Dual-miRNAs in Living Cells

In this work, the orderly aggregated catalytic hairpin assembly (OA-CHA) was developed for synchronous ultrasensitive detection and high-efficiency colocalization imaging of dual-miRNAs by a carefully designed tetrahedral conjugated ladder DNA structure (TCLDS). Exactly, two diverse hairpin probes were fixed on tetrahedron conjugated DNA nanowires to form the TCLDS without fluorescence response, which triggered OA-CHA in the aid of output DNA 1 and output DNA 2 produced by targets miRNA-217 and miRNA-196a cycle to generate TCLDS with remarkable fluorescence response. Impressively, compared with the traditional CHA strategy, OA-CHA avoided the fluorescence group and quenching group from approaching again because of the spatial confinement effect to significantly enhance the fluorescence signal, resulting in the simultaneous ultrasensitive detection of dual-miRNAs with detection limits of 21 and 32 fM for miRNA-217 and miRNA-196a, respectively. Meanwhile, the TCLDS with lower diffusivity could achieve accurate localization imaging for reflecting the spatial distribution of dual-miRNAs in living cells. The strategy based on OA-CHA provided a flexible and programmable nucleic amplification method for the synchronous ultrasensitive detection and precise imaging of multiple biomarkers and had potential in disease diagnostics..

Metabolic Engineering of Escherichia coli for Efficient Production of Pseudouridine

Pseudouridine-incorporated mRNA vaccines can enhance protein expression and reduce immunogenicity, leading to a high demand for pseudouridine to be used in mRNA drug production. To achieve the low-cost production of pseudouridine, Escherichia coli was systematically modified to utilize inexpensive raw materials to efficiently produce pseudouridine. First, in the pyrimidine biosynthesis pathway, genes related to the precursor competing pathway and the negative regulator were deleted, which increased pseudouridine production. Second, two critical genes, pseudouridine-5'-phosphate glycosidase (psuG) and phosphatase genes from different bacteria, were screened and employed in various genetic constructs, and the pseudouridine yield of the optical strain increased to 599 mg/L. The accumulation of pseudouridine was further increased by the deletion of pseudouridine catabolism-related genes. Ultimately, the pseudouridine titer in a 5 L bioreactor reached 7.9 g/L, and the yield of pseudouridine on glucose was 0.15 g/g. Overall, a cell factory producing pseudouridine was successfully constructed and showed potential for industrial production.

Effects of fermented herbal tea residue on meat quality, rumen fermentation parameters and microbes of black goats

Herbal tea residue (HTR) is generally considered to be a reusable resource which has still retains considerable proportion of nutrients and active substances. This study aimed to investigate the effects of substitution of whole corn silage with fermented herbal tea residue (FHTR) on meat quality, serum indices, rumen fermentation, and microbes in Chuanzhong black goats. Twenty-two female Chuanzhong black goats (4 months old) with similar weight (9.55 ± 0.95 kg) were selected and randomly divided into two groups. FHTR was used to replace 0% (CON group) and 30% (FHTR group) of whole corn silage in the diets and fed as a total mixed ration (TMR) for Chuanzhong black goats. The adaptation feeding period was 7 days, and the experimental period was 35 days. Results illustrated that the FHTR group had higher value of a* and concentrations of DM and CP and lower rate of water loss (P < 0.05) than the CON group. For the serum indices, goats fed with 30% FHTR had higher (P < 0.05) concentration of CR on day 35. For rumen fermentation, the pH and ratio of acetic acid/propionic acid (AA/PA) in the FHTR group were significantly lower than those in the CON group (P < 0.05). In addition, we studied the goats's rumen microbial community composition and found that the dominant phyla were Firmicutes, Bacteroidetes,and Tenericutes; and the dominant genera were Quinella, Candidatus_Saccharimonas, and Saccharofermentans. There was a significant difference in the beta diversity of the rumen microbiota between groups (P < 0.05). To sum up, the addition of FHTR can affect the meat quality, serum indices, improved rumen fermentation by adjusted the diversity and function of the rumen microbiota.

Synthesis of Peptidyl-tRNA Mimics for Structural Biology Applications

Protein biosynthesis is a central process in all living cells that is catalyzed by a complex molecular machine─the ribosome. This process is termed translation because the language of nucleotides in mRNAs is translated into the language of amino acids in proteins. Transfer RNA (tRNA) molecules charged with amino acids serve as adaptors and recognize codons of mRNA in the decoding center while simultaneously the individual amino acids are assembled into a peptide chain in the peptidyl transferase center (PTC). As the nascent peptide emerges from the ribosome, it is threaded through a long tunnel referred to as a nascent peptide exit tunnel (NPET). The PTC and NPET are the sites targeted by many antibiotics and are thus of tremendous importance from a biomedical perspective and for drug development in the pharmaceutical industry.Researchers have achieved much progress in characterizing ribosomal translation at the molecular level; an impressive number of high-resolution structures of different functional and inhibited states of the ribosome are now available. These structures have significantly contributed to our understanding of how the ribosome interacts with its key substrates, namely, mRNA, tRNAs, and translation factors. In contrast, much less is known about the mechanisms of how small molecules, especially antibiotics, affect ribosomal protein synthesis. This mainly concerns the structural basis of small molecule-NPET interference with cotranslational protein folding and the regulation of protein synthesis. Growing biochemical evidence suggests that NPET plays an active role in the regulation of protein synthesis.Much-needed progress in this field is hampered by the fact that during the preparation of ribosome complexes for structural studies (i.e., X-ray crystallography, cryoelectron microscopy, and NMR spectroscopy) the aminoacyl- or peptidyl-tRNAs are unstable and become hydrolyzed. A solution to this problem is the application of hydrolysis-resistant mimics of aminoacyl- or peptidyl-tRNAs.In this Account, we present an overview of synthetic methods for the generation of peptidyl-tRNA analogs. Modular approaches have been developed that combine (i) RNA and peptide solid-phase synthesis on 3'-aminoacylamino-adenosine resins, (ii) native chemical ligations and Staudinger ligations, (iii) tailoring of tRNAs by the selective cleavage of natural native tRNAs with DNAzymes followed by reassembly with enzymatic ligation to synthetic peptidyl-RNA fragments, and (iv) enzymatic tailing and cysteine charging of the tRNA to obtain modified CCA termini of a tRNA that are chemically ligated to the peptide moiety of interest. With this arsenal of tools, in principle, any desired sequence of a stably linked peptidyl-tRNA mimic is accessible. To underline the significance of the synthetic conjugates, we briefly point to the most critical applications that have shed new light on the molecular mechanisms underlying the context-specific activity of ribosome-targeting antibiotics, ribosome-dependent incorporation of multiple consecutive proline residues, the incorporation of d-amino acids, and tRNA mischarging.Furthermore, we discuss new types of stably charged tRNA analogs, relying on triazole- and squarate (instead of amide)-linked conjugates. Those have pushed forward our mechanistic understanding of nonribosomal peptide synthesis, where aminoacyl-tRNA-dependent enzymes are critically involved in various cellular processes in primary and secondary metabolism and in bacterial cell wall synthesis.

Angelica polysaccharides relieve blood glucose levels in diabetic KKAy mice possibly by modulating gut microbiota: an integrated gut microbiota and metabolism analysis

BACKGROUND

Angelica polysaccharides (AP) have numerous benefits in relieving type 2 diabetes (T2D). However, the underlying mechanisms have yet to be fully understood. Recent many reports have suggested that altering gut microbiota can have adverse effects on the host metabolism and contribute to the development of T2D. Here, we successfully established the T2D model using the male KKAy mice with high-fat and high-sugar feed. Meanwhile, the male C57BL/6 mice were fed with a normal feed. T2D KKAy mice were fed either with or without AP supplementation. In each group, we measured the mice's fasting blood glucose, weight, and fasting serum insulin levels. We collected the cecum content of mice, the gut microbiota was analyzed by targeted full-length 16S rRNA metagenomic sequencing and metabolites were analyzed by untargeted-metabolomics.

RESULTS

We found AP effectively alleviated glycemic disorders of T2D KKAy mice, with the changes in gut microbiota composition and function. Many bacteria species and metabolites were markedly changed in T2D KKAy mice and reversed by AP. Additionally, 16 altered metabolic pathways affected by AP were figured out by combining metagenomic pathway enrichment analysis and metabolic pathway enrichment analysis. The key metabolites in 16 metabolic pathways were significantly associated with the gut microbial alteration. Together, our findings showed that AP supplementation could attenuate the diabetic phenotype. Significant gut microbiota and gut metabolite changes were observed in the T2D KKAy mice and AP intervention.

CONCLUSIONS

Administration of AP has been shown to improve the composition of intestinal microbiota in T2D KKAy mice, thus providing further evidence for the potential therapeutic application of AP in the treatment of T2D.

ELIXIR and Toxicology: a community in development

Toxicology has been an active research field for many decades, with academic, industrial and government involvement. Modern omics and computational approaches are changing the field, from merely disease-specific observational models into target-specific predictive models. Traditionally, toxicology has strong links with other fields such as biology, chemistry, pharmacology and medicine. With the rise of synthetic and new engineered materials, alongside ongoing prioritisation needs in chemical risk assessment for existing chemicals, early predictive evaluations are becoming of utmost importance to both scientific and regulatory purposes. ELIXIR is an intergovernmental organisation that brings together life science resources from across Europe. To coordinate the linkage of various life science efforts around modern predictive toxicology, the establishment of a new ELIXIR Community is seen as instrumental. In the past few years, joint efforts, building on incidental overlap, have been piloted in the context of ELIXIR. For example, the EU-ToxRisk, diXa, HeCaToS, transQST, and the nanotoxicology community have worked with the ELIXIR TeSS, Bioschemas, and Compute Platforms and activities. In 2018, a core group of interested parties wrote a proposal, outlining a sketch of what this new ELIXIR Toxicology Community would look like. A recent workshop (held September 30th to October 1st, 2020) extended this into an ELIXIR Toxicology roadmap and a shortlist of limited investment-high gain collaborations to give body to this new community. This Whitepaper outlines the results of these efforts and defines our vision of the ELIXIR Toxicology Community and how it complements other ELIXIR activities.

Genome-wide identification, comparative analysis and functional roles in flavonoid biosynthesis of cytochrome P450 superfamily in pear (Pyrus spp.)

BACKGROUND

The cytochrome P450 (CYP) superfamily is the largest enzyme metabolism family in plants identified to date, and it is involved in many biological processes, including secondary metabolite biosynthesis, hormone metabolism and stress resistance. However, the P450 gene superfamily has not been well studied in pear (Pyrus spp.).

RESULTS

Here, the comprehensive identification and a comparative analysis of P450 superfamily members were conducted in cultivated and wild pear genomes. In total, 338, 299 and 419 P450 genes were identified in Chinese white pear, European pear and the wild pear, respectively. Based on the phylogenetic analyses, pear P450 genes were divided into ten clans, comprising 48 families. The motif and gene structure analyses further supported this classification. The expansion of the pear P450 gene family was attributed to whole-genome and single-gene duplication events. Several P450 gene clusters were detected, which have resulted from tandem and proximal duplications. Purifying selection was the major force imposed on the long-term evolution of P450 genes. Gene dosage balance, subfunctionalization and neofunctionalization jointly drove the retention and functional diversification of P450 gene pairs. Based on the association analysis between transcriptome expression profiles and flavonoid content during fruit development, three candidate genes were identified as being closely associated with the flavonoid biosynthesis, and the expression of one gene was further verified using qRT-PCR and its function was validated through transient transformation in pear fruit.

CONCLUSIONS

The study results provide insights into the evolution and biological functions of P450 genes in pear.

Graph regularized non-negative matrix factorization with [Formula: see text] norm regularization terms for drug-target interactions prediction

BACKGROUND

Identifying drug-target interactions (DTIs) plays a key role in drug development. Traditional wet experiments to identify DTIs are costly and time consuming. Effective computational methods to predict DTIs are useful to speed up the process of drug discovery. A variety of non-negativity matrix factorization based methods are proposed to predict DTIs, but most of them overlooked the sparsity of feature matrices and the convergence of adopted matrix factorization algorithms, therefore their performances can be further improved.

RESULTS

In order to predict DTIs more accurately, we propose a novel method iPALM-DLMF. iPALM-DLMF models DTIs prediction as a problem of non-negative matrix factorization with graph dual regularization terms and [Formula: see text] norm regularization terms. The graph dual regularization terms are used to integrate the information from the drug similarity matrix and the target similarity matrix, and [Formula: see text] norm regularization terms are used to ensure the sparsity of the feature matrices obtained by non-negative matrix factorization. To solve the model, iPALM-DLMF adopts non-negative double singular value decomposition to initialize the nonnegative matrix factorization, and an inertial Proximal Alternating Linearized Minimization iterating process, which has been proved to converge to a KKT point, to obtain the final result of the matrix factorization. Extensive experimental results show that iPALM-DLMF has better performance than other state-of-the-art methods. In case studies, in 50 highest-scoring proteins targeted by the drug gabapentin predicted by iPALM-DLMF, 46 have been validated, and in 50 highest-scoring drugs targeting prostaglandin-endoperoxide synthase 2 predicted by iPALM-DLMF, 47 have been validated.

Identifying potential biomarkers for non-obstructive azoospermia using WGCNA and machine learning algorithms

Objective

The cause and mechanism of non-obstructive azoospermia (NOA) is complicated; therefore, an effective therapy strategy is yet to be developed. This study aimed to analyse the pathogenesis of NOA at the molecular biological level and to identify the core regulatory genes, which could be utilised as potential biomarkers.

Methods

Three NOA microarray datasets (GSE45885, GSE108886, and GSE145467) were collected from the GEO database and merged into training sets; a further dataset (GSE45887) was then defined as the validation set. Differential gene analysis, consensus cluster analysis, and WGCNA were used to identify preliminary signature genes; then, enrichment analysis was applied to these previously screened signature genes. Next, 4 machine learning algorithms (RF, SVM, GLM, and XGB) were used to detect potential biomarkers that are most closely associated with NOA. Finally, a diagnostic model was constructed from these potential biomarkers and visualised as a nomogram. The differential expression and predictive reliability of the biomarkers were confirmed using the validation set. Furthermore, the competing endogenous RNA network was constructed to identify the regulatory mechanisms of potential biomarkers; further, the CIBERSORT algorithm was used to calculate immune infiltration status among the samples.

Results

A total of 215 differentially expressed genes (DEGs) were identified between NOA and control groups (27 upregulated and 188 downregulated genes). The WGCNA results identified 1123 genes in the MEblue module as target genes that are highly correlated with NOA positivity. The NOA samples were divided into 2 clusters using consensus clustering; further, 1027 genes in the MEblue module, which were screened by WGCNA, were considered to be target genes that are highly correlated with NOA classification. The 129 overlapping genes were then established as signature genes. The XGB algorithm that had the maximum AUC value (AUC=0.946) and the minimum residual value was used to further screen the signature genes. IL20RB, C9orf117, HILS1, PAOX, and DZIP1 were identified as potential NOA biomarkers. This 5 biomarker model had the highest AUC value, of up to 0.982, compared to other single biomarker models; additionally, the results of this biomarker model were verified in the validation set.

Conclusions

As IL20RB, C9orf117, HILS1, PAOX, and DZIP1 have been determined to possess the strongest association with NOA, these five genes could be used as potential therapeutic targets for NOA patients. Furthermore, the model constructed using these five genes, which possessed the highest diagnostic accuracy, may be an effective biomarker model that warrants further experimental validation.

Context-dependent function of the transcriptional regulator Rap1 in gene silencing and activation in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, heterochromatin is formed through interactions between site-specific DNA-binding factors, including the transcriptional activator Repressor Activator Protein (Rap1), and Sir proteins. Despite an understanding of the establishment and maintenance of Sir-silenced chromatin, the mechanism of gene silencing by Sir proteins has remained a mystery. Utilizing high-resolution chromatin immunoprecipitation, we found that Rap1, the native activator of the bidirectional HMLα promoter, bound its recognition sequence in silenced chromatin, and its binding was enhanced by the presence of Sir proteins. In contrast to prior results, various components of transcription machinery were not able to access HMLα in the silenced state. These findings disproved the long-standing model of indiscriminate steric occlusion by Sir proteins and led to investigation of the role of the transcriptional activator Rap1 in Sir-silenced chromatin. Using a highly sensitive assay that monitors loss-of-silencing events, we identified a role for promoter-bound Rap1 in the maintenance of silent chromatin through interactions with the Sir complex. We also found that promoter-bound Rap1 activated HMLα when in an expressed state, and aided in the transition from transcription initiation to elongation. Highlighting the importance of epigenetic context in transcription factor function, these results point toward a model in which the duality of Rap1 function was mediated by local chromatin environment rather than binding-site availability.

pLM-BLAST: distant homology detection based on direct comparison of sequence representations from protein language models

MOTIVATION

The detection of homology through sequence comparison is a typical first step in the study of protein function and evolution. In this work, we explore the applicability of protein language models to this task.

RESULTS

We introduce pLM-BLAST, a tool inspired by BLAST, that detects distant homology by comparing single-sequence representations (embeddings) derived from a protein language model, ProtT5. Our benchmarks reveal that pLM-BLAST maintains a level of accuracy on par with HHsearch for both highly similar sequences (with >50% identity) and markedly divergent sequences (with <30% identity), while being significantly faster. Additionally, pLM-BLAST stands out among other embedding-based tools due to its ability to compute local alignments. We show that these local alignments, produced by pLM-BLAST, often connect highly divergent proteins, thereby highlighting its potential to uncover previously undiscovered homologous relationships and improve protein annotation.

AVAILABILITY AND IMPLEMENTATION

pLM-BLAST is accessible via the MPI Bioinformatics Toolkit as a web server for searching precomputed databases (https://toolkit.tuebingen.mpg.de/tools/plmblast). It is also available as a standalone tool for building custom databases and performing batch searches (https://github.com/labstructbioinf/pLM-BLAST).

Histone regulator KAT2A acts as a potential biomarker related to tumor microenvironment and prognosis of diffuse large B cell lymphoma

BACKGROUND

Recent studies have indicated that epigenetic alterations contribute significantly to lymphoma pathogenesis. A type of epigenetic regulation known as histone acetylation plays a crucial role in transcriptional regulation in eukaryotic cells. Specifically, a significant effect of histone acetylation modifications on the abnormal progression and microenvironment of diffuse large B-cell lymphoma (DLBCL) has been observed.

METHODS

To provide insight into the significance of histone acetylation-related genes, we developed a HAscore model for analyzing histone acetylation patterns in DLBCL samples. Furthermore, KAT2A, a regulator of histone acetylation, was knocked down in DLBCL cell lines to investigate its role in proliferation, cell cycle, and apoptosis.

RESULTS

The HAscore model has been demonstrated to provide insight into the significance of these patterns, showing that patients with a low HAscore have distinct tumor immune microenvironments and poorer prognoses. Besides, KAT2A was identified as a potential biomarker related to immune infiltration and malignant pathways in DLBCL.

CONCLUSION

According to these findings, it is evident that the histone acetylation pattern score model is helpful in describing the immune status of DLBCL and that KAT2A may be used as a biomarker for its treatment.

First field and laboratory evaluation of LAMP assay for malaria diagnosis in Cubal, Angola

BACKGROUND

Malaria is a globally distributed infectious disease. According to the World Health Organization, Angola is one of the six countries that account for over half the global malaria burden in terms of both malaria cases and deaths. Diagnosis of malaria still depends on microscopic examination of thin and thick blood smears and rapid diagnostic tests (RDTs), which often lack analytical and clinical sensitivity. Molecular methods could overcome these disadvantages. The aim of this study was to evaluate, for the first time to our knowledge, the performance of a loop-mediated isothermal amplification (LAMP) for the diagnosis of malaria in an endemic area in Cubal, Angola, and to assess the reproducibility at a reference laboratory.

METHODS

A total of 200 blood samples from patients attended at Hospital Nossa Senhora da Paz, Cubal, Angola, were analysed for Plasmodium spp. detection by microscopy, RDTs, and LAMP. LAMP assay was easily performed in a portable heating block, and the results were visualized by a simple colour change. Subsequently, the samples were sent to a reference laboratory in Spain to be reanalysed by the same colorimetric LAMP assay and also in real-time LAMP format.

RESULTS

In field tests, a total of 67/200 (33.5%) blood samples were microscopy-positive for Plasmodium spp., 98/200 RDT positive, and 112/200 (56%) LAMP positive. Using microscopy as reference standard, field LAMP detected more microscopy-positive samples than RDTs (66/67; 98% vs. 62/67; 92.5%). When samples were reanalysed at a reference laboratory in Spain using both colorimetric and real-time assays, the overall reproducibility achieved 84.5%.

CONCLUSIONS

This is the first study to our knowledge in which LAMP has been clinically evaluated on blood samples in a resource-poor malaria-endemic area. The colorimetric LAMP proved to be more sensitive than microscopy and RDTs for malaria diagnosis in field conditions. Furthermore, LAMP showed an acceptable level of reproducibility in a reference laboratory. The possibility to use LAMP in a real-time format in a portable device reinforces the reliability of the assay for molecular diagnosis of malaria in resource-poor laboratories in endemic areas.

Adipose Mesenchymal Stem Cell-Derived Exosomes Promote the Regeneration of Corneal Endothelium Through Ameliorating Senescence

Purpose

Human corneal endothelial cells (hCECs) have been considered unable to regenerate in vivo, resulting in corneal decompensation after significant loss of hCECs. adipose-derived mesenchymal stem cell (ASC)-derived exosomes can regenerate tissues and organs. In this study, we investigated whether ASC-derived exosomes could protect and regenerate CECs.

Methods

We performed cell viability and cell-cycle analyses to evaluate the effect of ASC-derived exosomes on the regeneration capacity of cultured hCECs. Transforming growth factor-β (TGF-β) and hydrogen peroxide (H2O2) were used to induce biological stress in CECs. The effect of ASC-derived exosomes on CECs was investigated in vivo. ASC-derived exosomes were introduced into rat CECs using electroporation, and rat corneas were injured using cryoinjury. Next-generation sequencing analysis was performed to compare the differentially expressed microRNAs (miRNAs) between ASC-derived and hCEC-derived exosomes.

Results

ASC-derived exosomes induced CEC proliferation and suppressed TGF-β- or H2O2-induced oxidative stress and senescence. ASC-derived exosomes protect hCECs against TGF-β- or H2O2-induced endothelial-mesenchymal transition and mitophagy. In an in vivo study, ASC-derived exosomes promoted wound healing of rat CECs and protected the corneal endothelium against cryoinjury-induced corneal endothelium damage. Next-generation sequencing analysis revealed differentially expressed miRNAs for ASC-derived and hCEC-derived exosomes. They are involved in lysine degradation, adherens junction, the TGF-β signaling pathway, the p53 signaling pathway, the Hippo signaling pathway, the forkhead box O (FoxO) signaling pathway, regulation of actin cytoskeleton, and RNA degradation based on Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis.

Conclusions

ASC-derived exosomes promoted wound healing and regeneration of endothelial cells by inducing a shift in the cell cycle and suppressing senescence and autophagy.

Broad functional profiling of fission yeast proteins using phenomics and machine learning

Many proteins remain poorly characterized even in well-studied organisms, presenting a bottleneck for research. We applied phenomics and machine-learning approaches with Schizosaccharomyces pombe for broad cues on protein functions. We assayed colony-growth phenotypes to measure the fitness of deletion mutants for 3509 non-essential genes in 131 conditions with different nutrients, drugs, and stresses. These analyses exposed phenotypes for 3492 mutants, including 124 mutants of 'priority unstudied' proteins conserved in humans, providing varied functional clues. For example, over 900 proteins were newly implicated in the resistance to oxidative stress. Phenotype-correlation networks suggested roles for poorly characterized proteins through 'guilt by association' with known proteins. For complementary functional insights, we predicted Gene Ontology (GO) terms using machine learning methods exploiting protein-network and protein-homology data (NET-FF). We obtained 56,594 high-scoring GO predictions, of which 22,060 also featured high information content. Our phenotype-correlation data and NET-FF predictions showed a strong concordance with existing PomBase GO annotations and protein networks, with integrated analyses revealing 1675 novel GO predictions for 783 genes, including 47 predictions for 23 priority unstudied proteins. Experimental validation identified new proteins involved in cellular aging, showing that these predictions and phenomics data provide a rich resource to uncover new protein functions.

A predictive model for preterm infants with bronchopulmonary dysplasia based on ferroptosis-related lncRNAs

BACKGROUND

Bronchopulmonary dysplasia (BPD) is the most challenging chronic lung disease for prematurity, with difficulties in early identification. Given lncRNA emerging as a novel biomarker and the regulator of ferroptosis, this study aims to develop a BPD predictive model based on ferroptosis-related lncRNAs (FRLs).

METHODS

Using a rat model, we firstly explored mRNA levels of ferroptosis-related genes and ferrous iron accumulation in BPD rat lungs. Subsequently, a microarray dataset of umbilical cord tissue from 20 preterm infants with BPD and 34 preterm infants without BPD were downloaded from the Gene Expression Omnibus databases. Random forest and LASSO regression were conducted to identify diagnostic FRLs. Nomogram was used to construct a predictive BPD model based on the FRLs. Finally, umbilical cord blood lymphocytes of preterm infants born before 32 weeks gestational age and term infants were collected and determined the expression level of diagnostic FRLs by RT-qPCR.

RESULTS

Increased iron accumulation and several dysregulated ferroptosis-associated genes were found in BPD rat lung tissues, indicating that ferroptosis was participating in the development of BPD. By exploring the microarray dataset of preterm infants with BPD, 6 FRLs, namely LINC00348, POT1-AS1, LINC01103, TTTY8, PACRG-AS1, LINC00691, were determined as diagnostic FRLs for modeling. The area under the receiver operator characteristic curve of the model was 0.932, showing good discrimination of BPD. In accordance with our analysis of microarray dataset, the mRNA levels of FRLs were significantly upregulated in umbilical cord blood lymphocytes from preterm infants who had high risk of BPD.

CONCLUSION

The incorporation of FRLs into a predictive model offers a non-invasive approach to show promise in improving early detection and management of this challenging chronic lung disease in premature infant, enabling timely intervention and personalized treatment strategies.

Rhizochalinin Exhibits Anticancer Activity and Synergizes with EGFR Inhibitors in Glioblastoma In Vitro Models

Rhizochalinin (Rhiz) is a recently discovered cytotoxic sphingolipid synthesized from the marine natural compound rhizochalin. Previously, Rhiz demonstrated high in vitro and in vivo efficacy in various cancer models. Here, we report Rhiz to be highly active in human glioblastoma cell lines as well as in patient-derived glioma-stem like neurosphere models. Rhiz counteracted glioblastoma cell proliferation by inducing apoptosis, G2/M-phase cell cycle arrest, and inhibition of autophagy. Proteomic profiling followed by bioinformatic analysis suggested suppression of the Akt pathway as one of the major biological effects of Rhiz. Suppression of Akt as well as IGF-1R and MEK1/2 kinase was confirmed in Rhiz-treated GBM cells. In addition, Rhiz pretreatment resulted in a more pronounced inhibitory effect of γ-irradiation on the growth of patient-derived glioma-spheres, an effect to which the Akt inhibition may also contribute decisively. In contrast, EGFR upregulation, observed in all GBM neurospheres under Rhiz treatment, was postulated to be a possible sign of incipient resistance. In line with this, combinational therapy with EGFR-targeted tyrosine kinase inhibitors synergistically increased the efficacy of Rhiz resulting in dramatic inhibition of GBM cell viability as well as a significant reduction of neurosphere size in the case of combination with lapatinib. Preliminary in vitro data generated using a parallel artificial membrane permeability (PAMPA) assay suggested that Rhiz cannot cross the blood brain barrier and therefore alternative drug delivery methods should be used in the further in vivo studies. In conclusion, Rhiz is a promising new candidate for the treatment of human glioblastoma, which should be further developed in combination with EGFR inhibitors.

Epigenetic Profiles of Triple-Negative Breast Cancers of African American and White Females

Importance

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype and appears to have disproportionately higher incidence and worse outcomes among younger African American females.

Objective

To investigate whether epigenetic differences exist in TNBCs of younger African American females that may explain clinical disparities seen in this patient group.

Design, Setting, and Participants

This cross-sectional study used clinical, demographic, DNA methylation (HumanMethylation450; Illumina), and gene expression (RNA sequencing) data for US patient populations from publicly available data repositories (The Cancer Genome Atlas [TCGA], 2006-2012, and Gene Expression Omnibus [GEO], 2004-2013) accessed on April 13, 2021. White and African American females with TNBC identified in TCGA (69 patients) and a validation cohort of 210 African American patients from GEO (GSE142102) were included. Patients without available race or age data were excluded. Data were analyzed from September 2022 through April 2023.

Main Outcomes and Measures

DNA methylation and gene expression profiles of TNBC tumors by race (self-reported) and age were assessed. Age was considered a dichotomous variable using age 50 years as the cutoff (younger [<50 years] vs older [≥50 years]).

Results

A total of 69 female patients (34 African American [49.3%] and 35 White [50.7%]; mean [SD; range] age, 55.7 [11.6; 29-82] years) with TNBC were included in the DNA methylation analysis; these patients and 210 patients in the validation cohort were included in the gene expression analysis (279 patients). There were 1115 differentially methylated sites among younger African American females. The DNA methylation landscape on TNBC tumors in this population had increased odds of enrichment of hormone (odds ratio [OR], 1.82; 95% CI, 1.21 to 2.67; P = .003), muscle (OR, 1.85; 95% CI, 1.44 to 2.36; P < .001), and proliferation (OR, 3.14; 95% CI, 2.71 to 3.64; P < .001) pathways vs other groups (older African American females and all White females). Alterations in regulators of these molecular features in TNBCs of younger African American females were identified involving hormone modulation (downregulation of androgen receptor: fold change [FC] = -2.93; 95% CI, -4.76 to -2.11; P < .001) and upregulation of estrogen-related receptor α (FC = 0.86; 95% CI, 0.34 to 1.38; P = .002), muscle metabolism (upregulation of FOXC1: FC = 1.33; 95% CI, 0.62 to 2.03; P < .001), and proliferation mediators (upregulation of NOTCH1: FC = 0.71; 95% CI, 0.23 to 1.19; P = .004 and MYC (FC = 0.81; 95% CI, 0.18 to 1.45; P = .01).

Conclusions and Relevance

These findings suggest that TNBC of younger African American females may represent a distinct epigenetic entity and offer novel insight into molecular alterations associated with TNBCs of this population. Understanding these epigenetic differences may lead to the development of more effective therapies for younger African American females, who have the highest incidence and worst outcomes from TNBC of any patient group.

Descriptor-Free Deep Learning QSAR Model for the Fraction Unbound in Human Plasma

Chemical-specific parameters are either measured in vitro or estimated using quantitative structure-activity relationship (QSAR) models. The existing body of QSAR work relies on extracting a set of descriptors or fingerprints, subset selection, and training a machine learning model. In this work, we used a state-of-the-art natural language processing model, Bidirectional Encoder Representations from Transformers, which allowed us to circumvent the need for calculation of these chemical descriptors. In this approach, simplified molecular-input line-entry system (SMILES) strings were embedded in a high-dimensional space using a two-stage training approach. The model was first pre-trained on a masked SMILES token task and then fine-tuned on a QSAR prediction task. The pre-training task learned meaningful high-dimensional embeddings based upon the relationships between the chemical tokens in the SMILES strings derived from the "in-stock" portion of the ZINC 15 dataset─a large dataset of commercially available chemicals. The fine-tuning task then perturbed the pre-trained embeddings to facilitate prediction of a specific QSAR endpoint of interest. The power of this model stems from the ability to reuse the pre-trained model for multiple different fine-tuning tasks, reducing the computational burden of developing multiple models for different endpoints. We used our framework to develop a predictive model for fraction unbound in human plasma (fu,p). This approach is flexible, requires minimum domain expertise, and can be generalized for other parameters of interest for rapid and accurate estimation of absorption, distribution, metabolism, excretion, and toxicity.

Targeting the melanoma-associated antigen CSPG4 with HLA-C*07:01-restricted T-cell receptors

Intorduction

Chondroitin sulfate proteoglycan 4 (CSPG4), also known as high molecular weight-melanoma associated antigen, is expressed in melanoma but also other tumor entities and constitutes an attractive target for immunotherapeutic approaches. While recent preclinical reports focused on anti-CSPG4 chimeric antigen receptors (CAR), we here explore T-cell receptor (TCR)-based approaches targeting CSPG4.

Methods

The TCRs of two CSPG4-reactive T-cell clones (11C/73 and 2C/165) restricted by the highly prevalent HLA-C*07:01 allele were isolated and the respective αβTCR pairs were retrovirally expressed in CRISPR/Cas9-edited TCR-knockout T cells for functional testing. We also combined alpha and beta TCR chains derived from 11C/73 and 2C/165 in a cross-over fashion to assess for hemichain dominance. CSPG4+ melanoma, glioblastoma and lung cancer cell lines were identified and, if negative, retrovirally transduced with HLA-C*07:01.

Results

Functional tests confirmed specific recognition of CSPG4+HLA-C*07:01+ target cells by the αβTCR retrieved from the parental T-cell clones and in part also by the cross-over TCR construct 2Cα-11Cβ. Despite high surface expression, the 11Cα-2Cβ combination, however, was not functional.

Discussion

Collectively, 11C/73- and 2C/165-expressing T cells specifically and efficiently recognized CSPG4+HLA-C*07:01+ cancer cells which warrants further preclinical and clinical evaluation of these TCRs.

Life-History and Ecological Correlates of Egg and Clutch Mass Variation in Sympatric Bird Species at High Altitude

The variation in egg and clutch mass in sympatric species at high altitudes is poorly understood, and the potential causes of variation are rarely investigated. This study aimed to describe the interspecific variation in avian egg and clutch mass among 22 sympatric bird species at an altitude of 3430 m. Our objective was to reduce potential confounding effects of biotic/abiotic factors and investigated hypotheses concerning allometry, clutch size, parental care, nest predation, and lifespan as possible correlates and explanations for the observed variation. Our findings indicated that both egg and clutch mass evolve with body mass across species. We found that egg mass variation was not explained by clutch size when controlling for allometric effects, which contrasts the "egg mass vs. clutch size trade-off" hypothesis. Additionally, we found that clutch mass was positively associated with parental care (reflected by development period) but negatively associated with predation rate. By substituting egg mass and clutch size into the models, we found that clutch size was significantly correlated with parental care, predation rate, and lifespan, while egg mass was only significantly associated with development period. Overall, these findings support life-history theories suggesting that reduced clutch size or mass is associated with a higher risk of predation, reduced parental care, but longer adult lifespan. Interestingly, our results indicate that clutch size has a greater influence on these factors compared to egg mass. This could be attributed to the fact that smaller clutch sizes result in a more notable decrease in energetic allocation, as they require a reduced effort in terms of offspring production, incubation, and feeding, as opposed to solely reducing egg size. These findings contribute to the growing evidence that life-history and ecological traits correlate with egg and clutch mass variation in sympatric species. However, further research is needed to explore the potential evolutionary causes underlying these patterns.

Time Spent Jogging/Running and Biological Aging in 4458 U.S. Adults: An NHANES Investigation

Telomere length is a good index of cellular aging. Longer telomeres are predictive of longer life, and healthy lifestyles are associated with longer telomeres. This study explored the relationship between time spent jogging or running each week and leukocyte telomere length (LTL) in 4458 randomly selected U.S. adults. The association was studied using data collected by the National Health and Nutrition Examination Survey (NHANES), and a cross-sectional design. Total weekly jog/run time was calculated from survey responses. From the minute totals, three categories were formed: <10 min/week, 10-74 min/week, and ≥75 min/week. Adults in the third category met the U.S. guidelines. Data were analyzed using one-way ANOVA. Partial correlation was used to adjust for differences in potential mediating factors, including demographic and lifestyle/medical factors. In the total sample, after adjusting for all the potential covariates, mean LTL significantly differed across the three jog/run categories (F = 4.1, p = 0.0272). Specifically, adults who met the guidelines via jogging and/or running had significantly longer telomeres than adults who performed no jogging/running. Adults in the middle category did not differ from the other two categories. A minimum of 75 min of jogging/running weekly is predictive of longer telomeres when compared to adults who do not jog or run regularly.

Metabolomics and machine learning technique revealed that germination enhances the multi-nutritional properties of pigmented rice

Enhancing the dietary properties of rice is crucial to contribute to alleviating hidden hunger and non-communicable diseases in rice-consuming countries. Germination is a bioprocessing approach to increase the bioavailability of nutrients in rice. However, there is a scarce information on how germination impacts the overall nutritional profile of pigmented rice sprouts (PRS). Herein, we demonstrated that germination resulted to increase levels of certain dietary compounds, such as free phenolics and micronutrients (Ca, Na, Fe, Zn, riboflavin, and biotin). Metabolomic analysis revealed the preferential accumulation of dipeptides, GABA, and flavonoids in the germination process. Genome-wide association studies of the PRS suggested the activation of specific genes such as CHS1 and UGT genes responsible for increasing certain flavonoid compounds. Haplotype analyses showed a significant difference (P < 0.05) between alleles associated with these genes. Genetic markers associated with these flavonoids were incorporated into the random forest model, improving the accuracy of prediction of multi-nutritional properties from 89.7% to 97.7%. Deploying this knowledge to breed rice with multi-nutritional properties will be timely to address double burden nutritional challenges.

Revision of Litoria rothii (Anura: Pelodryadidae) from northern Australia

Litoria rothii is a widespread pelodryadid frog with a charismatic laughing advertisement call, distributed across the Australian Monsoon Tropics and southern New Guinea. Given its large distribution spanning well-known biogeographic barriers, variation in male advertisement calls and the prevalence of unresolved species complexes in the Australian frog fauna, we examine the genetic, morphological and acoustic diversity in the species from across its range. Our analyses reveal the presence of a previously unrecognised species in western parts of the range of L. rothii sensu lato, which we describe herein as a new species. Litoria ridibunda sp. nov. is distinguished from L. rothii on the basis of paraphyly of nuclear gene trees with L. everetti from Indonesia, colour patterns on the posterior thigh and male advertisement calls. Compared to L. rothii, the new species has a less contrasting pattern on the posterior thigh and a male advertisement call with a greater number of notes per call and a greater call duration. In particular, the magnitude of call differences between the species is highest where the ranges of the two species are in proximity in north-western Queensland. Our study further emphasises the undiagnosed diversity that remains in Australian frogs, even in relatively large, charismatic, frequently encountered species that often share human dwellings.

Type I interferons drive MAIT cell functions against bacterial pneumonia

Mucosal-associated invariant T (MAIT) cells are abundant in the lung and contribute to host defense against infections. During bacterial infections, MAIT cell activation has been proposed to require T cell receptor (TCR)-mediated recognition of antigens derived from the riboflavin synthesis pathway presented by the antigen-presenting molecule MR1. MAIT cells can also be activated by cytokines in an MR1-independent manner, yet the contribution of MR1-dependent vs. -independent signals to MAIT cell functions in vivo remains unclear. Here, we use Klebsiella pneumoniae as a model of bacterial pneumonia and demonstrate that MAIT cell activation is independent of MR1 and primarily driven by type I interferons (IFNs). During Klebsiella infection, type I IFNs stimulate activation of murine and human MAIT cells, induce a Th1/cytotoxic transcriptional program, and modulate MAIT cell location within the lungs. Consequently, adoptive transfer or boosting of pulmonary MAIT cells protect mice from Klebsiella infection, with protection being dependent on direct type I IFN signaling on MAIT cells. These findings reveal type I IFNs as new molecular targets to manipulate MAIT cell functions during bacterial infections.

Anti-methicillin-resistant Staphylococcus aureus and antibiofilm activity of new peptides produced by a Brevibacillus strain

Background

Methicillin-resistant Staphylococcus aureus (MRSA) is listed as a highly prioritized pathogen by the World Health Organization (WHO) to search for effective antimicrobial agents. Previously, we isolated a soil Brevibacillus sp. strain SPR19 from a botanical garden, which showed anti-MRSA activity. However, the active substances were still unknown.

Methods

The cell-free supernatant of this bacterium was subjected to salt precipitation, cation exchange, and reversed-phase chromatography. The antimicrobial activity of pure substances was determined by broth microdilution assay. The peptide sequences and secondary structures were characterized by tandem mass spectroscopy and circular dichroism (CD), respectively. The most active anti-MRSA peptide underwent a stability study, and its mechanism was determined through scanning electron microscopy, cell permeability assay, time-killing kinetics, and biofilm inhibition and eradication. Hemolysis was used to evaluate the peptide toxicity.

Results

The pure substances (BrSPR19-P1 to BrSPR19-P5) were identified as new peptides. Their minimum inhibition concentration (MIC) and minimum bactericidal concentration (MBC) against S. aureus and MRSA isolates ranged from 2.00 to 32.00 and 2.00 to 64.00 µg/mL, respectively. The sequence analysis of anti-MRSA peptides revealed a length ranging from 12 to 16 residues accompanied by an amphipathic structure. The physicochemical properties of peptides were predicted such as pI (4.25 to 10.18), net charge at pH 7.4 (-3 to +4), and hydrophobicity (0.12 to 0.96). The CD spectra revealed that all peptides in the water mainly contained random coil structures. The increased proportion of α-helix structure was observed in P2-P5 when incubated with SDS. P2 (NH2-MFLVVKVLKYVV-COOH) showed the highest antimicrobial activity and high stability under stressed conditions such as temperatures up to 100 °C, solution of pH 3 to 10, and proteolytic enzymes. P2 disrupted the cell membrane and caused bacteriolysis, in which its action was dependent on the incubation time and peptide concentration. Antibiofilm activity of P2 was determined by which the half-maximal inhibition of biofilm formation was observed at 2.92 and 4.84 µg/mL for S. aureus TISTR 517 and MRSA isolate 2468, respectively. Biofilm eradication of tested pathogens was found at the P2 concentration of 128 µg/mL. Furthermore, P2 hemolytic activity was less than 10% at concentrations up to 64 µg/mL, which reflected the hemolysis index thresholds of 32.

Conclusion

Five novel anti-MRSA peptides were identified from SPR19. P2 was the most active peptide and was demonstrated to cause membrane disruption and cell lysis. The P2 activity was dependent on the peptide concentration and exposure time. This peptide had antibiofilm activity against tested pathogens and was compatible with human erythrocytes, supporting its potential use as an anti-MRSA agent in this post-antibiotic era.

FOXO1 orchestrates the intestinal homeostasis via neuronal signaling in group 3 innate lymphoid cells

The neuro-immune regulation is associated with homeostasis of the intestine. Intestinal group 3 innate lymphoid cells (ILC3s) are tissue-resident lymphocytes whose functions are affected by the intestine niche. However, how a gut neuronal signal coordinates the immune response of ILC3s is largely unknown. Here, we found that cyclic adenosine monophosphate (cAMP) signaling exacerbated the inflammatory response and attenuated the expression level of the transcription factor forkhead box O1 (FOXO1) in ILC3s. Deficiency of FOXO1 drove the hyperactivation of ILC3s and resulted in gut inflammation independently of T cells. Mechanistically, FOXO1 promoted the transcription of neuropeptide receptor VIPR2 and inhibited the transcription of adrenoceptor ADRA2A in ILC3s. FOXO1-related regulation of VIPR2 and ADRA2A signaling balanced the activation of ILC3s under steady condition or during colitis. Moreover, chronic stress elevated cAMP level and downregulated FOXO1 level, exacerbating intestinal inflammation. Our findings reveal that FOXO1 balances the activation of ILC3s via VIP and adrenergic signaling and regulates intestinal homeostasis.