Transfer learning enables identification of multiple types of RNA modifications using nanopore direct RNA sequencing

Nanopore direct RNA sequencing (DRS) has emerged as a powerful tool for RNA modification identification. However, concurrently detecting multiple types of modifications in a single DRS sample remains a challenge. Here, we develop TandemMod, a transferable deep learning framework capable of detecting multiple types of RNA modifications in single DRS data. To train high-performance TandemMod models, we generate in vitro epitranscriptome datasets from cDNA libraries, containing thousands of transcripts labeled with various types of RNA modifications. We validate the performance of TandemMod on both in vitro transcripts and in vivo human cell lines, confirming its high accuracy for profiling m6A and m5C modification sites. Furthermore, we perform transfer learning for identifying other modifications such as m7G, Ψ, and inosine, significantly reducing training data size and running time without compromising performance. Finally, we apply TandemMod to identify 3 types of RNA modifications in rice grown in different environments, demonstrating its applicability across species and conditions. In summary, we provide a resource with ground-truth labels that can serve as benchmark datasets for nanopore-based modification identification methods, and TandemMod for identifying diverse RNA modifications using a single DRS sample.

Molecular Diagnostics of Cryptococcus spp. and Immunomics of Cryptococcosis-Associated Immune Reconstitution Inflammatory Syndrome

Cryptococcal infection poses a significant global public health challenge, particularly in regions near the equator. In this review, we offer a succinct exploration of the Cryptococcus spp. genome and various molecular typing methods to assess the burden and genetic diversity of cryptococcal pathogens in the environment and clinical isolates. We delve into a detailed discussion on the molecular pathogenesis and diagnosis of immune reconstitution inflammatory syndrome (IRIS) associated with cryptococcosis, with a specific emphasis on cryptococcal meningitis IRIS (CM-IRIS). Our examination includes the recent literature on CM-IRIS, covering host cellulomics, proteomics, transcriptomics, and genomics.

Transcriptomic and metabolomic correlation analysis: effect of initial SO2 addition on higher alcohol synthesis in Saccharomyces cerevisiae and identification of key regulatory genes

Introduction

Higher alcohols are volatile compounds produced during alcoholic fermentation that affect the quality and safety of the final product. This study used a correlation analysis of transcriptomics and metabolomics to study the impact of the initial addition of SO2 (30, 60, and 90 mg/L) on the synthesis of higher alcohols in Saccharomyces cerevisiae EC1118a and to identify key genes and metabolic pathways involved in their metabolism.

Methods

Transcriptomics and metabolomics correlation analyses were performed and differentially expressed genes (DEGs) and differential metabolites were identified. Single-gene knockouts for targeting genes of important pathways were generated to study the roles of key genes involved in the regulation of higher alcohol production.

Results

We found that, as the SO2 concentration increased, the production of total higher alcohols showed an overall trend of first increasing and then decreasing. Multi-omics correlation analysis revealed that the addition of SO2 affected carbon metabolism (ko01200), pyruvate metabolism (ko00620), glycolysis/gluconeogenesis (ko00010), the pentose phosphate pathway (ko00030), and other metabolic pathways, thereby changing the precursor substances. The availability of SO2 indirectly affects the formation of higher alcohols. In addition, excessive SO2 affected the growth of the strain, leading to the emergence of a lag phase. We screened the ten most likely genes and constructed recombinant strains to evaluate the impact of each gene on the formation of higher alcohols. The results showed that ADH4, SER33, and GDH2 are important genes of alcohol metabolism in S. cerevisiae. The isoamyl alcohol content of the EC1118a-ADH4 strain decreased by 21.003%; The isobutanol content of the EC1118a-SER33 strain was reduced by 71.346%; and the 2-phenylethanol content of EC1118a-GDH2 strain was reduced by 25.198%.

Conclusion

This study lays a theoretical foundation for investigating the mechanism of initial addition of SO2 in the synthesis of higher alcohols in S. cerevisiae, uncovering DEGs and key metabolic pathways related to the synthesis of higher alcohols, and provides guidance for regulating these mechanisms.

Transcriptomic analysis of the effect of remote ischaemic conditioning in an animal model of necrotising enterocolitis

Necrotising enterocolitis (NEC) has a complex pathophysiology but the common end-point is ischaemia reperfusion injury (IRI) and intestinal necrosis. We have previously reported that RIC significantly reduces the intestinal injury in a rat model of NEC. Here we describe the changes in intestinal mRNA occurring in the intestine of animals exposed to IRI, both with and without RIC. Related rat-pups were randomly assigned to four groups: SHAM, IRI only, RIC only and RIC + IRI. IRI animals, underwent 40 min of intestinal ischaemia, and 90 min of reperfusion. Animals that underwent RIC had three cycles of 5 min of alternating ischaemia/reperfusion by means of a ligature applied to the hind limb. Samples from the terminal ileum were immediately stored in RNA-preserving media for later next generation sequencing and transciptome analysis using R v 3.6.1. Differential expression testing showed that 868 genes differentially expressed in animals exposed to RIC alone compared to SHAM and 135 in the IRI and RIC group compared to IRI alone. Comparison between these two sets showed that 25 genes were differentially expressed in both groups. Pro-inflammatory molecules: NF-ĸβ2, Cxcl1, SOD2 and Map3k8 all show reduced expression in response to RIC. Targeted gene analysis revealed increased expression in PI3K which is part of the so-called RISK-pathway which is a key part of the protective mechanisms of RIC in the heart. Overall, this transcriptomic analysis shows that RIC provides a protective effect to the intestine via anti-inflammatory pathways. This could be particularly relevant to treating and preventing NEC.

Integrated Metabolomic and Transcriptomic Analysis Reveals the Underlying Antibacterial Mechanisms of the Phytonutrient Quercetin-Induced Fatty Acids Alteration in Staphylococcus aureus ATCC 27217

The utilization of natural products in food preservation represents a promising strategy for the dual benefits of controlling foodborne pathogens and enhancing the nutritional properties of foods. Among the phytonutrients, flavonoids have been shown to exert antibacterial effects by disrupting bacterial cell membrane functionality; however, the underlying molecular mechanisms remain elusive. In this study, we investigated the effect of quercetin on the cell membrane permeability of Staphylococcus aureus ATCC 27217. A combined metabolomic and transcriptomic approach was adopted to examine the regulatory mechanism of quercetin with respect to the fatty acid composition and associated genes. Kinetic analysis and molecular docking simulations were conducted to assess quercetin's inhibition of β-ketoacyl-acyl carrier protein reductase (FabG), a potential target in the bacterial fatty acid biosynthesis pathway. Metabolomic and transcriptomic results showed that quercetin increased the ratio of unsaturated to saturated fatty acids and the levels of membrane phospholipids. The bacteria reacted to quercetin-induced stress by attempting to enhance fatty acid biosynthesis; however, quercetin directly inhibited FabG activity, thereby disrupting bacterial fatty acid biosynthesis. These findings provide new insights into the mechanism of quercetin's effects on bacterial cell membranes and suggest potential applications for quercetin in bacterial inhibition.

Molecular Mechanism of Exogenous ABA to Enhance UV-B Resistance in Rhododendron chrysanthum Pall. by Modulating Flavonoid Accumulation

With the depletion of the ozone layer, the intensity of ultraviolet B (UV-B) radiation reaching the Earth's surface increases, which in turn causes significant stress to plants and affects all aspects of plant growth and development. The aim of this study was to investigate the mechanism of response to UV-B radiation in the endemic species of Rhododendron chrysanthum Pall. (R. chrysanthum) in the Changbai Mountains and to study how exogenous ABA regulates the response of R. chrysanthum to UV-B stress. The results of chlorophyll fluorescence images and OJIP kinetic curves showed that UV-B radiation damaged the PSII photosystem of R. chrysanthum, and exogenous ABA could alleviate this damage to some extent. A total of 2148 metabolites were detected by metabolomics, of which flavonoids accounted for the highest number (487, or 22.67%). KEGG enrichment analysis of flavonoids that showed differential accumulation by UV-B radiation and exogenous ABA revealed that flavonoid biosynthesis and flavone and flavonol biosynthesis were significantly altered. GO analysis showed that most of the DEGs produced after UV-B radiation and exogenous ABA were distributed in the cellular process, cellular anatomical entity, and catalytic activity. Network analysis of key DFs and DEGs associated with flavonoid synthesis identified key flavonoids (isorhamnetin-3-O-gallate and dihydromyricetin) and genes (TRINITY_DN2213_c0_g1_i4-A1) that promote the resistance of R. chrysanthum to UV-B stress. In addition, multiple transcription factor families were found to be involved in the regulation of the flavonoid synthesis pathway under UV-B stress. Overall, R. chrysanthum actively responded to UV-B stress by regulating changes in flavonoids, especially flavones and flavonols, while exogenous ABA further enhanced its resistance to UV-B stress. The experimental results not only provide a new perspective for understanding the molecular mechanism of the response to UV-B stress in the R. chrysanthum, but also provide a valuable theoretical basis for future research and application in improving plant adversity tolerance.

Diversity and host interaction of the gut microbiota in specific pathogen-free pigs

Pigs are widely used as animal models in various studies related to humans. The interaction between the gut microbiota and the host has significant effects on the host's health and disease status. However, although there have been many studies investigating the pig gut microbiota, the findings have been inconsistent due to variations in rearing conditions. Interactions between the gut microbiota and host have not been fully explored in pigs. Specific pathogen-free (SPF) pigs are ideal non-primate large animals to study the interactions between the gut microbiota and the host. In this study, we performed high-throughput sequencing analysis of the gut microbiota and the gut tissue transcriptome of six SPF pigs to provide a systematic understanding of the composition, function, and spatial distribution of gut microbiota in SPF pigs. We identified significant differences in microbial diversity and functionality among different gastrointestinal tract sites. Metagenomics data analysis revealed significant differences in alpha diversity and beta diversity of microbiota in different gastrointestinal sites of SPF pigs. Additionally, transcriptomic data indicated significant differences in gene expression as well as KEGG and GO functional enrichment between the small intestine and large intestine. Furthermore, by combining microbial metagenomics and host transcriptomics analyses, specific correlations were found between gut microbiota and host genes. These included a negative correlation between the TCN1 gene and Prevotella dentalis, possibly related to bacterial metabolic pathways involving vitamin B12, and a positive correlation between the BDH1 gene and Roseburia hominis, possibly because both are involved in fatty acid metabolism. These findings lay the groundwork for further exploration of the co-evolution between the microbiota and the host, specifically in relation to nutrition, metabolism, and immunity. In conclusion, we have elucidated the diversity of the gut microbiota in SPF pigs and conducted a detailed investigation into the interactions between the gut microbiota and host gene expression. These results contribute to our understanding of the intricate dynamics between the gut microbiota and the host, offering important references for advancements in life science research, bioproduct production, and sustainable development in animal husbandry.

Pre-vaccination transcriptomic profiles of immune responders to the MUC1 peptide vaccine for colon cancer prevention

Self-antigens abnormally expressed on tumors, such as MUC1, have been targeted by therapeutic cancer vaccines. We recently assessed in two clinical trials in a preventative setting whether immunity induced with a MUC1 peptide vaccine could reduce high colon cancer risk in individuals with a history of premalignant colon adenomas. In both trials, there were immune responders and non-responders to the vaccine. Here we used PBMC pre-vaccination and 2 weeks after the first vaccine of responders and non-responders selected from both trials to identify early biomarkers of immune response involved in long-term memory generation and prevention of adenoma recurrence. We performed flow cytometry, phosflow, and differential gene expression analyses on PBMCs collected from MUC1 vaccine responders and non-responders pre-vaccination and two weeks after the first of three vaccine doses. MUC1 vaccine responders had higher frequencies of CD4 cells pre-vaccination, increased expression of CD40L on CD8 and CD4 T-cells, and a greater increase in ICOS expression on CD8 T-cells. Differential gene expression analysis revealed that iCOSL, PI3K AKT MTOR, and B-cell signaling pathways are activated early in response to the MUC1 vaccine. We identified six specific transcripts involved in elevated antigen presentation, B-cell activation, and NF-kB1 activation that were directly linked to finding antibody response at week 12. Finally, a model using these transcripts was able to predict non-responders with accuracy. These findings suggest that individuals who can be predicted to respond to the MUC1 vaccine, and potentially other vaccines, have greater readiness in all immune compartments to present and respond to antigens. Predictive biomarkers of MUC1 vaccine response may lead to more effective vaccines tailored to individuals with high risk for cancer but with varying immune fitness.

A common polymorphism in the Intelectin-1 gene influences mucus plugging in severe asthma

By incompletely understood mechanisms, type 2 (T2) inflammation present in the airways of severe asthmatics drives the formation of pathologic mucus which leads to airway mucus plugging. Here we investigate the molecular role and clinical significance of intelectin-1 (ITLN-1) in the development of pathologic airway mucus in asthma. Through analyses of human airway epithelial cells we find that ITLN1 gene expression is highly induced by interleukin-13 (IL-13) in a subset of metaplastic MUC5AC+ mucus secretory cells, and that ITLN-1 protein is a secreted component of IL-13-induced mucus. Additionally, we find ITLN-1 protein binds the C-terminus of the MUC5AC mucin and that its deletion in airway epithelial cells partially reverses IL-13-induced mucostasis. Through analysis of nasal airway epithelial brushings, we find that ITLN1 is highly expressed in T2-high asthmatics, when compared to T2-low children. Furthermore, we demonstrate that both ITLN-1 gene expression and protein levels are significantly reduced by a common genetic variant that is associated with protection from the formation of mucus plugs in T2-high asthma. This work identifies an important biomarker and targetable pathways for the treatment of mucus obstruction in asthma.

Molecular and Physiological Effects of 17α-methyltestosterone on Sex Differentiation of Black Rockfish, Sebastes schlegelii

It is widely known that all-female fish production holds economic value for aquaculture. Sebastes schlegelii, a preeminent economic species, exhibits a sex dimorphism, with females surpassing males in growth. In this regard, achieving all-female black rockfish production could significantly enhance breeding profitability. In this study, we utilized the widely used male sex-regulating hormone, 17α-methyltestosterone (MT) at three different concentrations (20, 40, and 60 ppm), to produce pseudomales of S. schlegelii for subsequent all-female offspring breeding. Long-term MT administration severely inhibits the growth of S. schlegelii, while short term had no significant impact. Histological analysis confirmed sex reversal at all MT concentrations; however, both medium and higher MT concentrations impaired testis development. MT also influenced sex steroid hormone levels in pseudomales, suppressing E2 while increasing T and 11-KT levels. In addition, a transcriptome analysis revealed that MT down-regulated ovarian-related genes (cyp19a1a and foxl2) while up-regulating male-related genes (amh) in pseudomales. Furthermore, MT modulated the TGF-β signaling and steroid hormone biosynthesis pathways, indicating its crucial role in S. schlegelii sex differentiation. Therefore, the current study provides a method for achieving sexual reversal using MT in S. schlegelii and offers an initial insight into the underlying mechanism of sexual reversal in this species.

Prediction of m6A and m5C at single-molecule resolution reveals a transcriptome-wide co-occurrence of RNA modifications

The epitranscriptome embodies many new and largely unexplored functions of RNA. A significant roadblock hindering progress in epitranscriptomics is the identification of more than one modification in individual transcript molecules. We address this with CHEUI (CH3 (methylation) Estimation Using Ionic current). CHEUI predicts N6-methyladenosine (m6A) and 5-methylcytosine (m5C) in individual molecules from the same sample, the stoichiometry at transcript reference sites, and differential methylation between any two conditions. CHEUI processes observed and expected nanopore direct RNA sequencing signals to achieve high single-molecule, transcript-site, and stoichiometry accuracies in multiple tests using synthetic RNA standards and cell line data. CHEUI's capability to identify two modification types in the same sample reveals a co-occurrence of m6A and m5C in individual mRNAs in cell line and tissue transcriptomes. CHEUI provides new avenues to discover and study the function of the epitranscriptome.

Transcriptome dynamics of Gossypium purpurascens in response to abiotic stresses by Iso-seq and RNA-seq data

Gossypium purpurascens is a member of the Malvaceae family, holds immense economic significance as a fiber crop worldwide. Abiotic stresses harm cotton crops, reduce yields, and cause economic losses. Generating high-quality reference genomes and large-scale transcriptomic datasets across diverse conditions can offer valuable insights into identifying preferred agronomic traits for crop breeding. The present research used leaf tissues to conduct PacBio Iso-seq and RNA-seq analysis. We carried out an in-depth analysis of DEGs using both correlations with cluster analysis and principal component analysis. Additionally, the study also involved the identification of both lncRNAs and CDS. We have prepared RNA-seq libraries from 75 RNA samples to study the effects of drought, salinity, alkali, and saline-alkali stress, as well as control conditions. A total of 454.06 Gigabytes of transcriptome data were effectively validated through the identification of differentially expressed genes and KEGG and GO analysis. Overwhelmingly, gene expression profiles and full-length transcripts from cotton tissues will aid in understanding the genetic mechanism of abiotic stress tolerance in G. purpurascens.

Transcriptomics reveals age-related changes in ion transport-related factors in yak lungs

Yaks inhabit high-altitude, low-oxygen regions, where ion transport functions play a crucial role in maintaining intracellular and extracellular ionic balance and regulating pulmonary vascular tension. These functions affect pulmonary ventilation and blood flow rate, aiding tissue development and enhancing oxygen transfer efficiency, thus facilitating better adaptation to hypoxic environments. To investigate the regulatory mechanisms of ion transport-related factors on the growth and development of yak lungs, we employed RNA sequencing (RNA-seq)for sequencing the transcriptome in the lung tissues of neonatal (1-day-old), juvenile (1-year-old), and adult (4-year-old) yaks. We also performed differential gene expression and functional analyses. The results yielded 26 genes associated with ion transport, mainly enriched in the salivary and pancreatic secretion pathways. Finally, we used several methods including quantitative polymerase chain reaction (qRT-PCR), and Western blotting (WB), immunohistochemical (IHC) and immunofluorescence (IF) staining to determine the distribution of the expression of the ion transport genes FOXI1, KCNMA1, and SLC12A2 in yak lung tissues. qRT-PCR and WB results indicated that mRNA and protein relative expression levels of FOXI1 and SLC12A2 were significantly higher in neonatal yaks than in juvenile and adult yaks (all p < 0.05), whereas those of KCNMA1 were significantly higher in adult yaks than in neonatal and juvenile yaks (all p < 0.05). IHC and IF results demonstrated that FOXI1, KCNMA1, and SLC12A2 were distributed among the epithelial mucosal layers (including ciliated, goblet, and Clara cells) of the yaks' bronchi and their branches in the lungs across different age groups of yak. Therefore, our results suggested that FOXI1, KCNMA1, and SLC12A2 may be strongly associated with the development and aging processes in yak lungs. These results provide insights into the molecular mechanisms underlying the yak's adaptation to high-altitude environments and valuable references for further research.

The landscape of cancer-rewired GPCR signaling axes

We explored the dysregulation of G-protein-coupled receptor (GPCR) ligand systems in cancer transcriptomics datasets to uncover new therapeutics opportunities in oncology. We derived an interaction network of receptors with ligands and their biosynthetic enzymes. Multiple GPCRs are differentially regulated together with their upstream partners across cancer subtypes and are associated to specific transcriptional programs and to patient survival patterns. The expression of both receptor-ligand (or enzymes) partners improved patient stratification, suggesting a synergistic role for the activation of GPCR networks in modulating cancer phenotypes. Remarkably, we identified many such axes across several cancer molecular subtypes, including many involving receptor-biosynthetic enzymes for neurotransmitters. We found that GPCRs from these actionable axes, including, e.g., muscarinic, adenosine, 5-hydroxytryptamine, and chemokine receptors, are the targets of multiple drugs displaying anti-growth effects in large-scale, cancer cell drug screens, which we further validated. We have made the results generated in this study freely available through a webapp (gpcrcanceraxes.bioinfolab.sns.it).

Multi-Omics Characterization of Colon Mucosa and Submucosa/Wall from Crohn's Disease Patients

Crohn's disease (CD) is a subtype of inflammatory bowel disease (IBD) characterized by transmural disease. The concept of transmural healing (TH) has been proposed as an indicator of deep clinical remission of CD and as a predictor of favorable treatment endpoints. Understanding the pathophysiology involved in transmural disease is critical to achieving these endpoints. However, most studies have focused on the intestinal mucosa, overlooking the contribution of the intestinal wall in Crohn's disease. Multi-omics approaches have provided new avenues for exploring the pathogenesis of Crohn's disease and identifying potential biomarkers. We aimed to use transcriptomic and proteomic technologies to compare immune and mesenchymal cell profiles and pathways in the mucosal and submucosa/wall compartments to better understand chronic refractory disease elements to achieve transmural healing. The results revealed similarities and differences in gene and protein expression profiles, metabolic mechanisms, and immune and non-immune pathways between these two compartments. Additionally, the identification of protein isoforms highlights the complex molecular mechanisms underlying this disease, such as decreased RTN4 isoforms (RTN4B2 and RTN4C) in the submucosa/wall, which may be related to the dysregulation of enteric neural processes. These findings have the potential to inform the development of novel therapeutic strategies to achieve TH.

Spatially revealed roles for lncRNAs in Drosophila spermatogenesis, Y chromosome function and evolution

Unlike coding genes, the number of lncRNA genes in organism genomes is relatively proportional to organism complexity. From plants to humans, the tissues with highest numbers and levels of lncRNA gene expression are the male reproductive organs. To learn why, we initiated a genome-wide analysis of Drosophila lncRNA spatial expression patterns in these tissues. The numbers of genes and levels of expression observed greatly exceed those previously reported, due largely to a preponderance of non-polyadenylated transcripts. In stark contrast to coding genes, the highest numbers of lncRNAs expressed are in post-meiotic spermatids. Correlations between expression levels, localization and previously performed genetic analyses indicate high levels of function and requirement. More focused analyses indicate that lncRNAs play major roles in evolution by controlling transposable element activities, Y chromosome gene expression and sperm construction. A new type of lncRNA-based particle found in seminal fluid may also contribute to reproductive outcomes.

tauFisher predicts circadian time from a single sample of bulk and single-cell pseudobulk transcriptomic data

As the circadian clock regulates fundamental biological processes, disrupted clocks are often observed in patients and diseased tissues. Determining the circadian time of the patient or the tissue of focus is essential in circadian medicine and research. Here we present tauFisher, a computational pipeline that accurately predicts circadian time from a single transcriptomic sample by finding correlations between rhythmic genes within the sample. We demonstrate tauFisher's performance in adding timestamps to both bulk and single-cell transcriptomic samples collected from multiple tissue types and experimental settings. Application of tauFisher at a cell-type level in a single-cell RNAseq dataset collected from mouse dermal skin implies that greater circadian phase heterogeneity may explain the dampened rhythm of collective core clock gene expression in dermal immune cells compared to dermal fibroblasts. Given its robustness and generalizability across assay platforms, experimental setups, and tissue types, as well as its potential application in single-cell RNAseq data analysis, tauFisher is a promising tool that facilitates circadian medicine and research.

The transcription factor OCT6 promotes the dissolution of the naïve pluripotent state by repressing Nanog and activating a formative state gene regulatory network

In the mouse embryo, the transition from the preimplantation to the postimplantation epiblast is governed by changes in the gene regulatory network (GRN) that lead to transcriptional, epigenetic, and functional changes. This transition can be faithfully recapitulated in vitro by the differentiation of mouse embryonic stem cells (mESCs) to epiblast-like cells (EpiLCs), that reside in naïve and formative states of pluripotency, respectively. However, the GRN that drives this conversion is not fully elucidated. Here we demonstrate that the transcription factor OCT6 is a key driver of this process. Firstly, we show that Oct6 is not expressed in mESCs but is rapidly induced as cells exit the naïve pluripotent state. By deleting Oct6 in mESCs, we find that knockout cells fail to acquire the typical morphological changes associated with the formative state when induced to differentiate. Additionally, the key naïve pluripotency TFs Nanog, Klf2, Nr5a2, Prdm14, and Esrrb were expressed at higher levels than in wild-type cells, indicating an incomplete dismantling of the naïve pluripotency GRN. Conversely, premature expression of Oct6 in naïve cells triggered a rapid morphological transformation mirroring differentiation, that was accompanied by the upregulation of the endogenous Oct6 as well as the formative genes Sox3, Zic2/3, Foxp1, Dnmt3A and FGF5. Strikingly, we found that OCT6 represses Nanog in a bistable manner and that this regulation is at the transcriptional level. Moreover, our findings also reveal that Oct6 is repressed by NANOG. Collectively, our results establish OCT6 as a key TF in the dissolution of the naïve pluripotent state and support a model where Oct6 and Nanog form a double negative feedback loop which could act as an important toggle mediating the transition to the formative state.

Comparative transcriptomics of seed nourishing tissues: uncovering conserved and divergent pathways in seed plants

The evolutionary and ecological success of spermatophytes is intrinsically linked to the seed habit, which provides a protective environment for the initial development of the new generation. This environment includes an ephemeral nourishing tissue that supports embryo growth. In gymnosperms this tissue originates from the asexual proliferation of the maternal megagametophyte, while in angiosperms it is a product of fertilization, and is called the endosperm. The emergence of these nourishing tissues is of profound evolutionary value, and they are also food staples for most of the world's population. Here, using Orthofinder to infer orthologue genes among newly generated and previously published datasets, we provide a comparative transcriptomic analysis of seed nourishing tissues from species of several angiosperm clades, including those of early diverging lineages, as well as of one gymnosperm. Our results show that, although the structure and composition of seed nourishing tissues has seen significant divergence along evolution, there are signatures that are conserved throughout the phylogeny. Conversely, we identified processes that are specific to species within the clades studied, and thus illustrate their functional divergence. With this, we aimed to provide a foundation for future studies on the evolutionary history of seed nourishing structures, as well as a resource for gene discovery in future functional studies.

PfHDAC1 is an essential regulator of P. falciparum asexual proliferation and host cell invasion genes with a dynamic genomic occupancy responsive to artemisinin stress

Plasmodium falciparum, the deadly protozoan parasite responsible for malaria, has a tightly regulated gene expression profile closely linked to its intraerythrocytic development cycle. Epigenetic modifiers of the histone acetylation code have been identified as key regulators of the parasite's transcriptome but require further investigation. In this study, we map the genomic distribution of Plasmodium falciparum histone deacetylase 1 (PfHDAC1) across the erythrocytic asexual development cycle and find it has a dynamic occupancy over a wide array of developmentally relevant genes. Overexpression of PfHDAC1 results in a progressive increment in parasite load over consecutive rounds of the asexual infection cycle and is associated with enhanced gene expression of multiple families of host cell invasion factors (merozoite surface proteins, rhoptry proteins, etc.) and with increased merozoite invasion efficiency. With the use of class-specific inhibitors, we demonstrate that PfHDAC1 activity in parasites is crucial for timely intraerythrocytic development. Interestingly, overexpression of PfHDAC1 results in decreased sensitivity to frontline-drug dihydroartemisinin in parasites. Furthermore, we identify that artemisinin exposure can interfere with PfHDAC1 abundance and chromatin occupancy, resulting in enrichment over genes implicated in response/resistance to artemisinin. Finally, we identify that dihydroartemisinin exposure can interrupt the in vitro catalytic deacetylase activity and post-translational phosphorylation of PfHDAC1, aspects that are crucial for its genomic function. Collectively, our results demonstrate PfHDAC1 to be a regulator of critical functions in asexual parasite development and host invasion, which is responsive to artemisinin exposure stress and deterministic of resistance to it.

IMPORTANCE

Malaria is a major public health problem, with the parasite Plasmodium falciparum causing most of the malaria-associated mortality. It is spread by the bite of infected mosquitoes and results in symptoms such as cyclic fever, chills, and headache. However, if left untreated, it can quickly progress to a more severe and life-threatening form. The World Health Organization currently recommends the use of artemisinin combination therapy, and it has worked as a gold standard for many years. Unfortunately, certain countries in southeast Asia and Africa, burdened with a high prevalence of malaria, have reported cases of drug-resistant infections. One of the major problems in controlling malaria is the emergence of artemisinin resistance. Population genomic studies have identified mutations in the Kelch13 gene as a molecular marker for artemisinin resistance. However, several reports thereafter indicated that Kelch13 is not the main mediator but rather hinted at transcriptional deregulation as a major determinant of drug resistance. Earlier, we identified PfGCN5 as a global regulator of stress-responsive genes, which are known to play a central role in artemisinin resistance generation. In this study, we have identified PfHDAC1, a histone deacetylase as a cell cycle regulator, playing an important role in artemisinin resistance generation. Taken together, our study identified key transcriptional regulators that play an important role in artemisinin resistance generation.

Dynamic molecular signatures of acute myocardial infarction based on transcriptomics and metabolomics

Acute myocardial infarction (AMI) commonly precedes ventricular remodeling, heart failure. Few dynamic molecular signatures have gained widespread acceptance in mainstream clinical testing despite the discovery of many potential candidates. These unmet needs with respect to biomarker and drug discovery of AMI necessitate a prioritization. We enrolled patients with AMI aged between 30 and 70. RNA-seq analysis was performed on the peripheral blood mononuclear cells collected from the patients at three time points: 1 day, 7 days, and 3 months after AMI. PLC/LC-MS analysis was conducted on the peripheral blood plasma collected from these patients at the same three time points. Differential genes and metabolites between groups were screened by bio-informatics methods to understand the dynamic changes of AMI in different periods. We obtained 15 transcriptional and 95 metabolite expression profiles at three time points after AMI through high-throughput sequencing. AMI-1d: enrichment analysis revealed the biological features of 1 day after AMI primarily included acute inflammatory response, elevated glycerophospholipid metabolism, and decreased protein synthesis capacity. Phosphatidylcholine (PC) and phosphatidylethanolamine (PE) might stand promising biomarkers to differentiate post-AMI stage. Anti-inflammatory therapy during the acute phase is an important direction for preventing related pathology. AMI-7d: the biological features of this stage primarily involved the initiation of cardiac fibrosis response and activation of platelet adhesion pathways. Accompanied by upregulated TGF-beta signaling pathway and ECM receptor interaction, GP5 help assess platelet activation, a potential therapeutic target to improve haemostasis. AMI-3m: the biological features of 3 months after AMI primarily showed a vascular regeneration response with VEGF signaling pathway, NOS3 and SHC2 widely activated, which holds promise for providing new therapeutic approaches for AMI. Our analysis highlights transcriptional and metabolomics signatures at different time points after MI, which deepens our understanding of the dynamic biological responses and associated molecular mechanisms that occur during cardiac repair.

Exploring CDKN1A Upregulation Mechanisms: Insights into Cell Cycle Arrest Induced by NC2603 Curcumin Analog in MCF-7 Breast Cancer Cells

Breast cancer stands out as one of the most prevalent malignancies worldwide, necessitating a nuanced understanding of its molecular underpinnings for effective treatment. Hormone receptors in breast cancer cells substantially influence treatment strategies, dictating therapeutic approaches in clinical settings, serving as a guide for drug development, and aiming to enhance treatment specificity and efficacy. Natural compounds, such as curcumin, offer a diverse array of chemical structures with promising therapeutic potential. Despite curcumin's benefits, challenges like poor solubility and rapid metabolism have spurred the exploration of analogs. Here, we evaluated the efficacy of the curcumin analog NC2603 to induce cell cycle arrest in MCF-7 breast cancer cells and explored its molecular mechanisms. Our findings reveal potent inhibition of cell viability (IC50 = 5.6 μM) and greater specificity than doxorubicin toward MCF-7 vs. non-cancer HaCaT cells. Transcriptome analysis identified 12,055 modulated genes, most notably upregulation of GADD45A and downregulation of ESR1, implicating CDKN1A-mediated regulation of proliferation and cell cycle genes. We hypothesize that the curcumin analog by inducing GADD45A expression and repressing ESR1, triggers the expression of CDKN1A, which in turn downregulates the expression of many important genes of proliferation and the cell cycle. These insights advance our understanding of curcumin analogs' therapeutic potential, highlighting not just their role in treatment, but also the molecular pathways involved in their activity toward breast cancer cells.

Transcriptomic responses to antibiotic exposure in Mycobacterium tuberculosis

Transcriptional responses in bacteria following antibiotic exposure offer insights into antibiotic mechanism of action, bacterial responses, and characterization of antimicrobial resistance. We aimed to define the transcriptional antibiotic response (TAR) in Mycobacterium tuberculosis (Mtb) isolates for clinically relevant drugs by pooling and analyzing Mtb microarray and RNA-seq data sets. We generated 99 antibiotic transcription profiles across 17 antibiotics, with 76% of profiles generated using 3-24 hours of antibiotic exposure and 49% within one doubling of the WHO antibiotic critical concentration. TAR genes were time-dependent, and largely specific to the antibiotic mechanism of action. TAR signatures performed well at predicting antibiotic exposure, with the area under the receiver operating curve (AUC) ranging from 0.84-1.00 (TAR <6 hours of antibiotic exposure) and 0.76-1.00 (>6 hours of antibiotic exposure) for upregulated genes and 0.57-0.90 and 0.87-1.00, respectfully, for downregulated genes. This work desmonstrates that transcriptomics allows for the assessment of antibiotic activity in Mtb within 6 hours of exposure.

Functional versatility of Zur in metal homeostasis, motility, biofilm formation, and stress resistance in Yersinia pseudotuberculosis

Zur (zinc uptake regulator) is a significant member of the Fur (ferric uptake regulator) superfamily, which is widely distributed in bacteria. Zur plays crucial roles in zinc homeostasis and influences cell development and environmental adaptation in various species. Yersinia pseudotuberculosis is a Gram-negative enteric that pathogen usually serves as a model organism in pathogenicity studies. The regulatory effects of Zur on the zinc transporter ZnuABC and the protein secretion system T6SS have been documented in Y. pseudotuberculosis. In this study, a comparative transcriptomics analysis between a ∆zur mutant and the wild-type (WT) strain of Y. pseudotuberculosis was conducted using RNA-seq. This analysis revealed global regulation by Zur across multiple functional categories, including membrane transport, cell motility, and molecular and energy metabolism. Additionally, Zur mediates the homeostasis not only of zinc but also ferric and magnesium in vivo. There was a notable decrease in 35 flagellar biosynthesis and assembly-related genes, leading to reduced swimming motility in the ∆zur mutant strain. Furthermore, Zur upregulated multiple simple sugar and oligopeptide transport system genes by directly binding to their promoters. The absence of Zur inhibited biofilm formation as well as reduced resistance to chloramphenicol and acidic stress. This study illustrates the comprehensive regulatory functions of Zur, emphasizing its importance in stress resistance and pathogenicity in Y. pseudotuberculosis.

IMPORTANCE

Bacteria encounter diverse stresses in the environment and possess essential regulators to modulate the expression of genes in responding to the stresses for better fitness and survival. Zur (zinc uptake regulator) plays a vital role in zinc homeostasis. Studies of Zur from multiple species reviewed that it influences cell development, stress resistance, and virulence of bacteria. Y. pseudotuberculosis is an enteric pathogen that serves a model organism in the study of pathogenicity, virulence factors, and mechanism of environmental adaptation. In this study, transcriptomics analysis of Zur's regulons was conducted in Y. pseudotuberculosis. The functions of Zur as a global regulator in metal homeostasis, motility, nutrient acquisition, glycan metabolism, and nucleotide metabolism, in turn, increasing the biofilm formation, stress resistance, and virulence were reviewed. The importance of Zur in environmental adaptation and pathogenicity of Y. pseudotuberculosis was emphasized.