Epigenetic loss of heterozygosity of Apc and an inflammation-associated mutational signature detected in Lrig1+/--driven murine colonic adenomas.

BACKGROUND: The loss of a single copy of adenomatous polyposis coli (Apc) in leucine-rich repeats and immunoglobulin-like domains 1 (Lrig1)-expressing colonic progenitor cells induces rapid growth of adenomas in mice with high penetrance and multiplicity. The tumors lack functional APC, and a genetic loss of heterozygosity of Apc was previously observed. METHODS: To identify genomic features of early tumorigenesis, and to profile intertumoral genetic heterogeneity, tumor exome DNA (n = 9 tumors) and mRNA (n = 5 tumors) sequences were compared with matched nontumoral colon tissue. Putative somatic mutations were called after stringent variant filtering. Somatic signatures of mutational processes were determined and splicing patterns were observed. RESULTS: The adenomas were found to be genetically heterogeneous and unexpectedly hypermutated, displaying a strong bias toward G:C > A:T mutations. A genetic loss of heterozygosity of Apc was not observed, however, an epigenetic loss of heterozygosity was apparent in the tumor transcriptomes. Complex splicing patterns characterized by a loss of intron retention were observed uniformly across tumors. CONCLUSION: This study demonstrates that early tumors originating from intestinal stem cells with reduced Lrig1 and Apc expression are highly mutated and genetically heterogeneous, with an inflammation-associated mutational signature and complex splicing patterns that are uniform across tumors.

High-specificity detection of rare alleles with Paired-End Low Error Sequencing (PELE-Seq).

BACKGROUND: Polymorphic loci exist throughout the genomes of a population and provide the raw genetic material needed for a species to adapt to changes in the environment. The minor allele frequencies of rare Single Nucleotide Polymorphisms (SNPs) within a population have been difficult to track with Next-Generation Sequencing (NGS), due to the high error rate of standard methods such as Illumina sequencing. RESULTS: We have developed a wet-lab protocol and variant-calling method that identifies both sequencing and PCR errors, called Paired-End Low Error Sequencing (PELE-Seq). To test the specificity and sensitivity of the PELE-Seq method, we sequenced control E. coli DNA libraries containing known rare alleles present at frequencies ranging from 0.2-0.4 % of the total reads. PELE-Seq had higher specificity and sensitivity than standard libraries. We then used PELE-Seq to characterize rare alleles in a Caenorhabditis remanei nematode worm population before and after laboratory adaptation, and found that minor and rare alleles can undergo large changes in frequency during lab-adaptation. CONCLUSION: We have developed a method of rare allele detection that mitigates both sequencing and PCR errors, called PELE-Seq. PELE-Seq was evaluated using control E. coli populations and was then used to compare a wild C. remanei population to a lab-adapted population. The PELE-Seq method is ideal for investigating the dynamics of rare alleles in a broad range of reduced-representation sequencing methods, including targeted amplicon sequencing, RAD-Seq, ddRAD, and GBS. PELE-Seq is also well-suited for whole genome sequencing of mitochondria and viruses, and for high-throughput rare mutation screens.

Genome-wide identification of hypoxia-induced enhancer regions.

Here we present a genome-wide method for de novo identification of enhancer regions. This approach enables massively parallel empirical investigation of DNA sequences that mediate transcriptional activation and provides a platform for discovery of regulatory modules capable of driving context-specific gene expression. The method links fragmented genomic DNA to the transcription of randomer molecule identifiers and measures the functional enhancer activity of the library by massively parallel sequencing. We transfected a Drosophila melanogaster library into S2 cells in normoxia and hypoxia, and assayed 4,599,881 genomic DNA fragments in parallel. The locations of the enhancer regions strongly correlate with genes up-regulated after hypoxia and previously described enhancers. Novel enhancer regions were identified and integrated with RNAseq data and transcription factor motifs to describe the hypoxic response on a genome-wide basis as a complex regulatory network involving multiple stress-response pathways. This work provides a novel method for high-throughput assay of enhancer activity and the genome-scale identification of 31 hypoxia-activated enhancers in Drosophila.

Local de novo assembly of RAD paired-end contigs using short sequencing reads.

Despite the power of massively parallel sequencing platforms, a drawback is the short length of the sequence reads produced. We demonstrate that short reads can be locally assembled into longer contigs using paired-end sequencing of restriction-site associated DNA (RAD-PE) fragments. We use this RAD-PE contig approach to identify single nucleotide polymorphisms (SNPs) and determine haplotype structure in threespine stickleback and to sequence E. coli and stickleback genomic DNA with overlapping contigs of several hundred nucleotides. We also demonstrate that adding a circularization step allows the local assembly of contigs up to 5 kilobases (kb) in length. The ease of assembly and accuracy of the individual contigs produced from each RAD site sequence suggests RAD-PE sequencing is a useful way to convert genome-wide short reads into individually-assembled sequences hundreds or thousands of nucleotides long.