Using the ACMG/AMP framework to capture evidence related to predicted and observed impact on splicing: Recommendations from the ClinGen SVI Splicing Subgroup

The American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1, PS3, PP3, BS3, BP4, and BP7. However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. We utilized empirically derived splicing evidence to (1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, (2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and (3) exemplify methodology to calibrate splice prediction tools. We propose repurposing the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely, BP7 may be used to capture RNA results demonstrating no splicing impact for intronic and synonymous variants. We propose that the PS3/BS3 codes are applied only for well-established assays that measure functional impact not directly captured by RNA-splicing assays. We recommend the application of PS1 based on similarity of predicted RNA-splicing effects for a variant under assessment in comparison with a known pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA-assay evidence described aim to help standardize variant pathogenicity classification processes when interpreting splicing-based evidence.

Copy number variants as modifiers of breast cancer risk for BRCA1/BRCA2 pathogenic variant carriers

The contribution of germline copy number variants (CNVs) to risk of developing cancer in individuals with pathogenic BRCA1 or BRCA2 variants remains relatively unknown. We conducted the largest genome-wide analysis of CNVs in 15,342 BRCA1 and 10,740 BRCA2 pathogenic variant carriers. We used these results to prioritise a candidate breast cancer risk-modifier gene for laboratory analysis and biological validation. Notably, the HR for deletions in BRCA1 suggested an elevated breast cancer risk estimate (hazard ratio (HR) = 1.21), 95% confidence interval (95% CI = 1.09-1.35) compared with non-CNV pathogenic variants. In contrast, deletions overlapping SULT1A1 suggested a decreased breast cancer risk (HR = 0.73, 95% CI 0.59-0.91) in BRCA1 pathogenic variant carriers. Functional analyses of SULT1A1 showed that reduced mRNA expression in pathogenic BRCA1 variant cells was associated with reduced cellular proliferation and reduced DNA damage after treatment with DNA damaging agents. These data provide evidence that deleterious variants in BRCA1 plus SULT1A1 deletions contribute to variable breast cancer risk in BRCA1 carriers.

P12.03.B Ascorbate alters the hypoxic pathway in glioblastoma cellsin vitro and associates with improved patient survival

Background

Glioblastomas are highly aggressive and hypoxic tumours. This environment activates the hypoxic pathway, driving glioma progression and treatment resistance. The hypoxic pathway is regulated by the hypoxia inducible factor (HIF) hydroxylases, which require oxygen as a substrate. Under normoxic conditions, the HIF hydroxylases are active, causing degradation and inhibition of HIF transcription factors. Under hypoxia, the activity of the hydroxylases reduces and HIF accumulates, activating the hypoxic response. HIF hydroxylases also require ascorbate as a cofactor for optimal function. The brain has one of the highest ascorbate levels in the human body, yet data on ascorbate levels in gliomas is scarce. Cellular ascorbate uptake occurs through solute carrier family 23 member 2 (SLC23A2). My aim is to understand the relationship between ascorbate, SLC23A2 and the hypoxic pathway in brain cancer using both in vitro cell culture and clinical samples.

Material and Methods

Ascorbate uptake was measured in human glioblastoma cell lines (T98G, U251MG, U87MG; ATCC) using reverse phase high performance liquid chromatography (HPLC-ECD). CRISPR-Cas was designed to knock-out SLC23A2. Clinical glioblastoma samples (n=37) and follow-up data were provided by the Cancer Society Tissue Bank and University of Otago Dunedin. Ethics and informed consent were obtained (H19/163, MEC/08/02/016). Ascorbate levels, measured by HPLC-ECD, and HIF-1α and downstream targets were measured using Western blotting or ELISA. A HIF score was calculated from HIF-1α and downstream target protein levels to estimate hypoxic pathway activity.

Results

In this study we have shown that T98G and U251 cells accumulate up to 15 nmol ascorbate/106 cells when exposed to 500 µM ascorbate for up to 24 hours, compared to U87MG cells with up to 3 nmol ascorbate/106 cells. Cancer Cell Line Encyclopaedia data shows that T98G and U251MG cells express higher levels of SLC23A2 compared to U87MG cells, aligning with our results. Clinical glioblastoma tissue contained a median of 7.6 µg ascorbate/100 mg tissue. Patients survival was significantly longer with above, vs below, median tumour ascorbate levels (Gehan-Breslow-Wilcoxon p = 0.027). The HIF score was negatively correlated with tumour ascorbate levels (Pearson r -0.327, p = 0.048). Patients with higher HIF-score had significantly shorter survival time compared to those with a lower HIF score (Gehan-Breslow-Wilcoxon p = 0.005).

Conclusion

Ascorbate uptake in glioblastoma cells varies between cell lines and appears reliant on the level of SLC23A2. Higher ascorbate content in clinical glioblastoma samples was associated with reduced hypoxic pathway activity and longer patient survival. Ongoing work, using SLC23A2 CRISPR-Cas knock-out cells, is investigating the effect of disrupting ascorbate uptake on hypoxic pathway signalling in glioblastoma cells.

P12.06.A Relationship between ascorbate and DNA methylation markers in clinical glioma tumours

Background

Members of the 2-oxoglutarate-dependent dioxygenase (OGDD) enzyme family play an important role in gliomas as they regulate epigenetic modifications and response to hypoxia. The OGDDs require 2-OG and O2 as substrates, and ferrous iron and ascorbate as cofactors. Both hypoxia and aberrant DNA methylation are prognostic indicators for gliomas. The ten-eleven translocase (TET) DNA demethylases are OGDDs that convert 5-methyl cytosine (5mC) to 5-hydroxymethylcytosine (5hmC), with 5hmC levels related to better prognosis. Despite this, there is limited data on the OGDD enzymes and their substrates/cofactors in glioma tissues. Our previous study showed an association between ascorbate content and markers of the hypoxic response in glioblastoma tissue. Here we determine whether there is an association between ascorbate and DNA methylation in glioma. In addition, we assess whether methylation of the methylguanine-DNA methyltransferase (DNA repair enzyme MGMT) promoter is associated with ascorbate content.

Materials and methods

Frozen clinical glioma samples from 37 patients (n=11 WHO grade I-III, n=26 glioblastoma) were obtained from the Cancer Society Tissue Bank (Ethics approval H19/163). Isocitrate dehydrogenase 1 (IDH1) mutation status was determined by sequencing. Samples were processed on dry ice in liquid nitrogen and analysed for ascorbate (high-performance liquid chromatography), global DNA methylation (mass spectrometry) and MGMT promoter analyses (methylation specific PCR).

Results

Many grade I-III tumours were IDH1 R132H mutant (6/11), and most glioblastomas were not (2/26). Glioblastoma had significantly lower ascorbate content than grade I-III tumours (p=0.026). Glioblastoma also had lower global 5hmC levels (p=0.0013). IDH1 R132H tumours tended to have a lower ascorbate content (p=0.09). Ascorbate and 5hmC levels were directly correlated (Spearman r= 0.466, p=0.004). However, cytosine and 5mC showed no association with grade or ascorbate. MGMT promoter methylation status was not associated with global methylation or ascorbate content (p=0.97, p=0.96, respectively).

Conclusion

Our data suggests that ascorbate supports TET activity in clinical glioma. It also appears that site-specific (promoter) methylation was not affected by ascorbate availability. These findings may have clinical implications, as higher 5hmC levels are associated with improved outcome, whilst continued MGMT suppression suggests chemotherapy responsiveness. However, evidence that raising tumour ascorbate leads to increased 5hmC levels, or an associated improvement in survival, requires intervention trials.

Increased Ascorbate Content of Glioblastoma Is Associated With a Suppressed Hypoxic Response and Improved Patient Survival

Glioblastoma multiforme is a challenging disease with limited treatment options and poor survival. Glioblastoma tumours are characterised by hypoxia that activates the hypoxia inducible factor (HIF) pathway and controls a myriad of genes that drive cancer progression. HIF transcription factors are regulated at the post-translation level via HIF-hydroxylases. These hydroxylases require oxygen and 2-oxoglutarate as substrates, and ferrous iron and ascorbate as cofactors. In this retrospective observational study, we aimed to determine whether ascorbate played a role in the hypoxic response of glioblastoma, and whether this affected patient outcome. We measured the ascorbate content and members of the HIF-pathway of clinical glioblastoma samples, and assessed their association with clinicopathological features and patient survival. In 37 samples (37 patients), median ascorbate content was 7.6 μg ascorbate/100 mg tissue, range 0.8 – 20.4 μg ascorbate/100 mg tissue. In tumours with above median ascorbate content, HIF-pathway activity as a whole was significantly suppressed (p = 0.005), and several members of the pathway showed decreased expression (carbonic anhydrase-9 and glucose transporter-1, both p < 0.01). Patients with either lower tumour HIF-pathway activity or higher tumour ascorbate content survived significantly longer than patients with higher HIF-pathway or lower ascorbate levels (p = 0.011, p = 0.043, respectively). Median survival for the low HIF-pathway score group was 362 days compared to 203 days for the high HIF-pathway score group, and median survival for the above median ascorbate group was 390 days, compared to the below median ascorbate group with 219 days. The apparent survival advantage associated with higher tumour ascorbate was more prominent for the first 8 months following surgery. These associations are promising, suggesting an important role for ascorbate-regulated HIF-pathway activity in glioblastoma that may impact on patient survival.

Increased gene expression variability in BRCA1-associated and basal-like breast tumours

Purpose

Inherited variants in the cancer susceptibility genes, BRCA1 and BRCA2 account for up to 5% of breast cancers. Multiple gene expression studies have analysed gene expression patterns that maybe associated with BRCA12 pathogenic variant status; however, results from these studies lack consensus. These studies have focused on the differences in population means to identified genes associated with BRCA1/2-carriers with little consideration for gene expression variability, which is also under genetic control and is a feature of cellular function.

Methods

We measured differential gene expression variability in three of the largest familial breast cancer datasets and a 2116 breast cancer meta-cohort. Additionally, we used RNA in situ hybridisation to confirm expression variability of EN1 in an independent cohort of more than 500 breast tumours.

Results

BRCA1-associated breast tumours exhibited a 22.8% (95% CI 22.3–23.2) increase in transcriptome-wide gene expression variability compared to BRCAx tumours. Additionally, 40 genes were associated with BRCA1-related breast cancers that had ChIP-seq data suggestive of enriched EZH2 binding. Of these, two genes (EN1 and IGF2BP3) were significantly variable in both BRCA1-associated and basal-like breast tumours. RNA in situ analysis of EN1 supported a significant (p = 6.3 × 10⁻⁰⁴) increase in expression variability in BRCA1-associated breast tumours.

Conclusion

Our novel results describe a state of increased gene expression variability in BRCA1-related and basal-like breast tumours. Furthermore, genes with increased variability may be driven by changes in DNA occupancy of epigenetic effectors. The variation in gene expression is replicable and led to the identification of novel associations between genes and disease phenotypes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10549-021-06328-y.

Differential Expression of BARD1 Isoforms in Melanoma

Melanoma comprises <5% of cutaneous malignancies, yet it causes a significant proportion of skin cancer-related deaths worldwide. While new therapies for melanoma have been developed, not all patients respond well. Thus, further research is required to better predict patient outcomes. Using long-range nanopore sequencing, RT-qPCR, and RNA sequencing analyses, we examined the transcription of BARD1 splice isoforms in melanoma cell lines and patient tissue samples. Seventy-six BARD1 mRNA variants were identified in total, with several previously characterised isoforms (γ, φ, δ, ε, and η) contributing to a large proportion of the expressed transcripts. In addition, we identified four novel splice events, namely, Δ(E3_E9), ▼(i8), IVS10+131▼46, and IVS10▼176, occurring in various combinations in multiple transcripts. We found that short-read RNA-Seq analyses were limited in their ability to predict isoforms containing multiple non-contiguous splicing events, as compared to long-range nanopore sequencing. These studies suggest that further investigations into the functional significance of the identified BARD1 splice variants in melanoma are warranted.

Comprehensive Assessment of BARD1 Messenger Ribonucleic Acid Splicing With Implications for Variant Classification

Introduction: Case-control analyses have shown BARD1 variants to be associated with up to >2-fold increase in risk of breast cancer, and potentially greater risk of triple negative breast cancer. BARD1 is included in several gene sequencing panels currently marketed for the prediction of risk of cancer, however there are no gene-specific guidelines for the classification of BARD1 variants. We present the most comprehensive assessment of BARD1 messenger RNA splicing, and demonstrate the application of these data for the classification of truncating and splice site variants according to American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) guidelines. Methods: Nanopore sequencing, short-read RNA-seq (whole transcriptome and targeted), and capillary electrophoresis analysis were performed by four laboratories to investigate alternative BARD1 splicing in blood, breast, and fimbriae/ovary related specimens from non-cancer affected tissues. Splicing data were also collated from published studies of nine different tissues. The impact of the findings for PVS1 annotation was assessed for truncating and splice site variants. Results: We identified 62 naturally occurring alternative spliced BARD1 splicing events, including 19 novel events found by next generation sequencing and/or reverse transcription PCR analysis performed for this study. Quantitative analysis showed that naturally occurring splicing events causing loss of clinically relevant domains or nonsense mediated decay can constitute up to 11.9% of overlapping natural junctions, suggesting that aberrant splicing can be tolerated up to this level. Nanopore sequencing of whole BARD1 transcripts characterized 16 alternative isoforms from healthy controls, revealing that the most complex transcripts combined only two alternative splicing events. Bioinformatic analysis of ClinVar submitted variants at or near BARD1 splice sites suggest that all consensus splice site variants in BARD1 should be considered likely pathogenic, with the possible exception of variants at the donor site of exon 5. Conclusions: No BARD1 candidate rescue transcripts were identified in this study, indicating that all premature translation-termination codons variants can be annotated as PVS1. Furthermore, our analysis suggests that all donor and acceptor (IVS+/-1,2) variants can be considered PVS1 or PVS1_strong, with the exception of variants targeting the exon 5 donor site, that we recommend considering as PVS1_moderate.

Correction: Evaluation of copy-number variants as modifiers of breast and ovarian cancer risk for BRCA1 pathogenic variant carriers

This Article was originally published under a CC BY-NC-SA 4.0 license, but has now been made available under a CC BY 4.0 license. The PDF and HTML versions of the Article have been modified accordingly.

Nanopore sequencing of full-length BRCA1 mRNA transcripts reveals co-occurrence of known exon skipping events

Background

Laboratory assays evaluating the effect of DNA sequence variants on BRCA1 mRNA splicing may contribute to classification by providing molecular evidence. However, our knowledge of normal and aberrant BRCA1 splicing events to date has been limited to data derived from assays targeting partial transcript sequences. This study explored the utility of nanopore sequencing to examine whole BRCA1 mRNA transcripts and to provide accurate categorisation of in-frame and out-of-frame splicing events.

Methods

The exon structure of BRCA1 transcripts from a previously studied control lymphoblastoid cell line were assessed using MinION nanopore sequencing of long-range reverse transcriptase-PCR amplicons.

Results

Our study identified and characterised 32 complete BRCA1 isoforms, including 18 novel isoforms which showed skipping of multiple contiguous and/or non-contiguous exons. Furthermore, we show that known BRCA1 exon skipping events, such as Δ(9,10) and Δ21, can co-occur in a single transcript, with some isoforms containing four or more alternative splice junctions. Fourteen novel isoforms were formed entirely from a combination of previously identified alternative splice junctions, suggesting that the total number of BRCA1 isoforms might be lower than the number of splicing events reported previously.

Conclusions

Our results highlight complexity in BRCA1 transcript structure that has not been described previously. This finding has key implications for predicting the translation frame of splicing transcripts, important for interpreting the clinical significance of spliceogenic variants. Future research is warranted to quantitatively assess full-length BRCA1 transcript levels, and to assess the application of nanopore sequencing for routine evaluation of potential spliceogenic variants.

Electronic supplementary material

The online version of this article (doi:10.1186/s13058-017-0919-1) contains supplementary material, which is available to authorized users.

Evaluation of copy-number variants as modifiers of breast and ovarian cancer risk for BRCA1 pathogenic variant carriers

Genome-wide studies of patients carrying pathogenic variants (mutations) in BRCA1 or BRCA2 have reported strong associations between single-nucleotide polymorphisms (SNPs) and cancer risk. To conduct the first genome-wide association analysis of copy-number variants (CNVs) with breast or ovarian cancer risk in a cohort of 2500 BRCA1 pathogenic variant carriers, CNV discovery was performed using multiple calling algorithms and Illumina 610k SNP array data from a previously published genome-wide association study. Our analysis, which focused on functionally disruptive genomic deletions overlapping gene regions, identified a number of loci associated with risk of breast or ovarian cancer for BRCA1 pathogenic variant carriers. Despite only including putative deletions called by at least two or more algorithms, detection of selected CNVs by ancillary molecular technologies only confirmed 40% of predicted common (>1% allele frequency) variants. These include four loci that were associated (unadjusted P<0.05) with breast cancer (GTF2H2, ZNF385B, NAALADL2 and PSG5), and two loci associated with ovarian cancer (CYP2A7 and OR2A1). An interesting finding from this study was an association of a validated CNV deletion at the CYP2A7 locus (19q13.2) with decreased ovarian cancer risk (relative risk=0.50, P=0.007). Genomic analysis found this deletion coincides with a region displaying strong regulatory potential in ovarian tissue, but not in breast epithelial cells. This study highlighted the need to verify CNVs in vitro, but also provides evidence that experimentally validated CNVs (with plausible biological consequences) can modify risk of breast or ovarian cancer in BRCA1 pathogenic variant carriers.

Increased genomic burden of germline copy number variants is associated with early onset breast cancer: Australian breast cancer family registry

BACKGROUND

Women with breast cancer who have multiple affected relatives are more likely to have inherited genetic risk factors for the disease. All the currently known genetic risk factors for breast cancer account for less than half of the average familial risk. Furthermore, the genetic factor(s) underlying an increased cancer risk for many women from multiple-case families remain unknown. Rare genomic duplications and deletions, known as copy number variants (CNVs), cover more than 10% of a human genome, are often not assessed in studies of genetic predisposition, and could account for some of the so-called "missing heritability".

METHODS

We carried out a hypothesis-generating case-control study of breast cancer diagnosed before age 40 years (200 cases, 293 controls) using population-based cases from the Australian Breast Cancer Family Study. Genome-wide scanning for CNVs was performed using the Human610-Quad BeadChip and fine-mapping was conducted using PennCNV.

RESULTS

We identified deletions overlapping two known cancer susceptibility genes, (BRCA1 and BLM), and a duplication overlapping SMARCB1, associated with risk. The number of deletions across the genome was 1.5-fold higher for cases than controls (P = 10^-16), and 2-fold higher when only rare deletions overlapping genes (frequency <1%) were assessed (P = 5 × 10^-4). Association tests of CNVs, followed by experimental validation of CNV calls, found deletions overlapping the OR4C11 and OR4P4 genes were associated with breast cancer (P = 0.02 and P = 0.03, respectively).

CONCLUSION

These results suggest rare CNVs might have a role in breast cancer susceptibility, at least for disease at a young age.

The Role of Constitutional Copy Number Variants in Breast Cancer

Constitutional copy number variants (CNVs) include inherited and de novo deviations from a diploid state at a defined genomic region. These variants contribute significantly to genetic variation and disease in humans, including breast cancer susceptibility. Identification of genetic risk factors for breast cancer in recent years has been dominated by the use of genome-wide technologies, such as single nucleotide polymorphism (SNP)-arrays, with a significant focus on single nucleotide variants. To date, these large datasets have been underutilised for generating genome-wide CNV profiles despite offering a massive resource for assessing the contribution of these structural variants to breast cancer risk. Technical challenges remain in determining the location and distribution of CNVs across the human genome due to the accuracy of computational prediction algorithms and resolution of the array data. Moreover, better methods are required for interpreting the functional effect of newly discovered CNVs. In this review, we explore current and future application of SNP array technology to assess rare and common CNVs in association with breast cancer risk in humans.

E-cadherin loss alters cytoskeletal organization and adhesion in non-malignant breast cells but is insufficient to induce an epithelial-mesenchymal transition

BACKGROUND

E-cadherin is an adherens junction protein that forms homophilic intercellular contacts in epithelial cells while also interacting with the intracellular cytoskeletal networks. It has roles including establishment and maintenance of cell polarity, differentiation, migration and signalling in cell proliferation pathways. Its downregulation is commonly observed in epithelial tumours and is a hallmark of the epithelial to mesenchymal transition (EMT).

METHODS

To improve our understanding of how E-cadherin loss contributes to tumorigenicity, we investigated the impact of its elimination from the non-tumorigenic breast cell line MCF10A. We performed cell-based assays and whole genome RNAseq to characterize an isogenic MCF10A cell line that is devoid of CDH1 expression due to an engineered homozygous 4 bp deletion in CDH1 exon 11.

RESULTS

The E-cadherin-deficient line, MCF10A CDH1-/- showed subtle morphological changes, weaker cell-substrate adhesion, delayed migration, but retained cell-cell contact, contact growth inhibition and anchorage-dependent growth. Within the cytoskeleton, the apical microtubule network in the CDH1-deficient cells lacked the radial pattern of organization present in the MCF10A cells and F-actin formed thicker, more numerous stress fibres in the basal part of the cell. Whole genome RNAseq identified compensatory changes in the genes involved in cell-cell adhesion while genes involved in cell-substrate adhesion, notably ITGA1, COL8A1, COL4A2 and COL12A1, were significantly downregulated. Key EMT markers including CDH2, FN1, VIM and VTN were not upregulated although increased expression of proteolytic matrix metalloprotease and kallikrein genes was observed.

CONCLUSIONS

Overall, our results demonstrated that E-cadherin loss alone was insufficient to induce an EMT or enhance transforming potential in the non-tumorigenic MCF10A cells but was associated with broad transcriptional changes associated with tissue remodelling.

E-cadherin loss alters cytoskeletal organization and adhesion in non-malignant breast cells but is insufficient to induce an epithelial-mesenchymal transition

Background

E-cadherin is an adherens junction protein that forms homophilic intercellular contacts in epithelial cells while also interacting with the intracellular cytoskeletal networks. It has roles including establishment and maintenance of cell polarity, differentiation, migration and signalling in cell proliferation pathways. Its downregulation is commonly observed in epithelial tumours and is a hallmark of the epithelial to mesenchymal transition (EMT).

Methods

To improve our understanding of how E-cadherin loss contributes to tumorigenicity, we investigated the impact of its elimination from the non-tumorigenic breast cell line MCF10A. We performed cell-based assays and whole genome RNAseq to characterize an isogenic MCF10A cell line that is devoid of CDH1 expression due to an engineered homozygous 4 bp deletion in CDH1 exon 11.

Results

The E-cadherin-deficient line, MCF10A CDH1-/- showed subtle morphological changes, weaker cell-substrate adhesion, delayed migration, but retained cell-cell contact, contact growth inhibition and anchorage-dependent growth. Within the cytoskeleton, the apical microtubule network in the CDH1-deficient cells lacked the radial pattern of organization present in the MCF10A cells and F-actin formed thicker, more numerous stress fibres in the basal part of the cell. Whole genome RNAseq identified compensatory changes in the genes involved in cell-cell adhesion while genes involved in cell-substrate adhesion, notably ITGA1, COL8A1, COL4A2 and COL12A1, were significantly downregulated. Key EMT markers including CDH2, FN1, VIM and VTN were not upregulated although increased expression of proteolytic matrix metalloprotease and kallikrein genes was observed.

Conclusions

Overall, our results demonstrated that E-cadherin loss alone was insufficient to induce an EMT or enhance transforming potential in the non-tumorigenic MCF10A cells but was associated with broad transcriptional changes associated with tissue remodelling.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2407-14-552) contains supplementary material, which is available to authorized users.

E-cadherin loss alters cytoskeletal organization and adhesion in non-malignant breast cells but is insufficient to induce an epithelial-mesenchymal transition

BACKGROUND

E-cadherin is an adherens junction protein that forms homophilic intercellular contacts in epithelial cells while also interacting with the intracellular cytoskeletal networks. It has roles including establishment and maintenance of cell polarity, differentiation, migration and signalling in cell proliferation pathways. Its downregulation is commonly observed in epithelial tumours and is a hallmark of the epithelial to mesenchymal transition (EMT).

METHODS

To improve our understanding of how E-cadherin loss contributes to tumorigenicity, we investigated the impact of its elimination from the non-tumorigenic breast cell line MCF10A. We performed cell-based assays and whole genome RNAseq to characterize an isogenic MCF10A cell line that is devoid of CDH1 expression due to an engineered homozygous 4 bp deletion in CDH1 exon 11.

RESULTS

The E-cadherin-deficient line, MCF10A CDH1-/- showed subtle morphological changes, weaker cell-substrate adhesion, delayed migration, but retained cell-cell contact, contact growth inhibition and anchorage-dependent growth. Within the cytoskeleton, the apical microtubule network in the CDH1-deficient cells lacked the radial pattern of organization present in the MCF10A cells and F-actin formed thicker, more numerous stress fibres in the basal part of the cell. Whole genome RNAseq identified compensatory changes in the genes involved in cell-cell adhesion while genes involved in cell-substrate adhesion, notably ITGA1, COL8A1, COL4A2 and COL12A1, were significantly downregulated. Key EMT markers including CDH2, FN1, VIM and VTN were not upregulated although increased expression of proteolytic matrix metalloprotease and kallikrein genes was observed.

CONCLUSIONS

Overall, our results demonstrated that E-cadherin loss alone was insufficient to induce an EMT or enhance transforming potential in the non-tumorigenic MCF10A cells but was associated with broad transcriptional changes associated with tissue remodelling.