Tbx15 controls skeletal muscle fibre-type determination and muscle metabolism

Multimodal cell atlas of the ageing human skeletal muscle

Muscle atrophy and functional decline (sarcopenia) are common manifestations of frailty and are critical contributors to morbidity and mortality in older people1. Deciphering the molecular mechanisms underlying sarcopenia has major implications for understanding human ageing2. Yet, progress has been slow, partly due to the difficulties of characterizing skeletal muscle niche heterogeneity (whereby myofibres are the most abundant) and obtaining well-characterized human samples3,4. Here we generate a single-cell/single-nucleus transcriptomic and chromatin accessibility map of human limb skeletal muscles encompassing over 387,000 cells/nuclei from individuals aged 15 to 99 years with distinct fitness and frailty levels. We describe how cell populations change during ageing, including the emergence of new populations in older people, and the cell-specific and multicellular network features (at the transcriptomic and epigenetic levels) associated with these changes. On the basis of cross-comparison with genetic data, we also identify key elements of chromatin architecture that mark susceptibility to sarcopenia. Our study provides a basis for identifying targets in the skeletal muscle that are amenable to medical, pharmacological and lifestyle interventions in late life.

Vitamin A regulates mitochondrial biogenesis and function through p38 MAPK-PGC-1α signaling pathway and alters the muscle fiber composition of sheep

BACKGROUND

Vitamin A (VA) and its metabolite, retinoic acid (RA), are of great interest for their wide range of physiological functions. However, the regulatory contribution of VA to mitochondrial and muscle fiber composition in sheep has not been reported.

METHOD

Lambs were injected with 0 (control) or 7,500 IU VA palmitate into the biceps femoris muscle on d 2 after birth. At the age of 3 and 32 weeks, longissimus dorsi (LD) muscle samples were obtained to explore the effect of VA on myofiber type composition. In vitro, we investigated the effects of RA on myofiber type composition and intrinsic mechanisms.

RESULTS

The proportion of type I myofiber was greatly increased in VA-treated sheep in LD muscle at harvest. VA greatly promoted mitochondrial biogenesis and function in LD muscle of sheep. Further exploration revealed that VA elevated PGC-1α mRNA and protein contents, and enhanced the level of p38 MAPK phosphorylation in LD muscle of sheep. In addition, the number of type I myofibers with RA treatment was significantly increased, and type IIx myofibers was significantly decreased in primary myoblasts. Consistent with in vivo experiment, RA significantly improved mitochondrial biogenesis and function in primary myoblasts of sheep. We then used si-PGC-1α to inhibit PGC-1α expression and found that si-PGC-1α significantly abrogated RA-induced the formation of type I myofibers, mitochondrial biogenesis, MitoTracker staining intensity, UQCRC1 and ATP5A1 expression, SDH activity, and enhanced the level of type IIx muscle fibers. These data suggested that RA improved mitochondrial biogenesis and function by promoting PGC-1α expression, and increased type I myofibers. In order to prove that the effect of RA on the level of PGC-1α is caused by p38 MAPK signaling, we inhibited the p38 MAPK signaling using a p38 MAPK inhibitor, which significantly reduced RA-induced PGC-1α and MyHC I levels.

CONCLUSION

VA promoted PGC-1α expression through the p38 MAPK signaling pathway, improved mitochondrial biogenesis, and altered the composition of muscle fiber type.

Effects of pregnancy and lactation on muscle-tendon morphology

Pregnancy and lactation hormones have been shown to mediate anatomical changes to the musculoskeletal system that generates animal movement. In this study, we characterize changes in the medial gastrocnemius muscle, its tendon and aponeuroses that are likely to have an effect on whole animal movement and energy expenditure, using the rat model system, Rattus norvegicus. We quantified muscle architecture (mass, cross-sectional area, and pennation angle), muscle fiber type and diameter, and Young's modulus of stiffness for the medial gastrocnemius aponeuroses as well as its contribution to Achilles tendon in three groups of three-month-old female rats: virgin, primiparous pregnant, and primiparous lactating animals. We found that muscle mass drops by 23% during lactation but does not change during pregnancy. We also found that during pregnancy muscle fibers switch from Type I to IIa and during lactation from Type IIb to Type I. The stiffness of connective tissues that has a demonstrated role in locomotion, the aponeurosis and tendon, also changed. Pregnant animals had a significantly less stiff aponeurosis. However, tendon stiffness was most affected during lactation, with a significant drop in stiffness and interindividual variation. We propose that the energetic demands of locomotion may have driven the evolution of these anatomical changes in muscle-tendon units during pregnancy and lactation to ensure more energy can be allocated to fetal development and lactation.

Metabolomics Analysis Provides Novel Insights into the Difference in Meat Quality between Different Pig Breeds

The Chuanzang black (CB) pig is a new crossbred between Chinese local breeds and modern breeds. Here, we investigated the growth performance, plasma indexes, carcass traits, and meat quality characteristics of conventional DLY (Duroc × Landrace × Yorkshire) crossbreed and CB pigs. The LC-MS/MS-based metabolomics of pork from DLY and CB pigs, as well as the relationship between the changes in the metabolic spectrum and meat quality, were analyzed. In this study, CB pigs presented lower final body weight, average daily gain, carcass weight, and eye muscle area than DLY pigs (p ˂ 0.05). Conversely, the ratio of feed to gain, marbling score, and meat color score of longissimus dorsi (LD) were higher in CB than DLY pigs (p ˂ 0.05). Moreover, psoas major (PM) showed a higher meat color score and a lower cooking loss in CB than DLY pigs (p ˂ 0.05). Interestingly, CB pigs showed lower myofiber diameter and area but higher myofiber density than DLY pigs (p ˂ 0.05). Furthermore, the mRNA expression levels of MyHC I, PPARδ, MEF2C, NFATC1, and AMPKα1 were higher in CB than DLY pigs (p ˂ 0.05). Importantly, a total of 753 metabolites were detected in the two tissues (e.g., psoas major and longissimus dorsi) of CB and DLY pigs, of which the difference in metabolite profiles in psoas major between crossbreeds was greater than that in longissimus dorsi. Specifically, palmitic acid, stearic acid, L-aspartic acid, corticosterone, and tetrahydrocorticosterone were the most relevant metabolites of psoas major meat quality, and tetrahydrocorticosterone, L-Palmitoylcarnitine, arachidic acid, erucic acid, and 13Z,16Z-docosadienoic acid in longissimus dorsi meat were positively correlated with meat quality. The most significantly enriched KEGG pathways in psoas major and longissimus dorsi pork were galactose metabolism and purine metabolism, respectively. These results not only indicated improved meat quality in CB pigs as compared to DLY pigs but may also assist in rational target selection for nutritional intervention or genetic breeding in the swine industry.

Large Maf transcription factor family is a major regulator of fast type IIb myofiber determination

Myofibers are broadly characterized as fatigue-resistant slow-twitch (type I) fibers and rapidly fatiguing fast-twitch (type IIa/IIx/IIb) fibers. However, the molecular regulation of myofiber type is not entirely understood; particularly, information on regulators of fast-twitch muscle is scarce. Here, we demonstrate that the large Maf transcription factor family dictates fast type IIb myofiber specification in mice. Remarkably, the ablation of three large Mafs leads to the drastic loss of type IIb myofibers, resulting in enhanced endurance capacity and the reduction of muscle force. Conversely, the overexpression of each large Maf in the type I soleus muscle induces type IIb myofibers. Mechanistically, a large Maf directly binds to the Maf recognition element on the promoter of myosin heavy chain 4, which encodes the type IIb myosin heavy chain, driving its expression. This work identifies the large Maf transcription factor family as a major regulator for fast type IIb muscle determination.

Lunar gravity prevents skeletal muscle atrophy but not myofiber type shift in mice

Skeletal muscle is sensitive to gravitational alterations. We recently developed a multiple artificial-gravity research system (MARS), which can generate gravity ranging from microgravity to Earth gravity (1 g) in space. Using the MARS, we studied the effects of three different gravitational levels (microgravity, lunar gravity [1/6 g], and 1 g) on the skeletal muscle mass and myofiber constitution in mice. All mice survived and returned to Earth, and skeletal muscle was collected two days after landing. We observed that microgravity-induced soleus muscle atrophy was prevented by lunar gravity. However, lunar gravity failed to prevent the slow-to-fast myofiber transition in the soleus muscle in space. These results suggest that lunar gravity is enough to maintain proteostasis, but a greater gravitational force is required to prevent the myofiber type transition. Our study proposes that different gravitational thresholds may be required for skeletal muscle adaptation.

Extracellular Matrix and Protein Phosphorylation Dysregulation Related to Diabetes-Induced Erectile Dysfunction

Diabetes can cause erectile dysfunction (ED) in more than half of male patients. However, the mechanisms underlying diabetes-induced erectile dysfunction (DED) remain unknown. This study is aimed at systematically analyzing the cellular and molecular mechanisms leading to DED using bioinformatic analysis and providing molecular targets for predicting and treating DED. In total, we identified 800 DEGs in the DED samples compared with those in the control group. The 407 upregulated DEGs were mainly enriched in glucose and lipid metabolism-related pathways, and the 393 downregulated DEGs were primarily enriched in tissue development and structure. Dysregulated extracellular matrix genes (especially collagen and elastin) may be closely related to damage to the erectile function of the corpus cavernosum. Sixteen hub genes and 24 modules were detected with hub genes and MCODE analysis. The consensus sequence AAA (G/C) AAA was observed at the promoter sites of most genes that were enriched in the “posttranslational protein phosphorylation” pathway. These genes had abundant phosphorylation sites. Furthermore, 20 TFs targeting DEGs were identified using ChEA3 tool. In conclusion, our research comprehensively and systematically describes the molecular characteristics of DED and suggests that dysregulated extracellular matrix genes and protein phosphorylation may play critical roles in DED. Therefore, they may be potential markers for diagnosing and treating DED.

KDM8 epigenetically controls cardiac metabolism to prevent initiation of dilated cardiomyopathy

Cardiac metabolism is deranged in heart failure, but underlying mechanisms remain unclear. Here, we show that lysine demethylase 8 (Kdm8) maintains an active mitochondrial gene network by repressing Tbx15, thus preventing dilated cardiomyopathy leading to lethal heart failure. Deletion of Kdm8 in mouse cardiomyocytes increased H3K36me2 with activation of Tbx15 and repression of target genes in the NAD+ pathway before dilated cardiomyopathy initiated. NAD+ supplementation prevented dilated cardiomyopathy in Kdm8 mutant mice, and TBX15 overexpression blunted NAD+-activated cardiomyocyte respiration. Furthermore, KDM8 was downregulated in human hearts affected by dilated cardiomyopathy, and higher TBX15 expression defines a subgroup of affected hearts with the strongest downregulation of genes encoding mitochondrial proteins. Thus, KDM8 represses TBX15 to maintain cardiac metabolism. Our results suggest that epigenetic dysregulation of metabolic gene networks initiates myocardium deterioration toward heart failure and could underlie heterogeneity of dilated cardiomyopathy.

Skeletal muscle structure, physiology, and function

Contractions of skeletal muscles provide the stability and power for all body movements. Consequently, any impairment in skeletal muscle function results in some degree of instability or immobility. Factors that influence skeletal muscle structure and function are therefore of great interest scientifically and clinically. Injury, neuromuscular disease, and old age are among the factors that commonly contribute to impairments in skeletal muscle function. The goal of this chapter is to summarize the fundamentals of skeletal muscle structure and function to provide foundational knowledge for this Handbook volume. We examine the molecular interactions that provide the basis for the generation of force and movement, discuss mechanisms of the regulation of contraction at the level of myofibers, and introduce concepts of the activation and control of muscle function in vivo. Where appropriate, the chapter updates the emerging science that will increase understanding of muscle function.

Spatially resolved transcriptomics reveals innervation-responsive functional clusters in skeletal muscle

Striated muscle is a highly organized structure composed of well-defined anatomical domains with integrated but distinct assignments. So far, the lack of a direct correlation between tissue architecture and gene expression has limited our understanding of how each unit responds to physio-pathologic contexts. Here, we show how the combined use of spatially resolved transcriptomics and immunofluorescence can bridge this gap by enabling the unbiased identification of such domains and the characterization of their response to external perturbations. Using a spatiotemporal analysis, we follow changes in the transcriptome of specific domains in muscle in a model of denervation. Furthermore, our approach enables us to identify the spatial distribution and nerve dependence of atrophic signaling pathway and polyamine metabolism to glycolytic fibers. Indeed, we demonstrate that perturbations of polyamine pathway can affect muscle function. Our dataset serves as a resource for future studies of the mechanisms underlying skeletal muscle homeostasis and innervation.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

Promoter-Adjacent DNA Hypermethylation Can Downmodulate Gene Expression: TBX15 in the Muscle Lineage

TBX15, which encodes a differentiation-related transcription factor, displays promoter-adjacent DNA hypermethylation in myoblasts and skeletal muscle (psoas) that is absent from non-expressing cells in other lineages. By whole-genome bisulfite sequencing (WGBS) and enzymatic methyl-seq (EM-seq), these hypermethylated regions were found to border both sides of a constitutively unmethylated promoter. To understand the functionality of this DNA hypermethylation, we cloned the differentially methylated sequences (DMRs) in CpG-free reporter vectors and tested them for promoter or enhancer activity upon transient transfection. These cloned regions exhibited strong promoter activity and, when placed upstream of a weak promoter, strong enhancer activity specifically in myoblast host cells. In vitro CpG methylation targeted to the DMR sequences in the plasmids resulted in 86−100% loss of promoter or enhancer activity, depending on the insert sequence. These results as well as chromatin epigenetic and transcription profiles for this gene in various cell types support the hypothesis that DNA hypermethylation immediately upstream and downstream of the unmethylated promoter region suppresses enhancer/extended promoter activity, thereby downmodulating, but not silencing, expression in myoblasts and certain kinds of skeletal muscle. This promoter-border hypermethylation was not found in cell types with a silent TBX15 gene, and these cells, instead, exhibit repressive chromatin in and around the promoter. TBX18, TBX2, TBX3 and TBX1 display TBX15-like hypermethylated DMRs at their promoter borders and preferential expression in myoblasts. Therefore, promoter-adjacent DNA hypermethylation for downmodulating transcription to prevent overexpression may be used more frequently for transcription regulation than currently appreciated.

Integrative analysis deciphers the heterogeneity of cancer-associated fibroblast and implications on clinical outcomes in ovarian cancers

Accumulating evidence has recognized that cancer-associated fibroblasts (CAFs) are major players in the desmoplastic stroma of ovarian cancer, modulating tumor progression and therapeutic response. However, it is unclear regarding the signatures of CAFs could be utilized to predict the clinical outcomes of ovarian cancer patients. To fill in this gap, we explored the intratumoral compartment of ovarian cancer by analyzing the single-cell RNA-sequencing (scRNA-seq) datasets of ovarian carcinoma samples, and identified two distinct CAFs (tumor-promoting CAF_c1 subtype and myofibroblasts-like CAF_c2 subtype). The clinical significance of CAF subtypes was further validated in The Cancer Genomics Atlas (TCGA) database and other independent immunotherapy response datasets, and the results revealed that the patients with a higher expression of CAF_c1 signatures had a worse prognosis and showed a tendency of resistance to immunotherapy. This work uncovered the signatures of the CAF_c1 subtype that could serve as a novel prognostic indicator and predictive marker for immunotherapy.

Population structure and genomic footprints of selection in five major Iranian horse breeds

The genetic structure and characteristics of Iranian native breeds are yet to be comprehensibly investigated and studied. Therefore, we employed genomic information of 364 Iranian native horses representing the Asil (n = 109), Caspian (n = 40), Dareshuri (n = 44), Kurdish (n = 95), and Turkoman (n = 76) breeds to reveal the genetic structure and characteristics. For these and 19 other horse breeds, principal component analysis, Bayesian model-based, Neighbor-Net, and bootstrap-based TreeMix approaches were applied to investigate and compare their genetic structure. Additionally, three haplotype-based methods including haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length were applied to trace genomic footprints of selection of Asil, Caspian, Dareshuri, Kurdish, and Turkoman groups. Then, the Mahalanobis distance based on the negative-log₁₀ rank-based P-values was estimated based on the haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length values. Asil, Caspian, Dareshuri, Kurdish, and Turkoman can be categorized into five different genetic clusters. Based on the top 1% of Mahalanobis distance based on the negative-log₁₀ rank-based P-values of SNPs, we identified 24 SNPs formerly reported to be associated with different traits and >100 genes undergoing selection pressures in Asil, Caspian, Dareshuri, Kurdish, and Turkoman. The detected QTL undergoing selection pressures were associated with withers height, equine metabolic syndrome, overall body size, insect bite hypersensitivity, guttural pouch tympany, white markings, Rhodococcus equi infection, jumping test score, alternate gaits, and body weight traits. Our findings will aid to have a better perspective of the genetic characteristics and population structure of Asil, Caspian, Dareshuri, Kurdish, and Turkoman horses as Iranian native horse breeds.

Transcriptome study digs out BMP2 involved in adipogenesis in sheep tails

Background

Hu sheep and Tibetan sheep in China are characterized by fat tails and thin tails, respectively. Several transcriptomes have been conducted in different sheep breeds to identify the differentially expressed genes (DEGs) underlying this trait. However, these studies identified different DEGs in different sheep breeds.

Results

Hence, RNA sequencing was performed on Hu sheep and Tibetan sheep. We obtained a total of 45.57 and 43.82 million sequencing reads, respectively. Two libraries mapped reads from 36.93 and 38.55 million reads after alignment to the reference sequences. 2108 DEGs were identified, including 1247 downregulated and 861 upregulated DEGs. GO and KEGG analyses of all DEGs demonstrated that pathways were enriched in the regulation of lipolysis in adipocytes and terms related to the chemokine signalling pathway, lysosomes, and glycosaminoglycan degradation. Eight genes were selected for validation by RT–qPCR. In addition, the transfection of BMP2 overexpression into preadipocytes resulted in increased PPAR-γ expression and expression. BMP2 potentially induces adipogenesis through LOX in preadipocytes. The number of lipid drops in BMP2 overexpression detected by oil red O staining was also greater than that in the negative control.

Conclusion

In summary, these results showed that significant genes (BMP2, HOXA11, PPP1CC and LPIN1) are involved in the regulation of adipogenesis metabolism and suggested novel insights into metabolic molecules in sheep fat tails.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08657-8.

Regulation of the evolutionarily conserved muscle myofibrillar matrix by cell type dependent and independent mechanisms

Skeletal muscles play a central role in human movement through forces transmitted by contraction of the sarcomere. We recently showed that mammalian sarcomeres are connected through frequent branches forming a singular, mesh-like myofibrillar matrix. However, the extent to which myofibrillar connectivity is evolutionarily conserved as well as mechanisms which regulate the specific architecture of sarcomere branching remain unclear. Here, we demonstrate the presence of a myofibrillar matrix in the tubular, but not indirect flight (IF) muscles within Drosophila melanogaster. Moreover, we find that loss of transcription factor H15 increases sarcomere branching frequency in the tubular jump muscles, and we show that sarcomere branching can be turned on in IF muscles by salm-mediated conversion to tubular muscles. Finally, we demonstrate that neurochondrin misexpression results in myofibrillar connectivity in IF muscles without conversion to tubular muscles. These data indicate an evolutionarily conserved myofibrillar matrix regulated by both cell-type dependent and independent mechanisms.

Harnessing tissue-specific genetic variation to dissect putative causal pathways between body mass index and cardiometabolic phenotypes

Body mass index (BMI) is a complex disease risk factor known to be influenced by genes acting via both metabolic pathways and appetite regulation. In this study, we aimed to gain insight into the phenotypic consequences of BMI-associated genetic variants, which may be mediated by their expression in different tissues. First, we harnessed meta-analyzed gene expression datasets derived from subcutaneous adipose (n = 1257) and brain (n = 1194) tissue to identify 86 and 140 loci, respectively, which provided evidence of genetic colocalization with BMI. These two sets of tissue-partitioned loci had differential effects with respect to waist-to-hip ratio, suggesting that the way they influence fat distribution might vary despite their having very similar average magnitudes of effect on BMI itself (adipose = 0.0148 and brain = 0.0149 standard deviation change in BMI per effect allele). For instance, BMI-associated variants colocalized with TBX15 expression in adipose tissue (posterior probability [PPA] = 0.97), but not when we used TBX15 expression data derived from brain tissue (PPA = 0.04) This gene putatively influences BMI via its role in skeletal development. Conversely, there were loci where BMI-associated variants provided evidence of colocalization with gene expression in brain tissue (e.g., NEGR1, PPA = 0.93), but not when we used data derived from adipose tissue, suggesting that these genes might be more likely to influence BMI via energy balance. Leveraging these tissue-partitioned variant sets through a multivariable Mendelian randomization framework provided strong evidence that the brain-tissue-derived variants are predominantly responsible for driving the genetically predicted effects of BMI on cardiovascular-disease endpoints (e.g., coronary artery disease: odds ratio = 1.05, 95% confidence interval = 1.04-1.07, p = 4.67 × 10^-14). In contrast, our analyses suggested that the adipose tissue variants might predominantly be responsible for the underlying relationship between BMI and measures of cardiac function, such as left ventricular stroke volume (beta = 0.21, 95% confidence interval = 0.09-0.32, p = 6.43 × 10^-4).

Zinc deficiency and supplementation in autism spectrum disorder and Phelan-McDermid syndrome

Approximately 1 in 36 children are diagnosed with autism spectrum disorder (ASD). The disorder is four times more common in males than in females. Zinc deficiency and mutations in SHANK2 and SHANK3 (members of a family of excitatory postsynaptic scaffolding proteins) are all risk factors that may contribute to the pathophysiology of the disease. The presence of shankopathies (loss of one copy of the SHANK3 gene) can lead to the development of Phelan-McDermid syndrome (PMDS)-a rare genetic disorder characterized by developmental delay, intellectual disability, poor motor tone, and ASD-like symptoms. We reviewed the relationship between zinc, ASD, and PMDS as well as the effect of zinc supplementation in improving symptoms of ASD and PMDS based on 22 studies published within 6 years (2015-2020). Zinc deficiency (assessed by either dietary intake, blood, hair, or tooth matrix) was shown to be highly prevalent in ASD and PMDS patients as well as in preclinical models of ASD and PMDS. Zinc supplements improved the behavioral deficits in animal models of ASD and PMDS. Clinical trials are still needed to validate the beneficial therapeutic effects of zinc supplements in ASD and PMDS patients.

Genome Wide Association Study of Beef Traits in Local Alpine Breed Reveals the Diversity of the Pathways Involved and the Role of Time Stratification

Knowledge of the genetic architecture of key growth and beef traits in livestock species has greatly improved worldwide thanks to genome-wide association studies (GWAS), which allow to link target phenotypes to Single Nucleotide Polymorphisms (SNPs) across the genome. Local dual-purpose breeds have rarely been the focus of such studies; recently, however, their value as a possible alternative to intensively farmed breeds has become clear, especially for their greater adaptability to environmental change and potential for survival in less productive areas. We performed single-step GWAS and post-GWAS analysis for body weight (BW), average daily gain (ADG), carcass fleshiness (CF) and dressing percentage (DP) in 1,690 individuals of local alpine cattle breed, Rendena. This breed is typical of alpine pastures, with a marked dual-purpose attitude and good genetic diversity. Moreover, we considered two of the target phenotypes (BW and ADG) at different times in the individuals' life, a potentially important aspect in the study of the traits' genetic architecture. We identified 8 significant and 47 suggestively associated SNPs, located in 14 autosomal chromosomes (BTA). Among the strongest signals, 3 significant and 16 suggestive SNPs were associated with ADG and were located on BTA10 (50-60 Mb), while the hotspot associated with CF and DP was on BTA18 (55-62 MB). Among the significant SNPs some were mapped within genes, such as SLC12A1, CGNL1, PRTG (ADG), LOC513941 (CF), NLRP2 (CF and DP), CDC155 (DP). Pathway analysis showed great diversity in the biological pathways linked to the different traits; several were associated with neurogenesis and synaptic transmission, but actin-related and transmembrane transport pathways were also represented. Time-stratification highlighted how the genetic architectures of the same traits were markedly different between different ages. The results from our GWAS of beef traits in Rendena led to the detection of a variety of genes both well-known and novel. We argue that our results show that expanding genomic research to local breeds can reveal hitherto undetected genetic architectures in livestock worldwide. This could greatly help efforts to map genomic complexity of the traits of interest and to make appropriate breeding decisions.

Skeletal Muscle Metabolic Alternation Develops Sarcopenia

Sarcopenia is a new type of senile syndrome with progressive skeletal muscle mass loss with age, accompanied by decreased muscle strength and/or muscle function. Sarcopenia poses a serious threat to the health of the elderly and increases the burden of family and society. The underlying pathophysiological mechanisms of sarcopenia are still unclear. Recent studies have shown that changes of skeletal muscle metabolism are the risk factors for sarcopenia. Furthermore, the importance of the skeletal muscle metabolic microenvironment in regulating satellite cells (SCs) is gaining significant attention. Skeletal muscle metabolism has intrinsic relationship with the regulation of skeletal muscle mass and regeneration. This review is to discuss recent findings regarding skeletal muscle metabolic alternation and the development of sarcopenia, hoping to contribute better understanding and treatment of sarcopenia.

Skin aging: Dermal adipocytes metabolically reprogram dermal fibroblasts

Emerging data connects the aging process in dermal fibroblasts with metabolic reprogramming, provided by enhanced fatty acid oxidation and reduced glycolysis. This switch may be caused by a significant expansion of the dermal white adipose tissue (dWAT) layer in aged, hair-covered skin. Dermal adipocytes cycle through de-differentiation and re-differentiation. As a result, there is a strongly enhanced release of free fatty acids into the extracellular space during the de-differentiation of dermal adipocytes in the catagen phase of the hair follicle cycle. Both caveolin-1 and adiponectin are critical factors influencing these processes. Controlling the expression levels of these two factors also offers the ability to manipulate the metabolic preferences of the different cell types within the microenvironment of the skin, including dermal fibroblasts. Differential expression of adiponectin and caveolin-1 in the various cell types may also be responsible for the cellular metabolic heterogeneity within the cells of the skin.

The effects of Tbx15 and Pax1 on facial and other physical morphology in mice

DNA variants in or close to the human TBX15 and PAX1 genes have been repeatedly associated with facial morphology in independent genome-wide association studies, while their functional roles in determining facial morphology remain to be understood. We generated Tbx15 knockout (Tbx15 -/-) and Pax1 knockout (Pax1 -/-) mice by applying the one-step CRISPR/Cas9 method. A total of 75 adult mice were used for subsequent phenotype analysis, including 38 Tbx15 mice (10 homozygous Tbx15 -/-, 18 heterozygous Tbx15 +/-, 10 wild-type Tbx15 +/+ WT littermates) and 37 Pax1 mice (12 homozygous Pax1 -/-, 15 heterozygous Pax1 +/-, 10 Pax1 +/+ WT littermates). Facial and other physical morphological phenotypes were obtained from three-dimensional (3D) images acquired with the HandySCAN BLACK scanner. Compared to WT littermates, the Tbx15 -/- mutant mice had significantly shorter faces (p = 1.08E-8, R2 = 0.61) and their ears were in a significantly lower position (p = 3.54E-8, R2 = 0.62) manifesting a "droopy ear" characteristic. Besides these face alternations, Tbx15 -/- mutant mice displayed significantly lower weight as well as shorter body and limb length. Pax1 -/- mutant mice showed significantly longer noses (p = 1.14E-5, R2 = 0.46) relative to WT littermates, but otherwise displayed less obvious morphological alterations than Tbx15 -/- mutant mice did. We provide the first direct functional evidence that two well-known and replicated human face genes, Tbx15 and Pax1, impact facial and other body morphology in mice. The general agreement between our findings in knock-out mice with those from previous GWASs suggests that the functional evidence we established here in mice may also be relevant in humans.

Effect of dietary leucine supplementation on skeletal muscle fiber type transformation in weaning piglets

To investigate the effects of dietary leucine supplementation on muscle fiber type transformation in weaning piglets, 54 21-day-old male DLY (Duroc × Landrace × Yorkshire) weaned piglets were randomly divided into control, 0.25% and 0.5% leucine groups. The experiment lasted for 42 d. The results showed that dietary supplementation of 0.25% leucine significantly increased the protein expressions of slow MyHC, myoglobin and Troponin I-SS and the mRNA expressions of MyHC I, MyHC IIa, Tnni1, Tnnc1, Tnnt1 and myoglobin, while decreased the protein level of fast MyHC and the mRNA level of MyHC IIb in longissimus dorsi (LD) muscle. Furthermore, 0.25% leucine significantly increased succinic dehydrogenase (SDH) activity and decreased lactate dehydrogenase (LDH) activity. In addition, our data found that 0.25% leucine significantly increased serum adiponectin (AdipoQ) concentration, and the protein levels of AdipoQ, adiponectin receptor 1 (AdipoR1), phosphorylated AMP-activated protein kinase (p-AMPK) and PPAR-γ coactivator-1α (PGC-1α) and the mRNA levels of AdipoQ, AdipoR1 and AMPKα2. Together, our findings indicate that leucine promotes porcine skeletal muscle fiber type transformation from fast-twitch to slow-twitch, and the effect may be mediated by AdipoQ-AMPK-PGC-1α signaling pathway.

Association Study among Comethylation Modules, Genetic Polymorphisms and Clinical Features in Mexican Teenagers with Eating Disorders: Preliminary Results

Eating disorders are psychiatric disorders characterized by disturbed eating behaviors. They have a complex etiology in which genetic and environmental factors interact. Analyzing gene-environment interactions could help us to identify the mechanisms involved in the etiology of such conditions. For example, comethylation module analysis could detect the small effects of epigenetic interactions, reflecting the influence of environmental factors. We used MethylationEPIC and Psycharray microarrays to determine DNA methylation levels and genotype from 63 teenagers with eating disorders. We identified 11 comethylation modules in WGCNA (Weighted Gene Correlation Network Analysis) and correlated them with single nucleotide polymorphisms (SNP) and clinical features in our subjects. Two comethylation modules correlated with clinical features (BMI and height) in our sample and with SNPs associated with these phenotypes. One of these comethylation modules (yellow) correlated with BMI and rs10494217 polymorphism (associated with waist-hip ratio). Another module (black) was correlated with height, rs9349206, rs11761528, and rs17726787 SNPs; these polymorphisms were associated with height in previous GWAS. Our data suggest that genetic variations could alter epigenetics, and that these perturbations could be reflected as variations in clinical features.

PROKR1 delivery by cell-derived vesicles restores the myogenic potential of Prokr1-deficient C2C12 myoblasts

Cell-derived vesicles (CDVs) have been investigated as an alternative to exosomes. Here, we generated CDVs from Prokineticin receptor 1 (PROKR1) overexpressing HEK293T cells using micro-extrusion. More than 60 billion PROKR1-enriched CDV (PROKR1^Tg CDVs) particles with canonical exosome properties were recovered from 10⁷ cells. With 25 μg/mL of PROKR1^Tg CDVs, we observed delivery of PROKR1, significant reduction of apoptosis, and myotube formation in C2C12^Prokr1-/- myoblasts that have lost their myogenic potential but underwent apoptosis following myogenic commitment. Expression levels of early and late myogenic marker genes and glucose uptake capacity were restored to equivalent levels with wild-type control. Furthermore, PROKR1^Tg CDVs were accumulated in soleus muscle comparable to the liver without significant differences. Therefore, CDVs obtained from genetically engineered cells appear to be an effective method of PROKR1 protein delivery and offer promise as an alternative therapy for muscular dystrophy.

Histone methyltransferase MLL4 controls myofiber identity and muscle performance through MEF2 interaction

Skeletal muscle depends on the precise orchestration of contractile and metabolic gene expression programs to direct fiber-type specification and to ensure muscle performance. Exactly how such fiber type-specific patterns of gene expression are established and maintained remains unclear, however. Here, we demonstrate that histone monomethyl transferase MLL4 (KMT2D), an enhancer regulator enriched in slow myofibers, plays a critical role in controlling muscle fiber identity as well as muscle performance. Skeletal muscle-specific ablation of MLL4 in mice resulted in downregulation of the slow oxidative myofiber gene program, decreased numbers of type I myofibers, and diminished mitochondrial respiration, which caused reductions in muscle fatty acid utilization and endurance capacity during exercise. Genome-wide ChIP-Seq and mRNA-Seq analyses revealed that MLL4 directly binds to enhancers and functions as a coactivator of the myocyte enhancer factor 2 (MEF2) to activate transcription of slow oxidative myofiber genes. Importantly, we also found that the MLL4 regulatory circuit is associated with muscle fiber-type remodeling in humans. Thus, our results uncover a pivotal role for MLL4 in specifying structural and metabolic identities of myofibers that govern muscle performance. These findings provide therapeutic opportunities for enhancing muscle fitness to combat a variety of metabolic and muscular diseases.

A maternal high-fat/low-fiber diet impairs glucose tolerance and induces the formation of glycolytic muscle fibers in neonatal offspring

PURPOSE

In our previous study, the maternal high-fat/low-fiber (HF-LF) diet was suggested to induce metabolic disorders and placental dysfunction of the dam, but the effects of this diet on glucose metabolism of neonatal offspring remain largely unknown. Here, a neonatal pig model was used to evaluate the effects of maternal HF-LF diet during pregnancy on glucose tolerance, transition of skeletal muscle fiber types, and mitochondrial function in offspring.

METHODS

A total of 66 pregnant gilts (Guangdong Small-ear Spotted pig) at day 60 of gestation were randomly divided into two groups: control group (CON group; 2.86% crude fat, 9.37% crude fiber), and high-fat/low-fiber diet group (HF-LF group; 5.99% crude fat, 4.13% crude fiber).

RESULTS

The maternal HF-LF diet was shown to impair the glucose tolerance of neonatal offspring, downregulate the protein level of slow-twitch fiber myosin heavy chain I (MyHC I), and upregulate the protein levels of fast-twitch fiber myosin heavy chain IIb (MyHC IIb) and IIx (MyHC IIx) in soleus muscle. Additionally, compared with the CON group, the HF-LF offspring showed inhibition of insulin signaling pathway and decrease in mitochondrial function in liver and soleus muscle.

CONCLUSION

Maternal HF-LF diet during pregnancy impairs glucose tolerance, induces the formation of glycolytic muscle fibers, and decreases the hepatic and muscular mitochondrial function in neonatal piglets.

Cell autonomous requirement of neurofibromin (Nf1) for postnatal muscle hypertrophic growth and metabolic homeostasis

BACKGROUND

Neurofibromatosis type 1 (NF1) is a multi-organ disease caused by mutations in neurofibromin 1 (NF1). Amongst other features, NF1 patients frequently show reduced muscle mass and strength, impairing patients' mobility and increasing the risk of fall. The role of Nf1 in muscle and the cause for the NF1-associated myopathy are mostly unknown.

METHODS

To dissect the function of Nf1 in muscle, we created muscle-specific knockout mouse models for NF1, inactivating Nf1 in the prenatal myogenic lineage either under the Lbx1 promoter or under the Myf5 promoter. Mice were analysed during prenatal and postnatal myogenesis and muscle growth.

RESULTS

Nf1^Lbx1 and Nf1^Myf5 animals showed only mild defects in prenatal myogenesis. Nf1^Lbx1 animals were perinatally lethal, while Nf1^Myf5 animals survived only up to approximately 25 weeks. A comprehensive phenotypic characterization of Nf1^Myf5 animals showed decreased postnatal growth, reduced muscle size, and fast fibre atrophy. Proteome and transcriptome analyses of muscle tissue indicated decreased protein synthesis and increased proteasomal degradation, and decreased glycolytic and increased oxidative activity in muscle tissue. High-resolution respirometry confirmed enhanced oxidative metabolism in Nf1^Myf5 muscles, which was concomitant to a fibre type shift from type 2B to type 2A and type 1. Moreover, Nf1^Myf5 muscles showed hallmarks of decreased activation of mTORC1 and increased expression of atrogenes. Remarkably, loss of Nf1 promoted a robust activation of AMPK with a gene expression profile indicative of increased fatty acid catabolism. Additionally, we observed a strong induction of genes encoding catabolic cytokines in muscle Nf1^Myf5 animals, in line with a drastic reduction of white, but not brown adipose tissue.

CONCLUSIONS

Our results demonstrate a cell autonomous role for Nf1 in myogenic cells during postnatal muscle growth required for metabolic and proteostatic homeostasis. Furthermore, Nf1 deficiency in muscle drives cross-tissue communication and mobilization of lipid reserves.

Muscle-Specific Insulin Receptor Overexpression Protects Mice From Diet-Induced Glucose Intolerance but Leads to Postreceptor Insulin Resistance

Skeletal muscle insulin resistance is a prominent early feature in the pathogenesis of type 2 diabetes. In attempt to overcome this defect, we generated mice overexpressing insulin receptors (IR) specifically in skeletal muscle (IRMOE). On normal chow, IRMOE mice have body weight similar to that of controls but an increase in lean mass and glycolytic muscle fibers and reduced fat mass. IRMOE mice also show higher basal phosphorylation of IR, IRS-1, and Akt in muscle and improved glucose tolerance compared with controls. When challenged with high-fat diet (HFD), IRMOE mice are protected from diet-induced obesity. This is associated with reduced inflammation in fat and liver, improved glucose tolerance, and improved systemic insulin sensitivity. Surprisingly, however, in both chow and HFD-fed mice, insulin-stimulated Akt phosphorylation is significantly reduced in muscle of IRMOE mice, indicating postreceptor insulin resistance. RNA sequencing reveals downregulation of several postreceptor signaling proteins that contribute to this resistance. Thus, enhancing early insulin signaling in muscle by overexpression of the IR protects mice from diet-induced obesity and its effects on glucose metabolism. However, chronic overstimulation of this pathway leads to postreceptor desensitization, indicating the critical balance between normal signaling and hyperstimulation of the insulin signaling pathway.

Single-nucleus RNA-seq and FISH identify coordinated transcriptional activity in mammalian myofibers

Skeletal muscle fibers are large syncytia but it is currently unknown whether gene expression is coordinately regulated in their numerous nuclei. Here we show by snRNA-seq and snATAC-seq that slow, fast, myotendinous and neuromuscular junction myonuclei each have different transcriptional programs, associated with distinct chromatin states and combinations of transcription factors. In adult mice, identified myofiber types predominantly express either a slow or one of the three fast isoforms of Myosin heavy chain (MYH) proteins, while a small number of hybrid fibers can express more than one MYH. By snRNA-seq and FISH, we show that the majority of myonuclei within a myofiber are synchronized, coordinately expressing only one fast Myh isoform with a preferential panel of muscle-specific genes. Importantly, this coordination of expression occurs early during post-natal development and depends on innervation. These findings highlight a previously undefined mechanism of coordination of gene expression in a syncytium.

Whether skeletal muscle fibre gene expression is coordinated as a whole in different nuclei in the fibre is unclear. Here, the authors use single nucleus RNAseq and ATACseq to show the transcriptome heterogeneity of muscle nuclei in the adult mouse fibre, with correlations between the two datasets.

Effects of lycopene on skeletal muscle-fiber type and high-fat diet-induced oxidative stress

Increasing studies report that many natural products can participate in formation of muscle fibers. This study aimed to investigate the effect of lycopene on skeletal muscle-fiber type in vivo and in vitro. C2C12 myoblasts were used in vitro study, and the concentration of lycopene was 10 µM. In vivo, 8-week-old male C57/BL6 mice were used and divided into four groups (n=8): (1) ND: normal-fat diet; (2) ND+Lyc: normal-fat diet mixed with 0.33% w/w lycopene; (3) HFD: high-fat diet; and (4) HFD+Lyc: high-fat diet mixed with 0.33% w/w lycopene. The mice tissue samples were collected after 8 weeks feeding. We found that lycopene supplementation enhanced the protein expression of slow-twitch fiber, succinate dehydrogenase, and malic dehydrogenase enzyme activities, whereas lycopene reduced the protein expression of fast-twitch fibers, lactate dehydrogenase, pyruvate kinase enzyme activities. Moreover, lycopene can promote skeletal muscle triglyceride deposition, enhanced the mRNA expression of genes related to lipid synthesis, reduced the mRNA expression of genes related to lipolysis. And high-fat diet-induced dyslipidemia and oxidative stress were attenuated after lycopene supplementation. Additionally, lycopene supplementation reduced the glycolytic reserve but enhanced mitochondrial ATP production in C2C12 cells. These results demonstrated that lycopene affects the activities of metabolic enzymes in muscle fibers, promotes the expression of slow-twitch fibers, and enhanced mitochondrial respiratory capacity. We speculated that lycopene affects the muscle-fiber type through aerobic oxidation, suggesting that lycopene exerts potential beneficial effects on skeletal muscle metabolism.

What do warming waters mean for fish physiology and fisheries?

Environmental signals act primarily on physiological systems, which then influence higher-level functions such as movement patterns and population dynamics. Increases in average temperature and temperature variability associated with global climate change are likely to have strong effects on fish physiology and thereby on populations and fisheries. Here we review the principal mechanisms that transduce temperature signals and the physiological responses to those signals in fish. Temperature has a direct, thermodynamic effect on biochemical reaction rates. Nonetheless, plastic responses to longer-term thermal signals mean that fishes can modulate their acute thermal responses to compensate at least partially for thermodynamic effects. Energetics are particularly relevant for growth and movement, and therefore for fisheries, and temperature can have pronounced effects on energy metabolism. All energy (ATP) production is ultimately linked to mitochondria, and temperature has pronounced effects on mitochondrial efficiency and maximal capacities. Mitochondria are dependent on oxygen as the ultimate electron acceptor so that cardiovascular function and oxygen delivery link environmental inputs with energy metabolism. Growth efficiency, that is the conversion of food into tissue, changes with temperature, and there are indications that warmer water leads to decreased conversion efficiencies. Moreover, movement and migration of fish relies on muscle function, which is partially dependent on ATP production but also on intracellular calcium cycling within the myocyte. Neuroendocrine processes link environmental signals to regulated responses at the level of different tissues, including muscle. These physiological processes within individuals can scale up to population responses to climate change. A mechanistic understanding of thermal responses is essential to predict the vulnerability of species and populations to climate change.

Comparative population genomic analysis uncovers novel genomic footprints and genes associated with small body size in Chinese pony

Background

Body size is considered as one of the most fundamental properties of an organism. Due to intensive breeding and artificial selection throughout the domestication history, horses exhibit striking variations for heights at withers and body sizes. Debao pony (DBP), a famous Chinese horse, is known for its small body size and lives in Guangxi mountains of southern China. In this study, we employed comparative population genomics to study the genetic basis underlying the small body size of DBP breed based on the whole genome sequencing data. To detect genomic signatures of positive selection, we applied three methods based on population comparison, fixation index (F_ST), cross population composite likelihood ratio (XP-CLR) and nucleotide diversity (θπ), and further analyzed the results to find genomic regions under selection for body size-related traits.

Results

A number of protein-coding genes in windows with the top 1% values of F_ST (367 genes), XP-CLR (681 genes), and log₂ (θπ ratio) (332 genes) were identified. The most significant signal of positive selection was mapped to the NELL1 gene, probably underlies the body size and development traits, and may also have been selected for short stature in the DBP population. In addition, some other loci on different chromosomes were identified to be potentially involved in the development of body size.

Conclusions

Results of our study identified some positively selected genes across the horse genome, which are possibly involved in body size traits. These novel candidate genes may be useful targets for clarifying our understanding of the molecular basis of body size and as such they should be of great interest for future research into the genetic architecture of relevant traits in horse breeding program.

The importance of the ZBED6-IGF2 axis for metabolic regulation in mouse myoblast cells

The transcription factor ZBED6 acts as a repressor of Igf2 and affects directly or indirectly the transcriptional regulation of thousands of genes. Here, we use gene editing in mouse C2C12 myoblasts and show that ZBED6 regulates Igf2 exclusively through its binding site 5'-GGCTCG-3' in intron 1 of Igf2. Deletion of this motif (Igf2^ΔGGCT ) or complete ablation of Zbed6 leads to ~20-fold upregulation of the IGF2 protein. Quantitative proteomics revealed an activation of Ras signaling pathway in both Zbed6^-/- and Igf2^ΔGGCT myoblasts, and a significant enrichment of mitochondrial membrane proteins among proteins showing altered expression in Zbed6^-/- myoblasts. Both Zbed6^-/- and Igf2^ΔGGCT myoblasts showed a faster growth rate and developed myotube hypertrophy. These cells exhibited an increased O₂ consumption rate, due to IGF2 upregulation. Transcriptome analysis revealed ~30% overlap between differentially expressed genes in Zbed6^-/- and Igf2^ΔGGCT myotubes, with an enrichment of upregulated genes involved in muscle development. In contrast, ZBED6-overexpression in myoblasts led to cell apoptosis, cell cycle arrest, reduced mitochondrial activities, and ceased myoblast differentiation. The similarities in growth and differentiation phenotypes observed in Zbed6^-/- and Igf2^ΔGGCT myoblasts demonstrates that ZBED6 affects mitochondrial activity and myogenesis largely through its regulation of IGF2 expression. This study adds new insights how the ZBED6-Igf2 axis affects muscle metabolism.

Skeletal muscle enhancer interactions identify genes controlling whole-body metabolism

Obesity and type 2 diabetes (T2D) are metabolic disorders influenced by lifestyle and genetic factors that are characterized by insulin resistance in skeletal muscle, a prominent site of glucose disposal. Numerous genetic variants have been associated with obesity and T2D, of which the majority are located in non-coding DNA regions. This suggests that most variants mediate their effect by altering the activity of gene-regulatory elements, including enhancers. Here, we map skeletal muscle genomic enhancer elements that are dynamically regulated after exposure to the free fatty acid palmitate or the inflammatory cytokine TNFα. By overlapping enhancer positions with the location of disease-associated genetic variants, and resolving long-range chromatin interactions between enhancers and gene promoters, we identify target genes involved in metabolic dysfunction in skeletal muscle. The majority of these genes also associate with altered whole-body metabolic phenotypes in the murine BXD genetic reference population. Thus, our combined genomic investigations identified genes that are involved in skeletal muscle metabolism.

Obesity and type 2 diabetes (T2D) are metabolic disorders characterized by insulin resistance in skeletal muscle. Here, the authors map skeletal muscle enhancer elements dynamically regulated after exposure to free fatty acid palmitate or inflammatory cytokine TNFα and identify target genes involved in metabolic dysfunction in skeletal muscle.

Profiling and Functional Analysis of Circular RNAs in Porcine Fast and Slow Muscles

The different skeletal muscle fiber types exhibit distinctively different physiological and metabolic properties, and have been linked to both human metabolic diseases and meat quality traits in livestock. Circular RNAs (circRNAs) are a new class of endogenous RNA regulating gene expression, but regulatory mechanisms of skeletal muscle fibers involved in circRNAs remain poorly understood. Here, we constructed circRNA expression profiles of three fast-twitch biceps femoris (Bf) and three slow-twitch soleus (Sol) muscles in pigs using RNA-seq and identified 16,342 distinct circRNA candidates. Notably, 242 differentially expressed (DE) circRNAs between Bf and Sol muscles were identified, including 105 upregulated and 137 downregulated circRNAs, and are thus potential candidates for the regulation of skeletal muscle fiber conversion. Moreover, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of host genes of DE circRNAs revealed that host genes were mainly involved in skeletal muscle fiber-related GO terms (e.g., muscle contraction, contractile fiber part, and Z disk) and skeletal muscle fiber-related signaling pathways (e.g., AMPK and cGMP-PKG). We also constructed co-expression networks of DE circRNA-miRNA-mRNA using previously acquired high-throughput sequencing mRNA and miRNA data, from which 112 circRNA-miRNA and 95 miRNA-mRNA interactions were identified. Multiple circRNAs essentially serve as a sponge for miR-499-5p, which is preferentially expressed in slow-twitch muscle and reduces the severity of Duchenne muscular dystrophy (DMD). Taken together, a series of novel candidate circRNAs involved in the growth and development of porcine skeletal muscle was identified. Furthermore, they provide a comprehensive circRNA resource for further in-depth research on the regulatory mechanisms of circRNA in the formation of skeletal muscle fiber, and may provide insights into human skeletal muscle diseases.

Myosin Sequestration Regulates Sarcomere Function, Cardiomyocyte Energetics, and Metabolism, Informing the Pathogenesis of Hypertrophic Cardiomyopathy

BACKGROUND

Hypertrophic cardiomyopathy (HCM) is caused by pathogenic variants in sarcomere protein genes that evoke hypercontractility, poor relaxation, and increased energy consumption by the heart and increased patient risks for arrhythmias and heart failure. Recent studies show that pathogenic missense variants in myosin, the molecular motor of the sarcomere, are clustered in residues that participate in dynamic conformational states of sarcomere proteins. We hypothesized that these conformations are essential to adapt contractile output for energy conservation and that pathophysiology of HCM results from destabilization of these conformations.

METHODS

We assayed myosin ATP binding to define the proportion of myosins in the super relaxed state (SRX) conformation or the disordered relaxed state (DRX) conformation in healthy rodent and human hearts, at baseline and in response to reduced hemodynamic demands of hibernation or pathogenic HCM variants. To determine the relationships between myosin conformations, sarcomere function, and cell biology, we assessed contractility, relaxation, and cardiomyocyte morphology and metabolism, with and without an allosteric modulator of myosin ATPase activity. We then tested whether the positions of myosin variants of unknown clinical significance that were identified in patients with HCM, predicted functional consequences and associations with heart failure and arrhythmias.

RESULTS

Myosins undergo physiological shifts between the SRX conformation that maximizes energy conservation and the DRX conformation that enables cross-bridge formation with greater ATP consumption. Systemic hemodynamic requirements, pharmacological modulators of myosin, and pathogenic myosin missense mutations influenced the proportions of these conformations. Hibernation increased the proportion of myosins in the SRX conformation, whereas pathogenic variants destabilized these and increased the proportion of myosins in the DRX conformation, which enhanced cardiomyocyte contractility, but impaired relaxation and evoked hypertrophic remodeling with increased energetic stress. Using structural locations to stratify variants of unknown clinical significance, we showed that the variants that destabilized myosin conformations were associated with higher rates of heart failure and arrhythmias in patients with HCM.

CONCLUSIONS

Myosin conformations establish work-energy equipoise that is essential for life-long cellular homeostasis and heart function. Destabilization of myosin energy-conserving states promotes contractile abnormalities, morphological and metabolic remodeling, and adverse clinical outcomes in patients with HCM. Therapeutic restabilization corrects cellular contractile and metabolic phenotypes and may limit these adverse clinical outcomes in patients with HCM.

Ovarian cancer detection by DNA methylation in cervical scrapings

Background

Ovarian cancer (OC) is the most lethal gynecological cancer, worldwide, largely due to its vague and nonspecific early stage symptoms, resulting in most tumors being found at advanced stages. Moreover, due to its relative rarity, there are currently no satisfactory methods for OC screening, which remains a controversial and cost-prohibitive issue. Here, we demonstrate that Papanicolaou test (Pap test) cervical scrapings, instead of blood, can reveal genetic/epigenetic information for OC detection, using specific and sensitive DNA methylation biomarkers.

Results

We analyzed the methylomes of tissues (50 OC tissues versus 6 normal ovarian epithelia) and cervical scrapings (5 OC patients versus 10 normal controls), and integrated public methylomic datasets, including 79 OC tissues and 6 normal tubal epithelia. Differentially methylated genes were further classified by unsupervised hierarchical clustering, and each candidate biomarker gene was verified in both OC tissues and cervical scrapings by either quantitative methylation-specific polymerase chain reaction (qMSP) or bisulfite pyrosequencing. A risk-score by logistic regression was generated for clinical application.

One hundred fifty-one genes were classified into four clusters, and nine candidate hypermethylated genes from these four clusters were selected. Among these, four genes fulfilled our selection criteria and were validated in training and testing set, respectively. The OC detection accuracy was demonstrated by area under the receiver operating characteristic curves (AUCs) in 0.80–0.83 of AMPD3, 0.79–0.85 of AOX1, 0.78–0.88 of NRN1, and 0.82–0.85 of TBX15. From this, we found OC-risk score, equation generated by logistic regression in training set and validated an OC-associated panel comprising AMPD3, NRN1, and TBX15, reaching a sensitivity of 81%, specificity of 84%, and OC detection accuracy of 0.91 (95% CI, 0.82–1) in testing set.

Conclusions

Ovarian cancer detection from cervical scrapings is feasible, using particularly promising epigenetic biomarkers such as AMPD3/NRN1/TBX15. Further validation is warranted.

Dynamic enhancers control skeletal muscle identity and reprogramming

Skeletal muscles consist of fibers of differing metabolic activities and contractility, which become remodeled in response to chronic exercise, but the epigenomic basis for muscle identity and adaptation remains poorly understood. Here, we used chromatin immunoprecipitation sequencing of dimethylated histone 3 lysine 4 and acetylated histone 3 lysine 27 as well as transposase-accessible chromatin profiling to dissect cis-regulatory networks across muscle groups. We demonstrate that in vivo enhancers specify muscles in accordance with myofiber composition, show little resemblance to cultured myotube enhancers, and identify glycolytic and oxidative muscle-specific regulators. Moreover, we find that voluntary wheel running and muscle-specific peroxisome proliferator-activated receptor gamma coactivator-1 alpha (Pgc1a) transgenic (mTg) overexpression, which stimulate endurance performance in mice, result in markedly different changes to the epigenome. Exercise predominantly leads to enhancer hypoacetylation, whereas mTg causes hyperacetylation at different sites. Integrative analysis of regulatory regions and gene expression revealed that exercise and mTg are each associated with myocyte enhancer factor (MEF) 2 and estrogen-related receptor (ERR) signaling and transcription of genes promoting oxidative metabolism. However, exercise was additionally associated with regulation by retinoid X receptor (RXR), jun proto-oncogene (JUN), sine oculis homeobox factor (SIX), and other factors. Overall, our work defines the unique enhancer repertoires of skeletal muscles in vivo and reveals that divergent exercise-induced or PGC1α-driven epigenomic programs direct partially convergent transcriptional networks.

Tbx15 is required for adipocyte browning induced by adrenergic signaling pathway

OBJECTIVE

The T-box gene Tbx15 is abundantly expressed in adipose tissues, especially subcutaneous and brown fat. Although its expression is correlated with obesity, its precise biological role in adipose tissue is poorly understood in vivo. Here we investigated the function of Tbx15 in brown adipose thermogenesis and white adipose browning in vivo.

METHODS

In the present study, we generated adipose-specific Tbx15 knockout (AKO) mice by crossing Tbx15 floxed mice with adiponectin-Cre mice to delineate Tbx15 function in adipose tissues. We systematically investigated the influence of Tbx15 on brown adipose thermogenesis and white adipose browning in mice, as well as the possible underlying molecular mechanism.

RESULTS

Upon cold exposure, adipocyte browning in inguinal adipose tissue was significantly impaired in Tbx15 AKO mice. Furthermore, ablation of Tbx15 blocked adipocyte browning induced by β3 adrenergic agonist CL 316243, which did not appear to alter the expression of Tbx15. Analysis of DNA binding sites using chromatin-immunoprecipitation (ChIP) revealed that TBX15 bound directly to a key region in the Prdm16 promoter, indicating it regulates transcription of Prdm16, the master gene for adipocyte thermogenesis and browning. Compared to control mice, Tbx15 AKO mice displayed increased body weight gain and decreased whole body energy expenditure in response to high fat diets.

CONCLUSION

Taken together, these findings suggest that Tbx15 regulates adipocyte browning and might be a potential target for the treatment of obesity.

Semaphorin 3E/PlexinD1 signaling is required for cardiac ventricular compaction

Left ventricular noncompaction (LVNC) is one of the most common forms of genetic cardiomyopathy characterized by excessive trabeculation and impaired myocardial compaction during fetal development. Patients with LVNC are at higher risk of developing left/right ventricular failure or both. Although the key regulators for cardiac chamber development are well studied, the role of semaphorin (Sema)/plexin signaling in this process remains poorly understood. In this article, we demonstrate that genetic deletion of Plxnd1, a class-3 Sema receptor in endothelial cells, leads to severe cardiac chamber defects. They were characterized by excessive trabeculation and noncompaction similar to patients with LVNC. Loss of Plxnd1 results in decreased expression of extracellular matrix proteolytic genes, leading to excessive deposition of cardiac jelly. We demonstrate that Plxnd1 deficiency is associated with an increase in Notch1 expression and its downstream target genes. In addition, inhibition of the Notch signaling pathway partially rescues the excessive trabeculation and noncompaction phenotype present in Plxnd1 mutants. Furthermore, we demonstrate that Semaphorin 3E (Sema3E), one of PlexinD1's known ligands, is expressed in the developing heart and is required for myocardial compaction. Collectively, our study uncovers what we believe to be a previously undescribed role of the Sema3E/PlexinD1 signaling pathway in myocardial trabeculation and the compaction process.

Gene manipulation in liver ductal organoids by optimized recombinant adeno-associated virus vectors

Understanding the mechanism of how liver ductal cells (cholangiocytes) differentiate into hepatocytes would permit liver-regenerative medicine. Emerging liver ductal organoids provide an ex vivo system to investigate cholangiocyte-to-hepatocyte differentiation. However, as current gene manipulation methods require organoid dissociation into single cells and have only low efficiency, it is difficult to dissect specific gene functions in these organoids. Here we developed the adeno-associated virus (AAV) vector AAV-DJ as a powerful tool to transduce mouse and human liver ductal organoids. Via AAV-DJ-mediated up- or down-regulation of target genes, we successfully manipulated cholangiocyte-to-hepatocyte differentiation. We induced differentiation by overexpressing the hepatocyte-specifying regulator hepatocyte nuclear factor 4α (HNF4α) and blocked differentiation by stimulating Notch signaling or interfering with Smad signaling. Further screening for transcriptional factors critical for cholangiocyte-to-hepatocyte differentiation identified HOP homeobox (HOPX), T-box 15 (TBX15), and transcription factor CP2-like 1 (TFCP2L1) as master regulators. We conclude that this highly efficient and convenient gene manipulation system we developed could facilitate investigation into genes involved in cell lineage transitions and enable application of engineered organoids in regenerative medicine.

BDNF is a mediator of glycolytic fiber-type specification in mouse skeletal muscle

Brain-derived neurotrophic factor (BDNF) influences the differentiation, plasticity, and survival of central neurons and likewise, affects the development of the neuromuscular system. Besides its neuronal origin, BDNF is also a member of the myokine family. However, the role of skeletal muscle-derived BDNF in regulating neuromuscular physiology in vivo remains unclear. Using gain- and loss-of-function animal models, we show that muscle-specific ablation of BDNF shifts the proportion of muscle fibers from type IIB to IIX, concomitant with elevated slow muscle-type gene expression. Furthermore, BDNF deletion reduces motor end plate volume without affecting neuromuscular junction (NMJ) integrity. These morphological changes are associated with slow muscle function and a greater resistance to contraction-induced fatigue. Conversely, BDNF overexpression promotes a fast muscle-type gene program and elevates glycolytic fiber number. These findings indicate that BDNF is required for fiber-type specification and provide insights into its potential modulation as a therapeutic target in muscle diseases.

Methyltransferase-like 21e inhibits 26S proteasome activity to facilitate hypertrophy of type IIb myofibers

Skeletal muscles contain heterogeneous myofibers that are different in size and contractile speed, with type IIb myofiber being the largest and fastest. Here, we identify methyltransferase-like 21e (Mettl21e), a member of newly classified nonhistone methyltransferases, as a gene enriched in type IIb myofibers. The expression of Mettl21e was strikingly up-regulated in hypertrophic muscles and during myogenic differentiation in vitro and in vivo. Knockdown (KD) of Mettl21e led to atrophy of cultured myotubes, and targeted mutation of Mettl21e in mice reduced the size of IIb myofibers without affecting the composition of myofiber types. Mass spectrometry and methyltransferase assay revealed that Mettl21e methylated valosin-containing protein (Vcp/p97), a key component of the ubiquitin-proteasome system. KD or knockout of Mettl21e resulted in elevated 26S proteasome activity, and inhibition of proteasome activity prevented atrophy of Mettl21e KD myotubes. These results demonstrate that Mettl21e functions to maintain myofiber size through inhibiting proteasome-mediated protein degradation.-Wang, C., Zhang, B., Ratliff, A. C., Arrington, J., Chen, J., Xiong, Y., Yue, F., Nie, Y., Hu, K., Jin, W., Tao, W. A., Hrycyna, C. A., Sun, X., Kuang, S. Methyltransferase-like 21e inhibits 26S proteasome activity to facilitate hypertrophy of type IIb myofibers.

Geniposide Improves Glucose Homeostasis via Regulating FoxO1/PDK4 in Skeletal Muscle

It is well-known that imbalance state of glucose metabolism triggers many metabolic diseases and glucose uptake in skeletal muscle accounts for 90% of body weight. Geniposide is one of the major natural bioactive constituents of gardenia fruit, and the regulation of geniposide on glucose metabolism in skeletal muscle has not yet been investigated. Here, on the basis of microarray analysis, we discovered that geinposide decreased pyruvate dehydrogenase kinase 4 (PDK4) expression in skeletal muscle of mice and subsequently found that geniposide inhibited the expressions of forkhead box O1 (FoxO1), PDK4, and phosphorylated pyruvate dehydrogenase in vitro and in vivo. Moreover, geniposide promoted a switch of slow-to-fast myofiber type and glucose utilization, suggesting that geniposide improved glucose homeostasis. In addition, mechanistic studies revealed that geniposide played above roles by regulating FoxO1/PDK4, which controlled fuel selection via pyruvate dehydrogenase. Meanwhile, effects of geniposide mentioned above could be reversed by FoxO1 overexpression. Together, these results establish that geniposide confers controls on fuel usage and glucose homeostasis through FoxO1/PDK4 in skeletal muscle.

Deficiency in the secreted protein Semaphorin3d causes abnormal parathyroid development in mice

Primary hyperparathyroidism (PHPT) is a common endocrinopathy characterized by hypercalcemia and elevated levels of parathyroid hormone. The primary cause of PHPT is a benign overgrowth of parathyroid tissue causing excessive secretion of parathyroid hormone. However, the molecular etiology of PHPT is incompletely defined. Here, we demonstrate that semaphorin3d (Sema3d), a secreted glycoprotein, is expressed in the developing parathyroid gland in mice. We also observed that genetic deletion of Sema3d leads to parathyroid hyperplasia, causing PHPT. In vivo and in vitro experiments using histology, immunohistochemistry, biochemical, RT-qPCR, and immunoblotting assays revealed that Sema3d inhibits parathyroid cell proliferation by decreasing the epidermal growth factor receptor (EGFR)/Erb-B2 receptor tyrosine kinase (ERBB) signaling pathway. We further demonstrate that EGFR signaling is elevated in Sema3d ^-/- parathyroid glands and that pharmacological inhibition of EGFR signaling can partially rescue the parathyroid hyperplasia phenotype. We propose that because Sema3d is a secreted protein, it may be possible to use recombinant Sema3d or derived peptides to inhibit parathyroid cell proliferation causing hyperplasia and hyperparathyroidism. Collectively, these findings identify Sema3d as a negative regulator of parathyroid growth.

Deciphering White Adipose Tissue Heterogeneity

Adipose tissue not only stores energy, but also controls metabolism through secretion of hormones, cytokines, proteins, and microRNAs that affect the function of cells and tissues throughout the body. Adipose tissue is organized into discrete depots throughout the body, and these depots are differentially associated with insulin resistance and increased risk of metabolic disease. In addition to energy-dissipating brown and beige adipocytes, recent lineage tracing studies have demonstrated that individual adipose depots are composed of white adipocytes that are derived from distinct precursor populations, giving rise to distinct subpopulations of energy-storing white adipocytes. In this review, we discuss this developmental and functional heterogeneity of white adipocytes both between and within adipose depots. In particular, we will highlight findings from our recent manuscript in which we find and characterize three major subtypes of white adipocytes. We will discuss these data relating to the differences between subcutaneous and visceral white adipose tissue and in relationship to previous work deciphering adipocyte heterogeneity within adipose tissue depots. Finally, we will discuss the possible implications of adipocyte heterogeneity may have for the understanding of lipodystrophies.

Protein-coding variants implicate novel genes related to lipid homeostasis contributing to body-fat distribution

Body-fat distribution is a risk factor for adverse cardiovascular health consequences. We analyzed the association of body-fat distribution, assessed by waist-to-hip ratio adjusted for body mass index, with 228,985 predicted coding and splice site variants available on exome arrays in up to 344,369 individuals from five major ancestries (discovery) and 132,177 European-ancestry individuals (validation). We identified 15 common (minor allele frequency, MAF ≥5%) and nine low-frequency or rare (MAF <5%) coding novel variants. Pathway/gene set enrichment analyses identified lipid particle, adiponectin, abnormal white adipose tissue physiology and bone development and morphology as important contributors to fat distribution, while cross-trait associations highlight cardiometabolic traits. In functional follow-up analyses, specifically in Drosophila RNAi-knockdowns, we observed a significant increase in the total body triglyceride levels for two genes (DNAH10 and PLXND1). We implicate novel genes in fat distribution, stressing the importance of interrogating low-frequency and protein-coding variants.

Ca2+-Calcineurin-NFAT pathway mediates the effect of thymol on oxidative metabolism and fiber-type switch in skeletal muscle

Thymol is a major component of thyme, and it has been reported that thymol administration reduces body weight, plasma insulin and blood glucose in type-2 diabetes.

Developmental and functional heterogeneity of white adipocytes within a single fat depot

Recent studies suggest that, even within a single adipose depot, there may be distinct subpopulations of adipocytes. To investigate this cellular heterogeneity, we have developed multiple conditionally immortalized clonal preadipocyte lines from white adipose tissue of mice. Analysis of these clones reveals at least three white adipocyte subpopulations. These subpopulations have differences in metabolism and differentially respond to inflammatory cytokines, insulin, and growth hormones. These also have distinct gene expression profiles and can be tracked by differential expression of three marker genes: Wilms' tumor 1, transgelin, and myxovirus 1. Lineage tracing analysis with dual-fluorescent reporter mice indicates that these adipocyte subpopulations have differences in gene expression and metabolism that mirror those observed in the clonal cell lines. Furthermore, preadipocytes and adipocytes from these subpopulations differ in their abundance in different fat depots. Thus, white adipose tissue, even in a single depot, is comprised of distinct subpopulations of white adipocytes with different physiological phenotypes. These differences in adipocyte composition may contribute to the differences in metabolic behavior and physiology of different fat depots.