Clonal evolution of esophageal squamous cell carcinoma from normal mucosa to primary tumor and metastases.

The clonal evolution which drives esophageal squamous cell carcinoma (ESCC) from initiation in normal cell to primary carcinoma and metastases is poorly understood. In this study, multi-region whole-exome sequencing (WES) (284X) and whole-genome single nucleotide polymorphism genotyping were performed on a total of 109 samples of ESCC from 10 patients. This included 42 apparently normal samples of esophageal mucosa at increasing distances from the upper or lower boundaries of the primary tumor to the surgical margins of resection, 43 spatially separated tissue samples within primary tumor and 24 regional lymph node metastases. Phylogenetic analysis was performed to reconstruct ancestor-descendant relationships of clones and the clonal composition of multi-region samples. Mutations of cancer-related genes were validated by deep targeted sequencing (1,168X). Both inter- and intra-tumoral genetic heterogeneity were obvious across multi-region samples among ESCC patients. Clones varying in number from one to seven were discovered within each regional tumor or metastatic sample. Phylogenetic analysis demonstrated complex clonal evolution patterns. Regional lymph node metastases had characteristics of early initiation and polyclonal spreading, and could be derived from carcinoma in situ (CIS) directly. TP53 was the only gene harboring non-silent mutations identified across all multi-region tumor samples of all ten patients. Mutations of TP53 were also found in histologically normal mucosa in sites away from primary tumor.

Genome Evolution Analysis of Recurrent Testicular Malignant Mesothelioma by Whole-Genome Sequencing.

BACKGROUND/AIMS: Malignant mesothelioma of the tunica vaginalis testis is a rare and lethal disease. The genomic characteristics and genetic changes of tumor cells during the progression of this disease are unknown. METHODS: we performed whole-genome sequencing of four successive tumor samples derived from surgery and a blood sample in a single patient. RESULTS: All tumors were found to have significant C-to-T and T-to-C mutations, and amplification of copy number in chromosomes 1 and 12 were notified in all tumor samples. Subclone analysis revealed a parallel evolution of the tumor in this patient. We also identified some mutations in mesothelioma-associated genes such as KIF25, AHNAK, and PRDM2. CONCLUSIONS: The results showed a comprehensive genomic change in malignant mesothelioma of the tunica vaginalis testis and provide a better understanding of the clonal evolution during tumor recurrence and metastasis.

The Complete Chloroplast Genome Sequences of the Medicinal Plant Forsythia suspensa (Oleaceae).

Forsythia suspensa is an important medicinal plant and traditionally applied for the treatment of inflammation, pyrexia, gonorrhea, diabetes, and so on. However, there is limited sequence and genomic information available for F. suspensa. Here, we produced the complete chloroplast genomes of F. suspensa using Illumina sequencing technology. F. suspensa is the first sequenced member within the genus Forsythia (Oleaceae). The gene order and organization of the chloroplast genome of F. suspensa are similar to other Oleaceae chloroplast genomes. The F. suspensa chloroplast genome is 156,404 bp in length, exhibits a conserved quadripartite structure with a large single-copy (LSC; 87,159 bp) region, and a small single-copy (SSC; 17,811 bp) region interspersed between inverted repeat (IRa/b; 25,717 bp) regions. A total of 114 unique genes were annotated, including 80 protein-coding genes, 30 tRNA, and four rRNA. The low GC content (37.8%) and codon usage bias for A- or T-ending codons may largely affect gene codon usage. Sequence analysis identified a total of 26 forward repeats, 23 palindrome repeats with lengths >30 bp (identity > 90%), and 54 simple sequence repeats (SSRs) with an average rate of 0.35 SSRs/kb. We predicted 52 RNA editing sites in the chloroplast of F. suspensa, all for C-to-U transitions. IR expansion or contraction and the divergent regions were analyzed among several species including the reported F. suspensa in this study. Phylogenetic analysis based on whole-plastome revealed that F. suspensa, as a member of the Oleaceae family, diverged relatively early from Lamiales. This study will contribute to strengthening medicinal resource conservation, molecular phylogenetic, and genetic engineering research investigations of this species.

Genomic discordances and heterogeneous mutational burden, PD-L1 expression and immune infiltrates of non-small cell lung cancer metastasis.

AIMS: To investigate the genomic discordances and heterogeneous mutational burden, PD-L1 expression and immune cell (IC) infiltrates of non-small cell lung cancer (NSCLC) metastasis. METHODS: Surgical samples from 41 cases of NSCLC with metastatic tumours (MTs) and paired primary tumours (PTs) were collected. PD-L1 expression and ICs were quantified using image-based immunohistochemistry profiling. Whole exome sequencing was employed to explore discrepancies in genomic characteristics, tumour mutational burden (TMB) and tumour neoantigen burden (TNB) in 28 cases. RESULTS: Non-synonymous mutations in MTs were slightly more than in PTs, with only 42.34% of mutations shared between paired PTs and MTs. The heterogeneity of TMB showed no significant difference (p=0.785) between MTs and PTs, while TNB significantly increased in MTs (p=0.013). MTs generally exhibited a higher density of PD-L1+ cells and a higher tumour proportion score with a lower density of IC infiltrates. Subgroup analysis considering clinicopathological factors revealed that the heterogeneity of immune biomarkers was closely associated with the histology of lung adenocarcinoma, metastatic sites of extrapulmonary, time intervals and treatment history. Prognosis analysis indicated that a high density of CD8+ T cells was a low-risk factor, whereas a high density of PD-L1+ cells in MTs was a high-risk factor for cancer-related death in metastatic NSCLC. CONCLUSIONS: The mutational burden, PD-L1 expression and IC infiltrates undergo changes during NSCLC metastasis, which may impact the immunotherapeutic benefits in patients with NSCLC with metastatic progression and should be monitored according to clinical scenarios.

Single-cell sequencing reveals CD133+CD44--originating evolution and novel stemness related variants in human colorectal cancer.

BACKGROUND: Tumor heterogeneity of human colorectal cancer (CRC)-initiating cells (CRCICs) in cancer tissues often represents aggressive features of cancer progression. For high-resolution examination of CRCICs, we performed single-cell whole-exome sequencing (scWES) and bulk cell targeted exome sequencing (TES) of CRCICs to investigate stemness-specific somatic alterations or clonal evolution. METHODS: Single cells of three subpopulations of CRCICs (CD133+CD44+, CD133-CD44+, and CD133+CD44- cells), CRC cells (CRCCs), and control cells from one CRC tissue were sorted for scWES. Then, we set up a mutation panel from scWES data and TES was used to validate mutation distribution and clonal evolution in additional 96 samples (20 patients) those were also sorted into the same three groups of CRCICs and CRCCs. The knock-down experiments were used to analyze stemness-related mutant genes. Neoantigens of these mutant genes and their MHC binding affinity were also analyzed. FINDINGS: Clonal evolution analysis of scWES and TES showed that the CD133+CD44- CRCICs were the likely origin of CRC before evolving into other groups of CRCICs/CRCCs. We revealed that AHNAK2, PLIN4, HLA-B, ALK, CCDC92 and ALMS1 genes were specifically mutated in CRCICs followed by the validation of their functions. Furthermore, four predicted neoantigens of AHNAK2 were identified and validated, which might have applications in immunotherapy for CRC patients. INTERPRETATION: All the integrative analyses above revealed clonal evolution of CRC and new markers for CRCICs and demonstrate the important roles of CRCICs in tumorigenesis and progression of CRCs. FUNDING: A full list of funding bodies that contributed to this study can be found in the Acknowledgements section.

Genomic alterations caused by HPV integration in a cohort of Chinese endocervical adenocarcinomas.

The association between human papillomavirus (HPV) integration and relevant genomic changes in uterine cervical adenocarcinoma is poorly understood. This study is to depict the genomic mutational landscape in a cohort of 20 patients. HPV+ and HPV- groups were defined as patients with and without HPV integration in the host genome. The genetic changes between these two groups were described and compared by whole-genome sequencing (WGS) and whole-exome sequencing (WES). WGS identified 2916 copy number variations and 743 structural variations. WES identified 6113 somatic mutations, with a mutational burden of 2.4 mutations/Mb. Six genes were predicted as driver genes: PIK3CA, KRAS, TRAPPC12, NDN, GOLGA6L4 and BAIAP3. PIK3CA, NDN, GOLGA6L4, and BAIAP3 were recognized as significantly mutated genes (SMGs). HPV was detected in 95% (19/20) of patients with cervical adenocarcinoma, 7 of whom (36.8%) had HPV integration (HPV+ group). In total, 1036 genes with somatic mutations were confirmed in the HPV+ group, while 289 genes with somatic mutations were confirmed in the group without HPV integration (HPV- group); only 2.1% were shared between the two groups. In the HPV+ group, GOLGA6L4 and BAIAP3 were confirmed as SMGs, while PIK3CA, NDN, KRAS, FUT1, and GOLGA6L64 were identified in the HPV- group. ZDHHC3, PKD1P1, and TGIF2 showed copy number amplifications after HPV integration. In addition, the HPV+ group had significantly more neoantigens. HPV integration rather than HPV infection results in different genomic changes in cervical adenocarcinoma.

Mutations of METTL3 predict response to neoadjuvant chemotherapy in muscle-invasive bladder cancer.

BACKGROUND AND AIM: Neoadjuvant chemotherapy (NAC) followed by radical cystectomy is the current gold standard treatment for muscle-invasive urothelial bladder cancer (MIBC). Nonetheless, some MIBC patients showed limited pathological response after NAC. Herein, we used whole-exome sequencing (WES) to identify genetic mutations in MIBC that can predict NAC response. METHODS: Forty MIBC patients were enrolled in this study, in which 33 were successfully examined by WES and Sanger sequencing in the discovery cohort (n=13) and the validation cohort (n=20), respectively. ANNOVAR software was used to identify the potential mutations based on the data of WES. In addition, tumor-specific somatic mutations including single nucleotide variants and indels were called with the muTECT and Strelka software. The mutational analysis of specific genes was carried out based on the data from cBioPortal for Cancer Genomics. RESULTS: In the discovery cohort, the mutation frequencies of TP53, MED16, DRC7, CEND1, ATAD5, SETD8, and PIK3CA were significantly higher in 13 MIBC patients. Specifically, the presence of somatic mutations of APC, ATM, CDH9, CTNNB1, METTL3, NBEAL1, PTPRH, RNASEL, and FBXW7 in NAC responder signifies that these mutations were potential predictors of pathological response to NAC. Furthermore, somatic mutations of CCDC141, PIK3CA, CHD5, GPR149, MUC20, TSC1, and USP54 were exclusively identified in NAC nonresponders, suggesting that these mutations may participate in the process of NAC resistance. In the validation cohort, the somatic mutations of CDH9, METTL3, and PTPRH were significantly enriched in NAC responders while the somatic mutation of CCDC141 was significantly enriched in NAC nonresponders. Furthermore, survival analysis revealed that the patients expressing mutated METTL3 have a longer overall survival and disease- or progression-free survival than the patients acquiring wild-type METTL3. CONCLUSION: The somatic mutation of METTL3 can be a potential predictive biomarker of NAC response in MIBC patients. RELEVANCE FOR PATIENTS: MIBC patients bearing mutated METTL3 display a pathological response to NAC and have a significantly longer overall survival or disease/progression-free survival as compared to the patients bearing wild-type METTL3. Thus, the somatic mutation of METTL3 is a potential biomarker for predicting response to NAC in MIBC patients, assisting doctors in making the clinical decision.

Genomic profiling in amyloid light-chain amyloidosis reveals mutation profiles associated with overall survival.

Background: Amyloid light chain (AL) amyloidosis is characterized by tissue deposition of amyloid fibres derived from immunoglobulin that can lead to irreversible organ damage. Information about genomic profiles of AL amyloidosis is lacking.Methods: In this study, we adopted a two-step strategy to investigate the mutational profile of AL amyloidosis bone marrow plasma cells (PCs) and their clinical implications. In step one, whole-exome sequencing was performed in bone marrow PCs and paired with normal tissue from 10 AL amyloidosis patients, by which we identified 10 significantly mutated genes (SMGs). In step two, we constituted a targeted gene sequencing (TGS) panel covering the frequently mutated genes identified in step one, genes reported in prior AL amyloidosis studies, and known cancer driver mutations. Then, we analysed an expanded cohort of AL amyloidosis patients (N = 48) with this panel comprising 98 genes.Results: Four recurrent mutations were identified by TGS and verified by Sanger sequencing: ASB15 (c. 844 C > T), ASCC3 (c. 1595 A > G), HIST1H1E (c. 311 C > T) and KRAS (c. 35 G > A), among which the first three mutations were associated with inferior overall survival (OS). Additionally, we found that the number of mutations identified by the TGS panel of 98 genes could be a prognostic predictor for OS.Conclusions: In summary, we revealed genomic profiling in AL amyloidosis and found mutation profiles associated with OS.

Mutation landscape and intra-tumor heterogeneity of two MANECs of the esophagus revealed by multi-region sequencing.

Mixed adenoneuroendocrine carcinoma (MANEC) in the esophagus is an infrequent but highly malignant cancer with few known genomic alterations. We conducted whole-exome sequencing and whole-genome SNP genotyping for 4-6 tumor subregions and 5-6 adjacent normal tissue sites and 1-3 lymph node metastases in two esophageal MANECs to detect somatic mutations and copy number alterations, and to explore their spatial heterogeneity and underlying clonal structure. TP53 mutation, RB1 deletion or LOH, and PIK3CA, PTEN, KRAS, SOX2, DVL3, TP63 amplification appeared in all regions in both tumors. Mutations falling in known cancer genes tended to show higher variant allele frequencies than those not falling in these genes in at least one of the cases. Phylogenetic analyses of the samples and underlying subclones suggested extensive migration across different tumor regions and from some regions to the lymph nodes. Lymph node metastases appeared to have been seeded by both early founder cells as well as subsequent, locally emerging daughter clones. A phenotypically normal tissue site carried most of the mutations found in neighboring tumor samples, implying field cancerization. Understanding such complex genetic heterogeneity within each patient will be important for guiding clinical decisions.

Intraspecific and heteroplasmic variations, gene losses and inversions in the chloroplast genome of Astragalus membranaceus.

Astragalus membranaceus is an important medicinal plant in Asia. Several of its varieties have been used interchangeably as raw materials for commercial production. High resolution genetic markers are in urgent need to distinguish these varieties. Here, we sequenced and analyzed the chloroplast genome of A. membranaceus (Fisch.) Bunge var. mongholicus (Bunge) P.K. Hsiao using the next generation DNA sequencing technology. The genome was assembled using Abyss and then subjected to gene prediction using CPGAVAS and repeat analysis using MISA, Tandem Repeats Finder, and REPuter. Finally, the genome was subjected phylogenetic and comparative genomic analyses. The complete genome is 123,582 bp long, containing only one copy of the inverted repeat. Gene prediction revealed 110 genes encoding 76 proteins, 30 tRNAs, and four rRNAs. Five intra-specific hypermutation loci were identified, three of which are heteroplasmic. Furthermore, three gene losses and two large inversions were identified. Comparative genomic analyses demonstrated the dynamic nature of the Papilionoideae chloroplast genomes, which showed occurrence of numerous hypermutation loci, frequent gene losses, and fragment inversions. Results obtained herein elucidate the complex evolutionary history of chloroplast genomes and have laid the foundation for the identification of genetic markers to distinguish A. membranaceus varieties.

Complete plastid genome of Astragalus membranaceus (Fisch.) Bunge var. membranaceus.

The complete nucleotide sequence of the Astragalus membranaceus (Fisch.) Bunge var. membranaceus chloroplast genome was reported and characterized in this study. The chloroplast genome is a circular molecule of 123623 bp that belongs to the inverted repeat-lacking clade (IRLC). It comprises 110 genes, including 76 protein-coding genes, four unique rRNAs and 30 tRNAs. Similar to the plastomes of A. membranaceus (Fisch.) Bunge var. mongholicus (Bunge) P. K. Hsiao and other closely related species, rpl22 and rps16 are absent. The phylogenetic analysis of 67 proteins from 29 chloroplast genomes belonging to IRLC provided strong support for the non-monophyly of Galegeae. This genome has provided a wealth of information for distinguishing varieties of A. membranaceus.

Comprehensive analysis of alternative splicing in Digitalis purpurea by strand-specific RNA-Seq.

Digitalis purpurea (D. purpurea) is one of the most important medicinal plants and is well known in the treatment of heart failure because of the cardiac glycosides that are its main active compounds. However, in the absence of strand specific sequencing information, the post-transcriptional mechanism of gene regulation in D. purpurea thus far remains unknown. In this study, a strand-specific RNA-Seq library was constructed and sequenced using Illumina HiSeq platforms to characterize the transcriptome of D. purpurea with a focus on alternative splicing (AS) events and the effect of AS on protein domains. De novo RNA-Seq assembly resulted in 48,475 genes. Based on the assembled transcripts, we reported a list of 3,265 AS genes, including 5,408 AS events in D. purpurea. Interestingly, both glycosyltransferases and monooxygenase, which were involved in the biosynthesis of cardiac glycosides, are regulated by AS. A total of 2,422 AS events occurred in coding regions, and 959 AS events were located in the regions of 882 unique protein domains, which could affect protein function. This D. purpurea transcriptome study substantially increased the expressed sequence resource and presented a better understanding of post-transcriptional regulation to further facilitate the medicinal applications of D. purpurea for human health.

Aberrant miRNA expression response to UV irradiation in human liver cancer cells.

microRNAs (miRNAs) are small, highly conserved, non-coding RNAs that regulate gene expression at the post­transcriptional level. The expression of these small RNA genes is tightly regulated during development, differentiation and apoptosis of normal cells, however, they are often deregulated in various types of cancer. miRNA expression is also affected by cellular stress, including radiation and chemotherapy. The present study monitored the expression levels of several miRNAs by using the stem-loop real-time polymerase chain reaction (PCR) in HepG2 cells following ultraviolet (UV) radiation. Our data demonstrated that UV irradiation is able to induce alterations in miRNA expression levels in HepG2 cells. Among them, miR-26a, miR-34a and miR-146a were significantly upregulated, while the expression of miR-21 was significantly downregulated. Bioinformatics analysis of these significantly regulated miRNAs was discussed. The results also indicated that miRNAs may be part of the innate response mechanism of the cells to radiation injury, which provides a rationale for miRNA replacement therapy by using specific miRNAs that may function as tumor suppressor genes in several types of cancer.

Preparation of heparin-immobilized PVA and its adsorption for low-density lipoprotein from hyperlipemia plasma.

In this study, heparin was covalently coupled by glutaraldehyde to Poly(vinyl alcohol) [PVA] in solid-liquid two-phase reaction system by two-step synthesis method to prepare a LDL-selective adsorbent. The parameters (the material ratio, reaction time and dosage of catalyzer) were investigated to evaluate their effect upon the immobilized amount of heparin onto the surface of PVA, IR was used to verify the covalent immobilization result and the heparin-modified PVA was also undergone the evaluation of its adsorption capability for low-density lipoprotein from hyperlipemia plasma, and its hemocompatibility was preliminarily evaluated by platelet adhesion test. Results showed: (1) under optimized reaction conditions the highest immobilization amount of heparin onto PVA surface within the experiments of this study has been obtained; (2) the optimized reaction conditions were: (i) at the refluxing temperature 78 degrees C; (ii) the material ratio of "PVA(g): 50% glutaraldehyde (ml)" was about "1:3"; (iii) the reaction time was about 5 h; and (iv) the amount of catalyzer (concentrated HCL) was about 1% of the 50% glutaraldehyde; (3) within the experiments of this study the highest immobilization amount would be up to 25 microg heparin on the surface of per g PVA granules; (4) the heparin-modified PVA granules showed significant adsorption for LDL under faintly alkaline environment (pH=7.2-9.5) ; (5) The result of platelet adhesion test showed no platelet adhered to its surface. Therefore, immobilization of heparin onto the surface of a support is one approach to prepare a kind of LDL adsorbent for blood purification.

[ICE expression in laryngeal carcinoma and it's clinical relationship].

OBJECTIVE: To study the role of ICE in Laryngeal Carcinoma. METHOD: Using polyclonal antibody for ICE gene protein, 73 cases of laryngeal carcinoma, 30 cases dysplasia and 5 cases normal were stained by SP immunohistochemistry. RESULT: ICE protein is situated mainly in the upper of epithelium and a little in the basal layer. Compared with the normal epithelium, the expression of ICE in the dysplasia is decreased significantly (P < 0.01). There is a relationship between expression of ICE in laryngeal carcinoma and the degree of differentiation. The better differentiation, the much expression of ICE (P < 0.05). CONCLUSION: There is an intimate relationship between innormal expression protein of ICE and dysplasia, and that of between laryngeal carcinoma and prognosis. The expression of ICE may be a monitoring marker for prognosis of laryngeal carcinoma.