Sequencing methods and datasets to improve functional interpretation of sleeping beauty mutagenesis screens

Background

Animal models of cancer are useful to generate complementary datasets for comparison to human tumor data. Insertional mutagenesis screens, such as those utilizing the Sleeping Beauty (SB) transposon system, provide a model that recapitulates the spontaneous development and progression of human disease. This approach has been widely used to model a variety of cancers in mice. Comprehensive mutation profiles are generated for individual tumors through amplification of transposon insertion sites followed by high-throughput sequencing. Subsequent statistical analyses identify common insertion sites (CISs), which are predicted to be functionally involved in tumorigenesis. Current methods utilized for SB insertion site analysis have some significant limitations. For one, they do not account for transposon footprints – a class of mutation generated following transposon remobilization. Existing methods also discard quantitative sequence data due to uncertainty regarding the extent to which it accurately reflects mutation abundance within a heterogeneous tumor. Additionally, computational analyses generally assume that all potential insertion sites have an equal probability of being detected under non-selective conditions, an assumption without sufficient relevant data. The goal of our study was to address these potential confounding factors in order to enhance functional interpretation of insertion site data from tumors.

Results

We describe here a novel method to detect footprints generated by transposon remobilization, which revealed minimal evidence of positive selection in tumors. We also present extensive characterization data demonstrating an ability to reproducibly assign semi-quantitative information to individual insertion sites within a tumor sample. Finally, we identify apparent biases for detection of inserted transposons in several genomic regions that may lead to the identification of false positive CISs.

Conclusion

The information we provide can be used to refine analyses of data from insertional mutagenesis screens, improving functional interpretation of results and facilitating the identification of genes important in cancer development and progression.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1150) contains supplementary material, which is available to authorized users.     

SASI-Seq: sample assurance Spike-Ins,and highly differentiating 384 barcoding for Illumina sequencing

Background

A minor but significant fraction of samples subjected to next-generation sequencing methods are either mixed-up or cross-contaminated. These events can lead to false or inconclusive results. We have therefore developed SASI-Seq; a process whereby a set of uniquely barcoded DNA fragments are added to samples destined for sequencing. From the final sequencing data, one can verify that all the reads derive from the original sample(s) and not from contaminants or other samples.

Results

By adding a mixture of three uniquely barcoded amplicons, of different sizes spanning the range of insert sizes one would normally use for Illumina sequencing, at a spike-in level of approximately 0.1%, we demonstrate that these fragments remain intimately associated with the sample. They can be detected following even the tightest size selection regimes or exome enrichment and can report the occurrence of sample mix-ups and cross-contamination.As a consequence of this work, we have designed a set of 384 eleven-base Illumina barcode sequences that are at least 5 changes apart from each other, allowing for single-error correction and very low levels of barcode misallocation due to sequencing error.

Conclusion

SASI-Seq is a simple, inexpensive and flexible tool that enables sample assurance, allows deconvolution of sample mix-ups and reports levels of cross-contamination between samples throughout NGS workflows.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-110) contains supplementary material, which is available to authorized users.     

The mutational burden of acral melanoma revealed by whole‐genome sequencing and comparative analysis

Acral melanoma is a subtype of melanoma with distinct epidemiological, clinical and mutational profiles. To define the genomic alterations in acral melanoma, we conducted whole‐genome sequencing and SNP array analysis of five metastatic tumours and their matched normal genomes. We identified the somatic mutations, copy number alterations and structural variants in these tumours and combined our data with published studies to identify recurrently mutated genes likely to be the drivers of acral melanomagenesis. We compared and contrasted the genomic landscapes of acral, mucosal, uveal and common cutaneous melanoma to reveal the distinctive mutational characteristics of each subtype.     

exonsampler: a computer program for genome‐wide and candidate gene exon sampling for targeted next‐generation sequencing

The computer program exonsampler automates the sampling of thousands of exon sequences from publicly available reference genome sequences and gene annotation databases. It was designed to provide exon sequences for the efficient, next‐generation gene sequencing method called exon capture. The exon sequences can be sampled by a list of gene name abbreviations (e.g. IFNG, TLR1), or by sampling exons from genes spaced evenly across chromosomes. It provides a list of genomic coordinates (a bed file), as well as a set of sequences in fasta format. User‐adjustable parameters for collecting exon sequences include a minimum and maximum acceptable exon length, maximum number of exonic base pairs (bp) to sample per gene, and maximum total bp for the entire collection. It allows for partial sampling of very large exons. It can preferentially sample upstream (5 prime) exons, downstream (3 prime) exons, both external exons, or all internal exons. It is written in the Python programming language using its free libraries. We describe the use of exonsampler to collect exon sequences from the domestic cow (Bos taurus) genome for the design of an exon‐capture microarray to sequence exons from related species, including the zebu cow and wild bison. We collected ~10% of the exome (~3 million bp), including 155 candidate genes, and ~16 000 exons evenly spaced genomewide. We prioritized the collection of 5 prime exons to facilitate discovery and genotyping of SNPs near upstream gene regulatory DNA sequences, which control gene expression and are often under natural selection.     

Global characterization and target identification of piRNAs and endo-siRNAs in mouse gametes and zygotes

A set of small RNAs known as rasRNAs (repeat-associated small RNAs) have been related to the down-regulation of Transposable Elements (TEs) to safeguard genome integrity. Two key members of the rasRNAs group are piRNAs and endo-siRNAs. We have performed a comparative analysis of piRNAs and endo-siRNAs present in mouse oocytes, spermatozoa and zygotes, identified by deep sequencing and bioinformatic analysis. The detection of piRNAs and endo-siRNAs in the spermatozoa and revealed also in zygotes, hints to their potential delivery to oocytes during fertilization. However, a comparative assessment of the three cell types indicates that both piRNAs and endo-siRNAs are mainly maternally inherited. Finally, we have assessed the role of the different rasRNA molecules in connection with amplification processes by way of the “ping-pong cycle”. Our results suggest that the ping-pong cycle can act on other rasRNAs, such as tRNA- and rRNA-derived fragments, thus not only being restricted to TEs during gametogenesis.     

Mechanisms of cohesin-mediated gene regulation and lessons learned from cohesinopathies

Cohesins are conserved and essential Structural Maintenance of Chromosomes (SMC) protein-containing complexes that physically interact with chromatin and modulate higher-order chromatin organization. Cohesins mediate sister chromatid cohesion and cellular long-distance chromatin interactions affecting genome maintenance and gene expression. Discoveries of mutations in cohesin's subunits and its regulator proteins in human developmental disorders, so-called “cohesinopathies,” reveal crucial roles for cohesins in development and cellular growth and differentiation. In this review, we discuss the latest findings concerning cohesin's functions in higher-order chromatin architecture organization and gene regulation and new insight gained from studies of cohesinopathies. This article is part of a Special Issue entitled: Chromatin and epigenetic regulation of animal development.     

应用亚硫酸氢盐修饰后PCR联合TA克隆测序对E-cadherin启动子区甲基化水平的检测

E-cadherin是一种细胞粘附因子,通过增强细胞之间的粘附而起到抑制肿瘤转移的作用.Ecadherin基因启动子区的高甲基化是导致其在众多肿瘤细胞中表达下调甚至缺失的主要原因之一.本实验首先抽提SGC-7901细胞(胃腺癌细胞)、A549细胞(肺腺癌细胞)、MCF-7细胞(乳腺癌细胞)等3个肿瘤细胞株的全基因组DNA,然后对抽提的DNA进行亚硫酸氢盐修饰和纯化回收,根据修饰后的DNA序列设计引物并对其进行PCR扩增.然后将PCR扩增产物与pUC-T TA载体连接并转化入感受态大肠杆菌DH5α中进行培养,对筛选出的含有阳性重组子的菌落进行测序.测序结果显示,3个肿瘤细胞株的E-cadherin基因启动子区的CpG岛都呈现了高度的甲基化,亚硫酸氢盐的修饰效率达到了99.2%.综上研究表明,亚硫酸氢盐修饰后PCR(BSP)联合TA克隆测序可以对肿瘤细胞某基因启动子区CpG岛的甲基化水平进行精确量化,研究所使用的3个肿瘤细胞株均可作为研究肿瘤细胞E-cadherin基因甲基化的细胞模型.