Suggestion of GLYAT gene underlying variation of bone size and body lean mass as revealed by a bivariate genome-wide association study

A systematic review and cross-database analysis of single nucleotide polymorphisms underlying hip morphology and osteoarthritis reveals shared mechanisms

OBJECTIVE

Understanding the mechanisms of hip disease, such as osteoarthritis (OA), is crucial to advance their treatment. Such hip diseases often involve specific morphological changes. Genetic variations, called single nucleotide polymorphisms (SNPs), influence various hip morphological parameters. This study investigated the biological relevance of SNPs correlated to hip morphology in genome-wide association studies (GWAS). The SNP-associated genes were compared to genes associated with OA in other joints, aiming to see if the same genes play a role in both hip development and the risk of OA in other lower limb joints.

METHODOLOGY

A systematic literature review was conducted to identify SNPs correlated with hip morphology, based on the Population, Intervention, Comparison, Outcome, and Study (PICOS) framework. Afterwards, Gene Ontology (GO) analysis was performed, using EnrichR, on the SNP-associated genes and compared with non-hip OA-associated genes, across different databases.

RESULTS

Reviewing 49 GWAS identified 436 SNPs associated with hip joint morphology, encompassing variance in bone size, structure and shape. Among the SNP-associated genes, SOX9 plays a pivotal role in size, GDF5 impacts bone structure, and BMP7 affects shape. Overall, skeletal system development, regulation of cell differentiation, and chondrocyte differentiation emerged as crucial processes influencing hip morphology. Eighteen percent of GWAS-identified genes related to hip morphology were also associated with non-hip OA.

CONCLUSION

Our findings indicate the existence of multiple shared genetic mechanisms across hip morphology and OA, highlighting the necessity for more extensive research in this area, as in contrast to the hip, the genetic background on knee or foot morphology remains largely understudied.

Pleiotropic effects on Sarcopenia subphenotypes point to potential molecular markers for the disease

Sarcopenia is a progressive age-related muscle disease characterized by low muscle strength, quantity and quality, and low physical performance. The clinical overlap between these subphenotypes (reduction in muscle strength, quantity and quality, and physical performance) was evidenced, but the genetic overlap is still poorly investigated. Herein, we investigated whether there is a genetic overlap amongst sarcopenia subphenotypes in the search for more effective molecular markers for this disease. For that, a Bioinformatics approach was used to identify and characterize pleiotropic effects at the genome, loci and gene levels using Genome-wide association study results. As a result, a high genetic correlation was identified between gait speed and muscle strength (rG=0.5358, p=3.39 × 10-8). Using a Pleiotropy-informed conditional and conjunctional false discovery rate method we identified two pleiotropic loci for muscle strength and gait speed, one of them was nearby the gene PHACTR1. Moreover, 11 pleiotropic loci and 25 genes were identified for muscle mass and muscle strength. Lastly, using a gene-based GWAS approach three candidate genes were identified in the overlap of the three Sarcopenia subphenotypes: FTO, RPS10 and CALCR. The current study provides evidence of genetic overlap and pleiotropy among sarcopenia subphenotypes and highlights novel candidate genes and molecular markers associated with the risk of sarcopenia.

Identification and Characterization of Genomic Predictors of Sarcopenia and Sarcopenic Obesity Using UK Biobank Data

The substantial decline in skeletal muscle mass, strength, and gait speed is a sign of severe sarcopenia, which may partly depend on genetic risk factors. So far, hundreds of genome-wide significant single nucleotide polymorphisms (SNPs) associated with handgrip strength, lean mass and walking pace have been identified in the UK Biobank cohort; however, their pleiotropic effects on all three phenotypes have not been investigated. By combining summary statistics of genome-wide association studies (GWAS) of handgrip strength, lean mass and walking pace, we have identified 78 independent SNPs (from 73 loci) associated with all three traits with consistent effect directions. Of the 78 SNPs, 55 polymorphisms were also associated with body fat percentage and 25 polymorphisms with type 2 diabetes (T2D), indicating that sarcopenia, obesity and T2D share many common risk alleles. Follow-up bioinformatic analysis revealed that sarcopenia risk alleles were associated with tiredness, falls in the last year, neuroticism, alcohol intake frequency, smoking, time spent watching television, higher salt, white bread, and processed meat intake; whereas protective alleles were positively associated with bone mineral density, serum testosterone, IGF1, and 25-hydroxyvitamin D levels, height, intelligence, cognitive performance, educational attainment, income, physical activity, ground coffee drinking and healthier diet (muesli, cereal, wholemeal or wholegrain bread, potassium, magnesium, cheese, oily fish, protein, water, fruit, and vegetable intake). Furthermore, the literature data suggest that single-bout resistance exercise may induce significant changes in the expression of 26 of the 73 implicated genes in m. vastus lateralis, which may partly explain beneficial effects of strength training in the prevention and treatment of sarcopenia. In conclusion, we have identified and characterized 78 SNPs associated with sarcopenia and 55 SNPs with sarcopenic obesity in European-ancestry individuals from the UK Biobank.

Multi-omics research in sarcopenia: Current progress and future prospects

Sarcopenia is a systemic disease with progressive and generalized skeletal muscle dysfunction defined by age-related low muscle mass, high content of muscle slow fibers, and low muscle function. Muscle phenotypes and sarcopenia risk are heritable; however, the genetic architecture and molecular mechanisms underlying sarcopenia remain largely unclear. In recent years, significant progress has been made in determining susceptibility loci using genome-wide association studies. In addition, recent advances in omics techniques, including genomics, epigenomics, transcriptomics, proteomics, and metabolomics, offer new opportunities to identify novel targets to help us understand the pathophysiology of sarcopenia. However, each individual technology cannot capture the entire view of the biological complexity of this disorder, while integrative multi-omics analyses may be able to reveal new insights. Here, we review the latest findings of multi-omics studies for sarcopenia and provide an in-depth summary of our current understanding of sarcopenia pathogenesis. Leveraging multi-omics data could give us a holistic understanding of sarcopenia etiology that may lead to new clinical applications. This review offers guidance and recommendations for fundamental research, innovative perspectives, and preventative and therapeutic interventions for sarcopenia.

Unveiling genetic variants for age-related sarcopenia by conducting a genome-wide association study on Korean cohorts

Sarcopenia is an age-related disorder characterised by a progressive decrease in skeletal muscle mass. As the genetic biomarkers for sarcopenia are not yet well characterised, this study aimed to investigate the genetic variations related to sarcopenia in a relatively aged cohort, using genome-wide association study (GWAS) meta-analyses of lean body mass (LBM) in 6961 subjects. Two Korean cohorts were analysed, and subgroup GWAS was conducted for appendicular skeletal muscle mass (ASM) and skeletal muscle index. The effects of significant single nucleotide polymorphisms (SNPs) on gene expression were also investigated using multiple expression quantitative trait loci datasets, differentially expressed gene analysis, and gene ontology analyses. Novel genetic biomarkers were identified for LBM (rs1187118; rs3768582) and ASM (rs6772958). Their related genes, including RPS10, NUDT3, NCF2, SMG7, and ARPC5, were differently expressed in skeletal muscle tissue, while GPD1L was not. Furthermore, the 'mRNA destabilisation' biological process was enriched for sarcopenia. Our study identified RPS10, NUDT3, and GPD1L as significant genetic biomarkers for sarcopenia. These genetic loci were related to lipid and energy metabolism, suggesting that genes involved in metabolic dysregulation may lead to the pathogenesis of age-related sarcopenia.

Influencing Factors and Molecular Pathogenesis of Sarcopenia and Osteosarcopenia in Chronic Liver Disease

The liver plays a pivotal role in nutrient/energy metabolism and storage, anabolic hormone regulation, ammonia detoxification, and cytokine production. Impaired liver function can cause malnutrition, hyperammonemia, and chronic inflammation, leading to an imbalance between muscle protein synthesis and proteolysis. Patients with chronic liver disease (CLD) have a high prevalence of sarcopenia, characterized by progressive loss of muscle mass and function, affecting health-related quality of life and prognosis. Recent reports have revealed that osteosarcopenia, defined as the concomitant occurrence of sarcopenia and osteoporosis, is also highly prevalent in patients with CLD. Since the differentiation and growth of muscles and bones are closely interrelated through mechanical and biochemical communication, sarcopenia and osteoporosis often progress concurrently and affect each other. Osteosarcopenia further exacerbates unfavorable health outcomes, such as vertebral fracture and frailty. Therefore, a comprehensive assessment of sarcopenia, osteoporosis, and osteosarcopenia, and an understanding of the pathogenic mechanisms involving the liver, bones, and muscles, are important for prevention and treatment. This review summarizes the molecular mechanisms of sarcopenia and osteosarcopenia elucidated to data in hopes of promoting advances in treating these musculoskeletal disorders in patients with CLD.

Identification of novel pleiotropic gene for bone mineral density and lean mass using the cFDR method

Bone mineral density (BMD) and whole-body lean mass (WBLM) are two important phenotypes of osteoporosis and sarcopenia. Previous studies have shown that BMD and lean mass were phenotypically and genetically correlated. To identify the novel common genetic factors shared between BMD and WBLM, we performed the conditional false discovery rate (cFDR) analysis using summary data of the genome-wide association study of femoral neck BMD (n = 53,236) and WBLM (n = 38,292) from the Genetic Factors for Osteoporosis Consortium (GEFOS). We identified eight pleiotropic Single Nucleotide Polymorphism (SNPs) (PLCL1 rs11684176 and rs2880389, JAZF1 rs198, ADAMTSL3 rs10906982, RFTN2/MARS2 rs7340470, SH3GL3 rs1896797, ST7L rs10776755, ANKRD44/SF3B1 rs11888760) significantly associated with femoral neck BMD and WBLM (ccFDR < 0.05). Bayesian fine-mapping analysis showed that rs11888760, rs198, and rs1896797 were the possible functional variants in the ANKRD44/SF3B1, JAZF1i, and SH3GL3 loci, respectively. Functional annotation suggested that rs11888760 was likely to comprise a DNA regulatory element and linked to the expression of RFTN2 and PLCL1. PLCL1 showed differential expression in laryngeal posterior cricoarytenoid muscle between rats of 6 months and 30 months of age. Our findings, together with PLCL1's potential functional relevance to bone and skeletal muscle function, suggested that rs11888760 was the possible pleiotropic functional variants appearing to coregulate both bone and muscle metabolism through regulating the expression of PLCL1. The findings enhanced our knowledge of genetic associations between BMD and lean mass and provide a rationale for subsequent functional studies of the implicated genes in the pathophysiology of diseases, such as osteoporosis and sarcopenia.

Twelve years of GWAS discoveries for osteoporosis and related traits: advances, challenges and applications

Osteoporosis is a common skeletal disease, affecting ~200 million people around the world. As a complex disease, osteoporosis is influenced by many factors, including diet (e.g. calcium and protein intake), physical activity, endocrine status, coexisting diseases and genetic factors. In this review, we first summarize the discovery from genome-wide association studies (GWASs) in the bone field in the last 12 years. To date, GWASs and meta-analyses have discovered hundreds of loci that are associated with bone mineral density (BMD), osteoporosis, and osteoporotic fractures. However, the GWAS approach has sometimes been criticized because of the small effect size of the discovered variants and the mystery of missing heritability, these two questions could be partially explained by the newly raised conceptual models, such as omnigenic model and natural selection. Finally, we introduce the clinical use of GWAS findings in the bone field, such as the identification of causal clinical risk factors, the development of drug targets and disease prediction. Despite the fruitful GWAS discoveries in the bone field, most of these GWAS participants were of European descent, and more genetic studies should be carried out in other ethnic populations to benefit disease prediction in the corresponding population.

Downregulation of GLYAT Facilitates Tumor Growth and Metastasis and Poor Clinical Outcomes Through the PI3K/AKT/Snail Pathway in Human Breast Cancer

Background

The Glycine N-acyltransferase (GLYAT) gene encodes a protein that catalyzes the transfer of acyl groups from acyl CoA to glycine, resulting in acyl glycine and coenzyme A. Aberrant GLYAT expression is associated with several malignant tumors, but its clinical importance in human breast cancer (BC), has yet to be fully addressed. This study aims to evaluate the clinical function of GLYAT in BC patients.

Methods

GLYAT expression was determined by immune blot and immunohistochemistry in three BC cell lines and primary cancer tissues. The MDA-MB 231 cell line was used for GLYAT gene knockdown experiments while the MCF7 cell line for overexpression experiments. Colony formation experiments, soft agar experiments, and transwell assays were utilized for further inspection of cell proliferation and migration capabilities. Immunofluorescence and western blot were used to detect markers of the epithelial-mesenchymal transition (EMT) and changes in the PI3K/AKT/Snail pathway. The role of GLYAT in tumor growth and metastasis was also assessed in nude mice in vivo. Also, a correlation analysis was performed between clinicopathological features and GLYAT expression in BC patients.

Results

GLYAT was decreased in human BC tissues and cell lines. Functional analysis showed that knockdown of GLYAT augmented BC cell proliferation in vitro and in vivo. However, this phenomenon was reversed when GLYAT was overexpressed in the transfected cells. Moreover, downregulation of GLYAT promoted the migratory properties of BC cells, likely through the activation of PI3K/AKT/Snail signaling, which subsequently induced the EMT. IHC analysis indicated that GLYAT was decreased in human BC tissues and lower GLYAT expression was correlated with histological grade, tumor TNM stage, Ki-67 status, and poorer survival in BC patients. Furthermore, lower GLYAT expression seemed as an independent risk factor related to poor prognosis in BC patients based on Cox regression analyses.

Conclusion

Our findings demonstrate that downregulation of GLYAT expression in human breast cancer is correlated with EMT via the PI3K/AKT/Snail pathway and is also associated with histological grade, tumor TNM stage, Ki-67 status, and poor survival in breast cancer patients.

Enhanced Identification of Novel Potential Variants for Appendicular Lean Mass by Leveraging Pleiotropy With Bone Mineral Density

Strong relationships have been found between appendicular lean mass (ALM) and bone mineral density (BMD). It may be due to a shared genetic basis, termed pleiotropy. By leveraging the pleiotropy with BMD, the aim of this study was to detect more potential genetic variants for ALM. Using the conditional false discovery rate (cFDR) methodology, a combined analysis of the summary statistics of two large independent genome wide association studies (GWAS) of ALM (n = 73,420) and BMD (n = 10,414) was conducted. Strong pleiotropic enrichment and 26 novel potential pleiotropic SNPs were found for ALM and BMD. We identified 156 SNPs for ALM (cFDR <0.05), of which 74 were replicates of previous GWASs and 82 were novel SNPs potentially-associated with ALM. Eleven genes annotated by 31 novel SNPs (13 pleiotropic and 18 ALM specific) were partially validated in a gene expression assay. Functional enrichment analysis indicated that genes corresponding to the novel potential SNPs were enriched in GO terms and/or KEGG pathways that played important roles in muscle development and/or BMD metabolism (adjP <0.05). In protein–protein interaction analysis, rich interactions were demonstrated among the proteins produced by the corresponding genes. In conclusion, the present study, as in other recent studies we have conducted, demonstrated superior efficiency and reliability of the cFDR methodology for enhanced detection of trait-associated genetic variants. Our findings shed novel insight into the genetic variability of ALM in addition to the shared genetic basis underlying ALM and BMD.

An Overview of the Molecular Mechanisms Contributing to Musculoskeletal Disorders in Chronic Liver Disease: Osteoporosis, Sarcopenia, and Osteoporotic Sarcopenia

The prevalence of osteoporosis and sarcopenia is significantly higher in patients with liver disease than in those without liver disease and osteoporosis and sarcopenia negatively influence morbidity and mortality in liver disease, yet these musculoskeletal disorders are frequently overlooked in clinical practice for patients with chronic liver disease. The objective of this review is to provide a comprehensive understanding of the molecular mechanisms of musculoskeletal disorders accompanying the pathogenesis of liver disease. The increased bone resorption through the receptor activator of nuclear factor kappa (RANK)-RANK ligand (RANKL)-osteoprotegerin (OPG) system and upregulation of inflammatory cytokines and decreased bone formation through increased bilirubin and sclerostin and lower insulin-like growth factor-1 are important mechanisms for osteoporosis in patients with liver disease. Sarcopenia is associated with insulin resistance and obesity in non-alcoholic fatty liver disease, whereas hyperammonemia, low amount of branched chain amino acids, and hypogonadism contributes to sarcopenia in liver cirrhosis. The bidirectional crosstalk between muscle and bone through myostatin, irisin, β-aminoisobutyric acid (BAIBA), osteocalcin, as well as the activation of the RANK and the Wnt/β-catenin pathways are associated with osteosarcopenia. The increased understandings for these musculoskeletal disorders would be contributes to the development of effective therapies targeting the pathophysiological mechanism involved.

Variants in NEB and RIF1 genes on chr2q23 are associated with skeletal muscle index in Koreans: genome-wide association study

Although skeletal muscle plays a crucial role in metabolism and influences aging and chronic diseases, little is known about the genetic variations with skeletal muscle, especially in the Asian population. We performed a genome-wide association study in 2,046 participants drawn from a population-based study. Appendicular skeletal muscle mass was estimated based on appendicular lean soft tissue measured with a multi-frequency bioelectrical impedance analyzer and divided by height squared to derive the skeletal muscle index (SMI). After conducting quality control and imputing the genotypes, we analyzed 6,391,983 autosomal SNPs. A genome-wide significant association was found for the intronic variant rs138684936 in the NEB and RIF1 genes (β = 0.217, p = 6.83 × 10^–9). These two genes are next to each other and are partially overlapped on chr2q23. We conducted extensive functional annotations to gain insight into the directional biological implication of significant genetic variants. A gene-based analysis identified the significant TNFSF9 gene and confirmed the suggestive association of the NEB gene. Pathway analyses showed the significant association of regulation of multicellular organism growth gene-set and the suggestive associations of pathways related to skeletal system development or skeleton morphogenesis with SMI. In conclusion, we identified a new genetic locus on chromosome 2 for SMI with genome-wide significance. These results enhance the biological understanding of skeletal muscle mass and provide specific leads for functional experiments.

Three pleiotropic loci associated with bone mineral density and lean body mass

Both bone mineral density (BMD) and lean body mass (LBM) are important physiological measures with strong genetic determination. Besides, BMD and LBM might have common genetic factors. Aiming to identify pleiotropic genomic loci underlying BMD and LBM, we performed bivariate genome-wide association study meta-analyses of femoral neck bone mineral density and LBM at arms and legs, and replicated in the large-scale UK Biobank cohort sample. Combining the results from discovery meta-analysis and replication sample, we identified three genomic loci at the genome-wide significance level (p < 5.0 × 10-8): 2p23.2 (lead SNP rs4477866, discovery p = 3.47 × 10-8, replication p = 1.03 × 10-4), 16q12.2 (rs1421085, discovery p = 2.04 × 10-9, replication p = 6.47 × 10-14) and 18q21.32 (rs11152213, discovery p = 3.47 × 10-8, replication p = 6.69 × 10-6). Our findings not only provide useful insights into lean mass and bone mass development, but also enhance our understanding of the potential genetic correlation between BMD and LBM.

Intersecting Genetics of Frailty and Cardiovascular Disease

OBJECTIVES

To determine the genetic correlates of physical frailty and sarcopenia, focusing on single nucleotide polymorphisms (SNPs) in genome-wide association studies (GWAS), and to explore the genetic overlap of frailty with cardiovascular disease (CVD) and its risk factors.

METHODS

PubMed was systematically searched for GWAS studies investigating the association between SNPs and objective measures of physical frailty or sarcopenia. SNPs were retained if they were associated with one of the phenotypes of interest by a p-value of 5.0x10-8 or less.

RESULTS

Ten studies were included, with a total of 237 SNPs in 181 genes being associated with physical frailty or sarcopenia; as measured by handgrip strength or lean (muscle) mass. These genes were cross-referenced in the GWAS Catalog, and many of them were found to be associated with CVD or metabolic syndrome.

CONCLUSIONS

Evidence from GWAS has shown that frailty is associated with common genetic polymorphisms. Many of these polymorphisms have been implicated in CVD, supporting the hypothesis of a shared pathophysiology between these entities. Future studies are eagerly anticipated to map out the mechanistic links and discover therapeutic targets and novel biomarkers for frailty.

Identification of critical genetic variants associated with metabolic phenotypes of the Japanese population

We performed a metabolome genome-wide association study for the Japanese population in the prospective cohort study of Tohoku Medical Megabank. By combining whole-genome sequencing and nontarget metabolome analyses, we identified a large number of novel associations between genetic variants and plasma metabolites. Of the identified metabolite-associated genes, approximately half have already been shown to be involved in various diseases. We identified metabolite-associated genes involved in the metabolism of xenobiotics, some of which are from intestinal microorganisms, indicating that the identified genetic variants also markedly influence the interaction between the host and symbiotic bacteria. We also identified five associations that appeared to be female-specific. A number of rare variants that influence metabolite levels were also found, and combinations of common and rare variants influenced the metabolite levels more profoundly. These results support our contention that metabolic phenotyping provides important insights into how genetic and environmental factors provoke human diseases.

Replication of FTO Gene associated with lean mass in a Meta-Analysis of Genome-Wide Association Studies

Sarcopenia is characterized by low skeletal muscle, a complex trait with high heritability. With the dramatically increasing prevalence of obesity, obesity and sarcopenia occur simultaneously, a condition known as sarcopenic obesity. Fat mass and obesity-associated (FTO) gene is a candidate gene of obesity. To identify associations between lean mass and FTO gene, we performed a genome-wide association study (GWAS) of lean mass index (LMI) in 2207 unrelated Caucasian subjects and replicated major findings in two replication samples including 6,004 unrelated Caucasian and 38,292 unrelated Caucasian. We found 29 single nucleotide polymorphisms (SNPs) in FTO significantly associated with sarcopenia (combined p-values ranging from 5.92 × 10⁻¹² to 1.69 × 10⁻⁹). Potential biological functions of SNPs were analyzed by HaploReg v4.1, RegulomeDB, GTEx, IMPC and STRING. Our results provide suggestive evidence that FTO gene is associated with lean mass.

Association of 3p27.1 Variants with Whole Body Lean Mass Identified by a Genome-wide Association Study

Whole body lean mass (WBLM) is a heritable trait predicting sarcopenia. To identify genomic locus underlying WBLM, we performed a genome-wide association study of fat-adjusted WBLM in the Framingham Heart Study (FHS, N = 6,004), and replicated in the Kansas City Osteoporosis Study (KCOS, N = 2,207). We identified a novel locus 3p27.1 that was associated with WBLM (lead SNP rs3732593 P = 7.19 × 10^-8) in the discovery FHS sample, and the lead SNP was successfully replicated in the KCOS sample (one-sided P = 0.04). Bioinformatics analysis found that this SNP and its adjacent SNPs had the function of regulating enhancer activity in skeletal muscle myoblasts cells, further confirming the regulation of WBLM by this locus. Our finding provides new insight into the genetics of WBLM and enhance our understanding of sarcopenia.

Genome-wide Associations Reveal Human-Mouse Genetic Convergence and Modifiers of Myogenesis, CPNE1 and STC2

Muscle bulk in adult healthy humans is highly variable even after height, age, and sex are accounted for. Low muscle mass, due to fewer and/or smaller constituent muscle fibers, would exacerbate the impact of muscle loss occurring in aging or disease. Genetic variability substantially influences muscle mass differences, but causative genes remain largely unknown. In a genome-wide association study (GWAS) on appendicular lean mass (ALM) in a population of 85,750 middle-aged (aged 38-49 years) individuals from the UK Biobank (UKB), we found 182 loci associated with ALM (p < 5 × 10-8). We replicated associations for 78% of these loci (p < 5 × 10-8) with ALM in a population of 181,862 elderly (aged 60-74 years) individuals from UKB. We also conducted a GWAS on hindlimb skeletal muscle mass of 1,867 mice from an advanced intercross between two inbred strains (LG/J and SM/J); this GWAS identified 23 quantitative trait loci. Thirty-eight positional candidates distributed across five loci overlapped between the two species. In vitro studies of positional candidates confirmed CPNE1 and STC2 as modifiers of myogenesis. Collectively, these findings shed light on the genetics of muscle mass variability in humans and identify targets for the development of interventions for treatment of muscle loss. The overlapping results between humans and the mouse model GWAS point to shared genetic mechanisms across species.

GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures

Bone area is one measure of bone size that is easily derived from dual-energy X-ray absorptiometry (DXA) scans. In a GWA study of DXA bone area of the hip and lumbar spine (N ≥ 28,954), we find thirteen independent association signals at twelve loci that replicate in samples of European and East Asian descent (N = 13,608 - 21,277). Eight DXA area loci associate with osteoarthritis, including rs143384 in GDF5 and a missense variant in COL11A1 (rs3753841). The strongest DXA area association is with rs11614913[T] in the microRNA MIR196A2 gene that associates with lumbar spine area (P = 2.3 × 10-42, β = -0.090) and confers risk of hip fracture (P = 1.0 × 10-8, OR = 1.11). We demonstrate that the risk allele is less efficient in repressing miR-196a-5p target genes. We also show that the DXA area measure contributes to the risk of hip fracture independent of bone density.

Genetics of Bone and Muscle Interactions in Humans

PURPOSE OF REVIEW

To summarize the evidence from recent studies on the shared genetics between bone and muscle in humans.

RECENT FINDINGS

Genome-wide association studies (GWAS) have successfully identified a multitude of loci influencing the variability of different bone or muscle parameters, with multiple loci overlapping between the traits. In addition, joint analyses of multiple correlated musculoskeletal traits (i.e., multivariate GWAS) have underscored several genes with possible pleiotropic effects on both bone and muscle including MEF2C and SREBF1. Notably, several of the proposed pleiotropic genes have been validated using human cells or animal models. It is clear that the study of pleiotropy may provide novel insights into disease pathophysiology potentially leading to the identification of new treatment strategies that simultaneously prevent or treat both osteoporosis and sarcopenia. However, the role of muscle factors (myokines) that stimulate bone metabolism, as well as osteokines that affect muscles, is in its earliest stage of understanding.

Assessing the genetic correlations between early growth parameters and bone mineral density: A polygenic risk score analysis

OBJECTIVE

The relationships between early growth parameters and bone mineral density (BMD) remain elusive now. In this study, we performed a large scale polygenic risk score (PRS) analysis to evaluate the potential impact of early growth parameters on the variations of BMD.

METHODS

We used 2286 Caucasian subjects as cohort 1 and 3404 Framingham Heart Study (FHS) subjects as cohort 2 in this study. BMD at ulna & radius, hip and spine were measured using dual energy X-ray absorptiometry. BMD values were adjusted for age, sex, height and weight as covariates. Genome-wide single-nucleotide polymorphism (SNP) genotyping of the 2286 Caucasian subjects was performed using Affymetrix Human SNP Array 6.0. The GWAS datasets of early growth parameters were driven from the Early Growth Genetics Consortium, including birth weight (BW), birth head circumference (BHC), childhood body mass index (CBMI), pubertal height growth related indexes and tanner stage. Polygenic Risk Score (PRSice) and linkage disequilibrium (LD) score regression analysis were conducted to assess the genetic correlation between early growth parameters and BMD.

RESULTS

We detected significant genetic correlations in cohort 1, such as total spine BMD vs. CBMI (p value = 1.51 × 10^-4, rg = 0.4525), right ulna and radius BMD vs. CBMI (p value = 1.51 × 10^-4, rg = 0.4399) and total body BMD vs. tanner stage (p value = 7.00 × 10^-4, rg = -0.0721). For cohort 2, significant correlations were observed for total spine BMD vs. height change standard deviation score (SDS) between 8 years and adult (denoted as PGF + PGM) (p value = 3.97 × 10^-4, rg = -0.1425), femoral neck BMD vs. the timing of peak height velocity by looking at the height change SDS between age 14 years and adult (denoted as PTF + PTM) (p value = 7.04 × 10^-4, rg = -0.2185), and total spine BMD vs. PTF + PTM (p value = 6.86 × 10^-4, rg = -0.2180).

CONCLUSION

Our study results suggest that some early growth parameters could affect the variations of BMD.

Burly1 is a mouse QTL for lean body mass that maps to a 0.8-Mb region of chromosome 2

To fine map a mouse QTL for lean body mass (Burly1), we used information from intercross, backcross, consomic, and congenic mice derived from the C57BL/6ByJ (host) and 129P3/J (donor) strains. The results from these mapping populations were concordant and showed that Burly1 is located between 151.9 and 152.7 Mb (rs33197365 to rs3700604) on mouse chromosome 2. The congenic region harboring Burly1 contains 26 protein-coding genes, 11 noncoding RNA elements (e.g., lncRNA), and 4 pseudogenes, with 1949 predicted functional variants. Of the protein-coding genes, 7 have missense variants, including genes that may contribute to lean body weight, such as Angpt41, Slc52c3, and Rem1. Lean body mass was increased by the B6-derived variant relative to the 129-derived allele. Burly1 influenced lean body weight at all ages but not food intake or locomotor activity. However, congenic mice with the B6 allele produced more heat per kilogram of lean body weight than did controls, pointing to a genotype effect on lean mass metabolism. These results show the value of integrating information from several mapping populations to refine the map location of body composition QTLs and to identify a short list of candidate genes.

Bivariate genome-wide association meta-analysis of pediatric musculoskeletal traits reveals pleiotropic effects at the SREBF1/TOM1L2 locus

Bone mineral density is known to be a heritable, polygenic trait whereas genetic variants contributing to lean mass variation remain largely unknown. We estimated the shared SNP heritability and performed a bivariate GWAS meta-analysis of total-body lean mass (TB-LM) and total-body less head bone mineral density (TBLH-BMD) regions in 10,414 children. The estimated SNP heritability is 43% (95% CI: 34–52%) for TBLH-BMD, and 39% (95% CI: 30–48%) for TB-LM, with a shared genetic component of 43% (95% CI: 29–56%). We identify variants with pleiotropic effects in eight loci, including seven established bone mineral density loci: WNT4, GALNT3, MEPE, CPED1/WNT16, TNFSF11, RIN3, and PPP6R3/LRP5. Variants in the TOM1L2/SREBF1 locus exert opposing effects TB-LM and TBLH-BMD, and have a stronger association with the former trait. We show that SREBF1 is expressed in murine and human osteoblasts, as well as in human muscle tissue. This is the first bivariate GWAS meta-analysis to demonstrate genetic factors with pleiotropic effects on bone mineral density and lean mass.

Bone mineral density and lean skeletal mass are heritable traits. Here, Medina-Gomez and colleagues perform bivariate GWAS analyses of total body lean mass and bone mass density in children, and show genetic loci with pleiotropic effects on both traits.

Large meta-analysis of genome-wide association studies identifies five loci for lean body mass

Lean body mass, consisting mostly of skeletal muscle, is important for healthy aging. We performed a genome-wide association study for whole body (20 cohorts of European ancestry with n = 38,292) and appendicular (arms and legs) lean body mass (n = 28,330) measured using dual energy X-ray absorptiometry or bioelectrical impedance analysis, adjusted for sex, age, height, and fat mass. Twenty-one single-nucleotide polymorphisms were significantly associated with lean body mass either genome wide (p < 5 × 10-8) or suggestively genome wide (p < 2.3 × 10-6). Replication in 63,475 (47,227 of European ancestry) individuals from 33 cohorts for whole body lean body mass and in 45,090 (42,360 of European ancestry) subjects from 25 cohorts for appendicular lean body mass was successful for five single-nucleotide polymorphisms in/near HSD17B11, VCAN, ADAMTSL3, IRS1, and FTO for total lean body mass and for three single-nucleotide polymorphisms in/near VCAN, ADAMTSL3, and IRS1 for appendicular lean body mass. Our findings provide new insight into the genetics of lean body mass.Lean body mass is a highly heritable trait and is associated with various health conditions. Here, Kiel and colleagues perform a meta-analysis of genome-wide association studies for whole body lean body mass and find five novel genetic loci to be significantly associated.

Evaluation of the relationship and genetic overlap between Kashin-Beck disease and body mass index

OBJECTIVES

Body mass index (BMI) is one of the major factors affecting the development of osteoarthritis (OA) but there is currently no information available regarding the relationship between BMI and Kashin-Beck disease (KBD). Our aim in this study was to investigate the relationship and genetic overlap between BMI and KBD.

METHOD

A total of 2050 Han Chinese subjects participated in this study. Using a cohort of 333 grade I KBD patients, logistic regression analysis was conducted to evaluate the correlation between BMI and KBD. Another independent sample of 1717 subjects was genotyped for a genome-wide association study (GWAS) using Affymetrix Human SNP 6.0 Arrays. Single nucleotide polymorphism (SNP) effect concordance analysis (SECA) was applied to the GWAS summaries of KBD and BMI for pleiotropy analysis. Genome-wide bivariate association analysis (GWBAA) of KBD and BMI was carried out to identify the genes with pleiotropic effects on KBD and BMI. The relevance of identified genes with KBD was validated by gene expression profiling and immunohistochemistry.

RESULTS

BMI correlated positively with knee movement disorder in KBD (coefficient β = 0.068, p = 0.045). SECA identified a significant pleiotropic effect (empirical p = 0.021) between KBD and BMI. In the GWBAA, the rs1893577 of the ADAMTS1 gene achieved the most significant association signal (p = 7.38 × 10^-9). ADAMTS1 was also up-regulated in KBD vs. normal (ratio = 2.64 ± 2.80) and KBD vs. OA (ratio = 2.31 ± 2.01). The rate of ADAMTS1-positive chondrocytes in KBD was significantly higher than that in OA (p < 0.05) and healthy controls (p < 0.05).

CONCLUSIONS

Our results suggest that ADAMTS1 is a novel susceptibility gene for KBD.

A bivariate genome-wide association study identifies ADAM12 as a novel susceptibility gene for Kashin-Beck disease

Kashin-Beck disease (KBD) is a chronic osteoarthropathy, which manifests as joint deformities and growth retardation. Only a few genetic studies of growth retardation associated with the KBD have been carried out by now. In this study, we conducted a two-stage bivariate genome-wide association study (BGWAS) of the KBD using joint deformities and body height as study phenotypes, totally involving 2,417 study subjects. Articular cartilage specimens from 8 subjects were collected for immunohistochemistry. In the BGWAS, ADAM12 gene achieved the most significant association (rs1278300 p-value = 9.25 × 10⁻⁹) with the KBD. Replication study observed significant association signal at rs1278300 (p-value = 0.007) and rs1710287 (p-value = 0.002) of ADAM12 after Bonferroni correction. Immunohistochemistry revealed significantly decreased expression level of ADAM12 protein in the KBD articular cartilage (average positive chondrocyte rate = 47.59 ± 7.79%) compared to healthy articular cartilage (average positive chondrocyte rate = 64.73 ± 5.05%). Our results suggest that ADAM12 gene is a novel susceptibility gene underlying both joint destruction and growth retardation of the KBD.

Contribution of Heritability and Epigenetic Factors to Skeletal Muscle Mass Variation in United Kingdom Twins

CONTEXT

Skeletal muscle mass (SMM) is one of the major components of human body composition, with deviations from normal values often leading to sarcopenia.

OBJECTIVE

Our major aim was to conduct a genome-wide DNA methylation study in an attempt to identify potential genomic regions associated with SMM.

DESIGN

This was a mixed cross-sectional and longitudinal study.

SETTING

Community-based study.

PARTICIPANTS

A total of 1550 middle-aged United Kingdom twins (monozygotic [MZ] and dizygotic [DZ]), 297 of which were repeatedly measured participated in the study.

MAIN OUTCOME MEASURE

Appendicular lean mass assessed using dual-energy X-ray absorptiometry technology, and methylated DNA immunoprecipitation sequencing DNA methylation profiling genome-wide were obtained from each individual.

RESULTS

Heritability estimate of SMM, with simultaneous adjustment for covariates obtained using variance decomposition analysis, was h(2) = 0.809 ± 0.050. After quality control and analysis of longitudinal stability, the DNA methylation data comprised of 723 029 genomic sites, with positive correlations between repeated measurements (Rrepeated = 0.114-0.905). Correlations between MZ and DZ twins were 0.51 and 0.38 at a genome-wide average, respectively, and clearly increased with Rrepeated. Testing for DNA methylation association with SMM in 50 discordant MZ twins revealed 36 081 nominally significant results, of which the top-ranked 134 signals (P < .01 and Rrepeated > 0.40) were subjected to replication in the sample of 1196 individuals. Seven SMM methylation association signals replicated at a false discovery rate less than 0.1, and these were located in or near genes DNAH12, CAND1, CYP4F29P, and ZFP64, which have previously been highlighted in muscle-related studies. Adjusting for age, smoking, and blood cell heterogeneity did not alter significance of these associations.

CONCLUSION

This epigenome-wide study, testing longitudinally stable methylation sites, discovered and replicated a number of associations between DNA methylation at CpG loci and SMM. Four replicated signals were related to genes with potential muscle functions, suggesting that the methylome of whole blood may be informative of SMM variation.

The genetics of bone mass and susceptibility to bone diseases

Osteoporosis is characterized by low bone mass and an increased risk of fracture. Genetic factors, environmental factors and gene-environment interactions all contribute to a person's lifetime risk of developing an osteoporotic fracture. This Review summarizes key advances in understanding of the genetics of bone traits and their role in osteoporosis. Candidate-gene approaches dominated this field 20 years ago, but clinical and preclinical genetic studies published in the past 5 years generally utilize more-sophisticated and better-powered genome-wide association studies (GWAS). High-throughput DNA sequencing, large genomic databases and improved methods of data analysis have greatly accelerated the gene-discovery process. Linkage analyses of single-gene traits that segregate in families with extreme phenotypes have led to the elucidation of critical pathways controlling bone mass. For example, components of the Wnt-β-catenin signalling pathway have been validated (in both GWAS and functional studies) as contributing to various bone phenotypes. These notable advances in gene discovery suggest that the next decade will witness cataloguing of the hundreds of genes that influence bone mass and osteoporosis, which in turn will provide a roadmap for the development of new drugs that target diseases of low bone mass, including osteoporosis.

Bone and muscle: Interactions beyond mechanical

The musculoskeletal system is significantly more complex than portrayed by traditional reductionist approaches that have focused on and studied the components of this system separately. While bone and skeletal muscle are the two largest tissues within this system, this system also includes tendons, ligaments, cartilage, joints and other connective tissues along with vascular and nervous tissues. Because the main function of this system is locomotion, the mechanical interaction among the major players of this system is essential for the many shapes and forms observed in vertebrates and even in invertebrates. Thus, it is logical that the mechanical coupling theories of musculoskeletal development exert a dominant influence on our understanding of the biology of the musculoskeletal system, because these relationships are relatively easy to observe, measure, and perturb. Certainly much less recognized is the molecular and biochemical interaction among the individual players of the musculoskeletal system. In this brief review article, we first introduce some of the key reasons why the mechanical coupling theory has dominated our view of bone-muscle interactions followed by summarizing evidence for the secretory nature of bones and muscles. Finally, a number of highly physiological questions that cannot be answered by the mechanical theories alone will be raised along with different lines of evidence that support both a genetic and a biochemical communication between bones and muscles. It is hoped that these discussions will stimulate new insights into this fertile and promising new way of defining the relationships between these closely related tissues. Understanding the cellular and molecular mechanisms responsible for biochemical communication between bone and muscle is important not only from a basic research perspective but also as a means to identify potential new therapies for bone and muscle diseases, especially for when they co-exist. This article is part of a Special Issue entitled "Muscle Bone Interactions".

Genomics and metabolomics of muscular mass in a community-based sample of UK females

The contribution of specific molecular-genetic factors to muscle mass variation and sarcopenia remains largely unknown. To identify endogenous molecules and specific genetic factors associated with appendicular lean mass (APLM) in the general population, cross-sectional data from the TwinsUK Adult Twin Registry were used. Non-targeted mass spec-based metabolomic profiling was performed on plasma of 3953 females (mostly dizygotic and monozygotic twins). APLM was measured using dual-energy X-ray absorptiometry (DXA) and genotyping was genome-wide (GWAS). Specific metabolites were used as intermediate phenotypes in the identification of single-nucleotide polymorphisms associated with APLM using GWAS. In all, 162 metabolites were found significantly correlated with APLM, and explained 17.4% of its variation. However, the top three of them (unidentified substance X12063, urate, and mannose) explained 11.1% (P ≤ 9.25 × 10(-26)) so each was subjected to GWAS. Each metabolite showed highly significant (P ≤ 9.28 × 10(-46)) associations with genetic variants in the corresponding genomic regions. Mendelian randomization using these SNPs found no evidence for a direct causal effect of these metabolites on APLM. However, using a new software platform for bivariate analysis we showed that shared genetic factors contribute significantly (P ≤ 4.31 × 10(-43)) to variance in both the metabolites and APLM--independent of the effect of the associated SNPs. There are several metabolites, having a clear pattern of genetic inheritance, which are highly significantly associated with APLM and may provide a cheap and readily accessible biomarker of muscle mass. However, the mechanism by which the genetic factor influences muscle mass remains to be discovered.

Recent genetic discoveries in osteoporosis, sarcopenia and obesity

Osteoporosis is a skeletal disorder characterized by low bone mineral density (BMD) and an increased susceptibility to fractures. Evidence from genetic studies indicates that BMD, a complex quantitative trait with a normal distribution, is genetically controlled. Genome-wide association studies (GWAS) as well as studies using candidate gene approaches have identified single-nucleotide polymorphisms (SNPs) that are associated with BMD, osteoporosis and osteoporotic fractures. These SNPs have been mapped close to or within genes including those encoding WNT/β-catenin signaling proteins. Understanding the genetics of osteoporosis will help to identify novel candidates for diagnostic and therapeutic targets. Genetic factors are also important for the development of sarcopenia, which is characterized by a loss of lean body mass, and obesity, which is characterized by high fat mass. Hence, in this review, we discuss the genetic factors, identified by genetic studies, which regulate the body components related to osteoporosis, sarcopenia, and obesity.

Large-scale analysis reveals a functional single-nucleotide polymorphism in the 5'-flanking region of PRDM16 gene associated with lean body mass

Genetic factors are important for the development of sarcopenia, a geriatric disorder characterized by low lean body mass. The aim of this study was to search for novel genes that regulate lean body mass in humans. We performed a large-scale search for 250K single-nucleotide polymorphisms (SNPs) associated with bone mineral density (BMD) using SNP arrays in 1081 Japanese postmenopausal women. We focused on an SNP (rs12409277) located in the 5′-flanking region of the PRDM16 (PRD1-BF-1-RIZ1 homologous domain containing protein 16) gene that showed a significant P value in our screening. We demonstrated that PRDM16 gene polymorphisms were significantly associated with total body BMD in 1081 postmenopausal Japanese women. The rs12409277 SNP affected the transcriptional activity of PRDM16. The subjects with one or two minor allele(s) had a higher lean body mass than the subjects with two major alleles. Genetic analyses uncovered the importance of the PRDM16 gene in the regulation of lean body mass.

METTL21C is a potential pleiotropic gene for osteoporosis and sarcopenia acting through the modulation of the NF-κB signaling pathway

Sarcopenia and osteoporosis are important public health problems that occur concurrently. A bivariate genome-wide association study (GWAS) identified METTL21c as a suggestive pleiotropic gene for both bone and muscle. The METTL21 family of proteins methylates chaperones involved in the etiology of both myopathy and inclusion body myositis with Paget's disease. To validate these GWAS results, Mettl21c mRNA expression was reduced with siRNA in a mouse myogenic C2C12 cell line and the mouse osteocyte-like cell line MLO-Y4. At day 3, as C2C12 myoblasts start to differentiate into myotubes, a significant reduction in the number of myocytes aligning/organizing for fusion was observed in the siRNA-treated cells. At day 5, both fewer and smaller myotubes were observed in the siRNA-treated cells as confirmed by histomorphometric analyses and immunostaining with myosin heavy chain (MHC) antibody, which only stains myocytes/myotubes but not myoblasts. Intracellular calcium (Ca(2+)) measurements of the siRNA-treated myotubes showed a decrease in maximal amplitude peak response to caffeine, suggesting that less Ca(2+) is available for release due to the partial silencing of Mettl21c, correlating with impaired myogenesis. In siRNA-treated MLO-Y4 cells, 48 hours after treatment with dexamethasone there was a significant increase in cell death, suggesting a role of Mettl21c in osteocyte survival. To investigate the molecular signaling machinery induced by the partial silencing of Mettl21c, we used a real-time PCR gene array to monitor the activity of 10 signaling pathways. We discovered that Mettl21c knockdown modulated only the NF-κB signaling pathway (ie, Birc3, Ccl5, and Tnf). These results suggest that Mettl21c might exert its bone-muscle pleiotropic function via the regulation of the NF-κB signaling pathway, which is critical for bone and muscle homeostasis. These studies also provide rationale for cellular and molecular validation of GWAS, and warrant additional in vitro and in vivo studies to advance our understanding of role of METTL21C in musculoskeletal biology.

Endocrine crosstalk between muscle and bone

The musculoskeletal system is a complex organ comprised of the skeletal bones, skeletal muscles, tendons, ligaments, cartilage, joints, and other connective tissue that physically and mechanically interact to provide animals and humans with the essential ability of locomotion. This mechanical interaction is undoubtedly essential for much of the diverse shape and forms observed in vertebrates and even in invertebrates with rudimentary musculoskeletal systems such as fish. It makes sense from a historical point of view that the mechanical theories of musculoskeletal development have had tremendous influence of our understanding of biology, because these relationships are clear and palpable. Less visible to the naked eye or even to the microscope is the biochemical interaction among the individual players of the musculoskeletal system. It was only in recent years that we have begun to appreciate that beyond this mechanical coupling of muscle and bones, these 2 tissues function at a higher level through crosstalk signaling mechanisms that are important for the function of the concomitant tissue. Our brief review attempts to present some of the key concepts of these new concepts and is outline to present muscles and bones as secretory/endocrine organs, the evidence for mutual genetic and tissue interactions, pathophysiological examples of crosstalk, and the exciting new directions for this promising field of research aimed at understanding the biochemical/molecular coupling of these 2 intimately associated tissues.

A new perspective on the importance of glycine conjugation in the metabolism of aromatic acids

A number of endogenous and xenobiotic organic acids are conjugated to glycine, in animals ranging from mosquitoes to humans. Glycine conjugation has generally been assumed to be a detoxification mechanism, increasing the water solubility of organic acids in order to facilitate urinary excretion. However, the recently proposed glycine deportation hypothesis states that the role of the amino acid conjugations, including glycine conjugation, is to regulate systemic levels of amino acids that are also utilized as neurotransmitters in the central nervous systems of animals. This hypothesis is based on the observation that, compared to glucuronidation, glycine conjugation does not significantly increase the water solubility of aromatic acids. In this review it will be argued that the major role of glycine conjugation is to dispose of the end products of phenylpropionate metabolism. Furthermore, glucuronidation, which occurs in the endoplasmic reticulum, would not be ideal for the detoxification of free benzoate, which has been shown to accumulate in the mitochondrial matrix. Glycine conjugation, however, prevents accumulation of benzoic acid in the mitochondrial matrix by forming hippurate, a less lipophilic conjugate that can be more readily transported out of the mitochondria. Finally, it will be explained that the glycine conjugation of benzoate, a commonly used preservative, exacerbates the dietary deficiency of glycine in humans. Because the resulting shortage of glycine can negatively influence brain neurochemistry and the synthesis of collagen, nucleic acids, porphyrins, and other important metabolites, the risks of using benzoate as a preservative should not be underestimated.

The role of muscle in bone repair: the cells, signals, and tissue responses to injury

Bone repair is a complicated process that includes many types of cells, signaling molecules, and growth factors. Fracture healing involves a temporally and spatially regulated biologic process that involves recruitment of stem cells to the injury site, tissue specific differentiation, angiogenesis, and remodeling. In light of its proximity to bone and abundant vascularity, muscle is an important potential source of cells and signals for bone healing. More complete understanding of the role of muscle in bone formation and repair will provide new therapeutic approaches to enhance fracture healing. Recent studies establish that muscle-derived stem cells are able to differentiate into cartilage and bone and can directly participate in fracture healing. The role of muscle-derived stem cells is particularly important in fractures associated with more severe injury to the periosteum. Sarcopenia is a serious consequence of aging, and studies show a strong association between bone mass and lean muscle mass. Muscle anabolic agents may improve function and reduce the incidence of fracture with aging.

Glycine conjugation: importance in metabolism, the role of glycine N-acyltransferase, and factors that influence interindividual variation

INTRODUCTION

Glycine conjugation of mitochondrial acyl-CoAs, catalyzed by glycine N-acyltransferase (GLYAT, E.C. 2.3.1.13), is an important metabolic pathway responsible for maintaining adequate levels of free coenzyme A (CoASH). However, because of the small number of pharmaceutical drugs that are conjugated to glycine, the pathway has not yet been characterized in detail. Here, we review the causes and possible consequences of interindividual variation in the glycine conjugation pathway.

AREAS COVERED

The authors review the importance of CoASH in metabolism, formation and toxicity of xenobiotic acyl-CoAs, and mechanisms for restoring levels of CoASH. They focus on GLYAT, glycine conjugation, how genetic variation in the GLYAT gene could influence glycine conjugation, and the emerging roles of glycine metabolism in cancer and musculoskeletal development.

EXPERT OPINION

The substrate selectivity of GLYAT and its variants needs to be further characterized, as organic acids can be toxic if the corresponding acyl-CoA is not a substrate for glycine conjugation. GLYAT activity affects mitochondrial ATP production, glycine availability, CoASH availability, and the toxicity of various organic acids. Therefore, variation in the glycine conjugation pathway could influence liver cancer, musculoskeletal development, and mitochondrial energy metabolism.