Inherited BRCA1 and RNF43 pathogenic variants in a familial colorectal cancer type X family.

Genetic susceptibility to familial colorectal cancer (CRC), including for individuals classified as Familial Colorectal Cancer Type X (FCCTX), remains poorly understood. We describe a multi-generation CRC-affected family segregating pathogenic variants in both BRCA1, a gene associated with breast and ovarian cancer and RNF43, a gene associated with Serrated Polyposis Syndrome (SPS). A single family out of 105 families meeting the criteria for FCCTX (Amsterdam I family history criteria with mismatch repair (MMR)-proficient CRCs) recruited to the Australasian Colorectal Cancer Family Registry (ACCFR; 1998-2008) that underwent whole exome sequencing (WES), was selected for further testing. CRC and polyp tissue from four carriers were molecularly characterized including a single CRC that underwent WES to determine tumor mutational signatures and loss of heterozygosity (LOH) events. Ten carriers of a germline pathogenic variant BRCA1:c.2681_2682delAA p.Lys894ThrfsTer8 and eight carriers of a germline pathogenic variant RNF43:c.988 C > T p.Arg330Ter were identified in this family. Seven members carried both variants, four of which developed CRC. A single carrier of the RNF43 variant met the 2019 World Health Organization (WHO2019) criteria for SPS, developing a BRAF p.V600 wildtype CRC. Loss of the wildtype allele for both BRCA1 and RNF43 variants was observed in three CRC tumors while a LOH event across chromosome 17q encompassing both genes was observed in a CRC. Tumor mutational signature analysis identified the homologous recombination deficiency (HRD)-associated COSMIC signatures SBS3 and ID6 in a CRC for a carrier of both variants. Our findings show digenic inheritance of pathogenic variants in BRCA1 and RNF43 segregating with CRC in a FCCTX family. LOH and evidence of BRCA1-associated HRD supports the importance of both these tumor suppressor genes in CRC tumorigenesis.

Genome-wide exploration reveals distinctive northern and southern variants of Clonorchis sinensis in the Far East.

Clonorchis sinensis is a carcinogenic liver fluke that causes clonorchiasis-a neglected tropical disease (NTD) affecting ~35 million people worldwide. No vaccine is available, and chemotherapy relies on one anthelmintic, praziquantel. This parasite has a complex life history and is known to infect a range of species of intermediate (freshwater snails and fish) and definitive (piscivorous) hosts. Despite this biological complexity and the impact of this biocarcinogenic pathogen, there has been no previous study of molecular variation in this parasite on a genome-wide scale. Here, we conducted the first extensive nuclear genomic exploration of C. sinensis individuals (n = 152) representing five distinct populations from mainland China, and one from Far East Russia, and revealed marked genetic variation within this species between "northern" and "southern" geographical regions. The discovery of this variation indicates the existence of biologically distinct variants within C. sinensis, which may have distinct epidemiology, pathogenicity and/or chemotherapic responsiveness. The detection of high heterozygosity within C. sinensis specimens suggests that this parasite has developed mechanisms to readily adapt to changing environments and/or host species during its life history/evolution. From an applied perspective, the identification of invariable genes could assist in finding new intervention targets in this parasite, given the major clinical relevance of clonorchiasis. From a technical perspective, the genomic-informatic workflow established herein will be readily applicable to a wide range of other parasites that cause NTDs.

Minimizing Sample Failure Rates for Challenging Clinical Tumor Samples.

Identification of somatic variants in cancer by high-throughput sequencing has become common clinical practice, largely because many of these variants may be predictive biomarkers for targeted therapies. However, there can be high sample quality control (QC) failure rates for some tests that prevent the return of results. Stem-loop inhibition mediated amplification (SLIMamp) is a patented technology that has been incorporated into commercially available cancer next-generation sequencing testing kits. The claimed advantage is that these kits can interrogate challenging formalin-fixed, paraffin-embedded tissue samples with low tumor purity, poor-quality DNA, and/or low-input DNA, resulting in a high sample QC pass rate. The study aimed to substantiate that claim using Pillar Biosciences oncoReveal Solid Tumor Panel. Forty-eight samples that had failed one or more preanalytical QC sample parameters for whole-exome sequencing from the Australian Translational Genomics Centre's ISO15189-accredited diagnostic genomics laboratory were acquired. XING Genomic Services performed an exploratory data analysis to characterize the samples and then tested the samples in their ISO15189-accredited laboratory. Clinical reports could be generated for 37 (77%) samples, of which 29 (60%) contained clinically actionable or significant variants that would not otherwise have been identified. Eleven samples were deemed unreportable, and the sequencing data were likely dominated by artifacts. A novel postsequencing QC metric was developed that can discriminate between clinically reportable and unreportable samples.

Subtype-Specific Analyses Reveal Infiltrative Basal Cell Carcinomas Are Highly Interactive with their Environment.

Little is known regarding the molecular differences between basal cell carcinoma (BCC) subtypes, despite clearly distinct phenotypes and clinical outcomes. In particular, infiltrative BCCs have poorer clinical outcomes in terms of response to therapy and propensity for dissemination. In this project, we aimed to use exome sequencing and RNA sequencing to identify somatic mutations and molecular pathways leading to infiltrative BCCs. Using whole-exome sequencing of 36 BCC samples (eight infiltrative) combined with previously reported exome data (58 samples), we determine that infiltrative BCCs do not contain a distinct somatic variant profile and carry classical UV-induced mutational signatures. RNA sequencing on both datasets revealed key differentially expressed genes, such as POSTN and WISP1, suggesting increased integrin and Wnt signaling. Immunostaining for periostin and WISP1 clearly distinguished infiltrative BCCs, and nuclear ß-catenin staining patterns further validated the resulting increase in Wnt signaling in infiltrative BCCs. Of significant interest, in BCCs with mixed morphology, infiltrative areas expressed WISP1, whereas nodular areas did not, supporting a continuum between subtypes. In conclusion, infiltrative BCCs do not differ in their genomic alteration in terms of initiating mutations. They display a specific type of interaction with the extracellular matrix environment regulating Wnt signaling.

Marked mitochondrial genetic variation in individuals and populations of the carcinogenic liver fluke Clonorchis sinensis.

Clonorchiasis is a neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, and is often associated with a malignant form of bile duct cancer (cholangiocarcinoma). Although some aspects of the epidemiology of clonorchiasis are understood, little is known about the genetics of C. sinensis populations. Here, we conducted a comprehensive genetic exploration of C. sinensis from endemic geographic regions using complete mitochondrial protein gene sets. Genomic DNA samples from C. sinensis individuals (n = 183) collected from cats and dogs in China (provinces of Guangdong, Guangxi, Hunan, Heilongjiang and Jilin) as well as from rats infected with metacercariae from cyprinid fish from the Russian Far East (Primorsky Krai region) were deep sequenced using the BGISEQ-500 platform. Informatic analyses of mitochondrial protein gene data sets revealed marked genetic variation within C. sinensis; significant variation was identified within and among individual worms from distinct geographical locations. No clear affiliation with a particular location or host species was evident, suggesting a high rate of dispersal of the parasite across endemic regions. The present work provides a foundation for future biological, epidemiological and ecological studies using mitochondrial protein gene data sets, which could aid in elucidating associations between particular C. sinensis genotypes/haplotypes and the pathogenesis or severity of clonorchiasis and its complications (including cholangiocarcinoma) in humans.

Regional Variation in Epidermal Susceptibility to UV-Induced Carcinogenesis Reflects Proliferative Activity of Epidermal Progenitors.

To better understand the influence of ultraviolet (UV) irradiation on the initial steps of skin carcinogenesis, we examine patches of labeled keratinocytes as a proxy for clones in the interfollicular epidermis (IFE) and measure their size variation upon UVB irradiation. Multicolor lineage tracing reveals that in chronically irradiated skin, patches near hair follicles (HFs) increase in size, whereas those far from follicles do not change. This is explained by proliferation of basal epidermal cells within 60 µm of HF openings. Upon interruption of UVB, patch size near HFs regresses significantly. These anatomical differences in proliferative behavior have significant consequences for the cell of origin of basal cell carcinomas (BCCs). Indeed, a UV-inducible murine BCC model shows that BCC patches are more frequent, larger, and more invasive near HFs. These findings have major implications for the prevention of field cancerization in the epidermis.

Erratum: PTEN deletion drives acute myeloid leukemia resistance to MEK inhibitors.

[This corrects the article DOI: 10.18632/oncotarget.27206.].

Chromosome arm aneuploidies shape tumour evolution and drug response.

Chromosome arm aneuploidies (CAAs) are pervasive in cancers. However, how they affect cancer development, prognosis and treatment remains largely unknown. Here, we analyse CAA profiles of 23,427 tumours, identifying aspects of tumour evolution including probable orders in which CAAs occur and CAAs predicting tissue-specific metastasis. Both haematological and solid cancers initially gain chromosome arms, while only solid cancers subsequently preferentially lose multiple arms. 72 CAAs and 88 synergistically co-occurring CAA pairs multivariately predict good or poor survival for 58% of 6977 patients, with negligible impact of whole-genome doubling. Additionally, machine learning identifies 31 CAAs that robustly alter response to 56 chemotherapeutic drugs across cell lines representing 17 cancer types. We also uncover 1024 potential synthetic lethal pharmacogenomic interactions. Notably, in predicting drug response, CAAs substantially outperform  mutations and focal deletions/amplifications combined. Thus, CAAs predict cancer prognosis, shape tumour evolution, metastasis and drug response, and may advance precision oncology.

PTEN deletion drives acute myeloid leukemia resistance to MEK inhibitors.

Kinases such as MEK are attractive targets for novel therapy in cancer, including acute myeloid leukaemia (AML). Acquired and inherent resistance to kinase inhibitors, however, is becoming an increasingly important challenge for the clinical success of such therapeutics, and often arises from mutations in the drug-binding domain of the target kinase. To identify possible causes of resistance to MEK inhibition, we generated a model of resistance by long-term treatment of AML cells with AZD6244 (selumetinib). Remarkably, resistance to MEK inhibition was due to acquired PTEN haploinsufficiency, rather than mutation of MEK. Resistance via this mechanism was confirmed using CRISPR/Cas9 technology targeting exon 5 of PTEN. While PTEN loss has been previously implicated in resistance to a number of other therapeutic agents, this is the first time that it has been shown directly and in AML.

Progression of Disease Within 24 Months in Follicular Lymphoma Is Associated With Reduced Intratumoral Immune Infiltration.

PURPOSE: Understanding the immunobiology of the 15% to 30% of patients with follicular lymphoma (FL) who experience progression of disease within 24 months (POD24) remains a priority. Solid tumors with low levels of intratumoral immune infiltration have inferior outcomes. It is unknown whether a similar relationship exists between POD24 in FL. PATIENTS AND METHODS: Digital gene expression using a custom code set-five immune effector, six immune checkpoint, one macrophage molecules-was applied to a discovery cohort of patients with early- and advanced-stage FL (n = 132). T-cell receptor repertoire analysis, flow cytometry, multispectral immunofluorescence, and next-generation sequencing were performed. The immune infiltration profile was validated in two independent cohorts of patients with advanced-stage FL requiring systemic treatment (n = 138, rituximab plus cyclophosphamide, vincristine, prednisone; n = 45, rituximab plus cyclophosphamide, doxorubicin, vincristine, and prednisone), with the latter selected to permit comparison of patients experiencing a POD24 event with those having no progression at 5 years or more. RESULTS: Immune molecules showed distinct clustering, characterized by either high or low expression regardless of categorization as an immune effector, immune checkpoint, or macrophage molecule. Low programmed death-ligand 2 (PD-L2) was the most sensitive/specific marker to segregate patients with adverse outcomes; therefore, PD-L2 expression was chosen to distinguish immune infiltrationHI (ie, high PD-L2) FL biopsies from immune infiltrationLO (ie, low PD-L2) tumors. Immune infiltrationHI tissues were highly infiltrated with macrophages and expanded populations of T-cell clones. Of note, the immune infiltrationLO subset of patients with FL was enriched for POD24 events (odds ratio [OR], 4.32; c-statistic, 0.81; P = .001), validated in the independent cohorts (rituximab plus cyclophosphamide, vincristine, prednisone: OR, 2.95; c-statistic, 0.75; P = .011; and rituximab plus cyclophosphamide, doxorubicin, vincristine, and prednisone: OR, 7.09; c-statistic, 0.88; P = .011). Mutations were equally proportioned across tissues, which indicated that degree of immune infiltration is capturing aspects of FL biology distinct from its mutational profile. CONCLUSION: Assessment of immune-infiltration by PD-L2 expression is a promising tool with which to help identify patients who are at risk for POD24.

Deconvolution of single-cell multi-omics layers reveals regulatory heterogeneity.

Integrative analysis of multi-omics layers at single cell level is critical for accurate dissection of cell-to-cell variation within certain cell populations. Here we report scCAT-seq, a technique for simultaneously assaying chromatin accessibility and the transcriptome within the same single cell. We show that the combined single cell signatures enable accurate construction of regulatory relationships between cis-regulatory elements and the target genes at single-cell resolution, providing a new dimension of features that helps direct discovery of regulatory patterns specific to distinct cell identities. Moreover, we generate the first single cell integrated map of chromatin accessibility and transcriptome in early embryos and demonstrate the robustness of scCAT-seq in the precise dissection of master transcription factors in cells of distinct states. The ability to obtain these two layers of omics data will help provide more accurate definitions of "single cell state" and enable the deconvolution of regulatory heterogeneity from complex cell populations.

Subclonal evolution in disease progression from MGUS/SMM to multiple myeloma is characterised by clonal stability.

Multiple myeloma (MM) is a largely incurable haematological malignancy defined by the clonal proliferation of malignant plasma cells (PCs) within the bone marrow. Clonal heterogeneity has recently been established as a feature in MM, however, the subclonal evolution associated with disease progression has not been described. Here, we performed whole-exome sequencing of serial samples from 10 patients, providing new insights into the progression from monoclonal gammopathy of undetermined significance (MGUS) and smouldering MM (SMM), to symptomatic MM. We confirm that intraclonal genetic heterogeneity is a common feature at diagnosis and that the driving events involved in disease progression are more subtle than previously reported. We reveal that MM evolution is mainly characterised by the phenomenon of clonal stability, where the transformed subclonal PC populations identified at MM are already present in the asymptomatic MGUS/SMM stages. Our findings highlight the possibility that PC extrinsic factors may play a role in subclonal evolution and MGUS/SMM to MM progression.

Single-cell RNA-seq reveals dynamic transcriptome profiling in human early neural differentiation.

Background: Investigating cell fate decision and subpopulation specification in the context of the neural lineage is fundamental to understanding neurogenesis and neurodegenerative diseases. The differentiation process of neural-tube-like rosettes in vitro is representative of neural tube structures, which are composed of radially organized, columnar epithelial cells and give rise to functional neural cells. However, the underlying regulatory network of cell fate commitment during early neural differentiation remains elusive. Results: In this study, we investigated the genome-wide transcriptome profile of single cells from six consecutive reprogramming and neural differentiation time points and identified cellular subpopulations present at each differentiation stage. Based on the inferred reconstructed trajectory and the characteristics of subpopulations contributing the most toward commitment to the central nervous system lineage at each stage during differentiation, we identified putative novel transcription factors in regulating neural differentiation. In addition, we dissected the dynamics of chromatin accessibility at the neural differentiation stages and revealed active cis-regulatory elements for transcription factors known to have a key role in neural differentiation as well as for those that we suggest are also involved. Further, communication network analysis demonstrated that cellular interactions most frequently occurred in the embryoid body stage and that each cell subpopulation possessed a distinctive spectrum of ligands and receptors associated with neural differentiation that could reflect the identity of each subpopulation. Conclusions: Our study provides a comprehensive and integrative study of the transcriptomics and epigenetics of human early neural differentiation, which paves the way for a deeper understanding of the regulatory mechanisms driving the differentiation of the neural lineage.

Erythrocyte microRNA sequencing reveals differential expression in relapsing-remitting multiple sclerosis.

BACKGROUND: There is a paucity of knowledge concerning erythrocytes in the aetiology of Multiple Sclerosis (MS) despite their potential to contribute to disease through impaired antioxidant capacity and altered haemorheological features. Several studies have identified an abundance of erythrocyte miRNAs and variable profiles associated with disease states, such as sickle cell disease and malaria. The aim of this study was to compare the erythrocyte miRNA profile of relapsing-remitting MS (RRMS) patients to healthy sex- and age-matched controls. METHODS: Erythrocytes were purified by density-gradient centrifugation and RNA was extracted. Following library preparation, samples were run on a HiSeq4000 Illumina instrument (paired-end 100 bp sequencing). Sequenced erythrocyte miRNA profiles (9 patients and 9 controls) were analysed by DESeq2. Differentially expressed miRNAs were validated by RT-qPCR using miR-152-3p as an endogenous control and replicated in a larger cohort (20 patients and 18 controls). After logarithmic transformation, differential expression was determined by two-tailed unpaired t-tests. Logistic regression analysis was carried out and receiver operating characteristic (ROC) curves were generated to determine biomarker potential. RESULTS: A total of 236 erythrocyte miRNAs were identified. Of twelve differentially expressed miRNAs in RRMS two showed increased expression (adj. p < 0.05). Only modest fold-changes were evident across differentially expressed miRNAs. RT-qPCR confirmed differential expression of miR-30b-5p (0.61 fold, p < 0.05) and miR-3200-3p (0.36 fold, p < 0.01) in RRMS compared to healthy controls. Relative expression of miR-3200-5p (0.66 fold, NS p = 0.096) also approached significance. MiR-3200-5p was positively correlated with cognition measured by audio-recorded cognitive screen (r = 0.60; p < 0.01). MiR-3200-3p showed greatest biomarker potential as a single miRNA (accuracy = 75.5%, p < 0.01, sensitivity = 72.7%, specificity = 84.0%). Combining miR-3200-3p, miR-3200-5p, and miR-30b-5p into a composite biomarker increased accuracy to 83.0% (p < 0.05), sensitivity to 77.3%, and specificity to 88.0%. CONCLUSIONS: This is the first study to report differences in erythrocyte miRNAs in RRMS. While the role of miRNAs in erythrocytes remains to be elucidated, differential expression of erythrocyte miRNAs may be exploited as biomarkers and their potential contribution to MS pathology and cognition should be further investigated.

Translocation Breakpoints Preferentially Occur in Euchromatin and Acrocentric Chromosomes.

Chromosomal translocations drive the development of many hematological and some solid cancers. Several factors have been identified to explain the non-random occurrence of translocation breakpoints in the genome. These include chromatin density, gene density and CCCTC-binding factor (CTCF)/cohesin binding site density. However, such factors are at least partially interdependent. Using 13,844 and 1563 karyotypes from human blood and solid cancers, respectively, our multiple regression analysis only identified chromatin density as the primary statistically significant predictor. Specifically, translocation breakpoints preferentially occur in open chromatin. Also, blood and solid tumors show markedly distinct translocation signatures. Strikingly, translocation breakpoints occur significantly more frequently in acrocentric chromosomes than in non-acrocentric chromosomes. Thus, translocations are probably often generated around nucleoli in the inner nucleoplasm, away from the nuclear envelope. Importantly, our findings remain true both in multivariate analyses and after removal of highly recurrent translocations. Finally, we applied pairwise probabilistic co-occurrence modeling. In addition to well-known highly prevalent translocations, such as those resulting in BCR-ABL1 (BCR-ABL) and RUNX1-RUNX1T1 (AML1-ETO) fusion genes, we identified significantly underrepresented translocations with putative fusion genes, which are probably subject to strong negative selection during tumor evolution. Taken together, our findings provide novel insights into the generation and selection of translocations during cancer development.

Targeted Next-Generation Sequencing for Detecting MLL Gene Fusions in Leukemia.

Mixed lineage leukemia (MLL) gene rearrangements characterize approximately 70% of infant and 10% of adult and therapy-related leukemia. Conventional clinical diagnostics, including cytogenetics and fluorescence in situ hybridization (FISH) fail to detect MLL translocation partner genes (TPG) in many patients. Long-distance inverse (LDI)-PCR, the "gold standard" technique that is used to characterize MLL breakpoints, is laborious and requires a large input of genomic DNA (gDNA). To overcome the limitations of current techniques, a targeted next-generation sequencing (NGS) approach that requires low RNA input was tested. Anchored multiplex PCR-based enrichment (AMP-E) was used to rapidly identify a broad range of MLL fusions in patient specimens. Libraries generated using Archer FusionPlex Heme and Myeloid panels were sequenced using the Illumina platform. Diagnostic specimens (n = 39) from pediatric leukemia patients were tested with AMP-E and validated by LDI-PCR. In concordance with LDI-PCR, the AMP-E method successfully identified TPGs without prior knowledge. AMP-E identified 10 different MLL fusions in the 39 samples. Only two specimens were discordant; AMP-E successfully identified a MLL-MLLT1 fusion where LDI-PCR had failed to determine the breakpoint, whereas a MLL-MLLT3 fusion was not detected by AMP-E due to low expression of the fusion transcript. Sensitivity assays demonstrated that AMP-E can detect MLL-AFF1 in MV4-11 cell dilutions of 10-7 and transcripts down to 0.005 copies/ng.Implications: This study demonstrates a NGS methodology with improved sensitivity compared with current diagnostic methods for MLL-rearranged leukemia. Furthermore, this assay rapidly and reliably identifies MLL partner genes and patient-specific fusion sequences that could be used for monitoring minimal residual disease. Mol Cancer Res; 16(2); 279-85. ©2017 AACR.

Corrigendum: Whole-genome landscape of pancreatic neuroendocrine tumours.

This corrects the article DOI: 10.1038/nature21063.

Cutting edge genomics reveal new insights into tumour development, disease progression and therapeutic impacts in multiple myeloma.

Multiple Myeloma (MM) is a haematological malignancy characterised by the clonal expansion of plasma cells (PCs) within the bone marrow. Despite advances in therapy, MM remains a largely incurable disease with a median survival of 6 years. In almost all cases, the development of MM is preceded by the benign PC condition Monoclonal Gammopathy of Undetermined Significance (MGUS). Recent studies show that the transformation of MGUS to MM is associated with complex genetic changes. Understanding how these changes contribute to evolution will present targets for clinical intervention. We discuss three models of MM evolution; the linear, the expansionist and the intraclonal heterogeneity models. Of particular interest is the intraclonal heterogeneity model. Here, distinct populations of MM PCs carry differing combinations of genetic mutations. Acquisition of additional mutations can contribute to subclonal lineages where "driver" mutations may influence selective pressure and dominance, and "passenger" mutations are neutral in their effects. Furthermore, studies show that clinical intervention introduces additional selective pressure on tumour cells and can influence subclone survival, leading to therapy resistance. This review discusses how Next Generation Sequencing approaches are revealing critical insights into the genetics of MM development, disease progression and treatment. MM disease progression will illuminate possible mechanisms underlying the tumour.