Germline genetic contribution to the immune landscape of cancer.

Understanding the contribution of the host's genetic background to cancer immunity may lead to improved stratification for immunotherapy and to the identification of novel therapeutic targets. We investigated the effect of common and rare germline variants on 139 well-defined immune traits in ∼9000 cancer patients enrolled in TCGA. High heritability was observed for estimates of NK cell and T cell subset infiltration and for interferon signaling. Common variants of IFIH1, TMEM173 (STING1), and TMEM108 were associated with differential interferon signaling and variants mapping to RBL1 correlated with T cell subset abundance. Pathogenic or likely pathogenic variants in BRCA1 and in genes involved in telomere stabilization and Wnt-ß-catenin also acted as immune modulators. Our findings provide evidence for the impact of germline genetics on the composition and functional orientation of the tumor immune microenvironment. The curated datasets, variants, and genes identified provide a resource toward further understanding of tumor-immune interactions.

Loss of epigenetic suppression of retrotransposons with oncogenic potential in aging mammary luminal epithelial cells.

A primary function of DNA methylation in mammalian genomes is to repress transposable elements (TEs). The widespread methylation loss that is commonly observed in cancer cells results in the loss of epigenetic repression of TEs. The aging process is similarly characterized by changes to the methylome. However, the impact of these epigenomic alterations on TE silencing and the functional consequences of this have remained unclear. To assess the epigenetic regulation of TEs in aging, we profiled DNA methylation in human mammary luminal epithelial cells (LEps)-a key cell lineage implicated in age-related breast cancers-from younger and older women. We report here that several TE subfamilies function as regulatory elements in normal LEps, and a subset of these display consistent methylation changes with age. Methylation changes at these TEs occurred at lineage-specific transcription factor binding sites, consistent with loss of lineage specificity. Whereas TEs mainly showed methylation loss, CpG islands (CGIs) that are targets of the Polycomb repressive complex 2 (PRC2) show a gain of methylation in aging cells. Many TEs with methylation loss in aging LEps have evidence of regulatory activity in breast cancer samples. We furthermore show that methylation changes at TEs impact the regulation of genes associated with luminal breast cancers. These results indicate that aging leads to DNA methylation changes at TEs that undermine the maintenance of lineage specificity, potentially increasing susceptibility to breast cancer.

Changes in Immune Cell Types with Age in Breast are Consistent with a Decline in Immune Surveillance and Increased Immunosuppression.

A majority of breast cancers (BC) are age-related and we seek to determine what cellular and molecular changes occur in breast tissue with age that make women more susceptible to cancer initiation. Immune-epithelial cell interactions are important during mammary gland development and the immune system plays an important role in BC progression. The composition of human immune cell populations is known to change in peripheral blood with age and in breast tissue during BC progression. Less is known about changes in immune populations in normal breast tissue and how their interactions with mammary epithelia change with age. We quantified densities of T cells, B cells, and macrophage subsets in pathologically normal breast tissue from 122 different women who ranged in age from 24 to 74 years old. Donor-matched peripheral blood from a subset of 20 donors was analyzed by flow cytometry. Tissue immune cell densities and localizations relative to the epithelium were quantified in situ with machine learning-based image analyses of multiplex immunohistochemistry-stained tissue sections. In situ results were corroborated with flow cytometry analyses of peri-epithelial immune cells from primary breast tissue preparations and transcriptome analyses of public data from bulk tissue reduction mammoplasties. Proportions of immune cell subsets in breast tissue and donor-matched peripheral blood were not correlated. Density (cells/mm2) of T and B lymphocytes in situ decreased with age. T cells and macrophages preferentially localized near or within epithelial bilayers, rather than the intralobular stroma. M2 macrophage density was higher than M1 macrophage density and this difference was due to higher density of M2 in the intralobular stroma. Transcriptional signature analyses suggested age-dependent decline in adaptive immune cell populations and functions and increased innate immune cell activity. T cells and macrophages are so intimately associated with the epithelia that they are embedded within the bilayer, suggesting an important role for immune-epithelial cell interactions. Age-associated decreased T cell density in peri-epithelial regions, and increased M2 macrophage density in intralobular stroma suggests the emergence of a tissue microenvironment that is simultaneously immune-senescent and immunosuppressive with age.

Discriminating external and internal causes for heading changes in freely flying Drosophila.

As animals move through the world in search of resources, they change course in reaction to both external sensory cues and internally-generated programs. Elucidating the functional logic of complex search algorithms is challenging because the observable actions of the animal cannot be unambiguously assigned to externally- or internally-triggered events. We present a technique that addresses this challenge by assessing quantitatively the contribution of external stimuli and internal processes. We apply this technique to the analysis of rapid turns ("saccades") of freely flying Drosophila melanogaster. We show that a single scalar feature computed from the visual stimulus experienced by the animal is sufficient to explain a majority (93%) of the turning decisions. We automatically estimate this scalar value from the observable trajectory, without any assumption regarding the sensory processing. A posteriori, we show that the estimated feature field is consistent with previous results measured in other experimental conditions. The remaining turning decisions, not explained by this feature of the visual input, may be attributed to a combination of deterministic processes based on unobservable internal states and purely stochastic behavior. We cannot distinguish these contributions using external observations alone, but we are able to provide a quantitative bound of their relative importance with respect to stimulus-triggered decisions. Our results suggest that comparatively few saccades in free-flying conditions are a result of an intrinsic spontaneous process, contrary to previous suggestions. We discuss how this technique could be generalized for use in other systems and employed as a tool for classifying effects into sensory, decision, and motor categories when used to analyze data from genetic behavioral screens.

Cell-free DNA Concentration as a Biomarker of Response and Recurrence in HER2-Negative Breast Cancer Receiving Neoadjuvant Chemotherapy.

PURPOSE: We previously demonstrated the clinical significance of circulating tumor DNA (ctDNA) in patients with HER2-negative breast cancer receiving neoadjuvant chemotherapy (NAC). Here, we compared its predictive and prognostic value with cell-free DNA (cfDNA) concentration measured in the same samples from the same patients. EXPERIMENTAL DESIGN: 145 patients with hormone receptor (HR)-positive/HER2-negative and 138 triple-negative breast cancer (TNBC) with ctDNA data from a previous study were included in the analysis. Associations of serial cfDNA concentration with residual cancer burden (RCB) and distant recurrence-free survival (DRFS) were examined. RESULTS: In TNBC, we observed a modest negative correlation between cfDNA concentration 3 weeks after treatment initiation and RCB, but none of the other timepoints showed significant correlation. In contrast, ctDNA was significantly positively correlated with RCB at all timepoints (all R > 0.3 and P < 0.05). In the HR-positive/HER2-negative group, cfDNA concentration did not associate with response to NAC, but survival analysis showed that high cfDNA shedders at pretreatment had a significantly worse DRFS than low shedders (hazard ratio, 2.12; P = 0.037). In TNBC, the difference in survival between high versus low cfDNA shedders at all timepoints was not statistically significant. In contrast, as previously reported, ctDNA at all timepoints was significantly correlated with DRFS in both subtypes. CONCLUSIONS: In TNBC, cfDNA concentrations during therapy were not strongly correlated with response or prognosis. In the HR-positive/HER2-negative group, pretreatment cfDNA concentration was prognostic for DRFS. Overall, the predictive and prognostic value of cfDNA concentration was more limited than that of ctDNA.

Clinical significance and biology of circulating tumor DNA in high-risk early-stage HER2-negative breast cancer receiving neoadjuvant chemotherapy.

Circulating tumor DNA (ctDNA) analysis may improve early-stage breast cancer treatment via non-invasive tumor burden assessment. To investigate subtype-specific differences in the clinical significance and biology of ctDNA shedding, we perform serial personalized ctDNA analysis in hormone receptor (HR)-positive/HER2-negative breast cancer and triple-negative breast cancer (TNBC) patients receiving neoadjuvant chemotherapy (NAC) in the I-SPY2 trial. ctDNA positivity rates before, during, and after NAC are higher in TNBC than in HR-positive/HER2-negative breast cancer patients. Early clearance of ctDNA 3 weeks after treatment initiation predicts a favorable response to NAC in TNBC only. Whereas ctDNA positivity associates with reduced distant recurrence-free survival in both subtypes. Conversely, ctDNA negativity after NAC correlates with improved outcomes, even in patients with extensive residual cancer. Pretreatment tumor mRNA profiling reveals associations between ctDNA shedding and cell cycle and immune-associated signaling. On the basis of these findings, the I-SPY2 trial will prospectively test ctDNA for utility in redirecting therapy to improve response and prognosis.

An integrated tumor, immune and microbiome atlas of colon cancer.

The lack of multi-omics cancer datasets with extensive follow-up information hinders the identification of accurate biomarkers of clinical outcome. In this cohort study, we performed comprehensive genomic analyses on fresh-frozen samples from 348 patients affected by primary colon cancer, encompassing RNA, whole-exome, deep T cell receptor and 16S bacterial rRNA gene sequencing on tumor and matched healthy colon tissue, complemented with tumor whole-genome sequencing for further microbiome characterization. A type 1 helper T cell, cytotoxic, gene expression signature, called Immunologic Constant of Rejection, captured the presence of clonally expanded, tumor-enriched T cell clones and outperformed conventional prognostic molecular biomarkers, such as the consensus molecular subtype and the microsatellite instability classifications. Quantification of genetic immunoediting, defined as a lower number of neoantigens than expected, further refined its prognostic value. We identified a microbiome signature, driven by Ruminococcus bromii, associated with a favorable outcome. By combining microbiome signature and Immunologic Constant of Rejection, we developed and validated a composite score (mICRoScore), which identifies a group of patients with excellent survival probability. The publicly available multi-omics dataset provides a resource for better understanding colon cancer biology that could facilitate the discovery of personalized therapeutic approaches.

Analysis of germline-driven ancestry-associated gene expression in cancers.

Differential mRNA expression between ancestry groups can be explained by both genetic and environmental factors. We outline a computational workflow to determine the extent to which germline genetic variation explains cancer-specific molecular differences across ancestry groups. Using multi-omics datasets from The Cancer Genome Atlas (TCGA), we enumerate ancestry-informative markers colocalized with cancer-type-specific expression quantitative trait loci (e-QTLs) at ancestry-associated genes. This approach is generalizable to other settings with paired germline genotyping and mRNA expression data for a multi-ethnic cohort. For complete details on the use and execution of this protocol, please refer to Carrot-Zhang et al. (2020), Robertson et al. (2021), and Sayaman et al. (2021).

Analytic pipelines to assess the relationship between immune response and germline genetics in human tumors.

Germline genetic variants modulate human immune response. We present analytical pipelines for assessing the contribution of hosts' genetic background to the immune landscape of solid tumors using harmonized data from more than 9,000 patients in The Cancer Genome Atlas (TCGA). These include protocols for heritability, genome-wide association studies (GWAS), colocalization, and rare variant analyses. These workflows are developed around the structure of TCGA but can be adapted to explore other repositories or in the context of cancer immunotherapy. For complete details on the use and execution of this protocol, please refer to Sayaman et al. (2021).

Redefining breast cancer subtypes to guide treatment prioritization and maximize response: Predictive biomarkers across 10 cancer therapies.

Using pre-treatment gene expression, protein/phosphoprotein, and clinical data from the I-SPY2 neoadjuvant platform trial (NCT01042379), we create alternative breast cancer subtypes incorporating tumor biology beyond clinical hormone receptor (HR) and human epidermal growth factor receptor-2 (HER2) status to better predict drug responses. We assess the predictive performance of mechanism-of-action biomarkers from ∼990 patients treated with 10 regimens targeting diverse biology. We explore >11 subtyping schemas and identify treatment-subtype pairs maximizing the pathologic complete response (pCR) rate over the population. The best performing schemas incorporate Immune, DNA repair, and HER2/Luminal phenotypes. Subsequent treatment allocation increases the overall pCR rate to 63% from 51% using HR/HER2-based treatment selection. pCR gains from reclassification and improved patient selection are highest in HR+ subsets (>15%). As new treatments are introduced, the subtyping schema determines the minimum response needed to show efficacy. This data platform provides an unprecedented resource and supports the usage of response-based subtypes to guide future treatment prioritization.

Volume-constrained microcontainers enable myoepithelial functional differentiation in highly parallel mammary organoid culture.

A long-standing constraint on organoid culture is the need to add exogenous substances to provide hydrogel matrix, which limits the study of fully human or fully native organoids. This paper introduces an approach to culture reconstituted mammary organoids without the impediment of exogenous matrix. We enclose organoids in nanoliter-scale, topologically enclosed, fluid compartments surrounded by agar. Organoids cultured in these "microcontainers" appear to secrete enough extracellular matrix to yield a self-sufficient microenvironment without exogenous supplements. In microcontainers, mammary organoids exhibit contractility and a high-level, physiological, myoepithelial (MEP) behavior that has not been previously reported in reconstituted organoids. The presence of contractility suggests that microcontainers elicit MEP functional differentiation, an important milestone. Microcontainers yield thousands of substantially identical and individually trackable organoids within a single culture vessel, enabling longitudinal studies and statistically powerful experiments, such as the evaluation of small effect sizes. Microcontainers open new doors for researchers who rely on organoid models.

Evidence for accelerated aging in mammary epithelia of women carrying germline BRCA1 or BRCA2 mutations.

During aging in the human mammary gland, luminal epithelial cells lose lineage fidelity by expressing markers normally expressed in myoepithelial cells. We hypothesize that loss of lineage fidelity is a general manifestation of epithelia that are susceptible to cancer initiation. In the present study, we show that histologically normal breast tissue from younger women who are susceptible to breast cancer, as a result of harboring a germline mutation in BRCA1, BRCA2 or PALB2 genes, exhibits hallmarks of accelerated aging. These include proportionately increased luminal epithelial cells that acquired myoepithelial markers, decreased proportions of myoepithelial cells and a basal differentiation bias or failure of differentiation of cKit+ progenitors. High-risk luminal and myoepithelial cells are transcriptionally enriched for genes of the opposite lineage, inflammatory- and cancer-related pathways. We have identified breast-aging hallmarks that reflect a convergent biology of cancer susceptibility, regardless of the specific underlying genetic or age-dependent risk or the associated breast cancer subtype.

Deep proteome profiling of human mammary epithelia at lineage and age resolution.

Age is the major risk factor in most carcinomas, yet little is known about how proteomes change with age in any human epithelium. We present comprehensive proteomes comprised of >9,000 total proteins and >15,000 phosphopeptides from normal primary human mammary epithelia at lineage resolution from ten women ranging in age from 19 to 68 years. Data were quality controlled and results were biologically validated with cell-based assays. Age-dependent protein signatures were identified using differential expression analyses and weighted protein co-expression network analyses. Upregulation of basal markers in luminal cells, including KRT14 and AXL, were a prominent consequence of aging. PEAK1 was identified as an age-dependent signaling kinase in luminal cells, which revealed a potential age-dependent vulnerability for targeted ablation. Correlation analyses between transcriptome and proteome revealed age-associated loss of proteostasis regulation. Age-dependent proteome changes in the breast epithelium identified heretofore unknown potential therapeutic targets for reducing breast cancer susceptibility.

Ancestry-associated transcriptomic profiles of breast cancer in patients of African, Arab, and European ancestry.

Breast cancer largely dominates the global cancer burden statistics; however, there are striking disparities in mortality rates across countries. While socioeconomic factors contribute to population-based differences in mortality, they do not fully explain disparity among women of African ancestry (AA) and Arab ancestry (ArA) compared to women of European ancestry (EA). In this study, we sought to identify molecular differences that could provide insight into the biology of ancestry-associated disparities in clinical outcomes. We applied a unique approach that combines the use of curated survival data from The Cancer Genome Atlas (TCGA) Pan-Cancer clinical data resource, improved single-nucleotide polymorphism-based inferred ancestry assignment, and a novel breast cancer subtype classification to interrogate the TCGA and a local Arab breast cancer dataset. We observed an enrichment of BasalMyo tumors in AA patients (38 vs 16.5% in EA, p = 1.30E - 10), associated with a significant worse overall (hazard ratio (HR) = 2.39, p = 0.02) and disease-specific survival (HR = 2.57, p = 0.03). Gene set enrichment analysis of BasalMyo AA and EA samples revealed differences in the abundance of T-regulatory and T-helper type 2 cells, and enrichment of cancer-related pathways with prognostic implications (AA: PI3K-Akt-mTOR and ErbB signaling; EA: EGF, estrogen-dependent and DNA repair signaling). Strikingly, AMPK signaling was associated with opposing prognostic connotation (AA: 10-year HR = 2.79, EA: 10-year HR = 0.34). Analysis of ArA patients suggests enrichment of BasalMyo tumors with a trend for differential enrichment of T-regulatory cells and AMPK signaling. Together, our findings suggest that the disparity in the clinical outcome of AA breast cancer patients is likely related to differences in cancer-related and microenvironmental features.

Breast-Specific Molecular Clocks Comprised of ELF5 Expression and Promoter Methylation Identify Individuals Susceptible to Cancer Initiation.

A robust breast cancer prevention strategy requires risk assessment biomarkers for early detection. We show that expression of ELF5, a transcription factor critical for normal mammary development, is downregulated in mammary luminal epithelia with age. DNA methylation of the ELF5 promoter is negatively correlated with expression in an age-dependent manner. Both ELF5 methylation and gene expression were used to build biological clocks to estimate chronological ages of mammary epithelia. ELF5 clock-based estimates of biological age in luminal epithelia from average-risk women were within three years of chronological age. Biological ages of breast epithelia from BRCA1 or BRCA2 mutation carriers, who were high risk for developing breast cancer, suggested they were accelerated by two decades relative to chronological age. The ELF5 DNA methylation clock had better performance at predicting biological age in luminal epithelial cells as compared with two other epigenetic clocks based on whole tissues. We propose that the changes in ELF5 expression or ELF5-proximal DNA methylation in luminal epithelia are emergent properties of at-risk breast tissue and constitute breast-specific biological clocks. PREVENTION RELEVANCE: ELF5 expression or DNA methylation level at the ELF5 promoter region can be used as breast-specific biological clocks to identify women at higher than average risk of breast cancer.

Comprehensive Analysis of Genetic Ancestry and Its Molecular Correlates in Cancer.

We evaluated ancestry effects on mutation rates, DNA methylation, and mRNA and miRNA expression among 10,678 patients across 33 cancer types from The Cancer Genome Atlas. We demonstrated that cancer subtypes and ancestry-related technical artifacts are important confounders that have been insufficiently accounted for. Once accounted for, ancestry-associated differences spanned all molecular features and hundreds of genes. Biologically significant differences were usually tissue specific but not specific to cancer. However, admixture and pathway analyses suggested some of these differences are causally related to cancer. Specific findings included increased FBXW7 mutations in patients of African origin, decreased VHL and PBRM1 mutations in renal cancer patients of African origin, and decreased immune activity in bladder cancer patients of East Asian origin.

Genetic Variation in CCL5 Signaling Genes and Triple Negative Breast Cancer: Susceptibility and Prognosis Implications.

Triple-negative breast cancer (TNBC) accounts for ~15-20% of breast cancer (BC) and has a higher rate of early relapse and mortality compared to other subtypes. The Chemokine (C-C motif) ligand 5 (CCL5) and its signaling pathway have been linked to TNBC. We aimed to investigate the susceptibility and prognostic implications of genetic variation in CCL5 signaling genes in TNBC in the present study. We characterized variants in CCL5 and that of six other CCL5 signaling genes (CCND1, ZMIZ1, CASP8, NOTCH2, MAP3K21, and HS6ST3) among 1,082 unrelated Tunisian subjects (544 BC patients, including 196 TNBC, and 538 healthy controls), assessed the association of the variants with BC-specific overall survival (OVS) and progression-free survival (PFS), and correlated CCL5 mRNA and serum levels with CCL5 genotypes. We found a highly significant association between the CCND1 rs614367-TT genotype (OR = 5.14; P = 0.004) and TNBC risk, and identified a significant association between the rs614367-T allele and decreased PFS in TNBC. A decreased risk of lymph node metastasis was associated with the MAP3K21 rs1294255-C allele, particularly in rs1294255-GC (OR = 0.47; P = 0.001). CCL5 variants (rs2107538 and rs2280789) were linked to CCL5 serum and mRNA levels. In the TCGA TNBC/Basal-like cohort the MAP3K21 rs1294255-G allele was associated with a decreased OVS. High expression of CCL5 in breast tumors was significantly associated with an increased OVS in all BC patients, but particularly in TNBC/Basal-like patients. In conclusion, genetic variation in CCL5 signaling genes may predict not only TNBC risk but also disease aggressiveness.

Tissue aging: the integration of collective and variant responses of cells to entropic forces over time.

Aging is driven by unavoidable entropic forces, physicochemical in nature, that damage the raw materials that constitute biological systems. Single cells experience and respond to stochastic physicochemical insults that occur either to the cells themselves or to their microenvironment, in a dynamic and reciprocal manner, leading to increased age-related cell-to-cell variation. We will discuss the biological mechanisms that integrate cell-to-cell variation across tissues resulting in stereotypical phenotypes of age.

Age-related gene expression in luminal epithelial cells is driven by a microenvironment made from myoepithelial cells.

Luminal epithelial cells in the breast gradually alter gene and protein expression with age, appearing to lose lineage-specificity by acquiring myoepithelial-like characteristics. We hypothesize that the luminal lineage is particularly sensitive to microenvironment changes, and age-related microenvironment changes cause altered luminal cell phenotypes. To evaluate the effects of different microenvironments on the fidelity of epigenetically regulated luminal and myoepithelial gene expression, we generated a set of lineage-specific probes for genes that are controlled through DNA methylation. Culturing primary luminal cells under conditions that favor myoepithelial propogation led to their reprogramming at the level of gene methylation, and to a more myoepithelial-like expression profile. Primary luminal cells' lineage-specific gene expression could be maintained when they were cultured as bilayers with primary myoepithelial cells. Isogenic stromal fibroblast co-cultures were unable to maintain the luminal phenotype. Mixed-age luminal-myoepithelial bilayers revealed that luminal cells adopt transcription and methylation patterns consistent with the chronological age of the myoepithelial cells. We provide evidence that the luminal epithelial phenotype is exquisitely sensitive to microenvironment conditions, and that states of aging are cell non-autonomously communicated through microenvironment cues over at least one cell diameter.

The aerodynamics of hovering flight in Drosophila.

Using 3D infrared high-speed video, we captured the continuous wing and body kinematics of free-flying fruit flies, Drosophila melanogaster, during hovering and slow forward flight. We then 'replayed' the wing kinematics on a dynamically scaled robotic model to measure the aerodynamic forces produced by the wings. Hovering animals generate a U-shaped wing trajectory, in which large drag forces during a downward plunge at the start of each stroke create peak vertical forces. Quasi-steady mechanisms could account for nearly all of the mean measured force required to hover, although temporal discrepancies between instantaneous measured forces and model predictions indicate that unsteady mechanisms also play a significant role. We analyzed the requirements for hovering from an analysis of the time history of forces and moments in all six degrees of freedom. The wing kinematics necessary to generate sufficient lift are highly constrained by the requirement to balance thrust and pitch torque over the stroke cycle. We also compare the wing motion and aerodynamic forces of free and tethered flies. Tethering causes a strong distortion of the stroke pattern that results in a reduction of translational forces and a prominent nose-down pitch moment. The stereotyped distortion under tethered conditions is most likely due to a disruption of sensory feedback. Finally, we calculated flight power based directly on the measurements of wing motion and aerodynamic forces, which yielded a higher estimate of muscle power during free hovering flight than prior estimates based on time-averaged parameters. This discrepancy is mostly due to a two- to threefold underestimate of the mean profile drag coefficient in prior studies. We also compared our values with the predictions of the same time-averaged models using more accurate kinematic and aerodynamic input parameters based on our high-speed videography measurements. In this case, the time-averaged models tended to overestimate flight costs.