Epigenomics

Abstract. Introduction: Preterm birth drives adverse perinatal maternal and infant health outcomes through heterogeneous symptoms, severity, and etiologies. Delivery prior to reaching 37 weeks of gestation may result from medically indicated intervention for pregnancy complications or spontaneously in the absence of prior symptoms. Placental tissue collected following preterm birth exhibits differential DNA methylation compared to full-term placentas and may indicate pregnancy health during gestation. Placental DNA currently has limited utility for assessing health of ongoing pregnancy, as sampling placental tissue during gestation increases the risk of infection and miscarriage. Risks associated with placental sampling during pregnancy limit the use of DNA methylation in clinical preterm birth prediction. Assessing preterm birth risk during gestation requires non-invasive methods for characterizing placental DNA methylation. Results: We quantified genome-wide DNA methylation patterns of hypermethylated cell-free DNA in pregnant (n = 99) and non-pregnant (n = 93) plasma using cell-free methylated DNA immunoprecipitation sequencing (cfMeDIP-seq). In each sample, we assessed DNA methylation status in 300-bp genomic windows, examining both sequencing read counts and calculated absolute molar DNA amount. Known hypermethylated placental regions, including RASSF1, STAT5A, and ERG promoters showed significantly increased odds of detection in pregnant samples, suggesting enrichment of cell-free placental DNA. Of the 536,444 300-bp windows examined, 173,071 (32%) showed significant enrichment in pregnant plasma. Linear modeling identified 107,505 differentially methylated regions (DMRs) associated with pregnancies later diagnosed with intrauterine growth restriction (IUGR) (n = 22). Alu elements showed increased representation in these DMRs than expected, while other repetitive elements exhibited underrepresentation. Discussion: These results demonstrate cfMeDIP-seq's ability to enrich for cell-free placental DNA and characterize cell-free DNA methylation signatures of pregnancies complicated by IUGR. Enrichment of cell-free placental DNA enables non-invasive profiling of placental DNA methylation from maternal plasma. Detectable epigenetic signatures in maternal plasma may identify pregnancies at elevated risk for preterm birth before clinical symptoms appear. Our findings further highlight the potential of cell-free placental DNA for monitoring pregnancy health.

Background. DNA methylation patterns in cord blood are robustly associated with birthweight in the general population. However, it remains unknown whether these associations extend to clinically relevant populations, such as preterm neonates or those born small for gestational age, and whether they directly reflect birthweight or are driven indirectly by genetic, familial, maternal, and obstetric factors. Methods. We calculated a birthweight methylation profile score (MPSBW) using weights of 835 CpGs previously associated with birthweight in the general population and evaluated its association with birthweight in 67 monochorionic (MC) twin pairs including 134 neonates (97% born preterm) from the Twinlife study. MC twin pairs are identical twins sharing a single placenta, often unequally, which can result in unequal resource distribution and differential fetal growth. Results. We examined the association between within-pair differences in birthweight and MPSBW, thereby estimating the association independent of factors shared equally by co-twins. A 500-gram increase in birthweight was associated with a 0.256 SD increase in MPSBW (p<0.005) in this population of preterm neonates. Adjustment for polygenic score for birthweight (PGSBW) confirmed that the observed epigenetic associations were not driven by common genetic variation underlying birthweight. Interestingly, a similar effect size (0.226 SD per 500 g birthweight increase; p<0.05) was observed in the within-pair analysis, which controls for all shared influences within a twin pair. Conclusion DNA methylation is associated with individual differences in birthweight in a high-risk clinical population of MC twins, independent of shared genetic, familial or maternal influences.

Preeclampsia (PE) is a hypertensive pregnancy-specific disorder and a leading cause of maternal and fetal mortality. A common feature of PE placentas and maternal plasma is dyslipidemia, or abnormal lipid levels, which can increase oxidative stress and endothelial dysfunction. However, the precise transcriptional, post-transcriptional, and epigenetic mechanisms underlying these abnormalities remain poorly characterized. Identifying such changes may clarify disease mechanisms and identify lipid-related PE biomarkers. We conducted a large-scale meta-analysis integrating public placental datasets from NCBI GEO, comprising four DNA methylation (DNAm) datasets (n = 172), three RNA-sequencing datasets (n = 92), and an independent RNA microarray validation cohort (n =146). We evaluated differential DNAm (limma), gene expression (DESeq2), transcript-level shifts (Swish), and alternative splicing (rMATS) in PE versus control placentas, with all analyses stratified by fetal sex via an interaction term model. We also performed placental cell-type deconvolution to quantify PE-associated cell-type proportion changes. Our results demonstrated that lipid-related regulation changes in PE placentas occur primarily at the gene and transcript level, with DNAm showing no changes. We also identified significant isoform switching in PE that were undetected by differential gene expression analysis, and primarily driven by alternative transcription initiation and termination sites rather than alternative splicing. A subset of these isoform switches mapped to pathways dysregulated in PE and were predicted to cause functional protein changes. An interaction term model identified several sex-specific differentially expressed genes (DEGs) in PE, including a subset of male-specific downregulated genes involved in oxidative metabolism. However, many of the remaining sex-specific DEGs across both sexes were previously uncharacterized in the literature. These findings suggest that transcriptional and isoform-level regulation play a role in PE-associated dyslipidemia, with certain regulatory pathways displaying fetal sex-specific patterns. Highlights- Preeclampsia-associated dyslipidemia manifests at the gene and transcript level - Reciprocal isoform switches were missed by standard gene-level analyses - Alternative transcript initiation and termination drove isoform switching - Sex-interaction modeling identified sex-specific transcriptional shifts in PE

IntroductionFetal development represents a critical window during which genetic and environmental influences shape lifelong health. DNA methylation (DNAm) is a candidate underlying mechanism. While individual prenatal exposures have been related to DNAm, no studies have investigated the broader prenatal exposome, nor incorporated genetics with the exposome. Here, we integrated the prenatal exposome and genetics as predictors of DNAm at birth. MethodsWe used data from the Dutch Generation R (n=2282) and English Avon Longitudinal Study of Parents and Children (ALSPAC; n=809) cohorts. We performed epigenome-wide elastic net regression, using Generation R for model development/internal validation and ALSPAC for external validation, to predict DNAm at each CpG site. We used three models: Model 1 included 42 prenatal exposures, Model 2 additionally included child sex, gestational age and birth weight, and Model 3 further included meQTLs. ResultsIn Model 1, the prenatal exposome explained on average 0.7% of DNAm variation across 347 validated CpGs (0.1% of tested CpGs). This increased to 40,044 CpGs (10.2%) with 1.3% of variation explained in Model 2, and 91,305 CpGs (23.2%) with 3.0% of variation explained in Model 3. In Model 1, prenatal smoking was the largest predictor, followed by delivery characteristics, among which meconium-stained amniotic fluid was a novel finding. In Model 3, typically both SNPs and multiple prenatal exposures were selected. DiscussionWe find that genomic associations with cord blood DNAm are stronger and more widespread than prenatal exposures, although typically, the prenatal exposome explains additional variation in DNAm beyond genetic influences.

Prostate cancer remains one of the most common cancers among men globally. While significant strides have been made in diagnosis and treatment, understanding the complex genetic and epigenetic underpinnings of the disease remains crucial for guiding intervention and developing more personalized and effective therapies. The importance of DNA methylation in prostate cancer has been known for some time, but important facets of the modulation of the epigenome during carcinogenesis remain obscure, partly because the bulk of cancer methylation data have been produced using microarray technologies. Here we utilise the TruSeq methyl capture method (EPICseq) to profile the, previously defined, UK Prostate ICGC cohort of well-annotated primary prostate cancers. To this we add methylation sequencing of benign tissue from the same men. These data allow us to identify differentially methylated regions distinguishing cancerous and non-cancerous prostate tissue, while identifying numerous genes whose methylation profiles can perform that task as well as distinguishing between classes of prostate cancer. We describe a describe a methylation-based control mechanism for prostate-cancer-associated SNPs, and show that this seems a likely mechanism of action for a SNP near the MMP7 gene. We describe three novel molecular signatures that arise from different aspects of the biology of prostate cancer revealed by sequencing. Each is shown to be an independent classifier of cancers into groups with different expected times to relapse. These consist of patterns in driver gene methylation, strand-specific methylation, and signal arising in mitochondrial reads. We show that these signatures, combined with existing molecular tools, provide a powerful predictor of time to recurrence. By substantially enhancing understanding of prostate cancer risk, detection, and prognosis, we pave the way for the development of clinical practices that will benefit patients and improve outcomes.

Mammalian preimplantation development requires precisely coordinated lineage decisions to establish the trophectoderm (TE), inner cell mass (ICM), epiblast (EPI), and primitive endoderm (PrE). Glucose metabolism and epigenetic regulation are increasingly recognised as key determinants of lineage specification during preimplantation development. However, how glucose-dependent metabolic cues interface with epigenetic mechanisms to regulate embryonic cell fate remains poorly understood. Here, we investigated the role of glucose in regulating protein methylation by protein arginine methyltransferases (PRMT) in bovine preimplantation development. PRMT1 and its associated histone mark H4R3me2a were detected throughout bovine oocyte maturation and embryo development. Pharmacological inhibition of Type I PRMTs using two structurally distinct inhibitors, GSK3368715 and MS023, markedly reduced global protein asymmetric dimethylarginine (ADMA) and H4R3me2a levels. PRMT inhibition impaired blastocyst cell proliferation, reduced total cell number, and disrupted both first and second lineage decisions, as demonstrated by decreased CDX2- and SOX2-positive TE and ICM cells and reduced NANOG- and GATA6-positive EPI and PrE cell allocation. Mechanistically, Type I PRMT inhibition downregulated key components of the Hippo-associated TE programme, including YAP, TEAD4, and TFAP2C. Consistent effects were observed in mouse embryos, where MS023 treatment reduced ADMA, CDX2, YAP, and TFAP2C expression and impaired TE and ICM allocation. Collectively, our findings identify Type I PRMT-mediated ADMA as an essential epigenetic regulator of early mammalian lineage specification and support a conserved ADMA-Hippo regulatory axis linking arginine methylation to embryonic cell fate decisions. In briefType I protein arginine methyltransferase (PRMT)-mediated asymmetric dimethylarginine (ADMA) is required for proper lineage specification during mammalian preimplantation development. ADMA depletion disrupts Hippo signalling, cell proliferation, and trophectoderm and inner cell mass allocation in bovine and mouse embryos.

Circular RNAs (circRNAs) are an emerging class of RNAs with biomarker potential, but their detection in liquid biopsies is challenging due to low abundance. We developed Ouro-seq, a novel long-read sequencing protocol optimized for full-length circRNA recovery. Applied to urine, cervico-vaginal self-samples from cervical cancer patients, and plasma from lung cancer patients and controls, Ouro-seq recovered 2-5 times more and substantially longer circRNA molecules than conventional methods. Plasma contained predominantly exonic circRNAs, while urine and cervico-vaginal samples were dominated by previously undercharacterized intergenic circRNAs. We also identified extensive alternative circularization and splicing events. Functional analysis revealed distinct specialization patterns: exonic circRNAs showed enhanced miRNA sponging potential, while circRNAs from unplaced genomic scaffolds demonstrated greater peptide-coding capacity. This study establishes Ouro-seq as a valuable tool for comprehensive circRNA characterization in low-yield clinical samples and advances circRNA biology understanding with potential biomarker discovery and disease monitoring applications. MotivationWhile circular RNAs (circRNAs) constitute a minor fraction of total RNA, they may play critical roles in cancer development. CircRNA concentrations are typically too low for detection by Oxford Nanopore Long-Read Sequencing (LRS), particularly in samples with limited RNA content, such as liquid biopsies. Consequently, LRS-based circRNA analysis from liquid biopsies remains unexplored. To overcome these technical limitations, we developed an optimized circRNA enrichment method utilizing short-amplicon suppression, enabling circRNA profiling from urine, plasma, and cervico-vaginal samples.

Many factors including ionizing radiation and iodine deficiency are known to increase thyroid carcinogenesis risk. Our dataset analysis of The Cancer Genome Atlas (TCGA) showed that lower mRNA expression of NK2 homeobox 1 (NKX2-1) transcription factor, a master regulator of genesis, homeostasis, and function of thyroid, is linked to poor prognosis of papillary thyroid cancer patients. Here we provide the findings that thyroid-specific Nkx2-1 conditional knockout (Nkx2-1{Delta}T) mice develop thyroid adenoma and carcinoma in higher frequency with combined exposure to radiation and iodine deficiency than control Nkx2-1fl/fl mice. Iodine deficiency caused oxidative stress, which subsequently resulted in DNA damage, leading to transformation of thyroid follicular cells. RNA-seq gene set enrichment analysis indicated higher production of reactive oxygen species (ROS) in the thyroids of Nkx2-1{Delta}T as compared to Nkx2-1fl/fl mice with combined exposure to radiation and iodine deficiency. This was accompanied by a feedback induction of SOD3 (superoxide dismutase 3) and GPX2 (glutathione peroxidase 2). These antioxidants were naturally expressed at higher levels in the thyroids of Nkx2-1{Delta}T than Nkx2-1fl/fl mice without iodine deficiency or radiation. Nkx2-1{Delta}T thyroids exhibited abnormal follicle architecture and up-regulation of Acox2 (encoding acyl-CoA oxidase 2), which produces hydrogen peroxide. These results suggest that loss of NKX2-1 may contribute to excess ROS production, which elevates basal oxidative stress resulting in the promotion of ROS-induced carcinogenesis. We propose a role for NKX2-1 as a regulator of ROS production homeostasis in the thyroid. Its disturbance would dispose thyroid follicular cells more vulnerable to the ROS-producing carcinogens.

Environmental exposure to endocrine disrupting chemicals, like bisphenol-A (BPA), can impart detrimental effects on developing feto-placental unit, during pregnancy. Placenta remains a central player maintaining this feto-placental homeostasis for sustenance of a healthy pregnancy. Thus, the bisphenol-A mediated endocrine disruption affects the healthy functioning of placenta by altering key processes, such as tissue remodelling, angiogenesis, and metabolism. However, the underlying mechanism of BPA-altered ECM remodelling remains elusive. Therefore, in this study we investigated the BPA mediated changes in placental tissue remodelling using a bisphenol-A exposed murine model during pregnancy. The results reveal that, the phenotypic changes in feto-placental interface correlates with perturbed placental proteome in response to BPA. Further investigation highlights a S100a10-Annexin A2 axis mediated upregulation of tissue plasminogen activator (tPA), which drives altered extracellular matrix (ECM) degradation in placental decidua. This culminates into functional dysregulation in feto-placental axis, leading to reduced size of fetus and placenta. Therefore, this study provides novel insights of a S100a10-Annexin A2 axis associated mechanism for alteration of ECM remodelling in placental decidua due to BPA exposure, which may lead to toxicity related adverse pregnancy outcome.

BackgroundPrenatal growth restriction has been associated with adverse neonatal and long- term health outcomes, yet the epigenetic mechanisms by which an adverse intrauterine environment shapes fetal immune development remain incompletely understood. Monozygotic dichorionic-diamniotic twins with selective fetal growth restriction (sFGR) provide a unique human model for investigating environmentally driven developmental programming independent of genetic variation and inter-twin placental vascular anastomoses. MethodsUmbilical cord blood buffy coat samples were collected from three sFGR and two gestational age-matched concordant control twin pairs. Bulk RNA sequencing and genome-wide DNA methylation analysis were performed, followed by differential expression, pathway enrichment, hematopoietic and immune module analyses, differential methylation, and integrative transcriptomic-epigenomic analyses. ResultsCompared with concordant control twin pairs, discordant twins exhibited broad attenuation of immune and inflammatory transcriptional programs alongside enrichment of erythroid- and hypoxia-related pathways, consistent with adaptive hematopoietic responses to intrauterine stress. Within discordant twin pairs, the growth-restricted co-twins displayed marked transcriptional asymmetry characterized by selective enrichment of cytotoxic lymphoid signatures despite global suppression of myeloid and antigen-presenting cell-associated programs. Integrated transcriptomic and epigenomic analyses further revealed coordinated epigenetic remodeling, with hypomethylated regions in growth-restricted twins enriched for immune regulatory pathways, including T cell differentiation and leukocyte activation. At selected loci, concordant hypomethylation and increased gene expression suggested a potential epigenetic basis for the observed immune remodeling. ConclusionsThese findings suggest that intrauterine growth restriction is associated with coordinated hematopoietic and immune reprogramming at birth, consistent with both compositional and cell-intrinsic alterations. In genetically identical twins, relative growth divergence was associated with polarized transcriptional states, highlighting how intrauterine environmental differences may shape early immune development independent of genetic background.

Background: DNA methylation (DNAm) may capture cumulative genetic, environmental, and lifestyle influences on cardiovascular health. Composite DNAm score based on the American Heart Association Life's Essential 8 (LE8) framework have been linked to clinical events, but their association with early vascular changes and intergenerational effects is unclear. Methods: We studied up to 1432 participants from the multigenerational Young Finns Study (YFS-3G), including parents (G0) and adult offspring (G1). DNAm was measured using Illumina EPIC arrays in 2011 and/or 2018, and carotid intima--media thickness (cIMT) was assessed in 2018. The LE8 DNAm score was calculated as a weighted sum of methylation levels. Associations with cIMT were evaluated in intergenerational, prospective, and cross-sectional settings, adjusting for demographic, technical, and biological covariates and conventional cardiovascular risk factors. Results: Higher parental LE8 DNAm score was associated with lower offspring cIMT ({beta} = -0.022 mm/SD; p-value = 0.02), although the association was attenuated after adjustment for parental cardiovascular risk factors. In G1, a higher baseline DNAm score was associated with lower cIMT measured seven years later ({beta} = -0.030 mm/SD; p-value = 1.1 x 10-5). This association remained significant after adjustment for follow-up cardiovascular risk factors (p-value=0.009) but not after additional adjustment for prior cIMT. Cross-sectionally, higher DNAm score was associated with lower cIMT in both generations, with attenuation after risk factor adjustment in G1 but not G0. Associations with carotid plaque were not significant. Genes associated with the DNAm score were enriched for immune and inflammatory pathways. Conclusions: An LE8-derived DNAm score was associated with lower cIMT across the life course and, to a lesser extent, across generations. These findings suggest that blood DNAm reflects cumulative cardiovascular health and vascular burden and may complement conventional cardiovascular risk assessment.

Aberrant DNA methylation is a hallmark of cancer, but its clinical interpretation remains debated. UHRF1, a key epigenetic adaptor for DNA methylation maintenance and chromatin bivalency regulation in embryonic stem cells, is frequently overexpressed yet shows context-dependent prognostic behaviour. By integrating bulk and single-cell transcriptomics, CpG-resolution methylation, developmental chromatin states, immune profiling and clinical outcomes across gastric (STAD), clear-cell renal (KIRC) and adrenal (ACC) carcinomas, we identified a four-class UHRF1-embryonic morphogenesis (UHRF1-EM) framework resolving this paradox. This axis revealed an inverse prognostic pattern: whilst across all three tumours EM-low and EM-high states mark better or worse prognosis, respectively, UHRF1-high levels associate with favourable outcome in STAD (UH-EML), and unfavourable in KIRC and ACC (UH-EMH). The classification proved reproducible and independently prognostic after adjustment for stage and molecular subtypes, outperforming existing classifiers and exceeding pathological stage in KIRC and ACC. Multivariable models incorporating UHRF1-EM yielded uniformly positive {Delta}C-indices. Hypermethylation associated with the UHRF1-EM axis was enriched at ESC bivalent developmental loci (EM and oncofoetal genes), but not at housekeeping cell-cycle sites. In STAD, this pattern was related to oncofoetal gene downregulation and best prognosis, whereas in KIRC and ACC it matched with gene-body/enhancer methylation, higher EM expression, immunosuppressive microenvironments and worst prognosis. Together, these findings establish the UHRF1-EM axis as a clinically robust molecular classifier and support a mechanistic model in which tumour-specific epigenetic engagement of developmental loci may contribute to the prognostic inversion, providing a foundation for further mechanistic experimental validation.

Melittin exhibits antitumor activity in cervical cancer models, yet the long non-coding RNA (lncRNA) response and associated regulatory networks remain poorly understood. Here, strand-specific RNA-seq data from melittin-treated and untreated U14 murine cervical cancer cells were analyzed to characterize melittin-responsive lncRNAs and explore their potential functional associations. A total of 28,162 lncRNAs were identified, including 27,307 known and 855 novel transcripts. Differential expression analysis revealed 404 differentially expressed lncRNAs (DElncRNAs), comprising 191 upregulated and 213 downregulated lncRNAs, w most of which were predicted to localize to the cytoplasm or nucleus. Cis-target analysis identified 52 neighboring mRNAs as putative targets of 46 DElncRNAs. Functional enrichment highlighted mitochondrial electron transfer and redox-related processes, including the mitochondrial electron transfer flavoprotein complex, electron-transferring-flavoprotein dehydrogenase activity, ubiquinone binding, and quinone binding. In parallel, melittin induced mitochondrial membrane depolarization and increased intracellular reactive oxygen species accumulation in U14 cells. Co-expression analysis further identified 138 lncRNAs co-expressed with 161 mRNAs, which were enriched in chromatin remodeling, DNA replication, and DNA repair. EdU incorporation decreased with increasing melittin concentrations, indicating suppression of DNA synthesis and proliferative activity. RT-qPCR analysis confirmed the expression trends of selected DElncRNAs. Collectively, these findings demonstrate extensive remodeling of the lncRNA landscape in melittin-treated U14 cells and suggest that melittin-responsive lncRNA-mRNA networks are associated with mitochondrial redox disruption and impaired DNA synthesis. This study provides a transcriptomic framework for identifying candidate lncRNA-mRNA regulatory axes underlying the antitumor response to melittin.

Pediatric acute myeloid leukemia (AML) is driven by aberrant transcriptional programs sustained by poorly defined cis-regulatory mechanisms. To systematically identify functional enhancer dependencies, we developed an integrative enhancer discovery strategy that combines H3K27ac CUT&Tag profiling, enhancer-associated transcription, and CRISPR interference (CRISPRi) screening. By leveraging enhancer-associated transcription to prioritize candidate regulatory elements, we identified 321 leukemia-associated enhancers for functional interrogation. This approach uncovered the hematopoietic MYB enhancer (H-ME) within the HBS1L-MYB-AHI1 locus as a critical regulator of leukemic growth. H-ME repression reduced chromatin accessibility and active histone marks at the MYB promoter, suppressed MYB expression, and induced differentiation-associated transcriptional programs. In contrast, selective depletion of the enhancer-associated transcript had no effect on MYB expression or leukemic proliferation, demonstrating that enhancer activity resides within the underlying regulatory DNA element rather than its mature RNA product. H-ME exhibited preferential activity in megakaryocytic leukemia, and its perturbation impaired leukemic growth in primary patient-derived models in vitro and in vivo. Together, our findings establish an integrative framework for the systematic discovery of functional enhancer dependencies and identified H-ME as an RNA-independent regulator of MYB in pediatric AML.

BackgroundOverweight and obese individuals often exhibit lower mortality rates or better prognoses than lean or normal-weight individuals with stroke and other diseases, a phenomenon called the "obesity paradox". Carotid atherosclerosis is the primary cause of ischemic stroke, and body weight serves as a reliable surrogate for adiposity in mice. MethodsPhenotypic and genetic connections of carotid atherosclerosis with body weight were evaluated in 299 F2 mice derived from BALB/cJ and LP/J Apoe knockout (Apoe-/-) mice. F2 mice were fed a Western diet for 12 weeks. Atherosclerotic lesion sizes in left carotid arteries, body weight, coat color, plasma lipids, glucose, small dense LDL ApoB, and malondialdehyde were measured, and 11,000 single nucleotide polymorphism (SNP) markers were genotyped. ResultsCarotid lesion sizes inversely correlated with body weight in both sexes. Genome-wide scans identified two significant quantitative trait loci (QTLs) for carotid atherosclerosis on chromosomes (Chr) 6 and 15 in an additive sex model, and five QTLs on Chr 6, 7, 12, 13, and 15 in an interactive sex model. Adjusting for body weight variation downgraded Chr 15 QTL (Cath5) in both models, whereas other QTLs upgraded in the additive sex model and downgraded in the interactive sex model. Human syntenic region of Cath5 associated with carotid intima-medial thickness (cIMT) and waist-to-hip ratio (WHR). ConclusionsThese findings indicate that the obesity paradox in carotid atherosclerosis is partially driven by shared genetic components that exert opposing effects on adiposity and plaque development and act through sex-dependent mechanisms.

Background: The placenta has a unique transcriptomic profile, including microRNAs that are secreted into maternal circulation throughout pregnancy. MicroRNAs are small, non-coding RNA that post-transcriptionally regulate gene expression. Spontaneous preterm birth (sPTB) is associated with substantial differences in both placental pathophysiology and placental gene expression compared to term birth. We aimed to generate microRNA signatures of sPTB and map them to target genes using a microRNA-mRNA network. Methods: This study was conducted within the Conditions Affecting Neurocognitive Development and Learning in Early childhood (CANDLE) study. Placental samples were collected at delivery, and RNA was isolated for mRNA and microRNA sequencing. To investigate sPTB, this study excluded placental samples of participants with iatrogenic indications for PTB or induced labor. We examined differences in microRNA expression in participants who delivered before 37 weeks (N=35) compared to term participants (N=404) in a series of covariate-adjusted linear regression models. We used paired placental microRNA and mRNA expression data from this cohort to validate associations between computationally predicted microRNA-mRNA pairs and establish a microRNA-mRNA network. Results: Expression of 7 microRNAs were increased in sPTB (FDR<0.05) and were inversely correlated with sPTB-associated genes involved in immune signaling. Expression of 12 microRNAs were decreased in sPTB, including 4 members of the maternally expressed chromosome 14 microRNA cluster (miR-376a-3p, miR-376c-3p, miR-377-3p, and miR-381-3p). These microRNAs were predicted to negatively regulate oxidative phosphorylation genes that were increased in sPTB. The associations between miR-376c-3p and miR-377-3p and oxidative phosphorylation were confirmed in microRNA knockdown experiments. Conclusions: This study highlights potential biological mechanisms by which placental microRNA dysfunction might contribute to sPTB and highlights putative sPTB biomarkers that may be detectable in maternal circulation.

Non-invasive prenatal testing (NIPT) was initially developed to detect chromosomal abnormalities in fetuses through the analysis of cell-free fetal DNA in maternal blood. Recent advancements have expanded NIPT's applications to include the detection of viral infections during pregnancy. However, interpreting pathogen-derived cell-free DNA (cf-DNA) remains clinically complex. This study explores the clinical relevance of hepatitis B virus (HBV) cf-DNA using a dataset of approximately 500,000 NIPT visits and an independent validation cohort of 582 pregnant women (40 HBV-infected), aligned with HBV epidemiology from both population and individual perspectives. Our analysis reveals that HBV cf-DNA is a strong biomarker of high viral infectivity rather than a general marker of infection, suggesting its potential to identify pregnant women at heightened risk of vertical transmission by the end of the first trimester. Additionally, HBV-positive women showed a small but consistent reduction in fetal fraction relative to HBV-negative women across gestational weeks 9 - 17, an association compatible with an early effect of HBV on the placental contribution to cell-free DNA, although the observational design and unmeasured maternal covariates preclude causal inference.

STUDY QUESTION Is genetic liability to endometriosis associated with iron homeostasis, and is this relationship potentially causal? SUMMARY ANSWER Genetic evidence indicates that reduced systemic iron status is associated with increased risk of endometriosis, with evidence of 8 shared genome-wide significant loci and suggestive but inconsistent evidence for causal bidirectional effects. WHAT IS KNOWN ALREADY Endometriosis is a chronic inflammatory condition associated with local iron accumulation within ectopic lesions and peritoneal cavity, resulting from retrograde menstruation and altered iron homeostasis. Epidemiological studies have suggested that women with endometriosis may exhibit reduced systemic iron stores compared to women without endometriosis, reflected by lower circulating ferritin concentrations, although findings have been inconsistent and may be confounded by menstrual blood loss and inflammation. As observational studies cannot distinguish causal relationships from secondary effects or residual confounding, the potential genetic basis linking iron homeostasis and endometriosis risk remains unclear. STUDY DESIGN, SIZE, DURATION We performed genetic analyses using summary statistics from large-scale genome-wide association studies (GWAS) of endometriosis (overall and stage III/IV disease) and five iron biomarkers (serum iron, ferritin, total iron-binding capacity (TIBC), transferrin saturation, and hepcidin). Analyses included genome-wide genetic correlation using linkage disequilibrium score regression (LDSC), identification of shared genetic variants using multi-trait GWAS (MTAG) and bidirectional Mendelian randomisation to evaluate potential causal relationships. PARTICIPANTS/MATERIALS, SETTING, METHODS Iron biomarker summary statistics came from a six-cohort GWAS meta-analysis (HUNT, MGI, SardiNIA, deCODE, Interval, DBDS; N up to 257,953) of blood-derived serum iron, ferritin, transferrin saturation and TIBC (Moksnes et al., 2022). Endometriosis summary statistics came from a 24-study GWAS meta-analysis (60,674 cases, 701,926 controls; European and East Asian ancestry), 12 of which had surgically confirmed cases (Rahmioglu et al., 2023). Genome-wide genetic correlations between iron biomarkers and endometriosis (overall and stage III/IV disease) were estimated using linkage disequilibrium score regression (LDSC), based on summary statistics aligned to the GRCh37 reference genome and restricted to HapMap3 variants. Multi-trait GWAS (MTAG) was applied to each iron biomarker jointly with endometriosis to enhance discovery of genetic loci. Shared loci were functionally annotated using reproductive and iron related tissues from GTEx v8 and blood from eQTLGen expression quantitative trait loci (eQTL) data. Bidirectional Mendelian randomisation (MR) analyses were performed using genome-wide significant variants across multiple clumping thresholds, with inverse-variance weighting (IVW) as the primary method and sensitivity analyses including weighted median, MR-Egger and MR-PRESSO. MAIN RESULTS AND THE ROLE OF CHANCE Genetic correlation analyses suggested that a genetic predisposition to endometriosis is associated with a profile of lower systemic iron availability. Specifically, genetic liability to endometriosis was associated with higher total iron-binding capacity (TIBC; rg=0.16, p=4x10-4), together with lower transferrin saturation (rg=-0.16, p=0.006) and lower ferritin levels (rg=-0.10, p=0.022), findings that are consistent with reduced iron stores. MTAG identified eight additional genome-wide significant loci for endometriosis and eight loci shared with iron biomarkers, including regions implicating coagulation (F5), reproductive biology (WNT4), and immune and vascular pathways (e.g. ABO, STAT6). Mendelian randomisation analyses provided limited and inconsistent evidence for a causal relationship between iron status and endometriosis. Although the inverse-variance weighted (IVW) model showed nominal associations between higher ferritin levels and a lower risk of endometriosis (OR = 0.85, 95% CI 0.76-0.94; p = 0.002), and between genetic liability to endometriosis and higher TIBC (OR = 1.02, 95% CI 1.01-1.04; p = 0.006), these findings were not consistently supported by sensitivity analyses. MR-PRESSO identified a small number of pleiotropic variants, but their removal did not materially alter the results. LIMITATIONS, REASONS FOR CAUTION Iron biomarker GWAS included males and females, potentially obscuring female-specific effects. Dataset availability restricted analyses to European ancestry, limiting applicability to other populations, and to overall and stage III/IV endometriosis, precluding assessment of other disease subtypes. Heterogeneity across SNP instruments, reflected by Cochran's Q statistics, reduced the precision of Mendelian randomisation estimates. Moreover, the genetic instruments explained only between approximately 1.0% and 18.8% of variance in the iron biomarkers, depending on the clumping threshold, which may have limited power to detect causal effects. WIDER IMPLICATIONS OF THE FINDINGS These findings suggest that endometriosis is genetically associated with reduced systemic iron availability and altered iron homeostasis. Thus, lower systemic iron status observed in women with endometriosis may not be explained solely by menstrual blood loss or dietary factors, but reflect an underlying genetic predisposition. Shared genetic loci implicate coagulation, ABO biology, and immune pathways as potential mechanisms linking iron metabolism and endometriosis. Although Mendelian randomisation did not provide consistent evidence for causality, these findings support a shared genetic architecture and warrant further investigation using female-specific GWAS, refined disease subtypes, and multi-omic approaches. Clinically, these findings suggest that low systemic iron status in women with endometriosis may reflect factors beyond established causes of iron deficiency, including an underlying genetic predisposition.

Background: Inflammation and the tumor immune microenvironment contribute to lung adenocarcinoma (LUAD) progression, but the relationship among inflammation-linked transcriptional heterogeneity, patient survival, and immune-state variation remains incompletely defined. Objective: We aimed to identify inflammation-associated LUAD subtypes, derive a parsimonious survival-stratification signature, and characterize its immune and pathway context across public transcriptomic cohorts. Methods: Expression profiles and clinical data were obtained from TCGA-LUAD, GTEx normal lung, and GEO datasets GSE11969, GSE30219, GSE31210, and GSE40791. A curated set of 596 inflammation-related genes was used for consensus clustering. Differential-expression analysis, functional enrichment, univariate Cox regression, and LASSO-Cox modeling were integrated to construct a gene-expression risk score. The prognostic dataset comprised 730 cases and was randomly divided into training (n=502) and internal-validation (n=228) sets; 85 GSE30219 cases formed an external-validation cohort. Immune-cell enrichment, gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), and pan-cancer analyses were used for biological contextualization. Results: The LUAD-versus-control comparison identified 1,305 differentially expressed genes, including 498 upregulated and 807 downregulated genes. Consensus clustering resolved two inflammation-associated subtypes and 67 subtype-associated genes, of which 64 were higher and 3 were lower in Cluster 1 relative to Cluster 2. Thirty-three genes overlapped between the tumor-control and subtype contrasts. LASSO-Cox regression selected CHRDL1, FDCSP, CXCL13, CYP4B1, and S100P. The 1-, 3-, and 5-year areas under the time-dependent receiver operating characteristic curve were 0.6625, 0.6581, and 0.6658 in the training set; 0.7422, 0.6537, and 0.6761 in internal validation; and 0.6560, 0.6387, and 0.6753 in external validation. Risk groups differed across multiple T-cell, B-cell, natural-killer-cell, myeloid, dendritic-cell, macrophage, and granulocyte signatures. Positive GSEA signals included cell cycle (normalized enrichment score [NES]=2.67; adjusted P=1.42 x 10-), DNA replication (NES=2.52; adjusted P=2.52 x 10-), and mismatch repair (NES=2.20; adjusted P=1.77 x 10-). Conclusions: The five-gene expression score separated LUAD survival groups and captured coordinated proliferative and immune transcriptional states. Its moderate discrimination supports further biological and clinical validation rather than immediate clinical application.

Copy-number variants (CNVs) represent an important source of genetic variation that can influence complex traits and disease risk by altering gene dosage, disrupting coding sequence, or modifying regulatory elements. Existing CNV association studies have been limited in scale and have largely focused on European-ancestry populations. We present a CNV genome-wide association study of 13 anthropometric and cardiometabolic traits in 94,730 adults from the China Kadoorie Biobank, a large East Asian study. We identify 19 independent locus-phenotype associations across 15 unique loci. Novel associations include random plasma glucose at 8p23.1 ({beta} = -0.29 SD, P = 5.40x10-) and 14q11.2 ({beta} = +0.43 SD, P = 8.41x10-), diastolic blood pressure at 7p21.1 ({beta} = +0.75 SD, P = 5.25x10-), and duplication-associated reductions in body fat percentage at 12p12.1 ({beta} = -0.74 SD, P = 8.11x10-) and 17q12 ({beta} = -0.56 SD, P = 7.36x10-). We also replicated established dosage-sensitive regions, most prominently at two distinct intervals within 16p11.2 (BP2-BP3 and BP4-BP5), where CNVs show large bidirectional dosage effects across 5 adiposity traits including body mass index ({beta} = -0.84 SD per copy, P = 1.77x10-). These findings identify structural variants contributing to cardiometabolic and anthropometric trait variation in Chinese adults and expand the ancestry diversity of CNV association studies.