Distinct 5-methylcytosine profiles in poly(A) RNA from mouse embryonic stem cells and brain

An Epigenome-Wide Association Study of DNA Methylation and Proliferative Retinopathy over 28 Years in Type 1 Diabetes

Purpose

To perform a prospective epigenome-wide association study of DNA methylation (DNAm) and 28-year proliferative diabetic retinopathy (PDR) incidence in type 1 diabetes (T1D).

Design

Prospective observational cohort study.

Participants

The Pittsburgh Epidemiology of Diabetes Complications (EDC) study of childhood-onset (< 17 years) T1D.

Methods

Stereoscopic fundus photographs were taken in fields 1, 2, and 4 at baseline, 2, 4, 6, 8, 16, 23, and 28 years after DNAm measurements. The photos were graded using the modified Airlie House System. In those free of PDR at baseline (n = 265; mean T1D duration of 18 years at baseline), whole blood DNAm (EPIC array) at 683 597 CpGs was analyzed in Cox models for time to event. Associations between significant CpGs and clinical risk factors were assessed; genetic variants associated with DNAm were identified (methylation quantitative trait loci [meQTLs]). Mendelian randomization was used to examine evidence of causal associations between DNAm and PDR. Post hoc regional and functional analyses were performed.

Main Outcome Measures

Proliferative diabetic retinopathy was defined as the first instance of a grade of ≥ 60 in at least 1 eye or pan-retinal photocoagulation for PDR. Follow-up time was calculated from the study visit at which DNAm data were available (baseline) until PDR incidence or censoring (December 31, 2018 or last follow-up).

Results

PDR incidence was 53% over 28-years' follow-up. Greater DNAm of cg27512687 (KIF16B) was associated with reduced PDR incidence (P = 6.3 × 10-9; false discovery rate [FDR]: < 0.01); 113 cis-meQTLs (P < 5 × 10-8) were identified. Mendelian randomization analysis using the sentinel meQTL as the instrumental variable supported a potentially causal association between cg27512687 and PDR. Cg27512687 was also associated with lower pulse rate and albumin excretion rate and higher estimated glomerular filtration rate, but its association with PDR remained independently significant after adjustment for those factors. In regional analyses, DNAm of FUT4, FKBP1A, and RIN2 was also associated with PDR incidence.

Conclusions

DNA methylation of KIF16B, FUT4, FKBP1A, and RIN2 was associated with PDR incidence, supporting roles for epigenetic regulation of iron clearance, developmental pathways, and autophagy in PDR pathogenesis. Further study of those loci may provide insight into novel targets for interventions to prevent or delay PDR in T1D.

Financial Disclosures

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Splice_sim: a nucleotide conversion-enabled RNA-seq simulation and evaluation framework

Nucleotide conversion RNA sequencing techniques interrogate chemical RNA modifications in cellular transcripts, resulting in mismatch-containing reads. Biases in mapping the resulting reads to reference genomes remain poorly understood. We present splice_sim, a splice-aware RNA-seq simulation and evaluation pipeline that introduces user-defined nucleotide conversions at set frequencies, creates mixture models of converted and unconverted reads, and calculates mapping accuracies per genomic annotation. By simulating nucleotide conversion RNA-seq datasets under realistic experimental conditions, including metabolic RNA labeling and RNA bisulfite sequencing, we measure mapping accuracies of state-of-the-art spliced-read mappers for mouse and human transcripts and derive strategies to prevent biases in the data interpretation.

Regulations of m6A and other RNA modifications and their roles in cancer

RNA modification is an essential component of the epitranscriptome, regulating RNA metabolism and cellular functions. Several types of RNA modifications have been identified to date; they include N6-methyladenosine (m6A), N1-methyladenosine (m1A), 5-methylcytosine (m5C), N7-methylguanosine (m7G), N6,2'-O-dimethyladenosine (m6Am), N4-acetylcytidine (ac4C), etc. RNA modifications, mediated by regulators including writers, erasers, and readers, are associated with carcinogenesis, tumor microenvironment, metabolic reprogramming, immunosuppression, immunotherapy, chemotherapy, etc. A novel perspective indicates that regulatory subunits and post-translational modifications (PTMs) are involved in the regulation of writer, eraser, and reader functions in mediating RNA modifications, tumorigenesis, and anticancer therapy. In this review, we summarize the advances made in the knowledge of different RNA modifications (especially m6A) and focus on RNA modification regulators with functions modulated by a series of factors in cancer, including regulatory subunits (proteins, noncoding RNA or peptides encoded by long noncoding RNA) and PTMs (acetylation, SUMOylation, lactylation, phosphorylation, etc.). We also delineate the relationship between RNA modification regulator functions and carcinogenesis or cancer progression. Additionally, inhibitors that target RNA modification regulators for anticancer therapy and their synergistic effect combined with immunotherapy or chemotherapy are discussed.

The Macrophage migration inhibitory factor is a vital player in Pan-Cancer by functioning as a M0 Macrophage biomarker

BACKGROUND

The role of the macrophage migration inhibitory factor (MIF) has recently attracted considerable attention in cancer research; nonetheless, the insights provided by current investigations remain constrained. Our main objective was to investigate its role and the latent mechanisms within the pan-cancer realm.

METHODS

We used comprehensive pan-cancer bulk sequencing data and online network tools to investigate the association between MIF expression and patient prognosis, genomic instability, cancer cell stemness, DNA damage repair, and immune infiltration. Furthermore, we validated the relationship between MIF expression and M0 macrophages using single-cell datasets, the SpatialDB database, and fluorescence staining. Additionally, we assessed the therapeutic response using the ROC plotter tool.

RESULTS

We observed the upregulation of MIF expression across numerous cancer types. Notably, elevated MIF levels were associated with a decline in genomic stability. We found a significant correlation between increased MIF expression and increased expression of mismatch repair genes, stemness features, and homologous recombination genes across diverse malignancies. Subsequently, through an analysis using ESTIMATE and cytokine results, we revealed the involvement of MIF in immune suppression. Then, we validated MIF as a hallmark of the M0 macrophages involved in tumor immunity. Our study suggests an association with other immune-inhibitory cellular populations and restraint of CD8 + T cells. In addition, we conducted a comparative analysis of MIF expression before and after treatment in three distinct sets of therapy responders and non-responders. Intriguingly, we identified notable disparities in MIF expression patterns in bladder urothelial carcinoma and ovarian cancer following particular therapeutic interventions.

CONCLUSION

Comprehensive pan-cancer analysis revealed notable enrichment of MIF within M0 macrophages, exerting a profound influence on tumor-associated immunosuppression and the intricate machinery of DNA repair.

m5C methylated lncRncr3-MeCP2 interaction restricts miR124a-initiated neurogenesis

Coordination of neuronal differentiation with expansion of the neuroepithelial/neural progenitor cell (NEPC/NPC) pool is essential in early brain development. Our in vitro and in vivo studies identify independent and opposing roles for two neural-specific and differentially expressed non-coding RNAs derived from the same locus: the evolutionarily conserved lncRNA Rncr3 and the embedded microRNA miR124a-1. Rncr3 regulates NEPC/NPC proliferation and controls the biogenesis of miR124a, which determines neuronal differentiation. Rncr3 conserved exons 2/3 are cytosine methylated and bound by methyl-CpG binding protein MeCP2, which restricts expression of miR124a embedded in exon 4 to prevent premature neuronal differentiation, and to orchestrate proper brain growth. MeCP2 directly binds cytosine-methylated Rncr3 through previously unrecognized lysine residues and suppresses miR124a processing by recruiting PTBP1 to block access of DROSHA-DGCR8. Thus, miRNA processing is controlled by lncRNA m5C methylation along with the defined m5C epitranscriptomic RNA reader protein MeCP2 to coordinate brain development.

Detecting genetic gain and loss events in terms of protein domain: Method and implementation

Continuous gain and loss of genes are the primary driving forces of bacterial evolution and environmental adaptation. Studying bacterial evolution in terms of protein domain, which is the fundamental function and evolutionary unit of proteins, can provide a more comprehensive understanding of bacterial differentiation and phenotypic adaptation processes. Therefore, we proposed a phylogenetic tree-based method for detecting genetic gain and loss events in terms of protein domains. Specifically, the method focuses on a single domain to trace its evolution process or on multiple domains to investigate their co-evolution principles. This novel method was validated using 122 Shigella isolates. We found that the loss of a significant number of domains was likely the main driving force behind the evolution of Shigella, which could reduce energy expenditure and preserve only the most essential functions. Additionally, we observed that simultaneously gained and lost domains were often functionally related, which can facilitate and accelerate phenotypic evolutionary adaptation to the environment. All results obtained using our method agree with those of previous studies, which validates our proposed method.

The transcriptional landscape of cancer stem-like cell functionality in breast cancer

BACKGROUND

Cancer stem-like cells (CSCs) have been extensively researched as the primary drivers of therapy resistance and tumor relapse in patients with breast cancer. However, due to lack of specific molecular markers, increased phenotypic plasticity and no clear clinicopathological features, the assessment of CSCs presence and functionality in solid tumors is challenging. While several potential markers, such as CD24/CD44, have been proposed, the extent to which they truly represent the stem cell potential of tumors or merely provide static snapshots is still a subject of controversy. Recent studies have highlighted the crucial role of the tumor microenvironment (TME) in influencing the CSC phenotype in breast cancer. The interplay between the tumor and TME induces significant changes in the cancer cell phenotype, leading to the acquisition of CSC characteristics, therapeutic resistance, and metastatic spread. Simultaneously, CSCs actively shape their microenvironment by evading immune surveillance and attracting stromal cells that support tumor progression.

METHODS

In this study, we associated in vitro mammosphere formation assays with bulk tumor microarray profiling and deconvolution algorithms to map CSC functionality and the microenvironmental landscape in a large cohort of 125 breast tumors.

RESULTS

We found that the TME score was a significant factor associated with CSC functionality. CSC-rich tumors were characterized by an immune-suppressed TME, while tumors devoid of CSC potential exhibited high immune infiltration and activation of pathways involved in the immune response. Gene expression analysis revealed IFNG, CXCR5, CD40LG, TBX21 and IL2RG to be associated with the CSC phenotype and also displayed prognostic value for patients with breast cancer.

CONCLUSION

These results suggest that the characterization of CSCs content and functionality in tumors can be used as an attractive strategy to fine-tune treatments and guide clinical decisions to improve patients therapy response.

Chemical Synthesis of Modified RNA

Ribonucleic acids (RNAs) play a vital role in living organisms. Many of their cellular functions depend critically on chemical modification. Methods to modify RNA in a controlled manner-both in vitro and in vivo-are thus essential to evaluate and understand RNA biology at the molecular and mechanistic levels. The diversity of modifications, combined with the size and uniformity of RNA (made up of only 4 nucleotides) makes its site-specific modification a challenging task that needs to be addressed by complementary approaches. One such approach is solid-phase RNA synthesis. We discuss recent developments in this field, starting with new protection concepts in the ongoing effort to overcome current size limitations. We continue with selected modifications that have posed significant challenges for their incorporation into RNA. These include deazapurine bases required for atomic mutagenesis to elucidate mechanistic aspects of catalytic RNAs, and RNA containing xanthosine, N4-acetylcytidine, 5-hydroxymethylcytidine, 3-methylcytidine, 2'-OCF3, and 2'-N3 ribose modifications. We also discuss the all-chemical synthesis of 5'-capped mRNAs and the enzymatic ligation of chemically synthesized oligoribonucleotides to obtain long RNA with multiple distinct modifications, such as those needed for single-molecule FRET studies. Finally, we highlight promising developments in RNA-catalyzed RNA modification using cofactors that transfer bioorthogonal functionalities.

HATCHet2: clone- and haplotype-specific copy number inference from bulk tumor sequencing data

Bulk DNA sequencing of multiple samples from the same tumor is becoming common, yet most methods to infer copy-number aberrations (CNAs) from this data analyze individual samples independently. We introduce HATCHet2, an algorithm to identify haplotype- and clone-specific CNAs simultaneously from multiple bulk samples. HATCHet2 extends the earlier HATCHet method by improving identification of focal CNAs and introducing a novel statistic, the minor haplotype B-allele frequency (mhBAF), that enables identification of mirrored-subclonal CNAs. We demonstrate HATCHet2's improved accuracy using simulations and a single-cell sequencing dataset. HATCHet2 analysis of 10 prostate cancer patients reveals previously unreported mirrored-subclonal CNAs affecting cancer genes.

NSUN2 promotes lung adenocarcinoma progression through stabilizing PIK3R2 mRNA in an m5C-dependent manner

It is well known that 5-methylcytosine (m5C) is involved in variety of crucial biological processes in cancers. However, its biological roles in lung adenocarcinoma (LAUD) remain to be determined. The LUAD samples were used to assess the clinical value of NOP2/Sun RNA Methyltransferase 2 (NSUN2). Dot blot was used to determine global m5C levels. ChIP and dual-luciferase assays were performed to investigate the MYC-associated zinc finger protein (MAZ)-binding sites in NSUN2 promoter. RNA-seq was used to explore the downstream molecular mechanisms of NSUN2. Dual luciferase reporter assay, m5C-RIP-qPCR, and mRNA stability assay were conducted to explore the effect of NSUN2-depletion on target genes. Cell viability, transwell, and xenograft mouse model were designed to demonstrate the characteristic of NSUN2 in promoting LUAD progression. The m5C methyltransferase NSUN2 was highly expressed and caused elevated m5C methylation in LUAD samples. Mechanistically, MAZ positively regulated the transcription of NSUN2 and was related to poor survival of LUAD patients. Silencing NSUN2 decreased the global m5C levels, suppressed proliferation, migration and invasion, and inhibited activation of PI3K-AKT signaling in A549 and SPAC-1 cells. Phosphoinositide-3-Kinase Regulatory Subunit 2 (PIK3R2) was upregulated by NSUN2-mediated m5C methylation by enhancing its mRNA stabilization and activated the phosphorylation of the PI3K-AKT signaling. The present study explored the underlying mechanism and biological function of NSUN2-meditated m5C RNA methylation in LUAD. NSUN2 was discovered to facilitate the malignancy progression of LUAD through regulating m5C modifications to stabilize PIK3R2 activating the PI3K-AKT signaling, suggesting that NSUN2 could be a novel biomarker and promising therapeutic target for LUAD patients.

Ca2+ homeostasis and male fertility: a target for a new male contraceptive system

Ca2+ is a key secondary messenger that determines sperm motility patterns. Mammalian sperm undergo capacitation, a process to acquire fertilizing ability, in the female reproductive tract. Capacitated sperm change their flagellar waveform to develop hyperactivated motility, which is crucial for successful sperm navigation to the eggs and fertilization. The sperm-specific channel, CATSPER, and an ATPase transporter, PMCA4, serve as major paths for Ca2+ influx and efflux, respectively, in sperm. The ionic paths coordinate Ca2+ homeostasis in the sperm, and their loss-of-function impairs sperm motility, to cause male infertility. In this review, we summarize the physiological significance of these two Ca2+ gates and suggest their potential applications in novel male contraceptives.

The malarial blood transcriptome: translational applications

The blood transcriptome of malaria patients has been used extensively to elucidate the pathophysiological mechanisms and host immune responses to disease, identify candidate diagnostic and prognostic biomarkers, and reveal new therapeutic targets for drug discovery. This review gives a high-level overview of the three main translational applications of these studies (diagnostics, prognostics, and therapeutics) by summarising recent literature and outlining the main limitations and future directions of each application. It highlights the need for consistent and accurate definitions of disease states and subject groups and discusses how prognostic studies must distinguish clearly between analyses that attempt to predict future disease states and those which attempt to discriminate between current disease states (classification). Lastly it examines how many promising therapeutics fail due to the choice of imperfect animal models for pre-clinical testing and lack of appropriate validation studies in humans, and how future transcriptional studies may be utilised to overcome some of these limitations.

RNA m5C methylation modification: a potential therapeutic target for SARS-CoV-2-associated myocarditis

The Corona Virus Disease (COVID-19), caused by the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2), has quickly spread worldwide and resulted in significant morbidity and mortality. Although most infections are mild, some patients can also develop severe and fatal myocarditis. In eukaryotic RNAs, 5-methylcytosine (m5C) is a common kind of post-transcriptional modification, which is involved in regulating various biological processes (such as RNA export, translation, and stability maintenance). With the rapid development of m5C modification detection technology, studies related to viral m5C modification are ever-increasing. These studies have revealed that m5C modification plays an important role in various stages of viral replication, including transcription and translation. According to recent studies, m5C methylation modification can regulate SARS-CoV-2 infection by modulating innate immune signaling pathways. However, the specific role of m5C modification in SARS-CoV-2-induced myocarditis remains unclear. Therefore, this review aims to provide insights into the molecular mechanisms of m5C methylation in SARS-CoV-2 infection. Moreover, the regulatory role of NSUN2 in viral infection and host innate immune response was also highlighted. This review may provide new directions for developing therapeutic strategies for SARS-CoV-2-associated myocarditis.

From Data to Cure: A Comprehensive Exploration of Multi-omics Data Analysis for Targeted Therapies

In the dynamic landscape of targeted therapeutics, drug discovery has pivoted towards understanding underlying disease mechanisms, placing a strong emphasis on molecular perturbations and target identification. This paradigm shift, crucial for drug discovery, is underpinned by big data, a transformative force in the current era. Omics data, characterized by its heterogeneity and enormity, has ushered biological and biomedical research into the big data domain. Acknowledging the significance of integrating diverse omics data strata, known as multi-omics studies, researchers delve into the intricate interrelationships among various omics layers. This review navigates the expansive omics landscape, showcasing tailored assays for each molecular layer through genomes to metabolomes. The sheer volume of data generated necessitates sophisticated informatics techniques, with machine-learning (ML) algorithms emerging as robust tools. These datasets not only refine disease classification but also enhance diagnostics and foster the development of targeted therapeutic strategies. Through the integration of high-throughput data, the review focuses on targeting and modeling multiple disease-regulated networks, validating interactions with multiple targets, and enhancing therapeutic potential using network pharmacology approaches. Ultimately, this exploration aims to illuminate the transformative impact of multi-omics in the big data era, shaping the future of biological research.

NSUN2-Mediated m5C Methylation Impairs Endometrial Receptivity

Impaired endometrial decidualization is the primary cause of recurrent implantation failure (RIF). RNA methylation modification, especially NSUN family mediated m5C, is crucial for various physiological events, such as maternal-to-zygotic transition, gametogenesis, embryonic development, organismal lifespan, and cell cycle. However, the regulatory mechanisms between NSUN family mediated m5C modification and RIF remain unknown. We acquired NSUN2 expression data of 15 human endometrium samples at proliferative and secretory stages from reproductive cell atlas. The overall pattern of m5C sites and genes was elucidated through m5C-BS-seq, whereas the overall m5C levels in different groups were revealed by dot blot assay. BrdU and western blotting assays were carried out to evaluate the role of NSUN2 in proliferation and autophagy. The effects of NSUN2-mediated m5C modification on embryo attachment were evaluated by an in vitro model of a confluent monolayer of Ishikawa cells cocultured with BeWo spheroids, and its downstream targets were evaluated by real-time reverse-transcription PCR and western blotting in Ishikawa cells. The molecular mechanism for NSUN2 regulating its downstream targets' expression was determined by Cut&Tag and coimmunoprecipitation assays. NSUN2 was increased in SOX9+ cells and widespread in epithelial cell type at the proliferative stage by previous single-cell RNA sequencing data. NSUN2 overexpression (NSUN2OE) in the Ishikawa cell line elevated m5C levels and promoted cell proliferation and autophagy. NSUN2OE reduced attachment efficiency of BeWo cell spheres. Overexpressed NSUN2 was found to increase STAT1 and MMP14 mRNA expressions by inducing exon skipping. NSUN2 interacted with CLDN4 through m5C modification, and NSUN2OE or NSUN2 knockdown resulted in a similar variation tendency of CLDN4. Overexpression of NSUN2 increased CLDN4 H3K9ac modification by downregulating SIRT4 expression at the protein level, leading to the upregulation of CLDN4 mRNA expression. Our results uncovered a novel intricate regulatory mechanism between NSUN2-mediated m5C and RIF and suggested a potential new therapeutic strategy for RIF.

eSVD-DE: cohort-wide differential expression in single-cell RNA-seq data using exponential-family embeddings

BACKGROUND

Single-cell RNA-sequencing (scRNA) datasets are becoming increasingly popular in clinical and cohort studies, but there is a lack of methods to investigate differentially expressed (DE) genes among such datasets with numerous individuals. While numerous methods exist to find DE genes for scRNA data from limited individuals, differential-expression testing for large cohorts of case and control individuals using scRNA data poses unique challenges due to substantial effects of human variation, i.e., individual-level confounding covariates that are difficult to account for in the presence of sparsely-observed genes.

RESULTS

We develop the eSVD-DE, a matrix factorization that pools information across genes and removes confounding covariate effects, followed by a novel two-sample test in mean expression between case and control individuals. In general, differential testing after dimension reduction yields an inflation of Type-1 errors. However, we overcome this by testing for differences between the case and control individuals' posterior mean distributions via a hierarchical model. In previously published datasets of various biological systems, eSVD-DE has more accuracy and power compared to other DE methods typically repurposed for analyzing cohort-wide differential expression.

CONCLUSIONS

eSVD-DE proposes a novel and powerful way to test for DE genes among cohorts after performing a dimension reduction. Accurate identification of differential expression on the individual level, instead of the cell level, is important for linking scRNA-seq studies to our understanding of the human population.

RNA m5C modification upregulates E2F1 expression in a manner dependent on YBX1 phase separation and promotes tumor progression in ovarian cancer

5-Methylcytosine (m5C) is a common RNA modification that modulates gene expression at the posttranscriptional level, but the crosstalk between m5C RNA modification and biomolecule condensation, as well as transcription factor-mediated transcriptional regulation, in ovarian cancer, is poorly understood. In this study, we revealed that the RNA methyltransferase NSUN2 facilitates mRNA m5C modification and forms a positive feedback regulatory loop with the transcription factor E2F1 in ovarian cancer. Specifically, NSUN2 promotes m5C modification of E2F1 mRNA and increases its stability, and E2F1 binds to the NSUN2 promoter, subsequently reciprocally activating NSUN2 transcription. The RNA binding protein YBX1 functions as the m5C reader and is involved in NSUN2-mediated E2F1 regulation. m5C modification promotes YBX1 phase separation, which upregulates E2F1 expression. In ovarian cancer, NSUN2 and YBX1 are amplified and upregulated, and higher expression of NSUN2 and YBX1 predicts a worse prognosis for ovarian cancer patients. Moreover, E2F1 transcriptionally regulates the expression of the oncogenes MYBL2 and RAD54L, driving ovarian cancer progression. Thus, our study delineates a NSUN2-E2F1-NSUN2 loop regulated by m5C modification in a manner dependent on YBX1 phase separation, and this previously unidentified pathway could be a promising target for ovarian cancer treatment.

Identification of quantitative trait loci (QTLs) regulating leaf SPAD value and trichome density in mungbean (Vigna radiata L.) using genotyping-by-sequencing (GBS) approach

Quantitative trait loci (QTL) mapping is used for the precise localization of genomic regions regulating various traits in plants. Two major QTLs regulating Soil Plant Analysis Development (SPAD) value (qSPAD-7-1) and trichome density (qTric-7-2) in mungbean were identified using recombinant inbred line (RIL) populations (PMR-1×Pusa Baisakhi) on chromosome 7. Functional analysis of QTL region identified 35 candidate genes for SPAD value (16 No) and trichome (19 No) traits. The candidate genes regulating trichome density on the dorsal leaf surface of the mungbean include VRADI07G24840, VRADI07G17780, and VRADI07G15650, which encodes for ZFP6, TFs bHLH DNA-binding superfamily protein, and MYB102, respectively. Also, candidate genes having vital roles in chlorophyll biosynthesis are VRADIO7G29860, VRADIO7G29450, and VRADIO7G28520, which encodes for s-adenosyl-L-methionine, FTSHI1 protein, and CRS2-associated factor, respectively. The findings unfolded the opportunity for the development of customized genotypes having high SPAD value and high trichome density having a possible role in yield and mungbean yellow vein mosaic India virus (MYMIV) resistance in mungbean.

Accounting for isoform expression increases power to identify genetic regulation of gene expression

A core problem in genetics is molecular quantitative trait locus (QTL) mapping, in which genetic variants associated with changes in the molecular phenotypes are identified. One of the most-studied molecular QTL mapping problems is expression QTL (eQTL) mapping, in which the molecular phenotype is gene expression. It is common in eQTL mapping to compute gene expression by aggregating the expression levels of individual isoforms from the same gene and then performing linear regression between SNPs and this aggregated gene expression level. However, SNPs may regulate isoforms from the same gene in different directions due to alternative splicing, or only regulate the expression level of one isoform, causing this approach to lose power. Here, we examine a broader question: which genes have at least one isoform whose expression level is regulated by genetic variants? In this study, we propose and evaluate several approaches to answering this question, demonstrating that "isoform-aware" methods-those that account for the expression levels of individual isoforms-have substantially greater power to answer this question than standard "gene-level" eQTL mapping methods. We identify settings in which different approaches yield an inflated number of false discoveries or lose power. In particular, we show that calling an eGene if there is a significant association between a SNP and any isoform fails to control False Discovery Rate, even when applying standard False Discovery Rate correction. We show that similar trends are observed in real data from the GEUVADIS and GTEx studies, suggesting the possibility that similar effects are present in these consortia.

Identification and characterization of CLEC11A and its derived immune signature in gastric cancer

Introduction

C-type lectin domain family 11 member A (CLEC11A) was characterized as a growth factor that mainly regulates hematopoietic function and differentiation of bone cells. However, the involvement of CLEC11A in gastric cancer (GC) is not well understood.

Methods

Transcriptomic data and clinical information pertaining to GC were obtained and analyzed from publicly available databases. The relationships between CLEC11A and prognoses, genetic alterations, tumor microenvironment (TME), and therapeutic responses in GC patients were analyzed by bioinformatics methods. A CLEC11A-derived immune signature was developed and validated, and its mutational landscapes, immunological characteristics as well as drug sensitivities were explored. A nomogram was established by combining CLEC11A-derived immune signature and clinical factors. The expression and carcinogenic effects of CLEC11A in GC were verified by qRT-PCR, cell migration, invasion, cell cycle analysis, and in vivo model analysis. Myeloid-derived suppressor cells (MDSCs), regulatory T cells (Tregs), M2 macrophages, and T cells in tumor samples extracted from mice were analyzed utilizing flow cytometry analysis.

Results

CLEC11A was over-expressed in GC, and the elevated CLEC11A expression indicated an unfavorable prognosis in GC patients. CLEC11A was involved in genomic alterations and associated with the TME in GC. Moreover, elevated CLEC11A was found to reduce the benefit of immunotherapy according to immunophenoscore (IPS) and the tumor immune dysfunction, exclusion (TIDE). After validation, the CLEC11A-derived immune signature demonstrated a consistent ability to predict the survival outcomes in GC patients. A nomogram that quantifies survival probability was constructed to improve the accuracy of prognosis prediction in GC patients. Using shRNA to suppress the expression of CLEC11A led to significant inhibitions of cell cycle progression, migration, and invasion, as well as a marked reduction of in vivo tumor growth. Moreover, the flow cytometry assay showed that the knock-down of CLEC11A increased the infiltration of cytotoxic CD8+ T cells and helper CD4+ T into tumors while decreasing the percentage of M2 macrophages, MDSCs, and Tregs.

Conclusion

Collectively, our findings revealed that CLEC11A could be a prognostic and immunological biomarker in GC, and CLEC11A-derived immune signature might serve as a new option for clinicians to predict outcomes and formulate personalized treatment plans for GC patients.

Epigenetic addition of m5C to HBV transcripts promotes viral replication and evasion of innate antiviral responses

Eukaryotic five-methylcytosine (m5C) is an important regulator of viral RNA splicing, stability, and translation. However, its role in HBV replication remains largely unknown. In this study, functional m5C sites are identified in hepatitis B virus (HBV) mRNA. The m5C modification at nt 1291 is not only indispensable for Aly/REF export factor (ALYREF) recognition to promote viral mRNA export and HBx translation but also for the inhibition of RIG-I binding to suppress interferon-β (IFN-β) production. Moreover, NOP2/Sun RNA methyltransferase 2 (NSUN2) catalyzes the addition of m5C to HBV mRNA and is transcriptionally downregulated by the viral protein HBx, which suppresses the binding of EGR1 to the NSUN2 promoter. Additionally, NSUN2 expression correlates with m5C modification of type I IFN mRNA in host cells, thus, positively regulating IFN expression. Hence, the delicate regulation of NSUN2 expression induces m5C modification of HBV mRNA while decreasing the levels of m5C in host IFN mRNA, making it a vital component of the HBV life cycle. These findings provide new molecular insights into the mechanism of HBV-mediated IFN inhibition and may inform the development of new IFN-α based therapies.

Roast: a tool for reference-free optimization of supertranscriptome assemblies

BACKGROUND

Transcriptomic studies involving organisms for which reference genomes are not available typically start by generating de novo transcriptome or supertranscriptome assembly from the raw RNA-seq reads. Assembling a supertranscriptome is, however, a challenging task due to significantly varying abundance of mRNA transcripts, alternative splicing, and sequencing errors. As a result, popular de novo supertranscriptome assembly tools generate assemblies containing contigs that are partially-assembled, fragmented, false chimeras or have local mis-assemblies leading to decreased assembly accuracy. Commonly available tools for assembly improvement rely primarily on running BLAST using closely related species making their accuracy and reliability conditioned on the availability of the data for closely related organisms.

RESULTS

We present ROAST, a tool for optimization of supertranscriptome assemblies that uses paired-end RNA-seq data from Illumina sequencing platform to iteratively identify and fix assembly errors solely using the error signatures generated by RNA-seq alignment tools including soft-clips, unexpected expression coverage, and reads with mates unmapped or mapped on a different contig to identify and fix various supertranscriptome assembly errors without performing BLAST searches against other organisms. Evaluation results using simulated as well as real datasets show that ROAST significantly improves assembly quality by identifying and fixing various assembly errors.

CONCLUSION

ROAST provides a reference-free approach to optimizing supertranscriptome assemblies highlighting its utility in refining de novo supertranscriptome assemblies of non-model organisms.

A chromosome-level genome assembly of the Rhus gall aphid Schlechtendalia chinensis provides insight into the endogenization of Parvovirus-like DNA sequences

The Rhus gall aphid, Schlechtendalia chinensis, feeds on its primary host plant Rhus chinensis to induce galls, which have economic importance in medicines and the food industry. Rhus gall aphids have a unique life cycle and are economically beneficial but there is huge gap in genomic information about this group of aphids. Schlechtendalia chinensis induces rich-tannin galls on its host plant and is emerging as a model organism for both commercial applications and applied research in the context of gall production by insects. Here, we generated a high-quality chromosome-level assembly for the S. chinensis genome, enabling the comparison between S. chinensis and non-galling aphids. The final genome assembly is 344.59 Mb with 91.71% of the assembled sequences anchored into 13 chromosomes. We predicted 15,013 genes, of which 14,582 (97.13%) coding genes were annotated, and 99% of the predicted genes were anchored to the 13 chromosomes. This assembly reveals the endogenization of parvovirus-related DNA sequences (PRDs) in the S. chinensis genome, which could play a role in environmental adaptations. We demonstrated the characterization and classification of cytochrome P450s in the genome assembly, which are functionally crucial for sap-feeding insects and have roles in detoxification and insecticide resistance. This genome assembly also revealed the whole genome duplication events in S. chinensis, which can be considered in comparative evolutionary analysis. Our work represents a reference genome for gall-forming aphids that could be used for comparative genomic studies between galling and non-galling aphids and provides the first insight into the endogenization of PRDs in the genome of galling aphids. It also provides novel genetic information for future research on gall-formation and insect-plant interactions.

ALYREF-mediated RNA 5-Methylcytosine modification Promotes Hepatocellular Carcinoma Progression Via Stabilizing EGFR mRNA and pSTAT3 activation

5-Methylcytosine (m5C) is one of the most ubiquitous modifications of mRNA and contributes to cancer pathogenesis. Aly/REF export factor (ALYREF), an m5C reader, is associated with the prognosis of liver hepatocellular carcinoma (LIHC). However, the effects of ALYREF on the progression of LIHC and the underlying molecular mechanisms remains elusive. Through an analysis of an online database and 3 independent LIHC cohorts, we found that ALYREF was markedly elevated in human liver cancer tissues and was significantly correlated with LIHC clinicopathological parameters, including Ki67+ cell rate, high-grade TNM stage, and poor prognosis. Several experiments were conducted to investigate the molecular basis and functional role of ALYREF-related progression in this study. ALYREF could enhance LIHC cell proliferation, migration, invasion, and epithelial-mesenchymal transition (EMT) in vitro and tumor formation in vivo. Mechanistically, ALYREF promoted the progression of human LIHC through EGFR pathways. Furthermore, ALYREF could directly bind to the m5C modification site of EGFR 3' untranslated region (3' UTR) to stabilize EGFR mRNA. Collectively, ALYREF played a crucial oncogenic role in LIHC via the stabilization of EGFR mRNA and subsequent activation of the STAT3 signaling pathway. Our results may help to elucidate the potential mechanisms of ALYREF-induced m5C modification in the progression of human LIHC.

Single-cell analysis of the epitranscriptome: RNA modifications under the microscope

The identification of mechanisms capable of modifying genetic information by the addition of covalent RNA modifications distinguishes a level of complexity in gene expression which challenges key long-standing concepts of RNA biology. One of the current challenges of molecular biology is to properly understand the molecular functions of these RNA modifications, with more than 170 different ones having been identified so far. However, it has not been possible to map specific RNA modifications at a single-cell resolution until very recently. This review will highlight the technological advances in single-cell methodologies aimed at assessing and testing the biological function of certain RNA modifications, focusing on m6A. These advances have allowed for the development of novel strategies that enable the study of the 'epitranscriptome'. Nevertheless, despite all these improvements, many challenges and difficulties still need fixing for these techniques to work efficiently.

Comprehensive characterization of early-programmed tumor microenvironment by tumor-associated macrophages reveals galectin-1 as an immune modulatory target in breast cancer

Background: In recent years, there has been considerable interest in the therapeutic targeting of tumor-associated macrophages (TAMs) to modulate the tumor microenvironment (TME), resulting in antitumoral phenotypes. However, key mediators suitable for TAM-mediated remodeling of the TME remain poorly understood. Methods: In this study, we used single-cell RNA sequencing analyses to analyze the landscape of the TME modulated by TAMs in terms of a protumor microenvironment during early tumor development. Results: Our data revealed that the depletion of TAMs leads to a decreased epithelial-to-mesenchymal transition (EMT) signature in cancer cells and a distinct transcriptional state characterized by CD8+ T cell activation. Moreover, notable alterations in gene expression were observed upon the depletion of TAMs, identifying Galectin-1 (Gal-1) as a crucial molecular factor responsible for the observed effect. Gal-1 inhibition reversed immune suppression via the reinvigoration of CD8+ T cells, impairing tumor growth and potentiating immune checkpoint inhibitors in breast tumor models. Conclusion: These results provide comprehensive insights into TAM-mediated early tumor microenvironments and reveal immune evasion mechanisms that can be targeted by Gal-1 to induce antitumor immune responses.

Genomic hypomethylation in cell-free DNA predicts responses to checkpoint blockade in lung and breast cancer

Genomic hypomethylation has recently been identified as a determinant of therapeutic responses to immune checkpoint blockade (ICB). However, it remains unclear whether this approach can be applied to cell-free DNA (cfDNA) and whether it can address the issue of low tumor purity encountered in tissue-based methylation profiling. In this study, we developed an assay named iMethyl, designed to estimate the genomic hypomethylation status from cfDNA. This was achieved through deep targeted sequencing of young LINE-1 elements with > 400,000 reads per sample. iMethyl was applied to a total of 653 ICB samples encompassing lung cancer (cfDNA n = 167; tissue n = 137; cfDNA early during treatment n = 40), breast cancer (cfDNA n = 91; tissue n = 50; PBMC n = 50; cfDNA at progression n = 44), and ovarian cancer (tissue n = 74). iMethyl-liquid predicted ICB responses accurately regardless of the tumor purity of tissue samples. iMethyl-liquid was also able to monitor therapeutic responses early during treatment (3 or 6 weeks after initiation of ICB) and detect progressive hypomethylation accompanying tumor progression. iMethyl-tissue had better predictive power than tumor mutation burden and PD-L1 expression. In conclusion, our iMethyl-liquid method allows for reliable noninvasive prediction, early evaluation, and monitoring of clinical responses to ICB therapy.

SRSF2 plays an unexpected role as reader of m5C on mRNA, linking epitranscriptomics to cancer

A common mRNA modification is 5-methylcytosine (m5C), whose role in gene-transcript processing and cancer remains unclear. Here, we identify serine/arginine-rich splicing factor 2 (SRSF2) as a reader of m5C and impaired SRSF2 m5C binding as a potential contributor to leukemogenesis. Structurally, we identify residues involved in m5C recognition and the impact of the prevalent leukemia-associated mutation SRSF2P95H. We show that SRSF2 binding and m5C colocalize within transcripts. Furthermore, knocking down the m5C writer NSUN2 decreases mRNA m5C, reduces SRSF2 binding, and alters RNA splicing. We also show that the SRSF2P95H mutation impairs the ability of the protein to read m5C-marked mRNA, notably reducing its binding to key leukemia-related transcripts in leukemic cells. In leukemia patients, low NSUN2 expression leads to mRNA m5C hypomethylation and, combined with SRSF2P95H, predicts poor outcomes. Altogether, we highlight an unrecognized mechanistic link between epitranscriptomics and a key oncogenesis driver.

The effects of N6-methyladenosine RNA methylation on the nervous system

Epitranscriptomics, also known as "RNA epigenetics", is a type of chemical modification that regulates RNA. RNA methylation is a significant discovery after DNA and histone methylation. The dynamic reversible process of m6A involves methyltransferases (writers), m6A binding proteins (readers), as well as demethylases (erasers). We summarized the current research status of m6A RNA methylation in the neural stem cells' growth, synaptic and axonal function, brain development, learning and memory, neurodegenerative diseases, and glioblastoma. This review aims to provide a theoretical basis for studying the mechanism of m6A methylation and finding its potential therapeutic targets in nervous system diseases.

scReadSim: a single-cell RNA-seq and ATAC-seq read simulator

Benchmarking single-cell RNA-seq (scRNA-seq) and single-cell Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq) computational tools demands simulators to generate realistic sequencing reads. However, none of the few read simulators aim to mimic real data. To fill this gap, we introduce scReadSim, a single-cell RNA-seq and ATAC-seq read simulator that allows user-specified ground truths and generates synthetic sequencing reads (in a FASTQ or BAM file) by mimicking real data. At both read-sequence and read-count levels, scReadSim mimics real scRNA-seq and scATAC-seq data. Moreover, scReadSim provides ground truths, including unique molecular identifier (UMI) counts for scRNA-seq and open chromatin regions for scATAC-seq. In particular, scReadSim allows users to design cell-type-specific ground-truth open chromatin regions for scATAC-seq data generation. In benchmark applications of scReadSim, we show that UMI-tools achieves the top accuracy in scRNA-seq UMI deduplication, and HMMRATAC and MACS3 achieve the top performance in scATAC-seq peak calling.

MLACNN: an attention mechanism-based CNN architecture for predicting genome-wide DNA methylation

Methylation is an important epigenetic regulation of methylation genes that plays a crucial role in regulating biological processes. While traditional methods for detecting methylation in biological experiments are constantly improving, the development of artificial intelligence has led to the emergence of deep learning and machine learning methods as a new trend. However, traditional machine learning-based methods rely heavily on manual feature extraction, and most deep learning methods for studying methylation extract fewer features due to their simple network structures. To address this, we propose a bottomneck network based on an attention mechanism and use new methods to ensure that the deep network can learn more effective features while minimizing overfitting. This approach enables the model to learn more features from nucleotide sequences and make better predictions of methylation. The model uses three coding methods to encode the original DNA sequence and then applies feature fusion based on attention mechanisms to obtain the best fusion method. Our results demonstrate that MLACNN outperforms previous methods and achieves more satisfactory performance.

The effects of Korean Red Ginseng on stress-related neurotransmitters and gene expression: A randomized, double-blind, placebo-controlled trial

Background

Korean Red Ginseng (KRG) is an effective anti-stress treatment. In this study, we investigated the therapeutic potential effects of KRG on relieving stress in a general population using transcriptome analysis.

Methods

We conducted an 8-week clinical pilot study on 90 healthy men who reported stress. The study was completed by 43 participants in the KRG group and 44 participants in the placebo group. Participants were randomized 1:1 to the KRG and placebo groups. We evaluated the stress by stress response inventory (SRI) at baseline and 8 weeks. The main outcomes were changes in the levels of neurotransmitters (NTs) and NT-related gene expression. NTs were analyzed using automated (GC) content, and levels of gene expression were measured by reads per kilobase of transcript per million mapped reads (RPKM).

Results

The KRG group showed significantly preserved epinephrine decrease compared with placebo group at 8 weeks (changes in epinephrine, KRG vs. placebo; -1623.2 ± 46101.5 vs. -35116.3 ± 86288.2, p = 0012). Among subjects who higher SRI score, meaning stress increased compared to baseline, the KRG group showed a smaller decrease in serotonin than the placebo group (changes in serotonin, KRG vs. placebo; -2627.5 ± 5859.1 vs, -8087.4 ± 7162.4, p = 0.005) and a smaller increase in cortisol than the placebo group (changes in cortisol, KRG vs. placebo; 1912.7 ± 10097.75 vs. 8046.2 ± 8050.6 , p = 0.019) in subgroup analysis. Transcriptome findings indicated that KRG intake affects gene expression related with metabolism of choline, adrenalin, and monoamine.

Conclusion

These findings suggest that KRG has beneficial effects on the amelioration of stress response in NTs, and this effect is more prominent in stressful situations. Further clinical studies are required to confirm the anti-stress effect of KRG.

Comprehensive dissection of meiotic DNA double-strand breaks and crossovers in cucumber

Meiotic recombination drives genetic diversity and crop genome optimization. In plant breeding, parents with favorable traits are crossed to create elite varieties. Different hybridizations produce diverse types of segment reshuffling between homologous chromosomes. However, little is known about the factors that cause hybrid-specific changes in crossovers (COs). Here, we constructed 2 F2 populations from crosses between a semiwild and 2 domesticated cucumber (Cucumis sativus) accessions and examined CO events. COs mainly occurred around genes and differed unevenly along chromosomes between the 2 hybrids. Fine-scale CO distributions were suppressed in regions of heterozygous structural variations (SVs) and were accelerated by high sequence polymorphism. C. sativus RADiation sensitive 51A (CsRAD51A) binding, histone H3 lysine 4 trimethylation (H3K4me3) modification, chromatin accessibility, and hypomethylation were positively associated with global CO landscapes and in local DNA double-strand break (DSB) hotspots and genes. The frequency and suppression of COs could be roughly predicted based on multiomic information. Differences in CO events between hybrids could be partially traced to distinct genetic and epigenetic features and were significantly associated with specific DSB hotspots and heterozygous SVs. Our findings identify the genomic and epigenetic features that contribute to CO formation and hybrid-specific divergence in cucumber and provide theoretical support for selecting parental combinations and manipulating recombination events at target genomic regions during plant breeding.

Four functional profiles for fibre and mucin metabolism in the human gut microbiome

BACKGROUND

With the emergence of metagenomic data, multiple links between the gut microbiome and the host health have been shown. Deciphering these complex interactions require evolved analysis methods focusing on the microbial ecosystem functions. Despite the fact that host or diet-derived fibres are the most abundant nutrients available in the gut, the presence of distinct functional traits regarding fibre and mucin hydrolysis, fermentation and hydrogenotrophic processes has never been investigated.

RESULTS

After manually selecting 91 KEGG orthologies and 33 glycoside hydrolases further aggregated in 101 functional descriptors representative of fibre and mucin degradation pathways in the gut microbiome, we used nonnegative matrix factorization to mine metagenomic datasets. Four distinct metabolic profiles were further identified on a training set of 1153 samples, thoroughly validated on a large database of 2571 unseen samples from 5 external metagenomic cohorts and confirmed with metatranscriptomic data. Profiles 1 and 2 are the main contributors to the fibre-degradation-related metagenome: they present contrasted involvement in fibre degradation and sugar metabolism and are differentially linked to dysbiosis, metabolic disease and inflammation. Profile 1 takes over Profile 2 in healthy samples, and unbalance of these profiles characterize dysbiotic samples. Furthermore, high fibre diet favours a healthy balance between profiles 1 and profile 2. Profile 3 takes over profile 2 during Crohn's disease, inducing functional reorientations towards unusual metabolism such as fucose and H2S degradation or propionate, acetone and butanediol production. Profile 4 gathers under-represented functions, like methanogenesis. Two taxonomic makes up of the profiles were investigated, using either the covariation of 203 prevalent genomes or metagenomic species, both providing consistent results in line with their functional characteristics. This taxonomic characterization showed that profiles 1 and 2 were respectively mainly composed of bacteria from the phyla Bacteroidetes and Firmicutes while profile 3 is representative of Proteobacteria and profile 4 of methanogens.

CONCLUSIONS

Integrating anaerobic microbiology knowledge with statistical learning can narrow down the metagenomic analysis to investigate functional profiles. Applying this approach to fibre degradation in the gut ended with 4 distinct functional profiles that can be easily monitored as markers of diet, dysbiosis, inflammation and disease. Video Abstract.

Epigenetic modulations of immune cells: from normal development to tumor progression

The dysfunction of immune cell development often impairs immunological homeostasis, thus causing various human diseases. Accumulating evidence shows that the development of different immune cells from hematopoietic stem cells are highly fine-tuned by different epigenetic mechanisms including DNA methylation, histone modifications, chromatin remodeling and RNA-related regulations. Understanding how epigenetic regulators modulate normal development of immune cells contributes to the identification of new strategies for various diseases. Here, we review recent advances suggesting that epigenetic modulations can orchestrate immune cell development and functions through their impact on critical gene expression. We also discuss the aberrations of epigenetic modulations in immune cells that influence tumor progression, and the fact that underlying mechanisms affect how epigenetic drugs interfere with tumor progression in the clinic.

Selective neuronal vulnerability to deficits in RNA processing

Emerging evidence indicates that errors in RNA processing can causally drive neurodegeneration. Given that RNA produced from expressed genes of all cell types undergoes processing (splicing, polyadenylation, 5' capping, etc.), the particular vulnerability of neurons to deficits in RNA processing calls for careful consideration. The activity-dependent transcriptome remodeling associated with synaptic plasticity in neurons requires rapid, multilevel post-transcriptional RNA processing events that provide additional opportunities for dysregulation and consequent introduction or persistence of errors in RNA transcripts. Here we review the accumulating evidence that neurons have an enhanced propensity for errors in RNA processing alongside grossly insufficient defenses to clear misprocessed RNA compared to other cell types. Additionally, we explore how tau, a microtubule-associated protein implicated in Alzheimer's disease and related tauopathies, contributes to deficits in RNA processing and clearance.

Comparative mucomic analysis of three functionally distinct Cornu aspersum Secretions

Every animal secretes mucus, placing them among the most diverse biological materials. Mucus hydrogels are complex mixtures of water, ions, carbohydrates, and proteins. Uncertainty surrounding their composition and how interactions between components contribute to mucus function complicates efforts to exploit their properties. There is substantial interest in commercializing mucus from the garden snail, Cornu aspersum, for skincare, drug delivery, tissue engineering, and composite materials. C. aspersum secretes three mucus-one shielding the animal from environmental threats, one adhesive mucus from the pedal surface of the foot, and another pedal mucus that is lubricating. It remains a mystery how compositional differences account for their substantially different properties. Here, we characterize mucus proteins, glycosylation, ion content, and mechanical properties that could be used to provide insight into structure-function relationships through an integrative "mucomics" approach. We identify macromolecular components of these hydrogels, including a previously unreported protein class termed Conserved Anterior Mollusk Proteins (CAMPs). Revealing differences between C. aspersum mucus shows how considering structure at all levels can inform the design of mucus-inspired materials.

Biallelic variants in NSUN6 cause an autosomal recessive neurodevelopmental disorder

PURPOSE

5-methylcytosine RNA modifications are driven by NSUN methyltransferases. Although variants in NSUN2 and NSUN3 were associated with neurodevelopmental diseases, the physiological role of NSUN6 modifications on transfer RNAs and messenger RNAs remained elusive.

METHODS

We combined exome sequencing of consanguineous families with functional characterization to identify a new neurodevelopmental disorder gene.

RESULTS

We identified 3 unrelated consanguineous families with deleterious homozygous variants in NSUN6. Two of these variants are predicted to be loss-of-function. One maps to the first exon and is predicted to lead to the absence of NSUN6 via nonsense-mediated decay, whereas we showed that the other maps to the last exon and encodes a protein that does not fold correctly. Likewise, we demonstrated that the missense variant identified in the third family has lost its enzymatic activity and is unable to bind the methyl donor S-adenosyl-L-methionine. The affected individuals present with developmental delay, intellectual disability, motor delay, and behavioral anomalies. Homozygous ablation of the NSUN6 ortholog in Drosophila led to locomotion and learning impairment.

CONCLUSION

Our data provide evidence that biallelic pathogenic variants in NSUN6 cause one form of autosomal recessive intellectual disability, establishing another link between RNA modification and cognition.

A novel signature of the ligand and receptor genes associated with disulfidoptosis for prediction of prognosis, immunologic therapy responses in hepatocellular carcinoma

Backgroud

We aimed to explore the prognostic features of ligand and receptor genes associated with disulfidoptosis in hepatocellular carcinoma (HCC) and establish a risk signature utilizing these genes to predict the prognosis of HCC patients.

Methods

We used scRNA-seq data from GSE166635 to differentiate malignant cells from normal cells using "copykat".The study thoroughly examined the disparities in disulfidoptosis scores and the associated gene expressions between malignant and normal cells.We identified key ligand and receptor genes that are specific to HCC cells.Subsequently, Correlation analysis was conducted to ascertain the ligand and receptor genes associated with disulfidoptosis.We performed univariate Cox regression analysis to identify prognostic ligand and receptor genes associated with disulfidoptosis.We employed LASSO to construct a risk signature using prognostic ligand and receptor genes associated with disulfidoptosis.Lastly, we developed a nomogram model that integrates the risk signature with clinicopathological characteristics.

Results

Malignant cells displayed a marked increase in disulfidoptosis scores and the expression of associated genes compared to normal cells.We identified 110 receptor and ligand genes significantly associated with disulfidoptosis, and narrowed them down to create a risk signature comprising eight genes.Multivariate analysis confirmed the risk signature as an independent prognostic factor for HCC and validated its predictive value for immunotherapy outcomes.A novel nomogram was developed, incorporating stage information and the risk signature derived from disulfidoptosis-related receptor and ligand genes, demonstrating excellent predictive accuracy and reliability in HCC prognosis prediction.

Conclusion

Risk signatures based on disulfidoptosis-associated ligand and receptor genes can effectively predict HCC prognosis and may inform immunotherapy strategies.

A comprehensive review on genomic resources in medicinally and industrially important major spices for future breeding programs: Status, utility and challenges

In the global market, spices possess a high-value but low-volume commodities of commerce. The food industry depends largely on spices for taste, flavor, and therapeutic properties in replacement of cheap synthetic ones. The estimated growth rate for spices demand in the world is ∼3.19%. Since spices grow in limited geographical regions, India is one of the leading producer of spices, contributing 25-30 percent of total world trade. Hitherto, there has been no comprehensive review of the genomic resources of industrially important major medicinal spices to overcome major impediments in varietal improvement and management. This review focuses on currently available genomic resources of 24 commercially significant spices, namely, Ajwain, Allspice, Asafoetida, Black pepper, Cardamom large, Cardamom small, Celery, Chillies, Cinnamon, Clove, Coriander, Cumin, Curry leaf, Dill seed, Fennel, Fenugreek, Garlic, Ginger, Mint, Nutmeg, Saffron, Tamarind, Turmeric and Vanilla. The advent of low-cost sequencing machines has contributed immensely to the voluminous data generation of these spices, cracking the complex genomic architecture, marker discovery, and understanding comparative and functional genomics. This review of spice genomics resources concludes the perspective and way forward to provide footprints by uncovering genome assemblies, sequencing and re-sequencing projects, transcriptome-based studies, non-coding RNA-mediated regulation, organelles-based resources, developed molecular markers, web resources, databases and AI-directed resources in candidate spices for enhanced breeding potential in them. Further, their integration with molecular breeding could be of immense use in formulating a strategy to protect and expand the production of the spices due to increased global demand.

CRISPR/dCas9 DNA methylation editing is heritable during human hematopoiesis and shapes immune progeny

Aging is associated with an abnormal increase in DNA methylation (DNAm) in human gene promoters, including in bone marrow stem cells. DNAm patterns are further perturbed in hematological malignancies such as acute myeloid leukemia but the physiological significance of such epigenetic changes is unknown. Using epigenetic editing of human stem/progenitor cells (HSPCs), we show that p15 methylation affects hematopoiesis in vivo. We edited the CDKN2B (p15) promoter and ARF (p14) using dCas9-3A3L and observed DNAm spreading beyond the gRNA location. We find that despite a transient delivery system, DNAm is maintained during myeloid differentiation in vitro, and hypermethylation of the p15 promoter reduces gene expression. In vivo, edited human HSPCs can engraft the bone marrow of mice and targeted DNAm is maintained in HSPCs long term. Moreover, epigenetic changes are conserved and inherited in both myeloid and lymphoid lineages. Although the proportion of myeloid (CD33+) and lymphoid (CD19+) cells is unaffected, monocyte (CD14+) populations decreased and granulocytes (CD66b+) increased in mice engrafted with p15 hypermethylated HSPCs. Monocytes derived from p15 hypermethylated HSPCs appear to be activated and show increased inflammatory transcriptional programs. We believe these findings have clinical relevance since we found p15 promoter methylation in the peripheral blood of patients with clonal hematopoiesis. Our study shows DNAm can be targeted and maintained in human HSPCs and demonstrated functional relevance of aberrant DNAm on the p15 locus. As such, other aging-associated aberrant DNAm may impact hematopoiesis in vivo.

Genomic insights into biased allele loss and increased gene numbers after genome duplication in autotetraploid Cyclocarya paliurus

BACKGROUND

Autopolyploidy is a valuable model for studying whole-genome duplication (WGD) without hybridization, yet little is known about the genomic structural and functional changes that occur in autopolyploids after WGD. Cyclocarya paliurus (Juglandaceae) is a natural diploid-autotetraploid species. We generated an allele-aware autotetraploid genome, a chimeric chromosome-level diploid genome, and whole-genome resequencing data for 106 autotetraploid individuals at an average depth of 60 × per individual, along with 12 diploid individuals at an average depth of 90 × per individual.

RESULTS

Autotetraploid C. paliurus had 64 chromosomes clustered into 16 homologous groups, and the majority of homologous chromosomes demonstrated similar chromosome length, gene numbers, and expression. The regions of synteny, structural variation and nonalignment to the diploid genome accounted for 81.3%, 8.8% and 9.9% of the autotetraploid genome, respectively. Our analyses identified 20,626 genes (69.18%) with four alleles and 9191 genes (30.82%) with one, two, or three alleles, suggesting post-polyploid allelic loss. Genes with allelic loss were found to occur more often in proximity to or within structural variations and exhibited a marked overlap with transposable elements. Additionally, such genes showed a reduced tendency to interact with other genes. We also found 102 genes with more than four copies in the autotetraploid genome, and their expression levels were significantly higher than their diploid counterparts. These genes were enriched in enzymes involved in stress response and plant defense, potentially contributing to the evolutionary success of autotetraploids. Our population genomic analyses suggested a single origin of autotetraploids and recent divergence (~ 0.57 Mya) from diploids, with minimal interploidy admixture.

CONCLUSIONS

Our results indicate the potential for genomic and functional reorganization, which may contribute to evolutionary success in autotetraploid C. paliurus.

Towards in silico CLIP-seq: predicting protein-RNA interaction via sequence-to-signal learning

We present RBPNet, a novel deep learning method, which predicts CLIP-seq crosslink count distribution from RNA sequence at single-nucleotide resolution. By training on up to a million regions, RBPNet achieves high generalization on eCLIP, iCLIP and miCLIP assays, outperforming state-of-the-art classifiers. RBPNet performs bias correction by modeling the raw signal as a mixture of the protein-specific and background signal. Through model interrogation via Integrated Gradients, RBPNet identifies predictive sub-sequences that correspond to known and novel binding motifs and enables variant-impact scoring via in silico mutagenesis. Together, RBPNet improves imputation of protein-RNA interactions, as well as mechanistic interpretation of predictions.

RNA modifications in cancer

Currently, more than 170 modifications have been identified on RNA. Among these RNA modifications, various methylations account for two-thirds of total cases and exist on almost all RNAs. Roles of RNA modifications in cancer are garnering increasing interest. The research on m6A RNA methylation in cancer is in full swing at present. However, there are still many other popular RNA modifications involved in the regulation of gene expression post-transcriptionally besides m6A RNA methylation. In this review, we focus on several important RNA modifications including m1A, m5C, m7G, 2'-O-Me, Ψ and A-to-I editing in cancer, which will provide a new perspective on tumourigenesis by peeking into the complex regulatory network of epigenetic RNA modifications, transcript processing, and protein translation.

A Chinese indicine pangenome reveals a wealth of novel structural variants introgressed from other Bos species

Chinese indicine cattle harbor a much higher genetic diversity compared with other domestic cattle, but their genome architecture remains uninvestigated. Using PacBio HiFi sequencing data from 10 Chinese indicine cattle across southern China, we assembled 20 high-quality partially phased genomes and integrated them into a multiassembly graph containing 148.5 Mb (5.6%) of novel sequence. We identified 156,009 high-confidence nonredundant structural variants (SVs) and 206 SV hotspots spanning ∼195 Mb of gene-rich sequence. We detected 34,249 archaic introgressed fragments in Chinese indicine cattle covering 1.93 Gb (73.3%) of the genome. We inferred an average of 3.8%, 3.2%, 1.4%, and 0.5% of introgressed sequence originating, respectively, from banteng-like, kouprey-like, gayal-like, and gaur-like Bos species, as well as 0.6% of unknown origin. Introgression from multiple donors might have contributed to the genetic diversity of Chinese indicine cattle. Altogether, this study highlights the contribution of interspecies introgression to the genomic architecture of an important livestock population and shows how exotic genomic elements can contribute to the genetic variation available for selection.

Molecular Characterisation of Blood Microbiome in Patients with Ankylosing Spondylitis and Healthy Controls

Background

In human and animal studies, ankylosing spondylitis (AS) has been increasingly linked to changes in the microbial inhabitants in the human body (microbiome). These studies have primarily now concentrated on the microbial communities that live in the gastrointestinal tract. However, evidence suggests that various molecular techniques can be used to detect microbial DNA in blood circulation. This DNA might be an unknown reservoir of biomarkers with the potential to track alterations in the microbiomes of remote locations, such as the gut. To this end, we compared the presence and identity of microbial DNA in blood samples taken from ankylosing spondylitis patients to healthy control subjects by amplifying and sequencing the bacterial 16S rRNA variable region four.

Methods

The study's design is a case study based on the presence and identity of bacterial DNA in the blood of Ankylosing spondylitis (AS) patients (n = 10) and healthy control subjects (n = 10) was investigated by amplifying and sequencing the bacterial 16S rRNA gene. Blood concentrations of the cytokines TNF alpha, IL-17A, and IL-23 were determined by the Human Magnetic Luminex Screening, and data were analysed using an Unpaired T-test.

Results

Using PCR amplification, 8 of 10 AS patients (80%) and 8 of 10 healthy control samples (80%) had microbial 16S rRNA in their blood. At the phylum level, Proteobacteria (Control = 48.5%, AS = 52%), Firmicutes (Control = 27.8%, AS = 26.1%), Actinobacteria (Control = 15.4%, AS = 10.7%), and Bacteroidetes (Control = 6.5%, AS = 10%) dominated the blood microbiome. A two-tailed Mann-Whitney test found that Ankylosing Spondylitis was associated with significantly elevated Bacteroides (P < 0.05), Prevotella (P < 0.001), and Micrococcus (P < 0.01), and significantly reduced levels of Corynebacterium 1 (P < 0.001), Gemella (P < 0.01), and Alloprevotella (P < 0.05), compared to healthy controls. Additionally, it was shown that the presence of the Prevotella genus was highly positively correlated with higher levels of TNF-alpha (P < 0.05; r = 0.8) in AS patients' blood.

Conclusion

This article reveals that a blood microbiome exists in healthy individuals and identifies particular taxa modulated in disease. These blood-derived signatures indicate that this field needs more research and may be helpful as disease biomarkers.

Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes

Fruiting bodies (sporocarps, sporophores or basidiomata) of mushroom-forming fungi (Agaricomycetes) are among the most complex structures produced by fungi. Unlike vegetative hyphae, fruiting bodies grow determinately and follow a genetically encoded developmental program that orchestrates their growth, tissue differentiation and sexual sporulation. In spite of more than a century of research, our understanding of the molecular details of fruiting body morphogenesis is still limited and a general synthesis on the genetics of this complex process is lacking. In this paper, we aim at a comprehensive identification of conserved genes related to fruiting body morphogenesis and distil novel functional hypotheses for functionally poorly characterised ones. As a result of this analysis, we report 921 conserved developmentally expressed gene families, only a few dozens of which have previously been reported to be involved in fruiting body development. Based on literature data, conserved expression patterns and functional annotations, we provide hypotheses on the potential role of these gene families in fruiting body development, yielding the most complete description of molecular processes in fruiting body morphogenesis to date. We discuss genes related to the initiation of fruiting, differentiation, growth, cell surface and cell wall, defence, transcriptional regulation as well as signal transduction. Based on these data we derive a general model of fruiting body development, which includes an early, proliferative phase that is mostly concerned with laying out the mushroom body plan (via cell division and differentiation), and a second phase of growth via cell expansion as well as meiotic events and sporulation. Altogether, our discussions cover 1 480 genes of Coprinopsis cinerea, and their orthologs in Agaricus bisporus, Cyclocybe aegerita, Armillaria ostoyae, Auriculariopsis ampla, Laccaria bicolor, Lentinula edodes, Lentinus tigrinus, Mycena kentingensis, Phanerochaete chrysosporium, Pleurotus ostreatus, and Schizophyllum commune, providing functional hypotheses for ~10 % of genes in the genomes of these species. Although experimental evidence for the role of these genes will need to be established in the future, our data provide a roadmap for guiding functional analyses of fruiting related genes in the Agaricomycetes. We anticipate that the gene compendium presented here, combined with developments in functional genomics approaches will contribute to uncovering the genetic bases of one of the most spectacular multicellular developmental processes in fungi. Citation: Nagy LG, Vonk PJ, Künzler M, Földi C, Virágh M, Ohm RA, Hennicke F, Bálint B, Csernetics Á, Hegedüs B, Hou Z, Liu XB, Nan S, M. Pareek M, Sahu N, Szathmári B, Varga T, Wu W, Yang X, Merényi Z (2023). Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes. Studies in Mycology 104: 1-85. doi: 10.3114/sim.2022.104.01.

Selective Functionalisation of 5-Methylcytosine by Organic Photoredox Catalysis

The epigenetic modification 5-methylcytosine plays a vital role in development, cell specific gene expression and disease states. The selective chemical modification of the 5-methylcytosine methyl group is challenging. Currently, no such chemistry exists. Direct functionalisation of 5-methylcytosine would improve the detection and study of this epigenetic feature. We report a xanthone-photosensitised process that introduces a 4-pyridine modification at a C(sp3)-H bond in the methyl group of 5-methylcytosine. We propose a reaction mechanism for this type of reaction based on density functional calculations and apply transition state analysis to rationalise differences in observed reaction efficiencies between cyanopyridine derivatives. The reaction is initiated by single electron oxidation of 5-methylcytosine followed by deprotonation to generate the methyl group radical. Cross coupling of the methyl radical with 4-cyanopyridine installs a 4-pyridine label at 5-methylcytosine. We demonstrate use of the pyridination reaction to enrich 5-methylcytosine-containing ribonucleic acid.

Selective Functionalisation of 5-Methylcytosine by Organic Photoredox Catalysis

The epigenetic modification 5-methylcytosine plays a vital role in development, cell specific gene expression and disease states. The selective chemical modification of the 5-methylcytosine methyl group is challenging. Currently, no such chemistry exists. Direct functionalisation of 5-methylcytosine would improve the detection and study of this epigenetic feature. We report a xanthone-photosensitised process that introduces a 4-pyridine modification at a C(sp3 )-H bond in the methyl group of 5-methylcytosine. We propose a reaction mechanism for this type of reaction based on density functional calculations and apply transition state analysis to rationalise differences in observed reaction efficiencies between cyanopyridine derivatives. The reaction is initiated by single electron oxidation of 5-methylcytosine followed by deprotonation to generate the methyl group radical. Cross coupling of the methyl radical with 4-cyanopyridine installs a 4-pyridine label at 5-methylcytosine. We demonstrate use of the pyridination reaction to enrich 5-methylcytosine-containing ribonucleic acid.

mRNA Regulation by RNA Modifications

Chemical modifications on mRNA represent a critical layer of gene expression regulation. Research in this area has continued to accelerate over the last decade, as more modifications are being characterized with increasing depth and breadth. mRNA modifications have been demonstrated to influence nearly every step from the early phases of transcript synthesis in the nucleus through to their decay in the cytoplasm, but in many cases, the molecular mechanisms involved in these processes remain mysterious. Here, we highlight recent work that has elucidated the roles of mRNA modifications throughout the mRNA life cycle, describe gaps in our understanding and remaining open questions, and offer some forward-looking perspective on future directions in the field.