Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “transcriptome evolution”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

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

At least 37 records · Page 2Linked to original sources
PubMed2024

Molecular mechanisms of the specialization of human synapses in the neocortex.

Synapses of the neocortex specialized during human evolution to develop over extended timescales, process vast amounts of information and increase connectivity, which is thought to underlie our advanced social and cognitive abilities. These features reflect species-specific regulations of neuron and synapse cell biology. However, despite growing understanding of the human genome and the brain transcriptome at the single-cell level, linking human-specific genetic changes to the specialization of human synapses has remained experimentally challenging. In this review, we describe recent progress in characterizing divergent morphofunctional and developmental properties of human synapses, and we discuss new insights into the underlying molecular mechanisms. We also highlight intersections between evolutionary innovations and disorder-related dysfunctions at the synapse.

Humans↗
PubMed2025

Evolutionary fingerprints of epithelial-to-mesenchymal transition.

Mesenchymal plasticity has been extensively described in advanced epithelial cancers; however, its functional role in malignant progression is controversial1-5. The function of epithelial-to-mesenchymal transition (EMT) and cell plasticity in tumour heterogeneity and clonal evolution is poorly understood. Here we clarify the contribution of EMT to malignant progression in pancreatic cancer. We used somatic mosaic genome engineering technologies to trace and ablate malignant mesenchymal lineages along the EMT continuum. The experimental evidence clarifies the essential contribution of mesenchymal lineages to pancreatic cancer evolution. Spatial genomic analysis, single-cell transcriptomic and epigenomic profiling of EMT clarifies its contribution to the emergence of genomic instability, including events of chromothripsis. Genetic ablation of mesenchymal lineages robustly abolished these mutational processes and evolutionary patterns, as confirmed by cross-species analysis of pancreatic and other human solid tumours. Mechanistically, we identified that malignant cells with mesenchymal features display increased chromatin accessibility, particularly in the pericentromeric and centromeric regions, in turn resulting in delayed mitosis and catastrophic cell division. Thus, EMT favours the emergence of genomic-unstable, highly fit tumour cells, which strongly supports the concept of cell-state-restricted patterns of evolution, whereby cancer cell speciation is propagated to progeny within restricted functional compartments. Restraining the evolutionary routes through ablation of clones capable of mesenchymal plasticity, and extinction of the derived lineages, halts the malignant potential of one of the most aggressive forms of human cancer.

Animals↗
PubMed2026

Genomic characterization of aggressiveness in pituitary neuroendocrine tumors.

BACKGROUND: Aggressive evolution of PitNETs is rare; metastatic spread is even more. Defining aggressiveness and malignancy is challenging, subsequently hard to predict, and to understand. The aim was to provide a molecular definition of aggressiveness using genomic approaches. METHODS: PitNETs from 206 patients were included. Associations between 9 clinicopathological features of aggressiveness and PitNETs' omics were explored. Omics included transcriptome, DNA methylation, chromosomal alterations, and mutations. Clonal tumor evolution was monitored in 7 patients. RESULTS: Among the 9 clinicopathological features of aggressiveness, only rapid progression, progression after radiotherapy, Ki67/MIB1 proliferation index ≥10%, temozolomide treatment, metastases, and specific death were associated with specific omics signatures, while tumour maximal diameter ≥40 mm, cavernous, and sphenoid invasion were not. The omic signatures associated with these features of aggressiveness overlapped but remained distinct between corticotroph and mammo-somato-thyrotroph lineages. For each lineage, a common signature of aggressiveness was identified, associating a proliferative transcriptome signature and DNA hypermethylation. Alterations in specific genes were associated with aggressive features, including a novel PitNET gene, LRP1B, and known cancer genes (TP53, CDKN2A), while USP8 and GNAS alterations were not. Integration of gene alterations with methylome and transcriptome signatures isolated a subset of molecularly aggressive PitNETs. Molecular signatures were stable during the course of the disease, despite evolution toward aggressiveness and potential clonal divergence. CONCLUSION: This systematic analysis of clinicopathological features of aggressiveness using an integrated multiomic approach establishes a histomolecular definition of aggressiveness in PitNETs. Prospective cohort studies are needed to validate these molecular signatures and establish their prognostic value.

Humans↗
PubMed2026

An integrated single-cell and spatial transcriptomic atlas of thyroid cancer progression identifies prognostic fibroblast subpopulations.

Although well-differentiated thyroid carcinoma (WDTC) is characterized by a robust treatment response, aggressive subtypes, such as anaplastic thyroid carcinoma (ATC), remain highly lethal. To understand thyroid cancer evolution in both children and adults, we analyzed single-cell transcriptomes of 423,733 cells from 81 samples and spatially resolved key tumor and microenvironment populations across 28 tumors with spatial transcriptomics, including rare and unique composite WDTC/ATC tumors and pediatric diffuse sclerosing thyroid carcinomas. Additionally, we identified gene signatures of stromal cell populations in 5 large thyroid cancer bulk RNA-sequencing cohorts. Through this multi-institutional effort, we defined a population of POSTN+ myofibroblast cancer-associated fibroblasts (myCAFs) that are intimately associated with invasive tumor cells and correlate with poor prognosis, lymph node metastasis, and disease progression in thyroid carcinoma. We also revealed a population of inflammatory CAFs that are distant to tumor cells and are found in the inflammatory stromal microenvironment of autoimmune thyroiditis. Together, our study provides spatial profiling of thyroid cancer evolution in samples with mixed WDTC/ATC histopathology and identifies a prognostic myCAF subtype with potential clinical utility in predicting aggressive disease in both children and adults.

Humans↗
PubMed2025

Transcriptome-wide analysis reveals sequence selection to avoid mRNA aggregation in E. coli.

The stability of RNA base pairing and its limited four-letter code create an intrinsic potential for promiscuous RNA-RNA interactions. In vitro, such interactions drive RNA to self-assemble into aggregates. This raises a fundamental unanswered question: within a confined cellular volume at physiological mRNA abundances, how much aggregation would arise from sequence-encoded chemistry alone? Here, we establish this baseline with large-scale kinetic simulations of the E. coli transcriptome. Our simulations reveal that sequence-encoded base-pairing energetics is sufficient to generate a dynamic network of large aggregates, organized by long, multivalent mRNA hubs. Strikingly, evolutionary analysis shows that native E. coli sequences exhibit clear signatures of selection to counteract this propensity: they fold more stably, minimize unstructured regions, and form weaker intermolecular contacts than dinucleotide-preserving controls. These findings demonstrate that maintaining transcriptome solubility has been a significant, previously unrecognized constraint shaping genome evolution, and provide a new lens to interpret cellular RNA management.

Biological Sciences (Biophysics and Computational ↗
PubMed2026

Spatial transcriptomics of Ciona adult brains reveals functional zonalization and insights into neural gland function.

The ascidian Ciona is a pivotal chordate model for illuminating the evolutionary origins of the vertebrate brain. Here, spatial transcriptomics of the adult Ciona neural complex, combined with image-based computational super-resolution mapping, resolved distinct tissue domains including the cerebral ganglion, neural gland, ciliated funnel, neural gland duct/dorsal strand, and body wall muscle. Within the cerebral ganglion, high-resolution mapping revealed clear molecular zonalization separating the cortex and medulla, alongside regional specialization within the cortex itself. The neural gland exhibited localized enrichment of genes associated with extracellular matrix and cell-cell interactions. These spatial features suggest that the neural gland functions as a homeostatic and signaling interface, reminiscent of primitive vertebrate meninges or choroid plexus. Overall, this spatially defined gene expression map provides a foundational framework for understanding functional regionalization in the tunicate brain and its evolutionary relationship to vertebrate nervous systems.

Ciona↗
PubMed2026

Transposable elements as modulators of homoeologous gene expression in bread wheat: lessons from the pan-transcriptome era.

Bread wheat (Triticum aestivum L.) is an allohexaploid (AABBDD) whose three ancestral subgenomes generate complex patterns of gene regulation. Most genes exist as homoeologous triads, and the relative expression balance among copies, homoeolog expression bias, is central to polyploid evolution and adaptation. Recent high-quality assemblies, long-read transcriptomics, and pan-transcriptome resources have uncovered extensive cultivar-specific transcriptional diversity. Because transposable elements (TEs) compose over 80% of the wheat genome, they are prime candidates for shaping subgenome asymmetry. We synthesize recent pan-genomic and transcriptomic evidence, including genome-wide associations between TE insertions and genome-specific expression, and propose a unifying framework in which TEs modulate homoeolog expression by donating cis-regulatory sequences, altering chromatin states, producing small RNAs, and driving structural variation. We discuss experimental and computational challenges for establishing causality, and outline future functional and translational strategies to leverage TE-associated regulatory diversity in wheat breeding.

Triticum↗
PubMed2026

Evolutionary Reorganization of Transcriptomic Architecture Across a UVB Tolerance Gradient in Fish.

Environmental stressors such as ultraviolet radiation impose strong selective pressures on organisms, yet how adaptation to such stressors shapes transcriptomic responses at the network level remains poorly understood. Although stratospheric ozone is recovering globally, substantial regional variation in UV exposure persists, particularly in high-altitude environments where extreme UV levels can occur. Here, we compared three fish models representing distinct biological responses to UVB exposure: wild-type zebrafish (Danio rerio), a melanin-deficient zebrafish mutant (nacre) lacking a major protective mechanism against UVB damage, and the high-altitude Andean killifish Orestias ascotanensis, a species naturally exposed to extreme UVB radiation. Together, these models define a gradient spanning physiological protection, impaired protection, and evolutionary adaptation to UVB stress. Using RNA-seq and protein-protein interaction networks, we show that transcriptomic responses differ markedly across this gradient. Wild-type and nacre zebrafish exhibited relatively limited transcriptomic changes (∼2%-2.4% of genes changing), whereas O. ascotanensis displayed a large-scale and highly coordinated response (∼21.6% of genes changing) characterized by functionally specialized networks enriched in DNA repair pathways. These differences involved not only transcriptomic magnitude but also marked reorganization of transcriptomic architecture. Integration with positive selection analyses revealed that positively selected genes were concentrated within highly interconnected regions of transcriptomic networks, consistent with adaptation involving network reorganization. Furthermore, ortholog-based analyses suggest that adaptive responses involve differential reorganization of a conserved functional background. Together, our results support a model in which adaptation to environmental stress is associated with the reorganization of conserved transcriptomic networks across physiological and evolutionary contexts, providing a systems-level perspective on the molecular basis of adaptation.

UVB radiation↗
PubMed2025

Advances in tumor subclone formation and mechanisms of growth and invasion.

Tumor subclones refer to distinct cell populations within the same tumor that possess different genetic characteristics. They play a crucial role in understanding tumor heterogeneity, evolution, and therapeutic resistance. The formation of tumor subclones is driven by several key mechanisms, including the inherent genetic instability of tumor cells, which facilitates the accumulation of novel mutations; selective pressures from the tumor microenvironment and therapeutic interventions, which promote the expansion of certain subclones; and epigenetic modifications, such as DNA methylation and histone modifications, which alter gene expression patterns. Major methodologies for studying tumor subclones include single-cell sequencing, liquid biopsy, and spatial transcriptomics, which provide insights into clonal architecture and dynamic evolution. Beyond their direct involvement in tumor growth and invasion, subclones significantly contribute to tumor heterogeneity, immune evasion, and treatment resistance. Thus, an in-depth investigation of tumor subclones not only aids in guiding personalized precision therapy, overcoming drug resistance, and identifying novel therapeutic targets, but also enhances our ability to predict recurrence and metastasis risks while elucidating the mechanisms underlying tumor heterogeneity. The integration of artificial intelligence, big data analytics, and multi-omics technologies is expected to further advance research in tumor subclones, paving the way for novel strategies in cancer diagnosis and treatment. This review aims to provide a comprehensive overview of tumor subclone formation mechanisms, evolutionary models, analytical methods, and clinical implications, offering insights into precision oncology and future translational research.

Humans↗
PubMed2026

Comparative Transcriptomic Analyses Identify Candidate Genes for Convergent Reproductive Shifts in a Bimodal Viviparous Amphibian.

Shifts in reproductive mode represent key evolutionary innovations that shape species' life histories and evolutionary trajectories. Species showing bimodal reproductive strategies with multiple independent origins offer a rare opportunity to gain insights into the adaptive processes and mechanisms underlying convergent traits. The fire salamander, Salamandra salamandra, is the only amphibian exhibiting intraspecific variation in reproductive mode across multiple independent reproductive shifts, enabling investigation of the transition between larviparity (females give birth to aquatic larvae) and pueriparity (females give birth to fully developed terrestrial juveniles) within a single species and across different timescales. Pueriparity is an adaptive innovation that skips the aquatic larval stage, allowing individuals to exploit habitats with no available water bodies. The fire salamander is larviparous across most of its range, but pueriparity has evolved independently at least three times: once in the early Pleistocene within S. s. bernardezi in the mountains of northern Spain, and more recently on two land-bridge islands (NW Spain) inhabited by S. s. gallaica. To identify candidate genes associated with these distinct reproductive modes, we compared gene expression profiles of the uterus and oviduct of pregnant females across two independent evolutionary transitions using RNA-sequencing. We detected shared changes in maternal gene expression among pueriparous S. s. bernardezi and S. s. gallaica relative to their larviparous counterparts, in addition to differences unique to each independent evolutionary transition. Functional enrichment analyses indicated that differentially expressed genes were associated with reproductive timing, angiogenesis, and maternal signalling, consistent with the phenotypic differences observed in the uterine environment and embryonic development between the two reproductive modes. This study represents an important first step towards understanding the genomic basis of the evolution of pueriparity in a remarkable bimodal reproductive system, and provides transcriptomic resources and candidate genes for future research into the genomic architecture underlying this poorly understood adaptive trait.

Animals↗
PubMed2026

Dynamic evolution of chaperone-mediated autophagy is associated with tumor microenvironment remodeling and prognostic stratification in lung adenocarcinoma: insights from single-cell transcriptomics, ensemble machine learning, and experimental validation.

BACKGROUND: Lung adenocarcinoma (LUAD) shows prognostic heterogeneity, and tumor-node-metastasis (TNM) staging is limited for individualized management. Chaperone-mediated autophagy (CMA) maintains proteostasis, but its role during adenocarcinoma in situ (AIS)-minimally invasive adenocarcinoma (MIA)-invasive adenocarcinoma (IAC) progression remains unclear. METHODS: Single-cell RNA sequencing (scRNA-seq) data from GSE189357 and bulk transcriptomes from The Cancer Genome Atlas (TCGA)-LUAD and Gene Expression Omnibus (GEO) cohorts were integrated. CMA activity, cell-cell communication, weighted gene co-expression network analysis (WGCNA), tumor-normal differential expression, machine-learning survival modeling, tumor microenvironment (TME) features, drug sensitivity, and EPC1 function were analyzed. RESULTS: CMA-high tumor epithelial cells increased from AIS (58.1%) to MIA (65.7%) but declined in IAC (44.4%; p < 0.001). CMA-low cells preferentially received fibroblast-derived extracellular matrix cues. A CMA-negatively correlated module identified 69 core genes. Random survival forest (RSF) performed best among 117 machine-learning combinations (mean concordance index > 0.873). High-risk patients had worse survival across cohorts, and the risk score was independently associated with overall survival (hazard ratio = 16.013, 95% confidence interval: 9.579-26.768, p < 0.001). High-risk tumors showed proliferative activation and M0 macrophage enrichment, whereas low-risk tumors showed stronger immune-related signaling. EPC1 overexpression suppressed malignant phenotypes in A549 cells. CONCLUSION: CMA dynamics are associated with stromal and immune remodeling during LUAD progression. A CMA-based model provides robust prognostic stratification and may offer a basis for future TME-guided studies.

Chaperone-mediated autophagy↗
PubMed2025

Deconvoluting clonal and cellular architecture in IDH-mutant acute myeloid leukemia.

Isocitrate dehydrogenase 1/2 (IDH) mutations are early initiating events in acute myeloid leukemia (AML). The complex clonal architecture and cellular heterogeneity in IDH-mutant AML underlies the heterogeneous clinical presentation and outcomes. Integrating single-cell genotyping and transcriptomics, we demonstrate a stem-like and inflammatory phenotype of IDH-mutant AML and identify clone-specific programs associated with NPM1, NRAS, and SRSF2 co-mutations. Furthermore, these clones had distinct responses to treatment with combination IDH inhibitors and chemotherapy, including elimination, reconstitution of myeloid differentiation, or retention within progenitor populations. At relapse after IDH inhibitor monotherapy, we identify upregulated stemness, inflammation, mitochondrial metabolism, and anti-apoptotic factors, as well as downregulated major histocompatibility complex (MHC) class II antigen presentation. At the pre-leukemic stage, we observe upregulation of IDH2-associated pathways, including inflammation. We deliver a detailed phenotyping of IDH-mutant AML and a framework for dissecting contributions of recurrently mutated genes in AML at diagnosis and following therapy, with implications for precision medicine.

Leukemia, Myeloid, Acute↗
PubMed2026

T2T Genome Assembly and Multi-Omics Data Reveal Terrestrial Adaptation and Mucus Biosynthesis in Tropical Leatherleaf Slug (Laevicaulis alte).

Laevichaulis alte is a slug in the order Systellommatophora that evolved from aquatic ancestors and now faces strong challenges from desiccation, respiration on land, and novel pathogens. Its mucus is essential for water retention, locomotion, and defense. To link terrestrial adaptation with mucus biosynthesis, we generated a gap-free genome assembly of L. alte using PacBio HiFi reads, Oxford Nanopore ultra-long reads, and Hi-C data. The genome shows low heterozygosity and holocentromeric chromosomes. Functional metabolomics revealed marked metabolic shifts between L. alte and the closely related aquatic species Peronia verruculata. In L. alte, differential metabolites were enriched in lipid metabolism, immune regulation, and stress response pathways, consistent with life in a dry and microbe-rich terrestrial environment. Comparative genomics and transcriptomics identified candidate genes linked to mucus secretion and physiological adaptation, including VEGF, ASGR2, and COL6A6. Further analyses highlighted the vascular endothelial growth factor (VEGF) gene family as a key regulator connecting angiogenesis, tissue remodeling, and mucus production pathways in L. alte. Together, this gap-free genome and multi-omics dataset establish a molecular framework that links genomic innovation, mucus biology, and terrestrial adaptation in Systellommatophora, and they offer a basis for understanding ecological niche specialization in land molluscs.

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↗
PubMed2025

The chromosome-level genome of Stylosanthes guianensis provides insights into genome evolution and environmental adaptation.

Stylosanthes guianensis is a leguminous forage crop of significant economic importance, primarily distributed in tropical and subtropical regions. It exhibits strong adaptability to various stresses, yet the genetic basis underlying this trait remains unclear. In this study, we constructed the first chromosome-scale reference genome of S. guianensis using a combination of Nanopore and Hi-C sequencing technologies. The assembled genome size is 1254&#x2009;Mb, with 10 pseudochromosomes. Using Nanopore full-length transcriptome data, we generated high-quality transcript-level gene annotations, identifying 36&#x2009;585 gene models and 110&#x2009;601 transcripts. The repetitive sequences in S. guianensis account for 79.16% of the genome, with the extensive expansion of Gypsy elements in long terminal repeats contributing to its genome size enlargement. Comparative genomic and transcriptomic analyses revealed that flavonoid metabolism plays a pivotal role in stress adaptation, providing new insights into the genetic basis of stress tolerance. Additionally, we generated whole-genome methylation profiles under cold treatment and control conditions, offering valuable data for future epigenomic research. These findings provide essential molecular resources for understanding stress resilience in S. guianensis and advancing its molecular breeding.

Genome, Plant↗
PubMed2026

Evolutionary architecture and lineage-specific diversification of Forkhead box transcription factors in Perna viridis.

The Forkhead box (Fox) transcription factors are evolutionarily conserved regulators of development, cell cycle, and apoptosis across metazoans. This study provides the first comprehensive genome-wide analysis of the Fox gene family in the Asian green mussel (Perna viridis). We identified 28 Fox genes distributed across 10 chromosomes. Comparative analysis reveals the absence of the FoxI, FoxQ1, FoxR and FoxS subfamily, consistent with other bivalves and indicative of lineage-specific gene loss during molluscan evolution. Notably, gene duplications in the FoxAB, FoxD, FoxH, FoxN1-4, FoxQ2 and FoxQD subfamilies may reflect functional diversification associated with environmental adaptation. Exon-intron structural variability, including intron loss in several paralogues, suggests structural diversification and potential regulatory variation. Phylogenetic reconstruction confirmed the monophyly of core Fox classes while highlighting divergent expansion patterns in lophotrochozoans. Selection analyses showed strong purifying selection across duplicated Fox paralogs, supporting functional conservation after lineage-specific expansion. Gene Ontology enrichment linked Fox genes to stress response, apoptosis, and transcriptional regulation. By integrating phylogenetic, structural, and transcriptomic analyses, this study provides a genomic framework for understanding Fox gene organisation, evolution, and tissue-associated expression patterns in Perna viridis and establishes a comparative resource for future functional studies in bivalves.

Animals↗
PubMed2026

highSpaClone enables copy number alteration inference and tumor subclone analysis for high-resolution spatial transcriptomics.

High-resolution spatially resolved transcriptomics (SRT) offers unprecedented opportunities to investigate tumor heterogeneity but poses substantial computational and analytical challenges. Here, we present highSpaClone, a computational framework for copy number alteration (CNA) inference and tumor subclone identification from high-resolution SRT data across multiple spatial scales. By integrating spatial constraints into CNA estimation and clonal clustering, highSpaClone enables neighboring spatial locations to share information, thereby improving the robustness of genomic signals and the accuracy of subclone delineation. Across multiple Xenium and Visium HD datasets, highSpaClone revealed unique transcriptional programs, clonal evolutionary trajectories, and distinct tumor-microenvironment interactions. Furthermore, in human colorectal cancer samples, highSpaClone detected CNA events in histologically normal epithelial regions, highlighting early genomic alterations associated with field cancerization. These findings establish highSpaClone as a scalable framework for studying clonal architecture and tumor evolution.

CP: cancer biology↗
PubMed2024

Comparative transcriptomics of Venus flytrap (Dionaea muscipula) across stages of prey capture and digestion.

The Venus flytrap, Dionaea muscipula, is perhaps the world's best-known botanical carnivore. The act of prey capture and digestion along with its rapidly closing, charismatic traps make this species a compelling model for studying the evolution and fundamental biology of carnivorous plants. There is a growing body of research on the genome, transcriptome, and digestome of Dionaea muscipula, but surprisingly limited information on changes in trap transcript abundance over time since feeding. Here we present the results of a comparative transcriptomics project exploring the transcriptomic changes across seven timepoints in a 72-hour time series of prey digestion and three timepoints directly comparing triggered traps with and without prey items. We document a dynamic response to prey capture including changes in abundance of transcripts with Gene Ontology (GO) annotations related to digestion and nutrient uptake. Comparisons of traps with and without prey documented 174 significantly differentially expressed genes at 1 hour after triggering and 151 genes with significantly different abundances at 24 hours. Approximately 50% of annotated protein-coding genes in Venus flytrap genome exhibit change (10041 of 21135) in transcript abundance following prey capture. Whereas peak abundance for most of these genes was observed within 3 hours, an expression cluster of 3009 genes exhibited continuously increasing abundance over the 72-hour sampling period, and transcript for these genes with GO annotation terms including both catabolism and nutrient transport may continue to accumulate beyond 72 hours.

Droseraceae↗