PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools

Cardiac tissue-resident vesicles differentially modulate anti-fibrotic phenotype by age and sex through synergistic miRNA effects

Aging is a risk factor for cardiovascular disease, the leading cause of death worldwide. Cardiac fibrosis is a harmful result of repeated myocardial infarction that increases risk of morbidity and future injury. Interestingly, both rates and outcomes of cardiac fibrosis differ between young and aged individuals, as well as men and women. Here, for the first time, we identify and isolate matrix-bound extracellular vesicles from the left ventricles (LVs) of young or aged males and females in both human and murine models. These LV vesicles (LVVs) show differences in morphology and content between these four cohorts in both humans and mice. LVV effects on fibrosis were also investigated in vitro, and aged male LVVs were pro-fibrotic while other LVVs were anti-fibrotic. From these LVVs, we could identify therapeutic miRNAs to promote anti-fibrotic effects. Four miRNAs were identified and together, but not individually, demonstrated significant cardioprotective effects when transfected. This suggests that miRNA synergy can regulate cell response, not just individual miRNAs, and also indicates that biological agent-associated therapeutic effects may be recapitulated using non-immunologically active agents. Furthermore, that chronic changes in LVV miRNA content may be a major factor in sex- and age-dependent differences in clinical outcomes of cardiac fibrosis.

Long non-coding RNA variability in porcine skeletal muscle

Recently, numerous studies including various tissues have been carried out on long non-coding RNAs (lncRNAs), but still, its variability has not yet been fully understood. In this study, we characterised the inter-individual variability of lncRNAs in pigs, in the context of number, length and expression. Transcriptomes collected from muscle tissue belonging to six Polish Landrace boars (PL1-PL6), including half-brothers (PL1-PL3), were investigated using bioinformatics (lncRNA identification and functional analysis) and statistical (lncRNA variability) methods. The number of lncRNA ranged from 1289 to 3500 per animal, and the total number of common lncRNAs among all boars was 232. The number, length and expression of lncRNAs significantly varied between individuals, and no consistent pattern has been found between pairs of half-brothers. In detail, PL5 exhibits lower expression than the others, while PL4 has significantly higher expression than PL2-PL3 and PL5-PL6. Noteworthy, comparing the inter-individual variability of lncRNA and mRNA expression, they exhibited concordant patterns. The enrichment analysis for common lncRNA target genes determined a variety of biological processes that play fundamental roles in cell biology, and they were mostly related to whole-body homeostasis maintenance, energy and protein synthesis as well as dynamics of multiple nucleoprotein complexes. The high variability of lncRNA landscape in the porcine genome has been revealed in this study. The inter-individual differences have been found in the context of three aspects: the number, length and expression of lncRNAs, which contribute to a better understanding of its complex nature.

Atg8/LC3 controls systemic nutrient surplus signaling in flies and humans

Organisms experience constant nutritional flux. Mechanisms at the interface of opposing nutritional states-scarcity and surplus-enable organismal energy homeostasis. Contingent on nutritional stores, adipocytes secrete adipokines, such as the fat hormone leptin, to signal nutrient status to the central brain. Increased leptin secretion underlies metabolic dysregulation during common obesity, but the molecular mechanisms regulating leptin secretion from human adipocytes are poorly understood. Here, we report that Atg8/LC3 family proteins, best known for their role in autophagy during nutrient scarcity, play an evolutionarily conserved role during nutrient surplus by promoting adipokine secretion. We show that in a well-fed state, Atg8/LC3 promotes the secretion of the Drosophila functional leptin ortholog unpaired 2 (Upd2) and leptin from human adipocytes. Proteomic analyses reveal that LC3 directs leptin to a secretory pathway in human cells. We identified LC3-dependent extracellular vesicle (EV) loading and secretion (LDELS) as a required step for leptin release, highlighting a unique secretory route adopted by leptin in human adipocytes. In Drosophila, mutations to Upd2's Atg8 interaction motif (AIM) result in constitutive adipokine retention. Atg8-mediated Upd2 retention alters lipid storage and hunger response and rewires the bulk organismal transcriptome in a manner conducive to starvation survival. Thus, Atg8/LC3's bidirectional role in nutrient sensing-conveying nutrient surplus and responding to nutrient deprivation-enables organisms to manage nutrient flux effectively. We posit that decoding how bidirectional molecular switches-such as Atg8/LC3-operate at the nexus of nutritional scarcity and surplus will inform therapeutic strategies to tackle chronic metabolic disorders.

Molecular Changes Implicate Angiogenesis and Arterial Remodeling in Systemic Sclerosis-Associated and Idiopathic Pulmonary Hypertension

BACKGROUND

Pulmonary hypertension (PH) is a common complication of systemic sclerosis (SSc) and a leading cause of mortality among patients with this disease. PH can also occur as an idiopathic condition (idiopathic pulmonary arterial hypertension). Investigation of transcriptomic alterations in vascular populations is critical to elucidating cellular mechanisms underlying pathobiology of SSc-associated and idiopathic PH.

METHODS

We analyzed single-cell RNA sequencing profiles of endothelial and perivascular mesenchymal populations from explanted lung tissue of patients with SSc-associated PH (n=16), idiopathic pulmonary arterial hypertension (n=3), and healthy controls (n=15). Findings were validated by immunofluorescence staining of explanted human lung tissue.

RESULTS

Three disease-associated endothelial populations emerged. Two angiogenic endothelial cell (EC) subtypes markedly expanded in SSc-associated PH lungs: tip ECs expressing canonical tip markers PGF and APLN and phalanx ECs expressing genes associated with vascular development, endothelial barrier integrity, and Notch signaling. Gene regulatory network analysis suggested enrichment of Smad1 (SMAD family member 1) and PPAR-γ (peroxisome proliferator-activated receptor-γ) regulon activities in these 2 populations, respectively. Mapping of potential ligand-receptor interactions highlighted Notch, apelin-APJ (apelin receptor), and angiopoietin-Tie (tyrosine kinase with immunoglobulin-like and EGF-like domains 1) signaling pathways between angiogenic ECs and perivascular cells. Transitional cells, expressing both endothelial and pericyte/smooth muscle cell markers, provided evidence for the presence of endothelial-to-mesenchymal transition. Transcriptional programs associated with arterial endothelial dysfunction implicated VEGF-A (vascular endothelial growth factor-A), TGF-β1 (transforming growth factor beta-1), angiotensin, and TNFSF12 (tumor necrosis factor ligand superfamily member 12)/TWEAK (TNF-related weak inducer of apoptosis) in the injury/remodeling phenotype of PH arterial ECs.

CONCLUSIONS

These data provide high-resolution insights into the complexity and plasticity of the pulmonary endothelium in SSc-associated PH and idiopathic pulmonary arterial hypertension and provide direct molecular insights into soluble mediators and transcription factors driving PH vasculopathy.

Hypothalamic transcriptome response to simulated diel earthen pond hypoxia cycles in channel catfish (Ictalurus punctatus)

Commercial culture of channel catfish (Ictalurus punctatus) occurs in earthen ponds that are characterized by diel swings in dissolved oxygen concentration that can fall to severe levels of hypoxia, which can suppress appetite and lead to suboptimal growth. Given the significance of the hypothalamus in regulating these processes in other fishes, an investigation into the hypothalamus transcriptome was conducted to identify specific genes and expression patterns responding to hypoxia. Channel catfish in normoxic water were compared with catfish subjected to 12 h of hypoxia (20% oxygen saturation; 1.8 mg O2/L; 27°C) followed by 12 h of recovery in normoxia to mimic 24 h in a catfish aquaculture pond. Fish were sampled at 0-, 6-, 12-, 18-, and 24-h timepoints, with the 6- and 12-h samplings occurring during hypoxia. A total of 190 genes were differentially expressed during the experiment, with most occurring during hypoxia and returning to baseline values within 6 h of normoxia. Differentially expressed genes were sorted by function into Gene Ontology biological processes and revealed that most were categorized as "response to hypoxia," "sprouting angiogenesis," and "cellular response to xenobiotic stimulus." The patterns of gene expression reported here suggest that transcriptome responses to hypoxia are broad and quickly reversibly with the onset of normoxia. Although no genes commonly reported to modulate appetite were found to be differentially expressed in this experiment, several candidates were identified for future studies investigating the interplay between hypoxia and appetite in channel catfish, including adm, igfbp1a, igfbp7, and stc2b.NEW & NOTEWORTHY Channel catfish are an economically important species that experience diel episodic periods of hypoxia that can reduce appetite. This is the first study to investigate their transcriptome from the hypothalamus in a simulated 24-h span in a commercial catfish pond, with 12 h of hypoxia and 12 h of normoxia. The research revealed functional groups of genes relating to hypoxia, angiogenesis, and glycolysis as well as individual target genes possibly involved in appetite regulation.

Spatial variation in genomic signatures of local adaptation during the cane toad invasion of Australia

Adaptive evolution can facilitate species' range expansions across environmentally heterogeneous landscapes. However, serial founder effects can limit the efficacy of selection, and the evolution of increased dispersal during range expansions may result in gene flow swamping local adaptation. Here, we study how genetic drift, gene flow and selection interact during the cane toad's (Rhinella marina) invasion across the heterogeneous landscape of Australia. Following its introduction in 1935, the cane toad colonised eastern Australia and established several stable range edges. The ongoing, more rapid range expansion in north-central Australia has occurred concomitant with an evolved increase in dispersal capacity. Using reduced representation genomic data of Australian cane toads from the expansion front and from two areas of their established range, we test the hypothesis that high gene flow constrains local adaptation at the expansion front relative to established areas. Genetic analyses indicate the three study areas are genetically distinct but show similar levels of allelic richness, heterozygosity and inbreeding. Markedly higher gene flow or recency of colonisation at the expansion front have likely hindered local adaptation at the time of sampling, as indicated by reduced slopes of genetic-environment associations (GEAs) estimated using a novel application of geographically weighted regression that accounts for allele surfing; GEA slopes are significantly steeper in established parts of the range. Our work bolsters evidence supporting adaptation of invasive species post-introduction and adds novel evidence for differing strengths of evolutionary forces among geographic areas with different invasion histories.

Thymic stromal lymphopoietin induces IL-4/IL-13 from T cells to promote sebum secretion and adipose loss

BACKGROUND

The cytokine TSLP promotes type 2 immune responses and can induce adipose loss by stimulating lipid loss from the skin through sebum secretion by sebaceous glands, which enhances the skin barrier. However, the mechanism by which TSLP upregulates sebaceous gland function is unknown.

OBJECTIVES

This study investigated the mechanism by which TSLP stimulates sebum secretion and adipose loss.

METHODS

RNA-sequencing analysis was performed on sebaceous glands isolated by laser capture microdissection and single-cell RNA-sequencing analysis was performed on sorted skin T cells. Sebocyte function was analyzed by histological analysis and sebum secretion in vivo and by measuring lipogenesis and proliferation in vitro.

RESULTS

This study found that TSLP sequentially stimulated the expression of lipogenesis genes followed by cell death genes in sebaceous glands to induce holocrine secretion of sebum. TSLP did not affect sebaceous gland activity directly. Rather, single-cell RNA-sequencing revealed that TSLP recruited distinct T-cell clusters that produce IL-4 and IL-13, which were necessary for TSLP-induced adipose loss and sebum secretion. Moreover, IL-13 was sufficient to cause sebum secretion and adipose loss in vivo and to induce lipogenesis and proliferation of a human sebocyte cell line in vitro.

CONCLUSIONS

This study proposes that TSLP stimulates T cells to deliver IL-4 and IL-13 to sebaceous glands, which enhances sebaceous gland function, turnover, and subsequent adipose loss.

A three-dimensional valve-on-chip microphysiological system implicates cell cycle progression, cholesterol metabolism and protein homeostasis in early calcific aortic valve disease progression

BACKGROUND

Calcific aortic valve disease (CAVD) is one of the most common forms of valvulopathy, with a 50 % elevated risk of a fatal cardiovascular event, and greater than 15,000 annual deaths in North America alone. The treatment standard is valve replacement as early diagnostic, mitigation, and drug strategies remain underdeveloped. The development of early diagnostic and therapeutic strategies requires the fabrication of effective in vitro valve mimetic models to elucidate early CAVD mechanisms.

METHODS

In this study, we developed a multilayered physiologically relevant 3D valve-on-chip (VOC) system that incorporated aortic valve mimetic extracellular matrix (ECM), porcine aortic valve interstitial cell (VIC) and endothelial cell (VEC) co-culture and dynamic mechanical stimuli. Collagen and glycosaminoglycan (GAG) based hydrogels were assembled in a bilayer to mimic healthy or diseased compositions of the native fibrosa and spongiosa. Multiphoton imaging and proteomic analysis of healthy and diseased VOCs were performed.

RESULTS

Collagen-based bilayered hydrogel maintained the phenotype of the VICs. Proteins related to cellular processes like cell cycle progression, cholesterol biosynthesis, and protein homeostasis were found to be significantly altered and correlated with changes in cell metabolism in diseased VOCs. This study suggested that diseased VOCs may represent an early, adaptive disease initiation stage, which was corroborated by human aortic valve proteomic assessment.

CONCLUSIONS

In this study, we developed a collagen-based bilayered hydrogel to mimic healthy or diseased compositions of the native fibrosa and spongiosa layers. When the gels were assembled in a VOC with VECs and VICs, the diseased VOCs revealed key insights about the CAVD initiation process.

STATEMENT OF SIGNIFICANCE

Calcific aortic valve disease (CAVD) elevates the risk of death due to cardiovascular pathophysiology by 50 %, however, prevention and mitigation strategies are lacking, clinically. Developing tools to assess early disease would significantly aid in the prevention of disease and in the development of therapeutics. Previously, studies have utilized collagen and glycosaminoglycan-based hydrogels for valve cell co-cultures, valve cell co-cultures in dynamic environments, and inorganic polymer-based multilayered hydrogels; however, these approaches have not been combined to make a physiologically relevant model for CAVD studies. We fabricated a bi-layered hydrogel that closely mimics the aortic valve and used it for valve cell co-culture in a dynamic platform to gain mechanistic insights into the CAVD initiation process using proteomic and multiphoton imaging assessment.

Transcriptional drift in aging cells: A global decontroller

As cells age, they undergo a remarkable global change: In transcriptional drift, hundreds of genes become overexpressed while hundreds of others become underexpressed. Using archetype modeling and Gene Ontology analysis on data from aging Caenorhabditis elegans worms, we find that the up-regulated genes code for sensory proteins upstream of stress responses and down-regulated genes are growth- and metabolism-related. We observe similar trends within human fibroblasts, suggesting that this process is conserved in higher organisms. We propose a simple mechanistic model for how such global coordination of multiprotein expression levels may be achieved by the binding of a single factor that concentrates with age in C. elegans. A key implication is that a cell's own responses are part of its aging process, so unlike wear-and-tear processes, intervention might be able to modulate these effects.

Label-free liquid chromatography-mass spectrometry comparison of the breast muscle proteome profiles in two fast-growing broilers

Poultry meat-production is increasing worldwide; leading to the selection of chickens for meat-production that show a fast growth. A label-free quantitative proteomic-approach and Western-blot were applied to investigate the dynamics of muscle protein under rapid growth conditions in two common fast-growing broiler genetic-lines (Ross 508 and AZ Extra Heavy Red-chicken). Muscle exudate from chicken Pectoralis major was used as substrate to unveil the proteome of these genetic-lines. Six-hundred forty-five proteins were identified in total from all samples, and after statistical-analysis 172 proteins were found to be differentially-expressed, clearly distinguishing the two chicken genetic-lines. Several of these differentially-expressed proteins were involved with the proteasome and glycolysis/gluconeogenesis-pathways. Changes in meat-quality traits were also observed, which were reflected in the proteomic-profile. Proteins involved in the ubiquitin-proteasome system were associated with the bigger muscle mass of Ross 508, while phosphoglucomutase 1 was associated with a possible higher capability of AZ Extra Heavy Red-chickens to cope with stressors. This pilot proteomic-approach applied on muscle exudate samples provided key evidence about the pathways and processes underlying these two chicken genetic-lines and their meat-quality parameters. We also identified potential biomarkers that could determine the peculiar production potentials (e.g. breast-growth) of these broilers-lines, which arise from differences in their genetic-backgrounds.

Interneuron diversity in the human dorsal striatum

Deciphering the striatal interneuron diversity is key to understanding the basal ganglia circuit and to untangling the complex neurological and psychiatric diseases affecting this brain structure. We performed snRNA-seq and spatial transcriptomics of postmortem human caudate nucleus and putamen samples to elucidate the diversity and abundance of interneuron populations and their inherent transcriptional structure in the human dorsal striatum. We propose a comprehensive taxonomy of striatal interneurons with eight main classes and fourteen subclasses, providing their full transcriptomic identity and spatial expression profile as well as additional quantitative FISH validation for specific populations. We have also delineated the correspondence of our taxonomy with previous standardized classifications and shown the main transcriptomic and class abundance differences between caudate nucleus and putamen. Notably, based on key functional genes such as ion channels and synaptic receptors, we found matching known mouse interneuron populations for the most abundant populations, the recently described PTHLH and TAC3 interneurons. Finally, we were able to integrate other published datasets with ours, supporting the generalizability of this harmonized taxonomy.

Novel ectopic expression of zona pellucida 3 glycoprotein in lung cancer promotes tumor growth

Zona pellucida 3 (ZP3) expression is classically found in the ZP-layer of the oocytes, lately shown in ovarian and prostate cancer. A successful ZP3 ovarian cancer immunotherapy in transgenic mice suggested its use as an attractive therapeutic target. The biological role of ZP3 in cancer growth and progression is still unknown. We found that ~88% of the analyzed adenocarcinoma, squamous and small cell lung carcinomas to express ZP3. Knockout of ZP3 in a ZP3-expressing lung adenocarcinoma cell line, significantly decreased cell viability, proliferation, and migration rates in vitro. Zona pellucida 3 knock out (ZP3-KO) cell tumors inoculated in vivo in immunodeficient non-obese diabetic, severe combined immunodeficient mice showed significant inhibition of tumor growth and mitigation of the malignant phenotype. RNA sequencing revealed the deregulation of cell migration/adhesion signaling pathways in ZP3-KO cells. This novel functional relevance of ZP3 in lung cancer emphasized the suitability of ZP3 as a target in cancer immunotherapy and as a potential cancer biomarker.

Inhibiting retinoic acid signaling in dendritic cells suppresses respiratory syncytial virus infection through enhanced antiviral immunity

Retinoic acid (RA), controls the immunoregulatory functions of many immune cells, including dendritic cells (DCs), and is important for mucosal immunity. In DCs, RA regulates the expression of pattern recognition receptors and stimulates interferon production. Here, we investigated the role of RA in DCs in mounting immunity to respiratory syncytial virus (RSV). To abolish RA signaling in DCs, we used mice expressing a dominant negative form of retinoic acid receptor-α (RARα) under the CD11c promoter (CD11c-dnRARα). Paradoxically, upon RSV challenge, these animals had lower viral burden, reduced pathology, and greater Th1 polarized immunity than wild-type (WT) mice. Moreover, CD11c-dnRARα DCs infected with RSV showed enhancement in innate and adaptive immunity genes, while genes associated with viral replication were downregulated. These findings suggest that the absence of RA signaling in DCs enhances innate immunity against RSV infection leading to decreased viral load and reduced pathogenicity.

Systems modeling of oncogenic G-protein and GPCR signaling reveals unexpected differences in downstream pathway activation

Mathematical models of biochemical reaction networks are an important and emerging tool for the study of cell signaling networks involved in disease processes. One promising potential application of such mathematical models is the study of how disease-causing mutations promote the signaling phenotype that contributes to the disease. It is commonly assumed that one must have a thorough characterization of the network readily available for mathematical modeling to be useful, but we hypothesized that mathematical modeling could be useful when there is incomplete knowledge and that it could be a tool for discovery that opens new areas for further exploration. In the present study, we first develop a mechanistic mathematical model of a G-protein coupled receptor signaling network that is mutated in almost all cases of uveal melanoma and use model-driven explorations to uncover and explore multiple new areas for investigating this disease. Modeling the two major, mutually-exclusive, oncogenic mutations (Gαq/11 and CysLT2R) revealed the potential for previously unknown qualitative differences between seemingly interchangeable disease-promoting mutations, and our experiments confirmed oncogenic CysLT2R was impaired at activating the FAK/YAP/TAZ pathway relative to Gαq/11. This led us to hypothesize that CYSLTR2 mutations in UM must co-occur with other mutations to activate FAK/YAP/TAZ signaling, and our bioinformatic analysis uncovers a role for co-occurring mutations involving the plexin/semaphorin pathway, which has been shown capable of activating this pathway. Overall, this work highlights the power of mechanism-based computational systems biology as a discovery tool that can leverage available information to open new research areas.

In situ Patch-seq analysis of microglia reveals a lack of stress genes as found in FACS-isolated microglia

We applied the patch-seq technique to harvest transcripts from individual microglial cells from cortex, hippocampus and corpus callosum of acute brain slices from adult mice. After recording membrane currents with the patch-clamp technique, the cytoplasm was collected via the pipette and underwent adapted SMART-seq2 preparation with subsequent sequencing. On average, 4138 genes were detected in 113 cells from hippocampus, corpus callosum and cortex, including microglia markers such as Tmem119, P2ry12 and Siglec-H. Comparing our dataset to previously published single cell mRNA sequencing data from FACS-isolated microglia indicated that two clusters of cells were absent in our patch-seq dataset. Pathway analysis of marker genes in FACS-specific clusters revealed association with microglial activation and stress response. This indicates that under normal conditions microglia in situ lack transcripts associated with a stress-response, and that the microglia-isolation procedure by mechanical dissociation and FACS triggers the expression of genes related to activation and stress.

DrugRepoBank: a comprehensive database and discovery platform for accelerating drug repositioning

In recent years, drug repositioning has emerged as a promising alternative to the time-consuming, expensive and risky process of developing new drugs for diseases. However, the current database for drug repositioning faces several issues, including insufficient data volume, restricted data types, algorithm inaccuracies resulting from the neglect of multidimensional or heterogeneous data, a lack of systematic organization of literature data associated with drug repositioning, limited analytical capabilities and user-unfriendly webpage interfaces. Hence, we have established the first all-encompassing database called DrugRepoBank, consisting of two main modules: the 'Literature' module and the 'Prediction' module. The 'Literature' module serves as the largest repository of literature-supported drug repositioning data with experimental evidence, encompassing 169 repositioned drugs from 134 articles from 1 January 2000 to 1 July 2023. The 'Prediction' module employs 18 efficient algorithms, including similarity-based, artificial-intelligence-based, signature-based and network-based methods to predict repositioned drug candidates. The DrugRepoBank features an interactive and user-friendly web interface and offers comprehensive functionalities such as bioinformatics analysis of disease signatures. When users provide information about a drug, target or disease of interest, DrugRepoBank offers new indications and targets for the drug, proposes new drugs that bind to the target or suggests potential drugs for the queried disease. Additionally, it provides basic information about drugs, targets or diseases, along with supporting literature. We utilize three case studies to demonstrate the feasibility and effectiveness of predictively repositioned drugs within DrugRepoBank. The establishment of the DrugRepoBank database will significantly accelerate the pace of drug repositioning. Database URL:  https://awi.cuhk.edu.cn/DrugRepoBank.

Comprehensive genome-wide identification of Snf2 gene family and their expression profile under salt stress in six Brassica species of U's triangle model

MAIN CONCLUSION

We comprehensively identified and analyzed the Snf2 gene family. Some Snf2 genes were involved in responding to salt stress based on the RNA-seq and qRT-PCR analysis. Sucrose nonfermenting 2 (Snf2) proteins are core components of chromatin remodeling complexes that not only alter DNA accessibility using the energy of ATP hydrolysis, but also play a critical regulatory role in growth, development, and stress response in eukaryotes. However, the comparative study of Snf2 gene family in the six Brassica species in U's triangle model remains unclear. Here, a total of 405 Snf2 genes were identified, comprising 53, 50, and 46 in the diploid progenitors: Brassica rapa (AA, 2n = 20), Brassica nigra (BB, 2n = 16), and Brassica oleracea (CC, 2n = 18), and 93, 91, and 72 in the allotetraploid: Brassica juncea (AABB, 2n = 36), Brassica napus (AACC, 2n = 38), and Brassica carinata (BBCC, 2n = 34), respectively. These genes were classified into six clades and further divided into 18 subfamilies based on their conserved motifs and domains. Intriguingly, these genes showed highly conserved chromosomal distributions and gene structures, indicating that few dynamic changes occurred during the polyploidization. The duplication modes of the six Brassica species were diverse, and the expansion of most Snf2 in Brassica occurred primarily through dispersed duplication (DSD) events. Additionally, the majority of Snf2 genes were under purifying selection during polyploidization, and some Snf2 genes were associated with various abiotic stresses. Both RNA-seq and qRT-PCR analysis showed that the expression of BnaSnf2 genes was significantly induced under salt stress, implying their involvement in salt tolerance response in Brassica species. The results provide a comprehensive understanding of the Snf2 genes in U's triangle model species, which will facilitate further functional analysis of the Snf2 genes in Brassica plants.

Distinct roles of H3K27me3 and H3K36me3 in vernalization response, maintenance, and resetting in winter wheat

Winter plants rely on vernalization, a crucial process for adapting to cold conditions and ensuring successful reproduction. However, understanding the role of histone modifications in guiding the vernalization process in winter wheat remains limited. In this study, we investigated the transcriptome and chromatin dynamics in the shoot apex throughout the life cycle of winter wheat in the field. Two core histone modifications, H3K27me3 and H3K36me3, exhibited opposite patterns on the key vernalization gene VERNALIZATION1 (VRN1), correlating with its induction during cold exposure. Moreover, the H3K36me3 level remained high at VRN1 after cold exposure, which may maintain its active state. Mutations in FERTILIZATION-INDEPENDENT ENDOSPERM (TaFIE) and SET DOMAIN GROUP 8/EARLY FLOWERING IN SHORT DAYS (TaSDG8/TaEFS), components of the writer complex for H3K27me3 and H3K36me3, respectively, affected flowering time. Intriguingly, VRN1 lost its high expression after the cold exposure memory in the absence of H3K36me3. During embryo development, VRN1 was silenced with the removal of active histone modifications in both winter and spring wheat, with selective restoration of H3K27me3 in winter wheat. The mutant of Tafie-cr-87, a component of H3K27me3 "writer" complex, did not influence the silence of VRN1 during embryo development, but rather attenuated the cold exposure requirement of winter wheat. Integrating gene expression with H3K27me3 and H3K36me3 patterns identified potential regulators of flowering. This study unveils distinct roles of H3K27me3 and H3K36me3 in controlling vernalization response, maintenance, and resetting in winter wheat.

Rapid evolution of promoters from germline-specifically expressed genes including transposon silencing factors

BACKGROUND

The piRNA pathway in animal gonads functions as an 'RNA-based immune system', serving to silence transposable elements and prevent inheritance of novel invaders. In Drosophila, this pathway relies on three gonad-specific Argonaute proteins (Argonaute-3, Aubergine and Piwi) that associate with 23-28 nucleotide piRNAs, directing the silencing of transposon-derived transcripts. Transposons constitute a primary driver of genome evolution, yet the evolution of piRNA pathway factors has not received in-depth exploration. Specifically, channel nuclear pore proteins, which impact piRNA processing, exhibit regions of rapid evolution in their promoters. Consequently, the question arises whether such a mode of evolution is a general feature of transposon silencing pathways.

RESULTS

By employing genomic analysis of coding and promoter regions within genes that function in transposon silencing in Drosophila, we demonstrate that the promoters of germ cell-specific piRNA factors are undergoing rapid evolution. Our findings indicate that rapid promoter evolution is a common trait among piRNA factors engaged in germline silencing across insect species, potentially contributing to gene expression divergence in closely related taxa. Furthermore, we observe that the promoters of genes exclusively expressed in germ cells generally exhibit rapid evolution, with some divergence in gene expression.

CONCLUSION

Our results suggest that increased germline promoter evolution, in partnership with other factors, could contribute to transposon silencing and evolution of species through differential expression of genes driven by invading transposons.

CSF proteomic profiles of neurodegeneration biomarkers in Alzheimer's disease

INTRODUCTION

We aimed to unravel the underlying pathophysiology of the neurodegeneration (N) markers neurogranin (Ng), neurofilament light (NfL), and hippocampal volume (HCV), in Alzheimer's disease (AD) using cerebrospinal fluid (CSF) proteomics.

METHODS

Individuals without dementia were classified as A+ (CSF amyloid beta [Aβ]42), T+ (CSF phosphorylated tau181), and N+ or N- based on Ng, NfL, or HCV separately. CSF proteomics were generated and compared between groups using analysis of covariance.

RESULTS

Only a few individuals were A+T+Ng-. A+T+Ng+ and A+T+NfL+ showed different proteomic profiles compared to A+T+Ng- and A+T+NfL-, respectively. Both Ng+ and NfL+ were associated with neuroplasticity, though in opposite directions. Compared to A+T+HCV-, A+T+HCV+ showed few proteomic changes, associated with oxidative stress.

DISCUSSION

Different N markers are associated with distinct neurodegenerative processes and should not be equated. N markers may differentially complement disease staging beyond amyloid and tau. Our findings suggest that Ng may not be an optimal N marker, given its low incongruency with tau pathophysiology.

HIGHLIGHTS

In Alzheimer's disease, neurogranin (Ng)+, neurofilament light (NfL)+, and hippocampal volume (HCV)+ showed differential protein expression in cerebrospinal fluid. Ng+ and NfL+ were associated with neuroplasticity, although in opposite directions. HCV+ showed few proteomic changes, related to oxidative stress. Neurodegeneration (N) markers may differentially refine disease staging beyond amyloid and tau. Ng might not be an optimal N marker, as it relates more closely to tau.

Exceptional longevity of mammalian ovarian and oocyte macromolecules throughout the reproductive lifespan

The mechanisms contributing to age-related deterioration of the female reproductive system are complex, but aberrant protein homeostasis is a major contributor. We elucidated the exceptionally stable proteins, structures, and macromolecules that persist in mammalian ovaries and gametes across the reproductive lifespan. Ovaries exhibit localized structural and cell-type specific enrichment of stable macromolecules in both the follicular and extrafollicular environments. Moreover, both ovaries and oocytes harbor a panel of exceptionally long-lived proteins, including cytoskeletal components, mitochondrial, and oocyte-derived proteins. The exceptional persistence of these long-lived molecules might play a critical role in both lifelong maintenance and age-dependent deterioration of reproductive tissues.

Gene Set Summarization Using Large Language Models

Molecular biologists frequently interpret gene lists derived from high-throughput experiments and computational analysis. This is typically done as a statistical enrichment analysis that measures the over- or under-representation of biological function terms associated with genes or their properties, based on curated assertions from a knowledge base (KB) such as the Gene Ontology (GO). Interpreting gene lists can also be framed as a textual summarization task, enabling Large Language Models (LLMs) to use scientific texts directly and avoid reliance on a KB. TALISMAN (Terminological ArtificiaL Intelligence SuMmarization of Annotation and Narratives) uses generative AI to perform gene set function summarization as a complement to standard enrichment analysis. This method can use different sources of gene functional information: (1) structured text derived from curated ontological KB annotations, (2) ontology-free narrative gene summaries, or (3) direct retrieval from the model. We demonstrate that these methods are able to generate plausible and biologically valid summary GO term lists for an input gene set. However, LLM-based approaches are unable to deliver reliable scores or p-values and often return terms that are not statistically significant. Crucially, in our experiments these methods were rarely able to recapitulate the most precise and informative term from standard enrichment analysis. We also observe minor differences depending on prompt input information, with GO term descriptions leading to higher recall but lower precision. However, newer LLM models perform statistically significantly better than the oldest model across all performance metrics, suggesting that future models may lead to further improvements. Overall, the results are nondeterministic, with minor variations in prompt resulting in radically different term lists, true to the stochastic nature of LLMs. Our results show that at this point, LLM-based methods are unsuitable as a replacement for standard term enrichment analysis, however they may provide summarization benefits for implicit knowledge integration across extant but unstandardized knowledge, for large sets of features, and where the amount of information is difficult for humans to process.

Insights into the pathophysiology and response of persistent spinal pain syndrome type 2 to spinal cord stimulation: a human genome-wide association study

BACKGROUND

Spinal cord stimulation (SCS) provides pain relief for some patients with persistent spinal pain syndrome type 2 (PSPS 2), but the precise mechanisms of action and prognostic factors for a favorable pain response remain obscure. This in vivo human genome-wide association study provides some pathophysiological clues.

METHODS

We performed a high-density oligonucleotide microarray analysis of serum obtained from both PSPS 2 cases and pain-free controls who had undergone lower back spinal surgery at the study site. Using multivariate discriminant analysis, we tried to identify different expressions between mRNA transcripts from PSPS 2 patients relative to controls, SCS responders to non-responders, or SCS responders to themselves before starting SCS. Gene ontology enrichment analysis was used to identify the biological processes that best discriminate between the groups of clinical interest.

RESULTS

Thirty PSPS 2 patients, of whom 23 responded to SCS, were evaluated together with 15 pain-free controls. We identified 11 significantly downregulated genes in serum of PSPS 2 patients compared with pain-free controls and two significantly downregulated genes once the SCS response became apparent. All were suggestive of enhanced inflammation, tissue repair mechanisms and proliferative responses among the former. We could not identify any gene differentiating patients who responded to SCS from those who did not respond.

CONCLUSIONS

This study points out various biological processes that may underlie PSPS 2 pain and SCS therapeutic effects, including the modulation of neuroimmune response, inflammation and restorative processes.

Nutri-cereal tissue-specific transcriptome atlas during development: Functional integration of gene expression to identify mineral uptake pathways in little millet (Panicum sumatrense)

Little millet (Panicum sumatrense Roth ex Roem. & Schult.) is an essential minor millet of southeast Asia and Africa's temperate and subtropical regions. The plant is stress-tolerant, has a short life cycle, and has a mineral-rich nutritional profile associated with unique health benefits. We report the developmental gene expression atlas of little millet (genotype JK-8) from ten tissues representing different stages of its life cycle, starting from seed germination and vegetative growth to panicle maturation. The developmental transcriptome atlas led to the identification of 342 827 transcripts. The BUSCO analysis and comparison with the transcriptomes of related species confirm that this study presents high-quality, in-depth coverage of the little millet transcriptome. In addition, the eFP browser generated here has a user-friendly interface, allowing interactive visualizations of tissue-specific gene expression. Using these data, we identified transcripts, the orthologs of which in Arabidopsis and rice are involved in nutrient acquisition, transport, and response pathways. The comparative analysis of the expression levels of these transcripts holds great potential for enhancing the mineral content in crops, particularly zinc and iron, to address the issue of "hidden hunger" and to attain nutritional security, making it a valuable asset for translational research.

Arabidopsis PROTODERMAL FACTOR2 binds lysophosphatidylcholines and transcriptionally regulates phospholipid metabolism

Plant homeodomain leucine zipper IV (HD-Zip IV) transcription factors (TFs) contain an evolutionarily conserved steroidogenic acute regulatory protein (StAR)-related lipid transfer (START) domain. While the START domain is required for TF activity, its presumed role as a lipid sensor is not clear. Here we used tandem affinity purification from Arabidopsis cell cultures to demonstrate that PROTODERMAL FACTOR2 (PDF2), a representative member that controls epidermal differentiation, recruits lysophosphatidylcholines (LysoPCs) in a START-dependent manner. Microscale thermophoresis assays confirmed that a missense mutation in a predicted ligand contact site reduces lysophospholipid binding. We additionally found that PDF2 acts as a transcriptional regulator of phospholipid- and phosphate (Pi) starvation-related genes and binds to a palindromic octamer with consensus to a Pi response element. Phospholipid homeostasis and elongation growth were altered in pdf2 mutants according to Pi availability. Cycloheximide chase experiments revealed a role for START in maintaining protein levels, and Pi starvation resulted in enhanced protein destabilization, suggesting a mechanism by which lipid binding controls TF activity. We propose that the START domain serves as a molecular sensor for membrane phospholipid status in the epidermis. Our data provide insights toward understanding how the lipid metabolome integrates Pi availability with gene expression.

Poglut2/3 double knockout in mice results in neonatal lethality with reduced levels of fibrillin in lung tissues

Fibrillin microfibrils play a critical role in the formation of elastic fibers, tissue/organ development, and cardiopulmonary function. These microfibrils not only provide structural support and flexibility to tissues, but they also regulate growth factor signaling through a plethora of microfibril-binding proteins in the extracellular space. Mutations in fibrillins are associated with human diseases affecting cardiovascular, pulmonary, skeletal, and ocular systems. Fibrillins consist of up to 47 epidermal growth factor-like repeats, of which more than half are modified by protein O-glucosyltransferase 2 (POGLUT2) and/or POGLUT3. Loss of these modifications reduces secretion of N-terminal fibrillin constructs overexpressed in vitro. Here, we investigated the role of POGLUT2 and POGLUT3 in vivo using a Poglut2/3 double knockout (DKO) mouse model. Blocking O-glucosylation caused neonatal death with skeletal, pulmonary, and eye defects reminiscent of fibrillin/elastin mutations. Proteomic analyses of DKO dermal fibroblast medium and extracellular matrix provided evidence that fibrillins were more sensitive to loss of O-glucose compared to other POGLUT2/3 substrates. This conclusion was supported by immunofluorescent analyses of late gestation DKO lungs where FBN levels were reduced and microfibrils appeared fragmented in the pulmonary arteries and veins, bronchioles, and developing saccules. Defects in fibrillin microfibrils likely contributed to impaired elastic fiber formation and histological changes observed in DKO lung blood vessels, bronchioles, and saccules. Collectively, these results highlight the importance of POGLUT2/3-mediated O-glucosylation in vivo and open the possibility that O-glucose modifications on fibrillin influence microfibril assembly and or protein interactions in the ECM environment.

A catalogue of recombination coldspots in interspecific tomato hybrids

Increasing natural resistance and resilience in plants is key for ensuring food security within a changing climate. Breeders improve these traits by crossing cultivars with their wild relatives and introgressing specific alleles through meiotic recombination. However, some genomic regions are devoid of recombination especially in crosses between divergent genomes, limiting the combinations of desirable alleles. Here, we used pooled-pollen sequencing to build a map of recombinant and non-recombinant regions between tomato and five wild relatives commonly used for introgressive tomato breeding. We detected hybrid-specific recombination coldspots that underscore the role of structural variations in modifying recombination patterns and maintaining genetic linkage in interspecific crosses. Crossover regions and coldspots show strong association with specific TE superfamilies exhibiting differentially accessible chromatin between somatic and meiotic cells. About two-thirds of the genome are conserved coldspots, located mostly in the pericentromeres and enriched with retrotransposons. The coldspots also harbor genes associated with agronomic traits and stress resistance, revealing undesired consequences of linkage drag and possible barriers to breeding. We presented examples of linkage drag that can potentially be resolved by pairing tomato with other wild species. Overall, this catalogue will help breeders better understand crossover localization and make informed decisions on generating new tomato varieties.

Tolerance to Haemophilus influenzae infection in human epithelial cells: Insights from a primary cell-based model

Haemophilus influenzae is a human respiratory pathogen and inhabits the human respiratory tract as its only niche. Despite this, the molecular mechanisms that allow H. influenzae to establish persistent infections of human epithelia are not well understood. Here, we have investigated how H. influenzae adapts to the host environment and triggers the host immune response using a human primary cell-based infection model that closely resembles human nasal epithelia (NHNE). Physiological assays combined with dualRNAseq revealed that NHNE from five healthy donors all responded to H. influenzae infection with an initial, 'unproductive' inflammatory response that included a strong hypoxia signature but did not produce pro-inflammatory cytokines. Subsequently, an apparent tolerance to large extracellular and intraepithelial burdens of H. influenzae developed, with NHNE transcriptional profiles resembling the pre-infection state. This occurred in parallel with the development of intraepithelial bacterial populations, and appears to involve interruption of NFκB signalling. This is the first time that large-scale, persistence-promoting immunomodulatory effects of H. influenzae during infection have been observed, and we were able to demonstrate that only infections with live, but not heat-killed H. influenzae led to immunomodulation and reduced expression of NFκB-controlled cytokines such as IL-1β, IL-36γ and TNFα. Interestingly, NHNE were able to re-activate pro-inflammatory responses towards the end of the 14-day infection, resulting in release of IL-8 and TNFα. In addition to providing first molecular insights into mechanisms enabling persistence of H. influenzae in the host, our data further indicate the presence of infection stage-specific gene expression modules, highlighting fundamental similarities between immune responses in NHNE and canonical immune cells, which merit further investigation.

Nobiletin regulates intracellular Ca2+ levels via IP3R and ameliorates neuroinflammation in Aβ42-induced astrocytes

Astrocytes are the major glial cells in the human brain and provide crucial metabolic and trophic support to neurons. The amyloid-β peptide (Aβ) alter the morphological and functional properties of astrocytes and induce inflammation and calcium dysregulation, contributing to Alzheimer's disease (AD) pathology. Recent studies highlight the role of Toll-like receptor (TLR) 4/nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) signaling in inflammation. Reactive oxygen species (ROS) generated due to Aβ, induce apoptosis in the brain cells worsening AD progression. Astrocytic cell surface receptors, such as purinergic receptors (P2Y1 and P2Y2), metabotropic glutamate receptor (mGLUR)5, α7 nicotinic acetylcholine receptor (α7nAChR), and N-methyl-d-aspartate receptors (NMDARs), have been suggested to interact with inositol trisphosphate receptor (IP3R) on the endoplasmic reticulum (ER) to induce Ca2+ movement from ER to cytoplasm, causing Ca2+ dysregulation. We found that the citrus flavonoid nobiletin (NOB) protected primary astrocytes from Aβ42-induced cytotoxicity and inhibited TLR4/NF-κB signaling in Aβ42-induced primary rat astrocytes. NOB was found to regulate Aβ42-induced ROS levels through Keap1-Nrf2 pathway. The receptors P2Y1, P2Y2, mGLUR5, α7nAChR, and NMDARs induced intracellular Ca2+ levels by activating IP3R and NOB regulated them, thereby regulating intracellular Ca2+ levels. Molecular docking analysis revealed a possible interaction between NOB and IP3R in IP3R regulation. Furthermore, RNA sequencing revealed various NOB-mediated biological signaling pathways, such as the AD-presenilin, AD-amyloid secretase, and Wnt signaling pathway, suggesting possible neuroprotective roles of NOB. To conclude, NOB is a promising therapeutic agent for AD and works by modulating AD pathology at various levels in Aβ42-induced primary rat astrocytes.

Ribosome rescue factor PELOTA modulates translation start site choice for C/EBPα protein isoforms

Translation initiation at alternative start sites can dynamically control the synthesis of two or more functionally distinct protein isoforms from a single mRNA. Alternate isoforms of the developmental transcription factor CCAAT/enhancer-binding protein α (C/EBPα) produced from different start sites exert opposing effects during myeloid cell development. This choice between alternative start sites depends on sequence features of the CEBPA transcript, including a regulatory uORF, but the molecular basis is not fully understood. Here, we identify the factors that affect C/EBPα isoform choice using a sensitive and quantitative two-color fluorescent reporter coupled with CRISPRi screening. Our screen uncovered a role of the ribosome rescue factor PELOTA (PELO) in promoting the expression of the longer C/EBPα isoform by directly removing inhibitory unrecycled ribosomes and through indirect effects mediated by the mechanistic target of rapamycin kinase. Our work uncovers further links between ribosome recycling and translation reinitiation that regulate a key transcription factor, with implications for normal hematopoiesis and leukemogenesis.

Exploring the diversity of galls on Artemisia indica induced by Rhopalomyia species through morphological and transcriptome analyses

Plant galls generated by insects have highly organized structures, providing nutrients and shelter to the insects living within them. Most research on the physiological and molecular mechanisms of gall development has focused on single galls. To understand the diversity of gall development, we examined five galls with different morphologies generated by distinct species of Rhopalomyia (gall midge; Diptera: Cecidomyiidae) on a single host plant of Artemisia indica var. maximowiczii (Asteraceae). Vasculature developed de novo within the galls, indicating active transport of nutrients between galls and the host plant. Each gall had a different pattern of vasculature and lignification, probably due to differences in the site of gall generation and the gall midge species. Transcriptome analysis indicated that photosynthetic and cell wall-related genes were down-regulated in leaf and stem galls, respectively, compared with control leaf and stem tissues, whereas genes involved in floral organ development were up-regulated in all types of galls, indicating that transformation from source to sink organs occurs during gall development. Our results help to understand the diversity of galls on a single herbaceous host plant.

The reference genome sequence of Artemisia argyi provides insights into secondary metabolism biosynthesis

Artemisia argyi, a perennial herb of the genus Artemisia in the family Asteraceae, holds significant importance in Chinese traditional medicine, referred to as "Aicao". Here, we report a high-quality reference genome of Artemisia argyi L. cv. beiai, with a genome size up to 4.15 Gb and a contig N50 of 508.96 Kb, produced with third-generation Nanopore sequencing technology. We predicted 147,248 protein-coding genes, with approximately 68.86% of the assembled sequences comprising repetitive elements, primarily long terminal repeat retrotransposons(LTRs). Comparative genomics analysis shows that A. argyi has the highest number of specific gene families with 5121, and much more families with four or more members than the other 6 plant species, which is consistent with its more expanded gene families and fewer contracted gene families. Furthermore, through transcriptome sequencing of A. argyi in response to exogenous MeJA treatment, we have elucidated acquired regulatory insights into MeJA's impact on the phenylpropanoid, flavonoid, and terpenoid biosynthesis pathways of A. argyi. The whole-genome information obtained in this study serves as a valuable resource for delving deeper into the cultivation and molecular breeding of A. argyi. Moreover, it holds promise for enhancing genome assemblies across other members of the Asteraceae family. The identification of key genes establishes a solid groundwork for developing new varieties of Artemisia with elevated concentrations of active compounds.

FDX1 as a novel biomarker and treatment target for stomach adenocarcinoma

BACKGROUND

Stomach adenocarcinoma (STAD) is one of the main reasons for cancer-related deaths worldwide. This investigation aimed to define the connection between STAD and Cuproptosis-related genes (CRGs). Cuproptosis is a newly identified form of mitochondrial cell death triggered by copper.

AIM

To explore the identification of potential biomarkers for STAD disease based on cuproptosis.

METHODS

A predictive model using Gene Ontology (GO), Least Absolute Shrinkage and Selection Operator (LASSO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Set Variation Analysis (GSVA), and Gene Set Enrichment Analysis analyzed gene interconnections, focusing on 3 copper-related genes and their expression in The Cancer Genome Atlas-STAD. Networks for mRNA-miRNA and mRNA-transcription factor interactions were constructed. The prognostic significance of CRG scores was evaluated using time-receiver operating characteristic, Kaplan-Meier curves, and COX regression analysis. Validation was conducted with datasets GSE26942, GSE54129, and GSE66229. Expression of copper-related differentially expressed genes was also analyzed in various human tissues and gastric cancer subpopulations using the human protein atlas.

RESULTS

Three significant genes (FDX1, LIAS, MTF1) were identified and selected via LASSO analysis to predict and classify individuals with STAD into high and low CRG score subgroups. These genes were down-regulated in both risk categories. GO and KEGG analyses highlighted their involvement mainly in the electron transport chain. After validating their differential expression, FDX1 emerged as the most accurate diagnostic marker for gastric cancer. Additionally, the RCircos package localized FDX1 on chromosome 11.

CONCLUSION

Our study revealed that FDX1 could be a potential biomarker and treatment target for gastric malignancy, providing new ideas for further scientific research.

Quantitative proteomics analysis reveals possible anticancer mechanisms of 5'-deoxy-5'-methylthioadenosine in cholangiocarcinoma cells

Cholangiocarcinoma (CCA) is an aggressive cancer originating from bile duct epithelium, particularly prevalent in Asian countries with liver fluke infections. Current chemotherapy for CCA often fails due to drug resistance, necessitating novel anticancer agents. This study investigates the potential of 5'-deoxy-5'-methylthioadenosine (MTA), a naturally occurring nucleoside, against CCA. While MTA has shown promise against various cancers, its effects on CCA remain unexplored. We evaluated MTA's anticancer activity in CCA cell lines and drug-resistant sub-lines, assessing cell viability, migration, invasion, and apoptosis. The potential anticancer mechanisms of MTA were explored through proteomic analysis using LC-MS/MS and bioinformatic analysis. The results show a dose-dependent reduction in CCA cell viability, with enhanced effects on cancer cells compared to normal cells. Moreover, MTA inhibits growth, induces apoptosis, and suppresses cell migration and invasion. Additionally, MTA enhanced the anticancer effects of gemcitabine on drug-resistant CCA cells. Proteomics revealed the down-regulation of multiple proteins by MTA, affecting various molecular functions, biological processes, and cellular components. Network analysis highlighted MTA's role in inhibiting proteins related to mitochondrial function and energy derivation, crucial for cell growth and survival. Additionally, MTA suppressed proteins involved in cell morphology and cytoskeleton organization, important for cancer cell motility and metastasis. Six candidate genes, including ZNF860, KLC1, GRAMD1C, MAMSTR, TANC1, and TTC13, were selected from the top 10 most down-regulated proteins identified in the proteomics results and were subsequently verified through RT-qPCR. Further, KLC1 protein suppression by MTA treatment was confirmed through Western blotting. Additionally, based on TCGA data, KLC1 mRNA was found to be upregulated in the tissue of CCA patients compared to that of normal adjacent tissues. In summary, MTA shows promising anticancer potential against CCA by inhibiting growth, inducing apoptosis, and suppressing migration and invasion, while enhancing gemcitabine's effects. Proteomic analysis elucidates possible molecular mechanisms underlying MTA's anticancer activity, laying the groundwork for future research and development of MTA as a treatment for advanced CCA.

scRNA-seq reveals the diversity of the developing cardiac cell lineage and molecular players in heart rhythm regulation

We utilized scRNA-seq to delineate the diversity of cell types in the zebrafish heart. Transcriptome profiling of over 50,000 cells at 48 and 72 hpf defined at least 18 discrete cell lineages of the developing heart. Utilizing well-established gene signatures, we identified a population of cells likely to be the primary pacemaker and characterized the transcriptome profile defining this critical cell type. Two previously uncharacterized genes, atp1b3b and colec10, were found to be enriched in the sinoatrial cardiomyocytes. CRISPR/Cas9-mediated knockout of these two genes significantly reduced heart rate, implicating their role in cardiac development and conduction. Additionally, we describe other cardiac cell lineages, including the endothelial and neural cells, providing their expression profiles as a resource. Our results established a detailed atlas of the developing heart, providing valuable insights into cellular and molecular mechanisms, and pinpointed potential new players in heart rhythm regulation.

Engineering an inducible leukemia-associated fusion protein enables large-scale ex vivo production of functional human phagocytes

Ex vivo expansion of human CD34+ hematopoietic stem and progenitor cells remains a challenge due to rapid differentiation after detachment from the bone marrow niche. In this study, we assessed the capacity of an inducible fusion protein to enable sustained ex vivo proliferation of hematopoietic precursors and their capacity to differentiate into functional phagocytes. We fused the coding sequences of an FK506-Binding Protein 12 (FKBP12)-derived destabilization domain (DD) to the myeloid/lymphoid lineage leukemia/eleven nineteen leukemia (MLL-ENL) fusion gene to generate the fusion protein DD-MLL-ENL and retrovirally expressed the protein switch in human CD34+ progenitors. Using Shield1, a chemical inhibitor of DD fusion protein degradation, we established large-scale and long-term expansion of late monocytic precursors. Upon Shield1 removal, the cells lost self-renewal capacity and spontaneously differentiated, even after 2.5 y of continuous ex vivo expansion. In the absence of Shield1, stimulation with IFN-γ, LPS, and GM-CSF triggered terminal differentiation. Gene expression analysis of the obtained phagocytes revealed marked similarity with naïve monocytes. In functional assays, the novel phagocytes migrated toward CCL2, attached to VCAM-1 under shear stress, produced reactive oxygen species, and engulfed bacterial particles, cellular particles, and apoptotic cells. Finally, we demonstrated Fcγ receptor recognition and phagocytosis of opsonized lymphoma cells in an antibody-dependent manner. Overall, we have established an engineered protein that, as a single factor, is useful for large-scale ex vivo production of human phagocytes. Such adjustable proteins have the potential to be applied as molecular tools to produce functional immune cells for experimental cell-based approaches.

A library of reporters of the global regulators of gene expression in Escherichia coli

The topology of the transcription factor network (TFN) of Escherichia coli is far from uniform, with 22 global regulator (GR) proteins controlling one-third of all genes. So far, their production rates cannot be tracked by comparable fluorescent proteins. We developed a library of fluorescent reporters for 16 GRs for this purpose. Each consists of a single-copy plasmid coding for green fluorescent protein (GFP) fused to the full-length copy of the native promoter. We tracked their activity in exponential and stationary growth, as well as under weak and strong stresses. We show that the reporters have high sensitivity and specificity to all stresses tested and detect single-cell variability in transcription rates. Given the influence of GRs on the TFN, we expect that the new library will contribute to dissecting global transcriptional stress-response programs of E. coli. Moreover, the library can be invaluable in bioindustrial applications that tune those programs to, instead of cell growth, favor productivity while reducing energy consumption.IMPORTANCECells contain thousands of genes. Many genes are involved in the control of cellular activities. Some activities require a few hundred genes to run largely synchronous transcriptional programs. To achieve this, cells have evolved global regulator (GR) proteins that can influence hundreds of genes simultaneously. We have engineered a library of Escherichia coli strains to track the levels over time of these, phenotypically critical, GRs. Each strain has a single-copy plasmid coding for a fast-maturing green fluorescent protein whose transcription is controlled by a copy of the natural GR promoter. By allowing the tracking of GR levels, with sensitivity and specificity, this library should become of wide use in scientific research on bacterial gene expression (from molecular to synthetic biology) and, later, be used in applications in therapeutics and bioindustries.

Secretome profiling of human epithelial cells exposed to cigarette smoke extract and their effect on human lung microvascular endothelial cells

Cigarette smoke (CS) is one of the leading causes of pulmonary diseases and can induce lung secretome alteration. CS exposure-induced damages to human pulmonary epithelial cells and microvascular endothelial cells have been extensively demonstrated; however, the effects of the secretome of lung epithelial cells exposed to CS extracts (CSE) on lung microvascular endothelial cells are not fully understood. In this study, we aimed to determine the effects of the secretome of lung epithelial cells exposed to CSE on lung microvascular endothelial cells. Human lung epithelial cells, A549, were exposed to CSE, and the secretome was collected. Human lung microvascular endothelial cells, HULEC-5a, were used to evaluate the effect of the secretome of A549 exposed to CSE. Secretome profile, endothelial cell death, inflammation, and permeability markers were determined. CSE altered the secretome expression of A549 cells, and secretome derived from CSE-exposed A549 cells caused respiratory endothelial cell death, inflammation, and moderately enhanced endothelial permeability. This study demonstrates the potential role of cellular interaction between endothelial and epithelial cells during exposure to CSE and provides novel therapeutic targets or beneficial biomarkers using secretome analysis for CSE-related respiratory diseases.

Novel Serum Proteomes Expressed from Benzene Exposure Among Gasoline Station Attendants

Background

Research on the proteomes impact of benzene exposure in fuel station employees remains sparse, underscoring the need for detailed health impact assessments focusing on biomarker evaluation.

Objectives

This investigation aimed to analyze the differences in blood parameters and serum proteomes resulting from benzene exposure between gasoline station attendants (B-GSA) and a control group.

Design and methods

A cross-sectional analytical study was conducted with 96 participants, comprising 54 in the B-GSA group and 42 in the control group. The methodology employed included an interview questionnaire alongside urine and blood sample collections. The urine samples were analyzed for trans,trans-muconic acid (t,t-MA) levels, while the blood samples underwent complete blood count analysis and proteome profiling.

Results

Post-shift analysis indicated that the B-GSA group exhibited significantly higher levels of t,t-MA and monocytes compared to the control group (P < .05). Proteome quantification identified 1448 proteins differentially expressed between the B-GSA and control groups. Among these, 20 proteins correlated with the levels of t,t-MA in urine. Notably, 4 proteins demonstrated more than a 2-fold down-regulation in the B-GSA group: HBS1-like, non-structural maintenance of chromosomes element 1 homolog, proprotein convertase subtilisin/kexin type 4, and zinc finger protein 658. The KEGG pathway analysis revealed associations with apoptosis, cancer pathways, p53 signaling, and the TNF signaling pathway.

Conclusion

The changes in these 4 significant proteins may elucidate the molecular mechanisms underlying benzene toxicity and suggest their potential as biomarkers for benzene poisoning in future assessments.

Antibacterial Activities of Phenolic Compounds in Miang Extract: Growth Inhibition and Change in Protein Expression of Extensively Drug-Resistant Klebsiella pneumoniae

The rising incidence of extensively drug-resistant (XDR) Klebsiella pneumoniae, including carbapenem- and colistin-resistant strains, leads to the limitation of available effective antibiotics. Miang, known as chewing tea, is produced from Camellia sinensis var. assamica or Assam tea leaves fermentation. Previous studies revealed that the extract of Miang contains various phenolic and flavonoid compounds with numerous biological activities including antibacterial activity. However, the antibacterial activity of Miang against XDR bacteria especially colistin-resistant strains had not been investigated. In this study, the compositions of phenolic and flavonoid compounds in fresh, steamed, and fermented Assam tea leaves were examined by HPLC, and their antibacterial activities were evaluated by the determination of the MIC and MBC. Pyrogallol was detected only in the extract from Miang and showed the highest activities with an MIC of 0.25 mg/mL and an MBC of 0.25-0.5 mg/mL against methicillin-susceptible Staphylococcus aureus, methicillin-resistant S. aureus, Escherichia coli ATCC 25922, colistin-resistant E. coli, and colistin-resistant K. pneumoniae. The effects on morphology and proteomic changes in K. pneumoniae NH54 treated with Miang extract were characterized by SEM and label-free quantitative shotgun proteomics analysis. The results revealed that Miang extract caused the decrease in bacterial cell wall integrity and cell lysis. The up- and downregulated expression with approximately a 2 to >5-fold change in proteins involved in peptidoglycan synthesis and outer membrane, carbohydrate, and amino acid metabolism were identified. These findings suggested that Miang containing pyrogallol and other secondary metabolites from fermentation has potential as an alternative candidate with an antibacterial agent or natural active pharmaceutical ingredient against XDR bacteria including colistin-resistant bacteria.

Pre-hypertrophic chondrogenic enhancer landscape of limb and axial skeleton development

Chondrocyte differentiation controls skeleton development and stature. Here we provide a comprehensive map of chondrocyte-specific enhancers and show that they provide a mechanistic framework through which non-coding genetic variants can influence skeletal development and human stature. Working with fetal chondrocytes isolated from mice bearing a Col2a1 fluorescent regulatory sensor, we identify 780 genes and 2'704 putative enhancers specifically active in chondrocytes using a combination of RNA-seq, ATAC-seq and H3K27ac ChIP-seq. Most of these enhancers (74%) show pan-chondrogenic activity, with smaller populations being restricted to limb (18%) or trunk (8%) chondrocytes only. Notably, genetic variations overlapping these enhancers better explain height differences than those overlapping non-chondrogenic enhancers. Finally, targeted deletions of identified enhancers at the Fgfr3, Col2a1, Hhip and, Nkx3-2 loci confirm their role in regulating cognate genes. This enhancer map provides a framework for understanding how genes and non-coding variations influence bone development and diseases.

Neurodevelopmental and synaptic defects in DNAJC6 parkinsonism, amenable to gene therapy

DNAJC6 encodes auxilin, a co-chaperone protein involved in clathrin-mediated endocytosis (CME) at the presynaptic terminal. Biallelic mutations in DNAJC6 cause a complex, early-onset neurodegenerative disorder characterized by rapidly progressive parkinsonism-dystonia in childhood. The disease is commonly associated with additional neurodevelopmental, neurological and neuropsychiatric features. Currently, there are no disease-modifying treatments for this condition, resulting in significant morbidity and risk of premature mortality. To investigate the underlying disease mechanisms in childhood-onset DNAJC6 parkinsonism, we generated induced pluripotent stem cells (iPSC) from three patients harbouring pathogenic loss-of-function DNAJC6 mutations and subsequently developed a midbrain dopaminergic neuronal model of disease. When compared to age-matched and CRISPR-corrected isogenic controls, the neuronal cell model revealed disease-specific auxilin deficiency as well as disturbance of synaptic vesicle recycling and homeostasis. We also observed neurodevelopmental dysregulation affecting ventral midbrain patterning and neuronal maturation. To explore the feasibility of a viral vector-mediated gene therapy approach, iPSC-derived neuronal cultures were treated with lentiviral DNAJC6 gene transfer, which restored auxilin expression and rescued CME. Our patient-derived neuronal model provides deeper insights into the molecular mechanisms of auxilin deficiency as well as a robust platform for the development of targeted precision therapy approaches.

SNARE chaperone Sly1 directly mediates close-range vesicle tethering

The essential Golgi protein Sly1 is a member of the Sec1/mammalian Unc-18 (SM) family of SNARE chaperones. Sly1 was originally identified through remarkable gain-of-function alleles that bypass requirements for diverse vesicle tethering factors. Employing genetic analyses and chemically defined reconstitutions of ER-Golgi fusion, we discovered that a loop conserved among Sly1 family members is not only autoinhibitory but also acts as a positive effector. An amphipathic lipid packing sensor (ALPS)-like helix within the loop directly binds high-curvature membranes. Membrane binding is required for relief of Sly1 autoinhibition and also allows Sly1 to directly tether incoming vesicles to the Qa-SNARE on the target organelle. The SLY1-20 mutation bypasses requirements for diverse tethering factors but loses this ability if the tethering activity is impaired. We propose that long-range tethers, including Golgins and multisubunit tethering complexes, hand off vesicles to Sly1, which then tethers at close range to initiate trans-SNARE complex assembly and fusion in the early secretory pathway.

Transcriptomic profile in carbendazim-induced developmental defects in zebrafish (Danio rerio) embryos/larvae

Carbendazim is a widely used fungicide to protect agricultural and horticultural crops against a wide array of fungal species. Published reports have shown that the wide usage of carbendazim resulted in reprotoxicity, carcinogenicity, immunotoxicity, and developmental toxicity in mammalian models. However, studies related to the developmental toxicity of carbendazim in aquatic organisms are not clear. To address this gap, an attempt was made by exposing zebrafish embryos to carbendazim (800 μg/L) and assessing the phenotypic and transcriptomic profile at different developmental stages [24 hour post fertilization (hpf), 48 hpf, 72 hpf and 96 hpf). At 48 hpf, phenotypic abnormalities such as delay in hatching rate, deformed spinal axial curvature, and pericardial edema were observed in zebrafish larvae over its respective controls. At 72 hpf, exposure of zebrafish embryos exposed to carbendazim resulted in scoliosis; however, unexposed larvae did not exhibit signs of scoliosis. Interestingly, the transcriptomic analysis revealed a total of 1253 DEGs were observed at selected time points, while unique genes at 24 hpf, 48 hpf, 72 hpf and 96 hpf was found to be 76.54 %, 61.14 %, 92.98 %, and 68.28 %, respectively. Functional profiling of downregulated genes revealed altered transcriptomic markers associated with phototransduction (24 hpf and 72 hpf), immune system (48 hpf), and SNARE interactions in the vesicular pathway (96 hpf). Whereas functional profiling of upregulated genes revealed altered transcriptomic markers associated with riboflavin metabolism (24 hpf), basal transcription factors (48 hpf), insulin signaling pathway (72 hpf), and primary bile acid biosynthesis (96 hpf). Taken together, carbendazim-induced developmental toxicity could be ascribed to pleiotropic responses at the molecular level, which in turn might reflect phenotypic abnormalities.

MulStack: An ensemble learning prediction model of multilabel mRNA subcellular localization

Subcellular localization of mRNA is related to protein synthesis, cell polarity, cell movement and other biological regulation mechanisms. The distribution of mRNAs in subcellulars is similar to that of proteins, and most mRNAs are distributed in multiple subcellulars. Recently, some computational methods have been designed to predict the subcellular localization of mRNA. However, these methods only employed a sin-gle level of mRNA features and did not employ the position encoding of nucleotides in mRNA. In this paper, an ensemble learning prediction model is proposed, named MulStack, which is based on random forest and deep learning for multilabel mRNA subcellular localization. The proposed method employs two levels of mRNA features, including sequence-level and residue-level features, and position encoding is employed for the first time in the field of subcellular localization of mRNA. Random forest is employed to learn mRNA sequence-level feature, deep learning is employed to learn mRNA sequence-level feature and mRNA residue-level combined with position encoding. And the outputs of random forest and deep learning model will be weighted sum as the prediction probability. Compared with existing methods, the results show that MulStack is the best in the localization of the nucleus, cytosol and exosome. In addition, position weight matrices (PWMs) are extracted by convolutional neural networks (CNNs) that can be matched with known RNA binding protein motifs. Gene ontology (GO) enrichment analysis shows biological processes, molecular functions and cellular components of mRNA genes. The prediction web server of MulStack is freely accessible at http://bliulab.net/MulStack.

DNA Electrotransfer Regulates Molecular Functions in Skeletal Muscle

Background

Tissues, such as skeletal muscle, have been targeted for the delivery of plasmid DNA (pDNA) encoding vaccines and therapeutics. The application of electric pulses (electroporation or electrotransfer) increases cell membrane permeability to enhance plasmid delivery and expression. However, the molecular effects of DNA electrotransfer on the muscle tissue are poorly characterized.

Materials and Methods

Four hours after intramuscular plasmid electrotransfer, we evaluated gene expression changes by RNA sequencing. Differentially expressed genes were analyzed by gene ontology (GO) pathway enrichment analysis.

Results

GO analysis highlighted many enriched molecular functions. The terms regulated by pulse application were related to muscle stress, the cytoskeleton and inflammation. The terms regulated by pDNA injection were related to a DNA-directed response and its control. Several terms regulated by pDNA electrotransfer were similar to those regulated by pulse application. However, the terms related to pDNA injection differed, focusing on entry of the plasmid into the cells and intracellular trafficking.

Conclusion

Each muscle stimulus resulted in specific regulated molecular functions. Identifying the unique intrinsic molecular changes driven by intramuscular DNA electrotransfer will aid in the design of preventative and therapeutic gene therapies.

TBX3 is essential for establishment of the posterior boundary of anterior genes and upregulation of posterior genes together with HAND2 during the onset of limb bud development

During limb bud formation, axis polarities are established as evidenced by the spatially restricted expression of key regulator genes. In particular, the mutually antagonistic interaction between the GLI3 repressor and HAND2 results in distinct and non-overlapping anterior-distal Gli3 and posterior Hand2 expression domains. This is a hallmark of the establishment of antero-posterior limb axis polarity, together with spatially restricted expression of homeodomain and other transcriptional regulators. Here, we show that TBX3 is required for establishment of the posterior expression boundary of anterior genes in mouse limb buds. ChIP-seq and differential gene expression analysis of wild-type and mutant limb buds identifies TBX3-specific and shared TBX3-HAND2 target genes. High sensitivity fluorescent whole-mount in situ hybridisation shows that the posterior expression boundaries of anterior genes are positioned by TBX3-mediated repression, which excludes anterior genes such as Gli3, Alx4, Hand1 and Irx3/5 from the posterior limb bud mesenchyme. This exclusion delineates the posterior mesenchymal territory competent to establish the Shh-expressing limb bud organiser. In turn, HAND2 is required for Shh activation and cooperates with TBX3 to upregulate shared posterior identity target genes in early limb buds.

The influence of biological sex in human skeletal muscle transcriptome during ageing

Sex is a crucial biological variable, and influence of biological sex on the change of gene expression in ageing skeletal muscle has not yet been fully revealed. In this study, the mRNA expression profiles were obtained from the Gene Expression Omnibus database. Key genes were identified by differential expression analysis and weighted gene co-expression network analysis. The gene set enrichment analysis software and Molecular Signatures Database were used for functional and enrichment analysis. A protein-protein interaction network was constructed using STRING and visualized in Cytoscape. The results were compared between female and male subgroups. Differentially expressed genes and enriched pathways in different sex subgroups shared only limited similarities. The pathways enriched in the female subgroup were more similar to the pathways enriched in the older groups without taking sex difference into consideration. The pathways enriched in the female subgroup were more similar to the pathways enriched in the older groups without taking sex difference into consideration. The muscle myosin filament pathways were downregulated in the both aged female and male samples whereas transforming growth factor beta pathway and extracellular matrix-related pathways were upregulated. With muscle ageing, the metabolism-related pathways, protein synthesis and degradation pathways, results of predicted immune cell infiltration, and gene cluster associated with slow-type myofibers drastically different between the female and male subgroups. This finding may indicate that changes in muscle type with ageing may differ between the sexes in vastus lateralis muscle.

Physiochemical interaction between osmotic stress and a bacterial exometabolite promotes plant disease

Various microbes isolated from healthy plants are detrimental under laboratory conditions, indicating the existence of molecular mechanisms preventing disease in nature. Here, we demonstrated that application of sodium chloride (NaCl) in natural and gnotobiotic soil systems is sufficient to induce plant disease caused by an otherwise non-pathogenic root-derived Pseudomonas brassicacearum isolate (R401). Disease caused by combinatorial treatment of NaCl and R401 triggered extensive, root-specific transcriptional reprogramming that did not involve down-regulation of host innate immune genes, nor dampening of ROS-mediated immunity. Instead, we identified and structurally characterized the R401 lipopeptide brassicapeptin A as necessary and sufficient to promote disease on salt-treated plants. Brassicapeptin A production is salt-inducible, promotes root colonization and transitions R401 from being beneficial to being detrimental on salt-treated plants by disturbing host ion homeostasis, thereby bolstering susceptibility to osmolytes. We conclude that the interaction between a global change stressor and a single exometabolite from a member of the root microbiome promotes plant disease in complex soil systems.

Proteomics and metabolic burden analysis to understand the impact of recombinant protein production in E. coli

The impact of recombinant protein production (RPP) on host cells and the metabolic burden associated with it undermine the efficiency of the production system. This study utilized proteomics to investigate the dynamics of parent and recombinant cells induced at different time points for RPP. The results revealed significant changes in both transcriptional and translational machinery that may have impacted the metabolic burden, growth rate of the culture and the RPP. The timing of protein synthesis induction also played a critical role in the fate of the recombinant protein within the host cell, affecting protein and product yield. The study identified significant differences in the expression of proteins involved in fatty acid and lipid biosynthesis pathways between two E. coli host strains (M15 and DH5⍺), with the E. coli M15 strain demonstrating superior expression characteristics for the recombinant protein. Overall, these findings contribute to the knowledge base for rational strain engineering for optimized recombinant protein production.