NSUN2/YBX1 promotes the progression of breast cancer by enhancing HGH1 mRNA stability through m5C methylation.

BACKGROUND: RNA m5C methylation has been extensively implicated in the occurrence and development of tumors. As the main methyltransferase, NSUN2 plays a crucial regulatory role across diverse tumor types. However, the precise impact of NSUN2-mediated m5C modification on breast cancer (BC) remains unclear. Our study aims to elucidate the molecular mechanism underlying how NSUN2 regulates the target gene HGH1 (also known as FAM203) through m5C modification, thereby promoting BC progression. Additionally, this study targets at preliminarily clarifying the biological roles of NSUN2 and HGH1 in BC. METHODS: Tumor and adjacent tissues from 5 BC patients were collected, and the m5C modification target HGH1 in BC was screened through RNA sequencing (RNA-seq) and single-base resolution m5C methylation sequencing (RNA-BisSeq). Methylation RNA immunoprecipitation-qPCR (MeRIP-qPCR) and RNA-binding protein immunoprecipitation-qPCR (RIP-qPCR) confirmed that the methylation molecules NSUN2 and YBX1 specifically recognized and bound to HGH1 through m5C modification. In addition, proteomics, co-immunoprecipitation (co-IP), and Ribosome sequencing (Ribo-Seq) were used to explore the biological role of HGH1 in BC. RESULTS: As the main m5C methylation molecule, NSUN2 is abnormally overexpressed in BC and increases the overall level of RNA m5C. Knocking down NSUN2 can inhibit BC progression in vitro or in vivo. Combined RNA-seq and RNA-BisSeq analysis identified HGH1 as a potential target of abnormal m5C modifications. We clarified the mechanism by which NSUN2 regulates HGH1 expression through m5C modification, a process that involves interactions with the YBX1 protein, which collectively impacts mRNA stability and protein synthesis. Furthermore, this study is the first to reveal the binding interaction between HGH1 and the translation elongation factor EEF2, providing a comprehensive understanding of its ability to regulate transcript translation efficiency and protein synthesis in BC cells. CONCLUSIONS: This study preliminarily clarifies the regulatory role of the NSUN2-YBX1-m5C-HGH1 axis from post-transcriptional modification to protein translation, revealing the key role of abnormal RNA m5C modification in BC and suggesting that HGH1 may be a new epigenetic biomarker and potential therapeutic target for BC.

Dynamic RNA methylation modifications and their regulatory role in mammalian development and diseases.

Among over 170 different types of chemical modifications on RNA nucleobases identified so far, RNA methylation is the major type of epitranscriptomic modifications existing on almost all types of RNAs, and has been demonstrated to participate in the entire process of RNA metabolism, including transcription, pre-mRNA alternative splicing and maturation, mRNA nucleus export, mRNA degradation and stabilization, mRNA translation. Attributing to the development of high-throughput detection technologies and the identification of both dynamic regulators and recognition proteins, mechanisms of RNA methylation modification in regulating the normal development of the organism as well as various disease occurrence and developmental abnormalities upon RNA methylation dysregulation have become increasingly clear. Here, we particularly focus on three types of RNA methylations: N6-methylcytosine (m6A), 5-methylcytosine (m5C), and N7-methyladenosine (m7G). We summarize the elements related to their dynamic installment and removal, specific binding proteins, and the development of high-throughput detection technologies. Then, for a comprehensive understanding of their biological significance, we also overview the latest knowledge on the underlying mechanisms and key roles of these three mRNA methylation modifications in gametogenesis, embryonic development, immune system development, as well as disease and tumor progression.

DDX21 mediates co-transcriptional RNA m6A modification to promote transcription termination and genome stability.

N6-methyladenosine (m6A) is a crucial RNA modification that regulates diverse biological processes in human cells, but its co-transcriptional deposition and functions remain poorly understood. Here, we identified the RNA helicase DDX21 with a previously unrecognized role in directing m6A modification on nascent RNA for co-transcriptional regulation. DDX21 interacts with METTL3 for co-recruitment to chromatin through its recognition of R-loops, which can be formed co-transcriptionally as nascent transcripts hybridize onto the template DNA strand. Moreover, DDX21's helicase activity is needed for METTL3-mediated m6A deposition onto nascent RNA following recruitment. At transcription termination regions, this nexus of actions promotes XRN2-mediated termination of RNAPII transcription. Disruption of any of these steps, including the loss of DDX21, METTL3, or their enzymatic activities, leads to defective termination that can induce DNA damage. Therefore, we propose that the R-loop-DDX21-METTL3 nexus forges the missing link for co-transcriptional modification of m6A, coordinating transcription termination and genome stability.

DNA methylation clocks for estimating biological age in Chinese cohorts.

Epigenetic clocks are accurate predictors of human chronological age based on the analysis of DNA methylation at specific CpG sites. However, available DNA methylation (DNAm) age predictors are based on datasets with limited ethnic representation. Moreover, a systematic comparison between DNAm data and other omics datasets has not yet been performed. To address these knowledge gaps, we generated and analyzed DNA methylation datasets from two independent Chinese cohorts, revealing age-related DNAm changes. Additionally, a DNA methylation (DNAm) aging clock (iCAS-DNAmAge) and a group of DNAm-based multi-modal clocks for Chinese individuals were developed, with most of them demonstrating strong predictive capabilities for chronological age. The clocks were further employed to predict factors influencing aging rates. The DNAm aging clock, derived from multi-modal aging features (compositeAge-DNAmAge), exhibited a close association with multi-omics changes, lifestyles, and disease status, underscoring its robust potential for precise biological age assessment. Our findings offer novel insights into the regulatory mechanism of age-related DNAm changes and extend the application of the DNAm clock for measuring biological age and aging pace, providing basis for evaluating aging intervention strategies.

HALL: a comprehensive database for human aging and longevity studies.

Diverse individuals age at different rates and display variable susceptibilities to tissue aging, functional decline and aging-related diseases. Centenarians, exemplifying extreme longevity, serve as models for healthy aging. The field of human aging and longevity research is rapidly advancing, garnering significant attention and accumulating substantial data in recent years. Omics technologies, encompassing phenomics, genomics, transcriptomics, proteomics, metabolomics and microbiomics, have provided multidimensional insights and revolutionized cohort-based investigations into human aging and longevity. Accumulated data, covering diverse cells, tissues and cohorts across the lifespan necessitates the establishment of an open and integrated database. Addressing this, we established the Human Aging and Longevity Landscape (HALL), a comprehensive multi-omics repository encompassing a diverse spectrum of human cohorts, spanning from young adults to centenarians. The core objective of HALL is to foster healthy aging by offering an extensive repository of information on biomarkers that gauge the trajectory of human aging. Moreover, the database facilitates the development of diagnostic tools for aging-related conditions and empowers targeted interventions to enhance longevity. HALL is publicly available at https://ngdc.cncb.ac.cn/hall/index.

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.

Emerging Roles of RNA Methylation in Development.

More than 170 different types of chemical modifications have been identified on diverse types of RNA, collectively known as the epitranscriptome. Among them, N6-methyladenine (m6A), 5-methylcytosine (m5C), N1-methyladenine (m1A), and N7-methylguanosine (m7G) as the ubiquitous post-transcriptional modification are widely involved in regulating the metabolic processes such as RNA degradation, translation, stability, and export, mediating important physiological and pathological processes such as stress regulation, immune response, development, and tumorigenesis. Recently, the regulatory role of RNA modification during developmental processes is getting more attention. Therefore, the development of low-input even single-cell and high-resolution sequencing technologies is crucial for the exploration of the regulatory roles of RNA modifications in these important biological events of trace samples.This account focuses on the roles of RNA modifications in various developmental processes. We describe the distribution characteristics of various RNA modifications, catalytic enzymes, binding proteins, and the development of sequencing technologies. RNA modification is dynamically reversible, which can be catalyzed by methyltransferases and eliminated by demethylases. RNA m6A is the most abundant post-transcriptional modification on eukaryote mRNA, which is mainly concentrated near the stop codon, and involves in RNA metabolism regulation. RNA m5C, another most studied RNA modification, has been identified in a various of organisms and RNA species, mainly enriched in the regions downstream of translation initiation sites and broadly distributes across the whole coding sequence (CDS) in mammalian mRNAs. RNA m1A, with a lower abundance than m6A, is widely distributed in various RNA types, mainly locates in the 5' untranslated region (5'UTR) of mRNA and regulates translation. RNA m7G, one of the most common RNA modifications in eukaryotes, has been identified at cap regions and internal positions of RNAs and recently gained considerable attention.Thanks to the development of sequencing technology, m6A has been found to regulate the tumorigenic process, including tumor proliferation, invasion, and metastasis by modulating oncogenes and tumor suppressor genes, and affect oocyte maturation and embryonic development through regulating maternal and zygotic genes. m5C related proteins have been identified to participate in embryonic development, plant growth, and neural stem cell differentiation in a m5C dependent manner. m1A also has been revealed to be involved in these developmental processes. m7G dysregulation mainly involves in neurodevelopmental disorders and neurodegenerative diseases.Collectively, we summarized the gradually exhibited roles of RNA methylation during development, and discussed the possibility of RNA modifications as candidate biomarkers and potential therapeutic targets. The technological development is anticipated as the major driving force to expand our knowledge in this field.

Differential Transcriptomic Landscapes of SARS-CoV-2 Variants in Multiple Organs from Infected Rhesus Macaques.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused the persistent coronavirus disease 2019 (COVID-19) pandemic, which has resulted in millions of deaths worldwide and brought an enormous public health and global economic burden. The recurring global wave of infections has been exacerbated by growing variants of SARS-CoV-2. In this study, the virological characteristics of the original SARS-CoV-2 strain and its variants of concern (VOCs; including Alpha, Beta, and Delta) in vitro, as well as differential transcriptomic landscapes in multiple organs (lung, right ventricle, blood, cerebral cortex, and cerebellum) from the infected rhesus macaques, were elucidated. The original strain of SARS-CoV-2 caused a stronger innate immune response in host cells, and its VOCs markedly increased the levels of subgenomic RNAs, such as N, Orf9b, Orf6, and Orf7ab, which are known as the innate immune antagonists and the inhibitors of antiviral factors. Intriguingly, the original SARS-CoV-2 strain and Alpha variant induced larger alteration of RNA abundance in tissues of rhesus monkeys than Beta and Delta variants did. Moreover, a hyperinflammatory state and active immune response were shown in the right ventricles of rhesus monkeys by the up-regulation of inflammation- and immune-related RNAs. Furthermore, peripheral blood may mediate signaling transmission among tissues to coordinate the molecular changes in the infected individuals. Collectively, these data provide insights into the pathogenesis of COVID-19 at the early stage of infection by the original SARS-CoV-2 strain and its VOCs.

Determining a multimodal aging clock in a cohort of Chinese women.

BACKGROUND: Translating aging rejuvenation strategies into clinical practice has the potential to address the unmet needs of the global aging population. However, to successfully do so requires precise quantification of aging and its reversal in a way that encompasses the complexity and variation of aging. METHODS: Here, in a cohort of 113 healthy women, tiled in age from young to old, we identified a repertoire of known and previously unknown markers associated with age based on multimodal measurements, including transcripts, proteins, metabolites, microbes, and clinical laboratory values, based on which an integrative aging clock and a suite of customized aging clocks were developed. FINDINGS: A unified analysis of aging-associated traits defined four aging modalities with distinct biological functions (chronic inflammation, lipid metabolism, hormone regulation, and tissue fitness), and depicted waves of changes in distinct biological pathways peak around the third and fifth decades of life. We also demonstrated that the developed aging clocks could measure biological age and assess partial aging deceleration by hormone replacement therapy, a prevalent treatment designed to correct hormonal imbalances. CONCLUSIONS: We established aging metrics that capture systemic physiological dysregulation, a valuable framework for monitoring the aging process and informing clinical development of aging rejuvenation strategies. FUNDING: This work was supported by the National Natural Science Foundation of China (32121001), the National Key Research and Development Program of China (2022YFA1103700 and 2020YFA0804000), the National Natural Science Foundation of China (81502304), and the Quzhou Technology Projects (2022K46).

Spatiotemporal transcriptomic atlas reveals the dynamic characteristics and key regulators of planarian regeneration.

Whole-body regeneration of planarians is a natural wonder but how it occurs remains elusive. It requires coordinated responses from each cell in the remaining tissue with spatial awareness to regenerate new cells and missing body parts. While previous studies identified new genes essential to regeneration, a more efficient screening approach that can identify regeneration-associated genes in the spatial context is needed. Here, we present a comprehensive three-dimensional spatiotemporal transcriptomic landscape of planarian regeneration. We describe a pluripotent neoblast subtype, and show that depletion of its marker gene makes planarians more susceptible to sub-lethal radiation. Furthermore, we identified spatial gene expression modules essential for tissue development. Functional analysis of hub genes in spatial modules, such as plk1, shows their important roles in regeneration. Our three-dimensional transcriptomic atlas provides a powerful tool for deciphering regeneration and identifying homeostasis-related genes, and provides a publicly available online spatiotemporal analysis resource for planarian regeneration research.

A spatiotemporal gene expression and cell atlases of the developing rat ovary.

Normal ovarian development is necessary for the production of healthy oocytes. However, the characteristics of oocytes development at different stages and the regulatory relationship between oocytes and somatic cells remain to be fully explained. Here, we combined scRNA-seq and spatial transcriptomic sequencing to profile the transcriptomic atlas of developing ovarian of the rat. We identified four components from developing granulosa cells including cumulus, primitive, mural, and luteal cells, and constructed their differential transcriptional regulatory networks. Several novel growth signals from oocytes to cumulus cells were identified, such as JAG1-NOTCH2 and FGF9-FGFR2. Moreover, we observed three cumulus sequential phases during follicle development determined by the key transcriptional factors in each cumulus phase (Bckaf1, Gata6, Cebpb, etc.), as well as the potential pinpointed roles of macrophages in luteal regression. Altogether, the single-cell spatial transcriptomic profile of the ovary provides not only a new research dimension for temporal and spatial analysis of ovary development, but also valuable data resources and a research basis for in-depth excavation of the mechanisms of mammalian ovary development.

Aberrant m5C hypermethylation mediates intrinsic resistance to gefitinib through NSUN2/YBX1/QSOX1 axis in EGFR-mutant non-small-cell lung cancer.

BACKGROUND: RNA 5-methylcytosine (m5C) modification plays critical roles in the pathogenesis of various tumors. However, the function and molecular mechanism of RNA m5C modification in tumor drug resistance remain unclear. METHODS: The correlation between RNA m5C methylation, m5C writer NOP2/Sun RNA methyltransferase family member 2 (NSUN2) and EGFR-TKIs resistance was determined in non-small-cell lung cancer (NSCLC) cell lines and patient samples. The effects of NSUN2 on EGFR-TKIs resistance were investigated by gain- and loss-of-function assays in vitro and in vivo. RNA-sequencing (RNA-seq), RNA bisulfite sequencing (RNA-BisSeq) and m5C methylated RNA immunoprecipitation-qPCR (MeRIP-qPCR) were performed to identify the target gene of NSUN2 involved in EGFR-TKIs resistance. Furthermore, the regulatory mechanism of NSUN2 modulating the target gene expression was investigated by functional rescue and puromycin incorporation assays. RESULTS: RNA m5C hypermethylation and NSUN2 were significantly correlated with intrinsic resistance to EGFR-TKIs. Overexpression of NSUN2 resulted in gefitinib resistance and tumor recurrence, while genetic inhibition of NSUN2 led to tumor regression and overcame intrinsic resistance to gefitinib in vitro and in vivo. Integrated RNA-seq and m5C-BisSeq analyses identified quiescin sulfhydryl oxidase 1 (QSOX1) as a potential target of aberrant m5C modification. NSUN2 methylated QSOX1 coding sequence region, leading to enhanced QSOX1 translation through m5C reader Y-box binding protein 1 (YBX1). CONCLUSIONS: Our study reveals a critical function of aberrant RNA m5C modification via the NSUN2-YBX1-QSOX1 axis in mediating intrinsic resistance to gefitinib in EGFR-mutant NSCLC.

O-GlcNAcylation promotes the cytosolic localization of the m6A reader YTHDF1 and colorectal cancer tumorigenesis.

O-linked GlcNAc (O-GlcNAc) is an emerging post-translation modification that couples metabolism with cellular signal transduction by crosstalk with phosphorylation and ubiquitination to orchestrate various biological processes. The mechanisms underlying the involvement of O-GlcNAc modifications in N6-methyladenosine (m6A) regulation are not fully characterized. Herein, we show that O-GlcNAc modifies the m6A mRNA reader YTH domain family 1 (YTHDF1) and fine-tunes its nuclear translocation by the exportin protein Crm1. First, we present evidence that YTHDF1 interacts with the sole O-GlcNAc transferase (OGT). Second, we verified Ser196/Ser197/Ser198 as the YTHDF1 O-GlcNAcylation sites, as described in numerous chemoproteomic studies. Then we constructed the O-GlcNAc-deficient YTHDF1-S196A/S197F/S198A (AFA) mutant, which significantly attenuated O-GlcNAc signals. Moreover, we revealed that YTHDF1 is a nucleocytoplasmic protein, whose nuclear export is mediated by Crm1. Furthermore, O-GlcNAcylation increases the cytosolic portion of YTHDF1 by enhancing binding with Crm1, thus upregulating downstream target (e.g. c-Myc) expression. Molecular dynamics simulations suggest that O-GlcNAcylation at S197 promotes the binding between the nuclear export signal motif and Crm1 through increasing hydrogen bonding. Mouse xenograft assays further demonstrate that YTHDF1-AFA mutants decreased the colon cancer mass and size via decreasing c-Myc expression. In sum, we found that YTHDF1 is a nucleocytoplasmic protein, whose cytosolic localization is dependent on O-GlcNAc modification. We propose that the OGT-YTHDF1-c-Myc axis underlies colorectal cancer tumorigenesis.

m6 A promotes planarian regeneration.

Regeneration is the regrowth of damaged tissues or organs, a vital process in response to damages from primitive organisms to higher mammals. Planarian possesses active whole-body regenerative capability owing to its vast reservoir of adult stem cells, neoblasts, providing an ideal model to delineate the underlying mechanisms for regeneration. RNA N6 -methyladenosine (m6 A) modification participates in many biological processes, including stem cell self-renewal and differentiation, in particular the regeneration of haematopoietic stem cells and axons. However, how m6 A controls regeneration at the whole-organism level remains largely unknown. Here, we demonstrate that the depletion of m6 A methyltransferase regulatory subunit wtap abolishes planarian regeneration, potentially through regulating genes related to cell-cell communication and cell cycle. Single-cell RNA-seq (scRNA-seq) analysis unveils that the wtap knockdown induces a unique type of neural progenitor-like cells (NP-like cells), characterized by specific expression of the cell-cell communication ligand grn. Intriguingly, the depletion of m6 A-modified transcripts grn, cdk9 or cdk7 partially rescues the defective regeneration of planarian caused by wtap knockdown. Overall, our study reveals an indispensable role of m6 A modification in regulating whole-organism regeneration.

Omics Views of Mechanisms for Cell Fate Determination in Early Mammalian Development.

During mammalian preimplantation development, a totipotent zygote undergoes several cell cleavages and two rounds of cell fate determination, ultimately forming a mature blastocyst. Along with compaction, the establishment of apicobasal cell polarity breaks the symmetry of an embryo and guides subsequent cell fate choice. Although the lineage segregation of the inner cell mass (ICM) and trophectoderm (TE) is the first symbol of cell differentiation, several molecules have been shown to bias the early cell fate through their inter-cellular variations at much earlier stages, including the 2- and 4-cell stages. The underlying mechanisms of early cell fate determination have long been an important research topic. In this review, we summarize the molecular events that occur during early embryogenesis, as well as the current understanding of their regulatory roles in cell fate decisions. Moreover, as powerful tools for early embryogenesis research, single-cell omics techniques have been applied to both mouse and human preimplantation embryos and have contributed to the discovery of cell fate regulators. Here, we summarize their applications in the research of preimplantation embryos, and provide new insights and perspectives on cell fate regulation.

m6A epitranscriptomic regulation of tissue homeostasis during primate aging.

How N6-methyladenosine (m6A), the most abundant mRNA modification, contributes to primate tissue homeostasis and physiological aging remains elusive. Here, we characterize the m6A epitranscriptome across the liver, heart and skeletal muscle in young and old nonhuman primates. Our data reveal a positive correlation between m6A modifications and gene expression homeostasis across tissues as well as tissue-type-specific aging-associated m6A dynamics. Among these tissues, skeletal muscle is the most susceptible to m6A loss in aging and shows a reduction in the m6A methyltransferase METTL3. We further show that METTL3 deficiency in human pluripotent stem cell-derived myotubes leads to senescence and apoptosis, and identify NPNT as a key element downstream of METTL3 involved in myotube homeostasis, whose expression and m6A levels are both decreased in senescent myotubes. Our study provides a resource for elucidating m6A-mediated mechanisms of tissue aging and reveals a METTL3-m6A-NPNT axis counteracting aging-associated skeletal muscle degeneration.

Cpeb1b-mediated cytoplasmic polyadenylation of shha mRNA modulates zebrafish definitive hematopoiesis.

During vertebrate embryogenesis, hematopoietic stem and progenitor cell (HSPC) production through endothelial-to-hematopoietic transition requires suitable developmental signals, but how these signals are accurately regulated remains incompletely understood. Cytoplasmic polyadenylation, which is one of the posttranscriptional regulations, plays a crucial role in RNA metabolism. Here, we report that Cpeb1b-mediated cytoplasmic polyadenylation is important for HSPC specification by translational control of Hedgehog (Hh) signaling during zebrafish early development. Cpeb1b is highly expressed in notochord and its deficiency results in defective HSPC production. Mechanistically, Cpeb1b regulates hemogenic endothelium specification by the Hedgehog-Vegf-Notch axis. We demonstrate that the cytoplasmic polyadenylation element motif-dependent interaction between Cpeb1b and shha messenger RNA (mRNA) in the liquid-like condensates, which are induced by Pabpc1b phase separation, is required for cytoplasmic polyadenylation of shha mRNA. Intriguingly, the cytoplasmic polyadenylation regulates translation but not stability of shha mRNA, which further enhances the Shha protein level and Hh signal transduction. Taken together, our findings uncover the role of Cpeb1b-mediated cytoplasmic polyadenylation in HSPC development and provide insights into how posttranscriptional regulation can direct developmental signals with high fidelity to translate them into cell fate transition.

Nsun2 coupling with RoRγt shapes the fate of Th17 cells and promotes colitis.

T helper 17 (Th17) cells are a subset of CD4+ T helper cells involved in the inflammatory response in autoimmunity. Th17 cells secrete Th17 specific cytokines, such as IL-17A and IL17-F, which are governed by the master transcription factor RoRγt. However, the epigenetic mechanism regulating Th17 cell function is still not fully understood. Here, we reveal that deletion of RNA 5-methylcytosine (m5C) methyltransferase Nsun2 in mouse CD4+ T cells specifically inhibits Th17 cell differentiation and alleviates Th17 cell-induced colitis pathogenesis. Mechanistically, RoRγt can recruit Nsun2 to chromatin regions of their targets, including Il17a and Il17f, leading to the transcription-coupled m5C formation and consequently enhanced mRNA stability. Our study demonstrates a m5C mediated cell intrinsic function in Th17 cells and suggests Nsun2 as a potential therapeutic target for autoimmune disease.

Serum folate mediates the associations of MTHFR rs1801133 polymorphism with blood glucose levels and gestational diabetes mellitus in Chinese Han pregnant women.

This study aimed to explore the mediation effects of one-carbon metabolism (OCM) related nutrients on the association between MTHFR rs1801133 polymorphism and gestational diabetes mellitus (GDM). Folate, vitamin B12 and homocysteine (Hcy) were measured in the serum of 1254 pregnant women. Linear and logistic regressions were used to estimate the associations of OCM nutrients and MTHFR rs1801133 polymorphism with blood glucose levels and GDM risk. Mediation analysis was applied to test the mediation effects of folate, vitamin B12 and Hcy on the association of MTHFR rs1801133 polymorphism with blood glucose concentrations and GDM. Pregnant women with MTHFR rs1801133 CC genotype had higher serum folate (10·75 v. 8·90 and 9·40 ng/ml) and lower serum Hcy (4·84 v. 4·93 and 5·20 µmol/l) than those with CT and TT genotypes. Folate concentrations were positively associated with fasting plasma glucose (FPG), 1-h plasma glucose (1-h PG), 2-h plasma glucose (2-h PG) and GDM risk. Vitamin B12 levels were negatively correlated with FPG and GDM. Although no direct association was found between MTHFR rs1801133 genotypes and GDM, there were significant indirect effects of MTHFR rs1801133 CC genotype on FPG (ß: 0·005; 95 % CI: 0·001, 0·013), 1-h PG (ß: 0·006; 95 % CI: 0·001, 0·014), 2-h PG (ß: 0·007; 95 % CI: 0·001, 0·015) and GDM (ß: 0·006; 95 % CI: 0·001, 0·014) via folate. In conclusion, serum folate mediates the effect of MTHFR rs1801133 on blood glucose levels and GDM. Our findings potentially provide a feasible GDM prevention strategy via individualised folate supplementation according to the MTHFR genotypes.

scm6A-seq reveals single-cell landscapes of the dynamic m6A during oocyte maturation and early embryonic development.

N6-methyladenosine (m6A) has been demonstrated to regulate RNA metabolism and various biological processes, including gametogenesis and embryogenesis. However, the landscape and function of m6A at single cell resolution have not been extensively studied in mammalian oocytes or during pre-implantation. In this study, we developed a single-cell m6A sequencing (scm6A-seq) method to simultaneously profile the m6A methylome and transcriptome in single oocytes/blastomeres of cleavage-stage embryos. We found that m6A deficiency leads to aberrant RNA clearance and consequent low quality of Mettl3Gdf9 conditional knockout (cKO) oocytes. We further revealed that m6A regulates the translation and stability of modified RNAs in metaphase II (MII) oocytes and during oocyte-to-embryo transition, respectively. Moreover, we observed m6A-dependent asymmetries in the epi-transcriptome between the blastomeres of two-cell embryo. scm6A-seq thus allows in-depth investigation into m6A characteristics and functions, and the findings provide invaluable single-cell resolution resources for delineating the underlying mechanism for gametogenesis and early embryonic development.