SEARCH · PubMed Health

Results for “Developmental transcriptome”

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 55 records · Page 3Linked to original sources

PubMed2026

De novo chromatin remodelling variants in sporadic Chiari 1 malformation.

Chiari 1 malformation (CM1) is the most common congenital malformation of the human hindbrain. Although prior studies have implicated chromatin-remodeling genes in CM1, the de novo genetic architecture and underlying neurodevelopmental mechanisms remain incompletely defined. To investigate the molecular genetics of a novel familial form of CM1 linked with syringomyelia and tethered cord and determine whether rare, damaging de novo variants (DNVs) contribute to sporadic CM1 risk with gene- and pathway-level resolution, we performed whole-exome sequencing in an ultra-rare multigenerational family with CM1 and associated spinal pathology, and in the largest assembled trio-based cohort to date, comprising 1,585 proband-parent trios with sporadic, idiopathic CM1 (2017-2025). The comparison cohort included 1,798 unaffected control siblings. Clinical phenotyping was by systematic medical record review. Structural domain mapping, in silico modeling, and integration with single-cell transcriptomic data from developing human cerebellum was conducted to assess biological plausibility. A heterozygous loss-of-function variant in CHD3 segregated with CM1 and syringomyelia in a multigenerational family. In the trio-based cohort, rare protein-altering DNVs were significantly enriched across multiple chromodomain helicase DNA-binding (CHD) genes, including CHD1, CHD3, CHD4, and CHD8, exceeding gene-specific mutation expectations (protein-damaging variants: P = 1.3 × 10-9; predicted loss-of-function variants: P = 8.6 × 10-5). CHD1 contained two pathogenic DNVs (p.A999D and p.E984K). CHD4 (p.D744N, p.T1813P, and p.I1102T) and CHD8 (p.R1402X, p.R1472X, and p.R2035X) each contained three new DNVs. Variants clustered within conserved ATPase, helicase, and chromodomain regions essential for chromatin remodeling, and these patients frequently had comorbid developmental delay and related neurodevelopmental features. Single-cell transcriptomic analyses demonstrated enrichment in Purkinje cells and inhibitory neurons of midgestational cerebellum, where CHD gene products form a coherent chromatin-regulatory network. Rare, large-effect DNVs that disrupt chromatin-remodeling programs contribute to sporadic CM1, implicating genetically encoded dysregulation of cerebellar development as a central disease mechanism. Exome sequencing may complement surgical evaluation of children with sporadic CM1, particularly when accompanied by neurodevelopmental concerns, informing prognosis and family counseling.

de novo variants↗

PubMed2026

Comparative Transcriptomics Reveals Shared Downstream Pathways in Craniofacial Pathology.

Treacher Collins syndrome and Nager syndrome are craniofacial developmental disorders caused by defects in ribosome biogenesis and RNA splicing, respectively, yet they exhibit overlapping abnormalities affecting neural crest cell-derived craniofacial structures. To investigate shared downstream pathogenic mechanisms, we performed a comparative transcriptomic analysis of zebrafish polr1c and sf3b4 mutant models from our previous studies. Comparative analysis identified 17 shared differentially expressed genes (DEGs) between polr1c and sf3b4 mutants, with the majority of shared genes dysregulated in the same direction, indicating a coordinated rather than random transcriptional response. Gene ontology analysis identified ATP-dependent protein folding chaperone activity as the only shared molecular function, driven in part by upregulation of hsp90aa1.2, indicating a common proteostasis response. Because chaperone activity is linked to extracellular matrix (ECM) protein processing, we cross-referenced DEGs from both mutants against the curated zebrafish matrisome. Three of the 17 shared DEGs (serpinh1b, il11a, and lepa) were matrisome-associated and upregulated in both mutants. Serpinh1b, a collagen-specific chaperone, was strongly expressed in craniofacial cartilage and mesenchymal populations during pharyngeal arch development and exhibited nearly identical fold changes in both mutants. Il11a is of particular interest because its receptor, IL11RA, is known to be associated with human craniosynostosis, suggesting potential relevance to craniofacial development. Together, it is possible to hypothesize that shared chaperone-associated transcriptional changes, together with altered ECM-related gene expression, may contribute to polr1c- and sf3b4-associated craniofacial disorders, warranting further functional validation.

Extracellular Matrix↗

PubMed2026

Primary Tumor Epigenetic and Transcriptomic Alterations Associated with Nodal Burden and Metastatic Risk in ER+/HER2- Breast Cancer.

De-escalation of axillary surgery has resulted in the loss of pathologic nodal information, yet the extent of lymph node involvement remains an important determinant of treatment decisions in estrogen receptor-positive (ER+)/HER2- disease. We examined whether primary tumors differed molecularly according to the extent of this regional dissemination. Genome-wide DNA methylation profiling of primary ER+/HER2- tumors from 47 patients with pN1 (n = 29) vs. >pN1 (n = 18) disease showed differences concentrated at promoters of developmental and cell-adhesion genes. By integrating methylomes with transcriptomes from the TCGA-BRCA cohort (n = 148) and clinical outcomes from KM Plotter (RFS, n = 1154; OS, n = 442; DMFS, n = 423), we identified four genes (ARL10, RIC3, CXCL14, KCNH2) showing concordant molecular and clinical associations, from which we derived the Lymph-node Involvement Outcome Numerator (LION) score. Lower LION scores were observed in metastatic lesions from the AURORA US cohort (n = 45). In SCAN-B (n = 3969), lower scores were associated with shorter distant recurrence-free intervals (HR = 0.38; 95% CI 0.23-0.62); this association persisted after adjustment for age, nodal and tumor category but was lost after adjustment for histological grade (HR = 0.83; 95% CI 0.48-1.44), indicating that the score and grade capture overlapping biology. These findings suggest that primary tumors already display coordinated epigenetic and transcriptional alterations associated with the extent of metastatic dissemination.

Humans↗

PubMed2025

Impaired natural killer cell maturation in lung adenocarcinoma driven by FABP4 and SPON2 downregulation through disrupted lipid metabolism.

BACKGROUND: Although natural killer (NK) cells play a crucial role in antitumor immunity, the metabolic changes driving their dysfunction in lung adenocarcinoma remain poorly understood. This study investigates how these metabolic modifications impact NK cell function within the lung adenocarcinoma microenvironment. METHODS: A total of 13 pairs of lung adenocarcinoma samples were obtained from The Cancer Genome Atlas. Differential gene expression, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, and single-cell metabolic quantification analyses were used to characterize the transcriptomic, pathway, and metabolic signatures of NK cells. The developmental trajectory was reconstructed via pseudotime analysis. The fatty acid-binding protein 4 (FABP4) and spondin2 (SPON2) expression was examined using immunofluorescence (IF) and immunohistochemistry (IHC) in patients with lung adenocarcinoma. In NK cells with FABP4 downregulation, FABP4 function was analyzed using antibody-independent cell-mediated cytotoxicity assays, flow cytometry (FCM), and liquid chromatography-mass spectrometry. RESULTS: The number of NK cells was significantly decreased in the lung adenocarcinoma microenvironment. FABP4 and SPON2 expression was significantly lower in NK cells within tumor tissues than in the adjacent tissues. FABP4 expression was significantly lower in tumor tissues than in the adjacent tissues, whereas no significant difference in SPON2 expression was observed. The cytotoxic function of NK cells with decreased FABP4 levels was impaired. Non-targeted lipid metabolism analysis indicated that differentially expressed lipids in NK cells with low FABP4 levels were functionally enriched in the glycerophospholipid metabolism pathway compared to those in normal NK cells. CONCLUSIONS: The study findings present new evidence showing that low FABP4 and SPON2 gene expression may impair NK cell maturity by affecting lipid metabolism in lung adenocarcinoma. These results provide a new perspective on restoring immune function in patients with lung cancer.

FABP4↗

PubMed2025

Combined somatic mutation and transcriptome analysis reveals region-specific differences in clonal architecture in human cortex.

The human cerebral cortex is specialized into regions, but little is known about how human cellular lineages shape cortical regional variation and neuronal cell-type distribution during development. Here, we map single-cell lineages of human cortical regions and neuronal subtypes using >1,000 somatic single-nucleotide variants (sSNVs) identified from deep bulk whole-genome sequencing and analyzed over 25 regions and >72,000 single cells. In the fronto-parietal cortex, sSNVs are rarely restricted, marking neuron-generating clones that disperse into neighboring regions. In contrast, the primary visual cortex harbors 30%-70% more sSNVs than the neighboring secondary visual cortex. Clones at this border exhibit more restricted dispersion, suggesting late developmental lineage segregation. Single-nucleus sSNV and whole-transcriptome analysis reveal glutamatergic neuron clones with modest regional restrictions that share low-mosaic sSNVs with some GABAergic neurons, suggesting a recent dorsal cortical progenitor. Our analysis reveals human-specific cortical lineage patterns, regional differences in clonal patterns, and late divergence of some glutamatergic/GABAergic lineages.

Humans↗

PubMed2026

Beyond the gene: isoform diversity as a key contributor to human brain disorders.

The human brain exhibits exceptional transcriptomic complexity, with alternative splicing, promoter usage, and polyadenylation generating extensive transcript-isoform diversity. Isoform dysregulation is increasingly implicated in neurodevelopmental and psychiatric disorders (NPDs), yet the landscape, function, and genetic regulation of brain isoforms remain poorly understood due to limitations of short-read RNA sequencing. Advances in long-read sequencing (LR-seq) enable scalable full-length transcriptome profiling with single-cell and spatial resolution across developmental stages. Here, we review recent progress in isoform discovery, quantification, functional annotation, and genetic regulation, highlighting emerging links to human neurodevelopment and disease. LR-seq studies have uncovered tens of thousands of previously unannotated brain isoforms, with neuronal maturation characterized by increased exon inclusion and progressive 3' untranslated region (3' UTR) lengthening. Isoform-resolved genetic mapping outperforms gene-level analyses for NPD gene discovery and mechanistic interpretation. We argue that a shift from gene-centric to isoform-centric frameworks is essential to fully capture regulatory complexity in human neurogenetics. Together, these advances establish isoform diversity as a fundamental yet underappreciated axis of brain gene regulation and a key entry point for dissecting NPD biology.

Humans↗

PubMed2026

A recurrent CCDC82 frameshift variant associated with syndromic neurodevelopmental disorder in a consanguineous Pakistani family.

BACKGROUND: Intellectual disabilities (IDs) are part of neurodevelopmental disorders (NDDs) and are genetically heterogeneous conditions characterized by impairments in cognition, learning, and adaptive functioning. Despite advances in gene discovery, many individuals, particularly those from understudied populations, remain without a molecular diagnosis. Recent reports implicate CCDC82 (HGNC: 26282) as an autosomal recessive ID gene, although the phenotypic spectrum and biological context remain incompletely defined. METHODS: Exome sequencing (ES) was performed in a consanguineous Pakistani family (PKMR06A) with four affected individuals presenting with moderate to severe ID. Variant segregation was confirmed by Sanger sequencing. In silico analyses, including pathogenicity prediction, protein structural modeling, and domain intolerance assessment, were used to evaluate the functional consequences of the identified variant. Spatiotemporal gene expression patterns were examined using bulk and single-cell human brain transcriptomic datasets. RESULTS: Clinically, affected individuals of family PKMR06A presented with early childhood global developmental delay, speech delay, hypotonia, gait abnormalities, spasticity, and mild facial dysmorphism. Genetic screening revealed a recurrent rare homozygous frameshift variant in CCDC82 (NM_024725.4): c.373del; p.(Asp125Ilefs*6), segregating with disease in all available affected individuals of the family. The identified c.373del variant was absent from the gnomAD database and was classified as pathogenic (PVS1, PM2, and PP1) based on ACMG/AMP criteria. The c.373del variant is predicted to introduce a premature termination codon, p.(Asp125Ilefs*6), leading to deletion of essential coiled-coil domains from the encoded protein, supporting a loss-of-function mechanism. In silico, transcriptomic analyses demonstrated preferential CCDC82 expression during prenatal human brain development, providing developmental context for the neurodevelopmental phenotype associated with the identified truncating variant. CONCLUSIONS: This study expands the mutational landscape of CCDC82 and provides additional clinical and molecular evidence supporting its role in autosomal recessive NDD. The findings reinforce the importance of CCDC82 in human neurodevelopment and highlight the value of genomic investigation in underrepresented populations.

Autosomal recessive↗

PubMed2026

Prominent Movement Disorders in RNU2-2-Related Spliceosomopathy.

Pediatric movement disorders often overlap with neurodevelopmental diseases, suggesting shared molecular mechanisms. Variants in small nuclear RNA (snRNA) genes encoding spliceosome components have recently been associated with neurodevelopmental disorders, termed "RNUopathies." We analyzed genome sequencing data from 14 patients with undiagnosed pediatric movement disorders for pathogenic variants in snRNA genes. We identified recurrent de novo RNU2-2 variants (n.35A > G and n.4G > A) in two patients with intellectual disability, epilepsy, and hyperkinetic movement disorders. RNA sequencing of fibroblasts in one patient showed no characteristic transcriptomic signature. Spliceosomopathies should be considered in neurodevelopmental disorders and developmental and epileptic encephalopathies with hyperkinetic features.

Humans↗

PubMed2026

Firemaster 550 differentially alters gene expression underlying synaptic function in amygdala of prairie voles after gestational or lactational exposure.

Neurodevelopmental disorders often share similar behavioral diagnostic criteria including socioemotional and cognitive deficits. The prairie vole is a uniquely suitable model to study these deficits because they demonstrate strong social affiliation, bi-parental care, and partner attachment. Previously, we have shown that developmental exposure to the flame-retardant mixture Firemaster 550 (FM 550) impairs socioemotional behavior in the prairie vole and alters underlying neuroanatomy and function. However, the mechanisms for impaired pair bonding in males and increased anxiety in females remain unknown, along with the specific critical window(s) of vulnerability. Herein, we exposed prairie vole dams to FM 550 during gestation or lactation, and performed bulk RNA-seq on the amygdala, a hub of socioemotional processing, in their adult offspring. Two mathematically orthogonal methods were utilized for analysis, a linear statistical method and an ensemble machine learning method, incorporating sex as a biological variable. Gene ontology (GO) pathway analysis was performed following both and results compared to identify potential mechanisms of toxicity. GO results indicated consistent expression changes in the Synapse cellular component in all conditions, and implicated glutamatergic signaling specifically. Additionally, gestational exposure (GE) altered genes underlying modulation of synaptic transmission and neural development, while lactational exposure (LE) impacted genes underlying synaptic plasticity, axon guidance, and mitophagy. Machine learning identified disruption of endocrine system development, regulation of biosynthetic processes in GE animals, and suppression of various neuroinflammatory genes across multiple groups. Finally, we performed RNA expression analysis using Nanostring and demonstrated stronger correlation with the differentially expressed genes (DEG) of interest in females than males. Overall, this study demonstrates both the intersecting and distinct impacts of FM 550 exposure on amygdalar gene expression depending on sex and timing of exposure.

Animals↗

PubMed2026

Social evolution and diminished olfactory function in larval honey bees.

Social evolution made larval honey bees dependent on adult colony members for feeding; they are confined to cells in waxen honeycombs and visited about 100 times per day by adult "nurse" bees. Based on organismal resource conservation theory, we predicted larvae have diminished olfactory capabilities at both the molecular and behavioral levels. Consistent with theory, larvae expressed very low levels of Orco, an essential gene for olfactory receptor (OR) function. By contrast, they showed higher expression of Ir25a, essential for other forms of sensory perception including gustation. Also consistent with theory, behavioral assays demonstrated that larvae cannot find food via olfaction, suggesting they use taste for feeding. By contrast, it is known that adult honey bees use OR-based olfaction extensively for a variety of behavioral functions, and the honey bee genome contains many OR-encoding genes. Comparative transcriptomic analyses of social and nonsocial insects suggest that this developmentally regulated suppression of olfactory function is related to social evolution, especially systems of offspring care.

Animals↗

PubMed2026

Co-option of stomata in the convergent evolution of fern nectaries.

Understanding the origin of new structures is a central goal of evolutionary biology. In many instances, novel phenotypes arise through heterotopy: the expression of a structure in a new location. Using bracken fern (Pteridium aquilinum) as a model, we combine genomics, transcriptomics and metabolomics to begin to explore the origin and developmental routes in the convergent evolution of ant-enticing nectaries. We observe that P. aquilinum does not exclusively express flowering plant 'nectary genes' during nectary development. Rather, this fern builds nectaries through co-option of stomata. Specifically, P. aquilinum heterotopically expresses canonical angiosperm stomatal regulatory genes, leading to stomatal development in novel positions along the petiole. These non-laminar stomata were co-opted for nectar secretion through the expression of putative sugar transport genes, forming secretory nectarostomata. This work provides two advances in our understanding of nectary evolution and the origin of complex structures. First, heterotopic expression of stomata, and later exaptation, represents one realized developmental mechanism for the evolution of nectar glands. Second, while there are many routes to nectary evolution, nectarostomata development is a repeatable path that has evolved in ferns and flowering plants, representing an impressive case of convergent evolution through the same developmental mechanism, despite over 400 million years of divergent history.

Plant Stomata↗

PubMed2024

A time-resolved single-cell roadmap of the logic driving anterior neural crest diversification from neural border to migration stages.

Neural crest cells exemplify cellular diversification from a multipotent progenitor population. However, the full sequence of early molecular choices orchestrating the emergence of neural crest heterogeneity from the embryonic ectoderm remains elusive. Gene-regulatory-networks (GRN) govern early development and cell specification toward definitive neural crest. Here, we combine ultradense single-cell transcriptomes with machine-learning and large-scale transcriptomic and epigenomic experimental validation of selected trajectories, to provide the general principles and highlight specific features of the GRN underlying neural crest fate diversification from induction to early migration stages using Xenopus frog embryos as a model. During gastrulation, a transient neural border zone state precedes the choice between neural crest and placodes which includes multiple converging gene programs. During neurulation, transcription factor connectome, and bifurcation analyses demonstrate the early emergence of neural crest fates at the neural plate stage, alongside an unbiased multipotent-like lineage persisting until epithelial-mesenchymal transition stage. We also decipher circuits driving cranial and vagal neural crest formation and provide a broadly applicable high-throughput validation strategy for investigating single-cell transcriptomes in vertebrate GRNs in development, evolution, and disease.

Animals↗

PubMed2026

Natural leaf shape variation reveals diverse transcriptional targets of GmJAG1 during soybean leaf development.

The JAGGED transcription factor family regulates lateral organ development across angiosperms. In soybean (Glycine max Merr.), a D9H mutation in the EAR repression motif of GmJAG1 causes a narrow-leaflet phenotype and explains over 70% of phenotypic variance in leaf shape. Because this mutation does not affect the zinc finger DNA-binding domain, both alleles bind identical targets but differ in repressor recruitment. Previous studies mapped GmJAG1 binding sites, but the functional targets controlling leaf morphology are uncharacterized. Here, we used comparative transcriptomics across four soybean genotypes with contrasting leaf shape, spanning a developmental time series from shoot apex to mature leaf, and identified 1567 putative candidate target genes. GmJAG1 expression was confined to the shoot apex, yet 99.1% of candidate targets maintained differential expression throughout development. We found that neither Kip-Related Protein (KRP) cell cycle inhibitors nor Cyclin-Dependent Kinases (CDKs) showed differential expression despite binding evidence in Arabidopsis. However, D-type cyclins were upregulated in narrow-leaf genotypes suggesting that soybean GmJAG1 acts through cyclin-mediated rather than KRP-mediated cell cycle regulation described in Arabidopsis- a divergence in regulatory logic between the two species. Pathway analysis revealed enrichment of auxin (1.8-fold, P = 0.02) and salicylic acid (fourfold, P = 0.016) genes among JAG1D9H targets. Filtering by differential expression, binding data, phenotype correlation, and co-expression network membership identified 79 high-confidence targets, including orthologs of NPH3 (phototropin-mediated leaf flattening), MIK2 (cell wall integrity sensing), RD22 (ABA-responsive stress signaling), and SCL23 (GRAS transcription factor in bundle sheath development). These candidates provide targets for functional validation and breeding in legumes.

Glycine max↗

PubMed2026

Studies on functional differentiation of xpr1a and xpr1b genes in zebrafish.

Xenotropic and polytropic retrovirus receptor 1 (XPR1) is known to be involved in various biological processes, including phosphate homeostasis, cellular signaling, brain and vascular mineralization, whereas its specific contribution to bone development remains incompletely characterized. Due to genome duplication in teleosts, zebrafish Danio rerio possess two paralogous genes of XPR1 namely xpr1a and xpr1b, whose functional divergence remains unclear. The amino acid sequence similarity between zebrafish xpr1a and xpr1b was 83.26%. In situ hybridization demonstrated overlapping localization in the head and spinal cord at 24-48 hpf, while diverged by 72 hpf, with xpr1a becoming restricted to the head while xpr1b persisted in both regions. CRISPR/Cas9 was used to generate xpr1a and xpr1b mutants. The xpr1a mutants are comparatively healthy, viable but with mild growth reduction, whereas the xpr1b mutants display high mortality, reduced body length and severe vertebral deformities. Interestingly, all the double mutants died at the embryonic stage. Moreover, to further investigate the molecular and regulatory mechanisms, we conducted comparative transcriptome analysis on bone and brain tissues from xpr1b+/+ and xpr1b-/- zebrafish. In bone tissue, 6749 DEGs were identified, comprising 3846 upregulated and 2903 downregulated genes. These DEGs were mainly enriched in the MAPK signaling pathway, Wnt signaling pathway, cysteine and methionine metabolism, and ECM-receptor interaction. RT-qPCR validated results showed that seven osteogenesis-related genes (col1a1a, sp7, runx2b, col1a2, col1a1b, alp1 and entpd5), and two phosphate homeostasis related genes (slc20a2 and pdgfba), which are essential for skeletal mineralization and phosphate homeostasis, exhibited significantly downregulated expression in bone tissue of xpr1b mutant zebrafish. These results highlight the pivotal role of xpr1b in regulating skeletal mineralization and phosphate metabolism, thereby elucidating the functional specialization of XPR1 paralogs while providing a theoretical basis for understanding bone developmental mechanism in teleost vertebrates.

Animals↗

PubMed2026

SpliceHarmonization: an integrated method for identifying RNA splicing events in therapeutics for splicing modulation.

MOTIVATION: Splicing, a critical co-transcriptional process in eukaryotes, enhances transcriptome diversity by generating isoforms specific to cell types, tissues, or developmental stages. Recent advancements in splicing modulators have opened new avenues for targeting previously undruggable genes by inducing significant perturbations in splicing events. These developments underscore the need for comprehensive methods to accurately identify and compare splicing events. While several tools have been developed to detect local splice variants, inconsistencies across methods remain a significant challenge. To address this, we present SpliceHarmonization, an integrated approach that combines the strengths of rMATS, LeafCutter, and MAJIQ, enabling robust and reliable splicing analysis with event type annotations. RESULTS: In a comprehensive evaluation using diverse simulated datasets, SpliceHarmonization streamlined and standardized the outputs from three detection methods into a unified format, thereby improving splicing detection with event type annotation and outperforming individual methods. By integrating the outputs from rMATS, LeafCutter, and MAJIQ, our approach not only enhanced identification of a wide range of splicing events but also effectively mitigated method-specific discrepancies. This integration led to an accuracy exceeding 0.8 and a recall of up to 0.5, with an observed increase in AUC of up to 10%. Furthermore, SpliceHarmonization demonstrated high sensitivity in detecting low-abundance and complex splicing events, providing annotations including genomic coordinates and event type. AVAILABILITY AND IMPLEMENTATION: SpliceHarmonization is available at https://github.com/interactivereport/SpliceHarmonization.

RNA Splicing↗

PubMed2026

Natural variants of CsSHN1 orchestrate a temporal regulatory cascade driving fruit skin netting in cucumber.

Fruit skin netting (russeting, Rs) forms when epidermal microcracks are sealed by a suberized periderm, reducing marketability. We previously identified the Rs locus (CsSHN1), which encodes an AP2/ERF transcription factor, as a major determinant of cucumber skin netting, but how fruit growth is temporally coupled to periderm formation remains unclear. Here, we integrated population genomics, time-series multiomics, DNA affinity purification sequencing (DAP-seq), and transgenic assays to decode the CsSHN1-mediated regulatory network. Six functionally relevant CsSHN1 variants were identified across 325 cucumber accessions. Allele distribution and selective sweep analyses revealed breeding-driven selection for smooth fruit skin. Overexpression of a netted allele in a smooth background induced epidermal fissures, altered cell geometry, and increased fruit size, demonstrating a dosage-sensitive effect. Time-series transcriptomics and metabolomics of near-isogenic lines (NILs) defined 3 developmental phases of netting: early suppression of lignin and trehalose genes preceding cracks, growth-driven fissuring accompanied by cell-wall remodeling and defense activation, and maturation-stage cell-wall degradation with strong induction of ligno-suberin biosynthesis. Across the cucumber genome, DAP-seq identified approximately 8,000 in vitro CsSHN1 binding sites. These binding sites were significantly enriched for the GCC-box motif and included genes involved in cutin and suberin biosynthesis. Together, these results show that CsSHN1 orchestrates fruit skin netting through a growth-coupled temporal regulatory cascade, providing a mechanistic framework for manipulating fruit epidermal properties.

Cucumis sativus↗

PubMed2026

Genome-Wide Characterization of the ZIP Transporter Family in Sea Island Cotton (Gossypium barbadense L.) and Expression Profiling Under Heavy Metal and Pathogen Stresses.

G. barbadense represents an indispensable germplasm resource for high-quality textile fiber and disease resistance; nevertheless, systematic information regarding its ZRT/IRT-like protein (ZIP) gene family remains limited. Here, a total of 46 GbZIP genes were identified across the G. barbadense genome. Comprehensive bioinformatic investigations revealed uneven chromosomal distribution and confirmed that segmental/whole-genome duplications, supplemented by localized tandem duplications, drove family expansion. Members clustered within the same phylogenetic clades shared conserved motif organization and gene architecture, while promoter regions harbored abundant cis-acting elements associated with phytohormone and stress signaling. Transcriptome profiling indicated distinct expression patterns across vegetative/reproductive tissues, fiber and ovule developmental stages, and diverse abiotic stress conditions (cold, hot, drought, and salt). Quantitative Real-Time PCR (qRT-PCR) further validated that several GbZIP candidates exhibited temporal expression variations upon exposure to cadmium toxicity, V. dahliae infection, and combined Cd-V. dahliae stress. Specifically, GbZIP13, GbZIP18, GbZIP27, and GbZIP36 displayed prominent broad-spectrum responses to all three stress conditions, whereas GbZIP16, GbZIP29, and GbZIP30 showed stress-specific regulatory divergence. Overall, this study aims to systematically analyze the evolutionary characteristics and expression patterns of the GbZIP family, and to specifically evaluate the response differences under Cd stress, V. dahliae stress, and combined stress, in order to identify potential key candidate genes.

Gossypium barbadense↗

PubMed2026

A conserved partner-switching system controls terminal differentiation in multicellular cyanobacteria.

UNLABELLED: Canonical partner-switching systems (PSSs) regulate sigma factor activity through reversible phosphorylation, but their established roles have been largely limited to stress responses and sporulation in Firmicutes. Whether this regulatory mechanism also controls developmental cell fate decisions in other bacterial phyla has remained unknown. Here, we identify a canonical PSS that governs heterocyst differentiation in the multicellular cyanobacterium Anabaena sp. PCC 7120. This system comprises the anti-sigma factor All2284 (NfsS) and the anti-anti-sigma factor All2283 (NfsR). Structural predictions and biochemical assays showed that NfsS phosphorylates NfsR on a conserved serine residue, whereas bacterial two-hybrid and co-purification assays demonstrated that NfsS binds the developmental sigma factors SigC and SigE. Deletion of nfsR abolished heterocyst formation and diazotrophic growth, and transcriptomic analysis revealed broad failure to induce late heterocyst genes, including nitrogen fixation functions such as nifHDK and fdxH. Phylum-wide comparative genomics further showed that PSS genes and putative functional clusters are strongly enriched in filamentous and heterocyst-forming taxa, indicating an association between the expansion of these signaling modules and the emergence of multicellularity and developmental specialization. Together, these findings establish a PSS as a direct regulator of terminal cell differentiation in a gram-negative bacterium and reveal partner switching as a conserved regulatory principle linking environmental signaling to developmental fate in a major bacterial phylum. IMPORTANCE: While partner-switching systems are classically associated with stress responses and sporulation control in Firmicutes, whether this regulatory logic governs developmental decisions in other bacterial phyla has remained unknown. Here, we establish that a related partner-switching mechanism operates in a distinct bacterial lineage, the cyanobacteria, where it controls a major developmental transition involving terminal cell differentiation. Specifically, we show that a phosphorylation-dependent checkpoint involving the anti-sigma factor NfsS and the anti-anti-sigma factor NfsR directly regulates heterocyst formation. Disruption of this switch abolishes cell differentiation and diazotrophic growth, revealing that this system is an obligate gatekeeper for terminal differentiation. Conceptually, these findings substantially extend the known functional repertoire of partner-switching circuits: rather than controlling stress adaptation or spore dormancy, this module has been co-opted to govern a complex, multicellular developmental program in an organism that underpins global carbon and nitrogen cycles. This work, therefore, establishes a new paradigm for phosphorylation-based control of developmental sigma factors and provides a tractable model for dissecting how conserved signaling modules are rewired to drive lineage-specific innovations across the bacterial domain.

cell differentiation↗