Alternative Poly(A) Tails Meet miRNA Targeting in Caenorhabditis elegans

Transcriptome changes during the initiation and progression of Duchenne muscular dystrophy in Caenorhabditis elegans

Duchenne muscular dystrophy (DMD) is a lethal, X-linked disease characterized by progressive muscle degeneration. The condition is driven by nonsense and missense mutations in the dystrophin gene, leading to instability of the sarcolemma and skeletal muscle necrosis and atrophy. Resulting changes in muscle-specific gene expression that take place in dystrophin's absence remain largely uncharacterized, as they are potentially obscured by the chronic inflammation elicited by muscle damage in humans. Caenorhabditis elegans possess a mild inflammatory response that is not active in the muscle, and lack a satellite cell equivalent. This allows for the characterization of the transcriptome rearrangements affecting disease progression independently of inflammation and regeneration. In effort to better understand these dynamics, we have isolated and sequenced body muscle-specific transcriptomes from C. elegans lacking functional dystrophin at distinct stages of disease progression. We have identified an upregulation of genes involved in mitochondrial function early in disease progression, and an upregulation of genes related to muscle repair in later stages. Our results suggest that in C. elegans, dystrophin may have a signaling role early in development, and its absence may activate compensatory mechanisms that counteract muscle degradation caused by loss of dystrophin. We have also developed a temperature-based screening method for synthetic paralysis that can be used to rapidly identify genetic partners of dystrophin. Our results allow for the comprehensive identification of transcriptome changes that potentially serve as independent drivers of disease progression and may in turn allow for the identification of new therapeutic targets for the treatment of DMD.

Gene expression dynamics are a proxy for selective pressures on alternatively polyadenylated isoforms

Alternative polyadenylation (APA) produces isoforms with distinct 3′-ends, yet their functional differences remain largely unknown. Here, we introduce the APA-seq method to detect the expression levels of APA isoforms from 3′-end RNA-Seq data by exploiting both paired-end reads for gene isoform identification and quantification. We detected the expression levels of APA isoforms in individual Caenorhabditis elegans embryos at different stages throughout embryogenesis. Examining the correlation between the temporal profiles of isoforms led us to distinguish two classes of genes: those with highly correlated isoforms (HCI) and those with lowly correlated isoforms (LCI) across time. We hypothesized that variants with similar expression profiles may be the product of biological noise, while the LCI variants may be under tighter selection and consequently their distinct 3′ UTR isoforms are more likely to have functional consequences. Supporting this notion, we found that LCI genes have significantly more miRNA binding sites, more correlated expression profiles with those of their targeting miRNAs and a relative lack of correspondence between their transcription and protein abundances. Collectively, our results suggest that a lack of coherence among the regulation of 3′ UTR isoforms is a proxy for selective pressures acting upon APA usage and consequently for their functional relevance.