Loss of pex5 sensitizes zebrafish to fasting due to deregulated mitochondria, mTOR, and autophagy.

Animal models have been utilized to understand the pathogenesis of Zellweger spectrum disorders (ZSDs); however, the link between clinical manifestations and molecular pathways has not yet been clearly established. We generated peroxin 5 homozygous mutant zebrafish (pex5-/-) to gain insight into the molecular pathogenesis of peroxisome dysfunction. pex5-/- display hallmarks of ZSD in humans and die within one month after birth. Fasting rapidly depletes lipids and glycogen in pex5-/- livers and expedites their mortality. Mechanistically, deregulated mitochondria and mechanistic target of rapamycin (mTOR) signaling act together to induce metabolic alterations that deplete hepatic nutrients and accumulate damaged mitochondria. Accordingly, chemical interventions blocking either the mitochondrial function or mTOR complex 1 (mTORC1) or a combination of both improve the metabolic imbalance shown in the fasted pex5-/- livers and extend the survival of animals. In addition, the suppression of oxidative stress by N-acetyl L-cysteine (NAC) treatment rescued the apoptotic cell death and early mortality observed in pex5-/-. Furthermore, an autophagy activator effectively ameliorated the early mortality of fasted pex5-/-. These results suggest that fasting may be detrimental to patients with peroxisome dysfunction, and that modulating the mitochondria, mTORC1, autophagy activities, or oxidative stress may provide a therapeutic option to alleviate the symptoms of peroxisomal diseases associated with metabolic dysfunction.

Slc25a17 acts as a peroxisomal coenzyme A transporter and regulates multiorgan development in zebrafish.

Slc25a17 is known as a peroxisomal solute carrier, but the in vivo role of the protein has not been demonstrated. We found that the zebrafish genome contains two slc25a17 genes that function redundantly, but additively. Notably, peroxisome function in slc25a17 knockdown embryos is severely compromised, resulting in an altered lipid composition. Along the defects found in peroxisome-associated phenotypic presentations, we highlighted that development of the swim bladder is also highly dependent on Slc25a17 function. As Slc25a17 showed substrate specificity towards coenzyme A (CoA), injecting CoA, but not NAD+ , rescued the defective swim bladder induced by slc25a17 knockdown. These results indicated that Slc25a17 acts as a CoA transporter, involved in the maintenance of functional peroxisomes that are essential for the development of multiple organs during zebrafish embryogenesis. Given high homology in protein sequences, the role of zebrafish Slc25a17 may also be applicable to the mammalian system.

Loss of abcd4 in zebrafish leads to vitamin B12-deficiency anemia.

ABCD4, a member of the ATP-binding cassette transporter superfamily, is associated with the transport of vitamin B12 which is crucial for the development of red blood cells (RBCs) and may also be involved in its metabolism. However, the molecular function of ABCD4 during RBC development in zebrafish is mostly unknown. Using a morpholino-based knockdown approach, we found that abcd4-knockdown resulted in abnormal RBCs of irregular shapes and various sizes. o-Dianisidine staining, as an indicator of hemoglobin in RBCs, further confirmed that abcd4 morphants possessed fewer hemoglobinized cells and impaired blood circulation. Multiple protein sequence alignment revealed that the amino acid sequence for residues 13-292, which is the domain of vitamin B12 transport, of the zebrafish Abcd4 was highly conserved compared to that of other species. Accordingly, the abcd4 morphants can be rescued with human ABCD4, demonstrating a conserved role of ABCD4 in vertebrates. Notably, the vitamin B12-deficient phenotype in abcd4 morphants, which causes anemia, was recapitulated in the newly-established abcd4 mutant, indicating the possibility that the abcd4 mutant could be used as a disease model of vitamin B12-deficiency anemia. Our study provides an insight that the analysis of the newly-established abcd4 mutant may contribute to understanding its roles in ABCD4-related vitamin B12-deficiency anemia and the associated pathogeneses in humans.

Down-Regulation of EPAS1 Transcription and Genetic Adaptation of Tibetans to High-Altitude Hypoxia.

Tibetans are well adapted to the hypoxic environments at high altitude, yet the molecular mechanism of this adaptation remains elusive. We reported comprehensive genetic and functional analyses of EPAS1, a gene encoding hypoxia inducible factor 2α (HIF-2α) with the strongest signal of selection in previous genome-wide scans of Tibetans. We showed that the Tibetan-enriched EPAS1 variants down-regulate expression in human umbilical endothelial cells and placentas. Heterozygous EPAS1 knockout mice display blunted physiological responses to chronic hypoxia, mirroring the situation in Tibetans. Furthermore, we found that the Tibetan version of EPAS1 is not only associated with the relatively low hemoglobin level as a polycythemia protectant, but also is associated with a low pulmonary vasoconstriction response in Tibetans. We propose that the down-regulation of EPAS1 contributes to the molecular basis of Tibetans' adaption to high-altitude hypoxia.

PEX13 prevents pexophagy by regulating ubiquitinated PEX5 and peroxisomal ROS.

Peroxisomes are rapidly degraded during amino acid and oxygen deprivation by a type of selective autophagy called pexophagy. However, how damaged peroxisomes are detected and removed from the cell is poorly understood. Recent studies suggest that the peroxisomal matrix protein import machinery may serve double duty as a quality control machinery, where they are directly involved in activating pexophagy. Here, we explored whether any matrix import factors are required to prevent pexophagy, such that their loss designates peroxisomes for degradation. Using gene editing and quantitative fluorescence microscopy on culture cells and a zebrafish model system, we found that PEX13, a component of the peroxisomal matrix import system, is required to prevent the degradation of otherwise healthy peroxisomes. The loss of PEX13 caused an accumulation of ubiquitinated PEX5 on peroxisomes and an increase in peroxisome-dependent reactive oxygen species that coalesce to induce pexophagy. We also found that PEX13 protein level is downregulated to aid in the induction of pexophagy during amino acid starvation. Together, our study points to PEX13 as a novel pexophagy regulator that is modulated to maintain peroxisome homeostasis.Abbreviations: AAA ATPases: ATPases associated with diverse cellular activities; ABCD3: ATP binding cassette subfamily D member; 3ACOX1: acyl-CoA oxidase; 1ACTA1: actin alpha 1, skeletal muscle; ACTB: actin beta; ATG5: autophagy related 5; ATG7: autophagy related 7; ATG12: autophagy related 12; ATG16L1: autophagy related 16 like 1; CAT: catalase; CQ: chloroquine; Dpf: days post fertilization: FBS: fetal bovine serum; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; GFP: green fluorescent protein; H2O2: hydrogen peroxide; HA - human influenza hemagglutinin; HBSS: Hanks' Balanced Salt Solution; HCQ; hydroxychloroquine; KANL: lysine alanine asparagine leucine; KO: knockout; MAP1LC3B: microtubule associated protein 1 light chain 3 beta; MEF: mouse embryonic fibroblast; MTOR: mechanistic target of rapamycin kinase; MTORC1: mechanistic target of rapamycin kinase complex 1; MTORC2: mechanistic target of rapamycin kinase complex 2; MYC: MYC proto-oncogene, bHLH transcription factor; MZ: maternal and zygotic; NAC: N-acetyl cysteine; NBR1 - NBR1 autophagy cargo receptor; PBD: peroxisome biogenesis disorder; PBS: phosphate-buffered saline; PEX: peroxisomal biogenesis factor; PTS1: peroxisome targeting sequence 1; RFP: red fluorescent protein; ROS: reactive oxygen speciess; iRNA: short interfering RNA; SKL: serine lysine leucine; SLC25A17/PMP34: solute carrier family 25 member 17; Ub: ubiquitin; USP30: ubiquitin specific peptidase 30.

nudt7 gene depletion causes transcriptomic change in early development of zebrafish.

The Nudt family has been identified as enzymes performing Coenzyme A to 3'5'-ADP + 4'-phospho pantetheine catalysis. The members of this family have been shown to be particularly involved in lipid metabolism, while their involvement in gene regulation through regulating transcription or mRNA metabolism has also been suggested. Here, we focused on peroxisomal NUDT7, possessing enzymatic activity similar to that of its paralog, peroxisomal NUDT19, which is involved in mRNA degradation. No reports have been published about the Nudt family in zebrafish. Our transcriptomic data showed that the Nudt family members are highly expressed around zygotic gene activation (ZGA) in developing zebrafish embryos. Therefore, we confirmed the computational prediction that the products of the nudt7 gene in zebrafish were localized in the peroxisome and highly expressed in early embryogenesis. The depletion of nudt7 genes by the CRISPR/Cas9 system did not affect development; however, it decreased the rate of transcription in ZGA. In addition, H3K27ac ChIP-seq analysis demonstrated that this decrease in transcription was correlated with the genome-wide decrease of H3K27ac level. This study suggests that peroxisomal Nudt7 functions in regulating transcription in ZGA via formation of the H3K27ac domain in active chromatin.

Population history and genome wide association studies of birth weight in a native high altitude Ladakhi population.

Pathological low birth weight due to fetal growth restriction (FGR) is an important predictor of adverse obstetric and neonatal outcomes. It is more common amongst native lowlanders when gestating in the hypoxic environment of high altitude, whilst populations who have resided at high altitude for many generations are relatively protected. Genetic study of pregnant populations at high altitude permits exploration of the role of hypoxia in FGR pathogenesis, and perhaps of FGR pathogenesis more broadly. We studied the umbilical cord blood DNA of 316 neonates born to pregnant women managed at the Sonam Norboo Memorial Hospital, Ladakh (altitude 3540m) between February 2017 and January 2019. Principal component, admixture and genome wide association studies (GWAS) were applied to dense single nucleotide polymorphism (SNP) genetic data, to explore ancestry and genetic predictors of low birth weight. Our findings support Tibetan ancestry in the Ladakhi population, with subsequent admixture with neighboring Indo-Aryan populations. Fetal growth protection was evident in Ladakhi neonates. Although no variants achieved genome wide significance, we observed nominal association of seven variants across genes (ZBTB38, ZFP36L2, HMGA2, CDKAL1, PLCG1) previously associated with birthweight.

A prospective surveillance of drug prescribing and dispensing in a teaching hospital in western Nepal.

OBJECTIVE: To evaluate the drug dispensing practices and patients' knowledge on drug use among the outpatients and to identify and analyze the problems in drug prescribing and dispensing. METHODS: A prospective cross-sectional descriptive study was conducted using World Health Organization (WHO) core drug use indicators from July 13, 2008 to August 15, 2008 in Manipal Teaching Hospital, Pokhara, Nepal. RESULTS: A total of 4231 prescriptions were encountered with the total of 10591 drugs prescribed. The average number of drug per prescription was 2.5. Only 13% (n=10591) of drugs were prescribed by generic name. Percentage of drug prescribed from WHO model list of Essential drugs, Essential drug list of Nepal and Nepalese National Formulary was 21.7%, 32.8% and 42.3% respectively. Antibiotics and injections encountered were 28.3% and 3.1% respectively. Average cost per prescription was found to be Nepalese Rupees (NRs) 285.99 (US $ 3.73). Patient knowledge on correct use of drugs and appropriate labeling was found to be 81% and 1.4% respectively. Average dispensing time per prescription was 52 seconds. CONCLUSION: The finding from current study shows a trend towards irrational prescribing and dispensing. Hence, there is a need for effective intervention programme to encourage the physicians and dispensing pharmacists in promoting more rational drug use.

Population History and Altitude-Related Adaptation in the Sherpa.

The first ascent of Mount Everest by Tenzing Norgay and Sir Edmund Hillary in 1953 brought global attention to the Sherpa people and human performance at altitude. The Sherpa inhabit the Khumbu Valley of Nepal, and are descendants of a population that has resided continuously on the Tibetan plateau for the past ∼25,000 to 40,000 years. The long exposure of the Sherpa to an inhospitable environment has driven genetic selection and produced distinct adaptive phenotypes. This review summarizes the population history of the Sherpa and their physiological and genetic adaptation to hypoxia. Genomic studies have identified robust signals of positive selection across EPAS1, EGLN1, and PPARA, that are associated with hemoglobin levels, which likely protect the Sherpa from altitude sickness. However, the biological underpinnings of other adaptive phenotypes such as birth weight and the increased reproductive success of Sherpa women are unknown. Further studies are required to identify additional signatures of selection and refine existing Sherpa-specific adaptive phenotypes to understand how genetic factors have underpinned adaptation in this population. By correlating known and emerging signals of genetic selection with adaptive phenotypes, we can further reveal hypoxia-related biological mechanisms of adaptation. Ultimately this work could provide valuable information regarding treatments of hypoxia-related illnesses including stroke, heart failure, lung disease and cancer.

Sherpas share genetic variations with Tibetans for high-altitude adaptation.

BACKGROUND: Sherpas, a highlander population living in Khumbu region of Nepal, are well known for their superior climbing ability in Himalayas. However, the genetic basis of their adaptation to high-altitude environments remains elusive. METHODS: We collected DNA samples of 582 Sherpas from Nepal and Tibetan Autonomous Region of China, and we measured their hemoglobin levels and degrees of blood oxygen saturation. We genotyped 29 EPAS1 SNPs, two EGLN1 SNPs and the TED polymorphism (3.4 kb deletion) in Sherpas. We also performed genetic association analysis among these sequence variants with phenotypic data. RESULTS: We found similar allele frequencies on the tested 32 variants of these genes in Sherpas and Tibetans. Sherpa individuals carrying the derived alleles of EPAS1 (rs113305133, rs116611511 and rs12467821), EGLN1 (rs186996510 and rs12097901) and TED have lower hemoglobin levels when compared with those wild-type allele carriers. Most of the EPAS1 variants showing significant association with hemoglobin levels in Tibetans were replicated in Sherpas. CONCLUSION: The shared sequence variants and hemoglobin trait between Sherpas and Tibetans indicate a shared genetic basis for high-altitude adaptation, consistent with the proposal that Sherpas are in fact a recently derived population from Tibetans and they inherited adaptive variants for high-altitude adaptation from their Tibetan ancestors.

The fatty acid chain elongase, Elovl1, is required for kidney and swim bladder development during zebrafish embryogenesis.

Very long chain fatty acids are required for sphingolipid synthesis, lipid homeostasis, myelin formation, epidermal permeability, and retinal function. Seven different enzymes are known to be involved in the elongation cycle of fatty acids, with different chain-length specificities. Elovl1 is one of those enzymes whose function has been linked mainly to the synthesis of sphingolipids and the epidermal barrier. However, the role of Elovl1 in organogenesis is not clear. In zebrafish, 2 Elovl1 genes, elovl1a and elovl1b, are highly expressed in the swim bladder, and elovl1b is also expressed in the kidney. We found that both elovl1 knockdown embryos contain increased levels of long chain fatty acids from carbon number 14 to 20 as compared to control embryos. Oil-Red-O staining shows that yolk lipid consumption is greatly reduced, whereas lipid droplets accumulate within the swim bladder. Notably, knockdown of either elovl1a or elovl1b affects the expression of genes involved in swim bladder development and impairs inflation of the swim bladder. Consistent with its expression in the pronephros, knockdown of elovl1b alone affects the expression of genes required for kidney development and reduces renal clearance. Our findings strongly suggest that both elovl1 genes are a key determinant of swim bladder and kidney development in zebrafish, which may be comparatively applicable to lung and kidney development in humans.

Genetic evidence of a recent Tibetan ancestry to Sherpas in the Himalayan region.

Sherpas living around the Himalayas are renowned as high-altitude mountain climbers but when and where the Sherpa people originated from remains contentious. In this study, we collected DNA samples from 582 Sherpas living in Nepal and Tibet Autonomous Region of China to study the genetic diversity of both their maternal (mitochondrial DNA) and paternal (Y chromosome) lineages. Analysis showed that Sherpas share most of their paternal and maternal lineages with indigenous Tibetans, representing a recently derived sub-lineage. The estimated ages of two Sherpa-specific mtDNA sub-haplogroups (C4a3b1 and A15c1) indicate a shallow genetic divergence between Sherpas and Tibetans less than 1,500 years ago. These findings reject the previous theory that Sherpa and Han Chinese served as dual ancestral populations of Tibetans, and conversely suggest that Tibetans are the ancestral populations of the Sherpas, whose adaptive traits for high altitude were recently inherited from their ancestors in Tibet.

Cartilage development requires the function of Estrogen-related receptor alpha that directly regulates sox9 expression in zebrafish.

Estrogen-related receptor alpha (ESRRa) regulates a number of cellular processes including development of bone and muscles. However, direct evidence regarding its involvement in cartilage development remains elusive. In this report, we establish an in vivo role of Esrra in cartilage development during embryogenesis in zebrafish. Gene expression analysis indicates that esrra is expressed in developing pharyngeal arches where genes necessary for cartilage development are also expressed. Loss of function analysis shows that knockdown of esrra impairs expression of genes including sox9, col2a1, sox5, sox6, runx2 and col10a1 thus induces abnormally formed cartilage in pharyngeal arches. Importantly, we identify putative ESRRa binding elements in upstream regions of sox9 to which ESRRa can directly bind, indicating that Esrra may directly regulate sox9 expression. Accordingly, ectopic expression of sox9 rescues defective formation of cartilage induced by the knockdown of esrra. Taken together, our results indicate for the first time that ESRRa is essential for cartilage development by regulating sox9 expression during vertebrate development.

Developmental roles of D-bifunctional protein-A zebrafish model of peroxisome dysfunction.

The peroxisome is an intracellular organelle that responds dynamically to environmental changes. Various model organisms have been used to study the roles of peroxisomal proteins in maintaining cellular homeostasis. By taking advantage of the zebrafish model whose early stage of embryogenesis is dependent on yolk components, we examined the developmental roles of the D-bifunctional protein (Dbp), an essential enzyme in the peroxisomal ß-oxidation. The knockdown of dbp in zebrafish phenocopied clinical manifestations of its deficiency in human, including defective craniofacial morphogenesis, growth retardation, and abnormal neuronal development. Overexpression of murine Dbp rescued the morphological phenotypes induced by dbp knockdown, indicative of conserved roles of Dbp during zebrafish and mammalian development. Knockdown of dbp impaired normal development of blood, blood vessels, and most strikingly, endoderm-derived organs including the liver and pancreas - a phenotype not reported elsewhere in connection with peroxisome dysfunction. Taken together, our results demonstrate for the first time that zebrafish might be a useful model animal to study the role of peroxisomes during vertebrate development.