Beyond tobacco: genomic disparities in lung cancer between smokers and never-smokers.

BACKGROUND: Tobacco use is one of the main risk factors for Lung Cancer (LC) development. However, about 10-20% of those diagnosed with the disease are never-smokers. For Non-Small Cell Lung Cancer (NSCLC) there are clear differences in both the clinical presentation and the tumor genomic profiles between smokers and never-smokers. For example, the Lung Adenocarcinoma (LUAD) histological subtype in never-smokers is predominately found in young women of European, North American, and Asian descent. While the clinical presentation and tumor genomic profiles of smokers have been widely examined, never-smokers are usually underrepresented, especially those of a Latin American (LA) background. In this work, we characterize, for the first time, the difference in the genomic profiles between smokers and never-smokers LC patients from Chile. METHODS: We conduct a comparison by smoking status in the frequencies of genomic alterations (GAs) including somatic mutations and structural variants (fusions) in a total of 10 clinically relevant genes, including the eight most common actionable genes for LC (EGFR, KRAS, ALK, MET, BRAF, RET, ERBB2, and ROS1) and two established driver genes for malignancies other than LC (PIK3CA and MAP2K1). Study participants were grouped as either smokers (current and former, n = 473) or never-smokers (n = 200) according to self-report tobacco use at enrollment. RESULTS: Our findings indicate a higher overall GA frequency for never-smokers compared to smokers (58 vs. 45.7, p-value < 0.01) with the genes EGFR, KRAS, and PIK3CA displaying the highest prevalence while ERBB2, RET, and ROS1 the lowest. Never-smokers present higher frequencies in seven out of the 10 genes; however, smokers harbor a more complex genomic profile. The clearest differences between groups are seen for EGFR (15.6 vs. 21.5, p-value: < 0.01), PIK3CA (6.8 vs 9.5) and ALK (3.2 vs 7.5) in favor of never-smokers, and KRAS (16.3 vs. 11.5) and MAP2K1 (6.6 vs. 3.5) in favor of smokers. Alterations in these genes are comprised almost exclusively by somatic mutations in EGFR and mainly by fusions in ALK, and only by mutations in PIK3CA, KRAS and MAP2K1. CONCLUSIONS: We found clear differences in the genomic landscape by smoking status in LUAD patients from Chile, with potential implications for clinical management in these limited-resource settings.

Distinct Driver Pathway Enrichments and a High Prevalence of TSC2 Mutations in Right Colon Cancer in Chile: A Preliminary Comparative Analysis.

Colorectal cancer (CRC) is the second leading cause of cancer deaths globally. While ethnic differences in driver gene mutations have been documented, the South American population remains understudied at the genomic level, despite facing a rising burden of CRC. We analyzed tumors of 40 Chilean CRC patients (Chp) using next-generation sequencing and compared them to data from mainly Caucasian cohorts (TCGA and MSK-IMPACT). We identified 388 mutations in 96 out of 135 genes, with TP53 (45%), KRAS (30%), PIK3CA (22.5%), ATM (20%), and POLE (20%) being the most frequently mutated. TSC2 mutations were associated with right colon cancer (44.44% in RCRC vs. 6.45% in LCRC, p-value = 0.016), and overall frequency was higher compared to TCGA (p-value = 1.847 × 10-5) and MSK-IMPACT cohorts (p-value = 3.062 × 10-2). Limited sample size restricts definitive conclusions, but our data suggest potential differences in driver mutations for Chilean patients, being that the RTK-RAS oncogenic pathway is less affected and the PI3K pathway is more altered in Chp compared to TCGA (45% vs. 25.56%, respectively). The prevalence of actionable pathways and driver mutations can guide therapeutic choices, but can also impact treatment effectiveness. Thus, these findings warrant further investigation in larger Chilean cohorts to confirm these initial observations. Understanding population-specific driver mutations can guide the development of precision medicine programs for CRC patients.

ADAR1-mediated RNA-editing of 3'UTRs in breast cancer.

BACKGROUND: Whole transcriptome RNA variant analyses have shown that adenosine deaminases acting on RNA (ADAR) enzymes modify a large proportion of cellular RNAs, contributing to transcriptome diversity and cancer evolution. Despite the advances in the understanding of ADAR function in breast cancer, ADAR RNA editing functional consequences are not fully addressed. RESULTS: We characterized A to G(I) mRNA editing in 81 breast cell lines, showing increased editing at 3'UTR and exonic regions in breast cancer cells compared to immortalized non-malignant cell lines. In addition, tumors from the BRCA TCGA cohort show a 24% increase in editing over normal breast samples when looking at 571 well-characterized UTRs targeted by ADAR1. Basal-like subtype breast cancer patients with high level of ADAR1 mRNA expression shows a worse clinical outcome and increased editing in their 3'UTRs. Interestingly, editing was particularly increased in the 3'UTRs of ATM, GINS4 and POLH transcripts in tumors, which correlated with their mRNA expression. We confirmed the role of ADAR1 in this regulation using a shRNA in a breast cancer cell line (ZR-75-1). CONCLUSIONS: Altogether, these results revealed a significant association between the mRNA editing in genes related to cancer-relevant pathways and clinical outcomes, suggesting an important role of ADAR1 expression and function in breast cancer.

Concordance Analysis of ALK Gene Fusion Detection Methods in Patients with Non-Small-Cell Lung Cancer from Chile, Brazil, and Peru.

About 4% to 7% of the non-small-cell lung cancer patients have anaplastic lymphoma kinase (ALK) rearrangements, and specific targeted therapies improve patients' outcomes significantly. ALK gene fusions are detected by immunohistochemistry or fluorescent in situ hybridization as gold standards in South America. Next-generation sequencing-based assays are a reliable alternative, able to perform simultaneous detection of multiple events from a single sample. We analyzed 4240 non-small-cell lung cancer samples collected in 37 hospitals from Chile, Brazil, and Peru, where ALK rearrangements were determined as part of their standard of care (SofC) using either immunohistochemistry or fluorescent in situ hybridization. A subset of 1450 samples was sequenced with the Oncomine Focus Assay (OFA), and the concordance with the SofC tests was measured. An orthogonal analysis was performed using a real-time quantitative PCR echinoderm microtubule-associated protein-like 4-ALK fusion detection kit. ALK fusion prevalence is similar for Chile (3.67%; N = 2142), Brazil (4.05%; N = 1013), and Peru (4.59%; N = 675). Although a comparison between OFA and SofC assays showed similar sensitivity, OFA had significantly higher specificity and higher positive predictive value, which opens new opportunities for a more specific determination of ALK gene rearrangements.

webCEMiTool: Co-expression Modular Analysis Made Easy.

Co-expression analysis has been widely used to elucidate the functional architecture of genes under different biological processes. Such analysis, however, requires substantial knowledge about programming languages and/or bioinformatics skills. We present webCEMiTool, a unique online tool that performs comprehensive modular analyses in a fully automated manner. The webCEMiTool not only identifies co-expression gene modules but also performs several functional analyses on them. In addition, webCEMiTool integrates transcriptomic data with interactome information (i.e., protein-protein interactions) and identifies potential hubs on each network. The tool generates user-friendly html reports that allow users to search for specific genes in each module, as well as check if a module contains genes overrepresented in specific pathways or altered in a specific sample phenotype. We used webCEMiTool to perform a modular analysis of single-cell RNA-seq data of human cells infected with either Zika virus or dengue virus.

High Proportion of Potential Candidates for Immunotherapy in a Chilean Cohort of Gastric Cancer Patients: Results of the FORCE1 Study.

Gastric cancer (GC) is a heterogeneous disease. This heterogeneity applies not only to morphological and phenotypic features but also to geographical variations in incidence and mortality rates. As Chile has one of the highest mortality rates within South America, we sought to define a molecular profile of Chilean GCs (ClinicalTrials.gov identifier: NCT03158571/(FORCE1)). Solid tumor samples and clinical data were obtained from 224 patients, with subsets analyzed by tissue microarray (TMA; n = 90) and next generation sequencing (NGS; n = 101). Most demographic and clinical data were in line with previous reports. TMA data indicated that 60% of patients displayed potentially actionable alterations. Furthermore, 20.5% were categorized as having a high tumor mutational burden, and 13% possessed micro-satellite instability (MSI). Results also confirmed previous studies reporting high Epstein-Barr virus (EBV) positivity (13%) in Chilean-derived GC samples suggesting a high proportion of patients could benefit from immunotherapy. As expected, TP53 and PIK3CA were the most frequently altered genes. However, NGS demonstrated the presence of TP53, NRAS, and BRAF variants previously unreported in current GC databases. Finally, using the Kendall method, we report a significant correlation between EBV+ status and programmed death ligand-1 (PDL1)+ and an inverse correlation between p53 mutational status and MSI. Our results suggest that in this Chilean cohort, a high proportion of patients are potential candidates for immunotherapy treatment. To the best of our knowledge, this study is the first in South America to assess the prevalence of actionable targets and to examine a molecular profile of GC patients.

Long Non-Coding RNAs Responsive to Salt and Boron Stress in the Hyper-Arid Lluteño Maize from Atacama Desert.

Long non-coding RNAs (lncRNAs) have been defined as transcripts longer than 200 nucleotides, which lack significant protein coding potential and possess critical roles in diverse cellular processes. Long non-coding RNAs have recently been functionally characterized in plant stress-response mechanisms. In the present study, we perform a comprehensive identification of lncRNAs in response to combined stress induced by salinity and excess of boron in the Lluteño maize, a tolerant maize landrace from Atacama Desert, Chile. We use deep RNA sequencing to identify a set of 48,345 different lncRNAs, of which 28,012 (58.1%) are conserved with other maize (B73, Mo17 or Palomero), with the remaining 41.9% belonging to potentially Lluteño exclusive lncRNA transcripts. According to B73 maize reference genome sequence, most Lluteño lncRNAs correspond to intergenic transcripts. Interestingly, Lluteño lncRNAs presents an unusual overall higher expression compared to protein coding genes under exposure to stressed conditions. In total, we identified 1710 putatively responsive to the combined stressed conditions of salt and boron exposure. We also identified a set of 848 stress responsive potential trans natural antisense transcripts (trans-NAT) lncRNAs, which seems to be regulating genes associated with regulation of transcription, response to stress, response to abiotic stimulus and participating of the nicotianamine metabolic process. Reverse transcription-quantitative PCR (RT-qPCR) experiments were performed in a subset of lncRNAs, validating their existence and expression patterns. Our results suggest that a diverse set of maize lncRNAs from leaves and roots is responsive to combined salt and boron stress, being the first effort to identify lncRNAs from a maize landrace adapted to extreme conditions such as the Atacama Desert. The information generated is a starting point to understand the genomic adaptabilities suffered by this maize to surpass this extremely stressed environment.

ADAR1-mediated RNA-editing of 3'UTRs in breast cancer

BACKGROUND: Whole transcriptome RNA variant analyses have shown that adenosine deaminases acting on RNA ( ADAR ) enzymes modify a large proportion of cellular RNAs, contributing to transcriptome diversity and cancer evolution. Despite the advances in the understanding of ADAR function in breast cancer, ADAR RNA editing functional consequences are not fully addressed. RESULTS: We characterized A to G(I) mRNA editing in 81 breast cell lines, showing increased editing at 3'UTR and exonic regions in breast cancer cells compared to immortalized non-malignant cell lines. In addition, tumors from the BRCA TCGA cohort show a 24% increase in editing over normal breast samples when looking at 571 well-characterized UTRs targeted by ADAR1. Basal-like subtype breast cancer patients with high level of ADAR1 mRNA expression shows a worse clinical outcome and increased editing in their 3'UTRs. Interestingly, editing was particularly increased in the 3'UTRs of ATM, GINS4 and POLH transcripts in tumors, which correlated with their mRNA expression. We confirmed the role of ADAR1 in this regulation using a shRNA in a breast cancer cell line (ZR-75-1). CONCLUSIONS: Altogether, these results revealed a significant association between the mRNA editing in genes related to cancer-relevant pathways and clinical outcomes, suggesting an important role of ADAR1 expression and function in breast cancer.