Paired analysis of tumor mutation burden calculated by targeted deep sequencing panel and whole exome sequencing in non-small cell lung cancer

Owing to rapid advancements in NGS (next generation sequen-cing), genomic alteration is now considered an essential pre-dictive biomarkers that impact the treatment decision in many cases of cancer. Among the various predictive biomarkers, tumor mutation burden (TMB) was identified by NGS and was con-sidered to be useful in predicting a clinical response in cancer cases treated by immunotherapy. In this study, we directly com-pared the lab-developed-test (LDT) results by target sequencing panel, K-MASTER panel v3.0 and whole-exome sequencing (WES) to evaluate the concordance of TMB. As an initial step, the reference materials (n = 3) with known TMB status were used as an exploratory test. To validate and evaluate TMB, we used one hundred samples that were acquired from surgically resected tissues of non-small cell lung cancer (NSCLC) patients. The TMB of each sample was tested by using both LDT and WES methods, which extracted the DNA from samples at the same time. In addition, we evaluated the impact of capture re-gion, which might lead to different values of TMB; the evalu-ation of capture region was based on the size of NGS and target sequencing panels. In this pilot study, TMB was evalu-ated by LDT and WES by using duplicated reference samples; the results of TMB showed high concordance rate (R² = 0.887). This was also reflected in clinical samples (n = 100), which showed R² of 0.71. The difference between the coding sequence ratio (3.49%) and the ratio of mutations (4.8%) indicated that the LDT panel identified a relatively higher number of mutations. It was feasible to calculate TMB with LDT panel, which can be useful in clinical practice. Furthermore, a customized approach must be developed for calculating TMB, which differs according to cancer types and specific clinical settings.     

Comparison of solution-based exome capture methods for next generation sequencing

Background

Techniques enabling targeted re-sequencing of the protein coding sequences of the human genome on next generation sequencing instruments are of great interest. We conducted a systematic comparison of the solution-based exome capture kits provided by Agilent and Roche NimbleGen. A control DNA sample was captured with all four capture methods and prepared for Illumina GAII sequencing. Sequence data from additional samples prepared with the same protocols were also used in the comparison.

Results

We developed a bioinformatics pipeline for quality control, short read alignment, variant identification and annotation of the sequence data. In our analysis, a larger percentage of the high quality reads from the NimbleGen captures than from the Agilent captures aligned to the capture target regions. High GC content of the target sequence was associated with poor capture success in all exome enrichment methods. Comparison of mean allele balances for heterozygous variants indicated a tendency to have more reference bases than variant bases in the heterozygous variant positions within the target regions in all methods. There was virtually no difference in the genotype concordance compared to genotypes derived from SNP arrays. A minimum of 11× coverage was required to make a heterozygote genotype call with 99% accuracy when compared to common SNPs on genome-wide association arrays.

Conclusions

Libraries captured with NimbleGen kits aligned more accurately to the target regions. The updated NimbleGen kit most efficiently covered the exome with a minimum coverage of 20×, yet none of the kits captured all the Consensus Coding Sequence annotated exons.     

Molecular defects identified by whole exome sequencing in a child with Fanconi anemia

Fanconi anemia is a rare genetic disease characterized by bone marrow failure, multiple congenital malformations, and an increased susceptibility to malignancy. At least 15 genes have been identified that are involved in the pathogenesis of Fanconi anemia. However, it is still a challenge to assign the complementation group and to characterize the molecular defects in patients with Fanconi anemia. In the current study, whole exome sequencing was used to identify the affected gene(s) in a boy with Fanconi anemia. A recurring, non-synonymous mutation was found (c.3971C>T, p.P1324L) as well as a novel frameshift mutation (c.989_995del, p.H330LfsX2) in FANCA gene. Our results indicate that whole exome sequencing may be useful in clinical settings for rapid identification of disease-causing mutations in rare genetic disorders such as Fanconi anemia.     

Rare germline copy number variants in colorectal cancer predisposition characterized by exome sequencing analysis

正Colorectal cancer(CRC)is one of the most common neoplasms and an important cause of mortality worldwide(http://globocan.iarc.fr/).Approximately 35%of the variation in CRC susceptibility is likely due to heritable factors(Lichtenstein et al.,2000).Genetic variations in the human genome include single nucleotide variants(SNVs),short insertions and deletions,and larger structural vari-     

Efficient mutation identification in zebrafish by microarray capturing and next generation sequencing

Fish models like medaka, stickleback or zebrafish provide a valuable resource to study vertebrate genes. However, finding genetic variants e.g. mutations in the genome is still arduous. Here we used a combination of microarray capturing and next generation sequencing to identify the affected gene in the mozartkugel^p11cv (mzl^p11cv) mutant zebrafish. We discovered a 31-bp deletion in macf1 demonstrating the potential of this technique to efficiently isolate mutations in a vertebrate genome.     

Genetic mutation analysis at early stages of cell line development using next generation sequencing

A central goal for most biopharmaceutical companies is to reduce the development timeline to reach clinical proof of concept. This objective requires the development of tools that ensure the quality of biotherapeutic material destined for the clinic. Recent advances in high throughput protein analytics provide confidence in our ability to assess productivity and product quality attributes at early stages of cell line development. However, one quality attribute has, until recently, been absent from the standard battery of analytical tests facilitating informed choices early in cell line selection: genetic sequence confirmation. Techniques historically used for mutation analysis, such as detailed mass spectrometry, have limitations on the sample number and turnaround times making it less attractive at early stages. Thus, we explored the utility of Next‐Generation Sequencing (NGS) as a solution to address these limitations. Amplicon sequencing is one such NGS technique that is robust, rapid, sensitive, and amenable to multiplexing, all of which are essential attributes for our purposes. Here we report a NGS method based upon amplicon sequencing that has been successfully incorporated into our cell line development workflow alongside other high‐throughput protein analytical assays. The NGS method has demonstrated its value by identifying at least one Chinese hamster ovary (CHO) clone expressing a variant form of the biotherapeutic in each of the four clinical programs in which it has been utilized. We believe this sequence confirmation method is essential to safely accelerating the time to clinical proof of concept of biotherapeutics, and guard against delays related to sequence mutations. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:813–817, 2016     

Computational and statistical approaches to analyzing variants identified by exome sequencing

New sequencing technology has enabled the identification of thousands of single nucleotide polymorphisms in the exome, and many computational and statistical approaches to identify disease-association signals have emerged.     

Novel and conserved micrornas in Dalian purple urchin (Strongylocentrotus nudus) identified by next generation sequencing

MicroRNAs are regulators in regulation of broad range of phenotypes. The purple urchin, Strongylocentrotus nudus, is one of the most important marine economic animals that widely distributed in the cold seas along the coasts of eastern pacific area. To date, only 45 microRNAs have been identified in a related species, Strongylocentrotus purpurtus, and there is no report on S. nudus microRNAs. Herein, solexa sequencing technology was used to high throughput sequencing analysis of microRNAs in small RNA library isolated from five tissues of S. nudus. Totally, 8,966,865 reads were yielded, 131,015 of which were related to 415 unique microRNAs including 345 deuterostoma conserved and 70 urchin specific microRNAs, as well as 5 microRNA* sequences. The miRNA features including length distribution, end variations and genomic locations were characterized. Annotation of targets revealed a broad range of biological processes and signal transduction pathways that regulated by urchin miRNAs, of which signal transduction mechanisms was the subgroup containing the maximum targets. In addition, the expression of 100 miRNAs in female gonad was confirmed using microRNA microarray analysis. This study provides a first large scale cloning and characterization of S.nudus miRNAs and their potential targets, providing the foundation for further characterization for their role in the regulation of diversity of physiological processes.     

The statistics of bulk segregant analysis using next generation sequencing

We describe a statistical framework for QTL mapping using bulk segregant analysis (BSA) based on high throughput, short-read sequencing. Our proposed approach is based on a smoothed version of the standard G statistic, and takes into account variation in allele frequency estimates due to sampling of segregants to form bulks as well as variation introduced during the sequencing of bulks. Using simulation, we explore the impact of key experimental variables such as bulk size and sequencing coverage on the ability to detect QTLs. Counterintuitively, we find that relatively large bulks maximize the power to detect QTLs even though this implies weaker selection and less extreme allele frequency differences. Our simulation studies suggest that with large bulks and sufficient sequencing depth, the methods we propose can be used to detect even weak effect QTLs and we demonstrate the utility of this framework by application to a BSA experiment in the budding yeast Saccharomyces cerevisiae.     

全基因组测序与生物信息学分析在细菌耐药性研究中的应用

随着耐药细菌的大量出现及广泛传播,细菌耐药性成为全球备受关注的问题。耐药细菌的特征如耐药基因、毒力因子、质粒分型等以及不同菌株间亲缘关系对于细菌耐药性流行病学及分子生物学的研究有着十分重要的意义。但是传统的技术手段如聚合酶链式反应(Polymerase chain reaction,PCR)和脉冲场凝胶电泳(Pulsed field gel electrophoresis,PFGE)等得到的结果不够全面且精确度低,对于现有的研究存在很大的局限性。全基因组测序技术(Whole genome sequencing,WGS)和生物信息学分析(Bioinformatics analysis)由于能够快速详尽地得到耐药细菌的特征,也能更加精细地判断不同菌株间的进化关系,逐渐成为更加有效的技术手段,为耐药性研究提供了有效的帮助。因此,文中系统地介绍全基因组测序分析流程中的各个步骤,主要包括文库构建、平台测序以及后期数据分析三大方面的不同方法和其相应的特点,期望相关研究人员对此能够有更全面的了解,并得到一定的帮助。     

Screening the human exome: a comparison of whole genome and whole transcriptome sequencing

Background  

There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important.     

Characterization of bovine miRNAs by sequencing and bioinformatics analysis

Background  

MicroRNAs (miRNAs) are a family of ~22 nucleotide small RNA molecules which regulate gene expression by fully or partially binding to their complementary sequences in mRNAs or promoters. A large number of miRNAs and their expression patterns have been reported in human, mouse and rat. However, miRNAs and their expression patterns in live stock species such as beef cattle are not well studied.     

Novel multi-nucleotide polymorphisms in the human genome characterized by whole genome and exome sequencing

Genomic sequence comparisons between individuals are usually restricted to the analysis of single nucleotide polymorphisms (SNPs). While the interrogation of SNPs is efficient, they are not the only form of divergence between genomes. In this report, we expand the scope of polymorphism detection by investigating the occurrence of double nucleotide polymorphisms (DNPs) and triple nucleotide polymorphisms (TNPs), in which two or three consecutive nucleotides are altered compared to the reference sequence. We have found such DNPs and TNPs throughout two complete genomes and eight exomes. Within exons, these novel polymorphisms are over-represented amongst protein-altering variants; nearly all DNPs and TNPs result in a change in amino acid sequence and, in some cases, two adjacent amino acids are changed. DNPs and TNPs represent a potentially important new source of genetic variation which may underlie human disease and they should be included in future medical genetics studies. As a confirmation of the damaging nature of xNPs, we have identified changes in the exome of a glioblastoma cell line that are important in glioblastoma pathogenesis. We have found a TNP causing a single amino acid change in LAMC2 and a TNP causing a truncation of HUWE1.     

Genome-wide microRNA expression analysis of clear cell renal cell carcinoma by next generation deep sequencing

MicroRNAs (miRNAs), non-coding RNAs regulating gene expression, are frequently aberrantly expressed in human cancers. Next-generation deep sequencing technology enables genome-wide expression profiling of known miRNAs and discovery of novel miRNAs at unprecedented quantitative and qualitative accuracy. Deep sequencing was performed on 11 fresh frozen clear cell renal cell carcinoma (ccRCC) and adjacent non-tumoral renal cortex (NRC) pairs, 11 additional frozen ccRCC tissues, and 2 ccRCC cell lines (n?=?35). The 22 ccRCCs patients belonged to 3 prognostic sub-groups, i.e. those without disease recurrence, with recurrence and with metastatic disease at diagnosis. Thirty-two consecutive samples (16 ccRCC/NRC pairs) were used for stem-loop PCR validation. Novel miRNAs were predicted using 2 distinct bioinformatic pipelines. In total, 463 known miRNAs (expression frequency 1-150,000/million) were identified. We found that 100 miRNA were significantly differentially expressed between ccRCC and NRC. Differential expression of 5 miRNAs was confirmed by stem-loop PCR in the 32 ccRCC/NRC samples. With respect to RCC subgroups, 5 miRNAs discriminated between non-recurrent versus recurrent and metastatic disease, whereas 12 uniquely distinguished non-recurrent versus metastatic disease. Blocking overexpressed miR-210 or miR-27a in cell line SKCR-7 by transfecting specific antagomirs did not result in significant changes in proliferation or apoptosis. Twenty-three previously unknown miRNAs were predicted in silico. Quantitative genome-wide miRNA profiling accurately separated ccRCC from (benign) NRC. Individual differentially expressed miRNAs may potentially serve as diagnostic or prognostic markers or future therapeutic targets in ccRCC. The biological relevance of candidate novel miRNAs is unknown at present.     

A novel mutation identified in PKHD1 by targeted exome sequencing: Guiding prenatal diagnosis for an ARPKD family

Autosomal recessive polycystic kidney disease (ARPKD) is a rare hereditary renal cystic disease involving multiple organs, mainly the kidney and liver. Parents who had an affected child with ARPKD are in strong demand for an early and reliable prenatal diagnosis to guide the future pregnancies. Here we provide an example of prenatal diagnosis of an ARPKD family where traditional antenatal ultrasound examinations failed to produce conclusive results till 26th week of gestation. Compound heterozygous mutations c.274C>T (p.Arg92Trp) and c.9059T>C (p.Leu3020Pro) were identified using targeted exome sequencing in the patient and confirmed by Sanger sequencing. Further, the mother and father were revealed to be carriers of heterozygous c.274C>T and c.9059T>C mutations, respectively. Molecular prenatal diagnosis was performed for the current pregnancy by direct sequencing plus linkage analysis. Two mutations identified in the patient were both found in the fetus. In conclusion, compound heterozygous PKHD1 mutations were elucidated to be the molecular basis of the patient with ARPKD. The newly identified c.9059T>C mutation in the patient expands mutation spectrum in PKHD1 gene. For those ultrasound failed to provide clear diagnosis, we propose the new prenatal diagnosis procedure: first, screening underlying mutations in PKHD1 gene in the proband by targeted exome sequencing; then detecting causative mutations by direct sequencing in the fetal DNA and confirming results by linkage analysis.     

Statistical guidance for experimental design and data analysis of mutation detection in rare monogenic mendelian diseases by exome sequencing

Recently, whole-genome sequencing, especially exome sequencing, has successfully led to the identification of causal mutations for rare monogenic Mendelian diseases. However, it is unclear whether this approach can be generalized and effectively applied to other Mendelian diseases with high locus heterogeneity. Moreover, the current exome sequencing approach has limitations such as false positive and false negative rates of mutation detection due to sequencing errors and other artifacts, but the impact of these limitations on experimental design has not been systematically analyzed. To address these questions, we present a statistical modeling framework to calculate the power, the probability of identifying truly disease-causing genes, under various inheritance models and experimental conditions, providing guidance for both proper experimental design and data analysis. Based on our model, we found that the exome sequencing approach is well-powered for mutation detection in recessive, but not dominant, Mendelian diseases with high locus heterogeneity. A disease gene responsible for as low as 5% of the disease population can be readily identified by sequencing just 200 unrelated patients. Based on these results, for identifying rare Mendelian disease genes, we propose that a viable approach is to combine, sequence, and analyze patients with the same disease together, leveraging the statistical framework presented in this work.     

Treatment with the PARP inhibitor,niraparib, sensitizes colorectal cancer cell lines to irinotecan regardless of MSI/MSS status

Background

Cells with homologous recombination (HR) deficiency, most notably caused by mutations in the BRCA1 or BRCA2 genes, are sensitive to PARP inhibition. Microsatellite instability (MSI) accounts for 10-15% of colorectal cancer (CRC) and is hypothesized to lead to HR defects due to altered expression of Mre11, a protein required for double strand break (DSB) repair. Indeed, others have reported that PARP inhibition is efficacious in MSI CRC.

Methods

Here we examine the response to niraparib, a potent PARP-1/PARP-2 inhibitor currently under clinical evaluation, in MSI versus microsatellite stable (MSS) CRC cell lines in vitro and in vivo. We compiled a large panel of MSI and MSS CRC cell lines and evaluated the anti-proliferative activity of niraparib. In addition to testing single agent cytotoxic activity of niraparib, we also tested irinotecan (or SN-38, the active metabolite of irinotecan) activity alone and in combination with niraparib in vitro and in vivo.

Results

In contrast to earlier reports, MSI CRC cell lines were not more sensitive to niraparib than MSS CRC cell lines¸ suggesting that the MSI phenotype does not sensitize CRC cell lines to PARP inhibition. Moreover, even the most sensitive MSI cell lines had niraparib EC50s greater than 10 fold higher than BRCA-deficient cell lines. However, MSI lines were more sensitive to SN-38 than MSS lines, consistent with previous findings. We have also demonstrated that combination of niraparib and irinotecan was more efficacious than either agent alone in both MSI and MSS cell lines both in vitro and in vivo, and that niraparib potentiates the effect of irinotecan regardless of MSI status.

Conclusions

Our results support the clinical evaluation of this combination in all CRC patients, regardless of MSI status.