SEARCH · PubMed Health

Results for “Copy Number Variations”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources

PubMed2025

Genetic insight into lung neuroendocrine tumors: Notch and Wnt signaling pathways as potential targets.

BACKGROUND: The molecular landscape of lung neuroendocrine neoplasms is still poorly characterized, making it difficult to develop a molecular classification and personalized therapeutic approaches. Significant clinical heterogeneity of these malignancies has been highlighted among poorly differentiated histotypes and within the subgroup of well-differentiated neuroendocrine tumors (NET). Currently, the main prognostic factors of lung NET include stage, histotype, grade, peripheral location, and demographic parameters. To gain deeper insights into the genomic underpinnings of lung NETs, we conducted a pilot investigation to uncover potential genetic mutations and copy number variations (CNVs) implicated in their pathogenesis. METHODS: Formalin-fixed, paraffin-embedded intraoperative tumor biopsies and matched peripheral blood mononuclear cell samples were collected from six consecutive patients with lung NETs. The whole exome sequencing (WES) was performed to profile germline and somatic mutations, identify novel genetic alterations, and detect CNVs. Clinical and pathological data were systematically documented at diagnosis and during follow-up. RESULTS: The WES analysis identified a subset of mutations shared between germline and somatic; some were of particular clinical interest as they were associated with tumor proliferation and potential therapeutic targets such as the genes KDM5C, ATR, COL7A1, NOTCH4, PTPRS, SMO, SPEN, SPTA1, TAF1. These mutations were predominantly linked to chromatin remodeling and were involved in critical oncogenic pathways such as Notch and Wnt signaling. CONCLUSIONS: This pilot study highlights the potential role of NGS analysis on solid biopsy in the assessment of the mutational profile of lung NET. A comparison of germline and somatic mutations is critical to identifying putative tumor driver mutations. In perspective, the enrichment of a subpopulation of cancer cells in the blood, with one or more specific mutations, is information of enormous clinical relevance, either for prognosis or therapeutic decisions. Translational studies on large prospective series are required to establish the role of liquid biopsy in lung NET.

Humans↗

PubMed2025

Mapping Cerebellar Morphology in 15q11.2 CNV Carriers Using Normative Modeling.

Copy number variations (CNVs) at the 15q11.2 locus of the human genome have been associated with altered brain structure and increased risk for neurodevelopmental and neuropsychiatric disorders. The cerebellum is increasingly seen as a crucial brain region for neurodevelopmental conditions, yet the effects of 15q11.2 CNVs on cerebellar morphology remain largely unclear. Importantly, 15q11.2 CNVs shows reduced or incomplete penetrance (meaning that not all CNV carriers are affected) and variable expressivity (meaning that symptoms may differ between individuals with the same genetic alteration). Thus, there is a need to not only assess group differences, but also to quantify anatomical variability at the individual level. Here, we address these issues using normative models of brain anatomy trained on large datasets (n > 52k, age range: 3-85) to assess both group and individual-level deviations in cerebellar anatomy in carriers of 15q11.2 deletions (n = 120, mean [SD] age= 64.95 [7.58]) and duplications (n = 149, mean [SD] age=64.31 [7.21]), compared to non-carriers (n = 19,028, mean [SD] age=64.31 [7.58]). Group-level case-control analyses revealed significantly smaller total and regional cerebellar volumes in both deletion and duplication carriers, though with small effect sizes. Individual-level deviation analyses, capturing pronounced alterations in specific individuals, revealed a heterogeneous pattern among carriers. Overall, our findings suggest that CNVs at the 15q11.2 locus exert modest and highly individualized effects on cerebellar morphology.

15q11.2↗

PubMed2025

Novel genetic determinants contribute to hearing loss in a central European cohort with enlarged vestibular aqueduct.

BACKGROUND: The enlarged vestibular aqueduct (EVA) is the most commonly detected inner ear malformation. Biallelic pathogenic variants in the SLC26A4 gene, coding for the anion exchanger pendrin, are frequently involved in determining Pendred syndrome and nonsyndromic autosomal recessive hearing loss DFNB4 in EVA patients. In Caucasian cohorts, the genetic determinants of EVA remain unknown in approximately 50% of cases. We have recruited a cohort of 32 Austrian patients with hearing loss and EVA to define the prevalence and type of pathogenic sequence alterations in SLC26A4 and discover novel EVA-associated genes. METHODS: Sanger sequencing, single nucleotide polymorphism (SNP) assays, copy number variation (CNV) testing, and Exome Sequencing (ES) were employed for gene analysis. Cell-based functional and molecular assays were used to discriminate between gene variants with and without impact on protein function. RESULTS: SLC26A4 biallelic variants were detected in 5/32 patients (16%) and monoallelic variants in 5/32 patients (16%). The pathogenicity of the uncharacterized SLC26A4 protein variants was assigned or excluded based on their ion transport function and cellular abundance. The monoallelic or biallelic Caucasian EVA haplotype was detected in 7/32 (22%) patients, but its pathogenicity could not be confirmed. X-linked pathogenic variants in POU3F4 (2/32, 6%) and biallelic pathogenic variants in GJB2 (2/32, 6%) were also found. No CNV of SLC26A4 and STRC genes was detected. ES of eleven undiagnosed patients with bilateral EVA detected rare sequence variants in six EVA-unrelated genes (monoallelic variants in SCD5, REST, EDNRB, TJP2, TMC1, and two variants in CDH23) in five patients (5/11, 45%). Cell-based assays showed that the TJP2 variant leads to a mislocalized protein product forming dimers with the wild-type, supporting autosomal dominant pathogenicity. The genetic causes of hearing loss and EVA remained unidentified in (14/32) 44% of patients. CONCLUSIONS: The present investigation confirms the role of SLC26A4 in determining hearing loss with EVA, identifies novel genes in this pathophysiological context, highlights the importance of functional testing to exclude or assign pathogenicity of a given gene variant, proposes a possible diagnostic workflow, suggests a novel pathomechanism of disease for TJP2, and highlights voids of knowledge that deserve further investigation.

Humans↗

PubMed2026

A Novel Splice Variant in the COL1A1 Gene Leads to Exon 46 Skipping and Osteogenesis Imperfecta.

BACKGROUND: Osteogenesis imperfecta (OI) is a clinical and genetic disorder characterised by bone fragility, growth deficiency and skeletal deformity. Ninety per cent of OI cases are attributable to autosomal dominant variants in the COL1A1 and COL1A2 genes. METHODS: Candidate variants were identified and verified through trio whole-exome sequencing (trio-WES), copy number variation sequencing (CNV-seq) and Sanger sequencing. Minigene splicing assays were performed in HeLa and HEK293T cells with pcDNA3.1 and pcMINI-C vectors to investigate the function of the candidate variants. A systematic review of COL1A1 splicing variants and the corresponding genotype-phenotype spectrum was performed. RESULTS: Trio-WES revealed a novel heterozygous variant in the C-terminal region of the COL1A1 gene: NM_000088.4:c.3423+5G>A. Sanger sequencing confirmed the variant in both the proband (II-2) and her foetus (III-1) who were clinically suspected of having OI. The c.3423+5G>A variant causes complete skipping of Exon 46, as demonstrated by a minigene splicing assay. We retrieved 419 COL1A1 splicing variants from PubMed, excluded 15 without phenotypic data and 2 linked to Ehlers-Danlos syndrome and stratified the remaining 402 variants into three types on the basis of splice site location: (1) Variants at canonical splicing sites (77.8%, 313/402) mostly cause mild phenotypes, whereas a minority may be severe. (2) Intron variants in other locations, such as splice region variants (17.9%, 72/402), usually cause mild clinical phenotypes, and deep intronic splice variants (0.4%, 2/402) that may result in severe phenotypes. (3) Other variants (3.7%, 15/402), such as exon variants or fragment loss, are extremely rare. We also preliminarily discuss the mechanisms underlying phenotypic variability and the characteristics of C-terminal variants. CONCLUSIONS: This intron variant in COL1A1 was classified as likely pathogenic and was confirmed to disrupt COL1A1 expression. The summary analysis results also revealed a correlation among splicing variants, C-terminal region variants and disease, suggesting that variant location provides a useful framework for prognosis prediction.

Female↗

PubMed2026

Methylome Profiling of Cartilage Tumors: A Promising New Diagnostic Tool?

DNA methylation and copy number variation (CNV) profiling has emerged as a promising tool for the classification of bone and soft tissue tumors. We evaluated its utility in cartilage tumors, where distinguishing low-grade from high-grade conventional central chondrosarcomas (CSs) and atypical cartilaginous tumors (ACTs) from enchondromas (ECs) is a frequent diagnostic challenge, particularly on biopsy material. We analyzed 214 chondrogenic tumors, including ECs, ACTs, conventional CSs, dedifferentiated chondrosarcomas (DDCSs), and clear cell CSs, and determined their IDH1/2 mutation status. Unsupervised dimensionality reduction of genome-wide DNA methylation patterns revealed 4 clusters among IDH-mutant (MUT) tumors (IDH-MUT-1: mostly ECs and ACTs and some high-grade CSs; IDH-MUT-2: predominantly high-grade CSs; IDH-MUT-3: largely DDCSs; and IDH-MUT-SB: distinct skull base group with a markedly different methylation pattern) and 2 clusters among IDH-wild-type (WT) tumors (IDH-WT-1 and IDH-WT-2: both primarily high-grade CSs, with IDH-WT-2 showing higher tumor grade and more extensive CNVs). Clear cell CSs formed a separate cluster. The amount of CNVs, including loss of CDKN2A, increased with tumor grade, reflecting increased genomic instability during chondrosarcoma progression. Supervised classifiers trained separately, both on methylation and CNV data, and distinguished low-grade and high-grade cartilaginous tumors with area under the curve values of 0.87 to 0.97 and 85% to 90% accuracy. Furthermore, we tested whether DDCSs can be distinguished from metastatic carcinomas and other high-grade sarcomas of the bone. Across 246 reference samples, a supervised classifier achieved 97.2% accuracy (area under the curve, 99.8%) and correctly identified 30 of 32 DDCSs (93.8%). These results indicate that DNA methylation and CNV data analysis provide a valuable tool for distinguishing most low- and high-grade CSs, with additional utility also in differentiating DDCS from morphologic mimics.

cartilaginous tumors↗

PubMed2026

Age-associated epigenomic heterogeneity in papillary tumors of the pineal region: a multicenter YoungNOA investigation.

BACKGROUND: Papillary tumors of the pineal region (PTPR) are rare CNS neoplasms with adult and pediatric presentations, but whether age defines distinct molecular biology is unclear. METHODS: We assembled a multicenter retrospective cohort of 86 histologically confirmed PTPR with genome-wide DNA methylation data, comprising 62 adult and 24 pediatric tumors. Molecular subgroup, array platform, sex, and tumor purity were incorporated into multivariable models. Analyses included DNA methylation class assignment, differential methylation, copy-number variation (CNV), epigenetic mitotic-clock scores, methylation-based tumor microenvironment deconvolution, and descriptive survival evaluation. RESULTS: Adult and pediatric tumors mapped within the established PTPR-A and PTPR-B methylation framework rather than forming age-defined methylation classes. Pediatric tumors were enriched for PTPR-B (22 of 24 tumors [91.7%]) compared with adult tumors (39 of 62 [62.9%]). After adjustment for methylation-based subgroup as well as technical and biological covariates, 2,923 CpG probes were associated with age at a false discovery rate (FDR) threshold below 0.05, and 530 also met the prespecified effect-size threshold. Global methylation summaries were similar between age groups. CNV patterns were dominated by molecular subgroup; adjusted genomic CNV load was not independently associated with pediatric age. In contrast, epiTOC2 intrinsic rate score and the methylation signature represented by the first principal component (PC1) showed age-associated effects independent of molecular subgroup. Methylation-based deconvolution suggested a limited microenvironmental signal, with neutrophil fraction showing the most consistent adjusted association. CONCLUSIONS: Adult and pediatric PTPR share the established PTPR-A/PTPR-B framework. Pediatric tumors, particularly within PTPR-B, showed age-associated DNA methylation differences and higher epigenetic mitotic-clock (epiTOC2) scores in this retrospective cohort. These tissue-level associations do not establish clinical risk or treatment implications and require prospective clinical annotation and orthogonal validation.

Humans↗

PubMed2026

Exploring precision risk in pediatric vesicoureteral reflux: Innate immune gene variations and reflux outcomes in the RIVUR cohort.

INTRODUCTION: Children with vesicoureteral reflux (VUR) are at increased risk for morbidity from recurrent urinary tract infections (UTIs), yet the factors influencing spontaneous VUR resolution remain poorly defined. This study evaluates whether genetic variations in key urinary innate immune effectors (DEFA1A3, DMBT1, and RNASE7) influences VUR resolution and interacts with prophylaxis to alter clinical response. METHODS: We conducted a secondary analysis of 303 RIVUR participants with available DEFA1A3 and DMBT1 copy number variation (CNV) data and RNASE7 rs1263872 genotype. Primary outcomes were (1) VUR improvement (decrease in grade) and (2) VUR resolution at study exit. Multivariable logistic regression models included genotype, treatment, and their interactions, adjusting for age, sex, baseline grade (high vs low), laterality, bowel/bladder dysfunction, and any UTI. Internal validation used 2000-sample bootstrap with bias-corrected and accelerated confidence intervals and influence diagnostics. RESULTS: Clinical covariates did not significantly predict VUR improvement. Children with DEFA1A3 CNV >5 had higher odds of improvement (OR 2.36, 95% CI 1.12-4.96, p = 0.023), an effect that remained significant in bootstrap analyses. High-grade VUR was associated with lower odds of resolution (OR 0.34, 95% CI 0.12-0.94, p = 0.038). A significant interaction was observed between prophylaxis and high DMBT1 copy number for VUR resolution (interaction OR 2.99, 95% CI 1.11-8.04, p = 0.031); no interaction was seen for improvement. RNASE7 rs1263872 was not associated with either outcome. CONCLUSION: Innate immune gene variation may contribute to heterogeneity in VUR outcomes. High DEFA1A3 copy number was associated with reflux improvement and a DMBT1-prophylaxis interaction was associated with reflux resolution. The results of this study is hypothesis-generating and prompt further evaluation to assess whether a subset of children may experience structural benefit from prophylaxis or have a more favorable natural history based on their innate immune genotype.

Humans↗

PubMed2026

A Novel Complete F8 Tandem Duplication Causing Elevated Factor VIII Activity and Associated with Venous Thromboembolism.

Background Coagulation factor VIII (FVIII) is a critical component of the intrinsic coagulation pathway. While elevated FVIII levels are an established risk factor for venous thromboembolism (VTE), genetic variants in the F8 gene directly causing such elevations remain scarce. Here, we report a novel complete F8 tandem duplication identified in a female patient with splanchnic venous thrombosis (SVT). Methods We performed genetic testing using a thrombophilia panel targeting 35 genes involved in thrombosis and haemostasis to detect both point variants and copy number variations (CNVs). Family co-segregation analysis and phenotypic assays for FVIII and von Willebrand factor (VWF) were conducted. The structural basis of the identified F8 copy number gain was elucidated using optical genome mapping (OGM). Full-length F8 mRNA amplification, quantitative PCR, plasma FVIII Western blotting, and X-chromosome inactivation analysis were performed to assess the functional consequences of the duplication. Thrombin generation test (TGT) was employed to assess the hypercoagulable state. Results Genetic testing identified three copies of all 26 exons of the F8 gene in the proband, which was also detected in her mother (CNVs = 3) and son (CNVs = 2). One-stage clotting and chromogenic assays confirmed persistently elevated FVIII activity in the proband and her mother, accompanied by increased FVIII antigen levels. The OGM analysis confirmed a 229 kb tandem duplication including the F8 gene on one of the proband's X chromosomes. The junction regions exhibited high sequence homology and were rich in repetitive sequences, which precluded precise breakpoint mapping. Full-length F8 mRNA amplification revealed no aberrant transcripts, whereas quantitative PCR showed increased F8 mRNA expression in all carriers. Plasma FVIII Western blotting indicated FVIII heavy and light chains of expected molecular weights with increased band intensity in carriers. X-chromosome inactivation analysis in female carriers showed no significant skewing. TGT in two available carriers showed increased thrombin generation compared with a normal control at both low (1 pM) and high (5 pM) tissue factor concentrations. Conclusion We identified a novel complete F8 tandem duplication associated with increased FVIII expression and a hypercoagulable phenotype in a female patient with SVT. These findings support F8 gene dosage gain as a rare gain-of-function mechanism contributing to elevated FVIII levels and thrombophilia, while variation in VWF levels and acquired risk factors may modify thrombotic penetrance.

coagulation factor VIII↗

PubMed2026

Metab8D: a metabolic regulome network from multiomics and machine learning.

To explore multiomic regulation of the metabolome, we used machine learning to predict metabolomic variation across ~1000 different cancer cell lines with matched omics data from eight biomolecular classes: genomic copy number variation, mutations, DNA methylation, histone post-translational modifications (PTMs), transcriptomics and RNA splice variants, non-coding transcriptomics (miRNA and lncRNA), proteomics, and phosphoproteomics. Overall, the metabolome is tightly associated with the transcriptome, with coding and non-coding RNAs emerging as top predictors. Peripheral metabolites are predictable via levels of corresponding enzymes, while those in central metabolism require combinatorial predictors in signaling and redox pathways, and may not reflect corresponding pathway expression. We reconstruct multiomic interaction subnetworks for highly predictable metabolites, and YAP1 signaling emerged as a top global predictor across four omic layers. We prioritize predictive multiomic features for single-cell and spatial metabolomics assays. Top predictors were enriched for synthetic-lethal interactions and synergistic combination therapies that target compensatory metabolic modulators.

Machine Learning↗

PubMed2026

Comparative genomics reveals lineage-associated structural variation and diversification in a barley fungal pathogen.

Leaf rust, caused by Puccinia hordei, is a major barley disease worldwide. Despite repeated shifts in virulence, contrasting reproductive histories, and emerging fungicide insensitivity, the genomic basis of its diversification and adaptation remains poorly understood. In this study, we generated haplotype-resolved, chromosome-level genome assemblies for two isolates with contrasting virulence and analyzed 41 Australian isolates collected over 54 yr (1966-2020), integrating comparative and population genomics, mating-type gene phylogenies, chromosome-specific k-mer profiling, genome-wide copy-number variation (CNV) analysis, and gene-expression analysis. We identified a structurally dynamic chromosome characterized by repeat-associated rearrangements, structural variation, and lineage-associated CNV, representing the first evidence in a rust fungus of chromosome-scale structural diversification of this extent. Population analyses distinguished clonally expanded lineages from recombination-associated lineages, with mating-type gene phylogenies providing further support for lineage differentiation. More recently collected isolates showed increased duplication-associated variation, and CNV boundaries were associated with structural-variant breakpoints. We also identified lineage-associated amplification of Cyp51, with increased copy number associated with higher transcript abundance, supporting a potential role in fungicide adaptation. Overall, our findings highlight structural variation, contrasting reproductive histories, and lineage-associated CNV as important contributors to diversification in P. hordei, providing insights for future rust pathogen surveillance and management strategies.

Cyp51 gene↗

PubMed2026

Case Report: Early infantile drug-resistant epilepsy and gastrointestinal dysmotility associated with a de novo GNAO1 variant and a 16p13.11 microdeletion.

BACKGROUND: The GNAO1 gene is located on the long arm of chromosome 16 (16q13) and is associated with Developmental and Epileptic Encephalopathy 17. It encodes a protein involved in regulating various neurotransmitters and neuronal growth and development. The 16p13.11 microdeletion is a genomic deletion in the p13.11 region of the short arm of chromosome 16, previously reported to predispose individuals to neurodevelopmental disorders. These two regions do not overlap on the chromosome, and the 16p13.11 microdeletion does not involve the GNAO1 gene locus. CASE PRESENTATION: We reports a 1-month-old infant presenting with concurrent GNAO1 (p. G203R) gene variation and a 16p13.11 microdeletion. The patient initially presented with recurrent abdominal distension in the neonatal period, which gradually progressed to include focal seizures and intractable epileptic spasms. Interictal EEG showed multifocal discharges and burst-suppression patterns. Brain MRI revealed no abnormalities. A full gastrointestinal series suggested intestinal obstruction and gastro-esophageal reflux. Improved whole-exome sequencing identified: 1) A single nucleotide variantion: GNAO1, NM_020988.3:exon6: c.607G > A (p.Gly203Arg). 2) A copy number variation: A pathogenic 1.453 Mb deletion in the 16p13.11 region. CONCLUSIONS: We describe the developmental trajectory of an infant with overlapping features of gastrointestinal symptoms and drug-resistant epilepsy, carrying a GNAO1 variant and a 16p13.11 microdeletion. This case provides insights into the contribution of both single nucleotide variant and copy number variant to the phenotype.

16p13.11 microdeletion↗

PubMed1991

Evolutionary transience of hypervariable minisatellites in man and the primates.

Using PCR, two minisatellite loci showing extreme repeat-unit copy-number variation in humans have been characterized in great apes and monkeys. In contrast to humans, minisatellite locus MS32 is monomorphic with only 3-4 diverged repeat units in great apes, Old World and New World monkeys, this organization presumably representing the relatively stable ancestral precursor state of the human hypervariable locus. Similarly, minisatellite MS1 shows extreme repeat-copy-number variability in man compared with low copy number and minimal variability in great apes. Analysis of variant repeat units shows that the 5' and 3' regions of MS1 are relatively stable in great apes and man, and that variability in man is confined to the central region of the minisatellite. In contrast to the great apes, MS1 is highly variable in Old World monkeys. These results, as well as computer simulations of minisatellite evolution based on known mutation rates, show that short minisatellites are stable within the genome, and that the degree of polymorphism at a given locus can change dramatically over a short period of evolutionary time. The ability of hypervariable minisatellites to detect highly informative loci by cross-species hybridization is therefore largely unpredictable.

Animals↗

PubMed2026

Genus-Wide Pan-Genome Analysis of Populus bZIP Transcription Factors with Reanalysis of Public Salt-Stress Transcriptomes.

Basic leucine zipper (bZIP) transcription factors regulate plant development and stress responses, but their genus-wide diversity in Populus remains unclear. We analyzed 19 Populus genomes and retained 1764 bZIP proteins, including 21 independent new loci and four annotation corrections. Of these, 1762 were assigned to 79 orthologous gene groups (OGGs), comprising 43 core, 20 soft-core, 15 shell and one cloud OGG, of which 59 showed copy-number variation. Phylogenetic analysis assigned 74 representative pangenes to 13 subfamilies, with five remaining unclassified and motif patterns differing among subfamilies. Whole-genome duplication (WGD)/segmental duplication accounted for 81.0% of OGG-assigned proteins and contributed predominantly to the conserved component. Although 72.2% of bZIP proteins overlapped a transposable element within the gene body or 2-kb flanks, this proportion was modestly lower than in matched non-bZIP genes, and copy-number-variable OGGs showed no greater TE coverage than invariant OGGs. Among retained homologous comparisons, 97.6% had Ka/Ks ≤ 1, supporting predominant purifying selection. Across the heterogeneous public salt-stress RNA-seq datasets analyzed, no OGG showed a significant, directionally concordant response in at least two Populus taxa. These results reveal a conserved bZIP framework shaped mainly by ancient duplication alongside variable genomic contexts and transcriptional responses.

Populus↗

PubMed1990

Structure and chromosomal mapping of a highly polymorphic repetitive DNA sequence from the pseudoautosomal region of the mouse sex chromosomes.

The pseudoautosomal region of the Mov15 mouse strain is marked by a Moloney murine leukemia provirus. The sequences flanking the Mov15 provirus were molecularly cloned and shown to consist of a tandemly repeated sequence of 31 nucleotides. Copy number variation of this repeat most likely accounts for the polymorphism in the mouse pseudoautosomal region detected with a probe from the flanking sequences. In situ hybridization to metaphase chromosomes showed heavy labeling of the pairing region of the X and Y chromosomes. The repetitive sequence was also found at the subtelomeric region of three autosomes. A similar level of amplification as the one seen on the sex chromosomes seems to be present on chromosomes 9 and 13. Lower copy number appear to be present on chromosome 4.

Animals↗

PubMed2025

PARTAGE: Parallel analysis of replication timing and gene expression.

The human genome is partitioned into functional compartments that replicate at specific times during the S-phase. This temporal program, referred to as replication timing (RT), is co-regulated with the 3D genome organization, is cell type-specific, and changes during development in coordination with gene expression. Moreover, RT alterations are linked to abnormal gene expression, genome instability, and structural variation in multiple diseases, including cancer. However, mechanistic links between RT, large-scale 3D genome architecture, and transcriptional regulation remain poorly understood. A major limitation is that current approaches require the separate profiling of RT and transcriptomes from independent batches of samples, obscuring the complex co-regulation between the epigenome and transcriptome. Here, we developed PARTAGE, a multiomics approach that enables joint profiling of copy number variation (CNV), RT, and gene expression from the same sample, providing a more accurate integrative view of the complex relationships between RT and gene regulation.

Journal Article↗

PubMed2026

The evolution of Ca2+-ATPases across plants with profiles in Rhododendron and the function of key members in alleviating high calcium stress.

Ca2+-ATPase (CAP) is a key Ca2+ efflux protein in plants. Our previous research suggests that CAPs may play a crucial role in the adaptation of rhododendrons to high calcium environments. However, the evolution, variation, characteristic expression, and subfunctionalization of this gene family in Rhododendron remain unknown. Through the analysis of pan-genomes and pan-transcriptomes, we elucidated the systematic evolution of CAPs in plants, as well as their characteristic expression patterns in Rhododendron. During the evolutionary process from lower to higher plants, CAPs can be divided into six clades and exhibit structural conservation. CAPs have emerged and differentiated in lower plants such as algae, and they have undergone significant amplification in Eudicots plants like rhododendrons. Three Rhododendron species (Rhododendron bailiense, R. delavayi, and R. irroratum) located in the karst province of Guizhou in Southwest China exhibit the highest copies of CAPs, suggesting a strong association between CAP copy number variation and habitat, particularly in high calcium environments. Through multiple transcriptome analyses, we revealed that CAPs are induced under various environmental/developmental conditions (e.g. karst environments, high altitude, early flower development, hormones, etc.). Co-expression network analysis highlighted key members of calcineurin B-like protein (CBL) and CBL-interacting protein kinases (CIPK) that are associated with the high expression of CAPs. Experimental validation demonstrated that CAPb1 and CAPd1 significantly alleviate high calcium stress, and the CAPb1-CIPK1-CBL1 and CAPd1-CIPK2-CBL1 modules can further enhance the alleviation. These findings provide new insights into the evolution, characteristic expression, and function of CAPs, as well as new perspectives on the high calcium adaptability of rhododendrons.

Journal Article↗

PubMed2025

Experimental evolution of a pathogen confronted with innate immune memory increases variation in virulence.

Understanding the drivers and mechanisms of virulence evolution is still a major goal of evolutionary biologists and epidemiologists. Theory predicts that the way virulence evolves depends on the balance between the benefits and costs it provides to pathogen fitness. Additionally, host responses to infections, such as resistance or tolerance, play a critical role in shaping virulence evolution. But, while the evolution of pathogens has been traditionally studied under the selection pressure of host adaptive immunity, less is known about their evolution when confronted to simpler and less effective forms of immunity such as immune priming. In this study, we used a well-established insect model for immune priming - red flour beetles and their bacterial pathogen Bacillus thuringiensis tenebrionis - to test how this form of innate immune memory drives the pathogen evolution. Through controlled experimental evolution of the pathogen in primed versus non-primed hosts, we found no change in average virulence after eight selection cycles in primed host. Nonetheless, we observed a notable rise in the variability of virulence, defined as the ability to kill hosts, among independent pathogen lines that evolved in primed hosts, and the bacteria were unable to develop resistance to host priming. Whole genome sequencing revealed increased activity in the bacterial mobilome (prophages and plasmids). Expression of the Cry toxin - a well-known virulence factor - was linked to evolved differences in copy number variation of the cry-carrying plasmid, though this did not correlate directly with virulence. These findings highlight that innate immune memory can drive variability in pathogen traits, which may favor adaptation to variable environments. This underscores the need to consider pathogen evolution in response to innate immune memory when applying these mechanisms in medicine, aquaculture, pest control, and insect mass production.

Animals↗

PubMed2026

Spatial Total RNA Sequencing of Formalin-Fixed Paraffin-Embedded Tissue by spRandom-seq.

The molecular pathogenesis of infectious diseases and cancer is orchestrated by nanoscale of host and microbial RNA transcripts within the tissue microenvironment. Nevertheless, spatially resolving the comprehensive transcriptional landscape within complex clinical tissues, like formalin-fixed paraffin-embedded (FFPE) specimens, still poses a formidable challenge. Here, we present spRandom-seq, a random primer-based spatial total RNA sequencing technology designed to spatially resolve complete transcriptomes from host, bacteria, and even nanoscale viruses in FFPE tissues. Capitalizing on the random primer design, our technology not only facilitated the discovery of specific lncRNAs and alternative splicing events in mouse brain and olfactory bulb, but also delineated pronounced spatial heterogeneity in clinical FFPE sections-across distinct tumor regions in breast cancer and microbial infection sites in Klebsiella pneumoniae-infected tissues. Importantly, integrated analysis of host and viral RNAs in FFPE samples from hepatitis B virus (HBV)‑positive hepatocellular carcinoma (HCC) demonstrated that complement and coagulation pathways were specifically activated across expansive HBV‑infected tumor areas, which also exhibited an increased burden of copy number variations (CNVs). Owing to its compatibility with existing spatial transcriptomics platforms and minimal operational complexity, spRandom-seq represents a practical and scalable approach for clinical pathology applications and infection diagnostics.

Paraffin Embedding↗