Citations

Polycystic kidney disease (PKD) arises from mutations in cilia-associated genes, such as PKD1 and PKD2 , expressed in renal epithelial cells, leading to progressive kidney dysfunction and end-stage kidney disease (ESKD). PKD patients exhibit significant heterogeneity in disease progression, largely due to genetic and environmental modifiers. Like patients, mouse models of PKD also exhibit significant heterogeneity with regards to the gene mutated, age of disease onset, and rate of disease progression. To elucidate the cellular and molecular consequences of these variables, we constructed an integrated single-cell RNA sequencing atlas across mouse models of PKD, mapping changes in cell type composition, gene expression, and intercellular signaling networks across the whole atlas and within individual models. Across models, single cell RNA sequencing (scRNAseq) data revealed increased Spp1 (osteopontin) expression and signaling from PKD-enriched clusters. Global deletion of Spp1 in Pkd1 RC/RC mice resulted in a modest reduction in cyst severity and improved kidney function. From these studies, we created a freely available, searchable website (https://bmblx.bmi.osumc.edu/scPKD/) that can be used to identify cross- and intra-model changes in gene expression, guiding researchers to new therapeutic targets for treating PKD.

Summary Long-read nanopore sequencing enables simultaneous detection of germline variation and native DNA base modifications on individual DNA molecules, providing a unique opportunity to investigate allele-specific epigenetic regulation. Here, we performed whole-genome nanopore sequencing on normal and tumor prostate tissues to characterize differential methylation, methylation entropy, and allele-specific methylation (ASM) associated with noncoding genetic variants. Genome-wide analysis identified extensive cancer-associated differentially methylated regions (DMRs), with hypermethylated DMRs significantly enriched near transcription start sites and transcriptional regulatory regions. Integration with transcriptomic datasets revealed strong inverse relationships between promoter methylation and gene expression, while 5-hydroxymethylcytosine (5hmC) levels positively correlated with transcriptional activity across gene bodies. Using fragment-level methylation patterns enabled by long-read sequencing, we further quantified methylation entropy incorporating both 5mCG and 5hmCG states. Cancer-hypermethylated DMRs exhibited markedly reduced entropy, consistent with clonal fixation of methylation states during tumor progression. Entropy profiling across chromatin annotations demonstrated maximal epigenetic heterogeneity at partially modified enhancer-associated regions. To investigate cis-regulatory genetic effects, we developed a simple ASM framework (nanoASM) that can partition sequencing reads by allelic state and identifies allele-specific DMRs directly from long-read data. Compared with conventional population-level mQTL analysis, ASM demonstrated substantially improved statistical efficiency by leveraging within-individual contrasts and reducing sample-level heterogeneity. Although germline single nucleotide polymorphisms (SNPs) were largely shared between normal and tumor tissues, ASM patterns differed substantially, with tumor-associated ASM regions displaying significantly larger genomic span and stronger allelic methylation differences. Comparative analysis with TCGA prostate mQTL and GTEx prostate eQTL datasets demonstrated substantial concordance between ASM directionality and downstream transcriptional effects, particularly for variants located within DMRs and near transcription start sites. At the IRX4 prostate cancer risk locus, ASM identified an androgen-responsive regulatory domain overlapping AR ChIP-seq and H3K27ac peaks, nominating rs6885084 as a candidate functional variant. At the PSCA locus, ASM anchored by rs4736369 was associated with allele-specific methylation, chromatin activation, transcript abundance, and isoform usage. Together, these findings establish nanopore-based ASM analysis as a powerful approach for resolving functional noncoding variants and their regulatory domains they control in prostate cancer.

We report that moderate hyperoxia during a critical window of mucociliary differentiation disrupts terminal maturation in neonatal airway epithelial cells cultured in a three-dimensional model. Hyperoxia induces mitochondrial bioenergetic dysfunction and metabolic reprogramming, with proteomic analysis identifying complex I disruption as a key driver of impaired differentiation. Overall, these findings reveal a previously underrecognized link between mitochondrial bioenergetics and airway epithelial development, positioning metabolic dysfunction as an early trigger of hyperoxia-induced neonatal airway injury.

ABSTRACT Premature infants often require supplemental oxygen therapy, however, exposure to supraphysiological oxygen (hyperoxia) can disrupt normal lung development and contribute to bronchopulmonary dysplasia (BPD). Mitochondrial dysfunction is increasingly recognized as a contributor to hyperoxia-induced BPD. However, the effects of hyperoxia on mitochondrial function and mucociliary differentiation in the developing upper airway epithelium remain poorly understood. This study tested the hypothesis that hyperoxia impairs neonatal airway mucociliary differentiation by disrupting mitochondrial bioenergetic function. Neonatal tracheal airway epithelial cells (nTAECs) from term infants (n=5) were cultured in a 3D air-liquid interface (ALI) model and exposed to 60% O₂ during the mid-phase of differentiation (ALI day 7-14). Cellular phenotype was assessed using immunofluorescence staining and gene expression analyses. Mitochondrial function was evaluated through Seahorse metabolic flux analysis, and global protein changes were characterized by quantitative proteomics. Hyperoxia exposure significantly impaired terminal epithelial differentiation, characterized by reductions in ciliated and goblet cells. Seahorse assay revealed a decrease in baseline oxygen consumption and mitochondrial ATP production, accompanied by a compensatory increase in glycolytic ATP production. Quantitative proteomics identified disruption of mitochondrial Complex I as a central feature of the hyperoxic response. Downstream proteomic pathway analyses further confirmed the metabolic shift from mitochondrial to glycolytic ATP production and demonstrated altered epithelial differentiation pathways, including NOTCH and TGF-β signaling. These findings reveal that moderate hyperoxia impairs mitochondrial bioenergetics and alters metabolic programming, leading to disrupted mucociliary differentiation. Future in vivo studies should evaluate mitochondrial oxidative fitness as a therapeutic target in neonatal lung disease. NEW & NOTEWORTHY We report that moderate hyperoxia during a critical window of mucociliary differentiation disrupts terminal maturation in neonatal airway epithelial cells cultured in a 3D model. Hyperoxia induced mitochondrial bioenergetic dysfunction and metabolic reprogramming, with proteomic analysis identifying Complex I disruption as a key driver of impaired differentiation. Overall, these findings reveal a previously underrecognized link between mitochondrial bioenergetics and airway epithelial development, positioning metabolic dysfunction as an early trigger of hyperoxia-induced neonatal airway injury.

Abstract Background The mechanisms underlying metastatic latency in renal cell carcinoma (RCC) remain poorly understood. Methods This study evaluated 2 large independent cohorts for differences in tumor biology between patients who developed metastases early (≤1 year after nephrectomy) and those with late onset (>3 years). Results In the discovery cohort (n = 161), late metastatic RCC was associated with clear cell histology (88.9% vs 78.7%), lower pathological stage (pT1-2; 40.3% vs 18.0%), and favorable histopathological features including low grade (40.0% vs 2.3%), less sarcomatoid (5.6% vs 21.8%), and reduced necrosis (37.7% vs 78.3%; all P < .02). Late metastatic RCC tumors exhibited increased angiogenesis (63.5% vs 19.4%) and reduced inflammation (78.8% vs 50.0%; all P < .02) profiles. Genomic driver analyses revealed comparable rates of PBRM1 and SETD2 loss in late and early metastatic RCC, while BAP1 loss was significantly less common in late metastatic RCC (7.5% vs 27.1%; P < .02). In multivariable models, BAP1/PBRM1/SETD2 status and tumor necrosis emerged as key discriminators of late metastatic RCCs. These findings were confirmed in the second cohort (n = 307). Late metastatic RCC was enriched for fatty acid oxidation and angiogenesis pathways, supporting a less aggressive phenotype. This was further evidenced by a lower engraftment rate in murine models (0% vs 36.5%; P < .001) and significantly longer overall survival from the time of metastasis (median survival doubled, P < .001). Interestingly, late metastatic RCC shared genomic and phenotypic features with RCC that metastasizes to the pancreas, suggesting a common underlying biology influencing both metastatic latency and pancreatic tropism. Conclusions Overall, these findings advocate for recognition of late metastatic RCC because of its distinct biology and improved prognosis.