Resolving candidate genes of mouse skeletal muscle QTL via RNA-Seq and expression network analyses

Abstract

BackgroundWe have recently identified a number of Quantitative Trait Loci (QTL) contributing to the 2-fold muscle weight difference between the LG/J and SM/J mouse strains and refined their confidence intervals. To facilitate nomination of the candidate genes responsible for these differences we examined the transcriptome of the tibialis anterior (TA) muscle of each strain by RNA-Seq.Results13,726 genes were expressed in mouse skeletal muscle. Intersection of a set of 1061 differentially expressed transcripts with a mouse muscle Bayesian Network identified a coherent set of differentially expressed genes that we term the LG/J and SM/J Regulatory Network (LSRN). The integration of the QTL, transcriptome and the network analyses identified eight key drivers of the LSRN (Kdr, Plbd1, Mgp, Fah, Prss23, 2310014F06Rik, Grtp1, Stk10) residing within five QTL regions, which were either polymorphic or differentially expressed between the two strains and are strong candidates for quantitative trait genes (QTGs) underlying muscle mass. The insight gained from network analysis including the ability to make testable predictions is illustrated by annotating the LSRN with knowledge-based signatures and showing that the SM/J state of the network corresponds to a more oxidative state. We validated this prediction by NADH tetrazolium reductase staining in the TA muscle revealing higher oxidative potential of the SM/J compared to the LG/J strain (p<0.03).ConclusionThus, integration of fine resolution QTL mapping, RNA-Seq transcriptome information and mouse muscle Bayesian Network analysis provides a novel and unbiased strategy for nomination of muscle QTGs.

RNA-seq-based transcriptome profiling reveals differential gene expression in the lungs of Sprague–Dawley rats during early-phase acute hypobaric hypoxia

“…RNA-seq technology has been applied to different model as well as non-model organisms for a variety of studies covering allele specific gene expression, new splice variants (Saminathan et al 2015), fusion transcript variants (Norton et al 2013;Ma et al 2014), quantitative trait loci regions (Lionikas et al 2012) and tissue specific putative biomarkers (Metsalu et al 2014). We believe that there may be enhanced expression of those early-phase response genes, which may trigger signaling events that would prepare the lung tissue of tolerant animals to undergo changes that help fight or protect against or else repair the injury caused due to hypobaric hypoxia exposure, while gene expression pattern in lungs of susceptible rats might indicate an opposite pattern.…”

Transcriptomics of salinity tolerance capacity in Arctic charr (Salvelinus alpinus): a comparison of gene expression profiles between divergent QTL genotypes

“…An alternative approach is to couple RNA-Seq data with knowledge of QTL locations. By studying patterns from comparisons of gene expression data with the locations of mapped QTL, one can infer whether differential expression is cis-or trans-induced (35,49,71). For example, colocalization between a differentially expressed gene and a mapped QTL would suggest that variation in expression is due to a polymorphism in the cis-regulatory region.…”

Replication and discovery of musculoskeletal QTLs in LG/J and SM/J advanced intercross lines

“…) and the transcriptional profiles of skeletal muscle are available (Lionikas et al. ), which facilitates identification of genes underlying the QTL.…”