Ryan M. Layer

LUMPY: a probabilistic framework for structural variant discovery

Comprehensive discovery of structural variation (SV) from whole genome sequencing data requires multiple detection signals including read-pair, split-read, read-depth and prior knowledge. Owing to technical challenges, extant SV discovery algorithms either use one signal in isolation, or at best use two sequentially. We present LUMPY, a novel SV discovery framework that naturally integrates multiple SV signals jointly across multiple samples. We show that LUMPY yields improved sensitivity, especially when SV signal is reduced owing to either low coverage data or low intra-sample variant allele frequency. We also report a set of 4,564 validated breakpoints from the NA12878 human genome. https://github.com/arq5x/lumpy-sv.

Breakpoint profiling of 64 cancer genomes reveals numerous complex rearrangements spawned by homology-independent mechanisms

Tumor genomes are generally thought to evolve through a gradual accumulation of mutations, but the observation that extraordinarily complex rearrangements can arise through single mutational events suggests that evolution may be accelerated by punctuated changes in genome architecture. To assess the prevalence and origins of complex genomic rearrangements (CGRs), we mapped 6179 somatic structural variation breakpoints in 64 cancer genomes from seven tumor types and screened for clusters of three or more interconnected breakpoints. We find that complex breakpoint clusters are extremely common: 154 clusters comprise 25% of all somatic breakpoints, and 75% of tumors exhibit at least one complex cluster. Based on copy number state profiling, 63% of breakpoint clusters are consistent with being CGRs that arose through a single mutational event. CGRs have diverse architectures including focal breakpoint clusters, large-scale rearrangements joining clusters from one or more chromosomes, and staggeringly complex chromothripsis events. Notably, chromothripsis has a significantly higher incidence in glioblastoma samples (39%) relative to other tumor types (9%). Chromothripsis breakpoints also show significantly elevated intra-tumor allele frequencies relative to simple SVs, which indicates that they arise early during tumorigenesis or confer selective advantage. Finally, assembly and analysis of 4002 somatic and 6982 germline breakpoint sequences reveal that somatic breakpoints show significantly less microhomology and fewer templated insertions than germline breakpoints, and this effect is stronger at CGRs than at simple variants. These results are inconsistent with replication-based models of CGR genesis and strongly argue that nonhomologous repair of concurrently arising DNA double-strand breaks is the predominant mechanism underlying complex cancer genome rearrangements.

MICRORNA-378 targets the myogenic repressor myor during myoblast differentiation

MicroRNAs play important roles in many cell processes, including the differentiation process in several different lineages. For example, microRNAs can promote differentiation by repressing negative regulators of transcriptional activity. These regulated transcription factors can further up-regulate levels of the microRNA in a feed-forward mechanism. Here we show that MyoD up-regulates miR-378 during myogenic differentiation in C2C12 cells. ChIP and high throughput sequencing analysis shows that MyoD binds in close proximity to the miR-378 gene and causes both transactivation and chromatin remodeling. Overexpression of miR-378 increases the transcriptional activity of MyoD, in part by repressing an antagonist, MyoR. The 3′ untranslated region of MyoR contains a direct binding site for miR-378. The presence of this binding site significantly reduces the ability of MyoR to prevent the MyoD-driven transdifferentiation of fibroblasts. MyoR and miR-378 were anticorrelated during cardiotoxin-induced adult muscle regeneration in mice. Taken together, this shows a feed-forward loop where MyoD indirectly down-regulates MyoR via miR-378.

SpeedSeq: ultra-fast personal genome analysis and interpretation

SpeedSeq is an open-source genome analysis platform that accomplishes alignment, variant detection and functional annotation of a 50× human genome in 13 h on a low-cost server and alleviates a bioinformatics bottleneck that typically demands weeks of computation with extensive hands-on expert involvement. SpeedSeq offers performance competitive with or superior to current methods for detecting germline and somatic single-nucleotide variants, structural variants, insertions and deletions, and it includes novel functionality for streamlined interpretation.

Extrachromosomal MicroDNAs and Chromosomal Microdeletions in Normal Tissues

MicroDNA, Microdeletion? Extrachromosomal circular DNA (eccDNA) is widespread in eukaryotes. Shibata et al. (p. 82, published online 8 March) have found a new form of eccDNA—microDNA—in embryonic mouse brain, heart, and liver in adult mouse brain and in mouse and human cell lines. Unlike eccDNA, which is derived from repetitive sequences, transposable elements, or viral genomes, circular microDNA is derived from nonrepetitive genomic sequences and, very often, regions associated with genes. Length profiles—at 200 to 400 base pairs long—suggest that microDNA preferentially arises from regions of DNA occupied by nucleosomes. Furthermore, microDNA has very short regions of micro-homology at the beginning and ends of the circles. Analysis of several very common microDNAs indicates that their formation seems to result in microdeletions in somatic and germline genomes. The formation of circular microDNAs can result in somatic and germline microdeletions in mammalian cells. We have identified tens of thousands of short extrachromosomal circular DNAs (microDNA) in mouse tissues as well as mouse and human cell lines. These microDNAs are 200 to 400 base pairs long, are derived from unique nonrepetitive sequence, and are enriched in the 5′-untranslated regions of genes, exons, and CpG islands. Chromosomal loci that are enriched sources of microDNA in the adult brain are somatically mosaic for microdeletions that appear to arise from the excision of microDNAs. Germline microdeletions identified by the “Thousand Genomes” project may also arise from the excision of microDNAs in the germline lineage. We have thus identified a previously unknown DNA entity in mammalian cells and provide evidence that their generation leaves behind deletions in different genomic loci.

MUNC, a Long Noncoding RNA That Facilitates the Function of MyoD in Skeletal Myogenesis

ABSTRACT An in silico screen for myogenic long noncoding RNAs (lncRNAs) revealed nine lncRNAs that are upregulated more than 10-fold in myotubes versus levels in myoblasts. One of these lncRNAs, MyoD upstream noncoding (MUNC, also known as DRReRNA), is encoded 5 kb upstream of the transcription start site of MyoD, a myogenic transcription factor gene. MUNC is specifically expressed in skeletal muscle and exists as in unspliced and spliced isoforms, and its 5′ end overlaps with the cis-acting distal regulatory region (DRR) of MyoD. Small interfering RNA (siRNA) of MUNC reduced myoblast differentiation and specifically reduced the association of MyoD to the DRR enhancer and myogenin promoter but not to another MyoD-dependent enhancer. Stable overexpression of MUNC from a heterologous promoter increased endogenous MyoD, Myogenin, and Myh3 (myosin heavy chain, [MHC] gene) mRNAs but not the cognate proteins, suggesting that MUNC can act in trans to promote gene expression but that this activity does not require an induction of MyoD protein. MUNC also stimulates the transcription of other genes that are not recognized as MyoD-inducible genes. Knockdown of MUNC in vivo impaired murine muscle regeneration, implicating MUNC in primary satellite cell differentiation in the animal. We also discovered a human MUNC that is induced during differentiation of myoblasts and whose knockdown decreases differentiation, suggesting an evolutionarily conserved role of MUNC lncRNA in myogenesis. Although MUNC overlaps with the DRR enhancer, our results suggest that MUNC is not a classic cis-acting enhancer RNA (e-RNA) acting exclusively by stimulating the neighboring MyoD gene but more like a promyogenic lncRNA that acts directly or indirectly on multiple promoters to increase myogenic gene expression.

Nifty assignments

Creating assignments is a difficult and time consuming part of teaching Computer Science. Nifty Assignments is a forum, operating at a very practical level, to promote the sharing of assignment ideas and assignment materials.Each presenter will introduce their assignment, give a quick demo, and describe its niche in the curriculum and its strengths and weaknesses. The presentations (and the descriptions below) merely introduce each assignment. For more detail, each assignment has its own web page with more detailed information and assignment materials such as handouts and data files to aid the adoption of the assignment. Information on participating in Nifty Assignments as well as all the assignment pages are available from our central page… http://cse.stanford.edu/nifty/

Regulation of several androgen-induced genes through the repression of the miR-99a/let-7c/miR-125b-2 miRNA cluster in prostate cancer cells

The androgen receptor (AR) stimulates and represses gene expression to promote the initiation and progression of prostate cancer. Here, we report that androgen represses the miR-99a/let7c/125b-2 cluster through AR and anti-androgen drugs block the androgen-repression of the miRNA cluster. AR directly binds to the host gene of the miR-99a/let7c/125b-2 cluster, LINC00478. Expression of the cluster is repressed or activated by chromatin remodelers EZH2 or JMJD3 in the presence or absence of androgen, respectively. Bioinformatics analysis reveals a significant enrichment of targets of miR-99a, let-7c and miR-125b in androgen-induced gene sets, suggesting that downregulation of the miR-99a/let7c/125b-2 cluster by androgen protects many of their target mRNAs from degradation and indirectly assists in the gene induction. We validated the hypothesis with 12 potential targets of the miR-99a/let7c/125b-2 cluster induced by androgen: 9 out of the 12 mRNAs are downregulated by the microRNA cluster. To ascertain the biological significance of this hypothesis, we focused on IGF1R, a known prostate cancer growth factor that is induced by androgen and directly targeted by the miR-99a/let7c/125b-2 cluster. The androgen-induced cell proliferation is ameliorated to a similar extent as anti-androgen drugs by preventing the repression of the microRNAs or induction of IGF1R in androgen-dependent prostate cancer cells. Expression of a microRNA-resistant form of IGF1R protects these cells from inhibition by the miR-99a/let7c/125b-2 cluster. These results indicate that a thorough understanding of how androgen stimulates prostate cancer growth requires not only an understanding of genes directly induced/repressed by AR, but also of genes indirectly induced by AR through the repression of key microRNAs.

Notch3 and Mef2c Proteins Are Mutually Antagonistic via Mkp1 Protein and miR-1/206 MicroRNAs in Differentiating Myoblasts

Background: Notch3 is expressed in myogenic precursors, but its function is not well known. Results: Notch3 represses the activity of Mef2c and is in turn inhibited by the microRNAs-1 and -206. Conclusion: Notch3 serves as a regulator for preventing premature myogenic differentiation. Significance: Understanding how precocious differentiation is prevented is critical for designing therapy for skeletal muscle regeneration. The Notch signaling pathway is a well known regulator of skeletal muscle stem cells known as satellite cells. Loss of Notch1 signaling leads to spontaneous myogenic differentiation. Notch1, normally expressed in satellite cells, is targeted for proteasomal degradation by Numb during differentiation. A homolog of Notch1, Notch3, is also expressed in these cells but is not inhibited by Numb. We find that Notch3 is paradoxically up-regulated during the early stages of differentiation by an enhancer that requires both MyoD and activated Notch1. Notch3 itself strongly inhibits the myogenic transcription factor Mef2c, most likely by increasing the p38 phosphatase Mkp1, which inhibits the Mef2c activator p38 MAP kinase. Active Notch3 decreases differentiation. Mef2c, however, induces microRNAs miR-1 and miR-206, which directly down-regulate Notch3 and allow differentiation to proceed. Thus, the myogenic differentiation-induced microRNAs miR-1 and -206 are important for differentiation at least partly because they turn off Notch3. We suggest that the transient expression of Notch3 early in differentiation generates a temporal lag between myoblast activation by MyoD and terminal differentiation into myotubes directed by Mef2c.

Binary Interval Search

MOTIVATION The comparison of diverse genomic datasets is fundamental to understand genome biology. Researchers must explore many large datasets of genome intervals (e.g. genes, sequence alignments) to place their experimental results in a broader context and to make new discoveries. Relationships between genomic datasets are typically measured by identifying intervals that intersect, that is, they overlap and thus share a common genome interval. Given the continued advances in DNA sequencing technologies, efficient methods for measuring statistically significant relationships between many sets of genomic features are crucial for future discovery. RESULTS We introduce the Binary Interval Search (BITS) algorithm, a novel and scalable approach to interval set intersection. We demonstrate that BITS outperforms existing methods at counting interval intersections. Moreover, we show that BITS is intrinsically suited to parallel computing architectures, such as graphics processing units by illustrating its utility for efficient Monte Carlo simulations measuring the significance of relationships between sets of genomic intervals. AVAILABILITY https://github.com/arq5x/bits.