DUX4 Transcript Knockdown with Antisense 2′-O-Methoxyethyl Gapmers for the Treatment of Facioscapulohumeral Muscular Dystrophy

In comparative high-throughput sequencing assays, a fundamental task is the analysis of count data, such as read counts per gene in RNA-seq, for evidence of systematic changes across experimental conditions. Small replicate numbers, discreteness, large dynamic range and the presence of outliers require a suitable statistical approach. We present DESeq2, a method for differential analysis of count data, using shrinkage estimation for dispersions and fold changes to improve stability and interpretability of estimates. This enables a more quantitative analysis focused on the strength rather than the mere presence of differential expression. The DESeq2 package is available at http://www.bioconductor.org/packages/release/bioc/html/DESeq2.html. Electronic supplementary material The online version of this article (doi:10.1186/s13059-014-0550-8) contains supplementary material, which is available to authorized users.

We present kallisto, an RNA-seq quantification program that is two orders of magnitude faster than previous approaches and achieves similar accuracy. Kallisto pseudoaligns reads to a reference, producing a list of transcripts that are compatible with each read while avoiding alignment of individual bases. We use kallisto to analyze 30 million unaligned paired-end RNA-seq reads in <10 min on a standard laptop computer. This removes a major computational bottleneck in RNA-seq analysis.

Facioscapulohumeral muscular dystrophy (FSHD) is a common form of muscular dystrophy in adults that is foremost characterized by progressive wasting of muscles in the upper body. FSHD is associated with contraction of D4Z4 macrosatellite repeats on chromosome 4q35, but this contraction is pathogenic only in certain "permissive" chromosomal backgrounds. Here, we show that FSHD patients carry specific single-nucleotide polymorphisms in the chromosomal region distal to the last D4Z4 repeat. This FSHD-predisposing configuration creates a canonical polyadenylation signal for transcripts derived from DUX4, a double homeobox gene of unknown function that straddles the last repeat unit and the adjacent sequence. Transfection studies revealed that DUX4 transcripts are efficiently polyadenylated and are more stable when expressed from permissive chromosomes. These findings suggest that FSHD arises through a toxic gain of function attributable to the stabilized distal DUX4 transcript.