Major regulators of the multi-step metastatic process are potential therapeutic targets for breast cancer management.

Metastasis is a multi-step process that leads to the dissemination of tumor cells to new sites and, consequently, to multi-organ neoplasia. Although most lethal breast cancer cases are related to metastasis occurrence, little is known about the dysregulation of each step, and clinicians still lack reliable therapeutic targets for metastasis impairment. To fill these gaps, we constructed and analyzed gene regulatory networks for each metastasis step (cell adhesion loss, epithelial-to-mesenchymal transition, and angiogenesis). Through topological analysis, we identified E2F1, EGR1, EZH2, JUN, TP63, and miR-200c-3p as general hub-regulators, FLI1 for cell-adhesion loss specifically, and TRIM28, TCF3, and miR-429 for angiogenesis. Applying the FANMOD algorithm, we identified 60 coherent feed-forward loops regulating metastasis-related genes associated with distant metastasis-free survival prediction. miR-139-5p, miR-200c-3p, miR-454-3p, and miR-1301-3p, among others, were the FFL's mediators. The expression of the regulators and mediators was observed to impact overall survival and to go along with metastasis occurrence. Lastly, we selected 12 key regulators and observed that they are potential therapeutic targets for canonical and candidate antineoplastics and immunomodulatory drugs, like trastuzumab, goserelin, and calcitriol. Our results highlight the relevance of miRNAs in mediating feed-forward loops and regulating the expression of metastasis-related genes. Altogether, our results contribute to understanding the multi-step metastasis complexity and identifying novel therapeutic targets and drugs for breast cancer management.

Novel lncRNAs Co-Expression Networks Identifies LINC00504 with Oncogenic Role in Luminal A Breast Cancer Cells.

Long non-coding RNAs (lncRNAs) are functional transcripts with more than 200 nucleotides. These molecules exhibit great regulatory capacity and may act at different levels of gene expression regulation. Despite this regulatory versatility, the biology of these molecules is still poorly understood. Computational approaches are being increasingly used to elucidate biological mechanisms in which these lncRNAs may be involved. Co-expression networks can serve as great allies in elucidating the possible regulatory contexts in which these molecules are involved. Herein, we propose the use of the pipeline deposited in the RTN package to build lncRNAs co-expression networks using TCGA breast cancer (BC) cohort data. Worldwide, BC is the most common cancer in women and has great molecular heterogeneity. We identified an enriched co-expression network for the validation of relevant cell processes in the context of BC, including LINC00504. This lncRNA has increased expression in luminal subtype A samples, and is associated with prognosis in basal-like subtype. Silencing this lncRNA in luminal A cell lines resulted in decreased cell viability and colony formation. These results highlight the relevance of the proposed method for the identification of lncRNAs in specific biological contexts.

Association between SNP rs527616 in lncRNA AQP4-AS1 and susceptibility to breast cancer in a southern Brazilian population.

Breast cancer (BC) is the leading cause of death by this disease in women worldwide. Among the factors involved in tumorigenesis, long non-coding RNAs (lncRNAs) and their differential expression have been associated. Differences in gene expression may be triggered by variations in DNA sequence, including single nucleotide polymorphisms (SNPs). In the present study, we analyzed the rs527616 (C>G), located in the lncRNA AQP4-AS1, using PCR-SSP in 306 BC patients and 312 controls, from a Brazilian population. In the BC group, the frequency found for CG heterozygotes was above the expected and the overdominant model is the best one to explain our results (OR: 1.70, IC 95%: 1.23-2.34, P<0.001). Furthermore, the SNP were associated with age at BC diagnosis and the risk genotype more frequent in the older age group. According to TCGA data, AQP4-AS1 is down-regulated in BC tissue, and the overexpression is associated with better prognoses, including Luminal A, HER2-, stage 1 of disease and smaller tumor. In conclusion, the CG genotype is associated with increased susceptibility in the southern Brazilian population. This SNP is mapped in the lncRNA AQP4-AS1, showing differential expression in BC samples. Based on these results, we emphasize the potential of the role of AQP4-AS1 in cancer.

Clinical and functional analysis of the germline TP53 p.K164E acetylation site variant.

TP53 plays a critical role as a tumor suppressor by controlling cell cycle progression, DNA repair, and apoptosis. Post-translational modifications such as acetylation of specific lysine residues in the DNA binding and carboxy-terminus regulatory domains modulate its tumor suppressor activities. In this study, we addressed the functional consequences of the germline TP53 p.K164E (NM_000546.5: c.490A>G) variant identified in a patient with early-onset breast cancer and a significant family history of cancer. K164 is a conserved residue located in the L2 loop of the p53 DNA binding domain that is post-translationally modified by acetylation. In silico, in vitro, and in vivo analyses demonstrated that the glutamate substitution at K164 marginally destabilizes the p53 protein structure but significantly impairs sequence-specific DNA binding, transactivation, and tumor cell growth inhibition. Although p.K164E is currently considered a variant of unknown significance by different clinical genetic testing laboratories, the clinical and laboratory-based findings presented here provide strong evidence to reclassify TP53 p.K164E as a likely pathogenic variant.

MiR-150-5p Overexpression in Triple-Negative Breast Cancer Contributes to the In Vitro Aggressiveness of This Breast Cancer Subtype.

MiR-150-5p is frequently deregulated in cancer, with expression and mode of action varying according to the tumor type. Here, we investigated the expression levels and role of miR-150-5p in the aggressive breast cancer subtype triple-negative breast cancer (TNBC). MiR-150-5p expression levels were analyzed in tissue samples from 113 patients with invasive breast cancer (56 TNBC and 57 non-TNBC) and 41 adjacent non-tumor tissues (ANT). Overexpression of miR-150-5p was observed in tumor tissues compared with ANT tissues and in TNBC compared with non-TNBC tissues. MiR-150-5p expression levels were significantly associated with high tumor grades and the Caucasian ethnicity. Interestingly, high miR-150-5p levels were associated with prolonged overall survival. Manipulation of miR-150-5p expression in TNBC cells modulated cell proliferation, clonogenicity, migration, and drug resistance. Manipulation of miR-150-5p expression also resulted in altered expression of its mRNA targets, including epithelial-to-mesenchymal transition markers, MYB, and members of the SRC pathway. These findings suggest that miR-150-5p is overexpressed in TNBC and contributes to the aggressiveness of TNBC cells in vitro.

Transcribed Ultraconserved Regions Are Associated with Clinicopathological Features in Breast Cancer.

Ultraconserved regions (UCRs) are 481 genome segments, with length longer than 200 bp, that are 100% conserved among humans, mice, and rats. The majority of UCRs are transcriptionally active (T-UCRs) as many of them produce non-coding RNAs. In a previous study, we evaluated the expression level of T-UCRs in breast cancer (BC) patients and found that 63% of transcripts correlated with some clinical and/or molecular parameter of BC. In this study, we delved into the expression levels of 12 T-UCRs and correlated them with clinicopathological parameters, immunohistochemical markers, and overall survival in two breast cancer cohorts: TCGA and Brazilian patients. We found that uc.268 is more expressed in TCGA patients under 40 years of age, associated with progesterone receptor (PR) and estrogen receptor (ER), and its high expression is found in luminal A. Lower uc.84 and uc.376 were respectively observed in metastatic and stage IV tumors associated with good prognostic in luminal B. Moreover, uc.84 was only related to the HER2+, while uc.376 was related to ER+ and PR+, and HER2+. A panel composed of uc.147, uc.271, and uc.427 distinguished luminal A from triple negative patients with an AUC of 0.9531 (sensitivity 92.19% and specificity 86.76%). These results highlight the potential role of T-UCRs in BC and provide insights into the potential application of T-UCRs as biomarkers.

Genome interaction of the virus and the host genes and non-coding RNAs in SARS-CoV-2 infection.

In this review, we highlight the interaction of SARS-CoV-2 virus and host genomes, reporting the current studies on the sequence analysis of SARS-CoV-2 isolates and host genomes from diverse world populations. The main genetic variants that are present in both the virus and host genomes were particularly focused on the ACE2 and TMPRSS2 genes, and their impact on the patients' susceptibility to the virus infection and severity of the disease. Finally, the interaction of the virus and host non-coding RNAs is described in relation to their regulatory roles in target genes and/or signaling pathways critically associated with SARS-CoV-2 infection. Altogether, these studies provide a significant contribution to the knowledge of SARS-CoV-2 mechanisms of infection and COVID-19 pathogenesis. The described genetic variants and molecular factors involved in host/virus genome interactions have significantly contributed to defining patient risk groups, beyond those based on patients' age and comorbidities, and they are promising candidates to be potentially targeted in treatment strategies for COVID-19 and other viral infectious diseases.

Simultaneous CXCL12 and ESR1 CpG island hypermethylation correlates with poor prognosis in sporadic breast cancer.

BACKGROUND: CXCL12 is a chemokine that is constitutively expressed in many organs and tissues. CXCL12 promoter hypermethylation has been detected in primary breast tumours and contributes to their metastatic potential. It has been shown that the oestrogen receptor alpha (ESR1) gene can also be silenced by DNA methylation. In this study, we used methylation-specific PCR (MSP) to analyse the methylation status in two regions of the CXCL12 promoter and ESR1 in tumour cell lines and in primary breast tumour samples, and correlated our results with clinicopathological data. METHODS: First, we analysed CXCL12 expression in breast tumour cell lines by RT-PCR. We also used 5-aza-2'-deoxycytidine (5-aza-CdR) treatment and DNA bisulphite sequencing to study the promoter methylation for a specific region of CXCL12 in breast tumour cell lines. We evaluated CXCL12 and ESR1 methylation in primary tumour samples by methylation-specific PCR (MSP). Finally, promoter hypermethylation of these genes was analysed using Fisher's exact test and correlated with clinicopathological data using the Chi square test, Kaplan-Meier survival analysis and Cox regression analysis. RESULTS: CXCL12 promoter hypermethylation in the first region (island 2) and second region (island 4) was correlated with lack of expression of the gene in tumour cell lines. In the primary tumours, island 2 was hypermethylated in 14.5% of the samples and island 4 was hypermethylated in 54% of the samples. The ESR1 promoter was hypermethylated in 41% of breast tumour samples. In addition, the levels of ER alpha protein expression diminished with increased frequency of ESR1 methylation (p < 0.0001). This study also demonstrated that CXCL12 island 4 and ESR1 methylation occur simultaneously at a high frequency (p = 0.0220). CONCLUSIONS: This is the first study showing a simultaneous involvement of epigenetic regulation for both CXCL12 and ESR1 genes in Brazilian women. The methylation status of both genes was significantly correlated with histologically advanced disease, the presence of metastases and death. Therefore, the methylation pattern of these genes could be used as a molecular marker for the prediction of breast cancer outcome.

Identification of miRNAs Enriched in Extracellular Vesicles Derived from Serum Samples of Breast Cancer Patients.

MicroRNAs derived from extracellular vesicles (EV-miRNAs) are circulating miRNAs considered as potential new diagnostic markers for cancer that can be easily detected in liquid biopsies. In this study, we performed RNA sequencing analysis as a screening strategy to identify EV-miRNAs derived from serum of clinically well-annotated breast cancer (BC) patients from the south of Brazil. EVs from three groups of samples (healthy controls (CT), luminal A (LA), and triple-negative (TNBC)) were isolated from serum using a precipitation method and analyzed by RNA-seq (screening phase). Subsequently, four EV-miRNAs (miR-142-5p, miR-150-5p, miR-320a, and miR-4433b-5p) were selected to be quantified by quantitative real-time PCR (RT-qPCR) in individual samples (test phase). A panel composed of miR-142-5p, miR-320a, and miR-4433b-5p distinguished BC patients from CT with an area under the curve (AUC) of 0.8387 (93.33% sensitivity, 68.75% specificity). The combination of miR-142-5p and miR-320a distinguished LA patients from CT with an AUC of 0.9410 (100% sensitivity, 93.80% specificity). Interestingly, decreased expression of miR-142-5p and miR-150-5p were significantly associated with more advanced tumor grades (grade III), while the decreased expression of miR-142-5p and miR-320a was associated with a larger tumor size. These results provide insights into the potential application of EVs-miRNAs from serum as novel specific markers for early diagnosis of BC.

XAF1 as a modifier of p53 function and cancer susceptibility.

Cancer risk is highly variable in carriers of the common TP53-R337H founder allele, possibly due to the influence of modifier genes. Whole-genome sequencing identified a variant in the tumor suppressor XAF1 (E134*/Glu134Ter/rs146752602) in a subset of R337H carriers. Haplotype-defining variants were verified in 203 patients with cancer, 582 relatives, and 42,438 newborns. The compound mutant haplotype was enriched in patients with cancer, conferring risk for sarcoma (P = 0.003) and subsequent malignancies (P = 0.006). Functional analyses demonstrated that wild-type XAF1 enhances transactivation of wild-type and hypomorphic TP53 variants, whereas XAF1-E134* is markedly attenuated in this activity. We propose that cosegregation of XAF1-E134* and TP53-R337H mutations leads to a more aggressive cancer phenotype than TP53-R337H alone, with implications for genetic counseling and clinical management of hypomorphic TP53 mutant carriers.

Integrated copy number and miRNA expression analysis in triple negative breast cancer of Latin American patients.

Triple negative breast cancer (TNBC), a clinically aggressive breast cancer subtype, affects 15-35% of women from Latin America. Using an approach of direct integration of copy number and global miRNA profiling data, performed simultaneously in the same tumor specimens, we identified a panel of 17 miRNAs specifically associated with TNBC of ancestrally characterized patients from Latin America, Brazil. This panel was differentially expressed between the TNBC and non-TNBC subtypes studied (p ≤ 0.05, FDR ≤ 0.25), with their expression levels concordant with the patterns of copy number alterations (CNAs), present mostly frequent at 8q21.3-q24.3, 3q24-29, 6p25.3-p12.2, 1q21.1-q44, 5q11.1-q22.1, 11p13-p11.2, 13q12.11-q14.3, 17q24.2-q25.3 and Xp22.33-p11.21. The combined 17 miRNAs presented a high power (AUC = 0.953 (0.78-0.99);95% CI) in discriminating between the TNBC and non-TNBC subtypes of the patients studied. In addition, the expression of 14 and 15 of the 17miRNAs was significantly associated with tumor subtype when adjusted for tumor stage and grade, respectively. In conclusion, the panel of miRNAs identified demonstrated the impact of CNAs in miRNA expression levels and identified miRNA target genes potentially affected by both CNAs and miRNA deregulation. These targets, involved in critical signaling pathways and biological functions associated specifically with the TNBC transcriptome of Latina patients, can provide biological insights into the observed differences in the TNBC clinical outcome among racial/ethnic groups, taking into consideration their genetic ancestry.

Cost-effective and fast KIR gene-content genotyping by multiplex melting curve analysis.

Killer cell immunoglobulin-like receptor (KIR) genes encode cell surface molecules that recognize HLA molecules and modulate the activity of natural killer (NK) cells. KIR genes exhibit presence and absence polymorphism, which generates a variety of gene-content haplotypes in worldwide populations. KIR gene-content variation is implicated in many diseases and is also important for placentation and transplantation. Because of the complexity of KIR polymorphism, variation in this family is still mostly studied at the gene-content level, even with the advent of next-generation sequencing (NGS) methods. Gene-content determination is generally expensive and/or time-consuming. To overcome these difficulties, we developed a method based on multiplex polymerase chain reaction with specific sequence primers (PCR-SSP) followed by melting curve analysis that allows cost-effective, precise and fast generation of results. Our method was 100% concordant with a gel-based method and 99.9% concordant with presence and absence determination by NGS. The limit of detection for accurate typing was 30 ng of DNA (0.42 µM) with 260/230 and 260/280 ratios as low as 0.19 and of 0.44. In addition, we developed a user-friendly Java-based computational application called killerPeak that interprets the raw data generated by Viia7 or QuantStudio 7 quantitative PCR machines and reliably exports the final genotyping results in spreadsheet file format. The combination of a reliable method that requires low amount of DNA with an automated interpretation of results allows scaling the KIR genotyping in large cohorts with reduced turnaround time.

Long Non-Coding RNA TUG1 Expression Is Associated with Different Subtypes in Human Breast Cancer.

Taurine upregulated 1 gene (TUG1) is a long non-coding RNA associated with several types of cancer. Recently, differential expression of TUG1 was found in cancerous breast tissues and associated with breast cancer malignancy features. Although this is evidence of a potential role in breast cancer, TUG1 expression could not be associated with different subtypes, possibly due to the small number of samples analyzed. Breast cancer is a heterogeneous disease and, based on molecular signatures, may be classified into different subtypes with prognostic implications. In the present study, we include analysis of TUG1 expression in 796 invasive breast carcinoma and 105 normal samples of RNA sequencing (RNA-seq) datasets from The Cancer Genome Atlas (TCGA) and describe that TUG1 expression is increased in HER2-enriched and basal-like subtypes compared to luminal A. Additionally, TUG1 expression is associated with survival in HER2-enriched patients. These results reinforce the importance of TUG1 in breast cancer and outline its potential impact on specific subtypes.

Copy number and expression analysis of FOSL1, GSTP1, NTSR1, FADD and CCND1 genes in primary breast tumors with axillary lymph node metastasis.

In breast cancer, lymph node (LN) metastasis is one of the strongest prognostic factors at diagnosis. Therefore the identification of molecular markers with metastatic potential that promote the development of LN metastasis is of critical clinical relevance. In this study, we evaluated the copy number status of the FOSL1, GSTP1, NTSR1, FADD and CCND1 genes by TaqMan assays in 137 breast cancer patients, 84 with LN metastasis (LN+) and 53 with no LN metastasis (LN-). The copy number data for four of these genes (FOSL1, GSTP1, FADD and CCND1) were integrated with their mRNA expression levels in 31 patients. In both groups of patients, gains were the most frequent copy number alteration (CNA) observed, involving mainly the CCND1, NTSR1 and FADD genes; mRNA overexpression was more commonly observed for the CCND1 and FADD genes. For the FADD gene in the LN+ group, gene expression was shown to be dependent on CNAs; for the other genes no association was found. In conclusion, increase copy number and mRNA overexpression of FOSL1, GSTP1, FADD, NTSR1 and CCND1 genes are frequently observed in primary breast tumors, and except for the FADD gene, they occur independently and irrespectively of the patients' LN axillary metastatic status.

Worse prognosis in breast cancer patients can be predicted by immunohistochemical analysis of positive MMP-2 and negative estrogen and progesterone receptors.

INTRODUCTION:: Breast cancer is the most cause of death, and approximately 90% of these deaths are due to metastases. Matrix metalloproteinase-2 (MMP-2) gelatinase activity is able to degrade a major constituent of the tumor microenvironment, type IV collagen. Two well-established proteins used as markers in clinical practice for breast cancer are the receptors for estrogen (ER) and progesterone (PR). Although the presence of these receptors has been associated with a better prognosis, loss of these proteins can occur during tumor progression, with subsequent resistance to hormone therapy. OBJECTIVE:: To study the correlation among MMP-2, ER, and PR, as well as the establishment of the metastatic process in primary breast tumors. METHOD:: Breast cancer samples (n=44) were analyzed by immunohistochemistry for MMP-2, ER, and PR. RESULTS:: We observed that 90% of patients who had metastases and died showed positive staining for MMP-2 (p=0.0082 for both). Using Kaplan-Meier analysis, we found that negative ER patients who were also positive for MMP-2 had even worse disease-free survival (DFS) and overall survival (OS) (p= 0.012 and p=0.005, respectively). Similar results were found in PR-negative patients for DFS (a trend p=0.077) and OS (p=0.038). CONCLUSION:: Regardless of our small sample size (n=44), the data obtained strongly suggest that MMP-2 in combination with already well-established markers could help to predict the emergence of metastases and death in patients with breast cancer.

Comparative Proteomics of Tumor and Paired Normal Breast Tissue Highlights Potential Biomarkers in Breast Cancer.

BACKGROUND/AIM: Breast cancer is the most common type of cancer among women worldwide, and about 57,000 new cases are expected for the Brazilian population in 2015. Elucidation of protein expression and modification is essential for the biological understanding, early diagnosis and therapeutics of breast cancer. The main objectives of the study are comparison between the proteome of tumor and paired non-tumor breast cancer tissues, describing all identified proteins, highlighting the ones most differentially expressed and comparing the data with existing literature. MATERIALS AND METHODS: The five paired samples from patients with invasive ductal carcinoma were analyzed by 2-DE and MS. RESULTS: We collected 161 identified spots corresponding to 110 distinct proteins. Forty-three differentially-expressed spots were common to at least two samples, and the ten proteins with the highest-fold changes were CASPE, ENOG, TPM1, CAPG, VIME, TPM3, TRFE, PDIA6, WDR61 and PDIA3. Metabolic enzymes and proteins with binding functions were the most representative functional classes of proteins with increased and decreased expression in tumor tissue respectively. CONCLUSION: Taking the fold change as a parameter, we point to future targets to be studied by functional methods in a search for biomarkers for initiation and progress of breast cancer.

Prevalence of an inherited cancer predisposition syndrome associated with the germ line TP53 R337H mutation in Paraguay.

The tumor suppressor gene TP53 is the most frequently mutated gene in human cancer, and the germline TP53 R337H mutation is the most common mutation reported to date. However, this mutation is associated with a lower cumulative lifetime cancer risk than other mutations in the p53 DNA-binding domain. A detailed statistical analysis of 171,500 DNA tests in Brazilian neonates found that 0.27% of the general population is positive for this mutation, and some of the estimated 200,000 Brazilian R337H carriers in southern and southeastern Brazil have already developed cancer. The present study was designed to estimate R337H prevalence in neighboring Paraguay. To address this question, 10,000 dried blood samples stored in Guthrie cards since 2008 were randomly selected from the Paraguayan municipalities located at the border with Brazil. These samples were tested for R337H mutation using the PCR-restriction fragment length polymorphism assay. This germline mutation was detected in five samples (5/10,000), indicating that the total number of R337H carriers in Paraguay may be as high as 3500. Previous studies have shown that other countries (i.e., Portugal, Spain, and Germany) presented one family with this mutation, leading us to conclude that, besides Brazil and Paraguay, other countries may have multiple families carrying this mutation, which is an inherited syndrome that is difficult to control.

Two different karyotypes with 1q abnormalities in a patient with Fanconi anemia.

We report a case of Fanconi anemia in which cytogenetic analysis of bone marrow (BM) samples revealed two distinct karyotypes: 46,XY,dup(1)(q21q42), in the first sample and 46,XY,del(1)(q32) in the second, aspirated 7 months later after acute myeloid leukemia (AML) developed. We discuss the cytogenetic clonal fluctuation common in Fanconi anemia, with the Fanconi's anemia (FA) reports available in the literature. Interestingly, we have identified that del(1)(q32) has not been reported before in FA.

Cytogenetic findings in phyllodes tumor and fibroadenomas of the breast.

The cytogenetic data on fibroadenomas and cystosarcoma phyllodes tumor of the breast, which are both biphasic breast tumors composed of epithelial and stromal components, are quite limited. In this study, we report on clonal chromosomal alterations in three fibroadenomas and one cystosarcoma phyllodes analyzed by GTG banding. The fibroadenomas presented mostly numerical abnormalities involving chromosomes 16, 18, and 21. One case presented a deletion on 17p. The cystosarcoma phyllodes presented numerous numerical abnormalities, mostly chromosome gains, and several marker chromosomes.

Comprehensive cytogenetic evaluation of a mature ovarian teratoma case.

Mature ovarian teratomas are benign ovarian germ cell tumors that usually present with a normal karyotype. There are very few reports describing chromosomal abnormalities in these tumors, none of which are recurrent. In this study we report on a mature teratoma case with clonal chromosomal alterations which include monosomies of chromosomes 6, 14, 16, and 21; trisomies of chromosomes 14 and 21; and deletions of Xq, 5p, 16p, and 17p. Comparative genomic hybridization evaluation of the sample revealed a normal profile. These findings are discussed together with the cytogenetic reports on other cases of ovarian teratomas described in the literature.