Publications

2026

Freeman, Samantha R, Joshua P Hayden, Kailey Hooper, Travis Sullivan, Carmen Perrino, Kimberly Rieger-Christ, and Alex J Vanni. (2026) 2026. “Differential Genetic Expression in the Tunica Albuginea Tissue of Patients With Peyronie’s Disease.”. The Journal of Sexual Medicine 23 (1). https://doi.org/10.1093/jsxmed/qdaf352.

BACKGROUND: Peyronie's disease is characterized by penile curvature, sexual dysfunction, and pain with sexual activity, and its pathogenesis is thought to involve several molecular mechanisms.

AIM: The aim of the present study was to identify differential gene expression in the tunica albuginea (TA) tissue and corpora cavernosa (CC) tissue of a group of patients with Peyronie's Disease relative to a control group of patients with erectile dysfunction.

METHODS: Tissue from the TA and CC was collected from patients at the time of surgery. RNA sequencing analysis was performed on total RNA isolated from penile tissue. Gene ontology analyses were then performed to identify which molecular and biologic processes were associated with the differentially expressed genes. Clinical demographic factors were also collected from chart review of patients.

OUTCOMES: Analysis revealed a distinct gene expression profile in the TA of Peyronie's patients.

RESULTS: A total of 58 patients were included in the study; 57 had CC for RNA isolation and 30 had TA for RNA isolation. Among those with tunica, 15 had Peyronie's and 15 did not. This analysis identified 22 significantly differentially expressed genes between the groups. Bioinformatics analyses revealed that differentially expressed genes were implicated in immune and apoptotic processes. For the corpora population only two genes were found to be differentially expressed.

CLINICAL IMPLICATIONS: Further understanding of the genes implicated in Peyronie's could lead to additional treatments targeted at the prevention and progression of this disease.

STRENGTHS & LIMITATIONS: This study uses next-generation RNA sequencing to analyze whole-gene expression in a relatively large cohort, identifying novel genes involved in Peyronie's disease pathogenesis. Limitations include the use of formalin-fixed tissue and potential confounding due to differences in erectile function between Peyronie's and control groups.

CONCLUSION: Distinct genes were differentially expressed in the TA tissue of patients with Peyronie's disease when compared with controls and were largely implicated in immune and apoptotic processes.

2025

Matheny, Michael E, Amy M Perkins, Kimberly Rieger-Christ, Joseph S Ross, Jialin Mao, Art Sedrakyan, Nilay D Shah, et al. (2025) 2025. “Capacity Assessment for EHR-Based Medical Device Post-Market Surveillance for Synthetic Mid-Urethral Slings Among Women With Stress Urinary Incontinence: A NEST Consortium Study.”. BMJ Surgery, Interventions, & Health Technologies 7 (1): e000193. https://doi.org/10.1136/bmjsit-2023-000193.

OBJECTIVES: To evaluate the feasibility for use of electronic health record (EHR) data in conducting adverse event surveillance among women who received mid-urethral slings (MUS) to treat stress urinary incontinence (SUI) in five health systems.

DESIGN: Retrospective observational study using EHR data from 2010 through 2021. Women with a history of MUS were identified using common data models; a common analytic code was executed at each site. A manual chart review was conducted in a per-site random patient subset to establish a reference standard. Automated text processing (Text Processed Integrated (TPI)) was developed and evaluated at each site to determine the surgical approach and synthetic mesh implantation. Patients were characterized and surgical outcomes were ascertained over 730 subsequent days.

SETTING: Five large tertiary care academic medical centers.

PARTICIPANTS: Across five health systems, 9,906 eligible patients (mean age 57-60 per site) were identified.

MAIN OUTCOME MEASURES: Determination of surgical approach, synthetic mesh implantation, and assessment of the duration of surveillance for mortality and reoperation rates following MUS implantation.

RESULTS: In the TPI cohort analysis, 3,331 patients were identified. Surgical approach per site was retropubic (42% to 77%), transobturator (6% to 44%), single incision (0% to 24%), and adjustable sling (0% to <4%). Concordance rates for TPI using chart review were 71%-90% at each site for the surgical approach and 28%-85% for synthetic mesh implantation. Patient follow-up observation rates for mortality and reoperation ranged from 22% to 36% at 90 days, 15% to 30% at 365 days, and 8% to 19% at 730 days.

CONCLUSION: Using EHR data alone, identification of medical devices and surgical approaches was feasible among women with MUS surgery for SUI, but long-term follow-up ascertainment rates were low. Medical device surveillance using EHR data should be evaluated in the context of the clinical use case, as applicability may vary.

Burks, Eric J, Travis B Sullivan, and Kimberly M Rieger-Christ. (2025) 2025. “Pulmonary Adenocarcinoma of Low Malignant Potential Defines Indolent NSCLC Associated With Overdiagnosis in the National Lung Screening Trial.”. Cancer Biomarkers : Section A of Disease Markers 42 (1): CBM230452. https://doi.org/10.3233/CBM-230452.

BackgroundThe national lung screening trial (NLST) demonstrated a reduction in lung cancer mortality with lowdose CT (LDCT) compared to chest x-ray (CXR) screening. Overdiagnosis was high (79%) among bronchoalveolar carcinoma (BAC) currently replaced by adenocarcinoma in situ (AIS), minimally invasive adenocarcinoma (MIA) and adenocarcinoma of low malignant potential (LMP) exhibiting 100% disease specific survival (DSS).ObjectiveCompare the outcomes and proportions of BAC, AIS, MIA, and LMP among NLST screendetected stage IA NSCLC with overdiagnosis rate.MethodsWhole slide images were reviewed by a thoracic pathologist from 174 of 409 NLST screen-detected stage IA LUAD. Overdiagnosis rates were calculated from follow-up cancer incidence rates.ResultsMost BAC were reclassified as AIS/MIA/LMP (20/35 = 57%). The 7-year DSS was 100% for AIS/MIA/LMP and 94% for BAC. Excluding AIS/MIA/LMP, BAC behaved similarly to NSCLC (7-year DSS: 86% vs. 83%, P =  0.85) The overdiagnosis rate of LDCT stage IA NSCLC was 16.6% at 11.3-years, matching the proportion of AIS/MIA/LMP (16.2%) but not AIS/MIA (3.5%) or BAC (22.8%).ConclusionsAIS/MIA/LMP proportionally matches the overdiagnosis rate among stage IA NSCLC in the NLST, exhibiting 100% 7-year DSS. Biomarkers designed to recognize AIS/MIA/LMP preoperatively, would be useful to prevent overtreatment of indolent screen-detected cancers.

2024

Lamb, Carla R, Kimberly M Rieger-Christ, Chakravarthy Reddy, Jing Huang, Jie Ding, Marla Johnson, Sean Walsh, et al. (2024) 2024. “A Nasal Swab Classifier to Evaluate The Probability of Lung Cancer in Patients With Pulmonary Nodules.”. Chest 165 (4): 1009-19. https://doi.org/10.1016/j.chest.2023.11.036.

BACKGROUND: Accurate assessment of the probability of lung cancer (pCA) is critical in patients with pulmonary nodules (PNs) to help guide decision-making. We sought to validate a clinical-genomic classifier developed using whole-transcriptome sequencing of nasal epithelial cells from patients with a PN ≤ 30 mm who smoke or have previously smoked.

RESEARCH QUESTION: Can the pCA in individuals with a PN and a history of smoking be predicted by a classifier that uses clinical factors and genomic data from nasal epithelial cells obtained by cytologic brushing?

STUDY DESIGN AND METHODS: Machine learning was used to train a classifier using genomic and clinical features on 1,120 patients with PNs labeled as benign or malignant established by a final diagnosis or a minimum of 12 months of radiographic surveillance. The classifier was designed to yield low-, intermediate-, and high-risk categories. The classifier was validated in an independent set of 312 patients, including 63 patients with a prior history of cancer (other than lung cancer), comparing the classifier prediction with the known clinical outcome.

RESULTS: In the primary validation set, sensitivity and specificity for low-risk classification were 96% and 42%, whereas sensitivity and specificity for high-risk classification was 58% and 90%, respectively. Sensitivity was similar across stages of non-small cell lung cancer, independent of subtype. Performance compared favorably with clinical-only risk models. Analysis of 63 patients with prior cancer showed similar performance as did subanalyses of patients with light vs heavy smoking burden and those eligible for lung cancer screening vs those who were not.

INTERPRETATION: The nasal classifier provides an accurate assessment of pCA in individuals with a PN ≤ 30 mm who smoke or have previously smoked. Classifier-guided decision-making could lead to fewer diagnostic procedures in patients without cancer and more timely treatment in patients with lung cancer.

McNall, Shannon, Kailey Hooper, Travis Sullivan, Kimberly Rieger-Christ, and Matthew Clements. (2024) 2024. “Treatment Modalities for Non-Muscle Invasive Bladder Cancer: An Updated Review.”. Cancers 16 (10). https://doi.org/10.3390/cancers16101843.

The landscape of treatment for non-muscle invasive bladder cancer is rapidly changing. A complete and careful transurethral resection is the mainstay of initial treatment and is followed by intravesical therapy in intermediate or high-risk cases. The standard of care is intravesical BCG. Many alternative or additive approaches to this are being explored. We divided this review into three relevant spaces to consider these novel treatment approaches: (1) low-risk disease, for which intravesical therapy is not usually considered, (2) BCG-naïve disease (i.e., considering alternatives to the standard therapy), and (3) BCG-unresponsive disease. We performed a review of published literature and summarized ongoing trials in the United States. Novel approaches that we explored include surgical techniques for resection, alterations in dwell time for intravesical therapy, delivery method and schedule of intravesical therapies, new intravesical therapy agents, and systemic therapies (especially immunotherapy). These are thoroughly outlined throughout this review article, and the numerous modalities being studied demonstrate significant promise for the future treatment of the expanding space of NMIBC.

Hayden, Joshua P, Adam Wiggins, Travis Sullivan, Thomas Kalantzakos, Kailey Hooper, Alireza Moinzadeh, and Kimberly Rieger-Christ. (2024) 2024. “Use of Droplet Digital Polymerase Chain Reaction to Identify Biomarkers for Differentiation of Benign and Malignant Renal Masses.”. Cancers 16 (4). https://doi.org/10.3390/cancers16040787.

Several microRNAs (miRNAs) have been identified as cell-free biomarkers for detecting renal cell carcinoma (RCC). Droplet digital polymerase chain reaction (ddPCR) is a unique technology for nucleic acid quantification. It has the potential for superior precision, reproducibility, and diagnostic performance in identifying circulating miRNA biomarkers compared to conventional quantitative real-time PCR (qRT-PCR). This study aims to evaluate the performance of ddPCR compared to qRT- PCR in identifying miRNA biomarkers that differentiate malignant from benign renal masses. Potential biomarkers of RCC were identified from a literature review. RNA was extracted from the plasma of 56 patients. All the samples underwent analysis via ddPCR as well as qRT-PCR, and expression levels were recorded for the following miRNAs: miR-93, -144, -210, -221, and -222. Tumors were grouped into low-grade ccRCC, high-grade ccRCC, papillary RCC, and benign masses (primarily angiomyolipoma). The miRNA miR-210 (p = 0.034) and the combination of miRs-210 and miR-222 (p = 0.003) were expressed at significantly higher rates among those with RCC than those with benign masses, as measured by ddPCR. Using the combination of miR-210 and miR-222, ddPCR identified significant differences between the subgroups: papillary RCC versus benign (p = 0.03), low-grade ccRCC versus benign (p = 0.026), and high-grade ccRCC versus benign (p = 0.002). The only significant difference between these subgroups using qRT-PCR was between high-grade ccRCC and benign (p = 0.045). All the AUCs were significant when comparing each RCC subgroup with benign for both PCR technologies. Using a combination of miR-210 and miR-222, ddPCR identified significant differences between benign and malignant renal masses that were not identified as significant by conventional qRT-PCR.

Steiner, Dylan, Lila Sultan, Travis Sullivan, Hanqiao Liu, Sherry Zhang, Ashley LeClerc, Yuriy O Alekseyev, et al. (2024) 2024. “Identification of a Gene Expression Signature of Vascular Invasion and Recurrence in Stage I Lung Adenocarcinoma via Bulk and Spatial Transcriptomics.”. BioRxiv : The Preprint Server for Biology. https://doi.org/10.1101/2024.06.07.597993.

Microscopic vascular invasion (VI) is predictive of recurrence and benefit from lobectomy in stage I lung adenocarcinoma (LUAD) but is difficult to assess in resection specimens and cannot be accurately predicted prior to surgery. Thus, new biomarkers are needed to identify this aggressive subset of stage I LUAD tumors. To assess molecular and microenvironment features associated with angioinvasive LUAD we profiled 162 resected stage I tumors with and without VI by RNA-seq and explored spatial patterns of gene expression in a subset of 15 samples by high-resolution spatial transcriptomics (stRNA-seq). Despite the small size of invaded blood vessels, we identified a gene expression signature of VI from the bulk RNA-seq discovery cohort (n=103) and found that it was associated with VI foci, desmoplastic stroma, and high-grade patterns in our stRNA-seq data. We observed a stronger association with high-grade patterns from VI+ compared with VI- tumors. Using the discovery cohort, we developed a transcriptomic predictor of VI, that in an independent validation cohort (n=60) was associated with VI (AUROC=0.86; p=5.42×10-6) and predictive of recurrence-free survival (HR=1.98; p=0.024), even in VI- LUAD (HR=2.76; p=0.003). To determine our VI predictor's robustness to intra-tumor heterogeneity we used RNA-seq data from multi-region sampling of stage I LUAD cases in TRACERx, where the predictor scores showed high correlation (R=0.87, p<2.2×10-16) between two randomly sampled regions of the same tumor. Our study suggests that VI-associated gene expression changes are detectable beyond the site of intravasation and can be used to predict the presence of VI. This may enable the prediction of angioinvasive LUAD from biopsy specimens, allowing for more tailored medical and surgical management of stage I LUAD.

2023

Suaiti, Lubna, Travis B Sullivan, Kimberly M Rieger-Christ, Elliot L Servais, Kei Suzuki, and Eric J Burks. (2023) 2023. “Vascular Invasion Predicts Recurrence in Stage IA2-IB Lung Adenocarcinoma But Not Squamous Cell Carcinoma.”. Clinical Lung Cancer 24 (3): e126-e133. https://doi.org/10.1016/j.cllc.2022.12.006.

BACKGROUND: Lymphovascular invasion (LVI) is an adverse prognostic feature in resected stage I non-small cell lung cancer (NSCLC); however, it is unclear if the prognostic significance applies to both lung adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC).

MATERIALS AND METHODS: A retrospective review of H&E-stained slides from surgically resected AJCC 8th ed. stage IA2-IB LUAD (n = 344) and LUSC (n = 102) from two institutions was performed. LVI was defined as either lymphatic (LI) or vascular (VI) invasion. Outcomes were assessed by 5-year recurrence-free survival (RFS) estimates using the Kaplan-Meier method.

RESULTS: The cohorts of LUAD and LUSC showed no significant differences in 5-year RFS (81% each), stage, age, race, or surgical procedure. The presence of LVI, VI, and LI was predictive of 5-year RFS for LUAD (LVI + 71% vs. LVI - 92%, P < 0.001; VI + 64% vs. VI - 90%, P < 0.001; LI + 75% vs. LI - 84%, P = 0.030) but not LUSC (LVI + 84% vs. LVI - 79%, P = 0.740; VI + 83% vs. VI- 80%, P = 0.852; LI + 84% vs. LI - 81%, P = 0.757). Among LUAD with LVI, VI was a stronger predictor of 5-year RFS than the remaining subset of VI-LI + tumors (64% vs. 87%, P = 004). Subset analysis of LI among LUAD stratified by VI showed no significant prognostic advantage to adding LI for risk stratification (VI-LI + 87% vs. VI-LI - 92%, P = 0.347 & VI+LI + 62% vs. VI + LI- 66%, P = 0.422). VI was present in 36% of LUAD.

CONCLUSION: Vascular invasion is a strong predictor of recurrence in stage IA2-IB LUAD but not in LUSC. Adjuvant therapy trials should be directed at this subgroup.

Das, Sanjna, Joshua Hayden, Travis Sullivan, and Kimberly Rieger-Christ. (2023) 2023. “The Roles of MiRNAs in Predicting Bladder Cancer Recurrence and Resistance to Treatment.”. International Journal of Molecular Sciences 24 (2). https://doi.org/10.3390/ijms24020964.

Bladder cancer (BCa) is associated with significant morbidity, with development linked to environmental, lifestyle, and genetic causes. Recurrence presents a significant issue and is managed in the clinical setting with intravesical chemotherapy or immunotherapy. In order to address challenges such as a limited supply of BCG and identifying cases likely to recur, it would be advantageous to use molecular biomarkers to determine likelihood of recurrence and treatment response. Here, we review microRNAs (miRNAs) that have shown promise as predictors of BCa recurrence. MiRNAs are also discussed in the context of predicting resistance or susceptibility to BCa treatment.

Jamil, Marcus L, Aaron Perecman, Amanda Sherman, Travis Sullivan, Kimberly Christ, Alexandra Hansma, Eric Burks, and Alex J Vanni. (2023) 2023. “Urinary Microbiome Differences Between Lichen Sclerosus Induced and Non-Lichen Sclerosus Induced Urethral Stricture Disease.”. World Journal of Urology 41 (9): 2495-2501. https://doi.org/10.1007/s00345-023-04490-0.

OBJECTIVE: To describe differences in the urinary microbiome of patients with pathologically confirmed lichen sclerosus (LS) urethral stricture disease (USD) vs non-lichen sclerosus (non-LS) USD pre- and post-operatively.

METHODS: Patients were pre-operatively identified and prospectively followed, all underwent surgical repair and had tissue samples obtained to make a pathological diagnosis of LS. Pre- and post-operative urine samples were collected. Bacterial genomic DNA was extracted. Alpha and beta diversity measurements were calculated and compared. A zero-inflated negative binomial model was utilized to compare taxa abundances between disease status and surgery status.

RESULTS: Urine samples were obtained from both cohorts, 69 samples in total: 36 samples were obtained pre-operatively and 33 samples were obtained post-operatively. Ten patients provided both a pre-operative and post-operative urine sample. Twenty-six patients had pathological evidence of LS and 33 patients did not. There was a statistically significant difference in alpha diversity between the pre-operative urine samples of patients with non-LS USD and LS USD, (p = 0.01). There was no significant difference in alpha diversity within post-operative urine samples between patients with non-LS USD and LS USD, (p = 0.1). A significant difference was observed in Weighed UniFrac distances with respect to disease and operative status, (p = 0.001 and 0.002).

CONCLUSIONS: LS USD have significant alterations in diversity and differential abundance of urine microbiota compared to non-LS USD controls. These findings could be used to guide further investigations into the role of the urinary microbiome in LS USD pathogenesis, severity of presentation, and stricture recurrence.