Segmentation of human functional tissue units in support of a Human Reference Atlas.

The Human BioMolecular Atlas Program (HuBMAP) aims to compile a Human Reference Atlas (HRA) for the healthy adult body at the cellular level. Functional tissue units (FTUs), relevant for HRA construction, are of pathobiological significance. Manual segmentation of FTUs does not scale; highly accurate and performant, open-source machine-learning algorithms are needed. We designed and hosted a Kaggle competition that focused on development of such algorithms and 1200 teams from 60 countries participated. We present the competition outcomes and an expanded analysis of the winning algorithms on additional kidney and colon tissue data, and conduct a pilot study to understand spatial location and density of FTUs across the kidney. The top algorithm from the competition, Tom, outperforms other algorithms in the expanded study, while using fewer computational resources. Tom was added to the HuBMAP infrastructure to run kidney FTU segmentation at scale-showcasing the value of Kaggle competitions for advancing research.

Organization of the human intestine at single-cell resolution.

The intestine is a complex organ that promotes digestion, extracts nutrients, participates in immune surveillance, maintains critical symbiotic relationships with microbiota and affects overall health1. The intesting has a length of over nine metres, along which there are differences in structure and function2. The localization of individual cell types, cell type development trajectories and detailed cell transcriptional programs probably drive these differences in function. Here, to better understand these differences, we evaluated the organization of single cells using multiplexed imaging and single-nucleus RNA and open chromatin assays across eight different intestinal sites from nine donors. Through systematic analyses, we find cell compositions that differ substantially across regions of the intestine and demonstrate the complexity of epithelial subtypes, and find that the same cell types are organized into distinct neighbourhoods and communities, highlighting distinct immunological niches that are present in the intestine. We also map gene regulatory differences in these cells that are suggestive of a regulatory differentiation cascade, and associate intestinal disease heritability with specific cell types. These results describe the complexity of the cell composition, regulation and organization for this organ, and serve as an important reference map for understanding human biology and disease.

Single-cell analyses define a continuum of cell state and composition changes in the malignant transformation of polyps to colorectal cancer.

To chart cell composition and cell state changes that occur during the transformation of healthy colon to precancerous adenomas to colorectal cancer (CRC), we generated single-cell chromatin accessibility profiles and single-cell transcriptomes from 1,000 to 10,000 cells per sample for 48 polyps, 27 normal tissues and 6 CRCs collected from patients with or without germline APC mutations. A large fraction of polyp and CRC cells exhibit a stem-like phenotype, and we define a continuum of epigenetic and transcriptional changes occurring in these stem-like cells as they progress from homeostasis to CRC. Advanced polyps contain increasing numbers of stem-like cells, regulatory T cells and a subtype of pre-cancer-associated fibroblasts. In the cancerous state, we observe T cell exhaustion, RUNX1-regulated cancer-associated fibroblasts and increasing accessibility associated with HNF4A motifs in epithelia. DNA methylation changes in sporadic CRC are strongly anti-correlated with accessibility changes along this continuum, further identifying regulatory markers for molecular staging of polyps.

Global, distinctive, and personal changes in molecular and microbial profiles by specific fibers in humans.

Dietary fibers act through the microbiome to improve cardiovascular health and prevent metabolic disorders and cancer. To understand the health benefits of dietary fiber supplementation, we investigated two popular purified fibers, arabinoxylan (AX) and long-chain inulin (LCI), and a mixture of five fibers. We present multiomic signatures of metabolomics, lipidomics, proteomics, metagenomics, a cytokine panel, and clinical measurements on healthy and insulin-resistant participants. Each fiber is associated with fiber-dependent biochemical and microbial responses. AX consumption associates with a significant reduction in LDL and an increase in bile acids, contributing to its observed cholesterol reduction. LCI is associated with an increase in Bifidobacterium. However, at the highest LCI dose, there is increased inflammation and elevation in the liver enzyme alanine aminotransferase. This study yields insights into the effects of fiber supplementation and the mechanisms behind fiber-induced cholesterol reduction, and it shows effects of individual, purified fibers on the microbiome.

Age-associated genes in human mammary gland drive human breast cancer progression.

BACKGROUND: Aging is a comorbidity of breast cancer suggesting that aging-associated transcriptome changes may promote breast cancer progression. However, the mechanism underlying the age effect on breast cancer remains poorly understood. METHOD: We analyzed transcriptomics of the matched normal breast tissues from the 82 breast cancer patients in The Cancer Genome Atlas (TCGA) dataset with linear regression for genes with age-associated expression that are not associated with menopause. We also analyzed differentially expressed genes between the paired tumor and non-tumor breast tissues in TCGA for the identification of age and breast cancer (ABC)-associated genes. A few of these genes were selected for further investigation of their malignancy-regulating activities with in vitro and in vivo assays. RESULTS: We identified 148 upregulated and 189 downregulated genes during aging. Overlapping of tumor-associated genes between normal and tumor tissues with age-dependent genes resulted in 14 upregulated and 24 downregulated genes that were both age and breast cancer associated. These genes are predictive in relapse-free survival, indicative of their potential tumor promoting or suppressive functions, respectively. Knockdown of two upregulated genes (DYNLT3 and P4HA3) or overexpression of the downregulated ALX4 significantly reduced breast cancer cell proliferation, migration, and clonogenicity. Moreover, knockdown of P4HA3 reduced growth and metastasis whereas overexpression of ALX4 inhibited the growth of xenografted breast cancer cells in mice. CONCLUSION: Our study suggests that transcriptome alterations during aging may contribute to breast tumorigenesis. DYNLT3, P4HA3, and ALX4 play significant roles in breast cancer progression.

Spatial EGFR Dynamics and Metastatic Phenotypes Modulated by Upregulated EphB2 and Src Pathways in Advanced Prostate Cancer.

Advanced prostate cancer is a very heterogeneous disease reflecting in diverse regulations of oncogenic signaling pathways. Aberrant spatial dynamics of epidermal growth factor receptor (EGFR) promote their dimerization and clustering, leading to constitutive activation in oncogenesis. The EphB2 and Src signaling pathways are associated with the reorganization of the cytoskeleton leading to malignancy, but their roles in regulating EGFR dynamics and activation are scarcely reported. Using single-particle tracking techniques, we found that highly phosphorylated EGFR in the advanced prostate cancer cell line, PC3, was associated with higher EGFR diffusivity, as compared with LNCaP and less aggressive DU145. The increased EGFR activation and biophysical dynamics were consistent with high proliferation, migration, and invasion. After performing single-cell RNA-seq on prostate cancer cell lines and circulating tumor cells from patients, we identified that upregulated gene expression in the EphB2 and Src pathways are associated with advanced malignancy. After dasatinib treatment or siRNA knockdowns of EphB2 or Src, the PC3 cells exhibited significantly lower EGFR dynamics, cell motility, and invasion. Partial inhibitory effects were also found in DU145 cells. The upregulation of parts of the EphB2 and Src pathways also predicts poor prognosis in the prostate cancer patient cohort of The Cancer Genome Atlas. Our results provide evidence that overexpression of the EphB2 and Src signaling pathways regulate EGFR dynamics and cellular aggressiveness in some advanced prostate cancer cells.

BRCA1 Interacting Protein COBRA1 Facilitates Adaptation to Castrate-Resistant Growth Conditions.

COBRA1 (co-factor of BRCA1) is one of the four subunits of the negative elongation factor originally identified as a BRCA1-interacting protein. Here, we provide first-time evidence for the oncogenic role of COBRA1 in prostate pathogenesis. COBRA1 is aberrantly expressed in prostate tumors. It positively influences androgen receptor (AR) target gene expression and promoter activity. Depletion of COBRA1 leads to decreased cell viability, proliferation, and anchorage-independent growth in prostate cancer cell lines. Conversely, overexpression of COBRA1 significantly increases cell viability, proliferation, and anchorage-independent growth over the higher basal levels. Remarkably, AR-positive androgen dependent (LNCaP) cells overexpressing COBRA1 survive under androgen-deprivation conditions. Remarkably, treatment of prostate cancer cells with well-studied antitumorigenic agent, 2-methoxyestradiol (2-ME2), caused significant DNA methylation changes in 3255 genes including COBRA1. Furthermore, treatment of prostate cancer cells with 2-ME2 downregulates COBRA1 and inhibition of prostate tumors in TRAMP (transgenic adenocarcinomas of mouse prostate) animals with 2-ME2 was also associated with decreased COBRA1 levels. These observations implicate a novel role for COBRA1 in progression to CRPC and suggest that COBRA1 downregulation has therapeutic potential.

Single-Cell RNA-seq Reveals a Subpopulation of Prostate Cancer Cells with Enhanced Cell-Cycle-Related Transcription and Attenuated Androgen Response.

Increasing evidence suggests the presence of minor cell subpopulations in prostate cancer that are androgen independent and poised for selection as dominant clones after androgen deprivation therapy. In this study, we investigated this phenomenon by stratifying cell subpopulations based on transcriptome profiling of 144 single LNCaP prostate cancer cells treated or untreated with androgen after cell-cycle synchronization. Model-based clustering of 397 differentially expressed genes identified eight potential subpopulations of LNCaP cells, revealing a previously unappreciable level of cellular heterogeneity to androgen stimulation. One subpopulation displayed stem-like features with a slower cell doubling rate, increased sphere formation capability, and resistance to G2-M arrest induced by a mitosis inhibitor. Advanced growth of this subpopulation was associated with enhanced expression of 10 cell-cycle-related genes (CCNB2, DLGAP5, CENPF, CENPE, MKI67, PTTG1, CDC20, PLK1, HMMR, and CCNB1) and decreased dependence upon androgen receptor signaling. In silico analysis of RNA-seq data from The Cancer Genome Atlas further demonstrated that concordant upregulation of these genes was linked to recurrent prostate cancers. Analysis of receiver operating characteristic curves implicates aberrant expression of these genes and could be useful for early identification of tumors that subsequently develop biochemical recurrence. Moreover, this single-cell approach provides a better understanding of how prostate cancer cells respond heterogeneously to androgen deprivation therapies and reveals characteristics of subpopulations resistant to this treatment.Significance: Illustrating the challenge in treating cancers with targeted drugs, which by selecting for drug resistance can drive metastatic progression, this study characterized the plasticity and heterogeneity of prostate cancer cells with regard to androgen dependence, defining the character or minor subpopulations of androgen-independent cells that are poised for clonal selection after androgen-deprivation therapy. Cancer Res; 78(4); 853-64. ©2017 AACR.

Co-targeting ALK and EGFR parallel signaling in oral squamous cell carcinoma.

Squamous cell carcinoma (SCC) comprises 90% of all head and neck cancers and has a poor survival rate due to late-stage disease that is refractive to traditional therapies. Epidermal growth factor receptor (EGFR) is over-expressed in greater than 80% of head and neck SCC (HNSCC). However, EGFR targeted therapies yielded little to no efficacy in clinical trials. This study investigated the efficacy of co-targeting EGFR and the anaplastic lymphoma kinase (ALK) whose promoter is hypomethylated in late-stage oral SCC (OSCC). We observed increased ALK activity in late-stage human OSCC tumors and invasive OSCC cell lines. We also found that while ALK inhibition alone had little effect on proliferation, co-targeting ALK and EGFR significantly reduced OSCC cell proliferation in vitro. Further analysis showed significant efficacy of combined treatment in HSC3-derived xenografts resulting in a 30% decrease in tumor volumes by 14days (p<0.001). Western blot analysis showed that co-targeting ALK and EGFR significantly reduced EGFR phosphorylation (Y1148) in HSC3 cells but not Cal27 cells. ALK and EGFR downstream signaling interactions are also demonstrated by Western blot analysis in which lone EGFR and ALK inhibitors attenuated AKT activity whereas co-targeting ALK and EGFR completely abolished AKT activation. No effects were observed on ERK1/2 activation. STAT3 activity was significantly induced by lone ALK inhibition in HSC3 cells and to a lower extent in Cal27 cells. Together, these data illustrate that ALK inhibitors enhance anti-tumor activity of EGFR inhibitors in susceptible tumors that display increased ALK expression, most likely through abolition of AKT activation.

Integrative analysis identifies targetable CREB1/FoxA1 transcriptional co-regulation as a predictor of prostate cancer recurrence.

Identifying prostate cancer-driving transcription factors (TFs) in addition to the androgen receptor promises to improve our ability to effectively diagnose and treat this disease. We employed an integrative genomics analysis of master TFs CREB1 and FoxA1 in androgen-dependent prostate cancer (ADPC) and castration-resistant prostate cancer (CRPC) cell lines, primary prostate cancer tissues and circulating tumor cells (CTCs) to investigate their role in defining prostate cancer gene expression profiles. Combining genome-wide binding site and gene expression profiles we define CREB1 as a critical driver of pro-survival, cell cycle and metabolic transcription programs. We show that CREB1 and FoxA1 co-localize and mutually influence each other's binding to define disease-driving transcription profiles associated with advanced prostate cancer. Gene expression analysis in human prostate cancer samples found that CREB1/FoxA1 target gene panels predict prostate cancer recurrence. Finally, we showed that this signaling pathway is sensitive to compounds that inhibit the transcription co-regulatory factor MED1. These findings not only reveal a novel, global transcriptional co-regulatory function of CREB1 and FoxA1, but also suggest CREB1/FoxA1 signaling is a targetable driver of prostate cancer progression and serves as a biomarker of poor clinical outcomes.

DNA methylation screening of primary prostate tumors identifies SRD5A2 and CYP11A1 as candidate markers for assessing risk of biochemical recurrence.

BACKGROUND: Altered DNA methylation in CpG islands of gene promoters has been implicated in prostate cancer (PCa) progression and can be used to predict disease outcome. In this study, we determine whether methylation changes of androgen biosynthesis pathway (ABP)-related genes in patients' plasma cell-free DNA (cfDNA) can serve as prognostic markers for biochemical recurrence (BCR). METHODS: Methyl-binding domain capture sequencing (MBDCap-seq) was used to identify differentially methylated regions (DMRs) in primary tumors of patients who subsequently developed BCR or not, respectively. Methylation pyrosequencing of candidate loci was validated in cfDNA samples of 86 PCa patients taken at and/or post-radical prostatectomy (RP) using univariate and multivariate prediction analyses. RESULTS: Putative DMRs in 13 of 30 ABP-related genes were found between tumors of BCR (n = 12) versus no evidence of disease (NED) (n = 15). In silico analysis of The Cancer Genome Atlas data confirmed increased DNA methylation of two loci-SRD5A2 and CYP11A1, which also correlated with their decreased expression, in tumors with subsequent BCR development. Their aberrant cfDNA methylation was also associated with detectable levels of PSA taken after patients' post-RP. Multivariate analysis of the change in cfDNA methylation at all of CpG sites measured along with patient's treatment history predicted if a patient will develop BCR with 77.5% overall accuracy. CONCLUSIONS: Overall, increased DNA methylation of SRD5A2 and CYP11A1 related to androgen biosynthesis functions may play a role in BCR after patients' RP. The correlation between aberrant cfDNA methylation and detectable PSA in post-RP further suggests their utility as predictive markers for PCa recurrence. .

Nanomechanical biomarkers of single circulating tumor cells for detection of castration resistant prostate cancer.

BACKGROUND: Emerging evidence shows that nanomechanical phenotypes of circulating tumor cells (CTC) could become potential biomarkers for metastatic castration resistant prostate cancer (mCRPC). METHODS: To determine the nanomechanical phenotypes of CTCs we applied atomic force microscopy (AFM) employing the PeakForce quantitative nanomechanical (QNM) imaging. We assessed biophysical parameters (elasticity, deformation, and adhesion) of 130 CTCs isolated from blood samples from five castration sensitive (CS) and 12 castration resistant prostate cancer (CRPCa) patients. RESULTS: We found that CTCs from CRPCa patients are three times softer, three times more deformable, and seven times more adhesive than counterparts from CSPCa patients. Both nonsupervised hierarchical clustering and principle component analysis show that three combined nanomechanical parameters could constitute a valuable set to distinguish between CSPCa and CRPCa. CONCLUSIONS: [corrected] Our study indicates that nanomechanical phenotypes of CTCs may serve as novel and effective biomarkers for mCRPC.