Prognosis recovery score of apical surgery Guided Bone Regeneration using cone beam computed tomography and digital bioinformatics.

Objective: Guided Bone Regeneration (GBR) is a dental surgical procedure that uses barrier membranes to prevent soft tissue invasion and conduct new bone growth. This study aimed to define a Prognosis Recovery score (PR score) to objectively classify post-surgery responders from non-responder patients who underwent GBR using Cone Beam Computed Tomography (CBCT). Methods: This prospective-observational-longitudinal-cohort study recruited 250 individuals who were assigned to: Conventional-Apical-Surgery (CAS, n = 39), Apical-Surgery using human fascia lata Membrane placement (ASM, n = 42), and Apical-Surgery using human fascia lata Membrane placement and lyophilized allograft Bone powder (ASMB, n = 39); and Apical-Surgery using collagen membrane Porcine origin and Bovine Bone-matrix (ASPBB, n = 130), an independent external validation cohort. Surgery was performed, and evolution was monitored by CBCTs at 0, 6-, 12-, 18-, and 24 months post-surgery. Results: Normalized lesion volumes were calculated, and non-linear time evolution morphology curves were characterized. The three-time evolution bone growth patterns were: a linear tendency (PR0), "S'' shaped log-logistic (PR1), and "C" cellular growth (PR2). The treatment success rates were PR2-46 %, PR2-88 %, and PR2-95 %/PR1-5% for CAS, ASM, and ASMB groups. The xenograft ASPBB counterpart achieved PR2-92 % and PR1-8%. The score PR had a sensitivity, specificity, and accuracy of 100 %. Conclusions: Patients' treatment success can be quantitatively, objectively, and precisely predicted with the Prognosis Recovery score (using only two CBCTs), according to their biological response to allograft or xenograft materials (time-evolution bone growth curves), reducing cost and radiation exposure. Clinical significance: Through digital imaging and bioinformatic analysis of bone regeneration observed in CBCTs, we defined a Prognosis Recovery (PR) score using only two CBCT volume assessments (0 and 6 months). The PR score allowed us to define three time-evolution curves depending on the biomaterials used and to classify patients in a quantitative, objective, and accurate way.

Transvaginal ultrasound-guided oocyte retrieval in cattle: State-of-the-art and its impact on the in vitro fertilization embryo production outcome.

Transvaginal ultrasound-guided oocyte retrieval (commonly called OPU) and in vitro embryo production (IVP) in cattle has shown significant progress in recent years, in part, as a result of a better understanding of the full potential of these tools by end users. The combination of OPU and IVP (OPU-IVP) has been successfully and widely commercially used worldwide. The main advantages are a greater number of embryos and pregnancies per unit of time, faster genetic progress due to donor quick turn around and more elite sires mating combinations, larger spectrum of female age (calves, prepuberal, heifer, cow) and condition (open, pregnant) from which to retrieve oocytes, a reduced number of sperm (even sexed) required to fertilize the oocytes, among other benefits. OPU-IVP requires significant less donor preparation in comparison to conventional embryo transfer (<50% of usual FSH injections needed) to the extent of no stimulating hormones (FSH) are necessary. Donor synchronization, stimulation, OPU technique, oocyte competence, embryo performance, and its impact on cryopreservation and pregnancy are discussed.

Clinical implementation of a cervical cancer screening program via co-testing at a university hospital.

The Human Papillomavirus (HPV) test is a crucial technology for cervical cancer prevention because it enables programs to identify women with high-risk HPV infection who are at risk of developing cervical cancer. Current U.S. Preventive Services Task Force recommendations include cervical cancer screening every three years with cervical cytology alone or every five years with either high-risk HPV testing alone or high-risk HPV testing combined with cytology (co-testing). In Argentina, 7,548 new cervical cancer cases are diagnosed each year with 3,932 deaths attributed to this cause. Our study aims to show the clinical implementation of a cervical cancer screening program by concurrent HPV testing and cervical cytology (co-testing); and to evaluate the possible cervical cancer screening scenarios for Latin America, focusing on their performance and average cost. A cervical cancer screening five year program via co-testing algorithm (Hybrid-2-Capture/cytology) was performed on women aged 30-65 years old at a university hospital. Statistical analysis included a multinomial logistic regression, and two cancer screening classification alternatives were tested (cytology-reflex and HPV-reflex). A total of 2,273 women were included, 91.11% of the participants were double-negative, 2.55% double-positive, 5.90% positive-Hybrid-2-Capture-/negative-cytology, and 0.44% negative-Hybrid-2-Capture/positive-cytology. A thorough follow-up was performed in the positive-Hybrid-2-Capture group. Despite our efforts, 21 (10.93%) were lost, mainly because of changes on their health insurance coverage which excluded them from our screening algorithm. Of the 171 women with positive-Hybrid-2-Capture results and follow-up, 68 (39.77%) cleared the virus infection, 64 (37.43%) showed viral persistence, and 39 (22.81%) were adequately treated after detection via colposcopy/biopsy of histological HSIL (High-Grade Squamous Intraepithelial Lesion). The prevalence of high-risk HPV in this population was 192 women (8.45%), with HSIL histology detection rates of 17 per 1,000 screened women. A multinomial logistic regression analysis was performed over the women with positive-Hybrid-2-Capture considering the follow up (clearance, persistence and HSIL) as dependent variable, and the cytology test results (positive- or negative-cytology and Atypical Squamous Cells of Undetermined Significance, ASC-US) as independent variable. The model supported a direct association between cytology test results and follow up: negative-cytology/clearance, ASC-US/persistence, and positive-cytology/HSIL with the following probabilities of occurrence for these pairs 0.5, 0.647 and 0.647, respectively. Cytology could be considered a prognostic-factor in women with a positive-Hybrid-2-Capture. These findings suggest that the introduction of co-testing could diminish the burden of cervical cancer in low-and middle-income-countries, acting as a tool against inequity in healthcare.

Liquid biopsy for KRAS, NRAS and BRAF mutation testing in advanced colorectal cancer patients: the Argentinean experience.

Objective: To analyze the frequency of KRAS, NRAS and BRAF hotspot mutations in circulating tumor DNA (ctDNA) from patients with metastatic colorectal cancer (mCRC). Methods: Observational, descriptive and retrospective study in mCRC patients with available ctDNA-based genotype of KRAS, NRAS and BRAF. Results: The frequencies of plasma mutations for KRAS, NRAS and BRAF were 34% (± 7), 4% (± 3) and 4% (± 3), respectively. Median overall survival of plasma-tested RAS/BRAF-mutated patients was 26.6 months (95% CI: 14.4-not estimable [NE]), while RAS/BRAF wild-type patients did not reach the median survival during follow-up. Median progression-free survival for RAS/BRAF wild-type and RAS/BRAF-mutated patients was 12 (95% CI: 7-NE) and 4 months (95% CI: 4-NE), respectively. Conclusion: Our work supports the utility of KRAS, NRAS and BRAF analysis in liquid biopsy from mCRC patients.

Correction: Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights.

[This corrects the article DOI: 10.1371/journal.pone.0249773.].

Chromatin topology defines estradiol-primed progesterone receptor and PAX2 binding in endometrial cancer cells.

Estrogen (E2) and Progesterone (Pg), via their specific receptors (ERalpha and PR), are major determinants in the development and progression of endometrial carcinomas, However, their precise mechanism of action and the role of other transcription factors involved are not entirely clear. Using Ishikawa endometrial cancer cells, we report that E2 treatment exposes a set of progestin-dependent PR binding sites which include both E2 and progestin target genes. ChIP-seq results from hormone-treated cells revealed a non-random distribution of PAX2 binding in the vicinity of these estrogen-promoted PR sites. Altered expression of hormone regulated genes in PAX2 knockdown cells suggests a role for PAX2 in fine-tuning ERalpha and PR interplay in transcriptional regulation. Analysis of long-range interactions by Hi-C coupled with ATAC-seq data showed that these regions, that we call 'progestin control regions' (PgCRs), exhibited an open chromatin state even before hormone exposure and were non-randomly associated with regulated genes. Nearly 20% of genes potentially influenced by PgCRs were found to be altered during progression of endometrial cancer. Our findings suggest that endometrial response to progestins in differentiated endometrial tumor cells results in part from binding of PR together with PAX2 to accessible chromatin regions. What maintains these regions open remains to be studied.

Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights.

There has been limited study of Native American whole genome diversity to date, which impairs effective implementation of personalized medicine and a detailed description of its demographic history. Here we report high coverage whole genome sequencing of 76 unrelated individuals, from 27 indigenous groups across Mexico, with more than 97% average Native American ancestry. On average, each individual has 3.26 million Single Nucleotide Variants and short indels, that together comprise a catalog of 9,737,152 variants, 44,118 of which are novel. We report 497 common Single Nucleotide Variants (with allele frequency > 5%) mapped to drug responses and 316,577 in enhancer or promoter elements; interestingly we found some of these enhancer variants in PPARG, a nuclear receptor involved in highly prevalent health problems in Mexican population, such as obesity, diabetes, and insulin resistance. By detecting signals of positive selection we report 24 enriched key pathways under selection, most of them related to immune mechanisms. No missense variants in ACE2, the receptor responsible for the entry of the SARS CoV-2 virus, were found in any individual. Population genomics and phylogenetic analyses demonstrated stratification in a Northern-Central-Southern axis, with major substructure in the Central region. The Seri, a northern group with the most genetic divergence in our study, showed a distinctive genomic context with the most novel variants, and the most population specific genotypes. Genome-wide analysis showed that the average haplotype blocks are longer in Native Mexicans than in other world populations. With this dataset we describe previously undetected population level variation in Native Mexicans, helping to reduce the gap in genomic data representation of such groups.

Metropolitan age-specific mortality trends at borough and neighborhood level: The case of Mexico City.

Understanding the spatial and temporal patterns of mortality rates in a highly heterogeneous metropolis, is a matter of public policy interest. In this context, there is no, to the best of our knowledge, previous studies that correlate both spatio-temporal and age-specific mortality rates in Mexico City. Spatio-temporal Kriging modeling was used over five age-specific mortality rates (from the years 2000 to 2016 in Mexico City), to gain both spatial (borough and neighborhood) and temporal (year and trimester) data level description. Mortality age-specific patterns have been modeled using multilevel modeling for longitudinal data. Posterior tests were carried out to compare mortality averages between geo-spatial locations. Mortality correlation extends in all study groups for as long as 12 years and as far as 13.27 km. The highest mortality rate takes place in the Cuauhtémoc borough, the commercial, touristic and cultural core downtown of Mexico City. On the contrary, Tlalpan borough is the one with the lowest mortality rates in all the study groups. Post-productive mortality is the first age-specific cause of death, followed by infant, productive, pre-school and scholar groups. The combinations of spatio-temporal Kriging estimation and time-evolution linear mixed-effect models, allowed us to unveil relevant time and location trends that may be useful for public policy planning in Mexico City.

Candidate pharmacological treatments for substance use disorder and suicide identified by gene co-expression network-based drug repositioning.

Patients with substance use disorders (SUD) are at high risk to die by suicide. So far, the neurobiology of the suicide-SUD association has not been elucidated. This study aimed to identify potential pharmacological targets among hub genes from brain gene co-expression networks of individuals with SUD in a suicidal and non-suicidal context. Post-mortem samples from the prefrontal cortex of 79 individuals were analyzed. Individuals were classified into the following groups: suicides with SUD (n = 28), suicides without SUD (n = 23), nonsuicides with SUD (n = 9), nonsuicides without SUD (n = 19). Gene expression profiles were evaluated with the Illumina HumanHT-12 v4 array. Co-expression networks were constructed in WGCNA using the differentially expressed genes found in the comparisons: (a) suicides with and without SUD and (b) nonsuicides with and without SUD. Hub genes were selected for drug-gene interaction testing in the DGIdb database. Among drugs interacting with hub genes in suicides we found MAOA inhibitors and dextromethorphan. In the nonsuicide individuals, we found interactions with eglumegad and antipsychotics (olanzapine, clozapine, loxapine). Modafinil was found to interact with genes in both suicides and nonsuicides. These drugs represent possible candidate treatments for patients with SUD with and without suicidal behavior and their study in each context is encouraged.

Gene Co-expression Is Distance-Dependent in Breast Cancer.

Breast carcinomas are characterized by anomalous gene regulatory programs. As is well-known, gene expression programs are able to shape phenotypes. Hence, the understanding of gene co-expression may shed light on the underlying mechanisms behind the transcriptional regulatory programs affecting tumor development and evolution. For instance, in breast cancer, there is a clear loss of inter-chromosomal (trans-) co-expression, compared with healthy tissue. At the same time cis- (intra-chromosomal) interactions are favored in breast tumors. In order to have a deeper understanding of regulatory phenomena in cancer, here, we constructed Gene Co-expression Networks by using TCGA-derived RNA-seq whole-genome samples corresponding to the four breast cancer molecular subtypes, as well as healthy tissue. We quantify the cis-/trans- co-expression imbalance in all phenotypes. Additionally, we measured the association between co-expression and physical distance between genes, and characterized the ratio of intra/inter-cytoband interactions per phenotype. We confirmed loss of trans- co-expression in all molecular subtypes. We also observed that gene cis- co-expression decays abruptly with distance in all tumors in contrast with healthy tissue. We observed co-expressed gene hotspots, that tend to be connected at cytoband regions, and coincide accurately with already known copy number altered regions, such as Chr17q12, or Chr8q24.3 for all subtypes. Our methodology recovered different alterations already reported for specific breast cancer subtypes, showing how co-expression network approaches might help to capture distinct events that modify the cell regulatory program.

Gender imbalance in executive management positions at the Mexican National Institutes of Health.

BACKGROUND: Around the world, there is a significant difference in the proportion of women with access to leadership in healthcare with respect to men. This article studies gender imbalance and wage gap in managerial, executive, and directive job positions at the Mexican National Institutes of Health. METHODS: Cohort data were described using a visual circular representation and modeled using a generalized linear model. Analysis of variance was used to assess model significance, and posterior Fisher's least significant differences were analyzed when appropriate. RESULTS: This study demonstrated that there is a gender imbalance distribution among the hierarchical position at the Mexican National Health Institutes and also exposed that the wage gap exists mainly in the (highest or lowest) ranks in hierarchical order. CONCLUSIONS: Since the majority of the healthcare workforce is female, Mexican women are still underrepresented in executive and directive management positions at national healthcare organizations.

Sex differences in brain gene expression among suicide completers.

BACKGROUND: Suicide rates vary substantially by sex. Suicides committed by males significantly outnumber female suicides. Disparities in community and social factors provide a partial explanation for this phenomenon. Thus, the evaluation of sex differences at a biological level might contribute to the elucidation of the factors involved in this imbalance. The aim of the present study was to evaluate sex-specific gene expression patterns in the suicidal brain. METHODS: postmortem samples from the dorsolateral prefrontal cortex (DLPFC) of 75 Latino individuals were analyzed. We considered the following groups: i) male suicides (n = 38), ii) female suicides (n = 10), iii) male controls (n = 20), and iv) female controls (n = 7). Gene expression profiles were evaluated by microarrays. Differentially expressed genes among the groups were identified with a linear model. Similarities and differences in the gene sets between the sexes were identified. RESULTS: Differentially expressed genes were identified between suicides and controls of each sex: 1,729 genes in females and 1,997 genes in males. Female-exclusive suicide genes were related to cell proliferation and immune response. Meanwhile, male-exclusive suicide genes were associated to DNA binding and ribonucleic protein complex. Sex-independent suicide genes showed enrichment in mitochondrial and vesicular functions. LIMITATIONS: Relatively small sample size. Our diagnosis approach was limited to information found on coroner's records. The analysis was limited to a single brain area (DLPFC) and we used microarrays. CONCLUSION: Previously unexplored sex differences in the brain gene expression of suicide completers were identified, providing valuable foundation for the evaluation of sex-specific factors in suicide.

Brain Gene Expression Profiling of Individuals With Dual Diagnosis Who Died by Suicide.

Objective: Dual diagnosis (DD) is the co-occurrence of at least one substance use disorder and one or more mental disorders in a given individual. Despite this comorbidity being highly prevalent and associated with adverse clinical outcomes, its neurobiology remains unclear. Furthermore, patients with DD are at higher risk for suicidal behavior in comparison with single disorder patients. Our objective was to evaluate brain gene expression patterns in individuals with DD who died by suicide. Methods: We compared the gene expression profile in the dorsolateral prefrontal cortex of suicides with DD (n = 10) to the transcriptome of suicides with substance use disorder alone (n = 10), suicides with mood disorders (MD) alone (n = 13), and suicides without mental comorbidities (n = 5). Gene expression profiles were assessed by microarrays. In addition, we performed a brain cell type enrichment to evaluate whether the gene expression profiles could reflect differences in cell type compositions among the groups. Results: When comparing the transcriptome of suicides with DD to suicides with substance use disorder alone and suicides with MD alone, we identified 255 and 172 differentially expressed genes (DEG), respectively. The overlap of DEG between both comparisons (112 genes) highlighted the presence of common disrupted pathways in substance use disorder and MD. When comparing suicides with DD to suicides without mental comorbidities, we identified 330 DEG, mainly enriched in neurogenesis. Cell type enrichment indicated higher levels of glial markers in suicides with DD compared to the other groups. Conclusions: Suicides with DD exhibited a gene expression profile distinct from that of suicides with a single disorder, being substance use disorder or MD, and suicides without mental disorders. Our results suggest alteration in the expression of genes involved in glial specific markers, glutamatergic and GABAergic neurotransmission in suicides with DD compared to suicides with a single disorder and suicides without mental comorbidities. Alterations in the expression of synaptic genes at different levels were found in substance use disorder and MD.

Spatio-Temporal Representativeness of Air Quality Monitoring Stations in Mexico City: Implications for Public Health.

Assessment of the air quality in metropolitan areas is a major challenge in environmental sciences. Issues related include the distribution of monitoring stations, their spatial range, or missing information. In Mexico City, stations have been located spanning the entire Metropolitan zone for pollutants, such as CO, NO2, O3, SO2, PM2.5, PM10, NO, NO x , and PM CO . A fundamental question is whether the number and location of such stations are adequate to optimally cover the city. By analyzing spatio-temporal correlations for pollutant measurements, we evaluated the distribution and performance of monitoring stations in Mexico City from 2009 to 2018. Based on our analysis, air quality evaluation of those contaminants is adequate to cover the 16 boroughs of Mexico City, with the exception of SO2, since its spatial range is shorter than the one needed to cover the whole surface of the city. We observed that NO and NO x concentrations must be taken into account since their long-range dispersion may have relevant consequences for public health. With this approach, we may be able to propose policy based on systematic criteria to locate new monitoring stations.

Brain Gene Expression Pattern of Subjects with Completed Suicide and Comorbid Substance Use Disorder.

BACKGROUND/AIM: Although individuals with substance use disorder (SUD) are at high risk of committing suicide, most studies of postmortem gene expression exclude subjects with SUD due to the potential confounding effect of drugs in the transcriptome. Thus, little is known about the gene expression profile in suicides with SUD. The identification of altered biological processes in suicides with SUD is crucial in the comprehension of the interaction between both pathologies. METHODS: We evaluated the gene expression profile in the dorsolateral prefrontal area of suicides and nonsuicides with and without SUD by microarrays. RESULTS: We identified 222 differentially expressed genes, predominately enriched in cell proliferation in the comparison between suicides with and without SUD. When comparing the transcriptome of suicides with SUD to nonsuicides with SUD, we identified 550 differentially expressed genes, mainly enriched in oxidative phosphorylation. Differentially expressed genes (1,417) between suicides and nonsuicides without SUD were detected. Most of them were related to mitochondrial function. CONCLUSION: Interaction between suicide and SUD seems to influence the expression of genes involved in glial proliferation and glutamatergic neurotransmission. These results highlight, for the first time, that suicides with SUD have a gene expression profile distinct from that of subjects with only one of these disorders.

IL-8, GRO and MCP-1 produced by hepatocellular carcinoma microenvironment determine the migratory capacity of human bone marrow-derived mesenchymal stromal cells without affecting tumor aggressiveness.

New therapies are needed for advanced hepatocellular carcinoma (HCC) and the use of mesenchymal stromal cells (MSCs) carrying therapeutic genes is a promising strategy. HCC produce cytokines recruiting MSCs to the tumor milieu and modifying its biological properties. Our aim was to study changes generated on human MSCs exposed to conditioned media (CM) derived from human HCC fresh samples and xenografts. All CM shared similar cytokines expression pattern including CXCL1-2-3/GRO, CCL2/MCP-1 and CXCL8/IL-8 being the latter with the highest concentration. Neutralizing and knockdown experiments of CCL2/MCP-1, CXCL8/IL-8, CXCR1 and CXCR2 reduced in vitro MSC migration of ≥20%. Simultaneous CXCR1 and CXCR2 neutralization resulted in 50% of MSC migration inhibition. MSC stimulated with CM (sMSC) from HuH7 or HC-PT-5 showed a 2-fold increase of migration towards the CM compared with unstimulated MSC (usMSC). Gene expression profile of sMSC showed ~500 genes differentially expressed compared with usMSC, being 46 genes related with cell migration and invasion. sMSC increased fibroblasts and endothelial cells chemotaxis. Finally, sMSC with HuH7 CM and then inoculated in HCC tumor bearing-mice did not modify tumor growth. In this work we characterized factors produced by HCC responsible for the changes in MSC chemotactic capacity with would have an impact on therapeutic use of MSCs for human HCC.

RNA-Seq based genome-wide analysis reveals loss of inter-chromosomal regulation in breast cancer.

Breast cancer is a complex heterogeneous disease. Common hallmark features of cancer can be found. Their origin may be traced back to their intricate relationships governing regulatory programs during the development of this disease. To unveil distinctive features of the transcriptional regulation program in breast cancer, a pipeline for RNA-seq analysis in 780 breast cancer and 101 healthy breast samples, at gene expression and network level, was implemented. Inter-chromosomal relationships between genes resulted strikingly scarce in a cancer network, in comparison to its healthy counterpart. We suggest that inter-chromosomal regulation loss may be a novel feature in breast cancer. Additional evidence was obtained by independent validation in microarray and Hi-C data as well as supplementary computational analyses. Functional analysis showed upregulation in processes related to cell cycle and division; while migration, adhesion and cell-to-cell communication, were downregulated. Both the BRCA1 DNA repairing signalling and the Estrogen-mediated G1/S phase entry pathways were found upregulated. In addition, a synergistic underexpression of the γ-protocadherin complex, located at Chr5q31 is also shown. This region has previously been reported to be hypermethylated in breast cancer. These findings altogether provide further evidence for the central role of transcriptional regulatory programs in shaping malignant phenotypes.

TarSeqQC: Quality control on targeted sequencing experiments in R.

Targeted sequencing (TS) is growing as a screening methodology used in research and medical genetics to identify genomic alterations causing human diseases. In general, a list of possible genomic variants is derived from mapped reads through a variant calling step. This processing step is usually based on variant coverage, although it may be affected by several factors. Therefore, undercovered relevant clinical variants may not be reported, affecting pathology diagnosis or treatment. Thus, a prior quality control of the experiment is critical to determine variant detection accuracy and to avoid erroneous medical conclusions. There are several quality control tools, but they are focused on issues related to whole-genome sequencing. However, in TS, quality control should assess experiment, gene, and genomic region performances based on achieved coverages. Here, we propose TarSeqQC R package for quality control in TS experiments. The tool is freely available at Bioconductor repository. TarSeqQC was used to analyze two datasets; low-performance primer pools and features were detected, enhancing the quality of experiment results. Read count profiles were also explored, showing TarSeqQC's effectiveness as an exploration tool. Our proposal may be a valuable bioinformatic tool for routinely TS experiments in both research and medical genetics.

A novel non-parametric method for uncertainty evaluation of correlation-based molecular signatures: its application on PAM50 algorithm.

Motivation: The PAM50 classifier is used to assign patients to the highest correlated breast cancer subtype irrespectively of the obtained value. Nonetheless, all subtype correlations are required to build the risk of recurrence (ROR) score, currently used in therapeutic decisions. Present subtype uncertainty estimations are not accurate, seldom considered or require a population-based approach for this context. Results: Here we present a novel single-subject non-parametric uncertainty estimation based on PAM50's gene label permutations. Simulations results ( n = 5228) showed that only 61% subjects can be reliably 'Assigned' to the PAM50 subtype, whereas 33% should be 'Not Assigned' (NA), leaving the rest to tight 'Ambiguous' correlations between subtypes. The NA subjects exclusion from the analysis improved survival subtype curves discrimination yielding a higher proportion of low and high ROR values. Conversely, all NA subjects showed similar survival behaviour regardless of the original PAM50 assignment. We propose to incorporate our PAM50 uncertainty estimation to support therapeutic decisions. Availability and Implementation: Source code can be found in 'pbcmc' R package at Bioconductor. Contacts: cristobalfresno@gmail.com or efernandez@bdmg.com.ar. Supplementary information: Supplementary data are available at Bioinformatics online.

Predictive Outcomes for HER2-enriched Cancer Using Growth and Metastasis Signatures Driven By SPARC.

Understanding the mechanism of metastatic dissemination is crucial for the rational design of novel therapeutics. The secreted protein acidic and rich in cysteine (SPARC) is a matricellular glycoprotein which has been extensively associated with human breast cancer aggressiveness although the underlying mechanisms are still unclear. Here, shRNA-mediated SPARC knockdown greatly reduced primary tumor growth and completely abolished lung colonization of murine 4T1 and LM3 breast malignant cells implanted in syngeneic BALB/c mice. A comprehensive study including global transcriptomic analysis followed by biological validations confirmed that SPARC induces primary tumor growth by enhancing cell cycle and by promoting a COX-2-mediated expansion of myeloid-derived suppressor cells (MDSC). The role of SPARC in metastasis involved a COX-2-independent enhancement of cell disengagement from the primary tumor and adherence to the lungs that fostered metastasis implantation. Interestingly, SPARC-driven gene expression signatures obtained from these murine models predicted the clinical outcome of patients with HER2-enriched breast cancer subtypes. In total, the results reveal that SPARC and its downstream effectors are attractive targets for antimetastatic therapies in breast cancer.Implications: These findings shed light on the prometastatic role of SPARC, a key protein expressed by breast cancer cells and surrounding stroma, with important consequences for disease outcome. Mol Cancer Res; 15(3); 304-16. ©2016 AACR.