Evolutionary origins of retinoid active short-chain dehydrogenases/reductases of SDR16C family

17β-hydroxysteroid dehydrogenases in the progression of nonalcoholic fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) has become a worldwide epidemic and a major public health problem, with a prevalence of approximately 25%. The pathogenesis of NAFLD is complex and may be affected by the environment and susceptible genetic factors, resulting in a highly variable disease course and no approved drugs in the clinic. Notably, 17β-hydroxysteroid dehydrogenase type 13 (HSD17B13), which belongs to the 17β-hydroxysteroid dehydrogenase superfamily (HSD17Bs), is closely related to the clinical outcome of liver disease. HSD17Bs consists of fifteen members, most related to steroid and lipid metabolism, and may have the same biological function as HSD17B13. In this review, we highlight recent advances in basic research on the functional activities, major substrates, and key roles of HSD17Bs in the progression of NAFLD to develop innovative anti-NAFLD drugs targeting HSD17Bs.

Novel acetylation-related gene signatures for predicting the prognosis of patients with colorectal cancer

Histone acetylation may affect the tumorigenesis and prognosis of colorectal cancer (CRC). However, there is still a lack of studies exploring the effect of acetylation-related genes on the prognosis of CRC. To explore the role of acetylation-related genes in CRC prognosis using bioinformatics strategies, the expression data and survival information of CRC patients were collected from the Gene Expression Omnibus. The Molecular Signatures Database was used to select acetylation-related genes. Univariate and least absolute shrinkage and selection operator regression analyses were used to screen prognostic genes. Kaplan-Meier curves were plotted for survival analysis. Cibersort and pRRophetics were used to analyze immune infiltration and predict drug sensitivity, respectively. By implementing independent prognostic factors, a nomogram model was constructed. The result showed that a total of 48 prognostic genes which screened from the acetylation-related gene set were mainly enriched in ABC transporters and acetylation/deacetylation-related pathways. Three gene signatures (SDR16C5, MEAF6, and SOX4) were further defined, and a prognostic model was constructed that showed high sensitivity and specificity for predicting CRC prognosis in both training and validation cohorts. Patients with different prognostic risks also presented differential expression of gene signatures, infiltration of activated CD4 memory T cells, and drug sensitivity to bicalutamide, gefitinib, Lenalidomide, and imatinib. The nomogram suggested the potential of a risk score-based model in predicting 1- and 2-year survival in patients with CRC. In conclusion, we proposed three gene signatures from an acetylation-related gene set as potential targets for epigenetic therapy and constructed a prognostic model for CRC.

Mechanisms of Feedback Regulation of Vitamin A Metabolism

Vitamin A is an essential nutrient required throughout life. Through its various metabolites, vitamin A sustains fetal development, immunity, vision, and the maintenance, regulation, and repair of adult tissues. Abnormal tissue levels of the vitamin A metabolite, retinoic acid, can result in detrimental effects which can include congenital defects, immune deficiencies, proliferative defects, and toxicity. For this reason, intricate feedback mechanisms have evolved to allow tissues to generate appropriate levels of active retinoid metabolites despite variations in the level and format, or in the absorption and conversion efficiency of dietary vitamin A precursors. Here, we review basic mechanisms that govern vitamin A signaling and metabolism, and we focus on retinoic acid-controlled feedback mechanisms that contribute to vitamin A homeostasis. Several approaches to investigate mechanistic details of the vitamin A homeostatic regulation using genomic, gene editing, and chromatin capture technologies are also discussed.

The Protection Conferred by HSD17B13 rs72613567 Polymorphism on Risk of Steatohepatitis and Fibrosis May Be Limited to Selected Subgroups of Patients With NAFLD

INTRODUCTION

Our study aimed to explore how PNPLA3 rs738409 or phenotypic risk factors may moderate the relationship between HSD17B13 rs72613567 and risk of steatohepatitis and fibrosis.

METHODS

This analysis consisted of 1,153 non-Hispanic whites with biopsy-proven nonalcoholic fatty liver disease enrolled in the nonalcoholic steatohepatitis Clinical Research Network studies. Nonalcoholic fatty liver disease severity was determined by liver histology scored centrally according to the nonalcoholic steatohepatitis Clinical Research Network criteria. Moderation and logistic regression analyses were performed to identify the influence of moderators (PNPLA3 rs738409, age, sex, body mass index, and diabetes) on the relationship between HSD17B13 rs72613567 and risk of steatohepatitis and fibrosis.

RESULTS

HSD17B13 rs72613567 genotype frequency was as follows: (-/-), 64%; (-/A), 30%; (A/A), 6%. Moderation analysis showed that the protective effect of HSD17B13 rs72613567 A-allele on risk of steatohepatitis remained only significant among patients with PNPLA3 rs738409 genotype CC (β coeff: -0.19, P = 0.019), women (β coeff: -0.18, P < 0.001), patients of age ≥ 45 years (β coeff: -0.18, P < 0.001), patients with body mass index ≥ 35 kg/m2 (β coeff: -0.17, P < 0.001), and patients with diabetes (β coeff: -0.18, P = 0.020). Among women, the protective effect of HSD17B131 rs72613567 A-allele on risk of steatohepatitis was stronger in those aged ≥ 51 years. Logistic regression-based sensitivity analysis including various important subgroups confirmed our observations.

DISCUSSION

The protection conferred by HSD17B13 rs72613567 A-allele on risk of steatohepatitis and fibrosis may be limited to selected subgroups of individuals who are aged ≥ 45 years, women and have class ≥ 2 obesity or diabetes, and those with PNPLA3 rs738409 CC genotype.

17-Beta Hydroxysteroid Dehydrogenase 13 Deficiency Does Not Protect Mice From Obesogenic Diet Injury

BACKGROUND AND AIMS

17-Beta hydroxysteroid dehydrogenase 13 (HSD17B13) is genetically associated with human nonalcoholic fatty liver disease (NAFLD). Inactivating mutations in HSD17B13 protect humans from NAFLD-associated and alcohol-associated liver injury, fibrosis, cirrhosis, and hepatocellular carcinoma, leading to clinical trials of anti-HSD17B13 therapeutic agents in humans. We aimed to study the in vivo function of HSD17B13 using a mouse model.

APPROACH AND RESULTS

Single-cell RNA-sequencing and quantitative RT-PCR data revealed that hepatocytes are the main HSD17B13-expressing cells in mice and humans. We compared Hsd17b13 whole-body knockout (KO) mice and wild-type (WT) littermate controls fed regular chow (RC), a high-fat diet (HFD), a Western diet (WD), or the National Institute on Alcohol Abuse and Alcoholism model of alcohol exposure. HFD and WD induced significant weight gain, hepatic steatosis, and inflammation. However, there was no difference between genotypes with regard to body weight, liver weight, hepatic triglycerides (TG), histological inflammatory scores, expression of inflammation-related and fibrosis-related genes, and hepatic retinoid levels. Compared to WT, KO mice on the HFD had hepatic enrichment of most cholesterol esters, monoglycerides, and certain sphingolipid species. Extended feeding with the WD for 10 months led to extensive liver injury, fibrosis, and hepatocellular carcinoma, with no difference between genotypes. Under alcohol exposure, KO and WT mice showed similar hepatic TG and liver enzyme levels. Interestingly, chow-fed KO mice showed significantly higher body and liver weights compared to WT mice, while KO mice on obesogenic diets had a shift toward larger lipid droplets.

CONCLUSIONS

Extensive evaluation of Hsd17b13 deficiency in mice under several fatty liver-inducing dietary conditions did not reproduce the protective role of HSD17B13 loss-of-function mutants in human NAFLD. Moreover, mouse Hsd17b13 deficiency induces weight gain under RC. It is crucial to understand interspecies differences prior to leveraging HSD17B13 therapies.

Short-Chain Dehydrogenase NcmD Is Responsible for the C-10 Oxidation of Nocamycin F in Nocamycin Biosynthesis

Nocamycins I and II, featured with a tetramic acid scaffold, were isolated from the broth of Saccharothrix syringae NRRL B-16468. The biosynthesis of nocamycin I require an intermediate bearing a hydroxyl group at the C-10 position. A short chain dehydrogenase/reductase NcmD was proposed to catalyze the conversion of the hydroxyl group to ketone at the C-10 position. By using the λ-RED recombination technology, we generated the NcmD deletion mutant strain S. syringae MoS-1005, which produced a new intermediate nocamycin F with a hydroxyl group at C-10 position. We then overexpressed NcmD in Escherichia coli BL21 (DE3), purified the His₆-tagged protein NcmD to homogeneity and conducted in vitro enzymatic assays. NcmD showed preference to the cofactor NAD⁺, and it effectively catalyzed the conversion from nocamyin F to nocamycin G, harboring a ketone group at C-10 position. However, NcmD showed no catalytic activity toward nocamyin II. NcmD achieved maximum catalytic activity at 45°C and pH 8.5. The kinetics of NcmD toward nocamycin F was investigated at 45°C, pH 8.5 in the presence of 2 mM NAD⁺. The K_m and k_cat values were 131 ± 13 μM and 65 ± 5 min⁻¹, respectively. In this study, we have characterized NcmD as a dehydrogenase, which is involved in forming the ketone group at the C-10 position of nocamycin F. The results provide new insights to the nocamycin biosynthetic pathway.

Characterization of essential domains in HSD17B13 for cellular localization and enzymatic activity

Human genetic studies recently identified an association of SNPs in the 17-β hydroxysteroid dehydrogenase 13 (HSD17B13) gene with alcoholic and nonalcoholic fatty liver disease development. Mutant HSD17B13 variants devoid of enzymatic function have been demonstrated to be protective from cirrhosis and liver cancer, supporting the development of HSD17B13 as a promising therapeutic target. Previous studies have demonstrated that HSD17B13 is a lipid droplet (LD)-associated protein. However, the critical domains that drive LD targeting or determine the enzymatic activity have yet to be defined. Here we used mutagenesis to generate multiple truncated and point-mutated proteins and were able to demonstrate in vitro that the N-terminal hydrophobic domain, PAT-like domain, and a putative α-helix/β-sheet/α-helix domain in HSD17B13 are all critical for LD targeting. Similarly, we characterized the predicted catalytic, substrate-binding, and homodimer interaction sites and found them to be essential for the enzymatic activity of HSD17B13, in addition to our previous identification of amino acid P260 and cofactor binding site. In conclusion, we identified critical domains and amino acid sites that are essential for the LD localization and protein function of HSD17B13, which may facilitate understanding of its function and targeting of this protein to treat chronic liver diseases.

Generation and isolation of recombinant retinoid oxidoreductase complex

All-trans-retinoic acid (RA) is a bioactive lipid that influences many processes in embryonic and adult tissues. Given its bioactive nature, cellular concentrations of this molecule are highly regulated. The oxidation of all-trans-retinol to all-trans-retinaldehyde represents the first and rate-limiting step of the RA synthesis pathway. As such, it is the target of mechanisms that fine-tune RA levels within the cell. RDH10 is one enzyme responsible for the oxidation of all-trans-retinol to all-trans-retinaldehyde, and together with the all-trans-retinaldehyde reductase DHRS3 forms an oligomeric protein complex. The resulting retinoid oxidoreductase complex (ROC) is bifunctional and has the capacity to regulate steady-state levels of the direct precursor of RA, all-trans-retinaldehyde. As ROC represents a major regulatory element within the RA synthesis pathway, it is essential that methods are in place that allow for the study of this complex. Here we describe the production and isolation of recombinant ROC using a baculovirus expression system. Recombinant proteins retain enzymatic activities in intact microsomes and can be affinity purified for analysis. These methods can be used to assist in the assessment of ROC properties and the regulation of this protein complex's functional attributes.

Generation of Retinaldehyde for Retinoic Acid Biosynthesis

The concentration of all-trans-retinoic acid, the bioactive derivative of vitamin A, is critically important for the optimal performance of numerous physiological processes. Either too little or too much of retinoic acid in developing or adult tissues is equally harmful. All-trans-retinoic acid is produced by the irreversible oxidation of all-trans-retinaldehyde. Thus, the concentration of retinaldehyde as the immediate precursor of retinoic acid has to be tightly controlled. However, the enzymes that produce all-trans-retinaldehyde for retinoic acid biosynthesis and the mechanisms responsible for the control of retinaldehyde levels have not yet been fully defined. The goal of this review is to summarize the current state of knowledge regarding the identities of physiologically relevant retinol dehydrogenases, their enzymatic properties, and tissue distribution, and to discuss potential mechanisms for the regulation of the flux from retinol to retinaldehyde.

Mice lacking the epidermal retinol dehydrogenases SDR16C5 and SDR16C6 display accelerated hair growth and enlarged meibomian glands

Retinol dehydrogenases catalyze the rate-limiting step in the biosynthesis of retinoic acid, a bioactive lipid molecule that regulates the expression of hundreds of genes by binding to nuclear transcription factors, the retinoic acid receptors. Several enzymes exhibit retinol dehydrogenase activities in vitro; however, their physiological relevance for retinoic acid biosynthesis in vivo remains unclear. Here, we present evidence that two murine epidermal retinol dehydrogenases, short-chain dehydrogenase/reductase family 16C member 5 (SDR16C5) and SDR16C6, contribute to retinoic acid biosynthesis in living cells and are also essential for the oxidation of retinol to retinaldehyde in vivo Mice with targeted knockout of the more catalytically active SDR16C6 enzyme have no obvious phenotype, possibly due to functional redundancy, because Sdr16c5 and Sdr16c6 exhibit an overlapping expression pattern during later developmental stages and in adulthood. Mice that lack both enzymes are viable and fertile but display accelerated hair growth after shaving and also enlarged meibomian glands, consistent with a nearly 80% reduction in the retinol dehydrogenase activities of skin membrane fractions from the Sdr16c5/Sdr16c6 double-knockout mice. The up-regulation of hair-follicle stem cell genes is consistent with reduced retinoic acid signaling in the skin of the double-knockout mice. These results indicate that the retinol dehydrogenase activities of murine SDR16C5 and SDR16C6 enzymes are not critical for survival but are responsible for most of the retinol dehydrogenase activity in skin, essential for the regulation of the hair-follicle cycle, and required for the maintenance of both sebaceous and meibomian glands.

17-Beta Hydroxysteroid Dehydrogenase 13 Is a Hepatic Retinol Dehydrogenase Associated With Histological Features of Nonalcoholic Fatty Liver Disease

Nonalcoholic fatty liver disease (NAFLD) is a common cause of chronic liver disease. A single-nucleotide polymorphism (SNP), rs6834314, was associated with serum liver enzymes in the general population, presumably reflecting liver fat or injury. We studied rs6834314 and its nearest gene, 17-beta hydroxysteroid dehydrogenase 13 (HSD17B13), to identify associations with histological features of NAFLD and to characterize the functional role of HSD17B13 in NAFLD pathogenesis. The minor allele of rs6834314 was significantly associated with increased steatosis but decreased inflammation, ballooning, Mallory-Denk bodies, and liver enzyme levels in 768 adult Caucasians with biopsy-proven NAFLD and with cirrhosis in the general population. We found two plausible causative variants in the HSD17B13 gene. rs72613567, a splice-site SNP in high linkage with rs6834314 (r² = 0.94) generates splice variants and shows a similar pattern of association with NAFLD histology. Its minor allele generates simultaneous expression of exon 6-skipping and G-nucleotide insertion variants. Another SNP, rs62305723 (encoding a P260S mutation), is significantly associated with decreased ballooning and inflammation. Hepatic expression of HSD17B13 is 5.9-fold higher (P = 0.003) in patients with NAFLD. HSD17B13 is targeted to lipid droplets, requiring the conserved amino acid 22-28 sequence and amino acid 71-106 region. The protein has retinol dehydrogenase (RDH) activity, with enzymatic activity dependent on lipid droplet targeting and cofactor binding site. The exon 6 deletion, G insertion, and naturally occurring P260S mutation all confer loss of enzymatic activity. Conclusion: We demonstrate the association of variants in HSD17B13 with specific features of NAFLD histology and identify the enzyme as a lipid droplet-associated RDH; our data suggest that HSD17B13 plays a role in NAFLD through its enzymatic activity.

De novo draft assembly of the Botrylloides leachii genome provides further insight into tunicate evolution

Tunicates are marine invertebrates that compose the closest phylogenetic group to the vertebrates. These chordates present a particularly diverse range of regenerative abilities and life-history strategies. Consequently, tunicates provide an extraordinary perspective into the emergence and diversity of these traits. Here we describe the genome sequencing, annotation and analysis of the Stolidobranchian Botrylloides leachii. We have produced a high-quality 159 Mb assembly, 82% of the predicted 194 Mb genome. Analysing genome size, gene number, repetitive elements, orthologs clustering and gene ontology terms show that B. leachii has a genomic architecture similar to that of most solitary tunicates, while other recently sequenced colonial ascidians have undergone genome expansion. In addition, ortholog clustering has identified groups of candidate genes for the study of colonialism and whole-body regeneration. By analysing the structure and composition of conserved gene linkages, we observed examples of cluster breaks and gene dispersions, suggesting that several lineage-specific genome rearrangements occurred during tunicate evolution. We also found lineage-specific gene gain and loss within conserved cell-signalling pathways. Such examples of genetic changes within conserved cell-signalling pathways commonly associated with regeneration and development that may underlie some of the diverse regenerative abilities observed in tunicates. Overall, these results provide a novel resource for the study of tunicates and of colonial ascidians.

Characterization of human short chain dehydrogenase/reductase SDR16C family members related to retinol dehydrogenase 10

All-trans-retinoic acid (RA) is a bioactive derivative of vitamin A that serves as an activating ligand for nuclear transcription factors, retinoic acid receptors. RA biosynthesis is initiated by the enzymes that oxidize retinol to retinaldehyde. It is well established that retinol dehydrogenase 10 (RDH10, SDR16C4), which belongs to the 16C family of the short chain dehydrogenase/reductase (SDR) superfamily of proteins, is the major enzyme responsible for the oxidation of retinol to retinaldehyde for RA biosynthesis during embryogenesis. However, several lines of evidence point towards the existence of additional retinol dehydrogenases that contribute to RA biosynthesis in vivo. In close proximity to RDH10 gene on human chromosome 8 are located two genes that are phylogenetically related to RDH10. The predicted protein products of these genes, retinol dehydrogenase epidermal 2 (RDHE2, SDR16C5) and retinol dehydrogenase epidermal 2-similar (RDHE2S, SDR16C6), share 59% and 56% sequence similarity with RDH10, respectively. Previously, we showed that the single ortholog of the human RDHE2 and RDHE2S in frogs, Xenopus laevis rdhe2, oxidizes retinol to retinaldehyde and is essential for frog embryonic development. In this study, we explored the potential of each of the two human proteins to contribute to RA biosynthesis. The results of this study demonstrate that human RDHE2 exhibits a relatively low but reproducible activity when expressed in either HepG2 or HEK293 cells. Expression of the native RDHE2 is downregulated in the presence of elevated levels of RA. On the other hand, the protein encoded by the human RDHE2S gene is unstable when expressed in HEK293 cells. RDHE2S protein produced in Sf9 cells is stable but has no detectable catalytic activity towards retinol. We conclude that the human RDHE2S does not contribute to RA biosynthesis, whereas the low-activity RA-sensitive human RDHE2 may have a role in adjusting the cellular levels of RA in accord with specific physiological conditions.

Coelimination and Survival in Gene Network Evolution: Dismantling the RA-Signaling in a Chordate

The bloom of genomics is revealing gene loss as a pervasive evolutionary force generating genetic diversity that shapes the evolution of species. Outside bacteria and yeast, however, the understanding of the process of gene loss remains elusive, especially in the evolution of animal species. Here, using the dismantling of the retinoic acid metabolic gene network (RA-MGN) in the chordate Oikopleura dioica as a case study, we combine approaches of comparative genomics, phylogenetics, biochemistry, and developmental biology to investigate the mutational robustness associated to biased patterns of gene loss. We demonstrate the absence of alternative pathways for RA-synthesis in O. dioica, which suggests that gene losses of RA-MGN were not compensated by mutational robustness, but occurred in a scenario of regressive evolution. In addition, the lack of drastic phenotypic changes associated to the loss of RA-signaling provides an example of the inverse paradox of Evo-Devo. This work illustrates how the identification of patterns of gene coelimination-in our case five losses (Rdh10, Rdh16, Bco1, Aldh1a, and Cyp26)-is a useful strategy to recognize gene network modules associated to distinct functions. Our work also illustrates how the identification of survival genes helps to recognize neofunctionalization events and ancestral functions. Thus, the survival and extensive duplication of Cco and RdhE2 in O. dioica correlated with the acquisition of complex compartmentalization of expression domains in the digestive system and a process of enzymatic neofunctionalization of the Cco, while the surviving Aldh8 could be related to its ancestral housekeeping role against toxic aldehydes.