Next generation sequencing in cancer: opportunities and challenges for precision cancer medicine

Over the past decade, testing the genes of patients and their specific cancer types has become standardized practice in medical oncology since somatic mutations, changes in gene expression and epigenetic modifications are all hallmarks of cancer. However, while cancer genetic assessment has been limited to single biomarkers to guide the use of therapies, improvements in nucleic acid sequencing technologies and implementation of different genome analysis tools have enabled clinicians to detect these genomic alterations and identify functional and disease-associated genomic variants. Next-generation sequencing (NGS) technologies have provided clues about therapeutic targets and genomic markers for novel clinical applications when standard therapy has failed. While Sanger sequencing, an accurate and sensitive approach, allows for the identification of potential novel variants, it is however limited by the single amplicon being interrogated. Similarly, quantitative and qualitative profiling of gene expression changes also represents a challenge for the cancer field. Both RT-PCR and microarrays are efficient approaches, but are limited to the genes present on the array or being assayed. This leaves vast swaths of the transcriptome, including non-coding RNAs and other features, unexplored. With the advent of the ability to collect and analyze genomic sequence data in a timely fashion and at an ever-decreasing cost, many of these limitations have been overcome and are being incorporated into cancer research and diagnostics giving patients and clinicians new hope for targeted and personalized treatment. Below we highlight the various applications of next-generation sequencing in precision cancer medicine

IsomiR Expression Profiles in Human Lymphoblastoid Cell Lines Exhibit Population and Gender Dependencies.

For many years it was believed that each mature microRNA (miRNA) existed as a single entity with fixed endpoints and a \u27static\u27 and unchangeable primary sequence. However, recent evidence suggests that mature miRNAs are more \u27dynamic\u27 and that each miRNA precursor arm gives rise to multiple isoforms, the isomiRs. Here we report on our identification of numerous and abundant isomiRs in the lymphoblastoid cell lines (LCLs) of 452 men and women from five different population groups. Unexpectedly, we find that these isomiRs exhibit an expression profile that is population-dependent and gender-dependent. This is important as it indicates that the LCLs of each gender/population combination have their own unique collection of mature miRNA transcripts. Moreover, each identified isomiR has its own characteristic abundance that remains consistent across biological replicates indicating that these are not degradation products. The primary sequences of identified isomiRs differ from the known miRBase miRNA either at their 5´-endpoint (leads to a different \u27seed\u27 sequence and suggests a different targetome), their 3´-endpoint, or both simultaneously. Our analysis of Argonaute PAR-CLIP data from LCLs supports the association of many of these newly identified isomiRs with the Argonaute silencing complex and thus their functional roles through participation in the RNA interference pathway

Knowledge about the presence or absence of miRNA isoforms (isomiRs) can successfully discriminate amongst 32 TCGA cancer types.

Isoforms of human miRNAs (isomiRs) are constitutively expressed with tissue- and disease-subtype-dependencies. We studied 10 271 tumor datasets from The Cancer Genome Atlas (TCGA) to evaluate whether isomiRs can distinguish amongst 32 TCGA cancers. Unlike previous approaches, we built a classifier that relied solely on \u27binarized\u27 isomiR profiles: each isomiR is simply labeled as \u27present\u27 or \u27absent\u27. The resulting classifier successfully labeled tumor datasets with an average sensitivity of 90% and a false discovery rate (FDR) of 3%, surpassing the performance of expression-based classification. The classifier maintained its power even after a 15× reduction in the number of isomiRs that were used for training. Notably, the classifier could correctly predict the cancer type in non-TCGA datasets from diverse platforms. Our analysis revealed that the most discriminatory isomiRs happen to also be differentially expressed between normal tissue and cancer. Even so, we find that these highly discriminating isomiRs have not been attracting the most research attention in the literature. Given their ability to successfully classify datasets from 32 cancers, isomiRs and our resulting \u27Pan-cancer Atlas\u27 of isomiR expression could serve as a suitable framework to explore novel cancer biomarkers

A transcription map of the 6p22.3 reading disability locus identifying candidate genes

BACKGROUND: Reading disability (RD) is a common syndrome with a large genetic component. Chromosome 6 has been identified in several linkage studies as playing a significant role. A more recent study identified a peak of transmission disequilibrium to marker JA04 (G72384) on chromosome 6p22.3, suggesting that a gene is located near this marker. RESULTS: In silico cloning was used to identify possible candidate genes located near the JA04 marker. The 2 million base pairs of sequence surrounding JA04 was downloaded and searched against the dbEST database to identify ESTs. In total, 623 ESTs from 80 different tissues were identified and assembled into 153 putative coding regions from 19 genes and 2 pseudogenes encoded near JA04. The identified genes were tested for their tissue specific expression by RT-PCR. CONCLUSIONS: In total, five possible candidate genes for RD and other diseases mapping to this region were identified

Churchill regulates cell movement and mesoderm specification by repressing Nodal signaling

<p>Abstract</p> <p>Background</p> <p>Cell movements are essential to the determination of cell fates during development. The zinc-finger transcription factor, Churchill (ChCh) has been proposed to regulate cell fate by regulating cell movements during gastrulation in the chick. However, the mechanism of action of ChCh is not understood.</p> <p>Results</p> <p>We demonstrate that ChCh acts to repress the response to Nodal-related signals in zebrafish. When ChCh function is abrogated the expression of mesodermal markers is enhanced while ectodermal markers are expressed at decreased levels. In cell transplant assays, we observed that ChCh-deficient cells are more motile than wild-type cells. When placed in wild-type hosts, ChCh-deficient cells often leave the epiblast, migrate to the germ ring and are later found in mesodermal structures. We demonstrate that both movement of ChCh-compromised cells to the germ ring and acquisition of mesodermal character depend on the ability of the donor cells to respond to Nodal signals. Blocking Nodal signaling in the donor cells at the levels of Oep, Alk receptors or Fast1 inhibited migration to the germ ring and mesodermal fate change in the donor cells. We also detect additional unusual movements of transplanted ChCh-deficient cells which suggests that movement and acquisition of mesodermal character can be uncoupled. Finally, we demonstrate that ChCh is required to limit the transcriptional response to Nodal.</p> <p>Conclusion</p> <p>These data establish a broad role for ChCh in regulating both cell movement and Nodal signaling during early zebrafish development. We show that <it>chch </it>is required to limit mesodermal gene expression, inhibit Nodal-dependant movement of presumptive ectodermal cells and repress the transcriptional response to Nodal signaling. These findings reveal a dynamic role for <it>chch </it>in regulating cell movement and fate during early development.</p

MINTbase v2.0: a comprehensive database for tRNA-derived fragments that includes nuclear and mitochondrial fragments from all The Cancer Genome Atlas projects.

MINTbase is a repository that comprises nuclear and mitochondrial tRNA-derived fragments (\u27tRFs\u27) found in multiple human tissues. The original version of MINTbase comprised tRFs obtained from 768 transcriptomic datasets. We used our deterministic and exhaustive tRF mining pipeline to process all of The Cancer Genome Atlas datasets (TCGA). We identified 23 413 tRFs with abundance of ≥ 1.0 reads-per-million (RPM). To facilitate further studies of tRFs by the community, we just released version 2.0 of MINTbase that contains information about 26 531 distinct human tRFs from 11 719 human datasets as of October 2017. Key new elements include: the ability to filter tRFs on-the-fly by minimum abundance thresholding; the ability to filter tRFs by tissue keywords; easy access to information about a tRF\u27s maximum abundance and the datasets that contain it; the ability to generate relative abundance plots for tRFs across cancer types and convert them into embeddable figures; MODOMICS information about modifications of the parental tRNA, etc. Version 2.0 of MINTbase contains 15x more datasets and nearly 4x more distinct tRFs than the original version, yet continues to offer fast, interactive access to its contents. Version 2.0 is available freely at http://cm.jefferson.edu/MINTbase/

Expression of Tryptophan 2,3-Dioxygenase in Metastatic Uveal Melanoma

Uveal melanoma (UM) is the most common primary eye malignancy in adults and up to 50% of patients subsequently develop systemic metastasis. Metastatic uveal melanoma (MUM) is highly resistant to immunotherapy. One of the mechanisms for resistance would be the immune-suppressive tumor microenvironment. Here, we have investigated the role of tryptophan 2,3-dioxygenase (TDO) in UM. Both TDO and indoleamine 2,3-dioxygenase (IDO) catalyze tryptophan and produce kynurenine, which could cause inhibition of T cell immune responses. We first studied the expression of TDO on tumor tissue specimens obtained from UM hepatic metastasis. High expression of TDO protein was confirmed in all hepatic metastasis. TDO was positive in both normal hepatocytes and the tumor cells with relatively higher expression in tumor cells. On the other hand, IDO protein remained undetectable in all of the MUM specimens. UM cell lines established from metastasis also expressed TDO protein and increasing kynurenine levels were detected in the supernatant of MUM cell culture. In TCGA database, higher TDO2 expression in primary UM significantly correlated to BAP1 mutation and monosomy 3. These results indicate that TDO might be one of the key mechanisms for resistance to immunotherapy in UM

Beyond the one-locus-one-miRNA paradigm: microRNA isoforms enable deeper insights into breast cancer heterogeneity.

Here we describe our study of miRNA isoforms (isomiRs) in breast cancer (BRCA) and normal breast data sets from the Cancer Genome Atlas (TCGA) repository. We report that the full isomiR profiles, from both known and novel human-specific miRNA loci, are particularly rich in information and can distinguish tumor from normal tissue much better than the archetype miRNAs. IsomiR expression is also dependent on the patient\u27s race, exemplified by miR-183-5p, several isomiRs of which are upregulated in triple negative BRCA in white but not black women. Additionally, we find that an isomiR\u27s 5\u27 endpoint and length, but not the genomic origin, are key determinants of the regulation of its expression. Overexpression of distinct miR-183-5p isomiRs in MDA-MB-231 cells followed by microarray analysis revealed that each isomiR has a distinct impact on the cellular transcriptome. Parallel integrative analysis of mRNA expression from BRCA data sets of the TCGA repository demonstrated that isomiRs can distinguish between the luminal A and luminal B subtypes and explain in more depth the molecular differences between them than the archetype molecules. In conclusion, our findings provide evidence that post-transcriptional studies of BRCA will benefit from transcending the one-locus-one-miRNA paradigm and taking into account all isoforms from each miRNA locus as well as the patient\u27s race

Recessive mutation in tetraspanin CD151 causes Kindler syndrome-like epidermolysis bullosa with multi-systemic manifestations including nephropathy

Epidermolysis bullosa (EB) is caused by mutations in as many as 19 distinct genes. We have developed a next-generation sequencing (NGS) panel targeting genes known to be mutated in skin fragility disorders, including tetraspanin CD151 expressed in keratinocytes at the dermal-epidermal junction. The NGS panel was applied to a cohort of 92 consanguineous families of unknown subtype of EB. In one family, a homozygous donor splice site mutation in CD151 (NM_139029; c.351 + 2T > C) at the exon 5/intron 5 border was identified, and RT-PCR and whole transcriptome analysis by RNA-seq confirmed deletion of the entire exon 5 encoding 25 amino acids. Immunofluorescence of proband's skin and Western blot of skin proteins with a monoclonal antibody revealed complete absence of CD151. Transmission electron microscopy showed intracellular disruption and cell-cell dysadhesion of keratinocytes in the lower epidermis. Clinical examination of the 33-year old proband, initially diagnosed as Kindler syndrome, revealed widespread blistering, particularly on pretibial areas, poikiloderma, nail dystrophy, loss of teeth, early onset alopecia, and esophageal webbing and strictures. The patient also had history of nephropathy with proteinuria. Collectively, the results suggest that biallelic loss-of-function mutations in CD151 underlie an autosomal recessive mechano-bullous disease with systemic features. Thus, CD151 should be considered as the 20th causative, EB-associated gene

Murine MPDZ-linked hydrocephalus is caused by hyperpermeability of the choroid plexus.

Though congenital hydrocephalus is heritable, it has been linked only to eight genes, one of which is MPDZ Humans and mice that carry a truncated version of MPDZ incur severe hydrocephalus resulting in acute morbidity and lethality. We show by magnetic resonance imaging that contrast medium penetrates into the brain ventricles of mice carrying a Mpdz loss-of-function mutation, whereas none is detected in the ventricles of normal mice, implying that the permeability of the choroid plexus epithelial cell monolayer is abnormally high. Comparative proteomic analysis of the cerebrospinal fluid of normal and hydrocephalic mice revealed up to a 53-fold increase in protein concentration, suggesting that transcytosis through the choroid plexus epithelial cells of Mpdz KO mice is substantially higher than in normal mice. These conclusions are supported by ultrastructural evidence, and by immunohistochemistry and cytology data. Our results provide a straightforward and concise explanation for the pathophysiology of Mpdz-linked hydrocephalus