A coding and non-coding transcriptomic perspective on the genomics of human metabolic disease

Identification of different myofiber types in pigs muscles and construction of regulatory networks

BACKGROUND

Skeletal muscle is composed of muscle fibers with different physiological characteristics, which plays an important role in regulating skeletal muscle metabolism, movement and body homeostasis. The type of skeletal muscle fiber directly affects meat quality. However, the transcriptome and gene interactions between different types of muscle fibers are not well understood.

RESULTS

In this paper, we selected 180-days-old Large White pigs and found that longissimus dorsi (LD) muscle was dominated by fast-fermenting myofibrils and soleus (SOL) muscle was dominated by slow-oxidizing myofibrils by frozen sections and related mRNA and protein assays. Here, we selected LD muscle and SOL muscle for transcriptomic sequencing, and identified 312 differentially expressed mRNA (DEmRs), 30 differentially expressed miRNA (DEmiRs), 183 differentially expressed lncRNA (DElRs), and 3417 differentially expressed circRNA (DEcRs). The ceRNA network included ssc-miR-378, ssc-miR-378b-3p, ssc-miR-24-3p, XR_308817, XR_308823, SMIM8, MAVS and FOS as multiple core nodes that play important roles in muscle development. Moreover, we found that different members of the miR-10 family expressed differently in oxidized and glycolytic muscle fibers, among which miR-10a-5p was highly expressed in glycolytic muscle fibers (LD) and could target MYBPH gene mRNA. Therefore, we speculate that miR-10a-5p may be involved in the transformation of muscle fiber types by targeting the MYHBP gene. In addition, PPI analysis of differentially expressed mRNA genes showed that ACTC1, ACTG2 and ACTN2 gene had the highest node degree, suggesting that this gene may play a key role in the regulatory network of muscle fiber type determination.

CONCLUSIONS

We can conclude that these genes play a key role in regulating muscle fiber type transformation. Our study provides transcriptomic profiles and ceRNA interaction networks for different muscle fiber types in pigs, providing reference for the transformation of pig muscle fiber types and the improvement of meat quality.

ECHDC3 Variant Regulates the Right Hippocampal Microstructural Integrity and Verbal Memory in Type 2 Diabetes Mellitus

ECHDC3 is a risk gene for white matter (WM) hyperintensity and is associated with insulin resistance. This study aimed to investigate whether ECHDC3 variants selectively regulate brain WM microstructures and episodic memory in patients with type 2 diabetes mellitus (T2DM). We enrolled 106 patients with T2DM and 111 healthy controls. A voxel-wise general linear model was employed to explore the interaction effect between ECHDC3 rs11257311 polymorphism and T2DM diagnosis on fractional anisotropy (FA). A linear modulated mediation analysis was conducted to examine the potential of FA value to mediate the influence of T2DM on episodic memory in an ECHDC3-dependent manner. We observed a noteworthy interaction between genotype and diagnosis on FA in the right inferior temporal WM, right anterior limb of the internal capsule, right frontal WM, and the right hippocampus. Modulated mediation analysis revealed a significant ECHDC3 modulation on the T2DM → right hippocampal FA → short-term memory pathway, with only rs11257311 G risk homozygote demonstrating significant mediation effect. Together, our findings provide evidence of ECHDC3 modulating the effect of T2DM on right hippocampal microstructural impairment and short-term memory decline, which might be a neuro-mechanism for T2DM related episodic memory impairment.

Mitochondrial Alpha-Keto Acid Dehydrogenase Complexes: Recent Developments on Structure and Function in Health and Disease

The present work delves into the enigmatic world of mitochondrial alpha-keto acid dehydrogenase complexes discussing their metabolic significance, enzymatic operation, moonlighting activities, and pathological relevance with links to underlying structural features. This ubiquitous family of related but diverse multienzyme complexes is involved in carbohydrate metabolism (pyruvate dehydrogenase complex), the citric acid cycle (α-ketoglutarate dehydrogenase complex), and amino acid catabolism (branched-chain α-keto acid dehydrogenase complex, α-ketoadipate dehydrogenase complex); the complexes all function at strategic points and also participate in regulation in these metabolic pathways. These systems are among the largest multienzyme complexes with at times more than 100 protein chains and weights ranging up to ~10 million Daltons. Our chapter offers a wealth of up-to-date information on these multienzyme complexes for a comprehensive understanding of their significance in health and disease.

MCCC2 is a novel mediator between mitochondria and telomere and functions as an oncogene in colorectal cancer

BACKGROUND

The mitochondrial gene MCCC2, a subunit of the heterodimer of 3-methylcrotonyl-CoA carboxylase, plays a pivotal role in catabolism of leucine and isovaleric acid. The molecular mechanisms and prognostic value still need to be explored in the context of specific cancers, including colorectal cancer (CRC).

METHODS

In vitro and in vivo cell-based assays were performed to explore the role of MCCC2 in CRC cell proliferation, invasion, and migration. Mitochondrial morphology, membrane potential, intracellular reactive oxygen species (ROS), telomerase activity, and telomere length were examined and analyzed accordingly. Protein complex formation was detected by co-immunoprecipitation (CO-IP). Mitochondrial morphology was observed by transmission electron microscopy (TEM). The Cancer Genome Atlas (TCGA) CRC cohort analysis, qRT-PCR, and immunohistochemistry (IHC) were used to examine the MCCC2 expression level. The association between MCCC2 expression and various clinical characteristics was analyzed by chi-square tests. CRC patients' overall survival (OS) was analyzed by Kaplan-Meier analysis.

RESULTS

Ectopic overexpression of MCCC2 promoted cell proliferation, invasion, and migration, while MCCC2 knockdown (KD) or knockout (KO) inhibited cell proliferation, invasion, and migration. MCCC2 KD or KO resulted in reduced mitochondria numbers, but did not affect the gross ATP production in the cells. Mitochondrial fusion markers MFN1, MFN2, and OPA1 were all upregulated in MCCC2 KD or KO cells, which is in line with a phenomenon of more prominent mitochondrial fusion. Interestingly, telomere lengths of MCCC2 KD or KO cells were reduced more than control cells. Furthermore, we found that MCCC2 could specifically form a complex with telomere binding protein TRF2, and MCCC2 KD or KO did not affect the expression or activity of telomerase reverse transcriptase (TERT). Finally, MCCC2 expression was heightened in CRC, and patients with higher MCCC2 expression had favorable prognosis.

CONCLUSIONS

Together, we identified MCCC2 as a novel mediator between mitochondria and telomeres, and provided an additional biomarker for CRC stratification.

Comprehensive interrogation of human skeletal muscle reveals a dissociation between insulin resistance and mitochondrial capacity

Insulin resistance and blunted mitochondrial capacity in skeletal muscle are often synonymous, however, this association remains controversial. The aim of this study was to perform an in-depth multifactorial comparison of skeletal muscle mitochondrial capacity between individuals who were lean and active (Active, n = 9), individuals with obesity (Obese, n = 9), and individuals with obesity, insulin resistance, and type 2 diabetes (T2D, n = 22). Mitochondrial capacity was assessed by ex vivo mitochondrial respiration with fatty-acid and glycolytic-supported protocols adjusted for mitochondrial content (mtDNA and citrate synthase activity). Supercomplex assembly was measured by Blue Native (BN)-PAGE and immunoblot. Tricarboxylic (TCA) cycle intermediates were assessed with targeted metabolomics. Exploratory transcriptomics and DNA methylation analyses were performed to uncover molecular differences affecting mitochondrial function among the three groups. We reveal no discernable differences in skeletal muscle mitochondrial content, mitochondrial capacity, supercomplex assembly, TCA cycle intermediates, and mitochondrial molecular profiles between obese individuals with and without T2D that had comparable levels of confounding factors (body mass index, age, and aerobic capacity). We highlight that lean, active individuals have greater mitochondrial content, mitochondrial capacity, supercomplex assembly, and TCA cycle intermediates. These phenotypical changes are reflected at the level of DNA methylation and gene transcription. The collective observation of comparable muscle mitochondrial capacity in individuals with obesity and T2D (vs. individuals without T2D) underscores a dissociation from skeletal muscle insulin resistance. Clinical trial number: NCT01911104.NEW & NOTEWORTHY Whether impaired mitochondrial capacity contributes to skeletal muscle insulin resistance is debated. Our multifactorial analysis shows no differences in skeletal muscle mitochondrial content, mitochondrial capacity, and mitochondrial molecular profiles between obese individuals with and without T2D that had comparable levels of confounding factors (BMI, age, aerobic capacity). We highlight that lean, active individuals have enhanced skeletal muscle mitochondrial capacity that is also reflected at the level of DNA methylation and gene transcription.

Dysregulation of early gene response to a mixed meal in skeletal muscle in obesity and type 2 diabetes

Obesity- and type 2 diabetes mellitus-induced changes in the expression of protein-coding genes in human skeletal muscle were extensively examined at baseline (after an overnight fast). We aimed to compare the early transcriptomic response to a typical single meal in skeletal muscle of metabolically healthy subjects and obese individuals without and with type 2 diabetes. Transcriptomic response (RNA-seq) to a mixed meal (nutritional drink, ∼25 kJ/kg of body mass) was examined in the vastus lateralis muscle (1 h after a meal) in 7 healthy subjects and 14 obese individuals without or with type 2 diabetes. In all obese individuals, the transcriptome response to a meal was dysregulated (suppressed and altered) and associated with different biological processes compared with healthy control. To search for potential transcription factors regulating transcriptomic response to a meal, the enrichment of transcription factor-binding sites in individual promoters of the human skeletal muscle was examined. In obese individuals, the transcriptomic response is associated with a different set of transcription factors than that in healthy subjects. In conclusion, metabolic disorders are associated with a defect in the regulation of mixed meal/insulin-mediated gene expression-insulin resistance in terms of gene expression. Importantly, this dysregulation occurs in obese individuals without type 2 diabetes, i.e., at the first stage of the development of metabolic disorders.NEW & NOTEWORTHY In skeletal muscle of metabolically healthy subjects, a typical single meal normalized to body mass induces activation of various transcription factors, expression of numerous receptor tyrosine kinases associated with the insulin signaling cascade, and transcription regulators. In skeletal muscle of obese individuals without and with type 2 diabetes, this signaling network is poorly regulated at the transcriptional level, indicating dysregulation of the early gene response to a mixed meal.

Leveraging genetic diversity to identify small molecules that reverse mouse skeletal muscle insulin resistance

Systems genetics has begun to tackle the complexity of insulin resistance by capitalising on computational advances to study high-diversity populations. 'Diversity Outbred in Australia (DOz)' is a population of genetically unique mice with profound metabolic heterogeneity. We leveraged this variance to explore skeletal muscle's contribution to whole-body insulin action through metabolic phenotyping and skeletal muscle proteomics of 215 DOz mice. Linear modelling identified 553 proteins that associated with whole-body insulin sensitivity (Matsuda Index) including regulators of endocytosis and muscle proteostasis. To enrich for causality, we refined this network by focusing on negatively associated, genetically regulated proteins, resulting in a 76-protein fingerprint of insulin resistance. We sought to perturb this network and restore insulin action with small molecules by integrating the Broad Institute Connectivity Map platform and in vitro assays of insulin action using the Prestwick chemical library. These complementary approaches identified the antibiotic thiostrepton as an insulin resistance reversal agent. Subsequent validation in ex vivo insulin-resistant mouse muscle and palmitate-induced insulin-resistant myotubes demonstrated potent insulin action restoration, potentially via upregulation of glycolysis. This work demonstrates the value of a drug-centric framework to validate systems-level analysis by identifying potential therapeutics for insulin resistance.

The Therapeutic Potential of Vitamins B1, B3 and B6 in Charcot-Marie-Tooth Disease with the Compromised Status of Vitamin-Dependent Processes

Understanding the molecular mechanisms of neurological disorders is necessary for the development of personalized medicine. When the diagnosis considers not only the disease symptoms, but also their molecular basis, treatments tailored to individual patients may be suggested. Vitamin-responsive neurological disorders are induced by deficiencies in vitamin-dependent processes. These deficiencies may occur due to genetic impairments of proteins whose functions are involved with the vitamins. This review considers the enzymes encoded by the DHTKD1, PDK3 and PDXK genes, whose mutations are observed in patients with Charcot-Marie-Tooth (CMT) disease. The enzymes bind or produce the coenzyme forms of vitamins B1 (thiamine diphosphate, ThDP) and B6 (pyridoxal-5'-phosphate, PLP). Alleviation of such disorders through administration of the lacking vitamin or its derivative calls for a better introduction of mechanistic knowledge to medical diagnostics and therapies. Recent data on lower levels of the vitamin B3 derivative, NAD+, in the blood of patients with CMT disease vs. control subjects are also considered in view of the NAD-dependent mechanisms of pathological axonal degeneration, suggesting the therapeutic potential of vitamin B3 in these patients. Thus, improved diagnostics of the underlying causes of CMT disease may allow patients with vitamin-responsive disease forms to benefit from the administration of the vitamins B1, B3, B6, their natural derivatives, or their pharmacological forms.

Transcriptomics for Clinical and Experimental Biology Research: Hang on a Seq

Sequencing the human genome empowers translational medicine, facilitating transcriptome-wide molecular diagnosis, pathway biology, and drug repositioning. Initially, microarrays are used to study the bulk transcriptome; but now short-read RNA sequencing (RNA-seq) predominates. Positioned as a superior technology, that makes the discovery of novel transcripts routine, most RNA-seq analyses are in fact modeled on the known transcriptome. Limitations of the RNA-seq methodology have emerged, while the design of, and the analysis strategies applied to, arrays have matured. An equitable comparison between these technologies is provided, highlighting advantages that modern arrays hold over RNA-seq. Array protocols more accurately quantify constitutively expressed protein coding genes across tissue replicates, and are more reliable for studying lower expressed genes. Arrays reveal long noncoding RNAs (lncRNA) are neither sparsely nor lower expressed than protein coding genes. Heterogeneous coverage of constitutively expressed genes observed with RNA-seq, undermines the validity and reproducibility of pathway analyses. The factors driving these observations, many of which are relevant to long-read or single-cell sequencing are discussed. As proposed herein, a reappreciation of bulk transcriptomic methods is required, including wider use of the modern high-density array data-to urgently revise existing anatomical RNA reference atlases and assist with more accurate study of lncRNAs.

GeneFriends: gene co-expression databases and tools for humans and model organisms

Gene co-expression analysis has emerged as a powerful method to provide insights into gene function and regulation. The rapid growth of publicly available RNA-sequencing (RNA-seq) data has created opportunities for researchers to employ this abundant data to help decipher the complexity and biology of genomes. Co-expression networks have proven effective for inferring the relationship between the genes, for gene prioritization and for assigning function to poorly annotated genes based on their co-expressed partners. To facilitate such analyses we created previously an online co-expression tool for humans and mice entitled GeneFriends. To continue providing a valuable tool to the scientific community, we have now updated the GeneFriends database and website. Here, we present the new version of GeneFriends, which includes gene and transcript co-expression networks based on RNA-seq data from 46 475 human and 34 322 mouse samples. The new database also encompasses tissue-specific gene co-expression networks for 20 human and 21 mouse tissues, dataset-specific gene co-expression maps based on TCGA and GTEx projects and gene co-expression networks for additional seven model organisms (fruit fly, zebrafish, worm, rat, yeast, cow and chicken). GeneFriends is freely available at http://www.genefriends.org/.

Signaling and Gene Expression in Skeletal Muscles in Type 2 Diabetes: Current Results and OMICS Perspectives

Skeletal muscles mainly contribute to the emergence of insulin resistance, impaired glucose tolerance and the development of type 2 diabetes. Molecular mechanisms that regulate glucose uptake are diverse, including the insulin-dependent as most important, and others as also significant. They involve a wide range of proteins that control intracellular traffic and exposure of glucose transporters on the cell surface to create an extensive regulatory network. Here, we highlight advantages of the omics approaches to explore the insulin-regulated proteins and genes in human skeletal muscle with varying degrees of metabolic disorders. We discuss methodological aspects of the assessment of metabolic dysregulation and molecular responses of human skeletal muscle to insulin. The known molecular mechanisms of glucose uptake regulation and the first results of phosphoproteomic and transcriptomic studies are reviewed, which unveiled a large-scale array of insulin targets in muscle cells. They demonstrate that a clear depiction of changes that occur during metabolic dysfunction requires systemic and combined analysis at different levels of regulation, including signaling pathways, transcription factors, and gene expression. Such analysis seems promising to explore yet undescribed regulatory mechanisms of glucose uptake by skeletal muscle and identify the key regulators as potential therapeutic targets.

Recent Advances on the Role of ATGL in Cancer

The hypoxic state of the tumor microenvironment leads to reprogramming lipid metabolism in tumor cells. Adipose triglyceride lipase, also known as patatin-like phospholipase= domain-containing protein 2 and Adipose triglyceride lipase (ATGL), as an essential lipid metabolism-regulating enzyme in cells, is regulated accordingly under hypoxia induction. However, studies revealed that ATGL exhibits both tumor-promoting and tumor-suppressing effects, which depend on the cancer cell type and the site of tumorigenesis. For example, elevated ATGL expression in breast cancer is accompanied by enhanced fatty acid oxidation (FAO), enhancing cancer cells’ metastatic ability. In prostate cancer, on the other hand, tumor activity tends to be negatively correlated with ATGL expression. This review outlined the regulation of ATGL-mediated lipid metabolism pathways in tumor cells, emphasizing the Hypoxia-inducible factors 1 (HIF-1)/Hypoxia-inducible lipid droplet-associated (HIG-2)/ATGL axis, peroxisome proliferator-activated receptor (PPAR)/G0/G1 switch gene 2 (G0S2)/ATGL axis, and fat-specific protein 27 (FSP-27)/Early growth response protein 1 (EGR-1)/ATGL axis. In the light of recent research on different cancer types, the role of ATGL on tumorigenesis, tumor proliferation, and tumor metastasis was systemically reviewed.

Altered mitochondrial microenvironment at the spotlight of musculoskeletal aging and Alzheimer's disease

Emerging evidence has linked Alzheimer's disease (AD) onset with musculoskeletal aging via a muscle-brain crosstalk mediated by dysregulation of the mitochondrial microenvironment. This study investigated gene expression profiles from skeletal muscle tissues of older healthy adults to identify potential gene biomarkers whose dysregulated expression and protein interactome were involved in AD. Screening of the literature resulted in 12 relevant microarray datasets (GSE25941, GSE28392, GSE28422, GSE47881, GSE47969, GSE59880) in musculoskeletal aging and (GSE4757, GSE5281, GSE16759, GSE28146, GSE48350, GSE84422) in AD. Retrieved differentially expressed genes (DEGs) were used to construct two unique protein-protein interaction networks and clustering gene modules were identified. Overlapping module DEGs in the musculoskeletal aging and AD networks were ranked based on 11 topological algorithms and the five highest-ranked ones were considered as hub genes. The analysis revealed that the dysregulated expression of the mitochondrial microenvironment genes, NDUFAB1, UQCRC1, UQCRFS1, NDUFS3, and MRPL15, overlapped between both musculoskeletal aging and AD networks. Thus, these genes may have a potential role as markers of AD occurrence in musculoskeletal aging. Human studies are warranted to evaluate the functional role and prognostic value of these genes in aging populations with sarcopenia and AD.

Delayed Impact of 2-Oxoadipate Dehydrogenase Inhibition on the Rat Brain Metabolism Is Linked to Protein Glutarylation

Background

The DHTKD1-encoded 2-oxoadipate dehydrogenase (OADH) oxidizes 2-oxoadipate—a common intermediate of the lysine and tryptophan catabolism. The mostly low and cell-specific flux through these pathways, and similar activities of OADH and ubiquitously expressed 2-oxoglutarate dehydrogenase (OGDH), agree with often asymptomatic phenotypes of heterozygous mutations in the DHTKD1 gene. Nevertheless, OADH/DHTKD1 are linked to impaired insulin sensitivity, cardiovascular disease risks, and Charcot-Marie-Tooth neuropathy. We hypothesize that systemic significance of OADH relies on its generation of glutaryl residues for protein glutarylation. Using pharmacological inhibition of OADH and the animal model of spinal cord injury (SCI), we explore this hypothesis.

Methods

The weight-drop model of SCI, a single intranasal administration of an OADH-directed inhibitor trimethyl adipoyl phosphonate (TMAP), and quantification of the associated metabolic changes in the rat brain employ established methods.

Results

The TMAP-induced metabolic changes in the brain of the control, laminectomized (LE) and SCI rats are long-term and (patho)physiology-dependent. Increased glutarylation of the brain proteins, proportional to OADH expression in the control and LE rats, represents a long-term consequence of the OADH inhibition. The proportionality suggests autoglutarylation of OADH, supported by our mass-spectrometric identification of glutarylated K155 and K818 in recombinant human OADH. In SCI rats, TMAP increases glutarylation of the brain proteins more than OADH expression, inducing a strong perturbation in the brain glutathione metabolism. The redox metabolism is not perturbed by TMAP in LE animals, where the inhibition of OADH increases expression of deglutarylase sirtuin 5. The results reveal the glutarylation-imposed control of the brain glutathione metabolism. Glutarylation of the ODP2 subunit of pyruvate dehydrogenase complex at K451 is detected in the rat brain, linking the OADH function to the brain glucose oxidation essential for the redox state. Short-term inhibition of OADH by TMAP administration manifests in increased levels of tryptophan and decreased levels of sirtuins 5 and 3 in the brain.

Conclusion

Pharmacological inhibition of OADH affects acylation system of the brain, causing long-term, (patho)physiology-dependent changes in the expression of OADH and sirtuin 5, protein glutarylation and glutathione metabolism. The identified glutarylation of ODP2 subunit of pyruvate dehydrogenase complex provides a molecular mechanism of the OADH association with diabetes.

Administration of Phosphonate Inhibitors of Dehydrogenases of 2-Oxoglutarate and 2-Oxoadipate to Rats Elicits Target-Specific Metabolic and Physiological Responses

In vitro and in cell cultures, succinyl phosphonate (SP) and adipoyl phosphonate (AP) selectively target dehydrogenases of 2-oxoglutarate (OGDH, encoded by OGDH/OGDHL) and 2-oxoadipate (OADH, encoded by DHTKD1), respectively. To assess the selectivity in animals, the effects of SP, AP, and their membrane-penetrating triethyl esters (TESP and TEAP) on the rat brain metabolism and animal physiology are compared. Opposite effects of the OGDH and OADH inhibitors on activities of OGDH, malate dehydrogenase, glutamine synthetase, and levels of glutamate, lysine, citrulline, and carnosine are shown to result in distinct physiological responses. ECG is changed by AP/TEAP, whereas anxiety is increased by SP/TESP. The potential role of the ester moiety in the uncharged precursors of the 2-oxo acid dehydrogenase inhibitors is estimated. TMAP is shown to be less efficient than TEAP, in agreement with lower lipophilicity of TMAP vs. TEAP. Non-monotonous metabolic and physiological impacts of increasing OADH inhibition are revealed. Compared to the non-treated animals, strong inhibition of OADH decreases levels of tryptophan and beta-aminoisobutyrate and activities of malate dehydrogenase and pyruvate dehydrogenase, increasing the R–R interval of ECG. Thus, both metabolic and physiological actions of the OADH-directed inhibitors AP/TEAP are different from those of the OGDH-directed inhibitors SP/TESP, with the ethyl ester being more efficient than methyl ester.

A human-based multi-gene signature enables quantitative drug repurposing for metabolic disease

Insulin resistance (IR) contributes to the pathophysiology of diabetes, dementia, viral infection, and cardiovascular disease. Drug repurposing (DR) may identify treatments for IR; however, barriers include uncertainty whether in vitro transcriptomic assays yield quantitative pharmacological data, or how to optimise assay design to best reflect in vivo human disease. We developed a clinical-based human tissue IR signature by combining lifestyle-mediated treatment responses (>500 human adipose and muscle biopsies) with biomarkers of disease status (fasting IR from >1200 biopsies). The assay identified a chemically diverse set of >130 positively acting compounds, highly enriched in true positives, that targeted 73 proteins regulating IR pathways. Our multi-gene RNA assay score reflected the quantitative pharmacological properties of a set of epidermal growth factor receptor-related tyrosine kinase inhibitors, providing insight into drug target specificity; an observation supported by deep learning-based genome-wide predicted pharmacology. Several drugs identified are suitable for evaluation in patients, particularly those with either acute or severe chronic IR.

Aberrant mitochondrial homeostasis at the crossroad of musculoskeletal ageing and non-small cell lung cancer

Cancer cachexia is accompanied by muscle atrophy, sharing multiple common catabolic pathways with sarcopenia, including mitochondrial dysfunction. This study investigated gene expression from skeletal muscle tissues of older healthy adults, who are at risk of age-related sarcopenia, to identify potential gene biomarkers whose dysregulated expression and protein interference were involved in non-small cell lung cancer (NSCLC). Screening of the literature resulted in 14 microarray datasets (GSE25941, GSE28392, GSE28422, GSE47881, GSE47969, GSE59880 in musculoskeletal ageing; GSE118370, GSE33532, GSE19804, GSE18842, GSE27262, GSE19188, GSE31210, GSE40791 in NSCLC). Differentially expressed genes (DEGs) were used to construct protein-protein interaction networks and retrieve clustering gene modules. Overlapping module DEGs were ranked based on 11 topological algorithms and were correlated with prognosis, tissue expression, and tumour purity in NSCLC. The analysis revealed that the dysregulated expression of the mammalian mitochondrial ribosomal proteins, Mitochondrial Ribosomal Protein S26 (MRPS26), Mitochondrial Ribosomal Protein S17 (MRPS17), Mitochondrial Ribosomal Protein L18 (MRPL18) and Mitochondrial Ribosomal Protein L51 (MRPL51) were linked to reduced survival and tumour purity in NSCLC while tissue expression of the same genes followed an opposite direction in healthy older adults. These results support a potential link between the mitochondrial ribosomal microenvironment in ageing muscle and NSCLC. Further studies comparing changes in sarcopenia and NSCLC associated cachexia are warranted.

Human and mouse muscle transcriptomic analyses identify insulin receptor mRNA downregulation in hyperinsulinemia-associated insulin resistance

Hyperinsulinemia is commonly viewed as a compensatory response to insulin resistance, yet studies have demonstrated that chronically elevated insulin may also drive insulin resistance. The molecular mechanisms underpinning this potentially cyclic process remain poorly defined, especially on a transcriptome-wide level. Transcriptomic meta-analysis in >450 human samples demonstrated that fasting insulin reliably and negatively correlated with INSR mRNA in skeletal muscle. To establish causality and study the direct effects of prolonged exposure to excess insulin in muscle cells, we incubated C2C12 myotubes with elevated insulin for 16 h, followed by 6 h of serum starvation, and established that acute AKT and ERK signaling were attenuated in this model of in vitro hyperinsulinemia. Global RNA-sequencing of cells both before and after nutrient withdrawal highlighted genes in the insulin receptor (INSR) signaling, FOXO signaling, and glucose metabolism pathways indicative of 'hyperinsulinemia' and 'starvation' programs. Consistently, we observed that hyperinsulinemia led to a substantial reduction in Insr gene expression, and subsequently a reduced surface INSR and total INSR protein, both in vitro and in vivo. Bioinformatic modeling combined with RNAi identified SIN3A as a negative regulator of Insr mRNA (and JUND, MAX, and MXI as positive regulators of Irs2 mRNA). Together, our analysis identifies mechanisms which may explain the cyclic processes underlying hyperinsulinemia-induced insulin resistance in muscle, a process directly relevant to the etiology and disease progression of type 2 diabetes.

Dysregulated paired related homeobox 1 impacts on hepatocellular carcinoma phenotypes

BACKGROUND

Hepatocellular carcinoma (HCC) is a major cause of cancer-related death. Paired related homeobox 1 (PRRX1) is a transcription factor that regulates cell growth and differentiation, but its importance in HCC is unclear.

METHODS

We examined the expression pattern of PRRX1 in nine microarray datasets of human HCC tumour samples (n > 1100) and analyzed its function in HCC cell lines. In addition, we performed gene set enrichment, Kaplan-Meier overall survival analysis, metabolomics and functional assays.

RESULTS

PRRX1 is frequently upregulated in human HCC. Pathway enrichment analysis predicted a direct correlation between PRRX1 and focal adhesion and epithelial-mesenchymal transition. High expression of PRRX1 and low ZEB1 or high ZEB2 significantly predicted better overall survival in HCC patients. In contrast, metabolic processes correlated inversely and transcriptional analyses revealed that glycolysis, TCA cycle and amino acid metabolism were affected. These findings were confirmed by metabolomics analysis. At the phenotypic level, PRRX1 knockdown accelerated proliferation and clonogenicity in HCC cell lines.

CONCLUSIONS

Our results suggest that PRRX1 controls metabolism, has a tumour suppressive role, and may function in cooperation with ZEB1/2. These findings have functional relevance in HCC, including in understanding transcriptional control of distinct cancer hallmarks.

Resistance Exercise, Aging, Disuse, and Muscle Protein Metabolism

Skeletal muscle is the organ of locomotion, its optimal function is critical for athletic performance, and is also important for health due to its contribution to resting metabolic rate and as a site for glucose uptake and storage. Numerous endogenous and exogenous factors influence muscle mass. Much of what is currently known regarding muscle protein turnover is owed to the development and use of stable isotope tracers. Skeletal muscle mass is determined by the meal- and contraction-induced alterations of muscle protein synthesis and muscle protein breakdown. Increased loading as resistance training is the most potent nonpharmacological strategy by which skeletal muscle mass can be increased. Conversely, aging (sarcopenia) and muscle disuse lead to the development of anabolic resistance and contribute to the loss of skeletal muscle mass. Nascent omics-based technologies have significantly improved our understanding surrounding the regulation of skeletal muscle mass at the gene, transcript, and protein levels. Despite significant advances surrounding the mechanistic intricacies that underpin changes in skeletal muscle mass, these processes are complex, and more work is certainly needed. In this article, we provide an overview of the importance of skeletal muscle, describe the influence that resistance training, aging, and disuse exert on muscle protein turnover and the molecular regulatory processes that contribute to changes in muscle protein abundance. © 2021 American Physiological Society. Compr Physiol 11:2249-2278, 2021.

Genome-Wide Association Study on Reproduction-Related Body-Shape Traits of Chinese Holstein Cows

Reproduction is an important production activity for dairy cows, and their reproductive performance can directly affect the level of farmers' income. To better understand the genomic regions and biological pathways of reproduction-related traits of dairy cows, in the present study, three body shape traits-Loin Strength (LS), Rump Angle (RA), and Pin Width (PW)-were selected as indicators of the reproductive ability of cows, and we conducted genome-wide association analyses on them. The heritability of these three traits was medium, ranging from 0.20 to 0.38. A total of 11 significant single-nucleotide polymorphisms (SNPs) were detected associated with these three traits. Bioinformatics analysis was performed on genes close to the significant SNPs (within 200 Kb) of LS, RA, and PW, and we found that these genes were totally enriched in 20 gene ontology terms and six KEGG signaling pathways. Finally, the five genes CDH12, TARP, PCDH9, DTHD1, and ARAP2 were selected as candidate genes that might affect LS. The six genes LOC781835, FSTL4, ATG4C, SH3BP4, DMP1, and DSPP were selected as candidate genes that might affect RA. The five genes USP6NL, CNTN3, LOC101907665, UPF2, and ECHDC3 were selected as candidate genes that might affect the PW of Chinese Holstein cows. Our results could provide useful biological information for the improvement of body shape traits and contribute to the genomic selection of Chinese Holstein cows.

Molecular Transducers of Human Skeletal Muscle Remodeling under Different Loading States

Loading of skeletal muscle changes the tissue phenotype reflecting altered metabolic and functional demands. In humans, heterogeneous adaptation to loading complicates the identification of the underpinning molecular regulators. A within-person differential loading and analysis strategy reduces heterogeneity for changes in muscle mass by ∼40% and uses a genome-wide transcriptome method that models each mRNA from coding exons and 3' and 5' untranslated regions (UTRs). Our strategy detects ∼3-4 times more regulated genes than similarly sized studies, including substantial UTR-selective regulation undetected by other methods. We discover a core of 141 genes correlated to muscle growth, which we validate from newly analyzed independent samples (n = 100). Further validating these identified genes via RNAi in primary muscle cells, we demonstrate that members of the core genes were regulators of protein synthesis. Using proteome-constrained networks and pathway analysis reveals notable relationships with the molecular characteristics of human muscle aging and insulin sensitivity, as well as potential drug therapies.

The transcriptome profile of human trisomy 21 blood cells

BACKGROUND

Trisomy 21 (T21) is a genetic alteration characterised by the presence of an extra full or partial human chromosome 21 (Hsa21) leading to Down syndrome (DS), the most common form of intellectual disability (ID). It is broadly agreed that the presence of extra genetic material in T21 gives origin to an altered expression of genes located on Hsa21 leading to DS phenotype. The aim of this study was to analyse T21 and normal control blood cell gene expression profiles obtained by total RNA sequencing (RNA-Seq).

RESULTS

The results were elaborated by the TRAM (Transcriptome Mapper) software which generated a differential transcriptome map between human T21 and normal control blood cells providing the gene expression ratios for 17,867 loci. The obtained gene expression profiles were validated through real-time reverse transcription polymerase chain reaction (RT-PCR) assay and compared with previously published data. A post-analysis through transcriptome mapping allowed the identification of the segmental (regional) variation of the expression level across the whole genome (segment-based analysis of expression). Interestingly, the most over-expressed genes encode for interferon-induced proteins, two of them (MX1 and MX2 genes) mapping on Hsa21 (21q22.3). The altered expression of genes involved in mitochondrial translation and energy production also emerged, followed by the altered expression of genes encoding for the folate cycle enzyme, GART, and the folate transporter, SLC19A1.

CONCLUSIONS

The alteration of these pathways might be linked and involved in the manifestation of ID in DS.

A spontaneous genetically induced epiallele at a retrotransposon shapes host genome function

Intracisternal A-particles (IAPs) are endogenous retroviruses (ERVs) responsible for most insertional mutations in the mouse. Full-length IAPs harbour genes flanked by long terminal repeats (LTRs). Here, we identify a solo LTR IAP variant (Iap5-1solo) recently formed in the inbred C57BL/6J mouse strain. In contrast to the C57BL/6J full-length IAP at this locus (Iap5-1full), Iap5-1solo lacks DNA methylation and H3K9 trimethylation. The distinct DNA methylation levels between the two alleles are established during preimplantation development, likely due to loss of KRAB zinc finger protein binding at the Iap5-1solo variant. Iap5-1solo methylation increases and becomes more variable in a hybrid genetic background yet is unresponsive to maternal dietary methyl supplementation. Differential epigenetic modification of the two variants is associated with metabolic differences and tissue-specific changes in adjacent gene expression. Our characterisation of Iap5-1 as a genetically induced epiallele with functional consequences establishes a new model to study transposable element repression and host-element co-evolution.

Isoforms of the DHTKD1-Encoded 2-Oxoadipate Dehydrogenase, Identified in Animal Tissues, Are not Observed upon the Human DHTKD1 Expression in Bacterial or Yeast Systems

Unlike the OGDH-encoded 2-oxoglutarate dehydrogenase (OGDH), which is an essential enzyme present in all animal tissues, expression of the DHTKD1-encoded isoenzyme, 2-oxoadipate dehydrogenase (OADH), depends on a number of factors, and mutant DHTKD1 phenotypes are rarely manifested. Physiological significance of OADH is also obscured by the fact that both isoenzymes transform 2-oxoglutarate and 2-oxoadipate. By analogy with other members of the 2-oxo acid dehydrogenases family, OADH is assumed to be a component of the multienzyme complex that catalyzes oxidative decarboxylation of 2-oxoadipate. This study aims at molecular characterization of OADH from animal tissues. Phylogenetic analysis of 2-oxo acid dehydrogenases reveals OADH only in animals and Dictyostelium discoideum slime mold, within a common branch with bacterial OGDH. Examination of partially purified animal OADH by immunoblotting and mass spectrometry identifies two OADH isoforms with molecular weights of about 130 and 70 kDa. These isoforms are not observed upon the expression of human DHTKD1 protein in either bacterial or yeast system, where the synthesized OADH is of expected molecular weight (about 100 kDa). Thus, the OADH isoforms present in animal tissues, may result from the animal-specific regulation of the DHTKD1 expression and/or posttranslational modifications of the encoded protein. Mapping of the peptides identified in the OADH preparations, onto the protein structure suggests that the 70-kDa isoform is truncated at the N-terminus, but retains the active site. Since the N-terminal domain of OGDH is required for the formation of the multienzyme complex, it is possible that the 70-kDa isoform catalyzes non-oxidative transformation of dicarboxylic 2-oxo acids that does not require the multienzyme structure. In this case, the ratio of the OADH isoforms in animal tissues may correspond to the ratio between the oxidative and non-oxidative decarboxylation of 2-oxoadipate.

Selective Inhibition of 2-Oxoglutarate and 2-Oxoadipate Dehydrogenases by the Phosphonate Analogs of Their 2-Oxo Acid Substrates

Phosphonate analogs of pyruvate and 2-oxoglutarate are established specific inhibitors of cognate 2-oxo acid dehydrogenases. The present work develops application of this class of compounds to specific in vivo inhibition of 2-oxoglutarate dehydrogenase (OGDH) and its isoenzyme, 2-oxoadipate dehydrogenase (OADH). The isoenzymes-enriched preparations from the rat tissues with different expression of OADH and OGDH are used to characterize their interaction with 2-oxoglutarate (OG), 2-oxoadipate (OA) and the phosphonate analogs. Despite a 100-fold difference in the isoenzymes ratio in the heart and liver, similar Michaelis saturations by OG are inherent in the enzyme preparations from these tissues (KmOG = 0.45 ± 0.06 and 0.27 ± 0.026 mM, respectively), indicating no significant contribution of OADH to the OGDH reaction, or similar affinities of the isoenzymes to OG. However, the preparations differ in the catalysis of OADH reaction. The heart preparation, where OADH/OGDH ratio is ≈ 0.01, possesses low-affinity sites to OA (KmOA = 0.55 ± 0.07 mM). The liver preparation, where OADH/OGDH ratio is ≈ 1.6, demonstrates a biphasic saturation with OA: the low-affinity sites (Km,2OA = 0.45 ± 0.12 mM) are similar to those of the heart preparation; the high-affinity sites (Km,1OA = 0.008 ± 0.001 mM), revealed in the liver preparation only, are attributed to OADH. Phosphonate analogs of C5-C7 dicarboxylic 2-oxo acids inhibit OGDH and OADH competitively to 2-oxo substrates in all sites. The high-affinity sites for OA are affected the least by the C5 analog (succinyl phosphonate) and the most by the C7 one (adipoyl phosphonate). The opposite reactivity is inherent in both the low-affinity OA-binding sites and OG-binding sites. The C6 analog (glutaryl phosphonate) does not exhibit a significant preference to either OADH or OGDH. Structural analysis of the phosphonates binding to OADH and OGDH reveals the substitution of a tyrosine residue in OGDH for a serine residue in OADH among structural determinants of the preferential binding of the bulkier ligands to OADH. The consistent kinetic and structural results expose adipoyl phosphonate as a valuable pharmacological tool for specific in vivo inhibition of the DHTKD1-encoded OADH, a new member of mammalian family of 2-oxo acid dehydrogenases, up-regulated in some cancers and associated with diabetes and obesity.

Skeletal muscle proteins important for work capacity are altered with type 2 diabetes - Effect of 10-20-30 training

The study examined whether men with type 2 diabetes exhibit lower expression of muscle proteins important for exercise capacity, and whether exercise training promotes adaptations in these proteins. In a cross-sectional and longitudinal study, conducted at the University of Copenhagen. Twelve men with type 2 diabetes (T2D) were compared to eleven nondiabetes counterparts (ND) matched for age and body composition (body fat percentage). T2D underwent 10 weeks of high-intensity interval exercise training (10-20-30 training). T2D had lower expression of SOD1 (-62%; p < 0.001) and ETC complex V (-34%; p = 0.003), along with higher expression of ETC complex IV (+66%; p = 0.007), MFN2 (+62%; p = 0.001), and DRP1 (+30%; p = 0.028) compared to ND. T2D had higher (p < 0.001) expression of Na⁺ /K⁺ α1 (+98%), α2 (+114%), and NHE1 (+144%) than ND. In T2D, training increased exercise capacity (+9%; p < 0.001) as well as expression of SOD2 (+44%; p = 0.029), ETC complex II (+25%; p = 0.035), III (+52%; p = 0.041), IV (+23%; p = 0.005), and V (+21%; p = 0.035), CS activity (+32%; p = 0.006) as well as Na⁺ /K⁺ α1 (+24%; p = 0.034), Kir6.2 (+36%; p = 0.029), and MCT1 (+20%; p = 0.007). Men with type 2 diabetes exhibited altered expression of a multitude of skeletal muscle proteins important for exercise capacity. Ten weeks of 10-20-30 training upregulated expression of muscle proteins regulating antioxidant defense, mitochondrial function, and ion handling while enhancing exercise capacity in men with type 2 diabetes.

SGCG rs679482 Associates With Weight Loss Success in Response to an Intensively Supervised Outpatient Program

Weight loss in response to energy restriction is highly variable, and identification of genetic contributors can provide insights into underlying biology. Leveraging 1000 Genomes imputed genotypes, we carried out genome-wide association study (GWAS) analysis in 551 unrelated obese subjects of European ancestry who participated in an intensively supervised weight loss program with replication of promising signals in an independent sample of 1,331 obese subjects who completed the program at a later date. By single nucleotide polymorphism-based and sib-pair analysis, we show that that weight loss is a heritable trait, with estimated heritability (h ² = 0.49) within the range reported for obesity. We find rs679482, intronic to SGCG (sarcoglycan γ), highly expressed in skeletal muscle, to concordantly associate with weight loss in discovery and replication samples reaching GWAS significance in the combined meta-analysis (β = -0.35, P = 1.7 × 10^-12). Located in a region of open chromatin, rs679482 is predicted to bind DMRT2, and allele-specific transcription factor binding analysis indicates preferential binding of DMRT2 to rs679482-A. Concordantly, rs679482-A impairs native repressor activity and increases basal and DMRT2-mediated enhancer activity. These findings confirm that weight loss is a heritable trait and provide evidence by which a novel variant in SGCG, rs679482, leads to impaired diet response.

(ADP-ribosyl)hydrolases: structure, function, and biology

ADP-ribosylation is an intricate and versatile posttranslational modification involved in the regulation of a vast variety of cellular processes in all kingdoms of life. Its complexity derives from the varied range of different chemical linkages, including to several amino acid side chains as well as nucleic acids termini and bases, it can adopt. In this review, we provide an overview of the different families of (ADP-ribosyl)hydrolases. We discuss their molecular functions, physiological roles, and influence on human health and disease. Together, the accumulated data support the increasingly compelling view that (ADP-ribosyl)hydrolases are a vital element within ADP-ribosyl signaling pathways and they hold the potential for novel therapeutic approaches as well as a deeper understanding of ADP-ribosylation as a whole.

Precision Medicine in Lifestyle Medicine: The Way of the Future?

Precision medicine has captured the imagination of the medical community with visions of therapies precisely targeted to the specific individual's genetic, biological, social, and environmental profile. However, in practice it has become synonymous with genomic medicine. As such its successes have been limited, with poor predictive or clinical value for the majority of people. It adds little to lifestyle medicine, other than in establishing why a healthy lifestyle is effective in combatting chronic disease. The challenge of lifestyle medicine remains getting people to actually adopt, sustain, and naturalize a healthy lifestyle, and this will require an approach that treats the patient as a person with individual needs and providing them with suitable types of support. The future of lifestyle medicine is holistic and person-centered rather than technological.

Recent advances in understanding resistance exercise training-induced skeletal muscle hypertrophy in humans

Skeletal muscle plays a pivotal role in the maintenance of physical and metabolic health and, critically, mobility. Accordingly, strategies focused on increasing the quality and quantity of skeletal muscle are relevant, and resistance exercise is foundational to the process of functional hypertrophy. Much of our current understanding of skeletal muscle hypertrophy can be attributed to the development and utilization of stable isotopically labeled tracers. We know that resistance exercise and sufficient protein intake act synergistically and provide the most effective stimuli to enhance skeletal muscle mass; however, the molecular intricacies that underpin the tremendous response variability to resistance exercise-induced hypertrophy are complex. The purpose of this review is to discuss recent studies with the aim of shedding light on key regulatory mechanisms that dictate hypertrophic gains in skeletal muscle mass. We also aim to provide a brief up-to-date summary of the recent advances in our understanding of skeletal muscle hypertrophy in response to resistance training in humans.

Synthetic analogues of 2-oxo acids discriminate metabolic contribution of the 2-oxoglutarate and 2-oxoadipate dehydrogenases in mammalian cells and tissues

The biological significance of the DHTKD1-encoded 2-oxoadipate dehydrogenase (OADH) remains obscure due to its catalytic redundancy with the ubiquitous OGDH-encoded 2-oxoglutarate dehydrogenase (OGDH). In this work, metabolic contributions of OADH and OGDH are discriminated by exposure of cells/tissues with different DHTKD1 expression to the synthesized phosphonate analogues of homologous 2-oxodicarboxylates. The saccharopine pathway intermediates and phosphorylated sugars are abundant when cellular expressions of DHTKD1 and OGDH are comparable, while nicotinate and non-phosphorylated sugars are when DHTKD1 expression is order(s) of magnitude lower than that of OGDH. Using succinyl, glutaryl and adipoyl phosphonates on the enzyme preparations from tissues with varied DHTKD1 expression reveals the contributions of OADH and OGDH to oxidation of 2-oxoadipate and 2-oxoglutarate in vitro. In the phosphonates-treated cells with the high and low DHTKD1 expression, adipate or glutarate, correspondingly, are the most affected metabolites. The marker of fatty acid β-oxidation, adipate, is mostly decreased by the shorter, OGDH-preferring, phosphonate, in agreement with the known OGDH dependence of β-oxidation. The longest, OADH-preferring, phosphonate mostly affects the glutarate level. Coupled decreases in sugars and nicotinate upon the OADH inhibition link the perturbation in glucose homeostasis, known in OADH mutants, to the nicotinate-dependent NAD metabolism.

Hibernation induces widespread transcriptional remodeling in metabolic tissues of the grizzly bear

Revealing the mechanisms underlying the reversible physiology of hibernation could have applications to both human and animal health as hibernation is often associated with disease-like states. The present study uses RNA-sequencing to reveal the tissue and seasonal transcriptional changes occurring in grizzly bears (Ursus arctos horribilis). Comparing hibernation to other seasons, bear adipose has a greater number of differentially expressed genes than liver and skeletal muscle. During hyperphagia, adipose has more than 900 differentially expressed genes compared to active season. Hibernation is characterized by reduced expression of genes associated with insulin signaling, muscle protein degradation, and urea production, and increased expression within muscle protein anabolic pathways. Across all three tissues we find a subset of shared differentially expressed genes, some of which are uncharacterized, that together may reflect a common regulatory mechanism. The identified gene families could be useful for developing novel therapeutics to treat human and animal diseases.

Longevity-related molecular pathways are subject to midlife "switch" in humans

Emerging evidence indicates that molecular aging may follow nonlinear or discontinuous trajectories. Whether this occurs in human neuromuscular tissue, particularly for the noncoding transcriptome, and independent of metabolic and aerobic capacities, is unknown. Applying our novel RNA method to quantify tissue coding and long noncoding RNA (lncRNA), we identified ~800 transcripts tracking with age up to ~60 years in human muscle and brain. In silico analysis demonstrated that this temporary linear "signature" was regulated by drugs, which reduce mortality or extend life span in model organisms, including 24 inhibitors of the IGF-1/PI3K/mTOR pathway that mimicked, and 5 activators that opposed, the signature. We profiled Rapamycin in nondividing primary human myotubes (n = 32 HTA 2.0 arrays) and determined the transcript signature for reactive oxygen species in neurons, confirming that our age signature was largely regulated in the "pro-longevity" direction. Quantitative network modeling demonstrated that age-regulated ncRNA equaled the contribution of protein-coding RNA within structures, but tended to have a lower heritability, implying lncRNA may better reflect environmental influences. Genes ECSIT, UNC13, and SKAP2 contributed to a network that did not respond to Rapamycin, and was associated with "neuron apoptotic processes" in protein-protein interaction analysis (FDR = 2.4%). ECSIT links inflammation with the continued age-related downwards trajectory of mitochondrial complex I gene expression (FDR < 0.01%), implying that sustained inhibition of ECSIT may be maladaptive. The present observations link, for the first time, model organism longevity programs with the endogenous but temporary genome-wide responses to aging in humans, revealing a pattern that may ultimately underpin personalized rates of health span.

Reference quantitative transcriptome dataset for adult Caenorhabditis elegans

Caenorhabditis elegans is a nematode widely used in biology and genomics as a model organism. We provide an integrated, quantitative reference map for the transcriptome of whole, wild type Bristol N2 strain C. elegans worms. The map has been obtained by meta-analysis of 110 gene expression profiles available in Gene Expression Omnibus (GEO) repository and integrated using the computational biology tool Transcriptome Mapper (TRAM). Following probe assignment to the relative locus and intra- and inter-sample normalization (in particular using the scaled quantile method), a mean, consensus reference value is provided for 45,932 transcripts, along with standard deviation. Expression values are all mapped in the context of genomic coordinates. The map provides easy access to relationships among expression values of different genes in this standard condition, highlights genomic segments with relatively high over-/under-expression and may serve as a reference to test for gene expression variation for both individual genes and the whole transcriptome in specific biological conditions (e.g. mutated strains or differently grown worms).