The discordance between genome size and the complexity of eukaryotes can partly be attributed to differences in repeat density. The Muller F element (∼5.2 Mb) is the smallest chromosome in Drosophila melanogaster, but it is substantially larger (>18.7 Mb) in D. ananassae To identify the major contributors to the expansion of the F element and to assess their impact, we improved the genome sequence and annotated the genes in a 1.4-Mb region of the D. ananassae F element, and a 1.7-Mb region from the D element for comparison. We find that transposons (particularly LTR and LINE retrotransposons) are major contributors to this expansion (78.6%), while Wolbachia sequences integrated into the D. ananassae genome are minor contributors (0.02%). Both D. melanogaster and D. ananassae F-element genes exhibit distinct characteristics compared to D-element genes (e.g., larger coding spans, larger introns, more coding exons, and lower codon bias), but these differences are exaggerated in D. ananassae Compared to D. melanogaster, the codon bias observed in D. ananassae F-element genes can primarily be attributed to mutational biases instead of selection. The 5' ends of F-element genes in both species are enriched in dimethylation of lysine 4 on histone 3 (H3K4me2), while the coding spans are enriched in H3K9me2. Despite differences in repeat density and gene characteristics, D. ananassae F-element genes show a similar range of expression levels compared to genes in euchromatic domains. This study improves our understanding of how transposons can affect genome size and how genes can function within highly repetitive domains.

Constitutively active estrogen receptor α (ER/ESR1) mutations have been identified in approximately one-third of ER+ metastatic breast cancers. Although these mutations are known as mediators of endocrine resistance, their potential role in promoting metastatic disease has not yet been mechanistically addressed. In this study, we show the presence of ESR1 mutations exclusively in distant but not local recurrences in five independent breast cancer cohorts. In concordance with transcriptomic profiling of ESR1-mutant tumors, genome-edited ESR1 Y537S and D538G-mutant cell models exhibited a reprogrammed cell adhesive gene network via alterations in desmosome/gap junction genes and the TIMP3/MMP axis, which functionally conferred enhanced cell-cell contacts while decreasing cell-extracellular matrix adhesion. In vivo studies showed ESR1-mutant cells were associated with larger multicellular circulating tumor cell (CTC) clusters with increased compactness compared with ESR1 wild-type CTCs. These preclinical findings translated to clinical observations, where CTC clusters were enriched in patients with ESR1-mutated metastatic breast cancer. Conversely, context-dependent migratory phenotypes revealed cotargeting of Wnt and ER as a vulnerability in a D538G cell model. Mechanistically, mutant ESR1 exhibited noncanonical regulation of several metastatic pathways, including secondary transcriptional regulation and de novo FOXA1-driven chromatin remodeling. Collectively, these data provide evidence for ESR1 mutation-modulated metastasis and suggest future therapeutic strategies for targeting ESR1-mutant breast cancer.

Invasive lobular carcinoma (ILC) is a histological subtype of breast cancer with distinct molecular and clinical features from the more common subtype invasive ductal carcinoma (IDC). ILC cells exhibit anchorage-independent growth in ultra-low attachment (ULA) suspension cultures, which is largely attributed to the loss of E-cadherin. In addition to anoikis resistance, herein we show that human ILC cell lines exhibit enhanced cell proliferation in ULA cultures as compared to IDC cells. Proteomic comparison of ILC and IDC cell lines identified induction of PI3K/Akt and p90-RSK pathways specifically in ULA culture in ILC cells. Further transcriptional profiling uncovered unique upregulation of the inhibitors of differentiation family transcription factors ID1 and ID3 in ILC ULA culture, the knockdown of which diminished the anchorage-independent growth of ILC cell lines through cell cycle arrest. We find that ID1 and ID3 expression is higher in human ILC tumors as compared to IDC, correlated with worse prognosis uniquely in patients with ILC and associated with upregulation of angiogenesis and matrisome-related genes. Altogether, our comprehensive study of anchorage independence in human ILC cell lines provides mechanistic insights and clinical implications for metastatic dissemination of ILC and implicates ID1 and ID3 as novel drivers and therapeutic targets for lobular breast cancer.

Invasive lobular carcinoma (ILC) is the second most common histological subtype of breast cancer following invasive ductal carcinoma (IDC). ILC differs from IDC in a number of histological and clinical features, such as single strand growth, difficulty in detection, and frequent late recurrences. To understand the molecular pathways involved in the clinical characteristics of ILC, we compared the gene expression profiles of luminal A ILC and luminal A IDC using data from TCGA and utilized samples from METABRIC as a validation data set. Top pathways that were significantly enriched in ILC were related to immune response. ILC exhibited a higher activity of almost all types of immune cells based on cell type-specific signatures compared to IDC. Conversely, pathways that were less enriched in ILC were related to protein translation and metabolism, which we functionally validated in cell lines. The higher immune activity uncovered in our study highlights the currently unexplored potential of a response to immunotherapy in a subset of patients with ILC. Furthermore, the lower rates of protein translation and metabolism - known features of tumor dormancy - may play a role in the late recurrences of ILC and lower detection rate in mammography and PET scanning.

Invasive lobular breast carcinoma (ILC) accounts for 10% to 15% of breast cancers diagnosed annually. Evidence suggests that some aspects of endocrine treatment response might differ between invasive ductal carcinoma (IDC) and ILC, and that patients with ILC have worse long-term survival. We analyzed The Cancer Genome Atlas dataset and observed lower levels of ESR1 mRNA (P = 0.002) and ERα protein (P = 0.038) in ER+ ILC (n = 137) compared to IDC (n = 554), and further confirmed the mRNA difference in a local UPMC cohort (ILC, n = 143; IDC, n = 877; P < 0.005). In both datasets, the correlation between ESR1 mRNA and ERα protein was weaker in ILC, suggesting differential post-transcriptional regulation of ERα. In vitro, 17β-estradiol (E2) decreased the rate of degradation and increased the half-life of ERα in ILC cell lines, whereas the opposite was observed in IDC cell lines. Further, E2 failed to induce robust ubiquitination of ERα in ILC cells. To determine the potential clinical relevance of these findings, we evaluated the effect of 2 selective estrogen receptor downregulators (SERDs), ICI 182,780 and AZD9496, on ERα turnover and cell growth. While ICI 182,780 and AZD9496 showed similar effects in IDC cells, in ILC cell lines, AZD9496 was not as effective as ICI 182,780 in decreasing ERα stability and E2-induced proliferation. Furthermore, AZD9496 exhibited partial agonist activity in growth assays in ILC cell lines. Our study provides evidence for a distinct ERα regulation by SERDs in ILC cell lines, and therefore it is important to include ILC models into preclinical and clinical testing of novel SERDs.

Mutations in the estrogen receptor alpha (ERα) 1 gene (ESR1) are frequently detected in ER+ metastatic breast cancer, and there is increasing evidence that these mutations confer endocrine resistance in breast cancer patients with advanced disease. However, their functional role is not well-understood, at least in part due to a lack of ESR1 mutant models. Here, we describe the generation and characterization of genome-edited T47D and MCF7 breast cancer cell lines with the two most common ESR1 mutations, Y537S and D538G.

Eukaryotic chromatin is composed of DNA and protein components-core histones-that act to compactly pack the DNA into nucleosomes, the fundamental building blocks of chromatin. These nucleosomes are connected to adjacent nucleosomes by linker histones. Nucleosomes are highly dynamic and, through various core histone post-translational modifications and incorporation of diverse histone variants, can serve as epigenetic marks to control processes such as gene expression and recombination. The Histone Sequence Database is a curated collection of sequences and structures of histones and non-histone proteins containing histone folds, assembled from major public databases. Here, we report a substantial increase in the number of sequences and taxonomic coverage for histone and histone fold-containing proteins available in the database. Additionally, the database now contains an expanded dataset that includes archaeal histone sequences. The database also provides comprehensive multiple sequence alignments for each of the four core histones (H2A, H2B, H3 and H4), the linker histones (H1/H5) and the archaeal histones. The database also includes current information on solved histone fold-containing structures. The Histone Sequence Database is an inclusive resource for the analysis of chromatin structure and function focused on histones and histone fold-containing proteins.

Flavin-dependent histone demethylases govern histone H3K4 methylation and act as important chromatin modulators that are extensively involved in regulation of DNA replication, gene transcription, DNA repair, and heterochromatin gene silencing. While the activities of lysine-specific demethylase 1 (LSD1/KDM1A) in facilitating breast cancer progression have been well characterized, the roles of its homolog LSD2 (KDM1B) in breast oncogenesis are relatively less understood. In this study, we showed that LSD2 protein level was significantly elevated in malignant breast cell lines compared with normal breast epithelial cell line. TCGA- Oncomine database showed that LSD2 expression is significantly higher in basal-like breast tumors compared to other breast cancer subtypes or normal breast tissue. Overexpression of LSD2 in MDA-MB-231 cells significantly altered the expression of key important epigenetic modifiers such as LSD1, HDAC1/2, and DNMT3B; promoted cellular proliferation; and augmented colony formation in soft agar; while attenuating motility and invasion. Conversely, siRNA-mediated depletion of endogenous LSD2 hindered growth of multiple breast cancer cell lines while shRNA-mediated LSD2 depletion augmented motility and invasion. Moreover, LSD2 overexpression in MDA-MB-231 cells facilitated mammosphere formation, enriched the subpopulation of CD49f+/EpCAM- and ALDHhigh, and induced the expression of pluripotent stem cell markers, NANOG and SOX2. In xenograft studies using immune-compromised mice, LSD2-overexpressing MDA-MB-231 cells displayed accelerated tumor growth but significantly fewer lung metastases than controls. Taken together, our findings provide novel insights into the critical and multifaceted roles of LSD2 in the regulation of breast cancer progression and cancer stem cell enrichment.

Increased evidence suggests that somatic mutations in the ligand-binding domain of estrogen receptor [ER (ERα/ESR1)] are critical mediators of endocrine-resistant breast cancer progression. Insulinlike growth factor-1 (IGF1) is an essential regulator of breast development and tumorigenesis and also has a role in endocrine resistance. A recent study showed enhanced crosstalk between IGF1 and ERα in ESR1 mutant cells, but detailed mechanisms are incompletely understood. Using genome-edited MCF-7 and T47D cell lines harboring Y537S and D538G ESR1 mutations, we characterized altered IGF1 signaling. RNA sequencing revealed upregulation of multiple genes in the IGF1 pathway, including insulin receptor substrate-1 (IRS1), consistent in both Y537S and D538G ESR1 mutant cell line models. Higher IRS1 expression was confirmed by quantitative reverse transcription polymerase chain reaction and immunoblotting. ESR1 mutant cells also showed increased levels of IGF-regulated genes, reflected by activation of an IGF signature. IGF1 showed increased sensitivity and potency in growth stimulation of ESR1 mutant cells. Analysis of downstream signaling revealed the phosphoinositide 3-kinase (PI3K)-Akt axis as a major pathway mediating the enhanced IGF1 response in ESR1 mutant cells. Decreasing IRS1 expression by small interfering RNA diminished the increased sensitivity to IGF1. Combination treatment with inhibitors against IGF1 receptor (IGF1R; OSI-906) and ER (fulvestrant) showed synergistic growth inhibition in ESR1 mutant cells, particularly at lower effective concentrations. Our study supports a critical role of enhanced IGF1 signaling in ESR1 mutant cell lines, pointing toward a potential for cotargeting IGF1R and ERα in endocrine-resistant breast tumors with mutant ESR1.

Invasive lobular breast carcinoma (ILC) is a histological subtype of breast cancer that is characterized by loss of E-cadherin and high expression of estrogen receptor alpha (ERα). In many cases, ILC is effectively treated with adjuvant aromatase inhibitors (AIs); however, acquired AI resistance remains a significant problem.

Invasive lobular carcinoma (ILC) is the second most common histological subtype of breast cancer following invasive ductal carcinoma (IDC). To identify potential genetic drivers of ILC progression, we used NanoString nCounter technology to investigate the DNA copy number (CN) in 70 well-curated primary ILC samples. We confirmed prior observations of frequent amplification of CCND1 (33%), and MYC (17%) in ILC, but additionally identified a substantial subset of ILCs with ESR1 and ERBB2 (19%) amplifications. Of interest, tumors with ESR1 CN gains (14%) and amplification (10%) were more likely to recur compared to those with normal CN. Finally, we observed that MDM4 (MDMX) was amplified in 17% of ILC samples. MDM4 knockdown in TP53 wild-type ILC cell lines caused increased apoptosis, decreased proliferation associated with cell cycle arrest, and concomitant activation of TP53 target genes. Similar effects were seen in TP53 mutant cells, indicting a TP53-independent role for MDM4 in ILC. To conclude, amplification of ESR1 and MDM4 are potential genetic drivers of ILC. These amplifications may represent actionable, targetable tumor dependencies, and thus have potential clinical implications and warrant further study.

Invasive lobular breast carcinoma (ILC), one of the major breast cancer histologic subtypes, exhibits unique features compared with the well-studied ductal cancer subtype (IDC). The pathognomonic feature of ILC is loss of E-cadherin, mainly caused by inactivating mutations, but the contribution of this genetic alteration to ILC-specific molecular characteristics remains largely understudied. To profile these features transcriptionally, we conducted single-cell RNA sequencing on a panel of IDC and ILC cell lines, and an IDC cell line (T47D) with CRISPR-Cas9-mediated E-cadherin knockout (KO). Inspection of intracell line heterogeneity illustrated genetically and transcriptionally distinct subpopulations in multiple cell lines and highlighted rare populations of MCF7 cells highly expressing an apoptosis-related signature, positively correlated with a preadaptation signature to estrogen deprivation. Investigation of E-cadherin KO-induced alterations showed transcriptomic membranous systems remodeling, elevated resemblance to ILCs in regulon activation, and increased sensitivity to IFNγ-mediated growth inhibition via activation of IRF1. This study reveals single-cell transcriptional heterogeneity in breast cancer cell lines and provides a resource to identify drivers of cancer progression and drug resistance. SIGNIFICANCE: This study represents a key step towards understanding heterogeneity in cancer cell lines and the role of E-cadherin depletion in contributing to the molecular features of invasive lobular breast carcinoma.

Invasive lobular carcinoma (ILC) is an understudied subtype of breast cancer that requires novel therapies in the advanced setting. To study acquired resistance to endocrine therapy in ILC, we have recently performed RNA-Sequencing on long-term estrogen deprived cell lines and identified FGFR4 overexpression as a top druggable target. Here, we show that FGFR4 expression also increases dramatically in endocrine-treated distant metastases, with an average fold change of 4.8 relative to the paired primary breast tumor for ILC, and 2.4-fold for invasive ductal carcinoma (IDC). In addition, we now report that FGFR4 hotspot mutations are enriched in metastatic breast cancer, with an additional enrichment for ILC, suggesting a multimodal selection of FGFR4 activation. These data collectively support the notion that FGFR4 is an important mediator of endocrine resistance in ILC, warranting future mechanistic studies on downstream signaling of overexpressed wild-type and mutant FGFR4.