Search PubMedSearch

SEARCH · Search PubMed

Results for “Gene Expression Profiling”

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.

862 records · Page 16Linked to original sources

Dynamics and virulence of Enterobacteriaceae reservoirs harboring blaCTX-M group 1 in community wastewater.

UNLABELLED: Extended-spectrum beta-lactamase (ESBL)-producing bacteria are ubiquitous and can cause serious infections. Here, we examined untreated community wastewater influent as a reservoir for blaCTX-M group 1 organisms and their virulence potential. Raw influent samples (n = 268) were collected from four wastewater treatment plants (WWTPs) representing dense urban populations. We observed that blaCTX-M group 1 levels were high at all WWTPs and only ~1-2 log10 lower and not correlated to common human-specific microbiome fecal markers, Lachno3 and HF183, indicating a lack of connection to human fecal inputs. Concentrations of blaCTX-M group 1 genes and markers for presumptive host organisms Escherichia coli and Klebsiella pneumoniae were influenced by travel time and season. Amplicon sequencing revealed high diversity of blaCTX-M group 1-9 genes, with 63% belonging to group 1. Selective culture and 16S rRNA gene sequencing showed blaCTX-M group 1 isolates were 26% E. coli, 26% K. pneumoniae, 40% other Enterobacteriaceae, and 8% Aeromonas. Overall, E. coli averaged 3.6E7 cells/L, with 3% of all E. coli found to contain blaCTX-M group 1. Whole-genome sequencing of blaCTX-M group 1 E. coli from wastewater revealed resistance and virulence gene profiles similar to clinical isolates and distinct from other wastewater ESBL-resistant and non-resistant E. coli. Interpretation of wastewater data needs to consider both the existence of environmental reservoirs that contain potentially pathogenic organisms and the strong influence the dynamics of the conveyance system can have on final concentrations measured at the WWTP. IMPORTANCE: The CTX-M enzyme family is highly abundant in nosocomial, community, and environmental settings and is leading to treatment of infections with carbapenem antibiotics, a last-line therapeutic option. The progressive increase of the clinically relevant blaCTX-M group 1 resistance genes in the human population warrants investigation, particularly to understand the establishment and dynamics of environmental reservoirs. This study utilized molecular and culture methods to gain insight into the possible origin, abundance, and dynamics of blaCTX-M group 1 genes in untreated wastewater influent samples. We found extremely high levels of these genes, with Escherichia coli as a major host organism that closely resembled clinical strains, suggesting they are seeded and propagate in sewer pipe systems. The significance of our research is in developing approaches to monitor antimicrobial resistance reservoirs in community wastewater, which could shed light on global burdens and potential transmission cycles and indicate increasing inputs of clinically relevant strains originating from human populations.

E. coli

Genetic overlap between estimated glomerular filtration rate and cardiovascular disease identifies potential targets for cardiorenal syndrome.

Heart and kidney diseases frequently coexist, but the genetic basis of this relationship remains unclear. We analyzed genetic data from large-scale studies to investigate how kidney function (estimated glomerular filtration rate, eGFR) and six common cardiovascular diseases share genetic risk factors. Using MiXeR method, and conjunctional false discovery rate (conjFDR) to identify overlapping genetic regions, we found 478 shared genomic loci between eGFR and cardiovascular diseases. These shared genes are involved in tissue development and structure. We also identified 29 genes that could be targeted by existing medications approved by the US Food and Drug Administration, such as PRKAG2, PDE1A, and IGF1R. Among these, genetically predicted higher level of IGF1R expression is associated with a higher eGFR, which reflects good kidney function and is protective against cardiorenal diseases, such as atrial fibrillation, and myocardial infarction. These findings reveal genetic overlap between kidney function and cardiovascular diseases, highlighting potential targets for understanding and treating cardiorenal syndrome.

Humans

A complete hlyCABD-like RTX operon marks a virulence-associated subset of trh-positive Vibrio parahaemolyticus from Hangzhou Bay, China.

Vibrio parahaemolyticus remains a major cause of seafood-associated gastroenteritis, yet routine surveillance still relies largely on the canonical hemolysin markers thermostable direct hemolysin (tdh) and tdh-related hemolysin (trh). To determine whether this framework overlooks accessory virulence determinants in trh-positive lineages, we analyzed 193 V. parahaemolyticus isolates collected between 2022 and 2025 from clinical, environmental, and seafood-associated sources in the Hangzhou Bay region of China. Serotyping identified 45 serotypes, with O10:K4 predominating among clinical isolates. Both clinical and non-clinical populations showed open pan-genomes, although the non-clinical group carried a larger accessory gene pool. We identified a complete hlyCABD-like RTX operon in 10 trh-positive isolates with T3SS2-associated virulence backgrounds. These RTX-positive isolates were distributed across seven sequence types and three of five phylogenetic groups. This distribution was lineage-restricted but non-clonal. In the representative hybrid-assembled genome, the operon occurred within a mosaic genomic region containing additional virulence- and mobility-associated genes, indicating a composite pathogenicity island-like element. In the tested subset, RTX-positive isolates showed significantly greater hemolytic activity than RTX-negative trh-positive isolates. This significant difference was consistently observed in both plate-based and liquid assays, and within the RTX-positive subset, hlyA expression correlated with hemolytic activity, whereas the trh gene and the tlh (thermolabile hemolysin) gene did not. A complete hlyCABD-like RTX operon therefore identifies a hemolysis-associated subset of trh-positive V. parahaemolyticus and supports its further evaluation as an additional target for food safety surveillance.

Vibrio parahaemolyticus

Integrated assessment of biocontrol potential and genome analysis of endophytic Bacillus velezensis MGL-B1 against mango stem-end rot.

Mango stem-end rot is a globally significant postharvest disease that severely threatens the mango industry, primarily caused by Botryosphaeria dothidea. However, information on biocontrol agents targeting this pathogen in mango remains limited. In this study, we isolated and identified a strain of Bacillus velezensis MGL-B1 from mango leaf tissues for the first time, which exhibited broad-spectrum antifungal activity. Both in vitro and in vivo assays demonstrated that MGL-B1 effectively inhibited the growth of B. dothidea, with an in vivo biocontrol efficacy reaching 83.72 ± 5.10%, comparable to that of the commonly used chemical fungicide thiabendazole. Further mechanistic analysis revealed that MGL-B1 acts by directly disrupting the integrity of the pathogen's mycelial cell membrane. In addition, its released volatile organic compounds (VOCs) also displayed significant antifungal activity, with components such as 2-nonanone, 2-nonanol, and phenylethyl alcohol being confirmed to exert antifungal effects in in vitro fumigation assays. qPCR analysis showed that MGL-B1 treatment significantly upregulated the transcriptional levels of genes involved in plant-pathogen interaction, phenylpropanoid biosynthesis, and antioxidant defense pathways in mango fruits, with upregulation folds of 16.32, 37.19, and 75.93, respectively; meanwhile, the expression of browning-related genes such as polyphenol oxidase (PPO) was markedly suppressed. Whole-genome sequencing further revealed 14 biosynthetic gene clusters for antimicrobial compounds, including five unknown gene clusters. Collectively, B. velezensis MGL-B1 represents a promising biocandidate strain with multiple antifungal mechanisms and excellent control efficacy, providing a valuable resource for green and sustainable management of mango diseases.

Mangifera

In vivo genome-wide CRISPR screens identify FOXR1 as a suppressor of CD8+ T cell antitumor immunity.

T cell dysfunction critically limits the efficacy of T cell-based immunotherapies in solid tumors, yet the intrinsic regulators of T cell dysfunction remain incompletely understood. Through an in vivo genome-wide CRISPR screen in tumor-infiltrating CD8+ T cells, we identified Forkhead Box R1 (FOXR1) as a potent transcriptional suppressor of CD8+ T cell effector functions. Genetic ablation of FOXR1 significantly enhanced cytokine production and cytotoxic capacity in both murine and human CD8+ T cells, whereas its overexpression impaired T cell activation and effector molecule expression. Mechanistically, multiomics integration of RNA-seq, CUT&Tag-seq, and ATAC-seq revealed that FOXR1 binds directly to promoter regions of key effector genes, including IL2, GZMB, and PRF1, and represses their expression. Importantly, FOXR1 deletion in human anti-CD19 CAR T cells improved their efficacy against solid tumors, demonstrating that FOXR1 is a checkpoint of T cell effector function and targeting FOXR1 is a promising strategy to enhance CAR T cell efficacy against solid tumors.

Animals

Discovery of NAT-6-321056 as a novel modulator of VEGFR2 signaling to suppress tumor angiogenesis.

Vascular endothelial growth factor receptor 2 (VEGFR2) is a master regulator of angiogenesis and cancer progression. However, current VEGFR2 modulators face significant challenges, including off-target toxicity and acquired resistance, underscoring the urgent need for novel therapeutic agents with improved efficacy and safety profiles. Here, we reported that virtual screening of 39,442 natural products from the ZINC natural products-derived library, coupled with molecular docking and molecular dynamics (MD) simulations to evaluate the binding stability of candidate compounds, identified NAT-6-321056 as a highly promising modulator of VEGFR2 signaling. Biological evaluations demonstrated that NAT-6-321056 exerted potent inhibition on the growth of a broad spectrum of cancer cells, including both solid tumors and hematological malignancies. In EA.hy 926 endothelial cells and SK-N-DZ neuroblast cells, the compound significantly suppressed proliferation, migration, and invasion. Microscale thermophoresis (MST) confirmed direct binding of NAT-6-321056 to VEGFR2 with favorable affinity. Kinase profiling against a panel of 33 kinases indicated that NAT-6-321056 exhibited a multi-kinase modulation profile. Mechanistic studies revealed that NAT-6-321056 suppressed the expression of hypoxia-inducible factor 1-alpha (HIF-1α) and was associated with reduced VEGFR2 phosphorylation and attenuation of the downstream ERK/JNK/AKT signaling pathways. Moreover, NAT-6-321056 exhibited robust in vivo anti-angiogenic effects in both the chick chorioallantoic membrane (CAM) assay and transgenic zebrafish vascular fluorescence imaging models. Computational absorption, distribution, metabolism, excretion, and toxicity (ADMET) prediction suggested acceptable drug-like properties. Collectively, these findings demonstrated that NAT-6-321056 is a promising modulator of VEGFR2 signaling with potent anti-angiogenic activity and represents a viable candidate for cancer therapy.

Vascular Endothelial Growth Factor Receptor-2

Circular RNAs in amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder characterized by the progressive loss of motor neurons, with most cases lacking a clear genetic basis. Emerging evidence highlights the involvement of non-coding RNAs, particularly circular RNAs (circRNAs), in disease onset and progression. Here, we investigated circRNAs implicated in ALS and related motor neuron diseases (MNDs). Here, we provide a general overview of circular RNA metabolism and cellular functions. We then present our systematic literature review that identified ALS-associated circRNAs, followed by in silico analyses of 15 circular RNA candidates that were selected based on the most compelling data regarding ALS. Our results revealed that several circular RNAs regulate ALS-related genes, such as unfolded protein response, oxidative stress, cell cycle regulation, and apoptosis. Protein-RNA interaction analysis further showed that ALS-related circRNAs can sponge 20 RNA-binding proteins. Additionally, molecular docking analysis demonstrated that ALS-associated FUS variants significantly alter its binding affinity to circular RNAs. RNA-seq data from ALS patients confirmed significant alterations in the expression of host genes of ALS-related circRNAs and hub proteins in ALS-affected CNS tissues. Collectively, our findings identify circRNAs as potential key contributors to ALS pathogenesis.

Amyotrophic Lateral Sclerosis

Tracking GAD-specific T-cell expansions in Type 1 diabetes by intradermal GAD-Alum challenge.

Identifying and monitoring autoreactive T cells that drive beta cell destruction remains a major obstacle to developing effective immunotherapies for type 1 diabetes (T1D). These cells are extremely rare in peripheral blood and cannot be accessed directly from the pancreas. We used intradermal injection of Glutamic Acid Decarboxylase (GAD)-Alum to recruit GAD-specific T cells to accessible sites in the skin and skin-draining lymph nodes (LNs), sampled by skin suction blisters and ultrasound-guided LN aspiration. Peripheral blood samples obtained before GAD injection were restimulated with GAD in vitro to detect reactive CD4+ T cells. Single-cell RNA sequencing (scRNAseq) followed by re-expression of selected T cell receptors (TCRs) confirmed antigen specificity. Up to 70% of T cells at the skin injection site were clonally-expanded and 4 of 14 (28%) re-expressed TCRs were GAD-reactive. In LNs 1 of 14 (4%) clonally-expanded TCRs was GAD-reactive, representing ~0.08% of all T-cells. GAD-reactive cells across compartments displayed Th1 and Th17-associated transcription signatures. These results demonstrate the intradermal autoantigen challenge and scRNAseq, enable direct identification and molecular profiling of autoreactive T cells in vivo. This minimally invasive approach provides a powerful platform for tracking antigen-specific T cells to monitor disease activity and evaluate immune interventions in T1D.

Autoimmunity

Early proteomic and metabolic signatures of liver and eye in OAT-deficient mice.

Ornithine aminotransferase (OAT) deficiency causes hyperornithinemia and gyrate atrophy (GA) of the choroid and retina, a rare inherited retinal degeneration. To understand the early molecular changes that make the eye susceptible to damage, we performed quantitative proteomic and metabolomic profiling of liver, retina, and retinal pigment epithelium and choroid (RPE/Cho) from OAT-deficient (Oatrhg) mice prior to detectable vision impairment. In addition to reduced OAT expression and elevated ornithine, methylation-related metabolites such as N(6)-methyl-lysine were altered in all examined tissues of Oatrhg mice. In the liver, excess ornithine was directed into urea cycle metabolism, together with altered expression of detoxification enzymes and histone H2B proteins. In contrast, the retina showed minimal proteomic changes but pronounced alterations in amino acid pathways that support glutamate homeostasis. The RPE/Cho demonstrated the most extensive proteomic changes, particularly in mitochondrial metabolism, cytoskeleton, and extracellular matrix, along with changes in metabolites involved in lysine metabolism, energy metabolism, and antioxidant capacity. Incubation with 13C lysine demonstrated that lysine was primarily degraded in RPE/Cho but not the retina, and ornithine enhanced lysine degradation in an OAT-dependent manner. Together, these findings highlight common and tissue-specific impacts of OAT on the liver and ocular tissues and provide insight into early molecular changes that contribute to the selective vulnerability of the eye in GA. Proteomics data are available via ProteomeXchange (PXD063614) and metabolomics data via MassIVE repository (MSV000101103).

Animals

Multiomic insights into fungal polylactic acid degradation: Metabolic adaptation and hydrolytic mechanisms of Sporobolomyces pararoseus.

Polylactic acid (PLA), a biodegradable polyester from renewable resources, is a sustainable alternative to petrochemical plastics. However, its environmental degradation is inefficient naturally, requiring specific microbial activities. While bacterial PLA-degrading mechanisms are well documented, fungal degrading systems-particularly their molecular mechanisms-are underexplored.We isolated Sporobolomyces pararoseus ZRQ01 from the gut microbiota of PLA-fed mealworms. This fungal strain noticeably degraded PLA in PLA-containing medium supplemented with 2% glucose. Biodegradation assays revealed 22.8% loss of the PLA film weight after 35 days of incubation, and scanning electron microscopy confirmed extensive surface erosion and pore formation. Integrated transcriptomic and proteomic analyses, together with the reference genome of S. pararoseus ZRQ01, revealed that S. pararoseus ZRQ01 upregulates hydrolytic enzymes at both transcript and protein levels to cleave PLA into lactic acid. After lactic acid is transferred into S. pararoseus ZRQ01 cells by monocarboxylate transporters with increased abundance, it is assimilated by pathways of pyruvate metabolism and the TCA cycle with increased protein abundance. Intriguingly, upregulation of genes in autophagy-related and MAPK signaling pathways underscores an adaptive stress response potentially supporting cellular homeostasis and degradation-related gene expression. Our results highlight S. pararoseus ZRQ01's metabolic potential for bioremediation and offer insights into fungal bioplastic degradation pathways.

Polyesters

Age-related differences in semen quality in Holstein-Friesian bulls: a paired within-bull comparison of early and mature reproductive stages.

Genomic selection has changed dairy cattle breeding by increasing the use of young bulls for artificial insemination and shortening the reproductive lifespan of sires. Under these conditions, semen quality at the beginning of commercial use has become an important practical issue. Semen samples from 39 fertile Holstein-Friesian bulls used for commercial AI were collected between 2013 and 2016, during the introduction of genomic selection in Poland. This paired within-bull study compared semen collected from the same bulls at an early reproductive stage (13-20 months; young bulls, YB) and at full maturity (5-6 years; mature bulls, MB). The evaluation included conventional ejaculate traits, CASA-derived motility and kinematic descriptors, mtDNA copy number, and mitochondrial content per sperm cell. Importantly, all ejaculates met the quality requirements for commercial insemination. Ejaculate volume, sperm concentration, mitochondrial DNA copy number, and mitochondrial content did not differ significantly between age groups. The CASA-derived sperm movement profile, in contrast, differed with age. Semen from young bulls showed a higher proportion of progressively motile spermatozoa, whereas semen from mature bulls showed higher velocity-related parameters, including VSL, VCL, and STR. These findings indicate that bull age mainly affected sperm movement characteristics rather than semen output or mitochondrial content. Overall, the results support the use of young bulls in artificial insemination programs and show that age-related differences in semen quality are expressed mainly through changes in the post-thaw sperm motility and kinematic profile.

Animals

Inducible flocculation in Komagataella phaffii enables enhanced biomass separation for biopharmaceutical production.

Biomass separation represents a critical bottleneck in Komagataella phaffii-based biopharmaceutical processes, as typically high cell densities of 40 - 50 % create significant operational, technical and economic challenges for harvest operations. Yeast cell aggregation (flocculation) provides a solution to accelerate cell sedimentation by increasing particle size, thus allowing to improve biomass-supernatant separation efficiency during both natural gravity settling and (continuous) centrifugation operations. This study demonstrates successful engineering of K. phaffii strains with an inducible flocculation phenotype using CRISPR/Cas9-based genome editing to integrate the Saccharomyces cerevisiae FLO1 (ScFLO1) gene under control of various regulatory elements, including methanol-inducible and derepressible promoters. Flocculation strength could be enhanced by implementing transcriptional positive feedback circuits based on the methanol-inducible AOX1 promoter. To address methanol-free production requirements, we developed alternative systems to retrofit PAOX1-based ScFLO1 expression and exploited the derepressible PDF promoter, offering broader compatibility with biopharmaceutical manufacturing facilities. Flocculating cells cultivated in a bioreactor demonstrated significantly improved sedimentation behavior, with considerably lower supernatant turbidity after short low-speed centrifugation or gravity sedimentation compared to non-flocculating controls. Crucially, cell flocculation had no negative impact on product amount and quality when expressing a multivalent NANOBODY® VHH molecule with pharmaceutical relevance. Thus, this work establishes the first genetically engineered flocculation system in K. phaffii compatible with recombinant protein production, providing the basis for an innovative approach to streamline harvest operations in biopharmaceutical processes.

Flocculation

An allograft inflammatory factor enhances sperm viability by modulating intracellular calcium in oyster Crassostrea gigas.

As an important aquaculture bivalve, the Pacific oyster Crassostrea gigas faces severe constraints in artificial reproduction, where low sperm motility often leads to fertilization failure and limits the sustainable development of the oyster aquaculture industry. In the present study, the variation of sperm from different oyster individuals was observed, and high-quality sperm possessed intact, elongated flagella with no structural abnormalities, while low-quality sperm showed shortened flagella with frequent tangling or coiling defects. Transcriptomic analysis comparing high- and low-quality sperm revealed significantly reduced expression of genes associated with sperm motility and release (CgAIF1, CgAchR, CgSEX), sperm quality and development (CgEP4, CgIFi2b), and cryoprotection (CgISPs) in low-quality sperm. Notably, an allograft inflammatory factor (designed as CgAIF1) encoding EF-hand domain, known as Ca2+ binding activity, was among the most significantly downregulated in low-motility sperm. CgAIF1 is highly expressed in haemocytes, ganglia, and gonads of oysters. Incubation with the recombinant AIF1 protein (rCgAIF1) significantly improved sperm curvilinear velocity, thereby enhancing the overall motility of C. gigas sperm. Furthermore, rCgAIF1 incubation increased intracellular Ca2+ levels (2.13-fold at 30 min, 2.71-fold at 60 min) and superoxide dismutase (SOD) activity (1.44-fold at 30 min, 1.24-fold at 60 min) in sperm, suggesting potential roles in calcium homeostasis regulation and antioxidant defense. In conclusion, this study demonstrates that CgAIF1 significantly enhances motility of oyster sperm, providing a scientific basis for artificial breeding and seed production in oyster aquaculture.

Animals

Genome-wide association study of estimated glomerular filtration rate using repeated measurements in the Taiwan Biobank.

BACKGROUND: Chronic kidney disease (CKD) is a major global public health issue, with genetic factors playing a significant role in kidney function. Although genome-wide association studies (GWAS) have identified numerous loci associated with estimated glomerular filtration rate (eGFR), most studies relied on a single time-point measurement, which limits the capacity to account for within-individual measurement variability. METHODS: We performed a repeated-measurement GWAS in the prospective Taiwan Biobank (Taiwanese ancestry; n = 25,004) using two repeated creatinine-based eGFR measurements. Repeated eGFR values were analyzed using a linear mixed-effects model with a subject-specific random intercept and time-varying covariates, providing a more precise estimate of eGFR level. Identified loci underwent functional annotation (expression quantitative trait locus, deleteriousness prediction, and epigenetic markers) and were compared with results from a single-measurement GWAS. RESULTS: Six loci associated with eGFR were identified, including four previously reported regions (1q22, 4q21.1, 11p14.1, and 17q21.2) and two additional loci (6p21.32 and 15q24.2). Functional annotation implicated several candidate genes-such as MUC1/EFNA1, SHROOM3, HLA-DQB1, MPPED2, NRG4, and PGAP3/FBXL20-in the regulation of kidney function. CONCLUSION: Incorporating repeated eGFR measurements into GWAS may improve phenotypic precision for identifying genetic associations with kidney function. This study identified eGFR-associated loci and biologically plausible candidate genes in a Taiwanese population, which require further replication and functional validation.

Chronic kidney disease

Lactylation-related immune-metabolic dysregulation defines prognostic and therapeutic stratification in lung adenocarcinoma.

BACKGROUND: Lactylation links lactate metabolism with inflammatory signaling and immune regulation in tumors. However, its cellular distribution and translational value in lung adenocarcinoma (LUAD) remain unclear. METHODS: Single-cell RNA-sequencing datasets GSE189357 and GSE171145 were integrated to characterize lactylation-related activity, intercellular communication, and malignant epithelial cell states in LUAD. Single-cell-derived lactylation-related differentially expressed genes were mapped to TCGA-LUAD and multiple GEO cohorts. Univariate Cox regression and machine learning algorithms were used to construct a lactylation-related prognostic signature (LRPS). The associations of LRPS with prognosis, immunotherapy response, drug sensitivity, genomic alterations, immune infiltration, and inflammation- and metabolism-related pathways were evaluated. KRT7 was further validated using virtual knockout analysis, spatial transcriptomics, and in vitro and in vivo experiments. RESULTS: lactylation-related transcriptional activity showed heterogeneous distribution across LUAD cell populations and was associated with altered cell-cell communication. In malignant epithelial cells, LRTS-high and LRTS-low states exhibited distinct metabolic, inflammatory, and tumor-related pathway activities. LRPS showed stable prognostic performance in TCGA-LUAD and multiple GEO cohorts and remained an independent prognostic factor. Low LRPS was associated with greater potential benefit from immunotherapy, whereas different LRPS groups displayed distinct drug sensitivity, genomic alteration, and immune microenvironment patterns. KRT7 was highly expressed in LUAD and associated with poor prognosis. KRT7 knockdown suppressed LUAD cell proliferation, migration, invasion, colony formation, and tumor growth in vivo. CONCLUSIONS: This study identifies lactylation-related immune-metabolic dysregulation as a clinically relevant feature of LUAD and develops a single-cell-guided LRPS for prognosis and therapeutic stratification. KRT7 emerged as an LRPS-related functional candidate with experimentally supported roles in malignant LUAD phenotypes.

Immunotherapy

Penalized Cumulative Probability Model for a Continuous Outcome Subject to Detection Limits.

Mixed-type outcome data occur when the outcome variable's distribution is a mixture of both continuous and discrete ordinal variables. Such mixed-type outcomes are common in biomedical, psychological, and the health sciences, particularly for variables having either a detection or quantitation limit. When interest lies in identifying a combination of genomic features associated with a mixed-type outcome, any method used would require a variable selection strategy for high-dimensional data. Unfortunately, few variable selection methods exist for modeling a mixed-type outcome when the covariate space is high dimensional. This study develops a high-dimensional penalized cumulative probability model (CPM), to allow for the identification of genomic features associated with mixed-type outcome of interest. We demonstrated how such model may be estimated using the iterative penalization procedure-the generalized monotone incremental forward stagewise (GMIFS) algorithm. The Model-X knockoffs procedure was combined with the estimation algorithm to control the false discovery rates (FDR) when performing variable selection. Through extensive simulation studies, our penalized CPM was shown to outperform alternative methods in terms of controlled variable selection performance by achieving high statistical power with the FDR being controlled at the target level. We demonstrate the utility of our method by applying it to predict estimated glomeruli filtration rate (eGFR) in kidney transplant recipients at 24 months post-transplant using baseline gene expression data as predictors. Our CPM model identified five genes associated with this mixed-type outcome which have important links to renal disease, which may provide prognostic guidance for kidney transplantation recipients.

Models, Statistical

Age- and sex-adjusted genomic differences between Korean and Beat AML cohorts.

Genomic profiling plays a central role in risk stratification and therapeutic decision-making in acute myeloid leukemia (AML), yet the clinical implications of population-specific genomic architectures remain incompletely defined. We conducted a prospective, multicenter study of 603 adults with newly diagnosed AML in Korea, integrating targeted sequencing of 83 recurrently mutated genes with comprehensive clinical annotation across treatment intensities, including allogeneic hematopoietic stem cell transplantation (allo-HSCT). For contextual comparison, genomic profiles were evaluated against the Beat AML cohort. The overall genomic landscape was broadly conserved, supporting shared core disease biology across populations. However, RUNX1::RUNX1T1, CEBPA, GATA2, KIT, and DDX41 mutations were more frequent in the Korean cohort, whereas FLT3 and NPM1 mutations were less common. These differences translated into a distinct distribution of European LeukemiaNet (ELN) 2022 risk categories, with implications for therapeutic stratification. Notably, most DDX41 alterations were germline (3.2%), highlighting the need for systematic germline evaluation with implications for genetic counseling and donor selection. Although unadjusted overall survival appeared longer in the Korean cohort, this difference was not significant after adjustment for key clinical variables. These findings indicate that population-specific genomic distributions reshape the clinical application of risk stratification and support population-aware precision medicine strategies in AML.

Journal Article

The Role of Small Segmental Duplications in Generating Identical Isoforms Through Alternative Splicing Sites.

Alternative splicing plays a crucial role in expanding proteomic diversity but can also generate identical isoforms under certain conditions. While mutually exclusive splicing of tandem exons has occasionally been reported to produce identical isoforms, the extent to which other splicing events contribute to this phenomenon remains unclear. In this study, we demonstrate that alternative 5' and 3' splice site selection can also lead to the formation of identical isoforms, providing an additional type of splicing event for functional redundancy in transcriptomes. To address this, we analyzed reference genome annotations from 15 plant species, including Arabidopsis thaliana and wheat (Triticum aestivum), obtained from the RefSeq database. Identical isoforms were computationally defined as transcripts with distinct exon-intron structures but identical coding sequences. Our analysis reveals that the majority of alternative 5' and 3' fragments originate from small segmental duplications, suggesting that sequence repetition within gene regions facilitates the emergence of such splicing patterns. We also observed differences in the annotated 5' UTRs of some identical isoforms. However, since the alternative splicing sites themselves were not located within UTRs, these differences may reflect annotation uncertainty rather than genuine AS-derived variation. Given that UTR predictions in reference databases are not always precise, such observations should be interpreted cautiously. Expression analysis using an isoform-specific k-mer approach confirmed that identical isoforms can be differentially regulated. These findings suggest that, beyond expanding protein diversity, alternative splicing can also generate redundant isoforms that are differentially expressed at the RNA level, indicating potential regulatory roles. By elucidating the structural and regulatory factors contributing to the formation and retention of identical isoforms, our study provides new insights into the evolutionary and functional significance of alternative splicing in plants.

Alternative Splicing