Targeted next generation sequencing identifies clinically actionable mutations in patients with melanoma

Somatic sequencing of cancers has produced new insight into tumorigenesis, tumor heterogeneity, and disease progression, but the vast majority of genetic events identified are of indeterminate clinical significance. Here, we describe a NextGen sequencing approach to fully analyzing 248 genes, including all those of known clinical significance in melanoma. This strategy features solution capture of DNA followed by multiplexed, high‐throughput sequencing and was evaluated in 31 melanoma cell lines and 18 tumor tissues from patients with metastatic melanoma. Mutations in melanoma cell lines correlated with their sensitivity to corresponding small molecule inhibitors, confirming, for example, lapatinib sensitivity in ERBB4 mutant lines and identifying a novel activating mutation of BRAF. The latter event would not have been identified by clinical sequencing and was associated with responsiveness to a BRAF kinase inhibitor. This approach identified focal copy number changes of PTEN not found by standard methods, such as comparative genomic hybridization (CGH). Actionable mutations were found in 89% of the tumor tissues analyzed, 56% of which would not be identified by standard‐of‐care approaches. This work shows that targeted sequencing is an attractive approach for clinical use in melanoma.     

Multiple mitochondrial introgression events and heteroplasmy in trypanosoma cruzi revealed by maxicircle MLST and next generation sequencing

Background

Mitochondrial DNA is a valuable taxonomic marker due to its relatively fast rate of evolution. In Trypanosoma cruzi, the causative agent of Chagas disease, the mitochondrial genome has a unique structural organization consisting of 20–50 maxicircles (∼20 kb) and thousands of minicircles (0.5–10 kb). T. cruzi is an early diverging protist displaying remarkable genetic heterogeneity and is recognized as a complex of six discrete typing units (DTUs). The majority of infected humans are asymptomatic for life while 30–35% develop potentially fatal cardiac and/or digestive syndromes. However, the relationship between specific clinical outcomes and T. cruzi genotype remains elusive. The availability of whole genome sequences has driven advances in high resolution genotyping techniques and re-invigorated interest in exploring the diversity present within the various DTUs.

Methodology/Principal Findings

To describe intra-DTU diversity, we developed a highly resolutive maxicircle multilocus sequence typing (mtMLST) scheme based on ten gene fragments. A panel of 32 TcI isolates was genotyped using the mtMLST scheme, GPI, mini-exon and 25 microsatellite loci. Comparison of nuclear and mitochondrial data revealed clearly incongruent phylogenetic histories among different geographical populations as well as major DTUs. In parallel, we exploited read depth data, generated by Illumina sequencing of the maxicircle genome from the TcI reference strain Sylvio X10/1, to provide the first evidence of mitochondrial heteroplasmy (heterogeneous mitochondrial genomes in an individual cell) in T. cruzi.

Conclusions/Significance

mtMLST provides a powerful approach to genotyping at the sub-DTU level. This strategy will facilitate attempts to resolve phenotypic variation in T. cruzi and to address epidemiologically important hypotheses in conjunction with intensive spatio-temporal sampling. The observations of both general and specific incidences of nuclear-mitochondrial phylogenetic incongruence indicate that genetic recombination is geographically widespread and continues to influence the natural population structure of TcI, a conclusion which challenges the traditional paradigm of clonality in T. cruzi.     

Intratumor heterogeneity and clonal evolution revealed in castration-resistant prostate cancer by longitudinal genomic analysis

Intratumor heterogeneity is a key driver for local relapse and treatment failure. Thus, using multifocal prostate cancer as a model to investigate tumor inter-clonal relationships and tumor evolution could aid in our understanding of drug resistance. Previous studies discovered genomic alterations by comparing hormone-sensitive prostate cancer (HSPC) with castration-resistant prostate cancer (CRPC) in large cohorts. However, most studies did not sequentially sample tumors from the same patient. In our study, we performed whole-exome sequencing (WES) on 14 specimens from five locally relapsed patients before and after androgen-deprivation therapy. We described the landscape of genomic alterations before and after treatment and identified critical driver events that could have contributed to the evolution of CRPC. In addition to confirming known cancer genes such as TP53 and CDK12, we also identified new candidate genes that may play a role in the progression of prostate cancer, including MYO15A, CHD6 and LZTR1. At copy number alteration (CNA) level, gain of 8q24.13-8q24.3 was observed in 60% of patients and was the most commonly altered locus in both HSPC and CRPC tumors. Finally, utilizing phylogenetic reconstruction, we explored the clonal progression pattern from HSPC to CRPC in each patient. Our findings highlight the complex and heterogeneous mechanisms underlying the development of drug resistance, and underscore the potential value of monitoring tumor clonal architectures during disease progression in a clinical setting.     

Dynamic evolution of clonal epialleles revealed by methclone

We describe methclone, a novel method to identify epigenetic loci that harbor large changes in the clonality of their epialleles (epigenetic alleles). Methclone efficiently analyzes genome-wide DNA methylation sequencing data. We quantify the changes using a composition entropy difference calculation and also introduce a new measure of global clonality shift, loci with epiallele shift per million loci covered, which enables comparisons between different samples to gauge overall epiallelic dynamics. Finally, we demonstrate the utility of methclone in capturing functional epiallele shifts in leukemia patients from diagnosis to relapse. Methclone is open-source and freely available at https://code.google.com/p/methclone.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0472-5) contains supplementary material, which is available to authorized users.     

Targeted next generation sequencing identifies novel NOTCH3 gene mutations in CADASIL diagnostics patients

Background

Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is a monogenic, hereditary, small vessel disease of the brain causing stroke and vascular dementia in adults. CADASIL has previously been shown to be caused by varying mutations in the NOTCH3 gene. The disorder is often misdiagnosed due to its significant clinical heterogeneic manifestation with familial hemiplegic migraine and several ataxia disorders as well as the location of the currently identified causative mutations. The aim of this study was to develop a new, comprehensive and efficient single assay strategy for complete molecular diagnosis of NOTCH3 mutations through the use of a custom next-generation sequencing (NGS) panel for improved routine clinical molecular diagnostic testing.

Results

Our custom NGS panel identified nine genetic variants in NOTCH3 (p.D139V, p.C183R, p.R332C, p.Y465C, p.C597W, p.R607H, p.E813E, p.C977G and p.Y1106C). Six mutations were stereotypical CADASIL mutations leading to an odd number of cysteine residues in one of the 34 NOTCH3 gene epidermal growth factor (EGF)-like repeats, including three new typical cysteine mutations identified in exon 11 (p.C597W; c.1791C>G); exon 18 (p.C977G; c.2929T>G) and exon 20 (p.Y1106C; c.3317A>G). Interestingly, a novel missense mutation in the CACNA1A gene was also identified in one CADASIL patient. All variants identified (novel and known) were further investigated using in silico bioinformatic analyses and confirmed through Sanger sequencing.

Conclusions

NGS provides an improved and effective methodology for the diagnosis of CADASIL. The NGS approach reduced time and cost for comprehensive genetic diagnosis, placing genetic diagnostic testing within reach of more patients.

     

新一代测序技术应用于microRNA检测

MicroRNAs(miRNAs)是一类在进化上高度保守的非编码小分子单链RNA(~22nt),在基因转录后调控中发挥至关重要的作用。越来越多的证据表明,miRNAs参与很多重要的生理和病理过程,例如发育、器官形成、调亡、细胞增殖、肿瘤发生等。近年来飞速发展的新一代测序技术在miRNA检测方面具有重要的应用。文章简要介绍了新一代测序技术3大平台的基本步骤和原理,测序数据的生物信息学分析方法以及新一代测序技术在miRNA方向的主要应用。相比于传统的miRNA检测方法,新一代测序技术具有通量高、对遗传物质检测完全且准确度高,可重复性好等优点,在探索新miRNA、miRNA互补链、miRNA编辑、miRNA异构体检测以及miRNA靶基因检测等方面具有巨大优势。随着新一代测序技术的不断发展,测序成本不断降低,在未来几年,新一代测序技术的使用率或将大大增加。新一代测序技术的不断应用将进一步促进人类对于miRNA在各种生理病理过程中的功能和调控的认识。     

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.     

Complete plastid genome sequence of Vaccinium macrocarpon: structure, gene content, and rearrangements revealed by next generation sequencing

The complete plastid genome sequence of the American cranberry (Vaccinium macrocarpon Ait.) was reconstructed using next-generation sequencing data by in silico procedures. We used Roche 454 shotgun sequence data to isolate cranberry plastid-specific sequences of “HyRed” via homology comparisons with complete sequences from several species available at the National Center for Biotechnology Information database. Eleven cranberry plastid contigs were selected for the construction of the plastid genome-based homologies and on raw reads flowing through contigs and connection information. We assembled and annotated a cranberry plastid genome (82,284 reads; 185x coverage) with a length of 176 kb and the typical structure found in plants, but with several structural rearrangements in the large single-copy region when compared to other plastid asterid genomes. To evaluate the reliability of the sequence data, phylogenetic analysis of 30 species outside the order Ericales (with 54 genes) showed Vaccinium inside the clade Asteridae, as reported in other studies using single genes. The cranberry plastid genome sequence will allow the accumulation of critical data useful for breeding and a suite of other genetic studies.     

Comparison of the experimental methods in haplotype sequencing via next generation sequencing

Although the diploid nature has been observed for over 50 years, phasing the diploid is still a laborious task. The speed and throughput of next generation sequencing have largely increased in the past decades. However, the short read-length remains one of the biggest challenges of haplotype analysis. For instance, reads as short as 150 bp span no more than one variant in most cases. Numerous experimental technologies have been developed to overcome this challenge. Distance, complexity and accuracy of the linkages obtained are the main factors to evaluate the efficiency of whole genome haplotyping methods. Here, we review these experimental technologies, evaluating their efficiency in linkages obtaining and system complexity. The technologies are organized into four categories based on its strategy: (i) chromosomes separation, (ii) dilution pools, (iii) crosslinking and proximity ligation, (ix) long-read technologies. Within each category, several subsections are listed to classify each technology. Innovative experimental strategies are expected to have high-quality performance, low cost and be labor-saving, which will be largely desired in the future.     

新一代测序技术在植物转录组研究中的应用

随着DNA测序技术的发展,新一代测序技术以其高通量、低成本的特点,成为越来越多的生物学研究者在开展工作时的首选。在所有的新一代测序技术中,454测序系统是最早实现商业化且发展相对成熟的一种,目前被广泛的应用于各个领域的生物学研究中。文章以454测序系统为例,综述了新一代测序系统的原理、优缺点,及其在植物转录组研究中的应用,并对其在植物研究领域中可能的发展应用方向进行了展望。     

Personalized pathway enrichment map of putative cancer genes from next generation sequencing data

BACKGROUND: Pathway analysis of a set of genes represents an important area in large-scale omic data analysis. However, the application of traditional pathway enrichment methods to next-generation sequencing (NGS) data is prone to several potential biases, including genomic/genetic factors (e.g., the particular disease and gene length) and environmental factors (e.g., personal life-style and frequency and dosage of exposure to mutagens). Therefore, novel methods are urgently needed for these new data types, especially for individual-specific genome data. METHODOLOGY: In this study, we proposed a novel method for the pathway analysis of NGS mutation data by explicitly taking into account the gene-wise mutation rate. We estimated the gene-wise mutation rate based on the individual-specific background mutation rate along with the gene length. Taking the mutation rate as a weight for each gene, our weighted resampling strategy builds the null distribution for each pathway while matching the gene length patterns. The empirical P value obtained then provides an adjusted statistical evaluation. PRINCIPAL FINDINGS/CONCLUSIONS: We demonstrated our weighted resampling method to a lung adenocarcinomas dataset and a glioblastoma dataset, and compared it to other widely applied methods. By explicitly adjusting gene-length, the weighted resampling method performs as well as the standard methods for significant pathways with strong evidence. Importantly, our method could effectively reject many marginally significant pathways detected by standard methods, including several long-gene-based, cancer-unrelated pathways. We further demonstrated that by reducing such biases, pathway crosstalk for each individual and pathway co-mutation map across multiple individuals can be objectively explored and evaluated. This method performs pathway analysis in a sample-centered fashion, and provides an alternative way for accurate analysis of cancer-personalized genomes. It can be extended to other types of genomic data (genotyping and methylation) that have similar bias problems.     

ProteinSeq: high-performance proteomic analyses by proximity ligation and next generation sequencing

Despite intense interest, methods that provide enhanced sensitivity and specificity in parallel measurements of candidate protein biomarkers in numerous samples have been lacking. We present herein a multiplex proximity ligation assay with readout via realtime PCR or DNA sequencing (ProteinSeq). We demonstrate improved sensitivity over conventional sandwich assays for simultaneous analysis of sets of 35 proteins in 5 μl of blood plasma. Importantly, we observe a minimal tendency to increased background with multiplexing, compared to a sandwich assay, suggesting that higher levels of multiplexing are possible. We used ProteinSeq to analyze proteins in plasma samples from cardiovascular disease (CVD) patient cohorts and matched controls. Three proteins, namely P-selectin, Cystatin-B and Kallikrein-6, were identified as putative diagnostic biomarkers for CVD. The latter two have not been previously reported in the literature and their potential roles must be validated in larger patient cohorts. We conclude that ProteinSeq is promising for screening large numbers of proteins and samples while the technology can provide a much-needed platform for validation of diagnostic markers in biobank samples and in clinical use.     

The reproductive biology of Polytrichum formosum: clonal structure and paternity revealed by microsatellites

Using highly polymorphic microsatellite markers, we assessed clonal structure and paternity in a population of the bryophyte species Polytrichum formosum. Identical multilocus genotypes of individual shoots were almost never observed in spatially separated cushions, but were found to be highly clustered within moss cushions. Therefore, asexual reproduction through dispersal of gametophyte fragments is not very important in P. formosum. However, asexual reproduction on a very localized scale through vegetative growth of genets (branching of gametophytes via clonal growth of rhizomes) is very extensive. The patchy spatial distribution of genets and the absence of intermingling among genets suggest that this species follows a 'phalanx' clonal growth strategy. Vegetative proliferation of genets will increase their size, and, consequently, will have considerable fitness consequences for individuals in terms of increased genet longevity and reproductive output. Although paternity analysis of sporophytes confirmed male genet size, i.e. gamete production, to be an important determinant of male reproductive fitness, it also showed that the spatial distance to female genets is the predominant factor that governs male reproductive success. Moreover, we showed that male gamete dispersal distances in P. formosum are much further than generally assumed, and are in the order of metres rather than centimetres. Combining the findings, we conclude that the high genotypic diversity observed for this facultatively clonal species is most likely explained by a preponderance of sexual reproduction over clonal reproduction.     

Characterization of a biogas-producing microbial community by short-read next generation DNA sequencing

ABSTRACT: BACKGROUND: Renewable energy production is currently a major issue worldwide. Biogas is a promising renewable energy carrier as the technology of its production combines the elimination of organic waste with the formation of a versatile energy carrier, methane. In consequence of the complexity of the microbial communities and metabolic pathways involved the biotechnology of the microbiological process leading to biogas production is poorly understood. Metagenomic approaches are suitable means of addressing related questions. In the present work a novel high-throughput technique was tested for its benefits in resolving the functional and taxonomical complexity of such microbial consortia. RESULTS: It was demonstrated that the extremely parallel SOLiDTM short-read DNA sequencing platform is capable of providing sufficient useful information to decipher the systematic and functional contexts within a biogas-producing community. Although this technology has not been employed to address such problems previously, the data obtained compare well with those from similar high-throughput approaches such as 454-pyrosequencing GS FLX or Titanium. The predominant microbes contributing to the decomposition of organic matter include members of the Eubacteria, class Clostridia, order Clostridiales, family Clostridiaceae. Bacteria belonging in other systematic groups contribute to the diversity of the microbial consortium. Archaea comprise a remarkably small minority in this community, given their crucial role in biogas production. Among the Archaea, the predominant order is the Methanomicrobiales and the most abundant species is Methanoculleus marisnigri. The Methanomicrobiales are hydrogenotrophic methanogens. Besides corroborating earlier findings on the significance of the contribution of the Clostridia to organic substrate decomposition, the results demonstrate the importance of the metabolism of hydrogen within the biogas producing microbial community. CONCLUSIONS: Both microbiological diversity and the regulatory role of the hydrogen metabolism appear to be the driving forces optimizing biogas-producing microbial communities. The findings may allow a rational design of these communities to promote greater efficacy in large-scale practical systems. The composition of an optimal biogas-producing consortium can be determined through the use of this approach, and this systematic methodology allows the design of the optimal microbial community structure for any biogas plant. In this way, metagenomic studies can contribute to significant progress in the efficacy and economic improvement of biogas production.     

Somatic mutation profiles of MSI and MSS colorectal cancer identified by whole exome next generation sequencing and bioinformatics analysis

Background

Colorectal cancer (CRC) is with approximately 1 million cases the third most common cancer worldwide. Extensive research is ongoing to decipher the underlying genetic patterns with the hope to improve early cancer diagnosis and treatment. In this direction, the recent progress in next generation sequencing technologies has revolutionized the field of cancer genomics. However, one caveat of these studies remains the large amount of genetic variations identified and their interpretation.

Methodology/Principal Findings

Here we present the first work on whole exome NGS of primary colon cancers. We performed 454 whole exome pyrosequencing of tumor as well as adjacent not affected normal colonic tissue from microsatellite stable (MSS) and microsatellite instable (MSI) colon cancer patients and identified more than 50,000 small nucleotide variations for each tissue. According to predictions based on MSS and MSI pathomechanisms we identified eight times more somatic non-synonymous variations in MSI cancers than in MSS and we were able to reproduce the result in four additional CRCs. Our bioinformatics filtering approach narrowed down the rate of most significant mutations to 359 for MSI and 45 for MSS CRCs with predicted altered protein functions. In both CRCs, MSI and MSS, we found somatic mutations in the intracellular kinase domain of bone morphogenetic protein receptor 1A, BMPR1A, a gene where so far germline mutations are associated with juvenile polyposis syndrome, and show that the mutations functionally impair the protein function.

Conclusions/Significance

We conclude that with deep sequencing of tumor exomes one may be able to predict the microsatellite status of CRC and in addition identify potentially clinically relevant mutations.