Frequency of carriers for rare recessive Mendelian diseases in a Brazilian cohort of 320 patients.

Several Mendelian disorders follow an autosomal recessive inheritance pattern. Epidemiological information on many inherited disorders may be useful to guide health policies for rare diseases, but it is often inadequate, particularly in developing countries. We aimed to calculate the carrier frequencies of rare autosomal recessive Mendelian diseases in a cohort of Brazilian patients using whole exome sequencing (WES). We reviewed the molecular findings of WES from 320 symptomatic patients who had carrier status for recessive diseases. Using the Hardy-Weinberg equation, we estimated recessive disease frequencies (q2 ) considering the respective carrier frequencies (2pq) observed in our study. We calculated the sensitivity of carrier screening tests based on lists of genes from five different clinical laboratories that offer them in Brazil. A total of 425 occurrences of 351 rare variants were reported in 278 different genes from 230 patients (71.9%). Almost half (48.8%) were carriers of at least one heterozygous pathogenic/likely pathogenic variant for rare metabolic disorders, while 25.9% of epilepsy, 18.1% of intellectual disabilities, 15.6% of skeletal disorders, 10.9% immune disorders, and 9.1% of hearing loss. We estimated that an average of 67% of the variants would not have been detected by carrier screening panels. The combined frequencies of autosomal recessive diseases were estimated to be 26.39/10,000 (or ~0.26%). This study shows the potential research utility of WES to determine carrier status, which may be a possible strategy to evaluate the clinical and social burden of recessive diseases at the population level and guide the optimization of carrier screening panels.

Exome sequencing and targeted gene panels: a simulated comparison of diagnostic yield using data from 158 patients with rare diseases.

Next-generation sequencing (NGS) has altered clinical genetic testing by widening the access to molecular diagnosis of genetically determined rare diseases. However, physicians may face difficulties selecting the best diagnostic approach. Our goal is to estimate the rate of possible molecular diagnoses missed by different targeted gene panels using data from a cohort of patients with rare genetic diseases diagnosed with exome sequencing (ES). For this purpose, we simulated a comparison between different targeted gene panels and ES: the list of genes harboring clinically relevant variants from 158 patients was used to estimate the theoretical rate of diagnoses missed by NGS panels from 53 different NGS panels from eight different laboratories. Panels presented a mean rate of missed diagnoses of 64% (range 14%-100%) compared to ES, representing an average predicted sensitivity of 36%. Metabolic abnormalities represented the group with highest mean of missed diagnoses (86%), while seizure represented the group with lowest mean (46%). Focused gene panels are restricted in covering select sets of genes implicated in specific diseases and they may miss molecular diagnoses of rare diseases compared to ES. However, their role in genetic diagnosis remains important especially for well-known genetic diseases with established genetic locus heterogeneity.

Diagnostic power and clinical impact of exome sequencing in a cohort of 500 patients with rare diseases.

Rare diseases comprise a diverse group of conditions, most of which involve genetic causes. We describe the variable spectrum of findings and clinical impacts of exome sequencing (ES) in a cohort of 500 patients with rare diseases. In total, 164 primary findings were reported in 158 patients, representing an overall diagnostic yield of 31.6%. Most of the findings (61.6%) corresponded to autosomal dominant conditions, followed by autosomal recessive (25.6%) and X-linked (12.8%) conditions. These patients harbored 195 variants, among which 43.6% are novel in the literature. The rate of molecular diagnosis was considerably higher for prenatal samples (67%; 4/6), younger children (44%; 24/55), consanguinity (50%; 3/6), gastrointestinal/liver disease (44%; 16/36) and syndromic/malformative conditions (41%; 72/175). For 15.6% of the cohort patients, we observed a direct potential for the redirection of care with targeted therapy, tumor screening, medication adjustment and monitoring for disease-specific complications. Secondary findings were reported in 37 patients (7.4%). Based on cost-effectiveness studies in the literature, we speculate that the reports of secondary findings may influence an increase of 123.2 years in the life expectancy for our cohort, or 0.246 years/cohort patient. ES is a powerful method to identify the molecular bases of monogenic disorders and redirect clinical care.

A comparison between SOLiD 5500XLand Ion Torrent PGM-derived miRNA expression profiles in two breast cell lines.

Next-generation sequencing (NGS) platforms allow the analysis of hundreds of millions of molecules in a single sequencing run, revolutionizing many research areas. NGS-based microRNA studies enable expression quantification in unprecedented scale without the limitations of closed-platforms. Yet, whereas a massive amount of data produced by these platforms is available, comparisons of quantification/discovery capabilities between platforms are still lacking. Here we compare two NGS-platforms: SOLiD and PGM, by evaluating their microRNA identification/quantification capabilities using two breast-derived cell-lines. A high expression correlation (R2 > 0.9) was achieved, encompassing 97% of the miRNAs, and the few discrepancies in miRNA counts were attributable to molecules that have very low expression. Quantification divergences indicative of artefactual representation were seen for 14 miRNAs (higher in SOLiD-reads) and another 10 miRNAs more abundant in PGM-data. An inspection of these revealed an increased and statistically significant count of uracyls and uracyl-stretches for PGM-enriched miRNAs, compared to SOLiD and to the miRBase. In parallel, adenines and adenine-stretches were enriched for SOLiDderived miRNA reads. We conclude that, whereas both platforms are overall consistent and can be used interchangeably for microRNA expression studies, particular sequence features appear to be indicative of specific platform bias, and their presence in microRNAs should be considered for database-analyses.

BRCA1 deficiency is a recurrent event in early-onset triple-negative breast cancer: a comprehensive analysis of germline mutations and somatic promoter methylation.

PURPOSE: BRCA1 germline mutation is closely associated with triple-negative breast cancer. BRCA deficiency leads to impaired DNA repair and tumor development, and understanding this deficiency, in both hereditary and sporadic scenarios, is of great clinical and biological interest. Here, we investigated germline or somatic events that might lead to BRCA1 impairment in triple-negative breast cancer. We also analyzed the clinical implications associated with BRCA deficiency. METHODS: Next-generation sequencing for the BRCA1/2 genes and multiplex ligation-dependent probe amplification (MLPA) for the BRCA1 gene were performed for mutation screening. A customized bisulfite next-generation sequencing approach was used for assessing BRCA1 promoter methylation status in tumor tissue. RESULTS: A total of 131 triple-negative cases were assessed, and germline pathogenic variants were detected in 13.0% of all cases and in 26% of cases diagnosed in young women. Most germline pathogenic variants (88.2%) occurred in the BRCA1 gene. BRCA1 promoter hypermethylation was detected in 20.6% of tumors; none of these tumors were in BRCA1/2 pathogenic variant carriers. BRCA1 impairment by either germline or somatic events was significantly more frequent in young women (55% in those ≤ 40 years; 33% in those 41-50 years; 22% in those > 50 years of age) and associated with better overall and disease-free survival rates in this group of patients. CONCLUSIONS: BRCA1 deficiency was recurrent in early-onset triple-negative breast cancer in Brazilian patients and associated with improved survival. With the new treatment modalities being investigated, including poly (ADP-ribose)-polymerase (PARP) inhibitor therapy, our results suggest that a significant proportion of young women with this subtype of tumor might benefit from PARP inhibitor treatment, which warrants further investigation.

A genomic case study of desmoplastic small round cell tumor: comprehensive analysis reveals insights into potential therapeutic targets and development of a monitoring tool for a rare and aggressive disease.

BACKGROUND: Genome-wide profiling of rare tumors is crucial for improvement of diagnosis, treatment, and, consequently, achieving better outcomes. Desmoplastic small round cell tumor (DSRCT) is a rare type of sarcoma arising from mesenchymal cells of abdominal peritoneum that usually develops in male adolescents and young adults. A specific translocation, t(11;22)(p13;q12), resulting in EWS and WT1 gene fusion is the only recurrent molecular hallmark and no other genetic factor has been associated to this aggressive tumor. Here, we present a comprehensive genomic profiling of one DSRCT affecting a 26-year-old male, who achieved an excellent outcome. METHODS: We investigated somatic and germline variants through whole-exome sequencing using a family based approach and, by array CGH, we explored the occurrence of genomic imbalances. Additionally, we performed mate-paired whole-genome sequencing for defining the specific breakpoint of the EWS-WT1 translocation, allowing us to develop a personalized tumor marker for monitoring the patient by liquid biopsy. RESULTS: We identified genetic variants leading to protein alterations including 12 somatic and 14 germline events (11 germline compound heterozygous mutations and 3 rare homozygous polymorphisms) affecting genes predominantly involved in mesenchymal cell differentiation pathways. Regarding copy number alterations (CNA) few events were detected, mainly restricted to gains in chromosomes 5 and 18 and losses at 11p, 13q, and 22q. The deletions at 11p and 22q indicated the presence of the classic translocation, t(11;22)(p13;q12). In addition, the mapping of the specific genomic breakpoint of the EWS-WT1 gene fusion allowed the design of a personalized biomarker for assessing circulating tumor DNA (ctDNA) in plasma during patient follow-up. This biomarker has been used in four post-treatment blood samples, 3 years after surgery, and no trace of EWS-WT1 gene fusion was detected, in accordance with imaging tests showing no evidence of disease and with the good general health status of the patient. CONCLUSIONS: Overall, our findings revealed genes with potential to be associated with risk assessment and tumorigenesis of this rare type of sarcoma. Additionally, we established a liquid biopsy approach for monitoring patient follow-up based on genomic information that can be similarly adopted for patients diagnosed with a rare tumor.

Upregulated genes at 2q24 gains as candidate oncogenes in hepatoblastomas.

AIM: Cytogenetic data of hepatoblastomas, a rare embryonal tumor of the liver, mostly consist of descriptions of whole-chromosome aneuploidies and large chromosome alterations. High-resolution cytogenetics may provide clues to hepatoblastoma tumorigenesis and indicate markers with clinical significance. PATIENTS & METHODS: We used array-CGH (180K) to screen for genomic imbalances in nine hepatoblastomas. Additionally, we investigated the expression pattern of selected genes exhibiting copy number changes. RESULTS: Analysis showed mainly whole-chromosome or chromosome-arm aneuploidies, but some focal aberrations were also mapped. Expression analysis of 48 genes mapped at one 10 Mb amplification at 2q24 revealed upregulation of DAPL1, ERMN, GALNT5, SCN1A and SCN3A in the set of tumors compared with differentiated livers. CONCLUSION: These genes appear as candidates for hepatoblastoma tumorigenesis.

Pheochromocytoma/Paraganglioma (PPGL): A Misdiagnosed Cause of Hypertension during Pregnancy.

Hypertension (HT) during pregnancy is not an infrequent obstetric problem, reaching a prevalence of 5-10%. This condition is highly associated with both maternal and fetal complications if not precisely diagnosed and managed. Even though primary HT, obesity, and preeclampsia are the main causes of HT in this period, other less familiar conditions must be considered during the investigation. Pheochromocytoma and paraganglioma (PPGL) are chromaffin cell tumors that produce, store, and secrete catecholamines, leading to HT and other adrenergic manifestations. Recognition of PPGL is crucial since misdiagnosis and improper management can lead to high morbidity and mortality, particularly during pregnancy. We report on two cases of PPGL diagnosed during pregnancy with different managements. Case 1 is a 25-year-old female at 31 weeks of first pregnancy, whose severe HT and life-threatening symptoms prompted an emergency delivery without previous confirmation or medical treatment of a suspected PPGL. After confirmation, a right adrenal PPGL was surgically resected 4 months later, following 15 days of medical therapy. Case 2 is a 22-year-old female at 18 weeks of pregnancy whose symptomatic PPGL was resected in the second trimester. A next-generation sequencing panel, including 23 PPGL-related genes, found no germline pathogenic variants (GPVs) in case 1 and an exon 1-4 germinative heterozygous deletion of the MAX gene in case 2. Despite the different medical approaches, both cases had satisfactory outcomes. Although uncommon, PPGL should be considered in the differential diagnosis of HT in pregnancy since missing the diagnosis and failing to introduce appropriate and timely treatment may lead to dramatic consequences for the mother and fetus. PPGL diagnosed during reproductive age is likely to result from GPV, prompting genetic investigation and counseling.

External beam radiation therapy for recurrent or residual thyroid cancer: What is the best treatment time and the best candidate for long-term local disease control?

INTRODUCTION: Cervical disease control might be challenging in advanced thyroid cancer (DTC). Indications for cervical external beam radiation therapy (EBRT) are controversial. PURPOSE: To identify clinical and molecular factors associated with control of cervical disease with EBRT. METHODS: Retrospective evaluation and molecular analysis of the primary tumor DTC patients who underwent cervical EBRT between 1995 and 2022 was performed. RESULTS: Eighty adults, median age of 61 years, were included. T4 disease was present in 43.7%, lymph node involvement in 42.5%, and distant metastasis in 47.5%. Those with cervical progression were older (62.5 vs. 57.3, p = 0.04) with more nodes affected (12.1 vs. 2.8, p = 0.04) and had EBRT performed later following surgery (76.6 vs. 64 months, p = 0.05). EBRT associated with multikinase inhibitors showed longer overall survival than EBRT alone (64.3 vs. 37.9, p = 0.018) and better local disease control. Performing EBRT before radioiodine (RAI) was associated with longer cervical progression-free survival (CPFS) than was RAI before (67.5 vs. 34.5, p < 0.01). EBRT ≥2 years after surgery was associated with worse CPFS (4.9 vs. 34, p = 0.04). The most common molecular alterations were ERBB2, BRAF, FAT1, RET and ROS1 and TERT mutation was predictive of worse disease control after EBRT (p = 0.04). CONCLUSION: Younger patients, with fewer affected nodes and treated earlier after surgery had better cervical disease control. Combination of EBRT with MKI improved OS. TERT mutation might indicate worse responders to EBRT; however, further studies are necessary to clarify the role of molecular testing in selecting candidates for cervical EBRT.

Impact of genetic background as a risk factor for atherosclerotic cardiovascular disease: A protocol for a nationwide genetic case-control (CV-GENES) study in Brazil.

Atherosclerotic Cardiovascular Disease (ASCVD) represents the leading cause of death worldwide, and individual screening should be based on behavioral, metabolic, and genetic profile derived from data collected in large population-based studies. Due to the polygenic nature of ASCVD, we aimed to assess the association of genomics with ASCVD risk and its impact on the occurrence of acute myocardial infarction, stroke, or peripheral artery thrombotic-ischemic events at population level. CardioVascular Genes (CV-GENES) is a nationwide, multicenter, 1:1 case-control study of 3,734 patients in Brazil. Inclusion criterion for cases is the first occurrence of one of the ASCVD events. Individuals without known ASCVD will be eligible as controls. A core lab will perform the genetic analyses through low-pass whole genome sequencing and whole exome sequencing. In order to estimate the independent association between genetic polymorphisms and ASCVD, a polygenic risk score (PRS) will be built through a hybrid approach including effect size of each Single Nucleotide Polymorphism (SNP), number of effect alleles observed, sample ploidy, total number of SNPs included in the PRS, and number of non-missing SNPs in the sample. In addition, the presence of pathogenic or likely pathogenic variants will be screened in 8 genes (ABCG5, ABCG8, APOB, APOE, LDLR, LDLRAP1, LIPA, PCSK9) associated with atherosclerosis. Multiple logistic regression will be applied to estimate adjusted odds ratios (OR) and 95% confidence intervals (CI), and population attributable risks will be calculated. Clinical trial registration: This study is registered in clinicaltrials.gov (NCT05515653).

Are both distinct epithelial and stromal cells molecular analysis from phyllodes tumors versus fibroadenoma components affected in breast fibroepithelial progression?

PURPOSE: To determine molecular events involved in the tumorigenesis of phyllodes tumors (PT) and the role of each stromal (SC) and epithelial (EC) cell. METHODS: Frozen breast samples enriched with epithelial and stromal cells from three fibroadenomas and 14 PT were retrieved and laser microdissected. Sanger and polymerase chain reaction-based sequencing of exon 2 MED12 and TERT promoter hotspot mutations were performed; 44K microarray platform was used to analyze gene expression. RESULTS: All three fibroadenomas (FAs) presented mutations in MED12, but not in TERT, whose mutation was observed in five of the 14 PTs. EC and SC of each affected tumor displayed identical alterations. Of the total differentially expressed genes (DEG) (EC = 1,543 and SC = 850), 984 were EC-eDEGs and 291 were SC-eDEGs. We found a high similarity of diseases and functions enriched by both cell types, but dissimilarity in the number of enriched canonical pathways. Three signaling canonical pathways overlapping with EC and SC were predicted to be activated in one cell type and inactivated in the other, while no overlap in eDEGs was assigned to them. We also identified 13 EC-eDEGs and five SC-eDEGs enriched networks, in which the SC-eDEGs were able to segregate FA from PT samples. CONCLUSIONS: Identical TERT mutations from both SC and ES origins might affect the PTs tumorigenesis. Gene expression differences suggest coordinated molecular processes between these components with determinant differences acquired by SC, able to fully distinguish PTs from FAs lesions.

Germline genetic variants in pheochromocytoma/paraganglioma: single-center experience.

Pheochromocytoma and paragangliomas (PPGLs) are rare neuroendocrine tumors carrying 25-40% pathogenic germline gene variants (PGVs). We evaluated clinical, laboratory, and germline molecular profile of 115 patients with pathologic (14 patients were relatives from 8 different families recruited for genetic survey) confirmed PPGL followed in our institution. Patients with classic MEN2A/MEN2B phenotypes and at-risk relatives underwent direct analysis of RET proto-oncogene, and the remaining had samples submitted to complete next-generation sequencing aiming 23 PPGL-related genes: ATM, ATR, CDKN2A, EGLN1, FH, HRAS, KIF1B, KMT2D, MAX, MDH2, MERTK, MET, NF1, PIK3CA, RET, SDHA, SDHAF2, SDHB, SDHC, SDHD, TMEM127, TP53, and VHL. We also developed a clinical judgment score (CJS) to determine the probability of patients having a potentially hereditary disease. The resulting genetic landscape showed that 67 patients (58.3%) had variants in at least one gene: 34 (50.7%) had exclusively pathogenic or likely pathogenic variants, 13 (19.4%) had pathogenic or likely pathogenic variants and variant of undetermined significance (VUS), and 20 (29.8%) carried only VUS. PGVs were found in RET (n = 18; 38.3%), VHL (n = 10; 21.3%), SDHB and NF1 (n = 8; 17% each), and MAX, SDHD, TMEM127, and TP53 (n = 1; 2.1% each). Direct genetic testing disclosed 91.3% sensitivity, 81.2% specificity, and 76.4% and 93.3% positive predictive value (PPV) and negative predictive values (NPV), respectively. The CJS to identify patients who would not benefit from genetic testing had 75% sensitivity, 96.4% specificity, and 60% and 98.2% PPV and NPV, respectively. In summary, the landscape of PPGL germline gene variants from 115 Brazilian patients resulted in slightly higher prevalent pathogenic and likely pathogenic variants, especially in the RET gene. We suggest a CJS to identify PPGL patients who would not require initial genetic evaluation, improving test specificity and reducing costs.

Mutational Portrait of Lung Adenocarcinoma in Brazilian Patients: Past, Present, and Future of Molecular Profiling in the Clinic.

Objectives: Approximately 60% of lung adenocarcinomas (LAs) carry mutations that can guide treatment with tyrosine-kinase inhibitors (TKI) and other targeted therapies. Data on activating mutations in EGFR and other tyrosine-kinase receptor (TKR) genes in highly admixed populations, such as that of Brazil, are scarce. In this study, we comprehensively analyzed the actionable alteration profile of LA in Brazilian patients. Materials and Methods: EGFR driver mutation data were collected from a large Brazilian LA cohort covering an 8-year period of molecular testing in a single institution. Tests were performed using three distinct methods, and demographic and histopathological data were analyzed. For a subset of patients, driver mutations in KRAS, NRAS, and BRAF and gene fusions involving TKR genes (before TKI treatment) and EGFR T790M (after TKI treatment) were assessed. Results: EGFR mutations were detected in 25% of 1,316 LAs evaluated, with exon 19 deletions and exon 21 L858R TKI sensitizing mutations representing 72.5% of all mutations. Mutation rates were higher in women and non-smokers (p < 0.001). Next-generation sequencing was very sensitive, with a lower rate of inconclusive results compared with Sanger sequencing and pyrosequencing. EGFR/RAS/BRAF hotspot gene panels were applied in 495 LA cases and detected oncogenic mutations in 51.3% of samples, most frequently in EGFR (22.4%) and KRAS (26.9%). In subgroups of 36 and 35 patients, gene fusions were detected in 11.1% of tumors and EGFR T790M resistance mutations were detected in 59% of plasma samples from patients previously treated with TKI, respectively. Conclusion: This report provides the first comprehensive actionable alteration portrait of LA in Brazil. The high rate of actionable alterations in EGFR and other driver genes in LA reinforces the need to incorporate TKI guided by molecular diagnostics into clinical routines for patients in both public and private healthcare systems.

Reciprocal changes in gene expression profiles of cocultured breast epithelial cells and primary fibroblasts.

The importance of epithelial-stroma interaction in normal breast development and tumor progression has been recognized. To identify genes that were regulated by these reciprocal interactions, we cocultured a nonmalignant (MCF10A) and a breast cancer derived (MDA-MB231) basal cell lines, with fibroblasts isolated from breast benign-disease adjacent tissues (NAF) or with carcinoma-associated fibroblasts (CAF), in a transwell system. Gene expression profiles of each coculture pair were compared with the correspondent monocultures, using a customized microarray. Contrariwise to large alterations in epithelial cells genomic profiles, fibroblasts were less affected. In MDA-MB231 highly represented genes downregulated by CAF derived factors coded for proteins important for the specificity of vectorial transport between ER and golgi, possibly affecting cell polarity whereas the response of MCF10A comprised an induction of genes coding for stress responsive proteins, representing a prosurvival effect. While NAF downregulated genes encoding proteins associated to glycolipid and fatty acid biosynthesis in MDA-MB231, potentially affecting membrane biogenesis, in MCF10A, genes critical for growth control and adhesion were altered. NAFs responded to coculture with MDA-MB231 by a decrease in the expression of genes induced by TGFbeta1 and associated to motility. However, there was little change in NAFs gene expression profile influenced by MCF10A. CAFs responded to the presence of both epithelial cells inducing genes implicated in cell proliferation. Our data indicate that interactions between breast fibroblasts and basal epithelial cells resulted in alterations in the genomic profiles of both cell types which may help to clarify some aspects of this heterotypic signaling.

Influence of BRCA1 Germline Mutations in the Somatic Mutational Burden of Triple-Negative Breast Cancer.

The majority of the hereditary triple-negative breast cancers (TNBCs) are associated with BRCA1 germline mutations. Nevertheless, the understanding of the role of BRCA1 deficiency in the TNBC tumorigenesis is poor. In this sense, we performed whole-exome sequencing of triplet samples (leucocyte, tumor, and normal-adjacent breast tissue) for 10 cases of early-onset TNBC, including 5 hereditary (with BRCA1 germline pathogenic mutation) and 5 sporadic (with no BRCA1 or BRCA2 germline pathogenic mutations), for assessing the somatic mutation repertoire. Protein-affecting somatic mutations were identified for both mammary tissues, and Ingenuity Pathway Analysis was used to investigate gene interactions. BRCA1 and RAD51C somatic promoter methylation in tumor samples was also investigated by bisulfite sequencing. Sporadic tumors had higher proportion of driver mutations (≥25% allele frequency) than BRCA1 hereditary tumors, whereas no difference was detected in the normal breast samples. Distinct gene networks were obtained from the driver genes in each group. The Cancer Genome Atlas data analysis of TNBC classified as hereditary and sporadic reinforced our findings. The data presented here indicate that in the absence of BRCA1 germline mutations, a higher number of driver mutations are required for tumor development and that different defective processes are operating in the tumorigenesis of hereditary and sporadic TNBC in young women.

Are both distinct epithelial and stromal cells molecular analysis from phyllodes tumors versus fibroadenoma components affected in breast fibroepithelial progression?

Purpose: To determine molecular events involved in the tumorigenesis of phyllodes tumors (PT) and the role of each stromal (SC) and epithelial (EC) cell. Methods: Frozen breast samples enriched with epithelial and stromal cells from three fibroadenomas and 14 PT were retrieved and laser microdissected. Sanger and polymerase chain reaction-based sequencing of exon 2 MED12 and TERT promoter hotspot mutations were performed; 44K microarray platform was used to analyze gene expression. Results: All three fibroadenomas (FAs) presented mutations in MED12, but not in TERT, whose mutation was observed in five of the 14 PTs. EC and SC of each affected tumor displayed identical alterations. Of the total differentially expressed genes (DEG) (EC = 1,543 and SC = 850), 984 were EC-eDEGs and 291 were SC-eDEGs. We found a high similarity of diseases and functions enriched by both cell types, but dissimilarity in the number of enriched canonical pathways. Three signaling canonical pathways overlapping with EC and SC were predicted to be activated in one cell type and inactivated in the other, while no overlap in eDEGs was assigned to them. We also identified 13 EC-eDEGs and five SC-eDEGs enriched networks, in which the SC-eDEGs were able to segregate FA from PT samples. Conclusions: Identical TERT mutations from both SC and ES origins might affect the PTs tumorigenesis. Gene expression differences suggest coordinated molecular processes between these components with determinant differences acquired by SC, able to fully distinguish PTs from FAs lesions.

Alternative splicing: a bioinformatics perspective.

The degree of diversity at the transcriptome and proteome levels generated by alternative splicing is astonishing. In this review, we discuss several issues related to alternative splicing with a special emphasis on identification strategies based on bioinformatics.

Epithelial cells captured from ductal carcinoma in situ reveal a gene expression signature associated with progression to invasive breast cancer.

Breast cancer biomarkers that can precisely predict the risk of progression of non-invasive ductal carcinoma in situ (DCIS) lesions to invasive disease are lacking. The identification of molecular alterations that occur during the invasion process is crucial for the discovery of drivers of transition to invasive disease and, consequently, biomarkers with clinical utility. In this study, we explored differences in gene expression in mammary epithelial cells before and after the morphological manifestation of invasion, i.e., early and late stages, respectively. In the early stage, epithelial cells were captured from both pre-invasive lesions with distinct malignant potential [pure DCIS as well as the in situ component that co-exists with invasive breast carcinoma lesions (DCIS-IBC)]; in the late stage, epithelial cells were captured from the two distinct morphological components of the same sample (in situ and invasive components). Candidate genes were identified using cDNA microarray and rapid subtractive hybridization (RaSH) cDNA libraries and validated by RT-qPCR assay using new samples from each group. These analyses revealed 26 genes, including 20 from the early and 6 from the late stage. The expression profile based on the 20 genes, marked by a preferential decrease in expression level towards invasive phenotype, discriminated the majority of DCIS samples. Thus, this study revealed a gene expression signature with the potential to predict DCIS progression and, consequently, provides opportunities to tailor treatments for DCIS patients.

Linear mRNA amplification approach for RNAseq from limited amount of RNA.

Whole-transcriptome evaluation by next-generation sequencing (NGS) has been widely applied in the investigation of diverse transcriptional scenarios. In many clinical situations, including needle biopsy samples or laser microdissected cells, limited amounts of RNA are usually available for the assessment of the whole transcriptome. Here, we describe an mRNA amplification protocol based on in vitro T7 transcription for transcriptome evaluation by NGS. Initially, we performed RNAseq from two human mammary epithelial cell lines and evaluated several aspects of the transcriptomes generated by linear amplification of Poly (A)(+) mRNA species, including transcript representation, variability and abundance. Our protocol showed to be efficient with respect to full-length transcript coverage and quantitative expression levels. We then evaluated the applicability of using this protocol in a more realistic research scenario, analyzing tumor tissue samples microdissected by laser capture. In order to increase the quantification power of the libraries only the 3' end of transcripts were sequenced. We found highly reproducible RNAseq data among amplified tumor samples, with a median Spearman's correlation of 80%, strongly suggesting that the amplification step and library protocol preparation lead to a consistent transcriptional profile. Altogether, we established a robust protocol for assessing the polyadenylated transcriptome derived from limited amounts of total RNA that is applicable to all NGS platforms.

Gene expression patterns through oral squamous cell carcinoma development: PD-L1 expression in primary tumor and circulating tumor cells.

Oral squamous cell carcinoma (OSCC) is the most common tumor of the oral cavity and has been associated with poor prognosis. Scarce prognostic markers are available for guiding treatment and/or sub-classifying patients. This study aims to identify biomarkers by searching for genes whose expression is increased or decreased during tumor progression (through T1 to T4 stages). Thirty-six samples from all tumor size stages (from T1 to T4) were analyzed using cDNA microarrays. Selected targets were analyzed by immunohistochemistry and in circulating tumor cells by immunofluorescence and Nanostring. Correlation was shown between PD-L1 and tumor size and lymph node metastasis, HOXB9 and tumor size, BLNK and perineural invasion, and between ZNF813 and perineural invasion. PD-L1 positivity was an independent prognostic factor in this cohort (p = 0.044, HH = 0.426). In CTCs from patients with locally advanced OSCC, we found a strong cytoplasmatic expression of PD-L1. PD-L1 is a ligand of PD-1 and is believed to limit T cell activity in inflammatory responses and limit autoimmune diseases. We demonstrated an important role for PD-L1 in primary tumors according to tumor size, and in disease specific survival. Therefore, we could further determine individuals with PD-L1+ CTCs, and possibly follow treatment using CTCs.