Rahul Kanwar

Somatic mutation of CDKN1B in small intestine neuroendocrine tumors

The diagnosed incidence of small intestine neuroendocrine tumors (SI-NETs) is increasing, and the underlying genomic mechanisms have not yet been defined. Using exome- and genome-sequence analysis of SI-NETs, we identified recurrent somatic mutations and deletions in CDKN1B, the cyclin-dependent kinase inhibitor gene, which encodes p27. We observed frameshift mutations of CDKN1B in 14 of 180 SI-NETs, and we detected hemizygous deletions encompassing CDKN1B in 7 out of 50 SI-NETs, nominating p27 as a tumor suppressor and implicating cell cycle dysregulation in the etiology of SI-NETs.

The genomic landscape of small intestine neuroendocrine tumors

Small intestine neuroendocrine tumors (SI-NETs) are the most common malignancy of the small bowel. Several clinical trials target PI3K/Akt/mTOR signaling; however, it is unknown whether these or other genes are genetically altered in these tumors. To address the underlying genetics, we analyzed 48 SI-NETs by massively parallel exome sequencing. We detected an average of 0.1 somatic single nucleotide variants (SNVs) per 106 nucleotides (range, 0-0.59), mostly transitions (C>T and A>G), which suggests that SI-NETs are stable cancers. 197 protein-altering somatic SNVs affected a preponderance of cancer genes, including FGFR2, MEN1, HOOK3, EZH2, MLF1, CARD11, VHL, NONO, and SMAD1. Integrative analysis of SNVs and somatic copy number variations identified recurrently altered mechanisms of carcinogenesis: chromatin remodeling, DNA damage, apoptosis, RAS signaling, and axon guidance. Candidate therapeutically relevant alterations were found in 35 patients, including SRC, SMAD family genes, AURKA, EGFR, HSP90, and PDGFR. Mutually exclusive amplification of AKT1 or AKT2 was the most common event in the 16 patients with alterations of PI3K/Akt/mTOR signaling. We conclude that sequencing-based analysis may provide provisional grouping of SI-NETs by therapeutic targets or deregulated pathways.

Sequencing of Charcot-Marie-Tooth disease genes in a toxic polyneuropathy.

Mutations in Charcot–Marie–Tooth disease (CMT) genes are the cause of rare familial forms of polyneuropathy. Whether allelic variability in CMT genes is also associated with common forms of polyneuropathy—considered “acquired” in medical parlance—is unknown. Chemotherapy‐induced peripheral neuropathy (CIPN) occurs commonly in cancer patients and is individually unpredictable. We used CIPN as a clinical model to investigate the association of non‐CMT polyneuropathy with CMT genes.

SAAP-RRBS

Summary: Reduced representation bisulfite sequencing (RRBS) is a cost-effective approach for genome-wide methylation pattern profiling. Analyzing RRBS sequencing data is challenging and specialized alignment/mapping programs are needed. Although such programs have been developed, a comprehensive solution that provides researchers with good quality and analyzable data is still lacking. To address this need, we have developed a Streamlined Analysis and Annotation Pipeline for RRBS data (SAAP-RRBS) that integrates read quality assessment/clean-up, alignment, methylation data extraction, annotation, reporting and visualization. This package facilitates a rapid transition from sequencing reads to a fully annotated CpG methylation report to biological interpretation. Availability and implementation: SAAP-RRBS is freely available to non-commercial users at the web site http://ndc.mayo.edu/mayo/research/biostat/stand-alone-packages.cfm. Contact: baheti.saurabh@mayo.edu or sun.zhifu@mayo.edu Supplementary Information: Supplementary data are available at Bioinformatics online.

Exome-level comparison of primary well-differentiated neuroendocrine tumors and their cell lines

Neuroendocrine cancer cell lines are used to investigate therapeutic targets in neuroendocrine tumors (NET) and have been instrumental in the design of clinical trials targeting the PI3K/AKT/mTOR pathways, VEGF inhibitors, and somatostatin analogues. It remains unknown, however, whether the genomic makeup of NET cell lines reflect that of primary NET since comprehensive unbiased genome sequencing has not been performed on the cell lines. Four bronchopulmonary NET (BP-NET)-NCI-H720, NCI-H727, NCI-H835, and UMC11-and two pancreatic neuroendocrine tumors (panNET)-BON-1 and QGP1-were cultured. DNA was isolated, and exome sequencing was done. GATK and EXCAVATOR were used for bioinformatic analysis. We detected a total of 1,764 nonsynonymous single nucleotide variants at a rate of 8 per Mb in BP-NET and 4.3 per Mb in panNET cell lines, including 52 mutated COSMIC cancer genes in these cell lines, such as TP53, BRCA1, RB1, TSC2, NOTCH1, EP300, GNAS, KDR, STK11, and APC but not ATRX, DAXX, nor MEN1. Our data suggest that mutation rate, the pattern of copy number variations, and the mutational spectra in the BP-NET cell lines are more similar to the changes observed in small cell lung cancer than those found in primary BP-NET. Likewise, mutation rate and pattern including the absence of mutations in ATRX/DAXX, MEN1, and YY1 in the panNET cell lines BON1 and QGP1 suggest that these cell lines do not have the genetic signatures of a primary panNET. These results suggest that results from experiments with BP-NET and panNET cell lines need to be interpreted with caution.

SERE: Single-parameter quality control and sample comparison for RNA-Seq

BackgroundAssessing the reliability of experimental replicates (or global alterations corresponding to different experimental conditions) is a critical step in analyzing RNA-Seq data. Pearson’s correlation coefficient r has been widely used in the RNA-Seq field even though its statistical characteristics may be poorly suited to the task.ResultsHere we present a single-parameter test procedure for count data, the Simple Error Ratio Estimate (SERE), that can determine whether two RNA-Seq libraries are faithful replicates or globally different. Benchmarking shows that the interpretation of SERE is unambiguous regardless of the total read count or the range of expression differences among bins (exons or genes), a score of 1 indicating faithful replication (i.e., samples are affected only by Poisson variation of individual counts), a score of 0 indicating data duplication, and scores >1 corresponding to true global differences between RNA-Seq libraries. On the contrary the interpretation of Pearson’s r is generally ambiguous and highly dependent on sequencing depth and the range of expression levels inherent to the sample (difference between lowest and highest bin count). Cohen’s simple Kappa results are also ambiguous and are highly dependent on the choice of bins. For quantifying global sample differences SERE performs similarly to a measure based on the negative binomial distribution yet is simpler to compute.ConclusionsSERE can therefore serve as a straightforward and reliable statistical procedure for the global assessment of pairs or large groups of RNA-Seq datasets by a single statistical parameter.

Diametrically opposite methylome-transcriptome relationships in high- and low-CpG promoter genes in postmitotic neural rat tissue.

DNA methylation can control some CpG-poor genes but unbiased studies have not found a consistent genome-wide association with gene activity outside of CpG islands or shores possibly due to use of cell lines or limited bioinformatics analyses. We performed reduced representation bisulfite sequencing (RRBS) of rat dorsal root ganglia encompassing postmitotic primary sensory neurons (n = 5, r > 0.99; orthogonal validation p < 10−19). The rat genome suggested a dichotomy of genes previously reported in other mammals: low CpG content (< 3.2%) promoter (LCP) genes and high CpG content (≥ 3.2%) promoter (HCP) genes. A genome-wide integrated methylome-transcriptome analysis showed that LCP genes were markedly hypermethylated when repressed, and hypomethylated when active with a 40% difference in a broad region at the 5′ of the transcription start site (p < 10−87 for -6000 bp to -2000 bp, p < 10−73 for -2000 bp to +2000 bp, no difference in gene body p = 0.42). HCP genes had minimal TSS-associated methylation regardless of transcription status, but gene body methylation appeared to be lost in repressed HCP genes. Therefore, diametrically opposite methylome-transcriptome associations characterize LCP and HCP genes in postmitotic neural tissue in vivo.

Testing of candidate single nucleotide variants associated with paclitaxel neuropathy in the trial NCCTG N08C1 (Alliance)

Paclitaxel‐induced peripheral neuropathy (PIPN) cannot be predicted from clinical parameters and might have a pharmacogenomic basis. Previous studies identified single nucleotide variants (SNV) associated with PIPN. However, only a subset of findings has been confirmed to date in more than one study, suggesting a need for further re‐testing and validation in additional clinical cohorts. Candidate PIPN‐associated SNVs were identified from the literature. SNVs were retested in 119 patients selected by extreme phenotyping from 269 in NCCTG N08C1 (Alliance) as previously reported. SNV genotyping was performed by a combination of short‐read sequencing analysis and Taqman PCR. These 22 candidate PIPN SNVs were genotyped. Two of these, rs7349683 in the EPHA5 and rs3213619 in ABCB1 were found to be significantly associated with PIPN with an Odds ratios OR = 2.07 (P = 0.02) and OR = 0.12 (P = 0.03), respectively. In addition, three SNVs showed a trend toward a risk‐ or protective effect that was consistent with previous reports. The rs10509681 and rs11572080 in the gene CYP2C8*3 showed risk effect with an OR = 1.49 and rs1056836 in CYP1B1 showed a protective effect with an OR = 0.66. None of the other results supported the previously reported associations, including some SNVs displaying an opposite direction of effect from previous reports, including rs1058930 in CYP2C8, rs17222723 and rs8187710 in ABCC2, rs10771973 in FGD4, rs16916932 in CACNB2 and rs16948748 in PITPNA. Alliance N08C1 validated or supported a minority of previously reported SNV‐PIPN associations. Associations previously reported by multiple studies appeared to have a higher likelihood to be validated by Alliance N08C1.

Plasticity of DNA methylation in a nerve injury model of pain

The response of the peripheral nervous system (PNS) to injury may go together with alterations in epigenetics, a conjecture that has not been subjected to a comprehensive, genome-wide test. Using reduced representation bisulfite sequencing, we report widespread remodeling of DNA methylation in the rat dorsal root ganglion (DRG) occurring within 24 h of peripheral nerve ligation, a neuropathy model of allodynia. Significant (P < 10−4) cytosine hyper- and hypo-methylation was found at thousands of CpG sites. Remodeling occurred outside of CpG islands. Changes affected genes with known roles in the PNS, yet methylome remodeling also involved genes that were not linked to neuroplasticity by prior evidence. Consistent with emerging models relying on genome-wide methylation and RNA-seq analysis of promoter regions and gene bodies, variation of methylation was not tightly linked with variation of gene expression. Furthermore, approximately 44% of the dynamically changed CpGs were located outside of genes. We compared their positions with the intergenic, tissue-specific differentially methylated CpGs (tDMCs) of an independent experimental set consisting of liver, spleen, L4 control DRG, and muscle. Dynamic changes affected those intergenic CpGs that were different between tissues (P < 10−15) and almost never the invariant portion of the methylome (those CpGs that were identical across all tissues). Our findings—obtained in mixed tissue—show that peripheral nerve injury leads to methylome remodeling in the DRG. Future studies may address which of the cell types found in the DRG, such as specific groups of neurons or non-neuronal cells are affected by which aspect of the observed methylome remodeling.

Neurotoxicity to DRG neurons varies between rodent strains treated with cisplatin and bortezomib

Chemotherapy-induced peripheral neuropathy (CIPN) is a major dose limiting side effect that can lead to long-term morbidity. Approximately one-third of patients receiving chemotherapy with taxanes, vinca alkaloids, platinum compounds or proteasome inhibitors develop this toxic side effect. It is not possible to predict who will get CIPN, however, genetic susceptibility may play a role. We explored this hypothesis using an established in vitro dorsal root ganglia neurite outgrowth (DRG-NOG) assay to assess possible genetic influences for cisplatin- and bortezomib-induced neurotoxicity. Almost all previous in vitro studies have used rats or mice. We compared DRG-NOG between four genetically defined, inbred mouse strains (C57BL/6J, DBA/2J, BALB/cJ, and C3H/HeJ) and one rat strain (Sprague Dawley). Our studies found differences in cisplatin and bortezomib-induced neurotoxicity between mouse and rat strains and between the different mouse strains. C57BL/6J and Balb/cJ DRG-NOG was more sensitive to cisplatin than DBA/2J and C3H/HeJ DRG-NOG, and all mouse strains were more sensitive to cisplatin than rat. Bortezomib induced a biphasic dose response in DBA/2J and C3H/H3J mice. C57BL/6J DRG-NOG was most sensitive and Balb/cJ DRG-NOG was least sensitive to bortezomib. Our animal data supports the hypothesis that genetic background may play a role in CIPN and care must be taken when rodent models are used to better understand the contribution of genetics in patient susceptibility to CIPN.