SifiNet: a robust and accurate method to identify feature gene sets and annotate cells.

SifiNet is a robust and accurate computational pipeline for identifying distinct gene sets, extracting and annotating cellular subpopulations, and elucidating intrinsic relationships among these subpopulations. Uniquely, SifiNet bypasses the cell clustering stage, commonly integrated into other cellular annotation pipelines, thereby circumventing potential inaccuracies in clustering that may compromise subsequent analyses. Consequently, SifiNet has demonstrated superior performance in multiple experimental datasets compared with other state-of-the-art methods. SifiNet can analyze both single-cell RNA and ATAC sequencing data, thereby rendering comprehensive multi-omic cellular profiles. It is conveniently available as an open-source R package.

The TLR2/TLR6 ligand FSL-1 mitigates radiation-induced hematopoietic injury in mice and nonhuman primates.

Thrombocytopenia, hemorrhage, anemia, and infection are life-threatening issues following accidental or intentional radiation exposure. Since few therapeutics are available, safe and efficacious small molecules to mitigate radiation-induced injury need to be developed. Our previous study showed the synthetic TLR2/TLR6 ligand fibroblast stimulating lipopeptide (FSL-1) prolonged survival and provided MyD88-dependent mitigation of hematopoietic acute radiation syndrome (H-ARS) in mice. Although mice and humans differ in TLR number, expression, and function, nonhuman primate (NHP) TLRs are like those of humans; therefore, studying both animal models is critical for drug development. The objectives of this study were to determine the efficacy of FSL-1 on hematopoietic recovery in small and large animal models subjected to sublethal total body irradiation and investigate its mechanism of action. In mice, we demonstrate a lack of adverse effects, an easy route of delivery (subcutaneous) and efficacy in promoting hematopoietic progenitor cell proliferation by FSL-1. NHP given radiation, followed a day later with a single subcutaneous administration of FSL-1, displayed no adversity but showed elevated hematopoietic cells. Our analyses revealed that FSL-1 promoted red blood cell development and induced soluble effectors following radiation exposure. Cytologic analysis of bone marrow aspirates revealed a striking enhancement of mononuclear progenitor cells in FSL-1-treated NHP. Combining the efficacy of FSL-1 in promoting hematopoietic cell recovery with the lack of adverse effects induced by a single administration supports the application of FSL-1 as a viable countermeasure against H-ARS.

DYNATE: Localizing rare-variant association regions via multiple testing embedded in an aggregation tree.

Rare-variants (RVs) genetic association studies enable researchers to uncover the variation in phenotypic traits left unexplained by common variation. Traditional single-variant analysis lacks power; thus, researchers have developed various methods to aggregate the effects of RVs across genomic regions to study their collective impact. Some existing methods utilize a static delineation of genomic regions, often resulting in suboptimal effect aggregation, as neutral subregions within the test region will result in an attenuation of signal. Other methods use varying windows to search for signals but often result in long regions containing many neutral RVs. To pinpoint short genomic regions enriched for disease-associated RVs, we developed a novel method, DYNamic Aggregation TEsting (DYNATE). DYNATE dynamically and hierarchically aggregates smaller genomic regions into larger ones and performs multiple testing for disease associations with a controlled weighted false discovery rate. DYNATE's main advantage lies in its strong ability to identify short genomic regions highly enriched for disease-associated RVs. Extensive numerical simulations demonstrate the superior performance of DYNATE under various scenarios compared with existing methods. We applied DYNATE to an amyotrophic lateral sclerosis study and identified a new gene, EPG5, harboring possibly pathogenic mutations.

Common variation in a long non-coding RNA gene modulates variation of circulating TGF-ß2 levels in metastatic colorectal cancer patients (Alliance).

BACKGROUND: Herein, we report results from a genome-wide study conducted to identify protein quantitative trait loci (pQTL) for circulating angiogenic and inflammatory protein markers in patients with metastatic colorectal cancer (mCRC). The study was conducted using genotype, protein marker, and baseline clinical and demographic data from CALGB/SWOG 80405 (Alliance), a randomized phase III study designed to assess outcomes of adding VEGF or EGFR inhibitors to systemic chemotherapy in mCRC patients. Germline DNA derived from blood was genotyped on whole-genome array platforms. The abundance of protein markers was quantified using a multiplex enzyme-linked immunosorbent assay from plasma derived from peripheral venous blood collected at baseline. A robust rank-based method was used to assess the statistical significance of each variant and protein pair against a strict genome-wide level. A given pQTL was tested for validation in two external datasets of prostate (CALGB 90401) and pancreatic cancer (CALGB 80303) patients. Bioinformatics analyses were conducted to further establish biological bases for these findings. RESULTS: The final analysis was carried out based on data from 540,021 common typed genetic variants and 23 protein markers from 869 genetically estimated European patients with mCRC. Correcting for multiple testing, the analysis discovered a novel cis-pQTL in LINC02869, a long non-coding RNA gene, for circulating TGF-ß2 levels (rs11118119; AAF = 0.11; P-value < 1.4e-14). This finding was validated in a cohort of 538 prostate cancer patients from CALGB 90401 (AAF = 0.10, P-value < 3.3e-25). The analysis also validated a cis-pQTL we had previously reported for VEGF-A in advanced pancreatic cancer, and additionally identified trans-pQTLs for VEGF-R3, and cis-pQTLs for CD73. CONCLUSIONS: This study has provided evidence of a novel cis germline genetic variant that regulates circulating TGF-ß2 levels in plasma of patients with advanced mCRC and prostate cancer. Moreover, the validation of previously identified pQTLs for VEGF-A, CD73, and VEGF-R3, potentiates the validity of these associations.

Pharmacogenetic and clinical risk factors for bevacizumab-related gastrointestinal hemorrhage in prostate cancer patients treated on CALGB 90401 (Alliance).

The objective of this study was to discover clinical and pharmacogenetic factors associated with bevacizumab-related gastrointestinal hemorrhage in Cancer and Leukemia Group B (Alliance) 90401. Patients with metastatic castration-resistant prostate cancer received docetaxel and prednisone ± bevacizumab. Patients were genotyped using Illumina HumanHap610-Quad and assessed using cause-specific risk for association between single nucleotide polymorphisms (SNPs) and gastrointestinal hemorrhage. In 1008 patients, grade 2 or higher gastrointestinal hemorrhage occurred in 9.5% and 3.8% of bevacizumab (n = 503) and placebo (n = 505) treated patients, respectively. Bevacizumab (P < 0.001) and age (P = 0.002) were associated with gastrointestinal hemorrhage. In 616 genetically estimated Europeans (n = 314 bevacizumab and n = 302 placebo treated patients), grade 2 or higher gastrointestinal hemorrhage occurred in 9.6% and 2.0% of patients, respectively. One SNP (rs1478947; HR 6.26; 95% CI 3.19-12.28; P = 9.40 × 10-8) surpassed Bonferroni-corrected significance. Grade 2 or higher gastrointestinal hemorrhage rate was 33.3% and 6.2% in bevacizumab-treated patients with the AA/AG and GG genotypes, versus 2.9% and 1.9% in the placebo arm, respectively. Prospective validation of these findings and functional analyses are needed to better understand the genetic contribution to treatment-related gastrointestinal hemorrhage.

High-fat diet and obesity are associated with differential angiogenic gene expression in epithelial ovarian cancer.

OBJECTIVE: We sought to evaluate the association between diet and angiogenic biomarkers in KpB mice, and the association between these markers, body mass index (BMI), and overall survival (OS) in high-grade serous cancers (HGSC). METHODS: Tumors previously obtained from KpB mice subjected to high-fat diets (HFD, n = 10) or low-fat diets (LFD, n = 10) were evaluated for angiogenesis based on CD-31 microvessel density (MVD). Data from prior microarray analysis (Agilent 244 K arrays) conducted in 10 mice were utilized to assess associations between diet and angiogenetic biomarkers. Agilent (mouse) and Affymetrix Human Genome U133a probes were linked to 162 angiogenic-related genes. The associations between biomarkers, BMI, and OS were evaluated in an HGSC internal database (IDB) (n = 40). Genes with unadjusted p < 0.05 were evaluated for association with OS in the TCGA-OV database (n = 339). RESULTS: There was no association between CD-31 and diet in mice (p = 0.66). Sixteen angiogenic-related genes passed the p < 0.05 threshold for association with HFD vs. LFD. Transforming growth factor-alpha (TGFA) demonstrated 72% higher expression in HFD vs. LFD mice (p = 0.04). Similar to the mouse study, in our HGSC IDB, higher TGFA expression correlated with higher BMI (p = 0.01) and shorter survival (p = 0.001). In the TCGA-OV dataset, BMI data was not available and there was no association between TGFA and OS (p = 0.48). CONCLUSIONS: HFD and obesity may promote tumor progression via differential modulation of TGFA. We were unable to confirm this finding in the TCGA dataset. Further evaluation of TGFA is needed to determine if this is a target unique to obesity-driven HGSC.

Genome-wide association studies of survival in 1520 cancer patients treated with bevacizumab-containing regimens.

Germline variants might predict cancer progression. Bevacizumab improves overall survival (OS) in patients with advanced cancers. No biomarkers are available to identify patients that benefit from bevacizumab. A meta-analysis of genome-wide association studies (GWAS) was conducted in 1,520 patients from Phase III trials (CALGB 80303, 40503, 80405 and ICON7), where bevacizumab was randomized to treatment without bevacizumab. We aimed to identify genes and single nucleotide polymorphisms (SNPs) associated with survival independently of bevacizumab treatment or through interaction with bevacizumab. A cause-specific Cox model was used to test the SNP-OS association in both arms combined (prognostic), and the effect of SNPs-bevacizumab interaction on OS (predictive) in each study. The SNP effects across studies were combined using inverse variance. Findings were tested for replication in advanced colorectal and ovarian cancer patients from The Cancer Genome Atlas (TGCA). In the GWAS meta-analysis, patients with rs680949 in PRUNE2 experienced shorter OS compared to patients without it (P = 1.02 × 10-7 , hazard ratio [HR] = 1.57, 95% confidence interval [CI] 1.33-1.86), as well as in TCGA (P = .0219, HR = 1.58, 95% CI 1.07-2.35). In the GWAS meta-analysis, patients with rs16852804 in BARD1 experienced shorter OS compared to patients without it (P = 1.40 × 10-5 , HR = 1.51, 95% CI 1.25-1.82) as well as in TCGA (P = 1.39 × 10-4 , HR = 3.09, 95% CI 1.73-5.51). Patients with rs3795897 in AGAP1 experienced shorter OS in the bevacizumab arm compared to the nonbevacizumab arm (P = 1.43 × 10-5 ). The largest GWAS meta-analysis of bevacizumab treated patients identified PRUNE2 and BARD1 (tumor suppressor genes) as prognostic genes of colorectal and ovarian cancer, respectively, and AGAP1 as a potentially predictive gene that interacts with bevacizumab with respect to patient survival.

Bevacizumab-induced hypertension and proteinuria: a genome-wide study of more than 1000 patients.

BACKGROUND: Hypertension and proteinuria are common bevacizumab-induced toxicities. No validated biomarkers are available for identifying patients at risk of these toxicities. METHODS: A genome-wide association study (GWAS) meta-analysis was performed in 1039 bevacizumab-treated patients of European ancestry in four clinical trials (CALGB 40502, 40503, 80303, 90401). Grade ≥2 hypertension and proteinuria were recorded (CTCAE v.3.0). Single-nucleotide polymorphism (SNP)-toxicity associations were determined using a cause-specific Cox model adjusting for age and sex. RESULTS: The most significant SNP associated with hypertension with concordant effect in three out of the four studies (p-value <0.05 for each study) was rs6770663 (A > G) in KCNAB1, with the G allele increasing the risk of hypertension (p-value = 4.16 × 10-6). The effect of the G allele was replicated in ECOG-ACRIN E5103 in 582 patients (p-value = 0.005). The meta-analysis of all five studies for rs6770663 led to p-value = 7.73 × 10-8, close to genome-wide significance. The most significant SNP associated with proteinuria was rs339947 (C > A, between DNAH5 and TRIO), with the A allele increasing the risk of proteinuria (p-value = 1.58 × 10-7). CONCLUSIONS: The results from the largest study of bevacizumab toxicity provide new markers of drug safety for further evaluations. SNP in KCNAB1 validated in an independent dataset provides evidence toward its clinical applicability to predict bevacizumab-induced hypertension. ClinicalTrials.gov Identifier: NCT00785291 (CALGB 40502); NCT00601900 (CALGB 40503); NCT00088894 (CALGB 80303) and NCT00110214 (CALGB 90401).

Differential alternative RNA splicing and transcription events between tumors from African American and White patients in The Cancer Genome Atlas.

Individuals of African ancestry suffer disproportionally from higher incidence, aggressiveness, and mortality for particular cancers. This disparity likely results from an interplay among differences in multiple determinants of health, including differences in tumor biology. We used The Cancer Genome Atlas (TCGA) SpliceSeq and TCGA aggregate expression datasets and identified differential alternative RNA splicing and transcription events (ARS/T) in cancers between self-identified African American (AA) and White (W) patients. We found that retained intron events were enriched among race-related ARS/T. In addition, on average, 12% of the most highly ranked race-related ARS/T overlapped between any two analyzed cancers. Moreover, the genes undergoing race-related ARS/T functioned in cancer-promoting pathways, and a number of race-related ARS/T were associated with patient survival. We built a web-application, CanSplice, to mine genomic datasets by self-identified race. The race-related targets have the potential to aid in the development of new biomarkers and therapeutics to mitigate cancer disparity.

Loss of ß-catenin triggers oxidative stress and impairs hematopoietic regeneration.

Accidental or deliberate ionizing radiation exposure can be fatal due to widespread hematopoietic destruction. However, little is known about either the course of injury or the molecular pathways that regulate the subsequent regenerative response. Here we show that the Wnt signaling pathway is critically important for regeneration after radiation-induced injury. Using Wnt reporter mice, we show that radiation triggers activation of Wnt signaling in hematopoietic stem and progenitor cells. ß-Catenin-deficient mice, which lack the ability to activate canonical Wnt signaling, exhibited impaired hematopoietic stem cell regeneration and bone marrow recovery after radiation. We found that, as part of the mechanism, hematopoietic stem cells lacking ß-catenin fail to suppress the generation of reactive oxygen species and cannot resolve DNA double-strand breaks after radiation. Consistent with the impaired response to radiation, ß-catenin-deficient mice are also unable to recover effectively after chemotherapy. Collectively, these data indicate that regenerative responses to distinct hematopoietic injuries share a genetic dependence on ß-catenin and raise the possibility that modulation of Wnt signaling may be a path to improving bone marrow recovery after damage.

A prospective trial of abiraterone acetate plus prednisone in Black and White men with metastatic castrate-resistant prostate cancer.

BACKGROUND: Retrospective analyses of randomized trials suggest that Black men with metastatic castration-resistant prostate cancer (mCRPC) have longer survival than White men. The authors conducted a prospective study of abiraterone acetate plus prednisone to explore outcomes by race. METHODS: This race-stratified, multicenter study estimated radiographic progression-free survival (rPFS) in Black and White men with mCRPC. Secondary end points included prostate-specific antigen (PSA) kinetics, overall survival (OS), and safety. Exploratory analysis included genome-wide genotyping to identify single nucleotide polymorphisms associated with progression in a model incorporating genetic ancestry. One hundred patients self-identified as White (n = 50) or Black (n = 50) were enrolled. Eligibility criteria were modified to facilitate the enrollment of individual Black patients. RESULTS: The median rPFS for Black and White patients was 16.6 and 16.8 months, respectively; their times to PSA progression (TTP) were 16.6 and 11.5 months, respectively; and their OS was 35.9 and 35.7 months, respectively. Estimated rates of PSA decline by ≥50% in Black and White patients were 74% and 66%, respectively; and PSA declines to <0.2 ng/mL were 26% and 10%, respectively. Rates of grade 3 and 4 hypertension, hypokalemia, and hyperglycemia were higher in Black men. CONCLUSIONS: Multicenter prospective studies by race are feasible in men with mCRPC but require less restrictive eligibility. Despite higher comorbidity rates, Black patients demonstrated rPFS and OS similar to those of White patients and trended toward greater TTP and PSA declines, consistent with retrospective reports. Importantly, Black men may have higher side-effect rates than White men. This exploratory genome-wide analysis of TTP identified a possible candidate marker of ancestry-dependent treatment outcomes.

Pharmacogenetic study in gastric cancer patients treated with adjuvant fluorouracil/leucovorin or epirubicin/cisplatin/fluorouracil before and after chemoradiation on CALGB 80101 (Alliance).

There is a lack of pharmacogenetic predictors of outcome in gastric cancer patients. The aim of this study was to assess previously identified candidate genes associated with 5-fluorouracil (5-FU), cisplatin, or epirubicin toxicity or response in a cohort of resected gastric cancer patients treated on CALGB (Alliance) 80101. Gastric or gastroesophageal cancer patients randomized to adjuvant 5-FU/leucovorin or epirubicin/cisplatin/5-FU before and after 5-FU chemoradiation were genotyped for single nucleotide polymorphisms (SNPs) in GSTP1 (rs1695), ERCC1 (rs11615 and rs3212986), XRCC1 (rs25487), UGT2B7 (rs7439366) and the 28 base-pair tandem repeats in TYMS (rs34743033). Logistic regression and log rank tests were used to assess the association between each SNP and incidence of grade 3/4 neutropenia and leukopenia, overall (OS) and progression-free survival (PFS), respectively. Toxicity endpoint analyses were adjusted for the treatment arm, while OS and PFS were also adjusted for performance status, sex, age, lymph node involvement, and primary tumor site and size. Of 281 subjects with successful genotyping results and available clinical (toxicity and efficacy) data, 166 self-reported non-Hispanic White patients were included in the final analysis. There was a lack of evidence of an association among any SNPs tested with grade 3/4 neutropenia and leukopenia or OS and PFS. Age, lymph node involvement, and primary tumor size were significantly associated with OS and PFS. This study failed to confirm results of previous gastric cancer pharmacogenetic studies.

Clinical Results and Biomarker Analyses of Axitinib and TRC105 versus Axitinib Alone in Patients with Advanced or Metastatic Renal Cell Carcinoma (TRAXAR).

LESSONS LEARNED: The combination of carotuximab with axitinib did not provide a benefit over axitinib monotherapy in patients with metastatic clear cell renal cell carcinoma who had previously progressed on one or more vascular endothelial growth factor (VEGF)-targeted therapies. Exploratory evaluation of pretreatment circulating biomarkers suggested the combination might benefit patients who have low baseline VEGF levels. BACKGROUND: Endoglin is an angiogenic receptor expressed on proliferating tumor vessels and renal cell carcinoma (RCC) stem cells that is implicated as a mechanism of resistance to vascular endothelial growth factor receptor (VEGFR) inhibitors. This study evaluated an antiendoglin monoclonal antibody (carotuximab, TRC105) combined with axitinib in patients with advanced or metastatic clear cell renal cell carcinoma (mccRCC) who had progressed following one or more prior VEGF inhibitors. METHODS: TRAXAR was a multicenter, international randomized 1:1 (stratified by ECOG, 0 vs. 1), phase II study of carotuximab combined with axitinib versus axitinib alone in mccRCC patients who had progressed following one or more prior VEGF inhibitors. The primary endpoint was progression-free survival (PFS) assessed by independent central review (ICR) per RECIST 1.1 RESULTS: A total of 150 patients were randomized. The combination therapy resulted in shorter median PFS by RECIST 1.1 than axitinib monotherapy (6.7 vs. 11.4 months). The combination was tolerated similarly to axitinib monotherapy, and there were no treatment related deaths. Exploratory evaluation of pretreatment circulating biomarkers suggested the combination might benefit patients who have low baseline VEGF levels. CONCLUSION: The combination of carotuximab with axitinib did not demonstrate additional efficacy over single agent axitinib in patients with mccRCC who progressed following one or more prior VEGF inhibitor treatment.

Obesity and altered angiogenic-related gene expression in endometrial cancer.

OBJECTIVES: Evaluate association between obesity and angiogenic-related gene expression in endometrial cancer (EC). Evaluate interaction between diet and metformin on angiogenic-related gene expression. METHODS: We evaluated the association between 168 human angiogenic-related genes and body mass index (BMI) in the TCGA Uterine Corpus Endometrial Carcinoma cohort (endometrioid endometrial cancer (EEC) cohort n = 290, and copy number high cohort n = 55), an independent validation cohort from Gynecologic Cancer Center of Excellence (GYN-COE) (n = 62) and corresponding 185 homologous mouse genes in an LKB1fl/flp53fl/fl mouse model of EC (n = 20). Mice received 60% of calories from fat in a high-fat diet (HFD), mimicking diet-induced obesity, versus 10% of calories from fat in a low-fat diet (LFD). After tumor growth, HFD (n = 5) and LFD (n = 5) mice were treated with metformin (200 mg/kg/day) or control. Whole transcriptome analysis of mouse tumors was performed using RNA-Seq. RESULTS: At a false-discovery rate of 10%, twenty-one angiogenic-related genes were differentially expressed with respect to BMI when adjusting for grade in the TCGA EEC cohort. Evaluation of these genes in the mouse model control group revealed association between increased Edil3 expression in HFD versus LFD mice (2.5-fold change (FC); unadjusted p = 0.03). An interaction was observed for expression of Edil3 between diet and metformin treatment (unadjusted p = 0.009). Association between BMI and increased expression of EDIL3 was validated in one of four EDIL3 probesets in the GYN-COE cohort (p = 0.0011, adjusted p = 0.0342). CONCLUSIONS: Obesity may promote tumor progression via differential modulation of angiogenic pathways in EEC. Our exploratory findings demonstrated that EDIL3 may be a candidate gene of interest.

Long-Term Follow-up of CALGB (Alliance) 100001: Autologous Followed by Nonmyeloablative Allogeneic Transplant for Multiple Myeloma.

CALGB (Alliance) 100001 was a phase II study evaluating autologous stem cell transplant (ASCT) followed by nonmyeloablative allogeneic stem cell transplant (alloSCT) in patients with multiple myeloma who had received no more than 18 months of prior therapy and had experienced no more than 1 prior progression event. Conditioning for ASCT was with high-dose melphalan (200 mg/m2). The alloSCT reduced-intensity conditioning (RIC) regimen consisted of fludarabine (30 mg/m2/d i.v. on days -7 through -3) and cyclophosphamide (1 g/m2/d i.v. on days -4 through -3). The primary objective was to determine the 6-month post-alloSCT treatment-related mortality (TRM) rate. Additional objectives included determining the proportion of patients who could complete this tandem ASCT-alloSCT approach in a cooperative group setting, overall response rates, rates of donor chimerism, rates of graft-versus-host disease (GVHD), disease-free survival, and overall survival (OS). Sixty patients were enrolled, of whom 57 (95%) completed ASCT and 49 (82%) completed tandem ASCT-alloSCT. The TRM rate was 2% (1/49; 90% confidence interval, 0.10% to 9.3%). Moderate to severe (grades 2 to 3) acute GVHD was observed in 13 of 49 alloSCT patients (27%). One patient died due to GVHD within 9 months of alloSCT. Twenty-seven of the 49 patients (55%) who underwent alloSCT reported chronic GVHD as either limited (15/49; 31%) or extensive (12/49; 24%) in the first year post-alloSCT and prior to the start of nonprotocol therapy for progressive disease. With a median follow-up for survival of 11 years, the median OS time is 6.6 years and the median time to disease progression is 3.6 years. Similar to other studies, this study confirmed that tandem ASCT/alloSCT is associated with durable disease control in a subset of patients. This study demonstrated the feasibility of performing tandem ASCT/alloSCT in a cooperative group setting and determined that a fludarabine/cyclophosphamide RIC regimen is associated with a very low TRM rate.

Differential expression of immune related genes in high-grade ovarian serous carcinoma.

OBJECTIVE: To identify novel immunologic targets and biomarkers associated with overall survival (OS) in high-grade serous ovarian cancer (HGSC). METHODS: In this retrospective study, microarray data from 51 HGSC specimens were analyzed (Affymetrix HG-U133A). A panel of 183 immune/inflammatory response related genes linked to 279 probe sets was constructed a priori and screened. Associations between gene expression and OS were assessed using logrank tests. Multiple testing was addressed within the False Discovery Rate (FDR) framework. For external validation, TCGA Ovarian dataset and five GSE publicly available HGSC datasets were evaluated. RESULTS: In Duke data, 110 probe sets linked to 83 immunologic/inflammatory-related genes were differentially expressed in tumors from long versus short-term HGSC survivors (adjusted p < 0.05). In TCGA, concordant with the results from the Duke discovery cohort, high expression of one probe (IL6R) demonstrated a consistent significance and concordant association with higher expression in long-term HGSC survivors (Duke q-value = 0.022) and improved OS in the TCGA dataset (p-value = 0.015, HR = 0.8). Thirteen genes in GSE14764 (N = 4) and GSE26712 (N = 9) datasets had significant p-values and consistent concordant with Duke Data. Despite the significant associations of gene expression and OS in the individual GSE datasets, in the GSE meta-analysis no genes were consistently concordant and significantly associated with survival. CONCLUSIONS: Evaluation of IL6R expression may be warranted based on higher expression in long-term survivors and association with improved survival in advanced HGSC. The other candidate genes may also be of worthy of further exploration to enhance immuno-oncology drug discovery.

Predictive accuracy of markers or risk scores for interval censored survival data.

Methods for the evaluation of the predictive accuracy of biomarkers with respect to survival outcomes subject to right censoring have been discussed extensively in the literature. In cancer and other diseases, survival outcomes are commonly subject to interval censoring by design or due to the follow up schema. In this article, we present an estimator for the area under the time-dependent receiver operating characteristic (ROC) curve for interval censored data based on a nonparametric sieve maximum likelihood approach. We establish the asymptotic properties of the proposed estimator and illustrate its finite-sample properties using a simulation study. The application of our method is illustrated using data from a cancer clinical study. An open-source R package to implement the proposed method is available on Comprehensive R Archive Network.

Plasmonic nanobiosensors for detection of microRNA cancer biomarkers in clinical samples.

MicroRNAs (miRNAs) play an important role in the regulation of biological processes and have demonstrated great potential as biomarkers for the early detection of various diseases, including esophageal adenocarcinoma (EAC) and Barrett's esophagus (BE), the premalignant metaplasia associated with EAC. Herein, we demonstrate the direct detection of the esophageal cancer biomarker, miR-21, in RNA extracted from 17 endoscopic tissue biopsies using the nanophotonics technology our group has developed, termed the inverse molecular sentinel (iMS) nanobiosensor, with surface-enhanced Raman scattering (SERS) detection. The potential of this label-free, homogeneous biosensor for cancer diagnosis without the need for target amplification was demonstrated by discriminating esophageal cancer and Barrett's esophagus from normal tissue with notable diagnostic accuracy. This work establishes the potential of the iMS nanobiosensor for cancer diagnostics via miRNA detection in clinical samples without the need for target amplification, validating the potential of this assay as part of a new diagnostic strategy. Combining miRNA diagnostics with the nanophotonics technology will result in a paradigm shift in achieving a general molecular analysis tool that has widespread applicability for cancer research as well as detection of cancer. We anticipate further development of this technique for future use in point-of-care testing as an alternative to histopathological diagnosis as our method provides a quick result following RNA isolation, allowing for timely treatment.

fastJT: An R package for robust and efficient feature selection for machine learning and genome-wide association studies.

BACKGROUND: Parametric feature selection methods for machine learning and association studies based on genetic data are not robust with respect to outliers or influential observations. While rank-based, distribution-free statistics offer a robust alternative to parametric methods, their practical utility can be limited, as they demand significant computational resources when analyzing high-dimensional data. For genetic studies that seek to identify variants, the hypothesis is constrained, since it is typically assumed that the effect of the genotype on the phenotype is monotone (e.g., an additive genetic effect). Similarly, predictors for machine learning applications may have natural ordering constraints. Cross-validation for feature selection in these high-dimensional contexts necessitates highly efficient computational algorithms for the robust evaluation of many features. RESULTS: We have developed an R extension package, fastJT, for conducting genome-wide association studies and feature selection for machine learning using the Jonckheere-Terpstra statistic for constrained hypotheses. The kernel of the package features an efficient algorithm for calculating the statistics, replacing the pairwise comparison and counting processes with a data sorting and searching procedure, reducing computational complexity from O(n2) to O(n log(n)). The computational efficiency is demonstrated through extensive benchmarking, and example applications to real data are presented. CONCLUSIONS: fastJT is an open-source R extension package, applying the Jonckheere-Terpstra statistic for robust feature selection for machine learning and association studies. The package implements an efficient algorithm which leverages internal information among the samples to avoid unnecessary computations, and incorporates shared-memory parallel programming to further boost performance on multi-core machines.

Efficient estimation of grouped survival models.

BACKGROUND: Time- and dose-to-event phenotypes used in basic science and translational studies are commonly measured imprecisely or incompletely due to limitations of the experimental design or data collection schema. For example, drug-induced toxicities are not reported by the actual time or dose triggering the event, but rather are inferred from the cycle or dose to which the event is attributed. This exemplifies a prevalent type of imprecise measurement called grouped failure time, where times or doses are restricted to discrete increments. Failure to appropriately account for the grouped nature of the data, when present, may lead to biased analyses. RESULTS: We present groupedSurv, an R package which implements a statistically rigorous and computationally efficient approach for conducting genome-wide analyses based on grouped failure time phenotypes. Our approach accommodates adjustments for baseline covariates, and analysis at the variant or gene level. We illustrate the statistical properties of the approach and computational performance of the package by simulation. We present the results of a reanalysis of a published genome-wide study to identify common germline variants associated with the risk of taxane-induced peripheral neuropathy in breast cancer patients. CONCLUSIONS: groupedSurv enables fast and rigorous genome-wide analysis on the basis of grouped failure time phenotypes at the variant, gene or pathway level. The package is freely available under a public license through the Comprehensive R Archive Network.